abaqus系列教程-07线性动态分析.doc

1、7线性动态分析 如果你只对结构承受载荷后的长期响应感兴趣，静力分析（static analysis）是足 够的。 然而，如果加载时间很短（例如在地震中）或者如果载荷在性质上是动态的（例 如来自旋转机械的荷载），你就必须采用动态分析（dynamic analysis）。本章将讨论应 用 ABAQUS/Standard 进行线性动态分析； 关于应用ABAQUS/Explicit进行非线性动态 分析的讨论，请参阅第9 章“非线性显式动态分析”。 7.1引言 动态模拟是将惯性力包含在动力学平衡方程中： MuIP0 其中 M结构的质量。 u结构的加速度。 I在结构中的内力。 P所施加的外力。 在上面公式

2、中的表述是牛顿第二运动定律（F = ma ） 。 在静态和动态分析之间最主要的区别是在平衡方程中包含了惯性力（Mu）。在两 类模拟之间的另一个区别在于内力I的定义。在静态分析中，内力仅由结构的变形引 起；而在动态分析中，内力包括源于运动（例如阻尼）和结构的变形的贡献。 7.1.1固有频率和模态 最简单的动态问题是在弹簧上的质量自由振动，如图7-1 所示。 7-1 图 7 1质量弹簧系统 在弹簧中的内力给出为ku，所以它的动态运动方程为 mukuP0 这个质量弹簧系统的固有频率（natral frequency）（单位是弧度/秒（rad/s）给出为 k m 如果质量块被移动后再释放，它将以这个频

3、率振动。若以此频率施加一个动态外 力，位移的幅度将剧烈增加，这种现象即所谓的共振。 实际结构具有大量的固有频率。因此在设计结构时，非常重要的是避免使可能的 载荷频率过分接近于固有频率。通过考虑非加载结构（在动平衡方程中令P0）的 动态响应可以确定固有频率。则运动方程变为 MuI0 对于无阻尼系统，IKu ，因此有 MuKu0 这个方程的解具有形式为 u i e t 将此式代入运动方程，得到了特征值（ eigenvalue）问题 KM 2 其中 。 该系统具有 n 个特征值， 其中n是在有限元模型中的自由度数目。 记j是第j个 7-2 特征值； 它的平方根j是结构的第j阶模态的固有频率（natu

4、ral frequency ），而j是 相应的第j阶特征向量（eigenvector）。特征向量也就是所谓的模态（ mode shape）（也 称为振型），因为它是结构以第j阶模态振动的变形形状。 在 ABAQUS/Standard 中，应用频率的提取过程确定结构的振型和频率。这个过程 应用起来十分容易，你只要指出所需要的振型数目或所关心的最高频率即可。 7.1.2振型叠加 在线性问题中，可以应用结构的固有频率和振型来定性它在载荷作用下的动态响 应。采用振型叠加（ modal superposition）技术，通过结构的振型组合可以计算结构的 变形，每一阶模态乘以一个标量因子。在模型中的位移矢

5、量u定义为 u ii i 1 其中i是振型i的标量因子。 这一技术仅在模拟小变形、线弹性材料和无接触条件的 情况下是有效的，换句话说，即线性问题。 在结构的动力学问题中，结构的响应往往被相对较少的几阶振型控制，在计算这 类系统的响应时，应用振型叠加成为特别有效的方法。考虑一个含有10,000 个自由度 的模型， 对动态运动方程的直接积分将在每个时间点上同时需要联立求解10,000 个方 程。如果通过100 个振型来描述结构的响应，则在每个时间增量步上只需求解100个 方程。更重要的是，振型方程是解耦的，而原来的运动方程是耦合的。在计算振型和 频率的过程中， 开始时需要一点成本，但是， 在计算响

6、应时将会节省大量的计算花费。 如果在模拟中存在非线性，在分析中固有频率会发生明显的变化，因此振型叠加 法将不再适用。在这种情况下，只能要求对动力平衡方程直接积分，它所花费的时间 比振型分析昂贵得多。 必须具备下列特点的问题才适合于进行线性瞬态动力分析： 系统应该是线性的：线性材料行为，无接触条件，以及没有非线性几何效应。 响应应该只受相对少数的频率支配。当在响应中频率的成分增加时，诸如是打击和 碰撞的问题，振型叠加技术的效率将会降低。 7-3 载荷的主要频率应该在所提取的频率范围之内，以确保对载荷的描述足够精确。 应用特征模态，应该精确地描述由于任何突然加载所产生的初始加速度。 系统的阻尼不能

