équation de la chaleur en VBA

[invite6d271278]

Bonjour,

Cela fait maintenant plus d'un mois, plusieurs soir par semaines que je bloque sur un exercice qu'on m'a donné a faire.
L'exercice est en 3 parties, la 1ere partie j'ai réussis a la faire (après plusieurs heures de travail). Je ne fait pas la 3eme parti (mes collègues devrait s'en charger). Par contre, je blogue littéralement sur la 2éme partie !!

il faut résoudre l'équation de la chaleur en VBA.
Autrement dit, il faut être calé en info, et en physique.
J'ai fait appel a tout mon réseau pour m'aider sur cette exercice, personne n'y est arrivé...

J'essaye donc ma chance ici... ma dernière chance :/

J'ai tout de même réaliser un programme en VBA qui affiche des résultats,mais ils sont incohérents avec les phénomènes physique.
Je vous mets en pièce jointe l’énoncer du problème ainsi que mon programme VBA.

Si quelqu'un trouve mon erreur, je suis preneur !!
Merci d'avance

A bientôt.

Cordialement
-----

[JPL]

Poste le programme dans le corps d'un message en utilisant la balise Code (# dans la barre d'outils de l'éditeur en mode avancé). Merci.

[invite6d271278]

Voici le programme en question :

 sub Fourier()
Dim nc(3) As Double 'Coupure

'paramètres physiques
L = Cells(2, 2)

rho = Cells(3, 2)
cp = Cells(4, 2)
lambda = Cells(5, 2)
alpha = lambda / (rho * cp)
alpha = 10 / (rho * cp)
N = Cells(12, 2)
Pi = Cells(15, 2)

temps = 360

T0 = Cells(6, 2) 'Tg
Tl = Cells(7, 2) 'Td
Ti = Cells(9, 2) 'T initiale

'déclaration des tableaux et des matrices
ReDim x(N) As Double  'points
ReDim txt(N) As Double 'Temperature


nc(1) = Cells(13, 2)
nc(2) = Cells(13, 3)
nc(3) = Cells(13, 4)

'paramètres de calcul intérmédiaires
dx = L / (N - 1)

' ainsi que les paramètres de calcul
For j = 1 To N
    x(j) = (j - 1) * dx
Next j
 
For k = 1 To 3
    Cells(1, 10 + k) = nc(k)
    Cells(2, 10 + k) = T0
    For j = 2 To N 'Boucle en chaque point
        txt(j) = Ti
        For i = 1 To nc(k) 'Somme sur les i
            Ipi = i * Pi
            A = (-2 / Ipi) * (Ti * (Cos(Ipi) - 1) - Tl * Cos(Ipi) + T0)
        
            txt(j) = txt(j) + A * Exp(-alpha * ((Ipi / L) ^ 2) * temps)
        Next i
        txt(j) = txt(j) + T0 + (Tl - T0) / L * x(j)
        'Afichage des résultats
        Cells(j + 1, 10 + k) = txt(j)
    Next j
Next k

'Afichage des x
For j = 1 To N
    Cells(j + 1, 10) = x(j)
Next j
Cells(1, 10) = "x"

End Sub