J-C模型子程序计算结果没有任何应力

qdx天风海涛 · 发表于 2014-6-16 19:57:05

我上传了一下我的模型和相关的子程序文件，我的子程序在运行结束后，在云图里面显示没有任何应力。后来我尝试在程序中间加入write语句输出相关应力的值，发现所有的应力值均为0（包括Mises应力），后来我把程序中关于塑性的语句注释掉，再次运行发现计算结果的云图里面会显示相关的应力云图，不知道这是为什么，为什么只要把后面塑性的语句加上，前面弹性部分计算出来的Mises应力也全部为0，请高手指点

qdx天风海涛 · 发表于 2014-6-16 19:58:34

忘了说了，这个子程序式是很常见的卢剑锋庄茁* 张帆写的JC模型

surenjia · 发表于 2014-6-16 20:05:05

我有个问题:aba自己的J-C不能满足您的要求么?

qdx天风海涛 · 发表于 2014-6-16 20:08:16

surenjia 发表于 2014-6-16 20:05
我有个问题:aba自己的J-C不能满足您的要求么?

abaqus自带的JC模型只能在abaqus/Explicit求解下应用，这个被别人修改过的，能够在abaqus/standard里面用，还有就是我需要写一个其他本构的子程序，想先把别人写出来的这个弄明白了，再写自己的会好一点儿，关键是现在别人的都运行不出来，更别提自己写了。

surenjia · 发表于 2014-6-16 20:11:54

哦,没有玩过子程序,帮不了你.

cartographer · 发表于 2014-6-16 21:27:27

程序状态变量没有更新

另外状态变量应该是13个
========================

加载时间和幅值函数估计有问题去找个帖子看看吧

加载时间2e-4
幅值函数时间都到1了

另外幅值是指比值不是加200的力幅值就填200

======
最好输出一下状态变量sdv

qdx天风海涛 · 发表于 2014-6-16 21:31:02

这次传上来的子程序能够运行了，一开始我觉得这句话是没有用的IF(NPROPS.GT.5.AND.PROPS(4).GT.0.0)THEN，就给删掉了，因为IF后面的判断语句是肯定会成立的。后来加上之后发现能够正确运行了，想请教一个有没有人懂这句话的意思，为什么加上这句话就能运行

qdx天风海涛 · 发表于 2014-6-16 21:40:31

cartographer 发表于 2014-6-16 21:27
程序状态变量没有更新

另外状态变量应该是13个

您好，对于这个幅值函数我是真的不太明白，想问一下左边的Time/Frequency里面是指的step里面的时间，所填的范围也应该和step里面的时间范围相对应；右边的Amplitude是指的一个比例，应该填0-1之间的数，而这个时间施加的力的值就应该是这个比例值和先前加载时设置的Magnitude的乘积，不知道我这样理解对吗？

qdx天风海涛 · 发表于 2014-6-17 09:28:20

cartographer 发表于 2014-6-16 21:27
程序状态变量没有更新

另外状态变量应该是13个

你说的这两个问题我都修改了，最后还是不行，计算结果还是没有任何变化，和先前的结果一样

dashengyue · 发表于 2014-6-17 17:21:32

删除IF(NPROPS.GT.5.AND.PROPS(4).GT.0.0)THEN这句时，有没有相应的删除END IF？

qdx天风海涛 · 发表于 2014-6-19 09:14:39

dashengyue 发表于 2014-6-17 17:21
删除IF(NPROPS.GT.5.AND.PROPS(4).GT.0.0)THEN这句时，有没有相应的删除END IF？

删除了，我就是发现那句话没用，才删掉的

dashengyue · 发表于 2014-6-19 11:13:31

qdx天风海涛发表于 2014-6-19 09:14
删除了，我就是发现那句话没用，才删掉的

我删除了这句话，运行结果无变化。

qdx天风海涛 · 发表于 2014-6-20 09:19:27

dashengyue 发表于 2014-6-19 11:13
我删除了这句话，运行结果无变化。

结果是没变化，但是我找的这个例子，应该是比较经典的了，那个人的论文也发表出来了，我一个字不落的复制过来，发现就是不能运行

spencer0417 · 发表于 2014-6-20 09:56:13

你的A B m n 的数值都是在哪定义的程序里不写么？

dashengyue · 发表于 2014-6-20 17:17:43

qdx天风海涛发表于 2014-6-20 09:19
结果是没变化，但是我找的这个例子，应该是比较经典的了，那个人的论文也发表出来了，我一个字不落的复制 ...

我就用你的INP，计算的。。。

qdx天风海涛 · 发表于 2014-6-21 09:07:50

spencer0417 发表于 2014-6-20 09:56
你的A B m n 的数值都是在哪定义的程序里不写么？

那几个值都是在User materials里面输进去的

spencer0417 · 发表于 2014-6-23 10:09:24

qdx天风海涛发表于 2014-6-21 09:07
那几个值都是在User materials里面输进去的

然后再加上子程序就能运行JC模型了？

spencer0417 · 发表于 2014-6-23 10:54:36

qdx天风海涛发表于 2014-6-21 09:07
那几个值都是在User materials里面输进去的

就是这个子程序，只要把庄茁的程序原封不动调用就可以了么？一点不用改么？

spencer0417 · 发表于 2014-6-23 14:44:32

spencer0417 发表于 2014-6-23 10:09
然后再加上子程序就能运行JC模型了？

你好可以加个qq讨论一下么 381438467 我也在研究jc模型也有些问题需要讨教

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

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

[子程序] J-C模型子程序计算结果没有任何应力

本帖子中包含更多资源

评分

本帖子中包含更多资源

点评