1、9显式非线性动态分析 在前面的章节中，已经考察了显式动态程序的基本内容；在本章中，将对这个问 题进行更详细的讨论。显式动态程序对于求解广泛的、各种各样的非线性固体和结构 力 学 问 题 是 一 种 非 常 有 效 的 工 具 。 它 常 常 对 隐 式 求 解 器 是 一 个 补 充 ， 如 ABAQUS/Standard ；从用户的观点来看，显式与隐式方法的区别在于： 显式方法需要很小的时间增量步，它仅依赖于模型的最高固有频率，而与载荷的类 型和持续的时间无关。通常的模拟需要取10,000 至 1,000,000 个增量步，每个增量 步的计算成本相对较低。 隐式方法对时间增量步的大小没有内在
2、的限制；增量的大小通常取决于精度和收敛 情况。典型的隐式模拟所采用的增量步数目要比显式模拟小几个数量级。然而，由 于在每个增量步中必须求解一套全域的方程组，所以对于每一增量步的成本，隐式 方法远高于显式方法。 了解两个程序的这些特性，能够帮助你确定哪一种方法是更适合于你的问题。 9.1ABAQUS/Explicit适用的问题类型 在讨论显式动态程序如何工作之前，有必要了解ABAQUS/Explicit适合于求解哪 些 类 问 题 。 贯 穿 这 本 手 册 ， 我 们 已 经 提 供 了 贴 切 的 例 题 ， 它 们 一 般 是 应 用 ABAQUS/Explicit求解的如下类型问题： 高
3、速动力学（high-speed dynamic）事件 最 初 发 展 显 式 动 力 学 方 法 是 为 了 分 析 那 些 用 隐 式 方 法 （ 如 ABAQUS/Standard ）分析起来可能极端费时的高速动力学事件。作为此类模拟的 例子， 在第 10 章“材料” 中分析了一块钢板在短时爆炸载荷下的响应。因为迅速 施加的巨大载荷，结构的响应变化的非常快。对于捕获动力响应，精确地跟踪板 内的应力波是非常重要的。由于应力波与系统的最高阶频率相关联，因此为了得 到精确解答需要许多小的时间增量。 9-1 复杂的接触（contact ）问题 应用显式动力学方法建立接触条件的公式要比应用隐式方法容
4、易得多。结论 是 ABAQUS/Explicit能够比较容易地分析包括许多独立物体相互作用的复杂接触 问题。 ABAQUS/Explicit是特别适合于分析受冲击载荷并随后在结构内部发生复 杂相互接触作用的结构的瞬间动态响应问题。在第 12章“接触” 中展示的电路板 跌落试验就是这类问题的一个例子。在这个例子中，一块插入在泡沫封装中的电 路板从1m 的高度跌落到地板上。这个问题包括封装与地板之间的冲击，以及在 电路板和封装之间的接触条件的迅速变化。 复杂的后屈曲（postbuckling）问题 ABAQUS/Explicit能够比较容易地解决不稳定的后屈曲问题。在此类问题中， 随着载荷的施加，
5、结构的刚度会发生剧烈的变化。在后屈曲响应中常常包括接触 相互作用的影响。 高度非线性的准静态（quasi-static）的问题 由于各种原因，ABAQUS/Explicit常常能够有效的解决某些在本质上是静态 的问题。准静态过程模拟问题包括复杂的接触，如锻造、滚压和薄板成型等过程 一般地属于这类问题。薄板成型问题通常包含非常大的膜变形、褶皱和复杂的摩 擦接触条件。块体成型问题的特征有大扭曲、瞬间变形以及与模具之间的相互接 触。在第 13 章“ABAQUS/Explicit准静态分析”中，将展示一个准静态成型模拟 的例子。 材料退化（ degradation）和失效（ failure ） 在隐式
6、分析程序中，材料的退化和失效常常导致严重的收敛困难，但是 ABAQUS/Explicit能够很好地模拟这类材料。混凝土开裂的模型是一个材料退化 的例子，其拉伸裂缝导致了材料的刚度成为负值。金属的延性失效模型是一个材 料失效的例子，其材料刚度能够退化并且一直降低到零，在这段时间中，单元从 模型中被完全除掉。 这些类型分析的每一个问题都有可能包含温度和热传导的影响。 9-2 9.2动力学显式有限元方法 这一节包括ABAQUS/Explicit求解器的算法描述，在隐式和显式时间积分之间进 行比较，并讨论了显式方法的优越性。 9.2.1显式时间积分 ABAQUS/Explicit应用中心差分方法对运动
7、方程进行显示的时间积分，应用一个 增量步的动力学条件计算下一个增量步的动力学条件。在增量步开始时，程序求解动 力学平衡方程， 表示为用节点质量矩阵M 乘以节点加速度u等于节点的合力 （在所施 加的外力P 与单元内力I 之间的差值）： Mu = P - I 在当前增量步开始时（t 时刻），计算加速度为： 1 u |(M )(PI )| (t)(t) 由于显式算法总是采用一个对角的、或者集中的质量矩阵，所以求解加速度并不 复杂；不必同时求解联立方程。任何节点的加速度是完全取决于节点质量和作用在节 点上的合力，使得节点计算的成本非常低。 对加速度在时间上进行积分采用中心差分方法，在计算速度的变化时假
8、定加速度 为常数。应用这个速度的变化值加上前一个增量步中点的速度来确定当前增量步中点 的速度： ( t |t | ) (tt)(t) u |u |u | tt(t) (t)(t)2 22 速度对时间的积分并加上在增量步开始时的位移以确定增量步结束时的位移： u |ttu |tt|t tu | t ()( )() (t) 2 这样，在增量步开始时提供了满足动力学平衡条件的加速度。得到了加速度，在 时间上“显式地”前推速度和位移。所谓“显式”是指在增量步结束时的状态仅依赖 于该增量步开始时的位移、速度和加速度。这种方法精确地积分常值的加速度。为了 9-3 使该方法产生精确的结果，时间增量必须相当小
9、，这样在增量步中加速度几乎为常数。 由于时间增量步必须很小，一个典型的分析需要成千上万个增量步。幸运的是，因为 不必同时求解联立方程组，所以每一个增量步的计算成本很低。大部分的计算成本消 耗在单元的计算上，以此确定作用在节点上的单元内力。单元的计算包括确定单元应 变和应用材料本构关系（单元刚度）确定单元应力，从而进一步地计算内力。 这里给出了显式动力学方法的总结： 1. 节点计算 a. 动力学平衡方程 1 u(t)(M )(P(t)I(t) b. 对时间显式积分 ( tt) (tt )(t) uuu ttt (t)(t) 22 2 uutu (tt)(t)(tt)t (t) 2 2. 单元计算
10、 a. 根据应变速率，计算单元应变增量d b. 根据本构关系计算应力 (tt )f ( (t),d ) c. 集成节点内力I(tt) 3. 设置时间t 为tt，返回到步骤1。 9.2.2比较隐式和显式时间积分程序 对于隐式和显式时间积分程序，都是以所施加的外力P、单元内力I 和节点加速 度的形式定义平衡： 9-4 Mu = P - I 其中 M 是质量矩阵。 两个程序求解节点加速度，并应用同样的单元计算以获得单元内 力。两个程序之间最大的不同在于求解节点加速度的方式上。在隐式程序中，通过直 接求解的方法求解一组线性方程组，与应用显式方法节点计算的相对较低成本比较， 求解这组方程组的计算成本要高
11、得多。 在完全 Newton 迭代求解方法的基础上，ABAQUS/Standard 使用自动增量步。在 时刻tt增量步结束时，Newton 方法寻求满足动力学平衡方程，并计算出同一时刻 的位移。 由于隐式算法是无条件稳定的，所以时间增量t比应用于显式方法的时间增 量相对地大一些。对于非线性问题，每一个典型的增量步需要经过几次迭代才能获得 满足给定容许误差的解答。每次 Newton 迭代都会得到对于位移增量u的修正值c j。 j 每次迭代需要求解的一组瞬时方程为 K ? jcPIMu jjjj j 对于较大的模型，这是一个昂贵的计算过程。有效刚度矩阵K ?j是关于本次迭代的切 向刚度矩阵和质量矩
12、阵的线性组合。直到一些量满足了给定的容许误差才结束迭代， 如力残差、 位移修正值等。 对于一个光滑的非线性响应，Newton 方法以二次速率收敛， 描述如下： 迭代相对误差 11 -2 210 -4 310 . . . 然而，如果模型包含高度的非连续过程，如接触和滑动摩擦，则有可能失去二次 收敛，而是可能需要大量的迭代过程。为了满足平衡条件，减小时间增量的值可能是 必要的。在极端情况下，在隐式分析中的求解时间增量值可能与在显式分析中的典型 9-5 稳定时间增量值在同一量级上，但是仍然承担着隐式迭代的高昂求解成本。在某些情 况下，应用隐式方法甚至可能不会收敛。 在隐式分析中，每一次迭代都需要求解
13、大型的线性方程组，这一过程需要占用相 当数量的计算资源、磁盘空间和内存。对于大型问题，对这些方程求解器的需求优于 对单元和材料的计算的需求，对于在ABAQUS/Explicit中的分析这是类似的。随着问 题尺度的增加，对方程求解器的需求迅速地增加，因此在实践中，隐式分析的最大尺 度常常取决于给定计算机中的磁盘空间的大小和可用内存的数量，而不是取决于需要 的计算时间。 9.2.3显式时间积分方法的优越性 显式方法特别地适用于求解高速动力学事件，它需要许多小的时间增量来获得高 精度的解答。如果事件持续的时间是非常短，则可能得到高效率的解答。 在显式方法中可以很容易地模拟接触条件和其它一些极度不连续
14、的情况，并且能 够一个节点一个节点地求解而不必迭代。为了平衡在接触时的外力和内力，可以调整 节点加速度。 显式方法最显著的特点是没有在隐式方法中所需要的整体切线刚度矩阵。由于是 显式地前推模型的状态，所以不需要迭代和收敛准则。 9.3自动时间增量和稳定性 稳定性限制了ABAQUS/Explicit求解器所能采用的最大时间步长，这是应用 ABAQUS/Explicit进行计算的一个重要因素。下面一节将描述稳定性限制并讨论在 ABAQUS/Explicit中如何确定这个值，还将讨论影响稳定性限制的有关模型设计参数 的问题，这些模型参数包括模型的质量、材料和网格剖分。 9.3.1显式方法的条件稳定性
15、 应用显式方法，基于在增量步开始时刻t 的模型状态，通过时间增量t前推到当 前时刻的模型状态。这个使得状态能够前推并仍能够保持对问题的精确描述的时间是 非常短的。如果时间增量是大于这个最大的时间步长，则此时间增量已经超出了稳定 9-6 性限制（stability limite ）。超过稳定性限制的一个可能后果就是数值不稳定，它可能导 致解答不收敛。由于一般不可能精确地确定稳定性限制，因而采用保守的估计值。因 为稳定性限制对可靠性和精确性有很大的影响，所以必须一致性和保守地确定这个值。 为了提高计算的效率，ABAQUS/Explicit选择时间增量，使其尽可能地接近而且又不 超过稳定性限制。 9
16、.3.2稳定性限制的定义 以在系统中的最高频率（ max）的形式定义稳定性限制。无阻尼的稳定性限制由 下式定义 t stable 2 max 而有阻尼的稳定性限制由下面的表达式定义 t stable 22 () max 式中，是最高频率模态的临界阻尼部分。（回顾临界阻尼，它定义了在自由的和有阻 尼的振动关系中在有振荡运动与无振荡运动之间的限制。为了控制高频振荡， ABAQUS/Explicit总是以体积粘性的形式引入一个小量的阻尼。）这也许与工程上的直 觉相反，阻尼通常是减小稳定性限制的。 在系统中的实际最高频率是基于一组复杂的相互作用因素，而且是不大可能计算 出确切的值。代替的办法是应用一个
17、有效的和保守的简单估算。我们不是考虑模型整 体，而是估算在模型中每个个体单元的最高频率，它总是与膨胀模态有关。可以证明， 由逐个单元为基础确定的最高单元频率总是高于有限元组合模型的最高频率。 基于逐个单元的估算，稳定极限可以用单元长度 e L和材料波速 c重新定义： d t stable e L c d 9-7 因为没有明确如何确定单元的长度，对于大多数单元类型，例如一个扭曲的四边形单 元，上述方程只是关于实际的逐个单元稳定极限的估算。作为近似值，可以采用最短 的单元尺寸，但是估算的结果并不一定是保守的。单元长度越短，稳定极限越小。波 速是材料的一个特性。对于泊松比为零的线弹性材料 c d E
18、 其中，E 是杨氏模量，是密度。材料的刚度越大，波速越高，导致越小的稳定极限； 密度越高，波速越低，导致越大的稳定极限。 这种简单的稳定极限定义提供了某些直觉上的理解。稳定极限是当膨胀波通过由 单元特征长度定义的距离时所需要的时间。如果我们知道最小的单元尺寸和材料的波 速，我们就能够估算稳定极限。例如， 如果最小单元尺寸是5 mm，和膨胀波速是5000 m/s，稳定的时间增量就是在1 10 -6s的量级上。 9.3.3在 ABAQUS/Explicit中的完全自动时间增量与固定时间增量 在分析的过程中，ABAQUS/Explicit应用在前一节讨论过的那些方程调整时间增 量的值，使得基于模型的
19、当前状态的稳定极限永不越界。时间增量是自动的，并不需 用户干涉，甚至不需要建议初始的时间增量。稳定极限是从数值模型得来的一个数学 概念。因为有限元程序包含了所有的相关细节，所以能够确定出一个有效的和保守的 稳定极限。 然而，ABAQUS/Explicit容许用户不必顾及自动时间增量。在第 9.7 节“摘 要”中简要地讨论了人工时间增量控制。 在显式分析中所采用的时间增量必须小于中心差分算子的稳定极限。如果未能使 用足够小的时间增量则会导致不稳定的解答。当解答成为不稳定时，求解变量（如位 移）的时间历史响应一般会随着振幅的增加而振荡。总体的能量平衡也将发生显著的 变化。如果模型只包含一种材料，则
20、初始时间增量是直接与网格中的最小单元尺寸成 正比。如果网格中包含了均匀尺寸的单元但是却包含有多种材料，那么具有最大波速 的单元将决定初始的时间增量。 在具有大变形和/或非线性材料响应的非线性问题中，模型的最高频率将连续地变 化，并因而导致稳定极限的变化。对于时间增量的控制，ABAQUS/Explicit有两种方 案：完全的自动时间增量（程序中考虑了稳定极限的变化）和固定的时间增量。 9-8 应用两种估算方法确定稳定极限：逐个单元法和整体法。在分析开始时总是使用 逐个单元估算法，并在一定的条件下转变为整体估算法。 逐个单元估算法是保守的；与基于整体模型最高频率的真正的稳定极限相比较， 它将给出一
21、个更小的稳定时间增量。一般说来，约束（如边界条件）和动力学接触具 有压缩特征值响应谱的效果，而逐个单元估算法没有考虑这种效果。 另一方面，整体估算法应用当前的膨胀波波速确定整个模型的最高阶频率。这种 算法为了得到最高频率将连续地更新估算值。整体估算法一般地将允许时间增量超出 逐个单元估算法得到的值。 在 ABAQUS/Explicit中也提供了固定时间增量算法。确定固定时间增量的值或者 采用在分析步中初始的逐个单元稳定性估算法，或者采用由用户直接指定的时间增量。 当要求更精确地表达问题的高阶模态响应时，固定时间增量算法可能是更有用的。在 这种情况下，可能采用比逐个单元估算法更小的时间增量值。当
22、在分析步中应用了固 定时间增量，ABAQUS/Explicit将不再检查计算的响应是否稳定。通过仔细地检查能 量历史和其他的响应变量，用户应当确保得到了有效的响应。 9.3.4质量缩放以控制时间增量 由于质量密度影响稳定极限，在某些情况下，缩放质量密度能够潜在地提高分析 的效率。例如，许多模型需要复杂的离散，因此有些区域常常包含着控制稳定极限的 非常小或者形状极差的单元。这些控制单元常常数量很少并且可能只存在于局部区域。 通过仅增加这些控制单元的质量，就可以显著地增加稳定极限，而对模型的整体动力 学行为的影响是可以忽略的。 在 ABAQUS/Explicit中的自动质量缩放功能，可以阻止这些有
23、缺陷的单元不影响 稳定极限。质量缩放可以采用两种基本方法：直接地定义一个缩放因子或者给那些质 量需要缩放的单元逐个地定义所需要的稳定时间增量。这两种方法都容许对稳定极限 附加用户控制，详细介绍请参考ABAQUS 分析用户手册第7.15.1 节“Mass scaling” 。 然而，当采用质量缩放时也要小心，因为模型质量的显著变化可能会改变问题的物理 模型。 9-9 9.3.5材料对稳定极限的影响 材料模型通过它对膨胀制波波速的限制作用来影响稳定极限。在线性材料中，波 速是常数；所以，在分析过程中稳定极限的唯一变化来自于最小单元尺寸的变化。在 非线性材料中，例如产生塑性的金属材料，当材料屈服和材
24、料的刚度变化时波速发生 变化。在整个分析过程中，ABAQUS/Explicit监督在模型中材料的有效波速，并应用 在每个单元中的当前材料状态估算稳定性。在屈服之后刚度下降，减小了波速并因而 相应地增加了稳定极限。 9.3.6网格对稳定极限的影响 因为稳定极限大致与最短的单元尺寸成比例，所以应该优先使单元的尺寸尽可能 的大。遗憾的是，对于精确的分析采用一个细划的网格常常是必要的。为了在满足网 格精度水平要求的前提下，尽可能地获得最高的稳定极限，最好的方法是采用一个尽 可能均匀的网格。由于稳定极限是基于在模型中最小的单元尺寸，甚至一个单独的微 小单元或者形状极差的单元都能够迅速地降低稳定极限。为了
25、便于发现问题， ABAQUS/Explicit在状态文件（ .sta）中提供了网格中具有最低稳定极限的10 个单 元的清单。如果在模型中包含了一些稳定极限比网格中其它单元小得多的单元，将模 型网格重新划分使其更加均匀可能是值得的。 9.3.7数值不稳定性 在大多数情况下，ABAQUS/Explicit对于大多数单元保持了稳定。但是，如果定 义了弹簧和减振器单元，它们在分析过程中有可能成为不稳定。因此，能够在你的分 析过程中识别是否发生了数值不稳定性是非常有用的。如果确实发生了数值不稳定， 典型的情况是结果是无界的，没有物理意义的，而且常常解是振荡的。 9.4例题：在棒中的应力波传播 本例题展示
26、了在前面第2 章“ABAQUS 基础” 中所描述过的显式动态分析的一些 9-10 基本思想。它也描述了稳定极限，以及在求解时网格细划和材料的影响。 棒的尺寸如图9-1 所示。 主视图侧视图 图 9-1在棒中波传播的问题描述 为了使问题成为一个一维的应变问题，所有四个侧面均由滚轴支撑；这样，三维 模型模拟了一个一维问题。材料为钢材，其性质如图9-1 所示。棒的自由端承受一个 量级为 1.0 105Pa的爆炸载荷，如图9-2 所示，爆炸载荷的持续时间为3.88 10-5s。 幅 振 时间 图 9-2爆炸载荷的幅值时间曲线 9-11 9.4.1前处理用 ABAQUS/CAE创建模型 在这一节中，我们
27、将讨论如何应用ABAQUS/CAE创建这个模拟所用的模型。在 本手册的在线文档第A.7 节“Stress wave propagation inabar”提供了输入文件。当通 过 ABAQUS/CAE运行这个输入文件时，将创建关于该问题的完整的分析模型。根据 下面给出的指导如果你遇到困难，或者如果你希望检查你的工作，则可以运行这个输 入文件。在附录A“Example Files”中，给出了如何提取和运行输入文件的指导。 如果你没有进入ABAQUS/CAE或者其它的前处理器，可以人工创建关于这个问 题的输入文件， 关于这方面的讨论， 见 GettingStarted withABAQUS/Sta
28、ndard ： Keywords Version，第 3.4 节“Example：stresswave propagation inabar” 。 定义模型几何 在这个例子中， 应用可拉伸实体的基本特征，将创建一个三维的可变形物体。 首先画一个棒的二维轮廓图然后将它拉伸成型。 创建部件： 1在 CreatePart 对话框中，创建一个部件并命名为Bar ，接受三维的变形体和 可拉伸实体的基本特征的默认设置，对于模型采用近似的尺寸为0.50 。 2应用在图9-3 中给出的尺寸画棒的横截面。 图 9-3矩形 9-12 可以采用如下的步骤： a. 应用位于画图工具箱右上角的CreateLines:C
29、onnected工具创建一个 0.20m 高 0.20m 宽的矩形。 b. 当完成绘制轮廓图后，在提示区点击Done 。 显示 Edit Base Extrusion（编辑基础拉伸） 对话框。 为了完成部件定义， 你必须指定横截面拉伸的距离。 c. 在对话框中，输入拉伸深度1.0m。 3将模型保存到名为Bar.cae的模型数据库文件中。 定义材料和截面性质 创建一个单一线弹性材料，命名为Steel，采用密度7800kg/m 3，杨氏模 量为 207E9Pa和泊松比0.3 。 创建一个均匀的实体截面定义，命名为 BarSection，接受 Steel 作为材料， 接受 Plane stress/
30、strain thickness为 1。 将截面定义BarSection赋予整个部件。 创建装配件 进入 Assembly 模块，并创建一个部件Bar 的实体。模型按照默认方向放置， 整体的 3 轴位于棒的长度方向。 创建几何集合和面 创建几何集合TOP、BOT、FRONT、BACK、FIX 和 OUT，如图 9-4 所示（集 合 OUT 包含楞边， 在图 9-4 中如黑粗线所示） 。创建面命名LOAD，如图 9-5 所示。 这些区域将用于施加载荷和边界条件，以及定义需要的输出变量。 9-13 图 9-4集合 图 9-5面 定义分析步 创建一个单一的动态、显式分析步，命名为BlastLoad。
31、键入Apply pressureloadpulse作为分析步的描述，并设置Timeperiod为 2.0E-4s。 在 EditStep对话框中，点击Other页。为了保持应力波尽可能的尖锐，将 Quadrastic bulk viscosityparameter（二次体积粘性参数，将在第 9.5.1节“体 积粘性”中讨论）设置为0。 设置输出要求 编辑默认的场变量输出要求，这样在分析步BlastLoad 中，将预先选择的场 变量数据以四个相等的空间间隔写入输出数据库。 9-14 删除已存在的默认的历史变量输出请求，而创建一个新的历史变量输出请求 的集合。在CreateHistoryOutpu
32、t （创建历史变量输出）对话框中，接受默认 的名称H-Output-1和选择的分析步BlastLoad，点击Continue。点击在 Domain （范围）选项框旁边的箭头，选择 Set name （集合名称） ，然后选择OUT。 在 OutputVarables （输出变量）列表中，点击在Stresses左边的三角形，点 击在 S, Stresscomponents and invariants（应力分量与不变量）左边的三角 形，并选中S33 变量，它是在棒的轴向的应力分量，指定在每1.0E-6s 保存一 次输出。 定义边界条件 创建一个边界条件，命名为Fix right end，并在所有三
33、个方向上约束棒 的右端面（几何集合FIX ）（见图9-1）。创建其他的边界条件，在这些面的法线 方向约束顶面、底面、前面和后面（集合FRONT和 BACK 为 1 方向，集合TOP 和 BOT为 2 方向）。 定义载荷历史 -5s。 爆炸载荷将以它的最大值瞬时地施加并保持为常数，持续时间为3.88 10 然后载荷突然全部去除并保持为零值。创建一个幅值定义，命名为Blast，采用 在图 9-6 中所示的数据。在本问题中在任意给定时刻的压力载荷值是指压力载荷 的给定量级乘以由幅值曲线插值的值。 创建压力载荷，命名为Blastload ，并选择BlastLoad作为载荷施加的 分析步。 将载荷施加在
34、LOAD面上。选择 Uniform （均匀） 分布，指定值为1.0E5 Pa 作为载荷量级，并选择幅值为Blast 。 9-15 图 9-6定义爆炸载荷幅值的数据 创建网格 利用材料性质 （忽略了泊松比） ，我们可以应用前面介绍的公式计算材料的波 速 c d E 9 20710 MPa 3 7800kg /m 3 5.15 10 m/s 我们感兴趣的是随着时间应力沿着棒长度方向的传播；所以，需要一个足够 精细的网格来精确捕捉应力波。看起来使爆炸载荷发生在10个单元的跨度内是适 -5 合的。因为爆炸是持续了3.88 10 s，这意味着我们希望爆炸持续时间乘以波速 等于 10 个单元的长度： 5
35、L10el(3.8810s)cd 9-16 -4s时通过棒的固定端。 10 个单元的长度为0.2 m。 波以这个速度在1.94 10 因为棒的长度为1.0 m，这意味着我们要在长度方向上划分50 个单元。为了保持 网格均匀，在每个横向上也划分了10 个单元，使得网格为50 10 10，这个网 格如图 9-7 所示。 图 9-750 10 10 网格 我们使用整体单元尺寸0.02 为播撒种子的目标。 选择 C3D8R作为单元类型， 并剖分网格。 创建、运行和监控作业 创建一个作业，命名为Bar ，并键入Stresswave propagation ina bar(SI units)作为作业的描述
36、。提交作业，并监控分析结果。如果遇到了 任何错误，必须修改模型和重新运行模拟。必须调查任何警告信息的来源和采取 适当的措施；回顾到某些警告信息可以安全地忽略，而其他的警告信息需要采取 纠正的措施。 状态文件（ .sta） 你也可以观察状态文件Bar.sta来监控作业的进程，其中的信息包括关于 9-17 惯性矩， 接着是关注稳定极限的信息。按照顺序列出了10 个具有最低稳定时间极 限的单元。 在状态文件中继续给出求解过程的信息。下面的信息也显示在JobMonitor 中。 9-18 9.4.2后处理 在 JobManager ，点击 Results进入 ABAQUS/CAE的 Visulizat
37、ion 模块，并自 动地打开由这个作业创建的输出数据库（.odb ）文件。另一种方法，从位于工具栏下 面的 Module列表中， 选择 Visulization进入 Visulization 模块； 从主菜单栏中， 通过 选择 File -Open 打开.odb文件并双点击合适的文件。 沿路径（ path ）绘制应力 我们希望观察沿着棒长度方向的应力分布是如何随着时间变化的。为此，我 们将观察在整个分析过程中的三个不同时刻的应力分布。 对于输出数据库文件的前三个框图的每一个图，创建一条沿着棒的中心线3 方向应力（ S33）变化的曲线。为了创建这些绘图，你首先需要定义沿着棒的中 心的直线路径。
38、沿着棒的中心创建一条由点构成的路径（point list path） ： 1在主菜单栏中，选择Tools -Path -Create 。 显示 CreatePath（创建路径）对话框。 9-19 2 命名路径为Center，选择 Point list （点列）作为路径类型， 并点击 Continue。 显示 Edit Point List Path（编辑点列路径）对话框。 3在PointCoordinates（点坐标）列表中，输入棒两端中心的坐标。例如，如 果应用我们前述的方法生成了几何和网格，那么在列表输入中是0,0, 1和 0,0,0。（这个输入指定了从（0,0,1）到（ 0,0,0）的一
39、条路径，如在模型的 整体坐标系中所定义的。） 4当完成后，点击OK 关闭 Edit Point List Path对话框。 保存在三个不同时刻沿此路径的应力的X-Y 曲线图： 1在主菜单栏中，选择Tools -XY Data -Manager 。 2在 XY DataManager（XY 数据管理器）中，点击Create 。 显示 CreateXY Data（创建 XY 数据）对话框。 3选择 Path （路径）作为XY 数据的来源，并点击Continue 。 显示 XY Datafrom Path（从路径中获取XY 数据）对话框，以及你已经创建 的在路径列表中可以找到的路径。如果当前显示的是
40、未变形的模型形状，在视 图中高亮度显示你所选择的路径。 4在 Point loaction（点位置），选中Include intersection（包括交叉点）。 5在对话框的XValues （X 值）部分中，接受True distance（真实距离）作为 选择。 6在对话框的Y Values （Y 值）部分中，点击Field Output（场变量输出）以打 开 FieldOutput对话框。 7选择 S33应力分量，并点击OK 。 在 XY Datafrom Path对话框中的场输出变量发生变化，表示将创建在3 方向 的应力数据。 注意 ：ABAQUS/CAE可能警告你场输出变量将不会影响当
41、前的图像，保留绘 图模式为Asis ，并点击 OK 继续。 8在 XY Datafrom Path对话框中的Y Values部分，点击Step/Frame。 9在弹出的Step/Frame对话框中， 选择 frame 1，它是 5个记录框图的第2个图。 （列出的第1 个框图为 frame0，它是模型在分析步开始时的基本状态。）点击 OK。 在 XY Datafrom Path对话框中的Y Values部分发生改变，表示将从第一个 9-20 分析步的frame 1 创建数据。 10保存 XY 数据，点击Save as 。 显示 Save XY Dataas对话框。 11命名 XY 数据为 S33
42、_T1 ，并点击 OK。 在 XY Data Manager中，显示出S33_T1 。 12重复步骤8到步骤 10，创建 frame 2和 frame 3 的 XY 数据，并分别命名数 据集合为S33_T2 和 S33_T3 。 13关闭 XY Data from Path对话框，点击Cancel 。 绘制应力曲线： 1在 XY Data Manager对话框中，拖动光标高亮度显示所有3 组 XY 数据集。 2点击 Plot 。 ABAQUS/CAE绘制出沿着棒中心3 方向上对应于frame 1、2 和 3 的应力， 它们 -5s、1 10-4s 和 1.5 10-4s。 对 应于近似的模拟时
43、刻分别为5 10 设置 XY 曲线图： 1从主菜单栏中，选择Options -XY Plot 。 显示 XY Plot Options（XY 图选项）对话框。 2点击 TickMarks（刻度）页。 使 TickMarks选项可以工作。 3指定 Y 轴的主要刻度出现在20E3s增量（ increments）。 4对于 X 轴和 Y 轴的次要刻度选项，在每个主要刻度间隔之间指定次要刻度为0。 你也可以设置每个轴的标题。 5点击 Titles （标题）页，使标题选项可以工作。 6在 X-Axis（X 轴）域，选择标题来源为User-specified（自定义） 。在 Title Text （标题内
44、容）域，输入Distancealong bar(m)。 7在 Y-Axis 域，指定Stress-S33 (Pa)为 Y 轴的标题。 8点击 OK 确认用户选择的XY 绘图参数，并关闭XY Plot Options对话框。 设置在 XY 绘图中曲线的显示： 1从主菜单栏中，选择Options -XY Curve 。 显示 XY Curve Options（XY 曲线选项）对话框。 2在 XY Data 数据域中，选择S33_T2 。 9-21 3对于 S33_T2 曲线，选择点线类型，并点击Apply 。 S33_T2 曲线变成为点线。 4重复步骤2 与 3，使 S33_T3成为虚线。 5点击
45、 Dismiss关闭 XY Curve Options对话框。 所设置的绘图显示在图9-8 中。 图 9-8在 3 个不同时刻沿着棒的应力（S33） 我们能够看到，在三条曲线的每一条中应力波在棒的长度上的影响是近似为0.2 m。这个距离应该对应于爆炸波在作用时间内传播的距离，这可以通过简单的计算来 3m/s，那么波传播 0.2 m 所用的时间 验 证。如果波前的长度为0.2 m 和波速为 5.15 10 -5 为 3.88 10 m/s。正如所预料的，这就是我们所施加的爆炸荷载的作用时间。当应力 波沿着棒传播时它并不是严格的方波。特别是在应力突然改变之后有回复或者摆动。 在本章后面将要讨论的线
46、性体粘性，减缓了这种回复，因此并未对结果有负面的影响。 创建历史曲线图 另一种研究结果的方法是观察在棒中的三个不同点的应力的时间历史；例如，距 离棒的加载端为0.25m、0.50 m 和 0.75 m 的三个点。 为此， 我们必须首先确定位于这 9-22 些位置处的单元编号。确定这些单元编号的一种容易的方法是在包含沿着棒边界的单 元（集合OUT）的显示组中查询这些单元。 创建和绘出显示组并查询单元编号： 1从主菜单栏中，选择Tools -DisplayGroup -Create 。ABAQUS/CAE显示 Create Display Group（创建显示组）对话框。 2选择 Elements
47、作为 Item ，和 ElementsSets（单元集） 作为 SeletionMethod （选择方式）。从几何集合列表中，选择OUT，点击 Save SelectionAs（保 存选择内容为）。 3在 Save Selection As对话框中，命名显示组为History plot。点击 OK 。 4点击 Dismiss关闭 Create Display Group对话框。 5从主菜单栏中，选择Plot -Undeformed Shape，绘制未变形形状。 6从主菜单栏中，选择Tools -DisplayGroup -Plot -Historyplot ，绘制所 创建的显示组。 7从主菜单
48、栏中，选择Tools -Query 。 8在弹出的Query （查询）对话框中，选择Probe Values（探测），并点击OK。 显示 Probe Values对话框。 9点击在图9-9 中的阴影单元（在棒中的每第13 个单元）。单元的ID（编号）显 示在 Probe Values对话框中。标记这三个阴影单元的编号。 10点击 Cancel关闭 Probe Values对话框。 图 9-9History plot显示组 9-23 当提示你是否将结果写入到一个文件时，点击No。 绘制应力历史： 1从主菜单栏中，选择Result -History Output。 ABAQUS/CAE显示 His
49、tory Output对话框。 在 Output Varialbes域中包含了 在输出数据库的历史变量部分中的所有变量的列表；这些也是你能够绘制的所 有变量。为了观察变量选择的完整描述，拖动对话框的左边或右边框，增加 History Output对话框的宽度。 2应用 Ctrl+ 点击， 选择多组XY 数据集合，对于已经标识的三个单元（每第 13 个单元），选择在3 方向上的应力（S33）数据。 3在 History Output对话框的底部，点击Plot 。 ABAQUS/CAE绘制出在每个单元中的应力（纵轴）随时间变化的XY 图。 4点击 Dismiss关闭对话框。 如前所述，你可以设置图
50、的显示。 设置 X-Y 图： 1从主菜单栏中，选择Options - XY Plot 。 显示 XY Plot Options对话框。 2点击 Title 页。 标题选项可以工作。 3在 X-axis 域，指定X 轴标题为Total time(s)。 4点击 OK 确认所设置的XY 曲线图选项，并关闭对话框。 设置在 X-Y 图中曲线的显示： 1从主菜单栏中，选择Options - XY Curve 。 显示 XY Curve Options对话框。 2在 XY Data域中，选择对应于最接近于棒自由端的单元的临时的XY 数据编 号。（在这个集合中的单元最先受到应力波的影响。） 3选择 Use
