低压燃气管网水力计算程序

用该程序 可实现燃气管网的水力计算 计算前现将管网信息输入相关的文件

#include 
#include
#include
#define NUMBER  50 

void main（）
{
	double LAMD（int xint yint zfloat hdouble q）;
	
	FILE *fp1*fp2*fp3;	
    int C[NUMBER][5]A[NUMBER][NUMBER]M[NUMBER][NUMBER];
	int flagdc=0;
	int zrijlnmwTkp;//m 为节点数 n为管段数  k为已知压力节点号 p为接点压力
    float vmdQ[NUMBER][NUMBER]jz;
	 //  黏度  密度  节点流量   局部阻力
 	double ldsumq[NUMBER][NUMBER]G[NUMBER][NUMBER]H[NUMBER][NUMBER]N[NUMBER][NUMBER]F[NUMBER][NUMBER];
    //摩阻系数    管段流量           阻抗矩阵           A*G              导纳矩阵                系数矩阵  
	double x[NUMBER] tempt[NUMBER][NUMBER]W[NUMBER]DP[NUMBER]wc[NUMBER][NUMBER];
   	//    压力矩阵         管段流量暂存数组   速度数组   管段压降
	fp1=fopen（“a1.dat““r“）;
	
	fscanf（fp1“%d %d %d %d %d %f %f %d %f“&m&n&w&k&p&v&md&T&jz）;
     for（i=0;i	fscanf（fp1“%f “&Q[i][0]）;

	 /*for（i=0;i		for（j=0;j			wc[i][j]=0;*/  // 误差数组
 
	for（i=0;i	{		
		for（j=0;j<5;j++）
		fscanf（fp1“%d“&C[i][j]）;// 起点  终点  管长  管径  管材
	
	}
	fclose（fp1）;
	
	for（i=0;i		for（j=0;j			A[i][j]=0;
		for（i=0;i		{ 
			A[C[i][0]-1][i]=1;
			
			A[C[i][1]-1][i]=-1;
		}

/*for（i=0;i		{	for（j=0;j	      printf（“%d“A[i][j]）;
		  printf（“\n“）;
		}*/
			
		
		
		for（i=0;i			for（j=0;j				G[i][j]=0;	
		
			for （i=0;i	
		q[i][0]=（3.1416*C[i][3]*C[i][3]*w）*3600/（4*pow（106））;	
				
		   	/*for（i=0;i			{ 
			    printf（“%lf“q[i][0]）;
				printf（“\n“）; 
			}*/
			
				for （i=0;i			 t[i][0]=0;
			
				
again:		for（i=0;i   
	while（fabs（（q[i][0]-t[i][0]）/q[i][0]）>0.01）	
			
	{
	      for （i=0;i		   t[i][0]=q[i][0];

			for （i=0;i			{	
		        ld=LAMD（ C[i][3] C[i][4]wvt[i][0]）;
			   G[i][i]=1/（（6.26*（1e+7））*ld*（1+jz）*t[i][0]*md*T*C[i][2]/（pow（C[i][3]5）*273））;	

	           }	
		for （i=0;i			{
			 for（l=0;l			 {
				{ sum=0;
				for（j=0;j					sum=sum+A[i][j]*G[j][l];
				}
				H[i][l]=sum;
			 }
			
			}
	 			
	  	/*for （i=0;i		{
			sum=0;
			for（j=0;j				sum=sum+A[i][j]*t[j][0];
			
			Q[i][0]=sum;
		}*/
	

  /*	for（j=0;j	{
		printf（“%lf“-
			Q[j][0]）;
		  printf（“\n“）;
		
	}*/
		for（i=0;i			for（j=0;j				M[i][j]=A[j][i];
			
	  /*	for（i=0;i			{
				for（j=0;j			
	        	printf（“%5d“M[i][j]）;
	       	    printf（“\n“）;
			}  */

			
			for （i=0;i			{
				for（l=0;l				{
					{ sum=0;
					for（j=0;j						sum=sum+H[i][j]*M[j][l];
					}
					N[i][l]=sum;
				}
				
			}
		 
    	 	for（i=0;i			{
				for（j=0;j			
	        	printf（“%10.3lf“N[i][j]）;
	       	    printf（“\n“）;
			}


	//改写导纳矩阵//	 
      for（i=0;i				N[i][k-1]=-Q[i][0]-N[i][k-1]*p;
	 
			/*for（i=0;i			{
				for（j=0;j			
	        	printf（“%lf“N[i][

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       7401  2009-06-06 19:31  低压燃气管网计算\Cpp1.cpp

     文件        134  2009-04-20 14:46  低压燃气管网计算\a1.dat

     目录          0  2009-06-16 17:06  低压燃气管网计算

----------- ---------  ---------- -----  ----

                 7535                    3