ansys旋轉經(jīng)典命令流

1 旋轉摩擦12. 電磁三d命令流實例（論壇看到）113. 幫助感應加熱例子induction heating of a solid cylinder billet154. 感應加熱溫度場的數(shù)值模擬（論文）inducheat30命令流195. 如何施加恒定的角速度？Simwe仿真論壇246. 旋轉一個已經(jīng)生成好的物體277. 產(chǎn)生這樣的磁力線288. 旋轉摩擦生熱簡單例子（二維旋轉）328.1. 原版328.2. 部分gui操作359. VM229 Input Listing3910 軸承-耦合+接觸分析4711. 板的沖壓仿真521 旋轉摩擦FINISH/FILNAME,Exercise24 !定義隱式熱分析文件名/PREP7 !進入前處理器ET,1,SOLID5 !選擇單元類型MPTEMP, MPTEMP,1,0 MPDATA,DENS,1,7800 !定義材料1的密度MPTEMP, MPTEMP,1,0 MPDATA,C,1,460 !定義材料1的比熱MPTEMP, MPTEMP,1,0 MPDATA,KXX,1,66.6 !定義材料1的熱傳導系數(shù)MPTEMP, MPTEMP,1,0 UIMP,1,REFT,30 !定義材料1的熱膨脹系數(shù)的參考溫度MPDATA,ALPX,1,1.06e-5 !定義材料1的熱膨脹系數(shù)MPTEMP, MPTEMP,1,0 MPDATA,EX,1,206e9 !定義材料1的彈性模量 MPDATA,PRXY,1,0.3 !定義材料1的泊松比MPTEMP, MPTEMP,1,0 MPDATA,DENS,2,8900 !定義材料2的密度MPTEMP, MPTEMP,1,0 MPDATA,C,2,390 !定義材料2的比熱MPTEMP, MPTEMP,1,0 MPDATA,KXX,2,383 !定義材料2的熱傳導系數(shù)MPTEMP, MPTEMP,1,0 UIMP,2,REFT,30 !定義材料2的熱膨脹系數(shù)的參考溫度MPDATA,ALPX,2,1.75e-5 !定義材料2的熱膨脹系數(shù)MPTEMP, MPTEMP,1,0 MPDATA,EX,2,103e9 !定義材料2的彈性模量 MPDATA,PRXY,2,0.3 !定義材料2的泊松比CYLIND,0.7,0.5,0,0.2,0,180, !建立鋼環(huán)的幾何模型CYLIND,0.5,0.3,0,0.2,60,70, !建立銅滑塊FLST,5,2,4,ORDE,2 FITEM,5,8 FITEM,5,12 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,_Y LESIZE,_Y1, , ,4, , , , ,1 !選擇第8和第12號線并劃分為4段FLST,5,2,4,ORDE,2 FITEM,5,2 FITEM,5,7 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,_Y MSHAPE,1,3D MSHKEY,0FLST,5,2,4,ORDE,2 FITEM,5,2 FITEM,5,7 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,_Y LESIZE,_Y1, , ,36, , , , ,1 !選擇第2和第7號線并劃分為36段FLST,5,1,4,ORDE,1 FITEM,5,5 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,_Y LESIZE,_Y1, , ,36, , , , ,1 !選擇第5號線并劃分為36段 FLST,5,2,4,ORDE,2 FITEM,5,17 FITEM,5,19 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,_Y LESIZE,_Y1, , ,2, , , , ,1 ! 選擇第17和第19號線并劃分為2段FLST,5,2,4,ORDE,2 FITEM,5,20 FITEM,5,24 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,_Y LESIZE,_Y1, , ,4, , , , ,1 !選擇第20和第24號線并劃分為4段TYPE, 1 !給鋼環(huán)附單元類型MAT, 1 !給鋼環(huán)附材料REAL, !給鋼環(huán)附實常數(shù) ESYS, 0 SECNUM, CM,_Y,VOLU VSEL, , , , 1 CM,_Y1,VOLU CHKMSH,VOLU CMSEL,S,_Y VSWEEP,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 !給銅滑塊劃分單元TYPE, 1 !給銅滑塊附單元類型MAT, 2 !給銅滑塊附材料REAL, !給銅滑塊附實常數(shù) ESYS, 0 SECNUM, CM,_Y,VOLU VSEL, , , , 2 CM,_Y1,VOLU CHKMSH,VOLU CMSEL,S,_Y VSWEEP,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 !給銅滑塊附實常數(shù) /COM, CONTACT PAIR CREATION - START CM,_NODECM,NODE CM,_ELEMCM,ELEM CM,_KPCM,KP CM,_LINECM,LINE CM,_AREACM,AREA CM,_VOLUCM,VOLU /GSAV,cwz,gsav,tempMP,MU,1,0.2 MAT,1 R,3 REAL,3 ET,2,170ET,3,174KEYOPT,3,9,0R,3,RMORE, RMORE,0RMORE,0 KEYOPT,3,1,3! Generate the target surface ASEL,S,4 CM,_TARGET,AREA TYPE,2 NSLA,S,1ESLN,S,0ESURF,ALL CMSEL,S,_ELEMCM ! Generate the contact surface ASEL,S,9 CM,_CONTACT,AREATYPE,3 NSLA,S,1ESLN,S,0ESURF,ALL ALLSEL ESEL,ALLESEL,S,TYPE,2 ESEL,A,TYPE,3 ESEL,R,REAL,3 /PSYMB,ESYS,1 /PNUM,TYPE,1/NUM,1 EPLOT 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,_AREACM CMSEL,S,_VOLUCM CMDEL,_VOLUCM /GRES,cwz,gsav CMDEL,_TARGET CMDEL,_CONTACT /COM, CONTACT PAIR CREATION - END !用接觸向導建立銅滑塊與鋼環(huán)間的接觸/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,1 FITEM,2,10 /GO SFA,P51X,1,PRES, %PRE% !在銅滑塊上表面施加壓力載荷TUNIF,20, !施加初始溫度FLST,2,1,5,ORDE,1 FITEM,2,3 /GO DA,P51X,UX, !在鋼環(huán)的外表面施加UX向約束FLST,2,1,5,ORDE,1 FITEM,2,3 /GO DA,P51X,UY, !在鋼環(huán)的外表面施加UY向約束FLST,2,1,5,ORDE,1 FITEM,2,3 /GO DA,P51X,UZ, !在鋼環(huán)的外表面施加UZ向約束ESEL,S,MAT,2 !選擇銅滑塊單元 NSLE !選擇銅滑塊節(jié)點CSYS,1 !激活柱坐標系FLST,2,75,1,ORDE,2 FITEM,2,876 FITEM,2,-950NROTAT,P51X !旋轉節(jié)點坐標系FLST,2,75,1,ORDE,2 FITEM,2,876 FITEM,2,-950 /GO D,P51X, ,0.01, , , ,UY, , , , , ! 給銅滑塊加Y向位移載荷（即切向位移）ALLSEL !選擇所有的節(jié)點，單元，幾何體FINISH /SOL !進入求解器ANTYPE,4 !設置為瞬態(tài)分析 TRNOPT,FULL !選擇瞬態(tài)全分析方法LUMPM,0 !關閉集中質量矩陣計算格式NLGEOM,1 !定義大變形NSUBST,100,200,100 !定義載荷步OUTRES,ALL,20 !定義結果輸出子步間隔AUTOTS,1 !打開時間開關KBC,0 !設置為斜坡載荷 TIME,10 !定義求解時間SSTIF,0 !關閉應力剛化影響NROPT,FULL, , !選擇瞬態(tài)全分析方法EQSLV, , ,0,PRECISION,0 !定義雙精度求解MSAVE,0 TOFFST,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,0 MPDATA,DENS,1,7800 !定義材料1的密度MPTEMP, MPTEMP,1,0 MPDATA,C,1,460 !定義材料1的比熱MPTEMP, MPTEMP,1,0 MPDATA,KXX,1,66.6 !定義材料1的熱傳導系數(shù)MPTEMP, MPTEMP,1,0 UIMP,1,REFT,30 !定義材料1的熱膨脹系數(shù)的參考溫度MPDATA,ALPX,1,1.06e-5 !定義材料1的熱膨脹系數(shù)MPTEMP, MPTEMP,1,0 MPDATA,EX,1,206e9 !定義材料1的彈性模量 MPDATA,PRXY,1,0.3 !定義材料1的泊松比MPTEMP, MPTEMP,1,0 MPDATA,DENS,2,8900 !定義材料2的密度MPTEMP, MPTEMP,1,0 MPDATA,C,2,390 !定義材料2的比熱MPTEMP, MPTEMP,1,0 MPDATA,KXX,2,383 !定義材料2的熱傳導系數(shù)MPTEMP, MPTEMP,1,0 UIMP,2,REFT,30 !定義材料2的熱膨脹系數(shù)的參考溫度MPDATA,ALPX,2,1.75e-5 !定義材料2的熱膨脹系數(shù)MPTEMP, MPTEMP,1,0 MPDATA,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, 1 MAT, 1REAL, !給鋼固定塊附材料屬性ESYS, 0 SECNUM, ESIZE,0.25,0, !設置單元劃分尺寸MSHAPE,0,2D MSHKEY,1CM,_Y,AREA ASEL, , , , 1 CM,_Y1,AREA CHKMSH,AREA CMSEL,S,_Y AMESH,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 !給鋼固定塊劃分單元TYPE, 1 MAT, 2REAL, !給銅滑塊附材料屬性ESYS, 0 SECNUM, CM,_Y,AREA ASEL, , , , 2 CM,_Y1,AREA CHKMSH,AREA CMSEL,S,_Y AMESH,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 !給銅滑塊劃分單元CM,_NODECM,NODE CM,_ELEMCM,ELEM CM,_KPCM,KP CM,_LINECM,LINE CM,_AREACM,AREA CM,_VOLUCM,VOLU MP,MU,1,0.2 MAT,1 MP,EMIS,1, R,3 REAL,3 ET,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.0KEYOPT,3,3,0KEYOPT,3,4,0KEYOPT,3,5,0KEYOPT,3,7,0KEYOPT,3,8,0KEYOPT,3,9,1KEYOPT,3,10,1 KEYOPT,3,11,0 KEYOPT,3,12,0 KEYOPT,3,2,0KEYOPT,3,1,1LSEL,S,3 CM,_TARGET,LINE TYPE,2 NSLL,S,1ESLN,S,0ESURF,ALL CMSEL,S,_ELEMCM LSEL,S,5 CM,_CONTACT,LINETYPE,3 NSLL,S,1ESLN,S,0ESURF,ALL ALLSEL ESEL,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,_AREACM CMSEL,S,_VOLUCM CMDEL,_VOLUCM CMDEL,_TARGET CMDEL,_CONTACT FLST,2,1,5,ORDE,1 FITEM,2,1 !用接觸向導建立滑塊和固定塊間的接觸面/GO DA,P51X,UX, !將下面的鋼塊施加X向約束 FLST,2,1,5,ORDE,1 FITEM,2,1 /GO DA,P51X,UY, !將下面的鋼塊施加Y向約束 FLST,2,1,4,ORDE,1 FITEM,2,6 /GO DL,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,1 FITEM,2,7 /GO SFL,P51X,PRES, %PRE% !在銅滑塊上表面施加壓力載荷TUNIF,30, !給銅塊和鋼塊施加初始溫度ALLSEL !選擇所有的節(jié)點，單元，幾何體ANTYPE,TRANS !設置為瞬態(tài)分析 TIMINT,OFF,STRUC !關閉結構的動態(tài)慣性影響TINTP,1.0 !定義瞬態(tài)積分參數(shù)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,RESIZE WPSTYLE,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,1 MP,RSVX,1,17.2e-9 !MP,MURX,2,500 MP,MURX,3,1 BLOCK,-w,w,-h,h,-l,l ! ! BLOCK,-w-a,w+a,h+a,h+2*a,-l,l BLOCK,w+a,w+a+b,h+a,h+2*a,-l,l BLOCK,-w-a,-w-a-b,h+a,h+2*a,-l,lBLOCK,w+a,w+a+b,h,h+a,-l,l BLOCK,-w-a,-w-a-b,h,h+a,-l,l! BLOCK,-w-a,w+a,-(h+a),-(h+2*a),-l,l BLOCK,w+a,w+a+b,-(h+a),-(h+2*a),-l,lBLOCK,-w-a,-w-a-b,-(h+a),-(h+2*a),-l,l BLOCK,w+a,w+a+b,-h,-(h+a),-l,l BLOCK,-w-a,-w-a-b,-h,-(h+a),-l,l! BLOCK,-22e-6,22e-6,-15e-6,15e-6,-15e-6,15e-6VOVLAP,ALL NUMCMP,VOLU VSEL,S,VOLU,1 VATT,1VSEL,S,VOLU,2,11,1VATT,2 VSEL,S,VOLU,12 VATT,3 VSEL,S,MAT,1,2,1 ASLV,S LSLA,S LESIZE,ALL,3e-6 MSHKEY,1 MSHAPE,0,3D VMESH,ALL VSEL,S,MAT,3 ASLV,S LSLA,S LESIZE,ALL,3e-6 MSHKEY,0 MSHAPE,1,3D VMESH,ALL ALLSEL,ALL FINISH /SOLU ESEL,S,MAT,1 BFE,ALL,JS,1,-1e9ASEL,S,EXT DA,ALL,AZ,0 ALLSEL,ALL ANTYPE,0NROPT,AUTO,EQSLV,SPAR,0, PRECISION,0 MSAVE,0 PIVCHECK,1 /STATUS,SOLUSOLVE FINISH /POST1 ! elements /VSCALE,1,1,0 ! !* PLVECT,B, , , ,VECT,ELEM,ON,0 ESEL,U,MAT,3 /REPLOT /color/size/b/alignbcolor=#0000ff/color/balign=leftcode/BATCH,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 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 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,1VATT,1 ! conductorVSEL,S,VOLU,2,11,1 ! magnetic coreVATT,2VSEL,S,VOLU,12 ! air ambientVATT,3! Mesh volumesVSEL,S,MAT,1,2,1 ! Conductor and magnetic core with hexaelementsASLV,SLSLA,SLESIZE,ALL,3e-6MSHKEY,1 ! MAPPED VOLUME MESHMSHAPE,0,3D ! USING HEXVMESH,ALLVSEL,S,MAT,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,MAT,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,0MSAVE,0PIVCHECK,1/STATUS,SOLUSOLVE/POST1!* Vectorplot of magnetic flux (BSUM) for magnetic cores and conductorselements */VSCALE,1,1,0!*PLVECT,B, , , ,VECT,ELEM,ON,0ESEL,U,MAT,3/REPLOTFINISH! /EXIT,MODEL/code/alignalign=centersize=4color=blueb=錯因告訴你了，稍微改改就可以得到正確的解，來改著玩玩吧，獎分哦（+1技術分，+50貢獻分）=/b/color/size/align3. 幫助感應加熱例子induction heating of a solid cylinder billetbatch,list/filenam, induc/config,nres,100000 ！將載荷步最大設置為100000/prep7shpp,off/title, induction 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 permeabilitympdata,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-6mpdata,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 workpiecekesize,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 size ？lsel,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 radiationet,3,0 ! reset type 3 to null element ？nsel,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 ther