Mercurial > repos > public > sbplib

comparison +scheme/Elastic2dCurvilinearAnisotropic.m @ 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 34997aced843
children b40359c9faed
comparison
equal deleted inserted replaced
1318:6650b111742d 1319:8b1110385ee2
700 forcing_u_n = forcing_u_n - 1/2*tau_n_u*H_gamma; 700 forcing_u_n = forcing_u_n - 1/2*tau_n_u*H_gamma;
701 701
702 % Continuity of tangential displacement 702 % Continuity of tangential displacement
703 closure = closure + 1/2*tau_t_u*H_gamma*et_u'; 703 closure = closure + 1/2*tau_t_u*H_gamma*et_u';
704 penalty = penalty + 1/2*tau_t_u*H_gamma*et_v'; 704 penalty = penalty + 1/2*tau_t_u*H_gamma*et_v';
705 forcing_u_t = forcing_u_t - 1/2*tau_t_u*H_gamma;
705 % ------------------------------------------------------------------ 706 % ------------------------------------------------------------------
706 707
707 % --- Continuity of tractions ----------------------------- 708 % --- Continuity of tractions -----------------------------
708 709
709 % Continuity of normal traction 710 % Continuity of normal traction