LS-DYNA显式-隐式计算.docVIP

在进行显式动力分析之前，很多结构由于各种原因而存在初始应力，如果在工程上不能确定初始应力对动力分析的影响程度，一般需考虑初始应力的影响，如旋转机械、压力容器、螺栓联接及加工预载荷部件等。 在LS-DYNA中，有两种方式来施加初试应力： 一种是运ANSYS的隐式分析软件求解初始应力，通过文件的传送输入到LS-DYNA显式分析中，LS-DYNA显式求解器在开始显式分析前，即时间为零的时刻通过一个虚假的时间段，101个时间步施加初始应力。 一种就是运用LS-DYNA中本身载荷的动力松弛方式来施加，即允许LS-DYNA显式求解器对线性和非线性的静态或准静态问题进行近似的求解，通过增加阻尼使动能降低为零，求解器每250个循环（缺省）检查动能状态直到为零。 现就两种方法的具体实现进行说明： 一. ANSYS软件的初始应力分析： 用到ANSYS/LS-DYNA隐式到显式序列求解功能，有8个基本步骤: 1. 进行隐式求解 2. 为进行显式求解改变当前的作业名 3. 将隐式单元转换为与之对应的具有适当属性的显式单元 (关键选项,实常数，材料特性等) 4. 移进行隐式分析时所加的附加约束 5. 将来自隐式分析的节点结果写到动力松弛文件（drelax文件）中（drelax文件中包含节点的标号、XYZ位移、XYZ角度和温度） 6. 使用动力松弛文件初始化用于显式分析模型的几何形状 7. 给显式分析施加另外的载荷条件 8. 进行显式求解 当在ANSYS/LS-DYNA环境中进行SOLVE时，在后台调用的是： I=jobename.k m=drelax 命令，jobename.k中包含下面关键字：*CONTROL_DYNAMIC_RELAXATION250 0.100E-020.9999 0.127E+31 0.900 1 0.400E-01 2（用第二种方式求解时，2改为1） 二. LS-DYNA中本身载荷的动力松弛方式 在LS-DYNA中载荷中定义的曲线有三种应用方式： 仅用于瞬态分析 仅用于动力松弛 同时用于瞬态分析和动力松弛 相关的关键字为：*DEFINE_CURVE*CONTROL_DYNAMIC_RELAXATION250 0.100E-020.9999 0.127E+31 0.900 1 0.400E-01 1*LOAD_BODY_RX（RY、RZ） 运行后收敛的结果即为初始化应力，同时生成动力松弛文件（drdisp.sif文件） 在该方法中，若载荷仅进行动力松弛，求解时间一定要设为零，若同时用于瞬态分析和动力松弛。建议用两个载荷方式来进行.（仅用于瞬态分析和仅用于动力松弛） 序列求解——是一种将隐式（ANSYS）和显式（ANSYS/LS-DYNA）求解方法综合使用的分析技术。在需要序列求解的问题中，为了获得最后的结果，显式分析的结果被输入模型中（或相反）。 进行序列求解原因：有些工程问题过程十分复杂，包括动态和静态两个阶段（如：压力容器在跌落测试之前的初始环向应力或金属成形之后的线性回弹）。显式技术适合于求解非线性动态碰撞问题，不适于求解自然现象中的静力阶段。而隐式方法最适合于求解静态或准静态的问题。将两者的求解器联合使用是一个特别强大的工具。 显式－隐式序列求解是首先使用ANSYS/LS-DYAAN程序进行动力求解，然后将变形后的几何形状和应力输入ANSYS隐式分析中，通过给定的合适的边界条件进行后续求解。 ——显式－隐式求解技术当前仅能用于壳单元。——在分析显式部分，SHELL163被使用，然后显式结果被输入到SHELL181隐式单元。——为了阻止刚体运动在隐式分析中必须指定合适的约束条件。——显式－隐式序列求解的主要目的是模拟回弹。 回弹——在金属成形过程中，回弹被定义为变形部件（板料）离开模具后的空间变化，回弹是由于线弹性卸载而引起的。在成形时，当板料和模具接触时弹性能被储存在板材中，当成形压力被移走后，弹性能被释放，导致板料朝着它原始的方向变形或回弹。 执行一个显式－隐式序列求解需九个基本步骤：1，求解分析的显式部分。2，进入隐式求解器，改变当前的作业名。3，将显式单元转变成具有适当属性的隐式单元。4，关闭隐式单元形状检查。5，将隐式单元的几何形状修改为显式求解后的变形形状。6，不选择隐式求解所不需要的单元（仅保留非刚体的SHELL181单元）。7，重新定义边界条件。8，输入来自显式分析的应力（单元膜力）和厚度。9，求解分析的隐式部分。 步骤1：求解显式 ——在考虑回弹效应的成形模拟中，必须用SHELL163单元来模拟板料。——确保用于模拟板料的壳单元的厚度是真实的。——为了加速整个模拟世界增加冲头速度。——在进行隐式求解之前验证显式分析的结果。——在显式求解完成之后利用时间历程