[cig-commits] r15692 - short/3D/PyLith/branches/pylith-friction/libsrc/faults

brad at geodynamics.org brad at geodynamics.org
Fri Sep 25 17:17:59 PDT 2009


Author: brad
Date: 2009-09-25 17:17:59 -0700 (Fri, 25 Sep 2009)
New Revision: 15692

Modified:
   short/3D/PyLith/branches/pylith-friction/libsrc/faults/FaultCohesiveDynL.cc
   short/3D/PyLith/branches/pylith-friction/libsrc/faults/FaultCohesiveDynL.hh
Log:
Added some notes on friction.

Modified: short/3D/PyLith/branches/pylith-friction/libsrc/faults/FaultCohesiveDynL.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/libsrc/faults/FaultCohesiveDynL.cc	2009-09-25 21:42:51 UTC (rev 15691)
+++ short/3D/PyLith/branches/pylith-friction/libsrc/faults/FaultCohesiveDynL.cc	2009-09-26 00:17:59 UTC (rev 15692)
@@ -623,7 +623,22 @@
 				       const double t,
 				       topology::SolutionFields* const fields)
 { // constrainSolution
-  // :TODO: Friction stuff goes here
+
+  // Loop over Lagrange vertices
+
+  // Compute current traction from Lagrange multiplier and area.
+
+  // Compute the current slip from current displacements.
+
+  // Use fault constitutive model to compute traction associated with
+  // friction.
+  
+  // Adjust Lagrange multiplier values accordingly.
+
+  // Update the slip estimate based on adjustement to the Lagrange
+  // multiplier values.
+
+
 } // constrainSolution
 
 // ----------------------------------------------------------------------

Modified: short/3D/PyLith/branches/pylith-friction/libsrc/faults/FaultCohesiveDynL.hh
===================================================================
--- short/3D/PyLith/branches/pylith-friction/libsrc/faults/FaultCohesiveDynL.hh	2009-09-25 21:42:51 UTC (rev 15691)
+++ short/3D/PyLith/branches/pylith-friction/libsrc/faults/FaultCohesiveDynL.hh	2009-09-26 00:17:59 UTC (rev 15692)
@@ -179,7 +179,7 @@
 
   /** Constrain solution based on friction.
    *
-   * @param t Current time
+   * @param t Current time.
    * @param fields Solution fields
    */
   void constrainSolution(const double t,



More information about the CIG-COMMITS mailing list