compute.h

《管状换热器计算机辅助设计系统ExhCAD绘图系统(版本:1.01a Final)》为自由软件

				  break;					   
			  }
			else {
			  p=q;
			  i++;
			  q=RecordData.AirPhysics[i].t;			  
			}	
		 
		return r;
					  }

		    double  __stdcall  FindAirPhysicsU(double &airt,ExhCADData &RecordData){
	    double p,q,r;
		int i;
		
		p=RecordData.AirPhysics[0].t;
		q=RecordData.AirPhysics[1].t;    
   
	    for(i=0;i<25;)
			if (p==airt){
				r=RecordData.AirPhysics[i].p;
				break;
			}   
			else if(airt>p  && airt<q) {
                 r=RecordData.AirPhysics[i-1].u+ 
					 (RecordData.AirPhysics[i].u-RecordData.AirPhysics[i-1].u) *
					 (q-p)/
					 (airt-p);
				  break;					   
			  }
			else {
			  p=q;
			  i++;
			  q=RecordData.AirPhysics[i].t;			  
			}	
		 
		return r;
					  }

			    double  __stdcall  FindAirPhysicsV(double &airt,ExhCADData &RecordData){
	    double p,q,r;
		int i;
		
		p=RecordData.AirPhysics[0].t;
		q=RecordData.AirPhysics[1].t;    

	    for(i=0;i<25;)
			if (p==airt){
				r=RecordData.AirPhysics[i].p;
				break;
			}   
			else if(airt>p  && airt<q) {
                 r=RecordData.AirPhysics[i-1].v+ 
					 (RecordData.AirPhysics[i].v-RecordData.AirPhysics[i-1].v) *
					 (q-p)/
					 (airt-p);
				  break;					   
			  }
			else {
			  p=q;
			  i++;
			  q=RecordData.AirPhysics[i].t;			  
			}	
		 
		return r;
					  }

       double  __stdcall  FindAirPhysicsPr(double &airt,ExhCADData &RecordData){
	    double p,q,r;
		int i;
		
		p=RecordData.AirPhysics[0].t;
		q=RecordData.AirPhysics[1].t;    
 
	    for(i=0;i<25;)
			if (q==airt){
				r=RecordData.AirPhysics[i].p;
				break;
			}   
			else if(airt>p  && airt<q) {
                 r=RecordData.AirPhysics[i-1].pr+ 
					 (RecordData.AirPhysics[i].pr-RecordData.AirPhysics[i-1].pr) *
					 (q-p)/
					 (airt-p);
				  break;					   
			  }
			else {
			  p=q;
			  i++;
			  q=RecordData.AirPhysics[i].t;			  
			}	
		 
		return r;
					  }


	    double  __stdcall  FindFumePhysicsP(double &fumet,ExhCADData &RecordData){
	    double p,q,r;
		int i;
		
		p=RecordData.FumePhysics[0].t;
		q=RecordData.FumePhysics[1].t;    
        
	    for(i=0;i<12;)
			if (p==fumet){
				r=RecordData.FumePhysics[i].p;
				break;
			}   
			else if(fumet>p  && fumet<q) {
                 r=RecordData.FumePhysics[i-1].p+ 
					 (RecordData.FumePhysics[i].p-RecordData.FumePhysics[i-1].p) *
					 (q-p)/
					 (fumet-p);
				  break;					   
			  }
			else {
			  p=q;
			  i++;
			  q=RecordData.FumePhysics[i].t;			  
			}	
		 
		return r;
					  }

		    double  __stdcall  FindFumePhysicsCp(double &fumet,ExhCADData &RecordData){
	    double p,q,r;
		int i;
		
		p=RecordData.FumePhysics[0].t;
		q=RecordData.FumePhysics[1].t;    
        
	    for(i=0;i<12;)
			if (p==fumet){
				r=RecordData.FumePhysics[i].cp;
				break;
			}   
			else if(fumet>p  && fumet<q) {
                 r=RecordData.FumePhysics[i-1].cp+ 
					 (RecordData.FumePhysics[i].cp-RecordData.FumePhysics[i-1].cp) *
					 (q-p)/(fumet-p);
				  break;					   
			  }
			else {
			  p=q;
			  i++;
			  q=RecordData.FumePhysics[i].t;			  
			}	
		 
		return r;
					  }

				    double  __stdcall  FindFumePhysicsR(double &fumet,ExhCADData &RecordData){
	    double p,q,r;
		int i;
		
		p=RecordData.FumePhysics[0].t;
		q=RecordData.FumePhysics[1].t;    
        
	    for(i=0;i<12;)
			if (p==fumet){
				r=RecordData.FumePhysics[i].r;
				break;
			}   
			else if(fumet>p  && fumet<q) {
                 r=RecordData.FumePhysics[i-1].r+ 
					 (RecordData.FumePhysics[i].r-RecordData.FumePhysics[i-1].r) *
					 (q-p)/(fumet-p);
				  break;					   
			  }
			else {
			  p=q;
			  i++;
			  q=RecordData.FumePhysics[i].t;			  
			}	
		 
		return r;
					  }

		
   double  __stdcall  FindFumePhysicsU(double &fumet,ExhCADData &RecordData){
	    double p,q,r;
		int i;
		
		p=RecordData.FumePhysics[0].t;
		q=RecordData.FumePhysics[1].t;    
        
	    for(i=0;i<12;)
			if (p==fumet){
				r=RecordData.FumePhysics[i].u;
				break;
			}   
			else if(fumet>p  && fumet<q) {
                 r=RecordData.FumePhysics[i-1].u+ 
					 (RecordData.FumePhysics[i].u-RecordData.FumePhysics[i-1].u) *
					 (q-p)/(fumet-p);
				  break;					   
			  }
			else {
			  p=q;
			  i++;
			  q=RecordData.FumePhysics[i].t;			  
			}	
		 
		return r;
					  }

  double  __stdcall  FindFumePhysicsV(double &fumet,ExhCADData &RecordData){
	    double p,q,r;
		int i;
		
		p=RecordData.FumePhysics[0].t;
		q=RecordData.FumePhysics[1].t;    
        
	    for(i=0;i<12;)
			if (p==fumet){
				r=RecordData.FumePhysics[i].v;
				break;
			}   
			else if(fumet>p  && fumet<q) {
                 r=RecordData.FumePhysics[i-1].v+ 
					 (RecordData.FumePhysics[i].v-RecordData.FumePhysics[i-1].v) *
					 (q-p)/(fumet-p);
				  break;					   
			  }
			else {
			  p=q;
			  i++;
			  q=RecordData.FumePhysics[i].t;			  
			}	
		 
		return r;
					  }

  	 double  __stdcall  FindFumePhysicsPr(double &fumet,ExhCADData &RecordData){
	    double p,q,r;
		int i;
		
		p=RecordData.FumePhysics[0].t;
		q=RecordData.FumePhysics[1].t;    
        
	    for(i=0;i<12;)
			if (p==fumet){
				r=RecordData.FumePhysics[i].pr;
				break;
			}   
			else if(fumet>p  && fumet<q) {
                 r=RecordData.FumePhysics[i-1].pr+ 
					 (RecordData.FumePhysics[i].pr-RecordData.FumePhysics[i-1].pr) *
					 (q-p)/(fumet-p);
				  break;					   
			  }
			else {
			  p=q;
			  i++;
			  q=RecordData.FumePhysics[i].t;			  
			}	
		 
		return r;
					  }

	 double __stdcall FindLenModify(double &ld,ExhCADData &RecordData){	 

		double p,q,r=1.0;
		int i;
		
		p=RecordData.LenModify[0].ld; 
		q=RecordData.LenModify[1].ld;    
        
	    for(i=0;i<14;)
			if (p==ld){
				r=RecordData.LenModify[i].cl;
				break;
			}   
			else if(ld>p  && ld<q) {
                 r=RecordData.LenModify[i-1].cl+ 
					 (RecordData.LenModify[i].cl-RecordData.LenModify[i-1].cl) *
					 (q-p)/(ld-p);
				  break;					   
			  }
			else {
			  p=q;
			  i++;			 
			  q=RecordData.LenModify[i].ld;			  
			}	
			
			if(i>=14 && p<ld ) r=1.0; 

		return r;
	 }

	 double __stdcall FindRowModify(short &rn,short arrange,ExhCADData &RecordData){
		   int i,p;
		   double q;
		   for(i=0;i<10;i++){
			   p=RecordData.RowModify[i].rn ;
			   if (p ==rn){q=RecordData.RowModify [i].cn[arrange];break;}
		   }
		   if(i>=10 && p<rn ) q=1.0; 
		   return q;
	 }

					
	double __stdcall TransferHeatQuantity(double &inairt,double &outairt,
	                           double &airq,double &lose,ExhCADData &RecordData){                                   
						
                       double inaircp=0,outaircp=0,aircp=0,heatq=0;					
                         
						inaircp=FindAirPhysicsCp(inairt,RecordData);
						outaircp=FindAirPhysicsCp(outairt,RecordData);
                           
						aircp=(inaircp+outaircp)/2;

						heatq=airq*aircp*(outairt-inairt)*(1+lose)/2;                   
					
						 return heatq;

					  }

   double __stdcall FumeOutTemperature(double &infumet,double &fumeq,double &protect,
	                                   double &heatq,long &order,ExhCADData &RecordData){

					  double infumecp,outfumecp,outfumet,tmp;					 

					   infumecp=FindFumePhysicsCp(infumet,RecordData);
					   outfumet=infumet-350;
					   do{ 
						   tmp=outfumet;
						   outfumecp=FindFumePhysicsCp(tmp,RecordData);
						   outfumet=infumecp*infumet/outfumecp-heatq/protect/fumeq/outfumecp;
						   order+=1;
					   }while(fabs(tmp-outfumet)>Precision && order<Limit);

					   return outfumet;

							   }

    double __stdcall TemperatureDifference(double &inairt,double &outairt,double &infumet,
		                                   double &outfumet,double &fumeq,double &airq,
										   double &nroute,short &direction,ExhCADData &RecordData){
					   double ft,inaircp,infumecp,p,r,dtairmax,dtairmin,
						      dtfumemax,dtfumemin,dtop,dtpa,tmp,dt;

					   p = fabs((inairt - outairt) /(inairt - infumet));
					   inaircp=FindAirPhysicsCp(inairt,RecordData);
					   infumecp=FindFumePhysicsCp(infumet,RecordData);