在 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 里施加地震惯性力的方法。

首先，将三个方向的地震加速度放到一个文本文件里，如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即可对结构进行地震动力分析。

地震响应分析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种方法中就有只是我使用不正确，希望大家能一起来讨论，彻底解决这个问题。

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。

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-> MaterialModels菜单，在Material Model Number 1中添加Structural-> Linear-> Elastic->Isotropic属性，输入材料的弹性模量EX和泊松比PRXY分别为210E9和0.3。

(6)继续给Material Model Number 1添加Density属性，输入密度为7800。

ANSYS软件在结构地震反应分析中所用的方法徐旻洋 1110109132工程体系多自由度运动方程为:[M]{x’’}+[C]{x’}+[K]{x}={F(t)}(1)式中，[M]表示质量矩阵；[C]为阻尼矩阵；[K]为刚度矩阵；{x}为结构体系的位移向量；{F(t)}表示t时刻的载荷向量。

典型的无阻尼模态分析求解的基本方程就是上式(1)的特征值问题：[K]{Φi}=ωi2[M]{Φi} (2)式中，{Φi}是第i阶模态的振型向量(特征向量)，ωi是第i阶模态的固有频率。

ANSYS软件可以求解式(2)，计算结构的固有频率ωi，然后计算相应的振型向量{Φi}，即模态分析。

当式(1)右边{F(t)}是一个已知的谱(如地震反应谱)时，可以用ANSYS软件进行谱分析。

当{F(t)}是任意的随时间变化载荷时，ANSYS软件可进行瞬态动力分析。

ANSYS结构振型分解反应谱分析有如下内容：1）首先要定义好加速度反应谱。

这里需要注意的是，规范上给的是地震影响系数谱曲线，这个曲线的函数值是以地面加速度为单位的。

而我们在用这个软件算的时候就需要给出绝对的加速度值，这个绝对加速度值当然就是要在地震影响系数的基础上再乘上一个地面加速度。

而地面加速度也并不一定是9.8，这与我们使用的单位制有关，如果是N/M/S，就应该是9.8，如果是N/MM/S就应该是9800。

2）求振型。

一定要是相对质量矩阵进行归一化，使用modopt命令默认的方法就可以了。

这个式子是求振型参与系数的，显然这个式子里面不是完整的求振型参与系数的式子，它少了分母，但是，由于对振型相对质量矩阵进行了归一化，这个分母就等于1了，这就是为什么必须要对振型相对质量进行归一化的原因了。

在这一步中，可以这样理解，程序只进行了一次特征值求解，即只求出了周期和振型。

如果需要看某个振型的“内力/应力/反力”，就需要对其进行模态扩展。

3）求谱解。

其实在这一步中，程序只做了一件事，那就是求模态系数。

ANSYS地震分析算例算例描述：在一个城市里，设计师计划建造一座高楼大厦。

然而，由于该地区位于地震活动带，设计师决定对该高楼进行地震分析，以确保其在地震时的稳定性。

假设该高楼大厦的结构为钢筋混凝土框架结构。

设计师希望通过地震分析来确定该结构的最大位移、最大应力和最大应变等参数，并评估结构的稳定性。

地震分析步骤：1.确定模型和边界条件：-在ANSYS中创建高楼大厦的模型。

可以使用ANSYS提供的建模工具，通过绘制位于平面上的轮廓线进行模型创建。

-定义地震加载条件：确定地震波的性质和加载方向，并将其应用于模型的适当位置。

-设置边界条件：确定结构的支撑方式，如固定支座或固定边界条件。

2.材料和单元属性设置：-模型中的材料属性：定义使用的材料的弹性模量、泊松比和密度等特性。

-定义单元属性：选择适当的单元类型，并设置单元的尺寸和属性值。

3.进行静态分析：-应用静态载荷：为了模拟自重和其他永久载荷，将这些载荷应用到模型中。

-执行静态分析：运行ANSYS分析以计算静态应力和变形。

4.开展地震分析：-应用地震波加载：将地震波加载应用于模型的适当位置。

-执行动态分析：运行ANSYS分析以计算结构在地震荷载下的动态响应。

5.结果分析：-输出结果：分析完毕后，ANSYS将提供关于位移、应力、应变和变形等参数的结果。

-结果评估：根据结果评估结构的稳定性和安全性。

