UNIVERSIDAD TECNOLOGICA NACIONAL CURSO S3001 & 3051
FACULTAD REGIONAL BUENOS AIRES CICLO 2018
Estabilidad 2
Vamos a estudiar las tensiones que aparecen en discos circulares de espesor constante que giran a altas velocidades alrededor de un eje normal a su plano. El estado de cargas sobre el disco que gira es el producido por las fuerzas centrífugas. Considerando que el espesor es pequeño y constante, que hay simetría geométrica y de cargas podemos asumir considerando coordenadas polares que el estado de tensión, el estado de deformación y los corrimientos solo son función del radio y no dependen de la coordenada y la tensión normal longitudinal es nula en los puntos en todo su espesor. Esto es rigurosamente cierto en las caras laterales, en el espesor la TME prueba su existencia pero si es pequeño (como es el caso) y varía anulándose en las caras podemos despreciarla. Así llegamos a la conclusión de que en todos los puntos existe un estado doble de tensiones con simetría axial.
Vamos a estudiar el equilibrio de un elemento de volumen como muestra la figura.
En el elemento de la figura la fuerza másica elemental será
dP = dma = dm 2r = dv 2r = hdrd 2r 2 = 2hr 2drd c g g g Dónde , g , son el peso específico del material del disco, la aceleración de la gravedad y la velocidad angular constante respectivamente. Como estos valores son constantes lo indicamos como 2 = k y entonces la fuerza centrífuga elemental resulta dP = khr 2 drd . g Estudiaremos las proyecciones sobre el eje r y sobre el eje normal al eje r, en el plano baricéntrico paralelo a las caras laterales. UTN – FRBA – Apunte Teórico – Discos Giratorios – rev 003 UNIVERSIDAD TECNOLOGICA NACIONAL CURSO S3001 & 3051
FACULTAD REGIONAL BUENOS AIRES CICLO 2018
Estabilidad 2
( r + d r ) hd ( r + dr ) − r hrd − 2 hdrsen d2 + khr 2drd = 0
d hrd + hd dr + hrd d + drd hd − hrd − 2 hdrsen + khr 2drd = 0 r r r r r 2 d d Despreciando infinitésimos de orden superior, siendo el sen = por ser d un ángulo 2 2 pequeño y simplificando, resulta dr + rd − dr + kr 2dr = 0 , si dividimos miembro a r r miembro por dr d Será finalmente − + r r = −kr 2 = d r r − (1) ( ) r dr dr
La proyección sobre la dirección es idénticamente nula.
Procedemos a poner a la ecuación (1) en función del corrimiento radial del punto (u), para ello calcularemos los alargamientos específicos en las direcciones radial y circunferencial. De la figura
d ( r +u ) − rd u El alargamiento específico circunferencial es = = ( 2) rd r u + du −(u − dr ) − dr du y el alargamiento específico radial = r = ( 3) dr dr Si el material se comporta según la ley de Hooke será 1 ( ) 1 ( ) = r − , = − r si multiplicamos a (3) por y la sumamos a (2) es r E E 1− 2 + = , repetimos lo anterior pero multiplicando (2) por y sumamos a (3) r E 1− 2 entonces + = r llamando a r E = E queda finalmente: = E + r 1 1 ( ) E 1− 2 (4)
y r 1 ( = E r + ) (5) u du du u y sustituyendo en (4) y (5) resulta: = E + y = E + 1 r dr r 1 dr r Si tomamos la expresión (1) y reemplazamos operando adecuadamente es: UTN – FRBA – Apunte Teórico – Discos Giratorios – rev 003 UNIVERSIDAD TECNOLOGICA NACIONAL CURSO S3001 & 3051
FACULTAD REGIONAL BUENOS AIRES CICLO 2018
Estabilidad 2
( d r r ) − = −kr 2 dr u d r 2u r −u u u + r u u 2 r + − = E1 + + rE1 − E1 +u =−kr dr r r 2 r r 2 r r u u 2u u u u u 2u u u kr 2 + +r + − − − =r + − =− r r 2r r r r r 2r r r E1
Para resolver esta ecuación diferencial debe obtenerse la solución de la ecuación
2u u u homogénea asociada r + − =0. 2r r r Se multiplica miembro a miembro por r y se reemplaza U = r m , se resuelve la ecuación de segundo grado m ( m −1) + m − 1 = 0 → m = 1 La solución general de dicha ecuación es u = C r + C r −1 GH 1 2 Utilizando el método de variación de parámetros se obtiene una solución particular kr 3 u =− que sumada a la anterior permite obtener la solución general de la ecuación p 8E1 diferencial a resolver. kr 3 u =− + C r + C r −1 . G 8E1 1 2 La expresión del corrimiento se reemplaza en las expresiones de las tensiones radiales y circunferenciales y se obtienen las expresiones siguientes: = E C1(1+ )+C2 (1− ) − kr 2 1+3 1 ( ) r2 8 r = E1 C1(1+ )−C2 (1− ) kr 2 − ( 3+ ) r2 8 Las constantes C y C se determinan por condiciones de borde. 1 2
Los discos pueden ser
• macizos donde el eje está incorporado • huecos o con agujero central, en los cuales el disco se monta sobre un eje de rotación.
UNIVERSIDAD TECNOLOGICA NACIONAL CURSO S3001 & 3051
FACULTAD REGIONAL BUENOS AIRES CICLO 2018
Estabilidad 2
DISCO MACIZO El disco constituye con el eje una sola pieza. Las tensiones deben ser finitas, pero en la expresión si calculamos la tensión circunferencial o radial para radio nulo resultarían infinitas. Esto nos indica que la constante C = 0 y la otra constante se calculan sabiendo 2 que la tensión radial se anula si el radio es el radio exterior del disco: kre2 kre2 ( 3+ ) r ( ) r = re = E C1(1+ ) − 1 8 ( 3+ ) = 0 → C1 = 8E1(1+ ) Expresiones finales k = ( 3 + ) re2 − (1 + 3 ) r 2 8 3+
( 3 + ) ( re2 − r 2 ) k = r 8 Ambas funciones son parabólicas y alcanzan su valor máximo en el centro del disco. Este valor máximo es el mismo para ambas tensiones y siempre es de tracción. k = = ( 3 + ) re2 mx rmx 8
En el radio exterior las tensiones radiales son nulas y la tensión circunferencial es
k ( r = re ) = ( 3 + ) re2 − (1 + 3 ) r 2 = k 1 − r 2 e ( )e 8 3+ 4
UNIVERSIDAD TECNOLOGICA NACIONAL CURSO S3001 & 3051
FACULTAD REGIONAL BUENOS AIRES CICLO 2018
Estabilidad 2
DISCO CON AGUJERO CENTRAL
En este caso el disco tiene un agujero central de radio interno r . i La determinación de las constantes se hace considerando que la tensión radial es nula tanto en el borde interno como el externo. (1− ) C − kre2 3+ r ( ) re = 0 = E (1+ )C1 − 1 2 2 8 ( ) re 1− ) kri 2 ( ( ) r ri =0= E1 (1+ )C1 − C − 2 2 8 ( 3+ ) ri Este sistema resuelto nos lleva a
C (1+ ) = 1 k 8E1 ( ( 3+ ) ri2 + re2 ) C (1− ) = 2 kri2re2 8E1 ( 3+ ) y a las expresiones finales de las tensiones radial y circunferencial r 2 r 2 (1+3 ) 2 k = ( 3+ ) ri2 + re2 + i e − r 8 r 2 ( 3 + ) ri2 re2 2 k r = ( 3+ ) ri + re − 2 2 −r 8 r 2 La tensión circunferencial alcanza su valor máximo para el radio interno y resulta:
k = ( 3+ ) re2 + ri2 (1− ) mx 4 ( 3+ ) Si el radio es el externo: (1− ) ( ) k r = re = ( 3+ ) ri2 + re2 4 ( 3+ ) La tensión radial será nula en los bordes y por lo tanto alcanzará un máximo cuya ubicación y valor calculamos derivando su expresión respecto del r e igualando a 0. r ri2re2 =0=2 − 2r = r 2r 2 − r 4 → r = ri re r r 3 i e
rmx ( ) k ( r = ri re = ( 3+ ) re − ri 8 2 ) Siendo los gráficos de ambas tensiones en función del radio
UNIVERSIDAD TECNOLOGICA NACIONAL CURSO S3001 & 3051
FACULTAD REGIONAL BUENOS AIRES CICLO 2018
Estabilidad 2
En ambos casos la tensión circunferencial crece hacia el centro o hasta el borde interno, con lo cual lo más racional es no mantener el espesor constante sino que este aumente desde el borde externo al centro o al borde interno, siendo esto más necesario cuando el número de revoluciones es alto como en las turbinas hidraúlicas.
UNIVERSIDAD TECNOLOGICA NACIONAL CURSO S3001 & 3051
FACULTAD REGIONAL BUENOS AIRES CICLO 2018
Estabilidad 2
DISCO MACIZO O CON AGUJERO CENTRAL DE ESPESOR VARIABLE
Considerando el análisis del elemento de volumen del disco, pero con un espesor t = f ( r ) en función del radio y planteando la ecuación de equilibrio se llega a:
( r + d r ) (t + dt )( r + dr ) d − r trd − 2tdr sen d2 + g 2r 2tdrd = 0
Simplificando, despreciando infinitésimos de orden superior, reemplazando el seno por el arco y dividiendo por r se llega a: d r 1 dt r − 2 + + + r = 0 (1) dr r t dr r g Esta ecuación diferencial tiene tres funciones incógnitas la tensión circunferencial, la tensión radial y el espesor, todas en función del radio. Dado que nos interesa tener el máximo aprovechamiento del material, plantearemos el caso que cumpla = = = cte r adm Para esta situación la (1) resulta 1 dt 2 dt 2 + r =0→ =− rdr = −k rdr adm t dr g t g adm 0 la ecuación diferencial tiene como solución general : dt r2 = −k rdr → n t = −k +C t 0 0 2 1
UNIVERSIDAD TECNOLOGICA NACIONAL CURSO S3001 & 3051
FACULTAD REGIONAL BUENOS AIRES CICLO 2018
Estabilidad 2
Para determinar la constante es necesario conocer el espesor en algún valor de r, por
r2 t r 2 − k 0 2 ejemplo en r = 0 es t = t → C = nt n = −k →t =t e 0 1 0 t0 0 2 0 Si el disco tuviera agujero central, suponiendo que en el radio exterior se conozca el espesor será:
re2 re2 − r 2 (k0 )re2 − r 2 t r = r es t = t → C = nt + k n =k →t =t e 2 e e 1 e 0 2 te 0 2 e Como vemos en ambos casos el espesor crece exponencialmente desde el borde exterior hacia el centro. Este sólido se lo llama de igual resistencia porque en todos los puntos su tensión es la admisible.
