1 旋转摩擦 (1)2. 电磁三d命令流实例（论坛看到） (11)3. 帮助感应加热例子induction heating of a solid cylinder billet (15)4. 感应加热温度场的数值模拟（论文）inducheat30命令流 (19)5. 如何施加恒定的角速度？Simwe仿真论坛 (24)6. 旋转一个已经生成好的物体 (27)7. 产生这样的磁力线 (28)8. 旋转摩擦生热简单例子（二维旋转） (32)8.1. 原版 (32)8.2. 部分gui操作 (35)9. VM229 Input Listing (39)10 轴承---耦合+接触分析 (47)11. 板的冲压仿真 (52)1 旋转摩擦FINISH/FILNAME,Exercise24 !定义隐式热分析文件名/PREP7 !进入前处理器ET,1,SOLID5 !选择单元类型MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,1,,7800 !定义材料1的密度MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,C,1,,460 !定义材料1的比热MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,KXX,1,,66.6 !定义材料1的热传导系数MPTEMP,,,,,,,,MPTEMP,1,0UIMP,1,REFT,,,30 !定义材料1的热膨胀系数的参考温度MPDATA,ALPX,1,,1.06e-5 !定义材料1的热膨胀系数MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,206e9 !定义材料1的弹性模量MPDATA,PRXY,1,,0.3 !定义材料1的泊松比MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,2,,8900 !定义材料2的密度MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,C,2,,390 !定义材料2的比热MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,KXX,2,,383 !定义材料2的热传导系数MPTEMP,,,,,,,,MPTEMP,1,0UIMP,2,REFT,,,30 !定义材料2的热膨胀系数的参考温度MPDATA,ALPX,2,,1.75e-5 !定义材料2的热膨胀系数MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,2,,103e9 !定义材料2的弹性模量MPDATA,PRXY,2,,0.3 !定义材料2的泊松比CYLIND,0.7,0.5,0,0.2,0,180, !建立钢环的几何模型CYLIND,0.5,0.3,0,0.2,60,70, !建立铜滑块FLST,5,2,4,ORDE,2FITEM,5,8FITEM,5,12CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,,_YLESIZE,_Y1, , ,4, , , , ,1 !选择第8和第12号线并划分为4段FLST,5,2,4,ORDE,2FITEM,5,2FITEM,5,7CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,,_YMSHAPE,1,3DMSHKEY,0FLST,5,2,4,ORDE,2FITEM,5,2FITEM,5,7CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,,_YLESIZE,_Y1, , ,36, , , , ,1 !选择第2和第7号线并划分为36段FLST,5,1,4,ORDE,1FITEM,5,5CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINELESIZE,_Y1, , ,36, , , , ,1 !选择第5号线并划分为36段FLST,5,2,4,ORDE,2FITEM,5,17FITEM,5,19CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,,_YLESIZE,_Y1, , ,2, , , , ,1 ! 选择第17和第19号线并划分为2段FLST,5,2,4,ORDE,2FITEM,5,20FITEM,5,24CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,,_YLESIZE,_Y1, , ,4, , , , ,1 !选择第20和第24号线并划分为4段TYPE, 1 !给钢环附单元类型MA T, 1 !给钢环附材料REAL, !给钢环附实常数ESYS, 0SECNUM,CM,_Y,VOLUVSEL, , , , 1CM,_Y1,VOLUCHKMSH,'VOLU'CMSEL,S,_YVSWEEP,_Y1CMDELE,_YCMDELE,_Y1CMDELE,_Y2 !给铜滑块划分单元TYPE, 1 !给铜滑块附单元类型MA T, 2 !给铜滑块附材料REAL, !给铜滑块附实常数ESYS, 0SECNUM,CM,_Y,VOLUVSEL, , , , 2CM,_Y1,VOLUCHKMSH,'VOLU'CMSEL,S,_YVSWEEP,_Y1CMDELE,_YCMDELE,_Y2 !给铜滑块附实常数/COM, CONTACT PAIR CREATION - STARTCM,_NODECM,NODECM,_ELEMCM,ELEMCM,_KPCM,KPCM,_LINECM,LINECM,_AREACM,AREACM,_VOLUCM,VOLU/GSA V,cwz,gsav,,tempMP,MU,1,0.2MA T,1R,3REAL,3ET,2,170ET,3,174KEYOPT,3,9,0R,3,RMORE,RMORE,,0RMORE,0KEYOPT,3,1,3! Generate the target surfaceASEL,S,,,4CM,_TARGET,AREATYPE,2NSLA,S,1ESLN,S,0ESURF,ALLCMSEL,S,_ELEMCM! Generate the contact surfaceASEL,S,,,9CM,_CONTACT,AREATYPE,3NSLA,S,1ESLN,S,0ESURF,ALLALLSELESEL,ALLESEL,S,TYPE,,2ESEL,A,TYPE,,3ESEL,R,REAL,,3/PSYMB,ESYS,1/PNUM,TYPE,1/NUM,1EPLOTESEL,ALLESEL,S,TYPE,,2ESEL,A,TYPE,,3ESEL,R,REAL,,3CMSEL,A,_NODECMCMDEL,_NODECMCMSEL,A,_ELEMCMCMDEL,_ELEMCMCMSEL,S,_KPCMCMDEL,_KPCMCMSEL,S,_LINECMCMDEL,_LINECMCMSEL,S,_AREACMCMDEL,_AREACMCMSEL,S,_VOLUCMCMDEL,_VOLUCM/GRES,cwz,gsavCMDEL,_TARGETCMDEL,_CONTACT/COM, CONTACT PAIR CREATION - END !用接触向导建立铜滑块与钢环间的接触/GO*DIM,pre,TABLE,2,1,1,time,press, !定义压力载荷随时间变化的TABLE表*SET,PRE(1,0,1) , 0*SET,PRE(1,1,1) , 20e6*SET,PRE(2,0,1) , 10*SET,PRE(2,1,1) , 20e6FLST,2,1,5,ORDE,1FITEM,2,10/GOSFA,P51X,1,PRES, %PRE% !在铜滑块上表面施加压力载荷TUNIF,20, !施加初始温度FLST,2,1,5,ORDE,1FITEM,2,3/GODA,P51X,UX, !在钢环的外表面施加UX向约束FLST,2,1,5,ORDE,1FITEM,2,3/GODA,P51X,UY, !在钢环的外表面施加UY向约束FLST,2,1,5,ORDE,1FITEM,2,3/GODA,P51X,UZ, !在钢环的外表面施加UZ向约束ESEL,S,MA T,,2 !选择铜滑块单元NSLE !选择铜滑块节点CSYS,1 !激活柱坐标系FLST,2,75,1,ORDE,2FITEM,2,876FITEM,2,-950NROTA T,P51X !旋转节点坐标系FLST,2,75,1,ORDE,2FITEM,2,876FITEM,2,-950/GOD,P51X, ,0.01, , , ,UY, , , , , ! 给铜滑块加Y向位移载荷（即切向位移）ALLSEL !选择所有的节点，单元，几何体FINISH/SOL !进入求解器ANTYPE,4 !设置为瞬态分析TRNOPT,FULL !选择瞬态全分析方法LUMPM,0 !关闭集中质量矩阵计算格式NLGEOM,1 !定义大变形NSUBST,100,200,100 !定义载荷步OUTRES,ALL,20 !定义结果输出子步间隔AUTOTS,1 !打开时间开关KBC,0 !设置为斜坡载荷TIME,10 !定义求解时间SSTIF,0 !关闭应力刚化影响NROPT,FULL, , !选择瞬态全分析方法EQSLV, , ,0,PRECISION,0 !定义双精度求解MSA VE,0TOFFST,273, !定义温度偏移量SOLVE !求解/POST1 !进入后处理器SET,LAST !读取最后子步结果PLNSOL,TEMP !显示温度分布云图PLNSOL,S,EQV,2 !显示等效应力分布云图FINISH/EXIT,NOSAV !退出ANSYSFINISH/FILNAME,Exercise23 ! 定义隐式热分析文件名/PREP7 !进入前处理器ET,1,PLANE13,4 !选择单元类型MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,1,,7800 !定义材料1的密度MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,C,1,,460 !定义材料1的比热MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,KXX,1,,66.6 !定义材料1的热传导系数MPTEMP,,,,,,,,MPTEMP,1,0UIMP,1,REFT,,,30 !定义材料1的热膨胀系数的参考温度MPDATA,ALPX,1,,1.06e-5 !定义材料1的热膨胀系数MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,206e9 !定义材料1的弹性模量MPDATA,PRXY,1,,0.3 !定义材料1的泊松比MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,2,,8900 !定义材料2的密度MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,C,2,,390 !定义材料2的比热MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,KXX,2,,383 !定义材料2的热传导系数MPTEMP,,,,,,,,MPTEMP,1,0UIMP,2,REFT,,,30 !定义材料2的热膨胀系数的参考温度MPDATA,ALPX,2,,1.75e-5 !定义材料2的热膨胀系数MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,2,,103e9 !定义材料2的弹性模量MPDATA,PRXY,2,,0.3 !定义材料2的泊松比RECTNG,0,5,0,1.25, !建立钢固定块的几何模型RECTNG,0,1.25,1.25,2.5, !建立铜滑块的几何模型TYPE, 1MA T, 1REAL, !给钢固定块附材料属性ESYS, 0SECNUM,ESIZE,0.25,0, !设置单元划分尺寸MSHAPE,0,2DMSHKEY,1CM,_Y,AREAASEL, , , , 1CM,_Y1,AREACHKMSH,'AREA'CMSEL,S,_YAMESH,_Y1CMDELE,_YCMDELE,_Y1CMDELE,_Y2 !给钢固定块划分单元TYPE, 1MA T, 2REAL, !给铜滑块附材料属性ESYS, 0SECNUM,CM,_Y,AREAASEL, , , , 2CM,_Y1,AREACHKMSH,'AREA'CMSEL,S,_YAMESH,_Y1CMDELE,_YCMDELE,_Y1CMDELE,_Y2 !给铜滑块划分单元CM,_NODECM,NODECM,_ELEMCM,ELEMCM,_KPCM,KPCM,_LINECM,LINECM,_AREACM,AREACM,_VOLUCM,VOLUMP,MU,1,0.2MA T,1MP,EMIS,1,R,3REAL,3ET,2,169ET,3,172R,3,,,1.0,0.1,0,RMORE,,,1.0E20,0.0,1.0,RMORE,0.0,0,1.0,,1.0,0.5RMORE,0,1.0,1.0,0.0,,1.0 KEYOPT,3,3,0 KEYOPT,3,4,0 KEYOPT,3,5,0 KEYOPT,3,7,0 KEYOPT,3,8,0 KEYOPT,3,9,1 KEYOPT,3,10,1 KEYOPT,3,11,0 KEYOPT,3,12,0 KEYOPT,3,2,0 KEYOPT,3,1,1LSEL,S,,,3CM,_TARGET,LINE TYPE,2NSLL,S,1ESLN,S,0ESURF,ALLCMSEL,S,_ELEMCM LSEL,S,,,5CM,_CONTACT,LINE TYPE,3NSLL,S,1ESLN,S,0ESURF,ALLALLSELESEL,ALLESEL,S,TYPE,,2 ESEL,A,TYPE,,3 ESEL,R,REAL,,3 ESEL,ALLESEL,S,TYPE,,2 ESEL,A,TYPE,,3 ESEL,R,REAL,,3 CMSEL,A,_NODECM CMDEL,_NODECM CMSEL,A,_ELEMCM CMDEL,_ELEMCM CMSEL,S,_KPCM CMDEL,_KPCM CMSEL,S,_LINECM CMDEL,_LINECM CMSEL,S,_AREACM CMDEL,_AREACMCMSEL,S,_VOLUCMCMDEL,_VOLUCMCMDEL,_TARGETCMDEL,_CONTACTFLST,2,1,5,ORDE,1FITEM,2,1 !用接触向导建立滑块和固定块间的接触面/GODA,P51X,UX, !将下面的钢块施加X向约束FLST,2,1,5,ORDE,1FITEM,2,1/GODA,P51X,UY, !将下面的钢块施加Y向约束FLST,2,1,4,ORDE,1FITEM,2,6/GODL,P51X, ,UX,3.75 ! 给铜块右侧的线加X向位移载荷FINISH/SOL*DIM,pre,TABLE,2,1,1,time,press, !定义压力载荷随时间变化的TABLE表*SET,PRE(1,0,1) , 0*SET,PRE(1,1,1) , 10e6*SET,PRE(2,0,1) , 10*SET,PRE(2,1,1) , 10e6FLST,2,1,4,ORDE,1FITEM,2,7/GOSFL,P51X,PRES, %PRE% !在铜滑块上表面施加压力载荷TUNIF,30, !给铜块和钢块施加初始温度ALLSEL !选择所有的节点，单元，几何体ANTYPE,TRANS !设置为瞬态分析TIMINT,OFF,STRUC !关闭结构的动态惯性影响TINTP,,,,1.0 !定义瞬态积分参数NLGEOM,ON !定义大变形TIME,3.75E-3 !定义求解时间TOFFSET,273 !定义温度偏移量AUTO,ON !打开时间开关NSUB,100,10000,100 !定义载荷步OUTRES,ALL,-10 !定义结果输出子步间隔SOLVE !求解/POST1 !进入后处理器SET,LAST !读取最后子步结果PLNSOL,TEMP ! 显示温度分布云图PLNSOL,S,EQV,2 !显示等效应力分布云图FINISH/EXIT,NOSAV !退出ANSYS2. 电磁三d命令流实例（论坛看到）/BATCH/COM,ANSYS RELEASE 12.0.1 UP20090415 11:10:14 03/29/2011 /input,menust,tmp,'',,,,,,,,,,,,,,,,1/GRA,POWER/GST,ON/PLO,INFO,3/GRO,CURL,ON/CPLANE,1/REPLOT,RESIZEWPSTYLE,,,,,,,,0/TITLE, Test for EDGE-Formulation with SOLID117!! Parameters*SET,l,5e-6 ! length*SET,h,1e-6 ! height*SET,w,5e-6 ! width*SET,a,3e-6 ! 1st addition length of magnetic core*SET,b,5e-6 ! 2nd addition length of magnetic core!/PREP7 !!ET,1,117 !!MP,MURX,1,1MP,RSVX,1,17.2e-9 !MP,MURX,2,500MP,MURX,3,1BLOCK,-w,w,-h,h,-l,l !!BLOCK,-w-a,w+a,h+a,h+2*a,-l,lBLOCK,w+a,w+a+b,h+a,h+2*a,-l,lBLOCK,-w-a,-w-a-b,h+a,h+2*a,-l,lBLOCK,w+a,w+a+b,h,h+a,-l,lBLOCK,-w-a,-w-a-b,h,h+a,-l,l!BLOCK,-w-a,w+a,-(h+a),-(h+2*a),-l,lBLOCK,w+a,w+a+b,-(h+a),-(h+2*a),-l,lBLOCK,-w-a,-w-a-b,-(h+a),-(h+2*a),-l,lBLOCK,w+a,w+a+b,-h,-(h+a),-l,lBLOCK,-w-a,-w-a-b,-h,-(h+a),-l,l!BLOCK,-22e-6,22e-6,-15e-6,15e-6,-15e-6,15e-6 VOVLAP,ALLNUMCMP,VOLUVSEL,S,VOLU,,1V ATT,1VSEL,S,VOLU,,2,11,1V ATT,2VSEL,S,VOLU,,12V ATT,3VSEL,S,MA T,,1,2,1ASLV,SLSLA,SLESIZE,ALL,3e-6MSHKEY,1MSHAPE,0,3DVMESH,ALLVSEL,S,MA T,,3ASLV,SLSLA,SLESIZE,ALL,3e-6MSHKEY,0MSHAPE,1,3DVMESH,ALLALLSEL,ALLFINISH/SOLUESEL,S,MA T,,1BFE,ALL,JS,1,,,-1e9ASEL,S,EXTDA,ALL,AZ,0ALLSEL,ALLANTYPE,0NROPT,AUTO,,EQSLV,SPAR,,0,PRECISION,0MSA VE,0PIVCHECK,1/STA TUS,SOLUSOLVEFINISH/POST1! elements/VSCALE,1,1,0!!*PLVECT,B, , , ,VECT,ELEM,ON,0ESEL,U,MAT,,3/REPLOT[/color][/size][/b][/align][b][color=#0000ff][/color][/b][align=left][code]/BA TCH,LIST /TITLE, Test for EDGE-Formulation with SOLID117!! Parametersl=5e-6 ! lengthh=1e-6 ! heightw=5e-6 ! widtha=3e-6 ! 1st addition length of magnetic coreb=5e-6 ! 2nd addition length of magnetic core! *** Create Model!/PREP7 ! Enter ANSYS preprocessor to define problem!ET,1,117 ! Element type #1 is a magnetic edge element, 117.! (We use only one in this example)!MP,MURX,1,1 ! Define relative permeability of material #1.! Material #1 is the conductor made from copperMP,RSVX,1,17.2e-9 ! Define electric resistivity of material #1.!MP,MURX,2,500 ! Define relative permeability of material #2.! Material #2 is the magnetic core made from NiFe!MP,MURX,3,1 ! Define relative permeability of material #3.! Material #3 is the air ambient!! Define rectangular block of conductor volume! regionBLOCK,-w,w,-h,h,-l,l !!! Define rectangular block of upper magnetic core! volume regionBLOCK,-w-a,w+a,h+a,h+2*a,-l,lBLOCK,w+a,w+a+b,h+a,h+2*a,-l,lBLOCK,-w-a,-w-a-b,h+a,h+2*a,-l,lBLOCK,w+a,w+a+b,h,h+a,-l,lBLOCK,-w-a,-w-a-b,h,h+a,-l,l!! Define rectangular block of lower magnetic core ! volume region BLOCK,-w-a,w+a,-(h+a),-(h+2*a),-l,lBLOCK,w+a,w+a+b,-(h+a),-(h+2*a),-l,lBLOCK,-w-a,-w-a-b,-(h+a),-(h+2*a),-l,lBLOCK,w+a,w+a+b,-h,-(h+a),-l,lBLOCK,-w-a,-w-a-b,-h,-(h+a),-l,l!! Define rectangular block of air ambientBLOCK,-22e-6,22e-6,-15e-6,15e-6,-15e-6,15e-6! Overlap volumesVOVLAP,ALLNUMCMP,VOLU! Assign attributes to the volumesVSEL,S,VOLU,,1V ATT,1 ! conductorVSEL,S,VOLU,,2,11,1 ! magnetic coreV ATT,2VSEL,S,VOLU,,12 ! air ambientV ATT,3! Mesh volumesVSEL,S,MA T,,1,2,1 ! Conductor and magnetic core with hexaelements ASLV,SLSLA,SLESIZE,ALL,3e-6MSHKEY,1 ! MAPPED VOLUME MESHMSHAPE,0,3D ! USING HEXVMESH,ALLVSEL,S,MA T,,3 ! Air ambientASLV,SLSLA,SLESIZE,ALL,3e-6MSHKEY,0 ! FREE VOLUME MESHMSHAPE,1,3D ! USING TETVMESH,ALLALLSEL,ALLFINISH/SOLU!*** Apply current density on conductor ***ESEL,S,MA T,,1 ! Select coil elementBFE,ALL,JS,1,,,-1e9 ! Apply current density!*** Flux parallel boundary condition at all extern area of air ambient *** ASEL,S,EXTDA,ALL,AZ,0ALLSEL,ALL! Solution parametersANTYPE,0NROPT,AUTO,,EQSLV,SPAR,,0,PRECISION,0MSA VE,0PIVCHECK,1/STA TUS,SOLUSOLVE/POST1!*** Vectorplot of magnetic flux (BSUM) for magnetic core's and conductor'selements ***/VSCALE,1,1,0!!*PLVECT,B, , , ,VECT,ELEM,ON,0ESEL,U,MAT,,3/REPLOTFINISH! /EXIT,MODEL[/code][/align][align=center][size=4][color=blue][b]===错因告诉你了，稍微改改就可以得到正确的解，来改着玩玩吧，奖分哦（+1技术分，+50贡献分）==[/b][/color][/size][/align]3. 帮助感应加热例子induction heating of a solid cylinder billetbatch,list/filenam, induc/config,nres,100000 ！将载荷步最大设置为100000/prep7shpp,off/title, heating of a solid cylinder billet/com,et,1,13,,,1 ! PLANE13, axisymmetric, AZ dofet,2,13,,,1et,3,151,,,1,1,1 ! SURF151, thermal, radiationr,3,0 ! Real constant set for SURF151row=.015 ! outer radius of workpieceric=.0175 ! inner radius of coilroc=.0200 ! outer radius of coilro=.05 ! outer radius of modelt=.001 ! model thicknessfreq=150000 ! frequency (Hz.)pi=4*atan(1) ! picond=.392e7 ! maximum conductivitymuzero=4e-7*pi ! free-space permeabilitymur=200 ! maximum relative permeabilityskind=sqrt(1/(pi*freq*cond*muzero*mur)) ! skin depthftime=3 ! final timetinc=.05 ! time increment for harmonic analysistime=0 ! initialize timedelt=.01 ! maximum delta time stepemunit,mks ! set magnetic unitsmp,murx,1,1 ! air relative permeabilitymp,murx,3,1 ! coil relative permeabilitymptemp,1,25.5,160,291.5,477.6,635,698 ! temps for relative permeabilitymptemp,7,709,720.3,742,761,1000mpdata,murx,2,1,200,190,182,161,135,104 ! steel relative permeability mpdata,murx,2,7,84,35,17,1,1mptempmptemp,1,0,125,250,375,500,625 ! temps for resistivity mptemp,7,750,875,1000mpdata,rsvx,2,1,.184e-6,.272e-6,.384e-6,.512e-6,.656e-6,.824e-6 mpdata,rsvx,2,7,1.032e-6,1.152e-6,1.2e-6 ! steel resistivityrectng,0,row,0,t ! billetrectng,row,ric,0,t ! air-gaprectng,ric,roc,0,t ! coilrectng,roc,ro,0,t ! outer airaglue,allnumcmp,area ！对面重新编号ksel,s,loc,x,row ! select keypoints at outer radius of workpiece kesize,all,skind/2 ! set meshing size to 1/2 skin depth(为什么是半个集肤深度)ksel,s,loc,x,0 ! select keypoints at centerkesize,all,40*skind ! set meshing sizelsel,s,loc,y,t/2 ! select vertical lineslesize,all,,,1 ! set 1 division through thicknesslsel,allasel,s,area,,1aatt,2,1,1 ! set attributes for billet regionasel,s,area,,3aatt,3,1,2 ! set attributes for coil regionasel,s,area,,2,4,2aatt,1,1,2 ! set attributes for air regionasel,allmshape,0,2dmshk,1amesh,1 ! mesh billet arealsel,s,loc,y,0lsel,a,loc,y,t ！also selectlsel,u,loc,x,row/2lesize,all,.001lsel,allamesh,all ! mesh remaining areasn ! create space node for SURF151*get,nmax,node,,num,maxlsel,s,loc,x,rowtype,3real,3 ! REAL, NSET Sets the element real constant set attribute pointer.mat,2lmesh,all ! mesh billet outer radius with SURF151*get,emax,elem,,num,maxemodif,emax,3,nmax ! modify element to add space node for radiation et,3,0 ! reset type 3 to null elementnsel,s,loc,xd,all,az,0 ! apply flux-normal b.c. az向量磁位能nsel,allesel,s,mat,,3bfe,all,js,,,,15e6 ! apply current density to coil施加体力载荷esel,allfinish/soluantyp,harmharfrq,150000physics,write,emag ! write emag physics filefinish/prep7lsclear,all ! clear all b.c.'s and optionset,1,55,,,1 ! PLANE55 thermal element, axisymmetricet,2,0 ! null element type for coil and air region et,3,151,,,1,1,1 ! SURF151 element for radiationkeyopt,3,9,1r,3,1,5.67e-8 ! form factor, Stefan-Boltzmann constant mptempmptemp,1,0,730,930,1000 ! temps for conductivitympdata,kxx,2,1,60.64,29.5,28,28mptemp ! temps for enthalpymptemp,1,0,27,127,327,527,727mptemp,7,765,765.001,927mpdata,enth,2,1,0,91609056,453285756,1.2748e9,2.2519e9,3.3396e9 mpdata,enth,2,7,3.548547e9,3.548556e9,4.3520e9mp,emis,2,.68 ! emissivityfinish/soluantype,transtoffst,273tunif,100 ! initial uniform temperatured,nmax,temp,25 ! ambient temperaturecnvtol,heat,1 ! convergence tolerancekbc,1 ! step loadstrnopt,fullautos,on ! auto time-steppingdeltim,1e-5,1e-6,delt,on ! time step controloutres,basic,all ! save all load step information physics,write,thermal ! write thermal physics filefinish*do,i,1,ftime/tinc ! solution *do looptime=time+tinc ! increment timephysics,read,emag ! read emag physics file/solu*if,i,eq,1,thentunif,100 ! initial temperature*elseldread,temp,last,,,,,rth ! read thermal analysis temperatures*endifsolve ! solve harmonic analysisfinishphysics,read,thermal ! read thermal physics file/assign,esav,therm,esav ! redirect files for use in thermal restart/assign,emat,therm,emat/solu*if,i,gt,1,thenantype,trans,rest ! thermal restart*endiftime,time ! time at end of thermal runesel,s,mat,,2 ! select billet regionldread,hgen,,,,2,,rmg ! apply coupled joule heating load from emagesel,allsolvefinish/assign,esav ! reassign files to default/assign,emat*enddo ! end of solution loopingfinishsave ! save database/post26 ! time-history postprocessor/shownsol,2,1,temp,,tempcl ! store temperature at billet centerline nsol,3,2,temp,,tempsurf ! store temperature at billet outer diameterplvar,2,3 ! plot temperature rise over time prvar,2,3 ! print temperature rise over time finishn 定义节点。