二维数值传热学C++程序代码.docxVIP

#includestdio.h #includeiostream #includemath.h #includefstream //////////////////////////////////////////////////////////////////////////////// /******************************************************************************/ //////////////////////////////////////////////////////////////////////////////// 从此开始至下一个 /*/ 标记带处之间的代码可以写在头文件里。 //////////////////////////////////////////////////////////////////////////////// 计算参数宏定义 开始 //////////////////////////////////////////////////////////////////////////////// #define GRID_ROW 82 #define GRID_COLUMN 102 网格横纵数目，边界处多出一列处理 #define DX 0.05 #define DY 0.05 网格单元尺寸， m #define DEN_COPPER 8900 #define DEN_STEEL 7850 #define DEN_CONCRETE 2430 // 材料的密度， kg/(m^3) #define C_COPPER 390 #define C_STEEL 460 #define C_CONCRETE 970 // 材料比热容， J/(kg*K) #define COND_COPPER 400 #define COND_STEEL 20 #define COND_CONCRETE 1.63 材料导热系数， W/(m*K) #define CONV 5 // 空气与各种材料的对流换热系数， W/(m^2*K) #define INITL_TMPRT 298 初始温度， K #define T_C2K(x) ((x)+273)//C to K #define T_K2C(x) ((x)-273)//K to C 摄氏温度与绝对温度相互转换 #define T_AIR 298; 空气温度 ,K #define t_END 3600 计算时长， s #define t_STEP 1 时间步长， s #define COND_COP_CON 3.2467 #define COND_COP_STL 38.0952 #define COND_CON_STL 3.0143 // 不同材料的两个相邻结点，导热系数取调和平均数， W/(m*K) //////////////////////////////////////////////////////////////////////////////// 计算参数宏定义 结束 //////////////////////////////////////////////////////////////////////////////// enum Property{ copper, steel, concrete, iso834, adiabat, air }; 结点属性的枚举类型声明 struct Grid { double mNodeT[GRID_ROW][GRID_COLUMN];// 结点温度 Property mNodeProperty[GRID_ROW][GRID_COLUMN];// 结点材料属性 double mNodeC[GRID_ROW][GRID_COLUMN];//结点比热容 double mNodeDen[GRID_ROW][GRID_COLUMN];//结点密度 }; 结构体 Grid 定义，此为计算的主要承载数据结构，将在计算中存储温度数据 函数全局声明开始 void SetProperty(Grid );// void InitlTMPRT(Grid );// void SetNodePara(Grid );// void CalculateTMPRT(Grid );//  属性设置函数 温度初始化函数 结点物性参数设置函数 温度计算函数，几乎全部数值计算的承担者 void write(int,Grid );// 写文件函数，每隔 100 秒输出温度分布文件到当前目录 函数全局说明结束 void SetProperty(Grid A) {