Numerical Solution Using Linear Interpolation of Excitation

Tn= 1 ᶓ WnΔt= 0.031415 K= 10 K/in

Wn= 6.283 sin WdΔt= 0.5871343 mass= 0.2533 kip-sec2/in
ᶓ= 0.05 Cos WdΔt=0.8094895 2ᶓ /WdΔt= 0.159359
Δt= 0.1 sqrt(1-ᶓ2) 0.9987492 2ᶓ /WnΔt=0.1591596 Chart Title
Wd= 6.2751413 WdΔt= 0.6275141 c= 0.1592 kip-sec/in
2

Constants 1.5

A= 0.8129392 A' = -3.579354 1

B= 0.0906715 B' = 0.7559703 0.5
C= 0.0123548 C' = 0.17087 0
D= 0.0063513 D' = 0.18706 0 0.2 0.4 0.6 0.8 1 1.2
-0.5
-1
i ti pi ui u'i Aui Bu'i Cpi Dpi+1 A'ui B'u'i C'pi D'pi+1
-1.5
0 0 0 0 0 0 0 0 0.0318 0 0 0 0.9353
1 0.1 5 0.0318 0.9353 0.0259 0.0848 0.0618 0.055 -0.1138 0.7071 0.8544 1.6199
2 0.2 8.66 0.2275 3.0676 0.1849 0.2781 0.107 0.0635 -0.8143 2.319 1.4798 1.8706
3 0.3 10 0.6335 4.8551 0.515 0.4402 0.1235 0.055 -2.2675 3.6703 1.7087 1.6199
4 0.4 8.66 1.1337 4.7314 0.9216 0.429 0.107 0.0318 -4.0579 3.5768 1.4798 0.9353
5 0.5 5 1.4894 1.934 1.2108 0.1754 0.0618 0 -5.3311 1.462 0.8544 0
6 0.6 0 1.448 -3.0147 1.1771 -0.2733 0 0 -5.1829 -2.279 0 0
7 0.7 0 0.9038 -7.4619 0.7347 -0.6766 0 0 -3.235 -5.641 0 0
8 0.8 0 0.0581 -8.876 0.0472 -0.8048 0 0 -0.208 -6.71 0 0
9 0.9 0 -0.7576 -6.918 -0.6159 -0.6273 0 0 2.7117 -5.2298 0 0
10 1 0 -1.2432 -2.5181 -1.0106 -0.2283 0 0 4.4499 -1.9036 0 0
 Numerical Solution By Central Difference Method

Δt= 0.1 c= 0.1592

K= 10 K/in K~= 26.126
mass= 0.2533 kip-sec2/in a= 24.534
Wn= 6.283 b= -40.66
ᶓ= 0.05
Chart Title
2
i ti pi ui-1 ui pi~ ui+1
1.5
0 0 0 0 0 0 0
1 0.1 5 0 0 5 0.19138 1
2 0.2 8.66 0 0.19138 16.44152 0.629316
0.5
3 0.3 10 0.19138 0.629316 30.89268 1.18245
4 0.4 8.66 0.629316 1.18245 41.29876 1.580753 0
0 0.2 0.4 0.6 0.8
5 0.5 5 1.18245 1.580753 40.2632 1.541116
-0.5
6 0.6 0 1.580753 1.541116 23.87959 0.914016
7 0.7 0 1.541116 0.914016 -0.645847 -0.02472 -1
8 0.8 0 0.914016 -0.02472 -23.42961 -0.896793
-1.5
9 0.9 0 -0.02472 -0.896793 -35.8571 -1.372468
10 1 0 -0.896793 -1.372468 -33.80264 -1.293831 -2
 Chart Title

0.4 0.6 0.8 1 1.2

 Newmark method

γ= 0.5
β= 0.166667

Δt= 0.08 c= 0.4

K= 8 K/in K~ = 491.75
mass= 0.5 kip-sec2/in a= 38.7
Wn= 4 b= 1.516
ᶓ= 0.1

i ti pi del pi ui u'i u"i del p~i del ui del u'i del u"i
0 0 0 1 0 0 0 1 0.002034 0.076258 1.906457
1 0.08 1 3 0.002034 0.076258 1.906457 8.841383 0.017979 0.369195 5.416973
2 0.16 4 5 0.020013 0.445454 7.323429 33.34138 0.067801 0.913257 8.184571
3 0.24 9 0 0.087814 1.358711 15.508 76.09224 0.154738 1.106209 -3.36077
4 0.32 9 -3 0.242552 2.46492 12.14723 110.8076 0.225333 0.569346 -10.06081
5 0.4 6 -6 0.467885 3.034266 2.086422 114.5891 0.233023 -0.447889 -15.37006
6 0.48 0 0 0.700908 2.586378 -13.28364 79.95482 0.162592 -1.130572 -1.697021
 Combine Displacement Graph
1.6

1.4

1.2

1
Displacement

0.8

0.6

0.4

0.2

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7

Time interval

Numerical Solution Using Linear Interpolation of Excitation