It has been known for decades that fatigue crack propagation in elastic-plastic media is very sensitive to load history since the nonlinear behavior of the material can have a great influence on propagation rates. However, the raw computation of millions of fatigue cycles with nonlinear material behavior on tridimensional structures would lead to prohibitive calculation times. In this respect, we propose a global model reduction strategy, mixing both the a posteriori and a priori approaches in order to drastically decrease the computational cost of these types of problems. First, the small scale yielding hypothesis is assumed, and an a posteriori model reduction of the plastic behavior of the cracked structure is performed. This reduced model provides incrementally the plastic state in the vicinity of the crack front, from which the instantaneous crack growth rate is inferred. Then an additional a priori model reduction technique is used to accelerate even more the time to solution of the whole problem. This a priori approach consists in building incrementally and without any previous calculations a reduced basis specific to the considered test-case, by extracting information from the evolving displacement field of the structure. Then the displacement solutions of the updated crack geometries are sought as linear combinations of those few basis vectors. The numerical method chosen for this work is the finite element method. Hence, during the propagation the spatial discretization of the model has to be updated to be consistent with the evolving crack front. For this purpose, a specific mesh morphing technique is used, that enables to discretize the evolving model geometry with meshes of the same topology. This morphing method appears to be a key component of the model reduction strategy. Finally, the whole strategy introduced above is embedded inside an adaptive approach, in order to ensure the quality of the results with respect to a given accuracy. The accuracy and the efficiency of this global strategy have been shown through several examples; either in bidimensional and tridimensional cases for model crack propagation, including the industrial example of a helicopter structure.

Dans le domaine de la mise en forme industrielle, des outils de simulation comme le logiciel Pam-stamp 2G permettent entre autres le prototypage et l’optimisation numérique des produits, réduisant ainsi les coûts expérimentaux de mise au point. Les éléments finis de type coques en hypothèse d’état plan de contrainte demeurent les plus utilisés car ils permettent une prise en compte réaliste des déformations majeures de membrane et de flexion. Cependant, de par leur définition, la contrainte normale pouvant apparaitre en cas de compression du flan dans la direction de l’épaisseur ou encore de flexion extrême sur petit rayon est systématiquement omise. De plus, Il existe de nouveaux procédés de mise en forme de tôle, comme le pliage/sertissage/emboutissage avec laminage et/ou écrasement ainsi que l’hydroformage, qui ne peuvent pas être traités avec ces formulations d’éléments coques classiques. L’utilisation de couches d’éléments volumiques est souvent considérée comme une alternative non convenable aux simulations de ces procédés. Outre le très haut coût CPU, s’ajoutent le rendu parfois non réaliste ainsi que la complexité liée à la découpe du maillage du flan. Récemment, des éléments de type “solid-shell” ont été mis à contribution mais requièrent des améliorations quant à leurs lois de comportements. L’objectif d’ESI group consiste en l’élaboration, l’évaluation, l’implémentation et la validation industrielle d’une nouvelle formulation d’élément fini. Cet élément devra permettre la gestion d’une éventuelle variation d’épaisseur avec prise en compte réaliste de la contrainte normale, tout en assurant des résultats dignes de ceux d’une coque conventionnelle en flexion. Une nouvelle formulation de type coque 3D est ainsi d’abord proposée. Des éléments coques triangulaires et quadrangulaire en théorie de Mindlin et de Kirchhoff sont utilisés. Cette approche est d’abord développée dans un solveur quasi-statique implicite de l’INSA de Lyon pour validation numérique sur cas académiques linéaires et non linéaires de référence. Une validation expérimentale sur opération de pliage dépliage est également réalisée. Aux vues des non linéarités dues au contact avec frottement, grandes déformations et grands déplacements posant des soucis de convergence en implicite au cours de la simulation d’opérations de mise en forme, la méthode est ensuite développée dans le solveur explicite de l’INSA de Lyon. Seuls les éléments en théorie de Mindlin sont considérés. Les particularités liées à cette méthode de résolution dynamique comme la matrice de masse, le pas de temps critique et l’optimisation du CPU sont traitées. Une nouvelle méthode de contact dédiée aux opérations de mise en forme impliquant du laminage et/ou de l’amincissement est également proposée. Elle permet une transition automatique d’éléments standards vers des éléments coque "3D", palliant ainsi le remaillage. Des essais de mise en forme en U avec ou sans laminage apportent une validation expérimentale concernant le retour élastique. Une fois validées, les techniques et formulations les plus abouties sont implémentées dans le code industriel dédié à la mise en forme Pam-stamp 2G v2011. Après une vérification sur tests de référence, des applications sont enfin menées sur cas critiques inspirés de procédés industriels complexes et nécessitant essentiellement une loi de comportement 3D. / In the sheet metal forming industry, shell elements in plane-stress assumption are employed, as they perform quite well in simulating the major membrane and flexural large deformations involved. However, the normal stress, caused by compression along thickness direction of the blank or local high bending over very small radii, is hence systematically omitted. Besides, when it comes to unusual and challenging processes such as hydro-forming, thinning/thickening, forming with ironing, bottoming and so on, makeshift solutions such as layers of 3D solid hexahedrons or even recent “solid-shell” elements are no longer appropriate. An innovative 3D finite element formulation methodology overcoming the overcoming the plane-stress definition of classification shell elements, while keeping their very good bending assets is first proposed in this work. The method basically consists in adding a central node endowed with two degrees of freedom at the element center. These two extra translations normal to the element mid-plane give a new quadratic displacement field along the shell normal direction. A derivative normal strain can hence be expressed and a linear normal stress comes via a full 3D constructive law. A very pioneering contact technique, dedicated to forming processes with ironing, thinning/bottoming operations and allowing a usual-to-enhanced automatic element switch is developed as well. Once widely assessed, most interesting achievements are implemented in the dynamic explicit industrial code Pam-stamp 2 G v2011 and evaluated over critical industrial forming processes that require essentially a full 3D strain-stress behavior.

