汽车仿真分析基础.ppt

《汽车仿真分析基础.ppt》由会员分享，可在线阅读，更多相关《汽车仿真分析基础.ppt（74页珍藏版）》请在三一办公上搜索。

1、汽车仿真分析基础臧孟炎,汽车正面碰撞安全仿真,汽车正面碰撞安全仿真,可变性障碍壁汽车偏置碰撞安全仿真,可变性障碍壁汽车偏置碰撞安全仿真,汽车侧面碰撞安全仿真,安全气囊展开与前挡玻璃作用仿真,安全气囊展开与驾驶员作用仿真,行人头部保护安全仿真,教学目的和要求,基本掌握汽车碰撞安全仿真分析软件LS-DYNA的使用方法了解隐式、特别是显式有限元方法的相关基础知识,理论教学内容,LS-DYNA简介：LS-DYNA软件的发展历程、数据结构、基本使用方法。隐式、显式有限元方法：特别是显式有限元分析的一般方法、时间步长的确定、沙漏产生原因及控制方法。材料模型：包含弹塑性材料、刚体材料在内的碰撞仿真分析常用材

2、料模型。接触方法：点面接触、面面接触的基本方法和注意事项。分析实例：介绍汽车及零部件冲击问题有限元分析实例介绍。,实验教学内容,构造一个钢球撞击一块板的有限元分析模型并进行分析和结果处理。了解基于Hypermesh的有限元网格划分、模型建立、仿真计算和结果处理的基本方法。,平衡方程：不依赖于变形状态的线性方程；几何方程：应变和位移的关系呈线性；物性方程：应力和应变的关系呈线性；边界条件：力边界上的外力和位移边界上的位移是独立的，或线性依赖于变形状态。上述方程或边界条件中的任何一个不符合所述特点，就成为非线性问题。,非线性有限元,非线性有限元,非线性问题：材料非线性 塑性域的发生，高温条件下结构

3、的蠕变几何非线性 大变形（板壳结构的大挠度，屈曲）边界非线性 接触和碰撞,非线性问题分析步骤：,1.明确分析目的，决定模型化方针明确分析目的：确认“什么”评价标准（设计规范、安全标准）、评价部位（冲击部位、焊接部位）选择求解器：关键是静态分析还是动态分析？决定模型形状：形状省略、对称性利用选择单元类型材料物性选择确定使用条件选择单位,非线性问题分析步骤：,2.构造分析数据（前处理器）网格划分分析条件设定（载荷、边界、约束和接触等）材料定义（物性和单元断面定义等）分析控制参数确定（分析结束时间等）输出选择设定,非线性问题分析步骤：,3.执行计算（求解器）分析计算开始计算数据文件确认分析计算过程确

4、认分析计算结束,非线性问题分析步骤：,4.分析计算结果研究（后处理器）,显式有限元分析软件-LS-DYNA,John O.HallquistPresident and FounderLawrence Software Technology Corporation,John Hallquist earned his B.S.in Industrial Engineering,from Western Michigan University in 1970.He received an M.S.degree in Engineering Mechanics from MTU in 1972，and

5、 earned a Ph.D.degree in Mechanical Engineering and Engineering Mechanics in 1974.John joined the weapons laboratory at Lawrence Livermore National Laboratory(LLNL).,显式有限元分析软件-LSDYNA,LS-DYNA的输入数据结构,KEYWORD,数值卡,节点定义,立体单元定义,初速度定义,面的定义,定义顺序可任意变动,*KEYWORD,*END,LS-DYNA的输入数据结构,材料定义,单元类型-体单元,定义体单元关键字需要的参数：

6、*ELEMENT_SOLID体单元的ID号。PART的ID号8个节点的ID号单元公式的指定，通过*SECTION_SOLIDSECTION的ID号体单元计算公式。0专门用于蜂窝材料*MAT_MODIFIED_HONEYCOMB的体单元1-缺省的常应力单元2-全积分S/R体单元,*ELEMENT_SOLID$#eid pid n1 n2 n3 n4 n5 n6 n7 n8 2843 7 285 284 289 287 6183 6186 6185 6184 2844 7 284 268 267 289 6186 6188 6187 6185 2845 7 287 289 290 288 6184

