src/ODE/sixDOF/sixDOFbodies/sixDOFbodies.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 Class
  25     sixDOFbodies
  26
  27 Description
  28     6-DOF solver for multiple bodies
  29
  30 Author
  31     Dubravko Matijasevic, FSB Zagreb.  All rights reserved.
  32
  33 \*---------------------------------------------------------------------------*/
  34
  35 #include "objectRegistry.H"
  36 #include "sixDOFbodies.H"
  37
  38 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  39
  40 void Foam::sixDOFbodies::setBodies()
  41 {
  42     // Find if duplicate name existes
  43     forAll (names_, bodyI)
  44     {
  45         for
  46         (
  47             label otherBody = bodyI + 1;
  48             otherBody < names_.size();
  49             otherBody++
  50         )
  51         {
  52             if (names_[bodyI] == names_[otherBody])
  53             {
  54                 FatalErrorIn("sixDOFbodies::setBodies()")
  55                     << "Duplicate names of bodies: this is not allowed"
  56                     << exit(FatalError);
  57             }
  58         }
  59     }
  60
  61     odes_.setSize(names_.size());
  62     solvers_.setSize(names_.size());
  63
  64     forAll (names_, bodyI)
  65     {
  66         odes_.set
  67         (
  68             bodyI,
  69             new sixDOFqODE
  70             (
  71                 IOobject
  72                 (
  73                     names_[bodyI],
  74                     runTime_.timeName(),
  75                     runTime_,
  76                     IOobject::MUST_READ,
  77                     IOobject::AUTO_WRITE
  78                 )
  79             )
  80         );
  81
  82         solvers_.set
  83         (
  84             bodyI,
  85             ODESolver::New
  86             (
  87                 lookup("solver"),
  88                 odes_[bodyI]
  89             )
  90         );
  91     }
  92 }
  93
  94
  95 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  96
  97 Foam::sixDOFbodies::sixDOFbodies
  98 (
  99     const Time& runTime
 100 )
 101 :
 102     IOdictionary
 103     (
 104         IOobject
 105         (
 106             "sixDOFsolverDict",
 107             runTime.system(),
 108             runTime,
 109             IOobject::MUST_READ,
 110             IOobject::NO_WRITE
 111         )
 112     ),
 113     runTime_(runTime),
 114     odes_(),
 115     solvers_(),
 116     names_(lookup("bodies"))
 117 {
 118     setBodies();
 119 }
 120
 121
 122 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 123
 124 void Foam::sixDOFbodies::solve()
 125 {
 126     const scalar tol = readScalar(lookup("eps"));
 127
 128     forAll (odes_, bodyI)
 129     {
 130         Info << "Solving 6-DOF for " << names_[bodyI] << " in time"
 131          << tab << "T = " << runTime_.value() << " s" << endl;
 132
 133         solvers_[bodyI].solve
 134         (
 135             runTime_.value(),
 136             runTime_.value() + runTime_.deltaT().value(),
 137             tol,
 138             runTime_.deltaT().value()
 139         );
 140     }
 141 }
 142
 143
 144 const Foam::wordList& Foam::sixDOFbodies::names() const
 145 {
 146     return names_;
 147 }
 148
 149
 150 const Foam::PtrList<Foam::sixDOFqODE>& Foam::sixDOFbodies::operator()() const
 151 {
 152     return odes_;
 153 }
 154
 155
 156 // ************************************************************************* //