CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

6DOF Solver - Center of Mass and Center of Rotation

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

LinkBack Thread Tools Display Modes
Old   October 13, 2016, 14:51
Default 6DOF Solver - Center of Mass and Center of Rotation
New Member
Join Date: May 2016
Posts: 27
Rep Power: 3
Verse is on a distinguished road

I've tried defining both centreOfMass and centreOfRotation inside the dynamicMeshDict, as follows:

centreOfMass      (0.4 0 0.1)
centreOfRotation (-0.5 0 0.1)
After doing this, I recieve an error message as follows:

#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86-64-linux-gnu/"
I see that from the sixDoFRigidBodyMotionState.C code that the solver will default to centreOfMass if the centreOfRotation isn't present:

  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
    \\  /    A nd           | Copyright (C) 2011-2014 OpenFOAM Foundation
     \\/     M anipulation  |
    This file is part of OpenFOAM.

    OpenFOAM is free software: you can redistribute it and/or modify it
    under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
    for more details.

    You should have received a copy of the GNU General Public License
    along with OpenFOAM.  If not, see <>.


#include "sixDoFRigidBodyMotionState.H"

// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //


    const dictionary& dict
            dict.lookupOrDefault("centreOfMass", vector::zero)
    Q_(dict.lookupOrDefault("orientation", tensor::I)),
    v_(dict.lookupOrDefault("velocity", vector::zero)),
    a_(dict.lookupOrDefault("acceleration", vector::zero)),
    pi_(dict.lookupOrDefault("angularMomentum", vector::zero)),
    tau_(dict.lookupOrDefault("torque", vector::zero))
But surely there should be a way I could define both in the same test case? Including the effects of gravity is significantly important for my test case (but I also need rotation), so having these two variables as separate is very important.

Thank you.
Verse is offline   Reply With Quote


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
center of mass for naca airfoil frederic Felten Main CFD Forum 4 November 15, 2016 12:24
Initial Center of Mass king_411 STAR-CCM+ 0 May 20, 2016 20:41
Can you help me with a problem in ansys static structural solver? sourabh.porwal Structural Mechanics 0 March 27, 2016 17:07
6DOF Model Rocket Stability with Changing Mass LWhitson2 CFX 1 March 18, 2016 20:38
How to solve for components of Center of mass? 13msmemusman OpenFOAM Programming & Development 3 March 17, 2015 12:46

All times are GMT -4. The time now is 06:52.