7、 6185 6190 6189,单元类型-体单元,1-缺省的常应力单元公式，也称为1点积分单元积分点在单元的中心，需要控制沙漏。通常用于8节点6面体单元，但也可用于退化后的4节点、6节点单元，但精度有所下降。对于弯曲板件，厚度方向至少划分两个单元（因为该单元为常应力单元，表示拉压应力状态至少需要两个单元2-全积分S/R体单元单元内8点高斯积分，没有沙漏发生；但需要1点积分单元4倍的计算时间。没有体积锁死（由于采用了选择性缩减积分，且各积分点压力成分相同）大变形问题时不是很稳定,单元类型-壳单元,定义壳单元关键字所需要的参数：*ELEMENT_SHELL壳单元的ID号。PART的ID号4个节点的

8、ID号单元公式的指定，通过*SECTION_SHELLSECTION的ID号体单元计算公式。1-Hughe-Liu壳单元2-Belyschko-Tsay单元16-全积分壳单元（非常快）,*ELEMENT_SHELL$#eid pid n1 n2 n3 n4 3639 8 7781 7783 7785 7780 3640 8 7780 7785 7776 7777 3641 8 7769 7786 7784 7768,*SECTION_SHELL$HMNAME PROPS 1glass$#secid elform shrf nip propt qr/irid icomp setyp 1 2$#t

9、1 t2 t3 t4 nloc marea 2.000000 2.000000 2.000000 2.000000,单元类型-壳单元,节点的逆时针方向确定壳单元正面。板厚方向1-5点积分指定可能。单元中使用平面应力单元计算公式。应力应变的输出为积分点位置的值。,单元类型-壳单元,1-Hughes-Liu壳单元LS-DYNA中最早的壳单元，定义来自体单元的退缩形。面内一点积分，采用Jaumann应力更新，对翘曲的几何体有效。能对壳单元的参考面进行偏置（中、内、外）计算时间是BT壳单元的两倍,单元类型-壳单元,2-Belytschko-Tsay壳单元缺省的壳单元公式，面内单点积分，计算速度很快，对

10、大变形问题是最稳定有效的方法。采用Co-rotational应力更新、单元坐标系置于单元中心、基于平面单元假定，所以对于翘曲的几何体不适用。,单元类型-壳单元,16-全积分壳单元以BT壳单元为基础，面内4点积分。只比BT壳多2-3倍计算时间。可以处理翘曲的几何体问题,此时需要激活第种沙漏控制公式。,单元类型-梁单元,定义梁单元关键字既需要的参数：*ELEMENT_BEAM梁单元的ID号。PART的ID号3个节点的ID号单元公式的指定，通过*SECTION_BEAMSECTION的ID号梁单元计算公式。1-Hughe-Liu断面积分梁2-Belytschko-Schwer合力梁3-杆单元9-可变

11、形焊点梁单元,*ELEMENT_BEAM$#eid pid n1 n2 n3 3639 8 7781 7783 7785 3640 8 7780 7785 7776 3641 8 7769 7786 7784,*SECTION_BEAM(Material)$secid elform shrf qr/irid cst scoor 3 9 0.0000000 0.0000000 1.0000000 0.0000000$ts1 ts2 tt1 tt2 nsloc ntloc 6.5000000 6.5000000 0.0000000 0.0000000 0.0000000 0.0000000,单元类

