Discussion :

[Duch9]

Bonjour à toutes et à tous,

Pour mes études, je dois faire l'analyse thermique sur un modèle 1D d'une paroi constituée de plusieurs matériaux (voir .docx joint). Pour cela, j'ai appliqué la méthode nodale.
Cela étant fait, je dois faire tourner ce modèle sous MATLAB. L'utilisant depuis peu, je rencontre quelques petits problèmes. J'ai joint 2 scripts :

1-Prog_Paroi_V3 : c'est mon 1er programme. le résultat n'est pas réaliste. J'obtient des complexes dans mes matrices !! Après analyse, la seule réponse que j'ai trouvé et que mon maillage n'était pas régulier.

2-Paroi_Maillage : 2e programme. J'ai tenté de faire un maillage régulier, mais j'obtient des erreurs dans les dimensions de mes matrices (alors qu'elles sont corrects !)

Tournant en rond sur ce problème depuis un bout de temps, je sollicite donc votre aide.
Est-ce que quelqu'un a une idée pour résoudre mon problème ?

D'avance merci.

[Invité]

Bonsoir,

Je n'ai pas regardé en détails, les formules et tout, manque de motivation

1. Prog_Paroi_V3 : outre la variable D_IM non définie, les valeurs complexes proviennent initialement de la ligne :
   

   Code :

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

   GV01=1.78*(Text-T(1,k-1))^0.25;

   On a Text = 10 et T(1,k-1) = 20, (-10)^0.25 complexe.
2. Paroi_Maillage : L'erreur étant
   

   ??? Error using ==> minus
   Matrix dimensions must agree.
   
   Error in ==> Paroi_Maillage at 200
   Z=(C/dt)-(aa*A);     %résolution de CT=AT+B

   Elle vient du fait que A est de taille 1642 x 1641 alors que C est de taille 1641 x 1641.
   A qui prend cette taille suite à :
   

   Code :

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

   1
   2
   3
   4
   5
   6
    %Dans Poly
       for i=1623:1641
           A(i-1,i)=GC78;
           A(i,i)=-2*GC78;
           A(i+1,i)=GC78;
       end

[Duch9]

Bonjour Winjerome,

Merci pour ta réponse rapide.
je vais faire les modifs et les tests (demain). Je te tiens au courant.

@+

[Duch9]

Re-bonjour,

J'ai apporté les modifs sur mes 2 prog.

Dans le 1er, j'ai pris les valeurs absolue , ce qui supprime le problème de complexe. Par contre, les résultats obtenus sont les mêmes qu'avant. Pourquoi mon point 6 (dans PU) ne varie-t-il pas ?

Dans le 2e, la modif m'a permis d'obtenir des résultats. Mais ces derniers m'ont l'air encore moins cohérents !

Du coup, je suis un peu perdu : est-ce un problème de programmation ou est-ce que mon système thermique est bien posé (pourtant j'ai fait vérifié ma partie analytique) ?

D'avance merci pour vos réponses. En espérant arriver à bout de ce problème.