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 ...
¯K = ∫ L 0 BT D B d¯x 1 = ¯K = ∫ +1 −1 EA L B T D B L 2 dξ ⎡ 1 0 −1 0 ⎢ 0 0 0 0 ⎣ −1 0 1 0 0 0 0 0 Recordemos entonces que esta matriz representa el trabajo virtual interno a través de la expresión T.V.I. = (δū e ) T Para reescribir esta expresión en términos de desplazamientos respecto al sistema global de coordenadas debemos expresar los ū e y δū e en función de u e y δu e referidos al sistema global ¯K ū e ⎤ ⎥ ⎦ X 2 U 1 X 2 U 1 U 2 X 1 U 2 α X 1 Figura 3 barra de reticulado en coordenadas globales similarmente ū e = ⎡ ⎢ ⎣ ū 1 1 ū 1 2 ū 2 1 ū 2 2 La matriz Λ tiene la forma ⎤ ⎡ ⎥ ⎦ = ⎢ ⎣ cos α sin α 0 0 − sin α cos α 0 0 0 0 cos α sin α 0 0 − sin α cos α δū e = Λ δu e Λ = [ R 0 0 R ] ⎤ ⎡ ⎥ ⎢ ⎦ ⎣ u 1 1 u 1 2 u 2 1 u 2 2 ⎤ ⎥ ⎦ = Λ ue donde R es en este caso particular la matriz de rotación del sistema global al local y se cumple que Λ −1 = Λ T (en un caso general las matrices de transformación no son sencillamente matrices de cambio de coordenadas entre dos sistemas ortogonales). Luego el T.V.I. puede escribirse T.V.I. = (δu e ) T Λ} T {{ ¯K Λ} u e = (δu e ) T K K u e notar entonces que para transformar la matriz de rigidez la expresión correcta es con Λ T y no con Λ −1 , por ser la matriz de rigidez un tensor y no una aplicación lineal Similarmente para el trabajo virtual externo notar que: 66 T.V.E. = (δū e ) T ¯f = (δu e ) T Λ} {{ T ¯f } = (δu e ) T f f
Notar como es la matriz R, si definimos por t 1 al versor orientado del nudo 1 al nudo 2, y por t 2 al versor normal a t 1 de tal forma que t 3 = t 1 × t 2 , sea la normal saliente al plano, entonces es fácil ver que t 1 = [ cos α sin α que son precisamente las columnas de R ] R = [t 1 , t 2 ] t 2 = [ − sin α cos α En forma completamente similar puede obtenerse las matriz de un elemento de barra en 3 dimensiones. Ahora hay 3 componentes de desplazamiento por nudo y la matriz de rigidez en coordenadas locales resulta ¯K = EA L ⎡ ⎢ ⎣ 1 0 0 −1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 −1 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Llamando nuevamente por t 1 al versor orientado del nudo 1 al 2, y t 2 y t 3 a dos versores ortogonales a t 1 , tales que t 1 × t 2 = t 3 , entonces la matriz de rotación resulta R = [t 1 , t 2 , t 3 ] [ R 03×3 Λ = 0 3×3 R y finalmente la matriz de la barra en coordenadas globales es ] ⎤ ⎥ ⎦ ] K = Λ T ¯K Λ (4.2) Los elementos de barra articulada, además de utilizarse para modelos estructurales de enrejado, suelen utilizarse para modelar cables en general y tensores en particular. En estos casos debe notarse que los cables suelen presentar grandes desplazamientos, por lo cual el problema es no lineal y la matriz de rigidez debe actualizarse a la configuración deformada. 4.4.1. Ejercicio Calcular la matriz de rigidez de la barra articulada en coordenadas globales usando (4.2), mostrar que sólo depende de t 1 , es decir que no es necesario definir las direcciones t 2 y t 3 4.5. Elemento de viga en 3 dimensiones 4.5.1. Hipótesis más significativas (teoría clásica) Las secciones se mantienen planas Los desplazamientos son pequeños La geometría inicial y final son indistinguibles Las secciones son indeformables Las únicas tensiones relevantes son las que actúan en el plano de la sección (σ 11 , σ 12 , σ 13 ) Se desprecian las deformaciones asociadas al corte transversal. 67
- Page 13 and 14: Notar que K es no simétrica aun cu
- Page 15 and 16: integrando por partes el primer té
- 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 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 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
Notar como es la matriz R, si <strong>de</strong>finimos por t 1 al versor orientado <strong>de</strong>l nudo 1 al nudo 2, y por<br />
t 2 al versor normal a t 1 <strong>de</strong> tal forma que t 3 = t 1 × t 2 , sea la normal saliente al plano, entonces es<br />
fácil ver que<br />
t 1 =<br />
[<br />
cos α<br />
sin α<br />
que son precisamente las columnas <strong>de</strong> R<br />
]<br />
R = [t 1 , t 2 ]<br />
t 2 =<br />
[<br />
− sin α<br />
cos α<br />
En forma completamente similar pue<strong>de</strong> obtenerse las matriz <strong>de</strong> un elemento <strong>de</strong> barra en 3<br />
dimensiones. Ahora hay 3 componentes <strong>de</strong> <strong>de</strong>splazamiento por nudo y la matriz <strong>de</strong> rigi<strong>de</strong>z en<br />
coor<strong>de</strong>nadas locales resulta<br />
¯K = EA<br />
L<br />
⎡<br />
⎢<br />
⎣<br />
1 0 0 −1 0 0<br />
0 0 0 0 0 0<br />
0 0 0 0 0 0<br />
−1 0 0 1 0 0<br />
0 0 0 0 0 0<br />
0 0 0 0 0 0<br />
Llamando nuevamente por t 1 al versor orientado <strong>de</strong>l nudo 1 al 2, y t 2 y t 3 a dos versores<br />
ortogonales a t 1 , tales que t 1 × t 2 = t 3 , entonces la matriz <strong>de</strong> rotación resulta<br />
R = [t 1 , t 2 , t 3 ]<br />
[ R 03×3<br />
Λ =<br />
0 3×3 R<br />
y finalmente la matriz <strong>de</strong> la barra en coor<strong>de</strong>nadas globales es<br />
]<br />
⎤<br />
⎥<br />
⎦<br />
]<br />
K = Λ T ¯K Λ (4.2)<br />
Los elementos <strong>de</strong> barra articulada, a<strong>de</strong>más <strong>de</strong> utilizarse para mo<strong>de</strong>los estructurales <strong>de</strong> enrejado,<br />
suelen utilizarse para mo<strong>de</strong>lar cables en general y tensores en particular. En estos casos <strong>de</strong>be notarse<br />
que los cables suelen presentar gran<strong>de</strong>s <strong>de</strong>splazamientos, por lo cual el problema es no lineal y la<br />
matriz <strong>de</strong> rigi<strong>de</strong>z <strong>de</strong>be actualizarse a la configuración <strong>de</strong>formada.<br />
4.4.1. Ejercicio<br />
Calcular la matriz <strong>de</strong> rigi<strong>de</strong>z <strong>de</strong> la barra articulada en coor<strong>de</strong>nadas globales usando (4.2),<br />
mostrar que sólo <strong>de</strong>pen<strong>de</strong> <strong>de</strong> t 1 , es <strong>de</strong>cir que no es necesario <strong>de</strong>finir las direcciones t 2 y t 3<br />
4.5. Elemento <strong>de</strong> viga en 3 dimensiones<br />
4.5.1. Hipótesis más significativas (teoría clásica)<br />
Las secciones se mantienen planas<br />
Los <strong>de</strong>splazamientos son pequeños<br />
La geometría inicial y final son indistinguibles<br />
Las secciones son in<strong>de</strong>formables<br />
Las únicas tensiones relevantes son las que actúan en el plano <strong>de</strong> la sección (σ 11 , σ 12 , σ 13 )<br />
Se <strong>de</strong>sprecian las <strong>de</strong>formaciones asociadas al corte transversal.<br />
67