[Atze]

Hope it can help someone

/* in this example a cilynder with CG=(0,0) in loaded by its own wake. A spring of K=15 keeps it in position*/

#include "udf.h"
#include "dynamesh_tools.h"

DEFINE_CG_MOTION(oscillazione,dt,vel,omega,time,dt ime)
{

real x[ND_ND];
Thread *t;
face_t f;
real A[ND_ND];
real Ft,Fel, dv;
real Keq= 15;
real m = 0.25;
real s;
real force=0;
real v_prev=0;

t = DT_THREAD(dt);

s=DT_CG(dt)[1];

begin_f_loop(f,t)
{
F_AREA(A,f,t);
force += F_P(f,t) * A[1];
}
end_f_loop(f,t)

Fel=-Keq*s;

Ft=force+Fel;
dv=(Ft/m)*dtime;

vel[1] = DT_VEL_CG(dt)[1]+dv;

}