7、过大。 7.2阻尼 如果允许一个无阻尼结构做自由振动，则它的振幅会是一个常数。然而在实际中， 能量被结构的运动耗散，振动的幅度减小直至振动停止。这种能量耗散被称为阻尼 （damping）。通常假定阻尼为粘滞的，或者正比于速度。包含阻尼的动力平衡方程可 以重新写为 MuIP0 IKuCu 其中 C是结构的阻尼矩阵 u是结构的速度。 能量耗散来自于诸多因素，其中包括结构连接处的摩擦和局部材料的迟滞效应。 阻尼是一种很方便的方法，它包含了重要的能量吸收而又无需模拟具体的效果。 在 ABAQUS/Standard 中，特征模态的计算是关于无阻尼系统的。然而， 大多数工 程问题都包含某种阻尼，尽管阻尼可

8、能很小。对于每个模态，在有阻尼和无阻尼的固 有频率之间的关系是 d1 2 其中 d是阻尼特征值， c c 0 是临界阻尼比， c是该振型的阻尼， c是临界阻尼。 0 7-4 对于的较小值（0.1），有阻尼系统的特征频率非常接近于无阻尼系统的相 应值；当增大时，无阻尼系统的特征频率成为不太准确的；而当接近于 1 时，采 用无阻尼系统的特征频率就成为无效的。 如果结构是处于临界阻尼（1），在任何扰动后，结构不会有摆动而是尽可能 迅速地恢复到它的初始静止构形。（见图 7-2） 图 7 2阻尼 7.2.1在 ABAQUS/Standard中阻尼的定义 对于瞬时模态分析，在 ABAQUS/Standar

9、d 中可以定义一些不同类型的阻尼：直接 模态阻尼（direct modal damping），瑞利阻尼（Rayleigh damping）和复合模态阻尼 （composite modal damping ） 。 阻尼是针对模态动力学过程定义的，阻尼是分析步定义的一部分，每阶模态可以 定义不同量值的阻尼。 7-5 直接模态阻尼 应用直接模态阻尼可以定义与每阶模态相关的临界阻尼比，其典型的取值 范围是在临界阻尼的1%到 10%之间。直接模态阻尼允许用户精确地定义系统的 每阶模态的阻尼。 Rayleigh阻尼 在 Rayleigh 阻尼中，假设阻尼矩阵是质量和刚度矩阵的线性组合， CMK， 其中和是

10、由用户定义的常数。尽管阻尼是正比于质量和刚度矩阵的假设没有 严格的物理基础，实际上我们对于阻尼的分布知之甚少，也就不能保证其它更为 复杂的模型是正确的。一般的，这个模型对于大阻尼系统不可靠；即超过临界阻 尼的大约10%。相对于其它形式的阻尼，你可以精确地定义系统的每阶模态的 Rayleigh 阻尼。 对于一个给定模态i，临界阻尼值为 i，而 Rayleigh 阻尼值和的关系为 i2 2 i i 复合阻尼 在复合阻尼中，对于每种材料定义一个临界阻尼比，这样就得到了对应于整 体结构的复合阻尼值。当结构中有多种不同的材料时，这一选项是有用的。在本 指南中将不对复合阻尼做进一步的讨论。 7.2.2选择

11、阻尼值 在大多数线性动力学问题中，恰当地定义阻尼对于获得精确的结果是十分重要的。 但是，在某种意义上阻尼只是近似地模拟了结构吸收能量的特性，并非试图去模拟引 起这种效果的物理机制。因此，在模拟中确定所需要的阻尼数据是很困难的。偶尔， 你可以从动态试验中获得这些数据，但是，你不得不通过查阅参考资料或者经验获得 这些数据。在这些情况下，你必须十分谨慎地解释模拟结果，并通过参数分析研究来 评估模拟对于阻尼值的敏感性。 7-6 7.3单元选择 事实上，ABAQUS 的所有单元均可用于动态分析，选取单元的一般原则与静力分 析相同。但是，在模拟冲击和爆炸载荷时，应该选用一阶单元，因为它们具有集中质 量公式

