如何区分本构方程相对于fortran中的变量?

时间:2018-12-02 23:02:55

标签: fortran derivative abaqus differentiation

我正在使用fortran编写子例程来实现ABAQUS中TypeError: '>' not supported between instances of 'str' and 'float'给定的本构方程。通过输入上次迭代的应变,T_star和edot_log的值,我得到每个迭代的 yield 值。但是,我的代码也涉及三个区分

1)。关于应变

屈服(km)

2)。关于应变率

产量(公里)

3)。关于 T

产量(km)

目前,我正在手工找出导数,然后写下等式:

(A + B*(strain**xn))*(one - T_star**xm)* (one + C*edot_log)

这很好用,但是,有一种方法可以直接从收益方程 dyieldDeqps(km,1) = (B*xn*(strain**(xn-one)))* 1 (one - T_star**xm) * (one + C*edot_log) dyieldDeqps(km,2) = (A + B*(strain**xn))* 1 (one - T_star**xm) * (C/strainrate) dyieldDtemp(km) = (-xm)*(one/(T - T_trans))* 1 (T_star**(xm))*(A + B*(strain**xn)) 2 *(one + C*edot_log) 中找到导数。

任何帮助或建议都是最有帮助的。我不知道该怎么做?

Fortran子例程的相关部分:

(A + B*(strain**xn))*(one - T_star**xm)* (one + C*edot_log)

0 个答案:

没有答案