CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Defining centre of mass/rotation in rigidBodyMotion solver

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 9, 2017, 10:51
Default Defining centre of mass/rotation in rigidBodyMotion solver
  #1
New Member
 
Owen Jones
Join Date: Aug 2017
Location: Martinique
Posts: 2
Rep Power: 0
tok3rat0r is on a distinguished road
I'm currently running a simple model of a boat hull moving at constant velocity on an initially flat free surface with interDyMFoam.

This is my dynamicMeshDict for the sixDoFRigidBodyMotion solver:
Code:
dynamicFvMesh   dynamicMotionSolverFvMesh;

motionSolverLibs    ("libsixDoFRigidBodyMotion.so");

solver    sixDoFRigidBodyMotion;

sixDoFRigidBodyMotionCoeffs {
patches         (boat);
innerDistance   0.2;
outerDistance   1.2;

centreOfMass    (-3.2 0 0.7);
mass            11000;
momentOfInertia (4413 22295 20859);
rhoInf          1;
report          on;

//value uniform (0, 0, 0);

accelerationRelaxation 0.05;
accelerationDamping 0.9;

solver
{
    type Newmark;
}
}
Using this, the centre of mass and rotation are coincident as expected, and everything seems fine (I've no idea what the value uniform (0, 0, 0) field is for, and commenting it out doesn't seem to change anything). My question is, how do I translate this to a suitable dictionary for the 'new' solver rigidBodyMotion? My best guess (based on the DTCHull tutorial) is as follows, but I'm not sure about the centreOfMass and transform fields: does 'transform' affect both the centre of mass and the centre of rotation, moving the centre of mass from (0,0,0) to (-3.2,0,0.7) in this case, or does it only affect the centre of rotation (relevant to the inertia tensor)? I've also no idea what the joint structure is there for.
Code:
dynamicFvMesh   dynamicMotionSolverFvMesh;

motionSolverLibs   ("librigidBodyMeshMotion.so");

solver    rigidBodyMotion;

rigidBodyMotionCoeffs {
report          on;

solver
{
    type Newmark;
}

accelerationRelaxation 0.05;
accelerationDamping 0.9;

bodies
{
    boat
    {
        type            rigidBody;
        parent          root;

        centreOfMass    (0 0 0);
        mass            11000;
        inertia (4413 0 0 22295 0 20859);
        transform       (1 0 0 0 1 0 0 0 1) (-3.2 0 0.7);

        joint
        {
            type    composite;
            joints
            (
                {
                    type Pz;
                }
                {
                    type Ry;
                }
            );
        }

        patches         (boat);
        innerDistance   0.2;
        outerDistance   1.2;
    }
}
}
tok3rat0r is offline   Reply With Quote

Old   November 26, 2018, 10:44
Default its transformational constraint in a simpler form
  #2
New Member
 
Vignesh S P
Join Date: Jul 2018
Location: Coimbatore, Tamil Nadu, India
Posts: 17
Rep Power: 8
VigneshSP is on a distinguished road
i think thats a translational constraint in a combined form if u run the sixDoFRigidBodyMotion solver we have to specify the constraint under constraints funtion like highlighted below


Code:
motionSolverLibs    ("libsixDoFRigidBodyMotion.so");

dynamicFvMesh       dynamicOversetFvMesh;

dynamicOversetFvMeshCoeffs
{
}

solver          sixDoFRigidBodyMotion;

sixDoFRigidBodyMotionCoeffs
{

    patches         (floatingObject);
    innerDistance   100.0;
    outerDistance   101.0;

    centreOfMass    (0 0 0);

    // Cuboid dimensions
    

    // Density of the solid
    rhoSolid        780;

    // Cuboid mass
    mass            0.40872;

    // Cuboid moment of inertia about the centre of mass
    momentOfInertia (1.21327e-05 1.21327e-05 1.21327e-05);
    velocity (0 0 0);
    report          on;
    accelerationRelaxation 0.6;

    solver
    {
        type Newmark;
    }

    constraints
    {

        
        /*fixedPoint
        {
             sixDoFRigidBodyMotionConstraint point;
             centreOfRotation (0.5 0.45 0.1);
        }*/
        

        /*fixedLine
        {
            sixDoFRigidBodyMotionConstraint line;
            //centreOfRotation (0.5 0.45 0.1);
            direction (0 0 1);
        }*/

        
        /*fixedAxis
        {
            sixDoFRigidBodyMotionConstraint axis;
            axis (1 1 1);
        }*/
    /*arrestRotation
     {
           sixDoFRigidBodyMotionConstraint   orientation;
           
     }*/

        
    }
}

while running the above the output log generated



Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1806                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/


// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Selecting dynamicFvMesh dynamicOversetFvMesh
Selecting motion solver: sixDoFRigidBodyMotion
Selecting sixDoFSolver Newmark
Translational constraint tensor (1 0 0 0 1 0 0 0 1)
Rotational constraint tensor (1 0 0 0 1 0 0 0 1)

PIMPLE: no residual control data found. Calculations will employ 2 corrector loops

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Creating cellMask field to block out hole cells
VigneshSP is offline   Reply With Quote

Reply

Tags
centre of mass, dynamic mesh, interdymfoam, rigid body motion

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
fluent divergence for no reason sufjanst FLUENT 2 March 23, 2016 17:08
3d vof Smaras FLUENT 2 February 19, 2013 07:58
CFX 5.5 Roued CFX 1 October 2, 2001 17:49
Setting a B.C using UserFortran in 4.3 tokai CFX 10 July 17, 2001 17:25
Error during Solver cfd guy CFX 4 May 8, 2001 07:04


All times are GMT -4. The time now is 00:21.