]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Update to the step-42 documentation. 1723/head
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Tue, 6 Oct 2015 21:46:47 +0000 (16:46 -0500)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Tue, 6 Oct 2015 21:46:47 +0000 (16:46 -0500)
This implements a suggestion by Korosh Taebi.

examples/step-42/doc/intro.dox
examples/step-42/step-42.cc

index 0a2397871e1fa4c12c5013ec518430ae841a9365..93573cbff0a8058579c93e7bfc4da85ad3b8968c 100644 (file)
@@ -169,6 +169,8 @@ condition:
   \left\{ \mathbf u\in V: \mathbf n \cdot \mathbf u\leq g \text{ on } \Gamma_C \right\}.
 @f}
 
+In the actual code, we will use the abbreviation $\gamma=\dfrac{\gamma^{\text{iso}}}{2\mu + \gamma^{\text{iso}}}$.
+
 Given this formulation, we will apply two techniques:
 - Run a Newton method to iterate out the nonlinearity in the projector.
 - Run an active set method for the contact condition, in much the same
index 077e1ee2294e231b4acc58e3733f74dc61ae6e12..54ebb04827f739fe10d13d614076e8ce25a2bf1b 100644 (file)
@@ -176,12 +176,18 @@ namespace Step42
 
   // @sect4{ConstitutiveLaw::get_stress_strain_tensor}
 
-  // This is the principal component of the constitutive law. It projects the
-  // deviatoric part of the stresses in a quadrature point back to the yield
-  // stress (i.e., the original yield stress $\sigma_0$ plus the term that
-  // describes linear isotropic hardening).  We need this function to calculate
-  // the nonlinear residual in PlasticityContactProblem::residual_nl_system. The
-  // computations follow the formulas laid out in the introduction.
+  // This is the principal component of the constitutive law. It
+  // computes the fourth order symmetric tensor that relates the
+  // strain to the stress according to the projection given above,
+  // when evaluated at a particular strain point. We need this
+  // function to calculate the nonlinear residual in
+  // <code>PlasticityContactProblem::residual_nl_system()</code> where
+  // we multiply this tensor with the strain given in a quadrature
+  // point. The computations follow the formulas laid out in the
+  // introduction. In comparing the formulas there with the
+  // implementation below, recall that $C_\mu : \varepsilon = \tau_D$
+  // and that $C_\kappa : \varepsilon = \kappa
+  // \text{trace}(\varepsilon) I = \frac 13 \text{trace}(\tau) I$.
   //
   // The function returns whether the quadrature point is plastic to allow for
   // some statistics downstream on how many of the quadrature points are

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.