src/postProcessing/functionObjects/forces/pointPatchFields/derived/sixDoFRigidBodyDisplacement/sixDoFRigidBodyDisplacementPointPatchVectorField.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 "sixDoFRigidBodyDisplacementPointPatchVectorField.H"
  28 #include "pointPatchFields.H"
  29 #include "addToRunTimeSelectionTable.H"
  30 #include "Time.H"
  31 #include "fvMesh.H"
  32 #include "volFields.H"
  33 #include "uniformDimensionedFields.H"
  34 #include "forces.H"
  35
  36 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  37
  38 namespace Foam
  39 {
  40
  41 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  42
  43 sixDoFRigidBodyDisplacementPointPatchVectorField::
  44 sixDoFRigidBodyDisplacementPointPatchVectorField
  45 (
  46     const pointPatch& p,
  47     const DimensionedField<vector, pointMesh>& iF
  48 )
  49 :
  50     fixedValuePointPatchField<vector>(p, iF),
  51     motion_(),
  52     initialPoints_(p.localPoints()),
  53     rhoInf_(1.0)
  54 {}
  55
  56
  57 sixDoFRigidBodyDisplacementPointPatchVectorField::
  58 sixDoFRigidBodyDisplacementPointPatchVectorField
  59 (
  60     const pointPatch& p,
  61     const DimensionedField<vector, pointMesh>& iF,
  62     const dictionary& dict
  63 )
  64 :
  65     fixedValuePointPatchField<vector>(p, iF, dict),
  66     motion_(dict),
  67     rhoInf_(readScalar(dict.lookup("rhoInf")))
  68 {
  69     if (!dict.found("value"))
  70     {
  71         updateCoeffs();
  72     }
  73
  74     if (dict.found("initialPoints"))
  75     {
  76         initialPoints_ = vectorField("initialPoints", dict , p.size());
  77     }
  78     else
  79     {
  80         initialPoints_ = p.localPoints();
  81     }
  82 }
  83
  84
  85 sixDoFRigidBodyDisplacementPointPatchVectorField::
  86 sixDoFRigidBodyDisplacementPointPatchVectorField
  87 (
  88     const sixDoFRigidBodyDisplacementPointPatchVectorField& ptf,
  89     const pointPatch& p,
  90     const DimensionedField<vector, pointMesh>& iF,
  91     const pointPatchFieldMapper& mapper
  92 )
  93 :
  94     fixedValuePointPatchField<vector>(ptf, p, iF, mapper),
  95     motion_(ptf.motion_),
  96     initialPoints_(ptf.initialPoints_, mapper),
  97     rhoInf_(ptf.rhoInf_)
  98 {}
  99
 100
 101 sixDoFRigidBodyDisplacementPointPatchVectorField::
 102 sixDoFRigidBodyDisplacementPointPatchVectorField
 103 (
 104     const sixDoFRigidBodyDisplacementPointPatchVectorField& ptf,
 105     const DimensionedField<vector, pointMesh>& iF
 106 )
 107 :
 108     fixedValuePointPatchField<vector>(ptf, iF),
 109     motion_(ptf.motion_),
 110     initialPoints_(ptf.initialPoints_),
 111     rhoInf_(ptf.rhoInf_)
 112 {}
 113
 114
 115 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 116
 117 void sixDoFRigidBodyDisplacementPointPatchVectorField::autoMap
 118 (
 119     const pointPatchFieldMapper& m
 120 )
 121 {
 122     fixedValuePointPatchField<vector>::autoMap(m);
 123
 124     initialPoints_.autoMap(m);
 125 }
 126
 127
 128 void sixDoFRigidBodyDisplacementPointPatchVectorField::rmap
 129 (
 130     const pointPatchField<vector>& ptf,
 131     const labelList& addr
 132 )
 133 {
 134     const sixDoFRigidBodyDisplacementPointPatchVectorField& sDoFptf =
 135         refCast<const sixDoFRigidBodyDisplacementPointPatchVectorField>(ptf);
 136
 137     fixedValuePointPatchField<vector>::rmap(sDoFptf, addr);
 138
 139     initialPoints_.rmap(sDoFptf.initialPoints_, addr);
 140 }
 141
 142
 143 void sixDoFRigidBodyDisplacementPointPatchVectorField::updateCoeffs()
 144 {
 145     if (this->updated())
 146     {
 147         return;
 148     }
 149
 150     const polyMesh& mesh = this->dimensionedInternalField().mesh()();
 151     const Time& t = mesh.time();
 152     const pointPatch& ptPatch = this->patch();
 153
 154     // Patch force data is valid for the current positions, so
 155     // calculate the forces on the motion object from this data, then
 156     // update the positions
 157
 158     motion_.updatePosition(t.deltaTValue());
 159
 160     dictionary forcesDict;
 161
 162     forcesDict.add("patches", wordList(1, ptPatch.name()));
 163     forcesDict.add("rhoInf", rhoInf_);
 164     forcesDict.add("CofR", motion_.centreOfMass());
 165
 166     forces f("forces", db(), forcesDict);
 167
 168     forces::forcesMoments fm = f.calcForcesMoment();
 169
 170     // Get the forces on the patch faces at the current positions
 171
 172     vector gravity = vector::zero;
 173
 174     if (db().foundObject<uniformDimensionedVectorField>("g"))
 175     {
 176         uniformDimensionedVectorField g =
 177             db().lookupObject<uniformDimensionedVectorField>("g");
 178
 179         gravity = g.value();
 180     }
 181
 182     motion_.updateForce
 183     (
 184         fm.first().first() + fm.first().second() + gravity*motion_.mass(),
 185         fm.second().first() + fm.second().second(),
 186         t.deltaTValue()
 187     );
 188
 189     Field<vector>::operator=
 190     (
 191         motion_.currentPosition(initialPoints_) - initialPoints_
 192     );
 193
 194     fixedValuePointPatchField<vector>::updateCoeffs();
 195 }
 196
 197
 198 void sixDoFRigidBodyDisplacementPointPatchVectorField::write(Ostream& os) const
 199 {
 200     pointPatchField<vector>::write(os);
 201     motion_.write(os);
 202     os.writeKeyword("rhoInf")
 203         << rhoInf_ << token::END_STATEMENT << nl;
 204     initialPoints_.writeEntry("initialPoints", os);
 205     writeEntry("value", os);
 206 }
 207
 208
 209 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 210
 211 makePointPatchTypeField
 212 (
 213     pointPatchVectorField,
 214     sixDoFRigidBodyDisplacementPointPatchVectorField
 215 );
 216
 217 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 218
 219 } // End namespace Foam
 220
 221 // ************************************************************************* //