Les alliages d'aluminium utilisés dans les structures aéronautiques (fuselage, voilure) sont soumis à des chargements cycliques, faisant de la fatigue l'un des facteurs dimensionnant. Dans cette thèse, les mécanismes d'amorçage de ces fissures de fatigue – au niveau des particules intermétalliques – et de micropropagation sont étudiés expérimentalement et numériquement sur les alliages 2050-T8 et 7050-T7451. Les analyses des premiers chapitres portent sur la description des particules intermétalliques qui sont les plus susceptibles de donner lieu à une amorce de fissure dans la matrice d'aluminium. Les effets de la nature des particules et de leur taille sont quantifiés. La proximité entre les particules intermétalliques et les pores y est décrite. Cette analyse expérimentale fait intervenir des observations de surface en microscopie électronique à balayage (MEB), ainsi que des caractérisations tridimensionnelles (3D) réalisées à l'aide de la tomographie par rayonnement synchrotron. Dans les chapitres suivants, les analyses traitent de la prévision des chemins de fissuration à l'échelle d'un grain. Elles s'appuient sur des observations expérimentales de surface et des simulations de plasticité cristalline 3D pour comprendre les chemins de fissuration. Un modèle d'endommagement tenant compte de la cristallographie est alors proposé pour simuler la propagation de fissure par éléments finis. / Aluminium alloys used for aerospace structures (wing, fuselage) are subjected to cyclic loading. Fatigue properties of such alloys are therefore taken into account for the design of such parts. In this thesis, initiation mechanisms of fatigue cracks – near intermetallic particles – and micropropagation are studied experimentally and numerically on alloys 2050-T8 and 7050-T7451. In the first chapters, the analysis focuses on intermetallic particles which are most prone to initiate a fatigue crack in the aluminium matrix. The effects of the nature of particles as well as their size are quantified. The proximity between intermetallic particles and pores is described. This experimental analysis use surface observations obtained with a scanning electron microscope (SEM), and three-dimensional (3D) characterizations using synchrotron tomography. In the last chapters, analysis are oriented towards the prediction of crack paths at the grain size. They rely on surface experimental observations and 3D crystal plasticity modelling in order to understand crack paths. A damage model taking into account crystallography is proposed to simulate crack propagation using the finite element method.

