地震分析算例.docx

地震分析算例.docx

《地震分析算例.docx》由会员分享，可在线阅读，更多相关《地震分析算例.docx（44页珍藏版）》请在冰豆网上搜索。

地震分析算例

地震分析算例（ANSYS）

土木工程中除了常见的静力分析以外，动力分析，特别是结构在地震荷载作用下的受力分析，也是土木工程中经常遇到的问题。

结构的地震分析根据现行抗震规范要求，一般分为以下两类：

基于结构自振特性的地震反应谱分析和基于特定地震波的地震时程分析。

本算例将以一个4质点的弹簧－质点体系来说明如何使用有限元软件进行地震分析。

更复杂结构的分析其基本过程也与之类似。

关键知识点：

（a）模态分析

（b）谱分析

（c）地震反应谱输入

（d）地震时程输入

（e）时程动力分析

（1）在ANSYS窗口顶部静态菜单，进入Parameters菜单，选择ScalarParameters选项，在输入窗口中填入DAMPRATIO=0.02，即所有振型的阻尼比为2%

（2）ANSYS主菜单Preprocessor->Elementtype->Add/Edit/Delete，添加Beam188单元

（3）在ElementTypes窗口中，选择Beam188单元，选择Options，进入Beam188的选项窗口，将第7个和第8个选项，Stress/Strain（SectPoints）K7,Stress/Strain（SectNods）K8，从None改为MaxandMinOnly。

即要求Beam188单元输出积分点和节点上的最大、最小应力和应变

（4）在ElementTypes窗口中，继续添加Mass21集中质量单元

（5）下面输入材料参数，进入ANSYS主菜单Preprocessor->MaterialProps->MaterialModels菜单，在MaterialModelNumber1中添加Structural->Linear->Elastic->Isotropic属性，输入材料的弹性模量EX和泊松比PRXY分别为210E9和0.3。

（6）继续给MaterialModelNumber1添加Density属性，输入密度为7800。

（7）继续给MaterialModelNumber1添加Damping属性，采用参数化建模，输入阻尼类型为Constant，数值为DAMPRATIO

（8）接着建立梁单元的几何属性，和上一个例子一样，采用Sections建模，进入ANSYS主菜单Preprocessor->Sections->Beam->CommonSections，选择Sub-Type为工字型，截面尺寸W1=0.2，W2=0.2,W3=0.5,t1=0.01,t2=0.01,t3=0.008

（9）通过实参数输入集中质量单元的质量和转动惯量，在ANSYS主菜单中选择Preprocessor->RealConstants->Add/Edit/Delete菜单，在RealConstants窗口中选择Add，在ElementtypeforRealConstants选择Mass21，在RealConstantforNumber1窗口中输入1.6E2,1.6E2，如图所示。

即该质量单元在X和Y方向的质量都为160，由于本例子模型为平面问题，所以不必考虑Z方向的质量，同样也不考虑单元的转动惯量。

（10）继续添加第二类集中质量，过程和上面一样，但是输入的质量数值为1.2E2,1.2E2

（11）完成以上工作就完成了模型的基本数据准备，下面开始建立物理模型。

（12）在ANSYS主菜单中选择Preprocessor->Modeling->Create->Keypoints->InActiveCS，依次输入关键点编号和坐标：

K,1,0,0,0

K,2,0,3,0

K,3,0,6,0

K,4,0,9,0

K,5,0,12,0

K,6,0,0,10000

需要说明的是，关键点6为后面建立梁单元所需的截面方向控制点，在上一个例子中已经做过介绍。

（13）完成关键点输入后下面建立直线模型。

在ANSYS主菜单中选择Preprocessor->Modeling->Create->Lines->Lines->StraightLine，依次连接关键点1～5。

