L'actuel investigation numérique du Gurson-Tvergaard-Needleman (GTN) modèle avancé est une extension du travail de Ben Bettaieb et al. (2011). Le modèle a été implémenté à l'aide d'une sous routine (VUMAT) contenu dans le code commerciale d'éléments finis Abaqus/explicit. Le modèle d'endommagement améliore l'original en intégrant les trois mécanismes d'endommagement, la nucléation, la croissance, et la coalescence des cavités. Le modèle d'endommagement intègre les lois de nucléation et de croissance basés sur les phénomènes purement physiques. Ces nouvelles contributions incluant l'influence de l'écrouissage cinématique, ont été validées par les résultats de mesures expérimentales de tomographie à rayon X à haute résolution. Aussi, l'implémentation numérique de l'écrouissage cinématique dans le modèle modifié a contraint de proposer et de réarranger la définition de la triaxialité que l'on trouve habituellement dans la littérature. A coté de cela, un second critère d'initiation à la rupture basé sur l'ultime distance inter-cavités a été inclue afin de localiser et de quantifier avec plus de précision la distribution des déformations peu avant que le matériau ne casse complètement. L'actuel modèle d'endommagement a été appliqué dans des conditions industrielles pour prédire l'évolution de l'endommagement, l'état de contraintes, et l'initiation à la rupture pour différentes géométries de tôles et sur des essais d'emboutissage de tôles minces. / This numerical investigation of an advanced Gurson–Tvergaard–Needleman (GTN) model is an extension of the original work of Ben Bettaieb et al. (2011). The model has been implemented as a user-defined material model subroutine (VUMAT) in the Abaqus/explicit FE code. The current damage model extends the previous version by integrating the three damage mechanisms: nucleation, growth and coalescence of voids. Physically based void nucleation and growth laws are considered, including an effect of the kinematic hardening. These new contributions are based and validated on experimental results provided by high-resolution X-ray absorption tomography measurements. Also, the numerical implementation of the kinematic hardening in this damage extension has obliged to readapt the classical triaxiality definition. Besides, a secondary fracture initiation criterion based on the ultimate average inter-cavities distance has been integrated to localize and quantify with good accuracy the strain distribution just before the material fails apart. The current damage model is applied in industrial conditions to predict the damage evolution, the stress state and the fracture initiation in various tensile thin flat sheet geometries and the cross-die drawing tests.

Dans le cadre de la déposition géologique profonde des déchets radioactifs, l'argile de Boom est choisie comme une des formations hôtes potentielles dans le programme belge. Par sa formation géologique, cette argile est considérée comme un matériau isotrope transversal. En effet, l'anisotropie de ses propriétés hydraulique et thermique a été mise en évidence dans plusieurs études. Il apparait maintenant nécessaire de mener une étude complète sur son comportement anisotrope. Sur le plan expérimental, l'anisotropie des propriétés thermo-hydro-mécaniques de l'argile de Boom a été mise en évidence à l'aide des mesures de la conductivité thermique, de la conductivité hydraulique, et du module de cisaillement. Grâce à ces mesures sur des carottes fraichement forées, l'endommagement dû à l'excavation de la galerie a été identifié (dans la zone près de la galerie « Connecting » à Mol, Belgique) et modélisé par un modèle empirique dont la variable d'endommagement est basée sur le volume des macro-pores. Ensuite, l'anisotropie du comportement en condition K0 a été étudiée en analysant le comportement de trois types d'éprouvettes différents (0°, 45° et 90° au plan du litage) pendant et après la re-saturation dans des cellules oedométriques et triaxiales à plusieurs états de contraintes différents. Enfin, le comportement anisotrope en condition triaxiale a été étudié en réalisant des essais triaxiaux sur trois types d'éprouvettes. Au cours de ces essais, les mesures des vitesses d'ondes Vs et Vp ont été également effectuées à l'aide des éléments piézoélectriques encastrés dans deux embases supérieure et inférieure d'une cellule triaxiale. Ces résultats ont permis de déterminer les paramètres de l'élasticité anisotrope. Sur le plan de modélisation, un modèle anisotrope avec huit paramètres anisotropes (cinq paramètres élastiques et trois paramètres plastiques) a été développé, en utilisant la théorie de Bohler (Boehler et Sawczuk, 1977). Ce modèle est basé sur un modèle élasto-plastique isotrope à deux surfaces de charge validé précédemment pour l'argile de Boom (Hong, 2013). La validation de ce modèle anisotrope a été réalisée grâce aux résultats expérimentaux des essais oedométriques et triaxiaux obtenus dans cette étude / In the program of deep geological radioactive waste disposal in Belgium, Boom Clay has been chosen as one of the potential host rocks. Due to the geological stratification, this stiff clay has been regarded as a transverse isotropic material. The anisotropy of its hydraulic and thermal properties was shown in several studies. It seems necessary now to conduct a more in-depth study on the anisotropic behaviour of Boom Clay. In terms of experimental works, the anisotropy of the thermo-hydro-mechanical properties of Boom Clay was evidenced using measurements of thermal conductivity, hydraulic conductivity and small-strain shear modulus. Through these measurements on freshly cored Boom Clay samples, the damage due to excavation of the gallery was identified (in the zone near the Connecting gallery at Mol, Belgium) and modelled using an empirical model in which the damage variable is based on the volume of macro-pores. Afterwards, the anisotropy behaviour under K0 condition was studied by analysing the behaviour of three types of specimens (0°, 45°, and 90° to the bedding plane) during and after the re-saturation in the triaxial and oedometer cells under different stress states. Finally, the anisotropic behaviour under triaxial condition was investigated through several triaxial tests on three types of specimens. During these tests, the velocity measurements of seismic waves Vs et Vp were performed thanks to the bender elements installed in the upper and bottom bases of a triaxial cell. These results were used to determine the parameters of anisotropic elasticity. In terms of modelling works, an anisotropic model with eight anisotropic parameters (fives elastic parameters and three plastic parameters) was developed using the theory of Boehler (Boehler et Sawczuk, 1977). This model is based on an elasto-plastic isotropic model with two yield surfaces elaborated previously for Boom Clay (Hong, 2013). The validation of this anisotropic model was made based on the results obtained from oedometer and triaxial tests performed in this study