12、，这种公式模拟应力波的效果优于二阶单元采用的一致质量公式。 7.4动态问题的网格剖分 当你正在设计应用于动态模拟的网格时，你需要考虑在响应中将被激发的振型， 并且使所采用的网格能够充分地反映出这些振型。这意味着能够满足静态模拟的网格， 不一定能够计算由于加载激发的高频振型的动态响应。 例如，考虑图7-3 所示的板。一阶壳单元的网格对于板受均布载荷的静力分析是 适合的，并也适合于一阶振型的预测。但是，该网格是明显地过于粗糙以至于不能够 精确地模拟第六阶振型。 图 7 3板的粗网格 图 7-4 显示了同样的板采用了一阶单元的精细网格的模拟。现在，第六阶振型的 位移形状看起来明显变好，对于该阶振型所

13、预测的频率更加准确。如果作用在板上的 动态载荷会显著地激发该阶振型，则必须采用精细的网格；采用粗网格将得不到准确 的结果。 7-7 图 7 4 板的精细网格 7.5例题：货物吊车动态载荷 这个例子采用在第6.4 节“例题：货物吊车”中已分析过的同样的货物吊车，现 在要求研究的问题是当10kN 的载荷在 0.2秒的时间中落到吊车挂钩上所引起的响应。 在 A, B, C和 D 点（见图7-5）处的连接仅能够承受的最大拉力为100 kN。你必须决 定这些连接的任何一个是否会断裂。 图 7 5货物吊车 加载的持续时间很短意味着惯性效应可能是很重要的，基本上要进行动态分析。 这里没有提供关于结构的阻尼的

14、任何信息。由于在桁架和交叉支撑之间采用的是螺栓 连接，因此由摩擦效应引起的能量吸收可能是比较显著的。因此，基于经验可以对每 7-8 一阶振型选择5的临界阻尼。 施加载荷的值与时间的关系，如图7-6 所示。 图 7 6 载荷时间特性 在本手册的在线文档第A.5 节“Cargo crane dynamic loading ”提供了输入文件。 当通过 ABAQUS/CAE运行这个输入文件时，将创建关于该问题的完整的分析模型。 根据下面给出的指导如果你遇到困难，或者如果你希望检查你的工作，则可以运行这 个输入文件。在附录A“Example Files”中，给出了如何提取和运行输入文件的指导。 如果你没

15、有进入ABAQUS/CAE或者其它的前处理器，可以人工创建关于这个问 题的输入文件， 关于这方面的讨论， 见 GettingStarted withABAQUS/Standard ： Keywords Version，第 9.5 节“Example：Cargo cranedynamic loading ” 。 7.5.1修改模型 打开模型数据库文件Crane.cae，将Static模型复制成一个名为Dynamic 7-9 的模型。除了下面描述的修改之外，动态分析的模型基本上与静力分析的模型相同。 材料 在动态分析中，必须给定每种材料的密度，这样才能形成质量矩阵。在吊车 中钢的密度为7800 k

16、g/m 3。 在这个模型中，材料属性是作为截面特性定义的一部分给出的。所以需要在 Property 模块中编辑BracingSection和 MainMemberSection截面定义来指 定密度。在 Edit Beam Section（编辑梁截面） 对话框的 Section materialdensity （截面材料密度）域中，为每个截面输入密度值为7800 。 注意 ：如果材料数据的定义是独立于截面属性的，通过编辑材料定义可以将密度 包括在内，即在Edit Material对话框中，选择General -Density 。 分析步 应用于动态分析的分析步定义与静力分析的分析步定义具有本质上

17、的不同。 因此，两个新的分析步将取代前面所建立的静力分析步。 在动态分析中的第一个分析步用于计算结构的自振频率和振型。第二个分析 步则应用这些数据来计算吊车的瞬态（模型）动态响应。在这个分析中，我们假 定一切都是线性的。如果你想在分析中模拟任何的非线性，必须使用隐式动态 （implicitdynamic）过程对运动方程进行直接积分。关于进一步的细节请参阅第 7.9.2 节“非线性动态分析” 。 ABAQUS/Standard 提供了 Lanczos 和子空间迭代（subspaceiteration）的特征 值提取方法。对于具有很多自由度的系统，当要求大量的特征模态时，一般来说 Lanczos

18、方法的速度更快。当需要仅少数几个（少于20）特征模态时，则应用子 空间迭代法的速度可能更快。 在这个分析中， 我们使用Lanczos 特征值求解器并求解前30 个的特征值。 除 了指定所要提取模态的数目，也可以指定所感兴趣的最小和最大频率，因此，一 旦 ABAQUS/Standard 已经提取了在这个指定范围内的所有特征值，就会结束该分 析步。也可以指定一个变换点（shift point），距离这个变换点最近的特征值将被 提取。在默认情况下，不使用最小或最大的频率或变换点。如果没有约束结构的 刚体模态，必须设置变换值为一个小的负值，以避免由于刚体运动产生的数值问 题。 7-10 采用频率提取分

