How to find beta and Pf values using Monte Carlo

%USE MONTE CARLO METHOD TO CALCULATE PROBABILITY OF FAILURE OF WIND TURBINE
%The given limit state function for the given is 
%g = -3.689 + 0.0184*D^2 + 0.2938*V ;
%Given :
D = 90.0219343 ;
V = 6.672408424 ;
mean_D = 90 ;
mean_V = 6.6724 ;
sigma_D = 10.29023 ;
sigma_V = 0.67708 ;
%Now use normrnd function
n = 10e6;
D = normrnd(mean_D, sigma_D, [n,1]);
V = normrnd(mean_V, sigma_V, [n,1]);
%
%Create for loop
K = nnz(-3.689 + 0.0184*D.^2 + 0.2938*V < 0);
P_fail = K/n 

Answer 1

[Aditya]

Hi Ashwini,

From the code that you have shared I can see that you have already calculated the Pf (probability of failure) value. To calculate the beta value, you can utilize the "norminv" function in MATLAB. 

For more detailed information on how to use the "norminv" function, please refer to the documentation available at: https://in.mathworks.com/help/stats/norminv.html

Hope this helps!