为什么集肤效应算出的结果不一样呢？拷贝lilulilu版主的命令流

春风化雨

旧事重提，想看到导体的集肤效应，拷贝lilulilu版主的命令流
/prep7

curr=2000.0           ! current real
curi=1000.0           ! current imaginary
et,1,solid117,1
MP,RSVX,1,4.39e-6!定义电阻率
mp,murx,1,1.! coil relative permeability

et,2,solid117
mp,murx,2,1000
MP,RSVX,2,4.39e-6
!mp,murx,2,1!空气

EMUNIT,MKS
CYL4,0,0,0.01, , , ,0.1 !定义半径10CM，高30CM
cyl4,0,0,0.02,,,,0.1
wprota,,,90
vsbw,all,,dele  
wprota,,90  
vsbw,all,,dele
VOVLAP,all

vSEL,s,,,3,6
vaTT,2,,2

vSEL,s,,,1,2
vSEL,a,,,7,8
vaTT,1,,1

mshkey,2
mshape,0,3d
esize,0.005
allsel,all
vmesh,all

finish
/solu              
csys,1  
NSEL,EXT        
D,ALL,AZ,0   

NSEL,S,LOC,Z,0
nsel,r,mat,,1
D,ALL,VOLT,0

NSEL,S,LOC,Z,0.1
nsel,r,mat,,1
CP,1,VOLT,all
*get,n1,node,,num,min  ! get a node at z=0.1
F,n1,AMPS,curr
!
NSEL,ALL         
alls
FINISH  
/SOL
ANTYPE,3
HARFRQ,1e5
SOLVE
FINISH
把频率改为1e5，其他地方都没动，出来的结果
并且求解过程中出现两个WARNING：
1.Mid-node elements forced to have straight sides.
2.The item label MAT is not recognized by the NSEL command.               
  The NSEL command is ignored.  
奇怪的是，在同学机子上运行时没有出现这些WARNING，但结果是一样的
为什么得不到lilulilu版主的结果呢，请教各位版主。

[ 本帖最后由 perch88 于 2008-7-25 06:29 编辑 ]