ansys旋转经典命令流

1 旋转摩擦12. 电磁三d命令流实例（论坛看到）113. 帮助感应加热例子induction heating of a solid cylinder billet154. 感应加热温度场的数值模拟（论文）inducheat30命令流195. 如何施加恒定的角速度？Simwe仿真论坛246. 旋转一个已经生成好的物体277. 产生这样的磁力线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的热传导系数MPTEMP, MPTEMP,1,0 UIMP,1,REFT,30 !定义材料1的热膨胀系数的参考温度MPDATA,ALPX,1,1.06e-5 !定义材料1的热膨胀系数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的热传导系数MPTEMP, MPTEMP,1,0 UIMP,2,REFT,30 !定义材料2的热膨胀系数的参考温度MPDATA,ALPX,2,1.75e-5 !定义材料2的热膨胀系数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, !建立钢环的几何模型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 !给钢环附单元类型MAT, 1 !给钢环附材料REAL, !给钢环附实常数 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, !给铜滑块附实常数 ESYS, 0 SECNUM, CM,_Y,VOLU VSEL, , , , 2 CM,_Y1,VOLU CHKMSH,VOLU CMSEL,S,_Y VSWEEP,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 !给铜滑块附实常数 /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 !用接触向导建立铜滑块与钢环间的接触/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, !在钢环的外表面施加UX向约束FLST,2,1,5,ORDE,1 FITEM,2,3 /GO DA,P51X,UY, !在钢环的外表面施加UY向约束FLST,2,1,5,ORDE,1 FITEM,2,3 /GO DA,P51X,UZ, !在钢环的外表面施加UZ向约束ESEL,S,MAT,2 !选择铜滑块单元 NSLE !选择铜滑块节点CSYS,1 !激活柱坐标系FLST,2,75,1,ORDE,2 FITEM,2,876 FITEM,2,-950NROTAT,P51X !旋转节点坐标系FLST,2,75,1,ORDE,2 FITEM,2,876 FITEM,2,-950 /GO D,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 !定义双精度求解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的热传导系数MPTEMP, MPTEMP,1,0 UIMP,1,REFT,30 !定义材料1的热膨胀系数的参考温度MPDATA,ALPX,1,1.06e-5 !定义材料1的热膨胀系数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的热传导系数MPTEMP, MPTEMP,1,0 UIMP,2,REFT,30 !定义材料2的热膨胀系数的参考温度MPDATA,ALPX,2,1.75e-5 !定义材料2的热膨胀系数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 !选择所有的节点，单元，几何体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,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