ABAQUS中冲击动力学问题的求解方法.docxVIP

ABAQUS中冲击动力学问题的求解方法 冲击载荷随时间迅速变化。当物体的局部位置受到冲击时，所产生的扰动会逐渐传到未扰动的区域去，这种现象称为应力波的传播。当载荷作用时间短、变化快，且受力物体在加载方向的尺寸又足够大时，这种应力波的传播就显得特别重要[35]。 研究动力学问题最终将简化为求解动力学平衡方程式：节点质量矩阵M乘以节点加速度等于节点的合力（所施加的外力P与单元内力I之间的差值）： (2-1) 由于考虑了惯性力的影响，动力学平衡方程中出现了质量矩阵，最后得到的求解方程不是代数方程组，而是常微分方程组。 1 冲击动力学求解方法 如果加载时间过短或者是动态载荷，需要采用动态分析（dynamic analysis）。复合材料的低速冲击就属于动态分析问题。 动态分析又分为隐式分析和显式分析。在隐式分析中，结构的刚度矩阵需要进行多次生成和求逆，这使得分析求解成本大大增加，而且刚度退化和材料失效常常引起计算收敛问题。在显示分析中，能够避免计算收敛，较好地求解这一问题。 1.1 显式与隐式分析的区别 显式与隐式分析的区别在于[5]： 显式分析需要很小的时间增量步，它仅依赖于模型的最高固有频率，而与载荷的类型和持续的时间无关。通常的模拟需要10000~1000000个增量步，每个增量步的计算成本相对较低。它的求解方法是在时间域中以很小的时间增量步向前推出结果，而无需在每一个增量步求解耦合的方程系统，或者生成总体刚度矩阵。 隐式分析对时间增量步的大小没有内在的限制，增量的大小通常取决于精度和收敛情况。典型的隐式模拟所采用的增量步数目要比显式模拟小几个数量级。然而，由于在每个增量步中必须求解一套全域的方程组，所以对于每一增量步的成本，隐式方法远高于显式方法。 1.2 计算方法选择 复合材料层合板低速冲击损伤涉及到复杂的接触问题、材料刚度随着载荷发生变化的问题、材料的退化（degradation）和失效（failure）导致的严重的收敛问题，这些问题在隐式分析中都无法实现或者求解成本比较昂贵。综上所述，本文选取显式分析方法。 1.3 软件介绍 本文的数值模拟采用商业有限元软件ABAQUS，下面对软件进行简单介绍： ABAQUS是一套功能强大基于有限元方法的工程模拟软件，可以解决从相对简单的线性分析到复杂的非线性模拟等各种问题。ABAQUS具备十分丰富的单元库，可以模拟任意实际形状，也具有很丰富的材料模型库，可以模拟大多数典型工程材料的性能，包括金属、复合材料、橡胶、钢筋混凝土以及地质材料等。作为一种通用的模拟工具，ABAQUS不仅能解决结构分析（应力/位移）问题，而且能够模拟和研究热传导、电子元器件的热控制（热-电耦合）、声学、土壤力学（渗流-应力耦合分析）和压电分析等领域的问题。 ABAQUS为用户提供了广泛的功能，使用起来十分方便，在许多模拟中（包括高度非线性问题），用户只需提供结构的几何形状、材料性质、边界条件和载荷工况等工程数据。在非线性分析中，ABAQUS能自动选择合适载荷增量和收敛准则。并在分析中不断地调整这些参数值，确保获得精确的解答，用户几乎不必去定义任何参数就能控制问题的数值求解过程。ABAQUS由多个模块构成，包括前后处理模块ABAQUS/CAE、主求解器模块ABAQUS/ Standard和ABAQUS / Explicit等。 ABAQUS/CAE（前后处理） ABAQUS/CAE是ABAQUS有限元分析的前后处理模块，也是建模，分析和仿真的人机交互平台。该模块根据结构的几何图形生成网格，将材料和截面的属性分配到网格上，并施加载荷和边界条件。该模块可以进一步将生成的模型提交到后处理分析模块中运行，对运行情况进行监测，并对计算结果进行后处理。ABAQUS/CAE的后处理支持ABAQUS分析模块的所有功能，并且对计算结果的描述和解释提供了范围很广的选择，除了云图，等值线和动画显示之外，还可以用列表，曲线等其他常用工具的来完成工程显示。 ABAQUS/Standard（隐式求解器） ABAQUS/Standard是一个通用分析模块，它能够求解广泛领域的线性和非线性问题，包括静态分析、动态分析，以及复杂的非线性耦合物理场分析等。在每一个求解增量步（increment）中，ABAQUS/Standard隐式地求解方程组。 ABAQUS/Explicit（显式求解器） 使用ABAQUS/Explicit可以进行显式动态分析，它适于求解复杂非线性动力学问题和准静态问题，特别是用于模拟短暂、瞬时的动态事件，如冲击和爆炸问题。此外，它对处理接触条件变化的高度非线性问题也非常有效，例如模拟成型问题，它的求解方法是在时间域中以很小的时间增量步向前推出结果，而无需在每一