Les procédés de trempe sont largement répandus dans l'industrie en particulier dans le domaine de l'automobile, du nucléaire et de l'aérospatiale car ils ont un impact direct sur la microstructure, les propriétés mécaniques et les contraintes résiduelles de pièces critiques. La trempe est un processus fortement non-linéaire à cause des couplages forts entre la mécanique des fluides, les transferts thermiques aux différentes interfaces, les transformations de phase du solide et l'ébullition du milieu de trempe. Malgré les progrès effectués par la simulation numérique, ce procédé reste extrêmement difficile à modéliser.Dans ce travail, nous proposons le développement d'outils numériques permettant la simulation réaliste de ce procédé à l'échelle industrielle. La mécanique des fluides est simulée en utilisant une méthode d'Elements Finis stabilisés permettant de considérer des écoulements à haut nombre de Reynolds. Les transferts thermiques sont calculés directement sans l'utilisation de coefficients de transferts empiriques, en utilisant le couplage fort entre le fluide et le solide. Nous avons développé un modèle de changement de phase pour l'eau permettant de considérer les différents régimes d'ébullition. Une formulation unifiée des équations de Navier-Stokes, considérant une phase compressible et une phase incompressible a été développée afin de prendre en compte plus précisément la dynamique de la vapeur et de l'eau. Une procédure dynamique d'adaptation anisotrope de maillage, permettant une description plus fine des interfaces et une prise en compte plus précise des caractéristiques des écoulements est utilisée.Des exemples numériques exigeants ainsi qu'une validation expérimentale permettent d'évaluer la précision et la robustesse des outils proposés.Les outils développés permettent ainsi l'optimisation du mode opératoire du procédé, des ressources consommées et servent ainsi d'outils prospectifs pour la conception de produits. / Quenching processes of metals are widely adopted procedures in the industry, in particular automotive, nuclear and aerospace industries, since they have direct impacts on changing mechanical properties, controlling microstructure and releasing residual stresses of critical parts. Quenching is a highly nonlinear process because of the strong coupling between the fluid mechanics, heat transfer at the interface solid-fluid, phase transformation in the metal and boiling. In spite of the maturity and the popularity of numerical formulations, several involved mechanisms are still not well resolved.Therefore we propose a Direct Numerical Simulation of quenching processes at the industrial scale dealing with these phenomena. The fluid mechanics is simulated using a Finite Element Method adapted for high convective flows allowing the use of high stirring velocity in the quenching bath. Heat transfers are computed directly without the use of heat transfer coefficients but using the strong coupling between the fluid and the solid. We use a phase change model for the water that models all boiling regimes. A unified formulation of the Navier-Stokes equations, taking into account a compressible gas and an incompressible liquid is developed to model more accurately the vapor-water dynamics. A dynamic mesh adaptation procedure is used, increasing the resolution in the description of the interfaces and capturing more accurately the features of the flows.We assess the behavior and the accuracy of the proposed formulation in the simulation of time-dependent challenging numerical examples and experimental results.These recent developments enable the optimization of the process in terms of operating conditions, resources consumed and products conception.

