Transcript MEF
Elementos Finitos en un continuo elástico Objetivo: Determinar el campo de tensiones y deformaciones en un sólido contínuo Fases del método: 1. El continuo se divide, mediante líneas o superficies imaginarias, en un número de elementos finitos 2. Los elementos se encuentran conectados entre si mediante un número discreto de puntos que se denominan nodos, situados en su contorno. Los desplazamientos de estos nodos serán las incógnitas del problema. 3. Definición de funciones que definan de manera única el campo de desplazamientos dentro de cada elemento finito en función de los desplazamientos nodales. 4. Ecuaciones constitutivas del material que relaciona los desplazamientos y deformaciones con las tensiones. 5. Determinación de un sistema de fuerzas concentradas en los nodos, tal que equilibre las tensiones en el contorno y cualesquiera carga repartidas [email protected] Aplicación del método Pasos 1 y 2 División mediante líneas Creación de nodos y vj v(x,y) j uj u(x,y) k i x Paso 3.- Solución Aproximada p u ( x, y ) N i ( x, y )ui i 1 Funciones de interpolación p v( x, y ) N i ( x, y )vi i 1 u ( x, y ) N i v ( x, y ) 0 0 Nj 0 Ni 0 Nj ui v i u j vj Propiedades de las funciones de interpolación Elemento triangular u ( xi , yi ) ui N i ( xi , yi )ui N j ( xi , yi )u j N k ( xi , yi )uk vi N i ( xi , yi ) 1 N j ( xi , yi ) 0 N k ( xi , yi ) 0 u ( x j , y j ) u j N i ( x j , y j )ui N j ( x j , y j )u j N k ( x j , y j )uk Ni ( x j , y j ) 0 N j (x j , y j ) 1 Nk (x j , y j ) 0 u ( xk , yk ) uk N i ( xk , yk )ui N j ( xk , yk )u j N k ( xk , yk )uk N i ( xk , y k ) 0 N j ( xk , y k ) 0 N k ( xk , y k ) 1 i v(x,y) u(x,y) ui u( x, y) Ni ( x, y)ui N j ( x, y)u j Nk ( x, y)uk Elemento 2-D Ni(x,y) Nk(x,y) Nj(x,y) i i i k k k j j j ui i uj uk p u( x, y) Ni ui i 1 k j Paso 4.- Ecuaciones constitutivas {s}[D]({e} - {e0}) {s0} {s}.- Campo de tensiones [D].- Matriz de elasticidad {e}.- Campo de Deformaciones {e0}.- Deformaciones de origen térmico {s0}.- Tensiones residuales Caso 2-D Deformaciones ex e } e y xy x [L 0 y 0 y x u x v e } y u v y x x e } 0 y u ( x, y ) e } [L v ( x, y ) 0 u ( x, y ) y v( x, y ) x u ( x, y ) N i v ( x, y ) 0 u e e } [L[N e v N i 0 x N i [B 0 y N i N i y x 0 Nj 0 Ni 0 Nj ui v i u j v j u e e } [B e v N j x 0 0 N j N j y N j y x Tensiones Cálculo de las Tensiones s } [De } 0 1 E [D 1 0 2 1 - 0 0 (1 - ) 2 u e s } [D[B e v N i 0 x 1 s } E 2 1 0 0 1 - 0 0 (1 - ) 2 N i y 0 N i y N i x N j x 0 0 N j N j y N j y x e u e v Paso 5.- Sistema de fuerzas concentrada en los nodos Vi i Ui ui U i v V i i u j U j v V j j Trabajos Virtuales Al interior del elemento u} [N u } e e } [Bu } e u} [N u } e } [B u } e e e } s } Trabajo de las Tensiones u}T b} Trabajo de las fuerzas de cuerpo T En los nodos del elemento u } q } e T e Trabajo de las fuerzas nodales Igualando Trabajos e } s } u} b} } [B s } u } [N b} e T T e T T u u T T } q } u } ( [B s }dV - [N b}dV ) e T T e T e T V q } [B s }dV - [N b}dV T e V T V V {s}[D]({e} - {e0}) {s0} q } [B [D(e }- e }) s }}dV - [N b}dV T e T 0 0 V V {e} [Bue} q } [B [D([Bu }- e }) s }}dV - [N b}dV T e T e 0 V 0 V e T T T T q [B [D [B dV u - [B [D e 0 }dV [B s 0 }dV - [N b}dV V V V V } e } q } [K ]{u } { f } { f e e e [K [B [D[BdV T s } { fb } Matriz de Rigidez V fe } - [B [De }dV T 0 Fuerzas deformación inicial V fs } [B s }dV T 0 Fuerzas tensión residual V f b } - [N b}dV T V Fuerzas de cuerpo Aplicación Elemento Barra: x x1 A,E x2 N1 ( x) a1 b1 x 1 a1 b1 x1 N 2 ( x) a2 b2 x 0 a1 b1 x2 x2 - x N1 x2 - x1 x2 a1 x2 - x1 -1 b1 x2 - x1 x - x1 N2 x2 - x1 1 1 x1 x2 x1 x2 Ejemplo: Se tiene una barra como muestra la figura, se ha determinado por algún método el desplazamiento que tienen los nodos extremos, además se sabe que inicialmente se tiene una deformación de origen térmico de e 1.0x10-4 y una tensión residual de 10 x106 N/m2. La barra es de Acero E=2.0x1011 N/m2, una sección transversal de A= 0.001m2 Determine la tensión en el centro de la barra. 500 2 1 x u2=3 mm u1=1 mm Cálculo q } [K ]{u } { f } { f e e e s 0 } { fb } U1 q U 2 } e [B [L[N [N [N1 N1 N2 xj - x x j - xi -1 [B x j - xi u1 u u2 } e [ K ] [B [D[BdV T V [D [E d [L dx x - xi N2 x j - xi 1 x j - xi -1 x - x T j i [B 1 x j - xi -1 x - x -1 1 T 2 1 [B [D[B 1 [E x x x x 2 1 2 1 x2 - x1 1 - 1 E T [B [D[B x2 - x1 L 2 ( x2 - x1 ) - 1 1 E 1 - 1 T V [B [D[BdV L2 - 1 1 V dV EA 1 - 1 [K L - 1 1 Fuerzas nodales equivalente a deformación iinicial fe } - [B [De }dV T 0 V 1 - 1 fe } - [E e }dV L1 V 0 - 1 Ee 0 - 1 fe } dV EAe 0 L 1 V 1 Fuerzas nodales equivalente a tensiones residuales fs } [B T s }dV 0 V 1 - 1 fs } s }dV L1 V 0 - 1 fs } s 0 A 1 q } [K ]{u } { f } { f e e e s } U1 EA 1 - 1 u1 - 1 - 1 - EAe 0 As 0 1 1 U 2 L - 1 1 u2 Remplazando valores - 1 0.001 U1 8 1 3 - 1 3 - 1 4.0 x10 - 2.0 x10 1.0 x10 - 1 1 0.003 1 1 U 2 U1 - 799000 U 2 799000 Elemento Viga qi w( x) 1 2 x 3 x2 4 x3 wi qj L Condiciones de Borde Determinar la matriz de rigidez del elemento Viga wj w 1 2 x 3 x 4 x 2 3 Aplicando condiciones de borde se llega a: wi 1 qi 2 w j 1 2 L 3 L2 4 L3 q j 2 2 3 L 3 4 L2 Escrito en forma matricial se tiene: 1 0 1 0 0 1 wi 1 0 0 2 q i 2 3 L L L 3 w j 1 2 L 3L2 4 q j 0 0 wi 1 q i 2 3 2 3 1 - 2 wi - q i 2 w j - q j L L L 3 L 4 2 w 1 q - 2 w 1 q j j L3 i L2 L3 L2 Para la viga en estudio se tiene que: d 2w EI 2 Vi x - M i dx M ( x) Vi x - M i w( x) wi q i x (- 3 2 3 1 2 1 2 1 2 3 w q w q ) x ( w q w q ) x i i j j i i j j L2 L L2 L L3 L2 L3 L2 d 2w - 6 4 6 2 2 1 2 1 w q w q 6 ( w q w q j )x i i j j i j 2 2 2 3 2 i 3 2 dx L L L L L L L L Igualando términos Por equilibrio : Vi -V j M j Vi L - M i 12EI L3 6 EI 2 L - 12EI L3 6 EI L2 6 EI L2 4 EI L - 6 EI L2 2 EI L2 - 12EI L3 - 6 EI L2 12EI L3 - 6 EI L2 6 EI L2 w V 2 EI i i q 2 L i M i - 6 EI w j V j 2 L q j M j 4 EI L Se recomienda hacerla con los polinomios de interpolación: w( x) (1 - 3 2 2 3 2 2 1 3 3 2 2 3 1 2 1 3 x x ) w ( x x x ) q ( x x ) w ( x 2 x )q j i i j 2 3 2 2 3 L L L L L L L L 3 2 2 3 x 3x ) L2 L 2 1 N iq ( x ) ( x - x 2 2 x 3 ) L L 3 2 N jw ( x) ( 2 x 2 - 3 x 3 ) L L 1 1 N jq ( x ) ( - x 2 2 x 3 ) L L N iw ( x) (1 - d2 [L 2 dx [D EI [B [L[N [N [Niw [B - 26 123x L L Niq - N jw 4 6x L L2 N jq 6 12x L2 L3 - 2 6x L L2 [K [B [D[BdV T V - 6 12x L2 L3 4 6x L- 2 L L [EI - 26 123x 0 6 12x L L L2 L3 2 6x - 2 L L - 4 6x - 2 L L 12 EI - 6 12 x B1,1 EI 2 3 dx 3 L L L 0 L 2 6 12x - 3 2 L L - 2 6x 2 dx L L Se amplía la matriz a una viga que además esta sometida a tracción o compresión Mi Vi Mj Uj Ui EA L 0 0 EA L 0 0 0 0 12EI L3 6 EI L2 6 EI L2 4 EI L 0 0 - 12EI L3 6 EI L2 - 6 EI L2 2 EI L2 Vj - EA L 0 0 EA L 0 0 0 - 12EI L3 - 6 EI L2 0 12EI L3 - 6 EI L2 6 EI u Ui i 2 L v 2 EI i Vi M L2 q i i u U 0 j j v j V j - 6 EI M j q 2 j L 4 EI L 0 Para que sea válida en cualquier posición en el plano debe ser multiplicada por la matriz de rotación dada por: cos(q ) - sen(q ) sen(q ) cos(q ) 0 0 [R 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 cos(q ) - sen(q ) 0 0 sen(q ) cos(q ) 0 0 0 0 1 0 0 0 2 - x - x ) y -y y -y - x - x ) x -x y -y y -y x1 - x0)2 EA 1 2 y1 - y0)2 EI y1 - y0)2 EI 1 0 1 0 1 0 1 0 1 0 1 0 1 0 E A y - y ) 1 2 E I x - x ) 6 E I E A - 1 2 E A y - y ) - 1 2 E I x - x ) 6 E I 1 0 1 0 1 0 1 0 3 5 3 5 3 3 5 3 5 3 L L L L L L L L L L 2 2 2 2 y -y y -y ) x -x) x -x x -x y -y - y - y ) x -x) x -x - x1 - x0) 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 E A y - y ) 1 2 E I x - x ) E A 1 2 E I 6 E I E A y - y ) - 1 2 E I x - x ) E A - 1 2 E I 6 E I 1 0 1 0 1 0 1 0 3 5 3 5 3 3 5 3 5 3 L L L L L L L L L L y -y x -x y -y x -x 1 0 1 0 I 1 0 1 0 I 6 E I 6 E I 4 E ( -6) E I ( -6) E I 2 E 3 3 L 3 3 2 L L L L L T R K R 2 2 2 2 - x - x ) y y x x y y y y x x y y x x y y y y 1 0) EI 1 0) EA 1 2 1 0) EI 1 0) E A y - y 1 2 1 0 E I x - x (-6) 1 0 EI 1 0 1 0 1 0 1 0 E A - 1 2 E A y - y ) - 1 2 E I x - x ) ( -6) E I 1 0) 1 0) 1 0 1 0 3 5 3 5 3 3 5 3 5 3 L L L L L L L L L L 2 2 2 2 x x y y y y x x x x x x y y y y x x x x ) ) ) ) ) 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 E A y - y ) - 1 2 E I x - x ) E A - 1 2 E I ( -6) E I E A y - y ) 1 2 E I x - x ) E A 1 2 E I ( -6) E I 1 0 1 0 1 0 1 0 3 5 3 5 3 3 5 3 5 3 L L L L L L L L L L y -y x -x y -y x -x 1 0 1 0 I 1 0 1 0 I 6 E I 6 E I 2 E ( -6) E I ( -6) E I 4 E 3 3 2 3 3 L L L L L L Elemento triangular u( x, y) a0 a1 x a2 y vj v( x, y) a0 a1 x a2 y j vk uk k uj u ( xi , yi ) ui a0 a1 xi a2 yi vi u ( x j , y j ) u j a0 a1 x j a2 y j ui i 1 xi y y a0 ui 1 x y a u j j 1 j 1 xk yk a2 uk Desarrollamos para u(x,y) a0 a 1 a2 u ( xk , yk ) uk a0 a1 xk a2 yk x y - y x x y - y x -xi) yk yi xk j k j k i j i j u u u xj y k - y j xk - xi y k y i xk xi y j - y i xj i xj y k - y j xk - xi y k y i xk xi y j - y i xj j xj y k - y j xk - xi y k y i xk xi y j - y i xj k y -y -y ) y -yk) y j k i j i u u u i j k x y - y x - x y y x x y - y x x y - y x - x y y x x y - y x x y - y x - x y y x x y - y x j k j k i k i k i j i j j k j k i k i k i j i j j k j k i k i k i j i j x -x x -x -xk) xi k j j i u u u xj y k - y j xk - xi y k y i xk xi y j - y i xj i xj y k - y j xk - xi y k y i xk xi y j - y i xj j xj y k - y j xk - xi y k y i xk xi y j - y i xj k Campo de desplazamientos. u( x, y) [(x y - x y ) ( y - y )x (x - x ) yu [(x y - x y ) ( y - y )x (x - x ) yu [(x y - x y ) ( y - y )x (x - x ) yu j k k j j k k j i k i i k k i i k x j yk - y j xk - xi yk yi xk xi y j - yi x j j i j j i i j j k k u ( x, y ) u ( x, y ) [( x y j k [ - xk y j ) ( y j - yk ) x ( xk - x j ) y ui [( xk yi - xi yk ) ( yk - yi ) x ( xi - xk ) y u j ( xi y j - x j yi ) ( yi - y j ) x ( x j - xk ) y uk x j yk - y j xk - xi yk yi xk xi y j - yi x j (ai bi x ci y)ui (a j b j x c j y)u j (ak bk x ck y)uk x j yk - y j xk - xi yk yi xk xi y j - yi x j ai x j yk - xk y j 1 xi bi y j - yk ci xk - x j j yi xj yj 1 xj yj xk 1 xk y k yk - xi 1 yj 1 yk yi 1 xj 1 xk x j yk - y j xk - xi yk yi xk xi y j - yi x j i k u ( x, y ) v ( x, y ) (ai bi x ci y )ui (a j b j x c j y )u j (ak bk x ck y )uk (ai bi x ci y )vi (a j b j x c j y )v j (ak bk x ck y )vk Tensión deformación ex u x ey v y xy u v y x u 1 (bi ui b j u j bk uk ) x v 1 e y (ci vi c j v j ck vk ) y v u 1 xy (bi vi b j v j bk vk ci ui c j u j ck uk ) x y ex ui v ex bi 0 b j 0 bk 0 i bi 0 b j 0 bk 0 u 1 j 1 e y 0 ci 0 c j 0 ck [B } [B 0 ci 0 c j 0 ck c b c b c b v j ci bi ck b j ck bk xy i i k j k k u k vk ui v i u j } v j uk vk Tensión plana s x 1 0 ex E 1 0 e y [C e } s y 2 (1 - ) 0 0 1 - xy xy 2 ui v i s x 1 0 b 0 b 0 b 0 i j k u j E 1 0 . 0 ci 0 c j 0 ck [C [. B } s y 2 (1 - ) 0 0 1 - c b c b c b v j i i k j k k u xy 2 k vk Matriz de rigidez (1 - ) 2 Ejemplo: La estructura de la figura esta compuesta por una viga, una barra y una placa, en su extremo se aplica una carga vertical de 10000 N hacia abajo. Determinar los desplazamientos nodales, las fuerzas en cada uno de los elementos y la tensión máxima en cada uno de ellos. Barra circular d = 50 Placa e=10mm Viga cuadrada de a = 100 30º 10000 N 1000 1 1b 1p 2 3 1v