Débit d'une pompe

Débit normé.

Il faut d'abord créer une fonction DebitNorm.

function Qnorm=DebitNorm(teta,n) %calcul du debit norme d'une pompe % if (n==1) Qnorm=0.5*(sin(teta)+abs(sin(teta)) ); else tempo=0; for k=0:n-1 tempo=tempo+0.5*(sin(teta+(2*pi*k/n))+abs(sin(teta+(2*pi*k/n))) ); end; Qnorm=tempo; end;

puis le script du tracé du débit pour divers nombre de pistons.

%----------------------------------------------------------- % calcul du débit normé d'une pompe % à N pistons axiaux % (c) JPaul Molina %------------------------------------------------------------ clear all; echo on %-------------------------------------------------------------- % calcul du débit normé d'une pompe % à N pistons axiaux %--------------------------------------------------------------- echo off teta=0:0.1:2*pi; %teta varie de 0 à 2pi %etude comparative de 1 à 7 pistons Q1=DebitNorm(teta,1); Q2=DebitNorm(teta,2); Q3=DebitNorm(teta,3); Q4=DebitNorm(teta,4); Q5=DebitNorm(teta,5); Q6=DebitNorm(teta,6); Q7=DebitNorm(teta,7); clf; %tracé plot(teta,Q1,teta,Q2,teta,Q3,teta,Q4,teta,Q5,teta,Q6,teta,Q7); % mise en place des légendes title('Debit normé en fonction du nombre de pistons'); xlabel('teta en rd'); % légende abscisse ylabel('Q/Q°'); % légende ordonnée % écrit les textes relatifs à chaque courbe text( teta(5), Q1(5), '1' ); text( teta(5), Q2(6), '2' ); % décalage car une partie de Q1 est identique text( teta(5), Q3(5), '3' ); text( teta(5), Q4(5), '4' ); text( teta(5), Q5(5), '5' ); text( teta(5), Q6(5), '6' ); text( teta(5), Q7(5), '7' ); grid;

Dernière Modification : Lun 19 Mars 2007 13:28 Copyright © 1999-2010 Jean-Paul Molina Tous droits réservés.