Elementos Finitos Isoparametrico

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 4

ELEMENTOS FINITOS: ELEMENTO ISOPARAMETRICO

EJERCICIO

Demostrar que un elemento rectangular de cuatro nodos isoparimétricos como se muestra en


la figura (a). Se transforma en un elemento triangular de deformación constante figura (b).
Cuando se hacen coincidir los nodos 3 y 4.

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.

𝑋 = 𝑁₁(𝜉, 𝜂)х₁ + 𝑁₂(𝜉, 𝜂)х₂ + 𝑁₃(𝜉, 𝜂)х₃ + 𝑁₄(𝝃, 𝜼)х₄,


𝑌 = 𝑁₁(𝜉, 𝜂)у₁ + 𝑁₂(𝜉, 𝜂)у₂ + 𝑁₃(𝜉, 𝜂)у₃ + 𝑁₄(𝝃, 𝜼)у₄,
(9.51)
Ahora obtenemos que.
1 1 1 1
𝑋 = ∑𝑁ᵢ𝑋₁ = (1 − 𝜉)(1 − 𝜂)𝑋₁ + (1 + 𝜉)(1 − 𝜂)𝑋₂ + (1 + 𝜉)(1 + 𝜂)𝑋₃ + (1 − 𝜉)(1 + 𝜂)𝑋₄
4 4 4 4
1 1 1 1
𝑌 = ∑𝑁ᵢ𝑌₁ = (1 − 𝜉)(1 − 𝜂)𝑌₁ + (1 + 𝜉)(1 − 𝜂)𝑌₂ + (1 + 𝜉)(1 + 𝜂)𝑌₃ + (1 − 𝜉)(1 + 𝜂)𝑌₄
4 4 4 4
Al ser x₃=x₄, y₃=y₄, las expresiones se reducen a: (a)

1 1 1
𝑋 = (1 − 𝜉)(1 − 𝜂)𝑋₁ + (1 + 𝜉)(1 − 𝜂)𝑋₂ + (1 + 𝜉)(1 + 𝜂)𝑋₃
4 4 4
1 1 1
𝑌= (1 − 𝜉)(1 − 𝜂)𝑌₁ + (1 + 𝜉)(1 − 𝜂)𝑌₂ + (1 + 𝜉)(1 + 𝜂)𝑌₃
4 4 4

Estas expresiones significan que el elemento triangular se genera a partir de un elemento


cuadrilateral que tenga las funciones de interpolación.

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)

La matriz de deformación del elemento cuadrilateral degenerado en triangular,


según (9.69)

(𝐵𝑒) = (𝐼ᵅ)(𝑇ᵅ)(𝜕𝑁ᵅ). (9.69)


DONDE:
(𝐵𝑒) = (𝐼ᵅ)(𝑇ᵅ)(𝜕𝑁𝑒).

Donde (Iᵅ) es la matriz

1000
(Iᵅ)=(0 0 0 1) (𝐜)
0110

Para hallar la matriz inversa del jacobiano ampliada (𝑻ᵅ), calcularemos


previamente el jacobiano (Je). Sustituyendo las ecuaciones (b) en (9.57), Se
obtiene.
𝜕х 𝜕𝑦 1
− 𝑎(1 − 𝜂) 0
𝜕𝜂 𝜕𝜉
(𝑱𝒆) = | | = | 14 1
|
𝜕𝑥 𝜕𝑥
− − 𝑎(1 + 𝜉) 𝑏
𝜕𝜂 𝜕𝜉 4 2
De la que se deduce su determinante

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 − 𝜂)
𝜕𝜂 𝜕𝜉
𝐷𝑒𝑡(𝐽𝑒)

Con la que se compone la matriz inversa del jacobiano ampliada

2𝑏 0 0 0
2 𝑎(1 + 𝜉) a(1 − 𝜂) 𝑎(1 + 𝜉) a(1 − 𝜂)
(𝑇ᵅ) = | | (𝐝)
𝑎𝑏(1 − 𝜂) 0 0 2𝑏 0
0 0 𝑎(1 + 𝜉) a(1 − 𝜂)

Según (9.66), la matriz derivada de las funciones de interpretación ampliada del


elemento cuadrilateral degenerado en un elemento triangular es.

𝜕𝜂₁ 𝜕𝜂₂ 𝜕𝜂₃


0 0 0
𝜕𝜉 𝜕𝜉 𝜕𝜉
| |
𝜕𝜂₁ 𝜕𝜂₂ 𝜕𝜂₃
0 0 0
𝜕𝜂 𝜕𝜂 𝜕𝜂
(𝝏𝜼ᵅ) = | |
𝜕𝜂₁ 𝜕𝜂₂ 𝜕𝜂₃
0 0 0
𝜕𝜉 𝜕𝜉 𝜕𝜉
| 𝜕𝜂₁ 𝜕𝜂₂ 𝜕𝜂₃ |
0 0 0
𝜕𝜂 𝜕𝜂 𝜕𝜂

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

Multiplicando La matrices (c), (d) y (e) se obtiene la matriz de deformación del


elemento cuadrilateral degenerando en triangular.
−𝑏 0 𝑏 0 0 0
1
(𝑩𝒆) = | 0 −𝑎 0 0 0 𝑎 |
𝑎𝑏
−𝑎 −𝑏 0 𝑏 𝑎 0

Al ser constantes los terminos de la matriz de deformacion, el vector de


deformaciones (e)=(Be)(ᵹe) sera constante en cualquier punto del elemento.

También podría gustarte