CALCULO NUMERICO DE INTEGRAL DE DUHAMEL - SIN AMORTIGUAMIENTO

A. DUHAMEL - SIN AMORTIGUAMIENTO

DATOS: 1. Determinacion de la masa:
K(Regidez) 50000000 N/m m 12000 Kg
P(Peso) 120000 KN 2. Determinacion de frecuencia angular:
g(gravedad) 10 m/s^2 w 64.5497224 rad/s

3. Calculo de respuesta:

t(s) P(T)(N) ∆Pi ∆t Ѡt ∆A(t) A(t) ∆B(t) B(t)

0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
0.20000 50000.00000 50000.00000 0.20000 12.90994 257.41958 257.41958 -709.11551 -709.11551
0.40000 100000.00000 50000.00000 0.20000 25.81989 711.67118 969.09076 -450.44139 -1159.55690
0.60000 150000.00000 50000.00000 0.20000 38.72983 994.80684 1963.89761 16.12414 -1143.43276
0.80000 100000.00000 -50000.00000 0.20000 51.63978 -454.52764 1509.36997 885.04535 -258.38742
1.00000 50000.00000 -50000.00000 0.20000 64.54972 -732.63904 776.73092 415.46763 157.08021
1.20000 0.00000 -50000.00000 0.20000 77.45967 -746.75512 29.97580 -107.08147 49.99874
1.40000 0.00000 0.00000 0.20000 90.36961 0.00000 29.97580 0.00000 49.99874
1.60000 0.00000 0.00000 0.20000 103.27956 0.00000 29.97580 0.00000 49.99874
1.80000 0.00000 0.00000 0.20000 116.18950 0.00000 29.97580 0.00000 49.99874
2.00000 0.00000 0.00000 0.20000 129.09944 0.00000 29.97580 0.00000 49.99874
2.20000 0.00000 0.00000 0.20000 142.00939 0.00000 29.97580 0.00000 49.99874
2.40000 0.00000 0.00000 0.20000 154.91933 0.00000 29.97580 0.00000 49.99874
TO

U(t)
0.00000
0.00097
0.00195
0.00293
0.00198
0.00102
0.00006
0.00007
0.00007
0.00007
0.00005
0.00003
0.00000
 EVALUACION NUMERICA DE LA INTEGRAL DE DUHAMEL

B. DUHAMEL - CON AMORTIGUAMIENTO

DATOS: 1.- Determinacion de masa, frecuencia angular natural y frecuencia forzada.
Regidez(K) 10 KN/m Tipo Valor
Peso(W) 97.7738 KN m 0.2533 Kg ʓ 0.05 Ѡd=Ѡ*raiz(1-ʓ^2)
Gravedad(g) 386 m/s^2 Ѡn=Ѡ 6.283222007843 rad/s
C 4 Ѡd 6.2753630654204 62.3432748517048

2.- calculo de integrales

t(s) ʓѠ Ѡ*t Ѡd*t I1 I2 I3 I4 A B C D P Q I1

0 0.3142 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0253 0.3142 -6.2754 -0.0080 0.1590 0.6291 0.0080
0.02 0.3142 0.1257 0.1255 0.0200 0.0013 0.0032 0.0000 0.0255 1.0972 -6.1867 0.0120 0.1590 0.6252
0.04 0.3142 0.2513 0.2510 0.0005 0.0001 0.0000 0.0000 0.0257 1.8630 -6.0007 0.0320 0.1590 0.6213
0.06 0.3142 0.3770 0.3765 0.0005 0.0002 0.0000 0.0001 0.0258 2.5995 -5.7203 0.0520 0.1590 0.6174
0.08 0.3142 0.5027 0.5020 0.0005 0.0002 0.0000 0.0002 0.0260 3.2951 -5.3499 0.0720 0.1590 0.6135
0.1 0.3142 0.6283 0.6275 0.0004 0.0003 0.0000 0.0003 0.0261 3.9389 -4.8953 0.0920 0.1590 0.6097
1.2 0.3142 7.5399 7.5304 0.0031 0.0017 -0.0011 0.0067 0.0369 6.0496 -1.6973 1.1920 0.1590 0.4315
1.4 0.3142 8.7965 8.7855 -0.0022 0.0068 0.0088 0.0103 0.0393 3.4918 5.2236 1.3920 0.1590 0.4052
1.6 0.3142 10.0532 10.0406 -0.0076 0.0000 -0.0001 -0.0010 0.0419 -3.8812 4.9412 1.5920 0.1590 0.3806
1.8 0.3142 11.3098 11.2957 -0.0025 -0.0077 -0.0131 -0.0164 0.0446 -5.9021 -2.1551 1.7920 0.1590 0.3574
2 0.3142 12.5664 12.5507 0.0069 -0.0051 -0.0096 -0.0110 0.0475 0.2160 -6.2795 1.9920 0.1590 0.3356
2.2 0.3142 13.8231 13.8058 0.0075 0.0053 0.0113 0.0151 0.0506 6.0362 -1.7445 2.1920 0.1590 0.3152

3.- calculo de respuesta

