[cig-commits] r7267 - short/3D/PyLith/trunk/libsrc/feassemble
brad at geodynamics.org
brad at geodynamics.org
Fri Jun 15 13:37:38 PDT 2007
Author: brad
Date: 2007-06-15 13:37:38 -0700 (Fri, 15 Jun 2007)
New Revision: 7267
Modified:
short/3D/PyLith/trunk/libsrc/feassemble/ElasticityImplicit.cc
Log:
Added comment. No change to code.
Modified: short/3D/PyLith/trunk/libsrc/feassemble/ElasticityImplicit.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/feassemble/ElasticityImplicit.cc 2007-06-15 20:11:09 UTC (rev 7266)
+++ short/3D/PyLith/trunk/libsrc/feassemble/ElasticityImplicit.cc 2007-06-15 20:37:38 UTC (rev 7267)
@@ -404,6 +404,10 @@
// Compute Jacobian for consistent tangent matrix
for (int iQuad=0; iQuad < numQuadPts; ++iQuad) {
const double wt = quadWts[iQuad] * jacobianDet[iQuad];
+ // tau_ij = C_ijkl * e_kl
+ // = C_ijlk * 0.5 (u_k,l + u_l,k)
+ // = 0.5 * C_ijkl * (u_k,l + u_l,k)
+ // divide C_ijkl by 0.5 if k != l
const double C1111 = elasticConsts[iQuad][0];
const double C1122 = elasticConsts[iQuad][1];
const double C1112 = elasticConsts[iQuad][2]/2.0;
@@ -456,6 +460,10 @@
// Compute Jacobian for consistent tangent matrix
for (int iQuad=0; iQuad < numQuadPts; ++iQuad) {
const double wt = quadWts[iQuad] * jacobianDet[iQuad];
+ // tau_ij = C_ijkl * e_kl
+ // = C_ijlk * 0.5 (u_k,l + u_l,k)
+ // = 0.5 * C_ijkl * (u_k,l + u_l,k)
+ // divide C_ijkl by 0.5 if k != l
const double C1111 = elasticConsts[iQuad][ 0];
const double C1122 = elasticConsts[iQuad][ 1];
const double C1133 = elasticConsts[iQuad][ 2];
More information about the cig-commits
mailing list