Cette thèse est consacrée à l’étude théorique ainsi qu’au traitement numérique de fluides incompressibles à surface libre. La première partie concerne un système d’équations appelé le système de Green–Naghdi. Comme le système de Saint-Venant, il s’agit d’une approximation d’eaux peu-profondes du problème de Zakharov. La différence est que le système de Green–Naghdi est d’un degré plus élevé en ordre d’approximation. C’est pourquoi il contient tous les termes du système de Saint-Venant plus de termes d’ordre trois non-linéairement dispersives. Autrement dit, le système de Green–Naghdi peut être vu comme une perturbation dispersive du système de Saint-Venant. Ce dernier système étant hyperbolique, il entre dans le cadre classique développé pour des systèmes hyperboliques. En particulier, il est entropique (au sense de Lax) et symétrique. On peut donc lui appliquer les résultats d’existence et d’unicité bien connus pour des systèmes hyperboliques. Dans la première partie de ce travail, on généralise la notion de symétrie à une classe plus générale de systèmes contenant le système de Green–Naghdi. Ceci nous permet de symétriser les équations de Green–Naghdi et d’utiliser la symétrie obtenue pour déduire un résultat d’existence globale après avoir ajouté un terme dissipative d’ordre 2 au système. Ceci est fait en adaptant l’approche utilisée dans la littérature pour des systèmes hyperboliques. La deuxième partie de ce travail concerne le traitement numérique des équations de Navier–Stokes à surface libre avec un terme de tension de surface. Ici, la surface libre est modélisée en utilisant la formulation des lignes de niveaux. C’est pourquoi la condition cinématique (condition de l’évolution de surface libre) s’écrit sous la forme d’une équation d’advection satisfaite par la fonction de ligne de niveaux. Cette équation est résolue sur une domaine de calcul contenant strictement le domaine de fluide, sur de petits sous-intervalles du temps. Chaque itération de l’algorithme global correspond donc à l’advection du domaine du fluide sur le sous-intervalle du temps associé et ensuite de résoudre le système de Navier–Stokes discrétisé en temps sur le domaine du fluide. Cette discrétisation en temps est faite par la méthode des caractéristiques. L’outil clé qui nous permet de résoudre ce système uniquement sur le domaine du fluide est l’adaptation de maillage anisotrope. Plus précisément, à chaque itération le maillage est adapté au domaine du fluide tel que l’erreur d’approximation et l’erreur géométrique soient raisonnablement petites au voisinage du domaine du fluide. La résolution du problème discrétisé en temps sur le domaine du fluide est faite par l’algorithme d’Uzawa utilisé dans la cadre de la méthode des éléments finis. Par ailleurs, la condition de glissement de Navier est traité ici en ajoutant un terme de pénalisation à la formulation variationnelle associée. / This thesis is about theoretical study and numerical treatment of some problems raised in incompressible free-surface fluid dynamics. The first part concerns a model called the Green–Naghdi (GN) equations. Similarly to the non linear shallow water system (called also Saint-Venant system), the Green–Naghdi equations is a shallow water approximation of water waves problem. Indeed, GN equation is one order higher in approximation compared to Saint-Venant system. For this reason, it contains all the terms of Saint-Venant system in addition to some non linear third order dispersive terms. In other words, the GN equations is a dispersive perturbation of the Saint-Venant system. The latter system is hyperbolic and fits the general framework developed in the literature for hyperbolic systems. Particularly, it is entropic (in the sense of Lax) and symmertizable. Therefore, we can apply the well-posedness results known for symmetric hyperbolic system. During the first part of this work, we generalize the notion of symmetry to a more general type of equations including the GN system. This lets us to symmetrize the GN equation. Then, we use the suggested symmetric structure to obtain a global existence result for the system with a second order dissipative term by adapting the approach classically used for hyperbolic systems. The second part of this thesis concerns the numerical treatment of the free surface incompressible Navier–Stokes equation with surface tension. We use the level set formulation to represent the fluid free-surface. Thanks to this formulation, the kinematic boundary condition is treated by solving an advection equation satisfied by the level set function. This equation is solved on a computational domain containing the fluid domain over small time subintervals. Each iteration of the algorithm corresponds to the adevction of the fluid domain on a small time subinterval and to solve the time-discretized Navier–Stokes equations only on the fluid domain. The time discretization of the Navier–Stokes equation is done by the characteristic method. Then, the key tool which lets us solve this equation on the fluid domain is the anisotropic mesh adaptation. Indeed, at each iteration the mesh is adapted to the fluid domain such that we get convenient approximation and geometric errors in the vicinity of the fluid domain. This resolution is done using the Uzawa algorithm for a convenient finite element method. The slip boundary conditions are considered by adding a penalization term to the variational formulation associated to the problem.

