流动规则（镍钛诺材料模型）
通过采用对数应变定义，应变和应力张量的偏量和体积分量及其关系可以通过分离形式正确表示。
首先，考虑用下式表示总的塑性和弹性应变向量：
ε(bar)p = εul ξs( n(bar) + α*m(bar))
ε(bar)e(bar) = ε(bar) - ε(bar)p
可以通过下式计算基尔霍夫应力向量：
τ(bar) = p m(bar) + t(bar)
p = K (θ - 3 α εul ξs)
t = 2 G (e(bar) - εul ξsn(bar))
在上述公式中：
εul | 代表最大材料塑性应变变形的标量参数 [EUL] |
ξ秒 | 作为测量塑性应变的介于 0 和 1 之间的参数 |
θ | 体积应变 = ε11 + ε22 + ε33 |
e(bar) | 偏应变向量 |
t(bar) | 偏应力向量 |
n(bar) | 偏应力的范数 = t(bar)/ (sqrt(2) σ(bar))。 |
m(bar) | 向量形式的单位矩阵：{1,1,1,0,0,0}T |
K 和 G | 体积和抗剪弹性模量： K = E / [3(1-2ν)]，G = E / [2(1+ν)] |
增量形式的线性流动规则可以相应表示为：
装载: Δξs = ( 1.0 - ξs) ΔF / ( F - Rf1)
卸载： Δξs = ξs ΔF / ( F - Rf2)
在定义非零 β 时所使用的指数流动规则为：
装载: Δξs = β1( 1.0 - ξs) ΔF / ( F - Rf1)2
卸载： Δξs = β2ξs ΔF / ( F - Rf2)2
一般会发现形状记忆合金对速率效应不敏感。因此，在上述公式中，“时间”代表一个假定变量，其长度不影响解。
上述所有公式都是针对张力装载-卸载，因为对于压缩装载-卸载条件可以使用类似的表达式（使用压缩属性参数）。
此处的增量求解算法在计算应力和本构公式时使用返回-映射步骤作为一个求解步骤。因此，求解由两部分组成。首先，计算初始状态；然后，如果初始状态违反流动准则，则进行调整，将应力返回至流动表面。
