diff --git a/src/rigidBodyDynamics/rigidBodyModel/rigidBodyModelI.H b/src/rigidBodyDynamics/rigidBodyModel/rigidBodyModelI.H index 906790abb8..11e5eaf96e 100644 --- a/src/rigidBodyDynamics/rigidBodyModel/rigidBodyModelI.H +++ b/src/rigidBodyDynamics/rigidBodyModel/rigidBodyModelI.H @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org - \\ / A nd | Copyright (C) 2016-2022 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2016-2023 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -216,12 +216,15 @@ inline Foam::spatialVector Foam::RBD::rigidBodyModel::a { return ( - spatialTransform ( - X0_[master(bodyID)].E().T(), - masterPoint(bodyID, p) + spatialTransform + ( + X0_[master(bodyID)].E().T(), + masterPoint(bodyID, p) + ) + & a_[master(bodyID)] ) - & a_[master(bodyID)] + + spatialVector(Zero, g_) ); } diff --git a/src/rigidBodyDynamics/rigidBodyMotion/rigidBodyMotion.C b/src/rigidBodyDynamics/rigidBodyMotion/rigidBodyMotion.C index 90cbda02c7..b090cd31c8 100644 --- a/src/rigidBodyDynamics/rigidBodyMotion/rigidBodyMotion.C +++ b/src/rigidBodyDynamics/rigidBodyMotion/rigidBodyMotion.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org - \\ / A nd | Copyright (C) 2016-2021 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2016-2023 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -171,12 +171,15 @@ void Foam::RBD::rigidBodyMotion::status(const label bodyID) const { const spatialTransform CofR(X0(bodyID)); const spatialVector vCofR(v(bodyID, Zero)); + const spatialVector aCofR(a(bodyID, Zero)); Info<< "Rigid-body motion of the " << name(bodyID) << nl << " Centre of rotation: " << CofR.r() << nl << " Orientation: " << CofR.E() << 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; }