applications/solvers/solidMechanics/deprecatedSolvers/contactStressFoam/contactStressFoam.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | foam-extend: Open Source CFD
   4    \\    /   O peration     | Version:     3.2
   5     \\  /    A nd           | Web:         http://www.foam-extend.org
   6      \\/     M anipulation  | For copyright notice see file Copyright
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of foam-extend.
  10
  11     foam-extend is free software: you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by the
  13     Free Software Foundation, either version 3 of the License, or (at your
  14     option) any later version.
  15
  16     foam-extend is distributed in the hope that it will be useful, but
  17     WITHOUT ANY WARRANTY; without even the implied warranty of
  18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19     General Public License for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 Application
  25     stressedFoam
  26
  27 Description
  28     Transient/steady-state solver of linear-elastic, small-strain deformation
  29     of solid bodies in contact.
  30
  31     Simple linear elasticity structural analysis code.
  32     Solves for the displacement vector field U,
  33     also generating the stress tensor field sigma.
  34
  35     Code can be compiled for simple stress analysis,
  36     and dynamic stress analysis.
  37
  38     Additionally, special handling of very non-orthogonal meshes may be added.
  39
  40 \*---------------------------------------------------------------------------*/
  41
  42 #include "fvCFD.H"
  43 #include "fixedGradientFvPatchFields.H"
  44 #include "directionMixedFvPatchFields.H"
  45 #include "contactPatchPair.H"
  46
  47 #include "pointFields.H"
  48 #include "transformField.H"
  49
  50 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  51
  52 int main(int argc, char *argv[])
  53 {
  54
  55 #   include "setRootCase.H"
  56
  57 #   include "createTime.H"
  58 #   include "createMesh.H"
  59 #   include "createFields.H"
  60
  61 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  62
  63     Info<< "\nCalculating displacement field\n" << endl;
  64
  65     volTensorField gradU = fvc::grad(U);
  66
  67     for (runTime++; !runTime.end(); runTime++)
  68     {
  69         Info<< "Iteration: " << runTime.timeName() << nl << endl;
  70
  71 #       include "../stressedFoam/readStressedFoamControls.H"
  72
  73         int iCorr=0;
  74         scalar initialResidual=0;
  75
  76 #       include "contactBoundaries.H"
  77
  78         do
  79         {
  80 #           include "tractionBoundaries.H"
  81
  82             fvVectorMatrix UEqn
  83             (
  84 #               ifdef Dynamic
  85                 fvm::d2dt2(U)
  86               ==
  87 #               endif
  88
  89                 fvm::laplacian(2*mu + lambda, U, "laplacian(DU,U)")
  90
  91               + fvc::div
  92                 (
  93                     mu*gradU.T() + lambda*(I*tr(gradU)) - (mu + lambda)*gradU,
  94                     "div(sigma)"
  95                 )
  96             );
  97
  98             UEqn.setComponentReference(6, 0, vector::Z, 0);
  99
 100             initialResidual = UEqn.solve().initialResidual();
 101
 102             gradU = fvc::grad(U);
 103
 104         } while (initialResidual > convergenceTolerance && ++iCorr < nCorr);
 105
 106 #       include "calculateStress.H"
 107
 108         Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
 109             << "  ClockTime = " << runTime.elapsedClockTime() << " s"
 110             << nl << endl;
 111     }
 112
 113     Info<< "End\n" << endl;
 114
 115     return(0);
 116 }
 117
 118
 119 // ************************************************************************* //