logo资料库

稳态导热数值计算C代码.docx

第1页 / 共2页
第2页 / 共2页
资料共2页,全文预览结束
#include #include #define N 10 #define M 10 main() { char s; int i,j,l; float cha,x,y; float t[N][M],a[N][M]; /*打印出题目*/ printf("\t\t\t 二维稳态导热问题\t\t"); printf("\n\t\t\t\t\t\t----何鹏举\n"); printf("\n 题目:补充材料练习题二\n"); printf("\n 矩形区域,边长 L=W=1,假设区域内无内热源,导热系数为常熟,三个边温度为 T1=0,一个边 温度为 T2=1,求该矩形区域内的温度分布。\n"); printf("\n 是否要手动对温度场赋予初值?(Y/N):"); scanf("%c",&s); if(s=='y'||s=='Y') /*手动赋予温度初场*/ { printf("\n 请首先假定一个温度场的初始分布,即给出各节点的温度初值(一行一行进行):\n"); for(i=0;i
{ t[i][0]=0; t[i][N-1]=1; } /*步长计算*/ x=1.0/(N-1); y=1.0/(M-1); /*迭代循环*/ cha=1; while(cha>0.0001) { for(i=0;i=0;j--) for(i=0;i
分享到:
收藏