Exemple of calculation routine CAMIL (Amil gates)

David Dorchies — 2009-10-22T15:33:51Z

SUBROUTINE CAMIL(U,Y,PARA)
	PARAMETER(NBCM=$NBCM,NBUM=$NBUM,NSING=$NTO,NBIEF=$IB)
	PARAMETER(NPOIN=$NPOIN,NSECT=$NSC)
	CHARACTER CFIC1*3
	DIMENSION Y(*),U(*),PARA(*)
C-----index of the current controller (INURE),
C-----current number of calls of the controller n (INUAP(n))
	COMMON/P1/IFIC1,INURE,INUAP(NBCM),IFIC2
	COMMON/P2/INU,IFIC3(2),FIC1(3)
	COMMON/REGU/IFIC4,CFIC1(NBUM),IFIC5(NBUM,2),NSUG(NBUM)
1 ,IFIC6(NBUM),NKSU(NBUM),FIC2(NBUM,2)
	COMMON/REG2/NUD(NBCM),IFIC7(NBCM),BFIC1(NBCM)
	COMMON/NGEO/IFIC8(NBIEF,2),IP1(NSECT),IFIC9(NSECT,2)
	COMMON/RGEO/IFIC10,FIC3(NSECT),YGEO(NPOIN),FIC4(NPOIN,3)
1 ,FIC5(NSECT,3)
	COMMON/RSIN/FIC6(NSING),ZS0(NSING),FIC7(NSING,3)
C-----Equations de Goussard
C-----Affectations
C-----PARA(1) Modele (D parmis LISTD) de vanne AMIL
	D=PARA(1)
	Dec=PARA(2)
	R=PARA(3)
	HA=PARA(4)
C-----On met un filtre d'ordre 1 : S/E = (1-a)Z-1 / 1-aZ-1 (FILTRE = a)
	FILTRE=PARA(5)
C-----ZD : hauteur du radier % cote de fond
	ZD=ZS0(NKSU(NUD(INURE)))-YGEO(IP1(NSUG(NUD(INURE))))
C-----Calcul pour cette methode
	alfmax=asin(0.45*D/R)-asin(0.2205*D/R)
	alfb=HA/R
	alfb=min(alfb,1.)
	alfb=max(alfb,-1.)
	alfb=acos(alfb)
	alfc=alfb+alfmax
	UMAX=HA-R*cos(alfc)
C-----Y(1) est en tirant d'eau (mode TYS obligatoire sur le fichier .REG)
	Z=Y(1)-HA-ZD
C-----Initialisation de Zold au premier appel
	IF(INUAP(INURE).EQ.1) PARA(10)=Z
	Z=FILTRE*PARA(10)+(1.0-FILTRE)*Z
C-----On stocke Zold filtré dans PARA(10) pour le filtre
	PARA(10)=Z
C-----Regular setting
	if (Z.LE.-Dec) then
 alfa=0
 c=0
	elseif (z.GE.0.) then
 alfa=alfmax
 c=1
	else
 alfa=alfmax-asin(-Z/Dec*sin(alfmax))
 c=(1 - (0.45*D/R)**2)**0.5 * sin(alfa)
 1 / (0.2295*D/R) + 1.9608 * (1 - cos(alfa))
	endif
C-----Calcul de la commande
	U(1)=c*UMAX
C-----Si plusieurs U (SIMO) on duplique la meme ouverture
C-----(interessant dans le cas de plusieurs vannes identiques en parallele
C en mode V par exemple)
	IF(INU.GT.1) THEN
 DO 10 I=2,INU
 U(I)=U(1)
 10 CONTINUE
	ENDIF
C-----Calcul du debit
C En réglage classique Zmax=0, en réglage alternatif Zmax=Dec
C Kj=1 + 2.2*zmax/D - 1.6*(zmax/D)**2
C Kq=1 + 2.82*zmax/D - 34.8*(zmax/D)**2 + 164*(zmax/D)**3
C x=Y(1)-z2
C if (x.LT.0) then
C Q=0.
C elseif (x.LE.(0.45*(1 - 0.7211*c**0.3)*Kj*D)) then
C Q=0.7096*c*Kq*D**2.103*x**0.397
C else
C Q=0.3113*c*Kq*Kj**0.397*D**2.5
C endif
END

Download

The AMIL method

David Dorchies — 2009-10-22T15:33:50Z

These gates can be modelled directly in Hydraulic mode from SIC Version 4.0. It is better to use this new option. Nevertheless the previous option using a regulation module is still available.

For the AMIL gates, the variables D, d, HA and FILTRE must be given under the format (3(F8.2,1X),F8.2):

* D d HA FILTRE
PS= 220. .10 1.50 0.8

- D is the gate reference from the constructor list : 80, 90, 100, 110, 125, 140, 160, 180, 200, 220, 250, 280, 315, 355, 400, 450, 500, 560, 630, 710, 800,

- d is the decrement (equal to L/50 for AMIL, where L the width of the gate, Cf. Documentation of the constructor), given in meters (m),
- HA is the level of the gate rotation axis measured from the sill (in m).

- FILTER is the coefficient of the first order filter used to damp the fluctuations of the controlled water level. It should be in the range [0 .. 1]. The filter is stronger for higher coefficients. This filter is usually compulsory to get a stable gate.

Remark: The radius of the gate R is calculated from D and the constructor table (63, 63, 63, 63, 90, 90, 90, 125, 125, 125, 160, 160, 200, 200, 250, 250, 315, 315, 400, 400, 450),

Remark: D and R are given in the above tables in centimeters. In the software, these variables are converted into (m), when reading the .REG file. These converted variables (in meters) are the ones written on the .LST file. If the decrement d is not provided (i.e. equal to zero), we take d = 2% of D (d=0.05*D).

An example of regulation file (.REG) for the AMIL gate is given bellow. The calculation routine CAMIL for the AMIL gate is given in the following section.

* Controlleur 1

D=AMIL
* DUMin DUMax UMin UMax DT ('PU',1X,4(F5.2,1X),F6.0)
PU= 0.00 0.50 0.00 2.00 0
U=AWS(X, 2415.0, 3)
Y=TYS(R, +0, +0)
* D d HA FILTER (3(F8.2,1X),F8.2)
PS= 220 .10 1.50 0.8

SIC^2: Simulation and Integration of Control for Canals

Exemple of calculation routine CAMIL (Amil gates)

The AMIL method