Mercurial > repos > public > sbplib

comparison +scheme/Elastic2dVariable.m @ 967:368a2773f78b feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix bug in Elastic2dVariable that resulted in using tuning^2
author Martin Almquist <malmquist@stanford.edu>
date Sun, 23 Dec 2018 14:06:26 +0100
parents 4a85cff53107
children a4ad90b37998
comparison
equal deleted inserted replaced
966:4a85cff53107 967:368a2773f78b
481 alpha_ij = mu/(2*th_H) + mu/(4*th_R); 481 alpha_ij = mu/(2*th_H) + mu/(4*th_R);
482 end 482 end
483 483
484 [alpha_ii_u, alpha_ij_u] = computeAlpha(theta_R_u, theta_H_u, theta_M_u, lambda_u, mu_u); 484 [alpha_ii_u, alpha_ij_u] = computeAlpha(theta_R_u, theta_H_u, theta_M_u, lambda_u, mu_u);
485 [alpha_ii_v, alpha_ij_v] = computeAlpha(theta_R_v, theta_H_v, theta_M_v, lambda_v, mu_v); 485 [alpha_ii_v, alpha_ij_v] = computeAlpha(theta_R_v, theta_H_v, theta_M_v, lambda_v, mu_v);
486 sigma_ii = tuning*(alpha_ii_u + alpha_ii_v); 486 sigma_ii = alpha_ii_u + alpha_ii_v;
487 sigma_ij = tuning*(alpha_ij_u + alpha_ij_v); 487 sigma_ij = alpha_ij_u + alpha_ij_v;
488 488
489 d = @kroneckerDelta; % Kronecker delta 489 d = @kroneckerDelta; % Kronecker delta
490 db = @(i,j) 1-d(i,j); % Logical not of Kronecker delta 490 db = @(i,j) 1-d(i,j); % Logical not of Kronecker delta
491 sigma = @(i,j) tuning*(d(i,j)*sigma_ii + db(i,j)*sigma_ij); 491 sigma = @(i,j) tuning*(d(i,j)*sigma_ii + db(i,j)*sigma_ij);
492 492