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