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

《ABAQUS中冲击动力学问题的求解方法.docx》由会员分享，可在线阅读，更多相关《ABAQUS中冲击动力学问题的求解方法.docx（8页珍藏版）》请在优知文库上搜索。

1、ABAQUS中冲击动力学问题的求解方法冲击载荷随时间迅速变化。当物体的局部位置受到冲击时，所产生的扰动会逐渐传到未扰动的区域去，这种现象称为应力波的传播。当载荷作用时间短、变化快，且受力物体在加载方向的尺寸又足够大时，这种应力波的传播就显得特别重要网。研究动力学问题最终将简化为求解动力学平衡方程式：节点质量矩阵M乘以节点加速度ii等于节点的合力(所施加的外力P与单元内力/之间的差值)：Mii=P-I(2-1)由于考虑了惯性力的影响，动力学平衡方程中出现了质量矩阵，最后得到的求解方程不是代数方程组，而是常微分方程组。1冲击动力学求解方法如果加载时间过短或者是动态载荷，需要采用动态分析(dynam

2、icanalysis)0复合材料的低速冲击就属于动态分析问题。动态分析又分为隐式分析和显式分析。在隐式分析中，结构的刚度矩阵需要进行多次生成和求逆，这使得分析求解成本大大增加，而且刚度退化和材料失效常常引起计算收敛问题。在显示分析中，能够避免计算收敛，较好地求解这一问题。1.1 显式与隐式分析的区别显式与隐式分析的区别在于：显式分析需要很小的时间增量步，它仅依赖于模型的最高固有频率，而与载荷的类型和持续的时间无关。通常的模拟需要100007000000个增量步，每个增量步的计算成本相对较低。它的求解方法是在时间域中以很小的时间增量步向前推出结果，而无需在每一个增量步求解耦合的方程系统,或者生成

3、总体刚度矩阵。隐式分析对时间增量步的大小没有内在的限制，增量的大小通常取决于精度和收敛情况。典型的隐式模拟所采用的增量步数目要比显式模拟小几个数量级。然而，由于在每个增量步中必须求解一套全域的方程组，所以对于每一增量步的成本，隐式方法远高于显式方法。1.2 计算方法选择复合材料层合板低速冲击损伤涉及到复杂的接触问题、材料刚度随着载荷发生变化的问题、材料的退化(degradation)和失效(failure)导致的严重的收敛问题,这些问题在隐式分析中都无法实现或者求解成本比较昂贵。综上所述，本文选取显式分析方法。13软件介绍本文的数值模拟采用商业有限元软件ABAQUS,下面对软件进行简单介绍：A

4、BAQUS是一套功能强大基于有限元方法的工程模拟软件，可以解决从相对简单的线性分析到复杂的非线性模拟等各种问题。ABAQUS具备十分丰富的单元库，可以模拟任意实际形状，也具有很丰富的材料模型库，可以模拟大多数典型工程材料的性能，包括金属、复合材料、橡胶、钢筋混凝土以及地质材料等。作为一种通用的模拟工具，ABAQUS不仅能解决结构分析（应力/位移）问题,而且能够模拟和研究热传导、电子元器件的热控制（热-电耦合）、声学、土壤力学（渗流-应力耦合分析）和压电分析等领域的问题。ABAQUS为用户提供了广泛的功能，使用起来十分方便，在许多模拟中（包括高度非线性问题），用户只需提供结构的几何形状、材料性质

5、、边界条件和载荷工况等工程数据。在非线性分析中，ABAQUS能自动选择合适载荷增量和收敛准则。并在分析中不断地调整这些参数值，确保获得精确的解答，用户几乎不必去定义任何参数就能控制问题的数值求解过程。ABAQUS由多个模块构成，包括前后处理模块ABAQUS/CAE主求解器模块ABAQUS/Standard和ABAQUS/Explicit等。ABAQUS/CAE（前后处理）ABAQUS/CAE是ABAQUS有限元分析的前后处理模块，也是建模，分析和仿真的人机交互平台。该模块根据结构的几何图形生成网格，将材料和截面的属性分配到网格上，并施加载荷和边界条件。该模块可以进一步将生成的模型提交到后处理分

