I. INTRODUCTION
La recherche des solutions de problèmes de conduction par la méthode des transformées de LAPLACE s’effectue en trois étapes :
ü Ecriture de l’équation transformée du problème
ü Résolution de cette équation transformée
ü Retour à la solution physique par application du théorème d’inversion ou par utilisation de tables de transformée de Laplace.
L’exemple d’un mur soumis à une température constante Tp sur ses deux faces est considéré dans ce TP. La solution en x = L /2 dans le domaine de LAPLACE appelée fonction de transfert permet de déduire des fonctions de transfert approchées d’ordre 1 et 2 facilement inversibles. Ceci conduit à des solutions approchées du problème qui sont comparées à la solution analytique exacte.
II. POSITION DE PROBLEME
On considère l’équation de la chaleur dans un mur d’épaisseur e et de diffusivitéa :
Les deux faces du mur sont soumises à une température initiale f(t) de sorte que l’on peut écrire :
On supposera enfin que le mur est à une température initiale T(x,0) = 0
III. FORMULATION DANS LE DOMAINE DE LAPLACE
On utilise la transformation de LAPLACE pour résoudre ce problème. Désignons respectivement par q (x,p) et F(p), les transformations de Laplace de T(x,t) et de f(t)
Effectuons à présent la transformation de Laplace sur les deux membres de l’équation :
Les propriétés de linéarité de l’opérateur L (qui désigne la transformation de Laplace ) permettent d’écrire :
Or:
Sachant que T(x,0)=0 on a donc :
Soit:
Les conditions aux limites s’écrivent quant à elles :
La résolution de l’équation différentielle du second ordre précédente conduit à :
Où les deux constantes d’intégration C1 et C2 sont déterminées par les conditions aux limites. On a en effet :
Soit:
Ou encore:
Or:
C2 = F(p) - C1
Donc:
La transformée de la température s’écrit donc :
Soit:
En introduisant la fonction Sinus hyperbolique, ceci s’écrit :
IV. FONCTION DE TRANSFERT
On définit la fonction de transfert H(x,p) à une abscisse x, comme étant le rapport dans le domaine de Laplace de la réponse q (x,p) à la sollicitation F(p).
Elle s’écrit donc :
La transformation inverse de la fonction de transfert conduit à la réponse impulsionnelle du système linéaire considéré.
La réponse à un échelon on à toute autre excitation variable s’en déduit par intégration (théorème de DUHAMEL).
V. FONCTIONS DE TRANSFERT APPROCHEES
La fonction de transfert précédente est difficilement inversible. Nous allons donc chercher à l’approcher en utilisant le fait que la fonction sinh est développable en série entière toujours convergente :
En remplaçant dans l’expression de la fonction de transfert, le sinus hyperbolique par les deux premiers termes de son développement, on obtient la fonction de transfert approchée à l’ordre un suivante :
La fonction de transfert approchée à l’ordre un s’écrira donc au centre du mur ( x = e/2 ) :
En considérant les trois premiers termes du développement en série du sinus hyperbolique, on montre de la même manière que la fonction de transfert approchée à l’ordre deux au centre du mur (x=e/2) s’écrit :
VI. REPONSES IMPULSIONNELLES ET REPONSES INDICIELLES
Les fonctions de transfert approchées précédentes s’écrivent sous la forme générale :
Afin de pouvoir effectuer facilement le retour vers l’espace physique, nous allons effectuer une division polynomiale suivant les puissances croissantes de p et écrire :
La fonction de transfert s’écrira donc sous la forme plus classique :
où Rm est un polynôme du premier ou du second degré.
On montre alors que :
ce qui s’écrit aussi :
La réponse impulsionnelle à l’ordre un est donc simplement :
La réponse indicielle ou réponse à un échelon s’en déduit par simple intégration :
A l’ordre deux, la division polynomiale conduit à :
soit encore :
Les pôles de cette fonction de transfert d’ordre deux sont :
soit :
On peut donc écrire :
La réponse impulsionnelle à l’ordre deux s’en déduit :
La réponse à un échelon sera donc ici :
VII. DEROULEMENT DU TP
A. Développements analytiques
Ø Déterminer les fonctions de transfert approchées d’ordre un et deux en fonction de x.
Ø Calculer les pôles de ces fonctions en fonction de x.
v En déduire le domaine de validité du modèle d’ordre deux.
B. Expérience numériques
Ø Effectuer des simulations pour différents matériaux et différentes épaisseurs
Ø Effectuer une simulation en utilisant des grandeurs adimensionnelles
v En déduire le domaine temporel de validité de ces solutions approchées en utilisant le nombre de FOURIER
v Quelle épaisseur doit avoir une plaque de cuivre pour que la solution d’ordre deux soit correcte à partir de 4 secondes ?
v Quelle épaisseur doit avoir une plaque de béton pour que la solution d’ordre deux soit correcte à partir de 15 minutes ?