ANSYS Example02地震分析算例 ANSYSWord格式.docx
《ANSYS Example02地震分析算例 ANSYSWord格式.docx》由会员分享，可在线阅读，更多相关《ANSYS Example02地震分析算例 ANSYSWord格式.docx（16页珍藏版）》请在冰豆网上搜索。
（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，依次输入关键点编号和坐标：
关键点1：
坐标（0,0,0）
关键点2：
坐标（0,3,0）
关键点3：
坐标（0,6,0）
关键点4：
坐标（0,9,0）
关键点5：
坐标（0,12,0）
关键点6：
坐标（0,0,10000）
需要说明的是，关键点6为后面建立梁单元所需的截面方向控制点，在上一个例子中已经做过介绍。
（13）完成关键点输入后下面建立直线模型。
在ANSYS主菜单中选择Preprocessor->
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->
Mesh->
Lines，选择所有的直线，完成直线的网格划分。
（17）为了便于后面操作，将网格划分后的单元和节点编号进行适当的清理。
NumberingCtrls->
MergeItems，在MergeCoincidentorEquivalentlyDefinedItems窗口中选择All，清理所有重复的元素。
同样选择NumberingCtrls->
CompressNumber菜单，在CompressNumber中选择All，对节点和单元进行重新编号。
（18）下面建立集中质量单元，采取直接输入单元的方法建立。
Elements->
ElemAttributes，在ElementAttributes中设定单元类型编号为2MASS21，材料编号任意，实参数编号为1。
如图所示
（19）在ANSYS主菜单中选择Preprocessor->
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->
NewAnalysis，选择分析类型为Model
（27）在ANSYS主菜单中选择Solution->
AnalysisOptions，输入模态分析方法为子空间法（Subspace），求解8阶模态，同时需要作模态扩展，扩展的模态为8阶，并计算单元应力和应变，输入窗口如图
（28）再次选择ANSYS主菜单Solution->
CurrentLS，计算当前问题
（29）这时，如果需要看结果，可以进入后处理模块，即ANSYS主菜单GeneralPostproc，可以看到计算的各阶频率和振型。
（30）完成结构自振分析后，下面就可以进行反应谱分析
（31）进入ANSYS主菜单Solution->
NewAnalysis，选择分析类型为Spectrum。
（32）在ANSYS主菜单Solution->
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.45s。
值得注意的是，我国规范给的反应谱横坐标是周期，ANSYS定义的反应谱横坐标是频率，应该注意上述区别。
选择ANSYS主菜单Solution->
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->
SpectrValues，输入对应的反应谱数值依次如下：
0.154350625,
0.191590625,
0.210210625,
0.224175625,
0.250716714,
0.329514922,
0.484352764,
0.598723486,
0.786897371,
1.034212766,
0.3528
（37）最后选择ANSYS主菜单Solution->
ModeCombine，设定振型组合方式为SRSS法，如图所示
（38）选择ANSYS主菜单Solution->
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，行数为NT个，如图所示。
（44）下面从数据文件中读入地震时程记录。
选择ANSYS窗口顶部菜单Parameters->
Readfromfile，在弹出窗口中依次输入将数据读入的数组名称AC，可以用Browse选择目标数据文件，最后要给出数据读入的格式（F8.3），可以参考Fortran相应的文件输入输出要求。
（45）完成数据读入后就可以进行地震反应分析。
由于本次地震分析一共要进行1001步，超过了ANSYS默认的最大1000步限制，因此首先需要进行调整。
在ANSYS的命令输入窗口内，输入以下命令：
Finish
/CONFIG,NRES,20000
（46）进入ANSYS主菜单Solution->
NewAnalysis，指定分析类型为瞬态分析Transient，Solutionmethod可以选择Full。
（47）进入ANSYS主菜单Solution->
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，点击工具栏上第三个按钮，即绘制出节点的位移时程曲线如图
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
REAL,2
E,8
E,14
E,20
!
*SET,NT,1001
*SET,DT,0.02
*DIM,AC,,NT
*VREAD,AC
（1）,RECORD,TXT
（F8.3）
/SOLU
D,1,ALL,
ANTYPE,0
SOLVE
FINISH
/SOLU!
模态分析
ANTYPE,2
MODOPT,SUBSP,8
MXPAND,8,,,1
得到自振频率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
/POST1
/INP,,MCOM
升级会员 