12、型-梁单元,1-Hughe-Liu断面积分梁,N1,N2两节点确定r轴方向，第三个节点确定断面内惯性主轴。节点立体单元节点退缩为节点,单元类型-梁单元,1-Hughe-Liu断面积分梁,断面内积分点数和断面形状需要定义断面形状为矩形和圆管可在SECTION_SHELL下定义，工字梁等的其他断面，使用*INTEGRATION定义与轴垂直方向上的应力假定为平面应力状态轴向积分点为单元中心一点，并在此截面计算应力单元轴向弯矩一定应力输出为局部坐标系。,2-Belytschko-Schwer合力梁,不计算梁单元应力，仅计算节点上的力和力矩,3-杆单元,不需要第三个节点；断面内一点积分,单元类型-梁单元

13、,单元类型-梁单元,9-可变形点焊梁单元,需要*MAT_SPOTWELD材料定义，包括塑性和失效功能,单元类型-梁单元,9-可变形点焊梁单元,单元类型-梁单元,9-可变形点焊梁单元,单元类型-梁单元,9-可变形点焊梁单元,单元类型-弹簧阻尼单元,使用两个节点和离散的材料来定义无质量的弹簧阻尼单元能与其他所有显式单元连接具有平动和转动自由度能定义复杂的力-位移关系由于只能同时定义一个弹簧或阻尼，所以定义弹簧-阻尼集合时需要重叠定义两个单元,单元类型-弹簧阻尼单元,求解方法,静态线性问题特征值问题瞬态分析方法：隐式和显式,静态线性问题,Solve KD=F for DGauss eliminati

14、onLU decompositionEtc.,自由振动分析（特征值分析）,(Homogeneous equation,F=0),Assume,Let,K-li M fi=0,(Eigenvector：特征向量),(Roots of equation are the eigenvalues:特征值),特征值方程,自由振动分析（特征值分析）,Methods of solving eigenvalue equationJacobis methodGivens method and Householders methodThe bisection method(Sturm sequences)Inve

15、rse iterationQR methodSubspace iterationLanczos method,瞬态分析方法,Structure systems are very often subjected to transient excitation.A transient excitation is a highly dynamic time dependent force exerted on the structure,such as earthquake,impact,and shocks.The discrete governing equation system usuall

16、y requires a different solver from that of eigenvalue analysis.The widely used method is the so-called direct integration method.,瞬态分析方法,The direct integration method is basically using the finite difference method for time stepping.There are mainly two types of direct integration method;one is impl

17、icit and the other is explicit.Implicit method(e.g.Newmarks method)is more efficient for relatively slow phenomenaExplicit method(e.g.central differencing method)is more efficient for very fast phenomena,such as impact and explosion.,Assume that,Substitute into,Newmarks method(Implicit),常数,Newmarks

18、method(Implicit),where,Therefore,初始条件：,显式时间积分,考虑简单的单自由度线形弹簧阻尼系统，根据迖朗贝尔动力学原理：,可以用解析方法求解该常微分方程。若为非线性问题，如k为u的函数,此时用解析方法一般很难求解，只好应用数值方法来求解，常用的为有限差分法和有限元法。上面的公式具有普遍意义，对于有限元法而言，上述方程的矩阵形式为：,求解该方程使用两种方法：振型叠加法；逐步积分法。逐步积分法又可分为增量法、迭代法和混合法。隐式求解方法一般采用增量迭代法，需要转置换刚度矩阵，通过一系列线性逼近求解，对于存在内部接触这样的高度非线性动力学问题往往无法保证收敛。,有限中

19、心差分方程,显式时间积分,考虑简单的单自由度线形弹簧阻尼系统，根据迖朗贝尔动力学原理：,可以用解析方法求解该常微分方程。若为非线性问题，如k为u的函数,此时用解析方法一般很难求解，只好应用数值方法来求解，常用的为有限差分法和有限元法。上面的公式具有普遍意义，对于有限元法而言，上述方程的矩阵形式为：,求解该方程使用两种方法：振型叠加法；逐步积分法。逐步积分法又可分为增量法、迭代法和混合法。隐式求解方法一般采用增量迭代法，需要转置换刚度矩阵，通过一系列线性逼近求解，对于存在内部接触这样的高度非线性动力学问题往往无法保证收敛。,位移速度,时间,1,2,显式时间积分,LS-DYNA采用显式中心差分法来