Le comportement poussant fait référence au phénomène de grande déformation différée et souvent anisotrope observée lors de l'excavation du tunnel en terrain tectonisé. Il est à l'origine de difficultés d'avancement ce qui exige une adaptation de la méthode de creusement et de la conception des soutènements. Le présent travail vise à étudier le comportement des tunnels en terrain poussant en portant une attention particulière à l'anisotropie du massif rocheux par des approches à la fois analytique et numérique. Après un état de l'art sur le creusement des tunnels en terrain poussant, on interprète les données d'auscultation récoltées pendant l'excavation de la descenderie de Saint-Martin-La-Porte dans la cadre du projet Lyon-Turin. Des solutions analytiques pour tunnel creusés en milieu anisotrope sont ensuite développées en prenant en compte la complexité géométrique de la section, l'interaction entre deux tunnels parallèles, l'interaction terrain-soutènement et aussi les grandes déformations. Enfin, un modèle différé anisotrope qui comprend des joints rocheux avec une orientation fixe, imbriqués dans une matrice viscoplastique est proposé et implémenté dans FLAC3D. Les résultats de la simulation numérique réalisée avec ce modèle sont comparés aux mesures de convergence réalisées pendant l'excavation du tunnel. De plus, l'approche numérique a été étendue pour analyser le comportement d'un soutènement déformable en prenant en compte l'effet de l'anisotropie de la masse rocheuse / Squeezing behavior is characterized by large time-dependent and often anisotropic deformation during and well after excavation of tunnel and may lead to tremendous operational difficulties. The present thesis aims to deal with tunneling in squeezing ground with a special emphasis on the anisotropic behavior combining analytical and numerical approach. Following an overview on the squeezing behavior, the attention moves on the interpretation of the data monitored during the excavation of Saint-Martin-La-Porte within the Lyon-Turin railway project. The closed-formed solutions for tunnel excavated in anisotropic ground are then developed considering the complexity of tunnel cross-section, the interaction between two tunnels, the interaction ground/support and also the large strain calculation. Finally, an anisotropic creep model which includes weak planes of specific orientation embedded in a viscoplastic medium is proposed and implemented in FLac3D in order to back analyze the convergence data of Saint-Martin-La-Porte access gallery. The numerical model is also applied to the analysis of the behavior of the yield-control support systems taking account the effect of anisotropy of rock mass

