Documentos de Académico
Documentos de Profesional
Documentos de Cultura
Metodos Numericos Unidad Nº2
Metodos Numericos Unidad Nº2
Problema: Disear la losa aligerada que se muestra utilizando aproximaciones sucesivas:
Universidad Nacional de Cajamarca Facultad de Ingeniara
Curso: Mtodos Numricos 2013 Escuela de Ingeniera Civil
Ing Wilson E. Vargas Vargas 2013 Pgina 30
Problema 5:
Nudo 1 Fx=0 - F1*cos30+F2*cos60+ F1h = 0
Fy=0 - F1*sen30+F3*sen60+F1v = 0
Nudo 2 Fx=0 F2 +F1*cos30+ F2h + H2 = 0
Fy=0 F1*sen30 + F2v + R2 = 0
Nudo 3 Fx=0- F2 - F3*cos60+ F3h = 0
Fy=0 F3*sen30 + F3v + R3 = 0
Usando Matrices- (completando coeficientes se tiene):
- F1*cos30 0 F3*sen60 0 0 0 F1 0
F1*sen30 -F3*sen60 0 0 0 F2 -1000
Universidad Nacional de Cajamarca Facultad de Ingeniara
Curso: Mtodos Numricos 2013 Escuela de Ingeniera Civil
Ing Wilson E. Vargas Vargas 2013 Pgina 31
Problema N6
Los sistemas idealizados masa-resorte tienen diversas aplicaciones en la ingeniera. La figura muestra cuatro
resortes en serie cuando se comprimen con una fuerza de 2500 Kg.
En el equilibrio, se puede desarrollar las ecuaciones de equilibrio de fuerza definiendo las interrelaciones
entre los resortes:
K2(X2-X1) = K1X1
K3(X3-X2) = K2(X2-X1)
K4(X4-X3) = K3(X2-X2)
F= K4(X4-X3)
Dnde K son la constantes del resorte, si K1=150, K2=50, K3=75, K4=225 Kg/s
2
, calcular X
Acomodando las ecuaciones tenemos:
(K1+ K2)* X1 - K2X2 = 0 (I)
K2* X1- (K2+K3)X2+K3X3 = 0 (II)
K3* X2- (K3+K4)X3+K4X4 = 0 (III)
-K4* X3 +K4X4 = F (IV)
Sustituyendo datos y escribiendo matricialmente tenemos:
SISTEMAS MASA-RESORTE:
Ejemplo de aplicacin
Si se tiene el siguiente sistema
K
m1
m2
m1
m2
Universidad Nacional de Cajamarca Facultad de Ingeniara
Curso: Mtodos Numricos 2013 Escuela de Ingeniera Civil
Ing Wilson E. Vargas Vargas 2013 Pgina 32
Se solicita hallar los desplazamientos de las masas, si se sabe que m1=50kg, m2=35kg, m3=44kg, K=500N/m
Solucin:
Haciendo un Diagrama de cuerpo libre del sistema dado:
DCL-Para la masa 1
Las fuerzas que actan en la masa 1 son:
K(x
1
) =500(x
1
)
K(x2-x1) =500(x2-x1) K(x2-x1) =500(x2-x1)
W=m1g=500N
Se abe adems que F= ma , entonces se tiene
ma=2k(x2-x1)-m1g-kx1
ma=1000(x2-x1)-500-500*x m1g-kx1
Para m2:
K(x2-x1)=500(x2-x1) K(x2-x1)=500(x2-x1)
K(x3-x2)=500(x3-x2)
W2=m2g
ma=k(x3-x2)-m2g-2k(x2-x1)ma=500(x3-x2)-350-1000(x2-x1)
m3
m3
m1
Universidad Nacional de Cajamarca Facultad de Ingeniara
Curso: Mtodos Numricos 2013 Escuela de Ingeniera Civil
Ing Wilson E. Vargas Vargas 2013 Pgina 33
Para m3:
k(x
3
-x
2
)
w=m
3
g
ma= m
3
g- k(x
3
-x
2
) ma=440-500(x
3
-x
2
)
Si a=
=0, se tiene:
ORDENANDO LAS ECUACIONES Y VARIABLES se tiene:
-150(x1) + 100(X2) =50
100(x1) -150(x2) + 50(x3) =35
-50(x2) + 50(x3) = 44
El sistema de ecuaciones ser resuelto por dos mtodos
Mtodo de Gauss Seidel:
Ordenando los coeficientes de las ecuaciones para poder aplicar el mtodo, se tiene:
-150(x1) +100(x2) +0x3 =50 (1)
100(x1) -150(x2) +50(x3) =35 (2)
0x1 -50(x2) +50(x3) =44 (3)
Despejando:
m3
Universidad Nacional de Cajamarca Facultad de Ingeniara
Curso: Mtodos Numricos 2013 Escuela de Ingeniera Civil
Ing Wilson E. Vargas Vargas 2013 Pgina 34
De(1)
De(2)
De(3)
PROCESO ITERATIVO:
I
0 0 0 0
1 -0,333333333 -0,455555556 0,424444444
2 -0,637037037 -0,51654321 0,36345679
3 -0,677695473 -0,563978052 0,316021948
4 -0,709318701 -0,600871818 0,279128182
5 -0,733914546 -0,62956697 0,25043303
6 -0,753044647 -0,651885421 0,228114579
7 -0,767923614 -0,669244216 0,210755784
8 -0,779496144 -0,682745502 0,197254498
9 -0,788497001 -0,693246501 0,186753499
10 -0,795497667 -0,701413945 0,178586055
11 -0,80094263 -0,707766402 0,172233598
12 -0,805177601 -0,712707202 0,167292798
13 -0,808471468 -0,716550046 0,163449954
14 -0,811033364 -0,719538924 0,160461076
15 -0,81302595 -0,721863608 0,158136392
16 -0,814575739 -0,723671695 0,156328305
17 -0,81578113 -0,725077985 0,154922015
18 -0,816718657 -0,726171766 0,153828234
19 -0,817447844 -0,727022485 0,152977515
Universidad Nacional de Cajamarca Facultad de Ingeniara
Curso: Mtodos Numricos 2013 Escuela de Ingeniera Civil
Ing Wilson E. Vargas Vargas 2013 Pgina 35
20 -0,81801499 -0,727684155 0,152315845
21 -0,818456103 -0,728198787 0,151801213
22 -0,818799191 -0,728599057 0,151400943
23 -0,819066038 -0,728910377 0,151089623
24 -0,819273585 -0,729152516 0,150847484
25 -0,81943501 -0,729340846 0,150659154
26 -0,819560564 -0,729487324 0,150512676
27 -0,819658216 -0,729601252 0,150398748
28 -0,819734168 -0,729689863 0,150310137
29 -0,819793242 -0,729758782 0,150241218
30 -0,819839188 -0,729812386 0,150187614
31 -0,819874924 -0,729854078 0,150145922
32 -0,819902719 -0,729886505 0,150113495
33 -0,819924337 -0,729911726 0,150088274
34 -0,819941151 -0,729931343 0,150068657
35 -0,819954228 -0,7299466 0,1500534
36 -0,8199644 -0,729958467 0,150041533
37 -0,819972311 -0,729967696 0,150032304
38 -0,819978464 -0,729974875 0,150025125
39 -0,81998325 -0,729980458 0,150019542
40 -0,819986972 -0,729984801 0,150015199
41 -0,819989867 -0,729988178 0,150011822
42 -0,819992119 -0,729990805 0,150009195
43 -0,81999387 -0,729992849 0,150007151
44 -0,819995232 -0,729994438 0,150005562
45 -0,819996292 -0,729995674 0,150004326
46 -0,819997116 -0,729996635 0,150003365
47 -0,819997757 -0,729997383 0,150002617
48 -0,819998255 -0,729997965 0,150002035
Universidad Nacional de Cajamarca Facultad de Ingeniara
Curso: Mtodos Numricos 2013 Escuela de Ingeniera Civil
Ing Wilson E. Vargas Vargas 2013 Pgina 36
49 -0,819998643 -0,729998417 0,150001583
50 -0,819998945 -0,729998769 0,150001231
51 -0,819999179 -0,729999042 0,150000958
52 -0,819999362 -0,729999255 0,150000745
53 -0,819999503 -0,729999421 0,150000579
54 -0,819999614 -0,729999549 0,150000451
55 -0,8199997 -0,72999965 0,15000035
56 -0,819999766 -0,729999727 0,150000273
57 -0,819999818 -0,729999788 0,150000212
58 -0,819999859 -0,729999835 0,150000165
59 -0,81999989 -0,729999872 0,150000128
60 -0,819999914 -0,7299999 0,1500001
61 -0,819999933 -0,729999922 0,150000078
62 -0,819999948 -0,72999994 0,15000006
63 -0,81999996 -0,729999953 0,150000047
64 -0,819999969 -0,729999963 0,150000037
65 -0,819999976 -0,729999972 0,150000028
66 -0,819999981 -0,729999978 0,150000022
67 -0,819999985 -0,729999983 0,150000017
68 -0,819999989 -0,729999987 0,150000013
69 -0,819999991 -0,72999999 0,15000001
70 -0,819999993 -0,729999992 0,150000008
71 -0,819999995 -0,729999994 0,150000006
72 -0,819999996 -0,729999995 0,150000005
73 -0,819999997 -0,729999996 0,150000004
74 -0,819999997 -0,729999997 0,150000003
75 -0,819999998 -0,729999998 0,150000002
76 -0,819999998 -0,729999998 0,150000002
77 -0,819999999 -0,729999999 0,150000001
Universidad Nacional de Cajamarca Facultad de Ingeniara
Curso: Mtodos Numricos 2013 Escuela de Ingeniera Civil
Ing Wilson E. Vargas Vargas 2013 Pgina 37
78 -0,819999999 -0,729999999 0,150000001
79 -0,819999999 -0,729999999 0,150000001
80 -0,819999999 -0,729999999 0,150000001
81 -0,82 -0,729999999 0,150000001
82 -0,82 -0,73 0,15
83 -0,82 -0,73 0,15
De esta manera obtenemos los valores para:
X1= -0,82, X2 = -0,73, X3 = 0,15, que vienen a ser las respectivas longitudes de las deformaciones que
sufren los resortes
MTODO: USANDO LA HOJA DE CLCULO DE EXCEL:
Las ecuaciones y variables obtenidas del sistema son:
100(x2) -150(x1) =50
50(x3) -150(x2)+100(x1) =35
50(x3) -50(x2) =4
El sistema expresado matricialmente es:
Hallando la inversa de la matriz de los coeficientes:
-0,02 -0,02 0,02 50
-0,02 -0,03 0,03 35
-0,02 -0,03 0,05 44
Hallando los valores de X1,X2,X3:
x1 -0,82
Universidad Nacional de Cajamarca Facultad de Ingeniara
Curso: Mtodos Numricos 2013 Escuela de Ingeniera Civil
Ing Wilson E. Vargas Vargas 2013 Pgina 38
x2 -0,73
x3 0,15
Que vienen a ser los valores las longitudes de las deformaciones que sufren los resortes
Conclusiones
-Comprender y resolver problemas que incluyan sistemas de ecuaciones no lineales.
-Graficar el sistema de ecuaciones no lineales.
-Representar en forma matricial sistemas de ecuaciones.
-Comprender, analizar, construir sistemas de ecuaciones, como solucin de problemas de ingeniera.
-Aplicar la herramientas de Excel para la solucin de sistemas de ecuaciones lineales y no lineales