19、析步代替静态分析步: 1从主菜单栏中，选择Step -Replace -TipLoad 。在 ReplaceStep （替换分 析步）对话框中，从Linearperturbation（线性摄动）过程的列表中，选择 Frequency（频率）。 将删除不能转换的模型参数。在本例中删除了集中力，因为在频率提取分析中 不能应用它们。但是，频率提取分析步继承了与静态分析步相关的边界条件和 输出需求。 2在EditStep （编辑分析步）对话框的Basic （基础）页中，输入分析步描述 First 30modes；接受 Lanczos 特征求解器选项；并要求前30 阶特征值。 3将分析步重新命名为Ext

20、ract Frequencies。 在结构动态分析中，响应通常地是与低阶模态有关。但是，应该提取足够的 模态以便较好地表达结构的动态响应。检查是否已经提取了足够数量的特征值的 一种方法是查看在每个自由度上的全部有效质量，它表明了在所提取模态的每个 方向上激活了多少质量。在数据文件（.dat）的特征值输出中，给出了有效质量 的列表。在理想情况下，对于每个振型在每个方向上，有效质量的总和应当至少 占总质量的90。在第7.6节“模态数目的影响”中将给出进一步的讨论。 应用模型动态过程进行瞬时动态分析。瞬时响应是基于在第一个分析步中提 取全部的模态；在全部的30 阶模态中均采用了5%的临界阻尼。 创建

21、瞬时模型动态分析步: 1 从主菜单栏中， 选择 Step -Create 。从 Linearperturbation过程列表中选择 Modal dynamics，并命名分析步为Transientmodaldynamics。在上面 定义的频率提取分析步之后插入这个分析步。 2 在 Edit Step对话框的Basic页中，输入分析步的描述Simulationof Load Droppedon Crane ，并指定分析步的时间为0.5和时间增量 （timeincrement） 为 0.005。在动态分析中，时间是一个真实的物理量。 3 在 EditStep 对话框的 Damping页中，指定直接模

22、态阻尼（direct modal）， 并对第 1 阶至第 30 阶的模态输入临界阻尼比为0.05 。 如果使用了模态阻尼，必须指定在基于模态的动态过程中使用的特征模态。 ABAQUS/CAE默认自动地选择所有可能得到的特征模态。当然, 如果你希望 改变默认的选择，也可以应用KeywordsEditor （关键词编辑器）编辑 7-11 *SELECT EIGENMODES块。在这个问题中，接受默认的选择。 输出 应用 FieldOutput Request Manager（ 场变量输出管理器） ，对于 Extract Frequencies分析步，修改场变量输出设置，因此，选择了Preselec

23、ted defaults （预选默认值） 。在默认情况下， ABAQUS/Standard 将振型写入到输出数 据库文件（ .odb ），以便应用Visualization 模块绘制振型图。对于每阶振型的节 点位移都是经过单位化的，所以最大的位移为1单位。因此，这些结果和对应的 应力和应变是没有物理意义的：它们仅能够用于相互的比较。 完成动态分析通常比静态分析需要更多更多的增量步。做为结果，来自动态 分析的输出量可能是非常大，你应该控制输出要求以确保输出文件具有一个合理 的量。在本例中，要求在每第5个增量步结束时，向输出数据库文件中输出一次 位移形状。在分析步中有100个增量步（ 0.5/0.

24、005）；因此，有20 组场变量输出。 另外在每个增量步，将在模型加载端（例如，Tip-a集合）的位移和在固定 端（Attach集合）的约束反力作为历史数据写入到输出数据库文件中，以便从 这些数据中得到更好的解答。在动态分析中，我们也关心在模型中的能量分布以 及能量采用的形式。在模型中表现出的动能是质量运动的结果；表现出的应变能 是结构位移的结果；通过阻尼也耗散了能量。在默认情况下， 对于模型动态过程， 整个模型的能量将作为历史数据写入到.odb文件中。 对瞬时模型动态分析步中的输出请求： 1 从主菜单栏中，选择Output -FieldOutputRequests-Manager 。在标记

