Discussion :

[anass-07]

Bonjour à tous,

Je débute avec matlab et j'aurais voulu simuler le comportement d'un compresseur.

J'ai réussi à simuler le comportement de mon compresseur mais je veux introduire des tests qui peuvent me permettre de couper le compresseur si je dépasse une certaine pression seuil et de le redémarrer si la pression tombe en dessous d'une certaine valeur. Est-ce que quelqu'un peut m'aider à faire ce genre de truc?
Je vous remercie d'avance.

Anass

[thecrazydonut]

Bonjour,

Jette un oeil sur comment bien poser les questions.

On ne sait pas ce que tu as fait réellement, c'est pas vraiment une question à propos de Matlab mais plus à propos de ton algorithme complet, non ?
Pour mettre une condition, tu utilises "if" mais je présume que tu connais ça si tu as réussi à faire tout le travail en amont.

[anass-07]

Bonjour thecrazydonut,
Tout d abord merci pour la réponse rapide.
Mon travail consiste à simuler à l aide d un code Matlab le comportement d un compresseur dont je connais la fonction caractéristique p= f( q) où p est la pression et q est le débit . J ai donc utiliser un schéma d euler afin de résoudre l équation de conservation de la masse. Pour ce faire j ai créer une fonction qui fait ce calcul et qui me renvois la valeur de la pression a chaque pas de temps. Le but de mon post est de savoir comment introduire des conditions pour que cette fonction fasse le calcul jusqu à ce que la pression atteind une valeur Max. Jusque la j y arrive avec une boucle if ( j ai mis if p< pmax faire le calcul else faire le calcul avec un débit nul) mais le problème c est qu après avoir atteind la pmax le code utilise la seconde condition qui lui dit de prendre un début nul. Jusque la ca va . Mais moi je voudrais que la pression ne descende pas en dessous d une pmin autrement dit je voudrais redémarrer le compresseur lorsque la pression descende en dessous de pmin et ça j y arrive pas.
J espère que j ai été clair s il y a encore des lacunes sollicitez moi je rester à votre disposition .
Merci
Anass

[Jerome Briot]

Montre nous la code que tu as déjà écris

[anass-07]

Bonjour Dut,

voici le code comme vous l'avez demandé:

code pricipal :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
clear all
close all
clc
 
 
%-------------------------------------------------------------------------%
 
%                       Programme de simulations
 
%-------------------------------------------------------------------------%
 
 
% Entrée des données du compresseur 
 
% air
 
r_air = 287 ;   % J/K/kg
 
 
% à l'admission
 
T_1   = input('Température à l''admission  [°C]:                 ') + 273.15;        % K
p_1   = input('Pression à l''admission   [bars]:                 ') * 10^5 ;         % Pa
RHO_1 =  p_1 / ( r_air * T_1 ) ;
disp(['Masse volumique de l''air à l''admission [kg/m³] : ',num2str(RHO_1)]);
 
 
% au refoulement
 
T_2   = input('Température au refoulement [°C]:                 ') + 273.15;        % K
 
 
% Entrée des données des réservoirs
 
V_R = input( 'Volume du buffer           [m³]:                 ') ; % m³
 
% Entrée des données d'utilisation
 
q_u = input('Débit d''utilisation       [l/s]:                 ') ;        % l/s
 
 
%-------------------------------------------------------------------------%
% Résolution de l'équation différentielle avec et sanas fuites 
%-------------------------------------------------------------------------%
 
% Débit de fuites
 
% sans fuites
 
q_f = 0 ;
 
h    = 0.01;                           % pas de temps = 0.01 secondes
T    = 1:h:600;                        % calculer jusqu'à 10 minutes
p    = zeros(length(V_R),length(T)); 
 
 
for n=1:(length(T)-1)                 % bouclage
 
    p(1) = p_1 / 10^5;                % condition initiale p1 = 1 bar 
 
    k_1   = fc(h, T(n), p(n), q_u, q_f, RHO_1, V_R, T_2)  ;
 
    p(n+1) =  k_1 ;    % Euler
 
end

fonction fc:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
function F = fc(Delta_t,T, p , q_u, q_f, RHO_1, V_R, T_2)
 
% -------------------------------------------------------------------------
%                               Données
% -------------------------------------------------------------------------
 
r_air = 287 ; 
 
% -------------------------------------------------------------------------
%                               Résolution
% ------------------------------------------------------------------------- 
if p<8
        m_1     = p * (10^5) * V_R / (r_air *T_2 ) ;
 
        Delta   = ( 1703 / 30 )^2 - 4 * ( 24160 + p ) / 30 ;
 
        q_v_1   = ( - ( 1703 / 30 ) - Delta ^ 0.5 ) * ( - 15 ) ;  % débit volumique à l'entrée
 
        q_m_1   = RHO_1 * q_v_1 / 1000    ;
 
        q_m_out = RHO_1 * q_u   / 1000    ;
 
        q_m_f   = RHO_1 * q_f   / 1000    ;
 
        Delta_m = Delta_t*(q_m_1 - q_m_out - q_m_f) ;
 
        m_f     = Delta_m + m_1  ;
 
        F       = m_f * r_air * T_2 / (V_R * 10^5) ;
else
        m_1     = p * (10^5) * V_R / (r_air *T_2 ) ;
 
        q_m_1   = 0    ;
 
        q_m_out = RHO_1 * q_u   / 1000    ;
 
        q_m_f   = RHO_1 * q_f   / 1000    ;
 
        Delta_m = Delta_t*(q_m_1 - q_m_out - q_m_f) ;
 
        m_f     = Delta_m + m_1  ;
 
        F       = m_f * r_air * T_2 / (V_R * 10^5) ;
end