abaqus第九章-动力问题

《abaqus第九章-动力问题》由会员分享，可在线阅读，更多相关《abaqus第九章-动力问题（21页珍藏版）》请在装配图网上搜索。

1、第九章 动力问题如果只对结构加载荷后的长期响应感兴趣的话，静力分析就足够了。然而，如果加载时间很短，例如地震；或者载荷性质为动态，例如来自旋转机械的荷载，这时就必须采用动力分析。 9.1 引言 动态模拟是将惯性力包含在动力学平衡方程中： 其中M是结构的质量。 是结构的加速度。 是结构中的内力。 P是所施加的外力。公式的表述无非是牛顿的第二运动定律（F=ma）的表现。 动态分析和静态分析最主要的不同在于平衡方程中包含惯性力项（M）。两者的另一个不同之处在于内力的定义。在静态分析中，内力仅由结构的变形引起；而动态分析中的内力包括运动（例如阻尼）和结构变形的共同贡献。9.1.1 固有频率和模态 最简

2、单的动力问题是在弹簧上的质量振动，如图9-1所示。图91质量弹簧系统 弹簧的内力为，所以运动方程为这个质量弹簧系统的固有频率（单位是弧度/秒）为如果质量块被移动后再释放，它将以这个频率振动。假若以此频率施加一个动态外力，位移的幅度将剧烈增加即所谓的共振现象。实际的结构具有多个固有频率。因此，在设计结构时避免使各固有频率与可能的荷载频率过分接近就很重要。固有频率可以通过分析结构在无荷载（动力平衡方程中的）时的动态响应而得到。此时，运动方程变为对于无阻尼系统，则上式变为这个方程解的形式为将此式代入到运动方程中便得到了特征值问题方程其中。 该系统具有个特征值，此处是有限元模型的自由度数。记为第个特征

3、值。它的平方根是结构的第阶固有频率，并且是相应的第阶特征向量。特征向量也就是所谓的模态（也称为振型），因为它是结构在第阶振型下的变形状态。 在ABAQUS中，频率提取程序用来求解结构的振型和频率。这个程序使用起来十分简单，只要给出所需振型的数目和所关心的最高频率即可。9.1.2 振型叠加 在线性问题中，结构在荷载作用下的动力响应可以用固有频率和振型来表示，即结构的变形可以采用振型叠加的技术由各振型的组合得到，每一阶模态都要乘以一个标量因子。模型中位移矢量被定义为其中是振型的标量因子。这一技术只在模拟小变形、线弹性材料、无接触条件情况下是有效的，即必须是线性问题。在结构动力学问题中，结构的响应往

4、往取决于相对较少的几阶振型，这使得振型叠加方法在计算这类系统的响应时特别有效。考虑一个含有10,000个自由度的模型，则对运动方程的直接积分需要在每个时间点上求解10,000个联立方程组。但若结构的响应采用100阶振型来描述，那么在每个时间步上只需求解100个方程。更重要的是，振型方程是解耦的，而原来的运动方程是耦合的。虽然在计算振型和频率时需要花费一些时间作为代价，但在计算响应时将节省大量的时间。如果在模拟中存在非线性，在分析中固有频率会发生明显的变化，因此振型叠加法将不再适用。在这种情况下，需要对动力平衡方程直接积分，这将比振型分析要花费更多的时间。具有下列特点的问题才适于进行线性瞬态动力

5、学分析： 系统应该是线性的：线性材料特性，无接触条件，无非线性几何效应。 响应应该只受较少的频率支配。当响应中各频率成分增加时，例如撞击和冲击问题，振型叠加技术的有效性将大大降低。 载荷的主要频率应在所提取的频率范围内，以确保对载荷的描述足够精确。 由于任何突然加载所产生的初始加速度应该能用特征模态精确描述。 系统的阻尼不能过大。9.2 阻尼 如果一个无阻尼结构做自由振动，则它的振幅会保持恒定不变。然而，实际上由于结构运动而能量耗散，振幅将逐渐减小直至振动停止，这种能量耗散称为阻尼。通常假定阻尼为粘滞的或正比于速度。动力平衡方程可以重新写成包含阻尼的形式：其中C是结构的阻尼阵，是结构的速度。能

