汽车碰撞仿真ls-dyna关键字
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_关键字文件
$
AOPT
2.0
$
XP
YP
ZP
A1
A2
A3
1.000E+00 0.000E+00 0.000E+00
$
V1
V2
V3
D1
D2
D3
0.000E+00 1.000E+00 0.000E+00
*SECTION_SHELL
$PROPERTY NAME:blankpro
$
SID ELFORM
SHRF
NIP
PROPT QR/IRID
*KEYWORD
该分析项目的标题
*TITLE Numisheet 96: S-rail Stamping
关键字文件
关键字文件基本结构
板料 (Blank)
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
SOFT SOFSCL LCIDAB MAXPAR PENTOL
0
$ PENMAX THKOPT SHLTHK
SNLOG
1
MBOXID
SPR
MPR
1
1
PENCHK
BT
DT
00.000E+0001.00000E20
SFMT
FSF
VSF
DEPTH
BSORT FRCFRQ
凹模运动控制
*BOUNDARY_PRESCRIBED_MOTION_RIGID
*SET_PART_LIST
10-ls-dyna关键字命令介绍介绍
LIVERMORE SOFTWARE TECHNOLOGY CORPORATION
*NODE
Nodal ID x, y, z 方向在全局坐标系中的坐标值 平动和转动约束 • 可使用 *boundary_spc添加约束
LIVERMORE SOFTWARE TECHNOLOGY CORPORATION
出也可以在材料坐标系中输出 ( *mat_orthotropic_elastic 除外)
LIVERMORE SOFTWARE TECHNOLOGY CORPORATION
多的存储结果) • D3DRLF: 动力释放求解阶段的绘图状态 • INTFOR: 接触力和应力的绘图状态
LIVERMORE SOFTWARE TECHNOLOGY CORPORATION
输出提示
如果在一个已有先前分析结果的目录下运行一个新的 分析，选择下列之一: • 删除旧的‘二进制结果输出文件 • 在输出数据中选择另外的名字, 如在命令行中使用 g=filename • 创建一个新的目录运行要做的分析 (推荐) 请求足够多的输出状态以得到合理的时间历程结果 • 输出请求可能随分析的不同而不同，通常， 50 到 100 个状态
LIVERMORE SOFTWARE TECHNOLOGY CORPORATION
*SECTION_option
Option 包括 BEAM, DISCRETE, SHELL, SOLID 输入参数包含 ...
• • • •
梁单元，壳单元和体单元使用的计算公式 梁和壳的积分方法 壳的厚度 梁的截面性质
LIVERMORE SOFTWARE TECHNOLOGY CORPORATION
*ELEMENT_option
Option 指定单元的类型, 如 *element_shell 单元 ID Part ID (PID): 设定材料 ID, 性质 ID等到单元上 Node ID’s (单元连通性) Optional data can be provided to override shell thickness, orthotropic material angle, etc. on an element-by-element basis.
Ls-dyna Vehicle safety
Tongji University
George Shao
Ls-Dyna Vehicle Safety 关键字KEYWORD
HYPERMESH OASYS PRIMER FEMB ANSYS FEMAP LS-INGRID PATRAN
Tongji University
George Shao
Ls-Dyna Vehicle Safety 关键字KEYWORD
Tongji University
George Shao
Ls-Dyna Vehicle Safety 输出文件系统
二进制 d3plot 图形输出 dump 重启动文件 d3hsp 模型信息 message 求解信息 d3thdt 时间历程 runrsf 重启动
ASCII Glstat 总体统计能量 Matsum 基于component材料能量 Rcforce 接触面反作用力 Rbdout 刚体数据 Nodout 节点数据 Elout 单元结果
Tongji University
George Shao
Ls-Dyna Vehicle Safety 单元沙漏
Tongji University
George Shao
Ls-Dyna Vehicle Safety 材料MAT
• 弹性 MAT1 • 弹塑性 MAT24 • 正交各向异性/各向异性 • 超弹性（橡胶） • 泡沫 MAT57 • 复合材料 • 粘弹性 • 心脏/肺等生物材料
Tongji University
George Shao
Ls-Dyna Vehicle Safety
非线性动力分析
多物理场耦合分析
跌落分析
热分析
隐式求解
Tongji University
汽车碰撞分析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时，检查。
后面将逐一介绍碰撞分析中经常用到的控制卡片，并对每个卡片的作 用进行说明。
控制卡片使用规则
卡片相应的使用规则如下：
第一章 LS-DYNA关键字简介讲解
•1986 梁、壳、刚体
单一面接触
多平台支持
•1988
支持汽车行业应用
LS-DYNA3D
•1989
LSTC Full-Time 更多的选项和特性
ASME WAM’89 - San Francisco
-交通运输系统中的防撞性和乘
客安全
登录协同仿真时代
•1993 关键字格式
新的自动单一面接触
第一届LS-DYNA3D用户 年会
*COMPONENT假人组件 刚体假人. Integrated implicitly - separately from FE model. 由接触面控制相互作用. GEBOD 男性，女性，小孩 Hybrid III 5th, 50th and 95th 百分位成年人
登录协同仿真时代
*CONSTRAINED 结构体内和结构体之间
后处理： • »LS-POST, LS-TAURUS • »GLview, ANSYS,
HyperMesh, Fem-B
对于同一个物理系统总是有 多种方法构造有限元分析模 型 所有的模型不能完全重现物 理系统
登录协同仿真时代
输出文件
二进制文件 » d3dump01 » d3plot » d3plot01 » d3plot02 » …. » d3thdt » D3thdt01
登录协同仿真时代
*INCLUDE 把文件分成子文件 把子文件分成子-子文件 专门应用：冲压 part导入和转换（缩放、转动和平移）
登录协同仿真时代
*INITIAL 点火和动量 初始应力/应变 初始温度 初始汽车运动关系 初始速度
*INTERFACE • 研究大型结构中小部件的详细响应 •隔离关键部件 记录界面/边界的运动： »面 » 节点线 » 节点点
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关键字简介讲解
登录协同仿真时代
谢 谢！
Welcome to
ANSYS
China@
登录协同仿真时代
•1995 热分析
建立培训 LAB
交通运输系统中的防撞性 和乘 客安全
第三届LS-DYNA3D用户 年会
•1996 PC 版本（KBS2） 整合ANSYS •1997 •LS-DYNA（940版） •Eulerian/Lagrangian耦合 •自适应网格 •1999 •LS-POST •» LS-OPT
DYNA高级应用培训
1、关键字简介
发展历史
•1976 DYNA3D：Lawrence Livermore
国家实验室 John O. Hallquist 固体结构低速冲击 •1979 Cray-1 上运行 改进的滑动界面 Order of Magnitude faster •1981 新的材料模式 -炸药结构、土壤结构 射弹冲击
*CONTROL •终止：时间，能量或质量的变化 •时间步 •全局参数 » 沙漏类型 » 接触罚函数缩放因子 » 壳单元公式 » 阻尼 • 自适应网格， ALE， CFD，隐式
登录协同仿真时代
*DATABASE 控制输出 »ASCII 文件 » Binary文件
*DEFINE 与其他方式联动的输入: 曲线 限制几何范围的Boxes 局部坐标系 矢量 转换 (缩放，旋转，平移) » 必须在 *INCLUDE_TRANSFORM前 定义
*COMPONENT假人组件 刚体假人. Integrated implicitly - separately from FE model. 由接触面控制相互作用. GEBOD 男性，女性，小孩 Hybrid III 5th, 50th and 95th 百分位成年人
基于LS-DYNA的车门防撞梁碰撞仿真分析
引言据统计，汽车侧面碰撞事故导致的死亡率比较高（约占事故死亡总人数的34%）[1]。
主要原因为碰撞受力面积小，碰撞力集中，对乘员的伤害很大。
车门是侧面柱碰撞过程中主要的承力部件，提高车门防撞性的主要措施是在车门内部加装防撞梁[2-3]。
超高强钢在汽车领域逐渐得到应用，如22MnB5等。
国内如宝钢于2006年研发了热冲压用钢（冷轧B1500HS、热轧BRl500HS）[4]，可用于生产防撞部件、加强部件等。
本文采用有限元仿真软件LS-DYNA，对BR1500HS高强钢车门防撞梁进行碰撞仿真。
１ ＬＳ－ＤＹＮＡ的理论基础根据动力学问题的分析方法，离散化的结构动力方程为：()()()()()Mx t P t F t H t Cx t=−+−（1）式中，M为结构质量矩阵；C为阻尼矩阵；x(t)和x4(t)分别为节点的加速度向量和速度向量；P(t)、F(t)、H(t)分别为载荷向量、内力向量、沙漏阻力向量。
LS-DYNA采取变步长积分算法，积分步长必须小于某个临界值。
一般情况下，积分补偿取决于网格中的最小单元尺寸。
各单元类型的临界积分步长为：∆t e=α(l e/c) （2）式中，∆t e表示单元e的临界时间步长；α为时间步长因子，缺省值为0.9；l e为单元e的特征尺寸；c为纵波的波速。
LS-DYNA采用基于主、从表面的动态接触算法，其缺省接触算法是对称罚函数法，程序在每一时步对从节点和主节点分别进行穿透检查。
当前时间步若发生了穿透，则在从节点上施加法向接触力，按式（3）计算：f=k∆（3）i为接触刚度因子[5]。
单独对其碰撞变形情况进行仿真。
NACP中的侧面柱碰撞要求，在软件2mm，刚性柱的直径为254mm。
型导入ANSYS Workbench中进行前处理。
BR1500HS高强钢的密度为7850kg/m3，杨氏模量为2.1×105MPa，泊松比为0.3，屈服强度为1200MPa。
将圆柱体定义为刚性体。
汽车碰撞仿真LS-DYNA控制卡片关键字.答案
控制卡片使用规则
卡片相应的使用规则如下：
•������ 大部分的命令是由下划线分开的字符串, 如*control_hourglass
•������ 字符可以是大写或小写 •������ 在输入文件中，命令的顺序是不重要的(除了*keyword 和*define_table)
•������ 关键字命令必须左对齐，以*号开始
8. CONTROL_OUTPUT 用于设置输出参数。设置如下，按红色线框内设置，其余默认值。
控制卡片参数设置
8. CONTROL_PARALLEL 并行计算控制，最大可以使用4个CPU。
9. CONTROL_SHELL 壳单元控制
10. CONTROL_TERMINATION 计算终止控制卡片，控制计算终止时间，设置如下：
ISTUPD：单元厚度改变选项。该选项对所有壳单元变形有影响。
THEORY：壳单元使用的理论。（默认的是Belytschko-Tsay，面内单点积分，计算速度很快，采用Co-rotaional 应力更新，单元坐标系统置于单元中心，基于平面单元假定，建议在大多数分析中使用） BWC：针对Belytschko-Tsay单元的翘曲刚度。 MITER：平面应力塑性选项，默认为1。（运用于材料3，18，19和24）。 PROJ：在Belytschko-Tsay和Belytschko-Wong-Chiang单元中翘曲刚度投影方法。这个方法主要运用于显示分 析，如果是隐式分析，那此项无效 。 默认为0.
ORIEN：在初始化时可选择性的对接触面部分自动再定位。
控制卡片参数说明
ENMASS:接触单元被腐蚀的质量处理。0-节点被移除，1-体单元节点被保留，2-体单元壳单元节点被保留。 USRSTR：每个接触面分配的存储空间，针对用户提供的接触控制子程序。 USRFRC：每个接触面分配的存储空间，针对用户提供的接触摩擦子程序。
ls-dyna关键字命令介绍
LIVERMORE SOFTWARE TECHNOLOGY CORPORATION
*DATABASE_option 所谓的 ‘ASCII’ 输出文件
• 指定的时间历程输出 • 后续的幻灯片给出输出内容 • 输出格式的选择 (任何一种都可以用 LS-Prepost 处
理) • 真实的 ASCII ( 文本 ) 输出 • 压缩二进制输出 (binout)
在输出数据中选择另外的名字如在命令行中使用gfilename输出请求可能随分析的不同而不同通常50到100个状态的输出是足够的了对大模型来说使用d3thdt数据输出得到有选择的节点或单元的时间历程数据livermoresoftwaretechnologycorporationlivermoresoftwaretechnologycorporationdatabaseextentbinarydatabaseextentbinary对写入到二进制输出文件中的数据进行控制默认输出3个积分点midbottomtop的数据默认情况下lsdyna计算的弹性应变是不输出的各向异性正交异性的材料其应力可以在全局坐标系中输出也可以在材料坐标系中输出matorthotropicelastic除外livermoresoftwaretechnologycorporationlivermoresoftwaretechnologycorporationdatabaseextentbinarycontd控制每一个输出状态创建一个单独的d3plot文件选项ieverp梁的积分点输出数目beamip如果打开该选项可以对壳材料的沙漏能密度进行云图显livermoresoftwaretechnologycorporationlivermoresoftwaretechnologycorporationdatabaseoption所谓的ascii输出文件文本输出压缩二进制输出binoutlivermoresoftwaretechnologycorporationlivermoresoftwaretechnologycorporationascii输出文件airbagstatisticsdatabaseabstatinternalenergymassflowoutputmassflowrateboundarynodalforcesdatabasebndout对应于节点和刚体边界条件的节点力livermoresoftwaretechnologycorporationlivermoresoftwaretechnologycorporationascii输出文件discreteelementdatadatabasedeforc离散单元的力和力矩