（14）下面给建立完的几何模型赋予材料属性，在ANSYS菜单中选择Meshing->MeshAttributes->PickedLines，选中所有的直线，进入LineAttributes窗口，选择相关选项材料属性，实参数，单元类型和截面类型都为1，并选择勾上参考点方向选项，然后点击OK后输入关键点6作为截面方向控制点。

（15）在ANSYS主菜单中选择Preprocessor->Meshing->SizeCntrls->ManualSize->Lines->PickedLines，在ElementSizesonPickedLines窗口中设定NDIVNo.ofelementdivisions为3，即将每条直线分为3段

（16）在ANSYS主菜单中选择Preprocessor->Meshing->Mesh->Lines，选择所有的直线，完成直线的网格划分。

（17）为了便于后面操作，将网格划分后的单元和节点编号进行适当的清理。

在ANSYS主菜单中选择Preprocessor->NumberingCtrls->MergeItems，在MergeCoincidentorEquivalentlyDefinedItems窗口中选择All，清理所有重复的元素。

同样选择NumberingCtrls->CompressNumber菜单，在CompressNumber中选择All，对节点和单元进行重新编号。

（18）下面建立集中质量单元，采取直接输入单元的方法建立。

在ANSYS主菜单中选择Preprocessor->Modeling->Create->Elements->ElemAttributes，在ElementAttributes中设定单元类型编号为2MASS21，材料编号任意，实参数编号为1。

（19）在ANSYS主菜单中选择Preprocessor->Modeling->Create->Elements->AutoNumbered->ThruNodes，选择节点2，建立第一个集中质量。

（20）再次进入第18步ElementAttributes窗口，设定实参数（Realconstantsetnumber）为2。

（21）重复19步，选择节点8，14，20，建立其他的三个集中质量单元。

（22）到此完成所有建模工作，下面开始进行结构分析

（23）进入ANSYS主菜单中Solution功能模块，选择Solution->DefineLoads->Apply->Structural->Displacement->OnNodes，选择节点1，设定约束所有的自由度。

（24）首先做一次静力分析，选择ANSYS主菜单Solution->AnalysisType->NewAnalysis，设定分析类型为Static

（25）选择ANSYS主菜单Solution->Solve->CurrentLS选项，进行一次静力分析

（26）分析完后，下面进行模态分析，在ANSYS主菜单中选择Solution->AnalysisType->NewAnalysis，选择分析类型为Model

（27）在ANSYS主菜单中选择Solution->AnalysisType->AnalysisOptions，输入模态分析方法为子空间法（Subspace），求解8阶模态，同时需要作模态扩展，扩展的模态为8阶，并计算单元应力和应变即改为yes，输入窗口如图。

（28）再次选择ANSYS主菜单Solution->Solve->CurrentLS，计算当前问题

（29）这时，如果需要看结果，可以进入后处理模块，即ANSYS主菜单GeneralPostproc，可以看到计算的各阶频率和振型。

!

存储各模态频率

*DIM,FRE,,NMODE

*DO,I,1,NMODE

*GET,FRE（I）,MODE,I,FREQ!

OBTAINMODEFREQENCYFORMODEI

*ENDDO

FINISH

ALLSEL

（30）完成结构自振分析后，下面就可以进行反应谱分析

（31）进入ANSYS主菜单Solution->AnalysisType->NewAnalysis，选择分析类型为Spectrum。

（32）在ANSYS主菜单Solution->AnalysisType->AnalysisOptions中，选择谱分析的类型为单点输入（Sing-ptresp）

（33）在ANSYS主菜单中选择Preprocessor->Loads->LoadStepOpts->Time/Frequenc->Damping，输入所有的阻尼为DAMPRATIO，如图所示

（34）下面需要定义地震的反应谱。

我国规范给定的是基于加速度的反应谱。

在ANSYS主菜单中选择Solution->LoadStepOpts->Spectrum->SinglePoint->Settings，设定反应谱类型为地震加速度，放大系数为1，输入方向为X方向（1，0，0）。

（35）接下来开始输入地震反应谱。

