[batuning]

Dear friends ,


can you run my code .i want simply rotate a turbine with wind velocity .but my code no compile and error , my fluent and visual studio works and no problem for compiling codes .just this code no works .can you find a problem ?


my other question is how can achieve ab omega (w) angular velocity every time?



#include "udf.h"
#include "math.h"

DEFINE_SDOF_PROPERTIES(rotor, prop, dt, time, dtime)
{
prop[SDOF_MASS] = 0.218786; /*Mass of the rigid body in [kg] */
prop[SDOF_IZZ] = 0.000161; /*Mass moment of inertia about Z axis [kg/m^2]*/

/* Translational motion setting, use TRUE and FALSE as applicable */
prop[SDOF_ZERO_TRANS_X] = TRUE; /*Translation allowed in global X-Direction? */
prop[SDOF_ZERO_TRANS_Y] = TRUE; /*Translation allowed in global Y-Direction? */
prop[SDOF_ZERO_TRANS_Z] = TRUE; /*Translation allowed in global Z-Direction? */

/* Rotational motion setting, use TRUE and FALSE as applicable*/
prop[SDOF_ZERO_ROT_X] = TRUE; /*Rotation allowed about global X-Axis? */
prop[SDOF_ZERO_ROT_Y] = TRUE; /*Rotation allowed about global Y-Axis? */
prop[SDOF_ZERO_ROT_Z] = FALSE; /*Rotation allowed about global Z-Axis? */


/* DT_THETA(dt): orientation of body-fixed axis vector, DT_CG(dt): center of gravity vector, DT_VEL_CG(dt): cg velocity vector, DT_OMEGA_CG(t): angular velocity vector */
prop[SDOF_LOAD_M_Z] = -9.81 * sin(DT_THETA(dt)[2]) ;
Message("\n 2D: updated 6DOF properties DT_THETA_Z: %e, Mz: %e, Mass: %e \n",
DT_THETA(dt)[2], prop[SDOF_MASS]);
}

[AlexanderZ]

code is correct except last line

Message("\n 2D: updated 6DOF properties DT_THETA_Z: %e, Mz: %e, Mass: %e \n",DT_THETA(dt)[2], prop[SDOF_LOAD_M_Z], prop[SDOF_MASS]);

however, it should not be a problem,
try it anyway

[batuning]

thanks .and can you tell me how i use (extract) omega (angular velocity ) for this code?
udf.c
i want move rotor and gate as wind power not as omega .


i should change omega (w) .and try put 6dof angular velocity