ABAQUS瑞利阻尼
瑞利阻尼在Abaqus中的实现
瑞利阻尼在Abaqus中的实现作者：王飞王衔陈涛顾浩声来源：《计算机辅助工程》2018年第05期摘要：为充分体现结构阻尼对结构振动响应的影响，在结构动力响应有限元分析中加入经典的瑞利阻尼，并在Abaqus中进行匀质梁的弯曲和轴向简谐振动分析，然后与结构动力学理论解进行对比，研究瑞利阻尼在Abaqus中的实现方法。
结果表明：有限元动力响应分析结果与理论结果吻合良好，该方法可以准确实现在Abaqus中添加瑞利阻尼。
关键词：瑞利阻尼；结构阻尼；匀质梁；简谐振动；动力响应中图分类号：TU311.3文献标志码：B0 引言结构振动不可避免地会受到阻尼的作用[1]，能量在振动中不断耗散。
经典阻尼是结构所有部位都具有相似阻尼机制时的一种合理抽象，瑞利阻尼是根据振型阻尼比建立结构经典阻尼矩阵的重要方法。
[2]作为结构动力学分析的重要手段，有限元在复杂结构中应用广泛。
因此，本文探究如何在商业有限元软件中实现瑞利阻尼，并利用经典的结构动力学算例，证明该方法的可行性。
1 瑞利阻尼瑞利阻尼假定結构阻尼与质量矩阵和刚度矩阵的组合成一定比例，是对质量比例阻尼和刚度比例阻尼的一种改进[3]，其认为圆频率与阻尼比的关系式为瑞利阻尼圆频率与阻尼比的关系见图1。
在应用中，ωm通常取多自由度体系的基频，ωn在对结构振动有显著贡献的高阶振型中选取。
这样，具有高频率的振型反应将因其高阻尼比而被有效地消除。
[3]2 单自由度黏滞阻尼体系简谐振动单自由度黏滞阻尼体系简谐振动见图2。
根据结构动力学理论[2]，单自由度黏滞阻尼体系在简谐载荷下的运动方程为由此可以完全确定体系在简谐载荷下的振动情况。
3 有阻尼简支匀质梁弯曲简谐振动3.1 有阻尼简支匀质梁的弯曲振动微分方程简支匀质梁及其阻尼组成示意见图3。
假定该简支匀质梁内的阻尼由与绝对速度有关的外阻尼力fD和材料应变的黏滞阻尼力σD 这2种成分组成。
当梁各点以横向位移u（x，t）弯曲振动时，设黏滞阻尼系数为C（x），则外阻尼力可以使得匀质梁各阶振型相互独立，并且在广义坐标（将在第3.2节详述）下各阶阻尼如式（1）形式。
Abaqus中一种考虑材料阻尼的随机响应分析方法
Abaqus中一种考虑材料阻尼的随机响应分析方法作者：邓长喜来源：《科技视界》2016年第19期[摘要]介绍了Abaqus中随机响应分析的特点，简要论述了虚拟激励法的理论基础，介绍了利用Abaqus谐响应分析结果得到随机响应结果的步骤。
并以悬臂梁为例，将本文方法和Abaqus中自带随机振动求解器计算结果进行对比。
结果表明，两者误差很小，本文方法计算结果可信。
[关键词]随机响应；振动；虚拟激励法；大阻尼0引言在工程上，随机响应分析常用来预先分析设备抵抗随机载荷的能力。
随机振动将概率论与统计学中的方法应用到评估设备结构的稳定性、识别、响应以及可靠性上，形成了一个相关学科相结合的新产物。
Abaqus作为一个通用有限元求解器，在各个行业得到了广泛的应用。
随机分析的激励作为时间的函数。
具体载荷（力、速度等）都是未知的，载荷的本质是通过统计学的方法描述。
在Abaqus中进行随机响应分析时，直接输入激励的功率谱密度（PSD）曲线，然后有求解器直接计算出结构响应的功率谱密度曲线和对应的均方根（RMS）值。
随机响应分析是一个频域的线性摄动分析过程，通过结构的特征模态来计算得到。
阻尼可以通过模态阻尼、结构阻尼、瑞利阻尼或者复合阻尼等方式施加到结构模型中。
直接随机响应分析具有方便操作，计算速度快等特点。
但是，直接随机响应基于模态叠加法，使用的是全局的模态阻尼，只适用于小阻尼的情况，不能形成阻尼矩阵参与动力学分析，无法为材料添加阻尼，在含有大阻尼的橡胶材料构成的减震器模型中无法应用该分析。
本文简要介绍了虚拟激励法，并通过该方法处理Abaqus稳态动力学分析结果得到对应PSD激励下的随机响应结果，并且和常规随机响应分析结果做出对比。
1基本原理稳态动力学和随机振动分析在Abaqus中都是基于扫频分析，本质上是一样的。
但是，稳态动力学（正弦振动）分析的输入（激励）是各个频率点的加速度或位移峰值，响应为各个频率点的加速度或位移峰值。
abaqus 总体刚度阻尼质量矩阵提取python程序
abaqus 总体刚度阻尼质量矩阵提取python程序摘要：一、引言1.介绍Abaqus 软件在工程领域中的应用2.阐述总体刚度、阻尼、质量矩阵在Abaqus 中的重要性3.介绍Python 编程语言及其在Abaqus 脚本开发中的应用二、Abaqus 总体刚度阻尼质量矩阵提取Python 程序设计1.程序设计目标2.程序结构及主要功能模块3.所需Python 库及其功能介绍三、程序实现1.导入所需库及模块2.定义Abaqus 模型及总体刚度阻尼质量矩阵提取函数3.实现模型导入及参数设置4.提取并输出总体刚度阻尼质量矩阵四、程序测试与验证1.测试数据准备2.程序运行及结果分析3.结果与Abaqus 软件内置结果对比，验证程序准确性五、结论1.总结程序的主要功能及应用场景2.指出程序可能的改进方向3.对未来Abaqus 与Python 结合的技术发展进行展望正文：一、引言Abaqus 是一款广泛应用于工程领域的有限元分析软件，可以对结构、材料、热传导等多种物理现象进行建模与分析。
在Abaqus 中，总体刚度、阻尼、质量矩阵是描述模型动力特性的重要参数，对于分析模型的振动响应、结构性能等具有关键意义。
随着Python 编程语言的普及，越来越多的工程师开始使用Python 开发Abaqus 脚本，以提高工作效率和实现复杂功能的自动化。
本文将介绍一种用于提取Abaqus 模型总体刚度阻尼质量矩阵的Python 程序。
二、Abaqus 总体刚度阻尼质量矩阵提取Python 程序设计本程序旨在通过Python 编程语言实现对Abaqus 模型总体刚度、阻尼、质量矩阵的提取。
程序设计时考虑了易用性和通用性，可以适应不同类型和规模的Abaqus 模型。
程序主要包括以下功能模块：导入所需库及模块、定义Abaqus 模型及总体刚度阻尼质量矩阵提取函数、实现模型导入及参数设置、提取并输出总体刚度阻尼质量矩阵。
三、程序实现首先，程序通过import 语句导入所需的Python 库，如abqpy、numpy 等。
ABAQUS_关键字详解
ABAQUS 关键字详解2010-06-03 19:47:14 作者：huright*Boundary：指定边界条件(用来在节点定义边界条件或在子模型分析中指定被驱动的节点。
可选参数：amplitude：该参数仅在一些预设的变量有非零大小时使用。
设置该参数等于amplitude曲线名。
如果在standard中忽略该参数，则是线性ramp或是阶越型step。
位移只能是ramp型，而移动速度和转动速度只能是step型。
如果在explicit里忽略该参数，则参考的数量会在step开始时立刻应用，并保持常数。
在standard动态或模态分析中，应用与位移或速度的振幅曲线会被自动光滑处理。
而在explicit动态分析中，用户必须请求平滑处理才可以。
load case：该参数只用于standard分析，它只在直接法稳态动力学和屈曲分析中使用，在这两个过程中，该参数可以设置等于1(默认)或2。
如果用于直接法稳态动力学中，load case=1定义边界条件的实部，而load case=2定义了虚部。
如果用于屈曲分析，load case=1为应用载荷定义边界条件，而load case=2用来为屈曲模态定义反对称边界条件op：设置op=mod(默认)更改已存边界条件或为以前未被约束的自由度添加边界条件 op=new则如果所有当前起作用的边界条件都被移除，为了移除边界条件，使用op=new并重新指定所有要被处理的边界条件。
如果在standard的应力/位移分析中边界条件被移除，他们会被与在前一个step中计算产生的反力相等的集中力代替，如果该step是通用非线性分析步，则集中力会根据*step中的amplitude 参数来移除。
因此，默认幅值被使用，而集中力将在该静态分析step结束后被线性减少到零，然后立刻到动态分析。
type：用于应力/位移分析指定数值是位移历程形式、速度历程形式还是加速度历程形式。
在standard中，type=velocity是指定有限转动。
ABAQUS阻尼的确定
abaqus中各种阻尼的定义阻尼在结构计算中，本身就是一个很复杂的问题。
为了应对不同的问题，ABAQUS通过不同的命令来满足不同形式的阻尼。
首先，牵涉到阻尼的命令有：*Damping*Modal damping*Dynamic*Mass*Damping属于材料阻尼，该命令在材料模块*materia l中定义。
其参数Alph a和Beta的定义用于直接积分法(模型分析法时失效)。
而参数comp osite与*Modal damping一起用于模型分析法中复合阻尼的定义。
*Modal damping用于模型分析法。
可以定义直接阻尼，瑞雷阻尼和结构阻尼。
这两个命令都可以用参数AL PHA 和Beta来定义瑞雷阻尼。
但需要注意的是，如果两者同时使用，也就是在进行模型分析时*Damping中的Alph a和Beta值将被忽略。
而*Dynamic中参数Alp ha定义的是人工阻尼，它只是便于计算的一个算子。
与所说的材料阻尼不是一个性质。
当其值为零时就是所谓的纽马克-贝塔法。
至于*Mass中的参数Alpha定义的自然是瑞雷阻尼中考虑质量那一部分的影响因素。
当不考虑刚度影响，即Beta值的影响时，就可以直接在*Mass中定义质量的影响因子。
ABAQUS-显式非线性动态分析-Rayleig h阻尼1.显隐分析两个程序最大的不同在于求解节点加速度的方式上。
隐式分析在完全Newton法迭代求解方法的基础上，使用自动增量步。
对于光滑的非线性响应，以二次速率收敛。
但对于高度的非连续和非线性过程，可能失去二次收敛，并需要大量的迭代过程。
隐式分析的最大分析尺度常常取决于给定计算机中的磁盘空间的大小和可用内存的数量，而不是取决于需要多少计算时间。
显式分析是显式地前推模型状态，所以不需要迭代和收敛准则。
[土木] 在ABAQUS中对框架结构施加地震波(对初学者普及,同时向大虾们求教)
![[土木] 在ABAQUS中对框架结构施加地震波(对初学者普及,同时向大虾们求教)](https://img.taocdn.com/s3/m/9f8fba7ff242336c1eb95eb1.png)
[土木]在ABAQUS中对框架结构施加地震波（对初学者普及，同时向大虾们求教）初学ABAQUS没多久，由于课程需要，想用实体单元建了一个五层的框架结构，要对其施加地震波。
但是我只学了石老师《实例》的前面部分，只知道些比较基本的操作之类的，于是上网求助于论坛。
在找完很多帖子后，建模并计算，基本算是成功的。
所以首先是要衷心感谢各位在论坛上指导了我的楼主及他们的帖子们。
不过感觉网上各位大侠可能都没体会到民间疾苦，只提了个大概，我们自己得总结半天。
所以在此详细点写下自己查到的方法，回报下论坛上的各位，给其他一些初学的人一些帮助，也达到交流学习的目的。
但同时还有许多自己不知其所以然的，想要请教各位大侠（红色字体是引用别人说的，蓝色字体是我的疑问，望大家讨论或帮忙解答）。
在网上查了些方法：module选load，在tools-----amplitude-----creat默认的continue在Edit Amplitude里面输入时间和加速度，点OK。
点creat boundary condition，出现对话框creat boundary condition，选择acceleration/angular acceleration,continue---选择要施加的边界---done----出现对话框edit bondary condition对话框，在amplitude里选择你所定义的时间和加速度。
点ok就完工了。
这是在CAE里输入地震波的方式，我用的方法是直接在inp文件里加地震波的。
首先在CAE里建好模型，定义两个分析步。
第一个分析步是加自重，采用线性加载的方式。
(a)加载方式：ABAQUS在施加Gravity时，默认为Instantaneous（瞬时加载），如果把结构自重以瞬间加载方式加到结构上，相当于对结构施加了一个脉冲荷载，会引起结构在竖向的振动，在不考虑结构阻尼的情况，这种振动会一直持续下去。
abaqus阻尼系数
abaqus阻尼系数摘要：一、引言二、abaqus 软件介绍三、阻尼系数的概念及意义四、abaqus 阻尼系数的设置与调整五、阻尼系数对abaqus 仿真结果的影响六、总结正文：一、引言Abaqus 是一款广泛应用于工程领域的有限元分析软件，通过abaqus 可以进行各种工况下的仿真分析。
阻尼系数是abaqus 中的一个重要参数，对仿真结果具有显著影响。
本文将详细介绍abaqus 阻尼系数的相关知识。
二、abaqus 软件介绍Abaqus 是一款强大的有限元分析软件，由法国公司Dassault Systemes 的SIMULIA 品牌开发。
它具有丰富的功能，可以进行线性和非线性结构分析、热传导分析、热膨胀分析、动力学分析等。
广泛应用于航空航天、汽车制造、建筑结构、生物医学等各个领域。
三、阻尼系数的概念及意义阻尼系数是描述材料阻尼特性的一个参数，它反映了材料在受到振动时，能量消耗的能力。
在有限元分析中，阻尼系数用于模拟材料在振动过程中的能量损耗，对于准确预测结构的动态响应具有重要意义。
四、abaqus 阻尼系数的设置与调整在abaqus 中，阻尼系数的设置可以通过以下步骤完成：1.打开abaqus 软件，创建或打开一个模型。
2.在模型树中，找到需要设置阻尼系数的材料。
3.在材料属性对话框中，找到“阻尼”选项卡。
4.输入阻尼系数的值，并根据需要选择阻尼类型。
5.点击“确定”按钮，完成阻尼系数的设置。
五、阻尼系数对abaqus 仿真结果的影响阻尼系数对abaqus 仿真结果具有显著影响。
不同的阻尼系数设置会导致结构在受到外力作用时的振动幅度、振动衰减时间等特性发生变化。
因此，合理设置阻尼系数对于获得准确的仿真结果至关重要。
六、总结本文对abaqus 阻尼系数的相关知识进行了详细介绍，包括阻尼系数的概念及意义、abaqus 阻尼系数的设置与调整，以及阻尼系数对仿真结果的影响。
ABAQUS瑞利阻尼
关于ABAQUS中的质量比例阻尼总结论：ABAQUS中的质量比例阻尼是和绝对速度有关的，即质量比例阻尼产生的阻尼力由绝对速度引起。
以阻尼系数表达的阻尼，产生的阻尼力由相对速度引起。
Abaqus Analysis User's Manual—Mass proportional damping：The factor introduces damping forces caused by the absolute velocities of the model and so simulates the idea of the model moving through a viscous “ether” (a permeating, still fluid, so that any motion of any point in the model causes damping).帮助手册也说明了质量比例阻尼是和绝对速度有关。
问题：1、应用直接积分法进行时程分析，地震波一般以边界条件的形式加到支座处，结构阻尼只能使用Rayleigh阻尼，而这时产生的阻尼力是绝对速度产生的，而运动方程中的阻尼项产生的阻尼力是与相对速度有关。
2、SAP2000中施加地震波，支座处相对位移为0，绝对位移不为0，其相对位移相对哪一点来说的？算例：单自由度体系，如图（1），质量m=0.02533kg ，k=1N/m ，阻尼比ξ=0.05，对应的阻尼系数c=0.0159，若应用直接积分法进行时程分析，结构的阻尼需要转换成Rayleigh 阻尼，使用如下公式：22n n n βωαξω=+如果只使用质量比例阻尼（结构只有一阶振型），即2n n αξω=，容易得出α=0.6283。
.图（1）情况（1）： 在ABAQUS 中用spring 单元模拟竖向的直杆，水平刚度k=1N/m ，采用Rayleigh 阻尼，通过*mass,alpha=0.6283（质量比例阻尼）施加，地震波需用Elcentrol 波，以边界条件的形式加在支座处(竖向杆下端)。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
关于ABAQUS中的质量比例阻尼总结论：ABAQUS中的质量比例阻尼是和绝对速度有关的，即质量比例阻尼产生的阻尼力由绝对速度引起。
以阻尼系数表达的阻尼，产生的阻尼力由相对速度引起。
Abaqus Analysis User's Manual—Mass proportional damping：The factor introduces damping forces caused by the absolute velocities of the model and so simulates the idea of the model moving through a viscous “ether” (a permeating, still fluid, so that any motion of any point in the model causes damping).帮助手册也说明了质量比例阻尼是和绝对速度有关。
问题：1、应用直接积分法进行时程分析，地震波一般以边界条件的形式加到支座处，结构阻尼只能使用Rayleigh阻尼，而这时产生的阻尼力是绝对速度产生的，而运动方程中的阻尼项产生的阻尼力是与相对速度有关。
2、SAP2000中施加地震波，支座处相对位移为0，绝对位移不为0，其相对位移相对哪一点来说的？算例：单自由度体系，如图（1），质量m=，k=1N/m，阻尼比ξ=，对应的阻尼系数c=，若应用直接积分法进行时程分析，结构的阻尼需要转换成Rayleigh阻尼，使用如下公式：22nnnβωαξω=+如果只使用质量比例阻尼（结构只有一阶振型），即2nnαξω=，容易得出α=。
.图（1）情况（1）：在ABAQUS中用spring单元模拟竖向的直杆，水平刚度k=1N/m，采用Rayleigh阻尼，通过*mass,alpha=（质量比例阻尼）施加，地震波需用Elcentrol波，以边界条件的形式加在支座处(竖向杆下端)。
为了作对比，在SAP2000中的结构阻尼在分析工况中以质量比例阻尼的形式施加。
MATLAB中变成使用NewMark-beta方法。
三者质量的相对位移时程对比如下：图（2）质量点相对位移时程对比图（3）质量点绝对位移时程对比图（4）支座位移时程对比结论：质量点相对位移时程SAP2000与MATLAB重合很好，ABAQUS与两者差别较大；SAP2000与ABAQUS在支座处得位移时程重合很好。
情况（2）：在ABAQUS中用spring单元模拟竖向的直杆，水平刚度k=1N/m，采用Rayleigh阻尼，通过*mass,alpha=（质量比例阻尼）施加，地震波需用Elcentrol波，以惯性力的形式加质量点处。
SAP2000的参数设置同情况（1）。
三者质量的相对位移时程对比如下：图（5）质量点相对位移时程对比图（6）质量点绝对位移时程对比图（7）支座位移时程对比结论：质量点相对位移时程SAP2000与MATLAB及ABAQUS 重合很好，几乎完全一致。
情况（3）：在ABAQUS中用spring单元模拟竖向的直杆，水平刚度k=1N/m，采用阻尼器模拟结构阻尼，通过*dashpot施加阻尼系数c=，地震波需用Elcentrol波，以边界条件的形式加在支座处(竖向杆下端)。
SAP2000的参数设置：结构阻尼在连接单元的属性中施加阻尼阻，尼系数c=。
图（8）质量点相对位移时程对比图（9）质量点绝对位移时程对比图（10）支座位移时程对比结论：位移时程SAP2000与MATLAB及ABAQUS重合很好，几乎完全一致。
情况（4）：在ABAQUS中用spring单元模拟竖向的直杆，水平刚度k=1N/m，采用阻尼器模拟结构阻尼，通过*dashpot施加阻尼系数c=，地震波需用Elcentrol波，以惯性力的形式加质量点处。
SAP2000的参数设置：结构阻尼在连接单元的属性中施加阻尼阻，尼系数c=。
图（11）质量点相对位移时程对比图（12）质量点绝对位移时程对比图（13）支座位移时程对比结论：质量点相对位移时程SAP2000与MATLAB及ABAQUS 重合很好，几乎完全一致。
情况（5）：含阻尼器结构体系不变，结构阻尼采用质量比例阻尼，alpha=。
在质量点水平方向加入阻尼器，采用Maxwell模型，在ABAQUS中用spring 单元和dashpot单元模拟，如图（14）。
阻尼器参数为c d=，弹簧刚度k d =1。
地震波需用Elcentrol波，以惯性力的形式加质量点处。
图（14）在SAP2000中的结构阻尼在分析工况中以质量比例阻尼的形式施加。
图（15）质量点相对位移时程对比图（16）阻尼力滞回曲线对比结论：SAP2000与ABAQUS重合很好，几乎完全一致，MATLAB 与前两者略有差别。
情况（6）：含阻尼器结构体系不变，结构阻尼采用质量比例阻尼，alpha=。
在质量点水平方向加入阻尼器，采用Maxwell模型，在ABAQUS中用spring 单元和dashpot单元模拟，如图（14）。
阻尼器参数为c d=，弹簧刚度k d =1。
地震波需用Elcentrol波，以边界条件的形式加在支座处(竖向杆下端)，另外，水平阻尼器右端的支座固定。
图（17）质量点相对位移时程对比图（18）阻尼力滞回曲线对比结论：ABAQUS与其他两者差别较大。
情况（7）：含阻尼器结构体系不变，结构阻尼采用质量比例阻尼，alpha=。
在质量点水平方向加入阻尼器，采用Maxwell模型，在ABAQUS中用spring 单元和dashpot单元模拟，如图（14）。
阻尼器参数为c d=，弹簧刚度k d =1。
地震波需用Elcentrol波，以边界条件的形式加在支座处(竖向杆下端)，另外，水平阻尼器右端的支座也加地震波。
结论：三者有差别，没有情况（5）吻合的好ABAQUS响应的INPUT文件如下：情况（1）*Heading不含阻尼器，采用质量比例阻尼，地震波以边界条件施加*Node1,0,02,0,1*Nset,Nset=Nout1,2*Element,Type=Mass,Elset=PointMass1,2*Mass,Elset=PointMass,alpha=*Element,Type=SPRING2,Elset=spring2,1,2*SPRING,ELSET=spring1,11*Boundary1,1,6**Step1:Gravity*STEP,NAME=GravitySTEP1:Gravity*Static,1,1e-5,1*Dload,GRAV,0,*Output,FIELD, VARIABLE=PRESELECT*END STEP**Step2:Modal*Step,Name=ModalAnalysis,PERTURBATIONStep2:Modal anlysis*Frequency3*Output,Field,Variable=Preselect*End Step**Step3:Earthquake*Step,Name=TimeHistory,Inc=2000,NLGEOMStep3:Earthquake*Dynamic,HAFTOL=,ALPHA=0,30,0,*Amplitude,Name=Earthquake,Input=*Boundary,op=Mod,Type=Acceleration,Amplitude=Earthquake 1,1,1,1E-2*Output, field, variable=PRESELECT*Output, history,FREQUENCY=1*Node Output,Nset=NoutU1,U2,RF*Element Output,Elset=springS11,E11*End Step情况（2）*Heading不含阻尼器，采用质量比例阻尼，地震波以惯性力形式施加*Node1,0,02,0,1*Nset,Nset=Nout1,2*Element,Type=Mass,Elset=PointMass1,2*Mass,Elset=PointMass,alpha=*Element,Type=SPRING2,Elset=spring2,1,2*SPRING,ELSET=spring1,11*Boundary1,1,6**Step1:Gravity*STEP,NAME=GravitySTEP1:Gravity*Static,1,1e-5,1*Dload,GRAV,0,*Output,FIELD, VARIABLE=PRESELECT*END STEP**Step2:Modal*Step,Name=ModalAnalysis,PERTURBATIONStep2:Modal anlysis*Frequency3*Output,Field,Variable=Preselect*End Step**Step3:Earthquake*Step,Name=TimeHistory,Inc=2000,NLGEOMStep3:Earthquake*Dynamic,HAFTOL=,ALPHA=0,30,0,*Amplitude,Name=Earthquake,Input=*Cload,Amplitude=Earthquake2,1,*Output, field, variable=PRESELECT*Output, history,FREQUENCY=1*Node Output,Nset=NoutU1,U2,RF*Element Output,Elset=springS11,E11*End Step情况（3）*Heading不含阻尼器，采用阻尼器模拟结构阻尼，地震波以边界条件形式施加*Node1,0,02,0,1*Nset,Nset=Nout1,2*Element,Type=Mass,Elset=PointMass1,2*Mass,Elset=PointMass*Element,Type=SPRING2,Elset=spring2,1,2*SPRING,ELSET=spring1,11*Element,Type=DASHPOT2,Elset=dashpot3,1,2*DASHPOT, ELSET=dashpot1,1*Boundary1,1,6**Step1:Gravity*STEP,NAME=GravitySTEP1:Gravity*Static,1,1e-5,1*Dload,GRAV,0,*Output,FIELD, VARIABLE=PRESELECT*END STEP**Step2:Modal*Step,Name=ModalAnalysis,PERTURBATIONStep2:Modal anlysis*Frequency3*Output,Field,Variable=Preselect*End Step**Step3:Earthquake*Step,Name=TimeHistory,Inc=2000,NLGEOMStep3:Earthquake*Dynamic,HAFTOL=,ALPHA=0,30,0,*Amplitude,Name=Earthquake,Input=*Boundary,op=Mod,Type=Acceleration,Amplitude=Earthquake1,1,1,1E-2*Output, field, variable=PRESELECT*Output, history,FREQUENCY=1*Node Output,Nset=NoutU1,U2,RF*Element Output,Elset=springS11,E11*End Step情况（4）*Heading不含阻尼器，采用阻尼器模拟结构阻尼，地震波以惯性力形式施加*Node1,0,02,0,1*Nset,Nset=Nout1,2*Element,Type=Mass,Elset=PointMass1,2*Mass,Elset=PointMass*Element,Type=SPRING2,Elset=spring2,1,2*SPRING,ELSET=spring1,11*Element,Type=DASHPOT2,Elset=dashpot3,1,2*DASHPOT, ELSET=dashpot1,1*Boundary1,1,6**Step1:Gravity*STEP,NAME=GravitySTEP1:Gravity*Static,1,1e-5,1*Dload,GRAV,0,*Output,FIELD, VARIABLE=PRESELECT*END STEP**Step2:Modal*Step,Name=ModalAnalysis,PERTURBATIONStep2:Modal anlysis*Frequency3*Output,Field,Variable=Preselect*End Step**Step3:Earthquake*Step,Name=TimeHistory,Inc=2000,NLGEOMStep3:Earthquake*Dynamic,HAFTOL=,ALPHA=0,30,0,*Amplitude,Name=Earthquake,Input=*Cload,Amplitude=Earthquake2,1,*Output, field, variable=PRESELECT*Output, history,FREQUENCY=1*Node Output,Nset=NoutU1,U2,RF*Element Output,Elset=springS11,E11*End Step情况（5）*HeadingMaxwell模型阻尼器模拟,地震力以惯性力形式施加,采用质量比例阻尼*Node1,0,02,0,13,,14,1,1*Nset,Nset=Nout1,2,4*Element,Type=Mass,Elset=PointMass1,2*Mass,Elset=PointMass,alpha=*Element,Type=SPRING2,Elset=spring_12,1,2*SPRING,ELSET=spring_11,11**定义阻尼器的属性，用Spring和Dashpot单元*Element,Type=SPRINGA,Elset=spring_23,3,4*Spring,Elset=spring_21*Element,Type=DASHPOTA,Elset=dashpot 4,2,3*DASHPOT, ELSET=dashpot*Boundary1,1,64,1,6**Step1:Gravity*STEP,NAME=GravitySTEP1:Gravity*Static,1,1e-5,1*Dload,GRAV,0,*Output,FIELD, VARIABLE=PRESELECT*END STEP*Step,Name=ModalAnalysis,PERTURBATION Step2:Modal anlysis*Frequency3*Output,Field,Variable=Preselect*End Step**Step3:Earthquake*Step,Name=TimeHistory,Inc=2000,NLGEOM Step3:Earthquake*Dynamic,HAFTOL=,ALPHA=0,30,0,*Amplitude,Name=Earthquake,Input=*Cload,Amplitude=Earthquake2,1,*Output, field, variable=PRESELECT*Output, history,FREQUENCY=1*Node Output,Nset=NoutU1,U2,RF*Element Output,Elset=spring_1S11,E11*Element Output,Elset=spring_2S11,E11*Element Output,Elset=dashpotS11,E11,ER11*End Step情况（6）*HeadingMaxwell模型阻尼器模拟,地震力以边界形式施加,采用质量比例阻尼*Node1,0,02,0,13,,14,1,1*Nset,Nset=Nout1,2,4*Element,Type=Mass,Elset=PointMass1,2*Mass,Elset=PointMass,alpha=*Element,Type=SPRING2,Elset=spring_12,1,2*SPRING,ELSET=spring_11,11**定义阻尼器的属性，用Spring和Dashpot单元*Element,Type=SPRINGA,Elset=spring_23,3,4*Spring,Elset=spring_21*Element,Type=DASHPOTA,Elset=dashpot4,2,3*DASHPOT, ELSET=dashpot*Boundary1,1,64,1,6**Step1:Gravity*STEP,NAME=GravitySTEP1:Gravity*Static,1,1e-5,1*Dload,GRAV,0,*Output,FIELD, VARIABLE=PRESELECT*END STEP*Step,Name=ModalAnalysis,PERTURBATIONStep2:Modal anlysis*Frequency3*Output,Field,Variable=Preselect*End Step**Step3:Earthquake*Step,Name=TimeHistory,Inc=2000,NLGEOMStep3:Earthquake*Dynamic,HAFTOL=,ALPHA=0,30,0,*Amplitude,Name=Earthquake,Input=*Boundary,op=Mod,Type=Acceleration,Amplitude=Earthquake1,1,1,1E-2*Output, field, variable=PRESELECT*Output, history,FREQUENCY=1*Node Output,Nset=NoutU1,U2,RF*Element Output,Elset=spring_1S11,E11*Element Output,Elset=spring_2S11,E11*Element Output,Elset=dashpotS11,E11,ER11*End Step情况（7）*HeadingMaxwell模型阻尼器模拟,地震力以边界形式施加,采用质量比例阻尼,右支座也加地震波*Node1,0,02,0,13,,14,1,1*Nset,Nset=Nout1,2,4*Element,Type=Mass,Elset=PointMass1,2*Mass,Elset=PointMass,alpha=*Element,Type=SPRING2,Elset=spring_12,1,2*SPRING,ELSET=spring_11,11**定义阻尼器的属性，用Spring和Dashpot单元*Element,Type=SPRINGA,Elset=spring_23,3,4*Spring,Elset=spring_21*Element,Type=DASHPOTA,Elset=dashpot 4,2,3*DASHPOT, ELSET=dashpot*Boundary1,1,64,1,6**Step1:Gravity*STEP,NAME=GravitySTEP1:Gravity*Static,1,1e-5,1*Dload,GRAV,0,*Output,FIELD, VARIABLE=PRESELECT*END STEP*Step,Name=ModalAnalysis,PERTURBATION Step2:Modal anlysis*Frequency3*Output,Field,Variable=Preselect*End Step**Step3:Earthquake*Step,Name=TimeHistory,Inc=2000,NLGEOMStep3:Earthquake*Dynamic,HAFTOL=,ALPHA=0,30,0,*Amplitude,Name=Earthquake,Input=*Boundary,op=Mod,Type=Acceleration,Amplitude=Earthquake 1,1,1,1E-24,1,1,1E-2*Output, field, variable=PRESELECT*Output, history,FREQUENCY=1*Node Output,Nset=NoutU1,U2,RF*Element Output,Elset=spring_1S11,E11*Element Output,Elset=spring_2S11,E11*Element Output,Elset=dashpotS11,E11,ER11*End Step。