6、量耗散来自于诸多因素，其中包括结构结合处的摩擦和局部材料的迟滞效应。阻尼概念对于无需顾及能量吸收过程的细节表征而言是一个很方便的方法。 在ABAQUS中，是针对无阻尼系统计算其振型的，然而，大多数工程问题还是包含阻尼的，尽管阻尼可能很小。有阻尼的固有频率和无阻尼的固有频率的关系是其中 是阻尼特征值， 是临界阻尼比， 是该振型的阻尼， 是临界阻尼。 对较小的情形（），有阻尼系统的特征频率非常接近于无阻尼系统的相应值。当增大时，采用无阻尼系统的特征频率就不太准确，当接近于1时，就不能采用无阻尼系统的特征频率了。 当结构处于临界阻尼（）时，施加一个扰动后，结构不会有摆动而是很快地恢复到静止的初始形态

7、。（见图9-2）图92 阻尼9.2.1 ABAQUS中阻尼的定义 在ABAQUS中，为了进行瞬时模态分析，可定义不同类型的阻尼：直接模态阻尼，瑞利(Rayleigh)阻尼和复合模态阻尼。模拟动力学过程要定义阻尼。阻尼是分析步定义的一部分，每阶振型可以定义不同数量的阻尼。直接模态阻尼采用直接模态阻尼可以定义对应于每阶振型的临界阻尼比。的典型取值范围是从1%到10%。直接模态阻尼允许精确定义每阶振型的阻尼。Rayleigh阻尼 Rayleigh阻尼假设阻尼矩阵是质量矩阵和刚度矩阵的线性组合，其中和是用户定义的常数。尽管假设阻尼正比于质量和刚度没有严格的物理基础，但实际上我们对于阻尼的分布知之甚少，

8、也就不能保证使用更为复杂的阻尼模型是正确的。一般来讲，这个模型对于大阻尼系统也就是临界阻尼超过10%时，是失效的。相对于其它形式的阻尼，可以精确地定义系统每阶模态的Rayleigh阻尼。 复合阻尼 在复合阻尼中，可以定义每种材料的临界阻尼比，并且复合阻尼是对应于整体结构的阻尼。当结构中有许多不同种类材料时，这一选项是十分有用的。在这本指南中将不对复合阻尼做进一步的讨论。9.2.2 阻尼值的选择 在大多数线性动力学问题中，恰当地定义阻尼对于获得精确的结果是十分重要的。但是阻尼只是对结构吸收能量这种特性的近似描述，而不是去仿真造成这种效果的物理机制。所以，确定分析中所需要的阻尼数据是很困难的。有时

9、，可以从动力试验中获得这些数据，但是在多数情况下，不得不通过经验或参考资料获得数据。在这些情况下，要仔细地分析计算结果，应该通过参数分析来评价阻尼系数对于模拟的敏感性。9.3 单元选择 事实上，ABAQUS的所有单元均可用于动力分析。选取单元的一般原则与静力分析相同。但是，在模拟冲击和爆炸载荷时，应选用一次单元。因为它们具有集中质量公式，在模拟应力波效果方面优于采用二次单元的一致质量公式。9.4 动力问题的网格剖分 在动力分析中，剖分网格需要考虑响应中将被激发的振型，网格剖分应能充分反映那些振型。这就意味着，能满足静态模拟要求的网格，不一定能计算高频振型的动态响应。 例如，考虑图93所示的平板

10、。用一次壳单元剖分的网格对于受均布载荷的静力分析以及一阶振型的预测是适用的。但是，该网格对于精确模拟第六阶振型就显得太粗糙了。图93 平板的粗网格图94显示了同样的板用一次单元进行更精细的网格剖分后所模拟的结果。第六阶振型的位移看起来明显好一些，预测的频率也准确了许多。如果作用在板上的动载会显著地激发该阶振型，则应采用精细的网格；采用粗网格将得不到精确的结果。图 94 平板的精细网格9.5 例子：货物起重机动态载荷这个例子采用在第6.4节中分析过的同样的结构，现在要研究的问题是当10kN荷重跌落在起重机挂钩上0.2秒后的响应。A, B, C和D点（见图9-5）处的接头可以承受的最大拉力为100

