二維數(shù)值傳熱學(xué)C程序代碼

維數(shù)值傳熱學(xué)代碼檔編制序號LLETD298-POI08]維數(shù)值傳熱學(xué)代碼檔編制序號LLETD298-POI08]#include<>#include<iostream>#include<>#include<fstream>/////////////////////////////////////////////////////////////////////////////////""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""http://個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)個(gè)/////////////////////////////////////////////////////////////////////////////////〃從此開始至下一個(gè)/*/標(biāo)記帶處之間的代碼可以寫在頭文件里。////////////////////////////////////////////////////////////////////////////////計(jì)算參數(shù)宏定義開始////////////////////////////////////////////////////////////////////////////////#defineGRID_ROW82#defineGRID_COLUMN102〃網(wǎng)格橫縱數(shù)目，邊界處多出一列處理#defineDX#defineDY〃網(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#defineCOND_CONCRETE〃材料導(dǎo)熱系數(shù)，W/(m*K)#defineCONV5//空氣與各種材料的對流換熱系數(shù)，W/(m"2*K)#defineINITL_TMPRT298//初始溫度，K#defineT_C2K(x)((x)+273)//CtoK#defineT_K2C(x)((x)-273)//KtoC〃攝氏溫度與絕對溫度相互轉(zhuǎn)換#defineT_AIR298;//空氣溫度,K#definet_END3600//計(jì)算時(shí)長，s#definet_STEP1//時(shí)間步長，s#defineCOND_COP_CON#defineCOND_COP_STL#defineCOND_CON_STL〃不同材料的兩個(gè)相鄰結(jié)點(diǎn)，導(dǎo)熱系數(shù)取調(diào)和平均數(shù)，W/(m*K)////////////////////////////////////////////////////////////////////////////////計(jì)算參數(shù)宏定義結(jié)束////////////////////////////////////////////////////////////////////////////////enumProperty{copper,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ì)算的主要承載數(shù)據(jù)結(jié)構(gòu)，將在計(jì)算中存儲(chǔ)溫度數(shù)據(jù)〃函數(shù)全局聲明開始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ù)全局說明結(jié)束voidSetProperty(Grid&A){for(intj=0;j<GRID_COLUMN;++j)[0][j]=air;〃第一行為空氣邊界，計(jì)算區(qū)域外for(intj=0;j<GRID_COLUMN;++j)[GRID_ROW-1][j]=adiabat;//最后一行為絕熱邊界，計(jì)算區(qū)域外for(inti=0;i<GRID_ROW;++i)if((i>=0&&i<=40)||(i>=61&&i<GRID_ROW))[i][0]=air;else[i][0]=iso834;//第一列的空氣邊界和ISO834邊界，計(jì)算區(qū)域外for(inti=0;i<GRID_ROW;++i)[i][GRID_COLUMN-1]=air;//最后一列，空氣邊界，計(jì)算區(qū)域外for(inti=1;i<GRID_ROW-1;++i)for(intj=1;j<GRID_COLUMN-1;++j){if((i>=1&&i<=44&&j>=1&&j<=20)||(i>=61&&i<GRID_ROW-2&&j>=1&&j<=(100-(80-i)*(49/19))))[i][j]=concrete;〃混凝土區(qū)域elseif((i-20)*(i-20)+(j-25)*(j-25)<=20*20)[i][j]=air;//圓孔空氣區(qū)域elseif(i>=41&&i<=60&&j>=1&&j<=90&&((i-20)*(i-20)+(j-25)*(j-25)>=20*20))[i][j]=copper;//金屬銅區(qū)域的一部分elseif(i>=1&&i<+40&&j>=21&&j<=90)[i][j]=copper;〃金屬銅區(qū)域的另一部分else[i][j]=steel;〃剩下的即為金屬鐵的區(qū)域}}//屬性設(shè)置函數(shù)，將Grid結(jié)構(gòu)體A的所有網(wǎng)格賦予材料屬性(枚舉類型Property)voidInitlTMPRT(Grid&A){for(inti=0;i<GRID_ROW;++i)for(intj=0;j<GRID_COLUMN;++j)[i][j]=INITL_TMPRT;}//溫度初始化函數(shù)，將Grid結(jié)構(gòu)體A的全部結(jié)點(diǎn)溫度設(shè)置為初始溫度INITL_TMPRTvoidCalculateTMPRT(Grid&A){doubleQE=0,QW=0,QN=0,QS=0;//四個(gè)方向傳來的熱量for(inti=1;i<GRID_ROW-1;++i)for(intj=1;j<GRID_COLUMN-1;++j){switch[i][j]){{switch[i][j+1])//判斷右側(cè)屬性，并計(jì)算QE{casecopper:QE=[i][j+1]-[i][j])*COND_COPPER;break;caseair:QE=[i][j+1]-[i][j])*CONV*DY;break;casesteel:QE=[i][j+1]-[i][j])*COND_COP_STL;break;caseconcrete:QE=[i][j+1]-[i][j])*COND_COP_CON;break;caseiso834:QE=2*[i][j+1]-[i][j])*COND_COPPER;break;caseadiabat:QE=0;break;default:std::cout<<i<<”,"<<j+1<<"結(jié)點(diǎn)屬性未指定'n";getchar();}//std::cout<<QE<<"\n";switch[i][j-1])//判斷左側(cè)屬性，并計(jì)算QW{casecopper:QW=[i][j-1]-[i][j])*COND_COPPER;break;caseair:QW=[i][j-1]-[i][j])*CONV*DY;break;casesteel:QW=[i][j-1]-[i][j])*COND_COP_STL;break;caseconcrete:QW=[i][j-1]-[i][j])*COND_COP_CON;break;caseiso834:QW=2*[i][j-1]-[i][j])*COND_COPPER;break;caseadiabat:QW=0;break;default:{std::cout<<i<<","<<j-1<<"結(jié)點(diǎn)屬性未指定\n";getchar();}}switch[i-1][j])//判斷上側(cè)屬性，并計(jì)算QN{casecopper:QN=[i-1][j]-[i][j])*COND_COPPER;break;caseair:QN=[i-1][j]-[i][j])*CONV*DX;break;casesteel:QN=[i-1][j]-[i][j])*COND_COP_STL;break;caseconcrete:QN=[i-1][j]-[i][j])*COND_COP_CON;break;caseiso834:QN=2*[i-1][j]-[i][j])*COND_COPPER;break;caseadiabat:QN=0;break;default:std::cout<<i-1<<”,"<<j<<"結(jié)點(diǎn)屬性未指定\n”;getchar();}switch[i+1][j])//判斷下側(cè)屬性，并計(jì)算QS{casecopper:QS=[i+1][j]-[i][j])*COND_COPPER;break;caseair:QS=[i+1][j]-[i][j])*CONV*DX;break;casesteel:QS=[i+1][j]-[i][j])*COND_COP_STL;break;caseconcrete:QS=[i+1][j]-[i][j])*COND_COP_CON;break;caseiso834:QS=2*[i+1][j]-[i][j])*COND_COPPER;break;caseadiabat:QS=0;break;default:std::cout<<i+1<<”,"<<j<<"結(jié)點(diǎn)屬性未指定\n”;getchar();}break;}{switch[i][j+1])//判斷右側(cè)屬性，并計(jì)算QE{casecopper:QE=[i][j+1]-[i][j])*COND_COP_STL;break;caseair:QE=[i][j+1]-[i][j])*CONV*DY;break;casesteel:QE=[i][j+1]-[i][j])*COND_STEEL;break;caseconcrete:QE=[i][j+1]-[i][j])*COND_CON_STL;break;caseiso834:QE=2*[i][j+1]-[i][j])*COND_STEEL;break;caseadiabat:QE=0;break;default:std::cout<<i<<”,"<<j+1<<"結(jié)點(diǎn)屬性未指定\n”;getchar();}switch[i][j-1])//判斷左側(cè)屬性，并計(jì)算QW{casecopper:QW=[i][j-1]-[i][j])*COND_COP_STL;break;caseair:QW=[i][j-1]-[i][j])*CONV*DY;break;casesteel:QW=[i][j-1]-[i][j])*COND_STEEL;break;caseconcrete:QW=[i][j-1]-[i][j])*COND_CON_STL;break;caseiso834:QW=2*[i][j-1]-[i][j])*COND_STEEL;break;caseadiabat:QW=0;break;default:std::cout<<i<<”,"<<j-1<<"結(jié)點(diǎn)屬性未指定\n”;getchar();}switch[i-1][j])//判斷上側(cè)屬性，并計(jì)算QN{casecopper:QN=[i-1][j]-[i][j])*COND_COP_STL;break;caseair:QN=[i-1][j]-[i][j])*CONV*DX;break;casesteel:QN=[i-1][j]-[i][j])*COND_STEEL;break;caseconcrete:QN=[i-1][j]-[i][j])*COND_CON_STL;break;caseiso834:QN=2*[i-1][j]-[i][j])*COND_STEEL;break;caseadiabat:QN=0;break;default:std::cout<<i-1<<”,"<<j<<"結(jié)點(diǎn)屬性未指定\n”;getchar();}switch[i+1][j])//判斷下側(cè)屬性，并計(jì)算QS{casecopper:QS=[i+1][j]-[i][j])*COND_COP_STL;break;caseair:QS=[i+1][j]-[i][j])*CONV*DX;break;casesteel:QS=[i+1][j]-[i][j])*COND_STEEL;break;caseconcrete:QS=[i+1][j]-[i][j])*COND_CON_STL;break;caseiso834:QS=2*[i+1][j]-[i][j])*COND_STEEL;break;caseadiabat:QS=0;break;default:std::cout<<i+1<<”,"<<j<<"結(jié)點(diǎn)屬性未指定\n”;getchar();}break;}{switch[i][j+1])//判斷右側(cè)屬性，并計(jì)算QE{casecopper:QE=[i][j+1]-[i][j])*COND_COP_CON;break;caseair:QE=[i][j+1]-[i][j])*CONV*DY;break;casesteel:QE=[i][j+1]-[i][j])*COND_CON_STL;break;caseconcrete:QE=[i][j+1]-[i][j])*COND_CONCRETE;break;caseiso834:break;caseadiabat:QE=0;break;default:std::cout<<i<<”,"<<j+1<<"結(jié)點(diǎn)屬性未指定'n";getchar();}switch[i][j-1])//判斷左側(cè)屬性，并計(jì)算QW{casecopper:break;caseair:QW=[i][j-1]-[i][j])*CONV*DY;break;casesteel:break;caseconcrete:QW=[i][j-1]-[i][j])*COND_CONCRETE;break;caseiso834:break;caseadiabat:QW=0;break;default:std::cout<<i<<”,"<<j-1<<"結(jié)點(diǎn)屬性未指定\n”;getchar();}switch[i-1][j])//判斷上側(cè)屬性，并計(jì)算QN{casecopper:QN=[i-1][j]-[i][j])*COND_COP_CON;break;caseair:QN=[i-1][j]-[i][j])*CONV*DX;break;casesteel:QN=[i-1][j]-[i][j])*COND_CON_STL;break;caseconcrete:QN=[i-1][j]-[i][j])*COND_CONCRETE;break;caseiso834:break;caseadiabat:QN=0;break;default:std::cout<<i-1<<”,"<<j<<"結(jié)點(diǎn)屬性未指定\n”;getchar();}switch[i+1][j])//判斷下側(cè)屬性，并計(jì)算QS{casecopper:QS=[i+1][j]-[i][j])*COND_COP_CON;break;caseair:QS=[i+1][j]-[i][j])*CONV*DX;break;casesteel:QS=[i+1][j]-[i][j])*COND_CON_STL;break;caseconcrete:QS=[i+1][j]-[i][j])*COND_CONCRETE;break;caseiso834:break;caseadiabat:QS=0;break;default:std::cout<<i+1<<”,"<<j<<"結(jié)點(diǎn)屬性未指定\n”;getchar();}break;}{[i][j]=T_AIR;break;}default:break;}[i][j]=(QE+QN+QS+QW)/[i][j]*[i][j]*DX*DY);}}〃函數(shù)會(huì)判斷每個(gè)結(jié)點(diǎn)相鄰4個(gè)結(jié)點(diǎn)的屬性，然后將四個(gè)結(jié)點(diǎn)傳來的熱量加和，〃最后根據(jù)能量守恒計(jì)算出中間結(jié)點(diǎn)的溫度值(時(shí)間全顯式)voidwrite(intt,Grid&A){charfilename[40];intx=1,y=1;sprintf_s(filename,"temperature%",t);std::ofstreamfout(filename);fout<<"X,Y,TEMPERATURE\n”;for(inti=GRID_ROW-2;i>0;--i){for(intj=1;j<GRID_COLUMN-1;++j){if(y>100)y=1;fout<<x<<"