Mercurial > repos > public > sbplib

comparison +scheme/Elastic2dCurvilinearAnisotropic.m @ 1211:3d7faa2ca312 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add physical derivatives Dx and Dy to ElasticCurvilinearAnisotropic
author Martin Almquist <malmquist@stanford.edu>
date Fri, 27 Sep 2019 13:18:09 -0700
parents 67eee83fd9c9
children 43f1cd11e8e8
comparison
equal deleted inserted replaced
1209:67eee83fd9c9 1211:3d7faa2ca312
20 J, Ji 20 J, Ji
21 RHO % Density 21 RHO % Density
22 C % Elastic stiffness tensor 22 C % Elastic stiffness tensor
23 23
24 D % Total operator 24 D % Total operator
25
26 Dx, Dy % Physical derivatives
25 27
26 % Boundary operators in cell format, used for BC 28 % Boundary operators in cell format, used for BC
27 T_w, T_e, T_s, T_n 29 T_w, T_e, T_s, T_n
28 30
29 % Traction operators 31 % Traction operators
152 K = cell(dim, dim); 154 K = cell(dim, dim);
153 K{1,1} = y_v./J; 155 K{1,1} = y_v./J;
154 K{1,2} = -y_u./J; 156 K{1,2} = -y_u./J;
155 K{2,1} = -x_v./J; 157 K{2,1} = -x_v./J;
156 K{2,2} = x_u./J; 158 K{2,2} = x_u./J;
159
160 % Physical derivatives
161 obj.Dx = spdiag( y_v./J)*Du + spdiag(-y_u./J)*Dv;
162 obj.Dy = spdiag(-x_v./J)*Du + spdiag( x_u./J)*Dv;
157 163
158 % Wrap around Aniosotropic Cartesian 164 % Wrap around Aniosotropic Cartesian
159 rho_tilde = J.*rho; 165 rho_tilde = J.*rho;
160 166
161 PHI = cell(dim,dim,dim,dim); 167 PHI = cell(dim,dim,dim,dim);