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

¯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

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

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

Saved successfully!

Ooh no, something went wrong!