11、kN，现在需要判断这些接头是否会断裂。图95 货物起重机载荷持续时间很短意味着惯性效应可能是很重要的，所以实质上为动态分析。该题没有提供结构阻尼的任何信息。由于桁架与交叉支撑之间采用的是螺栓连接，因此摩擦造成的能量损失是主要的。根据经验，对于每一阶振型可以选择5的临界阻尼。 作用载荷的大小与时间的关系见图9-6。图96 载荷时间特性在A.8节“Cargo crane dynamic loading，”中提供了一个命令执行文件。在ABAQUS/CAE中运行这个文件，可以生成该问题的完整分析模型。如果按照以下的指令步骤去做有困难，或者希望检查自己的工作，可以运行这个文件。在附录A中介绍了如何获取和

12、运行这个文件的内容。如果你无权使用ABAQUS/CAE 或其它的前处理程序，可以手工生成该问题的输入文件，在Getting Started with ABAQUS/Standard: Keywords Version 9.5节中的“例子：货物起重机动力载荷”中讨论了这个内容。 9.5.1 修改模型打开模型数据库文件Crane.cae，以Dynamic为名拷贝这个Static模型。除了以下所述的一些修改之外，动力分析的模型基本上和静力分析的模型是一样的。材料在动力分析中，每种材料的密度必须给定，这样才能形成质量矩阵。在起重机中钢的密度为7800 kg/m3。在这个模型中，材料属性是作为截面定义的

13、一部分予以配置的。所以需要在Property模块中编辑BracingSection和MainMemberSection截面定义来指定密度。在Edit Beam Section对话框的Section material density域中为每个截面输入值7800.0。注意: 如果材料数据的定义与截面属性无关，密度可在对材料定义进行编辑时予以指定，这可以通过在Edit Material 对话框中选择GeneralDensity 来给出。分析步 动力分析中的分析步定义和静力分析中的分析步定义完全不一样。因此，需将以前建立的静力分析步删除，并定义两个新的分析步。从Step模块中的主菜单条中，选择Step

14、Delete 删除静力分析步Tip load。与这个分析步相应的所有载荷，边界条件以及输出要求也将被删除。动力分析中的第一个分析步用于计算结构的自振频率和振型。第二个分析步则用这些数据来计算提升机的瞬态动力响应（模态）。我们假定这是一个线性的分析过程。如果想在分析中模拟非线性，必须使用隐式动力程序对运动方程进行直接积分。ABAQUS提供Lanczos 法和子空间迭代法来提取特征值。当系统的自由度很多，要求提取大量的模态时，一般来说Lanczos方法的速度更快。当所要求提取的模态很少（少于20阶）时，则用子空间迭代法更快。在该分析中，使用Lanczos法求解前30阶特征值。可以不指定所要提取的模

15、态的数目，而代之以指定所感兴趣的最小和最大的频率，当ABAQUS找到了在这个指定的范围内的所有特征值之后，该分析步的分析便结束。也可以指定一个截止（shift）点，在这个截止点附近的特征值将被提取。缺省情况下，不使用后两种方法。如果结构的刚体位移模式未被约束，应该把截止值设定为一个小的负值以避免由于刚体位移所带来的数值上的困难。生成频率提取分析步:1. 从主菜单条中，选择StepCreate。从有效的Linear perturbation程序表中选择Frequency ，并将分析步命名为Extract Frequencies。 2. 在Edit Step对话框的 Basic 选项卡中，输入分析

16、步的描述First 30 modes； 接受Lanczos 法选项，提取前30阶特征值。在结构动力分析中，响应通常和低阶模态有关。当然，应该提取足够的模态来表达结构的动力响应。检查是否提取了足够的模态的一个方法是查看每个自由度的全部有效质量，有效质量表明了在所提取模态的每个方向上被激活的质量的多少。数据文件中，在特征值数据之后是有效质量的表格形式。在理想情况下，对应于每个振型的每个方向，模型有效质量的总和至少应占总质量的90。在9.6节中对此作了进一步的讨论。下面用模态动力程序进行瞬态动力分析。瞬态响应将以第一个分析步中提取出来的所有模态为依据；在所有的30阶模态中均用5% 的临界阻尼比。建立