Il est connu depuis des décennies que la propagation des fissures de fatigue dans les matériaux élastoplastiques est très sensible à l’histoire du chargement car le comportement non-linéaire du matériau peut avoir une grande influence sur les vitesses de propagation. Cependant, le calcul brut de millions de cycles de fatigue avec des comportements matériaux non-linéaires sur des structures tridimensionnelles réalistes conduirait à des temps de calcul prohibitifs. Ainsi, nous proposons de coupler deux approches de réduction de modèle a priori et a posteriori, afin de diminuer considérablement le coût de calcul de ce type de problèmes. Tout d’abord, considérant l’hypothèse de plasticité confinée, une stratégie de réduction de modèle a posteriori du comportement plastique de la structure fissurée est proposée. Le modèle réduit ainsi obtenu fournit incrémentalement l’état plastique autour du front de fissure, duquel est déduite la vitesse instantanée de la fissure. De plus, une seconde approche de réduction de modèle, a priori cette fois, est aussi mise en place afin d’accélérer encore plus les temps de résolution du problème global. Cette approche a priori consiste à construire incrémentalement —et sans calculs préalables— une base réduite spécifique à chaque cas-test, en extrayant de l’information des champs de déplacement de la structure au cours du temps et pendant la propagation éventuelle de la fissure. Ainsi, les champs de déplacement solutions de la géométrie fissurée réactualisée sont vus comme une combinaison linéaire de cette base réduite de vecteurs. La méthode numérique considérée ici est la méthode des éléments finis. De fait, pendant la propagation de la fissure, la discrétisation spatiale du modèle doit être réactualisée afin d’être conforme avec le front de la fissure. Dans ce but, une technique spécifique de déformation de maillage est utilisée, et permet de discrétiser la géométrie variable du modèle avec des maillages de même topologie. Cette technique de déformation de maillage apparaît comme une étape clé de la stratégie de réduction de modèle. Finalement, une approche adaptative est construite autour de cette stratégie. Elle permet de garantir la qualité des résultats obtenus par rapport à un critère de précision donné. La précision et l’efficacité de cette stratégie globale sont démontrées à travers de nombreux exemples bidimensionnels et tridimensionnels dans le cadre de propagation de fissure en model, de même que pour un exemple industriel d’une pièce fissurée d’hélicoptère. / It has been known for decades that fatigue crack propagation in elastic-plastic media is very sensitive to load history since the nonlinear behavior of the material can have a great influence on propagation rates. However, the raw computation of millions of fatigue cycles with nonlinear material behavior on tridimensional structures would lead to prohibitive calculation times. In this respect, we propose a global model reduction strategy, mixing both the a posteriori and a priori approaches in order to drastically decrease the computational cost of these types of problems. First, the small scale yielding hypothesis is assumed, and an a posteriori model reduction of the plastic behavior of the cracked structure is performed. This reduced model provides incrementally the plastic state in the vicinity of the crack front, from which the instantaneous crack growth rate is inferred. Then an additional a priori model reduction technique is used to accelerate even more the time to solution of the whole problem. This a priori approach consists in building incrementally and without any previous calculations a reduced basis specific to the considered test-case, by extracting information from the evolving displacement field of the structure. Then the displacement solutions of the updated crack geometries are sought as linear combinations of those few basis vectors. The numerical method chosen for this work is the finite element method. Hence, during the propagation the spatial discretization of the model has to be updated to be consistent with the evolving crack front. For this purpose, a specific mesh morphing technique is used, that enables to discretize the evolving model geometry with meshes of the same topology. This morphing method appears to be a key component of the model reduction strategy. Finally, the whole strategy introduced above is embedded inside an adaptive approach, in order to ensure the quality of the results with respect to a given accuracy. The accuracy and the efficiency of this global strategy have been shown through several examples; either in bidimensional and tridimensional cases for model crack propagation, including the industrial example of a helicopter structure.

