汽车碰撞分析LS DYNA控制卡片设置

1、2011-4-14汽车碰撞分析ls-dyna 控制卡片的设置作者:张远岭 碰撞分析控制卡片包括求解控制和结果输出控制，其中keyword、control_termination、 database_binary_d3plot是必不可少的。其他一些控制卡片如沙漏能控制、时间步控制、接触控制等则对计算过程进行控制，以便在发现模型中存在错误时及时的终止程序。 后面将逐一介绍碰撞分析中经常用到的控制卡片，并对每个卡片的作用进行说明。 控制卡片大分的的命是是由下划分分的的字符, 如*control_hourglass的字可以是大写或大在输输入文中，的命的的序是不不重的(除了*keyword 和*defi

2、ne_table) 关的的命必必须对对，以*号分始一一的“$”表示该行是注释行输输的的参可以是是固定格写者用用号分分定写者0 的参用该的参的的默认卡片相应的用规则如由：控制卡片用规则用hypermesh在ls-dyna模板由，选择analysis面板点击 control cards,选择相应卡片。直接在key入文中输输控制卡片的建立控制卡片可通过以由方格建立：由面介绍在hypermesh中给出碰撞分析中经常用的卡片的的参设置1. keywordkey文件起始关键字。该卡片可不作任何设置。文件起始关键字。该卡片可不作任何设置。2. title输入标题名称输入标题名称xxx。3. control_

3、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_terminationendtim：强制终止计算时间，必选，默认0.0；endcyc：终止循环。在计算终止时间之前，程序达到指定循环次数即终止计算。dtmin：确定最小时间步长tsmin的因子。tsmin=dtmin *

6、dtstart，其中dtstart为程序自动确定的初始步长。当迭代步长小于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$

7、dt2msf 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_plasticity等材料模型的壳单元，不建议使用该选项，因为使用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全局控制壳单元参数卡片wrpang

10、：最大翘曲角度，默认20度；esort：程序自动把退化的四边形单元处理为c0三角形单元公式；irnxx：单元法向更新开关，该选项只对hughe_liu，belytschko-wong-chiang,belytschko-tsay等公式起作用 。（默认为-1，每次循环都重新计算单元方向）。 istupd：单元厚度改变选项。该选项对所有壳单元变形有影响。 theory：壳单元使用的理论。（默认的是belytschko-tsay，面内单点积分，计算速度很快，采用co-rotaional应力更新，单元坐标系统置于单元中心，基于平面单元假定，建议在大多数分析中使用）bwc：针对belytschko-ts

11、ay单元的翘曲刚度。 miter：平面应力塑性选项，默认为1。（运用于材料3，18，19和24）。 proj：在belytschko-tsay和belytschko-wong-chiang单元中翘曲刚度投影方法。这个方法主要运用于显示分 析，如果是隐式分析，那此项无效 。 默认为0.ms1st：限制第一步的质量缩放，根据时间步确定质量向量一次。默认为0。dt2msf：决定最小时间步长的初始时间步长缩减系数，如果使用，dt2ms=-dt2msf*t。dt2mslc：在显示分析中把dt2ms指定为时间的函数，使用load-curve定义。控制卡片的参说明rotascl：为旋转单元质量定义一个缩放系

12、数。（不太常用）。intgrd：通过厚度数值积分法则的默认壳单元。当积分点为1到2个的时候使用gauss积分，当积分点从3个到10的时候使用lobatto积分，积分点为2个时，lobatto法则非常不准，须用gauss积分。lamsht：薄壳理论开关 。0：不更新切应变修正；1：薄壳理论切应变修正cstyp6：第6种壳单元坐标系的选用。1：可变的局部坐标系（默认）；2：统一局部坐标系 （计算结果有偏差，但效率比较高）。tshell：允许热传导通过有厚度的壳单元。控制卡片的参说明*control_hourglassihq qh 5 0.1沙漏控制卡片ihq：总体附加刚度或黏性阻尼方式选项；qh：

13、沙漏能系数 ,超过0.15会导致计算不稳定。*control_bulk_viscosityq1 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 xpene ssthk ecdt tiedprj 0 0 0 0 0.0 0 0 0 s

14、fric 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和nts接触类型中，是否考虑壳单元厚度对接触过程的影响，为0时，不考虑。为1时，除刚体以外考虑。为2时，

15、全部考虑。在单面接触、自动面接触等类型中，壳厚度的偏移总是被考虑。penopt：对称刚度检查。如果两个接触物体的材料性质与单元大小的巨大差异，引起接触主面与从面之间接触应力不匹配，可能导致计算不稳定和计算结果不切实际，这时可以调整该选项克服。为0或为1时，接触主面与从节点中最小刚度值，为2时，主面刚度；为3时，从节点刚度。thkchg：在单面接触时考虑壳厚度的改变（默认时不考虑）。orien：在初始化时可选择性的对接触面部分自动再定位。 控制卡片的参说明enmass:接触单元被腐蚀的质量处理。0-节点被移除，1-体单元节点被保留，2-体单元壳单元节点被保留。usrstr：每个接触面分配的存储空