25、Transient modal dynamics的列中（可能需要拉大这列表格才能看见完整的 分析步名称），选择标有Created的单元格。 2 编辑场变量输出要求，使得仅将每第5 个增量步的节点的位移写入到.odb文 件中。 3 从主菜单栏中，选择Output -History Output Requests-Manager 。在标 记 Transient modal dynamics的分析步中创建两个新的输出要求。在第 一个中， 输出集合Tip-a在每个增量步结束时的位移；在第二个中， 输出集合 Attach在每个增量步结束时的约束反力。 7-12 载荷和边界条件 边界条件与在静力分析中的条

26、件相同。由于在分析步替换过程中保留了这些 条件，无需再定义新的边界条件。 在吊车的端部施加一个集中力，它的量级是与时间相关的，如图7-6 所示。 与时间相关的载荷可以应用幅值曲线（amplitude curve）进行定义。 通过幅值曲线 上的值乘以载荷的值（- 10,000 N），可以获得当时的任意点处施加荷载的实际值。 指定与时间相关的载荷: 1 首 先 定 义 幅 值 。 从 在Load模 块 的 主 菜 单 栏 中 ， 选 择Tools - Amplitude-Create ，命名幅值为Bounce ，并选择 Tabular （数据表）类型。 在 EditAmplitude（编辑幅值）对

27、话框中，输入在表7-1 中所示的数据。接受 默认的 Step time （分析步时间）的选择作为时间跨度，并输入0.25作为光滑 参数值。 注意 ：点击鼠标键3，进入表格选项。 表 7 1 幅值曲线数据 时间（秒）幅值 0.00.0 0.011.0 0.21.0 0.210.0 2现在定义载荷。从主菜单栏中，选择Load -Create 。在 Transient modal dynamics分析步中施加载荷， 命名为 Tip load ，并选择 Concentrated force （集中力）作为载荷类型。 施加载荷到集合Tip-b。在集合 Tip-a和 Tip-b之 间，前面定义的约束方程意

28、味着吊车的两半部分将平均地分担载荷。 3在 Edit Load对话框中，输入-1.E4作为 CF2（2 方向作用力）的值，并为幅 值选择 Bounce 。 在本例中，结构默认没有初始的速度或者加速度。然而，如果你希望定义初 始的速度，则你可以做。通过在主菜单栏中，选择Field -Create ，并在分析步 开始时，将初始平移速度设置到在模型中所选择的区域。为了引入初始条件，你 也需要编辑模型动态分析步的定义。 7-13 运行分析 在作业（ Job）模块中，创建一个名为DynCrane的作业，采用下面的描述： 3-D model of light-servicecargo crane-dyna

29、mic analysis。 将模型保存在模型数据库文件中，并提交作业进行分析和监控求解过程；纠 正发现的任何一个模拟错误，研究引起任何警告信息的原因，如果必要则采取相 应的措施。 7.5.2结果 在分析中对于每一个增量步，JobMonitor （作业监视器）给出了所采用的自动 时间增量步的简明总结。一旦该增量步结束就立刻写出相应的信息，这样你可以在作 业运行中监控分析的过程。对于大型、复杂的问题， 这个功能十分有用。 在 JobMonitor 中给出的信息与在状态文件（DynCrane.sta）中给出的信息相同。 查看 JobMonitor和打印的输出数据文件（DynCrane.dat）以评估

30、分析结果。 JobMonitor 在 JobMonitor中，第 1 列显示了分析步编号，第2 列给出了增量步编号。 在每个增量步中为了得到收敛的结果，第 6列显示了 ABAQUS/Standard 所需要的 迭代次数。观察JobMonitor的内容，可以发现在分析步1中与单一增量步相关 的时间增量非常小。因为时间是与频率提取过程无关的，所以这个分析步没有占 用时间。 在分析步 2 的输出显示，在整个分析步中时间增量的大小保持为常数，并且 每个增量步只需迭代一次。在图7-7 中显示了JobMonitor的结束部分。 7-14 图 7-7Job Monitor的结束部分：货物吊车动态分析 数据文

31、件 对于分析步1 的主要结果是提取的特征值（eigenvalue）、参与系数（ participation factor）和有效质量（effectivemass），如下所示： 7-15 所提取的最高频率为96Hz，与此频率对应的周期为0.0104 秒，可以将它与 固定的时间增量0.005 秒相比较。在所提取的振型中，其周期没有远小于时间增 量的。相反地，时间增量必须能够求解感兴趣的最高频率。 广义质量列（generalized mass）给出了对应于该阶振型的单自由度系统的质 量。 振型参与系数（participation factor）列表反映了在哪个自由度上该振型起主 导作用。例如，根据结

