MAT_ADD_EROSION 关键字参数解释(期待高手完善)
本帖最后由 alinheng 于 2011-4-3 13:45 编辑*MAT_ADD_EROSION
$# mid excl mxpres mneps effeps voleps numfip ncs
1 1234.0000 1.0000E+8-1.0000E-3 1.0000E-3-1.0000E-3
$#mnpres sigp1 sigvm mxeps epssh sigth impulse failtm
-2.5000E+7 1.0000E+7 1.0000E+7 1.0000E-3 1.0000E-3
以下是对有关项的解释。括号内为LS-PrePost 查看结果时对应的：Fcomp或 History中的选项。近几个版本的LS-DYNA 中 *MAT_ADD_EROSION 关键字
格式改变较大，这个应该是971的格式。
mid: 要加侵蚀关键字的材料号；
excl:后面与破坏有关选项取与此项相同的值时，将被绕过，即不会考虑这种破坏；
mxpres(Stress-Pressure):最大静水压力，应取为正值，代表受压。所谓静水压力是三个主应力的平均值。本选项用于控制受压破坏；
mneps(Infin-3rd-principal strain):第三主应变极限值，应取为负值，代表受拉。本选项用于控制受拉破坏；
effeps(Infin-effective strain(v-m)):有效应变，V-M应变；
voleps(Infin-mean strain):体积应变极限值。正负都可，分别代表体积压缩和膨胀。在三维模型中voleps是mean strain的三倍；
numfip:在删除单元前失效积分点的数量，默认为1；
ncs:设置满足几条失效条件时删除单元，默认值为1.也即只要有一条失效条件满足就删除单元；
mnpres(Stress-Pressure):与mxpres相反，用以控制受拉破坏，一般取负值，代表受拉；
sigp1(Stress-1st-principal stress):第一主应力极限值，用以控制受压破坏；
sigvm(Stress-von mises stress):V-M应力限值；
mxeps(Infin-1st-principal strain):与mneps相反，用以控制受压破坏，一般取正值，代表受压；
epssh(Infin-max shear strain):最大剪应力限值；
sigth, impulse, failtm:与有效作用时间有关的失效判断相关参数。
以上是每一项的典型值单位制为SI，仅供定性参考。这些数据都是根据只有一个尺寸为1*1*1的单元的模型材料参数不断变化试验得到的，不妥之处还望指正。
*MAT_ADD_EROSION
$# mid excl mxpres mneps effeps voleps numfip ncs
1 1234.0000 1.0000E+8-1.0000E-3 1.0000E-3-1.0000E-3
$#mnpres sigp1 sigvm ...
alinheng 发表于 2011-4-3 13:41 http://forum.simwe.com/images/common/back.gif
静水以压力是压为正，拉为负，而应力则相反。
那么楼主的解释中MNEPS\ mxeps\sigp1是不是反了
mneps(Infin-3rd-principal strain):应取为负值，代表受压。本选项用于控制受压破坏；
mxeps(Infin-1st-principal strain):与mneps相反，用以控制受拉破坏，一般取正值，代表受拉；
sigp1(Stress-1st-principal stress):第一主应力极限值，用以控制受拉破坏
望大侠们解惑 6# zenis
那个贴应该针对的是老版本的DYNA手册，971手册中MAT_ADD_EROSION关键字CARD2中已更改为MNPRES（Minimum pressure at failure，亦即Pmin），并在CARD1中增加了MXPRES（Maximum pressure at failure，亦即Pmax）。当P>Pmax或P<Pmin时材料均会失效。如果定义岩石类材料受拉失效，将Pmin定义为某负值即可；如果定义岩石类材料受压失效，将Pmax定义为某正值即可。（P压为正，positive in compression） 麻烦楼主更正一下 以免给大家带来不便 楼主关于拉压正负号 搞反了 我认为lz对应力和应变的正负号的说法有问题。压应力为负，拉应力为正，应变也是这样的。建议楼主做个简单的压缩载荷或拉伸载荷例子，取出应力，看正负号。 本帖最后由 zpyc01 于 2011-4-7 23:10 编辑
1# alinheng
“mxpres(Stress-Pressure):最大静水压力，应取为正值，代表受压。所谓静水压力是三个主应力的平均值。本选项用于控制受压破坏”这句和另一个贴说的有出入。地址：http://forum.simwe.com/viewthread.php?tid=961333贴中说：“关于PFAIL关键字的说明：此关键字表示物体的静水压破坏，即各个方向受到压力时的破坏准则，其中压为正，拉为负，一般材料尤其是混凝土材料都是拉伸破坏，故此参数一般定义为负数，对于大小比较的是代数值的大小，因此当低于此准则即拉伸力超过允许数值，材料即宣告破坏（类似抗压强度）。当实际的静水压力（其实应该是拉力）小于此值（代数大小），材料即宣告破坏”。
不知道楼主取正值针对的是什么材料。对于岩石类材料，是应该取正值还是负值？ 感谢LZ热心的工作，能否再将模拟的参数调整过程说明清楚一些，比如针对什么样的材料，是否有试验数据的验证？:handshake MAT_ADD_EROSION 在lsdyna971R5手册中变动比较大。
而且R4之前的只支持单点积分的体单元。These failure models apply only to solid elements with one point integration in 2 and 3 dimensions。
但是In the R4 release of 971, this option applies to
nonlinear element formulations including the 2D continuum, 3D solid elements, 3D shell
elements, and the thick shell elements types 1 and 2.
而且在我计算中，不能用单点积分，只能用实体单元中的算法2，否则会一下子删除很多的单元。 "sigvm(Stress-von mises stress):V-M应力限值；"
我看有些帖子里说是抗拉强度,不知道对不对? 本帖最后由 zpyc01 于 2011-4-7 23:10 编辑
1# alinheng
“mxpres(Stress-Pressure):最大静水压力，应取为正值，代表受压。所谓静水压力是三个主应力的平均值。本选项用于控制受压破坏”这句和另一个贴说的有出入。地址：http://forum.simwe.com/viewthread.php?tid=961333贴中说：“关于PFAIL关键字的说明：此关键字表示物体的静水压破坏，即各个方向受到压力时的破坏准则，其中压为正，拉为负，一般材料尤其是混凝土材料都是拉伸破坏，故此参数一般定义为负数，对于大小比较的是代数值的大小，因此当低于此准则即拉伸力超过允许数值，材料即宣告破坏（类似抗压强度）。当实际的静水压力（其实应该是拉力）小于此值（代数大小），材料即宣告破坏”。
不知道楼主取正值针对的是什么材料。对于岩石类材料，是应该取正值还是负值？ 7# zpyc01
明白了,thanks 我还是有很多的疑问呀。 请问，在计算中怎样能够不使用单点积分，只用实体单元中的算法2呢？在哪里定义呢？ 10# chaofan710
*SECTION_SOLID
1 2 11# zpyc01
我刚试算了一下，发现出错了，它说
好像是不能用第二种算法，只能用第一种。那怎么办呢？有很多单元被删掉的问题怎么解决啊 本帖最后由 zpyc01 于 2011-4-12 16:34 编辑
12# chaofan710
你是不是加了MAT_ADD_EROSION，该关键字目前还只能用于单点积分的2D、3D单元，MAT_ADD_EROSION还不适用于采用全积分的单元 对于这个关键字终于有了一个全面的认识，谢谢！ 13# zpyc01
是的，就是添加了MAT-ADD-EROSION。请问，如果不添加这个，能不能模拟出混凝土局部破坏，就是某些单元被删除？像图中那样
15# chaofan710
可以采用自身带失效模式的本构模型，如*MAT_PLASTIC_KINEMATIC 首先感谢楼主的热心解释
楼主的解释中 epssh--Shear strain at failure ,应为失效剪应变吧？ 18# laoqiang117
我觉得好像是反了。 楼主有关于混凝土失效的参数吗？谢谢