Práctica de Elementos Finitos
Práctica de Elementos Finitos
Práctica de Elementos Finitos
Mdulo:
Elementos Finitos
Docente:
Fecha:
Sucre Bolivia
SET 1 y 2
r3=0
si
Estticamente
determinada
externamente
m = 2j - r
no
r3<0
r3>0
si
Estticamente
determinada
internamente
si
Inestable o
hiposttica
no
si
Inestable o
hiposttica
m < 2j r
no
no
si
no
no
m > 2j r
si
Estticamente
indeterminad
a
externamente
cuyo grado
ser la
diferencia
entre ambos
trminos de la
inecuacin
Estticamente
indeterminad
a
internamente
cuyo grado
ser la
diferencia
entre ambos
trminos de la
inecuacin
FIN
FIN
)|
)|
(
|
Donde:
-
) |
) |
) |
) |
3. Indicar los pasos para pasar de la ecuacin 5.11 a la ecuacin 5.12, identificando las
partes correspondientes a las cargas distribuidas, puntuales y momentos.
) |
) |
Ya que:
en los apoyos
Pm = salto en fuerza cortante en xm
-Mk = salto en momento flector en xk
Sustituyendo:
H1
1,00
0,99
0,97
0,94
0,90
0,84
0,78
0,72
0,65
0,57
0,50
0,43
0,35
0,28
0,22
0,16
0,10
0,06
0,03
0,01
0,00
Grfica de la funcin 1:
-1,00
-0,90
-0,80
-0,70
-0,60
-0,50
-0,40
-0,30
-0,20
-0,10
0,00
0,10
0,20
0,30
0,40
0,50
0,60
0,70
0,80
0,90
1,00
1.20
1.00
VALORES DE H1
0.80
0.60
0.40
0.20
0.00
-1.50
-1.00
-0.50
0.00
VALORES DE
0.50
1.00
1.50
H2
0,00
0,09
0,16
0,22
0,26
0,28
0,29
0,30
0,29
0,27
0,25
0,22
0,19
0,16
0,13
0,09
0,06
0,04
0,02
0,00
0,00
Grfica de la funcin 2:
-1,00
-0,90
-0,80
-0,70
-0,60
-0,50
-0,40
-0,30
-0,20
-0,10
0,00
0,10
0,20
0,30
0,40
0,50
0,60
0,70
0,80
0,90
1,00
0.35
0.30
0.25
VALORES DE H2
0.20
0.15
0.10
0.05
0.00
-1.50
-1.00
-0.50
0.00
VALORES DE
0.50
1.00
1.50
H3
0,00
0,01
0,03
0,06
0,10
0,16
0,22
0,28
0,35
0,43
0,50
0,57
0,65
0,72
0,78
0,84
0,90
0,94
0,97
0,99
1,00
Grfica de la funcin 3:
-1,00
-0,90
-0,80
-0,70
-0,60
-0,50
-0,40
-0,30
-0,20
-0,10
0,00
0,10
0,20
0,30
0,40
0,50
0,60
0,70
0,80
0,90
1,00
1.20
1.00
0.80
VALORES DE H3
0.60
0.40
0.20
0.00
-1.50
-1.00
-0.50
0.00
VALORES DE
0.50
1.00
1.50
H4
0,00
0,00
-0,02
-0,04
-0,06
-0,09
-0,13
-0,16
-0,19
-0,22
-0,25
-0,27
-0,29
-0,30
-0,29
-0,28
-0,26
-0,22
-0,16
-0,09
0,00
Grfica de la funcin 4:
-1,00
-0,90
-0,80
-0,70
-0,60
-0,50
-0,40
-0,30
-0,20
-0,10
0,00
0,10
0,20
0,30
0,40
0,50
0,60
0,70
0,80
0,90
1,00
0.00
-1.50
-1.00
-0.50
-0.05
0.00
0.50
1.00
1.50
-0.10
VALORES DE H4
-0.15
-0.20
-0.25
-0.30
-0.35
VALORES DE
) (
Desarrollando se tiene:
)
(
Reemplazando 5.25 y
(
(
)
(
(
[
(
(
)
)
)
)
(
(
)
(
(
)
)
Dnde:
SET 3
Datos:
E = 2x1012 [KN/m2] = 200000 [KN/cm2]
P = 300 [KN]
A1= 270 [mm2] = 2.7 [cm2]
A2= 270 [mm2] = 2.7 [cm2]
A3= 450 [mm2] = 4.5 [cm2]
Solucin
a)
Tabla de conectividad
ELEMENTO N. I. N. F. A (cm2) L (cm)
1
a
b
2,7
17
2
b
c
2,7
17
3
c
d
4,5
34
a
b
1 -1
31764,706 -31764,706 a
=
K1= 31764,706
-1 1
-31764,706 31764,706 b
b
c
1 -1
31764,706 -31764,706 b
=
K2= 31764,706
-1 1
-31764,706 31764,706 c
c
d
1 -1
26470,588 -26470,588 c
=
K3= 26470,588
-1 1
-26470,588 26470,588 d
La matriz de rigidez global es (se eliminan las filas y columnas correspondientes a los
nudos prescritos):
a
b
a
b
31764,706 -31764,706
0
0
a
-31764,706 63529,412 -31764,706
0
b
KG =
0
-31764,706 58235,294 -26470,588 a
0
0
-26470,588 26470,588 b
Qa
Qb
Qc
Qd
63529,412 -31764,706
Qb
300
*
=
-31764,706 58235,294
Qc
0
Los desplazamientos en los nudos son los siguientes:
Qa
Qb
Qc
Qd
=
0
= 6,49E-03
cm
= 3,54E-03
=
0
2= 11764,70588
2=
2= 5882,352941
2=
-1
-11764,70588 11764,70588
-1
-11764,70588 11764,70588
-1
-5882,352941 5882,352941
Qa
Qb
0
= 76,389
6,49E-03
Qb
Qc
6,49E-03
= -34,722
3,54E-03
Qc
Qd
3,54E-03
= -20,833
0
R1= -206,25 KN
R4= -93,75 KN
b)
En caso que las condiciones de apoyo de la izquierda sean las mismas que las de la derecha,
los resultados no sufriran cambios ya que el nico grado de libertad estara restringido de
la misma manera. Como no hay otros grados de libertad no interesara el tipo de apoyo con
tal que ese nico grado de libertad este restringido.
Datos:
E = 30x106 [psi]
A1 = A3 = A4 = A6 = A7 = 20 [in2]
A2 = A5 = 15 [in2]
Solucin.a)
Tabla de incidencia nodal
Nudo coord.(x) coord.(y)
1
0
0
2
240
0
3
480
0
4
360
144
5
120
144
Tabla de conectividad
Elemento nudo i nudo j
1
1
5
2
1
2
3
5
2
4
2
4
5
2
3
6
4
3
7
5
4
le (in) Ae (in)
187,44
20
240
15
187,44
20
187,44
20
240
15
187,44
20
240
20
2
0,4928
0,5929
-0,4928
-0,5929
1
1
0
k2= 1,875x106 -1
0
4
0
0
0
0
2
0
0
0
0
3
-1
0
1
0
9
-0,4096
-0,4928
0,4096
0,4928
10
-0,4928
-0,5929
0,4928
0,5929
1
2
9
10
1
2
3
4
9
0,4096
-0,4928
k3= 3,201x106 -0,4096
0,4928
10
-0,4928
0,5929
0,4928
-0,5929
3
-0,4096
0,4928
0,4096
-0,4928
4
0,4928
-0,5929
-0,4928
0,5929
9
10
3
4
3
0,4096
0,4928
6
k4= 3,201x10 -0,4096
-0,4928
4
0,4928
0,5929
-0,4928
-0,5929
7
-0,4096
-0,4928
0,4096
0,4928
8
-0,4928
-0,5929
0,4928
0,5929
3
4
7
8
l
0,64
1
-0,64
0,64
1
-0,64
1
m
0,77
0
0,77
0,77
0
0,77
0
3
1
0
k5= 1,875x106 -1
0
4
0
0
0
0
5
-1
0
1
0
6
0
0
0
0
3
4
5
6
7
0,4096
-0,4928
6
k6= 3,201x10 -0,4096
0,4928
8
-0,4928
0,5929
0,4928
-0,5929
9
1
0
6
k7= 2,5x10 -1
0
8
0
0
0
0
10
0
0
0
0
7
-1
0
1
0
5
-0,4096
0,4928
0,4096
-0,4928
6
0,4928
-0,5929
-0,4928
0,5929
7
8
5
6
9
10
7
8
Matriz de rigidez global (Se eliminan las filas y columnas correspondientes a los grados de
libertad 4, 5 y 6 ya que se encuentran prescritos)
1
3,1861
1,5775
-1,875
1,5775
1,8979
-1,875
6,3723
-1,875
K= 10
10
-1,3111 -1,5775 1
-1,5775 -1,8979 2
3,1861
1,5775
-1,875
106
0
0
-1,3111
-1,5775
1,5775
1,8979
0
0
0
-1,5775
-1,8979
3,7958
Q1
0
-1,875
0
0 -1,3111 -1,5775
Q2
-200000
0
0
0 -1,5775 -1,8979
Q
0
6,3723 -1,3111 -1,5775 -1,3111 1,5775
3
0
-1,3111 5,1223
0
-2,5
0 x Q7 =
Q8
0
-1,5775
0 3,7958
0
0
Q9
0
-1,3111
-2,5
0 5,1223
0
Q
0
1,5775
0
0
0 3,7958
10
Luego
Q1 =
Q2 =
Q3 =
Q7 =
Q8 =
Q9 =
Q10 =
0,1773
-0,79
0,0887
-0,0825 in
0,0368
-0,2154
-0,3581
0,1773
-0,79
6
2 = 30x10 x -1 0 1 0 x 0,0887 = -11082,25108 psi
240
0
-0,2154
-0,3581
6
3 = 30x10 x 0,64 -0,77 -0,64 0,77 x 0,0887 = 12987,01299 psi
187,44
0
0,0887
0
6
4 = 30x10 x -0,64 -0,77 0,64 0,77 x -0,0825 = -12987,01299 psi
187,44
0,0368
0,0887
0
6
0
= -11082,25108 psi
5 = 30x10 x -1 0 1 0 x
240
0
-0,0825
0,0368
6
0
= -12987,01299 psi
6 = 30x10 x 0,64 -0,77 -0,64 0,77 x
187,44
0
-0,2154
-0,3581
6
7 = 30x10 x -1 0 1 0 x -0,0825 = 16623,37662 psi
240
0,0368
Las reacciones de los apoyos sern las siguientes
0,1773
-0,79
0,0887
0
R4
R5 =10
R6
0 0
6
3,7958
0 0 -1,875
0 0
-1,5775 1,8979
0
-0,0825
0,0368
-0,2154
-0,3581
R4
400000
0
R5 =
lb
R6
-200000
b)
En el caso de que la carga sea horizontal, la tensin en el elemento inclinado ser cero, tal
como se demuestra a continuacin:
1,5775
1,8979
0
0
0
-1,5775
-1,8979
Q1
-200000
-1,875
0
0 -1,3111 -1,5775
Q2
0
0
0
0 -1,5775 -1,8979
Q3
0
6,3723 -1,3111 -1,5775 -1,3111 1,5775
0
-1,3111 5,1223
0
-2,5
0 x Q7 =
Q8
0
-1,5775
0 3,7958
0
0
Q9
0
-1,3111
-2,5
0 5,1223
0
Q10
0
1,5775
0
0
0 3,7958
Luego:
Q1 =
Q2 =
Q3 =
Q7 =
Q8 =
Q9 =
Q10 =
-0,2133
0,1773
-0,1067
-0,0533 in
-0,0443
-0,0533
0,0443
A su vez:
(
)
(
[
[
SET 4
Datos.
E = 30x106 [psi]
I = 345 [in4]
P = 5000 [lb]
Q = 1500 [lb/ft]
Solucin
a)
Las matrices de rigidez de cada elemento sern:
1
2
12 360 -12 360
575000 1,7E+07
17250000 6,9E+08
K1= 47916,7 360 14400 -360 7200
=
-12 -360 12 -360
-575000 -2E+07
360 7200 -360 14400
17250000 3,5E+08
3
-575000
-2E+07
575000
-2E+07
4
1,7E+07
3,5E+08
-2E+07
6,9E+08
1
2
3
4
3
4
12 360 -12 360
575000 1,7E+07
17250000 6,9E+08
K2= 47916,7 360 14400 -360 7200
=
-12 -360 12 -360
-575000 -2E+07
360 7200 -360 14400
17250000 3,5E+08
5
-575000
-2E+07
575000
-2E+07
6
1,7E+07
3,5E+08
-2E+07
6,9E+08
3
4
5
6
5
12 720 -12 720
71875
4312500
K3= 5989,58 720 57600 -720 28800
=
-12 -720 12 -720
-71875
720 28800 -720 57600
4312500
6
4312500
3,5E+08
-4E+06
1,7E+08
7
-71875
-4E+06
71875
-4E+06
8
4312500
1,7E+08
-4E+06
3,5E+08
5
6
7
8
7
8
12 504 -12 504
209548,1 8801020
8801020,4 4,9E+08
K4= 17462,3 504 28224 -504 14112
=
-12 -504 12 -504
-209548,1 -9E+06
504 14112 -504 28224
8801020,4 2,5E+08
9
-209548
-9E+06
209548
-9E+06
10
8801020
2,5E+08
-9E+06
4,9E+08
7
8
9
10
10
575000
17250000
-575000
1,7E+07
1,7E+07
6,9E+08
-17250000 3,5E+08
Ta
-575000
-1,7E+07
-575000
1,7E+07
Tb
-1E+07
-71875
4312500
1E+09
-4312500
1,7E+08
Tc
-4E+06
281423,1
1,7E+07 3,45E+08
K=
1150000
0
-575000
17250000
-2E+07
646875
-9E+06
2,5E+08
Td
-209548,1
-9E+06
209548
-9E+06
8801020,4 2,5E+08
-9E+06
4,9E+08 10 Te
3,5E+08 -1,3E+07
-71875
2
3
4
6
8
10
6,9E+08 -1,7E+07 345000000
0
0
0
2
Ta
0
-2E+07 1150000
0
1,7E+07
0
0
3
b
-5000
3,5E+08
0
1,38E+09 3,5E+08
0
0
4 * Tb =
0
0
17250000 345000000 1E+09 1,73E+08
0
6
Tc
0
0
0
0
1,7E+08 8,38E+08 2,5E+08 8
Td
-73500
0
0
0
0
2,46E+08 4,9E+08 10
Te
73500
De donde los desplazamientos y rotaciones son los siguientes:
Ta
-0,0003411
b
-0,013177
Tb = 2,3413E-05
Tc
0,00024748
Td
-0,000214
Te
0,00025614
1
2
3
4
5
6
7
8
9
10
575000 17250000 -575000 1,7E+07
0
0
0
0
0
0
1
a
0
0
-575000 -2E+07 646875 -1E+07 -71875 4312500
0
0
5 * Ta
0
0
0
0
-71875 -4E+06 281423,1 4488520 -209548 8801020 7
b
0
0
0
0
0
0
-209548,1 -9E+06 209548 -9E+06 9
Tb
c
Tc
d
Td
e
Te
Ra
2096,131
0
Rc = 3048,2143 - 0
Rd
226,40306
-5250
Re
-370,7483
-5250
Ra
Rc
Rd
Re
= 2096,131 lb
= 3048,2143 lb
= 5476,4031 lb
= 4879,2517 lb
a-c = -3,41E-04 in
c-d = 3,71E-03 +
3,21E-03
c-d = 6,92E-03 in
d-e = -2,25E-03 + -2,69E-03
d-e = -4,94E-03 in
b)
Obtencin de las matrices de rigidez para cada elemento:
1
2
12 360 -12 360
575000 1,7E+07
K 1= 47916,7 360 14400 -360 7200 = 17250000 6,9E+08
-12 -360 12 -360
-575000 -2E+07
360 7200 -360 14400
17250000 3,5E+08
3
-575000
-2E+07
575000
-2E+07
4
1,7E+07
3,5E+08
-2E+07
6,9E+08
1
2
3
4
3
4
12 360 -12 360
575000 1,7E+07
K 2= 47916,7 360 14400 -360 7200 = 17250000 6,9E+08
-12 -360 12 -360
-575000 -2E+07
360 7200 -360 14400
17250000 3,5E+08
5
-575000
-2E+07
575000
-2E+07
6
1,7E+07
3,5E+08
-2E+07
6,9E+08
3
4
5
6
5
12 720 -12 720
71875
K 3= 5989,58 720 57600 -720 28800 = 4312500
-12 -720 12 -720
-71875
720 28800 -720 57600
4312500
6
4312500
3,5E+08
-4E+06
1,7E+08
7
-71875
-4E+06
71875
-4E+06
7
8
12 504 -12 504
209548,1 8801020
K 4= 17462,3 504 28224 -504 14112 = 8801020,4 4,9E+08
-12 -504 12 -504
-209548,1 -9E+06
504 14112 -504 28224
8801020,4 2,5E+08
8
4312500
1,7E+08
-4E+06
3,5E+08
9
-209548
-9E+06
209548
-9E+06
5
6
7
8
10
8801020
2,5E+08
-9E+06
4,9E+08
7
8
9
10
La matriz de rigidez global de la que se eliminan las columnas y filas correspondientes a los
grados de libertad prescritos, ser:
1
2
3
4
5
6
7
8
575000 17250000 -575000 1,7E+07
0
0
0
0
1,7E+07 6,9E+08 -17250000 3,5E+08
0
0
0
0
-575000 -1,7E+07 1150000
0
-575000 1,7E+07
0
0
0
1,4E+09 -1,7E+07 3,5E+08
0
0
K = 1,7E+07 3,45E+08
0
0
-575000 -2E+07 646875 -1E+07
-71875 4312500
0
0
17250000 3,5E+08 -1,3E+07 1E+09 -4312500 1,7E+08
0
0
0
0
-71875 -4E+06 281423,1 4488520
0
0
0
0
4312500 1,7E+08 4488520,4 8,4E+08
0
0
0
0
0
0
-209548,1 -9E+06
0
0
0
0
0
0
8801020,4 2,5E+08
9
10
0
0
1
0
0
2
0
0
3
0
0
4
0
0
5
0
0
6
-209548 8801020 7
-9E+06 2,5E+08 8
209548 -9E+06 9
-9E+06 4,9E+08 10
2
3
3
6
2 6,9E+08 -1,7E+07 345000000
0
3 -2E+07 1150000
0
1,7E+07
4 3,5E+08
0
1,38E+09 3,5E+08
6
0
17250000 345000000 1E+09
7
0
0
0
-4E+06
8
0
0
0
1,7E+08
10
0
0
0
0
7
0
0
0
-4312500
281423,1
4488520
8801020
8
10
0
0
Ta
0
0
0
b
-5000
0
0
Tb
0
1,7E+08
0
* Tc =
0
4488520 8801020,4
d
-5250
8,4E+08 246428571
Td
-73500
2,5E+08 492857143
Te
73500
Datos.
I1=I3=350 in4
I2=175 in4
A1=A3= 18 in2
A2=9in2
Q=1500 lb/ft
E= CTTE
Solucin
Conectividad
Nudo
A
B
C
D
Elem.
1
2
3
coord.(x) coord.(y)
0
0
ft
10
20
ft
30
20
ft
40
0
ft
nudo i
A
B
D
nudo j
B
C
C
le (in)
268,32
240
268,32
Ae (in)
18
9
18
I(in4)
350
175
350
l
0,447
1
-0,447
k'e=
EA/L
0
0
(-EA/L)
0
0
0
12EI/L3
6EI/L2
0
(-12EI/L3) (6EI/L2)
0
6EI/L2
4EI/L
0
(-6EI/L2)
2EI/L
(-EA/L)
0
0
EA/L
0
0
0
(-12EI/L3) (-6EI/L2)
0
12EI/L3 (-6EI/L2)
0
6EI/L2
2EI/L
0
(-6EI/L2)
4EI/L
Elemento 1
0,06708
0
k1= E *
0
-0,06708
0
0
0
0,000217
0,029168
0
-0,00022
0,029168
0
-0,06708
0
0
0,02917
0
-0,0002174 0,02916844
5,21765
0
-0,0291684 2,60882528
0
0,06708
0
0
-0,02917
0
0,00021742 -0,0291684
2,60883
0
-0,0291684 5,21765057
m
0,894
0
0,894
Elemento 2
0,0375
0
k2= E *
0
-0,0375
0
0
0
0,000152
0,018229
0
-0,00015
0,018229
0
-0,0375
0
0
0,01823
0
-0,0001519 0,01822917
2,91667
0
-0,0182292 1,45833333
0
0,0375
0
0
-0,01823
0
0,00015191 -0,0182292
1,45833
0
-0,0182292 2,91666667
Elemento 3
0,06708
0
3
E
*
k =
0
-0,06708
0
0
0
0,000217
0,029168
0
-0,00022
0,029168
0
-0,06708
0
0
0,02917
0
-0,0002174 0,02916844
5,21765
0
-0,0291684 2,60882528
0
0,06708
0
0
-0,02917
0
0,00021742 -0,0291684
2,60883
0
-0,0291684 5,21765057
Matriz de Transformacin L
l
(-m)
L= 0
0
0
0
m
l
0
0
0
0
0 0 0
0 0 0
1 0 0
0 l m
0 (-m) l
0 0 0
0
0
0
0
0
1
Elemento 1
0.447 0.894 0
0
0 0
-0.894 0.447 0
0
0 0
L=
0
0
1
0
0 0
0
0
0 0.447 0.894 0
0
0
0 -0.894 0.447 0
0
0
0
0
0 1
Elemento 2
1
0
L= 0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
Elemento 3
-0.447 0.894 0
0
0 0
-0.894 -0.447 0
0
0 0
L=
0
0
1
0
0 0
0
0
0 -0.447 0.894 0
0
0
0 -0.894 -0.447 0
0
0
0
0
0 1
Transpuesta de L
Elemento 1
0.447 -0.894 0
0
0 0
0.894 0.447 0
0
0 0
LT=
0
0 1
0
0 0
0
0 0 0.447 -0.894 0
0
0 0 0.894 0.447 0
0
0 0
0
0 1
Elemento 2
1
0
LT= 0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
Elemento 3
-0.447 -0.894 0
0
0 0
0.894 -0.447 0
0
0 0
LT=
0
0 1
0
0 0
0
0 0 -0.447 -0.894 0
0
0 0 0.894 -0.447 0
0
0 0
0
0 1
Elemento 1
0.013577769 0.026721122 -0.026076585 -0.013577769 -0.026721122 -0.026076585
0.026721122 0.053659452 0.013038293 -0.026721122 -0.053659452 0.013038293
k1=E -0.026076585 0.013038293 5.217650566 0.026076585 -0.013038293 2.608825283
-0.013577769 -0.026721122 0.026076585 0.013577769 0.026721122 0.026076585
-0.026721122 -0.053659452 -0.013038293 0.026721122 0.053659452 -0.013038293
-0.026076585 0.013038293 2.608825283 0.026076585 -0.013038293 5.217650566
Elemento 2
0.0375
0
0 -0.0375
0
0
0 0.00015191 0.018229167
0 -0.00015191 0.018229167
k2=E
0 0.018229167 2.916666667
0 -0.018229167 1.458333333
-0.0375
0
0 0.0375
0
0
0 -0.00015191 -0.018229167
0
0.00015191 -0.018229167
0 0.018229167 1.458333333
0 -0.018229167 2.916666667
Elemento 3
0.013577769
-0.026721122
k3=E* -0.026076585
-0.013577769
0.026721122
-0.026076585
0.026721122 0.026076585
-0.0375
0
0
0.053811362 0.005190874
0 -0.00015191 0.018229167
0.005190874 8.134317233
0 -0.018229167 1.458333333
0
0 0.051077769 -0.026721122 0.026076585
-0.00015191 -0.018229167 -0.026721122 0.053811362 -0.005190874
0.018229167 1.458333333 0.026076585 -0.005190874 8.134317233
VECTOR DE CARGA
15000 lb
15000 lb
600000 lb*in
600000 lb*in
0
-15000
F= -600000
0
-15000
600000
Lb
Lb
Lb*in
Lb
Lb
Lb*in
Resolviendo obtenemos:
138671.9216
-370741.4072
1
-91140.11273
Q=
*
E
-138671.9216
-370741.4072
91140.11273
in
in
rad
in
in
rad
B
B
B
BC
RAX'
RAY'
RBX'
RBY'
RA
RB
Componente Y:
Por simetra:
Conectividad
Nudo
A
B
C
D
Elem.
1
2
3
coord.(x) coord.(y)
0
0
ft
10
20
ft
30
20
ft
40
0
ft
nudo i
A
B
D
nudo j
B
C
C
le (in)
268,32
240
268,32
Ae (in)
18
9
18
I(in4)
350
175
350
l
0,447
1
-0,447
k'e=
EA/L
0
0
(-EA/L)
0
0
0
12EI/L3
6EI/L2
0
(-12EI/L3) (6EI/L2)
0
6EI/L2
4EI/L
0
(-6EI/L2)
2EI/L
(-EA/L)
0
0
EA/L
0
0
0
(-12EI/L3) (-6EI/L2)
0
12EI/L3 (-6EI/L2)
0
6EI/L2
2EI/L
0
(-6EI/L2)
4EI/L
m
0,894
0
0,894
Elemento 1
0,06708
0
k1= E *
0
-0,06708
0
0
0
0,000217
0,029168
0
-0,00022
0,029168
0
-0,06708
0
0
0,02917
0
-0,0002174 0,02916844
5,21765
0
-0,0291684 2,60882528
0
0,06708
0
0
-0,02917
0
0,00021742 -0,0291684
2,60883
0
-0,0291684 5,21765057
Elemento 2
0,0375
0
k2= E *
0
-0,0375
0
0
0
0,000152
0,018229
0
-0,00015
0,018229
0
-0,0375
0
0
0,01823
0
-0,0001519 0,01822917
2,91667
0
-0,0182292 1,45833333
0
0,0375
0
0
-0,01823
0
0,00015191 -0,0182292
1,45833
0
-0,0182292 2,91666667
Elemento 3
0,06708
0
k3= E *
0
-0,06708
0
0
0
0,000217
0,029168
0
-0,00022
0,029168
0
-0,06708
0
0
0,02917
0
-0,0002174 0,02916844
5,21765
0
-0,0291684 2,60882528
0
0,06708
0
0
-0,02917
0
0,00021742 -0,0291684
2,60883
0
-0,0291684 5,21765057
Matriz de Transformacin L
l
(-m)
L= 0
0
0
0
m
l
0
0
0
0
0 0 0
0 0 0
1 0 0
0 l m
0 (-m) l
0 0 0
0
0
0
0
0
1
Elemento 1
0.447 0.894 0
0
0 0
-0.894 0.447 0
0
0 0
L=
0
0
1
0
0 0
0
0
0 0.447 0.894 0
0
0
0 -0.894 0.447 0
0
0
0
0
0 1
Elemento 2
1
0
L= 0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
Elemento 3
-0.447 0.894 0
0
0 0
-0.894 -0.447 0
0
0 0
L=
0
0
1
0
0 0
0
0
0 -0.447 0.894 0
0
0
0 -0.894 -0.447 0
0
0
0
0
0 1
Transpuesta de L
Elemento 1
0.447 -0.894 0
0
0 0
0.894 0.447 0
0
0 0
LT=
0
0 1
0
0 0
0
0 0 0.447 -0.894 0
0
0 0 0.894 0.447 0
0
0 0
0
0 1
Elemento 2
1
0
LT= 0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
Elemento 3
-0.447 -0.894 0
0
0 0
0.894 -0.447 0
0
0 0
LT=
0
0 1
0
0 0
0
0 0 -0.447 -0.894 0
0
0 0 0.894 -0.447 0
0
0 0
0
0 1
ke= LT * ke * L
Elemento 1
0.013577769 0.026721122 -0.026076585 -0.013577769 -0.026721122 -0.026076585
0.026721122 0.053659452 0.013038293 -0.026721122 -0.053659452 0.013038293
k1=E -0.026076585 0.013038293 5.217650566 0.026076585 -0.013038293 2.608825283
-0.013577769 -0.026721122 0.026076585 0.013577769 0.026721122 0.026076585
-0.026721122 -0.053659452 -0.013038293 0.026721122 0.053659452 -0.013038293
-0.026076585 0.013038293 2.608825283 0.026076585 -0.013038293 5.217650566
Elemento 2
0.0375
0
0 -0.0375
0
0
0 0.00015191 0.018229167
0 -0.00015191 0.018229167
k2=E
0 0.018229167 2.916666667
0 -0.018229167 1.458333333
-0.0375
0
0 0.0375
0
0
0 -0.00015191 -0.018229167
0
0.00015191 -0.018229167
0 0.018229167 1.458333333
0 -0.018229167 2.916666667
Elemento 3
0.013577769
-0.026721122
k3=E* -0.026076585
-0.013577769
0.026721122
-0.026076585
0.026721122 0.026076585
-0.0375
0
0
0.053811362 0.005190874
0 -0.00015191 0.018229167
0.005190874 8.134317233
0 -0.018229167 1.458333333
0
0 0.051077769 -0.026721122 0.026076585
-0.00015191 -0.018229167 -0.026721122 0.053811362 -0.005190874
0.018229167 1.458333333 0.026076585 -0.005190874 8.134317233
VECTOR DE CARGA
Debido a que existen rtulas en los nudos B y C los momentos en dichos nudos son igual a
0:
15000 lb
15000 lb
0 lb*in
0 lb*in
0
-15000
0
F=
0
-15000
0
Lb
Lb
Lb*in
Lb
Lb
Lb*in
Resolviendo obtenemos:
99683,11685
-329431,4291
1
-1032,749453
Q=
*
E
-99683,11685
-329431,4291
1032,749453
in
in
rad
in
in
rad
C
B
B
BC
Identificar cada una de las componentes de la ecuacin 5.11 indicando a cual de cada una
de las componentes de la 5.10 corresponde.
Ecuacin 5.10
*
)|
)|
(
|
Donde:
-
) |
) |
) |
) |
SET 5
a)
Mediante la teora elemental de las vigas:
Calculo de tension:
( )
( )
Calculo de deformacin:
Si:
x=90 ; y=0
(
Si:
x=90 ;
(
Si:
x=0 ; y=ymax
b)
2
e=1
e=1
1
0,33
0
78554595444
0,33
1
0
0
0
0,335
78554595444 25923016496
25923016496 78554595444
0
0
0
26315789474
Elemento
1
2
1
3
k
2
4
4
2
Haciendo el producto
C*B1=
8,72829E+11 -5,76067E+11
0 5,76067E+11 -8,72829E+11
2,88034E+11 -1,74566E+12
0 1,74566E+12 -2,88034E+11
-5,848E+11
2,92398E+11 -5,84795E+11
0 -2,924E+11
C*B2=
-8,73E+11
5,76E+11
-5,76E+11
8,73E+11
0,00E+00
-2,88E+11
1,75E+12
-1,75E+12
2,88E+11
0,00E+00
5,85E+11
-2,92E+11
5,85E+11
0,00E+00
0,00E+00
2,92E+11
459544383,3 -261194029,9
263157894,7
129615082,5 -196386488,6
131578947
-261194030
K1=
263157894,7 -131578947,4
263157894,7
129615082,5 -785545954,4
-196386489
129615082,5
-261194030
K2=
131578947
785545954,4 -129615082,5
131578947,4
263157894,7
263157894,7
129615082,5 -785545954,4
129615082,5
131578947,4 -65789473,68
65789473,7
129615082,5 -196386488,6
131578947
129615082,5 -65789473,7
131578947
785545954,4 -129615082,5
0 -129615082,5
131578947,4
4,60E+08 -2,61E+08
-2,61E+08
196386488,6
263157894,7 -131578947,4
-196386489
0 -129615082,5
131578947,4 -65789473,68
459544383,3 -261194029,9
129615082,5 -65789473,7
196386488,6
65789473,7
1,30E+08
Q1
2,63E+08
Q2
O
=
2,63E+08 -1,32E+08
4,60E+08
0,00E+00
Q3
1,30E+08 -7,86E+08
0,00E+00
8,51E+08
Q4
-7000
Q1
-1,57E-05
Q2
-8,60E-05
Q3
-1,57E-05
Q4
-8,53E-05
Q5
0,00E+00
Q6
0,00E+00
luego:
N/m2
Para el elemento 2 :
N/m2
c) Los resultados obtenidos mejoraran con una mayor discretizacion ya que los resultados
de los esfuerzos y deformaciones serian obtenidos en ms puntos de esta manera se tendra
mayor precisin pero siempre es conveniente iniciar con pocos elementos e ir aumentando
el nmero segn la precisin deseada en especial en las zonas donde se tiene saltos bruscos
de tensiones y se necesitan mayor precisin de los resultados.
d) Las causas para la diferencia de resultados entre los dos mtodos seria que las tensiones
en un elemento cuadrilateral no son constantes dentro del elemento ya que son funciones de
x y h por lo que varan.
En este problema existirn cuatro variaciones de tensiones ya que se tendrn las
coordenadas naturales para cada nudo estos valores debern ser promediados para obtener
el valor general.