二维数值传热学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){for (int j = 0; j GRID_COLUMN;++j)A.mNodeProperty[0][j] = air;//第一行为空气边界，计算区