CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Pre-Processing (https://www.cfd-online.com/Forums/openfoam-pre-processing/)
-   -   Initial torque for 6DoF simulation_OpenFOAM v1706 (https://www.cfd-online.com/Forums/openfoam-pre-processing/193192-initial-torque-6dof-simulation_openfoam-v1706.html)

Fieschi September 19, 2017 03:04

Initial torque for 6DoF simulation_OpenFOAM v1706
 
1 Attachment(s)
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.


Attachment 58481
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


All times are GMT -4. The time now is 03:45.