-
Notifications
You must be signed in to change notification settings - Fork 1
/
info_Comp_Fragility_NonStructural_Drift.m
98 lines (63 loc) · 4.33 KB
/
info_Comp_Fragility_NonStructural_Drift.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
98
function [PDS_ij_EDP, xm_Cost, numCompPerStory] = info_Comp_Fragility_NonStructural_Drift(i_n, i_m, x_PSDR_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) / 6.0; % total area of each floor divided by 6m^2 (Hwang and Lignos, 2017)
if i_n == 0 % Drywall partition
xm_EDP=0.0039; beta_EDP=0.17; xm_Cost=0.% 90.;
F_DS_ij = normcdf((log(x_PSDR_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.0039; beta_EDP=0.17; xm_Cost=90.;
F_DS_i1 = normcdf((log(x_PSDR_pdf/xm_EDP))/beta_EDP); % compute fragility function using Eq. 1 and estimated parameters
xm_EDP=0.0085; beta_EDP=0.23;
F_DS_i2 = normcdf((log(x_PSDR_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=0.0085; beta_EDP=0.23; xm_Cost=530.;
F_DS_ij = normcdf((log(x_PSDR_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) / 6.0; % total area of each floor divided by 6m^2 (Hwang and Lignos, 2017)
if i_n == 0 % Drywall finish
xm_EDP=0.0039; beta_EDP=0.17; xm_Cost=0.% 90.;
F_DS_ij = normcdf((log(x_PSDR_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.0039; beta_EDP=0.17; xm_Cost=90.;
F_DS_i1 = normcdf((log(x_PSDR_pdf/xm_EDP))/beta_EDP); % compute fragility function using Eq. 1 and estimated parameters
xm_EDP=0.0085; beta_EDP=0.23;
F_DS_i2 = normcdf((log(x_PSDR_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=0.0085; beta_EDP=0.23; xm_Cost=250.;
F_DS_ij = normcdf((log(x_PSDR_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 = (9.144*3.0*4.0 + 9.144*8.0) * 3.6576 * 6.0 / 2.81278 /6.0; % For each story.. See SMRF paper by Hwang
if i_n == 0 % Exterior glazing
xm_EDP=0.04; beta_EDP=0.36; xm_Cost=0.% 440.;
F_DS_ij = normcdf((log(x_PSDR_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.04; beta_EDP=0.36; xm_Cost=440.;
F_DS_i1 = normcdf((log(x_PSDR_pdf/xm_EDP))/beta_EDP); % compute fragility function using Eq. 1 and estimated parameters
xm_EDP=0.046; beta_EDP=0.33;
F_DS_i2 = normcdf((log(x_PSDR_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=0.046; beta_EDP=0.33; xm_Cost=440.;
F_DS_ij = normcdf((log(x_PSDR_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