-
Notifications
You must be signed in to change notification settings - Fork 0
/
info_Comp_Fragility_NonStructural_Accel.m
97 lines (63 loc) · 4.13 KB
/
info_Comp_Fragility_NonStructural_Accel.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
function [PDS_ij_EDP, xm_Cost, numCompPerStory] = info_Comp_Fragility_NonStructural_Accel(i_n, i_m, x_Accel_pdf)
% This function file returns from story number IDs (for each story):
% n = number of damage states a component may experience (this is per m)
% i_m = ID of component
% % % % % % % % % % % % % % % %
if i_m == 1
% % % % % % % % % % % % % % % %
numCompPerStory = 21*(9.144*9.144) /232.0; % total area divided by 232 m^2
if i_n == 0 % Suspended Ceiling (corrected on 30Jan2021)
% xm_EDP=1.09; beta_EDP=0.30; xm_Cost=0.% 3542.;
xm_EDP=0.56; beta_EDP=0.25; xm_Cost=0.;
F_DS_ij = normcdf((log(x_Accel_pdf/xm_EDP))/beta_EDP); % compute fragility function using Eq. 1 and estimated parameters
PDS_ij_EDP = 1.0 - F_DS_ij; % j=0, i.e. no damage
elseif i_n == 1
% xm_EDP=1.09; beta_EDP=0.30; xm_Cost=3542.;
xm_EDP=0.56; beta_EDP=0.25; xm_Cost=3542.;
F_DS_i1 = normcdf((log(x_Accel_pdf/xm_EDP))/beta_EDP); % compute fragility function using Eq. 1 and estimated parameters
% xm_EDP=1.69; beta_EDP=0.30;
xm_EDP=1.08; beta_EDP=0.25;
F_DS_i2 = normcdf((log(x_Accel_pdf/xm_EDP))/beta_EDP); % compute fragility function using Eq. 1 and estimated parameters
PDS_ij_EDP = F_DS_i1 - F_DS_i2; % 1<=j<=n, i.e. some damage
elseif i_n == 2
% xm_EDP=1.69; beta_EDP=0.30; xm_Cost=29337.;
xm_EDP=1.08; beta_EDP=0.25; xm_Cost=29337.;
F_DS_i1 = normcdf((log(x_Accel_pdf/xm_EDP))/beta_EDP); % compute fragility function using Eq. 1 and estimated parameters
% xm_EDP=1.91; beta_EDP=0.30;
xm_EDP=1.31; beta_EDP=0.25;
F_DS_i2 = normcdf((log(x_Accel_pdf/xm_EDP))/beta_EDP); % compute fragility function using Eq. 1 and estimated parameters
PDS_ij_EDP = F_DS_i1 - F_DS_i2; % 1<=j<=n, i.e. some damage
elseif i_n == 3
% xm_EDP=1.91; beta_EDP=0.30; xm_Cost=55200.;
xm_EDP=1.31; beta_EDP=0.25; xm_Cost=55200.;
F_DS_ij = normcdf((log(x_Accel_pdf/xm_EDP))/beta_EDP); % compute fragility function using Eq. 1 and estimated parameters
PDS_ij_EDP = F_DS_ij; % j=n, i.e. biggest damage
end
% % % % % % % % % % % % % % % %
elseif i_m == 2
% % % % % % % % % % % % % % % %
numCompPerStory = 21*(9.144*9.144) / (3.66*3.66); % total area of each floor divided by 3.66^2m^2 (each sprinkler protects 3.66*3.66m^2)
if i_n == 0 % Automatic sprinklers
xm_EDP=32.0; beta_EDP=1.40; xm_Cost=0.% 900.;
F_DS_ij = normcdf((log(x_Accel_pdf/xm_EDP))/beta_EDP); % compute fragility function using Eq. 1 and estimated parameters
PDS_ij_EDP = 1.0 - F_DS_ij; % j=0, i.e. no damage
elseif i_n == 1
xm_EDP=32.0; beta_EDP=1.40; xm_Cost=900.;
F_DS_ij = normcdf((log(x_Accel_pdf/xm_EDP))/beta_EDP); % compute fragility function using Eq. 1 and estimated parameters
PDS_ij_EDP = F_DS_ij; % j=n, i.e. biggest damage
end
% % % % % % % % % % % % % % % %
elseif i_m == 3
% % % % % % % % % % % % % % % %
numCompPerStory = 2.0; % total number of elevators (same as: Hwang and Lignos, 2017)
if i_n == 0 % Drywall finish
xm_EDP=0.50; beta_EDP=0.28; xm_Cost=0.% 868.;
F_DS_ij = normcdf((log(x_Accel_pdf/xm_EDP))/beta_EDP); % compute fragility function using Eq. 1 and estimated parameters
PDS_ij_EDP = 1.0 - F_DS_ij; % j=0, i.e. no damage
elseif i_n == 1
xm_EDP=0.50; beta_EDP=0.28; xm_Cost=868.;
F_DS_ij = normcdf((log(x_Accel_pdf/xm_EDP))/beta_EDP); % compute fragility function using Eq. 1 and estimated parameters
PDS_ij_EDP = F_DS_ij; % j=n, i.e. biggest damage
end
end
end