用Opensees进行IDA分析桥墩模型命令流Word文档格式.docx
《用Opensees进行IDA分析桥墩模型命令流Word文档格式.docx》由会员分享,可在线阅读,更多相关《用Opensees进行IDA分析桥墩模型命令流Word文档格式.docx(17页珍藏版)》请在冰豆网上搜索。
#tagDXDYRZfix1111
#Concretetagfcec0fcuecu
#Coreconcrete(confined)
uniaxialMaterialConcrete011-25600.0-0.00219-17780.0-0.01#Coverconcrete(unconfined)
uniaxialMaterialConcrete012-23400.0-0.002-0.0-0.006
#STEEL
#Reinforcingsteelsetfy400000.0;
#YieldstresssetE200000000.0;
#Young'
smodulus
#tagfyE0buniaxialMaterialSteel023$fy$E0.0118.50.9250.15uniaxialMaterialElastic1129043600uniaxialMaterialElastic1212326600uniaxialMaterialElastic13587247596#Definecross-sectionfornonlinearcolumns
#setsomeparameters
setcolWidth8.18
setcolDepth4.28
setcover0.05
setAs0.00049
#somevariablesderivedfromtheparameters
sety1[expr$colDepth/2.0]
setz1[expr$colWidth/2.0]
sectionFiber1{
#Createtheconcretecorefibers
patchrect12030[expr$cover-$y1][expr$cover-$z1][expr$y1-$cover][expr$z1-$cover]
#Createtheconcretecoverfibers(top,bottom,left,right)
patchrect2205[expr-$y1][expr$z1-$cover]$y1$z1
patchrect2205[expr-$y1][expr-$z1]$y1[expr$cover-$z1]
patchrect2510[expr-$y1][expr$cover-$z1][expr$cover-$y1][expr$z1-$cover]
patchrect2510[expr$y1-$cover][expr$cover-$z1]$y1[expr$z1-$cover]
#Createthereinforcingfibers(left,middle,right)
layerstraight3175$As[expr$y1-$cover][expr$z1-$cover][expr$y1-$cover][expr$cover-$z1]layerstraight3175$As[expr$cover-$y1][expr$z1-$cover][expr$cover-$y1][expr$cover-$z1]layerstraight3115$As[expr$y1-$cover][expr$z1-$cover][expr$cover-$y1][expr$z1-$cover]layerstraight3115$As[expr$y1-$cover][expr$cover-$z1][expr$cover-$y1][expr$cover-$z1]}
#Definecolumnelements
#Geometryofcolumnelements
#tag
geomTransfLinear1
#Numberofintegrationpointsalonglengthofelements
setnp5
#CreatethecolumnsusingBeam-columnelements
#tagndIndJsecIDtransfTag
elementnonlinearBeamColumn223$np11
elementnonlinearBeamColumn334$np11
elementnonlinearBeamColumn445$np11
elementnonlinearBeamColumn556$np11
elementnonlinearBeamColumn667$np11
elementnonlinearBeamColumn778$np11
elementnonlinearBeamColumn889$np11
elementnonlinearBeamColumn9910$np11
elementnonlinearBeamColumn101011$np11
elementnonlinearBeamColumn111112$np11
equalDOF1212
elementzeroLength11112-mat13-dir3
setm[expr3355.7]
setn[expr175.1]#tagMXMYRZmass2$n$n1e-8mass3$n$n1e-8mass4$n$n1e-8mass5$n$n1e-8mass6$n$n1e-8mass7$n$n1e-8mass8$n$n1e-8mass9$n$n1e-8mass10$n$n1e-8mass11$n$n1e-8mass12$m$m1e-8
#SetaparameterfortheaxialloadsetP33557.0;
#ofaxialcapacityofcolumnssetQ1715.5;
#ofaxialcapacityofcolumns
#CreateaPlainloadpatternwithaLinearTimeSeriespatternPlain1"
Constant"
{
#Createnodalloadsatnodes9
#ndFXFYMZload120.0[expr-$P]0.0load20.0[expr-$Q]0.0load30.0[expr-$Q]0.0load40.0[expr-$Q]0.0load50.0[expr-$Q]0.0load60.0[expr-$Q]0.0load70.0[expr-$Q]0.0load80.0[expr-$Q]0.0load90.0[expr-$Q]0.0load100.0[expr-$Q]0.0load110.0[expr-$Q]0.0
}systemSparseGeneral-piv#CreatetheconstrainthandlerconstraintsTransformation
#Createthetimeintegrationscheme
#CreatetheDOFnumberernumbererRCMtestNormDispIncr1.0e-8305
#CreatethesolutionalgorithmalgorithmNewton
#createthetransientanalysisintegratorLoadControl1
#Endofanalysisgeneration
#analysisStatic
initializeanalyze1
#Endofstaticanalysis
#--建立文件
#第1条波---15
setdataDirDynamic-Output-15;
#nameofoutputfolder
filemkdir$dataDir;
#createoutputfolder
#定义峰值加速度(g)
foreachpga{0.050.10.150.20.250.30.350.40.450.5}{
puts[format"
ThePGAis%3.2fg"
$pga]
recorderNode-file[format"
$dataDir/zhuanjiao%4.3f.out"
$pga]-time-node2-dof3disprecorderNode-file[format"
$dataDir/new%4.3f.out"
$pga]-time-node12-dof1disprecorderplot[format"
$pga]Node_Disp6500710390-columns12recorderElement-file[format"
$dataDir/Curvature-%4.3f.out"
$pga]-time-ele2section1deformationrecorderElement-file[format"
$dataDir/Force-%4.3f.out"
$pga]-time-ele2section1forcewipeAnalysis
setTime0.0setg[expr6.289*9.81*$pga]
settabasFN"
Path-filePath15.txt-dt0.01-factor$g"
#Definetheexcitationusingthe960groun