L'Interaction Fluide-Structure (IFS) décrit une classe très générale de problème physique, ce qui explique la nécessité de développer une méthode numérique capable de simuler le problème FSI. Pour cette raison, un solveur IFS est développé qui peut traiter un écoulement de fluide incompressible en interaction avec des structures différente: élastique ou rigide. Dans cet aspect, le solveur peut couvrir une large gamme d'applications.La méthode proposée est développée dans le cadre d'une formulation monolithique dans un contexte Eulérien. Cette méthode consiste à considérer un seul maillage et résoudre un seul système d'équations avec des propriétés matérielles différentes. La fonction distance permet de définir la position et l'interface de tous les objets à l'intérieur du domaine et de fournir les propriétés physiques pour chaque sous-domaine. L'adaptation de maillage anisotrope basé sur la variation de la fonction distance est ensuite appliquée pour assurer une capture précise des discontinuités à l'interface fluide-solide.La formulation monolithique est assurée par l'ajout d'un tenseur supplémentaire dans les équations de Navier-Stokes. Ce tenseur provient de la présence de la structure dans le fluide. Le système est résolu en utilisant une méthode élément fini et stabilisé suivant la formulation variationnelle multiéchelle. Cette formulation consiste à décomposer les champs de vitesse et pression en grande et petite échelles. La particularité de l'approche proposée réside dans l'enrichissement du tenseur de l'extra contraint.La première application est la simulation IFS avec un corps rigide. Le corps rigide est décrit en imposant une valeur nul du tenseur des déformations, et le mouvement est obtenu par la résolution du mouvement de corps rigide. Nous évaluons le comportement et la précision de la formulation proposée dans la simulation des exemples 2D et 3D. Les résultats sont comparés avec la littérature et montrent que la méthode développée est stable et précise.La seconde application est la simulation IFS avec un corps élastique. Dans ce cas, une équation supplémentaire est ajoutée au système précédent qui permet de résoudre le champ de déplacement. Et la contrainte de rigidité est remplacée par la loi de comportement du corps élastique. La déformation et le mouvement du corps élastique sont réalisés en résolvant l'équation de convection de la Level-Set. Nous illustrons la flexibilité de la formulation proposée par des exemples 2D. / Numerical simulations of fluid-structure interaction (FSI) are of first interest in numerous industrial problems: aeronautics, heat treatments, aerodynamic, bioengineering... Because of the high complexity of such problems, analytical study is in general not sufficient to understand and solve them. FSI simulations are then nowadays the focus of numerous investigations, and various approaches are proposed to treat them. We propose in this thesis a novel monolithic approach to deal with the interaction between an incompressible fluid flow and rigid/ elastic material. This method consists in considering a single grid and solving one set of equations with different material properties. A distance function enables to define the position and the interface of any objects with complex shapes inside the volume and to provide heterogeneous physical properties for each subdomain. Different anisotropic mesh adaptation algorithms based on the variations of the distance function or on using error estimators are used to ensure an accurate capture of the discontinuities at the fluid-solid interface. The monolithic formulation is insured by adding an extra-stress tensor in the Navier-Stokes equations coming from the presence of the structure in the fluid. The system is then solved using a finite element Variational MultiScale (VMS) method, which consists of decomposition, for both the velocity and the pressure fields, into coarse/resolved scales and fine/unresolved scales. The distinctive feature of the proposed approach resides in the efficient enrichment of the extra constraint. In the first part of the thesis, we use the proposed approach to assess its accuracy and ability to deal with fluid-rigid interaction. The rigid body is prescribed under the constraint of imposing the nullity of the strain tensor, and its movement is achieved by solving the rigid body motion. Several test case, in 2D and 3D with simple and complex geometries are presented. Results are compared with existing ones in the literature showing good stability and accuracy on unstructured and adapted meshes. In the second, we present different routes and an extension of the approach to deal with elastic body. In this case, an additional equation is added to the previous system to solve the displacement field. And the rigidity constraint is replaced with a corresponding behaviour law of the material. The elastic deformation and motion are captured using a convected level-set method. We present several 2D numerical tests, which is considered as classical benchmarks in the literature, and discuss their results.

L’objectif de la thèse est de développer un dispositif de micro positionnement intégrant à la fois les actionneurs et les capteurs. Un dispositif a été conçu afin de réaliser des déplacements dans les plans sur une course millimétrique. Le dispositif compact ne nécessite pas de système de guidage additionnel et selon le mode d’utilisation de ces moteurs, il est capable de réaliser des translations dans le plan ou des rotations autour d’un axe perpendiculaire au plan. Le dispositif comprend quatre moteurs électromagnétiques linéaires fixés orthogonalement sur une structure en forme de croix. Chaque moteur consiste en une paire de bobines planes entrelacées fixe et une barre ’aimants mobile. Un capteur de déplacement intégré dans la structure en croix permettant de mesurer le déplacement de celle-ci a été conçu et fabriqué. Ce capteur est constitué d’une tête de mesure à fibres optiques placé face à un réseau en silicium réalisé par des techniques de microfabrication. Afin de minimiser les erreurs d’assemblage, la structure en croix a également été micro fabriquée. Le dispositif est capable de réaliser un déplacement de 10 mm et une rotation de ±11° autour de l’axe perpendiculaire au plan du dispositif. La résolution de déplacement du dispositif est de 1,4 µm avec une précision de 31 nm en boucle fermée. Le dispositif peut également atteindre une vitesse de déplacement de 12 mm/s. / The objective of thesis is to develop an integrated micro positioning system for micro applications. A unique micro positioning system design capable to deliver millimeter level strokes with pre-embedded auto guidance feature in micro application has been realized. The design integrates, a stack of orthogonally arranged four electromagnetic linear motors. Each linear motor consists of a fixed planar electric drive coil and mobile permanent magnet array. The optimal design of the system delivers a small footprint size. In addition, to measure and control the displacement, a high resolution compact optical displacement measurement sensor has been designed and fabricated in silicon material using microfabrication technology. Furthermore, a light weight silicon cross structure was fabricated using dry etching technology to reduce components assembly errors. The device is capable to deliver 10 mm displacement stroke with a rotation of ±11° about an axis perpendicular to the plane of the device. The displacement resolution of the device is 1.4 µm with a precision of 31 nm in closed loop control. The device can realize displacement with a speed of 12 mm/s.

