COMSOL_Multiphysics自定义方程办法PDE

第一章LS-DYNA简介1.1LS-DYNA发展概况DYNA程序系列最初是1976年在美国LawrenceLivermoreNationalLab.由J.O.Hallquist博士主持开发完成的，主要目的是为武器设计提供分析工具，后经1979、1981、1982、1986、1987、1988年版的功能扩充和改进，成为国际著名的非线性动力分析软件，在武器结构设计、内弹道和终点弹道、军用材料研制等方面得到了广泛的应用。1988年J.O.Hallquist创建LSTC公司，推出LS-DYNA程序系列，主要包括显式LS-DYNA2D、LS-DYNA3D、隐式LS-NIKE2D、LS-NIKE3D、热分析LS-TOPAZ2D、LS-TOPAZ3D、前后处理LS-MAZE、LS-ORION、LS-INGRID、LS-TAURUS等商用程序，进一步规范和完善DYNA的研究成果，陆续推出930版（1993年）、936版（1994年）、940版（1997年），增加了汽车安全性分析（汽车碰撞、气囊、安全带、假人）、薄板冲压成型过程模拟，以及流体与固体耦合（ALE和Euler算法）等新功能，使得LS-DYNA程序系统在国防和民用领域的应用范围进一步扩大，并建立了完备的质量保证体系。1997年LSTC公司将LS-DYNA2D、LS-DYNA3D、LS-TOPAZ2D、LS-TOPAZ3D等程序合成一个软件包，称为LS-DYNA，PC版的前后处理采用ETA公司的FEMB，新开发的后处理器为LS-POST。1996年LSTC与ANSYS公司合作推出ANSYS/LS-DYNA，大大增强了LS-DYNA的分析能力，用户可以充分利用ANSYS的前后处理和统一数据库的优点。2001年5月推出960版，它在950版基础上增加了不可压缩流体求解程序模块，并增加了一些新的材料模型和新的接触计算功能，从2001年到2003年初LSTC公司不断完善960版的新功能，2003年3月正式发布970版。并对LS-DYNA的通用后处理器LS-POST增加了前处理器的功能，2003年初在LS-POST的基础上发布了LS-PREPOST1.0版。1.2分析能力LS-DYNA程序960版是功能齐全的几何非线性（大位移、大转动和大应变）、材料非线性（140多种材料动态模型）和接触非线性（30多种接触类型）程序。它以Lagrange算法为主，兼有ALE和Euler算法；以显式求解为主，兼有隐式求解功能；以结构分析为主，兼有热分析、流体-结构耦合功能；以非线性动力分析为主，兼有静力分析功能（如动力分析前的预应力计算和薄板冲压成型后的回弹计算）；军用和民用相结合的通用结构分析非线性有限元程序。LS-DYNA具有很广泛的分析功能，可模拟许多二、三维结构的物理特性：◆非线性动力分析◆热分析◆失效分析◆裂纹扩展分析◆接触分析◆准静态分析◆欧拉场分析◆任意拉格朗日-欧拉（ALE）分析◆流体-结构相互作用分析◆实时声场分析◆多物理场耦合分析（结构、热、流体、声场等）1.3材料模型LS-DYNA程序目前有140多种金属和非金属材料模型可供选择，如弹性、弹塑性、超弹性、泡沫、玻璃、地质、土壤、混凝土、流体、复合材料、炸药及起爆燃烧、刚性及用户自定义材料，并可考虑材料失效、损伤、粘性、蠕变、与温度相关、与应变率相关等性质。1.4单元库LS-DYNA程序现有16种单元类型，有二维、三维单元，薄壳、厚壳、体、梁单元，ALE、Euler、Lagrange单元等。各类单元又有多种理论算法可供选择，具有大位移、大应变和大转动性能，单元积分采用沙漏粘性阻尼以克服零能模式，单元计算速度快，节省存储量，可以满足各种实体结构、薄壁结构和流体-固体耦合结构的有限元网格剖分的需要。960版本有SPH（SmoothedParticleHydrodynamics）算法，SPHE算法（光顺质点流体动力算法）是一种无网格Lagrange算法，最早用于模拟天体物理问题，后来发现解决其它物理问题也是非常有用的工具，如连续体结构的解体、碎裂、固体的层裂、脆性断裂等。SPH算法可以解决许多常用算法解决不了的问题，是一种非常简单方便的解决动力学问题的研究方法。由于它是无网格的，它可以用于研究变形很大的结构，具体说明参见单元章节。1.5接触分析功能LS-DYNA程序的全自动接触分析功能易于使用，功能强大。现有40多种接触类型可以求解下列接触问题：变形体对变形体的接触、变形体对刚体的接触、刚体对刚体的接触、板壳结构的单面接触(屈曲分析)、与刚性墙接触、表面与表面的固连、节点与表面的固连、壳边与壳面的固连、流体与固体的界面等，并可考虑接触表面的静动力摩擦(库伦摩擦、粘性摩擦和用户自定义摩擦模型)、热传导和固连失效等。这种技术成功地用于整车碰撞研究、乘员与柔性气囊或安全带接触的安全性分析、薄板与冲头和模具接触的金属成型、水下爆炸对结构的影响，高速弹丸对靶板的穿甲模拟计算等,关于各种接触的区别和说明具体参见第三章基本概念和第十章汽车碰撞仿真分析。1.6初始条件、载荷和约束●初始速度、初应力、初应变、初始动量（模拟脉冲载荷）；●高能炸药起爆；●节点载荷、压力载荷、体力载荷、热载荷、重力载荷；●循环约束、对称约束（带失效）、无反射边界；●给定节点运动（速度、加速度或位移）、节点约束；●铆接、焊接（点焊、对焊、角焊）；●二个刚性体之间的连接－球形连接、旋转连接、柱形连接、平面连接、万向连接、平移连接；●位移/转动之间的线性约束、壳单元边与固体单元之间的固连；●带失效的节点固连。1.7ALE和Euler算法LS-DYNA程序具有Lagrange算法和Euler算法，Lagrange算法的单元网格附着在材料上，随着材料的流动而产生单元网格的变形。但是在结构变形过于巨大时，有可能使有限元网格造成严重畸变，引起数值计算的困难，甚至程序终止运算。ALE算法和Euler算法可以克服单元严重畸变引起的数值计算困难，并实现流体-固体耦合的动态分析。ALE算法先执行一个或几个Lagrange时步计算，此时单元网格随材料流动而产生变形，然后执行ALE时步计算：（1）保持变形后的物体边界条件，对内部单元进行重分网格，网格的拓扑关系保持不变，称为SmoothStep；（2）将变形网格中的单元变量（密度、能量、应力张量等）和节点速度矢量输运到重分后的新网格中，称为AdvectionStep。用户可以选择ALE时步的开始和终止时间，以及其频率。Euler算法则是材料在一个固定的网格中流动，在LS-DYNA中只要将有关实体单元标志Euler算法，并选择输运(advection)算法就可以了。LS-DYNA还可将Euler网格与全Lagrange有限元网格方便地耦合，以处理流体与结构在各种复杂载荷条件下的相互作用问题具体说明参见第十一章流固耦合分析。物体跌落水面1.8不可压缩流场分析LS-DYNA不可压缩流求解器是960版新增加的功能，用于模拟分析瞬态、不可压、粘性流体动力学现象。求解器中采用了超级计算机的算法结构，在确保有限元算法优点的同时计算性能得到大幅度提高，从而在广泛的流体力学领域具有很强的适用性。LS-DYNA不可压缩流求解器是基于隐式时间积分、显式时间积分两种算法的。它的显式算法(一阶精度)，解耦动量守恒方程并减少了对内存的需求，但带来损失部分精度。另一方面，由于显式算法满足扩散和CFD稳定性条件，又可提高计算精度。因此，采用单点积分和沙漏稳定性的显式算法被证明在不可压缩流求解分析中是简便、高效的。LS-DYNA不可压缩流求解器的二阶精度算法，采用了恒定质量的预置算法和物质质量的校正算法，合理解耦了速度场和压力场，从而减少了计算Navier_Stokes方程对CPU和内存的需求。二阶精度算法，可用于分析流场中的涡流，而且很容易推广应用到流体力学领域中湍流现象的计算分析。水射流形成的温度场、速度场、压力场1.9隐式求解950版本后，LS-DYNA增加了隐式求解功能，用于非线性结构静动力分析，包括结构固有频率和振型计算。LS-DYNA中可以交替使用隐式求解和显式求解，进行薄板冲压成型的回弹计算、结构动力分析之前施加预应力、特征值分析等。隐式求解分析具体参见第十七章隐式求解功能。1.10热分析LS-DYNA程序有二维和三维热分析模块，进行稳态或瞬态的热分析，和热固耦合分析，可以处理热传导、对流和辐射各种热问题，在焊接、冲压、锻压及碰撞等过程中可方便的考虑热问题（如塑性能转化为热能的问题）及热应力问题等。具体应用参见第十三章热分析和热固耦合分析。1.11LS-DYNA应用领域1.12LS-DYNA软硬件要求可在各主要工作站UNIX，LINUX操作系统及各种PC机Windows98、WingdowsNT/2000等操作系统上运行，支持MPP分布式计算，为巨型问题的求解提供了很好的保证。各种平台上的数据库统一一致，不需要转换数据格式，为前后处理带来了很大的方便性。注：该章节引用的一些图片来自LSTC公司。更多内容见附件免责声明：本页面/内容部分素材来源于互联网公开信息，旨在传递更多信息，不代表本平台立场。版权归原作者或机构所有，如涉及侵权，请通过平台联系我们，我们将在核实后第一时间处理。本平台对转载内容的真实性、准确性不作任何保证，用户需自行判断并承担使用风险。