![]() |
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 23:22. |