La thèse présente une méthodologie pour construire des éléments finis de type « solide-coque » avec intégration réduite en vue des applications à la simulation de la mise en forme tel que l’emboutissage des tôles où ces éléments finis doivent présenter de bonnes aptitudes à modéliser la flexion mais également les situations de laminage de la tôle. A partir des éléments volumiques à 8 nœuds et 3 degrés de liberté par nœud (les 3 composantes du déplacement), un neuvième nœud est rajouté au centre de l’élément. Ce neuvième nœud n’est pourvu que d’un seul degré de liberté, le déplacement le long de la direction de l’épaisseur. Cette direction privilégiée a un nombre de points d’intégration supérieur ou égal à 3 mais l’intégration est réduite au centre de l’élément diminuant très sensiblement les temps CPU par rapport à une intégration complète. Un soin particulier a été pris pour contrôler tous les modes à énergie nulle dus à l’intégration réduite. Ce nœud supplémentaire permet une distribution linéaire de la déformation normale. Avec les lois de comportement complètement 3D ces nouveaux éléments solide-coque donnent des résultats similaires en flexion à ceux obtenus avec des éléments coques et état plan de contrainte. Le neuvième nœud joue le rôle d’un paramètre supplémentaire pour l’interpolation quadratique du déplacement dans la direction de l’épaisseur. Ce degré de liberté a une signification physique et un effort équivalent à une pression normale peut être prescrit. Dans les situations de pression normale et dans le cas du contact, la contrainte normale obtenue est physique ce qui n’est pas le cas de nombreux éléments solide-coque de la littérature. Le pincement ou le laminage des tôles est correctement modélisé. Pour valider ces éléments, un module d’emboutissage en U avec passage et laminage de la bande de tôle sur des rouleaux a été construit au laboratoire. La comparaison entre les efforts d’emboutissage calculés et mesurés est très bonne ainsi que la géométrie des bandes de tôle obtenue après retour élastique. / This thesis presents a methodology for developing under-integrated “solid-shell” finite elements for sheet forming simulations like deep drawing where these elements must offer a bending capability and sheet thinning conditions as well. Starting from 8-node elements endowing three degrees of freedom per node (three displacement components), a ninth node is added at the centre of the element. This extra node has just one degree of freedom: a displacement along the ‘thickness’ direction. Several integration points are distributed along this privileged direction (5 points, generally) but the in-plane reduced integration at the centre of the element decreases CPU costs compared to a full integration. A special care has been taken to control all zero-energy modes due to the reduced integration. This additional node allows a linear distribution of the normal strain. With fully-3D constitutive laws, these new solid-shell elements give similar bending results as those obtained with shell elements and a plane stress state hypothesis. This ninth node acts as an additional parameter for the quadratic interpolation of the displacement in the ‘thickness’ direction. The corresponding degree of freedom has a physical meaning and a force, equivalent to a normal pressure for instance, may be prescribed. In situations of a normal pressure and in the case of contact, the obtained normal stress is physically defined, which is not the case for many solid-shell elements found in the literature. The pinching (or the thinning) of sheets is properly modelled. To validate these elements, an apparatus for U-drawing tests with ironing or thinning on strip sheets has been built in the laboratory. The comparison between numerical and experimental punch force during sheet forming is pretty good as well as the geometry of blank after springback.

