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 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81
| function salle_test
global Tair_exterieur iter Tair_interieur
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Modélisation d'une salle comportant 2 baies
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% données
e=20;
%R_Cond =0.084;
Tair_interieur = 30 ;
Tair_exterieur= 42;
lamda=0.95;
S=25;
%S_tot=100;
R_ext=0.704;
R_int=0.704;
Rm_int=0.296;
Rm_ext=0.296;
hcve=10;
hcvi=10;
%P_surf=1
iter=100;
%Constantes thermo-physiques
Cp_air=1004.5;
rho_air=101325/(287*(273.15+Tair_interieur));
C1=Cp_air*rho_air;
C2=Cp_air*rho_air;
h1=12.2;
h2=0.9;
%Initialisation des vecteurs de températures
T1=zeros(iter,1);
T2=zeros(iter,1);
P1=zeros(iter,1);
P2=zeros(iter,1);
Tm=zeros(iter,1);
U=zeros(5,iter);
t=zeros(iter,1);
dt=15;
I=eye(2);
%Définition des vecteurs A, T, B et U
A=[-h1/C1 h1/C1 ; h2/C2 -h2/C2];
B=[1/C1 0 ((hcvi/C1)*R_int -(hcve/C1)* R_ext + (lamda *S/e*C1)*R_ext -(lamda *S/e*C1)*R_int) ((hcvi/C1)*Rm_int-hcvi/C1-(lamda*S/e*C1)*Rm_int) ((hcve/C1)+(hcve/C1)*Rm_ext+(lamda*S/e*C1)*Rm_ext);0 1/C2 ((hcvi/C2)*R_int +(hcve/C2)* R_ext + (lamda *S /e*C2)*R_ext -(lamda *S /e*C2)*R_int) ((hcvi/C2)*Rm_int-hcvi/C2-(lamda*S/e*C2)*Rm_int) ((hcve/C2)+(hcve/C2)*Rm_ext+(lamda*S/e*C2)*Rm_ext)];
T=[T1; T2];
%Boucle itérative
for i=1:100
t(i)=(i-1)*dt;
T1(i)=T1(1);
T2(i)=T1(2);
P1(i)=P1(1);
P2(i)=P2(1);
Tm(i)=Tm(1);
U(:,i)=[P1 ; P2 ; Tm ; Tair_interieur ; Tair_exterieur];
T=(I-dt*A)\(T+dt*B*U(:,i));
end
%Tracé des mesures et des résultats du modèle
hold on
xlabel('Temps (min)')
ylabel('Température (°C)')
plot(t,T1,'c')
plot(t,T2,'c:')
plot(t,P1,'r')
plot(t,P2,'b')
plot(t,Tm,'g')
legend('T1','T2','Tm','P1','P2') |
Partager