请问，二维轴对称模型可以用求得的电流密度来求磁场了吗？(已修改)

看了论文觉得可以。
但是自己做的时候，怎么也不能成功。
求得的磁场力为0。有人说是二维求得的电流密度方向为XOY平面，所以求磁场力必然为0。可是我的模型是轴对称的啊？到底能不能求呢？反正我的程序现在怎么着都是0。
/clear
/PREP7
ET,1,PLANE67,,,1
!First process.This process is thermal-electric coupled calculation.
!Material properity.Material 1 is the argon,material 3 is tungsten electrobe.
MPTEMP,1,273,2500,5000,7500,10000,12500
MPTEMP,7,15000,17500,20000
MPDATA,DENS,1,1,1.67,0.2,0.1,0.07,0.04,0.03
MPDATA,DENS,1,7,0.02,0.015,0.014
MPDATA,KXX,1,1,0,0.1,0.2,0.3,0.6,1.6
MPDATA,KXX,1,7,2.4,2.0,2.2
MPDATA,RSVX,1,1,1000,100,1,1.3e-3,3.3e-4,1.6e-4
MPDATA,RSVX,1,7,1.2e-4,9.8e-5,8.5e-5
MPDATA,C,1,1,521,521,521,600,1250,6000
MPDATA,C,1,7,10000,4000,2500
MPDATA,DENS,3,,19350
MPDATA,RSVX,3,,0.0000000528
MPDATA,KXX,3,,173
K,1,0,0,0,
k,2,0.01,0,0,
K,3,0.0002,0,0,
K,4,0.0016,0,0,
K,5,0,0.005,0,
K,6,0.0002,0.005,0,
K,7,0,0.005+0.0014*sqrt(3),0,
K,8,0.0016,0.005+0.0014*sqrt(3),0,
K,9,0.01,0.005+0.0014*sqrt(3),0,
K,10,0.0016,0.005,0,
K,11,0.01,0.005,0,
A,5,6,8,7
A,1,3,6,5
A,3,4,10,6
A,6,10,8
A,8,9,11,10
A,10,11,2,4
LGLUE,ALL
AGLUE,ALL
ASEL,S,AREA,,2,6
AATT,1,,1,0,
ASEL,S,AREA,,1
AATT,3,,1,0,
allsel，all
ESIZE,0.0002,0,
MSHKEY,1
AMESH,ALL
/PNUM,KP,1
/PNUM,AREA,1
APLOT
FINISH
/SOL
TUNIF,10000, !uniform temperature
DL,12,,TEMP,1000,0
LSEL,S,LINE,,1,3 !the top of the argon
DL,ALL,,TEMP,3000,0
LSEL,S,LOC,Y,0 !the surface of the tungsten electrobe.
DL,ALL,,TEMP,5000,0
DL,ALL,,VOLT,0 !the bottom of the argon(the top of the workpiece,here the workpiece is neglected)
LSEL,S,LINE,,3
NSLL,S,1 !the top of the tunsten electrobe
F,ALL,AMPS,-160 !Here is the difficulity, because I can't make the direction of current from top to bottom.
ALLSEL,ALL
CNVTOL,AMPS,,,,-1 !If not so, the calculation can't be convergened.
SOLVE
FINISH
PHYSICS,CLEAR
/PREP7
ET,1,PLANE53,,,1 !Second process: calculate magnetic force
MPTEMP,1,273,2500,5000,7500,10000,12500
MPTEMP,7,15000,17500,20000
MPDATA,DENS,1,1,1.67,0.2,0.1,0.07,0.04,0.03
MPDATA,DENS,1,7,0.02,0.015,0.014
MPDATA,KXX,1,1,0,0.1,0.2,0.3,0.6,1.6
MPDATA,KXX,1,7,2.4,2.0,2.2
MPDATA,RSVX,1,1,1000,100,1,1.3e-3,3.3e-4,1.6e-4
MPDATA,RSVX,1,7,1.2e-4,9.8e-5,8.5e-5
MPDATA,C,1,1,521,521,521,600,1250,6000
MPDATA,C,1,7,10000,4000,2500
MPDATA,DENS,3,,19350
MPDATA,RSVX,3,,0.0000000528
MPDATA,KXX,3,,173
MPDATA,MURX,1,,1
MPDATA,MURX,3,,1
FINISH
/SOL
TUNIF,10000, !Conditions of temperature are the same as the first process
DL,12,,TEMP,1000,0
LSEL,S,LINE,,1,3
DL,ALL,,TEMP,3000,0
LSEL,S,LOC,Y,0
DL,ALL,,TEMP,5000,0
ALLSEL,ALL
LDREAD,JS,,,,,,RST !Using the previous result of JS(current density)
SOLVE
FINISH
SAVE
/POST1
PLNSOL, FMAG,SUM, 0 !magnetic force