发一个粘弹损伤子程序，还有个问题要请教大家

hanfei_5759 · 发表于 2010-4-2 22:54:34

这是我编写的UMAT用户子程序，可以用来进行三维推进剂本构的计算，可以用来计算空穴和填充粒子的演化规律。拿来给大家分享，共同学习一下，子程序为lefmpp3d.for，inp文件为void.inp。
其中应变是用应变梯度张量来获取的，但是将子程序改成二维的lefmpp2d.for用来对一般的平板计算还可以，但是用来计算含裂纹的紧凑拉伸计算时就会出现不收敛，调了半天还是搞不定，inp文件名为ct.inp。不知道是不是基础知识太差还是别的什么原因？请大家帮忙看看，指导一下，不胜感激！

xxy_1986 · 发表于 2010-4-3 11:26:32

刚开始接触损伤，帮你顶一下~

warmwormdk · 发表于 2010-4-3 23:58:39

大概看了下，有些疑问
一般基于超弹性的本构模型都是用能量来表示的，而在楼主的umat中怎么没有发现能量的计算呢？
空穴率国内应该还不能测吧，楼主的式子正确吗？
看到你的形变张量进行了分解，但是没有看明白是采用的什么应变……

magicalarm · 发表于 2010-4-8 09:34:27

楼主的文件我都大致看了一下，有几个建议你可以参考一下：
1.裂纹处存在很大的应力集中，在很多时候都是会很难收敛的。裂尖处应进一步的细化网格。
2.在CT.inp中，你只在左上角的一点加了2方向的位移，也很容易产生应力集中。最好是采用耦合约束（coupling constraint）的方法添加载荷，具体方法见石亦平《ABAQUS有限元分析实例详解》。
3.楼主能否保证二维子程序的正确？是否可以先在无裂纹的平板中进行验算，再在CT试样中进行计算。

希望能给楼主一些帮助。共同探讨。

rongzj0830 · 发表于 2012-11-23 10:34:33

您好，子程序是用fortune语言编的吗

suyin · 发表于 2020-10-16 16:04:06

11111111111111

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

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

[断裂与失效] 发一个粘弹损伤子程序，还有个问题要请教大家

本帖子中包含更多资源

相关帖子