La méthode AVIS - AVIO

David Dorchies — 2009-10-22T15:06:47Z

Voir la remarque sur la modélisation des vannes AMIL, AVIS et AVIO dans les versions de SIC 4.0 et ultérieures.

Pour les vannes AVIS et AVIO les 6 paramètres sont donnés au format (5(F8.2,1X),F8.2) :

* ALFMAX d R HA CHUTE FILTRE
PS= 60.00 .10 2.00 1.50 0.03 0.8

- ALFMAX est l'angle d'ouverture maximum de la vanne mesuré par rapport à la position fermée (en degrés),

- d est le décrément (de l'ordre de r/20 = 0.05*r pour AVIS et AVIO, où r est le rayon du flotteur, ou également 0.028*R pour les AVIO et les AVIS basse-chute et 0.032*R pour les AVIS haute-chute, où R est le rayon de la vanne. Le décrément peut être ajusté entre 0.05*r et 0.10*r ; Cf. Documentation du constructeur),
- R est le rayon de la vanne (en m),

- HA est la cote de l'axe de rotation de la vanne par rapport au radier (en m),
- CHUTE est la hauteur de chute éventuelle au passage à l'ouvrage entre la cote du radier et la cote du fond à l'aval de l'ouvrage, elle est positive si le fond baisse vers l'aval, elle est donnée en (m),

- FILTRE est le coefficient du filtre d'ordre 1 utilisé pour filtrer le niveau régulé. Il doit être compris entre 0 et 1 et le filtre est d'autant plus fort que le coefficient est proche de 1. Ce filtre permet de garantir la stabilité de la vanne.


Coupe vanne AVIS

Remarque : Si le décrément d n'est pas rentré (i.e. mis à zéro) on prend d=0.028*R. Attention, cette valeur n'est pas la valeur recommandée pour les AVIS haute-chute (d=0.032*R). Nous ne pouvons pas faire ce choix automatiquement car nous n'avons pas l'information si c'est une vanne AVIS ou AVIO, ni si c'est une basse-chute ou une haute-chute. Pour certaines valeurs de R il n'y a pas d'ambiguïté (ex : si R=0.90 m alors ce ne peut être qu'une AVIS haute-chute), mais pour d'autres valeurs de R les divers cas sont possibles (ex : les 4 vannes existent pour R=2.5 m).

Un exemple de fichier de régulation pour des vannes AVIS est donné ci-dessous. La routine de calcul pour la vanne AVIS (ou AVIO) est donné ci-dessous. Vous pouvez vous inspirer de cette routine pour écrire un module de régulation CUSER si vous désirez tester une autre formule.

Attention la cote régulée Y doit être déclarée en mode "tirant d'eau" (mode TYS obligatoire sur le fichier .REG)

Voici un exemple de fichier .reg correspondant à un module de régulation AVIS. Voir également les pages précédentes pour la signification des divers paramètres.

* Contrôleur 1

D=AVIS
* 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, +1)
* ALFMAX d R HA CHUTE FILTRE (5(F8.2,1X),F8.2)
PS= 60.00 .10 2.00 1.50 0.03 0.8

Exemple de routine de calcul CGEC (Vannes Avis et Avio)

David Dorchies — 2009-10-22T15:06:47Z

Voilà le module de régulation CGEC utilisé pour les vannes AVIS et AVIO. Vous pouvez vous inspirer de ce module de régulation pour en écrire d'autres. Les équations programmées proviennent d'un document de M. Goussard.

SUBROUTINE CGEC(U,Y,PARA)
	PARAMETER(NBCM=$NBCM)
	DIMENSION Y(*),U(*),PARA(*)
C-----index of the current controller (INURE),
C-----current number of calls of the controller n (INUAP(n))
	COMMON/P1/ILNPS,INURE,INUAP(NBCM),IREGS
	COMMON/P2/INU,INY,INZ,DTU,DTY,DTZ
C-----Equations de Goussard
C-----Affectations
C-----PARA(1) en degres, ALFMAX en radiants
	ALFMAX=PARA(1)*0.01745
C-----Dec est negatif dans la formule de Goussard
	Dec=-PARA(2)
	R=PARA(3)
	HA=PARA(4)
C-----CHUTE > 0 si ca descend vers l'aval
	CHUTE=PARA(5)
C-----On met un filtre d'ordre 1 : S/E = (1-a)Z-1 / 1-aZ-1 (FILTRE = a)
	FILTRE=PARA(6)
C-----Calcul pour cette methode
C-----Y(1) est en tirant d'eau (mode TYS obligatoire sur le fichier .REG)
	Z=Y(1)-CHUTE-HA
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
	IF (Z.GE.0.) THEN
 ALFA=0.
	ELSEIF (Z.LE.Dec) THEN
 ALFA=ALFMAX
	ELSE
 ALFA=0.
 IF(Dec.NE.0.) ALFA=SIN(ALFMAX)*Z/Dec
 ALFA=MIN(ALFA,1.)
 ALFA=MAX(ALFA,-1.)
 ALFA=ASIN(ALFA)
	ENDIF
C-----Calcul de U
	ALFB=0.
	IF(R.NE.0.) ALFB=HA/R
	ALFB=MIN(ALFB,1.)
	ALFB=MAX(ALFB,-1.)
	ALFB=ACOS(ALFB)
C-----Calcul de la commande
	U(1)=HA-R*COS(ALFA+ALFB)
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
END

Télécharger

SIC^2 : Logiciel de Simulation Intégrée des Canaux et de leur Contrôle

La méthode AVIS - AVIO

Exemple de routine de calcul CGEC (Vannes Avis et Avio)