[Fieschi]

Hello,

I am trying to model a 2D simulation of an object with 6dof model in a fixed velocity domain in x-direction.

My object is a plane with a -5 deg origin angle on x axis.


Domaine_plaque2D.png
The rectangular bloc is my background (fixed block) and my object is inside a moving cylinder wich is my overset.

I would like to give it an initial torque on z axis and then let gravity, drag and so on act on my object.

This is my dynamicMeshDict file :

FoamFile
{
version 2.0;
format ascii;
class dictionary;
object dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


motionSolverLibs ("libsixDoFRigidBodyMotion.so");


dynamicFvMesh dynamicOversetFvMesh;


dynamicOversetFvMeshCoeffs
{
}


solver sixDoFRigidBodyMotion;


sixDoFRigidBodyMotionCoeffs
{
patches (object);
innerDistance 100.0;
outerDistance 101.0;


centreOfMass (0.0 0.0 0.125);


Lx 1.0;
Ly 0.1;
Lz 0.05;


rhoSolid 1000;


mass #calc "$rhoSolid*$Lx*$Ly*$Lz";
momentOfInertia #codeStream
{
codeInclude
#{
#include "diagTensor.H"
#};


code
#{
scalar sqrLx = sqr($Lx);
scalar sqrLy = sqr($Ly);
scalar sqrLz = sqr($Lz);
os <<
$mass
*diagTensor(sqrLy + sqrLz, sqrLx + sqrLz, sqrLx + sqrLy)/12.0;
#};
};


angularMomentum (0 0 -1);


rho rhoInf;
rhoInf 1;
report on;
reportToFile on;


accelerationRelaxation 0.1;


solver
{
type Newmark;
}


constraints
{

fixedPoint
{
sixDoFRigidBodyMotionConstraint point;
centreOfRotation (0.0 0.0 0.125);
}
}
}


// ************************************************** *********************** //

I would like to know if « angularMomemtum » will act only at the initial state (and then let my object in a free motion) or if it will be constant during the run ?

Best regards