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 ...
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
- Page 46 and 47: de tal modo que todo punto x en el
- Page 48 and 49: N 1 (ξ) = (ξ − ξ 2) (ξ 1 −
- Page 50 and 51: modo de ubicar un nodo en todo punt
- Page 52 and 53: Figura 5 Malla de N nudos y N-1 ele
- Page 54 and 55: 3. Condiciones naturales de Neumann
- Page 56 and 57: (c) Considere las condiciones (ii)
- Page 58 and 59: los polinomios de interpolación in
- Page 60 and 61: La derivada segunda de u respecto a
- Page 62 and 63: 4.3.3. Formulación a partir del Pr
- Page 64 and 65: ¯K = ∫ L 0 BT D B d¯x 1 = ¯K =
- Page 66 and 67: Las cargas normales al eje de la vi
- Page 68 and 69: 2. θ 1 (ξ) = ∑ 2 I=1 N I (ξ)
- Page 70 and 71: 4.5.7. Cambio de base La expresión
- Page 72 and 73: Figura 5 Viga con deformación de c
- Page 74 and 75: 4.6.1. Matriz de rigidez de una vig
- Page 76 and 77: Donde las funciones de forma son: N
- Page 78 and 79: (una para cada W J ) en función de
- Page 80 and 81: Ejercicio 1-Sea el problema de conv
- Page 82 and 83: o directamente las coordenadas noda
- Page 84 and 85: donde r es el residuo que se quiere
- Page 86: ⎡ K 3−4 = ⎢ ⎣ K 4−5 = ⎡
- Page 89 and 90: Figura 1 Conducción del calor en 2
- Page 91 and 92: 5.3. Forma variacional del problema
- Page 93 and 94: 5.5. Flujo en un medio poroso El fl
- Page 95: ortotropía coinciden con las direc
- Page 99 and 100: Como φ es conocida sobre el contor
- Page 101 and 102: La formulación débil que resulta
- Page 103 and 104: 5.9.2. Formulación débil usando r
- Page 105 and 106: ⎡ σ = ⎢ ⎣ ⎤ σ 11 σ 22 σ
- Page 107 and 108: Figura 7 Teoría de placas clásica
- Page 109 and 110: Los esfuerzos de corte transversal
- Page 111 and 112: Capítulo 6 Elementos finitos en do
- Page 113 and 114: el punto p (ξ, η), y calculamos s
- Page 115 and 116: Para el caso cuadrático (ξ, η) =
- Page 117 and 118: ∂u ∂x 2 = NN∑ I=1 [ ∂N I
- Page 119 and 120: Figura 5 Mapeamientos para funcione
- Page 121 and 122: Elemento (L ( 1 , L 2 , L 3 ) w l L
- Page 123 and 124: ahora dependiente sólo del valor d
- Page 125 and 126: ∫ v ∫ σ ij δε ij dv = v δε
- Page 127 and 128: 6.8.1. Funciones de interpolación,
- Page 129 and 130: Siendo todas las variables constant
- Page 131 and 132: (el I por ejemplo) con dicha fuente
- Page 133 and 134: Capítulo 7 Aspectos generales asoc
- Page 135 and 136: 7.2. Imposición de restricciones n
- Page 137 and 138: entonces ¯k ij = k ij + 1 ǫ a ia
- Page 139 and 140: Figura 2 viga flexible entre column
- Page 141 and 142: D Matriz diagonal (definida positiv
- Page 143 and 144: 7.5. Suavizado de Variables En el m
- Page 145 and 146: 138
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