32、果可以看出1 阶振型主要在3 方向上起作用。 7-16 有效质量（ effective mass）列表反映了对于任何一个模态在每个自由度上所 激活的质量的大小。结果表明，在方向2 上具有显著质量的第一个模态是第3 阶 模态。在该方向上总的模型有效质量为378.23 kg。 前面在数据文件中，给出了模型的总质量为414.34kg。 为了保证已经采用了足够的模态，在每个方向上的总有效质量必须占模型质 量的绝大部分（即90）。然而，模型中的某些质量是与约束节点相联系的，这 些约束的质量占与约束节点相连接的所有单元质量的大约1/4，在本例中，约为 7-17 28 kg。因此，在模型中能够运动的质量是3

33、85 kg。在 x-、y-和 z-方向上的有效质 量分别为可运动质量的6%，98%和 97%。在 2-和 3-方向上的总有效质量远远超 过了前面所建议的90%，在 1-方向上的总有效质量是低得多。然而，由于载荷是 作用在 2-方向上的，在1-方向上的响应是不明显的。 对于模型动态分析步，由于关闭了所有数据文件的输出要求，所以在数据文 件中没有包含任何结果。 7.5.3后处理 进入 Visualization 模块，并打开输出数据库文件DynCrane.odb。 绘制振型 通过绘制与该频率相应的振型可以观察与一个给定的频率相应的变形状态。 选择一个模态并绘制对应的振型: 1从主菜单栏中，选择Re

34、sult -Step/Frame。 弹出 Step/Frame （分析步 /画面）对话框。 2 从StepName （ 分 析 步 名 称 ） 表 中 ， 选 择 第 一 个 分 析 步 （ Extract Frequencies）。 从 Frame 列表中，选择Mode 1 。 3从主菜单栏中， 选择 Plot -DeformedShape；或者使用工具箱中的工具。 ABAQUS/CAE显示了关于第一阶振型的变形形态，如图7 8 所示。 7-18 图 7 8第一阶模态 4从 Step/Frame对话框中，选择第三阶模态。 5点击 OK 。 ABAQUS/CAE显示出第三阶振型，如图7-9 所

35、示，并且Step/Frame对话框 消失了。 7-19 图 7 9第三阶模态 结果的动画演示 用动画（ animate）演示分析的结果。首先创建一个第三阶特征模态的动画放 大系数，然后创建一个瞬时结果的时间历史动画（time-history animation ） 。 创建一个特征模态的动画放大系数: 1从主菜单栏中， 选择 Animate -Scale Factor；或者使用工具箱中的工具。 通过从 0 到 1 不同的变形放大系数，ABAQUS/CAE显示第三阶振型和步骤。 在提示区的左侧，ABAQUS/CAE也显示了电影播放控制器。 2在提示区中，点击停止动画。 创建瞬时结果的时间历史动画

36、: 1从主菜单栏中，选择Options -Animation观察动画选项。 ABAQUS/CAE显示 Animation Options（动画选项）对话框。 2点击 Time History（时间历史）页。 7-20 3选择第二个分析步（Transient modal dynamics）。 4点击 OK 接受以上的选择并关闭对话框。 5从主菜单栏中， 选择 Animate -TimeHistory；或者使用工具箱中的工具。 ABAQUS/CAE在提示区的左边显示出电影播放控制器，并开始播放第二个分析 步中的每一帧画面。状态块 （status block）在动画放映过程中显示了当前的分析 步和增

37、量步。在达到了该分析步的最后一个增量步后，动画便自动地重播。 6在动画的播放过程中，你可以根据需要改变变形图。 a.显示 DeformedShape Plot Options（变形形状图选项）对话框。 b.在 Deformation Scale Factor（变形放大系数） 域中，选择 Uniform （一 致性）。 c. 输入 15.0作为变形放大系数值。 d. 点击 Apply ，采用所作的修改。 现在，ABAQUS/CAE以变形放大系数为15.0播放第二个载荷步的每一帧 图片。 e. 在 Deformation Scale Factor域中，选择Auto-compute（自动计算）。 f