可以根据设计准则或标准来判断结构是否合格。

总结：地震分析是建筑和结构设计中至关重要的一步。

ANSYS提供了强大的分析工具，可以帮助工程师和建筑师评估建筑结构在地震荷载下的响应。

通过ANSYS地震分析，设计者可以确定结构的稳定性和安全性，并采取必要的措施来增加其抗震能力。

ANSYS软件在高层建筑地震反应分析中的应用一、概要随着全球地震灾害的频发，高层建筑在地震作用下的安全性受到了广泛关注。

地震反应分析是评估高层建筑在地震作用下结构响应及损坏情况的关键技术手段。

ANSYS软件作为一款广泛应用于工程领域的有限元分析软件，在高层建筑地震反应分析中发挥着重要作用。

1. 抗震设计的重要性随着全球地震活动的增加，高层建筑在其生命周期中的抗震性能显得尤为重要。

高层建筑由于其特殊的结构形式和高度，往往成为地震作用下的关键受力构件。

如果在地震作用下发生倒塌，将造成巨大的人员伤亡和财产损失。

抗震设计的核心目标是确保建筑物在可能发生的地震中能够保持足够的稳定性和完整性，从而保护人员安全并减少财产损失。

ANSYS 软件作为一种强大的有限元分析工具，在高层建筑地震反应分析中发挥着至关重要的作用。

通过ANSYS，工程师们可以模拟和分析建筑物在地震作用下的动态行为，包括应力和变形分布、结构的失效模式以及能量耗散等。

这有助于设计师在建筑设计阶段就识别出潜在的薄弱环节，并采取相应的加固措施来提高建筑的抗震性能。

ANSYS还可以用于验证设计的合理性，通过与其他软件或实验结果的对比，确保建筑物在实际地震中的表现符合预期。

这对于保证高层建筑在地震中的安全性至关重要。

抗震设计是高层建筑安全性的重要保障。

ANSYS软件的应用使得这一过程更加高效、准确，为设计师提供了强有力的工具来应对地震带来的挑战。

2. 高层建筑地震反应分析的挑战随着城市化的加速和土地资源的紧张，高层建筑越来越多。

高层建筑在地震作用下的地震反应分析是一个复杂且具有挑战性的问题。

在地震作用下，高层建筑会受到水平、竖向和扭转等多种振动模态的影响，使得地震反应分析变得非常复杂。

高层建筑结构的耦合作用使得地震反应分析更加困难。

高层建筑结构中，各构件之间存在复杂的相互作用，如梁柱、梁墙、墙柱等。