20、进行时间积分。在已知0,tn时间步阶的情况下，求解tn+1时间步的解，运动方程为：,：外力矢量,：内力矢量,：HOURGLASS（沙漏）阻力矢量,：粘性阻力矢量,位移速度,时间,1,2,3,4,5,6,7,显式时间积分,显式中心差分法的稳定是有条件的（库仑条件）：使用中心差分法进行时间积分时的时间增量，必须比在介质中传播的应力波穿过任何单元所需要的时间短一些。,不满足库仑条件，就不能保证解的数值稳定性。设构成模型的各单元应力波通过时间的最小值为,应力波的传播速度计算公式依单元形式而异：,显式时间积分,在运动方程中有一项沙漏阻力：该项阻力是人为加上的力，为了防止沙漏变形。沙漏变形是由于LS-DY

21、NA采用单点高斯积分进行时间积分而导致的零能变形模式。,如果使用一点积分，在计算应变成份时即使节点发生位移应变值为0，因而不产生应力的现象。这种现象如果任其泛滥，除变形失真外，单元代表长度变短，时间步长变小；单元形状发生畸变，计算精度降低。,显式时间积分,以立体单元为研究对象：,：节点a在i方向上的位移,显式时间积分,用中心一点积分进行应变评价时，下式成立：,由于对角线上的节点位移相等，应变为零,显式时间积分,显式时间积分,显式时间积分,接触问题,分析对象的一部分与另一部分接触而发生相互作用时，需要在解除部分定义接触面，不然的话，原本相接触的各物体就会毫无阻力地贯穿，而不发生相互作用。接触面的

22、定义：一面为主面（Master），另一面为从面（Slave）。每一接触面分别由三节点或四节点构成的segment组成，主面称为Master segment，从面称为slave segment。构成segment的节点分别称为slave node和master node。接触判断是以节点为单位进行。同一物体的不同部分相接触的接触面，成为单面接触，只要定义从面。在被定义的接触面上，在被判断为接触后的节点和segment上，将作用由单元刚度计算的接触反力。同时还可以考虑摩擦力。,接触问题,在LS-DYNA中，有四十多种接触处理方法。根据分析模型的不同，有多种处理主从面关系的方法。但归纳起来，有以下4

23、种。滑动接触（sliding only）：接触面不分离，无摩擦滑动。固结接触（tied interface）：接触面被刚接，没有相对位置变化。一般接触（slide&void）：接触面时而接触，时而分离。单面接触（single surface）：一个面不同部分的相互接触。还有将接触面上发生严重畸变的单元自动剔除的接触类型。,接触问题,接触面的处理方法，大致可分为如下三类。Kinematic Constraint Method：节点约束法，用于tied interface。Penalty Method：惩罚法。Distributed Parameter Method：分布参数法，用于滑动接触类型。

24、,接触问题,接触计算过程中花费时间最多的是接触点的搜索（slave search）。从输入的大量segment和节点中，每一循环都必须调查相接触的节点和segment。所以计算效率的提高非常重要。Ls-dyna中，根据接触类型的不同，有如下两种节点搜索方法。利用segment形状连续性的方法（connectivity）。利用调查空间节点位置的bucket sort方法。利用bucket sort的方法，segment的方向性不需考虑时，计算效率高。,接触类型分类,接触计算过程包含以下三步：接触计算用segment数据定义接触面搜索被定义构成接触面的节点，或segment各节点与哪个segmen

25、t接触？接触反力计算，消除接触面渗透,接触类型分类,接触类型可以按下述方法分类按接触点搜索方法分类：connectivity,bucket sort接触反力计算方法：penalty,constraint接触面的定义方法：面-面接触(symmetry)，点-面接触(non-symmetry)将接触点搜索方法，反力计算和接触面定义组合起来，形成接触类型：SURFACE_TO_SURFACE（C,P,S）NODE_TO_SURFACE（C,P,N）AUTOMATIC_SURFACE_TO_SURFACE（B,P,S）CONSTRAINT_SURFACE_TO_SURFACE（B,C,S）AUTOMA

