电磁场感应涡流问题

崔佳乐

下面曲折线圈加载余弦交流电，在磁场作用下感应涡流，但为什么会出现下面图2这种情况？按原理分析感应涡流不应该是图4吗？
这是我的命令流，求高手看看是哪里设置有问题？非常感谢啊
finish  
/clear  
/filename,EMAT_3D_A0,1                           !定义文件名
/title,EMAT_3D_A0                          !定义工作名称
/PREP7  
!定义单元
ET,1,SOLID117                               !钢板

ET,2,SOLID117
                              !线圈

ET,3,SOLID117                               !空气

ET,4,SOLID117                              !永磁体
ET,5,SOLID117
ET,6,SOLID117                           
KEYOPT,1,1,1
KEYOPT,2,1,1
EMUNIT,MKS

MP,MURX,2,1    !线圈
MP,RSVX,2,0.175e-7   

MP,MURX,1,1
MP,RSVX,1,0.283e-7                !电阻率
MP,MURX,3,1                      !空气

HC=895000

MP,MURX,4,11
MP,MURX,5,11

MP,MURX,6,1000

MP,MGXX,4,0
MP,MGYY,4,HC
MP,MGZZ,4,0

MP,MGXX,5,0
MP,MGYY,5,-HC
MP,MGZZ,5,0

BLOCK,,0.0005,,0.0005,,0.0355           !!!线圈模型
VGEN,12,ALL,,,0.0030,,,,0
VSEL,NONE
BLOCK,,0.0035,,0.0005,-0.0005,0
VGEN,6,ALL,,,0.0060,,,,0
VSEL,NONE
BLOCK,0.0030,0.0065,,0.0005,0.0355,0.0360
VGEN,5,ALL,,,0.0060,,,,0
VSEL,NONE
BLOCK,0.0035,0.0060,,0.0005,-0.0005,0
VGEN,5,ALL,,,0.0060,,,,0
VSEL,NONE
BLOCK,,0.0005,,0.0005,0.0355,0.0360           
VGEN,2,ALL,,,0.0330,,,,0

VSEL,NONE
BLOCK,0.0005,0.0030,,0.0005,,0.0355           
VGEN,11,ALL,,,0.0030,,,,0
VSEL,NONE
BLOCK,0.0005,0.0030,,0.0005,0.0355,0.0360
VGEN,6,ALL,,,0.0060,,,,0

ALLSEL
VGLUE,ALL
NUMCMP,VOLU

!磁铁
BLOCK,-0.0450,-0.0050,,0.0150,-0.0070,0.0430
BLOCK,0.0385,0.0785,,0.0150,-0.0070,0.0430
BLOCK,-0.0450,0.0785,0.0150,0.0350,-0.0070,0.0430

!钢板模型
BLOCK,-0.0600,0.0935,-0.0010,-0.0310,-0.0240,0.0600

!空气模型
BLOCK,-0.0050,0.0385,,0.0150,-0.0070,0.0430 !线圈上空气
BLOCK,-0.0600,0.0935,-0.0010,0,-0.0240,0.0600!钢板上空气
BLOCK,-0.25,0.25,-0.25,0.25,-0.25,0.25

VOVLAP,ALL                                  !布尔搭接操作
NUMCMP,VOLU  

VSEL,S,,,54
VSEL,INVE                              
VSEL,U,,,53
VSEL,U,,,51

!为实体分配单元属性
!线圈

VSEL,S,,,1,13,1
VSEL,A,,,25
VSEL,A,,,37,47,1
VATT,2,,2,0  

!永磁体
VSEL,S,,,48
VATT,4,,4,0  
VSEL,S,,,49
VATT,5,,5,0  
VSEL,S,,,51
VATT,6,,6,0
!钢板
VSEL,S,,,50
VATT,1,,1,0
!空气
VSEL,S,MAT,,1                          
VSEL,A,MAT,,2   
VSEL,A,MAT,,4
VSEL,A,MAT,,5     
VSEL,A,MAT,,6         
VSEL,INVE
VATT,3,,3,0
!网格划分
!钢板
LSEL,S,,,211
LSEL,A,,,215
LESIZE,ALL,,,40,,,,,0
LSEL,S,,,210
LESIZE,ALL,,,10,,,,,0
ALLSEL
MSHAPE,0,3D
MSHK,1
VMESH,50

VSEL,S,,,1,49,1
VSEL,A,,,51,53,1
MSHAPE,1,3D
MSHK,0
VMESH,ALL         


VSEL,S,,,54
SMRT,4
VMESH,ALL

                          !加电压耦合
ASEL,S,,,170
NSLA,S,1
CP,1,VOLT,ALL



!求解
FINISH  
/SOL                      !进入求解器
ALLSEL
ANTYPE,TRANS
OUTRES,ALL,,
TRNOPT,FULL
ALLSEL
TM_START=1E-6
TM_END=1.25E-2
TM_INCR=0.069995E-2
PHI=0
AUTOTS,ON
*DO,TM,TM_START,TM_END,TM_INCR
TIME,TM
PHI=PHI+1
I=4*COS(PHI)
AS=6e-7
JS=I/AS

allsel                 !选择所有
!*  
ANTYPE,4               !选择瞬态分析
LUMPM,0
!*
ASEL,S,EXT             !选择所有外表面
DA,all,AZ,0,0          !指定电磁为0边界条件


DA,176,VOLT,0         !线圈另一面指定电压为0

FK,230,AMPS,JS      !给线圈耦合面上的一点施加总电流

allsel
NSUBST,2,0,0           !每步指定2子步
OUTRES,ERASE
OUTRES,ALL,1           !输出所有
AUTOTS,0               !关闭自动荷载步
solve                  !求解
SAVE,EMAT_3D_A0_resu.db              
*enddo                 !结束循环，与*do结合使用
FINISH

kathyYu

唔。。。好像好深奥的样子，涡流的不太懂啊
不过“PHI=0
  PHI=PHI+1
I=4*COS(PHI)”
但是FK赋予的是电流峰值，所以还用cos有意义吗？而且，PHI都是整数，不明白你这里是什么意思

崔佳乐

哦，我就是要加载一个余弦电流，峰值为4A,我看到一个实例上是这么加的。如果不对，那应该怎么加呢？希望你多多指导。

kathyYu

指导不敢当，本人也是菜鸟一枚。只是我觉得直接把
FK,230,AMPS,JS      !给线圈耦合面上的一点施加总电流，改为
FK,230,AMPS,4
就可以了啊，我之前都是这么做的。余弦电流不是在实际中和正弦没差吗，直接加就好，不过你要是研究一个周期的变化就另当别论了。
你加载的那个地方好复杂，尤其是那个do循环，因为不太清楚你具体设定的意义，只是感觉不太对劲，你好好看看吧。

崔佳乐

FK,230,AMPS,4 如果这样设的话，不就是在节点230处加4A的电流，那不就是直流了吗？我需要加的是交流电。

zgtly0723

值得借鉴，虽有不妥之处！