Mai 2001

L'usage de calculatrices électroniques de poche à alimentation autonome, non imprimantes et sans document d'accompagnement, est autorisé pour toutes les épreuves d'admissibilité, sauf pour les épreuves de français et de langues. Cependant, une seule calculatrice à la fois est admise sur la table ou le poste de travail, et aucun échange n'est autorisé entre les candidats.

Le chauffage par induction est le procédé qui permet de transférer une énergie électrique à une pièce métallique dans laquelle les pertes Joule sont transformées en chaleur.

Ce chauffage permet de transmettre l'énergie directement à l'intérieur du matériau sans contact matériel. Il est largement répandu dans l'industrie grâce à divers avantages (puissance massique plus élevée que tout procédé traditionnel, échauffement plus rapide, grande précision de température, etc.). On le trouve ainsi en fusion et dans les traitements thermiques superficiels mais aussi pour les chauffages avant déformation et pour des opérations d'assemblage comme le soudage ou le brasage.

Le problème étudié ici porte sur le chauffage avant déformation. Ceci concerne le réchauffage avant laminage ou le chauffage avant formage.

Les produits considérés sont de longs parallélépipèdes d'acier dont les dimensions typiques peuvent être de mètres. Un bobinage inducteur (solénoïde) entoure la barre selon le schéma de principe de la figure Afig1, Dans certains systèmes de chauffage la barre est introduite complètement dans l'inducteur et chauffée entièrement. Dans d'autres elle pénètre progressivement et n'est chauffée que sur une partie de sa longueur seulement.

La maîtrise du chauffage est complexe car elle met en jeu au sein du milieu conducteur chauffé différents phénomènes physiques partiellement couplés :

Ce problème propose sous certaines hypothèses une modélisation des différents phénomènes impliqués. Il envisage pour chaque partie (électromagnétisme, thermique, mécanique) une étude basée sur la résolution analytique des équations physiques en cause.

Les 3 parties A, B, C sont indépendantes. Elles seront rédigées sur des copies indépendantes.

Notations
: champ magnétique
: induction magnétique
: densité de courant
: perméabilité magnétique
: conductivité électrique
f : fréquence
période électrique
pulsation ( )
On rappelle la formule d'analyse vectorielle donnant le rotationnel d'un vecteur ;

La distribution des champs électriques est régie par les équations de Maxwell (avec le courant de déplacement négligeable):

On ajoute à ces deux équations les lois de comportement suivantes au sein du milieu:

L'étude s'effectue en régime sinusoïdal où toutes les grandeurs vectorielles ont des variations harmoniques. On adopte la forme complexe:

Compte tenu des dimensions de la barre, on considère que le système est infini dans les directions y et et que tous les phénomènes sont invariants en et .

Dans cette étude préliminaire la pénétration des champs à l'intérieur du métal est faible ; l'épaisseur de la barre n'intervient pas. Nous assimilons alors le domaine conducteur à un demi-espace (correspondant à ).

On suppose que le champ magnétique est dirigé selon y et ne dépend que de x :

A la surface ( ), le champ magnétique est imposé par les courants inducteurs:

On note:

A.1.1.) Montrer à partir des équations de Maxwell que l'amplitude complexe vérifie:

A.1.2.) On introduit . Montrer que l'équation vérifiée par H est:

A.1.3.) Donner l'expression générale de la solution et préciser les conditions aux limites vérifiées en et .
A.1.4.) Déduire de A.1.2.) et A.1.3.) la solution du problème électromagnétique en fonction de et .
A.1.5.) Donner l'expression de la densité de courant .
A.1.6.) Quelle est la dimension de ? Justifier l'appellation retenue pour : «épaisseur de peau» ou «profondeur de pénétration».
A.1.7.) Déterminer la densité volumique de puissance dissipée :

A.1.8.) Quelle est la moyenne temporelle sur une période de (pertes Joule):

A.1.9.) Montrer que la puissance totale dissipée dans l'ensemble de la barre, par mètre carré de surface (en y et z), est:

Vérifier l’homogénéité de cette expression.

Nous supposons dans cette partie que le milieu conducteur est assimilé à une plaque, infinie en y et z , de largeur 2a. L'origine est prise au milieu de la plaque ( ). Sur chacune des deux faces le champ magnétique est imposé :

A.2.1.) Montrer en exploitant les symétries du problème que .

Dans les questions suivantes on utilisera la fonction , qui se manipule algébriquement comme et par suite on utilisera les fonctions et .
A.2.2.) Vérifier l'expression du champ magnétique :

A.2.3.) Déterminer l'expression de la densité de courant .
A.2.4.) Déterminer l'expression de la densité de puissance volumique .
A.2.5.) Quelle est la moyenne temporelle sur une période de :

A.2.6.) Donner un équivalent simple de lorsque .
A.2.7.) Donner un équivalent simple de lorsque .
A.2.8.) Déterminer la puissance totale dissipée dans la plaque, par unité de surface (en y et z), et montrer qu'elle peut s'écrire:

