Comments added about including mu in weighted volume and dilatation c…

…alculations
ORNL · Aug 17, 2023 · ab4ab6a · ab4ab6a
1 parent cdeca2b
commit ab4ab6a
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/src/CabanaPD_Force.hpp b/src/CabanaPD_Force.hpp
@@ -588,6 +588,7 @@ class Force<ExecutionSpace, ForceModel<LPS, Fracture>>
                 // Get the reference positions and displacements.
                 double xi, r, s;
                 getDistance( x, u, i, j, xi, r, s );
+                // mu is included to account for bond breaking.
                 double m_j = mu( i, n ) * model.influence_function( xi ) * xi *
                              xi * vol( j );
                 m( i ) += m_j;
@@ -628,6 +629,7 @@ class Force<ExecutionSpace, ForceModel<LPS, Fracture>>
                 // Get the bond distance, displacement, and stretch.
                 double xi, r, s;
                 getDistance( x, u, i, j, xi, r, s );
+                // mu is included to account for bond breaking.
                 double theta_i = mu( i, n ) * model.influence_function( xi ) *
                                  s * xi * xi * vol( j );
                 // Check if all bonds are broken (m=0) to avoid dividing by