分析解.cpp

二维稳态导热问题。设材料物性均为常数

#include <iostream>
#include<stdio.h>
#include<math.h>
#define NM  50              //未知数个数
void main()
{
	int n,i,j,N;            //N为x方向结点个数
	double x,y,t[NM][NM],pi;//t[NM][NM]存放温度
	float dx;               //定义步长
	FILE *outfile;
	pi=3.141592654;
	printf("输入步长(dx/L):\n");
	scanf("%f",&dx);
	N=int(1.0/dx)+2;        //计算x方向结点个数
	for(i=0;i<N;i++)        //假设初始温度场
		for(j=0;j<N;j++)
			t[i][j]=0.0;
	for(i=0;i<N;i++)        //根据公式计算温度
		for(j=0;j<N;j++)
		{
			x=dx*i;
			y=dx*j;
			for(n=1;n<200;n++)
			t[i][j]=t[i][j]+2/pi*(pow(-1,n+1)+1)/n*sin(n*pi*x)*sinh(n*pi*y)/sinh(n*pi);
		}
    outfile=fopen("分析解.txt","w");
	if(outfile==NULL)
	{
		printf("打不开数据文件:分析解.txt!");
		exit(1);
	}
	fprintf(outfile,"步长=%f\n",dx);
    for(i=N-1;i>=0;i--)
	{
		for(j=N-1;j>=0;j--)
		{	printf("%5.4f ",t[i][j]);
	     	fprintf(outfile,"%5.4f ",fabs(t[i][j]));
		}
		printf("\n");
		fprintf(outfile,"\n");
	}
		
	fclose(outfile);	


}