Simulation Multiphysique
Thermique, mecanique, electrique et interpretation
Objectifs du Cours
- Comprendre les principes fondamentaux de la simulation numerique appliquee aux phenomenes physiques couples
- Maitriser les lois constitutives essentielles en thermique, mecanique (RDM) et electrique, avec leurs unites SI
- Analyser et interpreter les resultats d'une simulation multiphysique pour valider un concept
- Relier la theorie a des applications industrielles concretes pour anticiper les comportements reels
Introduction
Dans l'industrie moderne, la conception et l'optimisation des produits (voitures, smartphones, eoliennes, batiments) reposent de moins en moins sur la seule fabrication de prototypes physiques, couteux et longs a realiser. La simulation numerique est devenue un outil indispensable.
La simulation multiphysique va plus loin : elle etudie les interactions entre plusieurs phenomenes physiques simultanes. Par exemple, la dilatation thermique d'un composant electronique (couplage thermique-mecanique) ou l'echauffement d'un moteur par effet Joule (couplage electrique-thermique).
Cette approche permet de predire avec precision les performances, d'identifier les points critiques et de reduire les risques, accelerant ainsi l'innovation tout en maitrisant les couts.
1. Simulation Thermique
La simulation thermique vise a predire la repartition des temperatures et les flux de chaleur dans un systeme.
Transferts thermiques
L'energie thermique se propage de trois manieres :
Conduction
Transfert au sein d'un materiau solide ou fluide immobile, par agitation moleculaire. Mode principal dans les solides.
Convection
Transfert entre une surface et un fluide en mouvement. Peut etre naturelle ou forcee (ventilateur).
Rayonnement
Transfert par ondes electromagnetiques (infrarouges). Ne necessite pas de milieu materiel.
Resistance thermique Rth et flux phi
Rth = e / (lambda . S) phi = DeltaT / Rth
- phi : Flux thermique ou puissance thermique (en watt, W)
- DeltaT : Difference de temperature (en kelvin, K ou degres C)
- e : Epaisseur de la paroi (en metre, m)
- lambda : Conductivite thermique du materiau (en W.m-1.K-1)
- S : Surface d'echange (en metre carre, m2)
Loi de Fourier
phi = - lambda . S . (dT/dx)
Ou dT/dx est le gradient de temperature (en K.m-1). Le signe negatif indique que la chaleur circule des zones chaudes vers les zones froides.
Exemple : Dissipateur thermique pour composant electronique
Problematique : Un microprocesseur dissipe P = 95 W. Il ne doit pas depasser T_jmax = 85 degres C. Temperature ambiante = 25 degres C.
- Calcul de la resistance thermique necessaire :
Rth_tot = (T_j - T_amb) / P = (85-25) / 95 = 0.63 K/W- La simulation determine la repartition de cette resistance entre la pate thermique, le dissipateur et la convection forcee par un ventilateur.
- Choix des materiaux : Cuivre (lambda = 400 W.m-1.K-1) ou Aluminium (lambda = 237 W.m-1.K-1)
2. Simulation Mecanique (RDM)
La Resistance Des Materiaux (RDM) etudie la reponse des structures aux efforts qui leur sont appliques.
Contrainte et deformation
Contrainte normale sigma
sigma = F / S
Unite : le pascal (Pa), 1 Pa = 1 N.m-2
Deformation epsilon
epsilon = DeltaL / L0
Sans unite (ou m/m)
Loi de Hooke (comportement elastique lineaire)
sigma = E . epsilon
Ou E est le module d'Young (en Pa). Il caracterise la rigidite du materiau.
Types de sollicitations elementaires
Traction/Compression
Forces axiales tendant a allonger ou ecraser la piece. Contrainte uniforme sigma = F/S
Flexion
Moments flechissants provoquant une courbure. Contrainte variable, maximale en fibre extreme.
Torsion
Moments de torsion autour de l'axe longitudinal. Contraintes de cisaillement.
Coefficient de securite
s = R_e / sigma_max_calculee
Ou R_e est la limite elastique. Un coefficient de 2 a 10 est courant selon les risques.
Exemple : Dimensionnement d'un longeron de chassis automobile
Etapes de simulation mecanique :
- Modelisation : Creation geometrique 3D de la poutre, application des materiaux (acier haute resistance)
- Maillage : Decoupage en petits elements finis pour le calcul
- Conditions aux limites : Fixation des points d'ancrage (blocage des deplacements)
- Chargement : Application des forces (poids de la carrosserie, moteur, passagers)
- Calcul et post-traitement : Champs de contraintes sigma, deformations, zones critiques
- Optimisation : Modification de la forme (profils en I, en C), epaisseur, ou materiaux composites
3. Simulation Electrique et Interpretation
Lois de Kirchhoff
Loi des noeuds
Somme(I_entrant) = Somme(I_sortant)
La somme algebrique des intensites en un noeud est nulle
Loi des mailles
Somme(U_k) = 0
La somme algebrique des tensions dans une maille fermee est nulle
Impedances en regime sinusoidal
En regime alternatif, la resistance R s'etend au concept d'impedance Z (en ohms, Ohm), qui oppose le passage du courant.
Resistor
Z_R = RBobine (inductance L)
Z_L = j.L.omegadephase I de +pi/2
Condensateur (capacite C)
Z_C = 1/(j.C.omega)dephase I de -pi/2
Ou omega = 2.pi.f est la pulsation (en rad/s) et j le nombre complexe tel que j2 = -1
Puissances en regime sinusoidal
| Puissance | Formule | Unite |
|---|---|---|
| Apparente S | S = U_eff . I_eff | VA (volt-ampere) |
| Active P | P = U_eff . I_eff . cos(phi) | W (watt) |
| Reactive Q | Q = U_eff . I_eff . sin(phi) | var (volt-ampere reactif) |
Analyse frequentielle et filtrage
Elle etudie la reponse d'un circuit (gain, phase) en fonction de la frequence.
Un filtre passe-bas RC attenue les hautes frequences.
f_c = 1 / (2.pi.R.C)
Frequence de coupure
Exemple : Dimensionnement d'une alimentation solaire autonome
Systeme : Panneaux solaires - Regulateur - Batterie - Onduleur - Charges (230V/50Hz)
- Source : Modele du panneau solaire (courant I_ph, tension V_oc) dependant de l'ensoleillement
- Stockage : Modele electrochimique de la batterie (capacite en Ah, tension, resistance interne)
- Conversion : Modele de l'onduleur (rendement eta, forme d'onde de sortie)
- Objectif : Determiner la puissance crete des panneaux (Pc en Wc) et la capacite de la batterie (en Ah) pour garantir 3 jours d'autonomie
Interpretation des resultats de simulation
- Verification de la coherence : Les ordres de grandeur sont-ils plausibles ? (Ex : 500 degres C sur du plastique = erreur !)
- Analyse des points critiques : Identifier les maxima/minima (contrainte max, temperature max, courant max)
- Comparaison avec les criteres de conception : Respect des limites materielles (T_max, R_e, I_max) et des coefficients de securite
- Sensibilite aux parametres : Comment le systeme reagit-il a une variation d'un parametre ?
- Prise de decision : La simulation valide-t-elle le concept ? Sinon, quels parametres modifier ?
Resume - Points Cles a Retenir
- 1La simulation multiphysique est un outil numerique indispensable pour modeliser les interactions entre phenomenes (thermique, mecanique, electrique)
- 2En thermique : loi de Fourier
phi = -lambda.S.dT/dx, resistance et capacite thermiques - 3En mecanique (RDM) : loi de Hooke
sigma = E.epsilon, contraintes, deformations, coefficient de securite - 4En electrique : lois de Kirchhoff, impedances, puissances (P, Q, S) pour dimensionner les systemes
- 5La force de la simulation est d'anticiper les comportements, identifier les points faibles et optimiser avant fabrication
- 6L'interpretation critique des resultats (coherence, extremes, criteres) est essentielle pour une decision eclairee
Mini-Quiz
1. Pour un dissipateur thermique, quelle formule donne le flux de chaleur phi par conduction ?
- a) phi = Rth . DeltaT
- b) phi = DeltaT / Rth
- c) phi = lambda . S . dT/dx
- d) phi = - lambda . S . dT/dx
Reponse : d) C'est la loi de Fourier (le signe - indique le sens du flux)
2. Un longeron en acier (E = 210 GPa) de section S = 10 cm2 subit un effort de traction F = 50 kN. Quelle est la contrainte sigma ?
- a) 5 MPa
- b) 50 MPa
- c) 500 MPa
- d) 5 GPa
Reponse : b) sigma = F/S = 50000 N / 0.001 m2 = 50 MPa
3. Dans un circuit RLC serie, que represente la puissance reactive Q ?
- a) L'energie convertie en chaleur dans la resistance
- b) L'energie echangee entre la source et les elements reactifs (bobine et condensateur)
- c) Le produit U_eff . I_eff
- d) La puissance disponible pour un moteur
Reponse : b) La puissance reactive est echangee entre source et elements reactifs
Pret a pratiquer ?
Testez vos connaissances avec nos exercices interactifs !
