Mercurial > repos > public > sbplib

comparison +scheme/Elastic2dVariable.m @ 966:4a85cff53107 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix bug in Elastic2dVariable that caused crash for m_x neq m_y
author Martin Almquist <malmquist@stanford.edu>
date Sat, 22 Dec 2018 18:08:51 +0100
parents db3411264b96
children 368a2773f78b
comparison
equal deleted inserted replaced
965:db3411264b96 966:4a85cff53107
262 for i = 1:dim 262 for i = 1:dim
263 tau_l{j}{i} = sparse(n_l, dim*m_tot); 263 tau_l{j}{i} = sparse(n_l, dim*m_tot);
264 tau_r{j}{i} = sparse(n_r, dim*m_tot); 264 tau_r{j}{i} = sparse(n_r, dim*m_tot);
265 for k = 1:dim 265 for k = 1:dim
266 T_l{j}{i,k} = ... 266 T_l{j}{i,k} = ...
267 -d(i,j)*LAMBDA_l*(d(i,k)*d1_l{k}' + db(i,k)*e_l{j}'*D1{k})... 267 -d(i,j)*LAMBDA_l*(d(i,k)*d1_l{j}' + db(i,k)*e_l{j}'*D1{k})...
268 -d(j,k)*MU_l*(d(i,j)*d1_l{i}' + db(i,j)*e_l{j}'*D1{i})... 268 -d(j,k)*MU_l*(d(i,j)*d1_l{j}' + db(i,j)*e_l{j}'*D1{i})...
269 -d(i,k)*MU_l*d1_l{j}'; 269 -d(i,k)*MU_l*d1_l{j}';
270 270
271 T_r{j}{i,k} = ... 271 T_r{j}{i,k} = ...
272 d(i,j)*LAMBDA_r*(d(i,k)*d1_r{k}' + db(i,k)*e_r{j}'*D1{k})... 272 d(i,j)*LAMBDA_r*(d(i,k)*d1_r{j}' + db(i,k)*e_r{j}'*D1{k})...
273 +d(j,k)*MU_r*(d(i,j)*d1_r{i}' + db(i,j)*e_r{j}'*D1{i})... 273 +d(j,k)*MU_r*(d(i,j)*d1_r{j}' + db(i,j)*e_r{j}'*D1{i})...
274 +d(i,k)*MU_r*d1_r{j}'; 274 +d(i,k)*MU_r*d1_r{j}';
275 275
276 tau_l{j}{i} = tau_l{j}{i} + T_l{j}{i,k}*E{k}'; 276 tau_l{j}{i} = tau_l{j}{i} + T_l{j}{i,k}*E{k}';
277 tau_r{j}{i} = tau_r{j}{i} + T_r{j}{i,k}*E{k}'; 277 tau_r{j}{i} = tau_r{j}{i} + T_r{j}{i,k}*E{k}';
278 end 278 end