/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)
/filnam,transient-2
/triad,rbot
/prep7
wpstyle,,,,,,,,1
wpstyle,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
fitem,2,1
fitem,2,3
aovlap,p51x
wpstyle,,,,,,,,0
aplot
et,1,plane55
keyopt,1,1,0
keyopt,1,3,0 !2d plane problem
keyopt,1,4,0
keyopt,1,8,0 295
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
mshkey,0
amesh,2
!boundary conditions
/psf,hflu, ,2,0
convection
flst,2,1,4,orde,1
fitem,2,4
sfl,p51x,conv,50, ,25,
finish
/solu
antype,trans
trnopt,full
nropt
移动温度荷载计算
最新推荐文章于 2025-02-01 00:05:09 发布