首页/文章/ 详情

LS-DYNA 分析报告

5月前浏览255


目 录

(一) LS-DYNA基本介绍

(二) 某电动玩具的跌落分析

(三) 汽车保险杠的碰撞分析

(四) 总结


(一)LS-DYNA 的基本介绍

LS-DYNA 是一个以显式求解为主、兼有隐式求解功能,以 Lagrange 算法为主、兼有ALE和Euler 算法,以结构分析为主、兼有热分析和流体-结构耦合功能,以非线性动力分析为主、兼有静力分析功能, 军用和民用相结合的通用非线性结构动力分析有限元程序,主要用于求解各种非线性结构的高速碰撞、爆炸和金属成型等结构非线性问题。

DYNA 程序系列最初是 1976年在美国Lawrence Livermore National Lab.由J.O.Hallquist 主持开发完成的,主要目的是为武器设计提供分析工具,后经1979、1981、1982、1986、1987、1988 年版的功能扩充和改进,成为国际著名的非线性动力分析软件,在武器结构设计、内弹道和终点弹道、军用材料研制等方面得到了广泛的应用。

1988 年 J.O.Hallquist 创建 LSTC 公司,自此开始DYNA程序的商业化开发,LSTC 陆续将DYNA 的显式、隐式、热分析等系列程序组合在一起,形成一个整体的LS-DYNA 软件包,并逐步增加汽车安全性分析、薄板冲压成型过程模拟、流体与固体耦合(ALE和欧拉算法)等功能,使LS-DYNA程序系统在国防和民用领域的应用范围不断扩大,并建立了完备的质量保证体系。因此LS-DYNA 一经推出,即在显式有限元分析领域引起轰动效应,大大拓展了LS-DYNA 的用户领域,在中国地区,LS-DYNA 的用户数在短短的几年时间内即超过了200家,远远领先于其它显式分析程序。

LS-DYNA 程序的主要强项在于:

  • 历史悠久、应用广泛。该软件是全世界范围内最知名的有限元显式求解程序。LS-DYNA 程序开发的最初目的是为北约组织的武器结构设计、防护结构设计服务,是该组织的PublicDomain程序,后来商业化后广泛传播到世界各地的研究机构。从理论和算法而言, LS-DYNA是目前所有的显式求解程序的鼻祖和理论基础,经过多年的发展,它已经成为功能最丰富、应用领域最广、全球用户最多的有限元显式求解程序,在高度非线性结构动力学分析领域具有无可争议的领导地位。

  • 功能丰富:LS-DYNA融合了显式和隐式两种求解技术,不但可以进行动力分析,还能进行一定程度的静力分析。在显式求解技术方面, 不但涵盖了传统的拉格朗日(Lagrange)、欧拉(Euler)和任意拉格朗日-欧拉耦合(ALE 流固耦合)算法,而且提供了目前数值技术的先进技术-无网格算法(Meshfree),典型的有 SPH和EFG算法。

  • 强大的结构分析功能:相比其他显式程序,LS-DYNA在结构动力分析的功能上可以说无出其右者。主要体现在两点:(1)最丰富的接触碰撞算法;LS-DYNA 提供了几乎所有结构相互作用所需的接触算法。从接触算法原理上,有罚函数法、约束法和分配参数法; 从计算空间上,有1D、2D和3D接触;从材料特性来看,有变形体间接触、变形体和刚体接触以及刚体间接触;从行业来看,有制造业行业的成形接触和拉延筋接触,汽车行业专用的气囊单面接触和焊点接触,军工行业常用的侵蚀接触和滑移接触(模拟炸药和结构间接触);从接触状态来看,有固连、断开接触、约束接触、过盈接触、边接触和单面接触。还可以考虑接触界面的传热并提供初始穿透及接触力的处理。(2)最全面的单元技术;LS-DYNA提供了所有显式程序中最全面的单元算法:Lagrange、Euler、ALE、SPH、EFG,所有这些单元算法,都可以用于结构响应分析,其流固耦合功能也是侧重研究结构所受的冲击响应。

  • 强大的并行计算和二次开发功能:LS-DYNA很早就开始高级并行算法的研究和应用,其SMP和MPP并行计算能力在各种平台上均具有极好的并行特性,且支持其所有计算功能。LS-DYNA 具备完善的二次开发功能,不仅提供了材料和摩擦类型的二次开发功能,还提供了用户定义的程序控制、接口控制等多种接口功能,使得我们可以灵活的用户化的发挥LS-DYNA的强大功能。

单元类型

