Análisis por Elementos Finitos
Clase 6
Ing. Eduardo Orcés P.
Junio 20/2011
Ing. Eduardo Orcés
2011-I
TEMAS
Objetivos de aprendizaje
Análisis de Pórticos en 2D y 3D
Programación en MATLAB / CALFEM
Ing. Eduardo Orcés 2011-I
Pórticos
Son estructuras reticulares formadas por
vigas unidas entre sí con diferentes grados de
rigidez y orientación arbitraria en el espacio.
Las vigas pueden soportar cargas axiales y
transversales que producen esfuerzos de
flexión, torsión y cortante.
Ing Eduardo Orcés 2011-I
Pórticos
Objetivos de aprendizaje:
o Entender el comportamiento estructural de los pórticos.
o Estar en capacidad de crear modelos estructurales
apropiados para pórticos en 2D y 3D.
o Aplicar el método de EF para analizar pórticos usando
Matlab/Calfem, Mastan2, VisualFEA.
Ing Eduardo Orcés 2011-I
Pórticos
Viga con orientación arbitraria (en 2D):
En coordenadas locales, se obtiene la matriz de rigidez
k’ del elemento, sumando los efectos axiales a los de
flexión y cortante. Esto es válido siempre que la
relación de esbeltez L/r sea pequeña (no hay pandeo).
Ing Eduardo Orcés 2011-I
• Transformando a coordenadas globales:
q ' R q, p ' R p
c s 0 0 0 0
s c 0 0 0 0
0 0 1 0 0 0
donde : R , s sen , c cos
0 0 0 c s 0
0 0 0 s c 0
0 0 0 0 0 1
p' k ' q' , R p k' R q p RT k ' R q k q
k RT k ' R
Ing Eduardo Orcés 2011-I
• La expresión explícita para k es la siguiente:
Ing Eduardo Orcés 2011-I
Ejemplo: Para el pórtico plano mostrado,
determinar los desplazamientos y fuerzas
nodales. Datos: E = 30 x 106 psi, I = 65 plg4,
A = 6.8 plg2.
Ing Eduardo Orcés 2011-I
• Solución: Cargas equivalentes nodales:
• En coordenadas locales, la matriz de rigidez para un
elemento de viga en 2D es:
Ing Eduardo Orcés 2011-I
• Elemento 1: R = I
• Elementos 2 y 3: En el sistema de coordenadas local, donde i=3, j=1
para el elemento 2, y i=4, j=2 para el elemento 3.
Ing Eduardo Orcés 2011-I
• En general, la matriz de rotación R para el elemento viga es:
l m 0 0 0 0
m m 0 0 0 0
0 0 1 0 0 0
R , l cos , m sen
0 0 0 m n 0
0 0 0 n m 0
0 0 0 0 0 1
• Para los elementos 2 y 3 tenemos que l=0. m=1, por lo que:
• Usando la transformación k = RT.k’.R para los elementos 2 y 3,
obtenemos las matrices de rigidez en el sistema de coordenadas
global:
Ing Eduardo Orcés 2011-I
• Condiciones de Borde y Carga:
Ing Eduardo Orcés 2011-I
• Ensamblando para obtener las ecuaciones globales de EF y
aplicando las C.B., obtenemos las ecuaciones para los GDL activos:
• Resolviendo:
Ing Eduardo Orcés 2011-I
• Para obtener las reacciones, substituimos los desplazamientos
obtenidos anteriormente en las ecuaciones globales (GDL inactivos)
de EF, para obtener:
• Chequeo:
Ing Eduardo Orcés 2011-I
Ejercicio 6.1: Hallar las fuerzas nodales en los
elementos del pórtico plano del ejemplo anterior.
Ing Eduardo Orcés 2011-I
Elemento viga en 3D:
• En coordenadas centroidales principales (locales):
Ing Eduardo Orcés 2011-I
• Superponiendo los efectos de axial, flexión y torsión, y
colocando en los grados de libertad respectivos, se
obtiene:
p' k ' q'
• donde:
Ing Eduardo Orcés 2011-I
Matriz de rotación de coordenadas:
Ing Eduardo Orcés 2011-I
• Los vectores (desplazamientos, fuerzas) están
relacionados por:
• p = R.Q, donde R es la matriz de rotación:
Ing Eduardo Orcés 2011-I
• Los cosenos directores li, mi, ni están dados por:
l1, m1, n1 = cosenos directores de eje x’ con ejes x, y, z.
l2, m2, n2 = cosenos directores de eje y’ con ejes x, y, z.
l3, m3, n3 = cosenos directores de eje z’ con ejes x, y, z.
l1 = (xj - xi)/L, m1 = (yj - yi)/L, n1 = (zj - zi)/L.
Notar que en la matriz de rotación R solo aparecen los
cosenos directores del eje x’ con respecto a los ejes
globales.
• Determinamos la matriz de rotación R1 , rotando los ejes
globales hasta hacerlos coincidir con los ejes principales
del elemento ij :
Ing Eduardo Orcés 2011-I
• Primero se rota un ángulo β alrededor del eje y. El eje x
se hace coincidir con el eje xβ.
cos 0 sen
Ry ( ) 0 1 0
- sen 0 cos
Ing Eduardo Orcés 2011-I
• Luego se rota un ángulo γ alrededor del eje zβ. El nuevo
eje xγ coincide con el eje local x . El eje zγ permanece
paralelo al plano xz.
cos sen 0
Rz ( ) - sen cos 0
0 0 1
• Estas dos rotaciones pueden combinarse en una sola
rotación:
R zy Rz( ) Ry( )
Ing Eduardo Orcés 2011-I
• Después de esta rotación, el eje zγ es paralelo al plano
xz , y el eje yγ tiene una proyección vertical en la
dirección positiva y. Notar que el plano xγ –yγ es
vertical, por lo que el miembro ij no ha sido rotado
alrededor de su eje longitudinal. En muchos casos,
esto es suficiente pues el plano xγ –yγ es un plano
principal de la sección.
• En ciertos casos, hay que rotar un ángulo α alrededor
del eje x , para alinear los ejes globales con los
principales de la sección.
• Consideremos, por ejemplo, el caso de un perfil I cuyo
eje principal y está inclinado un ángulo α con respecto
al eje vertical.
Ing Eduardo Orcés 2011-I
• El plano x -yγ es vertical y define el ángulo α que hay
que rotar hasta un plano principal:
1 0 0
Rx( ) 0 cos sen
0 - sen cos
• Finalmente se obtiene la matriz R1:
R1 R x (α ) R zy
• Para determinar el ángulo α , generalmente se
selecciona un nodo de referencia k (el cual puede ser
ímaginario ) en el plano x y , pero no sobre el eje x .
Ing Eduardo Orcés 2011-I
• Los cosenos directores de la matriz de rotación se
pueden calcular a partir de las coordenadas de los 2
nodos del elemento y del nodo k en el plano x y del
elemento.
• Cosenos directores del eje x :
l11 = (xB - xA)/L
l12 = (yB - yA)/L
l13 = (zB - zA)/L
• Cosenos directores de eje z :
Z X K
X ( xB x A )i ( y B y A ) j ( z B z A )k
K ( xK x A )i ( y K y A ) j ( z K z A )k
Zx ( yB y A )( zk z A ) ( zB z A )( yk yA )
Zy ( zB z A )(xk z A ) ( xB x A )( zk zA )
Zz ( xB x A )( yk z A ) ( yB y A )( xk xA )
Ing Eduardo Orcés 2011-I
Zx Zy Zz
l31 , l32 , l33 ,Z Z x2 Z y2 Z z2
Z Z Z
• Cosenos directores de l eje y :
Y Z X
Z l31i l32 j l33k
X l11i l12 j l13k
Yx l13l32 l12l33 , Yy l11l33 l13l31 ,
Yz l12l31 l11l32
Yx Yy Yz
l21 , l22 , l23
Y Y Y
Y Yx2 Yy2 Yz2
Ing Eduardo Orcés 2011-I
• Notar que la dirección del eje y es en el
sentido del nodo de referencia k .
• Una vez calculados los lij , se puede
hallar la matriz de rotación R:
R1 0 0 0
l11 l12 l13
0 R1 0 0
R , R1 l21 l22 l23
0 0 R1 0
l31 l32 l33
0 0 0 R1
Ing Eduardo Orcés 2011-I
• La matriz de rigidez del elemento en coordenadas
globales es:
(e) T
k R k' R
• Si existe una carga distribuida con componentes wy’, wz’
(fuerza/unidad de longitud) aplicada al elemento,
entonces las cargas equivalentes nodales son:
2 2
w y ´l w z ´l w z ´l 2 w y ´l w y ´l w z ´l w z ´l 2 w y ´l
fe ' 0 , 0 - , 0 , 0 -
2 2 12 12 2 2 12 12
• Estas cargas se las transforma a componentes globales:
T
fe R fe '
Ing Eduardo Orcés 2011-I
• Luego se ensambla la matriz de rigidez K y el vector de
fuerzas F, para formar el sistema de ecuaciones K.Q = F.
Después de aplicar las C.B.’s se obtienen los grados de
libertad ‘activos’ Qa .
• Finalmente, usando la conectividad, se pueden hallar las
fuerzas en los elementos:
f ' k ' q fe '
Ing Eduardo Orcés 2011-I
Ejemplo: Analizar el pórtico 3D mostrado. Todas
las vigas son secciones de acero de ala ancha
W8x24, E = 29 x 106 psi, G = 11.15 x 106 psi, Iy =
18.3 plg4, Iz = 82.8 plg4, J = 0.35 plg4, A = 7.08
plg2.
• Nodos de
referencia:
• Elemento 1:
nodo 5 (0,0,0)
• Elemento 2:
nodo 4 (180,0,0)
• Elemento 3:
nodo 5 (0,0,0)
Ing Eduardo Orcés 2011-I
• Cosenos directores:
• Elemento 1: L=180
l11=(0-180)/180=-1, l12=l13=0
Zx=(yB-yA)(zk-zA)-(zB-zA)(yk-yA)
=(0)(0-180)-(0)(0)=0
Zy= -1802, Zz= 0.
Z= 1802
Ing Eduardo Orcés 2011-I
• Cosenos directores:
• Elemento 1:
l31=0, l32=-1, l33=0
Yx=(l13l32-l13l32)=0-0=0, Yy=0,
10 0 Yz=-1
R1 0 0 -1 l21=l22=0, l33=-1
0 -1 0
Ing Eduardo Orcés 2011-I
• Cosenos directores:
• Elemento 2: 0 1 0
R2 0 0 -1
-1 0 0
• Elemento 3: 0 0 -1
R3 -1 0 0
0 1 0
Ing Eduardo Orcés 2011-I
• Hallamos la matriz de rigidez global para
cada elemento:
Ri 0 0 0
T
0 Ri 0 0
ki R k i ' R, R , i 1,2,3
0 0 Ri 0
0 0 0 Ri
• Puesto que el nodo 1 es el único con
GDL’s diferentes de cero, y puesto que es
el primer nodo de cada elemento, solo se
requiere la parte superior izquierda (6x6)
de la matriz ki de cada elemento:
Ing Eduardo Orcés 2011-I
1'140 667 0 0 0 0 0
. 1092 0 0 0 - 98 278
. . 4941 0 444 667 0
ki
. sim. . . 21 681 0 0
. . . . 53'360000 0
. . . . . 4'793333
• Similar para k2 y k3.
• Ensamblando las matrices anteriores se
obtienen las ecuaciones de equilibrio para
los GDL’s activos de la estructura
Ing Eduardo Orcés 2011-I
0 1'146 699 0 0 0 - 444 667 - 98 278 u1
0 . 1'142851 0 98 278 0 - 98 278 u2
0 . . 1'150 548 444 667 444 667 0 u3
0 . sim. . . 65'175010 0 0 u4
492000 . . . . 106'741700 0 u5
0 . . . . . 23'608350 u6
• Resolviendo para los GDL’s activos:
u1 0.001794"
u2 0
u3 - 0.001792"
u4 0.000012 rad
u5 0.004624 rad
u6 0.000007 rad
Ing Eduardo Orcés 2011-I
• Fuerzas en los elementos: p' k ' q' k ' R q
Elemento1 Elemento 2 Elemento 3
p1 - 2046 - 0.47 2044
p2 - 2047 3.42 2047
p3 0.73 - 1.23 1.2
p4 - 0.27 100.3 0.16
p5 - 88 88.3 144.2
p6 - 245 950 144.4 245 949
p7 2046 0.47 - 2044
p8 2047 - 3.42 - 2047
p9 - 0.73 1.23 - 1.2
p10 0.27 - 100.3 0.16
p11 - 44 132.27 - 72.05
p12 - 122 577 470.6 122 576
Ing Eduardo Orcés 2011-I
Ing Eduardo Orcés 2011-I