注意：

在进行反应谱前，若按7度多遇地震，即地震影响系数0.08第一组III场地TG＝0.45，可以通过如下方式从国标给定的地震谱频率中的结果中算出相应的谱值，不过其频率需要提前定义数组。

过程如下：

利用前面模态分析得到的频率。

（因为这些频率的取值附近容易引起共振，破坏最大）

!

计算反应谱数据（依据规范GB50011-2001第5.1.5条）

!

地震影响系数

GRAV=9.81!

重力加速度G

TG=0.45!

特性周期

AMAX=0.08!

水平地震影响系数最大值

C=0.05!

阻尼比

!

*DIM,A,,NMODE

*DIM,T,,NMODE

*DO,I,1,NMODE

T（I）=1.0/FRE（I）

*ENDDO

R=0.9+（0.05-C）/（0.5+5.0*C）

P1=0.02+（0.05-C）/8

P2=1+（0.05-C）/（0.06+1.7*C）

*DO,I,1,NMODE

*IF,T（I）,GE,0.0,AND,T（I）,LT,0.1,THEN

A（I）=（0.45+（10.0*P2-4.5）*T（I））*AMAX*GRAV

*ELSEIF,T（I）,GE,0.1,AND,T（I）,LE,TG

A（I）=P2*AMAX*GRAV

*ELSEIF,T（I）,GT,TG,AND,T（I）,LE,5*TG

A（I）=（TG/T（I））**R*P2*AMAX*GRAV

*ELSE

A（I）=（P2*0.2**R-P1*（T（I）-5*TG））*AMAX*GRAV

*ENDIF

*ENDDO

这里输入的反应谱按7度多遇地震，取地震影响系数为0.08，第一组，III类场地，卓越周期Tg=0.45s。

选择ANSYS主菜单Solution->LoadStepOpts->Spectrum->SinglePoint->FreqTable，输入频率反应谱Freq1~Freq12为0.167,0.25,0.333,0.44444,0.5,0.667,1,1.25,1.667,2.222,10,100000,如图

（36）选择ANSYS主菜单Solution->LoadStepOpts->Spectrum->SinglePoint->SpectrValues，输入对应的反应谱数值依次如下：

（37）最后选择ANSYS主菜单Solution->LoadStepOpts->Spectrum->SinglePoint->ModeCombine，设定振型组合方式为SRSS法，如图所示

（38）选择ANSYS主菜单Solution->Solve->CurrentLS，计算反应谱结果

（39）进入ANSYS主菜单后处理模块GeneralPostproc，在ANSYS窗口顶部菜单选择File->Readinputfrom，选择文件后缀名为*.mcom的文件

（40）进入ANSYS主菜单GeneralPostproc->PlotResults->DeformedShapes，选择绘制变形后形状和结构形状，得到地震反应谱分析的结构变形如图

（41）最后我们来进行地震时程分析，进行地震时程分析以前，首先要有一个地震时程记录，本例子给定的地震时程记录总长20秒，记录点间隔0.02s，共有1001个记录点。

该地震记录存放在RECORD.TXT文件中。

（42）首先建立两个变量，在ANSYS窗口顶部菜单选择Parameters->Scalarparameters，在窗口中输入NT=1001，即总共1001个记录点，DT=0.02，即记录点间隔0.02s

（43）然后还需要定义一个数组来存放地震记录，选择ANSYS窗口顶部菜单Parameters->Arrayparameters->Define/Edit，在弹出的ArrayParameters窗口中点击Add按钮，在AddNewArrayParameter窗口中输入数组的名称为AC，Parametertype类型选择默认的Array，再I，J，K，的行数设为NT个，其余默认。

如图所示。

（44）下面从数据文件中读入地震时程记录。