LS-DYNA程序的单元类型众多,有二维、三维单元,详细可分为薄壳、厚壳、体、梁单元,ALE、Euler、Lagrange单元等。单元积分采用沙漏粘性阻尼以克服零能模式,单元计算速度快、存储空间少,可以满足各种实体结构、薄壁结构和流固耦合结构的有限元网格划分的需要。其中960版本提供SPH(SmoothedParticleHydrodynamics)算法,该算法是一种无网格Lagrange算法,可用于分析连续体结构的解体、破裂、固体的层裂、脆性断裂等。各种

单元都有多种算法可供用户选择,如实体单元单点积分,薄壳单元面内单点积分、沿壳厚多点积分。线性位移函数和单点积分的显式单元能够很好的用于大变形和材料失效等非线性问题。

材料类型

LS—DⅥqA970版本有100多种金属和非金属材料模型可供选择,如弹性、弹塑性、超弹性、泡沫、玻璃土壤、 混凝土、流体、复合材料、炸药及起爆燃烧、刚性及用户自定义材料,并可考虑材料的失效、损伤、粘性、蠕变、与温度或应变率相关等特性。

接触分析功能

LS.DYNA程序的全自动接触分析易于使用,功能强大。由40多种选择可以求解不同的接触问题,如变形体与变形体的接触、变形体与刚体的接触、板壳结构的单面接触(屈曲分析)、与刚性墙接触、表面与表面的固连、节点与表谣的圆连、壳边与壳面的固连、流体与固体的接触等, 并可考虑接触表面的静、动摩擦力(库伦摩擦、粘性摩擦和用户自定义摩擦)、热传导和固连失效等。这些技术成功地运用于整车碰撞模拟研究、乘员与柔性气囊或安全带接触的安全性分析、金属成型中的薄板与冲头和模具接触,高速弹丸对靶板的穿甲模拟计算等。


(二)某电动玩具的跌落分析

下 图 几 何 模 型 由 某 公 司 提 供 , 几 何 模 型 包 括 手 枪 钻 各 部 件  :

 image.png

大部分部件采用壳单元离散,电池、马达,轴等划分为实体单元。

 5f35a46f-83f2-400e-a066-d5da8cf43ec3.png

 有限元网格图形,节点总数为 86126,实体单元为 18229(含地面模型),壳单元为 58779。

 image.png

•  载荷、初始条件

手 枪钻跌落过程为自由落体,高度 1.829 米,跌落方向为模型中坐标系 Y 正向或负向。为减少计算时间,计算从手 枪钻与地面即将接触开始,即初始条件中整个手 枪钻具有沿+Z 方向或-Z 方向初速度 5.98 米/秒。

•  定义联接和接触

手 枪钻中存在的联接主要为螺钉联接和卡槽联接,采用 LS-DYNA 的约束方程算法定义,在此模型中,网格划分过程中最大程度遵照实际结构。模型中任何正在接触及受力变形后可能发生接触的部位都定义了接触界面。此模型中共定义了 20 个接触界面。

材料模式和参数

主要使用了 2 种材料模式:弹塑性材料、刚性材料。手 枪钻的塑料部件采用弹塑性材料马达、地面等采用刚性材料。

整机实际质量为 2.7Kg,有限元模型的总质量可由 LS-DYNA 输出列表,为2.705Kg,质量误差为 0.2%。

根据实验条件,选取以下两种跌落作为计算方案。

方案一(如下图)

image.png 

方案二(如下图)

image.png36f99368-ce62-4253-9336-d5e145649a90.png

image.png

8c6d5705-f129-4c7f-8337-6d5002158e15.png

上图说明在冲击发生后,电池包上盖根部处由于受到瞬间的剪切作用,发生较大的塑性变形。

5f91c264-069d-488a-a315-666b72ef10f9.png

上图说明在冲击发生后,冲击力传递到变速箱壳,在螺钉连接处发生较大的塑性变形,可能发生破坏。

对于有限元仿真计算,系统的能量变化非常重要,因为对于一个封闭系统而言,总的能量必须保持平衡 。

7dc29ac1-1024-44df-b6be-6ec17d3851f0.png
4968b9d7-c2eb-48c9-b5e2-c0ec8f6096ed.png

上图说明手 枪钻对地面的冲击力峰值约为 25KN。

30b20855-5ed9-478d-9367-2764c7c883df.png

98826438-95de-43db-a31d-b6d60edc1b00.png

9b3bfba6-d794-43a2-9faa-ab37101d7798.png




a178feda-edce-45db-9cf9-bf3634860ad3.png

通过与方案一比较可知,在该方案下电池包上盖承受更大的冲击能量,,屈服更大,很有可能发生破坏。


image.png

更多内容见附件


免责声明:

本页面/内容部分素材来源于互联网公 开 信 息,旨在传递更多信息,不代表本平台立场。

