自编材料subroutine没应力应变问题的解决方法

huihui123 · 发表于 2011-11-6 11:45:25

本帖最后由 huihui123 于 2011-11-9 17:39 编辑

白金泽和赵海鸥老师书上的例子是在ls970或之前的版本上编译成功的，用ls971的朋友们只需把subroutine umat41 (cm,eps,sig,hisv,dt,capa,etype,time,crv)改成 subroutine umat41(cm,eps,sig,epsp,hisv,dt,capa,etype,time,temper,failel,crv)就可以了！
补充：我看了一些朋友的子程序，他们的问题出在：
1.if (etype.eq.'brick')改成if (etype.eq.'solid')，或者把这个if（etype.eq.'...') endif语句去掉。
2.如果有考虑应变率，除以时间增量dt的时候要变成除以（dt+1E-8），避免除0，1E-8可以根据自己精度的需要而定。

suplukn · 发表于 2011-11-6 15:00:46

本帖最后由 suplukn 于 2011-11-6 15:25 编辑

:lol太感谢了，请问大牛，如何知道的？我赶快试试

suplukn · 发表于 2011-11-7 09:08:42

反映一下，修改过后，还是没有加速度和力云图

candy_xue · 发表于 2011-12-21 20:09:10

我也照着改了也不行啊。请问楼主还有什么别的方法吗？谢谢~~

飞云走壁 · 发表于 2012-3-30 10:38:06

期待更多解决措施的出现..

preking · 发表于 2013-2-1 15:06:21

正迷惑为什么不行呢

491742118 · 发表于 2014-2-25 15:59:25

楼主你好，我是个用户子程序新手，用了白金泽老师书里的例子和相应的f文件，其中第一个例子成功了，可第二个例子就是各向同性强化双线性塑性材料的能生成求解器，但是一求解k文件就产生负体积不知道怎么回事，求助！

491742118 · 发表于 2014-2-25 16:02:00

而且我看ls971版本的41号是这种结构啊：
subroutine umat41(cm,eps,sig,epsp,hsv,dt1,capa,etype,tt,temper,failel,crv)和您上面的有出入啊，应该是哪个呢？

liuxianlong · 发表于 2014-3-11 17:02:16

正在试，不知道群里面的各位怎么这么神

liuxianlong · 发表于 2014-3-13 11:18:04

大家参照此贴http://forum.simwe.com/forum.php ... D156%26typeid%3D156

lgh790723 · 发表于 2014-9-22 07:34:28

产生负体积一般是由于材料参数不合适，仔细检查输入参数的单位是否统一，有可能解决负体积问题

一尘不变 · 发表于 2018-6-15 15:50:11

491742118 发表于 2014-2-25 15:59
楼主你好，我是个用户子程序新手，用了白金泽老师书里的例子和相应的f文件，其中第一个例子成功了，可第二 ...

把关键字中NHV的1改为2就好了

ligenforever · 发表于 2018-7-2 17:04:54

suplukn 发表于 2011-11-7 09:08
反映一下，修改过后，还是没有加速度和力云图

想问修改后,答主的问题解决了吗?我现在也是没有应力的问题

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

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

[二次开发及其理论] 自编材料subroutine没应力应变问题的解决方法

评分