A.2.9.) On considère que .

Quelles sont les valeurs de lorsque a vaut respectivement: ?
A.2.10.) Montrer que la puissance dissipée a deux équivalents simples correspondant respectivement à et .
A.2.11.) De façon générale, énoncer d’après ce qui précède un (ou des) critère(s) permettant de choisir une fréquence de travail si on souhaite un chauffage efficace en volume de la barre.

Notations

T : température
: température ambiante
: élévation de la température par rapport à la température ambiante:
k : conductivité thermique
: masse volumique
: chaleur spécifique
densité surfacique du flux de chaleur:
q: densité volumique des sources de chaleur.
On cherche à déterminer l'évolution de la température au sein de la barre en fonction du temps de chauffage. On considère dans tout ce qui suit, comme en A, que la plaque est infinie dans les directions y et et que la température n'est fonction que de et de . La distribution de est alors donnée par l’équation de la chaleur dans laquelle la densité des sources de chaleur est la densité de puissance électromagnétique induite au sein de la barre (pertes Joule):

avec

Les phénomènes électromagnétiques et thermiques possèdent deux échelles de temps différentes: dynamique « rapide » pour l'électromagnétisme (fréquence de quelques dizaines de Hertz) et dynamique «lente» pour la thermique (temps de chauffage de plusieurs dizaines de minutes). Il est alors possible de mener une analyse partiellement découplée : on étudie l'évolution thermique du problème en supposant qu’à chaque instant le régime périodique des phénomènes électriques est établi. On fait alors la moyenne des pertes Joule sur une période électrique pour obtenir la puissance volumique injectée. On cherchera ainsi à résoudre :

où est la moyenne temporelle de .
On suppose qu'au début du chauffage la barre est à température ambiante .
B. 1 Cas d'une plaque large devant la profondeur de pénétration a On se place dans la situation . On considère alors que la puissance est dissipée dans une couche superficielle et localisée en et .
B.1.1.) La puissance totale dissipée dans l'ensemble de la plaque est donnée par Aeq 15.

Montrer que le chauffage de la plaque peut se représenter par une densité de flux thermique entrant sur les faces et de valeur:

B.1.2.) Montrer, en exploitant la symétrie du problème que dans ces conditions est solution de l'équation:

où D est la diffusivité thermique ( ), avec les conditions aux limites :

et la condition initiale :

B.1.3.) Montrer que la fonction définie par:

vérifie les équations Beq 5 et Beq 6 de B.1.2.) .
B.1.4.) Afin de résoudre le problème défini en B.1.2.) on effectue un changement d'inconnue en posant: . Montrer que est solution du problème:

B.1.5.) On envisage une résolution du problème défini en B.1.4.) par variables séparées. On pose . Montrer que les fonctions u et v vérifient l'égalité :

où K est une constante réelle.
B.1.6.) On pose où est réel. Montrer en utilisant les conditions aux limites que est nécessairement de la forme:

où n est entier relatif
B.1.7.) Montrer que la solution peut s'écrire:

où les coefficients sont des réels.
B.1.8.) Calculer pour en exploitant Beq 13 les quantités:

En déduire les coefficients et donner l'expression générale de la solution du problème défini en B.1.2.) .
B.1.9.) On considère que .

Donner une estimation du temps au bout duquel on peut considérer que l'évolution temporelle de est voisine de .
B.1.10.) On suppose que t est suffisamment grand pour considérer que . Au bout de combien de temps la température à la surface de la barre atteint-elle 1000 K ? Quelle est alors l'écart de température entre la surface et le centre de la plaque?

On cherche à résoudre l'équation:

lorsque est donnée par :

Les conditions aux limites sont données par:

où h est une constante traduisant la convection thermique ;

On suppose que . On note la solution correspondante.
B.2.1.1 ): Montrer que la solution du problème est alors donnée par :

où .
B.2.1.2 ) Quelle est la limite de lorsque a tend vers 0 . Commenter.
B.2.1.3 ) Quelle est la limite de lorsque a tend vers . Commenter.

B.2.2.1 ) Afin de résoudre le problème on effectue le changement de fonction inconnue :

Montrer que est solution de l'équation :

avec les conditions aux limites:

et la condition initiale:

B.2.2.2 ) On envisage une résolution du problème défini en B.2.2.1 ) par variables séparées.

On pose . Montrer que les fonctions u et v vérifient l'égalité:

où K est une constante réelle.
B.2.2.3 ) On pose où est réel. Montrer en utilisant les conditions aux limites que les valeurs sont solutions de:

B.2.2.4 ) On note la suite des solutions positives de Beq 25. Donner une approximation de pour .
B.2.2.5 ) Montrer que la solution peut s 'écrire :

où les coefficients sont des réels.
B.2.2.6 ) Déduire l' expression générale de .
B.2.2.7 ) Afin d'obtenir une évaluation de la montée en température de la barre, on limite le développement précédent au premier terme ( ).

On considère : . En vous aidant éventuellement d'un graphique donner une estimation de dans les deux cas et .
B.2.2.8 ) On considère: ; ; .

