src/dynamicFvMesh/solidBodyMotionFvMesh/solidBodyMotionFunctions/SKA/SKA.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2008 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 "SKA.H"
  28 #include "addToRunTimeSelectionTable.H"
  29 #include "Tuple2.H"
  30 #include "IFstream.H"
  31 #include "interpolateXY.H"
  32 #include "mathematicalConstants.H"
  33
  34 using namespace Foam::mathematicalConstant;
  35
  36 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  37
  38 namespace Foam
  39 {
  40 namespace solidBodyMotionFunctions
  41 {
  42     defineTypeNameAndDebug(SKA, 0);
  43     addToRunTimeSelectionTable(solidBodyMotionFunction, SKA, dictionary);
  44 };
  45 };
  46
  47
  48 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  49
  50 Foam::solidBodyMotionFunctions::SKA::SKA
  51 (
  52     const dictionary& SBMFCoeffs,
  53     const Time& runTime
  54 )
  55 :
  56     solidBodyMotionFunction(SBMFCoeffs, runTime)
  57 {
  58     read(SBMFCoeffs);
  59 }
  60
  61
  62 // * * * * * * * * * * * * * * * * Destructors * * * * * * * * * * * * * * * //
  63
  64 Foam::solidBodyMotionFunctions::SKA::~SKA()
  65 {}
  66
  67
  68 // * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * * //
  69
  70 Foam::septernion Foam::solidBodyMotionFunctions::SKA::transformation() const
  71 {
  72     scalar t = time_.value();
  73
  74     if (t < times_[0])
  75     {
  76         FatalErrorIn
  77         (
  78             "solidBodyMotionFunctions::SKA::transformation()"
  79         )   << "current time (" << t
  80             << ") is less than the minimum in the data table ("
  81             << times_[0] << ')'
  82             << exit(FatalError);
  83     }
  84
  85     if (t > times_[times_.size()-1])
  86     {
  87         FatalErrorIn
  88         (
  89             "solidBodyMotionFunctions::SKA::transformation()"
  90         )   << "current time (" << t
  91             << ") is greater than the maximum in the data table ("
  92             << times_[times_.size()-1] << ')'
  93             << exit(FatalError);
  94     }
  95
  96     translationRotationVectors TRV = interpolateXY
  97     (
  98         t,
  99         times_,
 100         values_
 101     );
 102
 103     // Convert the rotational motion from deg to rad
 104     TRV[1] *= pi/180.0;
 105
 106     quaternion R(TRV[1].x(), TRV[1].y(), TRV[1].z());
 107     septernion TR(septernion(CofG_ + TRV[0])*R*septernion(-CofG_));
 108
 109     Info<< "solidBodyMotionFunctions::SKA::transformation(): "
 110         << "Time = " << t << " transformation: " << TR << endl;
 111
 112     return TR;
 113 }
 114
 115
 116 bool Foam::solidBodyMotionFunctions::SKA::read(const dictionary& SBMFCoeffs)
 117 {
 118     solidBodyMotionFunction::read(SBMFCoeffs);
 119
 120     // If the timeDataFileName has changed read the file
 121
 122     fileName newTimeDataFileName(SBMFCoeffs_.lookup("timeDataFileName"));
 123
 124     if (newTimeDataFileName != timeDataFileName_)
 125     {
 126         timeDataFileName_ = newTimeDataFileName;
 127
 128         IFstream dataStream(timeDataFileName_);
 129
 130         if (dataStream.good())
 131         {
 132             List<Tuple2<scalar, translationRotationVectors> > timeValues
 133             (
 134                 dataStream
 135             );
 136
 137             times_.setSize(timeValues.size());
 138             values_.setSize(timeValues.size());
 139
 140             forAll(timeValues, i)
 141             {
 142                 times_[i] = timeValues[i].first();
 143                 values_[i] = timeValues[i].second();
 144             }
 145         }
 146         else
 147         {
 148             FatalErrorIn
 149             (
 150                 "solidBodyMotionFunctions::SKA::read(const dictionary&)"
 151             )   << "Cannot open time data file " << timeDataFileName_
 152                 << exit(FatalError);
 153         }
 154     }
 155
 156     SBMFCoeffs_.lookup("CofG") >> CofG_;
 157
 158     return true;
 159 }
 160
 161
 162 // ************************************************************************* //