6、析模块中运行，对运行情况进行监测，并对计算结果进行后处理。ABAQUS/CAE的后处理支持ABAQUS分析模块的所有功能，并且对计算结果的描述和解释提供了范围很广的选择，除了云图，等值线和动画显示之外，还可以用列表，曲线等其他常用工具的来完成工程显示。ABAQUSZStandard（隐式求解器）ABAQUSZStandard是一个通用分析模块，它能够求解广泛领域的线性和非线性问题，包括静态分析、动态分析，以及复杂的非线性耦合物理场分析等。在每,个求解增量步（increment）中，ABAQUS/Standard隐式地求解方程组。ABAQUS/Explicit（显式求解器）使用ABAQUS/Ex

7、plicit可以进行显式动态分析,它适于求解复杂非线性动力学问题和准静态问题，特别是用于模拟短暂、瞬时的动态事件，如冲击和爆炸问题。此外，它对处理接触条件变化的高度非线性问题也非常有效，例如模拟成型问题，它的求解方法是在时间域中以很小的时间增量步向前推出结果，而无需在每个增量步求解耦合的方程系统，或者生成总体刚度矩阵。ABAQUS/Explicit不但支持应力/位移分析，而且还支持完全耦合的瞬态温度/位移分析、声固耦合分析。任意的拉格朗日一欧拉自适应网格功能可以有效地模拟大变形非线性问题。将ABAQUSZStandard和ABAQUS/Explicit结合使用，结合二者的隐式和显式求解技术，可

8、以求解更广泛的实际问题。综上所述，本文应用显示求解器ABAQUS/Explicit进行数值模拟分析。2动力学显式有限元方法ABAQUS/Explicit是基于显式算法的有限元程序。显式有限元程序可以分析材料由于力学性能退化而引起材料切线刚度矩阵出现非正定的问题。因此ABAQUS/Explicit程序可以分析纤维增强树脂基复合材料结构受冲击这高度非线性问题。使用显式有限元分析时，根据材料的本构模型来确定单元的应力、应变和应变率等。一般情况下，ABAQUS/Explicit求解引擎将每个积分点的应变增量和应变率传递给材料本构，然后再根据相应的材料积分点的信息确定当前状态的应力。对于线弹性材料的计算

9、很简单，当前的总应变与材料刚度矩阵的乘积就是当前状态的应力；而对于非线性材料本构(如塑性材料)的计算就很复杂，本构模型模块将计算得到的应力再传递给求解引擎，由其计算出节点的加速度。然后，根据系统的运动方程求解出节点速度和位移，接着确定下一时间步的应变和应变率。2.1显式时间积分ABAQUS/Explicit应用中心差分方法对运动方程进行显式的时间积分，应用一个增量步的动力学条件计算下一个增量步的动力学条件。在增量步开始时，程序求解动力学平衡方程式：节点质量矩阵M乘以节点加速度U等于节点的合力(所施加的外力P与单元内力/之间的差值)，即式(2-1)：M=P-/在当前增量步开始时(t时刻)，计算加

10、速度：ii=(M)xP-I)(2-2)由于显式算法总是采用一个对角的或者集中的质量矩阵，所以求解加速度并不复杂，不必同时求解联立方程。采用节点的加速度完全取决于节点质量和作用在节点上的合力，使得节点计算的成本非常低。对加速度在时间上进行积分采用中心差分法，在计算速度的变化时假定加速度为常数。应用这个速度的变化值加上前个增量步中点的速度来确定当前增量步中点的速度：(2-3).I.I(AflC+&)+&W.wl/=w,+I(/H-(/)2220速度对时间的积分加上在增量步开始时的位移可以确定增量步结束时的位移：uIs=+&swIa(2-4)2这样,在增量步开始时提供了满足动力学平衡条件的加速度。得

11、到了加速度,在时间上“显式地”前推速度和位移。所谓“显式”是指在增量步结束时的状态仅依赖于该增量步开始时的位移、速度和加速度。这种方法精确地积分常值的加速度。为了使该方法得到精确的结果，时间增量必须相非常小，这样在增量步中加速度几乎为常数。由于时间增量步必须很小，所以一个典型的分析需要成千上万个增量步。因为不必同时求解联立方程组，所以每一个增量步的计算成本很低。大部分计算成本消耗在单元计算上，以此确定作用在节点上的单元内力。单元的计算包括确定单元应变和应用材料本构关系(单元刚度)确定单元应力，从而进一步计算内力。下面给出显式动力学方法的总结：(1)节点计算1)动力学平衡方程%=(M)7(%-/

12、)(2-5)2)对时间显式积分w(r+=%+%(2-7)(2)单元计算1)根据应变速率C，计算单元应变增量d2)根据本构关系计算应力。Sr+)=/(S”，d)(2-8)3)集成节点内力/.+&)(3)设置时间，为f+r,返回到步骤(1)。2.2 比较隐式和显式时间积分程序对于隐式和显式积分程序，都是以所施加的外力p、单元内力/和节点加速度形式定义平衡：Mii=P-I(2-9)其中，M是质量矩阵。两个程序求解节点加速度，并应用同样的单元计算以获得单元内力。两个程序最大的区别在于求解节点加速度上。在隐式程序中，通过直接求解的方法求解一组线性方程组，与显式方法节点计算的成本比较，求解这组方程组的计算

