• 大小: 6.85MB
    文件类型: .zip
    金币: 1
    下载: 0 次
    发布日期: 2023-11-01
  • 语言: 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

评论

共有 条评论