26、TIC_SINGLE_SURFACE（B,P,N）使用Soft constraint，penalty接触类型使用constraint方法计算反力。,1.接触面搜索 接触处理是计算中最耗时的工作，而其中最主要的部分是接触面搜索。有接触可能的零件，segment在输入文件中已经定义，接触面搜索的工作在于确定哪个部分与哪个部分相接触。在程序内部，调查某一特定的slave节点与哪个master segment接触。为了选出相接触的segment，需要找到与slave节点最近的master节点（构成master segment的节点）。这个工作是单纯地计算slave节点和各master节点间的距离，选出

27、最近的节点，但这需要花费许多时间，因此，需设法节约计算时间。,接触面搜索方法特征,1.接触面搜索 方法之一是利用单元形状的连续性。将上一个时间步接触的master segment事先记忆下来，首先确认该segment是否接触。如果没有可能接触，就从与slave节点最靠近的、构成master segment的节点开始顺序搜索。方法之二是利用bucket sort方法。这个方法是将节点分成组，在每个组内计算节点距离，达到减少计算工作量的目的。无论使用哪种方法，第一步的工作都是要找到最近的master节点。,接触面搜索方法特征,2.利用单元形状的连续性 利用master segment形状的连续性计

28、算接触点。在接触范围限定的情况下，效率高，但当网格不连续时，可能带来接触不安定。,接触面搜索方法特征,3.利用三维节点的空间布置的方法 特征：接触范围广时，搜索高效，不需考虑接触面的方向。说明：将包含网格的空间分为许多立方体（MAX5000）；搜索包含master，slave的立方体。分析错误：立方体最大为5000，立方体大小依存于网格。如果有单元发散而产生大变形时，一个立方体内包含segment将超过2000，程序会给出错误信息（more than 2000 nodes in bucket)。Bucket sort并非每个时间步都做，single_surface（25），automatic_

29、single_surface(100)。这种搜索方法与segment指向无关，壳单元必须考虑板厚。,接触面搜索方法特征,3.接触面搜索方法不同的影响,接触面搜索方法特征,3.接触面搜索方法不同的影响,接触面搜索方法特征,由接触面搜索，确定与slave节点接触的segment后，下面计算slave节点与该segment的接触位置。,接触反力计算方法特征,确定接触位置后，给slave节点和master节点加上接触反力。接触反力的计算方法有两种。,接触面反力计算方法特征,penalty方法(2)Constraint 方法,k:接触弹簧刚度，l:接触面穿透深度，n:接触面法向矢量，fsi:penalt

30、y系数，K:体积弹性模量，A:接触面积，V:接触单元体积，dt:时间增量，ms:slave节点集中质量,(1)Penalty方法的特征：该方法的关键在于在接触节点和接触面上设定有垂直弹簧。弹簧刚度以接触面上的单元为单位计算，而且用户可以进行刚性参数的调节。最初，确认slave节点是否穿透主面，如果没穿透，当然不进行接触反力计算。如有穿透，接触反力被加载在从节点和主面节点之间，接触反力的大小与穿透量成反比。,接触反力计算方法特征,(1)Penalty方法的特征：,接触反力计算方法特征,If l0,If l0,立体单元,壳单元,li:壳单元对角线长，fsi:penalty系数（通常取0.1），Ki:体积弹性模量，Ai:接触segment面积，Vi:接触单元体积,接触面定义,面-面接触(symmetry)，点-面接触(non-symmetry),SURFACE_TO_SURFACENODE_TO_SURFACEAUTOMATIC_SURFACE_TO_SURFACECONSTRAINT_SURFACE_TO_SURFACEAUTOMATIC_SINGLE_SURFACE,