資源描述:
《二維數(shù)值傳熱學(xué)C++程序代碼》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫(kù)。
1、#include#include#include#include/////////////////////////////////////////////////////////////////////////////////******************************************************************************////////////////////////////////////////
2、///////////////////////////////////////////從此開(kāi)始至下一個(gè)/*/標(biāo)記帶處之間的代碼可以寫在頭文件里。////////////////////////////////////////////////////////////////////////////////計(jì)算參數(shù)宏定義開(kāi)始////////////////////////////////////////////////////////////////////////////////#defineGRID_ROW82#defineGR
3、ID_COLUMN102//網(wǎng)格橫縱數(shù)目,邊界處多出一列處理#defineDX0.05#defineDY0.05//網(wǎng)格單元尺寸,m#defineDEN_COPPER8900#defineDEN_STEEL7850#defineDEN_CONCRETE2430//材料的密度,kg/(m^3)#defineC_COPPER390#defineC_STEEL460#defineC_CONCRETE970//材料比熱容,J/(kg*K)#defineCOND_COPPER400#defineCOND_STEEL20#defineCON
4、D_CONCRETE1.63//材料導(dǎo)熱系數(shù),W/(m*K)#defineCONV5//空氣與各種材料的對(duì)流換熱系數(shù),W/(m^2*K)#defineINITL_TMPRT298//初始溫度,K#defineT_C2K(x)((x)+273)//CtoK#defineT_K2C(x)((x)-273)//KtoC//攝氏溫度與絕對(duì)溫度相互轉(zhuǎn)換#defineT_AIR298;//空氣溫度,K#definet_END3600//計(jì)算時(shí)長(zhǎng),s#definet_STEP1//時(shí)間步長(zhǎng),s#defineCOND_COP_CON3.2467
5、#defineCOND_COP_STL38.0952#defineCOND_CON_STL3.0143//不同材料的兩個(gè)相鄰結(jié)點(diǎn),導(dǎo)熱系數(shù)取調(diào)和平均數(shù),W/(m*K)////////////////////////////////////////////////////////////////////////////////計(jì)算參數(shù)宏定義結(jié)束////////////////////////////////////////////////////////////////////////////////enumProperty{co
6、pper,steel,concrete,iso834,adiabat,air};//結(jié)點(diǎn)屬性的枚舉類型聲明structGrid{doublemNodeT[GRID_ROW][GRID_COLUMN];//結(jié)點(diǎn)溫度PropertymNodeProperty[GRID_ROW][GRID_COLUMN];//結(jié)點(diǎn)材料屬性doublemNodeC[GRID_ROW][GRID_COLUMN];//結(jié)點(diǎn)比熱容doublemNodeDen[GRID_ROW][GRID_COLUMN];//結(jié)點(diǎn)密度};//結(jié)構(gòu)體Grid定義,此為計(jì)算的主要
7、承載數(shù)據(jù)結(jié)構(gòu),將在計(jì)算中存儲(chǔ)溫度數(shù)據(jù)//函數(shù)全局聲明開(kāi)始voidSetProperty(Grid&);//屬性設(shè)置函數(shù)voidInitlTMPRT(Grid&);//溫度初始化函數(shù)voidSetNodePara(Grid&);//結(jié)點(diǎn)物性參數(shù)設(shè)置函數(shù)voidCalculateTMPRT(Grid&);//溫度計(jì)算函數(shù),幾乎全部數(shù)值計(jì)算的承擔(dān)者voidwrite(int,Grid&);//寫文件函數(shù),每隔100秒輸出溫度分布文件到當(dāng)前目錄//函數(shù)全局說(shuō)明結(jié)束voidSetProperty(Grid&A){for(intj=0;j<
8、GRID_COLUMN;++j)A.mNodeProperty[0][j]=air;//第一行為空氣邊界,計(jì)算區(qū)域外for(intj=0;j