Ejemplo 4.

AE 10000

Matriz AF Elementos L AE/L Qp P

0 1 1 6 1666.667 10 20
0.8 0.6 2 10 1000 10 10
1 0 3 8 1250 10
0.8 -0.6 4 10 1000 10

Paso 1
Qres Matriz Ks Matriz K (Af'*Ks*Af) U V=Af*U Q=Ks*V
10 1666.667 0 0 0 2530 0 0.007905 0.00419 6.98
10 0 1000 0 0 0 2386.667 0.00419 0.008838 8.84
10 0 0 1250 0 0.007905 9.88
10 0 0 0 1000 0.00381 3.81

Paso 2 (Barra 3 está fluyendo)

Qres Matriz Ks Matriz K (Af'*Ks*Af) U V=Af*U Q=Ks*V
2.933 1666.667 0 0 0 1280 0 0.015625 0.00419 6.98
1.056 0 1000 0 0 0 2386.667 0.00419 0.015014 15.01
0.000 0 0 0.0 0 0.015625 0.00
6.144 0 0 0 1000 0.009986 9.99

Paso 3 (Barra 2 y 3 están fluyendo)

Qres Matriz Ks Matriz K (Af'*Ks*Af) U V=Af*U Q=Ks*V
2.442 1666.667 0 0 0 640 -480 0.0425 0.015 25.00
0.000 0 0 0 0 -480 2026.667 0.015 0.043 0.00
0.000 0 0 0.0 0 0.0425 0.00
5.442 0 0 0 1000 0.025 25.00
 Resumen
Lambda U1
0.000 0
1.012 0.0080
1.082 0.0091
1.180 0.0133

Q/Qres Lamba dU dV dQ U V Q
0.698 1.012 0.0080 0.00424 7.067 0.0080 0.0042 7.0670
0.884 0.0042 0.00894 8.944 0.0042 0.0089 8.9441
0.988 0.00800 10.000 0.0080 10.0000
0.381 0.00386 3.856 0.0039 3.8559

Q/Qres Lamba dU dV dQ U V Q
2.381 0.070326 0.0011 0.00029 0.491 0.0091 0.0045 7.5581
14.220 0.0003 0.00106 1.056 0.0045 0.0100 10.0000
#DIV/0! 0.00110 0.000 0.0091 10.0000
1.625 0.00070 0.702 0.0046 4.5581

Q/Qres Lamba dU dV dQ U V Q
10.238 0.097674 0.0042 0.00147 2.442 0.0133 0.0060 10.0000
#DIV/0! 0.0015 0.00420 0.000 0.0060 0.0142 10.0000
#DIV/0! 0.00415 0.000 0.0133 10.0000
4.594 0.00244 2.442 0.0070 7.0000
Ejemplo 4.3 Marco Paso a Paso

EI 30000 Qp P
Af 0.2 0 0 0 0 150 30
0.2 1 0 0 0 150 0
0 1 -0.25 0 0 120 -50
0 0 -0.25 1 0 120 0
0 0 0.25 1 0 120 0
0 0 0.25 0 1 120
0.2 0 0 0 0 150
0.2 0 0 0 1 150

Paso 1
Qres Matriz Ks K=Af'*Ks*Af
150 24000 12000 0 0 0 0 0 0 5760
150 12000 24000 0 0 0 0 0 0 7200
120 0 0 30000 15000 0 0 0 0 0
120 0 0 15000 30000 0 0 0 0 0
120 0 0 0 0 30000 15000 0 0 7200
120 0 0 0 0 15000 30000 0 0
150 0 0 0 0 0 0 24000 12000
150 0 0 0 0 0 0 12000 24000

Paso 2 (Q6 está fluyendo)

Qres Matriz Ks K=Af'*Ks*Af
103.299 24000 12000 0 0 0 0 0 0 5760
134.951 12000 24000 0 0 0 0 0 0 7200
104.951 0 0 30000 15000 0 0 0 0 0
10.273 0 0 15000 30000 0 0 0 0 0
10.273 0 0 0 0 22500 0 0 0 7200
0.000 0 0 0 0 0 0 0 0
35.775 0 0 0 0 0 0 24000 12000
30.000 0 0 0 0 0 0 12000 24000

Paso 3 (Q4 y Q6 está fluyendo)

Qres Matriz Ks K=Af'*Ks*Af
95.455 24000 12000 0 0 0 0 0 0 5760
134.727 12000 24000 0 0 0 0 0 0 7200
104.727 0 0 22500 0 0 0 0 0 0
0.000 0 0 0 0 0 0 0 0 0
0.000 0 0 0 0 22500 0 0 0 7200
0.000 0 0 0 0 0 0 0 0
27.818 0 0 0 0 0 0 24000 12000
30.000 0 0 0 0 0 0 12000 24000

Paso 4 (Q4, Q6 y Q7 está fluyendo)

Qres Matriz Ks K=Af'*Ks*Af
80.000 24000 12000 0 0 0 0 0 0 3600
110.000 12000 24000 0 0 0 0 0 0 7200
80.000 0 0 22500 0 0 0 0 0 0
0.000 0 0 0 0 0 0 0 0 0
0.000 0 0 0 0 22500 0 0 0 3600
0.000 0 0 0 0 0 0 0 0
0.000 0 0 0 0 0 0 0 0
30.000 0 0 0 0 0 0 0 18000
K=Af'*Ks*Af U V=Af*U Q Q/Qres dlambda dU
7200 0 0 7200 0.0085 0.0017 26.35 0.176 1.773 0.0151
54000 -11250 15000 0 -0.0029 -0.0012 -8.49 0.057 -0.0051
-11250 11250 0 11250 -0.0076 -0.0010 8.49 0.071 -0.0135
15000 0 60000 15000 0.0007 0.0026 61.90 0.516 0.0012
0 11250 15000 54000 0.0003 -0.0012 -61.90 0.516 0.0005
-0.0016 -67.70 0.564
0.0017 64.44 0.430
0.0020 67.70 0.451

K=Af'*Ks*Af U V=Af*U Q Q/Qres dlambda dU

7200 0 0 7200 0.0213 0.0043 75.54 0.731 0.104 0.0022
54000 -11250 15000 0 -0.0065 -0.0022 -2.16 0.016 -0.0007
-11250 7031.25 -5625 0 -0.0175 -0.0021 2.16 0.021 -0.0018
15000 -5625 52500 0 0.0000 0.0043 98.92 9.629 0.0000
0 0 0 24000 -0.0064 -0.0044 -98.92 9.629 -0.0007
-0.0108 0.00 #DIV/0!
0.0043 76.62 2.142
-0.0021 0.00 0.000

K=Af'*Ks*Af U V=Af*U Q Q/Qres dlambda dU

7200 0 0 7200 0.0625 0.0125 125.00 1.310 0.124 0.0077
46500 -5625 0 0 -0.0271 -0.0146 -200.00 1.484 -0.0033
-5625 2812.5 5625 0 -0.1439 0.0089 200.00 1.910 -0.0178
0 5625 22500 0 0.0360 0.0719 0.00 #DIV/0! 0.0044
0 0 0 24000 -0.0188 0.0000 0.00 #DIV/0! -0.0023
-0.0547 0.00 #DIV/0!
0.0125 225.00 8.088
-0.0063 0.00 0.000
K=Af'*Ks*Af U V=Af*U Q Q/Qres dlambda dU
7200 0 0 3600 0.1250 0.0250 350.00 4.375 0.229 0.0286
46500 -5625 0 0 -0.0458 -0.0208 -200.00 1.818 -0.0105
-5625 2812.5 5625 0 -0.2189 0.0089 200.00 2.500 -0.0500
0 5625 22500 0 0.0547 0.1094 0.00 #DIV/0! 0.0125
0 0 0 18000 -0.0250 0.0000 0.00 #DIV/0! -0.0057
-0.0797 0.00 #DIV/0!
0.0250 0.00 #DIV/0!
0.0000 0.00 0.000
dV dQ U V Q
0.0030 46.7 0.0151 0.0030 46.70
-0.0021 -15.0 -0.0051 -0.0021 -15.05
-0.0018 15.0 -0.0135 -0.0018 15.05
0.0045 109.7 0.0012 0.0045 109.73
-0.0022 -109.7 0.0005 -0.0022 -109.73
-0.0029 -120.0 -0.0029 -120.00
0.0030 114.2 0.0030 114.22
0.0035 120.0 0.0035 120.00

dV dQ U V Q
0.0004 7.8 0.0173 0.0035 54.55
-0.0002 -0.2 -0.0058 -0.0024 -15.27
-0.0002 0.2 -0.0153 -0.0020 15.27
0.0005 10.3 0.0012 0.0050 120.00
-0.0005 -10.3 -0.0002 -0.0027 -120.00
-0.0011 0.0 -0.0040 -120.00
0.0004 8.0 0.0035 122.18
-0.0002 0.0 0.0033 120.00

dV dQ U V Q
0.0015 15.5 0.0250 0.0050 70.00
-0.0018 -24.7 -0.0092 -0.0042 -40.00
0.0011 24.7 -0.0331 -0.0009 40.00
0.0089 0.0 0.0056 0.0139 120.00
0.0000 0.0 -0.0025 -0.0027 -120.00
-0.0068 0.0 -0.0108 -120.00
0.0015 27.8 0.0050 150.00
-0.0008 0.0 0.0025 120.00
dV dQ U V Q
0.0057 80.0 0.0536 0.0107 150.00
-0.0048 -45.7 -0.0196 -0.0089 -85.71
0.0020 45.7 -0.0831 0.0011 85.71
0.0250 0.0 0.0181 0.0389 120.00
0.0000 0.0 -0.0082 -0.0027 -120.00
-0.0182 0.0 -0.0290 -120.00
0.0057 0.0 0.0107 150.00
0.0000 0.0 0.0025 120.00
Ejemplo 4.4. Reticulado con material no lineal

alpha 0.5
sy 10
ey 0.001

Af L A Pref
0 1 6 1 20
0.8 0.6 10 1 10
1 0 8 1
0.8 -0.6 10 1

Incremento de Carga
Lambda 0.55
P U(0) Matriz Ks Matriz K Pr=Bf*Q R(0)=P-Pr(0)
11 0.0000 833.3 0 0 0 1265.000 0.000 0.000 11.000
5.5 0.0000 0 500.0 0 0 0.000 1193.333 0.000 5.500
0 0 625.0 0
0 0 0 500.0

Estado de los elementos

Barra v=Af*U e=v/L sigma Et q k
1 0.00 0.0000 0.00 5000.00 0.00 833.3
2 0.00 0.0000 0.00 5000.00 0.00 500.0
3 0.00 0.0000 0.00 5000.00 0.00 625.0
4 0.00 0.0000 0.00 5000.00 0.00 500.0
dU(0) U(1)=U(0)+dU(0)
0.0087 0.0087
0.0046 0.0046