38、. 点击 OK 采用所作的修改，并关闭Deformation Scale Factor对话框。 现在， ABAQUS/CAE以默认的变形放大系数为0.8播放第二个载荷步的 每一帧图片。 确定拉力的峰值 为了找出固定连接点处的拉力峰值，创建在固定连接点处在1方向的约束反 力（变量RF1）的 X Y曲线图。在曲线图中可以同时绘制多条曲线。 绘制多条曲线: 1从主菜单栏中，选择Result -History Output。 ABAQUS/CAE显示 History Output（历史输出）对话框。 2从 Variables （变量）选项页中的OutputVariables （输出变量）域中，选择 具

39、有以下形式的4 条曲线（用 Ctrl+ 点击 ）： ReactionForce: RF1 PI:TRUSS-1 Node xxxinNSET ATTACH 3点击 Plot （绘图）。 ABAQUS/CAE显示选择的曲线。 7-21 4从主菜单栏中，选择Viewport- Viewport Annotation Options。 ABAQUS/CAE显示 Viewport Annotation Options（图形窗标注选项） 对话框。 5点击 Legend （图例） 页，并选中 Show min/maxvalues（显示最小 /最大值） 。 6点击 OK 确认所作的修改，并关闭对话框。 AB

40、AQUS/CAE显示出最大值和最小值。 结果图显示在图7-10 中（用户可以修改） 。对于在每榀桁架的顶端的两个节 点（B 点和 C 点）的曲线几乎是在每榀桁架的底端的两个节点（A 点和 D 点）的 曲线的反射。 图 7 10在固定端连接点处约束反力的历史 在每个桁架结构的顶端的固定连接点处的峰值拉力约为80 kN，它低于连接 点 100 kN 的承载能力。 注意到在1 方向的约束反力为负值，意味着杆件被拉出墙 体。当施加载荷时，下面的连接点受压（正的约束反力），但是在卸载之后，约束 反力在拉力和压力之间振荡。峰值拉力是约为40 kN，远小于允许值。可通过观 察 X Y图发现这些值。 7-22

41、 查看XY图: 1从主菜单栏中，选择Tools -Query 。 显示 Query （查询）对话框。 2在 Visualization Queries（可视化查询） 域中， 选择 Probe values（查看值） 。 3单击 OK 。 显示 Probe Values对话框。 4选择在图7-10 所示的点。 该点的 Y-坐标值是 -40.3 kN，它对应于在1方向的约束反力值。 7.6模态数量的影响 对于这个模拟，采用了30 个模态来表征结构的动力特性。这30 个模态的总模态 有效质量占到在y-方向和 z-方向可运动的结构质量的90以上， 这表明已经充分地反 映了结构的动态特性。 图 7-11

42、显示的是在集合Tip-a中的节点在第2个自由度方向的位移时间曲线， 说明了使用少量的模态对结果质量的影响。 图 7 11不同模态数量对结果的影响 7-23 如果检查有效质量列表，你会发现在2 方向上起重要作用的第一个模态是第3阶 模态，可见仅当采用两个模态时的动态响应是不足的。分析该节点在自由度2 方向的 位移， 采用五个模态与三十个模态的结果在0.2秒之后是相似的；但是， 在 0.2 秒之前 的响应却是有区别的， 这表明在第5 至第 30 阶模态之间存在着对早期响应起重要作用 的模态。在采用五个模态时，在2 方向上总的模态有效质量仅占可运动质量的57%。 8.8阻尼的影响 在这个模拟中，对所

43、有的模态均采用5%的临界阻尼。这个值是根据经验选择的， 它基于这样一个事实：作为局部摩擦效应的结果，在桁架和交叉支撑之间的螺栓连接 可能吸收显著的能量。在这种难以得到准确数据的情况下，研究所选取的数据对结果 的影响是很重要的。 当使用 1%、5%和 10%的临界阻尼时，图7-12 比较了在顶部连接处的一点处（C 点）的约束反力的变化历史。 图 7 12阻尼对结果的影响 7-24 正如所预料的那样，在高阻尼水平时比在低阻尼水平时的振动衰减得快的多，并 且在采用低阻尼时在模型中力的峰值是更高一些。即使当阻尼低到1%时，拉力的峰 值为 85kN，它仍低于连接的强度（100kN）。因此，在此跌落载荷作

44、用下，货物吊车 依然能够保持完好的状态。 7.8与直接时间积分的比较 由于这是个瞬时动态分析，所以会很自然地想到将结果与对运动方程采用直接积 分得到的结果进行比较。进行直接积分可以或者采用隐式方法（ABAQUS/Standard ） 或者采用显式方式（ABAQUS/Explicit ）。这里我们采用显式动态过程以扩展该分析。 直接比较前面给出的结果是不可能的，因为在ABAQUS/Explicit中没有提供B33 单元类型和临界阻尼。因此，在ABAQUS/Explicit分析中，单元类型改换成为B31 和 采用了 Rayleigh 阻尼以代替临界阻尼。 将 Dynamic模型复制成一个名为exp

