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

transform in the new motion library rigidBodyMeshMotion

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 13, 2018, 23:46
Default transform in the new motion library rigidBodyMeshMotion
  #1
New Member
 
Jiadong Wang
Join Date: Mar 2018
Location: Australia
Posts: 23
Rep Power: 2
jiadongw is on a distinguished road
Hello Foamers,

Recently, Iím trying to conduct the tutorial wingMotion with the new motion library rigidBodyMeshMotion. This tutorial originally employs the sixDoFRigidBodyMotion.

Unfortunately, the new case with rigidBodyMeshMotion has some problems. One of quite weird things is the moment of linearAxialAngularSpring has a same direction with angle, which means the axis spring accelerate the rotation instead of restraining it.
Also, Iím quite confused about the parameter Ďtransformí in dynamicMeshDict of rigidBodyMeshMotion.

I have attached the modified dynamicMeshDict (rigidBodyMeshMotion) and original dynamicMeshDict (sixDoFRigidBodyMotion) to this post.
Does anyone have the experience with this new motion library rigidBodyMeshMotion?

Any suggestions will be appreciated.

Thank you.

dynamicMeshDict of rigidBodyMeshMotion:
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5                                     |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dynamicFvMesh      dynamicMotionSolverFvMesh;

motionSolverLibs ("libRY05rigidBodyMeshMotion.so");

motionSolver    rigidBodyMotion;

report          on;

solver
{
    type symplectic;
}

g               (0 -9.81 0);
rho             rhoInf;
rhoInf          1;

bodies
{
   wing
  {

    type            rigidBody;
    parent          root;
    centreOfMass    (0.4974612746 -0.01671895744 0.125);
    mass            22.9;
    inertia         (1.958864357 0 0 3.920839234 0 2.057121362);
    transform       (1 0 0 0 1 0 0 0 1) (0.25 0.007 0.125);
   
    joint
    {
       type   composite;
       joints
       (
            {
                 type     Py;
            }                      
            {
                 type     Rz;
            }
       );
    }

        patches         (wing);
        innerDistance   0.3;
        outerDistance   1;

  }
}

restraints
{
    verticalSpring
    {
        type            linearSpring;
        body            wing;
        anchor          (0.25 0.007 0.125);
        refAttachmentPt (0.25 0.007 0.125);
        stiffness       4000;
        damping         2;
        restLength      0;
    }

    axialSpring
    {
        type            linearAxialAngularSpring;
        body            wing;
        axis            (0 0 1);
        referenceOrientation (1 0 0 0 1 0 0 0 1);
        stiffness       700;
        damping         0.5;
    }
}


// ************************************************************************* //
dynamicMeshDict of sixDoFRigidBodyMotion:
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5                                     |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dynamicFvMesh      dynamicMotionSolverFvMesh;

motionSolverLibs ("libtest00.so");

motionSolver    sixDoFRigidBodyMotion;

patches         (wing);
innerDistance   0.3;
outerDistance   1;

mass            22.9;
centreOfMass    (0.4974612746 -0.01671895744 0.125);
momentOfInertia (1.958864357 3.920839234 2.057121362);
orientation
(
    0.9953705935 0.09611129781 0
    -0.09611129781 0.9953705935 0
    0 0 1
);
angularMomentum (0 0 -2);
g               (0 -9.81 0);
rho             rhoInf;
rhoInf          1;
report          on;

solver
{
    type symplectic;
}

constraints
{
    yLine
    {
        sixDoFRigidBodyMotionConstraint line;
        centreOfRotation    (0.25 0.007 0.125);
        direction           (0 1 0);
    }

    zAxis
    {
        sixDoFRigidBodyMotionConstraint axis;
        axis                (0 0 1);
    }
}

restraints
{
    verticalSpring
    {
        sixDoFRigidBodyMotionRestraint linearSpring;

        anchor          (0.25 0.007 0.125);
        refAttachmentPt (0.25 0.007 0.125);
        stiffness       4000;
        damping         2;
        restLength      0;
    }

    axialSpring
    {
        sixDoFRigidBodyMotionRestraint linearAxialAngularSpring;

        axis            (0 0 1);
        stiffness       700;
        damping         0.5;
        referenceOrientation $orientation;
    }
}


// ************************************************************************* //
jiadongw is offline   Reply With Quote

Old   July 13, 2018, 23:51
Default
  #2
New Member
 
Jiadong Wang
Join Date: Mar 2018
Location: Australia
Posts: 23
Rep Power: 2
jiadongw is on a distinguished road
Just ignore the name of library in both dynamicMeshDict. That’s the recompiled version without any changes except I add some output orders using Info.
jiadongw is offline   Reply With Quote

Old   July 15, 2018, 21:43
Default
  #3
New Member
 
Jiadong Wang
Join Date: Mar 2018
Location: Australia
Posts: 23
Rep Power: 2
jiadongw is on a distinguished road
Any ideas?
jiadongw is offline   Reply With Quote

Reply

Tags
rigid body motion, transform, wing motion

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
How to use "translation" in solidBodyMotionFunction in OpenFOAM rupesh_w OpenFOAM Running, Solving & CFD 5 August 16, 2016 04:27
sliding mesh problem in CFX Saima CFX 45 September 22, 2015 10:53
External management of solid motion (using dynamicmesh) maxou1993 Main CFD Forum 0 July 28, 2015 11:37
decomposePar is missing a library whk1992 OpenFOAM Pre-Processing 8 March 7, 2015 07:53
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 06:20


All times are GMT -4. The time now is 11:02.