Fem Portico Plano E1

2.
1 Nudos
Las coordenadas de los nudos considerados en la estructura
columna 1: coordenada "x" del nudo en el sistema global
columna 2: coordenada "y" del nudo en el sistema global
NODE :=
1
1
2
240
180
480
180
2.2 Propiedades de los elementos

Las propiedades de la seccin trasversal de los elementos y del material que est compuesto, cada fila
representa una propiedad distinta y las columnas son:
columna 1: rea de la seccin trasversal del elemento
columna 2: momento de inercia
columna 3: mdulo de elasticidad del material
PROP :=
1
1
2
12
3
600
2.9104
2.3 Elementos
Cada fila representa una barra, contiene la informacin de la conectividad del elemento en el sistema, cada
columna representa:
columna 1: nudo inicial del elemento
columna 2: nudo final del elemento
columna 3: nmero de propiedad del elemento
MEMB :=
1
2.4 Restricciones - Apoyos

Cada fila representa un apoyo de la estructura, las columnas informan el comportamiento para cada grado de
libertad, la convencin es:
"1" para los grados de libertad de desplazamiento restringido
"0" para los grados de libertad donde existe desplazamiento libre.
Cada columna representa:
columna 1: nmero del nudo donde existe el apoyo
columna 2: "ux?" informacin del desplazamiento en la direccin "x" global
columna 3: "uy?" informacin del desplazamiento en la direccin "y" global
columna 4: "Rz?" iinformacin de la rotacin en el eje "z" global
SUPP :=
1
2.5 Cargas
2.5.1 Carga distribuda en la direccin y'
Se debe ingresar considerando la orientacin local del sistema, cada columna representa:
columna 1: nmero del elemento
columna 2: carga distribuida sobre el elemento
mly :=
1
2
1
Obteniendo la matriz de rigidez respecto al sistema global para el elemento:
m := 1
3.1 Matriz de rigidez respecto al sistema local

