用Opensees进行IDA分析桥墩模型命令流Word文档格式.docx

用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