二维有限体积法计算热传导及源码//#include"stdafx.h"#include#include#include#include#includeusingnamespacestd;#defineq500#definek1000voidmain(){//inputthevalueyouwant:doubleL,dx,dy,T,Ax,Ay,d;intm,n,i,j,kk,mm;//charstr1[20];char*str1="conf.dat";//cout>str1;//stringstr1("data.");ifstreamfin(st...
//#include"stdafx.h"#include#include#include#include#includeusingnamespacestd;#defineq500#definek1000voidmain(){//inputthevalueyouwant:doubleL,dx,dy,T,Ax,Ay,d;intm,n,i,j,kk,mm;//charstr1[20];char*str1="conf.dat";//cout<<"请?输º?入¨?文?件t名?：êo"<>str1;//stringstr1("data.");ifstreamfin(str1);if(!fin){cout<>L>>m>>n>>T;cout<ae(mm),ap(mm),aw(mm),an(mm),as(mm),su(mm),sp(mm);Ax=L*dy;Ay=L*dx;//结¨¢点Ì?0,左Á¨®下?角?顶£¤点Ì?aw[0]=0.0;ae[0]=k*Ax/dx;an[0]=k*Ay/dy;as[0]=0.0;sp[0]=0.0;su[0]=q*Ax*k;ap[0]=aw[0]+ae[0]+an[0]+as[0]-sp[0];for(i=1;imm-n)//最Á?后¨®一°?列¢D{ae[i]=0;aw[i]=k*Ax/dx;an[i]=k*Ay/dy;as[i]=an[i];su[i]=0.0;sp[i]=0.0;ap[i]=aw[i]+ae[i]+as[i]+an[i]-sp[i];}else//内¨²点Ì?{ae[i]=k*Ax/dx;aw[i]=k*Ax/dx;an[i]=k*Ay/dy;as[i]=an[i];su[i]=0.0;sp[i]=0.0;ap[i]=aw[i]+ae[i]+as[i]+an[i]-sp[i];}}//......................///double*b=newdouble[mm];double**maxtrix;maxtrix=newdouble*[mm];for(i=0;i=mm-n){if(i==mm-1){maxtrix[i][i-1]=-as[i];maxtrix[i][i-n]=-aw[i];}elseif(i==mm-n){maxtrix[i][i+1]=-an[i];maxtrix[i][i-n]=-aw[i];}else{maxtrix[i][i+1]=-an[i];maxtrix[i][i-n]=-aw[i];maxtrix[i][i-1]=-as[i];}}elseif(i%n==0){maxtrix[i][i+1]=-an[i];maxtrix[i][i+n]=-ae[i];maxtrix[i][i-n]=-aw[i];}elseif((i+1)%n==0){maxtrix[i][i-1]=-as[i];maxtrix[i][i+n]=-ae[i];maxtrix[i][i-n]=-aw[i];}else{maxtrix[i][i+1]=-an[i];maxtrix[i][i+n]=-ae[i];maxtrix[i][i-n]=-aw[i];maxtrix[i][i-1]=-as[i];}}}cout<=0;kk--){d=0.0;for(j=kk+1;j
本文档为【二维有限体积法计算热传导及源码】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。
该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。
[版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
浙公网安备 33021202002483