718933105 发表于 2016-7-5 10:32:31

Mat-brittle-damage

整体式模型
Mat-brittle-damage
密度,杨氏模量,泊松比,拉伸极限,剪切极限,剪力传递系数,压缩屈服应力
配筋率,钢筋杨氏模量,钢筋屈服应力,钢筋硬化模量,钢筋失效应变
2500,30Gpa,0.2,2.5Mpa,14.5Mpa,0.03,28.97Mpa
0.8%,200Gpa,335Mpa,1.18Gpa,0.1

*MAT_ADD_EROSION
CARD1
MID(材料号)      EXCL(删除号)MXPRES(最大压力)MNPEPS(最小主应变)
CARD2
PFAIL(最小压力)SIGPL(最大主应力)SIGVM (平均最大应力) EPSPL(最大主应变)EPSSH(剪应变)SIGTH(开始应力)IMPULSE(应力脉冲)FAILTM(失效时间)

以上是Mat-brittle-damage本构模型,我用此本构模型模拟框架结构爆破倒塌,并且设置了,*MAT_ADD_EROSION中最大主应力 SIGPL(30Mpa)或者最大主应变EPSPL(0.015)失效准则
可是倒塌结果和现实相差太远,是不是参数取得不合理?而且,框架中有的单元最大主应力没有到达30Mpa就失效,有的超过还没失效,请问又是怎么回事?有知道的大神请解答。

empty4585 发表于 2016-7-6 09:42:53

718933105 发表于 2016-7-6 08:46
我再K文件*mat-erosion
*MAT_ADD_EROSION(两个卡片,总共12个参数)
CARD1


哥们:
1, sigmaY、sigma1、sigmaVM要分清楚,你定义的是sigma1,程序就只判断sigma1,跟那两个无关。
2,你是看的d3plot下的history吧,没有单独采用ascii的方式来记录你要关心的这几个单元,所以关键的问题是:采样频率不够!明白了不,所以你绘出的曲线过于离散,你看到的曲线最大值到后面删除归零之间的部分没有描绘出来。采用d3thdt或者ascii来重新捕捉下,会达到满意的效果。
3,你以后叙述问题能不能认真,谨慎一些,一会mpa,一会kpa,一会stress-y的,叙述清楚点,要不谁有那么多时间一点点看你写得这一堆,去判断你要说什么。

采用频率问题,明白了吧?好了,就这样,不用私聊我的。

wangyu341124 发表于 2016-7-5 16:49:53

你好 我也是在做这方面的模拟 ,有好多不懂不知道能否交流下 ,我的QQ:1527578674,

718933105 发表于 2016-7-6 08:40:06

麻烦你看一下我之前发的帖子,有图片和曲线图,就是没达到就失效了啊,而且有的超过了却不失效:(

718933105 发表于 2016-7-6 08:46:22

我再K文件*mat-erosion
*MAT_ADD_EROSION(两个卡片,总共12个参数)
CARD1
MID      EXCLMXPRESMNPEPS
CARD2
PFAILSIGPL(最大主应力)SIGVMEPSPL (最大主应变) EPSSHSIGTHIMPULSEFAILTM
部分中作出了如下设置了:


*MAT_PLASTIC_KINEMATIC
1,0.250E+04,0.300E+11,0.200000,2.0E+7,0.0,0.0

*MAT_BRITTLE_DAMAGE(96号整体式砼材料)
         2,2500,0.300E+11,0.200000,6.00E+06,1.40E+7,0.5E+06,0.03
0.724E+06,0.050,0.200E+12,3.35E+08,1.1E+10,0.1,30.0E+06
*MAT_BRITTLE_DAMAGE
         3,2500,0.300E+11,0.200000,6.00E+06,1.40E+7,0.5E+06,0.03
0.724E+06,0.050,0.200E+12,3.35E+08,1.1E+10,0.025,30.0E+06
$
*mat_add_erosion
2,567
567,3.0E+7(最大主应力),567,0.015(最大主应变),567,567,567,567
$
*mat_add_erosion
3,567
567,567,567,567,567,567,567,0.5(失效时间)


我之前也分别设置了最大主应力和最大主应变的值,材料失效是的值也对应不上,请问怎么回事?
请大神指教

718933105 发表于 2016-7-7 09:10:03

empty4585 发表于 2016-7-6 09:42
哥们:
1, sigmaY、sigma1、sigmaVM要分清楚,你定义的是sigma1,程序就只判断sigma1,跟那两个无关。
2 ...

绝对的大神,谢谢谢谢,是这样的,本来应该是每隔0.01s输出一个步,但是这样计算时间太长了,为了节约时间,我按每0.1s输出一个步,才导致这样的。能不能加一下QQ好友啥的,我还有其他少数问题没解决,可能只需要你几句话,就能醍醐灌顶,却需要我自己摸索十天半个月的,甚至更长时间,我的QQ:718933105

我吃大橙子 发表于 2024-5-9 20:41:59

您好!请问下您在MAT96号材料定义的断裂韧度取值0.5E+06的取值单位是什么,数据有没有相关资料的支撑,感谢您的帮助

718933105 发表于 6 天前

我吃大橙子 发表于 2024-5-9 20:41
您好!请问下您在MAT96号材料定义的断裂韧度取值0.5E+06的取值单位是什么,数据有没有相关资料的支撑,感谢 ...

推荐一本书。辛春亮的,里面就有参数的来源和单位
《有限元分析常用材料参数手册》
页: [1]
查看完整版本: Mat-brittle-damage