选择ANSYS窗口顶部菜单Parameters->Arrayparameters->Readfromfile，选择类型，Array，点击OK在弹出窗口中依次输入将数据读入的数组名称AC，可以用Browse选择目标数据文件RECORD.TXT，Label中选择IJK项，其中ROWS行的数改为NT，其余保持默认，并在最后要给出数据FormatSurroundedBy（）读入的格式（F8.3）既可。

并且可以参考Fortran相应的文件输入输出要求。

（45）完成数据读入后就可以进行地震反应分析。

由于本次地震分析一共要进行1001步，超过了ANSYS默认的最大1000步限制，因此首先需要进行调整。

在ANSYS的命令输入窗口内，输入以下命令：

Finish

/CONFIG,NRES,20000

更改后可以进入，List-〉Status-〉configuration，中查看。

（46）进入ANSYS主菜单Solution->AnalysisType->NewAnalysis，指定分析类型为瞬态分析Transient，Solutionmethod可以选择Full。

（47）进入ANSYS主菜单Solution->AnalysisType->Sol'nControls，在SolutionControls窗口里面选择Transient页面，设定瑞雷阻尼的数值。

由于阻尼的机理十分复杂，因此不同问题可能各不相同，这里仅介绍一种阻尼的取法。

输入质量阻尼系数为2*DAMPRATIO*FREQ1*2*3.1415926，刚度阻尼系数为2*DAMPRATIO/（FREQ1*2*3.1415926），FREQ1为结构的第一阶自振频率。

（48）下面将用一组循环控制语句来定义作用在结构上的加速度并进行计算。

读者可以将下面这部分语句用任意编辑器（比如windows的记事本）输入好以后复制ANSYS的输入栏中。

*DO,I,1,1001!

对变量I循环1001

ACEL,AC（I）,0,0!

对结构施加X方向加速度

TIME,I*0.02!

计算时间步长为0.02s

OUTRES,ALL,ALL!

输出所有结果

SOLVE!

求解

*ENDDO!

循环结束

（49）下面介绍适用ANSYS的时程后处理器TimeHistPostpro来处理计算结果，在ANSYS主菜单上选择TimeHistPostpro，这时会弹出时程变量窗口TimeHistoryVariables。

点击工具栏上第一个绿色加号按钮，添加时程变量。

在AddTime-HistoryVariable窗口中选择要添加的变量为节点的X方向位移。

（50）选择第20号节点，回到时程变量窗口TimeHistoryVariables，点击工具栏上第三个按钮，即绘制出节点的位移时程曲线如图

一个四层弹簧－质点模型的地震ANSYS分析

!

ExampleofseismicanalysisinANSYS

!

包括频率分析，谱分析和时程分析

!

Withmodel&frequencyanalysis,spectrumanalysisandtransientanalysis

FINI

/CLEAR

/UNITS,SI

/PREP7

DAMPRATIO=0.02!

振型阻尼系数都为0.02

ET,1,BEAM188

KEYOPT,1,7,1

KEYOPT,1,8,1

ET,2,MASS21,,,4

MP,EX,1,210E9

MP,NUXY,1,0.3

MP,DENS,1,7800

MPDATA,DAMP,1,,DAMPRATIO

SECTYPE,3,BEAM,I

SECDATA,0.2,0.2,0.5,0.01,0.01,0.008

R,1,1.6E2,1.6E2

R,2,1.2E2,1.2E2

K,1,0,0,0

K,2,0,3,0

K,3,0,6,0

K,4,0,9,0

K,5,0,12,0

K,6,0,0,10000

L,1,2

L,2,3

L,3,4

L,4,5

LSEL,,,,1,4

LATT,1,0,1,,6,,3

ALLS

LESIZE,ALL,,,3,,,,,1

LMESH,ALL

NUMMRG,ALL

NUMCMP,ALL

TYPE,2

REAL,1

E,2

TYPE,2

REAL,2

E,8

E,14

E,20

*SET,NT,1001

*SET,DT,0.02

*DIM,AC,,NT

