《汽车碰撞分析LS_DYNA控制卡片设置.》由会员分享，可在线阅读，更多相关《汽车碰撞分析LS_DYNA控制卡片设置.（24页珍藏版）》请在金锄头文库上搜索。
1、汽车碰撞分析LS-DYNA 控制卡片的设置 作者:张远岭 2011-4-14 碰撞分析控制卡片包括求解控制和结果输出控制，其中KEYWORD、 CONTROL_TERMINATION、 DATABASE_BINARY_D3PLOT是必不可少的。其他一 些控制卡片如沙漏能控制、时间步控制、接触控制等则对计算过程进行控 制，以便在发现模型中存在错误时及时的终止程序。 后面将逐一介绍碰撞分析中经常用到的控制卡片，并对每个卡片的作 用进行说明。 控制卡片 大部分的命令是由下划线分开的字符串, 如*control_hourglass 字符可以是大写或小写 在输入文件中，命令的顺序是不重要的(除了*key
2、word 和*define_table) 关键字命令必须左对齐，以*号开始 第一列的“$”表示该行是注释行 输入的参数可以是固定格式或者用逗号分开 空格或者0 参数使用该参数的默认值 卡片相应的使用规则如下： 控制卡片使用规则 用hypermesh在LS-DYNA模板下，选择Analysis面板点击 control cards,选择相应卡片。 直接在key文件中输入 控制卡片的建立 控制卡片可通过以下方式建立： 下面介绍在hypermesh中给出碰撞分析中经常使用的卡片的参 数设置 1. KEYWORD KEY文件起始关键字。该卡片可不作任何设置。 2. TITLE 输入标题名称xxx。 3.
3、 CONTROL_ACCURACY 提高计算精度的控制卡片。设置INN值为2，其余默认，不起作用。 控制卡片参数设置 4. CONTROL_BULK_VISCOSITY 体积黏性项q，人工施加压力，用于处理应力波传播突变引起的不连续。Q1,Q2设置为默认,IBQ设为-1。 5. CONTROL_CONTACT 用于接触面的计算修改。控制摩擦惩罚系数、接触初始穿透检查、接触厚度等，设置见红框，其余默认。 控制卡片参数设置 6. CONTROL_ENERGY 用于控制沙漏能、耗散能、滑移能等的计算与否。全部设置成2，即全部参与计算。 7. CONTROL_HOURGLASS 用于沙漏的控制,通过施
4、加附加力来阻止沙漏变形。IHQ设置为5，QH默认为0.1。 8. CONTROL_OUTPUT 用于设置输出参数。设置如下，按红色线框内设置，其余默认值。 控制卡片参数设置 8. CONTROL_PARALLEL 并行计算控制，最大可以使用4个CPU。 9. CONTROL_SHELL 壳单元控制 10. CONTROL_TERMINATION 设置如下： 控制卡片参数设置 11. CONTROL_TERMINATION 计算终止控制卡片，控制计算终止时间。 12. DATABASE_BINARY_D3DUMP 设置如下： 13. DATABASE_BINARY_D3PLOT 设置如下： 控制
5、卡片参数设置 14. DATABASE_BINARY_D3THDT 设置如下： 15. DATABASE_BINARY_RUNRSF 设置如下： 16. DATABASE_BINARY_RUNRSF 设置如下： 控制卡片参数设置 17. DATABASE_BINARY_RUNRSF 设置如下： 控制卡片参数设置 控制卡片参数说明 *CONTROL_TERMINATION ENDTIM：强制终止计算时间，必选，默认0.0； ENDCYC：终止循环。在计算终止时间之前，程序达到指定循环次数即终止计算。 DTMIN：确定最小时间步长TSMIN的因子。TSMIN=DTMIN * DTSTART，其中D
6、TSTART为程序自动确定的 初始步长。当迭代步长小于TSMIN时，程序终止。 ENDENG：能量改变百分比，超过设定值则终止计算。默认0.0，不起作用； ENDMAS：质量变化百分比，超过设定值则终止计算。仅用于质量缩放DT2MS被使用时。默认0.0 ，不起作 用。 $ ENDTIM ENDCYC DTMIN ENDENG ENDMAS 150 0 0.0 0.0 0.0 *CONTROL_TIMESTEP(时间步长控制卡片) $ DTINIT TSSFAC ISDO TSLIMT DT2MS LCTM ERODE MSIST 0.0 0.9 2 0.0 -0.001 0 1 1 $ DT2
7、MSF DT2MSLC 计算所需时间步长时，要检查所有的单元。出于稳定性原因，用0.9（缺省）来 减小时间步：t = 0.9 l/c ，特征长度l，和波的传播速度c，都与单元的类型有关。 DTINIT：初始时间步长，如为0.0，由DYNA自行决定初始步长； TSSFAC：时间步长缩放系数，用于确定新的时间步长。默认为0.9，当计算不稳定时，可以减小该值，但同时 增加计算时间； ISDO：计算4节点壳单元时间步长的（不同的值对应特征长度的不同算法，推荐使用2，因为此选项可以获得 最大的时间步长，但有三角形单元存在时会导致计算不稳定）； TSLIMT：壳单元最小时间步分配 ,使单元的时间步长控制在
8、最小时间步长之上；只适用于使用 *mat_plastic_kinematic,*mat_power_law_plasticity*mat_strain_rate_dependent_plasticity,*mat_piecewise_linear_pla sticity等材料模型的壳单元，不建议使用该选项，因为使用DT2MS选项更好。 DT2MS：因质量缩放计算得到的时间步长。当设置为一个负值时，初始时间将不会小于TSSFAC*|DT2MS|。质 量只是增加到时间步小于TSSAFC*|DT2MS|的单元上。当质量缩放可接受时，推荐用这种方法。用这种方法时 质量增加是有限的，过多的增加质量会导致
9、计算终止。当设置为正值时，初始时间步长不会小于DT2MS。单 元质量会增件或者减小以保证每一个单元的时间步都一样。这种方法尽管不会因为过多增加质量而导致计算终 止，但更难以作出合理的解释。默认为0.0，不进行质量缩放； LCTM：限制最大时间步长的Load-curve,该曲线定义最大允许时间步长和时间的关系（可选择） ； ERODE：当计算时间步长小于TSMIN（最小时间步长）时体单元和t-shell被自动删除。 控制卡片参数说明 *CONTROL_SHELL $ WRPANG ESORT IRNXX ISTUPD THEORY BWC MITER PROJ 20 1 -1 0 2 2 1 全
10、局控制壳单元参数卡片 WRPANG：最大翘曲角度，默认20度； ESORT：程序自动把退化的四边形单元处理为C0三角形单元公式； IRNXX：单元法向更新开关，该选项只对Hughe_Liu，Belytschko-Wong-Chiang,Belytschko-Tsay等公式起作用 。（默认为-1，每次循环都重新计算单元方向）。 ISTUPD：单元厚度改变选项。该选项对所有壳单元变形有影响。 THEORY：壳单元使用的理论。（默认的是Belytschko-Tsay，面内单点积分，计算速度很快，采用Co-rotaional 应力更新，单元坐标系统置于单元中心，基于平面单元假定，建议在大多数分析中使用
11、） BWC：针对Belytschko-Tsay单元的翘曲刚度。 MITER：平面应力塑性选项，默认为1。（运用于材料3，18，19和24）。 PROJ：在Belytschko-Tsay和Belytschko-Wong-Chiang单元中翘曲刚度投影方法。这个方法主要运用于显示分 析，如果是隐式分析，那此项无效 。 默认为0. MS1ST：限制第一步的质量缩放，根据时间步确定质量向量一次。默认为0。 DT2MSF：决定最小时间步长的初始时间步长缩减系数，如果使用，DT2MS=-DT2MSF*t。 DT2MSLC：在显示分析中把DT2MS指定为时间的函数，使用load-curve定义。 控制卡片参
12、数说明 ROTASCL：为旋转单元质量定义一个缩放系数。（不太常用）。 INTGRD：通过厚度数值积分法则的默认壳单元。当积分点为1到2个的时候使用Gauss积分，当积分点从3个到10的时 候使用Lobatto积分，积分点为2个时，Lobatto法则非常不准，须用Gauss积分。 LAMSHT：薄壳理论开关 。0：不更新切应变修正；1：薄壳理论切应变修正 CSTYP6：第6种壳单元坐标系的选用。1：可变的局部坐标系（默认）；2：统一局部坐标系 （计算结果有偏差，但效率比较高）。 TSHELL：允许热传导通过有厚度的壳单元。 控制卡片参数说明 *CONTROL_HOURGLASS IHQ QH
13、5 0.1 沙漏控制卡片 IHQ：总体附加刚度或黏性阻尼方式选项； QH：沙漏能系数 ,超过0.15会导致计算不稳定。 *CONTROL_BULK_VISCOSITY Q1 Q2 TYPE 0.0 0.0 -1 体积黏度控制卡片 Q1：二次黏性系数 ； Q2：线性黏性系数 。 TYPE：体积黏性项。（当壳单元是类型2，10和16时使用-1）。 控制卡片参数说明 *CONTROL_CONTACT SLSFAC RWPNAL ISLCHK SHLTHK PENOPT THKCHG ORIEN ENMASS 0.1 1.0 2 2 1 0 1 0 USRSTR USRFRC NSBCS INTERM
14、 XPENE SSTHK ECDT TIEDPRJ 0 0 0 0 0.0 0 0 0 SFRIC DFRIC EDC INTVFC TH TH_SF TIPEN_SF 0.0 0.0 0.0 0.0 0.0 0.9 0.0 IGNORE FRCENG 1 SLSFAC：滑动接触惩罚系数 ，默认为0.1。当发现穿透量过大时，可以调整该参数； RWPNAL： 刚体作用于固定刚性墙时，刚性墙罚函数因子系数，为0.0时，不考虑刚体与刚性墙的作用，0时 ，刚体作用于固定的刚性墙，建议选择1.0； ISLCHK：接触面初始穿透检查，为0或1（默认）时，不检查。为2时，检查。 SHLTHK：在STS和NT
15、S接触类型中，是否考虑壳单元厚度对接触过程的影响，为0时，不考虑。为1时，除刚 体以外考虑。为2时，全部考虑。在单面接触、自动面接触等类型中，壳厚度的偏移总是被考虑。 PENOPT：对称刚度检查。如果两个接触物体的材料性质与单元大小的巨大差异，引起接触主面与从面之间接 触应力不匹配，可能导致计算不稳定和计算结果不切实际，这时可以调整该选项克服。为0或为1时，接触主面 与从节点中最小刚度值，为2时，主面刚度；为3时，从节点刚度。 THKCHG：在单面接触时考虑壳厚度的改变（默认时不考虑）。 ORIEN：在初始化时可选择性的对接触面部分自动再定位。 控制卡片参数说明 ENMASS:接触单元被腐蚀的质量处理。0-节点被移除，1-体单元节点被保留，2-体单元壳单元节点被保留。 USRSTR：每个接触面分配的存储空间，针对用户提供的接触控制子程序。 USRFRC：每个接触面分配的存储空间，针对用户提供的接触摩擦子程序。 NSBCS：接触搜寻的循环数（使用三维B
