Fix mines

a-latyshev · Oct 13, 2024 · a38f5db · a38f5db
1 parent 83f4d54
commit a38f5db
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/doc/demo/demo_plasticity_mohr_coulomb.py b/doc/demo/demo_plasticity_mohr_coulomb.py
@@ -56,7 +56,7 @@
 # $$
 #     F(\boldsymbol{u}; \boldsymbol{v}) = \int\limits_\Omega
 #     \boldsymbol{\sigma}(\boldsymbol{u}) \cdot
-#     \boldsymbol{\varepsilon}(\boldsymbol{v}) \, \mathrm{d}\boldsymbol{x} +
+#     \boldsymbol{\varepsilon}(\boldsymbol{v}) \, \mathrm{d}\boldsymbol{x} -
 #     \int\limits_\Omega \boldsymbol{q} \cdot \boldsymbol{v} \, \mathrm{d}\boldsymbol{x} = \boldsymbol{0}, \quad
 #     \forall \boldsymbol{v} \in V,
 # $$

diff --git a/doc/demo/demo_plasticity_von_mises.py b/doc/demo/demo_plasticity_von_mises.py
@@ -444,7 +444,7 @@ def sigma_external(derivatives):
 # %%
 u = fem.Function(V, name="displacement")
 du = fem.Function(V, name="Newton_correction")
-external_operator_problem = LinearProblem(J_replaced, F_replaced, Du, bcs=bcs)
+external_operator_problem = LinearProblem(J_replaced, -F_replaced, Du, bcs=bcs)
 
 # %%
 # Defining a cell containing (Ri, 0) point, where we calculate a value of u
@@ -508,7 +508,7 @@ def sigma_external(derivatives):
     sigma_n.x.array[:] = sigma.ref_coefficient.x.array
 
     if len(points_on_process) > 0:
-        results[i + 1, :] = (-u.eval(points_on_process, cells)[0], loading.value / q_lim)
+        results[i + 1, :] = (u.eval(points_on_process, cells)[0], loading.value / q_lim)
 
 # %% [markdown]
 # ### Post-processing