Commit df3a560d authored by mattijs's avatar mattijs
Browse files

ENH: surfaceInertia: print info for use in RBD. See also #1014.

parent 0b95bb00
......@@ -3,7 +3,7 @@
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\/ M anipulation | Copyright (C) 2015 OpenCFD Ltd.
\\/ M anipulation | Copyright (C) 2015-2018 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
......@@ -384,6 +384,26 @@ int main(int argc, char *argv[])
<< " momentOfInertia " << eVal << token::END_STATEMENT << nl
<< " orientation " << eVec.T() << token::END_STATEMENT
<< endl;
// Get inertia as a symmetric tensor
const symmTensor st
(
eVal.xx(),
eVal.xy(),
eVal.xz(),
eVal.yy(),
eVal.yz(),
eVal.zz()
);
Info<< nl
<< "Entries for rigidBodyDynamics dynamicMeshDict:"
<< nl
<< " mass " << m << token::END_STATEMENT << nl
<< " centreOfMass " << cM << token::END_STATEMENT << nl
<< " inertia " << st << token::END_STATEMENT << nl
<< endl;
}
if (calcAroundRefPt)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment