Mercurial > repos > public > sbplib

changeset 1319:8b1110385ee2 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Bugfix interface forcing in ElasticCurveAnisotropic
author Martin Almquist <malmquist@stanford.edu>
date Tue, 28 Jul 2020 21:58:22 -0700
parents 6650b111742d
children cd0bfcaef0ba
files +scheme/Elastic2dCurvilinearAnisotropic.m
diffstat 1 files changed, 1 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
diff -r 6650b111742d -r 8b1110385ee2 +scheme/Elastic2dCurvilinearAnisotropic.m
--- a/+scheme/Elastic2dCurvilinearAnisotropic.m	Sun Jul 26 20:42:06 2020 -0700
+++ b/+scheme/Elastic2dCurvilinearAnisotropic.m	Tue Jul 28 21:58:22 2020 -0700
@@ -702,6 +702,7 @@
             % Continuity of tangential displacement
             closure = closure + 1/2*tau_t_u*H_gamma*et_u';
             penalty = penalty + 1/2*tau_t_u*H_gamma*et_v';
+            forcing_u_t = forcing_u_t - 1/2*tau_t_u*H_gamma;
             % ------------------------------------------------------------------
 
             % --- Continuity of tractions -----------------------------