2 Tipos de elemento ejercicio - SRK Consulting

Dr. Alejo O. Sfriso. Universidad de Buenos Aires materias.fi.uba.ar/6408 [email protected]. SRK Consulting (Argentina) latam.srk.com [email protected]. AOSA.
612KB Größe 14 Downloads 92 vistas
Tipos de elemento Ejercicio elemento cuatro nodos

Dr. Alejo O. Sfriso Universidad de Buenos Aires SRK Consulting (Argentina) AOSA

materias.fi.uba.ar/6408 latam.srk.com www.aosa.com.ar

[email protected] [email protected] [email protected]

Ejercicio elemento 4 nodos

Elemento 2D de 4 a 9 nodos

2

Para aumentar la precisión se pueden agregar más elementos o usar elementos con mas nodos Cuadrado • 4 nodos: interpolación bi-lineal

(Bathe 1996)

1

Ejemplo (mecánica continuo) ⎧⎪ 1.5 + 0.1xy ⎫⎪ u ⎡⎣ x, y ⎤⎦ = x − x = ⎨ ⎬ ⎪⎩ 0.5 + 0.1xy ⎪⎭ ⎧⎪ 1.5 + 0.1x0.5x0.4 ⎫⎪ ⎧⎪ 1.52 ⎫⎪ u ⎡⎣0.5,0.4 ⎤⎦ = ⎨ ⎬= ⎨ ⎬ ⎩⎪ 0.5 + 0.1x0.5x0.4 ⎭⎪ ⎩⎪ 0.52 ⎪⎭

Ejercicio elemento 4 nodos

end

ini

⎧ ∂u1 1 ⎛ ∂u1 ∂u2 ⎞ ⎪ + ∂x 2 ⎜⎝ ∂x2 ∂x1 ⎟⎠ ⎪ 1 ε ⎡⎣ x, y ⎤⎦ = ⎨ ∂u2 ⎪ 1 ⎛ ∂u1 ∂u2 ⎞ ⎪ 2 ⎜⎝ ∂x + ∂x ⎟⎠ ∂x 2 1 2 ⎩ ⎧⎪ 0.04 0.045 ⎫⎪ ε ⎡⎣0.5,0.4 ⎤⎦ = ⎨ ⎬ ⎪⎩ 0.045 0.05 ⎪⎭

∂u1 ∂ ⎡⎣1.5 + 0.1xy ⎤⎦ = = 0.1y ∂x1 ∂x ∂u1 ∂ ⎡⎣1.5 + 0.1xy ⎤⎦ = = 0.1x ∂x2 ∂y ∂u2 ∂ ⎡⎣0.5 + 0.1xy ⎤⎦ = = 0.1y ∂x1 ∂x ∂u2 ∂ ⎡⎣0.5 + 0.1xy ⎤⎦ = = 0.1x ∂x2 ∂y

⎫ ⎪ 0.1y 0.05( x + y ) ⎪ ⎧⎪ ⎬= ⎨ 0.1x ⎪ ⎪⎩ 0.05( x + y ) ⎪ ⎭

⎫ ⎪ ⎬ ⎪⎭

3

Ejercicio elemento 4 nodos

Ejemplo 2D: desplazamientos FEM

4

1

ue = f[Ue] = H x Ue

4

3

2

⎧ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ue = H ⋅U e = H ⋅ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎩

U1 ⎫ ⎪ ⎧ V1 ⎪ ⎪ ⎪ ⎪ U2 ⎪ ⎪ ⎪⎪ V2 ⎪⎪ ⎬ = H⋅⎨ U3 ⎪ ⎪ ⎪ ⎪ V3 ⎪ ⎪ ⎪ U4 ⎪ ⎪⎩ ⎪ V4 ⎪ ⎭

1.5 0.5 1.5 0.5 1.6 0.6 1.5 0.5

⎫ ⎪ ⎪ ⎪ ⎪⎪ ⎬ ⎪ ⎪ ⎪ ⎪ ⎪⎭

2

Ejercicio elemento 4 nodos

Ejemplo 2D: desplazamientos FEM

4

3

1

2

⎡ (1− x ) (1− y ) H=⎢ ⎢ 0 ⎣

0

(1− x )(1− y )

x (1− y ) 0

0

xy

0

x (1− y )

(1− x ) y

0

xy

0

⎤ ⎥ (1− x ) y ⎥⎦ 0

5

Ejercicio elemento 4 nodos

