Elementos Finitos Isoparametrico
Elementos Finitos Isoparametrico
Elementos Finitos Isoparametrico
EJERCICIO
a y
4 3 3=4
1 2 x 1 2 x
(a) (b)
SOLUCIÓN: Según (9.44), Las funciones de interpolación del elemento cuadrilateral de cuatro
nodos son:
1 1 1 1
𝑁₁ = (1 − 𝜉)(1 − 𝜂), 𝑁₂ = (1 + 𝜉)(1 − 𝜂), 𝑁₃ = (1 + 𝜉)(1 + 𝜂), 𝑁₄ = (1 − 𝜉)(1 + 𝜂)
4 4 4 4
(9.44)
Ahora las ecuaciones de trasformación de coordenadas del elemento rectangular
isoparimétricos, según (9.51), son.
1 1 1
𝑋 = (1 − 𝜉)(1 − 𝜂)𝑋₁ + (1 + 𝜉)(1 − 𝜂)𝑋₂ + (1 + 𝜉)(1 + 𝜂)𝑋₃
4 4 4
1 1 1
𝑌= (1 − 𝜉)(1 − 𝜂)𝑌₁ + (1 + 𝜉)(1 − 𝜂)𝑌₂ + (1 + 𝜉)(1 + 𝜂)𝑌₃
4 4 4
1 1 1
𝑁₁ = (1 − 𝜉)(1 − 𝜂), 𝑁₂ = (1 + 𝜉)(1 − 𝜂), 𝑁₃ = (1 + 𝜂), 𝑁₄ = 0
4 4 2
Sustituyendo en (a) las coordenadas nodales X₁=0, Y₁=0, X₂=a, Y₂=0, X₃=0, Y₃=b, se
obtiene.
1 1
𝑋 = 4 (1 − 𝜉)(1 − 𝜂)𝑋𝑎 𝑌 = 4 (1 − 𝜉)(1 − 𝜂)𝑏 (b)
1000
(Iᵅ)=(0 0 0 1) (𝐜)
0110
1 1 𝑎𝑏
det (𝐽𝑒) = 𝑎(1 − 𝜂) 𝑏 = (1 − 𝜂)
4 2 8
Por lo consiguiente, de acuerdo con (9.62), la matriz inversa del jacobiano es
𝜕х 𝜕𝑦
−
𝜕𝜂 𝜕𝜉 2
| | 2𝑏 0
(𝑻𝒆) = (𝑱𝒆) = 𝜕𝑥 𝜕𝑥 =| || |
| − | 𝑎𝑏(1 − 𝜂) 𝑎(1 + 𝜉) a(1 − 𝜂)
𝜕𝜂 𝜕𝜉
𝐷𝑒𝑡(𝐽𝑒)
2𝑏 0 0 0
2 𝑎(1 + 𝜉) a(1 − 𝜂) 𝑎(1 + 𝜉) a(1 − 𝜂)
(𝑇ᵅ) = | | (𝐝)
𝑎𝑏(1 − 𝜂) 0 0 2𝑏 0
0 0 𝑎(1 + 𝜉) a(1 − 𝜂)
O sea.
−1 + 𝜂 0 1−𝜂 0 0 0
1 1+𝜉 0 −1−𝜉 0 2 0
(𝝏𝜼ᵅ) = | | (e)
4 0 −1+𝜂 0 −1−𝜂 0 0
0 1+𝜉 0 −1−𝜉 0 2