UDF添加磁力，可为什么没有效果

tchandsome

各位大侠，我通过UDF加入一体积力（磁力），编译成功了，可为什么求解后的图中却没有效果，我刚刚接触UDF，希望各位帮帮忙

#include "udf.h"
#include "dpm.h"
#include "math.h"
#define M 500000.00 //介质磁化强度
#define A 0.5 //磁介质直径
#define B 0.000006 //颗粒直径
#define PI 3.1415 //圆周率
#define H 796178.00 //背景磁场强度
#define Miu 0.000001256 //真空磁导率
#define XP 0.025 //颗粒磁化率
DEFINE_DPM_BODY_FORCE(particle_body_force,p,i)
{real bforce;
real x,y;
x=P_POS(p)[0];
y=P_POS(p)[1];
if(i==0)
bforce=-4*PI*Miu*XP*M*A*A*B/2*B/2*B/2*((M*A*A*x-4*H*y*y*x-4*H*x*y*y)/
(2*pow(x*x +y*y,1.5))+H*x/(pow(x*x+y*y,0.5)))/(3*pow(x*x+y*y,1.5));
else if(i==1)
bforce=-4*PI*Miu*XP*M*A*A*B/2*B/2*B/2*((M*A*A*y-4*H*y*y*y+4*H*x*x*y)/
(2*pow(x*x+y*y,1.5))+H*y/(pow(x*x+y*y,0.5)))/(3*pow(x*x+y*y,1.5));
return (bforce/P_MASS(p)); }