Mecánica del sólido deformable - Método de los elementos finit.

...

Los momentos serán cargas distribuidas con un valor distinto para cada elemento. También le añadiremos el momento puntual al nodo correspondiente.

Al final del todo copiaré el código de Máxima para que se entienda todo mucho mejor. Y a continuación muestro los resultados del ejercicio.

[pic 86]

Nodo Giro

1 0.0000e+00 0.0000e+00 0.0000e+00

2 3.5344e-01 0.0000e+00 0.0000e+00

3 7.0885e-01 0.0000e+00 0.0000e+00

4 1.0656e+00 0.0000e+00 0.0000e+00

5 1.4232e+00 0.0000e+00 0.0000e+00

6 1.7810e+00 0.0000e+00 0.0000e+00

7 2.1388e+00 0.0000e+00 0.0000e+00

8 2.4960e+00 0.0000e+00 0.0000e+00

9 2.8525e+00 0.0000e+00 0.0000e+00

10 3.2078e+00 0.0000e+00 0.0000e+00

11 3.5617e+00 0.0000e+00 0.0000e+00

12 3.2473e+00 0.0000e+00 0.0000e+00

13 2.9526e+00 0.0000e+00 0.0000e+00

14 2.6760e+00 0.0000e+00 0.0000e+00

15 2.4161e+00 0.0000e+00 0.0000e+00

16 2.1717e+00 0.0000e+00 0.0000e+00

17 1.9415e+00 0.0000e+00 0.0000e+00

18 1.7245e+00 0.0000e+00 0.0000e+00

19 1.5197e+00 0.0000e+00 0.0000e+00

20 1.3263e+00 0.0000e+00 0.0000e+00

21 1.1433e+00 0.0000e+00 0.0000e+00

22 9.7002e-01 0.0000e+00 0.0000e+00

23 8.1681e-01 0.0000e+00 0.0000e+00

24 6.7991e-01 0.0000e+00 0.0000e+00

25 5.5643e-01 0.0000e+00 0.0000e+00

26 4.4408e-01 0.0000e+00 0.0000e+00

27 3.4106e-01 0.0000e+00 0.0000e+00

28 2.4588e-01 0.0000e+00 0.0000e+00

29 1.5734e-01 0.0000e+00 0.0000e+00

30 7.4401e-02 0.0000e+00 0.0000e+00

31 1.0000e-05 0.0000e+00 0.0000e+00

Debido a nuestra forma de representar la barra en MEFI el nodo 1 de MEFI se corresponde con el 4 del enunciado, el 11 con el 3, el 21 con el 2 y el 31 con el nodo 1. En radianes tenemos:

NODO

MEFI

MEF(papel)

1

1,00E-05

1,00E-05

2

1,3430

2,06

3

3,561

4,201

4

0

0

Podemos comprobar que los datos, sin ser iguales, si que siguen una misma tendencia. Esto se puede deber a que hemos utilizado muy pocos elementos a la hora de calcularlo nosotros a papel ya que solo hemos utilizado 2 elementos. También podemos tener ciertas variaciones debidas al método que usa MEFI para resolver sus diagramas. Para reducir este error deberíamos aumentar el número de elementos que utilizamos para aproximar.

ANEXO CÓDIGO MEFI

TITULO Trabajo 2

PUNTOS

$ pun X Y

1 0.0 0.0 $nudo4

2 -0.1 0.0

3 -0.2 0.0

4 -0.3 0.0

5 -0.4 0.0

6 -0.5 0.0

7 -0.6 0.0

8 -0.7 0.0

9 -0.8 0.0

10 -0.9 0.0

11 -1.0 0.0 $nudo3

12 -1.1 0.0

13 -1.2 0.0

14 -1.3 0.0

15 -1.4 0.0

16 -1.5 0.0

17 -1.6 0.0

18 -1.7 0.0

19 -1.8 0.0

20 -1.9 0.0

21 -2.0 0.0 $nudo2

22 -2.1 0.0

23 -2.2 0.0

24 -2.3 0.0

25 -2.4 0.0

26 -2.5 0.0

27 -2.6 0.0

28 -2.7 0.0

29 -2.8 0.0

30 -2.9 0.0

31 -3.0 0.0 $nudo1

LINEAS

$ lín tipo pun

1 POL 1 2

2 POL 2 3

3 POL 3 4

4 POL 4 5

5 POL 5 6

6 POL 6 7

7 POL 7 8

8 POL 8 9

9 POL 9 10

10 POL