La modélisation numérique des structures bois dans des conditions d’incendie nécessite la connaissance : de la variation des propriétés physiques du bois telles que la conductivité thermique, la chaleur spécifique et la densité en fonction de la température ; de la dégradation thermique du bois au cours des phases de séchage, de pyrolyse et de combustion. En particulier, nous nous sommes intéressés à l’étude du comportement thermomécanique du matériau bois. La loi thermique est décrite par l’équation de la chaleur. Le modèle choisi intègre les trois modes du transfert de chaleur : la conduction, le rayonnement et la convection. La loi mécanique est modélisée dans le cadre de la thermodynamique des processus irréversibles utilisant la notion des variables d’état. Elle tient compte du couplage entre le comportement élastique orthotrope, plastique anisotrope à écrouissage non linéaire isotrope et un endommagement isotrope. L’intégration numérique de la loi mécanique par un schéma implicite itératif combinant la technique du retour radial avec la réduction du nombre des équations est présentée. Le couplage thermomécanique est réalisé, selon l’approche réglementaire de l’Eurocode 5 relatif à la résistance au feu des structures en bois, en appliquant le facteur de réduction Kθ sur la résistance mécanique d’un résineux. Les aspects théoriques et les conditions aux limites associés au modèle thermomécanique sont abordés. L’identification des paramètres du modèle est réalisée sur des données expérimentales obtenues sur des tests réels d’incendie disponibles dans la littérature. À ce titre, plusieurs comparaisons avec différentes applications sont réalisées. Le modèle éléments finis reproduit avec précision la distribution du champ de température dans l’épaisseur des panneaux en bois, la formation du charbon ainsi que l’évolution de la résistance mécanique au cours de l’exposition au feu / Numerical modelling of timber structures in fire conditions requires the knowledge of the variation with temperature of the physical properties of the wood material (the thermal conductivity, the specific heat and the density) in order to take into account the thermal degradation of wood under high temperatures during the drying, pyrolysis and combustion phases, as well as the temperature profiles in the thickness of the surfaces exposed to fire. In particular, this work focusses on the thermomechanical behaviour of timber. The heat transfer analysis is described by the standard equations of heat conduction. It includes the three modes of heat transfer: conduction, radiation and convection. The structural response is modelled within the framework of thermodynamics of irreversible processes using the notion of state variables. It takes into account the coupling between the orthotropic elastic behaviour, the anisotropic plastic behaviour with isotropic nonlinear hardening, and isotropic damage. The numerical integration of the equilibrium equations is carried out with an iterative implicit scheme combining the technique of radial re- turn with the reduction of the number of equations. The thermomechanical coupling is carried out according to the approach recommended by Eurocode 5 for the fire resistance of timber structures by applying the reduction factor Kθ to the strength of a softwood. The theoretical aspects and boundary conditions associated with the thermomechanical model are also discussed. The parameters of the model are identified with experimental data obtained from actual fire tests available in the literature. Several comparative applications are carried out. The finite element model accurately reproduces the distribution of the temperature profile in the thickness of timber planks, the formation of the charred layer, and the evolution of the mechanical resistance during exposure to fire

