changeset 730:9f28cf266f86 feature/poroelastic

Bugfix Elastic2dVariable interface
author Martin Almquist <malmquist@stanford.edu>
date Sat, 21 Apr 2018 16:04:46 -0700
parents aa8cf3851de8
children c3b460a7d7ac
files +scheme/Elastic2dVariable.m
diffstat 1 files changed, 1 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
--- a/+scheme/Elastic2dVariable.m	Fri Apr 20 16:56:49 2018 -0700
+++ b/+scheme/Elastic2dVariable.m	Sat Apr 21 16:04:46 2018 -0700
@@ -400,7 +400,7 @@
                 a2 = ( 16*(th1 + th2)*lambda*mu ) / (th1*th2*th3);
                 alpha_ii = a1 + sqrt(a2 + a1^2);
 
-                alpha_ij = 2/h11 + 1/(phi*h11);
+                alpha_ij = mu*(2/h11 + 1/(phi*h11));
             end
 
             [alpha_ii_u, alpha_ij_u] = computeAlpha(phi_u,h_u,h11_u,gamma_u,lambda_u,mu_u);