Marco 3d Metodo Matricial
Marco 3d Metodo Matricial
Marco 3d Metodo Matricial
Facultad De Ingeniera
Programa De Maestra Y Doctorado En Ingeniera
Teora General de las Estructuras
Marco 3D
Obtener las fuerzas internas en un marco tridimensional como el que se muestra a continuacin, empleando el mtodo matricial. Proponga dimensiones
diferentes a las mostradas en la figura, tanto para la geometra general como para las secciones transversales de los elementos. Puede, incluso, proponer
otro modelo estructural, "ms interesante". Considere una condicin de carga lateral (sobre nudos) y lineal uniforme sobre las vigas. Adems:
a) Represente los resultados mediante vectores como se ilustra al final del ejemplo que se incluye en los apuntes Marco 3D.
Columnas: Vigas:
NB := 13
bc := 0.55 bv := 0.3
hc := 0.3 hv := 0.4
Longitudes:
Lc := 3 Lv7 := 7 Lv9 := 5
3
bv hv
I0 := = 0.0009
12
Se propone lo siguiente:
Ejes locales:
En la siguiente imagen se muestra la definicin de los ejes locales para cada uno de los elementos del marco 3D adems se muestra el eje global que se
utilizara en los clculos:
Calculo de momentos de inercia de las barras:
I :=
4 4
Ba rra B (m) h (m) Iy (m ) Iz (m )
1 0.55 0.3 0.00416 0.00124
2 0.55 0.3 0.00416 0.00124
3 0.55 0.3 0.00416 0.00124
4 0.55 0.3 0.00416 0.00124
5 0.55 0.3 0.00416 0.00124
6 0.55 0.3 0.00416 0.00124
7 0.3 0.4 0.00090 0.00160 I=
8 0.3 0.4 0.00090 0.00160
9 0.3 0.4 0.00160 0.00090
10 0.3 0.4 0.00160 0.00090
11 0.3 0.4 0.00160 0.00090
12 0.3 0.4 0.00090 0.00160
13 0.3 0.4 0.00090 0.00160
F I :=
Ba rra s
Fuerzas
1 2 3 4 5 6 7 8 9 10 11 12 13
FAX 0 0 0 0 0 0 0 0 0 0 0 0 0
FAY 0 0 0 0 0 0 24.5 12.25 0 0 0 24.5 12.25
FAZ 0 0 0 0 0 0 0 0 12.5 25 12.5 0 0
MAX 0 0 0 0 0 0 0 0 0 0 0 0 0
MAY 0 0 0 0 0 0 0 0 -15.625 -31.25 -15.625 0 0
MAZ 0 0 0 0 0 0 28.583 14.292 0 0 0 28.583 14.292
FBX 0 0 0 0 0 0 0 0 0 0 0 0 0
FBY 0 0 0 0 0 0 24.5 12.25 0 0 0 24.5 12.25
FBZ 0 0 0 0 0 0 0 0 12.5 25 12.5 0 0
MBX 0 0 0 0 0 0 0 0 0 0 0 0 0
MBY 0 0 0 0 0 0 0 0 15.625 31.25 15.625 0 0
MBZ 0 0 0 0 0 0 -28.583 -14.292 0 0 0 -28.583 -14.292
F fij := ( 0 0 37 15.625 - 28.583 0 0 0 61.75 31.25 14.291 0 0 0 24.75 15.625 14.292 0 0 0 37 - 15.625 - 28.583 0 0 0 61.75
Esquema de momentos:
UNIVERSIDAD NACIONAL AUTNOMA DE MXICO
PROGRAMA DE MAESTRA Y DOCTORADO EN INGENIERA
TEORIA GENERAL DE LAS ESTRUCTURAS
MARCO 3D
Ademas de las cargas que se aplican en las barras se contemplan fuerzas horizontales sobre los siguientes nudos
Nudo 1 : Nudo 3 :
Fx = 50 Fy = 80
Fy = 100
T
F ext := ( 50 100 0 0 0 0 0 0 0 0 0 0 0 80 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 )
F := - F fij + Fext
T
F = ( 50 100 - 37 - 15.625 28.583 0 0 0 - 61.75 - 31.25 - 14.291 0 0 80 - 24.75 -15.625 - 14.292 0 0 0 - 37 15.625 28.583 0
ux'2 := ux'1 uy'2 := uy'1 uz'2 := uz'1 ux'8 := ux'7 uy'8 := uy'7 uz'8 := uz'7
ux'3 := ux'1 uy'3 := uy'1 uz'3 := uz'1 ux'12 := ux'7 uy'12 := uy'7 uz'12 := uz'7
ux'4 := ux'1 uy'4 := uy'1 uz'4 := uz'1 ux'13 := ux'7 uy'13 := uy'7 uz'13 := uz'7
Barras 9, 10 y 11
0 1 0
ux'9 := - 1 uy'9 := 0 uz'9 := 0 ux'10 := ux'9 uy'10 := uy'9 uz'10 := uz'9
0 0 1
ux'11 := ux'9 uy'11 := uy'9 uz'11 := uz'9
UNIVERSIDAD NACIONAL AUTNOMA DE MXICO
PROGRAMA DE MAESTRA Y DOCTORADO EN INGENIERA
TEORIA GENERAL DE LAS ESTRUCTURAS
MARCO 3D
( (
A7_8 := Ensamblaje A ux'7 , uy'7 , uz'7 , Lv7 )) 0 0.143 0 0 0 1 0 - 0.143 0 0 0 0
Construyendo las submatrices de A: 0 0.286 0 0 0 1 0 - 0.286 0 0 0 1
0 0.143 0 0 0 0 0 - 0.143 0 0 0 1
(
A7A := submatrix A7_8 , 1 , 8 , 1 , 6 )
0 0 0.143 0 -1 0 0 0 -0.143 0 0
0
A7_8 =
A7B := submatrix( A7_8 , 1 , 8 , 7 , 12) 0 0 0.286 0 -1 0 0 0 -0.286 0 - 1 0
Las columnas restantes tienen las mismas caractersticas: 0 0 0.143 0 0 0 0 0 -0.143 0 - 1 0
- 1 0 0 0 0 0 1 0 0 0 0 0
A8A := A7A A12A := A7A A13A := A7A
0 0 0 -1 0 0 0 0 0 1 0 0
A8B := A7B A12B := A7B A13B := A7B
1B
A M0 M0 M0 M0 M0
M A2B M0 M0 M0 M0
0
M0 M0 A3B M0 M0 M0
M0 M0 M0 A4B M0 M0
M0 M0 M0 M0 A5B M0
M M0 M0 M0 M0 A6B
0
A := A7A A7B M0 M0 M0 M0 A := Ensamblaje( A )
M0 A8A A8B M0 M0 M0
A9B M0 M0 A9A M0 M0
M0 A10B M0 M0 A10A M0
M M0 A11B M0 M0 A11A
0
M M0 M0 A12A A12B M0
0
M0 M0 M0 M0 A13A A13B
UNIVERSIDAD NACIONAL AUTNOMA DE MXICO
PROGRAMA DE MAESTRA Y DOCTORADO EN INGENIERA
TEORIA GENERAL DE LAS ESTRUCTURAS
MARCO 3D
1 2 3 4 5 6 7 8
1 0.333 0 0 0 0 0 0 0
2 0.667 0 0 0 -1 0 0 0
3 0.333 0 0 0 -1 0 0 0
4 0 0.333 0 0 0 0 0 0
A= 5 0 0.667 0 1 0 0 0 0
6 0 0.333 0 1 0 0 0 0
7 0 0 1 0 0 0 0 0
8 0 0 0 0 0 1 0 0
9 0 0 0 0 0 0 0.333 0
10 0 0 0 0 0 0 0.667 ...
0.004159 0.001238
(
k1 := k E , I1 , y' , I1 , z' , Lc , Areac , G , Jc )
0.004159 0.001238
k2 := k1 k3 := k1 k4 := k1 k5 := k1 k6 := k1 0.004159 0.001238
0.004159 0.001238
(
k7 := k E , I7 , y' , I7 , z' , Lv7 , Areav , G , Jv ) 0.004159 0.001238
0.004159 0.001238
k8 := k7 k912:=:=.7k7
k k13 := k7
I = 0.0009 0.0016
0.0009 0.0016
(
k9 := k E , I9 , y' , I9 , z' , Lv9 , Areav , G , Jv )
0.0016 0.0009
0.0016 0.0009
k10 := k9 k11 := k9
0.0016 0.0009
0.0009 0.0016
0.002773 0.000257 0.00064
0.0009 0.0016
0.002773 0.000257 0.00064
0.002773 0.000257 0.00064
Matriz diagonal de rigidez para toda la estructura:
0.000825 0.000457 0.00036
k1 = k7 = k9 =
0.000825 0.000457 0.00036 ( ( ))
kd := diag stack k1 , k2 , k3 , k4 , k5 , k6 , k7 , k8 , k9 , k10 , k11 , k12 , k13
0.000825 0.000457 0.00036
0.055 0.017143 0.024
0.00024 0.000077 0.000108
UNIVERSIDAD NACIONAL AUTNOMA DE MXICO
PROGRAMA DE MAESTRA Y DOCTORADO EN INGENIERA
TEORIA GENERAL DE LAS ESTRUCTURAS
MARCO 3D
1 2 3 4 5 6 7 8
1 0.019 0 0 0 -0.003 -0 -0.017 0
2 0 0.025 0 0.001 0 0 0 -0
3 0 0 0.055 0 -0 0 0 0
4 0 0.001 0 0.003 0 0 0 0
5 -0.003 0 -0 0 0.007 0 0 0
6 -0 0 0 0 0 0.001 0 -0
7 -0.017 0 0 0 0 0 0.036 0
K= 8 0 -0 0 0 0 -0 0 0.025
9 0 0 -0 0 0 0 0 0
10 0 0 0 -0 0 0 0 0.001
11 0 0 -0 0 0 0 -0.003 0
12 0 0 0 0 0 0 -0 0
13 0 0 0 0 0 0 -0.017 0
14 0 0 0 0 0 0 0 -0
15 0 0 0 0 0 0 0 0
16 0 0 0 0 0 0 0 ...
A1 := A1_6
0
0
0 M1y'
7246.049 20.093
0 M2y'
0 1121.79 3.111
- 6124.259 M3y' - 16.982
0
43575.689 M1z' := diag k e = 35.95
e1 := A1 d
1
=
50212.869 M ( 1) 1 41.426
2z'
d2
d 6637.18 M3z' 5.476
3 - 202.027 - 11.111
d4 N
- 4700.514 M -1.128
T
d5
d6
Fuerzas reales:
Momentos flexionantes en "y" y "z": Momento Torsionante:
M Ay + MBy M Az + MBz
F Az := = 3.111 F Ay := = 41.426 N = - 11.111 F Ax := N F Bx := N
L0 L0
1 1
UNIVERSIDAD NACIONAL AUTNOMA DE MXICO
PROGRAMA DE MAESTRA Y DOCTORADO EN INGENIERA
TEORIA GENERAL DE LAS ESTRUCTURAS
MARCO 3D
F Bz := - F Az F By := -F Ay
F
Ax
FAy
Agrupando las fuerzas de las barras en un vector:
FAz
MAx
MAy
M
F Az
A F 1 :=
Fi = FBx
FB
FBy
FBz
MBx
M
By
MBz
Sistematizando el problema:
F II := for i 1 .. NB
Fz1 , i P ( 8 i) - 1
P ( 8i) - 4 + 2 P( 8 i) - 3 + P ( 8 i) - 2
Fz2 , i
L0
i
P ( 8 i) - 7 + 2 P ( 8 i) - 6 + P ( 8 i) - 5
Fz3 , i
L0
i
Fz4 , i P 8 i
Fz5 , i P ( 8 i) - 7 + P ( 8 i) - 6
Fz6 , i P ( 8 i) - 4 + P ( 8 i) - 3
Fz7 , i P ( 8 i) - 1
Fz8 , i - Fz2 , i
Fz9 , i - Fz3 , i
Fz10 , i P8 i
Fz11 , i P( 8 i) - 6 + P ( 8 i) - 5
Fz12 , i P( 8 i) - 3 + P ( 8 i) - 2
x augment( Fz)
x
1 2 3 4 5 6
1 -11.111 -59.646 -8.336 -58.642 -66.745 -42.52
2 41.426 -5.819 32.466 52.427 15.924 43.576
3 3.111 24.27 19.285 -10.962 7.642 6.654
4 -1.128 0.029 1.762 -1.085 0.034 1.726
5 23.203 42.03 36.097 -8.249 9.57 8.131
F II = 6 77.376 -1.281 61.604 87.741 20.336 72.193
7 -11.111 -59.646 -8.336 -58.642 -66.745 -42.52
8 -41.426 5.819 -32.466 -52.427 -15.924 -43.576
9 -3.111 -24.27 -19.285 10.962 -7.642 -6.654
10 -1.128 0.029 1.762 -1.085 0.034 1.726
11 -13.871 30.78 21.759 -24.638 13.357 11.831
12 46.901 -16.175 35.795 69.54 27.435 ...
SOLUCION: ESTADO I + ESTADO 2 Las fuerzas en las barras son las columnas de la matriz
Sol := F I + F II
Ba rra s
Fuerzas
1 2 3 4 5 6 7 8 9 10 11 12 13
FAX -11.111 -59.646 -8.3361 -58.642 -66.745 -42.52 -43.273 -17.58 -55.482 -10.905 -45.539 -14.579 -8.3596
FAY 41.4256 -5.8189 32.4664 52.4269 15.9237 43.5764 21.1351 10.39 -3.6166 -1.4229 1.70562 23.6179 12.7544
FAZ 3.11063 24.27 19.2855 -10.962 7.64225 6.65394 3.0919 -1.9945 -10.024 21.3913 -5.774 3.05563 -1.9628
MAX -1.1281 0.02889 1.7617 -1.0853 0.03448 1.72615 1.45685 -1.0207 -0.8057 -0.5856 -0.7026 2.36572 -1.9407
MAY 23.2033 42.0295 36.0971 -8.2488 9.56966 8.13063 10.2018 -7.909 -67.174 -31.741 -56.596 10.0948 -7.7761
MAZ 77.3756 -1.2813 61.6044 87.7411 20.3363 72.1925 14.6771 8.03646 -9.0095 -3.553 4.23739 23.8324 16.0646
FBX -11.111 -59.646 -8.3361 -58.642 -66.745 -42.52 -43.273 -17.58 -55.482 -10.905 -45.539 -14.579 -8.3596
FBY -41.426 5.81892 -32.466 -52.427 -15.924 -43.576 27.8649 14.11 3.61663 1.42288 -1.7056 25.3821 11.7456
FBZ -3.1106 -24.27 -19.285 10.9623 -7.6422 -6.6539 -3.0919 1.99448 35.0236 28.6087 30.774 -3.0556 1.9628
MBX -1.1281 0.02889 1.7617 -1.0853 0.03448 1.72615 1.45685 -1.0207 -0.8057 -0.5856 -0.7026 2.36572 -1.9407
MBY -13.871 30.7805 21.7593 -24.638 13.3571 11.8312 11.4415 -6.0524 -45.444 13.6979 -34.774 11.2946 -5.9635
MBZ 46.9013 -16.175 35.7948 69.5395 27.4348 58.5365 -38.231 -21.057 -9.0737 -3.5614 4.2907 -30.007 -12.534
Sol
UNIVERSIDAD NACIONAL AUTNOMA DE MXICO
PROGRAMA DE MAESTRA Y DOCTORADO EN INGENIERA
TEORIA GENERAL DE LAS ESTRUCTURAS
MARCO 3D
8.24
Barra 7 Barra 8
y' y'
21.13 3.09 27.86 1.99 10.39 14.11
10.20 11.44
38.23 21.05
x' 17.58 x'
1.45 43.27 3.09 1.45 1.02 17.58 1.02
43.27 11.44
7.9
14.67 8.03
Barra 9 Barra 10
z' z'
35.02 3.61 28.60 1.42
21.39
45.44 67.14 31.74
x' 55.48 x' 10.90
0.80 55.48 3.61 0.80 0.58 10.90 1.42 0.58
10.02
9.07 3.56 13.69
Barra 11
z'
1.70 30.77
4.23
4.29
34.77 56.59
x' 45.53
0.70 45.53 0.70
5.77
1.70
y'
Barra 12 Barra 13
y' y'
23.61 3.05 25.38 1.96 12.75 11.74
10.09 11.29
30.00 12.53
x' 8.35 x'
2.36 14.57 3.05 2.36 1.94 8.35 1.94
14.57 1.96
7.77
23.83 16.06
-14.58
VIGA
-55.49
-5 VIG
COLUMNA
-43.28 8.
63
A
COL11U.
VI
-10.91
GA
-
-8.36
A
VIG
1M
-6
COLUMNA
3 NA
6.7
4 VIG
A
-17.59
Z
COLU
-5
9.
65
Y
MNA
VI
GA
54
-42
NA
-45.
X .51
UM
GA
VI
COL
NA
-8
M
.3
5
COLU
-23.62
VIGA
-35.01
VIG
COLUMNA
-1 A
25.38
0.
96
-21.14
COL3U.1
VI
GA
-12.75
21.4
-28.6
27.86
A
VIG
M
COLUMNA
1 NA
7.6 VIG
4 A
-10.39
5
11.7
COLU2M Z
Y
4. NA
27
VI
GA
76
NA
-30.
X
GA
UM
6.6
VI 5
14.11
COL
NAM
COLU
19
.2
9
-3
. 61 5
3.0
VIGA
VIG
COLUMNA
43
A
52.
3.09
-1
C O LU
VI .42
41.
43 GA 96
A -1.
VIG
MNA
COLUMNA
.92 VIG
15 A
99
COLUMNA -1. Z
83 Y
-5. 1.7
VI
GA
.5NA
X
8
GA
U43M
VI
COL
32MNA
7
.4
COLU
-9 8
. .501.0
-691
VIGA
VIG
COLUMNA
9. A
69
160.8.1 06
-4 .28
-11-3
C O LU
VI .55
GA 79 .43 .77
87.
A -27 -7
VIG
MNA
COLUMNA
.43 3.
4 2 -11 56 VIG
77. A
.2 .9
16 -7 .34
Z 20 95
COLUMNA 5.
Y 4.2
1
.5
3
58
VI
NA -
GA
29
-1. X
GA
UM
VI
COL
0 4
6. -4
.2
8
23.
72
76
5.
NA -3
M
COLU
4
.6
61
-67.14
24
.64
VIGA
VIG
COLUMNA
15.95
A
45.4
-30.
-31.74
-16.06
C O LU
VI
17.23 -8
.2 -13
GA 5
.36
-38.23
A
VIG
MN2A3.2
COLUMNA
39.77
7.14
VIG
A
1
9.
-3 CO Z
57
0. LUMNA04
77
7
5
-56.
Y
7.38
-11
VI .83
42
GA
.
NA
X
GA
UM
VI
-21.05
COL
8.1
3
4
34 . 7
-2
1 .7
6
NAM
COLU
36
.1
1
2.37
-0.81
VIGA
VIG
COLUMNA
1.46 -1 A
.0
8
COL1.U
VI
GA
-0.59
-
A
13M
VIG
-1.94
COLUMNA
NA
0.0 VIG
3 A
-1.02
COLU
Y
0.
0M
3 NA
VI
-0.7
GA
NA
X
GA
UM
VI 1.7
COL
3
NAM
COLU
1.
76