Mercurial > repos > public > sbplib
diff +scheme/elasticShearVariable.m @ 677:eeaf9a00e304 feature/poroelastic
Correct factor of 2 in Dirichlet penalty parameter.
| author | Martin Almquist <malmquist@stanford.edu> |
|---|---|
| date | Fri, 26 Jan 2018 15:17:49 -0800 |
| parents | 9e1d2351f539 |
| children |
line wrap: on
line diff
--- a/+scheme/elasticShearVariable.m Thu Jan 18 14:36:59 2018 -0800 +++ b/+scheme/elasticShearVariable.m Fri Jan 26 15:17:49 2018 -0800 @@ -237,7 +237,7 @@ for i = 1:obj.dim H_gamma_i = obj.H_boundary{i}; sigma_ij = tuning*delta(i,j)*2/(gamma*h(j)) +... - tuning*delta_b(i,j)*(2/(H11*h(j)) + 2/(H11*h(j)*phi)); + tuning*delta_b(i,j)*(2/(H11*h(j)) + 1/(H11*h(j)*phi)); ci = E{i}*inv(RHO)*nj*Hi*... ( (e{j}*H_gamma*e{j}'*A*e{j}*d{j}')'*E{i}' + ...