matlab 模特卡罗计算节点可靠度

利用matlab编写装配式节点蒙特卡洛可靠度计算 包括opensees程序

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.24;‘）;
f

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-03-05 11:12  节点可靠度\
     文件       22033  2018-02-26 20:27  节点可靠度\csj.m
     文件       15194  2018-02-26 20:36  节点可靠度\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