ANSYS地震波的输入
ANSYS_地震分析算例
ANSYS_地震分析算例地震是指地球上因地壳运动而产生的震动现象。
在地震发生后，建筑物的结构稳定性和抗震性能至关重要，因为地震可以对建筑物造成严重破坏。
因此，在建设和设计建筑物时，地震分析变得非常重要。
在此我将介绍一种用ANSYS进行地震分析的算例。
在地震分析中，我们首先需要建立一个合适的模型。
在这个算例中，我们将使用ANSYS提供的有限元分析方法。
首先，我们需要创建一个建筑物的三维模型。
在建筑物的模型中，我们需要包括所有的结构细节，例如建筑物的基础、柱子、梁和地板等。
我们可以使用ANSYS的几何建模工具来创建这个模型。
接下来，我们需要为建筑物定义材料特性。
建筑物的材料特性会对地震分析的结果产生重要影响。
例如，不同种类的混凝土、钢铁和木材等材料在地震作用下的响应是不同的。
我们需要使用ANSYS的材料库来定义这些材料的特性。
完成模型和材料定义后，我们需要定义地震荷载。
地震荷载是指地震发生时作用在建筑结构上的力量。
地震荷载可以根据地震的震级和地震波的性质来确定。
我们可以使用ANSYS的预处理工具来定义这些地震荷载。
接下来，我们需要定义边界条件。
边界条件是指建筑物与外部环境之间的约束关系。
例如，建筑物的基础是固定的，地震发生时不能移动。
我们需要使用ANSYS的加载工具来定义这些边界条件。
完成了模型、材料、地震荷载和边界条件的定义后，我们可以进行地震分析。
地震分析是指通过模拟地震发生时结构的动力响应来评估建筑物的抗震性能。
在ANSYS中，我们可以使用动力分析工具来进行这个分析。
在地震分析过程中，我们可以观察到各个部位的应力和位移等响应。
这些响应可以帮助我们评估建筑物的破坏机制和结构的安全性能。
例如，我们可以观察到柱子是否出现弯曲、梁是否发生裂缝等。
根据地震分析的结果，我们可以对建筑物的设计和结构进行优化。
例如，我们可以调整柱子和梁的尺寸、材质和布置方式，以提高建筑物的抗震能力。
综上所述，通过ANSYS进行地震分析可以帮助我们评估和优化建筑物的抗震性能。
在 ansys 中如何 施加 地震波
三向输入简化后的单向输入首先，将三个方向的地震加速度放到一个文本文件里，如accexyz.txt，在这个数据文件里共放三列数据，每列为一个方向的地震加速度值，这里仅给出数据文件中前几行的数据：-0.227109E-02 -0.209046E+00 0.467072E+01-0.413893E-02 -0.168195E+00 0.261523E+01-0.574753E-02 -0.157890E+00 0.809014E-01-0.731227E-02 -0.152996E+00 0.119975E+01-0.876865E-02 -0.138102E+00 0.130902E+01-0.101067E-01 -0.131582E+00 0.143611E+00 .......................然后，再建一个文本文件用来存放三个方向的地震加速度时间点，如time.txt，在这个数据文件里仅一列数据，对应于加速度数据文件里每一行的时间点，这里给出数据文件中前几行数据：0.100000E-010.200000E-010.300000E-010.400000E-010.500000E-010.600000E-01.......................编写如下的命令流文件，并命名为acce.inp*dim,ACCEXYZ,TABLE,2000,3 !01行*vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000 !02行(3e16.6) !03行*vread,ACCEXYZ(1,0),time,txt !04行(e16.6) !05行ACCEXYZ(0,1)=1 !06行ACCEXYZ(0,2)=2 !07行，同上ACCEXYZ(0,3)=3 !08行，同上finish/SOLUANTYPE,transbtime=0.01 !定义计算起始时间etime=15.00 !定义计算结束时间dtime=0.01 !定义计算时间步长*DO,itime,btime,etime,dtimetime,itimeAUTOTS,0NSUBST,1, , ,1KBC,1acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3) !施加三个方向的地震加速度SOLVE*ENDDO最后，在命令窗口里输入/input,acce,inp即可对结构进行地震动力分析。
ANSYS地震分析实例
ANSYS地震分析实例土木工程中除了常见的静力分析以外，动力分析，特别是结构在地震荷载作用下的受力分析，也是土木工程中经常碰到的题目。
结构的地震分析根据现行抗震规范要求，一般分为以下两类：基于结构自振特性的地震反应谱分析和基于特定地震波的地震时程分析。
本算例将以一个4质点的弹簧－质点体系来说明如何使用有限元软件进行地震分析。
更复杂结构的分析其基本过程也与之类似。
关键知识点：(a) 模态分析(b) 谱分析(c) 地震反应谱输进(d) 地震时程输进(e) 时程动力分析(1) 在ANSYS窗口顶部静态菜单，进进Parameters菜单，选择Scalar Parameters选项，在输进窗口中填进DAMPRATIO=0.02，即所有振型的阻尼比为2%(2) ANSYS主菜单Preprocessor->Element type->Add/Edit/Delete，添加Beam 188单元(3) 在Element Types窗口中，选择Beam 188单元，选择Options，进进Beam 188的选项窗口，将第7个和第8个选项，Stress/Strain (Sect Points) K7, Stress/Strain (Sect Nods) K8，从None 改为Max and Min Only。
即要求Beam 188单元输出积分点和节点上的最大、最小应力和应变(4) 在Element Types 窗口中，继续添加Mass 21集中质量单元(5) 下面输进材料参数，进进ANSYS主菜单Preprocessor->Material Props-> Material Models菜单，在Material Model Number 1中添加Structural-> Linear-> Elastic->Isotropic 属性，输进材料的弹性模量EX和泊松比PRXY分别为210E9和0.3。
ANSYS地震时程分析
在ANSYS里做地震分析时，需要对结构施加地震惯性荷载，地震惯性力是通过加速度的方式输入进结构的，然后与结构的质量一起形成动力计算时的惯性荷载，下面说一下在ANSYS 里施加地震惯性力的方法。
首先，将三个方向的地震加速度放到一个文本文件里，如accexyz.txt，在这个数据文件里共放三列数据，每列为一个方向的地震加速度值，这里仅给出数据文件中前几行的数据：-0.227109E-02 -0.209046E+00 0.467072E+01-0.413893E-02 -0.168195E+00 0.261523E+01-0.574753E-02 -0.157890E+00 0.809014E-01-0.731227E-02 -0.152996E+00 0.119975E+01-0.876865E-02 -0.138102E+00 0.130902E+01-0.101067E-01 -0.131582E+00 0.143611E+00 .......................然后，再建一个文本文件用来存放三个方向的地震加速度时间点，如time.txt，在这个数据文件里仅一列数据，对应于加速度数据文件里每一行的时间点，这里给出数据文件中前几行数据：0.100000E-010.200000E-010.300000E-010.400000E-010.500000E-010.600000E-01.......................编写如下的命令流文件，并命名为acce.inp*dim,ACCEXYZ,TABLE,2000,3 !01行*vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000 !02行(3e16.6) !03行*vread,ACCEXYZ(1,0),time,txt !04行(e16.6) !05行ACCEXYZ(0,1)=1 !06行ACCEXYZ(0,2)=2 !07行，同上ACCEXYZ(0,3)=3 !08行，同上finish/SOLUANTYPE,transbtime=0.01 !定义计算起始时间etime=15.00 !定义计算结束时间dtime=0.01 !定义计算时间步长*DO,itime,btime,etime,dtimetime,itimeAUTOTS,0NSUBST,1, , ,1KBC,1acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3) !施加三个方向的地震加速度SOLVE*ENDDO最后，在命令窗口里输入/input,acce,inp即可对结构进行地震动力分析。
关于ANSYS地震响应分析的一些讨论剖析
地震响应分析1模态组合就是根据模态分析中的几阶振型（也可以少于这几阶，看你要求的精度）进行组合（类似于结构最不利组合），从而求出地震响应的最大值。
2组合各振型反应的最大值，求得结构地震响应的最大值。
这个问题在论坛上已经有很多人问过，也有各种各样的回答，但是至今没有令人满意的解答。
我自己试过很多种方法，加上论坛上其他人提到的方法，大致归类如下：1.先做静力恒载工况分析，打开预应力pstres开关；然后转到时程分析。
结果：恒载对后面的时程计算不起作用，时程计算依然从0开始。
2.直接在antype，trans中考虑恒载：先把timint，off加acel,,9.81,打开应力刚化，sstif，on，lswrite，1，然后timint，on开始时程计算。
结果：恒载9.81起作用了，但结果是错的，它被积分了。
3.不用什么前处理，直接把9.81加在地震波上acel，9.81+ac（i）。
结果，同2，9.81带入了积分，这个9.81相当于阶跃荷载，而不是产生恒载。
4.ansys帮助中施加初始加速度的方法（篇幅限制请自己看帮助）。
结果，同2、3，9.81还是带进时间积分。
5.这种是我受到别人的启发，通过结构受ramp荷载的特点施加的，可以近似的解决问题。
即1）求出结构的自振一阶频率w2）令tr＝1/w3 定义ramp荷载为从0到tr加到9.81，然后在整个时间积分中保持不变4）antype,trans中分几个荷载步将荷载从0加到9.815 在随后的荷载步中acel，，9.81+ac（i）这种做法虽然也是将9.81++加到地震波中，但是因为满足TR的要求，所以这个动力效应被削弱到了静力效应，它作用在结构上就像静载一样。
对于单自由度结构理论上跟静载是完全一样的，但是多自由度会子静力效应上下很小的范围内波动，所以可以认为相当于静载的作用，这样我们就可以达到考虑恒载的目的了。
第5种是我至今为止考虑恒载的做法，我也很想知道还有没有更简单精确的方法，或者在前4种方法中就有只是我使用不正确，希望大家能一起来讨论，彻底解决这个问题。
ANSYS地震分析算例
02 地震分析算例 (ANSYS)土木工程中除了常见的静力分析以外，动力分析，特别是结构在地震荷载作用下的受力分析，也是土木工程中经常遇到的问题。
结构的地震分析根据现行抗震规范要求，一般分为以下两类：基于结构自振特性的地震反应谱分析和基于特定地震波的地震时程分析。
本算例将以一个4质点的弹簧－质点体系来说明如何使用有限元软件进行地震分析。
更复杂结构的分析其基本过程也与之类似。
关键知识点：(a) 模态分析(b) 谱分析(c) 地震反应谱输入(d) 地震时程输入(e) 时程动力分析(1) 在ANSYS窗口顶部静态菜单，进入Parameters菜单，选择Scalar Parameters选项，在输入窗口中填入DAMPRATIO=0.02，即所有振型的阻尼比为2%(2) ANSYS主菜单Preprocessor->Element type->Add/Edit/Delete，添加Beam 188单元(3) 在Element Types窗口中，选择Beam 188单元，选择Options，进入Beam 188的选项窗口，将第7个和第8个选项，Stress/Strain (Sect Points) K7, Stress/Strain (Sect Nods) K8，从None改为Max and Min Only。
即要求Beam 188单元输出积分点和节点上的最大、最小应力和应变(4) 在Element Types 窗口中，继续添加Mass 21集中质量单元(5) 下面输入材料参数，进入ANSYS主菜单Preprocessor->Material Props-> Material Models菜单，在Material Model Number 1中添加Structural-> Linear-> Elastic->Isotropic属性，输入材料的弹性模量EX和泊松比PRXY分别为210E9和0.3。
Ansys-谱分析实例(地震位移谱分析)
二．地震位移谱分析如图所示为一板梁结构，试计算在Y方向地震位移谱作用下的构件响应情况。
板梁结构相关参数见下表所示。
板梁结构几何参数和材料参数板梁结构(模型图)进行题目2的分析。
第一步是建立实体模型（如图4），并选择梁单元和壳单元模拟梁和板进行求解。
建此模型并无特别的难处，只要定义关键点正确，还有就是在建模过程当中注意对全局坐标系的运用，很容易就能做出模型。
此题的难点在于对梁和板的分析求解。
进行求解，首先进行的就是模态分析，约束好六条梁，就可以进行模态的分析求解了。
模态分析后，相应的就进行频谱分析，在输入频率和位移后开始运算求解。
此后进行模态扩展分析，最后进行模态合并分析。
分析完后，再对结果进行查看。
通过命令Main Menu>General Postproc>List Results>Nodal Solution查看节点位移结果、节点等效应力结果（图5）及反作用力结果（图6）。
通过图片我们看清晰的看到梁和板的受力情况及变形情况，在板与梁的连接处，板所受的应力最大，这些地方较容易受到破坏，故可考虑对其进行加固。
而梁主要是中间两层变形较大，所以在设计时应充分考虑材料的选用及直径的大小。
1.指定分析标题1．选取菜单路径Utility Menu | File | Change Jobname，将弹出Change Jobname (修改文件名)对话框。
2．在Enter new jobname (输入新文件名)文本框中输入文字“CH”，为本分析实例的数据库文件名。
单击对话框中的“OK”按钮，完成文件名的修改。
3．选取菜单路径Utility Menu | File | Change Title，将弹出Change Title (修改标题)对话框。
4．在Enter new title (输入新标题)文本框中输入文字“response analysis of a beam-shell structure”，为本分析实例的标题名。
ANSYS地震波的输入
三向输入简化后的单向输入首先，将三个方向的地震加速度放到一个文本文件里，如accexyz.txt，在这个数据文件里共放三列数据，每列为一个方向的地震加速度值，这里仅给出数据文件中前几行的数据：-0.227109E-02 -0.209046E+00 0.467072E+01-0.413893E-02 -0.168195E+00 0.261523E+01-0.574753E-02 -0.157890E+00 0.809014E-01-0.731227E-02 -0.152996E+00 0.119975E+01-0.876865E-02 -0.138102E+00 0.130902E+01-0.101067E-01 -0.131582E+00 0.143611E+00.......................然后，再建一个文本文件用来存放三个方向的地震加速度时间点，如time.txt，在这个数据文件里仅一列数据，对应于加速度数据文件里每一行的时间点，这里给出数据文件中前几行数据：0.100000E-010.200000E-010.300000E-010.400000E-010.500000E-010.600000E-01.......................编写如下的命令流文件，并命名为acce.inp*dim,ACCEXYZ,TABLE,2000,3 !01行*vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000 !02行(3e16.6) !03行*vread,ACCEXYZ(1,0),time,txt !04行(e16.6) !05行ACCEXYZ(0,1)=1 !06行ACCEXYZ(0,2)=2 !07行，同上ACCEXYZ(0,3)=3 !08行，同上finish/SOLUANTYPE,transbtime=0.01 !定义计算起始时间etime=15.00 !定义计算结束时间dtime=0.01 !定义计算时间步长*DO,itime,btime,etime,dtimetime,itimeAUTOTS,0NSUBST,1, , ,1KBC,1acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3) !施加三个方向的地震加速度SOLVE*ENDDO最后，在命令窗口里输入/input,acce,inp即可对结构进行地震动力分析。
在ansys中如何施加地震波
三向输入简化后的单向输入首先，将三个方向的地震加速度放到一个文本文件里，如accex yz.txt，在这个数据文件里共放三列数据，每列为一个方向的地震加速度值，这里仅给出数据文件中前几行的数据：-0.227109E-02 -0.209046E+00 0.467072E+01-0.413893E-02 -0.168195E+00 0.261523E+01-0.574753E-02 -0.157890E+00 0.809014E-01-0.731227E-02 -0.152996E+00 0.119975E+01-0.876865E-02 -0.138102E+00 0.130902E+01-0.101067E-01 -0.131582E+00 0.143611E+00 .......................然后，再建一个文本文件用来存放三个方向的地震加速度时间点，如time.txt，在这个数据文件里仅一列数据，对应于加速度数据文件里每一行的时间点，这里给出数据文件中前几行数据：0.100000E-010.200000E-010.300000E-010.400000E-010.500000E-010.600000E-01.......................编写如下的命令流文件，并命名为ac ce.inp*dim,ACCEXY Z,TABLE,2000,3 !01行*vread,A CCEXY Z(1,1),accexy z,txt,,JIK,3,2000 !02行(3e16.6) !03行*vread,A CCEXY Z(1,0),time,txt !04行(e16.6) !05行ACCEXY Z(0,1)=1 !06行XYZ(0,2)=2 !07行，同上ACCEXY Z(0,3)=3 !08行，同上sh/SOLU9ANTYPE,transbtime=0.01 !定义计算起始时间e=15.00 !定义计算结束时间dtime=0.01 !定义计算时间步长*DO,itime,btime,etime,dtime,itimeTS,0NSUBST,1, , ,1KBC,1acel,ACCEXY Z(itime,1),ACCEXY Z(itime,2),ACCEXY Z(itime,3) !施加三个方向的地震加速度SOLVE*ENDDO最后，在命令窗口里输入/input,acce,inp即可对结构进行地震动力分析。
时程分析中地震波输入位置的讨论
时程分析中地震波输入位置的讨论摘要:时程分析法通过直接动力分析可得到结构相应随时间的变化关系，能真实地反应结构地震相应随时间变化的全过程，是抗震分析的一种重要方法[1]。
目前有限元软件可以实现结构的时程分析，但是在不同的软件中，其实现方式不同，主要区别在地震波的输入位置不同。
本文通过有限元软件ABAQUS采用不同的地震波输入位置对同一结构进行时程分析分析，对比结构相同位置的时程位移曲线，结果表明结构在采用不同地震波输入位置的时程分析中，结构的地震响应基本一致。
关键词：时程分析、有限元软件、钢筋混凝土剪力墙Abstract: The time history analysis method to analyze the available structure through direct power to the relationship between the corresponding changes over time, truly reflect the structure of earthquake corresponding to the whole process of change over time, is an important method of seismic analysis [1]. Finite element software can be time-history analysis of the structure, but in different software in different ways, the main difference between the different positions in the seismic wave input. In this paper the finite element software ABAQUS using different seismic wave input location on the same structure, process analysis analysis, contrast structure the same location of when the process displacement curve, the results show that the structure using different seismic waves enter the position time history analysis, the seismic response basically the same.Keywords: time history analysis, finite element software, reinforced concrete shear walls一、引言在时程分析等动力学问题中，地震力以加速度形式从基础固定处输入。
基于ANSYS的框架结构地震分析教程(静力分析+模态分析+反应谱分析+LS-DYNA时程分析)
!---!---! !
QQ: 1549221758 Website: phipsi.top
-------
!****************************************************************
!**************************************************************** !----------------------------初始化-----------------------------!**************************************************************** finish /clear /FILNAME,ANSYS_Seismic_analysis_Static /Titie,Seismic analysis of the frame structure /PREP7 /DSCALE,ALL,10.0 /eshape,1 /UIS,MSGPOP,3 !进入前处理器 !后处理变形缩放系数 10 倍 !显示单元真实形状 !忽略警告信息
ANSYS中动力边界条件和地震波输入的讨论
在ANSYS里面用输入加速度的方法进行地震波的输入比较简单，现在用的也比较多。
原因可能是大部分朋友分析的都是刚性基础上的结构，或者不考虑地基的情况。
如果考虑上部结构和地基的相互作用，以前的做法就是输入地震加速度，看过很多文章都是简单的这样说一句，但具体怎么输入加速度没有讲(或者简单，不需要讲)，个人认为输入加速度就是通过加速度与质量形成惯性力，而惯性力是体力，有质量的单元在给定的加速度的情况下都会产生惯性力。
因此，在ANSYS里的作法可以通过给定结构各个方向加速度的方法来输入地震波，但问题是如果考虑地基时，地基如果用无质量地基，地基截断边界上加固定边界，这样地基上就不会产生附加惯性力了，即最为古老的无质量地基模型，但现在很多文献都认为无质量地基不能考虑地基辐射阻尼(因为上部结构或自由地基辐射出来的外行散射波在无质量地基的截断固定边界上会发生反射，与实际情况地基无限大不相符)，所得结果欠佳，需要考虑地基的辐射阻尼。
因此，引入了在截断边界上加人工边界条件来模拟由于上部结构或地表自由面产生的散射波向无限域地基辐射，避免在固定边界上产生反射。
但是，当在截断边界上采用人工边界时，地震波的输入方式受人工边界形式的影响，采用透射人工边界时，在人工边界上是加的地震波的位移时程曲线；采用粘弹性人工边界时，是通过在人工边界上施加等效结点荷载的方式实现的，这时，要在人工边界上加位移和速度时程曲线。
最近看了刘晶波老师的一篇文章“粘弹性人工边界及地震动输入在通用有限元软件中的实现”，里面讲了在ANSYS中实现粘弹性人工边界，对此很感兴趣，里面讲了一致粘弹性边界可以通过在边界上加COMBIN14单元来模拟，由于粘弹性人工边界就是在边界上加并联的弹簧和阻尼器来实现的，而ANSYS中的COMBIN14单元就是并联的弹簧和阻尼器单元，因此，可以借助于此单元来实现粘弹性人工边界条件；对于他们提出的等效人工边界单元，可以直接用实体单元来模拟，但实体单元的的弹性参数和泊松比采用与弹簧+阻尼单元相等效的值，单元的密度给一个很小的值。
ANSYS加载命令流
一，在ANSYS里施加地震惯性力的方法在ANSYS里做地震分析时，需要对结构施加地震惯性荷载，地震惯性力是通过加速度的方式输入进结构的，然后与结构的质量一起形成动力计算时的惯性荷载，下面说一下在ANSYS 里施加地震惯性力的方法。
首先，将三个方向的地震加速度放到一个文本文件里，如accexyz.txt，在这个数据文件里共放三列数据，每列为一个方向的地震加速度值，这里仅给出数据文件中前几行的数据：-0.227109E-02 -0.209046E+00 0.467072E+01-0.413893E-02 -0.168195E+00 0.261523E+01-0.574753E-02 -0.157890E+00 0.809014E-01-0.731227E-02 -0.152996E+00 0.119975E+01-0.876865E-02 -0.138102E+00 0.130902E+01-0.101067E-01 -0.131582E+00 0.143611E+00 .......................然后，再建一个文本文件用来存放三个方向的地震加速度时间点，如time.txt，在这个数据文件里仅一列数据，对应于加速度数据文件里每一行的时间点，这里给出数据文件中前几行数据：0.100000E-010.200000E-010.300000E-010.400000E-010.500000E-010.600000E-01.......................编写如下的命令流文件，并命名为acce.inp*dim,ACCEXYZ,TABLE,2000,3 !01行*vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000 !02行(3e16.6) !03行*vread,ACCEXYZ(1,0),time,txt !04行(e16.6) !05行ACCEXYZ(0,1)=1 !06行ACCEXYZ(0,2)=2 !07行，同上ACCEXYZ(0,3)=3 !08行，同上finish/SOLUANTYPE,transbtime=0.01 !定义计算起始时间etime=15.00 !定义计算结束时间dtime=0.01 !定义计算时间步长*DO,itime,btime,etime,dtimetime,itimeAUTOTS,0NSUBST,1, , ,1KBC,1acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3) !施加三个方向的地震加速度SOLVE*ENDDO最后，在命令窗口里输入/input,acce,inp即可对结构进行地震动力分析。
二维黏弹性人工边界单元及地震波输入在ANSYS中的实现
水利与建筑工程学报
ＪｏｕｒｎａｌｏｆＷａｔｅｒＲｅｓｏｕｒｃｅｓａｎｄＡｒｃｈｉｔｅｃｔｕｒａｌＥｎｇｉｎｅｅｒｉｎｇ
ＤＯＩ：１０．３９６９／ｊ．ｉｓｓｎ．１６７２－１１４４．２０１９．０５．０１０
Ｖｏｌ．１７Ｎｏ．５ Ｏｃｔ．，２０１９
二维黏弹性人工边界单元及地震波输入 在 ＡＮＳＹＳ中的实现
对无限域或半无限域进行地震响应分析时，往 往采用有限元法。该方法通常采取的处理措施是在 无限域或半无限域空间中人为地截取一定范围的有 限域，进而可以将无限区域转化为有限区域，然后在 有限域的基础上，对其进行离散法处理，进而无限域
或半无限域最终被转化为有限个单元。同时为了达 到近似处理的效果，通常在截取的有限域的边界上 人为的施加约束，但是这一处理的缺点是在处理过 程中，不仅没有考虑远场无限地基的辐射阻尼效应， 同时半空间无限地基弹卷
截取的有限域的范围存在不确定性，随着有限域截 取范围的增加，自由度也随之增加。基于此离散范 围也更加宽广，进而网格划分的单元数及结点数成 倍增加，因此，在地震动的动力计算中工作量也进一 步被增加，对计算机的要求也更高，进而，在动力学 中非线性问题的求解难度进一步增大。目前，针对 该关键性问题处理的诸多措施中，较好的处理方法 是在截取的有限域边界处设置局部人工边界。目前 在工程中较常用的局部边界主要有黏性边界［１］、透 射边界［２］、黏弹性边界［３］等。黏性边界的优点是物 理概念清晰，在实际动力计算程序中实现起来相对 比较容易，对于 １阶动力响应精确度较高。基于此， 黏性边界单元被嵌入到了很多大型商用软件中以便 于工程的实际应用，如 ＬＳＤＹＮＡ、ＦＩＡＣ、ＡＢＡＱＵＳ［４］ 等。它的缺点是在施加人工边界过程中只着重考虑 了人工边界辐射阻尼效应，忽略了人工边界处半无 限空间介质的弹性恢复性能。因此在实际工程应用 中容易发生整体漂移的现象，进而使计算结果精度 不高。透射边界的优点是对于 ２阶动力响应计算结 果存在较高的精度，其缺点是动力计算分析程序的 代码编写较复杂、在实际工程应用中可能存在高频 失稳的现象。黏弹性边界不仅能有效的模拟半无限 空间介质的辐射阻尼效应，同时还能较好的模拟半 无限地基的弹性恢复性能。此外还具有良好的频率 稳定性和较高的精度。因此在实际工程中得到了广 泛的应用［５－１０］。
ANSYS软件在高层建筑地震反应分析中的应用
ANSYS软件在高层建筑地震反应分析中的应用一、概要随着全球地震灾害的频发，高层建筑在地震作用下的安全性受到了广泛关注。
地震反应分析是评估高层建筑在地震作用下结构响应及损坏情况的关键技术手段。
ANSYS软件作为一款广泛应用于工程领域的有限元分析软件，在高层建筑地震反应分析中发挥着重要作用。
1. 抗震设计的重要性随着全球地震活动的增加，高层建筑在其生命周期中的抗震性能显得尤为重要。
高层建筑由于其特殊的结构形式和高度，往往成为地震作用下的关键受力构件。
如果在地震作用下发生倒塌，将造成巨大的人员伤亡和财产损失。
抗震设计的核心目标是确保建筑物在可能发生的地震中能够保持足够的稳定性和完整性，从而保护人员安全并减少财产损失。
ANSYS 软件作为一种强大的有限元分析工具，在高层建筑地震反应分析中发挥着至关重要的作用。
通过ANSYS，工程师们可以模拟和分析建筑物在地震作用下的动态行为，包括应力和变形分布、结构的失效模式以及能量耗散等。
这有助于设计师在建筑设计阶段就识别出潜在的薄弱环节，并采取相应的加固措施来提高建筑的抗震性能。
ANSYS还可以用于验证设计的合理性，通过与其他软件或实验结果的对比，确保建筑物在实际地震中的表现符合预期。
这对于保证高层建筑在地震中的安全性至关重要。
抗震设计是高层建筑安全性的重要保障。
ANSYS软件的应用使得这一过程更加高效、准确，为设计师提供了强有力的工具来应对地震带来的挑战。
2. 高层建筑地震反应分析的挑战随着城市化的加速和土地资源的紧张，高层建筑越来越多。
高层建筑在地震作用下的地震反应分析是一个复杂且具有挑战性的问题。
在地震作用下，高层建筑会受到水平、竖向和扭转等多种振动模态的影响，使得地震反应分析变得非常复杂。
高层建筑结构的耦合作用使得地震反应分析更加困难。
高层建筑结构中，各构件之间存在复杂的相互作用，如梁柱、梁墙、墙柱等。
这些相互作用使得地震力在建筑物内的传递变得复杂，难以准确模拟实际的地震反应。
地震接触命令流ansys
地震接触命令流ansys本人现在在做一个地震波的时程分析，经过计算，是收敛的，有位移时程结果，但是没有加速度及速度的时程结果。
请各位高手指点。
谢谢！望大家一起交流，共同进步。
命令流如下，大概需要计算1个小时。
fini/clear/COM,Structural/prep7!单元定义!土及桩所用单元et,1,plane42 !平面单元keyopt,1,3,2 !打开plane42的平面应变选项r,1!边界所用弹簧单元et,2,14keyopt,2,3,2r,2,826!***场地土的材料模型**!***上层4.5m粉质砂土层mp,dens,1,2000mp,ex,1,2.5e8mp,nuxy,1,0.4mp,mu,1,0.3TB,DP,1TBdata,1,1.2e4,13!***下层24m基岩mp,ex,2,1.8e9mp,nuxy,2,0.25mp,mu,2,0.4TB,DP,2TBdata,1,8e4,32!***下部换土层mp,dens,3,2600mp,ex,3,1.8e9mp,nuxy,3,0.25mp,mu,3,0.4TB,DP,3TBdata,1,8e4,32!***四周换上层mp,dens,4,2000mp,ex,4,2.5e8mp,nuxy,4,0.4mp,mu,4,0.3TB,DP,4TBdata,1,1.2e4,13!***四周换下层mp,dens,5,2600mp,ex,5,1.8e9mp,nuxy,5,0.25mp,mu,5,0.4TB,DP,5TBdata,1,8e4,32!***混凝土桩、承台及墩****mp,ex,6,3e10mp,nuxy,6,0.2!***平面自由场地层几何模型** !***第一层土*** RECTNG,-9,-4.5,28.5,24 RECTNG,-23,-11,28.5,24!***第二层土*** RECTNG,-2,0,24,12 RECTNG,-4.5,-4,24,12 RECTNG,-9,-4.5,24,12 RECTNG,-23,-11,24,0 RECTNG,-2,0,12,8 RECTNG,-4,-2,12,8 RECTNG,-9,-4,12,8 RECTNG,-11,-9,12,0 RECTNG,-9,-6,8,6 RECTNG,-9,0,4,0 RECTNG,-23,0,0,-16!***下部换土层*** RECTNG,-6,0,8,6 RECTNG,-9,0,6,4!***四周换土层***!四周上部换土层RECTNG,-11,-9,28.5,24!四周下部换土层RECTNG,-11,-9,24,12!**桩*******rectng,-4,-2,24,12!**承台******* rectng,-4.5,-4,28.5,24 rectng,-4,-3,28.5,24k,100,-3,28.5,0k,101,-3,24,0k,102,-2,24,0k,103,-1.65,28.5,0a,100,101,102,103k,104,-1.65,28.5,0k,105,-2,24,0k,106,0,24,0k,107,0,28.5,0a,104,105,106,107!**墩*******rectng,-3,-1.65,38.5,28.5!生成对称面FLST,3,23,5,ORDE,2 FITEM,3,1FITEM,3,-23 ARSYM,X,P51X, , , ,0,0!布尔运算!土的布尔allselFLST,2,34,5,ORDE,4 FITEM,2,1FITEM,2,-17FITEM,2,24FITEM,2,-40AGLUE,P51X!桩、承台及墩的布尔allselFLST,2,12,5,ORDE,4 FITEM,2,18FITEM,2,-23FITEM,2,41FITEM,2,-46AGLUE,P51X!有限网格划分!划分网格lsel,s,loc,y,29,38 !墩lesize,all,,,10lsel,s,loc,y,25,28 !包括部分承台lesize,all,,,5lsel,s,loc,y,13,22lesize,all,,,12lsel,s,loc,y,8.5,11.5lesize,all,,,4lsel,s,loc,y,6.5,7.5lsel,a,loc,y,4.5,5.5lesize,all,,,2lsel,s,loc,y,0.5,3.5 lesize,all,,,4lsel,s,loc,x,-23 lsel,a,loc,x,23lsel,r,loc,y,0,24 lesize,all,,,24lsel,s,loc,x,-11 lsel,a,loc,x,11lsel,r,loc,y,0,12 lesize,all,,,12lsel,s,loc,x,-11.6,-22.5 lsel,a,loc,x,11.6,22.5 lesize,all,,,12lsel,s,loc,x,-9.6,-10.5 lsel,a,loc,x,9.6,10.5 lesize,all,,,2lsel,s,loc,y,12,28.5 lsel,r,loc,x,4.6,8.5 lesize,all,,,5lsel,s,loc,y,12,28.5 lsel,r,loc,x,-4.6,-8.5 lesize,all,,,5lsel,s,loc,y,12,28.5 lsel,r,loc,x,4.1,4.4 lesize,all,,,1lsel,s,loc,y,12,28.5 lsel,r,loc,x,-4.1,-4.4 lesize,all,,,1lsel,s,loc,y,24,28.5lesize,all,,,1 lsel,s,loc,y,24,28.5 lsel,r,loc,x,-3.1,-3.9 lesize,all,,,1 lsel,s,loc,y,24,38.5 lsel,r,loc,x,2.1,2.9 lesize,all,,,2 lsel,s,loc,y,24,38.5 lsel,r,loc,x,-2.1,-2.9 lesize,all,,,2 lsel,s,loc,y,8,28.5 lsel,r,loc,x,0.1,1.6 lesize,all,,,2 lsel,s,loc,y,8,28.5 lsel,r,loc,x,-0.1,-1.6 lesize,all,,,2 lsel,s,loc,y,8,12 lsel,r,loc,x,2.1,3.9 lesize,all,,,3 lsel,s,loc,y,8,12 lsel,r,loc,x,-2.1,-3.9 lesize,all,,,3 lsel,s,loc,y,8,11 lsel,r,loc,x,4.1,5.9 lesize,all,,,3 lsel,s,loc,y,8,11 lsel,r,loc,x,-4.1,-5.9 lesize,all,,,3 lsel,s,loc,y,6,8lesize,all,,,3lsel,s,loc,y,6,8lsel,r,loc,x,-6.1,-8.9lesize,all,,,3lsel,s,loc,y,6lsel,r,loc,x,0.1,5.9lesize,all,,,8lsel,s,loc,y,6lsel,r,loc,x,-0.1,-5.9lesize,all,,,8lsel,s,loc,y,0,4lsel,r,loc,x,0.1,5.9lesize,all,,,11lsel,s,loc,y,0,4lsel,r,loc,x,-0.1,-5.9lesize,all,,,11lsel,s,loc,y,-0.1,-15.4lesize,all,,,8lsel,s,loc,y,-16lesize,all,,,25!第一层土划分网格allselasel,s,,,1,2asel,a,,,24,25aatt,1,1,1,0 !aatt,mat,real,type,esys amesh,all!第二层土划分网格asel,s,,,3asel,a,,,47asel,a,,,49,55asel,a,,,57asel,a,,,58,60,2asel,a,,,61,63,2asel,a,,,64,68,2asel,a,,,69,71,2asel,a,,,72,74,2asel,a,,,75aatt,2,1,1,0 !aatt,mat,real,type,esys amesh,all!下部换土层划分网格allsel,allasel,s,,,67,73,3asel,a,,,65aatt,3,1,1,0 !aatt,mat,real,type,esys amesh,all!四周换土层划分网格!四周上层换土层划分网格allsel,allasel,s,,,48asel,a,,,56aatt,4,1,1,0 !aatt,mat,real,type,esys amesh,all!四周下层换土层划分网格asel,s,,,59asel,a,,,62aatt,5,1,1,0 !aatt,mat,real,type,esysamesh,all!桩、承台及墩划分网格allsel,allasel,s,,,4,15aatt,6,1,1,0 !aatt,mat,real,type,esysamesh,all!******接触分析******et,3,169et,4,171KEYOPT,4,12,0 !选择接触类型(CONTA171)为STANDARD r,3,1,0,0.039,0.1!RMORE,,,2.6E4 !设置实常数,最大接触摩擦26kPa!定义桩的刚性接触面（目标面）lsel,s,line,,76lsel,a,line,,15lsel,a,line,,72lsel,a,line,,70lsel,a,line,,69lsel,a,line,,23lsel,a,line,,26lsel,a,line,,162lsel,a,line,,164lsel,a,line,,161lsel,a,line,,31lsel,a,line,,168CM,target,LINElplottype,3real,3nsll,s,1 !选择所有在已选线上的节点esurf,all !在已选择的单元的自由面上产生新单元!esel,s,type,,3 !选择3＃类型的单元!eplot!定义土的柔性接触面（接触面）lsel,s,line,,2lsel,a,line,,194lsel,a,line,,14lsel,a,line,,212lsel,a,line,,12lsel,a,line,,11lsel,a,line,,193lsel,a,line,,104lsel,a,line,,219lsel,a,line,,106lsel,a,line,,207lsel,a,line,,94CM,contact,LINElplottype,4real,3nsll,s,1 !选择所有在已选线上的节点esurf,all !在已选择的单元的自由面上产生新单元!esel,s,type,,4 !选择4＃类型的单元!eplotALLSEL,ALLFINISH!****求解*******!**施加约束**/SOLUnsel,s,loc,x,-23nsel,a,loc,x,23D,ALL,all,0nsel,s,loc,y,-16D,ALL,uy,0ALLSEL/replot,fastfinisave!****以上几何建模完成*******!****以下开始求解******!***模态分析***为了计算基频！/soluantype,modal!pstres,offmodopt,lanb,20,0,0,,off !选择模态分析方法mxpand,20,,,yes !扩展振型outpr,basic,all !控制求解屏幕输出outres,all,all !输出所有项，每一步都输出allselsolve!***通用后处理***!/post1!set,list !显示各阶频率!set,last !第一阶模态!pldisp,1!fini!***地震时程分析***（只有x方向的振动）!***将地震数据存在ff文件中****fini/config,nres,20000 !定义最大子步（只在还没有进行求解时才能设置）/config,fsplit,750 !由于结果文件太大(>4G)，指定每个结果文件大小为750*4=3000M!/config,statu !显示修改后的值/soluallselantype,trans !指定为时程分析（采用荷载步法，而非连续法）dampratio=0.239pi=3.1415freq1=6.0590*2*pifreq2=6.8381*2*pi!得到阻尼系数（质量系数和刚度系数）alphad,2*dampratio*freq1*freq2/(freq1+freq2)betad,2*dampratio/(freq1+freq2)trnopt,full !指定为完全瞬态分析timint,off !关闭时间积分开关outres,basic,all !输出基本项，每一项都输出nlgeom,on !打开大变形开关!ncnv,0 !即使不收敛也不中断分析!sstif,offtime,1e-8 !设置一个十分小的时间步kbc,1 !设置加载方式为阶跃式nsubst,2 !设置两个子步acel,0,9.8,0 !施加重力加速度allslswr,1 !把这个写入第一子步!solve !求解!计算在重力作用下，结构的地震响应timint,on !打开时间积分效应!outres,basic,all !输出基本项，每一项都输出kbc,0 !设置加载方式为线性加载nsubst,10 !设置15个加载子步!定义和读入时程曲线nt=250 !时程曲线有1500个点dt=0.02 !时间步长*dim,acex,array,1,nt,0,,, !1t行nt列的Elcentro地震数据文件(1,nt)*creat,ff !创建宏文件ff*vread,acex(1),acex,txt,,nt(e11.3)*end!***将地震数据导入***/input,ff控制结果输出outpr,all,none !不输出屏幕打印信息!outres,all,all !不输出输出结果信息outres,nsol,all !仅输出节点自由度解!求解器选项nlgeom,onlnsrch,on !打开线性搜索sstif,on !激活在非线性分析中的主应力刚度效应,考虑大变形时默认为打开nropt,modi !自动选择牛顿－莱布利兹收敛lumpm,on !使用集中质量矩阵，容易收敛eqslv,SPARSE !可由程序自动选择（选择求解方式稀疏矩阵，一致矩阵，雅克比矩阵等）当可能出现病态时，用稀疏矩阵!***输入地震波求解*****DO,tm,1,nt,1TIME,dt*tm!***选择地震波加速度施加的位置****nsel,s,loc,y,-16 !选择土体底部的所有节点作为地震波的输入位置ACEL,acex(1,tm)allsel,alllswr,tm!SOLVE*ENDDO*do,i,1,ntlssolve,i*enddosave。
在ANSYS里施加地震惯性力的方法
[分享]在ANSYS里施加地震惯性力的方法[b][size=16px]在ANSYS里施加地震惯性力的方法[/size][/b]在ANSYS里做地震分析时，需要对结构施加地震惯性荷载，地震惯性力是通过加速度的方式输入进结构的，然后与结构的质量一起形成动力计算时的惯性荷载，下面说一下在ANSYS 里施加地震惯性力的方法。
首先，将三个方向的地震加速度放到一个文本文件里，如accexyz.txt，在这个数据文件里共放三列数据，每列为一个方向的地震加速度值，这里仅给出数据文件中前几行的数据：[font=宋体] -0.227109E-02 -0.209046E+00 0.467072E+01-0.413893E-02 -0.168195E+00 0.261523E+01-0.574753E-02 -0.157890E+00 0.809014E-01-0.731227E-02 -0.152996E+00 0.119975E+01-0.876865E-02 -0.138102E+00 0.130902E+01-0.101067E-01 -0.131582E+00 0.143611E+00[/font] .......................然后，再建一个文本文件用来存放三个方向的地震加速度时间点，如time.txt，在这个数据文件里仅一列数据，对应于加速度数据文件里每一行的时间点，这里给出数据文件中前几行数据：[font=宋体] 0.100000E-010.200000E-010.300000E-010.400000E-010.500000E-010.600000E-01[/font].......................编写如下的命令流文件，并命名为acce.inp[font=宋体]*dim,ACCEXYZ,TABLE,2000,3 !01行*vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000 !02行[/font][font=宋体](3e16.6) !03行[/font][font=宋体]*vread,ACCEXYZ(1,0),time,txt !04行(e16.6) !05行[/font][font=宋体]ACCEXYZ(0,1)=1 !06行ACCEXYZ(0,2)=2 !07行，同上ACCEXYZ(0,3)=3 !08行，同上[/font]finish[font=宋体]/SOLUANTYPE,transbtime=0.01 !定义计算起始时间etime=15.00 !定义计算结束时间dtime=0.01 !定义计算时间步长*DO,itime,btime,etime,dtimetime,itime[/font][font=宋体][font=宋体]AUTOTS,0NSUBST,1, , ,1KBC,1[/font]acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3) !施加三个方向的地震加速度SOLVE*ENDDO[/font]最后，在命令窗口里输入/input,acce,inp即可对结构进行地震动力分析。
ANSYS地震分析算例
ANSYS地震分析算例地震分析是通过模拟地震波在结构体系中传播和反应的过程，来评估结构的抗震性能。
ANSYS提供了丰富的工具和功能来支持地震分析，包括地震波输入、地震响应计算和结构的抗震设计。
接下来，我们将介绍一个ANSYS地震分析的算例，来说明如何使用ANSYS进行地震分析。
首先，我们需要定义地震波的输入。
在ANSYS中，可以通过加载事先记录的地震波时程数据来模拟地震波。
这些地震波数据可以从观测站或数字模拟中获取。
通过加载地震波数据，可以将地震波的荷载施加在相应的结构上。
其次，我们需要建立地震分析的数值模型。
在ANSYS中，可以使用各种元素和材料模型来表示结构。
对于地震分析，通常会使用3D有限元模型。
在建立数值模型时，需要根据实际情况定义结构的几何形态和材料特性。
建议使用精细的网格划分来确保模型的准确性和可靠性。
然后，我们需要设置地震分析的边界条件。
这包括定义结构的支撑条件、荷载施加方式以及结构的初始条件等。
在地震分析中，结构通常会受到来自地震波的水平和垂直方向两个方向上的振动力。
因此，需要设置适当的支撑条件和加载方式来模拟地震波对结构的影响。
接着，我们可以进行地震分析计算。
在ANSYS中，可以使用不同的求解方法来进行地震分析，包括静力分析、模态分析和时程历程分析。
静力分析适用于弹性结构，可以用来评估结构在地震荷载下的变形和应力分布。
模态分析可以计算结构的固有频率和振型，并用于评估结构的抗震性能。
时程历程分析是一种更为准确的地震分析方法，可以模拟地震波在结构中的传播和反应的过程。
最后，我们可以进行地震分析结果的后处理。
在ANSYS中，可以使用各种功能来对地震分析的结果进行可视化和分析。
可以绘制结构的变形图、应力分布图和振动模态图，以评估结构的抗震性能。
此外，还可以计算结构的位移响应和应力峰值，以更详细地评估结构的动力响应。
总结起来，ANSYS提供了一个完整的地震分析解决方案，可以用于评估结构的抗震性能。
ANSYS地震波的输入
对于地震波的输入，可以把荷载记录作成文件，利用apdl的读取功能读入数据库中。
下面的例子是自己编的一个小文件。
修改一下可以更简洁。
Fini/config,nres,1000*dim,aceX,TABLE,3000,1*dim,aceY,TABLE,3000,1*dim,aceZ,TABLE,3000,1*creat,ff*vread,aceX(1,1),acex,txt,,1(e16.6)*vread,aceX(1,0),acexTT,txt,,,1(e16.6)ACEX(0,1)=1*end/input,ff*creat,ff*vread,aceY(1,1),txt,,1(e16.6)*vread,aceY(1,1),ACETT,,,1(e17.6)ACEY(0,1)=1*end/input,ff*creat,ff*vread,aceZ(1,1),txt,,1(e16.6)*vread,aceZ(1,0),ACETT,,,1(e17.6)ACEZ(0,1)=1*end/input,ff!地震波时程记录分成了3个文件，每个文件是一列。
分别记录x，y，z方向的加速度。
Accett是时间记录。
这样就可以把加速度记录读取到ansys数据库中作为数组。
也可以把加速度记录作成一个文件，这样程序就简单多了。
下面是计算部分语句：/SOLUANTYPE,trans!求解其自己选了TM_START=0.01TM_END=15.00TM_INCR=0.01*DO,TM,TM_START,TM_END,TM_INCRTIME,tmalpha,BETAD,ACEL,acex(tm),acey(tm),acez(tm) SOLVE*ENDDOfini。
地震波动输入方法笔记
地下结构地震动输入方法：
1．将地震波动输入问题化为波源问题，即将输入地震动转化为作用于人工边界上的等效荷载，使人工边界产生与自由波场相同的位移、速度和应力。
（刘晶波）
2．沿用地面建筑物的地震动输入方式，即模型边界采用固定边界，仅从计算模型底部输入地震波。
3．直接加速度输入法（不适用于多点激励）；直接位移输入法；大质量法（large mass method, LMM）（适用于长周期结构计算），大刚度法（large stiffness method, LSM）
4．等效荷载法，直接求解法，相对运动法，大质量法。
（杨斌斌）
5．拉格朗日乘子法（于海丰）。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
三向输入简化后的单向输入
首先，将三个方向的地震加速度放到一个文本文件里，如accexyz.txt，在这个数据文件里共放三列数据，每列为一个方向的地震加速度值，这里仅给出数据文件中前几行的数据：-0.227109E-02 -0.209046E+00 0.467072E+01
-0.413893E-02 -0.168195E+00 0.261523E+01
-0.574753E-02 -0.157890E+00 0.809014E-01
-0.731227E-02 -0.152996E+00 0.119975E+01
-0.876865E-02 -0.138102E+00 0.130902E+01
-0.101067E-01 -0.131582E+00 0.143611E+00 .......................
然后，再建一个文本文件用来存放三个方向的地震加速度时间点，如time.txt，在这个数据文件里仅一列数据，对应于加速度数据文件里每一行的时间点，这里给出数据文件中前几行数据：
0.100000E-01
0.200000E-01
0.300000E-01
0.400000E-01
0.500000E-01
0.600000E-01
.......................
编写如下的命令流文件，并命名为acce.inp
*dim,ACCEXYZ,TABLE,2000,3 !01行
*vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000 !02行(3e16.6) !03行
*vread,ACCEXYZ(1,0),time,txt !04行
(e16.6) !05行
ACCEXYZ(0,1)=1 !06行
ACCEXYZ(0,2)=2 !07行，同上
ACCEXYZ(0,3)=3 !08行，同上
finish
/SOLU
ANTYPE,trans
btime=0.01 !定义计算起始时间
etime=15.00 !定义计算结束时间
dtime=0.01 !定义计算时间步长
*DO,itime,btime,etime,dtime
time,itime
AUTOTS,0
NSUBST,1, , ,1
KBC,1
acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3) !施加三个方向的地震加速度SOLVE
*ENDDO
最后，在命令窗口里输入/input,acce,inp即可对结构进行地震动力分析。
说明和讨论：
1、命令流中各行说明：
01行：定义2000行，3列的数组，(行数根据数据文件里加速度点数来定)ACCEXYZ用来存放三个方向的加速度值。
02行：从数据文件accexyz.txt里读加速度值到数组ACCEXYZ,2000为行数，可根据情况修改。
其中的JIK，3，2000非常重要，它决定着将从加速度数据文件中数据输入到ACCEXYZ数组时的存放格式，这里用到按JIK方式变化,即读进来的数据依次放入ACCEXYZ(I,J)（J从1到3，I从1到2000，K默认从1到1）。
根据这行命令下面(3e16.6)的格式，每次从数据文件里读一行三个数据，分别放入ACCEXYZ(1,1)、ACCEXYZ(1,2)、ACCEXYZ(1,3)，接着再读下一行的三个数据，分别放入ACCEXYZ(2,1)、ACCEXYZ(2,2)、ACCEXYZ(2,3)，依次类推。
当然依据数据文件的格式，也可以选用IJK, IKJ, JIK, JKI, KIJ, KJI等格式，其中IJK为默认。
03行：读数据的格式，每行三个数值。
由于数据文件中的数据是用科学记数年法表示的，因此，这里也用相应FORTRAN的科学记数法的格式。
如果数据文件里的数值是如“0.2876 1.2333 2.9938”这样的格式，此行的数据格式也就修改为“(3f10.4)”这样的格式。
04行：从数据文件time.txt里读时间值到数组ACCEXYZ的第零列。
06行：将数组ACCEXYZ的第零行赋值，如果不对行插值的话也可以不赋值
例2：
（转自/question/120580469.html?fr=qrl&cid=89&index=1&fr2=query）地震波为一个单列的txt文档（不要有时间列，时间间隔假定0.02s）
命令流（下面的命令必须用文本读入的方法实现，命令窗无法实现）：
*DIM,EQ,,6000,1,,, !首先定义一个EQ数组
*VREAD,EQ,earthquake,TXT,,,,,,0, !将earthquake.txt读入EQ数组里
(e16.0)
用的时候：
*DO,EQ,0.02,0.02*6000,0.02
TIME,TM
acel,EQ
!求解过程...。