applications/solvers/solidMechanics/deprecatedSolvers/solidEquilibriumDisplacementFoam/solidEquilibriumDisplacementFoam.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     solidEquilibriumDisplacementFoam
  26
  27 Description
  28     Steady-state segregated finite-volume solver of linear-elastic,
  29     small-strain deformation of a solid body, with optional thermal
  30     diffusion and thermal stresses.
  31
  32     Simple linear elasticity structural analysis code.
  33     Solves for the displacement vector field D, also generating the
  34     stress tensor field sigma.
  35
  36 \*---------------------------------------------------------------------------*/
  37
  38 #include "fvCFD.H"
  39 #include "Switch.H"
  40
  41 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  42
  43 int main(int argc, char *argv[])
  44 {
  45
  46 #   include "setRootCase.H"
  47
  48 #   include "createTime.H"
  49 #   include "createMesh.H"
  50 #   include "readMechanicalProperties.H"
  51 #   include "createFields.H"
  52
  53 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  54
  55     Info<< "\nCalculating displacement field\n" << endl;
  56
  57     while (runTime.loop())
  58     {
  59         Info<< "Iteration: " << runTime.value() << nl << endl;
  60
  61 #       include "readSteadyStressFoamControls.H"
  62
  63         solve
  64         (
  65             fvm::laplacian(2*mu + lambda, Dcorr, "laplacian(DD,Dcorr)")
  66           + fvc::div(sigmaExp + sigmaD)
  67         );
  68
  69         D += accFac*Dcorr;
  70
  71         {
  72             volTensorField gradDcorr = fvc::grad(Dcorr);
  73
  74             sigmaExp =
  75                 (lambda - mu)*gradDcorr + mu*gradDcorr.T()
  76               + (lambda*I)*tr(gradDcorr);
  77
  78             sigmaD += accFac*(mu*twoSymm(gradDcorr) + (lambda*I)*tr(gradDcorr));
  79         }
  80
  81 #       include "calculateStress.H"
  82 #       include "kineticEnergyLimiter.H"
  83
  84         Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
  85             << "  ClockTime = " << runTime.elapsedClockTime() << " s"
  86             << nl << endl;
  87     }
  88
  89     Info<< "\n end \n";
  90
  91     return 0;
  92 }
  93
  94
  95 // ************************************************************************* //