17、瞬态模态动力分析步:1. 从主菜单条中，选择StepCreate。从 Linear perturbation 程序表中选择Modal dynamics ，命名分析步为Transient modal dynamics。在前面所生成的频率提取分析步之后插入这个分析步。2. 在Edit Step对话框的Basic选项卡中，输入分析步的描述Simulation of Load Dropped on Crane 并指定时间周期为0.5，时间增量为0.005。在动力分析中时间是一个真实的物理量。3. 在Edit Step 对话框的Damping选项卡中指定直接模态阻尼并输入第1阶到第30阶模态的临界阻尼比

18、为0.05。如果使用了模态阻尼，必须指定在模态分析中使用的特征模态。缺省情况下，ABAQUS/CAE自动地选择所有可以得到的特征模态。当然也可以通过Keywords Editor，编辑*SELECT EIGENMODES块来改变缺省选择。在该问题中，接受缺省选择。输出结果 缺省情况下，ABAQUS 将振型写入输出数据库文件（.odb）中，以便能在Visualization 模块中绘制振型图。因此，可以接受和Extract Frequencies 分析步相关的缺省的场和历史输出请求。每阶振型的节点位移是经过规范化的，所以最大的位移为1。因此位移结果和相应的应力应变结果是没有物理意义的：它们只具有

19、相对的意义。完成动力分析通常比静力分析需要多得多的增量步。因此动力分析的输出量可能非常大，所以应该通过控制输出要求来保持结果文件大小的合理性。在本例中，要求每5个增量步结束时向输出数据库文件中输出位移。在该载荷步中有100个增量步（0.5/0.005）；因此共有20组场变量的输出。另外，在每个增量步将模型加载端（例如Tip-a集）的位移以及固定端（Attach集）的反力作为历史数据写入输出数据库中，以便利用这些数据作进一步的分析。在动态分析中，我们同样关心模型中的能量分布及其形式。模型中质点的运动产生动能；结构的位移产生应变能；能量由于阻尼而耗散。缺省情况下模态动力程序会将整个模型的能量作为历

20、史数据写入到.odb 文件中。瞬态模态动力分析步中的输出请求：1. 从主菜单条中，选择OutputField Output RequestsManager。 选择Transient modal dynamics栏（可能需要拉大表格栏才能看见完整的分析步名称）中的标识为Created的小块。2. 编辑场变量输出要求使得每隔5个增量步向.odb 文件中输出节点的位移一次。3. 从主菜单条中，选择OutputHistory Output RequestsManager。 在Transient modal dynamics 分析步中建立两个新的输出要求。在第一个输出要求中，输出Tip-a 集每个增量步

21、的位移，在第二个输出要求中，输出Attach集每个增量步的反力。载荷和边界条件边界条件和静力分析的相同。因为在删除静力分析步时将边界条件一并删除了，所以需在频率提取分析步中重新指定Attach集的边界条件。在结构端部施加一个集中力，集中力的大小是和时间相关的，如图96所示。与时间相关的载荷可以通过振幅曲线进行定义，在每个时间点上荷载的实际大小可以通过振幅曲线在该时间时的值乘上载荷的大小 (10,000N)而得。给出与时间相关的载荷:1. 先定义振幅。从Load模块的主菜单条中，选择ToolsAmplitudeCreate。命名振幅为Bounce，选择Tabular类型，在Edit Amplit

22、ude对话框中输入表91所示的数据。 接受缺省的Step time 选择作为时间步长，并输入0.25作为滤波参数值。注意：按下鼠标3键进入表格选项。表91 振幅曲线数据Time (sec)Amplitude0.00.00.011.00.21.00.210.02定义载荷。从主菜单条中，选择LoadCreate。在Transient modal dynamics分析步中施加载荷，命名为load Tip load，选择Concentrated force 作为载荷类型。给Tip-b 集施加载荷。先前定义的Tip-a 集与Tip-b 集之间的约束方程意味着载荷将被起重机的两部分均等地承担。3在Edit

23、 Load 对话框中，输入-1.E4 作为CF2 的值，选择amplitude 的类型为Bounce。在本例中，结构没有初始的速度和加速度，因此使用默认值。若想指定初始速度，则可以在主菜单条中选择FieldCreate ，并在初始分析步中给模型中所选择的区域指定初始速度，同时还需要编辑模态动力分析步的定义来引入初始条件。进行分析 在Job 模块中，以DynCrane 为名生成一个作业，描述为：3-D model of light-service cargo crane-dynamic analysis。将模型保存在模型数据库文件中，提交作业进行分析并监控求解过程；修改所发现的每一个错误；研究产