*VREAD,AC,RECORD,txt,,IJK,NT,1,1,,

（F8.3）

/SOLU

D,1,ALL,

ANTYPE,0

SOLVE

FINISH

/SOLU!

模态分析

ANTYPE,2

MODOPT,SUBSP,8

MXPAND,8,,,1

SOLVE

FINI

!

得到自振频率1

*GET,FREQ1,MODE,1,FREQ

/ESHAPE,1.0

/SOLU

ANTYPE,SPECTR!

谱分析

SPOPT,SPRS

MDAMP,1,DAMPRATIO,,,,,,

SED,1

SVTYPE,2,1

FREQ,0.167,0.25,0.333,0.44444,0.5,0.667,1,1.25,1.667

SV,,0.154350625,0.191590625,0.210210625,0.224175625,0.250716714,0.329514922,0.484352764,0.598723486,0.786897371

FREQ,2.222,10,100000

SV,,1.034212766,1.034212766,0.3528

!

按7度多遇地震，地震影响系数0.08第一组III场地TG＝0.45

SRSS,0.01,DISP

SOLVE

FINISH

/POST1

/INP,,MCOM

NSORT,U,SUM

*GET,DMAX_SPE,SORT,,MAX!

反应谱最大变形

!

开始地震计算

FINISH

/CONFIG,NRES,20000

/config,fsplit,750这样的话在生成结果文件的时候，单个结果文件最大控制在3G

/SOLU

ANTYPE,TRANS

TRNOPT,FULL

ALPHAD,2*DAMPRATIO*FREQ1*2*3.1415926

BETAD,2*DAMPRATIO/（FREQ1*2*3.1415926）

*DO,I,1,1001

ACEL,AC（I）,0,0

TIME,I*0.02

OUTRES,ALL,ALL

（注意：

若是输出东西格式太大，所以可以仔细利用输出项：

kbc,0  

nsubst,2,10,1　  !

设置子步数  

outres,all,none　!

输出控制  

outres,nsol,last  

outres,v,last  

outres,a,last  

outres,nload,last, 

outres,strs,last  

）

SOLVE

*ENDDO

FINISH

/POST26

FILE

NSOL,2,20,U,X,

PLVAR,2,,,,,,,,,,

STORE,MERGE

XVAR,1

PLVAR,2,

*GET,DMAX_TRANS,VARI,2,EXTREM,VMAX,,!

时程分析最大位移

数据文件

record.txt

-0.0014

-0.0108

-0.0101

-0.0088

-0.0095

-0.012

-0.0142

-0.0128

-0.011

-0.0085

-0.0085

-0.0131

-0.0176

-0.0194

-0.0162

-0.0144

-0.0108

-0.0082

-0.0042

-0.0066

-0.0131

-0.019

-0.0196

-0.0066

0.003

0.0141

-0.0049

-0.0128

-0.0144

-0.0203

-0.026

-0.0325

-0.0306

-0.0172

-0.0197

-0.0163

-0.0164

-0.0067

0.0025

0.015

0.0236

0.0252

0.0336

0.0463

0.0492

0.0419

0.0359

0.0271

0.0235

0.0339

0.0412

0.053

0.0639

0.0732

0.0652

0.0599

0.04

0.04

0.0063

-0.0515

-0.0787

-0.0603

-0.0484

-0.025

-0.0059

0.0134

0.0308

0.0499

0.071

0.0995

0.1219

0.1529

0.1449

0.1155

0.0935

0.0892

0.0926

0.0839

0.0901

0.0993

0.1209

0.0328

-0.1475

-0.2066

-0.1989

-0.2034

-0.1816

-0.1725

-0.1752

-0.1753

-0.1805

-0.163

-0.1347

-0.1087

-0.0782

-0.0429

-0.0017

0.036

0.0785

0.1164

0.1598

0.196

0.2412

0.2729

0.3036

0.32

0.3417

0.2821

0.2324

-0.1198

-0