CFD中文网

@大喵 请问问题解决了吗？

在代码里加入

volScalarField rho ( IOobject ( "rho", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh );

然后直接用就行

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

最近打算重写一个 :mihu:

我也没看懂，楼主写得太复杂了，别人没功夫仔细看你的公式。我建议你自己用matlab算好每个mesh上的数值，然后写成csv格式用swak4Foam导入。