24、生警告信息的原因，如果有必要采取相应的措施。9.5.2结果输出检查Job Monitor 或者状态文件(.sta) 和输出数据文件(.dat) 来评估分析的结果。作业监控器和状态文件查看Job Monitor 或者状态文件DynCrane.sta的内容。可以发现在分析步1中的每一时间增量非常小。因为时间与频率提取程序是无关的，所以这个分析步没有占用时间。状态文件中的内容如下所示：状态文件中的输出表明第二个分析步中的时间增量大小为常数，并且每个增量步只需迭代一次。这是因为模态动力分析为振型的线性叠加，无需迭代。 出于同样的原因，在信息文件文中没有关于平衡和残差的信息。数据文件分析步1的主要结果是

25、所提取的特征值、振型参与系数和有效质量，如下所示：所提取的最高频率为96Hz，与此对应的周期为0.0104秒，可以将它与固定的时间增量0.005秒相比较。对于所提取的振型的周期远小于时间增量没有异议。相反地，时间增量必须能够求解所感兴趣的最高频率。在广义质量列中给出了对应于该阶振型的单自由度系统的质量。振型参与系数表反映了振型主要在哪个自由度上起作用。例如，根据结果可以看出第一阶振型主要在第三个方向上起作用。有效质量表反映了任一个模态在每个自由度上所激活的质量的大小。从中可以看出，在方向2上具有显著质量的第一个模态是第3阶。 该方向上总的模态有效质量为378.23 kg。在数据文件中已经给出了

26、模型总的质量为 414.34kg。 为了保证采用了足够的振型，每个方向上的总的等效质量应占模型质量的绝大部分（即90）。然而，模型中的某些质量是与约束节点相联系的，这些被约束的质量大致为与约束节点相连接的所有单元质量的1/4左右，在本例中，约为28kg。因此模型中可运动的质量是385kg。x-、 y- 和z- 方向上的有效质量分别为可运动质量的6，98和97。2-和3-方向上的有效质量远超过了先前所建议的90%，1-方向的则低得多。然而，由于载荷是作用在2-方向的，1-方向的响应就意义不大。 由于数据文件输出要求被关闭，所以数据文件没有包含任何模态动力分析步的结果。9.5.3后处理进入Visu

27、alization模块，并打开输出数据库文件DynCrane.odb。绘制振型 通过绘制与频率相应的振型来观察其变形形态。选择一个模态并绘制其变形形态:1. 从主菜单条中，选择ResultStep/Frame。弹出Step/Frame 对话框。2. 从Step Name列表中选择第一个分析步(Extract Frequencies) 。从Frame 表中选择Mode 1 。3. 从主菜单条中，选择PlotDeformed Shape；或者使用工具箱中的工具。ABAQUS/CAE 显示第一阶振型的变形形态，如图97所示。图 97 第一阶模态4. 在Step/Frame对话框中选择第三阶模态。5.

28、 点击OK。ABAQUS/CAE 即显示第三阶振型，如图98所示，Step/Frame 对话框关闭。图 98 第三阶模态结果的动画演示可以用动画演示分析的结果。首先建立第三阶模态动画演示的比例系数，然后生成一个瞬态结果的时间历程动画。建立特征模态动画的比例系数:1. 从主菜单条中，选择AnimateScale Factor；或者使用工具箱中的 工具。ABAQUS/CAE 以不同的变形比例系数（从0变为1）显示第三阶振型。ABAQUS/CAE 还在提示区的左边显示了媒体播放控制器。 用户可以利用这些控件来进行开始、停止以及跳过动画的操作。从左到右执行以下功能：放映、停止、第一幅图像、前一幅图像、

29、下一幅图像和最后一幅图像。2. 在提示区中，点击停止放映动画。建立瞬态结果的时间历程动画:1. 从主菜单条中，选择OptionsAnimation 。Animation Options 对话框弹出。2. 点击Time History 选项卡。3. 选择第二个分析步(Transient modal dynamics)。4. 单击OK 接受选择并关闭对话框。5. 从主菜单条中，选择 AnimateTime History；或者使用工具箱中的工具。在提示区的左边出现媒体播放控制器。6. 点击开始播放动画。ABAQUS/CAE 即播放第二个分析步中的每一帧画面。其中的状态块显示了当前的分析步和所经历的

