Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Rs-Using Newmark's Linear Acceleration Method

Download as xls, pdf, or txt
Download as xls, pdf, or txt
You are on page 1of 14

RESPONSE SPE

PERIOD DEFORMATION VELOCITY ACCELER


0.05 0.0005572486 0.01895149 0.897015
0.06 0.0008309989 0.02948597 0.928942
0.07 0.0012222307 0.04533449 1.003802
0.08 0.0018027641 0.07085377 1.133573 Deformation response spectru
0.09 0.0028284834 0.1210778 1.405269
0.8
0.1 0.0054556392 0.25656567 2.195515

D in m
0.7
0.11 0.0075298773 0.35105509 2.50434
0.6
0.12 0.0084800307 0.37005179 2.369876 0.5
0.13 0.0096776575 0.43909699 2.304487 0.4
0.14 0.0094483516 0.39971761 1.939951 0.3
0.15 0.0146883197 0.61477329 2.627123 0.2
0.16 0.0150210552 0.5945671 2.3613 0.1
0
0.17 0.0155358945 0.53860991 2.163362 0 0.5 1 1.5 2 2.5
0.18 0.0188767699 0.65811614 2.344625 Period (sec)
0.19 0.0189281858 0.63446444 2.110049
0.2 0.018196347 0.58416225 1.830691
0.21 0.0232848531 0.69854172 2.124837
0.22 0.0279681922 0.80498713 2.325465 Psedou-Velocity response spect
V in m/sec

0.23 0.0322107145 0.87767287 2.450391 2.5


