umat中if语句前后输出同一变量值不同

akakac · 发表于 2018-5-2 15:19:52

本人采用庄茁老师书后面JC本构umat,算了一个例子，结果不收敛，检查发现是在应力达到屈服应力附近时收敛困难。检查umat发现，在判断屈服的if语句处出现问题，在调入子程序后（其中包含有屈服应力），我在if语句前后各输出了一次屈服应力，if语句之前还比较正常，但if语句之后就有很大问题，某些单元节点应力达到了屈服应力的三倍（明显错误），一个单元内的高斯积分点有的正常，有的为零（明显错误），求大佬们指教。

akakac · 发表于 2018-5-2 15:25:18

自己顶一下，求大佬指教。

akakac · 发表于 2018-5-2 17:07:54

来个大佬解释一下，灰常感谢

akakac · 发表于 2018-5-3 10:02:36

自己再顶一下，来个大佬解释一下，或者哪位斑竹解释一下，非常感谢

akakac · 发表于 2018-5-3 10:04:22

我现在基本上发现问题出在哪儿，就是高斯积分点的值在if语句前后算到不一样，然后导致单元节点值差的很远。

akakac · 发表于 2018-5-3 10:34:56

貌似想明白了，每个积分点处进行屈服判断，超过屈服极限的积分点在if语句后才会输出新的屈服应力，未达到屈服应力的积分点就不输出了，所以为零

akakac · 发表于 2018-5-4 16:50:53

终于发现问题，是庄茁老师书上umat有一行打错了，写书的人不负责任真是坑死人

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

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

[二次开发] umat中if语句前后输出同一变量值不同