CON.AMR. SIN. AMR.
t(s) P(T)(kgf/m) ∆Pi ∆t ∆Ad(t) Ad(t) ∆Bd(t) Bd(t) U(t)
0.0 0.0 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.02 120000 120000.0000 0.0200 245.3691 245.3691 19175.3148 19175.3148 -1501.5839 -150158.3853 0.0974
0.04 120000 0.0000 0.0200 60.2432 305.6122 11.4814 19186.7963 -2926.6120 -292661.1964 0.1951
0.06 0 -120000.0000 0.0200 -359.8854 -54.2732 8.9017 19195.6980 -4238.4586 -423845.8618 0.2934
0.08 0 0.0000 0.0200 0.0000 -54.2732 0.0000 19195.6980 -5329.4354 -532943.5352 0.1976
0.1 0 0.0000 0.0200 0.0000 -54.2732 0.0000 19195.6980 -6154.6913 -615469.1293 0.1022
1.2 0 0.0000 1.1000 0.0000 -54.2732 0.0000 19195.6980 -29192.5599 -2919255.9918 0.0065
1.4 0 0.0000 0.2000 0.0000 -54.2732 0.0000 19195.6980 85072.4023 8507240.2255 0.0074
1.6 0 0.0000 0.2000 0.0000 -54.2732 0.0000 19195.6980 99301.3864 9930138.6367 0.0074
1.8 0 0.0000 0.2000 0.0000 -54.2732 0.0000 19195.6980 -40003.3677 -4000336.7671 0.0066
2 0 0.0000 0.2000 0.0000 -54.2732 0.0000 19195.6980 -151729.4562 -15172945.6189 0.0051
2.2 0 0.0000 0.2000 0.0000 -54.2732 0.0000 19195.6980 -54754.8868 -5475488.6846 0.0029
0.0004
 I2 I3
-0.15895495 0.00252983
 EVALUACION NUMERICA DE LA INTEGRAL DE DUHA

B. DUHAMEL - CON AMORTIGUAMIENTO

DATOS: 1.- Determinacion de masa, frecuencia angular nat
Regidez(K) 100000 KN/m Tipo Valor
Peso(W) 98000 KN m 100
Gravedad(g) 980 m/s^2 Ѡn=Ѡ 31.622776601684
C Ѡd 30.983866769659

2.- calculo de integrales

2.718281828459
t(s) ʓѠ Ҽ^(-ʓѠt) Ѡd*t I1
0 6.32455532 1 0 0
0.02 6.32455532 0.881181986383 0.619677335 0.0199392398282
0.04 6.32455532 0.776481693126 1.239354671 0.0039536908282
0.06 6.32455532 0.684221680738 1.859032006 -0.00525188352
0.08 6.32455532 0.602923819759 2.478709342 -0.010604742911
0.1 6.32455532 0.531285609133 3.098386677 -0.011136425422
1.2 6.32455532 0.000505750859 37.18064012 0.0026410177963
1.4 6.32455532 0.000142755462 43.37741348 3.214493569E-06
1.6 6.32455532 4.029478457E-05 49.57418683 1.172992621E-06
1.8 6.32455532 1.137378313E-05 55.77096019 4.036084045E-07
2 6.32455532 3.210414053E-06 61.96773354 1.335423559E-07
2.2 6.32455532 9.061855915E-07 68.16450689 4.295044733E-08

3.- calculo de respuesta

t(s) P(T)(kgf/m) ∆Pi ∆t ∆Ad(t) Ad(t)

0.0 0.0 0 0 0 0
0.02 120000 120000.0 0.0 -1372.90772 -1372.907717604
0.04 120000 0 0.02 474.4428994 -898.4648182207
0.06 0 -120000.0 0.02 1817.335791 918.87097237067
0.08 0 0 0.02 0 918.87097237067
0.1 0 0.0 0.02 0 918.87097237067
1.2 0 0 1.1 0 918.87097237067
1.4 0 0.0 0.2 0 918.87097237067
1.6 0 0 0.2 0 918.87097237067
1.8 0 0.0 0.2 0 918.87097237067
2 0 0 0.2 0 918.87097237067
2.2 0 0.0 0.2 0 918.87097237067
LA INTEGRAL DE DUHAMEL

masa, frecuencia angular natural y frecuencia forzada. Ѡd=Ѡ*raiz(1-ʓ^2)

Kg ʓ 0.2
rad/s
0

I2 I3 I4
0 0 0 1
0.006528326833 0.0007070725228 -0.000228817953 1.13483936 0.0007070725 0.1262990685
-0.012990282813 -0.001022025443 -0.101715924442
-0.014349793939 -0.00061800231 -0.058258142422 0.0007070725 -0.000228818
-0.010500651353 -0.000169260869 -0.033727405876
-0.003922726418 0.0003897034362 -0.003150941313
0.016312812737 0.0202665456042 0.0877831010816
-8.90072457E-06 -0.019566338065 -0.017209712558
-2.32549883E-06 8.635377292E-06 -1.269937076E-05
-5.98763584E-07 2.608264347E-06 -4.810282856E-06
-1.5147816E-07 7.6362156E-07 -1.693611585E-06
-3.748931E-08 2.169520214E-07 -5.692745629E-07

CON.AMR. SIN. AMR.

∆Bd(t) Bd(t) U(t)
0 0 0 0.00000 0.00000
4242.435136649 4242.4351366494 -1.208976653166 -120.89767 0.09739
-1558.833937542 2683.6011991075 -0.431754794705 -43.17548 0.19509
344382.9287116 347066.52991068 21.981243618699 2198.12436 0.29336
0 347066.52991068 53.343856502173 5334.38565 0.19762
0 347066.52991068 59.463352263048 5946.33523 0.10216
0 347066.52991068 -0.049280811534 -4.92808 0.00645
0 347066.52991068 -0.01317758928 -1.31776 0.00738
0 347066.52991068 -0.003485012326 -0.34850 0.00744
0 347066.52991068 -0.000910148527 -0.09101 0.00664
0 347066.52991068 -0.000234225897 -0.02342 0.00506
0 347066.52991068 -5.92194618E-05 -0.00592 0.00288
0.00037
 #NAME?

Ѡ*raiz(1-ʓ^2)

#NAME?