45、licit的新模型， 必须对 explicit模型 进行如下的修改。 修改模型： 1 删除模型动态分析步。 2 用一个显式动态分析步（dynamic, explicit）替换保留下的频率提取分析步，并 指定分析步时间期限为0.5s，另外，为了应用几何线性（取消几何非线性 Nlgeom ），编辑分析步。这将导致一个线性分析。 3 将分析步改名为Transient dynamics。 4 创建两个新的历史变量输出要求。第一个要求输出集合Tip-a的位移历史； 第 二个要求输出集合Attach的约束反力历史。 5 在支撑的截面特性中，添加质量比例阻尼（在主菜单栏中，选择Section - Edit

46、-BracingSection；在截面编辑器中，点击Damping ）。 对 alpha 采用的值为15 ，而其它的值保持为0。 对于在结构的低阶和高阶频率上临界阻尼的值，这些值作出了一个合理的权衡。 对于三个最低的固有频率，的有效值是大于0.05，但是如图7-11所示，前两 阶模态对于响应没有做出显著的贡献。对于余下的模态，的有效值均小于 0.05。随固有频率的变化如图7-13 所示。 7-25 尼 阻 界 临 频率 图 7-13 阻尼对结果的影响 6 对于主要构件的截面特性，重复上述步骤。 7 重新定义在集合Tip-b的尖端载荷。设置CF2 = -10000，并使用幅值定义 Bounce

47、。 8 改变单元库为Explicit，对于模型的所有区域设置单元类型为B31。 9创建一个新作业，命名为expDynCrane，并将其提交分析。 当作业完成后，进入Visualization模块查看结果。特别是，比较在前面从 ABAQUS/Standard得到的与现在从ABAQUS/Explicit得到的尖端位移历史。如图 7-14 所示， 在动态响应方面它们的区别很小。这些区别是由于在模型动态分析中采 用了不同的单元和阻尼类型。实际上，如果修改ABAQUS/Standard的分析，使其 采用 B31 单元和质量比例阻尼，那么由两种分析方法得到的结果几乎没有区别（见 图 7-14），这确认了模

48、型动态方法的准确性。 7-26 图 7-14比较由 ABAQUS/Standard 和 ABAQUS/Explicit得到的尖端位移 7.9其它的动态过程 现在简要的回顾在ABAQUS 中的其它动态过程，即线性模型动态（linear modal dynamics）和非线性动态（nonlinear dynamics） 。 7.9.1线性模态法的动态分析 在 ABAQUES/Standard 中还有其它几种采用了振型叠加技术的线性、动态过程。 与模型动态过程在时域（time domain）上计算响应不同，这些过程在频域（frequency domain）上提供结果，这可以使我们从另外的角度来分析结

49、构的行为。 在 ABAQUS 分析用户手册的第6.3节“ Dynamic stress/displacement analysis”中 7-27 给出了这些过程的完整描述。 稳态动态（ Steady-state dynamics） 在用户指定频率范围内的谐波激励下，这个程序用于计算引起结构响应的振 幅和相位。以下是一些典型的例子： 汽车发动机支座在发动机运转速度范围内的响应。 在建筑物中的旋转机械。 飞机发动机的部件。 反应谱（ Response spectrum） 当结构承受在它的固定点处的动态运动时，这个程序提供了对峰值响应的评 估（位移、应力等） 。固定点处的运动是所谓的“基础运动”（b

50、ase motion）；地震 发生时引起的地面运动就是一个例子。当为了设计需要估计峰值响应时，这是一 种典型的方法。 随机响应（ Random response） 在承受随机连续的激励时，该程序用于预测系统的响应。激励是采用具有统 计意义的能量谱密度函数来表示的。以下是随机响应分析的例子： 飞机对扰动的响应。 结构对噪音的响应，例如喷气发动机产生的噪音。 7.9.2非线性动态分析 如前面所述，模型动态过程仅适用于线性问题。当对非线性动态响应感兴趣时， 必须对运动方程进行直接地积分。在 ABAQUS/Standard 中，完成对运动方程的直接积 分是采用了一个隐式动态过程。在应用这个过程时，在每