From cee8307209c7b357b306acd92f8f2b3c9974c350 Mon Sep 17 00:00:00 2001 From: labourasse <labourassee@gmail.com> Date: Wed, 1 Feb 2023 09:17:38 +0100 Subject: [PATCH] symmetrize the Cauchy Tensor --- src/scheme/HyperelasticSolver.cpp | 2 ++ 1 file changed, 2 insertions(+) diff --git a/src/scheme/HyperelasticSolver.cpp b/src/scheme/HyperelasticSolver.cpp index c1cf4fdfb..962e133fc 100644 --- a/src/scheme/HyperelasticSolver.cpp +++ b/src/scheme/HyperelasticSolver.cpp @@ -464,6 +464,8 @@ class HyperelasticSolverHandler::HyperelasticSolver final : public HyperelasticS new_u[j] += dt_over_Mj * momentum_fluxes; new_E[j] += dt_over_Mj * energy_fluxes; new_CG[j] += dt_over_Mj * cauchy_green_fluxes; + new_CG[j] += transpose(new_CG[j]); + new_CG[j] *= 0.5; }); CellValue<const double> new_Vj = MeshDataManager::instance().getMeshData(*new_mesh).Vj(); -- GitLab