30、增量。当达到了该分析步的最后一个增量步后，动画便自动地重放。7. 在动画的播放过程中可以根据用户自己的要求改变变形图：a. 显示Deformed Shape Plot Options 对话框。b. 从Deformation Scale Factor 域中选择Uniform。c. 输入15.0 作为变形比例系数值。d. 点击Apply 。ABAQUS/CAE 现在以15.0为变形比例系数播放第二个载荷步的每一帧图片。e. 从Deformation Scale Factor区中选择Auto-compute 。f. 点击OK 确认所作的选择并关闭Deformed Shape Plot Options

31、 对话框。ABAQUS/CAE现在以缺省的变形比例系数0.8播放第二个载荷步的每一帧图片。确定拉力的峰值 为了找出连接点上的拉力峰值，生成连接节点方向1的反力(变量 RF1)的XY 图。为此需要同时绘制多条曲线。绘制多条曲线:1. 从主菜单条中，选择ResultHistory Output。History Output 对话框弹出。2. 从Variables 选项卡中的Output Variables 域中，选择以下4条曲线（用Ctrl+Click）：Reaction Force: RF1 PI: TRUSS-1 Node xxx in NSET ATTACH3. 单击Plot。ABAQUS/

32、CAE即显示选择的曲线。4. 从主菜单条中，选择CanvasViewport Annotation Options。Viewport Annotation Options 对话框弹出。5. 单击Legend 选项卡，选择Show min/max values。6. 单击OK 确认所作的改变并关闭对话框。ABAQUS/CAE 显示最大值和最小值。图99中显示了这个结果图。顶端（B点和 C点）的两个节点的每根曲线几乎和桁架底端（A点和D点）的两个节点的曲线反对称。图 99 连接节点处反力的历程曲线每根桁架顶端连接点处的峰值拉力约为80kN，低于接头的承载能力100kN。记住方向1的反力为负值，这意

33、味着杆件被拉离墙面。当施加载荷时，下面的连接点受压（正的反力），但卸载之后反力在拉力和压力之间摆动。峰值拉力约为40 kN，远小于允许值。可通过查询XY图找到这些值。查询XY 图:1. 从主菜单条中，选择ToolsQuery。Query 对话框弹出。2. 在Visualization Queries域中选择Probe values。3. 单击OK。Probe Values 对话框弹出。4. 选择图99所示的点。该点的-坐标值是40.3 kN。这就是方向1的反力值。9.6模态阶数的影响在该例中采用了30阶模态来表征结构的动力特性。对应于所有模态的总的模态有效质量占到结构y方向和x方向可运动质量的

34、90以上，这表明当前的模态阶数充分反映出了结构的动力特性。图910显示的是Tip-a 集中的节点在第二个自由度方向的位移时间曲线，可以看出选取较少量模态后对结果的影响。 图 910 不同模态数对结果的影响如果检查有效质量表，会发现在方向2上起重要作用的第一个模态是第三阶模态，这就是只采用两个模态时无动力响应的原因。该节点自由度2方向的位移用5阶模态和用30阶模态分析的结果在0.2秒之后是相似的，当然在0.2秒之前的响应是有区别的，这表明在第5阶和第30阶模态之间存在着对早期响应起重要作用的模态。在使用5阶模态时，方向2总的模态有效质量仅占可动质量的57% 。9.7 阻尼的影响在该例中，对于所有

35、的振型均采用5的临界阻尼比。这个值是根据经验选择的。它基于这样一个事实：由于局部摩擦，连接桁架和交叉支撑的螺栓会吸收大量的能量。在这种情况下，难以得到准确的数据，所以分析所选数据对结果的影响是很重要的。图911比较了使用1%, 5%, and 10%的临界阻尼比时，顶部连接处的C点的反力随时间的变化。图911 阻尼对结果的影响正如预料的那样，低阻尼时的振动比高阻尼时衰减得慢，模型反力的峰值在低阻尼时要高些。即使当阻尼比很低，达到1%时，拉力的峰值也只有85kN，仍低于接头的强度（100kN）。因此，在此跌落载荷作用下，货物起重机依然能保持完好。9.8其它的动力程序现在简要的介绍一下ABAQUS

