DEFINE_CG_MOTION宏【注释版】

线速度是通过物体上的x方向的力平衡达到的。表达形式为：

此处v为速度，F为外力，m为质量。使用显示欧拉格式表达t时刻速度为：

源代码：

#include "udf.h"

 

static real v_prev = 0.0;     

DEFINE_CG_MOTION(moving_body, dt, vel, omega, time, dtime)

{

    Thread *t;

    face_t f;

    real NV_VEC(A);        //定义一个向量变量A

    real force, dv;

 

    NV_S(vel, =, 0.0);        //对变量vel赋值为0.0

    NV_S(omega, =, 0.0);    //对变量omega赋值为0.0

    if(!Data_Valid_P())

        return;

 

    t = DT_THREAD(dt);        //DT_THREAD 就是获取网格运动所在的 thread 指针

 

    force = 0.0;

    begin_f_loop(f, t)

    {

         F_AREA(A, f, t);    //获取面上单元的面积，并且存储在向量A当中

         force = force + F_P(f, t) * A[0];        //循环累加获取x方向的力的大小

    }

    end_f_loop(f, t)

    

    dv = dtime * force / 100.0;    

    v_prev = v_prev + dv;    //原来的速度加上速度的增量，得到现在的速度

    vel[0] = v_prev;    //将现在的速度返回给动网格宏函数

}