rigidBodyModel: Removed g from the body acceleration

This commit is contained in:
Henry Weller
2023-09-27 21:26:05 +01:00
parent b7e88c09ac
commit bb127a4523
2 changed files with 13 additions and 7 deletions

View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org \\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2016-2022 OpenFOAM Foundation \\ / A nd | Copyright (C) 2016-2023 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -216,12 +216,15 @@ inline Foam::spatialVector Foam::RBD::rigidBodyModel::a
{ {
return return
( (
spatialTransform
( (
X0_[master(bodyID)].E().T(), spatialTransform
masterPoint(bodyID, p) (
X0_[master(bodyID)].E().T(),
masterPoint(bodyID, p)
)
& a_[master(bodyID)]
) )
& a_[master(bodyID)] + spatialVector(Zero, g_)
); );
} }

View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org \\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2016-2021 OpenFOAM Foundation \\ / A nd | Copyright (C) 2016-2023 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -171,12 +171,15 @@ void Foam::RBD::rigidBodyMotion::status(const label bodyID) const
{ {
const spatialTransform CofR(X0(bodyID)); const spatialTransform CofR(X0(bodyID));
const spatialVector vCofR(v(bodyID, Zero)); const spatialVector vCofR(v(bodyID, Zero));
const spatialVector aCofR(a(bodyID, Zero));
Info<< "Rigid-body motion of the " << name(bodyID) << nl Info<< "Rigid-body motion of the " << name(bodyID) << nl
<< " Centre of rotation: " << CofR.r() << nl << " Centre of rotation: " << CofR.r() << nl
<< " Orientation: " << CofR.E() << nl << " Orientation: " << CofR.E() << nl
<< " Linear velocity: " << vCofR.l() << nl << " Linear velocity: " << vCofR.l() << nl
<< " Angular velocity: " << vCofR.w() << " Angular velocity: " << vCofR.w() << nl
<< " Linear acceleration: " << aCofR.l() << nl
<< " Angular acceleration: " << aCofR.w()
<< endl; << endl;
} }