表面效应单元-单元质量不合格的原因分析
本帖最后由 wxwieandy 于 2019-7-1 02:03 编辑使用surf151表面效应单元，含有外部附加节点，在模拟感应加热求解计算时提示单元质量错误，很疑惑，为什么会出现这样的错误提示呢？（单元质量不合格），命令流如下：思考了3天，始终没有找到表面效应单元质量不合格原因，希望大家指点迷津，谢谢！
!------------感应线圈截面参数20×20mm-10×10mm/filenam, induc/prep7!shpp,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,0row_in=0.0065 ! inner radius of copper conductorrow_out=0.014 ! out radius of copper conductorsleeve_out=0.015 ! out radius of sleeveric=0.019 ! inner radius of induction coilroc=0.039 ! outer radius of induction coilro=0.125 d_coil=0.014 h_coil=0.02 t=0.08 h_sleeve=0.01 freq=1500 ! frequency (Hz.)pi=4*atan(1) ! picond=5.71e7 muzero=4e-7*pi ! free-space permeability murx=1 skind=sqrt(1/(pi*freq*cond*muzero*murx)) ! skin depthftime=4 ! final time tinc=2 ! time increment for harmonic analysistime=0 ! initialize timedelt=0.05 ! maximum delta time stepemunit,mks ! set magnetic unitsmp,murx,1,1 ! air relative permeabilitymp,murx,3,1 ! coil relative permeabilitymp,murx,2,1 !copper relative permeabilitymptemp,1,20,100,200,250,500,600 ! temps for resistivitymptemp,7,750,900,1000mpdata,rsvx,2,1,1.78e-8,2.33e-8,2.92e-8,3.35e-8,5.386e-8,6.25e-8mpdata,rsvx,2,7,7.41e-8,8.75e-8,9.72e-8 ! copper resistivityrectng,0,row_in,0,t/2 rectng,row_in,row_out,0,t/2 ! 铜管rectng,row_out, ric,0,t/2 !线圈和铜管之间的空气场rectng,ric,roc,d_coil/2,d_coil/2+h_coil ! 实心线圈截面rectng,ric+h_coil/4,roc-h_coil/4,d_coil/2+h_coil/4,d_coil/2+0.75*h_coil !线圈中的空心截面rectng,0,ro,0,t/2 ! outer air
allsaptn,all !面布尔运算allsarsym,y,all,,,,0,0!创建对称面nummrg,kpnumcmp,kp lsel,s,length,,row_out-row_inasll,saatt,2,1,1 !为铜管分配材料编号2,实常数集1和单元类型1cm,copper,area lsel,s,length,,0.5*h_coilasll,s,1cm,air_1,area lsel,s,length,,h_coillsel,a,length,,0.5*h_coilasll,s,1asel,u,,,air_1aatt,3,1,2 cm,coil,area allsasel,u,,,copperasel,u,,,coilaatt,1,1,2 cm,air,area !创建空气面组件 allslesize,all,0.001mshape,0,2damesh,copperamesh,coilamesh,airnummrg,nodenumcmp,noden !为surf151创建空间节点
*get,nmax,node,,num,maxtype,3 !激活第三种单元类型real,3 !实常数集编号为3mat,2ksel,s,loc,x,row_outlslk,s,1lmesh,allesll,semodif,all,3,nmax
et,3,0 ! reset type 3 to null elementalls/solunsel,s,loc,xnsel,u,,,nmax d,all,az,0 ! apply flux-normal b.c.allsesel,s,mat,,3bfe,all,js,,,,33.3e6 ! apply current density to coil allsantyp,harmharfrq, freq !感应加热的频率physics,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 regionet,3,151,,,1,1,1 ! SURF151 element for radiationkeyopt,3,9,1r,3,1,5.67e-8 ! form factor, Stefan-Boltzmann constantmp,dens,2,8.9e3 mptempmptemp,1,0,20,100,300,400,600 mptemp,7,800,1000,1100mpdata,kxx,2,1,401,398,393,384,379,366mpdata,kxx,2,7,352,336,170 mptempmptemp,1,0,20,100,200,600,1000 mptemp,7,1100mpdata,c,2,1,0.34e3,0.38e3,0.39e3,0.4e3,0.44e3,0.48e3mpdata,c,2,7,0.5e3 mp,emis,2,.60 finish/soluallsantype,transtoffst,273 tunif,25 ! initial uniform temperatured,nmax,temp,25 ! ambient temperaturecnvtol,heat,1 ! convergence tolerancekbc,1 ! step loadstrnopt,fullautots,on ! auto time-steppingdeltim,1e-3,1e-4,delt,on ! time step controloutres,all,all ! save all load step informationphysics,write,thermal ! write thermal physics filefinish alls*do,i,1,ftime/tinc ! solution *do looptime=time+tinc ! increment timephysics,read,emag ! read emag physics file/solu*if,i,eq,1,then tunif,25 ! initial temperature*else ldread,temp,last,,,,,rth ! read thermal analysis temperatures*endifsolve ! solve harmonic analysisfinish !------------------------------------------physics,read,thermal ! read thermal physics file/assign,esav,therm,esav ! redirect files for use in thermal restart/assign,emat,therm,emat/soluparsav,scalar,parameter,sav !save parameters before multiframe restart*if,i,gt,1,then antype,trans,rest ! thermal restart*endifparres,new,parameter,sav !restore parameters after multiframe restarttime,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 centerlinensol,3,2,temp,,tempsurf !store temperature at billet outer diameterplvar,2,3 ! plot temperature rise over timeprvar,2,3
页:
[1]