Ejemplo 2D: desplazamientos FEM

⎡ (1− x ) (1− y ) H=⎢ ⎢ 0 ⎣

0

xy

0

x (1− y )

(1− x ) y

0

xy

0

⎧ ⎪ ⎪ ⎪ ⎡ 0.3 0 0.3 0 0.2 0 0.2 0 ⎤ ⎪⎪ ue = H ⋅U e = ⎢ ⎥⎨ ⎣ 0 0.3 0 0.3 0 0.2 0 0.2 ⎦ ⎪ ⎪ ⎪ ⎪ ue = f[Ue] = H x Ue ⎪⎩

1.5 0.5 1.5 0.5 1.6 0.6 1.5 0.5

0

(1− x )(1− y )

x (1− y ) 0

⎤ ⎥ (1− x ) y ⎥⎦ 0

⎫ ⎪ ⎪ ⎪ ⎪⎪ ⎧⎪ 1.52 ⎫⎪ ⎬= ⎨ ⎬ ⎪ ⎩⎪ 0.52 ⎭⎪ ⎪ ⎪ ⎪ ⎪⎭

6

3

Ejercicio elemento 4 nodos

Ejemplo 2D: deformaciones FEM ⎡ − (1− y ) 0 1− y ⎢ ∂H ⎢ B= = 0 − (1− x ) 0 ∂x ⎢ ⎢ − (1− x ) 2 − (1− y ) 2 − x 2 ⎣

0

y

0

−y

−x

0

x

0

(1− y )

2 x 2

y 2

(1− x )

⎧ ⎪ ⎪ ⎪ ⎡ −0.6 0 0.6 0 0.4 0 −0.4 0 ⎤⎪ ⎥⎪ ε = B ⋅ U e = ⎢⎢ 0 −0.5 0 −0.5 0 0.5 0 0.5 ⎥ ⎨ ⎢⎣ −0.25 −0.3 −0.25 0.3 0.25 0.2 0.25 −0.2 ⎥⎦ ⎪ ⎪ ⎪ ⎪ ue = H x Ue ⎪⎩

1.5 0.5 1.5 0.5 1.6 0.6 1.5 0.5

ε = B x Ue

⎤ ⎥ 1− x ⎥ ⎥ 2 −y 2 ⎥ ⎦ ⎫ ⎪ ⎪ ⎪ ⎧ ⎪⎪ ⎪ 0.04 ⎬ = ⎨ 0.05 ⎪ ⎪ 0.045 ⎪ ⎩ ⎪ ⎪ ⎪⎭ 0

⎫ ⎪ ⎬ ⎪ ⎭

7

Ejemplo 2D: tensiones (elasticidad)

Ejercicio elemento 4 nodos

ìs xx ü E ï ï ís yy ý = 2 ïs ï 1 -n î xy þ

8

0 ü ìe xx ü ì1 n ï ïï ï 0 ý íe yy ý ín 1 ï0 0 1 -n ï ïe ï î þ î xy þ

E = 1000 MPa ν = 0.20 ⎧ σ ⎪⎪ xx ⎨ σ yy ⎪ ⎪⎩ σ xy

ue = f[Ue]⎧⎪ σ xx ⎪ ε = f[ue] ⎨ σ yy ⎪ σ = f[ε] ⎪⎩ σ xy

⎫ ⎧ ε ⎪⎪ 1000 MPa ⎧⎪ 1 0.2 0 ⎫⎪ ⎪⎪ xx ⎬= ⎨ 0.2 1 0 ⎬ ⎨ ε yy 1− 0.22 ⎪ ⎪ 0 0.8 ⎪⎭ ⎪ ε xy ⎩ 0 ⎪⎭ ⎩⎪

⎫ ⎧ 1 0.2 0 ⎪⎪ 1000 MPa ⎪ = ⎬ ⎨ 0.2 1 0 1− 0.22 ⎪ ⎪ 0 0.8 ⎩ 0 ⎪⎭⎡0.5,0.4⎤ ⎣ ⎦

⎫ ⎪⎪ ⎬ ⎪ ⎪⎭

⎫ ⎧ 0.04 ⎫ ⎧ 52 ⎪⎪ ⎪ ⎪ ⎬ ⎨ 0.05 ⎬ = ⎨ 60 ⎪ ⎪ 0.045 ⎪ ⎪ 38 ⎭⎩ ⎭ ⎩

⎫ ⎪ ⎬ MPa ⎪ ⎭

4