ansys命令流源文件.doc

磁悬浮+球形载流线圈+2-14+2-15的ansys命令流文件磁悬浮仿真! 磁悬浮! 铝板有涡流! 谐波分析finish/clear! 定义参数，单位均采用国际制单位cr1=0.031 ! 盘状线圈内半径cr2=0.195 ! 盘状线圈外半径ch=0.0125 ! 盘状线圈高度n=250 ! 线圈匝数!lh=0.002 !lh=0.007 ! 铝板高度lh=0.014 ! 铝板高度lr1=0.02 ! 铝板内半径lr2=0.25 ! 铝板外半径 pi=2*asin(1) ! 3.1415926!xfh=0.005 ! 线圈悬浮高度xfh=0.01 ! 线圈悬浮高度w=lr2 ! 场域外空气范围相关尺寸h=lh+xfh+chim0=8.27 ! 线圈电流js0=im0*n*sqrt(2)/(cr2-cr1)*ch) ! 线圈截面上的电流密度（幅值） ! 前处理/prep7 ! 前处理et, 1, plane53, , , 1 ! 指定单元类型，轴对称场分析mp, murx, 1, 1 ! 指定1号材料（空气）的相对磁导率mp, murx, 2, 1 ! 指定2号材料（线圈）的相对磁导率mp, murx, 3, 1 ! 指定3号材料（铝板）的相对磁导率!mp, rsvx, 3, 2.62e-6 ! 指定3号材料（铝板）的电阻率mp, rsvx, 3, 2.62e-8 ! 指定3号材料（铝板）的电阻率! 建立几何模型rectng, 0, lr1, 0, lh ! 铝板rectng, 0, lr2, 0, lh rectng, 0, cr1, lh+xfh, lh+xfh+chrectng, 0, cr2, lh+xfh, lh+xfh+chrectng, 0, lr2, 0, lh+xfh+chrectng, 0, w+h, -h, 2*hrectng, 0, w+8*h, -5*h, 6*haovlap, all! 对几何模型（即，面）设置属性 ! 选择线圈所对应的面，根据位置来选择asel, s, loc, x, cr1, cr2asel, r, loc, y, lh+xfh, lh+xfh+chlsel, s, loc, x, cr1asll, raatt, 2, , 1, 0, ! 选择铝板asel, s, loc, x, lr1, lr2asel, r, loc, y, 0, lhlsel, s, loc, x, lr1asll, raatt, 3, , 1, 0, ! 选择空气allselasel, u, mat, , 2, 3aatt, 1, , 1, 0! 剖分，建立网格! 先划分铝板所在区域asel, s, mat, , 3esize, 0.003amesh, all! 划分线圈所在区域asel, s, mat, , 2 ! 根据材料号来选择线圈esize, 0.003 ! 单元分割数为1，即每个线圈截面就是一个单元amesh, all! 划分线圈外的空气区域lsel, s, loc, y, 0, hlsel, r, loc, x, 0asllcm, airin, areamshape, 1, 2damesh, allesize, 0.02! 划分线圈内的空气区域smrtsize, 6mshape, 1, 2d ! 三角形单元mshkey, 0 ! 自由剖分asel, s, mat, , 1cmsel, u, airinamesh, all! 加载线圈电流密度asel, s, mat, , 2eslabfe, all, js, , , , js0! 加载外边界磁力线平行边界条件allsellsel, s, ext ! 选择外边界处的线dl, all, , asym ! 磁力线平行allsel! 加载求力边界条件asel, s, mat, , 2eslacm, ccoil, elemfmagbc, ccoilallselsavefinish/soluantype, 3harfrq,50 ! 50 Hz.solvefinish/post1set, 1, , 1, 0 ! 读实部结果PLF2D,27,0,10,1 fmagsum, ccoilcsys, 0 ! 选择直角坐标系PATH,FIELD,2,10 ! DEFINE PATH WITH NAME = FIELDPPATH,1,0.02,0.014 ! DEFINE PATH POINTS BY NODEPPATH,2,0.25,0.014PDEF,By,B,yPRPATH,By ! PRINT By ALONG COIL AXIS/SHOW,GRPH,1PLPATH,By ! DISPLAY By ALONG COIL AXISSENERGY,0,1 !计算磁场能量，hamonicESEL,S,MAT,3POWERH! CALCULATE TIME-AVERAGE POWER LOSS in plane单匝线圈/PREP7!前处理ET,1,PLANE53,1Mp,murx,1,1Mp,murx,2,1RECTNG,0.05,0.051,-0.0005,0.0005, PCIRC,0,0.1,-90,90, PCIRC,0,1,-90,90, asel,allaovlap,allnumcmp,allaplot !显示面asel,s,area,1 !选择1号面AATT,2,1 !为该面赋材料属性，介质1asel,s,area,2,3 !选择2号面AATT,1,1 !为该面赋材料属性，介质1Allsel !选择所有元素SMRT,1 !设置剖分规格，数字越小，网格越小，计算越准确，但计算量增大MSHAPE,0,2D !QuadMapped剖分MSHKEY,0ASEL,ALL !选择所有面AMESH, ALL!剖分所有面finish!前处理结束/solu!进入求解器!ANTYPE,HARMONIC!HARFRQ,50ASEL,S,MAT,2 APLOT BFA,all,JS, , ,1000000,0 ALLSEL,ALL LSEL,S,EXT DL,all, ,Az,0 Allselsolve !求解FINISH !求解结束/POST1 !SET,1,LAST,1,0, , , PLF2D,27,0,10,1 PATH,l,2,30,20, PPATH,1,0,0,0,0,0, PPATH,2,0,0.1,0,0,0,PDEF, ,B,Y,AVG PLPATH,BY PRPATH,BY 2-14在ANSYS中输入命令流如下：*create,outputdata*cfopen,E:ansys_workoutput,dat*do,i,1,40*vwrite,i,VOLTG(I,0) ,VOLTG(I,1)(1x,f8.5,1x,f8.5,1x,f8.5)*enddo*END/PREP7 !前处理*DIM,VOLTG,TABLE,40,x !为函数边界条件作准备*DO,INC,1,40,1T=INC*0.1/40VOLTG(INC,0) = TVOLTG(INC,1) = 100*SIN(3.1415926/0.1*T)*ENDDO*USE,outputdataET,1,PLANE121 !定义单元类型1为PLANE121MP,PERX,1,1 !定义材料1的相对介电常数是1K,1,0,0 !定义6个关键点K,2,0.1,0K,3,0.1,0.1K,4,0.098,0.1K,5,0.002,0.1K,6,0,0.1L,1,2 !由关键点形成线段L,2,3L,3,4L,4,5L,5,6L,6,1AL,1,2,3,4,5,6 !由线段构成面aplot !显示面asel,s,area,1 !选择1号面AATT,1,1 !为该面赋材料属性，介质1Allsel !选择所有元素!SMRT,3 !设置剖分规格，数字越小，网格越小，计算越准确，但计算量增大!MSHAPE,0,2D !QuadMapped剖分!MSHKEY,0LESIZE,all, , ,20, , , , ,1 !手动剖分：所有线分为20段ASEL,ALL !选择所有面AMESH, ALL !剖分所有面finish !前处理结束/solu !进入求解器LSEL,S,LINE,1,2 ! SELECT LOW BOUNDARYLSEL,A,LINE,6DL,ALL,VOLT,0LSEL,S,LINE,4 ! SELECT UP BOUNDARY!DL,ALL,VOLT,100NSLL,S !选择线上的节点NplotD,ALL,VOLTG, %VOLTG% ! Heat flux of 100 on all selected nodes,Allselsolve !求解FINISH !求解结束/POST1 !进入后处理PLNSOL,VOLT, ,0, !画电压云图2-15/PREP7 !前处理 ET,1,PLANE121 !定义单元类型1为PLANE121 MP,PERX,1,9 !定义材料1的相对介电常数是9 MP,PERX,2 ,1 !定义材料2的相对介电常数是1 K,1,-1,-1 !定义4个关键点 K,2,1,-1 K,3,1,1 K,4,-1,1 L,1,2 !由关键点形成线段 L,2,3 L,3,4 L,4,1 AL,1,2,3,4 !由线段构成面 CYL4,0,0,0,0,0.1,360 !画圆 asel,all aovlap,all !布尔操作构成面 numcmp,all !压缩编号 aplot !显示面 asel,s,area,1 !选择1号面 AATT,1,1 !为该面赋材料属性，介质1 asel,s,area,2 !选择2号面 AATT,2,1 !为该面赋材料属性，介质1 Allsel !选择所有元素 SMRT,3 !设置剖分规格，数字越小，网格越小，计算越准确，但计算量增大 MSHAPE,0,2D !QuadMapped剖分 MSHKEY,0 ASEL,ALL !