[rubenparedes]

Hi,

I am simulating a marine turbine rotating due to incoming flow. I am trying to implement a break that applies a constant torque on a shaft using codeStream as a rigidBodyMotion restraint within a dynamicMeshDict. The sign of the torque depends only on the direction of rotation.

I tried to follow some codeStream examples without success.

restraints
{
    breakTorque
    {
        type        externalForce;
        
        body        Turbine;
        location    (1.0 0 0.11432977);
        
        force         #codeStream
        {
            codeInclude
            #{
//                #include "rigidBodyState.H"
                #include "rigidBodyMotion.H"
//                #include "rigidBodyMeshMotion.H"
            #};
            
            codeOptions
            #{
                -I$(LIB_SRC)/rigidBodyDynamics/lnInclude \
//                -I$(LIB_SRC)/rigidBodyMeshMotion/lnInclude
            #};

            codeLibs
            #{

                -lrigidBodyDynamics
 //               -lrigidBodyMeshMotion
            #};

            code
            #{                
               scalar omega = vCofR.w(v(0,Zero))[2]
               torque = 16.0 
                
                if (omega < 0)
                {
                    os  << (0 torque 0.0)
                }
                else
                {	
                    os  << (0 -torque 0.0)
                }
             #};
        };        
    }
}

Any suggestion will be appreciated.

Ruben