MAT_ADD_EROSION 关键字参数解释(期待高手完善)

alinheng · 发表于 2011-4-3 13:41:47

本帖最后由 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的单元的模型材料参数不断变化试验得到的，不妥之处还望指正。

laoqiang117 · 发表于 2011-4-21 20:36:49

*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

静水以压力是压为正，拉为负，而应力则相反。
那么楼主的解释中MNEPS\ mxeps\sigp1是不是反了
mneps(Infin-3rd-principal strain):应取为负值，代表受压。本选项用于控制受压破坏；
mxeps(Infin-1st-principal strain):与mneps相反，用以控制受拉破坏，一般取正值，代表受拉；
sigp1(Stress-1st-principal stress):第一主应力极限值，用以控制受拉破坏
望大侠们解惑

zpyc01 · 发表于 2011-4-7 23:18:38

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）

marlboro_25214 · 发表于 2014-7-22 04:53:19

麻烦楼主更正一下以免给大家带来不便楼主关于拉压正负号搞反了

雕栏玉砌 · 发表于 2012-12-24 09:25:53

我认为lz对应力和应变的正负号的说法有问题。压应力为负，拉应力为正，应变也是这样的。建议楼主做个简单的压缩载荷或拉伸载荷例子，取出应力，看正负号。

zenis · 发表于 2011-4-7 20:33:33

本帖最后由 zpyc01 于 2011-4-7 23:10 编辑

1# alinheng
“mxpres(Stress-Pressure):最大静水压力，应取为正值，代表受压。所谓静水压力是三个主应力的平均值。本选项用于控制受压破坏”这句和另一个贴说的有出入。地址：http://forum.simwe.com/viewthread.php?tid=961333贴中说：“关于PFAIL关键字的说明：此关键字表示物体的静水压破坏，即各个方向受到压力时的破坏准则，其中压为正，拉为负，一般材料尤其是混凝土材料都是拉伸破坏，故此参数一般定义为负数，对于大小比较的是代数值的大小，因此当低于此准则即拉伸力超过允许数值，材料即宣告破坏（类似抗压强度）。当实际的静水压力（其实应该是拉力）小于此值（代数大小），材料即宣告破坏”。
不知道楼主取正值针对的是什么材料。对于岩石类材料，是应该取正值还是负值？

ywluck · 发表于 2011-4-3 16:45:09

感谢LZ热心的工作，能否再将模拟的参数调整过程说明清楚一些，比如针对什么样的材料，是否有试验数据的验证？

lojade · 发表于 2011-4-3 18:24:08

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，否则会一下子删除很多的单元。

zenis · 发表于 2011-4-4 10:01:23

"sigvm(Stress-von mises stress):V-M应力限值；"
我看有些帖子里说是抗拉强度,不知道对不对?

zenis · 发表于 2011-4-7 20:33:37

本帖最后由 zpyc01 于 2011-4-7 23:10 编辑

1# alinheng

“mxpres(Stress-Pressure):最大静水压力，应取为正值，代表受压。所谓静水压力是三个主应力的平均值。本选项用于控制受压破坏”这句和另一个贴说的有出入。地址：http://forum.simwe.com/viewthread.php?tid=961333贴中说：“关于PFAIL关键字的说明：此关键字表示物体的静水压破坏，即各个方向受到压力时的破坏准则，其中压为正，拉为负，一般材料尤其是混凝土材料都是拉伸破坏，故此参数一般定义为负数，对于大小比较的是代数值的大小，因此当低于此准则即拉伸力超过允许数值，材料即宣告破坏（类似抗压强度）。当实际的静水压力（其实应该是拉力）小于此值（代数大小），材料即宣告破坏”。
不知道楼主取正值针对的是什么材料。对于岩石类材料，是应该取正值还是负值？

zenis · 发表于 2011-4-8 09:42:47

7# zpyc01
明白了,thanks

369770836 · 发表于 2011-4-8 16:17:56

我还是有很多的疑问呀。

chaofan710 · 发表于 2011-4-12 15:23:03

请问，在计算中怎样能够不使用单点积分，只用实体单元中的算法2呢？在哪里定义呢？

zpyc01 · 发表于 2011-4-12 15:29:29

10# chaofan710
*SECTION_SOLID
1 2

chaofan710 · 发表于 2011-4-12 16:07:28

11# zpyc01
我刚试算了一下，发现出错了，它说

好像是不能用第二种算法，只能用第一种。那怎么办呢？有很多单元被删掉的问题怎么解决啊

zpyc01 · 发表于 2011-4-12 16:31:00

本帖最后由 zpyc01 于 2011-4-12 16:34 编辑

12# chaofan710
你是不是加了MAT_ADD_EROSION，该关键字目前还只能用于单点积分的2D、3D单元，MAT_ADD_EROSION还不适用于采用全积分的单元

hantianyi · 发表于 2011-4-12 20:29:46

对于这个关键字终于有了一个全面的认识，谢谢！

chaofan710 · 发表于 2011-4-14 14:22:22

13# zpyc01
是的，就是添加了MAT-ADD-EROSION。请问，如果不添加这个，能不能模拟出混凝土局部破坏，就是某些单元被删除？像图中那样

zpyc01 · 发表于 2011-4-14 15:05:46

15# chaofan710
可以采用自身带失效模式的本构模型，如*MAT_PLASTIC_KINEMATIC

laoqiang117 · 发表于 2011-4-19 09:40:25

首先感谢楼主的热心解释
楼主的解释中 epssh--Shear strain at failure ,应为失效剪应变吧？

alinheng · 发表于 2011-4-22 13:40:59

18# laoqiang117
我觉得好像是反了。

NWLNLENWLNL0379 · 发表于 2011-4-25 13:16:26

楼主有关于混凝土失效的参数吗？谢谢

账号		自动登录	找回密码
密码			注册

账号		自动登录	找回密码
密码			立即注册

[材料模型] MAT_ADD_EROSION 关键字参数解释(期待高手完善)

评分

相关帖子

点评

评分

本帖子中包含更多资源

本帖子中包含更多资源