Rs-Using Newmark's Linear Acceleration Method
Rs-Using Newmark's Linear Acceleration Method
Rs-Using Newmark's Linear Acceleration Method
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
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
'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
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
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
MaxD = 0
MaxV = 0
MaxA = 0
For J = 1 To N
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