#include "udf.h"

?#include "math.h"

?DEFINE_SDOF_PROPERTIES(Moving_wall, prop, dt, time, dtime)

{

prop[SDOF_MASS] = 1.0;

?prop[SDOF_IZZ] = 0.1;

?prop[SDOF_ZERO_TRANS_X] = TRUE;

prop[SDOF_ZERO_TRANS_Y] = TRUE;

prop[SDOF_ZERO_TRANS_Z] = TRUE;

?prop[SDOF_ZERO_ROT_X] = TRUE;

prop[SDOF_ZERO_ROT_Y] = TRUE;

?prop[SDOF_ZERO_ROT_Z] = FALSE;

real ang = DT_THETA(dt)[2];

?real ang_vel = DT_OMEGA_CG(dt)[2];

?real moment_z = 0;

?if ( ang GOET 4.0 * 3.1415 / 9.0 - 3.1415 / 2.0)

{

moment_z = 0;

?}

else if ( ang LOET 3.1415 / 9.0 - 3.1415 / 2.0)

{

moment_z = 0.8;

?}

?else if ( (ang LT 4.0 * 3.1415 / 9.0 - 3.1415 / 2.0) && (ang GT 3.1415 / 9.0 - 3.1415 / 2.0) && (ang_vel GT 0) )

{

moment_z = 0.8;

?}

else

?{

moment_z = 0.0;

}

?prop[SDOF_LOAD_M_Z] = moment_z;

?}