Mercurial > repos > public > sbplib

comparison +scheme/Elastic2dVariable.m @ 730:9f28cf266f86 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Bugfix Elastic2dVariable interface
author Martin Almquist <malmquist@stanford.edu>
date Sat, 21 Apr 2018 16:04:46 -0700
parents aa8cf3851de8
children eebe24a636c7
comparison
equal deleted inserted replaced
729:aa8cf3851de8 730:9f28cf266f86
398 th3 = h*gamma; 398 th3 = h*gamma;
399 a1 = ( (th1 + th2)*th3*lambda + 4*th1*th2*mu ) / (2*th1*th2*th3); 399 a1 = ( (th1 + th2)*th3*lambda + 4*th1*th2*mu ) / (2*th1*th2*th3);
400 a2 = ( 16*(th1 + th2)*lambda*mu ) / (th1*th2*th3); 400 a2 = ( 16*(th1 + th2)*lambda*mu ) / (th1*th2*th3);
401 alpha_ii = a1 + sqrt(a2 + a1^2); 401 alpha_ii = a1 + sqrt(a2 + a1^2);
402 402
403 alpha_ij = 2/h11 + 1/(phi*h11); 403 alpha_ij = mu*(2/h11 + 1/(phi*h11));
404 end 404 end
405 405
406 [alpha_ii_u, alpha_ij_u] = computeAlpha(phi_u,h_u,h11_u,gamma_u,lambda_u,mu_u); 406 [alpha_ii_u, alpha_ij_u] = computeAlpha(phi_u,h_u,h11_u,gamma_u,lambda_u,mu_u);
407 [alpha_ii_v, alpha_ij_v] = computeAlpha(phi_v,h_v,h11_v,gamma_v,lambda_v,mu_v); 407 [alpha_ii_v, alpha_ij_v] = computeAlpha(phi_v,h_v,h11_v,gamma_v,lambda_v,mu_v);
408 sigma_ii = tuning*(alpha_ii_u + alpha_ii_v)/4; 408 sigma_ii = tuning*(alpha_ii_u + alpha_ii_v)/4;