Boucle Matlab

[invite3b879c03]

Bonjour,

J'ai un une courbe de distribution de température à réaliser en transfert de chaleur. Dans mon programme Matlab, je dois trouver une valeur de température optimisé que je donne au début arbitraire. Voici le code:

L=1.0; %length of the insulated warm-compost (m)
H=1.5; %height of the insulated warm-compost (m)
t=0.02; %plastic wall thickness (m)

b_ins=0.05;

T_compost=30+273;
T1=10+273; %inside wall temperature (°C)
T_outside=-20+273;

k_PET=0.2; %thermal conductivity of the polyethylene terephtalate (PET)(W/m.K)
k_PVC=0.14; %thermal conductivity of the polyvinylechloride (PVC) (W/m.K)

A=4*L*H+L*L; %external area of the compost (m^2)
V=L*L*H; %volume of the compost (m^3)

k_ins=k_PVC;

boltzman=5.67e-8;
epsilon=0.95;

h_conv=20;
T4app=-10+273;

while T4>T4app
h_rad=boltzman*epsilon*((T4app ^4)-(T_outside^4))/(T4app-T_outside);
U_4_ext=1/((1/h_conv)+(1/h_rad));
Q_4_ext=U_4_ext*A*(T4app-T_outside);

Q_1_4=Q_4_ext;
U_1_4=1/((2*t/k_PET)+(b_ins/k_ins));
T4=T1-(Q_1_4/(A*U_1_4));
T4app=T4;
end

Donc voila je donne T4app a -10°C par exemple. Donc cette valeur je veux qu'elle soit prise dans mes formules dnas la boucle et obtenir une nouvelle valeur T4 (plus précise en théorie). Je veux donc ensuite comparer cette valeur à celle que j'ai approximé. Si T4-T4app<1 par exemple, alors je veux obtenir cette dernière valeur de T4. Sinon T4 = mon nouveau T4app et on recommence une boucle. etc

Voilà j'espère que vous pourrez m'aider, merci.
-----

[invite3b879c03]

Le problème que j'ai, c'est que quand je veux faire tourner, il me dit que T4 n'est pas défini (normal T4 est dans l boucle), j'ai donc essayer de duppliquer le code pour faire un tour de boucle "manuellement" mais apres il rentre même plus dans la boucle. Je vois pas mon erreur...

[invited1c1a33e]

Bonjour, je pense qu'une boucle do...while serait plus appropriée (si ça existe en matlab), au pire initialisez T4 avec n'importe quelle valeur > T4app

par exemple:

T4app=-10+273;
T4=T4app+1;
while T4>T4app

mais le do ... while est plus rigoureux.

Si ça ne rentre pas dans la boucle c'est que T4 <= T4app dès la première itération (les afficher?)


A mon avis ce code est à corriger:
T4=T1-(Q_1_4/(A*U_1_4));
T4app=T4;

[invite3b879c03]

Do while n'existe pas dans matlab.

Pourquoi devrais-je modifier ces 2 lignes?

La première ligne correspond a la nouvelle temperature que j'obtiens avec celle approximée.

La seconde c'est juste pour remplacer mon T4app par on T4 calculé pour refaire une boucle jusqu'a remplir la condition.

[A voir en vidéo sur Futura]

[invited1c1a33e]

Pourquoi devrais-je modifier ces 2 lignes?

Parce que si T4app et T4 sont identiques comment à la prochaine itération on peut avoir la condition T4>T4app qui soit vraie ?

[invite3b879c03]

Ok c'est bon je pense avoir trouvé avec ce code-ci:

%-----------------1st loop to get a initial value of T4------------------%
T4app=-19+273;
h_rad=boltzman*epsilon*((T4app ^4)-(T_outside^4))/(T4app-T_outside);
U_4_ext=1/((1/h_conv)+(1/h_rad));
Q_4_ext=U_4_ext*A*(T4app-T_outside);

Q_1_4=Q_4_ext;
U_1_4=1/((2*t/k_PET)+(b_ins/k_ins));
T4=T1-(Q_1_4/(A*U_1_4));
T4-273

%-----------------While loop to get the real value of T4------------------%
while abs(T4-T4app)>0.5
h_rad=boltzman*epsilon*((T4app ^4)-(T_outside^4))/(T4app-T_outside);
U_4_ext=1/((1/h_conv)+(1/h_rad));
Q_4_ext=U_4_ext*A*(T4app-T_outside);

Q_1_4=Q_4_ext;
U_1_4=1/((2*t/k_PET)+(b_ins/k_ins));
T4=T1-(Q_1_4/(A*U_1_4));
if T4>T4app
T4app=T4app+0.1;
else
T4app=T4app-0.1;
end
end