Mercurial > repos > public > sbplib
changeset 919:e30aaa4a3e09 feature/poroelastic
Bugfix in Elastic2dVariable.get_boundary_ops
| author | Martin Almquist <malmquist@stanford.edu> |
|---|---|
| date | Wed, 28 Nov 2018 14:04:31 -0800 |
| parents | 7d4f57725192 |
| children | 72cd29107a9a |
| files | +scheme/Elastic2dVariable.m |
| diffstat | 1 files changed, 10 insertions(+), 10 deletions(-) [+] |
line wrap: on
line diff
--- a/+scheme/Elastic2dVariable.m Fri Nov 02 10:47:03 2018 -0700 +++ b/+scheme/Elastic2dVariable.m Wed Nov 28 14:04:31 2018 -0800 @@ -515,21 +515,21 @@ LAMBDA = obj.LAMBDA; MU = obj.MU; - phi = obj.phi{j}; - h = obj.h(j); - h11 = obj.H11{j}*h; - gamma = obj.gamma{j}; dim = obj.dim; + theta_R = obj.theta_R{j}; + theta_H = obj.theta_H{j}; + theta_M = obj.theta_M{j}; - a_lambda = dim/h11 + 1/(h11*phi); - a_mu_i = 2/(gamma*h); - a_mu_ij = 2/h11 + 1/(h11*phi); + a_lambda = dim/theta_H + 1/theta_R; + a_mu_i = 2/theta_M; + a_mu_ij = 2/theta_H + 1/theta_R; d = @kroneckerDelta; % Kronecker delta db = @(i,j) 1-d(i,j); % Logical not of Kronecker delta - alpha = @(i,k) d(i,k)*tuning*( d(i,j)* a_lambda*LAMBDA ... - + d(i,j)* a_mu_i*MU ... - + db(i,j)*a_mu_ij*MU ); + alpha = @(i,j) tuning*( d(i,j)* a_lambda*LAMBDA ... + + d(i,j)* a_mu_i*MU ... + + db(i,j)*a_mu_ij*MU ); + varargout{i} = alpha; otherwise error(['No such operator: operator = ' op{i}]);