L’instabilité des cavités souterraines (mines, carrières, tunnels,…) peut induire les mouvements de terrains d’amplitude suffisante pour endommager les bâtiments et les infrastructures en surface. Les méthodes traditionnelles, utilisées dans les pratiques d’ingénieur pour prévoir les déformations dans les structures, sont basées sur les caractéristiques des mouvements de terrain en condition de terrain vierge sans prendre en compte l’effet de la présence des structures en surface. L’objectif de cette thèse est de prédire les déformations des ouvrages en tenant compte de l’influence de l’interaction sol-structure, d’une part ; et d’évaluer la performance d’une solution de protection (tranchée périphérique), d’autre part. Cela a été achevé par la réalisation d’études paramétriques utilisant deux approches complémentaires : une approche expérimentale à l’aide d’un modèle réduit physique 3D sous gravité normale et une modélisation numérique 3D par la méthode des éléments finis. En particulier l’effet d’un certain nombre de paramètres géométriques et mécaniques a pu être investigué dans l’étude de l’interaction sol-structure : la position de la structure par rapport à la cuvette d’affaissement, le poids de la structure et la raideur relative entre le sol et la structure. Concernant l’étude de l’efficacité de tranchées périphériques, l’effet de la position de la structure, de la position de la tranchée vis-à-vis de la structure et de la rigidité de la tranchée a été analysé. Les résultats obtenus ont abouti à une meilleure compréhension du problème d’interaction sol-structure et ont montré l’importance de cet effet qui doit être pris en compte dans l’évaluation de la vulnérabilité du bâti. Le transfert des mouvements du sol à la structure est faible (moins de 2,5%), dans le cas modélisé : structure rigide et interface glissante. Les différents résultats ont permis par ailleurs de mettre en évidence l’efficacité de la tranchée périphérique pour réduire les sollicitations affectant les structures. La tranchée doit être remplie avec un matériau très déformable et surtout placée à une distance de l’ordre d’un mètre de la structure. / The instability of underground cavities due to mining or tunneling activities can induce surface ground movements which damage overlying buildings and infrastructures. The conventional design methods, used in engineering practices to assess the structure’s deformations, are based on greendfield ground movements without taking into account the effect of the surface structure. The objective of this thesis is to predict the building deformations by taking the in-fluence of the soil-structure interaction into consideration, and to evaluate the performance of a mitigation technique (peripheral trench). This has been achieved by performing parametric studies using both the experimental approach by means of a 3D small-scale physical model under earth gravity condition and the 3D finite element numerical modeling. In particular the effect of building position, building weight, and relative stiffness of building and underlying soil has been investigated in the soil-structure interaction study. Concerning the trench effectiveness analysis, the effect of building position, trench position relatively to the building, and the trench stiffness has been investigated. The results of these investigations have led to a better understanding of the soil-structure interaction problem and have pointed out the importance of this effect which should be taken into account for building damage assessment. The transfer of ground movements to the buildings is low (less than 2.5%) for the modeled case: stiff structure and sliding contact. In addition, the results have proved the effectiveness of the peripheral trench to reduce the ground movements affecting the buildings. The trench should be filled with a very deformable material and located at a distance of around one meter from the building.

