thermal stress in openfoam

[大喵]

最近在学习solidDisplacementFoam时遇到了一个问题，在计算热应力是，openfoam是这样写的：

       if (thermalStress)
        {
            const volScalarField& T = Tptr();
            sigma = sigma - I*(rho*threeKalpha*T);
        }

可弹性力学中，应该是用温度变化计算热应力，即：

 sigma = sigma - I*(rho*threeKalpha*(T-T0));

其中T0是初始的参考温度。请问有人了解这个问题吗？请不吝赐教

[lushifer]

你好！可以参考foam-extend版本，相对而言，extend版本的固体力学求解器选择更多。我看了下，extend版本的考虑温度影响的弹性求解器是用的温度差更新应力的。