13、成本要高得多。在完全Newton迭代求解方法的基础上,ABAQUSZStandard使用自动增量步。在时刻Z+增量步结束时，NeWtOn方法寻求满足动力学平衡方程，并计算出同一时刻的位移。由于隐式算法是无条件稳定的，所以时间增量&比应用显式方法的时间增量相对大些。对于非线性问题，每个典型的增量步需要经过几次迭代才能获得满足给定容许误差的解答。每次Newton迭代都会得到对于位移增量町的修正值与。每次迭代需要求解的一组瞬时方程为KjCj=Pj-Ij-MjUj(2-10)对于较大的模型，这是一个非常大的计算过程。有效刚度矩阵M是关于本次迭代的切向刚度矩阵和质量矩阵的线性组合。直到一些变量满足了给定

14、的容许误差才结束迭代，如位移修正值等。在隐式分析中，每一次迭代都需要求解大型的线性方程组，这一过程需要相当数量的计算资源、内存和磁盘空间。2.3 显式时间积分方法的优越性显式方法特别适用于求解高速动力学问题，它需要许多很小的时间增量来获得高精度的解答。如果事件持续的时间非常短，则可能得到高效率的解答。在显式方法中可以很容易的模拟接触条件和其他一些极度不连续的情况，并且能够一个节点、一个节点地求解而不必迭代。为了平衡接触时的内力和外力，可以调整节点加速度。显式方法最显著的特点是没有在隐式方法中所需要的整体切向刚度矩阵。由于是显式地前推模型的状态，所以不需要迭代和收敛准则。2.4 显式方法的条件稳

15、定性稳定性限制了ABAQUS/Explicit求解器所能采用的最大时间步长，这是ABAQuSZExplicit进行计算的一个重要因素。应用显式方法，基于在增量步开始时刻f的模型状态，通过时间增量/前推到当前时刻的模型状态。使得状态能够前推并仍能保持对问题的精确描述的时间非常短。如果时间增量大于这个最大的时间步长，则时间增量已经超出了稳定性限制(stabilitylimit)o超过稳定性限制的后果是数值不稳定，可能导致计算不收敛。由于一般不可能精确地确定稳定性限制，因而采用保守的估计值。为了提高计算效率，ABAQUS/Explicit选择时间增量，使其尽可能地接近而且又不超过稳定性限制。在系统中，以最高频率(wt)的形式定义稳定性限制。无阻尼的稳定极限由下式定义：Atablc=(2-11)有阻尼的max稳定极限由下面的表达式定义：e=(4-14)(2-12)式中，4是最高频率模态的临界阻尼部分。由此可以看出，阻尼通常是减小稳定性限制的。在系统中的实际最高频率基于一组复杂的相互作用因素，而且不可能计算出确切的值。代替的方法是应用一个有效的、保守的简单估算。不是考虑模型整体，而是估算在模型中每个个体单元的最高频率，它总是与膨胀波的波速有关。可以证明，以逐个单元为基础确定的最高单元频率总是高于有限元组合模型的最高频率。基于逐个单元的估算，稳定极限可以