CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Pre-Processing

Initial torque for 6DoF simulation_OpenFOAM v1706

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 19, 2017, 03:04
Default Initial torque for 6DoF simulation_OpenFOAM v1706
  #1
New Member
 
Laurent Fieschi
Join Date: Aug 2017
Posts: 6
Rep Power: 7
Fieschi is on a distinguished road
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

Last edited by Fieschi; September 19, 2017 at 05:48.
Fieschi is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Segmentation fault when using reactingFOAM for Fluids Tommy Floessner OpenFOAM Running, Solving & CFD 4 April 22, 2018 12:30
chtMultiRegionSimpleFoam turbulent case Aditya Patil OpenFOAM Running, Solving & CFD 6 April 24, 2017 22:13
simpleFoam error - "Floating point exception" mbcx4jc2 OpenFOAM Running, Solving & CFD 12 August 4, 2015 02:20
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 06:20
calculation stops after few time steps sivakumar OpenFOAM Running, Solving & CFD 7 March 17, 2013 06:37


All times are GMT -4. The time now is 14:37.