16、间，针对用户提供的接触控制子程序。usrfrc：每个接触面分配的存储空间，针对用户提供的接触摩擦子程序。nsbcs：接触搜寻的循环数（使用三维bucket分类搜索），推荐使用默认项。interm：间歇搜寻主面和从面接触次数。xpene：接触面穿透检查最大乘数，默认4.0。ssthk：在单面接触中是否使用真实壳单元厚度，默认0,不使用真实厚度。ecdt：时间步长内忽略腐蚀接触。*control_output$ npopt neecho nrefup iaccop opifs ipnint ikedit iflush 1 3 0 0 0.5 0 100 npopt：是否全部输出。（如果选1，那么坐

17、标系、单元链接、刚性墙定义和初始速度将不输出）。neecho：与npopt作用基本相同，只是可以屏蔽的输出选项不同。（如果选择3，则节点和单元都不输出到echo文件）。nrefup：beam单元的参考节点坐标是否更新，0不更新，1更新。iaccop:：从时间历程和节点速度得到平均加速度。时间历程文件“d3thdt”;速度文件“nodout”opifs：输出接触文件的时间间隔。ipnint：输出第一次循环所有单元的初始时间步长，默认0,输出100个时间步最小的单元。ikedit：在d3hsp输出间隔步的状态，如果输出glstat文件，忽略。iflush：针对i/o缓存的时间步间隔数，默认值500

18、0，如果缓存不是空的，计算非正常终止，输出文件将不完整。iprtf：在rbdout和matsum中默认输出。该选项是为了降低输出文件大小，排除一些不必要的输出。控制卡片的参说明*control_energy$ hgen rwen slnten rylen 2 2 2 2hgen:沙漏能计算选项。（默认选项1是不考虑沙漏能。选项2是计算沙漏能，并通过glstat和matsum两个ascii文件输出。）rwen：延迟能量耗散选项。slnten：接触滑移能耗散选项。（如果有接触那么这个选项设置成2）。rylen：阻尼能耗散选项。（默认选项1是不考虑阻尼能。选项2是计算阻尼能，并通过glstat文件输

19、出。）*control_accuracy$ osu inn pidosu 0 2 0osu:显示分析中旋转机械修正项 (second order objective stress updates)inn:壳和体单元的不变节点编号 pidosu:part id更新 *control_parallel$ ncpu numrhs const accu 4 0 1 1 ncpu：cpu使用数量。numrhs：内存数。const：并行计算。（默认设置为2，更快计算。）para：如果const=1，那么表示并行力集合。 控制卡片的参说明*database_abstat 输出安全带相关计算结果，如体积、压力

20、、质量、密度等。仅在安全气囊分析时使用。*database_avsflt 输出avs数据，与*database_extent_avs配合使用。*database_bndout 输出边界力及能量。与关键字*load_node_set（定义载荷边界）以及关键字*boundary_prescribed_motion_option（定义位移边界）等配合使用。*database_defgeo 输出变形几何信息。*database_deforc 输出离散单元作用力信息*database_elout 输出单元计算结果，如应力应变等。必须和*database_history_option配合使用，该命令定义输

21、出单元编号。*database_gceout 输出几何接触实体作用力。可以获得接触力和力矩。*database_glstat 输出模型整体信息，如动能、势能、沙漏能、阻尼能等计算结果。*database_h3out 输出假人相关计算结果。*database_option输出文本控制卡片dt：输出时间间隔；binary：输出文本格式。1为ascii，2为二进制，3两个都输出。$ dt binary 0.1 3控制卡片的参说明*database_jntforc 输出连接点信息。*database_matsum 输出与材料相关的信息，如动能、内能等。*database_movie 输出movie数据

22、，与*database_extent_movie配合使用。*database_mpgs 输出mpgs数据，与*database_extent_mpgs配合使用。*database_ncforc 输出节点界面力。定义接触时，必须将关键字*contact的spr和mpr参数设置成1。该命令将所有的节点力全部输出，因此文件会较大。*database_nodfor 输出节点力。该命令通过定义节点组输出节点力。必须和关键字*set_node_option及*database_nodal_force_group配合使用。该命令获得的节点力与ncforc文件中节点力完全一致，只是该文件可以由用户自定义输出节

23、点编号。*database_nodout 输出节点变形、速度、加速度等计算结果信息，和关键字 *database_history_node配合使用，该命令定义输出节点编号。*database_rbdout 输出刚体数据。*database_rcforc 输出合成界面力。该文件获得冲击力合力。*database_rwforc 输出刚性墙力。*database_sbtout 输出安全带相关计算结果。*database_secforc 输出cross类型（如焊点）的力，必须和关键字*database_cross_section_option配合使用。控制卡片的参说明*database_sleout

24、输出滑移能，和*control_energy配合使用。*database_spcforc 输出单点约束（spc）的反作用力。*database_sphout 输出sph数据。*database_ssstat 输出子系统数据。与*database_extent_ssstat配合使用。*database_swforc 输出节点约束反作用力（点焊和铆钉）（合成力）*database_tprint 输出热量。*database_trhist 输出轨迹线，与*database_tracer配合使用。控制卡片的参说明*database_binary_option 二进制文件的输出设置dt/cycl lcdt/nr beam npltc psetid istats tstart lavgdt/cycl：输出时间间隔/时间步长输出间隔。lcdt/nr：dump文件之间的时间间隔/重启文件数量（默认为1）beam：把弹簧阻尼单元的坐标和合力写进d3plot和d3part。npltc：dt=endtime/npltc，覆盖前面指定的dt。psetid：set_part id,只适用于d3part。istats：设置统计量级别