applications/solvers/solidMechanics/deprecatedSolvers/stressedFoam/stressedFoam.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 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 U, 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 "readThermalProperties.H"
  52 #   include "createFields.H"
  53
  54 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  55
  56     Info<< "\nCalculating displacement field\n" << endl;
  57
  58     for (runTime++; !runTime.end(); runTime++)
  59     {
  60         Info<< "Iteration: " << runTime.timeName() << nl << endl;
  61
  62 #       include "readStressedFoamControls.H"
  63
  64         int iCorr = 0;
  65         scalar initialResidual = 0;
  66
  67         do
  68         {
  69             volTensorField gradU = fvc::grad(U);
  70
  71             if (thermalStress)
  72             {
  73                 solve
  74                 (
  75                     fvm::ddt(T) == fvm::laplacian(DT, T)
  76                 );
  77             }
  78
  79             fvVectorMatrix UEqn
  80             (
  81                 fvm::d2dt2(U)
  82              ==
  83                 fvm::laplacian(2*mu + lambda, U, "laplacian(DU,U)")
  84
  85               + fvc::div
  86                 (
  87                     mu*gradU.T() + lambda*(I*tr(gradU)) - (mu + lambda)*gradU,
  88                     "div(sigma)"
  89                 )
  90             );
  91
  92             if (thermalStress)
  93             {
  94                 UEqn += threeKalpha*fvc::grad(T);
  95             }
  96
  97             //UEqn.setComponentReference(1, 0, vector::X, 0);
  98             //UEqn.setComponentReference(1, 0, vector::Z, 0);
  99
 100             initialResidual = UEqn.solve().initialResidual();
 101
 102         } while (initialResidual > convergenceTolerance && ++iCorr < nCorr);
 103
 104 #       include "calculateStress.H"
 105
 106         Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
 107             << "  ClockTime = " << runTime.elapsedClockTime() << " s"
 108             << nl << endl;
 109     }
 110
 111     deleteDemandDrivenData(Tptr);
 112
 113     Info<< "End\n" << endl;
 114
 115     return(0);
 116 }
 117
 118
 119 // ************************************************************************* //