我正在使用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)