src/dynamicMesh/meshMotion/solidBodyMotion/linearOscillation/linearOscillation.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 \*---------------------------------------------------------------------------*/
  25
  26 #include "linearOscillation.H"
  27 #include "addToRunTimeSelectionTable.H"
  28 #include "mathematicalConstants.H"
  29
  30 using namespace Foam::mathematicalConstant;
  31
  32 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  33
  34 namespace Foam
  35 {
  36 namespace solidBodyMotionFunctions
  37 {
  38     defineTypeNameAndDebug(linearOscillation, 0);
  39     addToRunTimeSelectionTable
  40     (
  41         solidBodyMotionFunction,
  42         linearOscillation,
  43         dictionary
  44     );
  45 };
  46 };
  47
  48
  49 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  50
  51 Foam::vector
  52 Foam::solidBodyMotionFunctions::linearOscillation::calcPosition
  53 (
  54     const scalar t
  55 ) const
  56 {
  57     return amplitude_*sin(2*pi*t/period_);
  58 }
  59
  60
  61 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  62
  63 Foam::solidBodyMotionFunctions::linearOscillation::linearOscillation
  64 (
  65     const dictionary& SBMFCoeffs,
  66     const Time& runTime
  67 )
  68 :
  69     solidBodyMotionFunction(SBMFCoeffs, runTime),
  70     amplitude_(SBMFCoeffs_.lookup("amplitude")),
  71     period_(readScalar(SBMFCoeffs_.lookup("period")))
  72 {}
  73
  74
  75 // * * * * * * * * * * * * * * * * Destructors * * * * * * * * * * * * * * * //
  76
  77 Foam::solidBodyMotionFunctions::linearOscillation::~linearOscillation()
  78 {}
  79
  80
  81 // * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * * //
  82
  83 Foam::septernion
  84 Foam::solidBodyMotionFunctions::linearOscillation::transformation() const
  85 {
  86     scalar t = time_.value();
  87
  88     septernion TR(calcPosition(t), quaternion::I);
  89
  90     Info<< "solidBodyMotionFunctions::linearOscillation::transformation(): "
  91         << "Time = " << t << " transformation: " << TR << endl;
  92
  93     return TR;
  94 }
  95
  96
  97 Foam::septernion
  98 Foam::solidBodyMotionFunctions::linearOscillation::velocity() const
  99 {
 100     scalar t = time_.value();
 101     scalar dt = time_.deltaT().value();
 102
 103     septernion TV
 104     (
 105         (calcPosition(t) - calcPosition(t - dt))/dt,
 106         quaternion::zero
 107     );
 108
 109     return TV;
 110 }
 111
 112
 113 bool Foam::solidBodyMotionFunctions::linearOscillation::read
 114 (
 115     const dictionary& SBMFCoeffs
 116 )
 117 {
 118     solidBodyMotionFunction::read(SBMFCoeffs);
 119
 120     SBMFCoeffs_.lookup("amplitude") >> amplitude_;
 121     SBMFCoeffs_.lookup("amplitude") >> period_;
 122
 123     return true;
 124 }
 125
 126
 127 // ************************************************************************* //