用Opensees进行IDA分析(桥墩模型命令流)

wipe#Openseesdandun##Units:kN,m,sec#-----------------#Startofmodelgeneration#-----------------#CreateModeBulider(withtwo-dimensionsand3DOF/node)modelbasic-ndm2-ndf3#-----------------#tagXYnode10.00.0node20.00.0node30.02.0node40.04.0node50.06.0node60.08.0node70.010.0node80.012.0node90.014.0node100.016.0node110.018.0node120.020.0#-----------------#Fixsupportsatbaseofcolumns#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#---------------------#setsomeparameterssetcolWidth8.18setcolDepth4.28setcover0.05setAs0.00049#somevariablesderivedfromtheparameterssety1[expr$colDepth/2.0]setz1[expr$colWidth/2.0]sectionFiber1{#Createtheconcretecorefiberspatchrect12030[expr$cover-$y1][expr$cover-$z1][expr$y1-$cover][expr$z1-$cover]#Createtheconcretecoverfibers(top,bottom,left,right)patchrect2205[expr-$y1][expr$z1-$cover]$y1$z1patchrect2205[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#taggeomTransfLinear1#Numberofintegrationpointsalonglengthofelementssetnp5#CreatethecolumnsusingBeam-columnelements#tagndIndJsecIDtransfTagelementnonlinearBeamColumn223$np11elementnonlinearBeamColumn334$np11elementnonlinearBeamColumn445$np11elementnonlinearBeamColumn556$np11elementnonlinearBeamColumn667$np11elementnonlinearBeamColumn778$np11elementnonlinearBeamColumn889$np11elementnonlinearBeamColumn9910$np11elementnonlinearBeamColumn101011$np11elementnonlinearBeamColumn111112$np11equalDOF1212elementzeroLength11112-mat13-dir3setm[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#CreateaPlainloadpatternwithaLinearTimeSeriespatternPlain1Constant{#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#-------------------analysisStaticinitializeanalyze1#Endofstaticanalysis#-------------------#--建立文件#第1条波---15setdataDirDynamic-Output-15;#nameofoutputfolderfilemkdir$dataDir;#createoutputfolder#------#定义峰值加速度（g）foreachpga{0.050.10.150.20.250.30.350.40.450.5}{puts[formatThePGAis%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$dataDir/new%4.3f.out$pga]Node_Disp6500710390-columns12recorderElement-file[format$dataDir/Curvature-%4.3f.out$pga]-time-ele2section1deformationrecorderElement-file[format$dataDir/Force-%4.3f.out$pga]-time-ele2section1forcewipeAnalysissetTime0.0setg[expr6.289*9.81*$pga]settabasFNPath-filePath15.txt-dt0.01-factor$g#Definetheexcitationusingthe960groundmotionrecords#tagdiraccelseriesargspatternUniformExcitation21-accel$tabasFN#addsomemassproportionaldamping#defineDAMPING#---------------------------------------------------------------------setxDamp0.05;#2%dampingratiosetlambda[eigen1];#eigenvaluemode1setomega[exprpow($lambda,0.5)];setalphaM0.;#M-prop.damping;D=alphaM*MsetbetaKcurr0.;#K-proportionaldamping;+beatKcurr*KCurrentsetbetaKcomm[expr2.*$xDamp/($omega)];#K-prop.dampingparameter;+betaKcomm*KlastCommittsetbetaKinit0.;#initial-stiffnessproportionaldamping+beatKinit*Kini#definedampingrayleigh$alphaM$betaKcurr$betaKinit$betaKcomm;#RAYLEIGHdamping#---------------------------------------------------------------------#CreatetheconvergencetesttestEnergyIncr1.0e-8305#CreatethesolutionalgorithmalgorithmNewton#CreatethesystemofequationstorageandsolversystemSparseGeneral-piv#Createtheconstrainthand