rea de la seccin transversal del elemento
A(m) = 12
Momento de inercia respecto al eje z
I(m) = 600
Mdulo de elasticidad del material
E(m) = 2.9 10
Longitud del elemento
(xf (m)
x i (m)
+ y f (m) y i (m)
( 240 0) + ( 180 0) = 300
=
l(m) = 300
reemplazando estos valores, resuta:
k e(m) =
1.16 10
7.733
1.16 10
1.16 10
2.32 10
5
1.16
1.16 10
0
1.16
1.16
7.733
1.16
0
10
103
10
3
105
3
3
1.16 10
1.16 10
7.733
1.16
1.16 10
1.16 10
7.733
1.16 103
2.32 10
10
3.2 matriz de transformacin local - global

Coseno del menor ngulo que forma el elemento respecto al eje global "x"
x (m) x (m) =
l(m)
( 240 0)
300
= 0.
x (m) = 0.8
8
Coseno del menor ngulo que forma el elemento respecto al eje global "y"
y (m) y (m) =
l(m)
( 180 0)
300
= 0.
6y (m) = 0.6
ordenando los resultados, se tiene:
0.8
0.6
0
0
0
0.6 0
0.8 0
0 0.8
0 0.6 0.8 0
0.6 0
3.3 Matriz de rigidez respecto al sistema global

La matriz de rigidez respecto al sistema global est dado por.
T
Ke := T(m) k e(m) T(m)
745.184
696
553.088
422.549
553.088
745.184 553.088
928
553.088 422.549
2.32 10
928
696
928
928
696
696
1.16 10
745.184
Ke = 553.088
553.088
696
745.184
553.088
696
422.549
928
553.088
422.549
928
gdl(m) =
928
928
696
1.16 10
2.32 10
696
3.4 Matriz de rigidez total ensamblado

Seguidamente se ensambla en una sola matriz segn los grados de libertad, lo que resulta.
2.195
553.088
696
696
745.184 553.088
1.45 10
10
553.088 422.549
553.088
696
884.5
5.22 10
928
553.088
928
1.16 105
0
745.184
928
696
1.16 10
696
15.104
1.813 10
1.813
103
0
1.45 10
696
696
553.08
884.5
696
745.184
437.654
0
0
0
2.32 10
5
928
1.45 103
1.45 103
15.104
10
1.813 10
0
15.104
1.813 10
103
1.813
1.45 105
103
1.813
2.9
4. VECTOR DE FUERZAS NODALES EQUIVALENTES
4.1 Obteniendo vector de fuerzas nodales equivalentes para el elemento
m := 2
Todo el anlisis ser para la carga distribuida sobre el elemento m = 2

la carga distribuida es:
w := mly
1,2
w = 0.25
4.1.1 Momentos de empotramiento perfecto en el sistema de orientacin local( la carga acta en el eje local y')
0 30
3
MEPmly(m , w)T =
0 30 1.2
1.2 10
10
4.1.2 Matriz de transformacin de desplazamientos, para el elemento m = 2
1
0
0
0
0
0
est dado por
0 0 0 0 0
1 0 0 0 0
0 1 0 0 0
0 0 1 0 0
0 0 0 1 0
0 0 0 0
4.1.3 Mometnos de empotramiento perfecto en el sistema de orientacin global, est dado por:
0
0
30
T(m ) MEPmly(m , w) =
1.2 10
MEPmlY(m , w) =
30
3
1.2 10
30
1.2
30
1.2
3
10
10
4.1.4 Carga nodal equivalente en el sistema de orientacin global, est dado por:
0
0
30
1.2 103
30
1.2 10
MEPmlY(m , w) =
CNEmlY(m , w) =
30
1.2 10
1.2 10
los grados de libertad asociado para la carga nodal equivalentes es;

0
30
2
1.2 10
CNEmlY(m , w) =
30
gdl(m) =
30
7
8

1.2 10
fuerzas nodales equivalentes para cada carga, segn los grados de libertad asociado para cada nudo en la
que actan, resultando.
0
30
3
1.2 10
0
cne
F
=
0
30
1.2 10
4.2 Vector de cargas puntuales en nudos

T
FNLF = ( 0 0 0 0 0 0 0 0 0 )
no existe cargas puntuales en nudos
4.3 Finalmente, el vector de carga nodal equivalente resulta
30
1.2
103
0
30
3
1.2 10
NLF
+ F
cne
0
0
0
0
30
1.2 10
F=
0
0
0
0
30
1.2 10
5 DESPLAZAMIENTOS EN NUDOS
5.1 Imponiendo las condiciones de contorno
La matriz de rigidez total(K) de la estructura obtenido anteriormente est libre en todos sus grados de
libertad(desplazamiento libre), pero como la estructura cuenta con apoyos en las cuales los desplazameintos
son nulos se modifica "K" para tal efecto.
696
696
3
2.195 103 553.088
745.184 553.088
1.45 10
Km =
437.654
553.088
696
745.184
553.088
696
884.5
553.088 422.549
0
884.5
5.22 105
553.088
928
422.549
928
696
696
553.088
1.16 10
928
696
1 1 307
1.16 105
0
5
928
1 1 307
553.088
928
928
1 0 1 307
696
1.813 103 1.45 10
0
0
0
0
0
0
15.104
1.813 10
1.45 10
5.2 Matriz aumentada
10
15.104 1.813 10
10
0
307
307
0
1.813 10
1 10
3
0
1.813
1.45 10
0
10
1.813
307
1 10
Teniendo los coeficientes de las incgnitas(Km) y los trminos independientes(F) de las ecuaciones
simultneas que se forman para cada grado de libertad en los nudos, existe muchas maneras para resolver la
ecuacin. Formando la matriz aumentada
augment(Km, F)
2.195 103 553.088
437.654
553.088
884.5
696
553.088
745.184
=
422.549
553.088
928
696
745.184
696
553.088
884.5
553.088
422.549
5.22 10
696
928
553.088
1 1 307
696
928
553.088
1.16 105
696
928
928
0
5
696
928
0
15.104
1.813 10
1.813 10
1.16 10
1 1 307
1.45 103
696
3 1.45 10
0
0
0
1 0 1 307
30
3
1.2 10
1.45 10
15.104 1.813 10
10
0
1.813
30
3
1.2 10
307
1 10
1.813
10
1.45 10
307
10
0
10
1.813
1 10
5.3 Resolviendo la ecuacicn

la ecuacin anterior en su forma escalonada reducida
00
rref (augment(Km, F)) =
00
00
0 0 0 0 0 0 0 0
0.025
0.095
1 0 0 0 0 0 0 0
0 0 0 0 0 0
3
0 1
2.17 10
0 0 1 0 0 0 0 0
0 0 0 1 0 0 0 0
0 0 0 0 1 0 0 0
0
0
0 0 0 0 0 1 0 0
0 0 0 0 0 0 1 0
0 0 0 0 0 0 0 1
5.4 Los desplazamientos

Ordenando
T
D =
0.095
2.17 10
0 0 0 0 0 0
0.024727 0.095411
0.00217
0.025
Do =
6 REACCIONES EN LOS APOYOS
2.195
10
553.088
696
745.184 553.088
696
1.45 103
437.654
884.5
553.088 422.549
928
15.104
553.088
553.088
696
745.184
553.088
696
0
1.813
3
0
3
10
0.025
0.095
30
K=
553.088
696
422.549
928
553.088
422.549
928
0
5
928
1.16 10
696
928
2.32 10
D=
1.45 10
1.45 10
F=
30
3
15.104 1.813 10
0
15.104
1.813 10
3
1.813 10
01.45
0 0
1.813
10 10
2.9 10
1.2 10
1.554 10 14
13
2.274 10
35.855
R=
24.625
145.986
35.855
Ordenado :
nudo
Ro =
1.688 10
Ry
35.854609
24.625498
145.98617
35.375 3
35.85460
9
Mz
Rx
35.374502
1687.604162
Obteniendo los resultados para el elemento
6.1 Matriz de rigidez respecto al sistema local
k e(m) =
15.104
15.104
1.45 10
0
3 1.45 10
1.813
0
10
1.813
10
0
1.45 10
2.9 10
1.813
1.45 10
1.813 10
1.45 10
0
10
10
0
15.104
0
1.813 10
m := 2
1.813 103
5
1.45 10
15.104
1.813 10
1.813
matriz de transformacin de desplazamientos
1
0
0
0
0
0
0 0 0 0 0
1 0 0 0 0
0 1 0 0 0
0 0 1 0 0
0 0 0 1 0
0 0 0 0
vector de desplazamiento en los nudos del elemento en el sistema global

0.025
0.095
De(m) =
2.17 10
3
gdl(m) =
7
8
desplazamiento en el sistema local
T(m) De(m) =
0.095
3
2.17 10
0.095
de(m) =
2.17 10
2.9 10
6.2 Desplazamientos en los nudos del elemento en el sistema local
0
0
6.3 Vector del momento de empotramiento perfecto

En el sistema de orientacin local, para todas las cargas que actan en el elemento.
MEPsl(m) =
30
1.2 103
0
30
1.2 10
6.4 Esfuerzos en los extremos del elemento
k e(m) de(m) + MEPsl(m)
35.855
24.625
397.724
= 35.855
35.375
1.688 10
............... de igual manera para cada elemento
Resumiendo para todos los elementos, se tiene:

1
MFOR =
Donde:
columna
columna
columna
columna
columna
1:
2:
3:
4:
5:
43.45898
-6
43.458986
35.85460
-9
35.854609
nmero de elemento
nmero de nudo
exfeurzo axial
esfuerzo cortante
momento flector
4
1.812367
1.81236
7
24.62549
8
35.37450
2
5
145.98617
397.7238
397.723
8
1687.60416
7.1 Grfico del desplazamiento en el eje "x" local

Se evaluar la variacin del desplazamiento de la directriz del elemento "m", para 11 puntos
m := 2
l(m) .. l(m)
x := 0 , 10
x =
ux(m , x) =
0
24
48
72
96
120
144
168
192
216
240
0.02472
7
0.02225
5
0.01978
2
0.01730
9
0.01483
6
0.01236
4
0.00989
1
0.00741
8
0.00494
5
0.00247
3
0
7.2 Grfico del desplazamiento en el eje "y" local
x =
uy(m , x) =
0
24
48
72
96
120
144
168
0.095411
0.134927
0.152155
0.151365
0.136827
- 0.11281
0.083585
-
216
0.053421
0.026589
-
240
0.007359
0
192

Fem Portico Plano E1

Cargado por

Copyright:

Formatos disponibles

Fem Portico Plano E1

Cargado por

Información del documento

Descripción original:

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

Fem Portico Plano E1

Cargado por

Copyright:

Formatos disponibles

2.

2.2 Propiedades de los elementos

2.4 Restricciones - Apoyos

Obteniendo la matriz de rigidez respecto al sistema global para el elemento:

3.1 Matriz de rigidez respecto al sistema local

Longitud del elemento

( 240 0) + ( 180 0) = 300

3.2 matriz de transformacin local - global

3.3 Matriz de rigidez respecto al sistema global

Ke := T(m) k e(m) T(m)

3.4 Matriz de rigidez total ensamblado

4. VECTOR DE FUERZAS NODALES EQUIVALENTES

4.1 Obteniendo vector de fuerzas nodales equivalentes para el elemento

Todo el anlisis ser para la carga distribuida sobre el elemento m = 2

4.1.2 Matriz de transformacin de desplazamientos, para el elemento m = 2

est dado por

los grados de libertad asociado para la carga nodal equivalentes es;

4.2 Vector de cargas puntuales en nudos

4.3 Finalmente, el vector de carga nodal equivalente resulta

1.813 103 1.45 10

5.2 Matriz aumentada

2.195 103 553.088

5.3 Resolviendo la ecuacicn

rref (augment(Km, F)) =

5.4 Los desplazamientos

6 REACCIONES EN LOS APOYOS

Obteniendo los resultados para el elemento

6.1 Matriz de rigidez respecto al sistema local

matriz de transformacin de desplazamientos

vector de desplazamiento en los nudos del elemento en el sistema global

desplazamiento en el sistema local

6.2 Desplazamientos en los nudos del elemento en el sistema local

6.3 Vector del momento de empotramiento perfecto

6.4 Esfuerzos en los extremos del elemento

k e(m) de(m) + MEPsl(m)

............... de igual manera para cada elemento

Resumiendo para todos los elementos, se tiene:

7.1 Grfico del desplazamiento en el eje "x" local

7.2 Grfico del desplazamiento en el eje "y" local

También podría gustarte