B.2. THE R-F METHOD FOR EXAMPLE 2.7 227
format short e
disp('reliability')
R=normcdf(beta0)
% Displaye iterative process and write it to Excel file
disp(ddp)
xlswrite('example2_6',ddp); % Write the iterative process
% into a Excel file
B.2 THE R-F METHOD FOR EXAMPLE 2.7
% The R-F method for Example 2.7
% The Limit State function: g(M,ka,se, kf, d)
clear; % Clean memory
% Input the distribution parameters dp1 and dp2
% Means or the first distribution parameter
dp1=[0.315,0.905,24.7,1.562,1.250];
% Standard deviation or the second distribution parameter
dp2=[0.142,0.0543,2.14,0.1250,0.00125];
kb=0.8507; % Size mofification factor
r=1; %Number of non-normal distribution
n=5; % Total five random variables
% Calculate the mean of lognorm M and the initial
% point x0(i)
x0(1)=exp(dp1(1)+dp2(1)^2/2);
for i=2:n-1
x0(i)=dp1(i);
end
% Use the surface of the limit state function
% to calculate x0(n)
x0(5)=(x0(4)*32*x0(1)/pi/(x0(2)*kb*x0(3)))^(1/3);
beta=0; %Set the beta =0
% Iterative process starting
for j=1:1000
%Calculate the equivalent mean and standard deviation
zteq=norminv(logncdf(x0(1),dp1(1),dp2(1)));
steq=normpdf(zteq)/lognpdf(x0(1),dp1(1),dp2(1));
mteq=x0(1)-zteq*steq;
% Mean and standard deviation matrix
meq(1)=mteq;