Les bassins de retenue des eaux pluviales sont utilisés pour préserver la qualité des eaux réceptrices par sédimentation pendant le temps de pluie. Cependant, les efficacités du bassins n'étaient pas satisfaisants en raison de la mal compréhension du processus de sédimentation. Afin de mieux comprendre ces processus dans des ouvrages in situ, cette thèse porte à la fois sur des expérimentations in situ et sur les modélisations de l'hydrodynamique et du transport particulaire dans les bassins de retenue pilotes et in situ. Cette recherche s’est appuyée en grande partie sur le bassin Django Reinhardt (BDR) à Chassieu (volume: 32000 m3, surface: 11000 m2) dans le cadre de l’OTHU et sur les données expérimentales obtenues par Dufresne (2008) et Vosswinkel et al. (2012). Les échantillons de sédiments ont été prélevés et leurs caractéristiques physiques ont été analysées en laboratoire dans le but de cerner leur distribution spatiale. Concernant la modélisation numérique, les simulations de l’hydrodynamique en régime permanent ont été réalisées à l'aide du logiciel CFD Fluent et ont été évaluées à partir de l’analyse de corrélation entre le comportement hydrodynamique du bassin et la distribution spatiale des caractéristiques physiques des sédiments. Les conditions limites sur le fond couramment utilisées et largement décrites dans la littérature ont été testées dans le but de représenter la distribution spatiale des sédiments et l’efficacité de décantation du BDR. Les conditions testées sont : i) contrainte de cisaillement critique ou bed shear stress – BSS et ii) énergie cinétique turbulente critique ou bed turbulent kinetic energy - BTKE. L’approche Euler-Lagrange dite « particle tracking » a été mise en œuvre. En raison de l'échec de prédiction des zones de dépôt à l’aide des conditions limites disponibles (BSS et BTKE), une nouvelle relation a été proposée pour estimer le seuil BTKE. La condition à la limite obtenue en utilisant cette nouvelle relation a été testée sur un bassin pilote (Dufresne, 2008) et sur le BDR en régime permanent. Les résultats obtenus n’étaient pas très satisfaisants concernant la prédiction des zones de dépôt et l’efficacité de décantation dans le bassin BDR, même en considérant une distribution granulométrique non uniforme. Afin de mieux prédire les zones de dépôt dans le BDR, une nouvelle méthode a été proposée en considérant le transport des particules, leur décantation et leur érosion en régime transitoire. Sur la base de la méthode proposée pour le transport des particules, la décantation et l'érosion en régime transitoire, plusieurs modélisations avec différentes conditions limites ont été réalisées dans un bassin de retenue pilote rectangulaire (Vosswinkel et al., 2012). Les prédictions des efficacités et des zones de dépôt en régime transitoire avec la méthode proposée sont satisfaisantes / Stormwater detention basins are used to preserve the quality of receiving waters by sedimentation during the wet weather. However, the removal efficiencies of basin were not satisfactory due to the not well understanding of the sedimentation processes. In order to further understand these processes in the real facilities, this thesis therefore focuses both on in situ experiments and modeling of the hydrodynamic and sediment transport in field detention basin and in small scale basin in laboratory. This research was supported by large part on the Django Reinhardt basin (DRB) in Chassieu within the OTHU program and the experimental data deriving from Dufresne (2008) and Vosswinkel et al. (2012). Samples of sediments accumulated in the basin were collected and their physical characteristics were analyzed in order to determine their spatial distribution. Concerning numerical modeling, the hydrodynamic simulations in steady state were performed using CFD software Fluent and were evaluated by the correlation analysis between the hydrodynamic behavior of DRB and the spatial distribution of the physical characteristics of sediments. The bed boundary conditions used in the literatures were tested in order to represent the spatial distribution of sediments and removal efficiency of DRB. The conditions tested were: i) critical bed shear stress - BSS and ii) critical bed turbulent kinetic energy - BTKE. Because of the failure prediction of DRB deposit zones with usual bed boundary conditions, a new relationship based on particle settling velocities has been proposed to estimate the BTKE threshold for the bed boundary condition. The proposed boundary condition was tested in a pilot basin (Dufresne, 2008) and the DRB using the Euler-Lagrange approach under steady flow conditions. The results were not very satisfactory regarding the DRB deposit zones, even considering non-uniform grain size. In order to better predict the deposit zones and settling efficiency in field detention basins, a new method has been proposed accounting for the sediment transport, settling and erosion under unsteady conditions. Based on this proposed method for representing the particle transport, settling and erosion processes under unsteady conditions, various simulations with different bed boundary conditions were carried out in a pilot rectangular basin (Vosswinkel et al., 2012). The predictions of removal efficiencies and deposition zones are satisfactory. Hence, taking into account transient effects on both hydrodynamics and sediment transport leads to drastically improve the spatial and temporal distributions of sediments in settling detention basins.

