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 波，以边界条件的形式加在支座处(竖向杆下端)。
abaqus材料阻尼比参数含义
【文章】abaqus材料阻尼比参数含义1. 引言在工程领域中，材料的阻尼比是一个非常重要的参数，它对结构的振动特性和动态响应有着重要影响。
在使用abaqus软件进行结构分析和仿真时，对材料的阻尼比参数进行准确理解和设定是至关重要的。
本文将从深度和广度两个方面全面评估abaqus材料阻尼比参数的含义，并根据此进行撰写有价值的文章。
2. abaqus材料阻尼比参数含义的深度探讨abaqus中的材料阻尼比参数指的是在结构振动分析中描述材料内部耗散能力的一个重要参数。
阻尼比越大，表示材料的能量耗散能力越强，振动衰减的速度越快。
在abaqus中，阻尼比参数通常用来描述材料在动态载荷作用下的振动特性，是与材料的内部结构和分子运动状态密切相关的一个参数。
3. abaqus材料阻尼比参数含义的广度探讨阻尼比参数在abaqus中有着广泛的应用，涉及到多个领域和行业。
在工程结构分析中，正确理解和设置材料阻尼比参数可以更准确地预测结构的动态响应和振动特性。
在地震工程中，材料阻尼比参数的设定对结构的抗震性能有着重要的影响。
在航空航天领域，材料阻尼比参数也是飞行器动力学分析和设计中不可或缺的重要参数。
4. 回顾与总结通过深度和广度的探讨，我们对abaqus材料阻尼比参数的含义有了更为全面和深入的理解。
在结构分析和仿真中，正确理解和设置材料阻尼比参数对于准确预测结构的振动特性和动态响应至关重要。
我们也认识到材料阻尼比参数在工程领域中有着广泛的应用，涉及到多个行业和领域的工程实践和研究。
深入理解abaqus材料阻尼比参数的含义对于提高工程设计和分析的准确性和可靠性具有重要意义。
5. 个人观点与理解作为一个工程师，我深知材料阻尼比参数的重要性。
在结构分析和设计中，我经常需要对abaqus中的材料阻尼比参数进行设置和调整。
通过学习和使用，我认识到正确理解并准确设置材料阻尼比参数对于工程实践具有重要意义。
我也意识到材料阻尼比参数的含义和作用并不局限于结构领域，而是涉及到更广泛的工程应用和研究领域。
ABAQUS详细阻尼说明
ABAQUS中阻尼的定义在ABAQUS中阻尼可以应用在下面的动力学分析中：◆非线性问题直接积分求解(显式分析或者隐式分析)；◆直接法或子空间法稳态动力学分析；◆模态动力学分析(线性)。
针对模态动力学分析，在ABAQUS/Standard中可定义几种不同类型的阻尼：直接模态阻尼(Direct Modal Damping)，瑞利阻尼(Rayleigh Damping)，复合模态阻尼(Composite Modal Damping)和结构阻尼(Structure Damping)。
ABAQUS模态动力学分析中用*MODAL DAMPING选项来定义阻尼。
阻尼是包含在分析步内定义的一部分，每阶模态可以定义不同量值的阻尼。
1、直接模态阻尼采用直接模态阻尼可以定义对应于每阶模态的阻尼比ξ。
其典型的取值范围是在临界阻尼的1％～10％之间。
直接模态阻尼允许用户精确定义系统的每阶模态的阻尼。
在分析步骤内定义直接模态阻尼。
如图1所示，激活直接模态阻尼选项(Direct modal)，并在数据行内输入数据。
图1直接模态阻尼定义对应的ABAQUS输入文件为：*MODAL DAMPING,MODAL=DIRECTm1,m2,ξa其中，*MODAL DAMPING选项中的MODAL=DIRECT参数表示被指定的直接模态阻尼，数据行输入的数据m1为起始模态序号，m2为截止模态序号，ξa为模态阻尼比。
例如，对于前10阶振型的阻尼定义为4%的临界模态阻尼，11～20阶振型的阻尼为5%的临界阻尼，在分析步骤中的定义如下：*MODAL DAMPING,MODAL=DIRECT1，10，0.0411，20，0.052、瑞利阻尼在瑞利阻尼中，假设阻尼矩阵可表示为质量矩阵和刚度矩阵的线性组合，即C=αM+βK(1)其中，α和β是用户根据材料特性定义的常数。
尽管假设阻尼正比于质量和刚度没有严格的物理基础，但是实际上我们对于阻尼分布的真实情况知之甚少，也就不能保证其它更为复杂的模型是正确的。
Abaqus地震分析的总结
Abaqus 地震分析的总结——时步、单元尺寸、滤波、等效非线性、无限元 笑看风云
1、自由场地震反应
经典的自由场地震反应用 shake91 或 proshake 等进行分析，在分析完可以导出各层土的等效线性参数，包 括阻尼（粘滞阻尼系数）和剪切模量，用剪切模量可以计算弹性模量，shake 中假定泊松比为常数，对地震 反应没影响。其实 FLAC 中有自带自由场边界，计算地震很方便。 如以下是 shake91 中自带的例子输出得到的参数 ITERATION NUMBER 8 VALUES IN TIME DOMAIN NO TYPE DEPTH UNIFRM. <---- DAMPING ----> <---- SHEAR MODULUS -----> G/Go (FT) STRAIN NEW USED ERROR NEW USED ERROR RATIO --- ---- ---- ------- ----- ------ ------------ ------- ------ ----1 2 2.5 .00077 .007 .007 .0 3851.5 3851.5 .0 .992 2 2 7.5 .00295 .014 .014 .0 3020.0 3020.0 .0 .960 3 2 15.0 .00634 .023 .023 .0 2803.8 2803.8 .0 .892 4 2 25.0 .00976 .028 .028 .0 2985.8 2985.8 .0 .852 5 1 35.0 .01099 .030 .030 .0 3621.7 3621.6 .0 .933 6 1 45.0 .01403 .035 .035 .0 3540.5 3540.4 .0 .912 7 1 55.0 .01362 .034 .034 .0 4296.0 4296.0 .0 .915 8 1 65.0 .01566 .037 .037 .0 4239.8 4239.8 .0 .903 9 2 75.0 .01356 .034 .034 .0 5402.7 5402.7 .0 .792 10 2 85.0 .01505 .037 .037 .0 5266.0 5266.0 .0 .772 11 2 95.0 .01336 .034 .034 .0 6288.2 6288.2 .0 .795 12 2 105.0 .01413 .035 .035 .0 6203.4 6203.4 .0 .784 13 2 115.0 .01233 .032 .032 .0 7357.2 7357.2 .0 .810 14 2 125.0 .01281 .033 .033 .0 7290.8 7290.8 .0 .803 15 2 135.0 .01115 .030 .030 .0 8570.4 8570.4 .0 .829 16 2 145.0 .00865 .026 .026 .0 11292.6 11292.6 .0 .863 Shake 中的 outcrop 指出露基岩，baserock 指土层底部的基岩，因此不考虑波的衰减的情况下，在 outcrop 处用加速度计测得的地震加速度幅值为 baserock 处的 2 倍。 如果在 abaqus 中的土层与基岩的交界面处输入 地震波，在跟 shake91 对比时，要是用 baserock 输入。在 abaqus 中建立的土层数要和 shake91 中的土层数 相等，使用同样的参数，注意是弹模，泊松比，密度和瑞利阻尼。注意 shake91 中必须是 8 列或 4 列数据， 如果使用两列数据，则无法读入加速度。Shake91 是频域有限元法，对于加速度幅值是先放大，再滤波。 自带的例子中滤掉大于 25Hz 的波。
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各材料瑞利阻尼参数瑞利阻尼参数是描述材料在振动过程中吸收和耗散能量的重要参数之一。
在Abaqus中，不同材料的瑞利阻尼参数可以通过材料定义和模型参数来进行设置。
本文将介绍Abaqus中常用材料的瑞利阻尼参数设置方法，并对其物理意义进行解释。
1. 线性弹性材料的瑞利阻尼参数（RAYLEIGH DAMPING）对于线性弹性材料，其瑞利阻尼参数可以通过RAYLEIGH DAMPING定义。
该参数是通过将材料的质量阻尼系数和刚度阻尼系数相加得到的。
质量阻尼系数描述了材料在振动过程中由于惯性而产生的能量耗散，刚度阻尼系数描述了材料在振动过程中由于刚度而产生的能量耗散。
在Abaqus中，可以通过在材料定义中设置RAYLEIGH DAMPING 参数来指定材料的瑞利阻尼参数。
这些参数的设置需要根据具体的材料特性和振动问题来确定。
2. 金属材料的瑞利阻尼参数（METAL DAMPING）对于金属材料，其瑞利阻尼参数可以通过METAL DAMPING定义。
金属材料的瑞利阻尼参数通常是非线性的，与材料的应变速率和温度有关。
在Abaqus中，可以通过在材料定义中设置METAL DAMPING参数来指定金属材料的瑞利阻尼参数。
这些参数的设置需要根据具体的金属材料特性和振动问题来确定。
3. 土壤材料的瑞利阻尼参数（SOIL DAMPING）对于土壤材料，其瑞利阻尼参数可以通过SOIL DAMPING定义。
土壤材料的瑞利阻尼参数通常是非线性的，与材料的应变速率、孔隙水压和应变路径有关。
在Abaqus中，可以通过在材料定义中设置SOIL DAMPING参数来指定土壤材料的瑞利阻尼参数。
这些参数的设置需要根据具体的土壤材料特性和振动问题来确定。
4. 混凝土材料的瑞利阻尼参数（CONCRETE DAMPING）对于混凝土材料，其瑞利阻尼参数可以通过CONCRETE DAMPING定义。
混凝土材料的瑞利阻尼参数通常是非线性的，与材料的应变速率、应力状态和应变路径有关。
abaqus中的动态分析方法
ABAQUS线性动态分析如果你只对结构承受载荷后的长期响应感兴趣，静力分析（StatiC analysis）是足够的。
然而，如果加载时间很短（例如在地震中）或者如果载荷在性质上是动态的（例如来自旋转机械的荷载），你就必须采用动态分析（dynamic analysis ）。
本章将讨论应用ABAQUS/Standard进行线性动态分析；关于应用ABAQUS/Explicit 进行非线性动态分析的讨论，请参阅第9章“非线性显式动态分析”。
7.1引言动态模拟是将惯性力包含在动力学平衡方程中：Mu T-P=O其中M 结构的质量。
U 结构的加速度。
I 在结构中的内力。
P 所施加的外力。
在上面公式中的表述是牛顿第二运动定律（ F = ma ）。
在静态和动态分析之间最主要的区别是在平衡方程中包含了惯性力（MU）。
在两类模拟之间的另一个区别在于内力I的定义。
在静态分析中，内力仅由结构的变形引起；而在动态分析中，内力包括源于运动（例如阻尼）和结构的变形的贡献。
7.1.1固有频率和模态最简单的动态问题是在弹簧上的质量自由振动，如图7-1所示。
图7- 质量一弹簧系统在弹簧中的内力给出为ku ,所以它的动态运动方程为mu ku - P = O这个质量—弹簧系统的固有频率（natral frequency ）（单位是弧度/秒（rad∕s））给出为如果质量块被移动后再释放，它将以这个频率振动。
若以此频率施加一个动态外力，位移的幅度将剧烈增加，这种现象即所谓的共振。
实际结构具有大量的固有频率。
因此在设计结构时，非常重要的是避免使可能的载荷频率过分接近于固有频率。
通过考虑非加载结构（在动平衡方程中令P=O ）的动态响应可以确定固有频率。
则运动方程变为MUT =0对于无阻尼系统，I=KU ,因此有MU Ku =0这个方程的解具有形式为U=Ψe^c t将此式代入运动方程，得到了特征值（eigenvalue）问题K=M其中，=∙2。
该系统具有n个特征值，其中n是在有限元模型中的自由度数目。
abaqus瑞利阻尼[详解]
![abaqus瑞利阻尼[详解]](https://img.taocdn.com/s3/m/1d0abb22f011f18583d049649b6648d7c1c70899.png)
关于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 movingth rough 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 波，以边界条件的形式加在支座处(竖向杆下端)。
阻尼
阻 尼系统结构特征值和模态的求解是在无阻尼情况下得到的，而在动力学问题中，任意结构都应存在或大或小的阻尼，阻尼的大小对系统的动力学响应将产生一定的影响。
本章主要讨论在动力学分析中怎样应用ABAQUS 定义系统的阻尼特性。
引言当系统作无阻尼自由振动时，由于没有能量输入与输出，系统机械能守恒，系统的振幅为常数。
然而在实际结构中，这种无阻尼自由振动并不存在。
结构运动时能量耗散，振幅将逐渐减小直至停止振动，这种能量耗散被称为阻尼(damping)。
通常假定阻尼为粘性的，其大小正比于速度，方向与速度相反。
有阻尼结构系统的动力学方程可以写为：0MuI P +−=&& (4.1.1) I Ku Cu=+& (4.1.2) 其中， C 为结构的阻尼矩阵，u&为结构的速度。
能量耗散来源于几个因素，其中包括结构连接处的摩擦和局部材料的迟滞效应。
阻尼对于表征结构吸收能量是一个很方便的方法，它包含了重要的能量吸收过程，而不需要模拟耗能的具体机制。
阻尼在ABAQUS/Standard 中，特征模态的计算是从无阻尼系统中提取出的。
然而，大多数工程问题都包含某种阻尼，尽管阻尼可能很小。
对于每一模态，有阻尼固有频率和无阻尼固有频率之间的关系是：d ω= (4.2.1)其中d ω是有阻尼的固有频率；cc c ξ=为临界阻尼；c 是该模态的阻尼，0c 是该模态的临界阻尼。
当临界阻尼ξ取较小值(< 0.1ξ)时，有阻尼系统的特征频率和特征向量与无阻尼系统非常接近；随着ξ的增加，采用无阻尼系统求得的特征频率就会变得不太准确，当ξ接近1时，无阻尼特征频率和特征向量就失效了。
但是，大多数用线性动力学分析的结构问题只有很小的阻尼，因而可以采用无阻尼特征频率。
当结构处于临界阻尼即1ξ=时，施加一个扰动后，结构不会振荡，而是尽可能迅速地恢复到它的初始静止构形，如图4.2.1所示。
图4.2.1 阻尼在ABAQUS中定义阻尼在ABAQUS中阻尼可以应用在下面的动力学分析中：z非线性问题直接积分求解(显式分析或者隐式分析)；z直接法或子空间法稳态动力学分析；z模态动力学分析(线性)。
abaqus计算不收敛阻尼设置
abaqus计算不收敛阻尼设置在使用ABAQUS进行数值模拟时，经常会遇到计算不收敛的问题。
不收敛意味着模拟过程中无法达到稳定状态，计算结果难以获得或不准确。
为了解决这个问题，我们可以通过设置阻尼来提高模拟的稳定性。
阻尼是指在模拟过程中引入的一种能量损耗机制，通过消耗一部分能量来减小模拟系统的振荡或不稳定性。
在ABAQUS中，可以通过设置各种类型的阻尼来调节模拟的收敛性。
在ABAQUS中，有多种类型的阻尼可供选择。
其中，常用的有模型阻尼、材料阻尼和数值阻尼等。
模型阻尼是通过改变模型的刚度、质量或阻尼系数来实现的。
材料阻尼是通过改变材料的物理性质或材料的损伤模型来实现的。
数值阻尼是通过改变数值求解器的参数或算法来实现的。
在设置阻尼之前，我们首先要确定模拟中存在的不收敛问题的原因。
不收敛问题可能由多种因素引起，如模型初始状态不合理、材料参数不准确、加载条件不合适、网格划分不合理等。
只有明确了问题的根源，才能有针对性地设置阻尼来解决不收敛问题。
针对不同的不收敛问题，我们可以采取不同的阻尼设置策略。
例如，对于模型初始状态不合理导致的不收敛问题，可以通过增加模型的初始速度或位移来提高收敛性。
对于材料参数不准确导致的不收敛问题，可以通过调整材料的阻尼系数或损伤模型来改善收敛性。
对于加载条件不合适导致的不收敛问题，可以通过调整加载的速度或幅值来增加模拟的稳定性。
对于网格划分不合理导致的不收敛问题，可以通过细化网格或采用更合适的网格类型来提高收敛性。
除了根据不同的不收敛问题采取不同的阻尼设置策略外，还可以通过调整数值求解器的参数来改善收敛性。
ABAQUS提供了丰富的数值求解器参数可供选择，如时间步长、收敛准则、迭代次数等。
通过调整这些参数，可以提高模拟的稳定性和收敛性。
在设置阻尼时，需要注意一些问题。
首先，阻尼设置不能过大或过小，过大会导致模拟结果不准确，过小会导致模拟不收敛。
其次，阻尼设置应该与模拟系统的实际情况相匹配，不能随意设置。
瑞利阻尼计算公式
瑞利阻尼计算公式
瑞利阻尼计算公式
瑞利阻尼计算公式是一种根据物体的动能和惯性计算它在特定情况下的阻尼的方法。
瑞利阻尼公式可以应用于不同的动态系统，如悬挂系统、离合器、减速器、液压转子等。
它可以用来测量物体动能和惯性之间的关系，以确定物体在特定情况下的阻尼。
瑞利阻尼计算公式的基本原理是，当物体的惯性被减少时，它的动能也会减少。
因此，用瑞利阻尼公式计算物体的阻尼，是根据物体的惯性和动能之间的关系来进行的。
它的公式如下：
F = m * a
其中，F表示物体的阻尼，m表示物体的质量，a表示物体的加速度。
瑞利阻尼计算公式可以用来计算物体在特定情况下的阻尼。
它可以用来测量物体动能和惯性之间的关系，以确定物体在特定情况下的阻尼。
它可以用于不同的动态系统，如悬挂系统、离合器、减速器、液压转子等。
瑞利阻尼计算公式是一种有效的方法，可以用来测量物体惯性和动能之间的关系，以确定物体在特定情况下的阻尼。
它可以用于不同的动态系统，是设计精确动态系统的重要方法。
abaqus系列之关键字编辑
总规则1、关键字必须以*号开头，且关键字前无空格2、**为注释行，它可以出现在文件中的任何地方3、当关键字后带有参数时，关键词后必须采用逗号隔开4、参数间都采用逗号隔开5、关键词可以采用简写的方式，只要程序能识别就可以了6、不需使用隔行符，如果参数比较多，一行放不下，可以另起一行，只要在上一行的末尾加逗号便可以*AMPLITUDE：定义幅值曲线这个选项允许任意的载荷、位移和其它指定变量的数值在一个分析步中随时间的变化(或者在ABAQUS/Standard分析中随着频率的变化)。
必需的参数：NAME：设置幅值曲线的名字可选参数：DEFINITION：设置definition=Tabular(默认)给出表格形式的幅值-时间(或幅值-频率)定义。
设置DEFINITION=EQUALLY SPACED/PERIODIC/MODULATED/DECAY/SMOOTH STEP/SOLUTION DEPENDENT或BUBBLE来定义其他形式的幅值曲线。
INPUT：设置该参数等于替换输入文件名字。
TIME：设置TIME=STEP TIME(默认)则表示分析步时间或频率。
TIME=TOTAL TIME表示总时间。
VALUE：设置VALUE=RELATIVE(默认)，定义相对幅值。
VALUE=ABSOLUTE表示绝对幅值，此时，数据行中载荷选项内的值将被省略，而且当温度是指定给已定义了温度TEMPERATURE=GRADIENTS(默认)梁上或壳单元上的节点，不能使用ABSOLUTE。
对于DEFINITION=TABULAR的可选参数：SMOOTH：设置该参数等于DEFINITION=TABULAR的数据行第一行1、时间或频率2、第一点的幅值(绝对或相对)3、时间或频率4、第二点的幅值(绝对或相对) 等等基本形式：*Amplitude，name=Amp-10.，0.，0.2，1.5，0.4，2.，1.，1.*BEAM SECTION：当需要数值积分时定义梁截面*BOND：定义绑定和绑定属性*BOUNDARY：定义边界条件用来在节点定义边界条件或在子模型分析中指定被驱动的节点。
- 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。