非序贯蒙特卡洛算法计算节点可靠度matlab程序

非序贯蒙特卡洛模拟法 通常被称为状态抽样法，它被广泛用在电力系统风险评估中。这个方法的依据是，一个系统状态是所有元件状态的组合，且每一元件状态可由对元件出现在该状态的概率进行抽样来确定。

clear
clc
mu=[29.1;33.33;565;200000;100];
sigma=[3.783;4.0;39.55;4000;30];%均值标准差向量
sLn=sqrt（log（1+（sigma（1:4）./mu（1:4））.^2））;%对数分布的正态化还原
mLn=log（mu（1:4））-sLn.^2/2;
aEv=sqrt（6）*sigma（5）/pi;
uEv=-psi（1）*aEv-mu（5）;
level1=0;
level2=0;
level3=0;
level4=0;

for i=1:5
CUnconfFc=-lognrnd（mLn（1）sLn（1））;
BUnconfFc=CUnconfFc;
PCUnconfFc=-lognrnd（mLn（2）sLn（2））;
PBUnconfFc=PCUnconfFc;
CFy=lognrnd（mLn（3）sLn（3））;
BFy=CFy;
CEs=lognrnd（mLn（4）sLn（4））;
BEs=CEs;
D=-evrnd（uEvaEv）;

fid=fopen（‘csj.tcl‘‘wt‘）;
fprintf（fid‘%s\n‘‘model BasicBuilder -ndm 2 -ndf 3‘）;
fprintf（fid‘%s‘‘set fName “file‘）;
fprintf（fid‘%d‘i）;
fprintf（fid‘%s\n‘‘“;‘）;
fprintf（fid‘%s\n‘‘source ProcMKPC2.tcl‘）;
fprintf（fid‘%s\n‘‘puts “a“‘）;
fprintf（fid‘%s\n‘‘source ProcRC2.tcl‘）;
fprintf（fid‘%s\n‘‘puts “a“‘）;
fprintf（fid‘%s\n‘‘set Strfactor 145; set Lenfactor [expr 1/25.4];‘）;
fprintf（fid‘%s\n‘‘set colY 600; set colZ 600;‘）;
fprintf（fid‘%s\n‘‘set bmY 950; set bmZ 600;‘）;
fprintf（fid‘%s\n‘‘set colCov 35; set bmCov1 30; set bmCov2 30; set bmCov $bmCov1;‘）;
fprintf（fid‘%s\n‘‘set JointWidth [expr $colY]; set JointHeight [expr $bmY]; set JointDepth $colZ ;‘）;
fprintf（fid‘%s\n‘‘set BeamLength1 500; set BeamLength2 1850; set BeamLength3 700; set BeamLength4 800; set ColumnLength1 300; set ColumnLength2 700; set ColumnLength3 600;‘）;
fprintf（fid‘%s\n‘‘set JointVolume [expr $JointWidth*$JointHeight*$JointDepth];‘）;

fprintf（fid‘%s‘‘set CUnconfFc ‘）;fprintf（fid‘%d\n‘CUnconfFc）; 

fprintf（fid‘%s\n‘‘set CUnconfEc -0.002;‘）;
fprintf（fid‘%s\n‘‘set CTSmispace 50; set CTSspace 200; set CTSlength 2120; set CTSFy 568.7; set CTSarea 50.24;‘）;
fprintf（fid‘%s‘‘set CFy ‘）;fprintf（fid‘%d\n‘CFy）; %
fprintf（fid‘%s‘‘set CEs ‘）;fprintf（fid‘%d\n‘CEs）;
fprintf（fid‘%s\n‘‘set CsHratio 0.004216; set CAs1 615.4; set CAs2 490.6‘）;

fprintf（fid‘%s\n‘‘procMKPC $CUnconfFc $CUnconfEc $colY $colZ $colCov $CTSmispace $CTSlength $CTSFy $CTSarea $Strfactor $Lenfactor‘）;
fprintf（fid‘%s\n‘‘set CUnconfFcu [lindex $concreteProp 2]; set CUnconfEcu [lindex $concreteProp 3];‘）;
fprintf（fid‘%s\n‘‘set CConfmiFc [lindex $concreteProp 4]; set CConfmiEc [lindex $concreteProp 5];‘）;
fprintf（fid‘%s\n‘‘set CConfmiFcu [lindex $concreteProp 6]; set CConfmiEcu [lindex $concreteProp 7];‘）;
fprintf（fid‘%s\n‘‘procMKPC $CUnconfFc $CUnconfEc $colY $colZ $colCov $CTSspace $CTSlength $CTSFy $CTSarea $Strfactor $Lenfactor‘）;
fprintf（fid‘%s\n‘‘set CUnconfFcu [lindex $concreteProp 2]; set CUnconfEcu [lindex $concreteProp 3];‘）;
fprintf（fid‘%s\n‘‘set CConfFc [lindex $concreteProp 4]; set CConfEc [lindex $concreteProp 5];‘）;
fprintf（fid‘%s\n‘‘set CConfFcu [lindex $concreteProp 6]; set CConfEcu [lindex $concreteProp 7];‘）;
fprintf（fid‘%s‘‘set BUnconfFc ‘）;fprintf（fid‘%d\n‘BUnconfFc）;
fprintf（fid‘%s\n‘‘set BUnconfEc -0.002;‘）;



fprintf（fid‘%s\n‘‘set BTSmispace 50; set BTSspace 100; set BTSlength 2860; set BTSFy 568.7; set BTSarea 50

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件       22043  2018-11-12 16:48  节点可靠度\csj.m
     文件       15194  2018-11-12 16:48  节点可靠度\csj.tcl
     文件    20073472  2018-02-26 19:51  节点可靠度\OpenSees.exe
     文件       16005  2018-02-26 19:52  节点可靠度\PR.tcl
     文件        2507  2018-02-26 19:52  节点可靠度\ProcMKPC2.tcl
     文件        4055  2018-02-26 19:52  节点可靠度\ProcRC2.tcl
     文件        1968  2018-02-26 19:52  节点可靠度\ProcUniaxialPinching.tcl