完整版用Opensees进行IDA分析桥墩模型命令流.docx

完整版用Opensees进行IDA分析桥墩模型命令流.docx

《完整版用Opensees进行IDA分析桥墩模型命令流.docx》由会员分享，可在线阅读，更多相关《完整版用Opensees进行IDA分析桥墩模型命令流.docx（26页珍藏版）》请在冰豆网上搜索。

完整版用Opensees进行IDA分析桥墩模型命令流

wipe

#Openseesdandun

#

#Units:

kN,m,sec

#-----------------

#Startofmodelgeneration

#-----------------

#CreateModeBulider（withtwo-dimensionsand3DOF/node）

modelbasic-ndm2-ndf3

#-----------------

#tagXY

node10.00.0

node20.00.0

node30.02.0

node40.04.0

node50.06.0

node60.08.0

node70.010.0

node80.012.0

node90.014.0

node100.016.0

node110.018.0

node120.020.0

#-----------------

#Fixsupportsatbaseofcolumns

#tagDXDYRZ

fix1111

#----------------

#Concretetagfcec0fcuecu

#Coreconcrete（confined）

uniaxialMaterialConcrete011-25600.0-0.00219-17780.0-0.01

#Coverconcrete（unconfined）

uniaxialMaterialConcrete012-23400.0-0.002-0.0-0.006

#STEEL

#Reinforcingsteel

setfy400000.0;#Yieldstress

setE200000000.0;#Young'smodulus

#tagfyE0b

uniaxialMaterialSteel023$fy$E0.0118.50.9250.15

uniaxialMaterialElastic1129043600

uniaxialMaterialElastic1212326600

uniaxialMaterialElastic13587247596

#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]

#tagMXMYRZ

mass2$n$n1e-8

mass3$n$n1e-8

mass4$n$n1e-8

mass5$n$n1e-8

mass6$n$n1e-8

mass7$n$n1e-8

mass8$n$n1e-8

mass9$n$n1e-8

mass10$n$n1e-8

mass11$n$n1e-8

mass12$m$m1e-8

#Setaparameterfortheaxialload

setP33557.0;#ofaxialcapacityofcolumns

setQ1715.5;#ofaxialcapacityofcolumns

#CreateaPlainloadpatternwithaLinearTimeSeries

patternPlain1"Constant"{

#Createnodalloadsatnodes9

#ndFXFYMZ

load120.0[expr-$P]0.0

load20.0[expr-$Q]0.0

load30.0[expr-$Q]0.0

load40.0[expr-$Q]0.0

load50.0[expr-$Q]0.0

load60.0[expr-$Q]0.0

load70.0[expr-$Q]0.0

load80.0[expr-$Q]0.0

load90.0[expr-$Q]0.0

load100.0[expr-$Q]0.0

load110.0[expr-$Q]0.0

}

systemSparseGeneral-piv

#Createtheconstrainthandler

constraintsTransformation

#Createthetimeintegrationscheme

#CreatetheDOFnumberer

numbererRCM

testNormDispIncr1.0e-8305

#Createthesolutionalgorithm

algorithmNewton

#createthetransientanalysis

integratorLoadControl1

#-----------------

#Endofanalysisgeneration

#-------------------

analysisStatic

initialize

analyze1

#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-dof3disp

recorderNode-file[format"$dataDir/new%4.3f.out"$pga]-time-node12-dof1disp

recorderplot[format"$dataDir/new%4.3f.out"$pga]Node_Disp6500710390-columns12

recorderElement-file[format"$dataDir/Curvature-%4.3f.out"$pga]-time-ele2section1deformation

recorderElement-file[format"$dataDir/Force-%4.3f.out"$pga]-time-ele2section1force

wipeAnalysis

setTime0.0

setg[expr6.289*9.81*$pga]

settabasFN"Path-filePath15.txt-dt0.01-factor$g"

#Definetheexcitationusingthe960groundmotionrecords

#tagdiraccelseriesargs

patternUniformExcitation21-accel$tabasFN

#addsomemassproportionaldamping

#defineDAMPING

#---------------------------------------------------------------------

setxDamp0.05;#2%dampingratio

setlambda[eigen1];#eigenvaluemode1

setomega[exprpow（$lambda,0.5）];

setalphaM0.;#M-prop.damping;D=alphaM*M

setbetaKcurr0.;#K-proportionaldamping;+beatKcurr*KCurrent

setbetaKcomm[expr2.*$xDamp/（$omega）];#K-prop.dampingparameter;+betaKcomm*KlastCommitt

setbetaKinit0.;#initial-stiffnessproportionaldamping+beatKinit*Kini

#definedamping

rayleigh$alphaM$betaKcurr$betaKinit$betaKcomm;#RAYLEIGHdamping

#---------------------------------------------------------------------

#Createtheconvergencetest

testEnergyIncr1.0e-8305

#Createthesolutionalgorithm

algorithmNewton

#Createthesystemofequationstorageandsolver

systemSparseGeneral-piv

#Createtheconstrainthandler

constraintsTransformation

#Createthetimeintegrationscheme

integratorNewmark0.50.25

#CreatetheDOFnumberer

numbererRCM

#createthetransientanalysis

