jsq.c

水电站流量水头计算器

#include "c8051f000.h"               
#include <intrins.h>     
#include <math.h>

#define unsigned int uint;


main()
{
	float num[9];
/*	float LL,ST,XL,GL;
	float I1,I2,I3,k;
	float LL_out,ST_out,XL_out;*/
	float D_LL,D_ST;
	float LL,ST,LL_out,ST_out;
	float I1,I2,k;
	k=9.3745;
	num[0]=10.00;
    num[1]=16.70;
    num[2]=150.00;
	num[3]=48.00;
	num[4]=0.0;
	num[5]=0.0;
    num[6]=9.80;
	num[7]=96;
    num[8]=126.00;
	/*num[0]=70.00;                                          //蜗壳压差变送器输出20mA对应的压差值
	num[1]=320.00;                                         //流量最大值，输出20mA对应的流量值，单位m3/s
	num[2]=70.00;                                          //水头压差变送器的量程，20mA对应的压差值，单位m
	num[3]=17.50;                                          //功率最大值 ，万千瓦
	num[4]=44.16;                                          //蜗壳进口断面积，米2
	num[5]=112.0;                                          //尾水管出口断面积，米2
	num[6]=9.8;                                            //当地重力加速度
	num[7]=98;                                             // 发电机效率
	num[8]=100;*/
/*	I1=6.64;
	I2=15.0;
	I3=8.75;
	k=num[3]/num[6]/num[7]/4095.0;
	LL=(I1*5.0/16.0-1.25)*2.4*4095.0/5.0/2.43;
	LL=pow(LL/4095,0.4998)*num[1];
	ST=(I2*5.0/16.0-1.25)*2.4*4095.0/5.0/2.43;
	ST=num[2]*ST/4095.0;
	XL=(I3*5.0/16.0-1.25)*2.4*4095.0/5.0/2.43;
	XL=XL*k*1.031381/LL/ST*100000000;
	GL=3.0*I3-12.0;
	LL_out=LL*4095.0*2.5/num[1]/4095.0;
	LL_out=(LL_out/2.0+10.0/16.0)*10.0/50.0;
	ST_out=ST*4095.0*2.5/num[2]/4095.0;
	ST_out=(ST_out/2.0+10.0/16.0)*10.0/50.0;
	XL_out=XL*4095.0*2.5/100.0/4095.0;
	XL_out=(LL_out/2.0+10.0/16.0)*10.0/50.0;*/
	D_LL=4095;
	D_ST=4095;
	LL=((D_LL*2.43*5.0/4095/2.43+1.25)*16/5-4)/(20-4);
	LL=9.3745*pow(num[0]*LL,0.4927);
	ST=((D_ST*2.43*5.0/4095/2.43+1.25)*16/5-4)/(20-4);
	ST=ST*num[2];
	LL_out=LL*4095.0/num[1];
	if((LL_out<=0)&&LL_out>=4095.0) LL_out=LL_out;
	else LL_out=4095.0;
	ST_out=ST*4095.0/num[8];
	if((ST_out<=0)&&LL_out>=4095.0) ST_out=LL_out;
	else ST_out=4095.0;
	I1=(20-4)/num[0]*pow(num[1]/k,1/0.4927)+4;
	I2=16*num[8]/num[2]+4;
	while(1);
}