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 ...
Figura 6 Elasticidad tridimensional Las ecuaciones anteriores desarrolladas para el problema tridimensional resultan: ⎡ ∂σ 11 + ∂σ 12 + ∂σ ⎤ 13 ∂x 1 ∂x 2 ∂x 3 ∂σ 21 ∇ · σ = + ∂σ 22 + ∂σ 23 = ∂σ ij t i = ∂σ ji t i ∂x 1 ∂x 2 ∂x 3 ∂x j ∂x j ⎢ ⎣ ∂σ 31 + ∂σ 32 + ∂σ ⎥ ⎦ 33 ∂x 3 ∂x 2 ∂x 3 ⎡ σ · n = σ T n = ⎢ ⎣ ∂σ ij ∂x j + ρ (b i − a i ) = 0 ⎤ σ 11 ν 1 + σ 12 ν 2 + σ 13 ν 3 σ 21 ν 1 + σ 22 ν 2 + σ 23 ν 3 ⎥ ⎦ = σ ijν j t i = σ ij ν i t j σ 31 ν 1 + σ 32 ν 2 + σ 33 ν 3 σ ij ν j = f i (∇u) ij = ∂u i ∂x j ε ij = (∇ s u) ij = 1 2 ( ∂uj + ∂u ) i ∂x i ∂x j Para un material isótropo, el tensor de elasticidad depende de sólo dos constantes y puede escribirse D =2µ I+λ 1 ⊗ 1 D ijkl = µ (δ ik δ jl + δ il δ jk ) + λ δ ij δ kl E Eν µ = λ = 2 (1 + ν) (1 + ν) (1 − 2ν) donde µ y λ son los parámetros de Lamé, E es el módulo de elasticidad de Young y ν es la relación de Poisson. I es el tensor identidad de cuarto orden, 1 es el tensor identidad de segundo orden y ⊗ denota el producto tensorial σ ij = D ijkl ε kl = [µ (δ ik δ jl + δ il δ jk ) + λ δ ij δ kl ] ε kl σ ij = 2µ ε ij + δ ij λ ε kk 97
5.9.2. Formulación débil usando residuos ponderados Apliquemos el método de residuos ponderados a la ecuación de equilibrio (5.10) con una función de ponderación w = (w 1 , w 2 , w 3 ) donde los w i son funciones independientes una de otra y ponderan cada componente de la ecuación de equilibrio ∫ Ω w· {∇ · σ + ρ (x) [b (x) − a (x)]} dΩ = 0 ∫ ∫ w· (∇ · σ) dΩ = − Ω Integremos por partes el primer miembro, para lo cual recordemos que: Ω ρ (x) w· [b (x) − a (x)] dΩ w · σ = w T σ = w i σ ij t j ∇ · (w · σ) = ∂ ∂x j (w i σ ij ) = ∂w i ∂σ ij σ ij + w i ∂x j ∂x j = ∇w : σ + w· (∇ · σ) la última expresión permite escribir el primer miembro del residuo como ∫ Ω ∫ w· (∇ · σ) dΩ = ∂Ω ∫ w·(σ · n) d∂Ω − ∇w : σ dΩ } {{ } Ω f(s) donde el operador “:” es el producto punto entre tensores de segundo orden, similar al de vectores (por ej.: σ : ε = σ ij ε ij = ∑ 3 ∑ 3 i=1 j=1 σ ij ε ij ). Notando además que debido a la simetría del tensor de tensiones ∇w : σ = ∇ s w : σ la integral ponderada del residuo puede escribirse: ∫ Ω ∫ ∇ s w : σ dΩ = Reemplazando la ecuación constitutiva tenemos: ∫ Ω ∫ ∇ s w : D : ∇ s u dΩ = Ω ∫ w·ρ (x) [b (x) − a (x)] dΩ + Ω ∂Ω ∫ w·ρ (x) [b (x) − a (x)] dΩ + w · f (s) d∂Ω ∂Ω w · f (s) d∂Ω Las condiciones sobre la solución u son: continuidad (compatibilidad), derivabilidad (∇u debe existir y poder ser calculado) y u = ū en ∂Ω u . Al usar Galerkin e integrar por partes, las condiciones sobre w resultan similares: continuidad y derivabilidad (∇w debe existir y poder ser calculado) y w = 0 en ∂Ω u . Esta última condición permite dividir la segunda integral del segundo miembro en dos partes, dividiendo el contorno en dos partes (∂Ω σ y ∂Ω u ) en la segunda parte la integral resulta entonces identicamente nula. Finalmente si llamamos ∫ Ω ∫ ¯ε ij D ijkl ε kl dΩ = Ω ¯ε = ∇ s w ∫ w i ρ (x) [b i (x) − a i (x)] dΩ + w i f i (s) d∂Ω s ∂Ω σ donde en el primer miembro se puede reemplazar el tensor D 98 ∫ Ω ∫ ¯ε ij D ijkl ε kl dΩ = Ω (2µ ¯ε ij ε ij + λ ¯ε kk ε ll ) dΩ
- 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 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: La formulación débil que resulta
- 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
- Page 147 and 148: Difícilmente se pueda encontrar un
- Page 149 and 150: que para el caso de placas esbeltas
- Page 151 and 152: cualquier variación de deformacion
Figura 6<br />
Elasticidad tridimensional<br />
Las ecuaciones anteriores <strong>de</strong>sarrolladas para el problema tridimensional resultan:<br />
⎡<br />
∂σ 11<br />
+ ∂σ 12<br />
+ ∂σ ⎤<br />
13<br />
∂x 1 ∂x 2 ∂x 3<br />
∂σ 21<br />
∇ · σ =<br />
+ ∂σ 22<br />
+ ∂σ 23<br />
= ∂σ ij<br />
t i = ∂σ ji<br />
t i<br />
∂x 1 ∂x 2 ∂x 3<br />
∂x j ∂x j<br />
⎢<br />
⎣ ∂σ 31<br />
+ ∂σ 32<br />
+ ∂σ ⎥<br />
⎦<br />
33<br />
∂x 3 ∂x 2 ∂x 3<br />
⎡<br />
σ · n = σ T n =<br />
⎢<br />
⎣<br />
∂σ ij<br />
∂x j<br />
+ ρ (b i − a i ) = 0<br />
⎤<br />
σ 11 ν 1 + σ 12 ν 2 + σ 13 ν 3<br />
σ 21 ν 1 + σ 22 ν 2 + σ 23 ν 3<br />
⎥<br />
⎦ = σ ijν j t i = σ ij ν i t j<br />
σ 31 ν 1 + σ 32 ν 2 + σ 33 ν 3<br />
σ ij ν j = f i<br />
(∇u) ij<br />
= ∂u i<br />
∂x j<br />
ε ij = (∇ s u) ij<br />
= 1 2<br />
( ∂uj<br />
+ ∂u )<br />
i<br />
∂x i ∂x j<br />
Para un material isótropo, el tensor <strong>de</strong> elasticidad <strong>de</strong>pen<strong>de</strong> <strong>de</strong> sólo dos constantes y pue<strong>de</strong><br />
escribirse<br />
D =2µ I+λ 1 ⊗ 1<br />
D ijkl = µ (δ ik δ jl + δ il δ jk ) + λ δ ij δ kl<br />
E<br />
Eν<br />
µ =<br />
λ =<br />
2 (1 + ν) (1 + ν) (1 − 2ν)<br />
don<strong>de</strong> µ y λ son los parámetros <strong>de</strong> Lamé, E es el módulo <strong>de</strong> elasticidad <strong>de</strong> Young y ν es la relación<br />
<strong>de</strong> Poisson. I es el tensor i<strong>de</strong>ntidad <strong>de</strong> cuarto or<strong>de</strong>n, 1 es el tensor i<strong>de</strong>ntidad <strong>de</strong> segundo or<strong>de</strong>n y<br />
⊗ <strong>de</strong>nota el producto tensorial<br />
σ ij = D ijkl ε kl = [µ (δ ik δ jl + δ il δ jk ) + λ δ ij δ kl ] ε kl<br />
σ ij = 2µ ε ij + δ ij λ ε kk<br />
97