Répondre à la discussion
Affichage des résultats 1 à 5 sur 5

Résolution d'un système d'équations différentielles d'ordre 1



  1. #1
    layougue

    Résolution d'un système d'équations différentielles d'ordre 1


    ------

    Bonjour,
    je suis confronté à un problème de thermique dans lequel
    je dois modéliser l'évolution en température de corps en interactions. le bilan des équations en stationnaire me donne le système suivant:
    P=[rr(1,4)*((Tk(4)^4)-(Tk(1)^4))-rr(1,7)*((Tk(1)^4)-(Tk(7)^4))-rcd(1)*(Tk(1)-Tk(2));
    ...;
    rr(7,1)*((Tk(1)^4)-(Tk(7)^4))-rr(7,3)*((Tk(7)^4)-(Tk(3)^4))+rr(7,4)*((Tk(4)^4)-(Tk(7)^4))+rr(7,6)*((Tk(6)^4)-(Tk(7)^4))-rcd(6)*(Tk(7)-Text)];
    7 inconnues Tk et des constantes comme rr(i,j) Text,etc.
    je résoud ce système par Newton-raphson et tout se passe bien.

    mais je suis completement perdu, quand il s'agit de passer en transitoire, càd à ce système se rajoute les termes inertiels de chaleur stockée par les corps : mC(dTk/dt)
    tout en conservant, ma méthode de résolution par newton-raphson.

    quelle procédure dois-je appliquer? aidez-moi svp!

    merci d'avance.

    -----

  2. Publicité
  3. #2
    canard

    Talking Re : Résolution d'un système d'équations différentielles d'ordre 1

    t'as pas plus long comme éuqation lol

  4. #3
    layougue

    Re : Résolution d'un système d'équations différentielles d'ordre 1

    c'est en fait un système d'équations (il y en a sept pour etre précis), mis sous forme matricielle qui peut etre résolu en permanent par newton raphson par:
    [Jacobienne]*[DT]k-k+1=[P]

    sinon, si tu veux j'ai des équation un peu plus longue, mais pour mon problème çà ne m'aide pas trop

    au plaisir

  5. #4
    edpiste

    Re : Résolution d'un système d'équations différentielles d'ordre 1

    Je crois que ce que veut dire Canard, c'est qu'il faut que tu fasses l'effort de présenter la question sous une forme synthétique, dans laquelle on voit apparaitre les données et les inconnues.
    Ce qui n'est toujours pas le cas de ton dernier message. Que sont ces termes "inertiels" dont tu parles ? Où est le vecteur inconnu dans ton système écrit sous forme matricielle ?

  6. A voir en vidéo sur Futura
  7. #5
    layougue

    Re : Résolution d'un système d'équations différentielles d'ordre 1

    excusez moi, je n'avais pas compris. désolé.

    je vais essayer d'être plus claire.
    mon problème est de résoudre un système d'équations polynomiale d'ordre 4 et multi variables (7 températures inconnues de 7 corps), correspondant aux echanges thermiques radiatifs en conductifs dans une enceinte.
    l'équation utilisé pour chacun des corps est l'équation de la chaleur:
    mC(dT/dt)=somme des flux entrants-somme des flux sortant + flux généré=R (T)
    ou mC(dT/dt) est mon terme inertiel de montée en température.
    en permanent dT/dt=0, donc je résoud itérativement mon système d'équations par newton raphson

    J*DT=-R
    R: équations bilans des corps (vecteur)
    J: jacob. de R
    DT:l'accroisement entre deux itérations de T (vecteur)
    il me revient à résoudre
    Tk+1=Tk+DT, avec DT=-J^(-1)*R
    çà, çà va (c'est bon tout le monde suit)

    maintenant, quand dT/dt~=0
    je voudrais faire le lien avec le calcul précédent et l'équation de la chaleur qui developpée donne:

    Tk+1(t+dt)=Tk(t)+dt*(R)/(mC)

    dt pas de temps

    j'espère que c'est plus claire...

    pas sur!

Discussions similaires

  1. Résolution numérique d'un système d'équations différentielles d'ordre 2 non linéaire
    Par mariejuana dans le forum Mathématiques du supérieur
    Réponses: 4
    Dernier message: 02/07/2008, 16h55
  2. Résolution algébrique d'un système d'équations sous Matlab
    Par jbollagnier dans le forum Mathématiques du supérieur
    Réponses: 6
    Dernier message: 22/10/2007, 10h45
  3. la résolution d'un systéme d'équations différentielles
    Par H.Dalila dans le forum Mathématiques du supérieur
    Réponses: 0
    Dernier message: 26/04/2007, 14h07
  4. Résolution d'équations différentielles
    Par redvivi dans le forum Mathématiques du supérieur
    Réponses: 6
    Dernier message: 14/01/2007, 12h59
  5. résolution d'équations différentielles partielles
    Par loki dans le forum Mathématiques du supérieur
    Réponses: 6
    Dernier message: 31/05/2006, 16h45