【www.myl5520.com--经典名言】
dyna材料模型失效准则确定
篇一:dyna中材料fail作用
有些材料类型中有关于失效准则的定义，但是也有些无失效准则的材料类型，这个时候需要额外的定义失效准则，与材料参数一块定义材料特性。用到*mat_add_erosion关键字，对于这个关键字有几个需要注意的地方。
1、材料的通用性准则：
材料通常为拉破坏或者剪切破坏，静水压是以压为正，拉为负，所以静水压破坏就是给出最小的承受压力，当然需要小于0（即压力），如果静水压小于该值，则材料破坏。相反，应力则是以压为负，拉为正，故最大主应力或最大等效应力或最大剪应力破坏等等都是给出最大的应力极限，当然大于0，如果拉应力大于该值，则材料破坏，无论是MAT_ADD_EROSION意：屈服不是失效。
2、单元失效模拟的功能与目的
单元删除功能本身是为了克服有限元本身的缺陷的一项方法，由于有限元本身就是基于连续介质力学的，而在连续介质理学中，所研究的物体需要是连续的，既物质域在空间中连续。在这样的理论假设框架下，单元本身是不会消失的。然而在实际情况下，由于损伤断裂的存在，势必会使得一些单元消失或者完全的失效，所以为了能够模拟这种情况，DYNA提供了单元失效功能。
破坏、失效、断裂，都是工程性的概念，它表示在达到某一准则后，结构、构件、或者构件中的某一部分，从结构中退出工作，不再影响整体结构的受力。而从有限元概念上说，对上述机制的模拟，基本手段都是一样的，就是当满足某一指标（比如某个应变大小）后，将一个单元或者一个积分点的质量、刚度和应力、应变都设为零（或者非常接近与零），这样它在整体结构计算中就不再发挥作用，进而实现了退出工作机制的模拟。所以，无论是把纤维模型中的某个纤维、或者分层壳模型中的某一层、或者实体模型中的某个积分点，或者结构中的某个单元，让其不再参与整体结构计算，都可以达到模拟破坏退出工作的目的。而所谓单元生死技术，是上述基本概念在有限元程序中的一个“打包”应用。它除了让单元不再参与计算外，一般还有一个重要的附加功能，就是对仅和“被杀死”单元相连的“孤岛”节点，让其自由度不再参与整体结构计算，以减少计算困难。而后来有限元程序的前后处理又不断改进，可以做到在后处理里面“看不到”已杀死的单元，这样就显得更加真实。但正因为这些包装，使得很多人反而忘记了所谓单元生死技术的基本概念。.
所以，不要被单元生死吓到，即便是有限元程序不提供“单元生死”功能，通过适当的设计单元质量、刚度和应力应变矩阵，也可以实现单元生死同样的效果。至于构件的部分或局部破坏（诸如钢筋的断裂），更是有多种实现方法，使用者可以灵活掌握。
3、关于关键字参数
这个参数有两行参数，第一行：MID（MID待失效的材料编号）excl（排除数字，任意假设）；第二行：PFAIL（失效压力）SIGPI（失效主应力）SIGVM（失效等效应力，一般指抗拉强度）EPSPI（失效主应变）EPSSH（失效剪应变）SIGTH（极限应力）IMPULSE（失效应力冲量）FAILTM（失效时间） 其中excl为排除数字，这个数字可以任意定义，如果第二行某个参数和这个数据相同，那么该参数定义的失效准则就被忽略。（第二行可以定义很多准则）。不选用其它失效准则不能留空，必须要填排除数字。
关于PFAIL关键字的说明：各个方此关键字表示物体的静水压破坏，即向受到压力时的破坏准则，其中压为正，拉为负，一般材料尤其是混凝土材料都是拉伸破坏，故此参数一般定义为负数，对于大小比较的是代数值的大小，因此当低于此准则即拉伸力超过允许数值，材料即宣告破坏（类似抗压强度）。当实际的静水压力（其实应该是拉力）小于此值（代数大小），材料即宣告破坏。
除最后一个是关于时间的破坏准则外，其余的六个破坏准则都是正数，表示拉力，当计算的数值大于此值时材料失效删除。
4、关于材料失效 压缩破坏在这个关键字中无法体现，要想施加压缩破坏准则，必须要自己定义关键字参数，即进行二次开发。另外，需要说明的是，动态破坏的基本特性是时率相关性和损伤积累性，损伤这一块，特别是微观上真实的损伤，而不是宏观上的唯象损伤，DYNA几乎是空白，所以就需要自定义材料了！ 另外，应力波的破坏形式有两种，即拉伸破坏和剪切破坏，很少有材料是压缩破坏的，因为还没有达到压缩破坏的阀值的时候由于可能由于泊松比导致的侧向拉力已经达到了极限，所以混凝土材料真正的压缩强度是多少没有人知道。
5、参数的使用范围
关键字的使用范围只是单点积分的2d和3d的实体单元
6、关于材料失效与裂纹
在DYNA中，材料一旦失效就被自动的删除，而结构之所以出现裂缝或者破碎，是因为结构单元中一部分单元失效，另一部分未失效，这些未失效的部分被孤立就形成了破碎。 裂纹的形成有两种方式，一种是定义单元失效准则（*mat_add_erosion关键字），这种在定义的时候裂纹部位的网格必须足够的密，否则大量单元失效对结算结果会有较大影响；另一种是定义节点约束失效形成裂纹，方法是单元之间不是通过共节点进行连接，而是相互独立的，通过定义失效约束进行连接，这种方法的问题在于建立模型的过程比较复杂。
ls_dyna中一些材料参数K文件
篇二:dyna中材料fail作用
1.紫铜(johnson_cook) EX=1.19 cm-g-us *MAT_JOHNSON_COOK 1 8.96000 0.46
0.900E-03 2.920E-03 0.310 0.250E-01 1.09 0.1356E+04 210 0.100E-05 0.383E-05 -9.00E+00 3.00 0.00 3.00 0.00 0.00 0.00 0.00
*EOS_GRUNEISEN
1 0.394 1.489 0.00 0.00 2.02 0.47 0.00 1.00
========================================================================================
2.高能炸药(LS-DYNA中无此材料模型，任取以模型代替，修改k文件) *MAT_HIGH_EXPLOSIVE_BURN
1 1.7870001 0.8390000 0.3400000 0.0000000 0.0000000 0.0000000 0.0000000 *EOS_JWL
1 5.8140002 6.8010E-2 4.1000000 1.0000000 0.3500000 0.0900000 1.0000000
1 M8--------------------------------------------------------------------------------------- 或者：
*MAT_HIGH_EXPLOSIVE_BURN
1 1.6400000 0.6930000 0.2700000 0.0000000 0.0000000 0.0000000 0.0000000 *EOS_JWL
1 3.7400000 3.2300e-2 4.1500001 0.9500000 0.3000000 0.0700000 1.0000000
======================================================================================== 3.高导无氧铜 (Steinberg):
*MAT_STEINBERG
2 8.93 0.477 0.120E-02 36.0 0.450 0.00 0.640E-02 2.83 2.83 0.377E-03 0.100E-02 63.5 0.179E+04 2.02 1.50 -9.00 3.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 *EOS_GRUNEISEN
2 0.394 1.49 0.00 0.00 2.02 0.470 0.00 1.00
-------------------------------------------------------------------------------------- (plastic_kinematic):
*MAT_PLASTIC_KINEMATIC
1 8.93 1.17 0.350000 0.400E-02 0.100E-02 1.00 0.00 0.00 0.00
========================================================================================dyna中材料fail作用。
4.钢(johnson_cook):cm-g-us EX=2.0
*MAT_JOHNSON_COOK 1 7.83000 0.77
7.920E-03 5.100E-03 0.260 0.140E-01 1.03 0.1793E+04 294 0.100E-05 0.477E-05 -9.00E+00 3.00 0.0 5.00 0.00 0.00 0.00 0.00
*EOS_GRUNEISEN
1 0.4569 1.490 0.00 0.00 2.17 0.46 0.00 1.00
--------------------------------------------------------------------------------------- *MAT_PLASTIC_KINEMATICdyna中材料fail作用。
1 7.83 2.07 0.300000 0.400E-02 5.00E-02 1.00 0.00 0.00 0.00
======================================================================================== 5.混凝土(solid_concrete) *MAT_SOIL_CONCRETE
2 2.500 1.250E-01 1.667E-01 1 2 3 4 -2.00E-05 0.0 0.3 1 *DEFINE_CURVE
1 0 0.0000000 0.0000000 0.0000000 0.0000000 0.0 0.0
0.0953 0.230e-01 0.1823 0.650e-01 0.2624 1.300e-01 0.3365 1.750e-01 0.4055 2.250e-01 0.4700 3.500e-01 *DEFINE_CURVE
2 0 0.0000000 0.0000000 0.0000000 0.0000000 0.0 0.0dyna中材料fail作用。
0.400e-02 0.240e-02 0.800e-02 0.550e-02 1.200e-02 0.950e-02 1.400e-02 1.260e-02 1.800e-02 1.550e-02 1.000e-01 1.550e-02 *DEFINE_CURVE
3 0 0.0000000 0.0000000 0.0000000 0.0000000 0.00 0.00 1.500e-04 0.00 1.000e-03 0.04 1.000e-02 0.49 1.000e-01 4.90 *DEFINE_CURVE
4 0 0.0000000 0.0000000 0.0000000 0.0000000 0.00 0.01 1.500e-04 0.01 1.000e-03 0.05 1.000e-02 0.50 1.000e-01 5.00
================================================================6.土壤(泡沫材料模型)
*MAT_SOIL_AND_FOAM
2 1.80 1.601E-04 1.328E+02 3.300E-03 1.310E-07 0.1232 0.0 0.0 0.0
0.0 0.050 0.090 0.11 0.15 0.19 0.21 0.22 0.25 0.30
0.0 3.420E-02 4.530E-02 6.760E-02 1.270E-01 2.080E-01 2.710E-01 3.920E-01 5.660E-01 1.230E+00
============================================== *MAT_PLASTIC_KINEMATIC 1 7.82 2.1 0.284 0.0100 0.04 1.0 0.000E+00 0.000E+00 0.000E+00 *MAT_PLASTIC_KINEMATIC
1 7.80 2.1 0.300 0.023 0. 0240. 1.00 0 0 0
PLASTIC_KINEMATIC
1 7.80 2.07 0.300000 0. 235 0.0700 1.00 0.00 0.00 0.00
=================================================== B4C陶瓷
*MAT_JOHNSON_HOLMQUIST_CERAMICS
2 2510e-6 197e+3 0.927 0.7 0.005 0.85 0.67 1.0 260 0.2 19e+3 8.71e+3 1.0 0.001 0.5 233e+3 -593e+3 28e+5 0.1*MAT_GEOLOGIC_CAP_MODEL
1,2700,1.4E10,1.1E10,2.7E7,0.11,8.0E6,1.4E-7 4.43,4.6E-10,0.42,1.1E8
$40号砼,单位：长度mm,时间msec,质量g,力N,压强Mpa *MAT_SCHWER_MURRAY_CAP_MODEL 1,2.3E-3,1.048E4,1.168E4,0,0,1
190,0,184.2,2.5E-3,0.999,0.7,2.5,2.5E3 5,100,1,0,0.999,0.3,0.94
5.0E-2,2.5E-4,3.5E-7,23,0,1,300,13.8
1,0,0,0
0.747,3.3E-4,0.17,5E-2,0.66,4E-4,0.16,5E-2 $Define Steel 4340 4340号钢 *MAT_JOHNSON_COOK
1, 7.8400E+00, 0.759E-00, 2.000E+00, 0.320E+00, 0.000E+00
0.793E-02, 0.510E-02, 2.60E-01, 1.400E-02, 1.030E+00, 1.793E+03, 2.930E+02, 1.000E-05
4.400E-06, 2.000E-02, 2.000E+00, 0.000E+00, 0.80E+00, 2.100E+00, -0.05E+00, 2.000E-03
0.610E+00 *EOS_GRUNEISEN
1, 4.578E-01, 1.330E+00, 0.000E+00, 0.000E+00, 1.670E+00, 0.430E+00, 0.000E+00 1.000E+00
$ SI unit :cm-g-microsecond $ unit conversion factor
$ 1 psi=6895 pa, 1 dyn=1e-5 n,1 bar=1e5 pa,1 g/cm**3=0.0361 lb/in**3,1 lb=0.454kg