偏微分方程組解法-文庫吧資料
2025-06-25 21:46本頁面
【正文】 %60、1180、2300、360、420s時(shí)刻溫度分布圖plot(x,sol(2,:,1)) 60s(t網(wǎng)格為0:60:2160,其時(shí)間間隔為0,60,120,180,……2160,第二點(diǎn)為60s)hold onplot(x,sol(3,:,1))plot(x,sol(4,:,1))plot(x,sol(5,:,1))plot(x,sol(6,:,1))plot(x,sol(7,:,1))plot(x,sol(8,:,1))%1080、141800、2160s時(shí)刻溫度分布圖plot(x,sol(19,:,1))hold onplot(x,sol(25,:,1))plot(x,sol(31,:,1))plot(x,sol(37,:,1))600s內(nèi)的溫度分布變化2160s內(nèi)的溫度分布變化不同時(shí)刻溫度分布圖，根據(jù)以上兩準(zhǔn)數(shù)可知：該傳熱過程內(nèi)外對(duì)流及導(dǎo)熱阻力相當(dāng)；當(dāng)加熱時(shí)間小于360s時(shí)為非正規(guī)階段，加熱時(shí)間大于360s后進(jìn)入正規(guī)階段，從上圖也可得到該結(jié)論。mesh(x,t,sol)%2160s內(nèi)的溫度分布變化t=[0:60:2160]。t=[0:60:600]。qr=。ql=1。%偏微分方程初始條件(一維動(dòng)態(tài)平壁兩側(cè)對(duì)流)function u0=icbun1(x)u0=20。f=dudx。 某厚度為20cm鋼板原溫度為20，現(xiàn)將其置于1000的爐中加熱，平壁導(dǎo)熱系數(shù)為，對(duì)流傳熱系數(shù)，導(dǎo)溫系數(shù)為；試分析溫度分布隨時(shí)間的變化及鋼板表面溫度達(dá)到500時(shí)所需的時(shí)間。sol=pdepe(0,pdefun,icfun,bcfun,x,t)。命令：x=linspace(0,10,20)*1e2。pr=ur20。%偏微分方程邊界條件(一維動(dòng)態(tài)傳熱)function [pl,ql,pr,qr]=bcfun(xl,ul,xr,ur,t)pl=ul120。s=0。 解： 模型轉(zhuǎn)化為標(biāo)準(zhǔn)形式：初始條件為：
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1