Quelle est la température atteinte au centre de la barre au bout de 3 minutes lorsque pour ? Comparer au résultat obtenu en B.1.10.) .

Dans cette partie on étudiera comment le barreau se dilate sans contrainte ou reste sans dilatation sous contrainte. L'extrémité de la barre d'abscisse sera considérée comme fixe par rapport au repère absolu.

La barre est de longueur L, de section carrée de côté b, Elle sera, ou non, chauffée à son extrémité gauche (dispositif électromagnétique placé en ). Nous proposerons un modèle qui traduit la déformation de cette barre sous l'effet d'actions mécaniques et thermiques, Ce modèle pourra être statique (pas d'accélération prise en compte) ou dynamique (accélération prise en compte).

On utilisera le repère absolu ( ), où est centre de gravité de la section d'abscisse de la barre et les axes sont définis sur la figure Cfigl.

Pour étudier cette barre, on utilise le modèle géométrique présenté sur la figure Cfig2.

La barre est représentée par sa ligne moyenne qui relie les centres de gravité des sections successives de la barre.

On note la projection suivant l'axe du déplacement du point , on suppose donc que le déplacement est le même pour tous les points d'une section d'abscisse x et qu'il est de la forme: .

Le torseur { } en des efforts intérieurs à la barre est défini comme le torseur des efforts exercés par la partie de la barre définie par les abscisses sur la partie de la barre définie par les abscisses . Ce torseur { } a pour éléments de réduction en une résultante et un moment .

On note la masse totale de la barre où p est la masse linéique (masse par unité de longueur) supposée constante.

On note l'élévation de température supposée ici connue (différence entre la température à l'instant t et la température ambiante à l'instant ) au point . On note et les deux constantes qui caractérisent le comportement de la barre.

On utilisera (sans démonstration) la modélisation suivante:

La première équation traduit, en projection sur l'axe , le principe fondamental de la dynamique.
Dans cette équation, représente la projection sur de l’accélération du point . La deuxième équation traduit la relation en tout point entre le déplacement. , l’effort interne et l'élévation de température .

Dans la suite, nous étudierons des problèmes dont les conditions aux limites seront définies comme suit:

Problème de Type I: si est la projection sur l'axe du déplacement imposé au point en et si est la projection sur l'axe de l'effort imposé au point en .

Problème de Type II: si est la projection sur l'axe du déplacement imposé au point en et si est la projection sur l'axe du déplacement imposé au point en .

Dans cette partie, on considère que la variation de température est nulle : il n'y " donc pas d'effet thermique".
C.2.1.) En utilisant les équations Ceq 1 et Ceq 2, écrire l’équation (Ceq3) qui relie et .
C.2.2.) Nous allons rechercher une interprétation de cette équation en considérant le système continu comme la limite d'un système discret décrit sur la figure Cfig3.

Nous adoptons la modélisation suivante: la barre est représentée comme un ensemble de N masses , identiques indicées par i. On notera le déplacement de chacune de ces masses sur l'axe . Elles sont reliées par ressorts identiques de raideur et de longueur à vide .

Écrire le principe fondamental de la dynamique pour chacune des masses intérieures avec i de 2 à (équation ).
C.2.3.) La distance 1 étant considérée comme petite, écrire le développement, au deuxième ordre, qui exprime (avec de 2 à ) en fonction de et de ses dérivées par rapport à (équation Ceq5).
C.2.4.) La distance 1 étant considérée comme petite, écrire le développement, au deuxième ordre, qui exprime (avec de 2 à ) en fonction de et de ses dérivées par rapport à (équation Ceq6).
C.2.5.) En vous inspirant de la forme de l'équation Ceq3, écrire l'équation Ceq7 à partir des équations Ceq4, Ceq5 et Ceq6. Exprimer les valeurs de et de k en fonction des données initiales du problème.
C.2.6.) On envisage de résoudre cette équation Ceq3 en utilisant la même méthode de séparation de variables que dans la partie Étude Thermique (B.1.5.). On pose . Montrer que les fonctions et vérifient l'égalité:

C.2.7.) On admet queK et on poseK où est réel. On adopte les conditions aux limites suivantes (problème de Type I): la barre est fixée au point d'abscisse , elle est libre de tout effort au point d'abscisse . Donner l'expression de .
C.2.8.) Déduire de ce qui précède que la solution générale de l’équation (Ceq3) est de la forme:

C.2.9.) Exprimer en fonction de .
C.2.10.) Donner l'expression finale de et sont alors déterminés, ce qui n'est pas demandé, en considérant les conditions initiales (à ) du mouvement. Cette expression est celle qui traduit les vibrations libres de la barre en traction/compression.

On se place dans le cas où il n'y a pas d'élévation de température et où les accélérations sont nulles.
C.3.1.) Écrire l'équation Ceq8 vérifiée par .
C.3.2.) On applique un effort à chaque extrémité de la barre en et en . La barre est-elle en équilibre? Trouver l'expression de dans la barre en fonction de .