版权归原作者或机构所有,如涉及侵权,请通过平台联系我们,我们将在核实后第一时间处理。

本平台对转载内容的真实性、准确性不作任何保证,用户需自行判断并承担使用风险。

附件

免费C4-LS-DYNA分析报告.pdf
Fluent
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-03-27
最近编辑:5月前
仿真支持爱好者
在仿真的路上越走越远
获赞 333粉丝 21文章 318课程 0
点赞
收藏
作者推荐

如何写input文件

一、输入文件的组成和结构:1.一个输入文件由模型数据和历史数据两部分组成.模型数据的作用:定义一个有限元模型.包括单元,节点,单元性质,定义材料等等有关说明模型自身的数据.模型数据可被组织到零件中(零件可以被组装成一个有意义的模型).历史数据的定义是模型发生了什么----事情的进展,模型响应的荷载,历史被分成一系列的时步层序.每一步就是一个响应(静态加载,动态响应等),时步的定义包括过程类型(比如静态应力分析,瞬时传热分析等)对于时间积分的控制参数或者非线性解过程,加载和输出要求.Ataminimumthemodelconsistsofthefollowinginformation:geometry,elementsectionproperties,materialdata,loadsandboundaryconditions,analysistype,andoutputrequests.2.ABAQUS输入文件的结构形式。1)必须有一个*HEADING开头。2)接下来就是模型数据部分,定义节点,单元,材料,初始条件等。模型数据的层次为:部件,组装,模型。必须的模型数据:(1)几何数据:模型的几何形状是用单元和节点来定义的,结构性单元的截面是必须定义的。比如梁单元。特殊的特征也可以用特殊的单元来定义,比如弹簧单元,阻尼器,点式群体等。(2)材料的定义:材料必须定义比如使用的是钢啊,岩石,土啊等材料。可选的模型数据:(1)零件和组合:一个模型可以用几个零件来定义有可以把几个零件组合成一个集来定义。(2)初始条件:比如初始应力,温度,或者速度等(3)边界条件:(4)运动约束(5)相互作用(6)振幅定义(7)输出控制(8)环境特性(9)用户子程序(10)分析附属部分3)接下来就是历史数据:定义分析的类型,荷载,输出要求等。分析的目的就是预测模型对某些外部荷载或者某些初始条件的反映。一个ABAQUS分析是建立在STEP的概念上的,(在历史数据中描述)在分析中可以定义多个STEPS。每个STEP用*STEP开始,用*ENDSTEP结束。*STEP是历史数据和模型数据的分界点,第一次出现*STEP的前面的是模型数据,后面的就是历史数据。必须的历史数据:响应类型:必须立刻出现在*STEP选项后面。ABAQUS中有两种响应步,一种是总体分析响应步,可以是线形和非线形的,另一种是线形扰动步。可选历史数据:(1)荷载:通常定义某种加的荷载类型和大小。荷载可以被描述成时间的函数。(2)边界条件输出控制(3)辅助控制(4)再生单元和曲面二、书写input文件的语法和规则:1.关键词行:1)必须以*开始,后面接的是选项的名字,然后随之定义选项的内容.如:*MATERIALNAME=STEEL注释行是以**开始的.2)如果有参数,则参数和关键词之间必须用“,”格开。3)在参数之间必须用“,”格开。4)关键词行中的空格可以忽略。5)每行的长度不能超过256个字符6)关键词和参数对大小写是不区分的。7)参数值通常对大小写也是不区分的,但是唯一的例外是文件名区分大小写。8)关键词和参数,以及大多数情况下的参数值是不需要全拼写出来的,只要他们之间可以相互区分就可以了。9)假如参数有响应的值,则付值号是“=”。10)关键词行可以延续,比如参数的名字很长,要在下一行继续这个关键词行的话就可以这样做,用“,”来连接。比如:*ELASTIC,TYPE=ISOTROPIC,DEPENDENCIES=111)有些选项允许INPUT和FILE的参数作为一个输入文件名,这样的文件名必须包括一个完整的路径名或者是一个相对路径名。2.数据行:(数据行如果和关键词相联系必须紧跟关键词行)1)一个数据行包括空格在内不能超过256个字符。2)所有的数据条目之间必须用“,”格开。3)一行中必须包括指定说明的数据条目的数字。4)每行结尾的空数据域可以省略。5)浮点数最多可以占用20个字符。6)整数可以是10个7)字符串可以是80个8)延续行可以被用到特定的情况。3.标签:什么是标签呢,比如集名,曲面名,钢筋名,他们是区分大小写的,长度可以有80个字符长。标签中的空格是可以省略的,除非用“”来标示。那就不能省略了。没有用“”来标示的标签必须用字母来开头。如果一个标签用“”来定义那么“”也是标签的一部分。标签的开始和结束不能用双重“_”下面是一个使用了引号和没有使用引号的例子:*ELEMENT,TYPE=SPRINGA,ELSET="Oneelement"1,1,2*SPRING,ELSET="Oneelement"1.0E-5,*NSET,ELSET="Oneelement",NSET=NODESET*BOUNDARYnodeset,1,24.数据行重复:数据行可以重复,也就是说每行数据可以有一行响应的变量也可以有几行。同样也可以有多行数据行,对应各自的变量行.如:*ELASTIC,TYPE=ISOTROPIC200.E3,0.3,20.150.E3,0.35,400.80.E3,0.42,700.定义了一个零件的材料性质,均质,线弹性,在不同应力下的扬式模量和泊松比.三、从外存储器中引入模型或者历史数据The*INCLUDEoptioncanbeusedtospecifyanexternalfilethatcontainsaportionoftheABAQUSinputfile.Thisfilecanincludemodelandhistorydefinitiondata,commentlines,andother*INCLUDEreferences.Whentheoptionisencountered,ABAQUSwillimmediatelyprocesstheinputdatawithinthefilespecifiedbytheINPUTparameter.Whentheend-of-fileisreached,ABAQUSwillreturntoprocessingtheoriginalfile.Amaximumoffivelevelsofnested*INCLUDEreferencescanbeused.TheoptionsupportsUNIXenvironmentvariables.Usage:*INCLUDE,INPUT=file_name四、举例:下面介绍两个简单的例子:例一:悬臂梁.输入文件的开始就是文件头,以HEADING开始,随后是模型的名字,如下所示:*HEADINGCANTILEVERBEAM然后是网格定义:现在就是模型数据的开始了.一般选择从网格的定义开始,网格包括(单元和节点)假如我们的悬臂梁有五个单元,六个节点,下面我们首先详细说明节点:*NODE,NSET=ENDS1,0.6,100.*NGEN1,6节点组集,NSET其值(名字)为ENDS.下面的就是这样理解的,第一个节点是从0开始的,第六个节点是在100结束的.同样我们来定义单元:*ELEMENT,TYPE=B21(单元类型)1,1,2(单元类型的参数)*ELGEN,ELSET=BEAM(产生单元集,及其名称)1,5(一个单元集,包括5个单元)现在定义单元的性质:*BEAMSECTION,SECTION=RECTANGULAR,ELSET=BEAM,MATERIAL=STEEL1.,2.梁截面,截面的形状是矩形,单元集的名称是梁单元,材料是钢。截面的尺寸是1*2。下面定义材料的性质:*MATERIAL,NAME=STEEL*ELASTIC30.E6,材料是钢,弹性,弹性模量是30E6。下面定义边界:*BOUNDARY6,ENCASTRE边界是在6节点,通过ENCASTRE来描述。边界也可以用下面的形式来定义:*BOUNDARY6,1,6ABAQUS对结构单元的中的节点的自由度使用常规的编号方式。1,2,3代表的是位移分量;4,5,6代表的是旋转分量。以上是模型数据的定义,下面开始历史数据的定义:(加载的次序,事件的发生,还有我们想看到的变量的响应)时步(步骤)的定义。*STEP,PERTURBATION(步骤的开始,扰动是其名称)*STATIC(静态分析)*CLOAD(集中荷载)1,2,-20000.(在节点1,y[2代表的是y方向,假如是1则代表的是x方向,3代表的是z方向]方向施加荷载,荷载的大小是-20000。也就是向下施加荷载。)*ENDSTEP(步骤的结束)下面来解释输出要求:*ELPRINT,POSITION=AVERAGEDATNODES,(节点的平均值)SUMMARY=YES在表的下部求和)S11,(积分点的应力分量在X方向)E11(在积分点的应变分量在X方向)SF,(在积分点的截面力)*NODEFILE,NSET=ENDSU,(节点的空间位移)CF,(节点的点荷载)RF节点的反作用力)*OUTPUT,FIELD,VARIABLE=PRESELECT*ELEMENTOUTPUTSF,*OUTPUT,HISTORY*NODEOUTPUT,NSET=ENDSU,CF,RF更多内容见附件免责声明:本页面/内容部分素材来源于互联网公开信息,旨在传递更多信息,不代表本平台立场。版权归原作者或机构所有,如涉及侵权,请通过平台联系我们,我们将在核实后第一时间处理。本平台对转载内容的真实性、准确性不作任何保证,用户需自行判断并承担使用风险。

有附件
未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