void unsIncrTotalLagrangianStress::updateTotalFields()
{
Info << "Update total fields" << endl;
D_ += DD_;
pointD_ += pointDD_;
gradD_ += gradDD_;
gradDf_ += gradDDf_;
sigma_ += DSigma_;
sigmaf_ += DSigmaf_;
epsilon_ += DEpsilon_;
epsilonf_ += DEpsilonf_;
U_ = fvc::ddt(D_);
if (rheology_.plasticityActive())
{
epsilonP_ += rheology_.DEpsilonP();
// epsilonPf_ += rheology_.DEpsilonPf();
rheology_.updateYieldStress();
}
}
这个文件夹下面也有这个updatetoalfields?