From ee09746c7e23a8cc34c81bfc5cfc794893c7ab9e Mon Sep 17 00:00:00 2001
From: labourasse <labourassee@gmail.com>
Date: Thu, 2 Feb 2023 08:39:26 +0100
Subject: [PATCH] hydroable version
---
src/scheme/HyperelasticSolver.cpp | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/scheme/HyperelasticSolver.cpp b/src/scheme/HyperelasticSolver.cpp
index b54a221a1..9a94c787f 100644
--- a/src/scheme/HyperelasticSolver.cpp
+++ b/src/scheme/HyperelasticSolver.cpp
@@ -235,7 +235,7 @@ class HyperelasticSolverHandler::HyperelasticSolver final : public HyperelasticS
for (size_t j = 0; j < node_to_cell.size(); ++j) {
const CellId J = node_to_cell[j];
const unsigned int R = node_local_number_in_its_cells[j];
- br += Ajr(J, R) * u[J] + sigma[J] * Cjr(J, R);
+ br += Ajr(J, R) * u[J] - sigma[J] * Cjr(J, R);
}
b[r] = br;
@@ -361,7 +361,7 @@ class HyperelasticSolverHandler::HyperelasticSolver final : public HyperelasticS
const auto& cell_nodes = cell_to_node_matrix[j];
for (size_t r = 0; r < cell_nodes.size(); ++r) {
- F(j, r) = Ajr(j, r) * (u[j] - ur[cell_nodes[r]]) + sigma[j] * Cjr(j, r);
+ F(j, r) = -Ajr(j, r) * (u[j] - ur[cell_nodes[r]]) + sigma[j] * Cjr(j, r);
}
});
--
GitLab