cuboidI.H
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration | Website: https://openfoam.org
5  \\ / A nd | Copyright (C) 2016-2019 OpenFOAM Foundation
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
13  the Free Software Foundation, either version 3 of the License, or
14  (at your 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, see <http://www.gnu.org/licenses/>.
23 
24 \*---------------------------------------------------------------------------*/
25 
26 // * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
27 
28 inline Foam::symmTensor Foam::RBD::cuboid::I
29 (
30  const scalar m,
31  const vector& L
32 ) const
33 {
34  const scalar mBy12 = m/12.0;
35  const scalar mSqrLx = mBy12*sqr(L.x());
36  const scalar mSqrLy = mBy12*sqr(L.y());
37  const scalar mSqrLz = mBy12*sqr(L.z());
38 
39  return symmTensor
40  (
41  mSqrLy + mSqrLz, 0, 0,
42  mSqrLx + mSqrLz, 0,
43  mSqrLx + mSqrLy
44  );
45 }
46 
47 
48 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
49 
51 (
52  const word& name,
53  const scalar m,
54  const vector& c,
55  const vector& L
56 )
57 :
58  rigidBody(name, m, c, I(m, L)),
59  L_(L)
60 {}
61 
62 
64 (
65  const word& name,
66  const dictionary& dict
67 )
68 :
69  rigidBody(name, rigidBodyInertia()),
70  L_(dict.lookup("L"))
71 {
72  const scalar m(dict.lookup<scalar>("mass"));
73  const vector c(dict.lookup("centreOfMass"));
74  rigidBodyInertia::operator=(rigidBodyInertia(m, c, I(m, L_)));
75 }
76 
77 
78 // * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * //
79 
80 inline const Foam::vector& Foam::RBD::cuboid::L() const
81 {
82  return L_;
83 }
84 
85 
86 // ************************************************************************* //
A list of keyword definitions, which are a keyword followed by any number of values (e...
Definition: dictionary.H:158
dimensionedSymmTensor sqr(const dimensionedVector &dv)
scalar m() const
Return the mass of the rigid-body.
Vector< scalar > vector
A scalar version of the templated Vector.
Definition: vector.H:49
rigidBodyInertia()
Null constructor, initializes to zero.
const vector & c() const
Return the centre of mass of the rigid-body.
SymmTensor< scalar > symmTensor
SymmTensor of scalars.
Definition: symmTensor.H:48
A class for handling words, derived from string.
Definition: word.H:59
rigidBody(const word &name, const scalar &m, const vector &c, const symmTensor &Ic)
Construct from mass, centre of mass and moment of inertia tensor.
Definition: rigidBodyI.H:29
const vector & L() const
Return the lengths of the sides of the cuboid.
Definition: cuboidI.H:80
cuboid(const word &name, const scalar m, const vector &c, const vector &L)
Construct from name, mass, centre of mass and lengths of the sides.
Definition: cuboidI.H:51
ITstream & lookup(const word &, bool recursive=false, bool patternMatch=true) const
Find and return an entry data stream.
Definition: dictionary.C:812