Practica MEF 2D Tension Plana

1.
Aplicar el Método de Elementos Finitos para

determinar:
a. El estado de Deformación
b. Los tensores principales
Para la estructura que se muestra en la figura:

Donde:
E = 175,000kg/cm2
v = 0.25
t = 12 mm
Tipo de problema: Esfuerzo plana
I. Discretización del continuo. – Adoptando elementos triangulares, tenemos
Coordenadas Elemento
i x [m] y [m] i 1 2 3
1 0.0 0.0 1 1 2 5
2 0.7 0.0 2 2 3 6
3 1.4 0.0 3 1 5 4
4 0.0 0.5 4 2 6 5
5 0.7 0.5
6 1.4 0.5
II. Matriz constitutiva. - Para todos los elementos
Nota: Esta matriz es constante para todos los elementos del modelo.
III. Cálculo de las matrices de rigidez

Elemento 1
X Y
1 0 0 t= 0.012 m
2 0.7 0 2A= 0.35
5 0.7 0.5 A= 0.175
Donde la matriz de deformación es:
-0.5 0 0.5 0 0 0
[B] = 1/|J| 0 0 0 -0.7 0 0.7
0 -0.5 -0.7 0.5 0.7 0
La matriz de rigidez está dada por:
1 2 3 4 9 10
8000 0 -8000 2800 0 -2800 1
0 3000 4200 -3000 -4200 0 2
-8000 4200 13880 -7000 -5880 2800 3
[K]= 2800 -3000 -7000 18680 4200 -15680 4
0 -4200 -5880 4200 5880 0 9
-2800 0 2800 -15680 0 15680 10
Elemento 2
X Y
2 0.7 0 t= 1.00 m
3 1.4 0 2A= 1.30
6 1.4 0.5 A= 0.65
-0.5 0 0.5 0 0 0
[B] = 1/|J| 0 0 0 -0.7 0 0.7
0 -0.5 -0.7 0.5 0.7 0
3 4 5 6 11 12
8000 0 -8000 2800 0 -2800 3
0 3000 4200 -3000 -4200 0 4
-8000 4200 13880 -7000 -5880 2800 5
[K]= 2800 -3000 -7000 18680 4200 -15680 6
0 -4200 -5880 4200 5880 0 11
-2800 0 2800 -15680 0 15680 12
Elemento 3
X Y
1 0 0 t= 0.012 m
5 0.7 0.5 2A= 0.35
4 0 0.5 A= 0.175
0 0 0.5 0 -0.5 0
[B] = 1/|J| 0 -0.7 0 0 0 0.7
-0.7 0 0 0.5 0.7 -0.5
1 2 9 10 7 8
5880 0 0 -4200 -5880 4200 1
0 15680 -2800 0 2800 -15680 2
0 -2800 8000 0 -8000 2800 9
[K]= -4200 0 0 3000 4200 -3000 10
-5880 2800 -8000 4200 13880 -7000 7
4200 -15680 2800 -3000 -7000 18680 8
Elemento 4:
X Y
2 0.7 0 t= 0.012 m
6 1.4 0.5 2A= 0.35
5 0.7 0.5 A= 0.175

0 0 0.5 0 -0.5 0
[B] = 1/|J| 0 -0.7 0 0 0 0.7
-0.7 0 0 0.5 0.7 -0.5

3 4 11 12 9 10
5880 0 0 -4200 -5880 4200 3
0 15680 -2800 0 2800 -15680 4
0 -2800 8000 0 -8000 2800 11
[K]= -4200 0 0 3000 4200 -3000 12
-5880 2800 -8000 4200 13880 -7000 9
4200 -15680 2800 -3000 -7000 18680 10
Ensamblando en la matriz de rigidez total, tenemos
IV. Matriz de rigidez y vector de fuerzas totales
13880 0 -8000 2800 -5880 4200 0 -7000 1

0 18680 4200 -3000. 2800 -15680 -7000 0 2
-8000 4200 27760 -7000. -8000. 2800 -11760 7000 0 -7000 3
2800 -3000 -7000 37360 4200 -3000. 7000 -31360 -7000 0 4
-8000 4200 13880 -7000. -5880 2800 5
2800 -3000. -7000. 18680 4200 -15680 6
[K]=
-5880 2800 13880 -7000. -8000 4200 7
4200 -15680 -7000. 18680 2800 -3000. 8
0 -7000 -11760 7000 -8000. 2800 27760 -7000. -8000. 4200 9
-7000 0 7000 -31360 4200 -3000. -7000. 37360 2800 -3000 10
0 -7000 -5880 4200 -8000. 2800 13880. 0 11
-7000 0 2800 -15680 4200 -3000. 0 18680 12
1 2 3 4 5 6 7 8 9 10 11 12
El vector de fuerzas nodales está dado por:
Rx1 1
Ry1 2
0 3
0 4
0 5
-4.2 6
[f]= Rx4 7
Ry4 8
0 9
-3.6 10
0 11
0 12
V. Condiciones de borde
Los nudos 1 y 4 están restringidos, es decir: U1 = V1 = U4 = V4 = 0
Por tanto, podemos eliminar las files y columnas: 1, 2, 7 y 8, por lo que el sistema reducido queda:
VI. Sistema de ecuaciones a resolver
Cuya solución es:
VII. Estado de deformación del modelo

