anys中求气隙磁密谐波的傅里叶分析中遇到的问题

环佩叮当

这是一段在三维模型里求永磁电机气隙磁密谐波的一段代码，是在别人的基础上改的，其中有错误，输出的结果不对，恳请各位大侠不吝赐教，指出错误，不甚感激！
我感觉有两个地方可能出错，一个就是tracdata(i,5)，tracdata是一个数组，我不明白的是后面这个5是什么意思；另一个是
VFILL,THETA(1),RAMP,0,js/4，这里面js是极对数，不知道这句什么意思。
另外还有一个问题，我的模型是三维的，求气隙磁密的话，路径是不是应该是一个面，而不是一条线？如果是面的话，具体怎么做啊

PATH,bm_plot,2,30,50, ！定义路径
WPOFF,0,0   
CSWPLA,1000,1
CSCIR,1000,1
ppath,1,,(dz_nj-zz_qx/2)*cos(10*pi/180),(dz_nj-zz_qx/2)*sin(10*pi/180),zz1+yct_k/2
ppath,2,,dz_nj-zz_qx/2,0.0000001,zz1+yct_k/2
PMAP,'UNI '
/PBC,PATH,1
PDEF, ,B,x,AVG
PAGET,TRACPOIN,POIN
PAGET,TRACDATA,TABL    ！将磁密B值附到数组TRACDATA里
PAGET,TRACLABL,LABE
*DIM,COEFF,,43
*DIM,MODE,TABLE,43
*DIM,ISYM,TABLE,43
*DIM,THETA,TABLE,50+1
*DIM,CURVEI,TABLE,50+1               
*DIM,CURVEO,TABLE,50+1
*VFILL,THETA(1),RAMP,0,js/4  
*do,i,1,50+1,1
CURVEI(i)=tracdata(i,5)   ！这句不懂
*enddo
MODE(1)=0
ISYM(1)=1
ISTART=2
*DO,i,ISTART,43,2
        MODE(i)=i/2
        ISYM(i)=1  
*enddo
ISTART=3
*DO,i,ISTART,43,2
        MODE(i)=(i/2)-.5
        ISYM(i)=-1
*enddo
*MFOURI,FIT,COEFF(1),MODE(1),ISYM(1),THETA(1),CURVEI(1)
*MFOURI,EVAL,COEFF(1),MODE(1),ISYM(1),THETA(1),CURVEO(1)

ct_factor=2.0*pi*(dz_nj-zz_qx/2)*yct_k*jh_xs/js/2
n_harm=1

*cfopen,Noload_Field_Result,txt,,append
*cfwrite,'n_harm_ct',SQRT(COEFF(2)*COEFF(2)+COEFF(3)*COEFF(3))*ct_factor!!!!!!!!!!基波磁通
*cfclose

*cfopen,Noload_Field_Result,txt,,append
*cfwrite,'n_harm', 0,COEFF(1)
*cfclose
*do,i,1,21,1
*cfopen,Noload_Field_Result,txt,,append
*cfwrite,'n_harm',COEFF(2*i),COEFF(2*i+1)!!!!!!!!!!!磁密的傅里叶分解
*cfclose  
*enddo

*do,i,1,50,1
*cfopen,Noload_Field_Result,txt,,append
*cfwrite,'gap_flux', i*360*2/50,tracdata(i,5) !!!!!!!!!!!磁密
*cfclose  
*enddo

环佩叮当

补充一下，数组COEFF（）里是电机磁密21次谐波峰值

环佩叮当

我已经解决问题了，谢谢关心！

沙隆巴斯

咋解决的，拿出来说说呗

cqlyj2008

回复 3# 环佩叮当
对，三维模型中的路径应该是一个面才对。但是，其实面上的节点Az值都不一样，这又如何处理呢

环佩叮当

回复 5# cqlyj2008
这个没办法，我在三维里面也是用一条线作为路径分析的，不知道怎么用面分析

环佩叮当

回复 4# 沙隆巴斯
这个程序是对的，只要改动一下这边VFILL,THETA(1),RAMP,0,js/4，这个数组里面存放的是与磁密相对应的相位