求助！！超磁致伸缩换能器的磁分析和热分析

sxqljj

本人毕业论文是做做超磁致伸缩换能器的磁分析和热分析，谐波分析，线圈加电流密度，频率18KHZ，模型为2D轴对称结构，目前有2个问题急需解决，希望各位大侠帮帮忙：
1.以下是2维建模的轴对称模型
.其中绿色标1是永磁体，2是线圈，3是超磁致伸缩材料(磁导率我设为常数），红色4是厄铁（纯铁）蓝色5为上下端盖（相当磁导率为1）
其中标3的磁致伸缩材料结构部分，为了减小涡流，采用了硅钢薄片叠压的方式，为了模拟这种模型，在模型中我对3部分进行了面切割，也就是3本来是个完整的矩形面，我采用切割将它分割成几个小矩形，每2个矩形面之间有间隙，看成是空气，利用这样来模拟，不知道是不是达到了真实硅钢片的效果？？？？还有对面3切割后，赋予材料属性和划分网格的时候是不是要有特殊的要求？？？？？
2.在磁场分析完后，得到了3部分的涡流损耗，我还想继续进行热分析，想得到3部分的温度分布，以及线圈的温度分布，接下来应该怎么来做呢？？？需要哪些命令呢？？？求指点啊！！多谢大家
以下为命令流：
fini
/clear,start
/TITLE,3D TRANSIENT ANALYSIS
KEYW,MAGNOD,1

/PREP7
ET,1,53,0,,1
ET,2,53,0,,1  
MP,MURX,1,1                         !定义空气的材料特性(A1)
MP,MURX,2,1                           !定义线圈的材料特性(A3面)
MP,RSVX,2,3E-8
MP,MURX,3,2000                        !定义厄铁的导磁率(A4面)
MP,RSVX,3,1E-7                        !定义厄电阻率
MP,MURX,4,1                           !定义上/下端盖，顶杆(A5面)
MP,MURX,5,9                           !定义稀土棒的材料(A2面)
MP,RSVX,5,6.6E-7
MP,MURX,6,1.1

/PNUM,AREA,1
CYL4,0,0,0,-90,0.20,90                 !定义空气半圆面   
RECTNG,0,0.002,0,0.086                !输出顶杆
RECTNG,0.002,0.0055,0.018,0.068        !定义稀土棒棒的面
RECTNG,0.002,0.0055,0.068,0.072        !定义上永磁体
RECTNG,0.002,0.0055,0.014,0.018        !定义下永磁体
RECTNG,0.0065,0.014,0.014,0.072        !定义线圈的面
RECTNG,0.002,0.014,0.072,0.076         !定义上厄铁的面
RECTNG,0.002,0.014,0.010,0.014         !定义下厄铁的面
RECTNG,0.014,0.016,0.010,0.076         !定义外壳
RECTNG,0.002,0.018,0.004,0.010         !定义下端盖
RECTNG,0.002,0.018,0.076,0.086         !定义上端盖
RECTNG,0.0055,0.0065,0.014,0.072         !定义骨架
AOVLAP,ALL       !搭接
NUMCMP,AREA
!*
WPSTYLE,,,,,,,,1         ! 打开 WP 显示 - Workplane > Display Working Plane
KWPAVE,      11          ! 打开 WP 工具条 － workPlane > Offset WP by Increments...
wpro,,,90.000000         ! 使用 WP 工具条，将 WP 绕 WY 轴旋转 90 度，使 WX-WY 平面垂直面 A1
!*
wpoff,,,0.001        ! 使用 WP 工具条，将 WP 沿 WZ 方向移动 0.001
ASBW,       1            ! 使用 WP 切割 A1，切割后待进一步切割部分为 A13
!*
wpoff,,,0.00025          ! 使用 WP 工具条，将 WP 沿 WZ 方向移动 0.00025 (间隙大小)
ASBW,      13            ! 使用 WP 切割 A13，切割后待进一步切割部分为 A15
!*
wpoff,,,0.001        ! 使用 WP 工具条，将 WP 沿 WZ 方向移动 0.001
ASBW,      15           ! 使用 WP 切割 A15，切割后待进一步切割部分为 A13  
!*
wpoff,,,0.00025          ! 使用 WP 工具条，将 WP 沿 WZ 方向移动 0.00025
ASBW,      13            ! 使用 WP 切割 A13，切割后待进一步切割部分为 A15  
!*  
  
!*  注意：由于切割时会重新使用以前操作空出的面号，有时同一个命令流在不同机器上运行会出错。
!*
NUMCMP,AREA
APLOT
WPSTYLE,,,,,,,,0           !关闭WP显示
save

ASEL,S,,,1
AATT,1,,1

ASEL,S,,,2
AATT,2,1,2

ASEL,S,,,3
AATT,6,,1

ASEL,S,,,4
AATT,6,,1
ASEL,S,,,5
AATT,1,,1
ASEL,S,,,6
AATT,3,,1

ASEL,S,,,7
AATT,3,,1
ASEL,S,,,8
AATT,3,,1
ASEL,S,,,9
AATT,4,,1

ASEL,S,,,10
AATT,4,,1

ASEL,S,,,11
AATT,4,,1

ASEL,S,,,12
AATT,1,,1

ASEL,S,,,13
AATT,5,,1

ASEL,S,,,14
AATT,5,,1

ASEL,S,,,15
AATT,5,,1

ASEL,S,,,16
AATT,1,,1

ASEL,S,AREA,,2
ESIZE,0.5/1000
MSHAPE,1,2D
MSHKEY,0
AMESH,all

ASEL,S,AREA,,3,4
ESIZE,0.5/1000
MSHAPE,1,2D
MSHKEY,0
AMESH,all

ASEL,S,AREA,,1
ESIZE,0.2/1000
MSHAPE,1,2D
MSHKEY,0
AMESH,all

ASEL,S,AREA,,16
ESIZE,0.2/1000
MSHAPE,1,2D
MSHKEY,0
AMESH,all

ASEL,S,AREA,,13,15
ESIZE,0.5/1000
MSHAPE,1,2D
MSHKEY,0
AMESH,all

ASEL,S,AREA,,6,8
asel,a,area,,9,11
ESIZE,0.5/1000
MSHAPE,1,2D
MSHKEY,0
AMESH,all

ASEL,S,AREA,,5
ESIZE,0.2/1000
MSHAPE,1,2D
MSHKEY,0
AMESH,all

ASEL,S,AREA,,12
ESIZE,5/1000
MSHAPE,1,2D
MSHKEY,0
AMESH,all

!ASEL,S,,,2
!ESEL,S,MAT,,2
!NSLE,S
!CP,1,AZ,ALL

allsel
CSYS,1
NSEL,S,LOC,X,0.2
D,ALL,AZ,0
CSYS,0
NSEL,S,LOC,X,0
D,ALL,AZ,0
ALLSEL
FINISH

/SOLU
ANTYPE,HARMIC
!ASEL,S,AREA,,2
!ESLA,S
esel,s,mat,,2
BFE,ALL,JS,,,,2E6,0
ESEL,ALL
HARFRQ,18000
SOLVE
FINISH
/POST1
set,1
ESEL,S,MAT,,1,6
PLNSOL,B,SUM
PLF2D,50
set,last,,,1
ESEL,S,MAT,,2,7
PLNSOL,B,SUM
PLF2D,30

sxqljj

希望各位大大给予帮助 谢谢了啊

sxqljj

有谁知道是怎么回事吗？？

sxqljj

急需帮助啊啊 啊

anonymouss

4# sxqljj
别指望谁给你解决问题，这个版面基本上不解决问题；
大家都热衷于学习，讨论过去的帖子，对你提出的新问题，基本上是围观。
我对你做的不了解，爱莫能助；有什么不懂，还是挖掘自身的潜力吧