36、/Standard中其它的动力程序即线性模态分析和非线性模态分析9.8.1线性模态动力分析在ABAQUES/Standard中还有其它几个采用振型叠加技术的线性动力程序。 与模态动力程序计算的是时域上的响应不同，这些程序计算频域上的结果，这可以使我们从另外一个角度来分析结构的特性。这些程序的完整描述见ABAQUS/Standard用户手册的6.3节。.稳态动力分析 这个程序用于计算结构响应的振幅和相位，其荷载是谐波激励形式，荷载的频率在用户给定的范围之内。以下是一些典型的例子： 汽车发动机底座在发动机运转速度范围内的响应。 建筑物中的旋转机械。 飞机发动机的部件反应谱 这个程序提供了一种结构承

37、受固定端运动的动力作用时的峰值响应（位移、应力等）的计算途径。 固定端的运动被称作为“基础运动”，地震引起地面运动就是一个例子。当设计目的是估计峰值响应时，这是一种被采用的典型方法。随机响应该程序用于预测系统在承受随机连续激励时的响应。激励是采用统计意义上的能量谱密度函数来表示的。随机响应分析的例子如下： 飞机对扰动的响应。 结构对噪音的响应，例如来自喷气飞机的噪音。9.8.2非线性动态分析前面提到过，模态动力程序只适用于线性问题。如果对非线性动力响应感兴趣，必须对运动方程进行直接积分。在ABAQUS/Standard中运动方程的直接积分由一个隐式动力程序来完成。在使用这个程序时，在每个时间点

38、上都要建立质量、阻尼和刚度矩阵并求解动力平衡方程。由于这些操作的计算量很大，因此直接积分的动力分析要比模态方法昂贵。由于ABAQUS/Standard中的非线性动力程序采用隐式的时间积分，所以适用于求解非线性结构动力问题。例如，某一突然事件（如冲击）激发的结构的动态响应，再例如由于塑性或粘性阻尼造成大量的能量耗散的结构响应问题。在这些问题中，高频响应由于模型的耗能机制而迅速地衰减，这一点对于振动初期是十分重要的。当应力波在模型中传播时，显式算法是一次一个单元地扩展计算结果。因此，这很适合于求解应力波的影响非常重要的问题以及所模拟的突发事件的持续时间很短（典型的是小于1秒）的问题。和ABAQUS

39、/Standard相比，显式算法的另一个优点是它能更容易地模拟非连续的非线性，例如接触和失效问题。而大型、高度非连续性的问题，即使响应是准静态的，采用ABAQUS/Explicit一般的会很容易地予以模拟。9.9 相关的ABAQUS的例子 ABAQUS 例题手册第2.2.2节，“印度Point核反应堆供水线的线性分析” ABAQUS 基准校验手册第1.3.3节，“爆炸载荷作用的圆柱面板”。 ABAQUS 基准校验手册第1.4.6节，“悬臂板的特征值分析”。9.10建议阅读的文献 Clough, R. W. and J. Penzien, Dynamics of Structures, McGr

40、aw-Hill, 1975. NAFEMS Ltd., A Finite Element Dynamics Primer, 1993. Spence, P. W. and C. J. Kenchington, The Role of Damping in Finite Element Analysis, Report R0021, NAFEMS Ltd., 1993. 9.11小结 动力分析包括了结构的惯性效应。 频率提取程序可提取结构的自振频率和振型。 通过振型叠加技术，可用振型确定线性系统的动力响应。这一方法尽管有效，但是不能用于非线性问题。 线性动力程序可以计算瞬态载荷下的瞬态响应、谐波载荷下的稳态响应、支座运动的峰值响应以及随机载荷的响应。 为了准确表现结构的动力特性，必须提取足够多的振型。运动方向上总的模态有效质量应占结构总的可运动质量的90以上。 用户可以定义直接模态阻尼、Rayleigh阻尼和复合模态阻尼。但是由于固有频率和振型的计算都是基于无阻尼的结构，所以只能分析低阻尼的结构。 模态技术不适用于非线性的动力分析。在这类分析中必须采用直接时间积分方法。 用振幅曲线可以定义任意的随时间变化的载荷和给定的边界条件。 振型和瞬态结果可以在ABAQUS/CAE的Visualization模块中用动画显示。这对于理解动力响应和非线性静力分析是十分有帮助的。921