Capítulo 1 Métodos de residuos ponderados Funciones de prueba ...

Capítulo 1 Métodos de residuos ponderados Funciones de prueba ... Capítulo 1 Métodos de residuos ponderados Funciones de prueba ...

efn.uncor.edu
from efn.uncor.edu More from this publisher
15.07.2014 Views

La ecuación de equilibrio en la dirección z es la que gobierna el alabeo ∂σ xz ∂x + ∂σ yz ∂y + ∂σ zz ∂z + F z = ∂ ∂x [G xz (−βy + w′ x)] + ∂ ∂y [G yz (βx + w′ y)] = 0 (5.2) G ∂w′ x xz ∂x + G ∂w′ y yz ∂y ∂ 2 w = G xz ∂x + G ∂ 2 w 2 yz ∂y = 0 (5.3) 2 Donde en la segunda expresión se ha supuesto que el material es homogéneo Si el material además es isótropo G xz = G yz = G se obtiene la ecuación de Laplace 5.6.4. Condiciones de contorno ( ) ∂ 2 w G ∂x + ∂2 w = G ∇ · ∇w = 0 2 ∂y 2 Figura 5 Condición de contorno en la torsión Las condiciones de contorno son exclusivamente naturales (Neumann). La tensión de corte normal al contorno debe ser cero τ ν = 0. Denominando con α al ángulo que forma la normal al contorno ν con el eje x, tenemos que (llamando s a lo longitud de arco sobre el contorno ν x = cos α = dy ds ν y = sin α = − dx ds la tensión de corte normal es reemplazando las expresiones 5.1a Si el material es isótropo dy τ ν = τ · ν = τ xz ds − τ dx yz ds τ ν = [G xz (−βy + w′ x)] dy ds − [G yz (βx + w′ y)] dx [ ds = −β G xz y dy ds + G yzx dx ] dy + G xz ds w′ x ds − G dx yzw′ y ds [ τ ν = −Gβ y dy ] ( ds + xdx + G ds } {{ } 1 dr 2 2 ds w′ x ) dy ds − dx w′ y ds (5.4) (5.5) 91

5.6.5. Forma débil de la ecuación de alabeo Aplicando residuos ponderados sobre la expresión 5.2, con v la función de ponderación ∫ { ∂ v A ∂x [G xz (−βy + w′ x)] + ∂ } ∂y [G yz (βx + w′ y)] dA = 0 ∫ [ ∂ v ∂x (Gxzw′ x) + ∂ ] [ ∂y (Gyzw′ y) + vβ − ∂ ∂x (G xzy) + ∂ ] ∂y (G yzx) dA = 0 A Notar que el segundo término es necesario sólo cuando el material no es homogéneo, es decir cuando hay una variación de la matriz constitutiva del material. Esta formulación permite entonces tratar el alabeo de secciones compuestas de distintos materiales. Notar además que el segundo término incluye sólo valores conocidos, por lo cual podríamos separar en dos miembros la ecuación, de tal forma que el segundo miembro es nulo para materiales homogéneos ∫ A v [ ∂ ∂x (G xzw′ x) + ∂ ∂y (G yzw′ y) ] ∫ dA = A vβ [ ∂ ∂x (G xzy) − ∂ ∂y (G yzx) ] dA Integrando por partes ambos miembros ∫ [ ∂ v A ∂x (Gxzw′ x) + ∂ ] ∫ ∫ ∂y (Gyzw′ y) dA = v [G xz w′ xν x + G yz w′ yν y ] ds − [G xz w′ xv′ x + G yz w′ yv′ y] dA S A ∫ [ ∂ vβ ∂x (G xzy) − ∂ ] ∫ ∫ ∂y (G yzx) dA = vβ [G xz yν x − G yz xν y ] ds − β [G xz yv′ x − G yz xv′ y] dA A S Reemplazando en la expresión anterior y notando que los términos en el contorno se anulan entre si (ver ecuación 5.5) ∫ ∫ v {G xz w′ xν x + G yz w′ yν y − β [G xz yν x − G yz xν y ]} ds = vτ ν ds = 0 resulta ∫ S ∫ [G xz w′ xv′ x + G yz w′ yv′ y] dA = A [ ] [ ] Gxz w′ x [v′ x, v′ y] dA A G yz y = ∫ ∫ A A β [G xz yv′ x − G yz xv′ y] dA [ ] [ Gxz y β [v′ x, v′ y] G yz −x s A ] dA Notar que para la solución del problema es necesario fijar en algún punto el valor de w, además habitualmente se resuelve el problema para un valor unitario de β. En secciones simétricas es suficiente con discretizar una de las porciones simétricas, en este caso deben imponerse condiciones de contorno sobre w (esenciales), w = 0 en las líneas de simetría. La formulación presentada hasta aquí sigue los lineamientos del método de los desplazamientos (o método de rigidez) consistente en resolver las ecuaciones de equilibrio del problema expresadas en función de las incógnitas de desplazamiento. Una vez obtenida la solución del problema la determinación de deformaciones y tensiones es directa. A continuación veremos una solución alternativa, consistente en resolver ecuaciones de compatibilidad, lo que se asocia habitualmente con el método de las fuerzas, contraparte de esta formulación 5.6.6. Función de tensión La formulación más renombrada para el análisis del problema de alabeo está asociada a la solución de una ecuación de compatibilidad. Esto es así porque, como veremos, para materiales isótropos y homogéneos , resulta como ecuación de gobierno la ecuación de Laplace con condiciones de contorno sólo esenciales y homogéneas (secciones simplemente conexas). Esto permite hacer analogías con otros problemas mecánicos, como por ejemplo la membrana traccionada que se describió antes (“analogía de la membrana” debida a Prandtl). 92

La ecuación <strong>de</strong> equilibrio en la dirección z es la que gobierna el alabeo<br />

∂σ xz<br />

∂x + ∂σ yz<br />

∂y<br />

+ ∂σ zz<br />

∂z + F z = ∂<br />

∂x [G xz (−βy + w′ x)] + ∂ ∂y [G yz (βx + w′ y)] = 0 (5.2)<br />

G ∂w′ x<br />

xz<br />

∂x + G ∂w′ y<br />

yz<br />

∂y<br />

∂ 2 w<br />

= G xz<br />

∂x + G ∂ 2 w<br />

2 yz<br />

∂y = 0 (5.3)<br />

2<br />

Don<strong>de</strong> en la segunda expresión se ha supuesto que el material es homogéneo<br />

Si el material a<strong>de</strong>más es isótropo G xz = G yz = G se obtiene la ecuación <strong>de</strong> Laplace<br />

5.6.4. Condiciones <strong>de</strong> contorno<br />

( )<br />

∂ 2 w<br />

G<br />

∂x + ∂2 w<br />

= G ∇ · ∇w = 0<br />

2 ∂y 2<br />

Figura 5<br />

Condición <strong>de</strong> contorno en la torsión<br />

Las condiciones <strong>de</strong> contorno son exclusivamente naturales (Neumann). La tensión <strong>de</strong> corte<br />

normal al contorno <strong>de</strong>be ser cero τ ν = 0. Denominando con α al ángulo que forma la normal al<br />

contorno ν con el eje x, tenemos que (llamando s a lo longitud <strong>de</strong> arco sobre el contorno<br />

ν x = cos α = dy<br />

ds<br />

ν y = sin α = − dx<br />

ds<br />

la tensión <strong>de</strong> corte normal es<br />

reemplazando las expresiones 5.1a<br />

Si el material es isótropo<br />

dy<br />

τ ν = τ · ν = τ xz<br />

ds − τ dx<br />

yz<br />

ds<br />

τ ν = [G xz (−βy + w′ x)] dy<br />

ds − [G yz (βx + w′ y)] dx<br />

[<br />

ds<br />

= −β G xz y dy<br />

ds + G yzx dx ]<br />

dy<br />

+ G xz<br />

ds<br />

w′ x<br />

ds − G dx<br />

yzw′ y<br />

ds<br />

[<br />

τ ν = −Gβ y dy ] (<br />

ds + xdx + G<br />

ds<br />

} {{ }<br />

1 dr 2<br />

2 ds<br />

w′ x<br />

)<br />

dy<br />

ds − dx<br />

w′ y<br />

ds<br />

(5.4)<br />

(5.5)<br />

91

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!