CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   Fluent UDF and Scheme Programming (http://www.cfd-online.com/Forums/fluent-udf/)
-   -   UDF code (http://www.cfd-online.com/Forums/fluent-udf/104800-udf-code.html)

yahya_azizi July 16, 2012 10:46

UDF code
 
I'm simulating a flapping wing aircraft and have designed a rectangular wing to model that. I have to use dynamic motion to model the flapping motion. However, I am not familiar with dynamic mesh. What I know is that I have to use a UDF code to define the motion. I am using the following UDF code which has been posted by a member on the forum to simulate the motion but everytime I compile the code I get an error. has anyone faced this problem before?


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

DEFINE_CG_MOTION(asymflap, dt, v_cg, omega, time, dtime)

{
Thread *t;

NV_S(v_cg, =, 0.0);
NV_S(omega, =, 0.0);

if (!Data_Valid_P())
{
Message("Wrong\n");
return;
}

t=DT_THREAD((Dynamic_Thread*)dt);

frequency = 3.0;

pi = 3.141592654;

/* define motion variables */

Ar = 30.0; /*roll amplitude*/
Ap = 5.0; /*pitch amplitude*/
Ay = 0.0; /*yaw amplitude*/

Aroll = Ar*pi/180.0; /*conversion to radians*/
Apitch = Ap*pi/180.0; /*conversion to radians*/

w=2.0*pi*frequency; /*omega (radians)*/
T=1.0/frequency;

roll = Aroll*sin(w*time + pi/2);
droll = w*Aroll*cos(w*time + pi/2);

pitch = Apitch*sin(w*time);
dpitch = w*Apitch*cos(w*time);

/*euler angles initial frame mod 1*/

omega[0] = droll;
omega[1] = dpitch*1.0*cos(roll);
omega[2] = dpitch*1.0*sin(roll);

Message("time:%f omegax:%f omegay:%f omegaz:%f \n", time, omega[0], omega[1],omega[2]);

flotus1 July 16, 2012 10:50

What kind of error do you get when trying to compile the UDF?

Did you try interpreting it instead?


All times are GMT -4. The time now is 02:19.