使用Umat和不使用Umat结果不一样？

zhangzheng_1979 · 发表于 2009-4-20 14:41:11

我最近用Umat做shell的分析，先通过一个最简单的拉伸来验证Umat子程序是否正确，刚开始就发现计算出来的应力与不用Umat子程序稍有差距，但误差还算可以接受。但模型复杂以后，怎么算都不对。后来重新对最简单的算例进行验证，通过比较刚度矩阵发现，主对角元相差较大，其他位置元素基本相同。不知道是什么原因？我将inp文件、umat子程序和输出的刚阵文件放在附件中，请大家鉴定。问题主要包括：1）我通过比较刚度矩阵K，发现使用与不使用UMAT主对角元有一定的差别，而其他项则基本完全相同，UMAT使用与否对计算结果是否有一定的影响引起误差？2）我编写的UMAT子程序过程是否正确？是否缺少或者漏掉某些东西？3）当计算扩展至非线性时，引入的误差是否比线性的更大？因为我们现在做的主要是涉及到几何非线性。谢谢大家。
我的文件名中注明：umat的说明是采用umat子程序的，nonumat的说明没有采用子程序。谢谢！

zhangzheng_1979 · 发表于 2009-4-20 19:39:08

请版主和大侠们指教！谢谢

zhangzheng_1979 · 发表于 2009-4-20 22:19:21

自己顶一下！

okyle · 发表于 2009-4-21 08:40:01

在线性区就有误差吗？

zhangzheng_1979 · 发表于 2009-4-21 11:26:24

是的，还没有引入非线性。

zhangzheng_1979 · 发表于 2009-4-21 20:28:47

高手救我！

zhangzheng_1979 · 发表于 2009-4-22 16:14:28

这个问题大家都没遇到过吗？用UMat之前都不验证一下？

zhangzheng_1979 · 发表于 2009-4-22 21:37:43

再次请教版主和大侠！

shawn2008 · 发表于 2009-4-22 23:00:18

8# zhangzheng_1979

平面应力状态下，你的ddsdde定义不对。

zhangzheng_1979 · 发表于 2009-4-23 11:44:14

shawn2008，谢谢您的回复。我看了您发的页面里的内容，其实和我的是一样的，只是写法不同。我直接用弹性模量和泊松比，你用的是lame常数，但计算出来的结果和我的Umat是一致的，但与不用Umat还是不一样。

zhangzheng_1979 · 发表于 2009-4-23 16:56:41

在线等待回答

zhangzheng_1979 · 发表于 2009-4-23 20:02:31

anybody can help me?!!

shawn2008 · 发表于 2009-4-24 00:11:32

10# zhangzheng_1979
我没仔细看，不好意思，这两个是一样的表达。

zhangzheng_1979 · 发表于 2009-4-24 09:55:45

呵呵，感谢您的关注。希望有高人出手！
13# shawn2008

zhangzheng_1979 · 发表于 2009-4-24 16:45:55

郁闷了好几天了，高手快出来吧

zhangzheng_1979 · 发表于 2009-4-24 20:03:14

自己顶起来！

zhangzheng_1979 · 发表于 2009-4-25 11:19:29

本帖最后由 zhangzheng_1979 于 2009-4-25 11:20 编辑

UMAT和不使用UMAT计算结果是否本身就有偏差？这个问题我现在还没得到权威的答复。昨天搜到一篇硕士论文其中对其进行了讨论，现在贴出来，请高手们辨析：里面对采用和不采用UMAT计算的结果进行了对比，发现二者之间并不相同，为什么会相差这么大呢？

zhangzheng_1979 · 发表于 2009-4-27 15:07:25

Still no answer
自己顶！

shawn2008 · 发表于 2009-4-27 18:21:47

你确定你的结果是在相同的时间比较的？

zhangzheng_1979 · 发表于 2009-4-27 21:59:04

我的inp文件两者不同之处就是将这行：
*Elastic
替换成：
*User Material, constants=2

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

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

[子程序] 使用Umat和不使用Umat结果不一样？

本帖子中包含更多资源

相关帖子

本帖子中包含更多资源

本帖子中包含更多资源