为什么我的hill48不是二次收敛呢

ETAYU · 发表于 2021-6-16 16:23:10

本帖最后由 ETAYU 于 2021-6-16 16:29 编辑

我写了一个hill48的子程序，运算结果应该没有问题，但是今天打开message file，然后发现好像我运算并不是二次收敛，而是一次的。我猜测可能是DDSDDE出了问题，但是反复检查公式又找不出问题，有点不知道该怎么办了。以下是进入塑性之后的某一个increment各个interation对应的最大残差以及平均力。
【Iteration】       【LARGEST RESIDUAL FORCE】       【TIME AVG. FORCE】       【Rα】
   1          -0.580                30.8       0.0188
   2          0.345                30.8       0.0112
   3          -0.183                30.9       0.0059
   4          0.105                30.9       0.0034
   5          -0.058                31.0       0.0019
   6          0.032                31.0       0.0010
   7          -1.80E-02       31.0       0.0006
而且每一次残差的正负号都变了，请问这个是不是能说明什么问题呢？
这是屈服准则

这是ddsdde的计算公式

其中

（发帖新人，这句话后面两张图片不知道为什么删除不了，请忽略）

ETAYU · 发表于 2021-6-16 20:53:26

希望有人能看到

ETAYU · 发表于 2021-6-21 15:18:37

再顶一次

abaqus109 · 发表于 2022-9-14 13:27:23

为什么不发一下for文件呢？也可能不是这几步的错误吧？

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

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

[子程序] 为什么我的hill48不是二次收敛呢

本帖子中包含更多资源