汽车安全性分析卡片诠释
汽车安全性分析卡片诠释
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CONTROL_TIMESTEP
$ dtinit tssfac isdo tslimt dt2ms lctm erode ms1st
0.0 0.90000 0 .2223E-5 -.1000E-5 0 0 0
tslimt:设置SHELL单元的最小时间步长，但是该选项仅仅适用于材料：*MAT_PLASTIC_KINDMATIC
*MAT_POWER_LAW_PLASTICITY *MAT_PIECEWISE_LINEAR_PLASTICITY等。
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
在PART定义中考虑摩擦系数（用于转向管柱，见图steer_system）
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PART_CONTACT_PRINT
Halter 1 3.0 mm
$ pid secid mid eosid hgid grav adpopt tmid
500011 1 1 0 0 0 0 0
$ fs fd dc vc optt sft ssf
.15000001 .08000000 1.0000000 0.000E+00 2.9000001 0.000E+00 1.0000000
$ prbf
0
该卡片的目的是碰撞分析中，只考虑局部零件的摩擦系数
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
汽车纵梁材料定义
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MAT_PIECEWISE_LINEAR_PLASTICITY_TITLE （24＃）
LINE_PLASTICITY
$ mid ro e pr sigy etan fail tdel
9 .7830E-8 207000.0 0.28000 200.000 0.0 0.0 0.0
$ c p lcss lcsr vp
0.0 0.0 0 0 0.0
$ eps1 eps2 eps3 eps4 eps5 eps6 eps7 eps8
0.0 0.03090 0.40900 0.05000 0.15100 0.30100 0.70100 0.91000
$ es1 es2 es3 es4 es5 es6 es7 es8
210.000 300.000 314.000 325.000 390.000 438.000 505.000 527.000
*MAT_PIECEWISE_LINEAR_PLASTICITY_TITLE
M0511213
$ mid ro e pr sigy etan fail tdel
511213 .78900E-8 210000.0 0.30000 400.000 0.0 0.0 0.0
$ c p lcss lcsr vp
80.000 4.500 2000008 0 1.000
$ eps1 eps2 eps3 eps4 eps5 eps6 eps7 eps8
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
$ es1 es2 es3 es4 es5 es6 es7 es8
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
分段线性塑性材料模型，LCSS：定义应力应变曲线（表），LCSR：应变率比例曲线。
c、p考虑应变率模型的参数。
eps1－eps8 定义至少8个点的等效塑性应变。
es1－es8， 定义相对应的屈服应力值。
*DEFINE_CURVE_TITLE
$ lcid sidr sfa sfo offa offo dattyp
2000008 0 1.000 1.000 0.0 0.0 0
$ abscissa (time) ordinate (value)
0.000000E+00 0.400000E+03
0.390000E-01 0.426000E+03
0.770000E-01 0.453000E+03
0.122000E+00 0.480000E+03
0.166000E+00 0.507000E+03
0.199000E+00 0.550000E+03
0.231000E+00 0.600000E+03
0.247000E+00 0.670000E+03
0.100000E+01 0.680000E+03
应力（MPa）
定义不同应变率下的应力应变曲线的材料模型
*MAT_PIECEWISE_LINEAR_PLASTICITY_TITLE
GENERAL-STEEL
$ mid ro e pr sigy etan fail tdel
1 .78500E-8 205000.0 0.28000 0.0 0.0 0.0 0.0
$ c p lcss lcsr vp
0.0 0.0 1400 0 0.0
$ eps1 eps2 eps3 eps4 eps5 eps6 eps7 eps8
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
$ es1 es2 es3 es4 es5 es6 es7 es8
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
*DEFINE_TABLE_TITLE
STRESS-STRAIN
$ tbid
1400
$ value
0.000000E+00
0.150000E-04
0.150000E+01
0.300000E+02
0.500000E+03
*DEFINE_CURVE_TITLE
C1-STREES-STRAIN
$ lcid sidr sfa sfo offa offo dattyp
1401 0 1.000 1.000 0.0 0.0 0
$ abscissa (time) ordinate (value)
0.000000E+00 0.265000E+03
0.100000E-01 0.280800E+03
0.250000E-01 0.299000E+03
0.550000E-01 0.320000E+03
0.100000E+00 0.344900E+03
0.117700E+00 0.354600E+03
0.130000E+00 0.359500E+03
0.140000E+00 0.364300E+03
0.150000E+00 0.367600E+03
0.160000E+00 0.369200E+03
0.170000E+00 0.370000E+03
0.190000E+00 0.370000E+03
0.220000E+00 0.370000E+03
0.250000E+00 0.370000E+03
0.450000E+00 0.370000E+03
*DEFINE_CURVE_TITLE
C2-STREES-STRAIN
$ lcid sidr sfa sfo offa offo dattyp
1402 0 1.000 1.000 0.0 0.0 0
$ abscissa (time) ordinate (value)
0.000000E+00 0.290000E+03
0.100000E-01 0.302700E+03
0.250000E-01 0.315000E+03
0.550000E-01 0.335000E+03
0.100000E+00 0.361900E+03
0.117700E+00 0.372400E+03
0.130000E+00 0.376500E+03
0.140000E+00 0.377300E+03
0.150000E+00 0.378900E+03
0.160000E+00 0.380000E+03
0.170000E+00 0.380000E+03
0.190000E+00 0.380000E+03
0.220000E+00 0.380000E+03
0.250000E+00 0.380000E+03
0.450000E+00 0.380000E+03
*DEFINE_CURVE_TITLE
C3-STREES-STRAIN
$ lcid sidr sfa sfo offa offo dattyp
1403 0 1.000 1.000 0.0 0.0 0
$ abscissa (time) ordinate (value)
0.000000E+00 0.321400E+03
0.100000E-01 0.330000E+03
0.250000E-01 0.345700E+03
0.550000E-01 0.370800E+03
0.100000E+00 0.400800E+03
0.117700E+00 0.410000E+03
0.130000E+00 0.414600E+03
0.140000E+00 0.417800E+03
0.150000E+00 0.418600E+03
0.160000E+00 0.420000E+03
0.170000E+00 0.420000E+03
0.190000E+00 0.420000E+03
0.220000E+00 0.420000E+03
0.250000E+00 0.420000E+03
0.450000E+00 0.420000E+03
*DEFINE_CURVE_TITLE
C4-STREES-STRAIN
$ lcid sidr sfa sfo offa offo dattyp
1404 0 1.000 1.000 0.0 0.0 0
$ abscissa (time) ordinate (value)
0.000000E+00 0.380000E+03
0.100000E-01 0.390000E+03
0.250000E-01 0.400000E+03
0.550000E-01 0.420000E+03
0.100000E+00 0.432400E+03
0.117700E+00 0.435000E+03
0.130000E+00 0.438000E+03
0.140000E+00 0.440000E+03
0.150000E+00 0.440000E+03
0.160000E+00 0.440000E+03
0.170000E+00 0.440000E+03
0.190000E+00 0.440000E+03
0.220000E+00 0.440000E+03
0.250000E+00 0.440000E+03
0.450000E+00 0.440000E+03
*DEFINE_CURVE_TITLE
C5-STREES-STRAIN
$ lcid sidr sfa sfo offa offo dattyp
1405 0 1.000 1.000 0.0 0.0 0
$ abscissa (time) ordinate (value)
0.000000E+00 0.475000E+03
0.100000E-01 0.485000E+03
0.250000E-01 0.495000E+03
0.550000E-01 0.510000E+03
0.100000E+00 0.518400E+03
0.117700E+00 0.520000E+03
0.140000E+00 0.520000E+03
0.150000E+00 0.520000E+03
0.160000E+00 0.520000E+03
0.170000E+00 0.520000E+03
0.190000E+00 0.520000E+03
0.220000E+00 0.520000E+03
0.250000E+00 0.520000E+03
0.450000E+00 0.520000E+03
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MAT_DAMPER_NONLINEAR_VISCOUS_TITLE
DAMPER
$ mid lcdr
30 60
非线性粘性阻尼材料，
*MAT_SPRING_NONLINEAR_ELASTIC_TITLE
SPRING
$ mid lcd lcr
31 29 0
非线性弹性弹簧材料，
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MAT_MOONEY-RIVLIN_RUBBER_TITLE （27号材料）
M_UTREAD
$ mid ro pr a b ref
210 .2040E-8 0.49900 5.541 0.98530 0.0
$ sgl sw st lcid
0.0 0.0 0.0 0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
吸能泡沫材料定义
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MAT_LOW_DENSITY_FOAM （57号材料 ）
5635349.1310E-11 30.6 20003011.0000E+20 0.01 0.0 0.0
8.0 0.0 0.0 0.0 0.0 0.0 0.0
高度可压缩性泡沫，主要用于坐垫、汽车前部吸能块、侧碰假人补丁。
*DEFINE_CURVE
2000301 0 1.0 1.0 0.0 0.0 0
0.0 0.0
0.05 0.50999999
0.1 0.75999999
0.15000001 0.89999998
0.2 0.97000003
0.25 1.03999996
0.30000001 1.12
0.34999999 1.20000005
0.40000001 1.30999994
0.44999999 1.48000002
0.5 1.72000003
0.55000001 2.01999998
0.60000002 2.38000011
0.64999998 2.9000001
0.89999998 5.5
0.99000001 185.5
应力应变曲线 横坐标是应变、纵坐标是应力。
*MAT_FU_CHANG_FOAM （ 83#材料）
1 8.800E-04 10.900000 1.0000000 0.000E+00 0.000E+00 0.000E+00 0.000E+00
0.000E+00
1.128E+03 1.0000000 3.0000000 1.0000000 0.5100000 0.5000000 1.9000000 1.2000000
800.00000 0.1490000 4.200E-03 3.0000000 2.0000000 1.430E-02 9.400E-02
MAT-FU-CHANG-FOAM是中低密度泡沫材料模型，考虑了应变率的影响
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
防止负体积措施
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
像*MAT_LOW_DENSITY_FOAM这种软材料，在高速碰撞下，变形很大，容易发生负体积。
*CONTACT_INTERIOR
$ psid
5
在泡沫单元内部定义接触，防止负体积产生。该算法自动生成实体表面，并形成接触。
把所有可能出现负体积的实体单元组成一个PART SET组。
*MAT_NULL_TITLE
M0563419
$ mid ro pc mu terod cerod ym pr
563419 .21040E-8 0.0 0.0 0.0 0.0 69000.00 0.30000
*SECTION_SHELL_TITLE
S0563419
$ secid elform shrf nip propt qr/irid icomp setyp
563419 2 1.000 3.000 0.0 0.0 0 1
$ t1 t2 t3 t4 nloc marea
0.50000 0.50000 0.50000 0.50000 0 0.0
在体单元的外表面构建SHELL单元，(更为可靠的做法是，在每个实体单元表面上构建SHELL单元)单元材料为MAT_NULL。
定义整车单面接触时，要包括新建的SHELL单元，不要包括实体单元。
时间步长可以调小，
!!!!!!!!!!!!!!!!!!!!!!!!!!!!
过小的实体单元（<5mm）在高速碰撞中很容易发生负体积，尽量避免过小单元。
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
挡风玻璃材料定义
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MAT_LAMINATED_GLASS_TITLE （32＃）
M0507121
$ mid ro eg prg syg etg efg ep
507121 .33190E-8 71000.00 0.21800 75.000 10000.00 0.00400 500.000
$ prp syp etp
0.45000 2.500 5.000
$ f1 f2 f3 f4 f5 f6 f7 f8
0.0 0.0 1.000 1.000 1.000 0.0 0.0 0.0
$ f1 f2 f3 f4 f5 f6 f7 f8
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
$ f1 f2 f3 f4 f5 f6 f7 f8
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
$ f1 f2 f3 f4 f5 f6 f7 f8
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
蜂窝材料定义（汽车散热风扇的蜂窝材料）
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MAT_HONEYCOMB_TITLE （26＃）
M0530404
$ mid ro e pr sigy vf mu bulk
530404 .14750E-9 2070.000 0.30000 140.000 0.20000 0.05000 0.0
$ lca lcb lcc lcs lcab lcbc lcca lcsr
2000012 2000012 2000012 2000012 0 0 0 0
$ eaau ebbu eccu gabu gbcu gcau aopt
20.700 20.700 20.700 2.100 2.100 2.100 0.0
$ xp yp zp a1 a2 a3
0.0 0.0 0.0 0.0 0.0 0.0
$ d1 d2 d3 tsef ssef
0.0 0.0 0.0 0.0 0.0
*DEFINE_CURVE_TITLE
$ lcid sidr sfa sfo offa offo dattyp
2000012 0 1.000 1.000 0.0 0.0 0
$ abscissa (time) ordinate (value)
0.160000E+00 0.200000E+01
0.200000E+00 0.897000E+00
0.800000E+00 0.896000E+00
0.100000E+01 0.000000E+00
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
玻璃和引擎盖的粘胶材料
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MAT_PLASTICITY_COMPRESSION_TENSION_TITLE （124＃）
M0507123
$ mid ro e pr c p fail tdel
507123 .48210E-9 350.000 0.39000 0.0 0.0 0.0 0.0
$ lcidc lcidt
2010005 2010006
$ pc pt
1.000 1.000
lcidc：材料压缩时，应力和应变的曲线
lcidt：材料拉伸时，应力和应变的曲线
*DEFINE_CURVE_TITLE
$ lcid sidr sfa sfo offa offo dattyp
2010005 0 1.000 1.000 0.0 0.0 0
$ abscissa (time) ordinate (value)
0.550000E-01 0.100000E+02
0.900000E-01 0.110000E+02
0.170000E+00 0.120000E+02
0.270000E+00 0.130000E+02
0.330000E+00 0.140000E+02
0.400000E+00 0.155000E+02
0.425000E+00 0.160000E+02
0.100000E+04 0.157000E+04
*DEFINE_CURVE_TITLE
$ lcid sidr sfa sfo offa offo dattyp
2010006 0 1.000 1.000 0.0 0.0 0
$ abscissa (time) ordinate (value)
0.000000E+00 0.160000E+02
0.500000E-01 0.100000E-03
0.100000E+01 0.100000E-03
该材料970版只适合Solid单元，是目前唯一能区分拉伸和压缩状态的材料模型。971会增加SHELL单元。
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
用接触来模拟粘胶
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CONTACT_TIED_SHELL_EDGE_TO_SURFACE_ID
$ cid
54
$ ssid msid sstyp mstyp sboxid mboxid spr mpr
6 7 2 2 0 0 0 0
$ fs fd dc vc vdc penchk bt dt
0.15000 0.08000 1.000 0.0 0.0 0 0.0 1.000
$ sfs sfm sst mst sfst sfmt fsf vsf
1.000 1.000 0.0 0.0 1.000 1.000 1.000 1.000
$ soft sofscl lcidab maxpar edge depth bsort frcfrq
1 0.10000 0 1.025 0.0 2.000 0 1
$ penmax thkopt shlthk snlog isym i2d3d sldthk sldstf
0.0 0 0 0 0 0 0.0 0.0
$ igap ignore
2 0
粘胶用实体单元（一层）模拟，为从接触
和粘胶接触的所有壳单元为主接触。
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
整车正碰定义
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*RIGIDWALL_GEOMETRIC_FLAT_ID
$ rwid
0FRONT NCAP FRONT END WALL
$ nsid nsidex boxid
0 0 0
$ xt yt zt xh yh zh fric
-960.283 -972.272 -290.000 -860.283 -972.272 -290.000 0.50000
$ xhev yhev zhev lenl lenm
-960.283 -872.272 -290.000 0.0 0.0
定义整车所有的节点为从节点，并与刚性墙发生接触，与刚性墙的接触的必须是NODE SET 组，
nsid＝0：表示所有的节点；
其实，可以定义BOX来限定矩形框内的节点和刚性墙发生接触，BOX的范围可以是车的前部，因为在整车正碰中，
与刚性墙发生接触的只是汽车前部的部分节点。
*RIGIDWALL_GEOMETRIC_FLAT_ID
$ rwid
0FRONT NCAP VEHICLE BOTTOM WALL
$ nsid nsidex boxid
0 0 4
$ xt yt zt xh yh zh fric
-1116.859 -1172.272 -290.000 -1116.859 -1172.272 -190.000 0.0
$ xhev yhev zhev lenl lenm
-1016.859 -1172.272 -290.000 0.0 0.0
定义box内所有的节点为从节点并与路面接触：
boxid=4：表示该矩形框内所有的节点有效。因为与路面接触的只有轮胎的部分节点，所有有必要定义BOX，目的
是定义只定义部分节点与路面发生接触，这样可以减少接触搜寻，提高计算效率。
*DEFINE_BOX_TITLE
$ boxid xmn xmx ymn ymx zmn zmx
RIGID WALL
4-.1117E+040.4153E+04-.1172E+040.1172E+04-.2900E+03-.9000E+02
定义矩形框X Y Z的最大最小值。
*INITIAL_VELOCITY
$ nsid nsidex boxid irigid
0 0 1 0
$ vx vy vz vxr vyr vzr
-15650.00 0.0 0.0 0.0 0.0 0.0
*DEFINE_BOX_TITLE
DEFINE_BOX: 1
$ boxid xmn xmx ymn ymx zmn zmx
1-.9200E+030.4000E+04-.1000E+040.1000E+04-.3000E+030.1400E+04
定义整车的初始速度，以及节点范围。
该卡片nsid＝0：表示所有的节点，
在整车碰撞分析中，最好定义nsid＝0，这样做非常好，因为碰撞分析定义后，如果修改了计算模型，可以不必再重新修改卡片。
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
定义焊点
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1，通过beam方式：
*MAT_SPOTWELD_TITLE
M0563407
$ mid ro e pr sigy et dt tfail
563407 .23530E-6 210000.0 0.30000 .9999E+9 200.000 .10500E-5 0.0
$ efail nrr nrs nrt mrr mss mtt nf
0.0 400.000 0.0 0.0 0.0 0.0 0.0 4.000
！！！100号材料。在材料中定义焊点的失效
nrr=400 opt=0定义轴向力达到400N时，焊点失效。
opt=1时，定义应力失效，opt=0时为合力失效。
nrs：定义剪切失效。
sigy:初始屈服应力。
et:硬化指数。
*SECTION_BEAM_TITLE
S0563407
$ secid elform shrf qr/irid cst scoor nsm
563407 9 1.000 1.000 1.000 0.0 0.0
$ ts1 ts2 tt1 tt2 print
5.100 5.100 0.0 0.0 0.0
！！！！第9种单元算法。
*CONTACT_TIED_SHELL_EDGE_TO_SURFACE_ID
$ cid
54
$ ssid msid sstyp mstyp sboxid mboxid spr mpr
6 7 2 2 0 0 0 0
$ fs fd dc vc vdc penchk bt dt
0.15000 0.08000 1.000 0.0 0.0 0 0.0 1.000
$ sfs sfm sst mst sfst sfmt fsf vsf
1.000 1.000 0.0 0.0 1.000 1.000 1.000 1.000
$ soft sofscl lcidab maxpar edge depth bsort frcfrq
1 0.10000 0 1.025 0.0 2.000 0 1
$ penmax thkopt shlthk snlog isym i2d3d sldthk sldstf
0.0 0 0 0 0 0 0.0 0.0
$ igap ignore
2 0
如果不是点对点，则通过定义beam单元与SHELL单元的接触来连接。
*ELEMENT_BEAM_PID
此卡片仅对焊点BEAM单元有效，如果没有激活PID，则节点向距离最近的SEGMENT投影。如果激活了PID，则向
卡片中提供的两个PART投影。
卡片中的第二个节点是可选的，如果只有第一个节点，则给节点分别向两边投影，因此，最好，该节点能在两个投影面的中间。
2， 通过SPOTWELD
*CONSTRAINED_SPOTWELD_ID
$VPG_SPOTWELD 591377
$ cnsid
363559
$ N1 N2 SN SS N M TF EP
341322 341616 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00
可以定义失效
3，通过刚性连接，不能考虑失效。
*CONSTRAINED_NODAL_RIGID_BODY
$VPG_NODAL_RB 591378 368575
$ pid cid nsid pnode iprt
4022
4，其他焊接
*CONSTRAINED_GENERALIZED_WELD
