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 ...
2. θ 1 (ξ) = ∑ 2 I=1 N I (ξ) θ I 1 3. u 2 (ξ) = φ 1 (ξ) u 1 2 + ϕ 1 (ξ) θ 1 3 + φ 2 (ξ) u 2 2 + ϕ 2 (ξ) θ 2 3 4. u 3 (ξ) = φ 1 (ξ) u 1 3 − ϕ1 (ξ) θ 1 2 + φ2 (ξ) u 2 3 − ϕ2 (ξ) θ 2 2 Donde puede verse que hemos usado para el comportamiento axial los polinomios de Lagrange de grado 1 y para el comportamiento flexional los polinomios de Hermite de orden 1. 4.5.5. Desarrollo Recordando que d = d dξ = d 2 dx 1 dξ dx 1 dξ L ( ) −1 dξ dx1 = = J −1 dx 1 dξ tenemos que las deformaciones generalizadas resultan ε = du 1 = N I 2 dx ′ ξ 1 L uI 1 = ( ) u 2 1 − u 1 1 1 L χ 1 = dθ 1 = N I 2 dx ′ ξ 1 L θI 1 = ( θ 2 1 − ) 1 θ1 1 L χ 3 = d2 u 2 dx 2 1 = [ 6ξ L ( ) = d2 u 2 dξ 2 dξ 2 = 4 [ ] φ 1 1 dx 1 L 2 ′ ξξu 1 2 + ϕ 1′ ξξθ 1 3 + φ 2′ ξξu 2 2 + ϕ 2′ ξξθ 2 3 ⎡ ⎤ −1 + 3ξ , , − 6ξ 2 L L , 1 + 3ξ ] u 1 2 ⎢ θ 1 3 ⎥ 2 L ⎣ u 2 ⎦ 2 θ 2 3 χ 2 = − d2 u 3 dx 2 1 [ 6ξ = − L ( ) 2 = − d2 u 3 dξ dξ 2 = − 4 [ φ 1 1 dx 1 L 2 ′ ξξ u1 3 − ϕ1′ ξξ θ1 2 + φ2′ ξξ u2 3 − ] ϕ2′ ξξ θ2 2 ⎡ ⎤ + 3ξ , −−1 , − 6ξ ] u 1 3 + 3ξ , −1 ⎢ θ 1 2 ⎥ 2 L L2 L ⎣ u 2 ⎦ 3 θ 2 2 ⎡ ⎢ ⎣ Que puede resumirse en la siguiente expresión matricial ⎤ ε χ 1 χ 3 χ 2 ⎥ ⎦ = 1 L ⎡ ⎢ ⎣ −1 +1 −1 +1 6ξ −1 + 3ξ − 6ξ L L −1 + 3ξ − 6ξ L 6ξ L +1 + 3ξ +1 + 3ξ ⎡ ⎤ ⎥ ⎦ ⎢ ⎣ u 1 1 u 1 2 u 1 3 θ 1 1 θ 1 2 θ 1 3 u 2 1 u 2 2 u 2 3 θ 2 1 θ 2 2 θ 2 3 ⎤ ⎥ ⎦ 70
}{{} ε =B (ξ) } {{ } }{{} a e 4×1 4×12 12×1 En tanto que la relaciones constitutivas pueden escribirse matricialmente como ⎡ ⎤ ⎡ ⎤ ⎡ ⎤ N EA ε ⎢ M t ⎥ ⎣ M 3 ⎦ = ⎢ αGI p ⎥ ⎢ χ 1 ⎥ ⎣ EI 3 ⎦ ⎣ χ 3 ⎦ M 2 EI 2 χ 2 σ }{{} 4×1 = D }{{} 4×4 ε }{{} 4×1 Finalmente la matriz de rigidez (en el sistema local) resulta de integrar a lo largo de la viga K }{{} L = 12×12 ∫ L 0 B T }{{} 12×4 D }{{} 4×4 }{{} B dx 1 = 4×12 ∫ 1 −1 B (ξ) T D B (ξ) L 2 dξ Notar que como máximo en B(ξ) hay polinomios de orden 1, luego en el producto B T D B hay como máximo polinomios de orden 2 en ξ, por lo cual bastan dos puntos de integración si se va a realizar una integración numérica. La integral puede hacerse en este caso en forma analítica sin ningún problema. La matriz de rigidez (cuya obtención se deja como ejercicio) es idéntica a la obtenida en los cursos de cálculo matricial de estructuras. Esto es así porque las funciones de interpolación utilizadas son capaces de reproducir la solución exacta de las ecuaciones diferenciales (homogéneas). 4.5.6. Término independiente (vector de cargas) Supongamos que las cargas están representadas localmente (es decir respecto al sistema local) y varían en forma lineal dentro del elemento ⎡ q (ξ) = ⎣ q ⎤ 1 2∑ q 2 ⎦ = N I (ξ) q I q 3 I=1 El trabajo virtual externo se puede escribir como ⎡ ∫ L ∫ L T.V.E. = δu T q (ξ) dx 1 = [δu 1 , δu 2 , δu 3 ] ⎣ q ⎤ 1 (ξ) q 2 (ξ) ⎦ dx 1 0 0 q 3 (ξ) ⎡ ⎤ ∫ L N 1 N 2 T = [δa e ] T ⎣ φ 1 ϕ 1 φ 2 ϕ 2 ⎦ 0 φ 1 0 −ϕ 1 φ 2 0 −ϕ 2 Donde F L × ⎡ = [δa e ] T } {{ } 1×12 ⎣ N 1 N 2 N 1 N 2 N 1 ∫ L 0 Φ T (ξ) } {{ } 12×3 N (ξ) } {{ } 3×6 dx 1 [ q 1 q 2 ] } {{ } 6×1 ⎡ ⎤ ⎦ dx 1 N 2 ⎢ ⎣ q 1 1 q 1 2 q 1 3 q 2 1 q 2 2 q 2 3 ⎤ ⎥ ⎦ = [δa e ] T F } {{ } }{{} L 1×12 12×1 es el término independiente referido al sistema de coordenadas locales. (4.3) 71
- Page 17 and 18: La sustitución de la aproximación
- Page 19 and 20: Figura 2 Ejercicio 21. a) Problema
- Page 21 and 22: Capítulo 2 El método de elementos
- Page 23 and 24: Figura 1 Aproximación de una funci
- Page 25 and 26: Figura 2 Comportamiento de tres fun
- Page 27 and 28: 2.4.1. Propiedades de la matriz K y
- Page 29 and 30: Figura 5 Descripción local y globa
- Page 31 and 32: 2.4.3.2. Cálculo del vector de car
- Page 33 and 34: de funciones. Una medida universalm
- Page 35: Figura 6 Curvas del error utilizand
- Page 38 and 39: eferiremos como ” ley de conserva
- Page 40 and 41: ambos miembros haciendo tender el p
- Page 42 and 43: 1. Condiciones de borde esenciales,
- Page 44 and 45: 2. Las condiciones de borde no pued
- 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 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 and 96: ortotropía coinciden con las direc
- Page 97 and 98: 5.6.5. Forma débil de la ecuación
- 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
2. θ 1 (ξ) = ∑ 2<br />
I=1 N I (ξ) θ I 1<br />
3. u 2 (ξ) = φ 1 (ξ) u 1 2 + ϕ 1 (ξ) θ 1 3 + φ 2 (ξ) u 2 2 + ϕ 2 (ξ) θ 2 3<br />
4. u 3 (ξ) = φ 1 (ξ) u 1 3 − ϕ1 (ξ) θ 1 2 + φ2 (ξ) u 2 3 − ϕ2 (ξ) θ 2 2<br />
Don<strong>de</strong> pue<strong>de</strong> verse que hemos usado para el comportamiento axial los polinomios <strong>de</strong> Lagrange<br />
<strong>de</strong> grado 1 y para el comportamiento flexional los polinomios <strong>de</strong> Hermite <strong>de</strong> or<strong>de</strong>n 1.<br />
4.5.5. Desarrollo<br />
Recordando que<br />
d<br />
= d dξ<br />
= d 2<br />
dx 1 dξ dx 1 dξ L<br />
( ) −1<br />
dξ dx1<br />
=<br />
= J −1<br />
dx 1 dξ<br />
tenemos que las <strong>de</strong>formaciones generalizadas resultan<br />
ε = du 1<br />
= N I 2<br />
dx ′ ξ<br />
1 L uI 1 = ( )<br />
u 2 1 − u 1 1<br />
1<br />
L<br />
χ 1 = dθ 1<br />
= N I 2<br />
dx ′ ξ<br />
1 L θI 1 = ( θ 2 1 − ) 1<br />
θ1 1<br />
L<br />
χ 3 = d2 u 2<br />
dx 2 1<br />
=<br />
[ 6ξ<br />
L<br />
( )<br />
= d2 u 2 dξ 2<br />
dξ 2 = 4 [ ]<br />
φ<br />
1<br />
1<br />
dx 1 L 2 ′ ξξu 1 2 + ϕ 1′ ξξθ 1 3 + φ 2′ ξξu 2 2 + ϕ 2′ ξξθ 2 3<br />
⎡ ⎤<br />
−1 + 3ξ<br />
, , − 6ξ<br />
2<br />
L L , 1 + 3ξ ]<br />
u 1 2<br />
⎢ θ 1 3 ⎥<br />
2 L ⎣ u 2 ⎦<br />
2<br />
θ 2 3<br />
χ 2 = − d2 u 3<br />
dx 2 1<br />
[ 6ξ<br />
= −<br />
L<br />
( ) 2<br />
= − d2 u 3 dξ<br />
dξ 2 = − 4 [<br />
φ<br />
1<br />
1<br />
dx 1 L 2 ′ ξξ u1 3 − ϕ1′ ξξ θ1 2 + φ2′ ξξ u2 3 − ]<br />
ϕ2′ ξξ θ2 2<br />
⎡ ⎤<br />
+ 3ξ<br />
, −−1 , − 6ξ ]<br />
u 1 3<br />
+ 3ξ<br />
, −1 ⎢ θ 1 2 ⎥<br />
2<br />
L L2 L ⎣ u 2 ⎦<br />
3<br />
θ 2 2<br />
⎡<br />
⎢<br />
⎣<br />
Que pue<strong>de</strong> resumirse en la siguiente expresión matricial<br />
⎤<br />
ε<br />
χ 1<br />
χ 3<br />
χ 2<br />
⎥<br />
⎦ = 1 L<br />
⎡<br />
⎢<br />
⎣<br />
−1 +1<br />
−1 +1<br />
6ξ<br />
−1 + 3ξ − 6ξ<br />
L L<br />
−1 + 3ξ<br />
− 6ξ<br />
L<br />
6ξ<br />
L<br />
+1 + 3ξ<br />
+1 + 3ξ<br />
⎡<br />
⎤<br />
⎥<br />
⎦<br />
⎢<br />
⎣<br />
u 1 1<br />
u 1 2<br />
u 1 3<br />
θ 1 1<br />
θ 1 2<br />
θ 1 3<br />
u 2 1<br />
u 2 2<br />
u 2 3<br />
θ 2 1<br />
θ 2 2<br />
θ 2 3<br />
⎤<br />
⎥<br />
⎦<br />
70