汽车碰撞仿真LS-DYNA控制卡片关键字.
汽车碰撞分析LS DYNA控制卡片设置共26页
谢谢你的阅读
❖ 知识就是财富 ❖ 丰富你的人生
71、既然我已经踏上这条道路，那么，任何东西都不应妨碍我沿着这条路走下去。——康德 72、家庭成为快乐的种子在外也不致成为障碍物但在旅行之际却是夜间的伴侣。——西塞罗 73、坚持意志伟大的事业需要始终不渝的精神。——伏尔泰 74、路漫漫其修道远，吾将上下而求索。——屈原 75、内外相应，言行相称。——韩非
汽车碰撞分析LS DYNA控制卡片设置
36、如果我们国家的法律中只有某种 神灵， 而不是 殚精竭 虑将神 灵揉进 宪法， 总体上 来说， 法律就 会更好 。—— 马克·吐 温 37、纲纪废弃之日，便是暴政兴起之 时。— —威·皮 物特
38、若是没有公众舆论的支持，法律 是丝毫 没有力 量的。 ——菲 力普斯 39、一个判例造出另一个判例，它们 迅速累 聚，进 而变成 法律。 ——朱 尼厄斯
LS-Dyna碰撞分析调试指南
LS-DYNA 碰撞分析调试LS-DYNA碰撞计算模型的主要检查、调试项目有：a、质量增加百分比小于5%；b、总沙漏能小于5%；c、滑移界面能；d、检查各部件之间的连接、接触关系是否定义正确，检查模型的完整性；e、检查数值输出的稳定性。
一、质量缩放Mass scale的检查：质量缩放——对于时间步长小于控制卡片中设置的最小时间步长的单元，我们通常采取增加单元材料密度的方法来增大其时间步长，以减短模型的计算时间。
关于LS-DYNA中单元时间步长的计算方法请参见附录一。
1、初步检查。
让模型在dyna中运行2个时间步，在Hyper view中调出glstat文件并检查mass scaling项（质量增加应该小于5%）；调出matsum文件并检查各部件的质量增加情况，对于质量增加过大以及有快速增长趋势的部件应检查此部件的网格质量和材料参数设置（质量增加一般是由于单元的特征长度太小或者是材料参数E、ρ设置错误，导致该单元的时间步长低于控制卡片中设置的最小时间步长，从而引起质量缩放）。
2、全过程检查。
调整模型使其符合初步检查的标准，计算模型至其正常结束。
再按[初步检查]的要求检查调试整个模型直至达到要求。
一个计算收敛的模型在其整个计算过程中，最大质量缩放应小于总质量的5% 。
二、沙漏能Hourglass energy的检查：沙漏能的出现是因为模型中采用了缩减积分引起的，我们常用的B-T单元采用的是面内单点积分，这种算法会引起沙漏效应（零能模式）。
具体介绍参见附录二。
检查：在dyna中计算模型至其正常结束。
在Hyper view中调出glstat文件并检查energy的total energy 、Hourglass energy两项，整个计算过程中沙漏能应小于总能量的5% 。
三、滑移界面能sliding interface energy的检查：滑移界面能是由摩擦和阻尼所引起的。
剧烈的滑动摩擦会引起大的正值的滑移界面能；未能检测到的穿透（undetected penetrations）常常会引起大的负值的滑移截面能。
Ls-dyna 高级教程 关键字讲解-_Lectures_1
Solver - Solution processing
DYNA3D developed at Lawrence Livermore National Laboratory by John Hallquist Low velocity impact of heavy, solid structures, military applications DYNA3D ported on Cray-1 Improved sliding interface Order of magnitude faster New material models - Explosive-structure, Soil-structure Impacts of penetration projectiles
The George Washington University
Pre-stress and Post-stress Interactive graphics Preprocessor - LS-PrePost Third party interfaces Postprocessor, LS-PrePost Other rigid body program coupling CAD d data t i interface t f
The George Washington University
Department of Civil and Environmental Engineering
Keyword Format Input File
*KEYWORD *TITLE SAMPLE INPUT FILE *CONTROL_TERMINATION 0.1000000 0 0.0000000 0 0.0000000 *DATABASE_BINARY_D3PLOT 1.00000-3 0 *DATABASE_BINARY_D3THDT 1.00000-3 *MAT_ELASTIC 1 7.89000-9 2.00000+5 0.3000000 *SECTION_SOLID 1 0 *SECTION_SHELL 1 2 1.0000000 1.0000000 1.0000000 1.0000000 0.0000000 *PART PART NAME 1 1 1 1 0 0 0 0 0 *NODE 1 0.000000000E+00 0.000000000E+00 0.000000000E+00 2 7.000000000E+00 0.000000000E+00 0.000000000E+00 3 0.000000000E+00 7.000000000E+00 0.000000000E+00 4 7.000000000E+00 7.000000000E+00 0.000000000E+00 5 0.000000000E+00 0.000000000E+00 7.000000000E+00 6 7.000000000E+00 0.000000000E+00 7.000000000E+00 7 0.000000000E+00 7.000000000E+00 7.000000000E+00 8 7.000000000E+00 7.000000000E+00 7.000000000E+00 *ELEMENT_SOLID 1 1 1 2 4 3 5 6 8 7 *PART PART NAME 2 2 2 2 0 0 0 0 0 *ELEMENT SHELL *ELEMENT_SHELL 1 2 1 2 4 3 *END
LS-Dyna碰撞分析调试指南
LS-DYNA 碰撞分析调试LS-DYNA碰撞计算模型的主要检查、调试项目有：a、质量增加百分比小于5%；b、总沙漏能小于5%；c、滑移界面能；d、检查各部件之间的连接、接触关系是否定义正确，检查模型的完整性；e、检查数值输出的稳定性。
一、质量缩放Mass scale的检查：质量缩放——对于时间步长小于控制卡片中设置的最小时间步长的单元，我们通常采取增加单元材料密度的方法来增大其时间步长，以减短模型的计算时间。
关于LS-DYNA中单元时间步长的计算方法请参见附录一。
1、初步检查。
让模型在dyna中运行2个时间步，在Hyper view中调出glstat 文件并检查mass scaling项（质量增加应该小于5%）；调出matsum文件并检查各部件的质量增加情况，对于质量增加过大以及有快速增长趋势的部件应检查此部件的网格质量和材料参数设置（质量增加一般是由于单元的特征长度太小或者是材料参数E、ρ设置错误，导致该单元的时间步长低于控制卡片中设置的最小时间步长，从而引起质量缩放）。
2、全过程检查。
调整模型使其符合初步检查的标准，计算模型至其正常结束。
再按[初步检查]的要求检查调试整个模型直至达到要求。
一个计算收敛的模型在其整个计算过程中，最大质量缩放应小于总质量的5% 。
二、沙漏能Hourglass energy的检查：沙漏能的出现是因为模型中采用了缩减积分引起的，我们常用的B-T单元采用的是面内单点积分，这种算法会引起沙漏效应（零能模式）。
具体介绍参见附录二。
检查：在dyna中计算模型至其正常结束。
在Hyper view中调出glstat文件并检查energy的total energy 、Hourglass energy两项，整个计算过程中沙漏能应小三、滑移界面能sliding interface energy的检查：滑移界面能是由摩擦和阻尼所引起的。
剧烈的滑动摩擦会引起大的正值的滑移界面能；未能检测到的穿透（undetected penetrations）常常会引起大的负值的滑移截面能。
汽车碰撞分析LS_DYNA控制卡片设置
控制卡片参数说明
*CONTROL_TIMESTEP(时间步长控制卡片) $ DTINIT TSSFAC ISDO TSLIMT DT2MS LCTM ERODE MSIST 0.0 0.9 2 0.0 -0.001 0 1 1 $ DT2MSF DT2MSLC 计算所需时间步长时，要检查所有的单元。出于稳定性原因，用0.9（缺省）来 减小时间步：Δt = 0.9 l/c ，特征长度l，和波的传播速度c，都与单元的类型有关。 DTINIT：初始时间步长，如为0.0，由DYNA自行决定初始步长； TSSFAC：时间步长缩放系数，用于确定新的时间步长。默认为0.9，当计算不稳定时，可以减小该值，但同时 增加计算时间； ISDO：计算4节点壳单元时间步长的（不同的值对应特征长度的不同算法，推荐使用2，因为此选项可以获得 最大的时间步长，但有三角形单元存在时会导致计算不稳定）； TSLIMT：壳单元最小时间步分配 ,使单元的时间步长控制在最小时间步长之上；只适用于使用 *mat_plastic_kinematic,*mat_power_law_plasticity*mat_strain_rate_dependent_plasticity,*mat_piecewise_linear_pla sticity等材料模型的壳单元，不建议使用该选项，因为使用DT2MS选项更好。 DT2MS：因质量缩放计算得到的时间步长。当设置为一个负值时，初始时间将不会小于TSSFAC*|DT2MS|。质 量只是增加到时间步小于TSSAFC*|DT2MS|的单元上。当质量缩放可接受时，推荐用这种方法。用这种方法时 质量增加是有限的，过多的增加质量会导致计算终止。当设置为正值时，初始时间步长不会小于DT2MS。单 元质量会增件或者减小以保证每一个单元的时间步都一样。这种方法尽管不会因为过多增加质量而导致计算终 止，但更难以作出合理的解释。默认为0.0，不进行质量缩放； LCTM：限制最大时间步长的Load-curve,该曲线定义最大允许时间步长和时间的关系（可选择） ； ERODE：当计算时间步长小于TSMIN（最小时间步长）时体单元和t-shell被自动删除。
LS-DYNA 关键字用户手册说明书
Copyright © 2001, Dr Brian Gladman < brg@>, Worcester, UK. All rights reserved.
The free distribution and use of this software in both source and binary form is allowed (with or without changes) provided that:
LS-DYNA®, LS-OPT® and LS-PrePost® are registered trademarks of Livermore Software Technology, an Ansys company in the United States. All other trademarks, product names and brand names belong to their respective owners.
LST reserves the right to modify the material contained within Байду номын сангаасhis manual without prior notice.
The information and examples included herein are for illustrative purposes only and are not intended to be exhaustive or all-inclusive. LSTC assumes no liability or responsibility whatsoever for any direct or indirect damages or inaccuracies of any type or nature that could be deemed to have resulted from the use of this manual.
汽车碰撞精确分析LSDYNA控制卡片设置
THKCHG：在单面接触时考虑壳厚度的改变（默认时不考虑）。
ORIEN：在初始化时可选择性的对接触面部分自动再定位。
控制卡片参数说明
ENMASS:接触单元被腐蚀的质量处理。0-节点被移除，1-体单元节点被保留，2-体单元壳单元节点被保留。 USRSTR：每个接触面分配的存储空间，针对用户提供的接触控制子程序。 USRFRC：每个接触面分配的存储空间，针对用户提供的接触摩擦子程序。 NSBCS：接触搜寻的循环数（使用三维Bucket分类搜索），推荐使用默认项。 INTERM：间歇搜寻主面和从面接触次数。 XPENE：接触面穿透检查最大乘数，默认4.0。 SSTHK：在单面接触中是否使用真实壳单元厚度，默认0,不使用真实厚度。 ECDT：时间步长内忽略腐蚀接触。
16. DATABASE_BINARY_RUNRSF 设置如下：
控制卡片参数设置
17. DATABASE_BINARY_RUNRSF 设置如下：
控制卡片参数说明
*CONTROL_TERMINATION
$ ENDTIM ENDCYC DTMIN ENDENG ENDMAS
150
0
0.0
0.0
0.0
SLSFAC：滑动接触惩罚系数 ，默认为0.1。当发现穿透量过大时，可以调整该参数；
RWPNAL： 刚体作用于固定刚性墙时，刚性墙罚函数因子系数，为0.0时，不考虑刚体与刚性墙的作用，>0时， 刚体作用于固定的刚性墙，建议选择1.0；
ISLCHK：接触面初始穿透检查，为0或1（默认）时，不检查。为2时，检查。
后面将逐一介绍碰撞分析中经常用到的控制卡片，并对每个卡片的作 用进行说明。
控制卡片使用规则
卡片相应的使用规则如下：
碰撞模拟中的控制卡片设置及意义.
1.指定输出文件【ABSTAT】——气囊统计表。
输出体积、压强、内能、气体质量流入率、气体质量流出率、质量、温度、密度。
【AVSFLT】——A VS数据【BNDOUT】——边界环境的力和能量。
输出三个方向的力。
【DEFGEO】——变形的几何体的文件【DEFORC】——离散单元。
输出三个方向的力。
【ELOUT】——单元数据。
（见DATABASE_HISTORY_OPTION）梁单元平面应力块平面应变轴向合力xx,yy,zz 应力xx,yy,zz 应力xx,yy,zz 应变S方向剪切合力xy,yz,zx 应力xy,yz,zx 应力xy,yz,zx 应变T方向剪切合力塑性应变有效应力下表面应变S方向合力矩屈服函数上表面应变T方向合力矩扭力合力【GCEOUT】——几何接触实体。
包含三个方向力和力矩。
【GLSTAT】——总体数据。
【JNTFORC】——运动副力文件【NATSUM】——材料能量。
GLSTAT JNTFORC MATSUM动能x,y,z三方向的力动能内能x,y,z三方向的力矩内能总能量沙漏能比率x,y,z三方向的动量刚性墙能量x,y,z三方向的刚体速度弹簧和阻尼能量总动能沙漏能总内能阻尼能总沙漏能滑移面能量外功x,y,z三方向速度时间步单元ID号控制的时间步【MOVIE】——【MPGS】——【NCFORC】——接触面节点力【NODFOR】——节点力组【NODOUT】——节点数据NCFORC NODOUT NODFORx方向力位移x,y,z三方向力y方向力速度z方向力加速度转动量角速度角加速度【RBDOUT】——刚体数据【RCFORC】——接触面合成力【RWFORC】——刚性墙所受的力RBDOUT RCFORC RWFORC三方向合位移三方向合力法向力三方向合速度三方向合力三方向合加速度【SBTOUT】——安全带输出文件【SECFORC】——横截面通过的力（见DA TABASE_CROSS_SECTION_OPTION）【SLEOUT】——滑移面的能量【SPCFORC】——单点约束的反作用力【SPHOUT】——SPH数据（见DATABASE_HISTORY_OPTION）【SSSTAT】——子系统数据【SWFORC】——节点约束反力（焊点和铆钉）SECFORC SLEOUT SPCFORC SWFORC x,y,z三方向力Slave能量x,y,z三方向力轴向力x,y,z三方向力矩Master能量x,y,z三方向力矩剪切力x,y,z三方向中心面积合力【TPRINT】——结构对的热量输出【TRHIST】——追踪质点时间历程信息2.时间步控制【DTNIT】——初始时间步长。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
7. CONTROL_HOURGLASS 用于沙漏的控制,通过施加附加力来阻止沙漏变形。IHQ设置为5，QH默认为0.1。
8. CONTROL_OUTPUT 用于设置输出参数。设置如下，按红色线框内设置，其余默认值。
•������ 第一列的“$”表示该行是注释行
•������ 输入的参数可以是固定格式或者用逗号分开 •������ 空格或者0 参数������ 使用该参数的默认值
控制卡片的建立
控制卡片可通过以下方式建立：
•用hypermesh在LS-DYNA模板下，选择Analysis面板点击 control cards,选择相应卡片。 •直接在key文件中输入
汽车碰撞分析LS-DYNA 控制卡片的设置
作者:张远岭
2011-4-14
控制卡片
碰撞分析控制卡片包括求解控制和结果输出控制，其中KEYWORD、 CONTROL_TERMINATION、 DATABASE_BINARY_D3PLOT是必不可少的。其他一 些控制卡片如沙漏能控制、时间步控制、接触控制等则对计算过程进行控 制，以便在发现模型中存在错误时及时的终止程序。 后面将逐一介绍碰撞分析中经常用到的控制卡片，并对每个卡片的作 用进行说明。
下面介绍在hypermesh中给出碰撞分析中经常使用的卡片的参
数设置
控制卡片参数设置
1. KEYWORD KEY文件起始关键字。该卡片可不作任何设置。
2. TITLE 输入标题名称xxx。
3. CONTROL_ACCURACY 提高计算精度的控制卡片。设置INN值为2，其余默认，不起作用。
控制卡片参数设置
4. CONTROL_BULK_VISCOSITY
体积黏性项q，人工施加压力，用于处理应力波传播突变引起的不连续。Q1,Q2设置为默认,IBQ设为-1。
5. CONTROL_CONTACT 用于接触面的计算修改。控制摩擦惩罚系数、接触初始穿透检查、接触厚度等，设置见红框，其余默认。
控制卡片参数设置
12. DATABASE_BINARY_D3DUMP 设置如下：
13. DATABASE_BINARY_D. DATABASE_BINARY_D3THDT 设置如下：
15. DATABASE_BINARY_RUNRSF 设置如下：
16. DATABASE_EXTENT_BINARY 设置如下：
控制卡片参数说明
*CONTROL_TIMESTEP(时间步长控制卡片) $ DTINIT TSSFAC ISDO TSLIMT DT2MS LCTM ERODE MSIST 0.0 0.9 2 0.0 -0.001 0 1 1 $ DT2MSF DT2MSLC 计算所需时间步长时，要检查所有的单元。出于稳定性原因，用0.9（缺省）来 减小时间步：Δt = 0.9 l/c ，特征长度l，和波的传播速度c，都与单元的类型有关。 DTINIT：初始时间步长，如为0.0，由DYNA自行决定初始步长； TSSFAC：时间步长缩放系数，用于确定新的时间步长。默认为0.9，当计算不稳定时，可以减小该值，但同时 增加计算时间； ISDO：计算4节点壳单元时间步长的（不同的值对应特征长度的不同算法，推荐使用2，因为此选项可以获得 最大的时间步长，但有三角形单元存在时会导致计算不稳定）； TSLIMT：壳单元最小时间步分配 ,使单元的时间步长控制在最小时间步长之上；只适用于使用 *mat_plastic_kinematic,*mat_power_law_plasticity*mat_strain_rate_dependent_plasticity,*mat_piecewise_linear_pla sticity等材料模型的壳单元，不建议使用该选项，因为使用DT2MS选项更好。 DT2MS：因质量缩放计算得到的时间步长。当设置为一个负值时，初始时间将不会小于TSSFAC*|DT2MS|。质 量只是增加到时间步小于TSSAFC*|DT2MS|的单元上。当质量缩放可接受时，推荐用这种方法。用这种方法时 质量增加是有限的，过多的增加质量会导致计算终止。当设置为正值时，初始时间步长不会小于DT2MS。单 元质量会增加或者减小以保证每一个单元的时间步都一样。这种方法尽管不会因为过多增加质量而导致计算终 止，但更难以作出合理的解释。默认为0.0，不进行质量缩放； LCTM：限制最大时间步长的Load-curve,该曲线定义最大允许时间步长和时间的关系（可选择） ； ERODE：当计算时间步长小于TSMIN（最小时间步长）时体单元和t-shell被自动删除。
控制卡片参数设置
17. DATABASE_BINARY_RUNRSF
设置如下：
控制卡片参数说明
*CONTROL_TERMINATION $ ENDTIM ENDCYC DTMIN 150 0 0.0
ENDENG 0.0
ENDMAS 0.0
ENDTIM：强制终止计算时间，必选，默认0.0； ENDCYC：终止循环。在计算终止时间之前，程序达到指定循环次数即终止计算。 DTMIN：确定最小时间步长TSMIN的因子。TSMIN=DTMIN * DTSTART，其中DTSTART为程序自动确定的 初始步长。当迭代步长小于TSMIN时，程序终止。 ENDENG：能量改变百分比，超过设定值则终止计算。默认0.0，不起作用； ENDMAS：质量变化百分比，超过设定值则终止计算。仅用于质量缩放DT2MS被使用时。默认0.0 ，不起作 用。
控制卡片参数设置
8. CONTROL_PARALLEL 并行计算控制，最大可以使用4个CPU。
9. CONTROL_SHELL
壳单元控制
10. CONTROL_TERMINATION 计算终止控制卡片，控制计算终止时间，设置如下：
控制卡片参数设置
11. CONTROL_TIMESTEP 计算时间步长相关参数
控制卡片使用规则
卡片相应的使用规则如下：
•������ 大部分的命令是由下划线分开的字符串, 如*control_hourglass •������ 字符可以是大写或小写 •������ 在输入文件中，命令的顺序是不重要的(除了*keyword 和*define_table) •������ 关键字命令必须左对齐，以*号开始