L’augmentation de la puissance électrique des composants électroniques pose le problème de la dissipation de la chaleur générée. Les boîtiers électriques doivent permettre la dissipation de cette chaleur en conservant une isolation électrique. La solution retenue pour évacuer la chaleur par transfert thermique consiste en matériaux composites dont les renforts par leur structuration vont améliorer la conductivité thermique. Des composites à matrice polymère ont été choisis pour leur aptitude de mise en forme. La conductivité thermique et l’isolation électrique sont assurées par des charges céramiques. Les méthodes d’homogénéisation donnent des pistes d’amélioration du comportement de composites en fonction des propriétés de leurs constituants, de leur géométrie et de leur distribution. Elles fournissent ainsi une formulation optimisée de matériaux répondant à certaines caractéristiques issues de cahiers des charges émanant du partenaire industriel (Institut de Soudure). La conductivité thermique attendue des composites impose une forte fraction volumique de charges pour compenser le caractère isolant de la matrice polymère. Des méthodes d’homogénéisation ont été développées pour prédire la conductivité thermique effective pour de forts taux de charges (supérieur à 20%) et des contrastes élevés de conductivité thermique. La présence d’une interphase provenant d’incompatibilités fortes entre les composants doit également être modélisée / The increase of electronic components in the integrated circuits and the required electrical power set the question of the dissipation of the heat generated. The electrical box must favor the heat dissipation while maintaining electrical insulation. The solution chosen to transfer the heat is to develop composite materials whose reinforcements by their structure will improve the thermal conductivity. Polymer-based composite materials were chosen for their building ability. Thermal conductivity and electrical insulation are insured by ceramic reinforcements. The homogenization methods allow to improve the composites’ design according to the properties of their constituents, their geometry and their distribution. They thus provide an optimized formulation of materials satisfying the characteristics emanating from the industrial partner (‘Institut de Soudure’). The expected thermal conductivity of the composites imposes a high volume fraction of reinforcements to counterbalance the insulating polymer matrix. Homogenization methods have been developed to provide predictions of effective thermal conductivity for high (greater than 20%) reinforcement rates and high thermal conductivity contrasts. The presence of an interphase resulting from strong physico-chemical incompatibilities between the components must also be modeled

Dans ce travail, les propriétés thermiques effectives du milieu fibreux sont déterminées en tenant compte du couplage conduction et rayonnement. Un échantillon numérique fibreux statistiquement homogène composé de deux phases a été généré en empilant des cylindres finis absorbant dans le vide. Ces cylindres sont dispersés selon des fonctions de distribution de la position de leur centre et de leur orientation. L'interpénétration des cylindres est permis. L'extinction, l'absorption et la diffusion sont caractérisées par des fonctions statistiques radiatives qui permettent de savoir si le milieu est Beerien (ou non). Elles sont déterminées précisément à l'aide d'une méthode de Monte Carlo. On montre que la phase gazeuse a un comportement Beerien et que le phase fibreuse a un comportement fortement non Beerien. Le champ de puissance radiative déposée dans le milieu fibreux est calculé en résolvant un modèle qui couple une Équation du Transfert Radiatif Généralisée (ETRG) et une Équation du Transfert radiatif Classique (ETR). Le modèle de conduction thermique est basé sur une méthode de marche aléatoire ne nécessitant aucun maillage. La simulation du mouvement Brownien de marcheurs dans les fibres permet de résoudre l'équation de l'énergie. L'idée de la méthode est de caractériser la température d'un volume élémentaire par une densité de marcheurs, qui peuvent parcourir le milieu. Le problème est gouverné par les conditions aux limites ; Une concentration constante de marcheurs (ou un flux constant) est associée à une température imposée (ou un flux). / In the present work, the effective heat transfer properties of fibrous medium are determined by taking into account a coupling of heat conduction and radiation. A virtual, statistically homogeneous, two-phase fibrous sample has been built by stacking finite absorbing cylinders in vaccum. These cylinders are dispersed according to prescribed distribution functions defining the cylinder positions and orientations. Cylinder overlappings are allowed. Extinction, absorption and scattering are characterised by radiative statistical functions which allow the Beerian behaviour of a medium to be assessed (or not). They are accurately determined with a Monte Carlo method. Whereas the gaseous phase exhibits a Beerian behaviour, the fibre phase is strongly non Beerian. The radiative power field deposited within the fibrous material is calculated by resolving a model which couples a Generalized Radiative Transfer Equation (GRTE) and a classic Radiative Transfer Equation (RTE). The model of conduction transfer is based on a random walk method without meshing. The simulation of Brownian motion of walkers in fibres allows the energy equation to be solved. The idea of the method is to characterize the temperature in an elementary volume by the density of walkers, which roam the medium. The problem is governed by boundary conditions ; A constant concentration of walkers (or a constant flux) is associated with a fixed temperature (or flux).