线性粘弹性本构模型（Prony形式）的二次开发问题

dbhj2 · 发表于 2012-9-5 22:52:07

本帖最后由 dbhj2 于 2012-9-8 17:57 编辑

问题描述：针对含损伤粘弹性本构模型 UMAT 二次开发，本人给出了算法并编制了相应的UMAT程序，尽管程序可以运行并得到了一些貌似可行的结果，但在进一步分析时发现了一些问题，于是将问题进一步退化到线性粘弹性本构问题，但问题仍然存在（后面附上了相关问题详细描述、算法及UMAT程序文件）。
所采用的线性粘弹性本构方程如下

其中

本人对程序存在问题的原因进行了一些初步检查分析，但目前仍未得到有效解。尽管问题简单却长时间苦于思索无果，寄希望于诸位同行高手帮忙指正，并恳请给出您宝贵的建议与意见，请各位赐教，万分感激！
我的邮箱dbhj0@sina.com

dbhj2 · 发表于 2012-9-5 22:56:27

本帖最后由 dbhj2 于 2012-9-12 09:50 编辑

下面是相关问题的详细描述及其对应的UMAT程序文件等附件

dbhj2 · 发表于 2012-9-6 13:36:30

本帖最后由 dbhj2 于 2012-9-8 17:50 编辑

各位同行，班主啊，期待给出您宝贵意见！

dbhj2 · 发表于 2012-9-9 12:04:32

本帖最后由 dbhj2 于 2012-9-9 12:43 编辑

尝试了几种算法，结果都基本一样的，后来发现程序有个错误:蠕变柔量系数应为正。似乎在于aba本身处理粘弹性问题：载荷处理，以及采用准静态法计算粘弹性时老提醒将采用弹性属性。

dbhj2 · 发表于 2012-9-12 09:48:43

本帖最后由 dbhj2 于 2012-9-12 09:48 编辑

对验证UMAT程序时，也做了采用ABAQUS中本身的粘弹性模型的分析设置，参照帮助文件进行了归一化处理，不过得到的结果好像却是弹性解，（与Marc的计算结果有出入），不知问题在哪里，怎么弄？好困惑啊，那位高手指点一下！

dbhj2 · 发表于 2012-9-12 09:51:17

dbhj2 发表于 2012-9-9 12:04
尝试了几种算法，结果都基本一样的，后来发现程序有个错误:蠕变柔量系数应为正。似乎在于aba本身处理粘弹性 ...

UMAT子程序已经更新

zcs220 · 发表于 2012-9-12 11:41:49

帮你顶下，UMAT不是我强项，下来看看。

dbhj2 · 发表于 2012-10-13 23:20:25

本帖最后由 dbhj2 于 2012-10-13 23:21 编辑

后来发现，上面的程序估计还没错，不过加压长时间计算确实有问题，偶尔再次读到UMAT帮助时发现，在解决近似不可压缩与不可压缩问题上，UMAT并不能很好解决，但是说可以采用UHYPER 或可以引入体积模量Kt进行解决，遗憾的是具体的也么有细说，估计这就是导致求解上述问题的UMAT的不足吧，近似不可压缩问题应该是UMAT的短板吧！

peterwei · 发表于 2015-10-19 05:14:10

楼主问题解决了吗？

轩落子 · 发表于 2016-1-27 11:19:21

加油，楼主

轩落子 · 发表于 2016-1-27 21:25:13

楼主编好了吗

danielzlf · 发表于 2017-8-15 13:33:09

楼主请问问题解决了吗？求回复，谢谢

danielzlf · 发表于 2017-8-16 16:44:24

请问版主可否联系，有些问题想请叫一下，谢谢QQ270296993

李1995 · 发表于 2018-11-1 09:36:35

值得学习111111111

好望矫 · 发表于 2019-5-29 20:33:46

值得学习，谢谢楼主

pcstudying · 发表于 2022-3-25 13:15:00

值得学习，谢谢楼主

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

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

[子程序] 线性粘弹性本构模型（Prony形式）的二次开发问题

本帖子中包含更多资源

本帖子中包含更多资源