0.24 0.033429449 0.84822987 2.335595
0.25 0.0327355406 0.78519368 2.107804 2
0.26 0.0361080827 0.70245389 2.149555 1.5
0.27 0.0363040455 0.7750892 2.004095
1
0.28 0.0408645786 0.82886003 2.097595
0.29 0.0451675822 0.92336361 2.161332 0.5
0.3 0.0484866388 0.9969195 2.168055 0
0.31 0.0501779315 1.0136732 2.101261 0 0.5 1 1.5 2 2.5
0.32 0.0491373408 0.95761546 1.93109 Period (sec)
0.33 0.0453422703 0.88435216 1.675584
0.34 0.0485586416 0.87407984 1.690439
0.35 0.0616438617 0.90418341 2.025091
0.36 0.0695641801 1.10927514 2.160088
0.37 0.0693498894 1.15806877 2.038605
0.38 0.0686381007 1.14528193 1.912885
0.39 0.0739821005 1.2222985 1.957439
0.4 0.0846547155 1.40617117 2.129227
0.41 0.0964846101 1.5578119 2.309836
0.42 0.1035887081 1.61072533 2.363222
0.43 0.1053555844 1.56096068 2.293039
0.44 0.1021214178 1.44755797 2.122766
0.45 0.0961908239 1.40778956 1.91161
0.46 0.0970834667 1.35664471 1.846377
0.47 0.0937313534 1.22562855 1.707575
0.48 0.0862862261 1.08377669 1.507127
0.49 0.0879349402 1.05086719 1.473873
0.5 0.0919726956 1.0275595 1.480504
0.51 0.0955109838 1.05244165 1.477759
0.52 0.0986578636 1.0774576 1.468303
0.53 0.1013406172 1.10971429 1.451853
0.54 0.1034621641 1.13845417 1.427857
0.55 0.1126441892 1.25887913 1.49856
0.56 0.131771089 1.42091295 1.690966
0.57 0.1489012676 1.55203606 1.844333
0.58 0.1623670379 1.63145081 1.942373
0.59 0.1711674749 1.64915812 1.978828
0.6 0.1751268367 1.61990074 1.957677
0.61 0.174059899 1.58606314 1.882477
0.62 0.1685988355 1.52344699 1.76507
0.63 0.1596585832 1.49432171 1.618832
0.64 0.148525376 1.43632754 1.459256
0.65 0.1357134805 1.34190848 1.292668
0.66 0.1263996466 1.27616279 1.167747
0.67 0.1458008227 1.2573287 1.307077
0.68 0.1665419648 1.41969416 1.449428
0.69 0.1845677156 1.59460615 1.560085
0.7 0.2000207005 1.73571505 1.642743
0.71 0.2138397438 1.84809159 1.707114
0.72 0.2273296475 1.94644679 1.764744
0.73 0.2435426174 2.03230338 1.839162
0.74 0.2581446646 2.11613184 1.897101
0.75 0.2687673663 2.17064353 1.922847
0.76 0.2738073839 2.23140148 1.907694
0.77 0.2717515982 2.24127903 1.844512
0.78 0.2636964696 2.19004534 1.744238
0.79 0.2505370026 2.08388135 1.615505
0.8 0.2354110353 1.94513582 1.480259
0.81 0.2207719532 1.81002457 1.354143
0.82 0.2091665341 1.697146 1.251858
0.83 0.2022239674 1.62050713 1.181319
0.84 0.2000972925 1.57848941 1.14123
0.85 0.2015094215 1.55854013 1.122401
0.86 0.2050009953 1.54797013 1.115449
0.87 0.2085943656 1.5723122 1.109059
0.88 0.211326427 1.59294658 1.098194
0.89 0.2129436305 1.60596114 1.081871
0.9 0.213760815 1.60899337 1.062023
0.91 0.2142675362 1.60771015 1.041272
0.92 0.2151575442 1.61091785 1.022991
0.93 0.217411341 1.62037463 1.011596
0.94 0.2207808856 1.63864441 1.005533
0.95 0.2252367193 1.66583159 1.004345
0.96 0.2302763665 1.69976898 1.005536
0.97 0.2349517911 1.73095017 1.004907
0.98 0.2384967388 1.75670916 0.999358
0.99 0.2406235762 1.76837344 0.988004
1 0.2467234318 1.76514735 0.99289
1.01 0.2505235348 1.74040609 0.988318
1.02 0.2518736876 1.72991239 0.974256
1.03 0.2509240019 1.76824029 0.951828
1.04 0.2510470557 1.78818508 0.934069
1.05 0.2550694512 1.78899336 0.931045
1.06 0.2573063776 1.77458864 0.921573
1.07 0.2578133765 1.74055576 0.90621
1.08 0.2566633171 1.69304292 0.885538
1.09 0.2539298361 1.66774896 0.860105
1.1 0.2496733133 1.66049789 0.830381
1.11 0.2488679141 1.64990269 0.812856
1.12 0.248924114 1.63412684 0.798586
1.13 0.2476919339 1.63603253 0.780631
1.14 0.2545370775 1.64772255 0.788192
1.15 0.2617446882 1.6575966 0.796476
1.16 0.2690877442 1.6652103 0.804764
1.17 0.2763413811 1.67189869 0.812391
1.18 0.2836965848 1.67674169 0.819938
1.19 0.290936518 1.67783485 0.82679
1.2 0.2977982539 1.67470921 0.832243
1.21 0.3046325144 1.66697208 0.837329
1.22 0.3109519588 1.65432858 0.840745
1.23 0.3168306815 1.63677185 0.842767
1.24 0.322301823 1.61717152 0.843549
1.25 0.3270765904 1.59236343 0.842404
1.26 0.3314610645 1.56248951 0.840199
1.27 0.3350879105 1.54827303 0.836069
1.28 0.3382354625 1.53902497 0.830787
1.29 0.3407482513 1.52715941 0.824034
1.3 0.3426030413 1.51334041 0.815822
1.31 0.3440799619 1.49644984 0.806877
1.32 0.3449117752 1.47675022 0.796619
1.33 0.3452099824 1.45454072 0.785364
1.34 0.3451327416 1.45900899 0.773512
1.35 0.3444263481 1.46364382 0.760536
1.36 0.3434289384 1.46177829 0.747222
1.37 0.3418897681 1.45521756 0.733054
1.38 0.3400708799 1.44185386 0.718624
1.39 0.3377985122 1.42234581 0.703589
1.4 0.3353785968 1.3973894 0.688605
1.41 0.3392756283 1.39026067 0.68676
1.42 0.3436395186 1.40045254 0.685831
1.43 0.3475586377 1.40807654 0.683985
1.44 0.3510309823 1.41340604 0.681257
1.45 0.3540652094 1.41672222 0.677701
1.46 0.3567736068 1.4275593 0.673562
1.47 0.3592445779 1.4907643 0.669031
1.48 0.3613594289 1.54872879 0.663906
1.49 0.3632483213 1.60257136 0.658448
1.5 0.3650471138 1.65231676 0.652916
1.51 0.3666009661 1.69984097 0.647039
1.52 0.3681413812 1.74684794 0.641236
1.53 0.3696379945 1.79179097 0.635455
1.54 0.371062108 1.83357414 0.629645
1.55 0.3813766053 1.87213834 0.638824
1.56 0.3931120801 1.90644484 0.650067
1.57 0.4045815929 1.93859262 0.660538
1.58 0.415619315 1.96745916 0.669996
1.59 0.4264094644 1.99245425 0.678771
1.6 0.4371021589 2.0136706 0.687122
1.61 0.4477182374 2.03325016 0.695094
1.62 0.4583522619 2.04983006 0.702846
1.63 0.4689199766 2.06279926 0.710255
1.64 0.4794127059 2.07217574 0.717319
1.65 0.4897779331 2.07796362 0.723972
1.66 0.4998632171 2.08275483 0.730005
1.67 0.5096698686 2.08400862 0.735439
1.68 0.5193060379 2.08161205 0.74045
1.69 0.5284668648 2.07554354 0.744621
1.7 0.5374548408 2.06578619 0.748402
1.71 0.5459800128 2.05475924 0.751407
1.72 0.5541106477 2.04014012 0.753755
1.73 0.5619222348 2.02179221 0.75557
1.74 0.5691659077 1.99974409 0.756539
1.75 0.5761213166 1.97492086 0.757057
1.76 0.5825808693 1.94796574 0.756871
1.77 0.5884753131 1.91743857 0.755914
1.78 0.5940425065 1.89067596 0.754516
1.79 0.5990906077 1.87314903 0.752449
1.8 0.6035628302 1.85305041 0.749667
1.81 0.6075694709 1.89150806 0.746328
1.82 0.6111333598 1.93270774 0.742479
1.83 0.6141015483 1.97152994 0.737953
1.84 0.6164659163 2.00786144 0.732764
1.85 0.6184538253 2.04160388 0.727201
1.86 0.6198475983 2.07267481 0.721024
1.87 0.6206030176 2.10100839 0.714203
1.88 0.6208924423 2.12655578 0.706954
1.89 0.6206421173 2.14928542 0.699211
1.9 0.6197131654 2.16918296 0.690835
1.91 0.6183728677 2.18625097 0.682141
1.92 0.6164360836 2.20050856 0.67294
1.93 0.6138866395 2.21271925 0.66323
1.94 0.610987415 2.22324572 0.65331
1.95 0.6074463728 2.23108462 0.642879
1.96 0.603542832 2.23630819 0.632247
1.97 0.5991810057 2.23900113 0.621321
1.98 0.5943260582 2.23925944 0.610078
1.99 0.5892705295 2.23718922 0.598824
2 0.5837830302 2.23290542 0.58733
2.01 0.5779825505 2.22653062 0.575723
2.02 0.5720140102 2.21819373 0.56415
2.03 0.5657292213 2.20802878 0.552468
2.04 0.5592487955 2.19617372 0.540798
2.05 0.5526328375 2.18276921 0.5292
2.06 0.5457410457 2.16795753 0.517539
2.07 0.5388489741 2.1518815 0.506078
2.08 0.5317212195 2.13468347 0.494593
2.09 0.5245880862 2.11650437 0.4833
2.1 0.5230955769 2.09748289 0.477346
2.11 0.5232193778 2.07775463 0.472944
2.12 0.523708909 2.05784839 0.468931
2.13 0.5244348451 2.03790199 0.465182
2.14 0.5255046207 2.01757317 0.461785
2.15 0.5269331954 1.99698047 0.458743
2.16 0.5284631945 1.97623671 0.455825
2.17 0.5299635101 1.9554487 0.452916
2.18 0.5313618093 1.93471686 0.449954
2.19 0.5327000411 1.91413507 0.446977
2.2 0.5341412076 1.89379044 0.444121
2.21 0.5357382157 1.87376323 0.441427
2.22 0.5375196164 1.8541268 0.438914
2.23 0.5394308375 1.83494761 0.436533
2.24 0.5412929001 1.81628523 0.434137
2.25 0.543069813 1.79819246 0.431699
2.26 0.544829094 1.78071546 0.429274
2.27 0.5464827759 1.76389388 0.426791
2.28 0.5480296055 1.74776108 0.424253
2.29 0.5495964222 1.7323443 0.421758
2.3 0.5511535912 1.71766497 0.419283
2.31 0.5526943995 1.70373891 0.416823
2.32 0.5542765637 1.69057662 0.414421
2.33 0.5560014993 1.67818357 0.41215
2.34 0.5578392986 1.6665605 0.409985
2.35 0.5598357775 1.65570372 0.407958
2.36 0.5620228843 1.64560542 0.406089
2.37 0.5644150119 1.63625395 0.404383
2.38 0.5670098011 1.62763418 0.402835
2.39 0.5697899675 1.61972777 0.40143
2.4 0.5727253024 1.61251349 0.400142
2.41 0.5757757345 1.6059675 0.398942
2.42 0.5789189419 1.60006365 0.397812
2.43 0.5822439247 1.59477378 0.39681
2.44 0.5855855764 1.59006794 0.395823
2.45 0.5891007526 1.58591471 0.394955
2.46 0.5926120356 1.58228143 0.394086
2.47 0.5962869553 1.57913443 0.393326
2.48 0.5999154147 1.57643926 0.392534
2.49 0.6037319602 1.57416094 0.391865
2.5 0.6074912683 1.57226409 0.391157
2.51 0.6113620049 1.57071321 0.390519
2.52 0.615221181 1.56947279 0.389871
2.53 0.6191072641 1.56850752 0.389238
2.54 0.6230285329 1.5677824 0.388625
2.55 0.6269053338 1.56726293 0.387983
2.56 0.6308457531 1.56691519 0.387377
2.57 0.6346994599 1.566706 0.386716
2.58 0.6386124642 1.56660299 0.38609
2.59 0.6424331586 1.56657472 0.385407
2.6 0.6462699956 1.56659075 0.384732
2.61 0.6500425705 1.56662174 0.384018
2.62 0.6537536127 1.56663944 0.383268
2.63 0.6574495306 1.56661685 0.382509
2.64 0.6609921769 1.56652817 0.381662
2.65 0.6645576347 1.56634891 0.38083
2.66 0.6679323992 1.56605586 0.379892
2.67 0.6712538223 1.56562715 0.378926
2.68 0.6743835036 1.56504226 0.377858
2.69 0.6774163228 1.56428199 0.37674
2.7 0.6802235537 1.56332854 0.375504
2.71 0.6829264233 1.56216541 0.374219
2.72 0.6853465392 1.56077746 0.372789
2.73 0.6876818876 1.5591509 0.371324
2.74 0.6896771024 1.5572732 0.369688
2.75 0.691608819 1.55513315 0.368032
2.76 0.6931871561 1.55272078 0.366204
2.77 0.6946679533 1.55002738 0.364341
2.78 0.6958032726 1.54704542 0.362316
2.79 0.6968535808 1.54376852 0.360266
2.8 0.6976076597 1.54019148 0.358085
2.81 0.6981850513 1.53631014 0.355835
2.82 0.6986061094 1.53212143 0.353529
2.83 0.698840777 1.52762328 0.351153
2.84 0.6989243241 1.52281458 0.348726
2.85 0.69886929 1.51769518 0.346256
2.86 0.698666679 1.51226579 0.343739
2.87 0.6983297691 1.50652797 0.341183
2.88 0.6978443266 1.50048407 0.338582
2.89 0.6971562211 1.50307407 0.335912
2.9 0.6962395707 1.50744252 0.33316
2.91 0.6951664243 1.5111071 0.330365
2.92 0.6938919729 1.51406723 0.327504
2.93 0.6923699517 1.51632392 0.324559
2.94 0.6906770873 1.52020058 0.321567
2.95 0.6888042293 1.52600087 0.318524
2.96 0.686707842 1.53107738 0.315413
2.97 0.684410753 1.53543185 0.312244
2.98 0.6820010417 1.53906746 0.30906
2.99 0.6794013555 1.54198877 0.305826
3 0.676623341 1.54453811 0.302549
ESPONSE SPECTRUM-WORKBOOK
N.B Select IE/ LA time stepping methods
NEWMARK`S LINEAR ACCELERATION METHOD
3

A /g
2.5

1.5

sponse spectrum
1

0.5

0
0 0.5 1 1.5 2 2. 5 3 3 .5
Peri od( sec)

5 2 2.5 3 3.5
eriod (sec)

response spectrum

5 2 2.5 3 3.5
eriod (sec)
THOD

3 3. 5
If Worksheets("sheet1").Cells(4, 12) = "METHOD BASED ON INTRPOLATION EXITATION" Then

For T = 0.05 To 3 Step 0.01


u(0) = 0
uu(0) = 0
p(0) = 0
Tno = (100 * T) - 1
'Cooficients
Wn = (2 * pi) / T
k = Wn ^ 2
WD = (Wn) * Sqr(1 - da ^ 2)
At = MAT(2, 1) - MAT(1, 1)
v = (2.718281828) ^ (-da * Wn * At)
si = Sin(WD * At)
co = Cos(WD * At)

'FINDING THE COOFICIENTS


a = v * ((da / Sqr(1 - da ^ 2)) * si + co)
b = v * ((1 / WD) * si)
E = 2 * da / (Wn * At)
F = (1 - 2 * da ^ 2) / (WD * At)
G = da / (Sqr(1 - da ^ 2))
H = 1 + 2 * da / (Wn * At)
c = (1 / k) * (E + v * ((F - G) * si - H * co))
D = (1 / k) * (1 - ((2 * da) / (Wn * At)) + v * (((-1 + 2 * da ^ 2) / (WD * At)) * si + ((2 * da) / (Wn * At)) * co

AA = -v * ((Wn * si) / Sqr(1 - da ^ 2))


BB = v * (co - ((da * si) / Sqr(1 - da ^ 2)))
EE = 1 / At
FF = Wn / Sqr(1 - da ^ 2)
GG = da / (At * (Sqr(1 - da ^ 2)))
CC = (1 / k) * ((-1 / At) + v * ((Wn / Sqr(1 - da ^ 2) + da / (At * Sqr(1 - da ^ 2))) * si + (co / At)))
DD = (1 / (k * At)) * (1 - v * (((da * si) / Sqr(1 - da ^ 2)) + co))

'The responses
For x = 1 To N

p(x) = MAT(x, 2)
u(x) = a * u(x - 1) + b * uu(x - 1) + c * p(x - 1) + D * p(x)
uu(x) = AA * u(x - 1) + BB * uu(x - 1) + CC * p(x - 1) + DD * p(x)

Next x

'FINDING THE MAXIMUM DISPLACEMENT AND VELOCITY


MaxD = 0
MaxV = 0
For J = 1 To N
If Abs(u(J)) >= MaxD Then
MaxD = Abs(u(J))
Else: MaxD = MaxD
End If

If Abs(uu(J)) >= MaxV Then


MaxV = Abs(uu(J))
Else: MaxV = MaxV
End If

Next J

Worksheets("sheet1").Cells(Tno, 1) = T
Worksheets("sheet1").Cells(Tno, 2) = MaxD
Worksheets("sheet1").Cells(Tno, 3) = MaxV
Worksheets("sheet1").Cells(Tno, 4) = (MaxD * Wn ^ 2 / 9.81) '/ (2 * pi)

Next T

Else

For T = 0.05 To 3 Step 0.01


u(0) = 0
uu(0) = 0
uuu(0) = 0
p(0) = 0

Tno = (100 * T) - 1
'Cooficients
Wn = (2 * pi) / T
k = Wn ^ 2
At = MAT(2, 1) - MAT(1, 1)
beta = 1 / 6
gama = 0.5
c = 2 * da * Wn
a = (1 / (beta * At)) + (gama / beta) * c
b = (1 / (2 * beta)) + At * c * ((gama / (2 * beta)) - 1)
kpe = k + ((gama / (beta * At)) * c) + (1 / (beta * (At) ^ 2))

'The responses
For x = 1 To N
p(x) = -(MAT(x + 1, 2) - MAT(x, 2)) + a * uu(x - 1) + b * uuu(x - 1)
deu = p(x) / kpe
deuu = (gama / (beta * At)) * deu - (gama / beta) * uu(x - 1) + At * (1 - (gama / (2 * beta))) * uuu(x -
deuuu = (1 / (beta * (At) ^ 2)) * deu - (1 / (beta * At)) * uu(x - 1) - (1 / (2 * beta)) * uuu(x - 1)
u(x) = u(x - 1) + deu
uu(x) = uu(x - 1) + deuu
uuu(x) = uuu(x - 1) + deuuu

Next x

'FINDING THE MAXIMUM DISPLACEMENT ,VELOCITY AND ACCELERATION

MaxD = 0
MaxV = 0
MaxA = 0
For J = 1 To N

If Abs(u(J)) >= MaxD Then


MaxD = Abs(u(J))
Else: MaxD = MaxD
End If

If Abs(uu(J)) >= MaxV Then


MaxV = Abs(uu(J))
Else: MaxV = MaxV
End If

If Abs(uuu(J)) >= MaxA Then


MaxA = Abs(uuu(J))
Else: MaxA = MaxA
End If

Next J

Worksheets("sheet1").Cells(Tno, 1) = T
Worksheets("sheet1").Cells(Tno, 2) = MaxD
Worksheets("sheet1").Cells(Tno, 3) = MaxV
Worksheets("sheet1").Cells(Tno, 4) = MaxA / 9.81

Next T

End If
Exit Sub
EXITATION" Then

i + ((2 * da) / (Wn * At)) * co))

))) * si + (co / At)))


ama / (2 * beta))) * uuu(x - 1)
beta)) * uuu(x - 1)

You might also like