Discussion :

[by.hugo82]

bonjour
dans le cadre de mon stage je dois faire une étude thermique d'un panneaux photovoltaïque hybride a récupération de chaleur.

j'ai effectué des bilans thermiques sur toutes les couches constituant le panneau.

j'ai un système de 8 équations différentielles d'ordre 1 a résoudre par la méthode de Runge kutta ordre 4 en scilab.

le but est d'avoir les évolutions des différentes températures des couches constituant le panneau en fonction du temps.

Ta est une constante fixes (température ambiantes)
A,B,C... : constantes fixes

dT1=A*(T2-T1)-B*(T1-Ta);

dT2=C*(T2-T1)-D*(T2-T3)-E;

dT3=F*(T2-T3)-G*(T3-T4);

0=G*(T3-T4)-H*(T4-T5)-I*(T4^4-T6^4);

dT5=H(T4-T5)-J*(T5-T6)-K*(T5-Ta);

0=I*(T4^4-T6^4)+J*(T5-T6)-L*(T6-T7);

dT(7)=L*(T6-T7)-M*(T7-T8);

0=M*(T7-T8)-N*(T8^4-Ta^4)-O*(T8-Ta);

je débute en programmation donc mon niveau est moyen
j'ai choisi scilab car c'est le language que je connais le plus

[by.hugo82]

j'ai essayé ceci mais il ne veux pas prendre les 3 équations qui = 0 (nœuds arithmétique)
veillez m'aidé svp

function dT=fct(t,T)

dT1=A*(T2-T1)-B*(T1-Ta);

dT2=C*(T2-T1)-D*(T2-T3)-E;

dT3=F*(T2-T3)-G*(T3-T4);

0=G*(T3-T4)-H*(T4-T5)-I*(T4^4-T6^4);

dT5=H(T4-T5)-J*(T5-T6)-K*(T5-Ta);

0=I*(T4^4-T6^4)+J*(T5-T6)-L*(T6-T7);

dT(7)=L*(T6-T7)-M*(T7-T8);

0=M*(T7-T8)-N*(T8^4-Ta^4)-O*(T8-Ta);

endfunction

t0=0;T0=[273;273;273;273;273;273;273;273];t=0:1:10;
z=ode('rk4',T0,t0,t,fct)

[FR119492]

Salut!
Ton système différentiel est d'ordre 5 et non d'ordre 8. Tu dois commencer par éliminer les variables T4, T6 et T8.
Jean-Marc Blanc