1、91第九章 动力问题如果只对结构加载荷后的长期响应感兴趣的话，静力分析就足够了。然而，如果加载时间很短，例如地震；或者载荷性质为动态，例如来自旋转机械的荷载，这时就必须采用动力分析。 9.1 引言动态模拟是将惯性力包含在动力学平衡方程中：0PIuM其中 M 是结构的质量。是结构的加速度。u是结构中的内力。IP 是所施加的外力。公式的表述无非是牛顿的第二运动定律（F=ma）的表现。动态分析和静态分析最主要的不同在于平衡方程中包含惯性力项（M ）。u两者的另一个不同之处在于内力 的定义。在静态分析中，内力仅由结构的变形I引起；而动态分析中的内力包括运动（例如阻尼）和结构变形的共同贡献。9.1.1
2、固有频率和模态最简单的动力问题是在弹簧上的质量振动，如图 9-1 所示。图 91 质量弹簧系统弹簧的内力为 ，所以运动方程为kumP0这个质量弹簧系统的固有频率（单位是弧度/秒）为 k如果质量块被移动后再释放，它将以这个频率振动。假若以此频率施加一个动态外力，位移的幅度将剧烈增加即所谓的共振现象。实际的结构具有多个固有频率。因此，在设计结构时避免使各固有频率与可92能的荷载频率过分接近就很重要。固有频率可以通过分析结构在无荷载（动力平衡方程中的 ）时的动态响应而得到。此时，运动方程变为 MuI0对于无阻尼系统， ，则上式变为 K这个方程解的形式为 tieu将此式代入到运动方程中便得到了特征值问
3、题方程 M其中 。 2该系统具有 个特征值，此处 是有限元模型的自由度数。记 为第 个特nnj征值。它的平方根 是结构的第 阶固有频率，并且 是相应的第 阶特征向量。jj j特征向量也就是所谓的模态（也称为振型），因为它是结构在第 阶振型下的变j形状态。在 ABAQUS 中，频率提取程序用来求解结构的振型和频率。这个程序使用起来十分简单，只要给出所需振型的数目和所关心的最高频率即可。9.1.2 振型叠加在线性问题中，结构在荷载作用下的动力响应可以用固有频率和振型来表示，即结构的变形可以采用振型叠加的技术由各振型的组合得到，每一阶模态都要乘以一个标量因子。模型中位移矢量 被定义为u1ii其中 是
4、振型 的标量因子。这一技术只在模拟小变形、线弹性材料、无接触条ii件情况下是有效的，即必须是线性问题。在结构动力学问题中，结构的响应往往取决于相对较少的几阶振型，这使得振型叠加方法在计算这类系统的响应时特别有效。考虑一个含有 10,000 个自由度的模型，则对运动方程的直接积分需要在每个时间点上求解 10,000 个联立方程组。但若结构的响应采用 100 阶振型来描述，那么在每个时间步上只需求解100 个方程。更重要的是，振型方程是解耦的，而原来的运动方程是耦合的。虽然在计算振型和频率时需要花费一些时间作为代价，但在计算响应时将节省大量的时间。如果在模拟中存在非线性，在分析中固有频率会发生明显
5、的变化，因此振型叠加法将不再适用。在这种情况下，需要对动力平衡方程直接积分，这将比振型分析要花费更多的时间。具有下列特点的问题才适于进行线性瞬态动力学分析： 93 系统应该是线性的：线性材料特性，无接触条件，无非线性几何效应。 响应应该只受较少的频率支配。当响应中各频率成分增加时，例如撞击和冲击问题，振型叠加技术的有效性将大大降低。 载荷的主要频率应在所提取的频率范围内，以确保对载荷的描述足够精确。 由于任何突然加载所产生的初始加速度应该能用特征模态精确描述。 系统的阻尼不能过大。9.2 阻尼如果一个无阻尼结构做自由振动，则它的振幅会保持恒定不变。然而，实际上由于结构运动而能量耗散，振幅将逐渐
6、减小直至振动停止，这种能量耗散称为阻尼。通常假定阻尼为粘滞的或正比于速度。动力平衡方程可以重新写成包含阻尼的形式： MuIPKC0其中 C 是结构的阻尼阵， 是结构的速度。能量耗散来自于诸多因素，其中包括结构结合处的摩擦和局部材料的迟滞效应。阻尼概念对于无需顾及能量吸收过程的细节表征而言是一个很方便的方法。在 ABAQUS 中，是针对无阻尼系统计算其振型的，然而，大多数工程问题还是包含阻尼的，尽管阻尼可能很小。有阻尼的固有频率和无阻尼的固有频率的关系是 d12其中是阻尼特征值，d是临界阻尼比，0c是该振型的阻尼，是临界阻尼。对 较小的情形（ ），有阻尼系统的特征频率非常接近于无阻尼系统1.0的
7、相应值。当 增大时，采用无阻尼系统的特征频率就不太准确，当 接近于 1时，就不能采用无阻尼系统的特征频率了。当结构处于临界阻尼（ ）时，施加一个扰动后，结构不会有摆动而是很快地恢复到静止的初始形态。（见图 9-2）94图 92 阻尼9.2.1 ABAQUS 中阻尼的定义在 ABAQUS 中，为了进行瞬时模态分析，可定义不同类型的阻尼：直接模态阻尼，瑞利(Rayleigh) 阻尼和复合模态阻尼。模拟动力学过程要定义阻尼。阻尼是分析步定义的一部分，每阶振型可以定义不同数量的阻尼。直接模态阻尼采用直接模态阻尼可以定义对应于每阶振型的临界阻尼比 。 的典型取值范围是从 1%到 10%。直接模态阻尼允许
8、精确定义每阶振型的阻尼。Rayleigh 阻尼Rayleigh 阻尼假设阻尼矩阵是质量矩阵和刚度矩阵的线性组合，KMC其中 和 是用户定义的常数。尽管假设阻尼正比于质量和刚度没有严格的物理基础，但实际上我们对于阻尼的分布知之甚少，也就不能保证使用更为复杂的阻尼模型是正确的。一般来讲，这个模型对于大阻尼系统也就是临界阻尼超过10%时，是失效的。相对于其它形式的阻尼，可以精确地定义系统每阶模态的Rayleigh 阻尼。95复合阻尼在复合阻尼中，可以定义每种材料的临界阻尼比，并且复合阻尼是对应于整体结构的阻尼。当结构中有许多不同种类材料时，这一选项是十分有用的。在这本指南中将不对复合阻尼做进一步的讨
9、论。9.2.2 阻尼值的选择在大多数线性动力学问题中，恰当地定义阻尼对于获得精确的结果是十分重要的。但是阻尼只是对结构吸收能量这种特性的近似描述，而不是去仿真造成这种效果的物理机制。所以，确定分析中所需要的阻尼数据是很困难的。有时，可以从动力试验中获得这些数据，但是在多数情况下，不得不通过经验或参考资料获得数据。在这些情况下，要仔细地分析计算结果，应该通过参数分析来评价阻尼系数对于模拟的敏感性。9.3 单元选择事实上，ABAQUS 的所有单元均可用于动力分析。选取单元的一般原则与静力分析相同。但是，在模拟冲击和爆炸载荷时，应选用一次单元。因为它们具有集中质量公式，在模拟应力波效果方面优于采用二
10、次单元的一致质量公式。9.4 动力问题的网格剖分在动力分析中，剖分网格需要考虑响应中将被激发的振型，网格剖分应能充分反映那些振型。这就意味着，能满足静态模拟要求的网格，不一定能计算高频振型的动态响应。例如，考虑图 93 所示的平板。用一次壳单元剖分的网格对于受均布载荷的静力分析以及一阶振型的预测是适用的。但是，该网格对于精确模拟第六阶振型就显得太粗糙了。图 93 平板的粗网格图 94 显示了同样的板用一次单元进行更精细的网格剖分后所模拟的结果。96第六阶振型的位移看起来明显好一些，预测的频率也准确了许多。如果作用在板上的动载会显著地激发该阶振型，则应采用精细的网格；采用粗网格将得不到精确的结果
11、。图 94 平板的精细网格9.5 例子：货物起重机动态载荷这个例子采用在第 6.4 节中分析过的同样的结构，现在要研究的问题是当10kN 荷重跌落在起重机挂钩上 0.2 秒后的响应。 A, B, C 和 D 点（见图 9-5）处的接头可以承受的最大拉力为 100kN，现在需要判断这些接头是否会断裂。图 95 货物起重机载荷持续时间很短意味着惯性效应可能是很重要的，所以实质上为动态分析。该题没有提供结构阻尼的任何信息。由于桁架与交叉支撑之间采用的是螺栓连接，因此摩擦造成的能量损失是主要的。根据经验，对于每一阶振型可以选择 5的临界阻尼。作用载荷的大小与时间的关系见图 9-6。97图 96 载荷时
12、间特性在 A.8 节“Cargo crane dynamic loading，”中提供了一个命令执行文件。在 ABAQUS/CAE 中运行这个文件，可以生成该问题的完整分析模型。如果按照以下的指令步骤去做有困难，或者希望检查自己的工作，可以运行这个文件。在附录 A 中介绍了如何获取和运行这个文件的内容。如果你无权使用 ABAQUS/CAE 或其它的前处理程序，可以手工生成该问题的输入文件，在 Getting Started with ABAQUS/Standard: Keywords Version 9.5 节中的“例子：货物起重机动力载荷”中讨论了这个内容。 9.5.1 修改模型打开模型数据
13、库文件 Crane.cae，以 Dynamic 为名拷贝这个 Static 模型。除了以下所述的一些修改之外，动力分析的模型基本上和静力分析的模型是一样的。材料在动力分析中，每种材料的密度必须给定，这样才能形成质量矩阵。在起重机中钢的密度为 7800 kg/m3。在这个模型中，材料属性是作为截面定义的一部分予以配置的。所以需要在Property 模块中编辑 BracingSection 和 MainMemberSection 截面定义来指定密度。在 Edit Beam Section 对话框的 Section material density 域中为每个截面输入值 7800.0。注意: 如果材
14、料数据的定义与截面属性无关，密度可在对材料定义进行编辑时予以指定，这可以通过在 Edit Material 对话框中选择 General Density 来给出。98分析步 动力分析中的分析步定义和静力分析中的分析步定义完全不一样。因此，需将以前建立的静力分析步删除，并定义两个新的分析步。从 Step 模块中的主菜单条中，选择 Step Delete 删除静力分析步 Tip load。与这个分析步相应的所有载荷，边界条件以及输出要求也将被删除。动力分析中的第一个分析步用于计算结构的自振频率和振型。第二个分析步则用这些数据来计算提升机的瞬态动力响应（模态）。我们假定这是一个线性的分析过程。如果想
15、在分析中模拟非线性，必须使用隐式动力程序对运动方程进行直接积分。ABAQUS 提供 Lanczos 法和子空间迭代法来提取特征值。当系统的自由度很多，要求提取大量的模态时，一般来说 Lanczos 方法的速度更快。当所要求提取的模态很少（少于 20 阶）时，则用子空间迭代法更快。在该分析中，使用 Lanczos 法求解前 30 阶特征值。可以不指定所要提取的模态的数目，而代之以指定所感兴趣的最小和最大的频率，当 ABAQUS 找到了在这个指定的范围内的所有特征值之后，该分析步的分析便结束。也可以指定一个截止（shift ）点，在这个截止点附近的特征值将被提取。缺省情况下，不使用后两种方法。如果
16、结构的刚体位移模式未被约束，应该把截止值设定为一个小的负值以避免由于刚体位移所带来的数值上的困难。生成频率提取分析步:1. 从主菜单条中，选择 Step Create。从有效的 Linear perturbation 程序表中选择 Frequency ，并将分析步命名为 Extract Frequencies。 2. 在 Edit Step 对话框的 Basic 选项卡中，输入分析步的描述 First 30 modes； 接受 Lanczos 法选项，提取前 30 阶特征值。在结构动力分析中，响应通常和低阶模态有关。当然，应该提取足够的模态来表达结构的动力响应。检查是否提取了足够的模态的一个方
17、法是查看每个自由度的全部有效质量，有效质量表明了在所提取模态的每个方向上被激活的质量的多少。数据文件中，在特征值数据之后是有效质量的表格形式。在理想情况下，对应于每个振型的每个方向，模型有效质量的总和至少应占总质量的 90。在 9.6节中对此作了进一步的讨论。下面用模态动力程序进行瞬态动力分析。瞬态响应将以第一个分析步中提取出来的所有模态为依据；在所有的 30 阶模态中均用 5% 的临界阻尼比。建立瞬态模态动力分析步:1. 从主菜单条中，选择 Step Create。从 Linear perturbation 程序表中选择 Modal dynamics ，命名分析步为 Transient mo
18、dal dynamics。在前面所生成的频率提取分析步之后插入这个分析步。2. 在 Edit Step 对话框的 Basic 选项卡中，输入分析步的描述 Simulation of Load Dropped on Crane 并指定时间周期为 0.5，时间增量为 0.005。在动力分析中时间是一个真实的物理量。993. 在 Edit Step 对话框的 Damping 选项卡中指定直接模态阻尼并输入第 1阶到第 30 阶模态的临界阻尼比为 0.05。如果使用了模态阻尼，必须指定在模态分析中使用的特征模态。缺省情况下，ABAQUS/CAE 自动地选择所有可以得到的特征模态。当然也可以通过Keyw
19、ords Editor，编辑 *SELECT EIGENMODES 块来改变缺省选择。在该问题中，接受缺省选择。输出结果 缺省情况下，ABAQUS 将振型写入输出数据库文件（.odb）中，以便能在Visualization 模块中绘制振型图。因此，可以接受和 Extract Frequencies 分析步相关的缺省的场和历史输出请求。每阶振型的节点位移是经过规范化的，所以最大的位移为 1。因此位移结果和相应的应力应变结果是没有物理意义的：它们只具有相对的意义。完成动力分析通常比静力分析需要多得多的增量步。因此动力分析的输出量可能非常大，所以应该通过控制输出要求来保持结果文件大小的合理性。在本例
20、中，要求每 5 个增量步结束时向输出数据库文件中输出位移。在该载荷步中有100 个增量步（ 0.5/0.005）；因此共有 20 组场变量的输出。另外，在每个增量步将模型加载端（例如 Tip-a 集）的位移以及固定端（Attach 集）的反力作为历史数据写入输出数据库中，以便利用这些数据作进一步的分析。在动态分析中，我们同样关心模型中的能量分布及其形式。模型中质点的运动产生动能；结构的位移产生应变能；能量由于阻尼而耗散。缺省情况下模态动力程序会将整个模型的能量作为历史数据写入到.odb 文件中。瞬态模态动力分析步中的输出请求：1. 从主菜单条中，选择 Output Field Output R
21、equests Manager。 选择 Transient modal dynamics 栏（可能需要拉大表格栏才能看见完整的分析步名称）中的标识为 Created 的小块。2. 编辑场变量输出要求使得每隔 5 个增量步向.odb 文件中输出节点的位移一次。3. 从主菜单条中，选择 Output History Output Requests Manager。 在 Transient modal dynamics 分析步中建立两个新的输出要求。在第一个输出要求中，输出 Tip-a 集每个增量步的位移，在第二个输出要求中，输出Attach 集每个增量步的反力。载荷和边界条件边界条件和静力分析的相
22、同。因为在删除静力分析步时将边界条件一并删除了，所以需在频率提取分析步中重新指定 Attach 集的边界条件。在结构端部施加一个集中力，集中力的大小是和时间相关的，如图 96 所示。与时间相关的载荷可以通过振幅曲线进行定义，在每个时间点上荷载的实际大小可以通过振幅曲线在该时间时的值乘上载荷的大小 (10,000 N)而得。910给出与时间相关的载荷:1. 先定义振幅。从 Load 模块的主菜单条中，选择Tools Amplitude Create。命名振幅为 Bounce，选择 Tabular 类型，在 Edit Amplitude 对话框中输入表 91 所示的数据。 接受缺省的 Step t
23、ime 选择作为时间步长，并输入 0.25 作为滤波参数值。注意：按下鼠标 3 键进入表格选项。表 91 振幅曲线数据Time (sec) Amplitude0.0 0.00.01 1.00.2 1.00.21 0.02定义载荷。从主菜单条中，选择 Load Create。在 Transient modal dynamics 分析步中施加载荷，命名为 load Tip load，选择 Concentrated force 作为载荷类型。给 Tip-b 集施加载荷。先前定义的 Tip-a 集与 Tip-b 集之间的约束方程意味着载荷将被起重机的两部分均等地承担。3在 Edit Load 对话框中
24、，输入-1.E4 作为 CF2 的值，选择 amplitude 的类型为 Bounce。在本例中，结构没有初始的速度和加速度，因此使用默认值。若想指定初始速度，则可以在主菜单条中选择 Field Create ，并在初始分析步中给模型中所选择的区域指定初始速度，同时还需要编辑模态动力分析步的定义来引入初始条件。进行分析 在 Job 模块中，以 DynCrane 为名生成一个作业，描述为： 3-D model of light-service cargo crane-dynamic analysis。将模型保存在模型数据库文件中，提交作业进行分析并监控求解过程；修改所发现的每一个错误；研究产生警
25、告信息的原因，如果有必要采取相应的措施。9.5.2 结果输出检查 Job Monitor 或者状态文件 (.sta) 和输出数据文件(.dat) 来评估分析的结果。作业监控器和状态文件查看 Job Monitor 或者状态文件 DynCrane.sta 的内容。可以发现在分析步 1中的每一时间增量非常小。因为时间与频率提取程序是无关的，所以这个分析步没有占用时间。状态文件中的内容如下所示：911状态文件中的输出表明第二个分析步中的时间增量大小为常数，并且每个增量步只需迭代一次。这是因为模态动力分析为振型的线性叠加，无需迭代。 出于同样的原因，在信息文件文中没有关于平衡和残差的信息。数据文件分析
26、步 1 的主要结果是所提取的特征值、振型参与系数和有效质量，如下所示：所提取的最高频率为 96 Hz，与此对应的周期为 0.0104 秒，可以将它与固定的时间增量 0.005 秒相比较。对于所提取的振型的周期远小于时间增量没有异议。相反地，时间增量必须能够求解所感兴趣的最高频率。在广义质量列中给出了对应于该阶振型的单自由度系统的质量。振型参与系数表反映了振型主要在哪个自由度上起作用。例如，根据结果可以看出第一阶振型主要在第三个方向上起作用。912有效质量表反映了任一个模态在每个自由度上所激活的质量的大小。从中可以看出，在方向 2 上具有显著质量的第一个模态是第 3 阶。 该方向上总的模态有效质
27、量为 378.23 kg。在数据文件中已经给出了模型总的质量为 414.34 kg。 为了保证采用了足够的振型，每个方向上的总的等效质量应占模型质量的绝大部分（即 90）。然而，模型中的某些质量是与约束节点相联系的，这些被约束的质量大致为与约束节点相连接的所有单元质量的 1/4 左右，在本例中，约为 28 kg。因此模型中可运动的质量是 385 kg。x- 、 y- 和 z- 方向上的有效质量分别为可运动质量的 6，98和 97。2-和 3-方向上的有效质量远超过了先前所建议的 90%，1- 方向的则低得多。然而，由于载荷是作用在 2-方向的，1-方向的响应就意义不大。 由于数据文件输出要求被
28、关闭，所以数据文件没有包含任何模态动力分析步的结果。9.5.3 后处理进入 Visualization 模块，并打开输出数据库文件 DynCrane.odb。913绘制振型 通过绘制与频率相应的振型来观察其变形形态。选择一个模态并绘制其变形形态:1. 从主菜单条中，选择 Result Step/Frame。弹出 Step/Frame 对话框。2. 从 Step Name 列表中选择第一个分析步(Extract Frequencies) 。从 Frame 表中选择 Mode 1 。3. 从主菜单条中，选择 Plot Deformed Shape；或者使用工具箱中的工具。ABAQUS/CAE 显示
29、第一阶振型的变形形态，如图 97 所示。图 97 第一阶模态4. 在 Step/Frame 对话框中选择第三阶模态。5. 点击 OK。ABAQUS/CAE 即显示第三阶振型，如图 98 所示，Step/Frame 对话框关闭。914图 98 第三阶模态结果的动画演示可以用动画演示分析的结果。首先建立第三阶模态动画演示的比例系数，然后生成一个瞬态结果的时间历程动画。建立特征模态动画的比例系数:1. 从主菜单条中，选择 Animate Scale Factor；或者使用工具箱中的 工具。ABAQUS/CAE 以不同的变形比例系数（从 0 变为 1）显示第三阶振型。ABAQUS/CAE 还在提示区的
30、左边显示了媒体播放控制器。 用户可以利用这些控件来进行开始、停止以及跳过动画的操作。从左到右执行以下功能：放映、停止、第一幅图像、前一幅图像、下一幅图像和最后一幅图像。2. 在提示区中，点击 停止放映动画。915建立瞬态结果的时间历程动画:1. 从主菜单条中，选择 Options Animation 。Animation Options 对话框弹出。2. 点击 Time History 选项卡。3. 选择第二个分析步(Transient modal dynamics)。4. 单击 OK 接受选择并关闭对话框。5. 从主菜单条中，选择 Animate Time History；或者使用工具箱中的
31、工具。在提示区的左边出现媒体播放控制器。6. 点击 开始播放动画。ABAQUS/CAE 即播放第二个分析步中的每一帧画面。其中的状态块显示了当前的分析步和所经历的增量。当达到了该分析步的最后一个增量步后，动画便自动地重放。7. 在动画的播放过程中可以根据用户自己的要求改变变形图：a. 显示 Deformed Shape Plot Options 对话框。b. 从 Deformation Scale Factor 域中选择 Uniform。c. 输入 15.0 作为变形比例系数值。d. 点击 Apply 。ABAQUS/CAE 现在以 15.0 为变形比例系数播放第二个载荷步的每一帧图片。e.
32、从 Deformation Scale Factor 区中选择 Auto-compute 。f. 点击 OK 确认所作的选择并关闭 Deformed Shape Plot Options 对话框。ABAQUS/CAE 现在以缺省的变形比例系数 0.8 播放第二个载荷步的每一帧图片。确定拉力的峰值 为了找出连接点上的拉力峰值，生成连接节点方向 1 的反力(变量 RF1)的XY 图。为此需要同时绘制多条曲线。绘制多条曲线:1. 从主菜单条中，选择 Result History Output。916History Output 对话框弹出。2. 从 Variables 选项卡中的 Output Va
33、riables 域中，选择以下 4 条曲线（用Ctrl+Click）：Reaction Force: RF1 PI: TRUSS-1 Node xxx in NSET ATTACH3. 单击 Plot。ABAQUS/CAE 即显示选择的曲线。4. 从主菜单条中，选择 Canvas Viewport Annotation Options。Viewport Annotation Options 对话框弹出。5. 单击 Legend 选项卡，选择 Show min/max values。6. 单击 OK 确认所作的改变并关闭对话框。ABAQUS/CAE 显示最大值和最小值。图 99 中显示了这个结果
34、图。顶端（ B 点和 C 点）的两个节点的每根曲线几乎和桁架底端（A 点和 D 点）的两个节点的曲线反对称。图 99 连接节点处反力的历程曲线每根桁架顶端连接点处的峰值拉力约为 80 kN，低于接头的承载能力100 kN。记住方向 1 的反力为负值，这意味着杆件被拉离墙面。当施加载荷时，下面的连接点受压（正的反力），但卸载之后反力在拉力和压力之间摆动。峰值拉力约为 40 kN，远小于允许值。可通过查询 XY 图找到这些值。917查询 XY 图:1. 从主菜单条中，选择 Tools Query。Query 对话框弹出。2. 在 Visualization Queries 域中选择 Probe v
35、alues。3. 单击 OK。Probe Values 对话框弹出。4. 选择图 9 9 所示的点。该点的 -坐标值是40.3 kN。这就是方向 1 的反力值。9.6 模态阶数的影响在该例中采用了 30 阶模态来表征结构的动力特性。对应于所有模态的总的模态有效质量占到结构 y方向和 x方向可运动质量的 90以上，这表明当前的模态阶数充分反映出了结构的动力特性。图 910 显示的是 Tip-a 集中的节点在第二个自由度方向的位移时间曲线，可以看出选取较少量模态后对结果的影响。 图 910 不同模态数对结果的影响如果检查有效质量表，会发现在方向 2 上起重要作用的第一个模态是第三阶模态，这就是只采
36、用两个模态时无动力响应的原因。该节点自由度 2 方向的位移用 5 阶模态和用 30 阶模态分析的结果在 0.2 秒之后是相似的，当然在 0.2 秒之918前的响应是有区别的，这表明在第 5 阶和第 30 阶模态之间存在着对早期响应起重要作用的模态。在使用 5 阶模态时，方向 2 总的模态有效质量仅占可动质量的 57% 。9.7 阻尼的影响在该例中，对于所有的振型均采用 5的临界阻尼比。这个值是根据经验选择的。它基于这样一个事实：由于局部摩擦，连接桁架和交叉支撑的螺栓会吸收大量的能量。在这种情况下，难以得到准确的数据，所以分析所选数据对结果的影响是很重要的。图 911 比较了使用 1%, 5%,
37、 and 10%的临界阻尼比时，顶部连接处的 C点的反力随时间的变化。图 911 阻尼对结果的影响正如预料的那样，低阻尼时的振动比高阻尼时衰减得慢，模型反力的峰值在低阻尼时要高些。即使当阻尼比很低，达到 1%时，拉力的峰值也只有 85 kN，仍低于接头的强度（100 kN）。因此，在此跌落载荷作用下，货物起重机依然能保持完好。9.8 其它的动力程序现在简要的介绍一下 ABAQUS/Standard 中其它的动力程序 即线性模态分析和非线性模态分析9199.8.1 线性模态动力分析在 ABAQUES/Standard 中还有其它几个采用振型叠加技术的线性动力程序。与模态动力程序计算的是时域上的响
38、应不同，这些程序计算频域上的结果，这可以使我们从另外一个角度来分析结构的特性。这些程序的完整描述见 ABAQUS/Standard 用户手册的 6.3 节。.稳态动力分析 这个程序用于计算结构响应的振幅和相位，其荷载是谐波激励形式，荷载的频率在用户给定的范围之内。以下是一些典型的例子： 汽车发动机底座在发动机运转速度范围内的响应。 建筑物中的旋转机械。 飞机发动机的部件反应谱 这个程序提供了一种结构承受固定端运动的动力作用时的峰值响应（位移、应力等）的计算途径。 固定端的运动被称作为“基础运动”，地震引起地面运动就是一个例子。当设计目的是估计峰值响应时，这是一种被采用的典型方法。随机响应该程序
39、用于预测系统在承受随机连续激励时的响应。激励是采用统计意义上的能量谱密度函数来表示的。随机响应分析的例子如下： 飞机对扰动的响应。 结构对噪音的响应，例如来自喷气飞机的噪音。9.8.2 非线性动态分析前面提到过，模态动力程序只适用于线性问题。如果对非线性动力响应感兴趣，必须对运动方程进行直接积分。在 ABAQUS/Standard 中运动方程的直接积分由一个隐式动力程序来完成。在使用这个程序时，在每个时间点上都要建立质量、阻尼和刚度矩阵并求解动力平衡方程。由于这些操作的计算量很大，因此直接积分的动力分析要比模态方法昂贵。由于 ABAQUS/Standard 中的非线性动力程序采用隐式的时间积分
40、，所以适用于求解非线性结构动力问题。例如，某一突然事件（如冲击）激发的结构的动态响应，再例如由于塑性或粘性阻尼造成大量的能量耗散的结构响应问题。在这些问题中，高频响应由于模型的耗能机制而迅速地衰减，这一点对于振动初期是十分重要的。当应力波在模型中传播时，显式算法是一次一个单元地扩展计算结果。因此，这很适合于求解应力波的影响非常重要的问题以及所模拟的突发事件的持续时间很短（典型的是小于 1 秒）的问题。920和 ABAQUS/Standard 相比，显式算法的另一个优点是它能更容易地模拟非连续的非线性，例如接触和失效问题。而大型、高度非连续性的问题，即使响应是准静态的，采用 ABAQUS/Exp
41、licit 一般的会很容易地予以模拟。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, McGraw-Hill, 1975. NAFEMS Ltd., A Finite Element Dynamics Primer, 1993. Spe
42、nce, 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
