3-2.txt

ANSYS电磁分析的命令流文件

!/BATCH,LIST
/TITLE,EDDY CURRENT LOSS IN THICK STEEL PLATE (NL HARMONIC)! 定义工作标题
/FILNAM,VOL_COIL_2D,1                        ! 定义工作文件名
KEYW,MAGNOD,1                               ! 指定磁场分析

/PREP7
DD=2.5E-3                                       ! 钢板1/2高度
HM=2644.1                                      ! 输入的磁场强度
SIGMA=5.0E6                                    ! 材料的电导率
FF=50                                           ! 频率(HZ)
DX=DD/6                                        ! 钢板模型宽度
DY=DD/6                                        ! 线圈高度
!
ET,1,PLANE53                                    ! 指定单元类型
MP,MURX,1,1                                     ! 定义线圈区域材料特性
TB,BH,2,,25                                       ! 定义钢板B-H曲线
TBPT,,100,0.46512
TBPT,,200,0.72993
TBPT,,300,0.90090
TBPT,,400,1.0204
TBPT,,500,1.1086
TBPT,,600,1.1765
TBPT,,700,1.2302
TBPT,,800,1.2739
TBPT,,900,1.3100
TBPT,,1000,1.3405
TBPT,,1400,1.4257
TBPT,,1800,1.4778
TBPT,,2200,1.5131
TBPT,,2600,1.5385
TBPT,,3000,1.5576
TBPT,,3400,1.5726
TBPT,,3800,1.5847
TBPT,,4200,1.5945
TBPT,,4600,1.6028
TBPT,,5000,1.6098
TBPT,,7000,1.6332
TBPT,,9000,1.6465
TBPT,,11000,1.6551
TBPT,,13000,1.6611 
TBPT,,15000,1.6656
MP,RSVX,2,1/SIGMA
!
/PNUM,AREA,1                                ! 打开面区域编号
RECT,0,DX,0,DD                                  ! 建立钢板模型
RECT,0,DX,DD,DD+DY                            ! 建立线圈模型
AGLUE,ALL                                      ! 布尔叠分操作
SAVE,VOL_COIL_2D_GEOM.DB
!
ASEL,S,LOC,Y,DD,DD+DY                         ! 保存几何模型到文件
AATT,1,0,1                                       ! 给线圈赋予材料特性
ASEL,S,LOC,Y,0,DD
AATT,2,0,1                                       ! 给钢板赋予材料特性
ASEL,ALL
AMESH,ALL                                     ! 自由面网格划分
SAVE,VOL_COIL_2D_MESH.DB                    ! 保存网格单元数据到文件
!
NSEL,S,LOC,Y,DD+DY
CP,1,AZ,ALL                                ! 将线圈节点上的磁矢势自由度进行耦合
NSEL,ALL
FINISH
!
/SOLU
ANTYPE,HARMIC                              ! 设定求解类型为谐性
NSEL,S,LOC,Y,0
D,ALL,AZ,0                                    ! 设定边界条件为磁通量与X轴平行
NSEL,ALL
ESEL,S,MAT,,1
BFE,ALL,JS,,0,0,HM/DY,0                        ! 给线圈施加电流密度载荷
ESEL,ALL
HMAGSOLV,FF,,1E-2                            ! 进行非线性谐性电磁场求解
SAVE,VOL_COIL_2D_RESU.DB                  ! 保存计算结果
FINISH
!
/POST1
SET,LAST,1,,0                                 ! 读取最后一步求解结果
ESEL,S,MAT,,2                                ! 选择材料号为2的所有单元，即钢板
POWERH                                     ! 计算涡流损耗 (W/M)
PAVG=PAVG*2/DX                             !转化涡流损耗(W/M**2)
*CFOPEN,VOL_COIL_2D,TXT,C:\
*VWRITE,PAVG
(/'COMPUTED EDDY CURRENT LOSS:',F9.4,'WATTS/M**2')
*CFCLOS
FINISH