VIII. Tensor de Deformaciones unitarias. – Se determina según: [𝜀] = [𝐵]{𝑎}
Elemento 1
ex -1.42857 0 1.428571 0 0 0 0 -0.0012181

ey = 0 0 0 -2 0 2 0 = 0.00016687
gxy 0 -1.42857 -2 1.428571 2 0 -0.00085 -0.0002407
-0.00251
0.00082
-0.00243
Elemento 2
ex -1.42857 0 1.428571 0 0 0 -0.00085 -0.00044

ey = 0 0 0 -2 0 2 -0.00251 = 0.000483
gxy 0 -1.42857 -2 1.428571 2 0 -0.00116 -0.00061
-0.00599
0.00101
-0.00575
Elemento 3
ex 0 0 1.428571 0 -1.42857 0 0 0.001176

ey = 0 -2 0 0 0 0.5 0 = 0
gxy -2 0 0 1.428571 2 0 0.000823 -0.00347
-0.00243
0
0
Elemento 4
ex 0 0 1.428571 0 -1.42857 0 -0.000853 0.000278

ey = 0 -2 0 0 0 2 -0.002515 = 0.000167
gxy -2 0 0 1.428571 2 -1.4285 0.001018 -0.00139
-0.005751
0.000823
-0.002432
Estos valores pueden representarse en el modelo analizado, de la siguiente manera:
IX. Tensor de esfuerzos. – Determinado según: [𝜎] = [𝐷]{𝜀}
Elemento 1
sx 1866667 466666.7 0 -0.00122 -2195.84 T/m2

sy = 466666.7 1866667 0 x 0.000167 = -256.931 T/m2
txy 0 0 700000 -0.00024 -168.459 T/m2
Elemento 2
sx 1866667 466666.7 0 -0.00044 -597.3263 T/m2

sy = 466666.7 1866667 0 x 0.000483 = 695.2425 T/m2
txy 0 0 700000 -0.00061 -426.6608 T/m2
Elemento 3
sx 1866667 466666.7 0 0.001176 2195.842 T/m2

sy = 466666.7 1866667 0 x 0 = 548.9606 T/m2
txy 0 0 700000 -0.003474 -2431.542 T/m2
Elemento 4
sx 1866667 466666.7 0 0.000278 597.3251 T/m2

sy = 466666.7 1866667 0 x 0.000167 = 441.3613 T/m2
txy 0 0 700000 -0.00139 -973.3393 T/m2
Cuya representación en el modelo es:
Potosí, Mayo 14 de 2020

Practica MEF 2D Tension Plana

Cargado por

Copyright:

Formatos disponibles

Practica MEF 2D Tension Plana

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Practica MEF 2D Tension Plana

Cargado por

Copyright:

Formatos disponibles

1.

Aplicar el Método de Elementos Finitos para

Para la estructura que se muestra en la figura:

Tipo de problema: Esfuerzo plana

I. Discretización del continuo. – Adoptando elementos triangulares, tenemos

II. Matriz constitutiva. - Para todos los elementos

III. Cálculo de las matrices de rigidez

La matriz de rigidez está dada por:

Donde la matriz de deformación es:

La matriz de rigidez está dada por:

Donde la matriz de deformación es:

La matriz de rigidez está dada por:

Donde la matriz de deformación es:

La matriz de rigidez está dada por:

Ensamblando en la matriz de rigidez total, tenemos

IV. Matriz de rigidez y vector de fuerzas totales

13880 0 -8000 2800 -5880 4200 0 -7000 1

El vector de fuerzas nodales está dado por:

Los nudos 1 y 4 están restringidos, es decir: U1 = V1 = U4 = V4 = 0

VI. Sistema de ecuaciones a resolver

Cuya solución es:

VII. Estado de deformación del modelo

ex -1.42857 0 1.428571 0 0 0 0 -0.0012181

ex -1.42857 0 1.428571 0 0 0 -0.00085 -0.00044

ex 0 0 1.428571 0 -1.42857 0 0 0.001176

ex 0 0 1.428571 0 -1.42857 0 -0.000853 0.000278

IX. Tensor de esfuerzos. – Determinado según: [𝜎] = [𝐷]{𝜀}

sx 1866667 466666.7 0 -0.00122 -2195.84 T/m2

sx 1866667 466666.7 0 -0.00044 -597.3263 T/m2

sx 1866667 466666.7 0 0.001176 2195.842 T/m2

sx 1866667 466666.7 0 0.000278 597.3251 T/m2

Potosí, Mayo 14 de 2020

También podría gustarte