这是一个加表面效应单元的例子,你看一下有没有帮助。 <BR>FINISH <BR>/CLEAR <BR>!!!gauss加载 <BR>/title, the simulate of laser cladding process <BR>*SET,P,2000, !the power of the laser <BR>*SET,V,0.008 !scanning speed <BR>!*SET, RADIUS,0.003 !radius of the laser beam <BR>*SET, R,0.003 <BR>*set,k,0.6 !the coefficient of the absorb <BR>*SET,L,0.05 !the length of the block <BR>/prep7 <BR>!熔覆层高5mm=0.005m <BR>!基体高10mm=0.01m,长x=50mm=0.05,宽=20mm=0.02m <BR>BLOCK,0,0.05,0, -0.0105,-0.01,0, <BR>Vsel,all <BR>vplot <BR>! <BR>/USER, 1 <BR>/VIEW, 1, -0.198026943995E-01, 0.382677277552 , 0.923669829831 <BR>/ANG, 1, -0.107453019382 <BR>/REPLO <BR>! <BR>ET,1,SOLID70 <BR>ET,2,Surf152 <BR>DOF,TEMP <BR>KEYOPT,2,3,0 <BR>KEYOPT,2,4,1 <BR>KEYOPT,2,5,0 <BR>KEYOPT,2,6,0 <BR>KEYOPT,2,7,0 <BR>KEYOPT,2,8,4 !根据固体表面温度计算对流系数 <BR>KEYOPT,2,9,0 <BR>!用温度表方式定义材料性能 <BR>Mp,kxx,1,82.9 <BR>Mp,dens,1,8900 <BR>Mp,c,1,471 <BR>mptemp,1,20,300,600,750,900,1200,1500,1800,2000 <BR>mpdata,hf,1,1,6,50,120,180,200,250,378,700,850 <BR>! <BR>! <BR>Mshkey,0 !自由网格化 <BR>Mshape,1,3d <BR>type,1 <BR>esize,0.002, <BR>Vmesh,all <BR>eplot <BR>! <BR>! <BR>asel,s,loc,y,0 !对流边界条件要施加在5个面的周围 <BR>asel,a,loc,z,-0.01 <BR>asel,a,loc,x,0 <BR>asel,a,loc,x,0.05 <BR>type,2 !定义属性 <BR>amesh,all <BR>esel,s,type,,2 !选择表面效应单元 <BR>sf,all,conv,-1,25 !确定对流边界条件,表格加载 <BR>!sf,all,conv,200,25 !确定对流边界条件 <BR>!!底面也施加轻微的对流 <BR>allsel,all <BR>asel,s,loc,y, -0.0105 <BR>type,2 !定义属性 <BR>amesh,all <BR>esel,s,type,,2 <BR>sf,all,conv,3,25 <BR>allsel,all <BR>! <BR>nsel,all <BR>/psf,hflux,,2,0 !以箭头方式显示载荷 <BR>! <BR>finish <BR>/SOLU <BR>!!!!! <BR>antype,transient,new <BR>TRNOPT,FULL <BR>Nropt,,,off !关闭自适应下降 <BR>LUMPM,0 <BR>Solcontrol,on <BR>Autots,on <BR>Kbc,1 <BR>Outres,all,last, <BR>Lnsrch,on <BR>neqit,50 !每一子步中方程的迭代次数限值 <BR>!load step1,initial conditions 25 <BR>time,0.001 <BR>deltim,0.001,0.001,0.001 <BR>allsel,all <BR>timint,off <BR>tunif,25 !均匀的节点初始温度25,只对分析的第一个子步有效 <BR>solve <BR>EPLOT <BR>! <BR>! load step2………,apply moving heat flux <BR>allsel <BR>nsel,all <BR>t=0 <BR>dt=0.5 <BR>local,11,1, 0.001+0.005,0,0,0,-90,0 <BR>CSYS,11 <BR>esel,s,type,,1 <BR>nsel,r,loc,x, 0,0.003 <BR>nsel,r,loc,y,0,180 <BR>nsel,r,loc,z,0 <BR>sf,,hflux,7000000 <BR>!sf,,hflux,%flux2% <BR>csys,0 <BR>t=t+dt <BR>time,t !每一载荷步结束时的时间 <BR>deltim,0.1,0.05,0.1,on <BR>kbc,1 !设定为阶越载荷 <BR>!neqit,25 !每个子步的默认次数 <BR>autots,on !打开自动时间步长 <BR>timint,on,therm !时间积分效果 <BR>TINTP,,,,1 <BR>Lnsrch,on <BR>!outres,all,all <BR>outres,all,last <BR>!outres,all,10 <BR>allsel <BR>solve <BR>save |