applications/solvers/stressAnalysis/solidEquilibriumDisplacementFoam/tractionDisplacementCorrection/tractionDisplacementCorrectionFvPatchVectorField.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2009 OpenCFD Ltd.
   6      \\/     M anipulation  |
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of OpenFOAM.
  10
  11     OpenFOAM 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 2 of the License, or (at your
  14     option) any later version.
  15
  16     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
  17     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  18     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  19     for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with OpenFOAM; if not, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 \*---------------------------------------------------------------------------*/
  26
  27 #include "tractionDisplacementCorrectionFvPatchVectorField.H"
  28 #include "addToRunTimeSelectionTable.H"
  29 #include "volFields.H"
  30
  31 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  32
  33 namespace Foam
  34 {
  35
  36 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  37
  38 tractionDisplacementCorrectionFvPatchVectorField::
  39 tractionDisplacementCorrectionFvPatchVectorField
  40 (
  41     const fvPatch& p,
  42     const DimensionedField<vector, volMesh>& iF
  43 )
  44 :
  45     fixedGradientFvPatchVectorField(p, iF),
  46     traction_(p.size(), vector::zero),
  47     pressure_(p.size(), 0.0)
  48 {
  49     fvPatchVectorField::operator=(patchInternalField());
  50     gradient() = vector::zero;
  51 }
  52
  53
  54 tractionDisplacementCorrectionFvPatchVectorField::
  55 tractionDisplacementCorrectionFvPatchVectorField
  56 (
  57     const tractionDisplacementCorrectionFvPatchVectorField& tdpvf,
  58     const fvPatch& p,
  59     const DimensionedField<vector, volMesh>& iF,
  60     const fvPatchFieldMapper& mapper
  61 )
  62 :
  63     fixedGradientFvPatchVectorField(tdpvf, p, iF, mapper),
  64     traction_(tdpvf.traction_, mapper),
  65     pressure_(tdpvf.pressure_, mapper)
  66 {}
  67
  68
  69 tractionDisplacementCorrectionFvPatchVectorField::
  70 tractionDisplacementCorrectionFvPatchVectorField
  71 (
  72     const fvPatch& p,
  73     const DimensionedField<vector, volMesh>& iF,
  74     const dictionary& dict
  75 )
  76 :
  77     fixedGradientFvPatchVectorField(p, iF),
  78     traction_("traction", dict, p.size()),
  79     pressure_("pressure", dict, p.size())
  80 {
  81     fvPatchVectorField::operator=(patchInternalField());
  82     gradient() = vector::zero;
  83 }
  84
  85
  86 tractionDisplacementCorrectionFvPatchVectorField::
  87 tractionDisplacementCorrectionFvPatchVectorField
  88 (
  89     const tractionDisplacementCorrectionFvPatchVectorField& tdpvf
  90 )
  91 :
  92     fixedGradientFvPatchVectorField(tdpvf),
  93     traction_(tdpvf.traction_),
  94     pressure_(tdpvf.pressure_)
  95 {}
  96
  97
  98 tractionDisplacementCorrectionFvPatchVectorField::
  99 tractionDisplacementCorrectionFvPatchVectorField
 100 (
 101     const tractionDisplacementCorrectionFvPatchVectorField& tdpvf,
 102     const DimensionedField<vector, volMesh>& iF
 103 )
 104 :
 105     fixedGradientFvPatchVectorField(tdpvf, iF),
 106     traction_(tdpvf.traction_),
 107     pressure_(tdpvf.pressure_)
 108 {}
 109
 110
 111 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 112
 113 void tractionDisplacementCorrectionFvPatchVectorField::autoMap
 114 (
 115     const fvPatchFieldMapper& m
 116 )
 117 {
 118     fixedGradientFvPatchVectorField::autoMap(m);
 119     traction_.autoMap(m);
 120     pressure_.autoMap(m);
 121 }
 122
 123
 124 // Reverse-map the given fvPatchField onto this fvPatchField
 125 void tractionDisplacementCorrectionFvPatchVectorField::rmap
 126 (
 127     const fvPatchVectorField& ptf,
 128     const labelList& addr
 129 )
 130 {
 131     fixedGradientFvPatchVectorField::rmap(ptf, addr);
 132
 133     const tractionDisplacementCorrectionFvPatchVectorField& dmptf =
 134         refCast<const tractionDisplacementCorrectionFvPatchVectorField>(ptf);
 135
 136     traction_.rmap(dmptf.traction_, addr);
 137     pressure_.rmap(dmptf.pressure_, addr);
 138 }
 139
 140
 141 // Update the coefficients associated with the patch field
 142 void tractionDisplacementCorrectionFvPatchVectorField::updateCoeffs()
 143 {
 144     if (updated())
 145     {
 146         return;
 147     }
 148
 149     const dictionary& mechanicalProperties = db().lookupObject<IOdictionary>
 150     (
 151         "mechanicalProperties"
 152     );
 153
 154     dimensionedScalar rho(mechanicalProperties.lookup("rho"));
 155     dimensionedScalar rhoE(mechanicalProperties.lookup("E"));
 156     dimensionedScalar nu(mechanicalProperties.lookup("nu"));
 157
 158     dimensionedScalar E = rhoE/rho;
 159     dimensionedScalar mu = E/(2.0*(1.0 + nu));
 160     dimensionedScalar lambda = nu*E/((1.0 + nu)*(1.0 - 2.0*nu));
 161
 162     Switch planeStress(mechanicalProperties.lookup("planeStress"));
 163
 164     if (planeStress)
 165     {
 166         lambda = nu*E/((1.0 + nu)*(1.0 - nu));
 167     }
 168
 169     vectorField n = patch().nf();
 170
 171     const fvPatchField<symmTensor>& sigmaD =
 172         patch().lookupPatchField<volSymmTensorField, symmTensor>("sigmaD");
 173
 174     const fvPatchField<tensor>& sigmaExp =
 175         patch().lookupPatchField<volTensorField, tensor>("sigmaExp");
 176
 177     gradient() =
 178     (
 179         (traction_ + pressure_*n)/rho.value() - (n & (sigmaD + sigmaExp))
 180     )/(2.0*mu + lambda).value();
 181
 182     fixedGradientFvPatchVectorField::updateCoeffs();
 183 }
 184
 185
 186 // Write
 187 void tractionDisplacementCorrectionFvPatchVectorField::write(Ostream& os) const
 188 {
 189     fvPatchVectorField::write(os);
 190     traction_.writeEntry("traction", os);
 191     pressure_.writeEntry("pressure", os);
 192     writeEntry("value", os);
 193 }
 194
 195
 196 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 197
 198 makePatchTypeField
 199 (
 200     fvPatchVectorField,
 201     tractionDisplacementCorrectionFvPatchVectorField
 202 );
 203
 204 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 205
 206 } // End namespace Foam
 207
 208 // ************************************************************************* //