applications/solvers/solidMechanics/deprecatedSolvers/contactStressFoam/createFields.H

   1     mesh.weights();
   2
   3     Info<< "Reading field U\n" << endl;
   4     volVectorField U
   5     (
   6         IOobject
   7         (
   8             "U",
   9             runTime.timeName(),
  10             mesh,
  11             IOobject::MUST_READ,
  12             IOobject::AUTO_WRITE
  13         ),
  14         mesh
  15     );
  16
  17     Info<< "Reading material properties\n" << endl;
  18
  19     IOdictionary mechanicalProperties
  20     (
  21         IOobject
  22         (
  23             "mechanicalProperties",
  24             "constant",
  25             mesh,
  26             IOobject::MUST_READ,
  27             IOobject::NO_WRITE
  28         )
  29     );
  30
  31     dimensionedScalar rho
  32     (
  33         mechanicalProperties.lookup("rho")
  34     );
  35
  36     dimensionedScalar rhoE
  37     (
  38         mechanicalProperties.lookup("E")
  39     );
  40
  41     dimensionedScalar nu
  42     (
  43         mechanicalProperties.lookup("nu")
  44     );
  45
  46     Info<< "Normalising E : E/rho\n" << endl;
  47     dimensionedScalar E = rhoE/rho;
  48
  49     Info<< "Reading boundary conditions\n" << endl;
  50
  51     IOdictionary boundaryConditions
  52     (
  53         IOobject
  54         (
  55             "boundaryConditions",
  56             "constant",
  57             mesh,
  58             IOobject::MUST_READ,
  59             IOobject::NO_WRITE
  60         )
  61     );
  62
  63     Info<< "Reading patchPressure\n" << endl;
  64     scalarField rhoPatchPressure
  65     (
  66         boundaryConditions.lookup("patchPressure")
  67     );
  68
  69     Info<< "Normalising patchPressure : patchPressure/rho\n" << endl;
  70     scalarField patchPressure = rhoPatchPressure/rho.value();
  71
  72     Info<< "Reading patchTraction\n" << endl;
  73     vectorField rhoPatchTraction
  74     (
  75         boundaryConditions.lookup("patchTraction")
  76     );
  77
  78     Info<< "Normalising patchTraction : patchTraction/rho\n" << endl;
  79     vectorField patchTraction = rhoPatchTraction/rho.value();
  80
  81     Info<< "Calculating Lame's coefficients\n" << endl;
  82
  83     dimensionedScalar mu = E/(2.0*(1.0 + nu));
  84     dimensionedScalar lambda = nu*E/((1.0 + nu)*(1.0 - 2.0*nu));
  85     dimensionedScalar threeK = E/(1.0 - 2.0*nu);
  86
  87     Switch planeStress
  88     (
  89         mechanicalProperties.lookup("planeStress")
  90     );
  91
  92     if (planeStress)
  93     {
  94         Info<< "Plane Stress\n" << endl;
  95
  96         //- change lambda and threeK for plane stress
  97         lambda = nu*E/((1.0 + nu)*(1.0 - nu));
  98     }
  99     else
 100     {
 101         Info<< "Plane Strain\n" << endl;
 102     }
 103
 104     Info<< "mu = " << mu.value() << " Pa/rho\n";
 105     Info<< "lambda = " << lambda.value() << " Pa/rho\n";
 106
 107 #   include "calculatePatchTraction.H"
 108
 109     Info<< "Reading contact patch IDs etc.\n" << endl;
 110     const dictionary& contactPatchDict
 111     (
 112      mesh.solutionDict().subDict("contactPatch")
 113     );
 114     label gradPatch(readLabel(contactPatchDict.lookup("gradPatchID")));
 115     label dirPatch(readLabel(contactPatchDict.lookup("dirPatchID")));
 116     scalar urf(readScalar(contactPatchDict.lookup("contactUnderrelaxation")));
 117     scalar touchTolerance
 118     (
 119         readScalar(contactPatchDict.lookup("touchTolerance"))
 120     );
 121
 122     Info << "Creating contact pairs" << endl;
 123     contactPatchPair contactPair(mesh, gradPatch, dirPatch, touchTolerance);
 124
 125     contactPatchPair reversePair(mesh, dirPatch, gradPatch, touchTolerance);
 126
 127     volScalarField contactArea
 128     (
 129         IOobject
 130         (
 131             "contactArea",
 132             runTime.timeName(),
 133             mesh,
 134             IOobject::NO_READ,
 135             IOobject::AUTO_WRITE
 136         ),
 137         mesh,
 138         dimensionedScalar(0)
 139     );
 140