mirror of
https://develop.openfoam.com/Development/openfoam.git
synced 2025-11-28 03:28:01 +00:00
rigidBodyDynamics/bodies/cuboid: New body shape
Calculate the inertia from the lengths of the sides
This commit is contained in:
Henry Weller
76
src/rigidBodyDynamics/bodies/cuboid/cuboid.C
Normal file
76
src/rigidBodyDynamics/bodies/cuboid/cuboid.C
Normal file
@ -0,0 +1,76 @@
|
||||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "cuboid.H"
|
||||
#include "addToRunTimeSelectionTable.H"
|
||||
|
||||
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
|
||||
|
||||
namespace Foam
|
||||
{
|
||||
namespace RBD
|
||||
{
|
||||
defineTypeNameAndDebug(cuboid, 0);
|
||||
|
||||
addToRunTimeSelectionTable
|
||||
(
|
||||
rigidBody,
|
||||
cuboid,
|
||||
dictionary
|
||||
);
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
||||
|
||||
Foam::autoPtr<Foam::RBD::rigidBody> Foam::RBD::cuboid::clone() const
|
||||
{
|
||||
return autoPtr<rigidBody>(new cuboid(*this));
|
||||
}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
|
||||
|
||||
Foam::RBD::cuboid::~cuboid()
|
||||
{}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * //
|
||||
|
||||
void Foam::RBD::cuboid::write(Ostream& os) const
|
||||
{
|
||||
os.writeKeyword("type")
|
||||
<< type() << token::END_STATEMENT << nl;
|
||||
|
||||
os.writeKeyword("mass")
|
||||
<< m() << token::END_STATEMENT << nl;
|
||||
|
||||
os.writeKeyword("L")
|
||||
<< L() << token::END_STATEMENT << nl;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
124
src/rigidBodyDynamics/bodies/cuboid/cuboid.H
Normal file
124
src/rigidBodyDynamics/bodies/cuboid/cuboid.H
Normal file
@ -0,0 +1,124 @@
|
||||
/*---------------------------------------------------------------------------*\3
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
Class
|
||||
Foam::cuboid
|
||||
|
||||
Description
|
||||
Specialization of rigidBody to construct a cuboid given the mass and
|
||||
lengths of the sides.
|
||||
|
||||
SourceFiles
|
||||
cuboidI.H
|
||||
cuboid.C
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#ifndef RBD_cuboid_H
|
||||
#define RBD_cuboid_H
|
||||
|
||||
#include "rigidBody.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
namespace Foam
|
||||
{
|
||||
namespace RBD
|
||||
{
|
||||
|
||||
/*---------------------------------------------------------------------------*\
|
||||
Class cuboid Declaration
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
class cuboid
|
||||
:
|
||||
public rigidBody
|
||||
{
|
||||
// Private member data
|
||||
|
||||
//- Lengths of the sides
|
||||
vector L_;
|
||||
|
||||
|
||||
// Private member functions
|
||||
|
||||
// Calculate and return the inertia tensor
|
||||
inline symmTensor I(const scalar m, const vector& L) const;
|
||||
|
||||
|
||||
public:
|
||||
|
||||
//- Runtime type information
|
||||
TypeName("cuboid");
|
||||
|
||||
|
||||
// Constructors
|
||||
|
||||
//- Construct from name, mass, centre of mass and lengths of the sides.
|
||||
inline cuboid
|
||||
(
|
||||
const word& name,
|
||||
const scalar m,
|
||||
const vector& c,
|
||||
const vector& L
|
||||
);
|
||||
|
||||
//- Construct from dictionary
|
||||
inline cuboid
|
||||
(
|
||||
const word& name,
|
||||
const dictionary& dict
|
||||
);
|
||||
|
||||
//- Return clone of this cuboid
|
||||
virtual autoPtr<rigidBody> clone() const;
|
||||
|
||||
|
||||
//- Destructor
|
||||
virtual ~cuboid();
|
||||
|
||||
|
||||
// Member Functions
|
||||
|
||||
//- Return the lengths of the sides of the cuboid
|
||||
inline const vector& L() const;
|
||||
|
||||
//- Write
|
||||
virtual void write(Ostream&) const;
|
||||
};
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
} // End namespace RBD
|
||||
} // End namespace Foam
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#include "cuboidI.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#endif
|
||||
|
||||
// ************************************************************************* //
|
||||
86
src/rigidBodyDynamics/bodies/cuboid/cuboidI.H
Normal file
86
src/rigidBodyDynamics/bodies/cuboid/cuboidI.H
Normal file
@ -0,0 +1,86 @@
|
||||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
|
||||
|
||||
inline Foam::symmTensor Foam::RBD::cuboid::I
|
||||
(
|
||||
const scalar m,
|
||||
const vector& L
|
||||
) const
|
||||
{
|
||||
const scalar mBy12 = m/12.0;
|
||||
const scalar mSqrLx = mBy12*sqr(L.x());
|
||||
const scalar mSqrLy = mBy12*sqr(L.y());
|
||||
const scalar mSqrLz = mBy12*sqr(L.z());
|
||||
|
||||
return symmTensor
|
||||
(
|
||||
mSqrLy + mSqrLz, 0, 0,
|
||||
mSqrLx + mSqrLz, 0,
|
||||
mSqrLx + mSqrLy
|
||||
);
|
||||
}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
||||
|
||||
inline Foam::RBD::cuboid::cuboid
|
||||
(
|
||||
const word& name,
|
||||
const scalar m,
|
||||
const vector& c,
|
||||
const vector& L
|
||||
)
|
||||
:
|
||||
rigidBody(name, m, c, I(m, L)),
|
||||
L_(L)
|
||||
{}
|
||||
|
||||
|
||||
inline Foam::RBD::cuboid::cuboid
|
||||
(
|
||||
const word& name,
|
||||
const dictionary& dict
|
||||
)
|
||||
:
|
||||
rigidBody(name, rigidBodyInertia()),
|
||||
L_(dict.lookup("L"))
|
||||
{
|
||||
const scalar m(readScalar(dict.lookup("mass")));
|
||||
const vector c(dict.lookup("centreOfMass"));
|
||||
rigidBodyInertia::operator=(rigidBodyInertia(m, c, I(m, L_)));
|
||||
}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * //
|
||||
|
||||
inline const Foam::vector& Foam::RBD::cuboid::L() const
|
||||
{
|
||||
return L_;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
Reference in New Issue
Block a user