analysisVariableTransient

recorderElement-file[format"$dataDir/ele1Concrete-%4.3f.out"$pga]-time-ele2section1fiber-2.05-4.01stressStrain

recorderElement-file[format"$dataDir/ele1sec1StressStraingangjin%4.3f.out"$pga]-time-ele2section1fiber-2.05-4.03stressStrain

recorderplot[format"$dataDir/ele1sec1StressStraingangjin%4.3f.out"$pga]strain-stress0300700400-columns32

recorderElement-file[format"$dataDir/ele1-%4.3f.out"$pga]-time-ele1localForce

recorderElement-file[format"$dataDir/ele2-%4.3f.out"$pga]-time-ele2localForce

#Performtheanalysis

#numStepsdt

analyze40000.0050.00000000010.00530

#------

setPI3.1415926

setlambda[eigen1]

setomega[exprpow（$lambda,0.5）]

setTperiod[expr2*$PI/$omega]

puts"T1:

$Tperiod"

removeloadPattern2

removerecorders;

reset;

}

注：

同文目录下需要编号为15.txt的地震记录文件，时间间隔为0.01s，PGA为0.159g，如下：

-0.0124

-0.00268

0.0082

0.0198

0.0308

0.0377

0.0378

0.0309

0.0201

0.00885

0.00074

-0.00136

0.0027

0.0105

0.0189

0.025

0.0274

0.0267

0.026

0.0282

0.0345

0.0436

0.0532

0.0595

0.0594

0.0521

0.0405

0.0278

0.0164

0.00667

-0.00233

-0.0117

-0.0221

-0.0331

-0.0444

-0.0552

-0.0628

-0.064

-0.0578

-0.0464

-0.0334

-0.0213

-0.0109

-0.00136

0.00811

0.0169

0.0234

0.0268

0.0275

0.0265

0.0251

0.023

0.0186

0.0106

0.000209

-0.00818

-0.0109

-0.00759

-0.000752

0.00608

0.00956

0.00827

0.0033

-0.00281

-0.00686

-0.00659

-0.00208

0.00454

0.0101

0.0121

0.0107

0.00815

0.00744

0.00959

0.0131

0.0152

0.0145

0.0113

0.00779

0.00629

0.00706

0.00692

0.00174

-0.00993

-0.0262

-0.0439

-0.0605

-0.0741

-0.0814

-0.0803

-0.0724

-0.061

-0.0496

-0.0398

-0.0313

-0.023

-0.014

-0.004

0.00632

0.0141

0.0152

0.00802

-0.00407

-0.0126

-0.0104

0.00346

0.024

0.0453

0.0637

0.077

0.0804

0.0699

0.0468

0.0167

-0.0125

-0.0303

-0.029

-0.0101

0.0178

0.0438

0.0612

0.0698

0.0742

0.0769

0.0757

0.0676

0.0525

0.033

0.0124

-0.00684

-0.0214

-0.0272

-0.0233

-0.0125

0.000534

0.0125

0.0224

0.0312

0.0397

0.0458

0.046

0.0391

0.0271

0.0134

0.00139

-0.00477

-0.00219

0.00819

0.0226

0.0361

0.0423

0.0376

0.0236

0.00522

-0.0128

-0.0283

-0.0418

-0.0549

-0.0678

-0.0777

-0.0823

-0.0818

-0.0784

-0.0748

-0.073

-0.0736

-0.0756

-0.0774

-0.0765

-0.071

-0.0611

-0.0485

-0.0357

-0.0238

-0.0131

-0.0026

0.00822

0.0196

0.0299

0.034

0.0279

0.0124

-0.00687

-0.0218

-0.0274

-0.0247

-0.0182

-0.0124

-0.00987

-0.0116

-0.0177

-0.027

-0.0378

-0.0483

-0.0582

-0.0672

-0.0724

-0.0688

-0.0548

-0.0331

-0.00849

0.015

0.0362

0.056

0.0759

0.097

0.119

0.137

0.145

0.141

0.128

0.113

0.0983

0.0875

0.0817

0.0814

0.0858

0.0923

0.0985

0.1

0.0938

0.0797

0.061

0.0417

0.024

0.00794

-0.00748

-0.0235

-0.0406

-0.0584

-0.0764

-0.0935

-0.106

-0.11

-0.104

-0.0925

-0.0788

-0.0664

-0.0563

-0.0482

-0.0441

-0.0465

-0.0553

-0.068

-0.0814

-0.0915

-0.0939

-0.0873

-0.0743

-0.059

-0.0449

-0.0328

-0.0219

-0.0107

0.00151

0.0148

0.0284

0.0419

0.0548

0.0673

0.0787

0.0869

0.0906

0.0905

0.0887

0.0865

0.0828

0.0746

0.0606

0.0421

0.0218

0.0021

-0.0163

-0.0335

-0.048

-0.0566

-0.0575

-0.0524

-0.0447

-0.0375

-0.032

-0.0275

-0.0231

-0.018

-0.0119

-0.0056

-0.000207

0.00316

0.00446

0.0051

0.0077

0.014

0.0236

0.0348

0.0447

0.0491

0.0452

0.0344

0.0208

0.0105