vi
niO
PA CAE
-U
TE
Cy
Introducción al MEF
Ejercicio 1
M
Ejercicio 1: Obtener las funciones de forma y la matriz de deformación en
M
coordenadas normalizadas (o naturales) en los siguientes casos:
a) Elemento lineal axial
b) Elemento cuadrático axial (nodo intermedio centrado)
E-
'
# ( #&
!" # = $
#" ( #&
CA
"%&
Introducción al MEF
Ejercicio 1
vi
Ø Desarrollamos el productorio: 1 2
#
niO
#) = (1 #* = +1
'
# ( #&
!" # = $
#" ( #&
-U
"%&
# ( #* # ( ,+1- #(1 1
!) # = = = = ,1 ( #-
#) ( #* (1 ( ,+1- (2 2
# ( #) # ( ,(1- #+1 1
!* # = = = = ,1 + #-
#* ( #) +1 ( ,(1- 2 2
TE
Cy
Introducción al MEF
Ejercicio 1
M
Ø La matriz de deformación del elemento en axial viene dada como:
M
0[! # ] 0 ! # 0#
./ = =
03 0# 03
Ø Necesitamos la relación entre las dos variables para calcular el
E-
Jacobiano del cambio:
1) Planteamos la interpolación de la geometría isoparamétrica:
CA
3= ! # /
4 / = !) # 3) + !* # 3*
3
1 1
3= 1 ( # 7 5 + 1 + # 7 6/ 1
#
2
2 2
6/ #) = (1 #* = +1
3= ,1 + #-
2 3) = 5 3* = 6/
Introducción al MEF
Ejercicio 1
Ø 2) Calculamos la derivada de x con respecto a #:
vi
6/ 03 6/
niO
3= ,1 + #- 8 =
/ =
2 0# 2
Ø Finalmente calculamos la matriz de deformación del elemento:
-U
0[! # ] 0 ! # 0# 2 11 1 1
./ = = ./ = ( = ( / /
03 0# 03 6/ 22 6 6
1 0!) # 1
!) # = ,1 ( #- =(
2 0#
TE
2
1 0!* # 1
!* # = ,1 + #- 0#
=
2
2
Cy
Introducción al MEF
Ejercicio 1
M
Ø Elemento cuadrático con nodo intermedio centrado:
M
1 2 3
# = (1 #=5 # = +1
E-
Ø Se aplica la fórmula de Lagrange en coordenadas naturales:
,# ( #* - ,# ( #9 - 1
!) # = = # ,# ( 1-
,#) ( #* - ,#) ( #9 - 2
CA
,# ( #) - ,# ( #9 -
!* # = = ,1 + #-,1 ( #-
,#* ( #) - ,#* ( #9 -
,# ( #) - ,# ( #* - 1
!9 # = = # ,# + 1-
,#9 ( #) - ,#9 ( #* - 2
Introducción al MEF
Ejercicio 1
vi
Ø Necesitamos al igual que antes el Jacobiano de la trasnformación:
1
3 = !) # 3) + !* # 3* + !9 # 39 = # #(1 75+
niO
2
6/ 1
+ 1+# 1(# 7 + # 1 + # 7 6/
2 2
1 ( # * 6/ # + # * 6/ 6/
3= + = ,#-
-U
2 2 2
Ø Se calculo la derivada de x con respecto a la variable # (Jacobiano):
03 6/
0#
=
2
TE
Cy
Introducción al MEF
Ejercicio 1
M
Ø Por último calculamos la matriz de deformación del elemento:
0[! # ] 0 ! # 0#
M
./ = =
03 0# 03
1
!) # = # ,# ( 1- 0!) # 1
E-
2 0#
=#(
2
0!* #
!* # = ,1 + #-,1 ( #- 0#
= (2#
1 0!9 # 1
CA
!9 # = # ,# + 1- 0#
=#+
2
2
2 1 1
./ = #( ( 2# #+
6/ 2 2
Introducción al MEF
Ejercicio 2
vi
Ø Realizar la transformación isoparamétrica de un elemento cuadrático
con el nodo central no centrado
6
;
niO
'
3 = < !" ? 3"
1 2 3 ">)
L
-U
Obtener asimismo el vector de deformación Be a partir de las funciones de
forma en coordenadas naturales:
0!,3-
./ =
03
TE
Cy
Introducción al MEF
Ejercicio 2
M
Ø Realizar la transformación isoparamétrica de un elemento cuadrático
con el nodo central no centrado
6
M
1 2 3
E-
Las funciones de forma para un elemento cuadrático normalizado:
'
1
CA
!) ,#- = ( # 1 ( # 3 = < !" ? 3"
2
!* # = 1 ( # 1 + # ">)
1
!9 # = #,1 + #- 0!,3-
2 ./ =
03
Introducción al MEF
Ejercicio 2
vi
• Transformación isoparamétrica
3)
6
3 = !) !* !9 3* = !) ? 3) + !* ? 3* + !9 ? 39 = !) ? 5 + !* ? + !9 ? 6 =
39 ;
niO
6 6 6 6# * 6# 6# * 6
= 1(# * + #+# * = ( + + = 2 + ;# + # *
; 2 ; ; 2 2 @
• Jacobiano de la transformación
-U
03 6
8= = ; + 2#
0# @
• Vector de deformación
0!,#- d# 0! #
./ = = 8A)
TE
0# 03 0#
Cy
Introducción al MEF
Ejercicio 2
M
• Matriz de deformación
0!) # 1
= ( 1 ( 2#
0# 2
M
!* # ;
= (2# ./ = (,1 ( 2#- (B# ,1 + 2#-
0# ; + 2# 6
!9 # 1
= ,1 + 2#-
0# 2
E-
CA
Introducción al MEF
Ejercicio 3
vi
Ejercicio 3: Calcular la matriz exacta de rigidez mediante integración
numérica en los siguientes casos:
niO
a) Elemento lineal axial
b) Elemento cuadrático axial
-U
TE
Cy
Introducción al MEF
Ejercicio 3
M
a) Elemento lineal axial:
M
E) 1 E)
( 1 1 6/ FG
6 FG (
/ 1 (1
C/ = D 0# = D 0#
1 6 6 2
/ / 2 6 (1 +1
/
A) A)
6/
E-
Ø Aplicamos integración numérica. Función a integrar: H # = 1
CA
- Con la fórmula de la cuadratura de Gauss-Legendre de orden “p” se puede
integrar exactamente un polinomio de grado n=2p-1
Tabla:
El integrando es cte -> n=0 -> p=1 #=5
I) = 2
Introducción al MEF
Ejercicio 3
vi
a) Elemento lineal axial:
niO
JK 1 (1 ) JK 1 (1 JK 1 (1
C/ = NO>) H #) 7 I) = M 72= M
* LM (1 +1 * L (1 +1 L (1 +1
-U
Como debe de ser, nos sale la misma matriz de rigidez que trabajando en
coordenadas cartesianas.
TE
Cy
Introducción al MEF
Ejercicio 3
M
b) Elemento lineal cuadrático:
M
E)
C / = D ./ P FG ./ 8/ 0#
A)
E-
1
E) #(
2 2 2 1 1 6/
C/ = D / (2# FG / # ( (2# #+ 0#
6 1 6 2 2 2
A)
#+
CA
- Los términos del integrando son cuadráticos (lineal por lineal). Para calcular de
manera exacta la integral despejamos el número de puntos necesario de la
expresión:
Q = 2R ( 1 S 2 = 2R ( 1 S R = 1TU S 2 RVQWXY
Introducción al MEF
Ejercicio 3
1
E) #(
2
vi
2 2 1 1 6/
C/ = D / (2# FG / # ( (2# #+ 0#
6 1 6 2 2 2
A)
#+
2
niO
*
1 1 1
#( #( ,(2#- #* (
2 2 B
E)
FG 1
C/ =2 / D Z\^ (2# *
(2# #+ 0#
-U
6 2
A)
*
1
Z\^ Z\^ #+
2
1
Tabla de Gauss-Legendre: #) = (
TE
;
_ I) = 1
1
#* = + _ I* = 1
;
Cy
Introducción al MEF
Ejercicio 3
M
*
1 1 1
#( #( ,(2#- #* (
2 2 B
Ø En nuestro caso cada integral será aplicar: E)
FG 1
C/ = 2 D Z\^ (2# * (2# #+ 0#
6/ 2
A)
*
1
Z\^ Z\^ #+
M
E) *
1 1
D H # 0# = < H #" 7 I" = H ( 71+H 71
`>)
; ;
A)
E-
2FG 1
1
*
1 1
*
b FG 2FG 1 1 1 FG
a)) = / 7 ( ( + + ( = a)9 = 7 27 ( =
6 ; 2 ; 2 ; 6/ 6/ ; B ; 6/
2FG 1 1@ FG
2FG 1
1 1 1 1 1 e FG a** = 7 27B7 =
6 ; ; 6/
CA
/
a)* = / 7 ( ( +2 + ( (2 =
6 ; 2 ; ; 2 c; ; 6/
2FG 1 1 1 1 1 1 e FG
a*9 = 7 ( + +2 + + (2 =(
6/ ; 2 ; ; 2 ; ; 6/
* *
2FG 1
1 1 1 b FG
a99 = / 7 ( + + + + =
6 ; 2 ; 2 ; 6/
vi
Elementos Bidimensionales
niO
-U
TE
Cy
Introducción al MEF
Ejemplo
M
Ejercicio 1. Calcular las funciones de forma para un elemento
bidimensional Lagrangiano de cuatro nodos en coordenadas
M
naturales.
Datos: funciones de forma para el elemento unidimensional
lagrangiano en coordenadas naturales:
E-
i# ( #$
%
!" ,#- = $ k m
#" ( #$
$>)i,$%"-
Obtención de las funciones de forma lagrangianas bidimensionales:
CA
fg ,!T "- = lgh ,!-ilgj ,"-
Introducción al MEF
Ejemplo
Ø Las funciones de forma bidimensionales se obtienen como
vi
producto de las funciones de forma unidimensionales
correspondientes.
niO
1
n)" # = ,1 + # #" -
2 1
!" #T o = n)" # n)" o = ,1 + # #" -,1 + o o" -
1 B
n)" o = ,1 + o o" -
2
-U
Nodo pq rq
1 -1 -1
2 1 -1 TE
3 1 1
4 -1 1
Cy
Introducción al MEF
Ejemplo
M
Ø Sustituimos:
M
E-
CA
vi
Transformación isoparamétrica
niO
-U
TE
Cy
Introducción al MEF
Ejemplo
M
Ejercicio 2: Realizar la transformación isoparamétrica del rectángulo de la figura.
M
Indicar el valor del Jacobinao y calcular a través del mismo el área del elemento. Para
un valor de 2a=3 cm y 2b=1 cm, ¿el área calculado coincide con el valor real del
mismo?
E-
y
3
2b 4
1 2
CA
2a
Introducción al MEF
Ejemplo
vi
Ø Se plantea la transformación isoparamétrica del elemento:
niO
4
x = å N i xi = N1 × 0 + N 2 × 2a + N 3 × 2a + N 4 × 0 = N 2 × 2a + N 3 × 2a
i =1
4
y = å Ni yi = N1 × 0 + N 2 × 0 + N 3 × 2b + N 4 × 2b = N 3 × 2b + N 4 × 2b
-U
i =1
Se calcula el Jacobiano de la
transformación:
TE
Cy
Introducción al MEF
Ejemplo
M
Ø Calculamos los términos del Jacobiano:
M
E-
Se plantea la relación entre
diferenciales
CA
Introducción al MEF
Ejemplo
vi
Ø Para el caso de 2a=3 y 2b=1:
niO
-U
TE
Cy
Introducción al MEF
Ejemplo
M
Ø Para el caso de 2a=3 y 2b=1:
M
E-
CA
Introducción al MEF
Ejemplo
vi
Ø Para el caso de 2a=3 y 2b=1:
niO
Ø Resolvemos la integral del área en coordenadas normalizadas:
-U
)
; ; E) E) ;
G/ =s 0#0o = # t o t = 7 2 7 2 = ;
A) B B A) A) B
TE
Coincidente con el área real
Cy
M
M
Transformación e integración numérica (2D)
E-
CA
Introducción al MEF
Ejemplo
vi
Ejercicio 3 (mayo13:1.5/10): Calcular el área del elemento de 4 nodos de la figura,
obteniendo la integral en el dominio (1 u ! 1; "1 ! # ! 1 mediante integración
niO
numérica de Gauss-Legendre con un punto de integración. Justificar si se obtiene la
solución exacta al utilizar en este caso un único punto de integración. Y
4
Datos: coordenadas de los nodos 3
Nodo 1 (0,1), nodo 2 (0.7,0.7), nodo 3 (1.5,1.5) y nodo 4 (0,2)
-U
1 2
X
TE
Cy
Introducción al MEF
Ejemplo
M
M
E-
CA
Introducción al MEF
Ejemplo
vi
niO
-U
TE
Cy
M
M
Elementos Bidimensionales
E-
CA
Introducción al MEF
Ejemplo
vi
niO
-U
TE
Cy
Introducción al MEF
Ejemplo
M
Ø Funciones de forma del elemento:
M
E-
CA
Introducción al MEF
Ejemplo
vi
Ø Se plantea la transformación isoparamétrica del elemento:
&!
niO
'!
&"
! 0 0 0 0 '
" =
! " # $
% &"
0 ! 0 " 0 # 0 $ #
'#
&$
-U
'$
Ø Sustituyendo los valores de x e y del elemento pedido:
! ! ! !
&= 1() 1(* 2 + 1 + ) 1 ( * 5.5 + 1+) 1+* 7 + 1() 1(* 4
$ $ TE $ $
! ! ! !
'= 1() 1(* 3 + 1+) 1(* 3 + 1+) 1+* 6 + 1() 1(* 6
$ $ $ $
Cy
Introducción al MEF
Ejemplo
M
! ! ! !
&= 1() 1(* 2 + 1 + ) 1 ( * 5.5 + 1+) 1+* 7 + 1() 1(* 4
$ $ $ $
M
! ! ! !
'= 1() 1(* 3 + 1+) 1(* 3 + 1+) 1+* 6 + 1() 1(* 6
$ $ $ $
E-
Ø Ambos apartados se resuelven de la misma manera:
1) En el primero sustituyendo en las ecuaciones del cambio de variable el punto de coordenadas
CA
normalizadas ), * = -1,1/ se obtienen las coordenadas cartesianas del elemento & , ' =
7,6 .
2) En el segundo el proceso es igual pero con las coordenadas ), * = (0.5, (0.5 . Es importante
en este segundo caso comprobar que el valor del punto en coordenadas cartesianas esté dentro
de la zona de dominio del elemento en estudio.