Cette thèse s’inscrit dans le contexte du développement de Bâtiments Basse Consommation. La conception de telles constructions les rend sensibles aux sollicitations internes. Aussi, les outils de thermique du bâtiment existants ne sont pas adaptés pour simuler assez fidèlement ce type de bâtiments, si bien qu’un modèle tridimensionnel et dynamique a été développé ici. Celui-ci présente plusieurs particularités : il s’appuie sur une discrétisation spatiale optimisée des parois, la tache solaire y est localisée et l’intégration des dynamiques des conditions environnementales est assurée par un solveur numérique à pas de temps adaptatif et un seul nœud d’air est considéré. La validation du modèle s’est suivant une confrontation avec des mesures en conditions réelles réalisées dans une cellule de BESTlab d’EDF R&D. Un suivi visuel de la tache solaire a permis de confirmer sa bonne localisation par notre modèle. Des mesures de température en surface complétées par des cartographies thermographiques ont été comparées aux champs de températures simulés, montrant une bonne concordance. Les comparaisons de températures d’air mesurées et simulées ont montré des résidus ne dépassant pas 1,5 ˚C, pour des erreurs moyennes de 0,5 ˚C. La pertinence des deux principales innovations du modèle a été ensuite démontrée : l’utilisation d’entrées échantillonnées à la minute associées à un solveur à pas de temps adaptatif permet de minimiser les erreurs de simulation : en mi-saison, les résidus maximaux sont respectivement de 1 ˚C et 2 ˚C pour des entrées à la minute et à l’heure. En hiver, les températures d’air simulées tendent à plus osciller autour de la consigne quand le pas d’échantillonnage des entrées s’allonge. Deux modèles unidimensionnels, représentatifs de modèles courants, M1D,sol diluant le rayonnement solaire sur le sol seul et M1D,parois le distribuant de façon homogène sur les parois au prorata de la taille de la tache solaire censée les frappées, ne dégradent que légèrement la précision des calculs de température d’air. Cependant, ces modèles 1D ne permettent pas de calcul des champs de températures sur les parois si bien qu’ils présentent des erreurs locales dépassant 20 ˚C aux endroits touchés par la tache solaire. Enfin en hiver, le modèle 3D permet de prédire des consommations de chauffage surestimées de 6,5 % quand M 1D,parois les surestime de 11 % et M1D,sol de 22 %. Les améliorations apportées par notre modèle ont été confirmées pour d’autres types de cellules. D’ailleurs des écarts plus importants entre M1D,sol et le modèle 3D ont été observés pour une cellule dont parois et sol ont des compositions très différentes, alors que l’orientation a aussi un impact. Ce travail confirme la nécessité de représenter plus finement les phénomènes physiques pour des locaux fortement isolés. Des améliorations sont à intégrer, comme la description de l’anisothermie de l’air. / Low energy building constructions become sensitive to internal gains : any internal heating source has an impact on the envelope. Therefore, it is important to evaluate the performance of current transient thermal models when adapted to low energy buildings. This work describes a numerical model to simulate a single room, using a refined spatial three-dimensional description of heat conduction in the envelope but a single air node is considered. The model has been developed for environmental conditions that vary over short time-steps and has integrated the projection of solar radiation through a window onto interior walls : the sun patch. The validation of the model has been done through a detailed comparison between model and measurements. The in situ experiment has been carried out in one of the BESTlab cells (EDF R&D). The sun patch has been followed by a camera to validate its calculated position and surface. Temperature measurements by thermocouples and by thermal cameras have been compared to the models outputs. Differences between air and surface temperatures measured and simulated were never above 1.5 ˚C and mean errors reached 0.5 ˚C. The two innovations of the model have then be proven. Using minute wise weather data and inputs associated to an adaptative solver, enabled to pull down simulation errors : in May maximal differences rised from 1 ˚C to 2 ˚C for respectivelly one minute and hourly wise inputs. More important errors are seen in summer whereas in winter, air temperatures simulated tend to more fluctuate around the set up temperature when the sampling step gets longer. Two one dimensional models, close to traditional taken simulation tools, were used. Model M 1D,sol supposed the incoming radiation to reach only the floor. A 1D model with sun patch movement, called 1D,parois , was also used. These two models evaluated the air temperature with an acceptable error. However, their surface temperatures were still subject to important errors. Thus, for temperature surfaces evaluation, both 1D model presented differences up to 20 ˚C for surfaces touched by the sun patch. In winter, the 3D model can predict heating energy consumptions overestimated by 6.5 % when M 1D,parois overestimated them by 11 % and M1D,sol by 22 %. The improvements brought by our model have been proven also for other cells with different thermal masses. For these cells, differences between M1D,sol and the 3D model could reach 4.5 ˚C. Differences seemed to be more important for low thermal mass cells, and the orientation of the building had a strong impact. This work has confirmed the necessity of representing more accuratelly the descriptions of the enveloppe for strongly insulated rooms. To improve the model, the anisothermal hypotheses of the air should be considered.