肋片散热分析—计算传热学课程设计(3)

{

double Kd,Slx,Bic;

Kd=0.2625/(2.026025+0.2625*y[0]); Slx=2*LGZ*sqrt(3.14*LGZ/V); Bic=0.08*LGZ/4.0; Nr=0.000191*LGZ/4.0; f[0]=y[1]; //f[0]=dy1/dx

f[1]=-Kd*y[1]*y[1]+Slx*(Bic*y[0]+Nr*(pow((y[0]+3.908667),4)-pow(3.868667,4))); return; }

void pqw1(double Y,double *P,double *Q,double *W) //左边界处的第三类边界条件(x=xa)

//P1*y1+Q1*y2=W1-用户应根据具体条件进行修改 {

*P=1.0; *Q=0.0; *W=1; }

void pqw2(double Y,double *P,double *Q,double *W) //右边界处的第三类边界条件(x=xb)

//P2*y1+Q2*y2=W2-用户应根据具体条件进行修改 { *P=0; *Q=1; *W=0; return; } 求最大传热量

10

//调用各个函数求最大传热量 for(g=0;L-g*bc>0;g++) {

//用肋高函数求肋高

LG[k]=LeiGao(L,bc,k);

LGZ=LG[k];

x=xa;

//利用试射法确定m

shoot(N,NS,x,h,M1,M2,Eps,y0,y);

//求热量

Q[k]=RL(V, y0); LG[g]=LeiGao(L,bc,g);

LGZ=LG[g];

x=xa;

//利用试射法确定m

shoot(N,NS,x,h,M1,M2,Eps,y0,y); Q[g]=RL(V, y0);

//输出任意长度及所对应的热量

printf(\输出长度 fprintf(fp,\

printf(\输出肋基出温度梯度 fprintf(fp,\

printf(\输出对应传热量 fprintf(fp,\ //保存最大传热量

if(Q[k]

}

//输出最大传热量，并输出对应圆肋的尺寸 max=k;

LG[max]=LeiGao(L,bc,max);

11

LGZ=LG[max]; A=V/LG[max]; r=sqrt(A/3.14); CJB=LG[max]/r;

printf(\最佳肋高LG[max]=%f\\n\ fprintf(fp,\最佳肋高LG[max]=%f\\n\ printf(\最佳面积A=%f\\n\ fprintf(fp,\最佳面积A=%f\\n\ printf(\最佳半径r=%f\\n\ fprintf(fp,\最佳半径r=%f\\n\ printf(\最佳长径比CJB=%f\\n\ fprintf(fp,\最佳长径比CJB=%f\\n\ x=xa;

//利用试射法确定m

shoot(N,NS,x,h,M1,M2,Eps,y0,y); Q[max]=RL(V, y0);//求最大传热量

printf(\最大热量Q[max]=%6.4f \ fprintf(fp,\最大热量Q[max]=%6.4f \ printf(\ fprintf(fp,\

//输出最大传热量时的温度分分布

printf(\输出最大传热量时的温度分布\\n\显示在屏幕上 fprintf(fp,\输出最大传热量时的温度分布\\n\保存到文件中 //输出表头

printf(\显示在屏幕上 fprintf(fp,\保存到文件中 //输出x=a时的结果 printf(\ fprintf(fp,\ for(i=0;i

12

{

printf(\ fprintf(fp,\ }

printf(\ fprintf(fp,\ x=xa;

//调用R-K方法计算并输出后续各点的值 for(i=0;i

rungek(N,&x,h,y); //根据求出的m解决问题 printf(\ fprintf(fp,\ for(i=0;i

printf(\ fprintf(fp,\ }

printf(\ fprintf(fp,\ } 验证程序

//令lmd=常数，无辐射，验证程序的正确性 if(kd==0.0&&Nr==0.0) {

printf(\令lmd=100，无辐射，验证程序理论温度\\n\ fprintf(fp,\令lmd=100，无辐射，验证程序理论温度\\n\ MM=sqrt(0.16/r);

13

MMH=MM*LGZ;

printf(\显示在屏幕上 fprintf(fp,\保存到文件中 x=0.0000;

for(i=0;h*i<=1.0;i++) { x=h*i;

printf(\ fprintf(fp,\

LLWD=(exp(MMH*(x-1))+exp(MMH*(1-x)))/(exp(MMH)+exp(-1.0*MMH));//温度分布

printf(\ fprintf(fp,\ printf(\ fprintf(fp,\ }

}

附录2 数学模型的无量纲化过程推导

针对式（1）进行无量纲化处理，为此定义

x?xL??，

T?TfTw?Tf （10）

其中Tf、Tw均为常数，假定λ=λοk(T)，则无量纲化过程如下。

2d?dT?dTd?dTdTd??dT??A??A?A ?A??A???? （11）22dx?dx?dxdxdxdxdx?dx?22

dTdx?d???Tw?T?dxLf??Tf?????Tw?Tfd?Ldx （12-a）

dTdx222d?dT?Tw?Tfd?????22dx?dx?Ldx (12-b)

14