混凝土损伤模拟在abaqus中的二次开发

敦诚

以下是本人自己利用umat子程序在abaqus下的二次开发计算混凝土损伤的算例，子程序的源代码因为涉及到版权问题，没法给出，但是其理论来源在文件中给出了。

shujiqqing

万丈高楼平地起！
想做损伤方面的研究，先试着做了个弹性损伤的，但是结果有点问题。即：
1、如果是只有一个单元可以收敛，结果还可以
2、单元个数很多（实际中肯定是这样的），如果材料特性都一样，可以收敛，但是如果其中有一个单元的材料特性和其他不一样（见附件中的单元1000），就收敛不了。
真诚的向楼主及其他大牛请教！

下面说下一我做的算例的思路：
单轴压缩试验
弹性损伤，采用应变控制损伤的发展（为了简单起见，设定当最大压应变超过某一固定值后发生损伤，且固定了损伤D=0.9)。

感觉是很简单的东西，但是只有从简单的开始，才能加深对问题的理解！

lizhg06

ding:}

bearfl333

thank you!

shujiqqing

感谢楼主！

shujiqqing

shujiqqing

input及umat子程序

流星晴空

问问楼主损伤模型的参数具体是什么意思

shujiqqing

[local]1[/local]

xiaobai

这个问题我也很纳闷，希望大家讨论讨论！

过眼烟云

楼主的umat没发完吧，塑性部分在哪里？希望能交流

shujiqqing

原帖由 xiaobai 于 2008-6-9 15:49 发表
这个问题我也很纳闷，希望大家讨论讨论！

可以一起讨论讨论
楼主已经答应帮我们看一看了，这个问题肯定就能解决了。
不过楼主最近身体不好，在家休养呢，祝愿好人早日康复！

这个问题解决后，我会将解决方法再发上来的，大家共同进步！
（另不好意思，由于网络不好，那个附件发了3次)

敦诚

您好：
我看了您的inp，一共是11的分析步，在地11个分析步计算停止收敛，原因就是损伤值太大，一下就是0.9.这样的话单元的平衡计算一定过不去，我试着将您的损伤值调小为0.01.计算还是可以进行的，您可以试着将损伤值逐渐平缓地累计。一直到0.9，这样才可以在每一次迭代的过程中计算通过。
附件是我修改的inp，我在里面加了一个后处理显示的子程序GETVRM，这样就可以在后处理结果中检查损伤值的大小了。您的for中的第9个变量是damage，GETVRM中输出了15个自定义变量，够用了。

敦诚

http://www.simwe.com/forum/viewt ... d=319285#pid1355517
其实我的模型中，塑性与损伤是耦合在一起的，如果从金属塑性的情况考虑应力应变曲线的话，它是不会出现软化段的，一般软化段的解释就是认为金属材料出现了颈缩现象。
但是我模拟的是混凝土材料，软化段是由于损伤导致有效介面积较小，而后承载力下降。这是我全文的思路。

敦诚

关于子程序的流程图还有介绍在这里传一下。附件是一个ppt文件。

wwgi555

下载了各位大牛的资料，要好好学习，谢谢！

zhgy1003

顶一下
好东西啊

capicious

下载了,好好学习,3Q

lantianren

谢谢！

zltju

谢谢各位大牛了！