这些相互作用使得地震力在建筑物内的传递变得复杂，难以准确模拟实际的地震反应。

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里施加地震惯性力的方法[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 软件进行地震分析，分析一个简单的二维框架结构在地震中的响应。

下面是详细的步骤和算例设置：1.几何建模：我们首先在ANSYS中进行几何建模，绘制一个二维的框架结构。

框架由4个节点组成，其中1号和4号节点是固定支座，2号和3号节点是自由节点。

我们可以设置框架的长度、宽度、高度等参数。

2.材料属性：我们需要为框架结构定义材料属性，包括弹性模量和泊松比等。

这些参数可以根据实际的材料特性进行设置。

3.边界条件：我们将1号和4号节点设置为固定支座，以防止结构物在地震中发生位移。

4.地震负荷：我们需要定义地震负荷，即地震的加速度记录。

这些加速度记录可以根据地震现场的实测数据来确定。

在ANSYS中，可以将地震加速度记录分为不同的时程，并将其作为负荷应用在结构上。

5.模型分析：在所有参数设置好后，我们可以进行模型分析。

在ANSYS中，可以选择静力分析或动力分析进行地震分析。

如果选择静力分析，将根据结构物的初始状态和地震负荷计算结构物的响应。

如果选择动力分析，则可以考虑结构物的动态特性和阻尼效应。

6.结果评估：结果评估可以包括结构物的最大位移、最大应力、最大应变等信息，以及结构物的破坏模式和安全性评估等。

在ANSYS中，可以通过可视化和结果输出等方式来进行结果评估。

总结：在这个地震分析算例中，我们使用ANSYS软件对一个二维框架结构进行了地震响应的模拟和分析。

通过设置几何模型、材料属性、边界条件和地震负荷等参数，进行模型分析并评估结构物的性能和安全性。

通过这个算例，我们可以更好地理解地震分析的过程和方法，并为实际工程项目提供参考和指导。

! 【算例】考虑地震作用的建筑物加速度瞬态分析C*** 采用等效实体单元进行建筑物的地震加速度瞬态分析c*** 进行中央电视台新址主楼的等效几何建模*** begin ***FINISH/CLEARc*** 设置参数L=150W=150H=235PI=3.1415926 ! 底座长度! 底座宽度! 整个建筑的高度A=PI*6/180L_HL=L/2T_HL=L/2-TAN(A)*H LO_H=40TO_H=75WW=40/PREP7! 内倾角! 底座长度的1/2 ，作为角点坐标! 顶部长度的1/2 ，作为焦点的坐标! 底座的高度! 顶部悬空部分的高度! “腿”部的宽度/TITLE, EX 9.2(12) by Zeng P, Lei L P, Fang GET,1,SOLID45MP,EX,1,1e10MP,PRXY,1,0.23 MP,DENS,1,1e3 MP,DAMP,1,0.05 ! 几何建模K,,-L_HL,-L_HL, K,,L_HL,-L_HL, K,,L_HL,L_HL, K,,-L_HL,L_HL, K,,-T_HL,-T_HL,H,K,,T_HL,-T_HL,H, K,,T_HL,T_HL,H, K,,-T_HL,T_HL,H,V,1,2,3,4,5,6,7,8 ! 定义单元和材料参数! 由于采用了的简化实体模型，材料参数也进行了等效! 定义8 个角点! 根据8 个角点生成方锥体V1VGEN,2,1, , ,WW,WW,LO_H,, ! V2 (THE COPY OF V1) VGEN,2,1, , ,-WW,-WW,-LO_H,, ! V3 (THE COPY OF V1)VSBV,1,2 VSBV,4,3! V4 (V4=V1-V2) 切掉一个角! NEW_V1 (NEW_V1=V4-V3) 切掉另一个角c*** 为了实现六面体网格划分，需要对实体进行分块处理CSKP,11,0,14,34,10,1,1, ! 通过 3 个点建立局部坐标系WPCSYS,-1,11 VSBW,ALL! 基于局部坐标系建立新的工作面( workplane ) ! 用当前的工作面剖分实体CSKP,12,0,18,30,34,1,1, ! 以下的命令流与前三条的作用相同WPCSYS,-1,12 $VSBW,ALL $CSKP,13,0,27,26,25,1,1,WPCSYS,-1,13 $VSBW,ALL $CSKP,14,0,32,33,34,1,1,WPCSYS,-1,14 $VSBW,ALL $CSKP,15,0,15,31,32,1,1,WPCSYS,-1,15 $VSBW,ALL $CSKP,16,0,29,28,25,1,1,WPCSYS,-1,16 $VSBW,ALLESIZE,L/10! 定义单元尺度 VMESH,ALL ! 对体划分网格CSDELE,11,16,1, ! 删除前面建立的局部坐标系ASEL,S,LOC,Z,,, ! 选择 z=0 处的面，DA,ALL,ALL,, ! 对所选择的面施加完全的约束，即将地面固定 ALLSELFINISHc*** 进行中央电视台新址主楼的等效几何建模*** end *** c*** 进行瞬态分析/SOLU*DIM,ACC,ARRAY,6,4, ! 定义数组，在后面的语句中分别为时间和三个方! 向的加速度赋值(天津波简化数据)ACC(1,1)=0.1 $ACC(1,2)= 0.03599 $ACC(1,3)= 0.24653 $ACC(1,4)= -0.24363ACC(2,1)=0.7 $ACC(2,2)= 0.0293 $ACC(2,3)= 0.14177 $ACC(2,4)= -0.31038ACC(3,1)=1.2 $ACC(3,2)= 0.03695 $ACC(3,3)= 0.71123 $ACC(3,4)= -0.48752ACC(4,1)=1.7 $ACC(4,2)= -0.00725 $ACC(4,3)= -0.2167 $ACC(4,4)= -0.10452ACC(5,1)=2.4 $ACC(5,2)= 0.00734 $ACC(5,3)= -0.81282 $ACC(5,4)= 0.66357ACC(6,1)=3.0 $ACC(6,2)= 0.01759 $ACC(6,3)= -0.00823 $ACC(6,4)= -0.0361ANTYPE,TRANS ! 设置瞬态分析方式TIME,0.001 ! 设置初始加载步，设一个很小的数，但不能为零ACEL,0,0,0 !三个方向都施加零加速度，作为初始计算 SOLVE*DO,TT,1,6,1 !根据不同的时间步，施加加速度给整个结构 TIME,ACC(TT,1)ACEL,ACC(TT,2),ACC(TT,3),ACC(TT,4) ! 三个方向都施加加速度SOLVE *ENDDOc*** 进入时间历程后处理器 /POST26 ! 进入时间 -历程后处理 /AXLAB,x,Time(s) ! /AXLAB,y,Displacement (m) ! NSOL,2,node(-T_HL,-T_HL,H),u,X NSOL,3,node(-T_HL,-T_HL,H),u,Y NSOL,4,node(-T_HL,-T_HL,H),u,Z PLVAR,2,3,4 ! 输出该点位移随时间的变化关系曲线/AXLAB,y,von Mises Equivalent Stress (Pa) ! 定义 y 轴标题ANSOL,5,node(-L_HL+LO_H,L_HL-LO_H*tan(A),LO_H),s,eqv !定义" 腿部 " 角点(受力最!大)等效应力为 y 轴的变量plvar,5 ! 输出该点等效应力随时间的变化关系曲线FINISH!%%%%% %%%% end %%%%%%%%%%%定义 x 轴标识定义 y 轴标识!将建筑物悬空部最高点的位移设置为 2 号变量欢迎您的下载，资料仅供参考!致力为企业和个人提供合同协议，策划案计划书，学习资料等等打造全网一站式需求。

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。

三向输入简‎化后的单向‎输入首先，将三个方向‎的地震加速‎度放到一个‎文本文件里‎，如acce‎x yz.txt，在这个数据‎文件里共放‎三列数据，每列为一个‎方向的地震‎加速度值，这里仅给出‎数据文件中‎前几行的数‎据：-0.22710‎9E-02 -0.20904‎6E+00 0.46707‎2E+01-0.41389‎3E-02 -0.16819‎5E+00 0.26152‎3E+01-0.57475‎3E-02 -0.15789‎0E+00 0.80901‎4E-01-0.73122‎7E-02 -0.15299‎6E+00 0.11997‎5E+01-0.87686‎5E-02 -0.13810‎2E+00 0.13090‎2E+01-0.10106‎7E-01 -0.13158‎2E+00 0.14361‎1E+00 .......................然后，再建一个文‎本文件用来‎存放三个方‎向的地震加‎速度时间点‎，如time‎.txt，在这个数据‎文件里仅一‎列数据，对应于加速‎度数据文件‎里每一行的‎时间点，这里给出数‎据文件中前‎几行数据：0.10000‎0E-010.20000‎0E-010.30000‎0E-010.40000‎0E-010.50000‎0E-010.60000‎0E-01.......................编写如下的‎命令流文件‎，并命名为a‎c ce.inp*dim,ACCEX‎Y Z,TABLE‎,2000,3 !01行*vread‎,A CCEX‎Y Z(1,1),accex‎y z,txt,,JIK,3,2000 !02行(3e16.6) !03行*vread‎,A CCEX‎Y Z(1,0),time,txt !04行(e16.6) !05行ACCEX‎Y Z(0,1)=1 !06行XYZ(0,2)=2 !07行，同上ACCEX‎Y Z(0,3)=3 !08行，同上sh/SOLU9‎ANTYP‎E,trans‎btime‎=0.01 !定义计算起‎始时间e=15.00 !定义计算结‎束时间dtime‎=0.01 !定义计算时‎间步长*DO,itime‎,btime‎,etime‎,dtime‎,itime‎TS,0NSUBS‎T,1, , ,1KBC,1acel,ACCEX‎Y Z(itime‎,1),ACCEX‎Y Z(itime‎,2),ACCEX‎Y Z(itime‎,3) !施加三个方‎向的地震加‎速度SOLVE‎*ENDDO‎最后，在命令窗口‎里输入/input‎,acce,inp即可‎对结构进行‎地震动力分‎析。

三向输入简化后的单向输入首先,将三个方向的地震加速度放到一个文本文件里,如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-01 0.300000E-01 0.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即可对结构进行地震动力分析。