ansys moving heat flux.zip_ansys heat_ansys moving heat_ansys 焊接

!com, units: SI, m - degree C - W - kg -sec - J ! ! Bi=h*dx/K = 60*0.002/51.9 = 0.0023 < 1 ! F0 = ( Kdt ) / (rho*C*dx^2) ! assume F0= beta = 0.4 ! dt = (beta*rho*C*dx^2)/K = 0.1179 (sec) ! /title,transient analysis 2 (moving heat flux) /filname,transient-2 ?????????????????? ! /triad,rbot ??????????????????????? ! /prep7 ????????????? ! wpstyl,,,,,,,,1 ???????????? wpstyl,0.01,0.01,0,0.1,0.003,0,1,,5 ???????????宏???? ! blc4,0,0.1,0.1,-0.1 !0.1*0.1 m ??????????????? cyl4,0.03,0.05,0.01 ???????????? asba, 1, 2 !drill a hole ! blc4,0.1,0.1,-0.01,-0.1 ???????????????? ! flst,2,2,5,orde,2 ???GUI?????????????? fitem,2,1 ???GUI???????????????? fitem,2,3 ???GUI???????????????? aovlap,p51x ???????????????? ! wpstyle,,,,,,,,0 ???????? aplot ????????????? ! et,1,plane55 ?????��??????plane55 !* keyopt,1,1,0 !????????????? keyopt,1,3,0 !2d plane problem keyopt,1,4,0 keyopt,1,8,0 keyopt,1,9,0 ! !aisi 1045 steel ! uimp,1,dens, , ,7870 , ????????????? !kg/m3 uimp,1,kxx, , , 51.9 , !w/m-c uimp,1,c, , , 486, !w-s / kg-c ! type,1 mat,1 ! esize,0.002, !element size = 0.002 m ! mshkey,1 ??????????????? amesh,1 ????????1 ! mshkey,0 ???????????????? amesh,2 ????????2 ! !boundary conditions ! /psf,hflux, ,2,0 ??????????????????? ! !convection ! flst,2,1,4,orde,1 fitem,2,4 sfl,p51x,conv,50, ,25, ! finish ! /solu ! antype,trans ???????????? trnopt,full ??????????????? nropt,auto, , ????????????????=???????? eqslv, ???????????????????? ! solcontrol,on ???????????????????? autots,on ?????????? ! kbc,1 !step outres,all,all, ????????快??? ! !load step 1, initial conditions 25?? ! time,0.001 ????????????? deltim,0.001,0.001,0.001 ????????????????? timint,off !time integration = off tunif,25, ??????梀???????? solve ????? ! !load step 2..... , apply moving heat flux ! timint,on !time integration = on ! j=1 ! *do,i,200,80,-5 ! time,j ????????????? deltim,0.11,0.11,0.11, ????????????????? ! sfedele,i+5,2,hflux !delete heat flux of previous step ! sfe,i,2,hflux, ,5/0.002/0.002, , , !apply heat flux, face 2 , 5 w ! eplot ! solve j=j+1 *enddo ! finish ! /post1 /efacet,1 ?????? set,last ???????????快??????? ! plnsol,temp, ,0, ???????????? finish