No description available.

Le bois est utilisé comme matériau absorbeur de chocs dans les conteneurs de transport du CEA. La simulation numérique des crashs tests de ces structures fait partie intégrante de leur développement et de leur validation. Il est donc impératif d’étudier expérimentalement les différentes essences de ces matériaux afin de mettre en évidence leurs caractéristiques et de déterminer leurs lois de comportement. En effet, ces données sont absolument nécessaires afin d’alimenter les modèles numériques retenus pour les simulations. A partir des modèles existants dans les codes de calcul ou par la programmation de lois de comportement spécifiques prenant en compte le paramètre vitesse de déformation, le comportement de ce matériau cellulaire a été simulé afin de participer à ce programme de validation produit. / Wood is usually used as material absorber of shocks in CEA’s containers. Numerical simulation of crash-test of these structures is an integral part of their development and their validation. It is thus imperative to experimentally study various type of wood to highlight their characteristics and determine their behavior. Indeed, these data are absolutely necessary to feed the numerical models chosen for the simulations. From the existing models in simulation tools or by programming specific laws, both taking into account the parameter strain rate, the behavior of this cellular material was modeled to participate in this product program.

Cette étude présente le développement d'une méthode expérimentale et numérique de caractérisation du comportement d'un matériau composite à renforts tissés, depuis l'échelle microscopique jusqu'au niveau macroscopique, en augmentant progressivement la complexité des géométries et des cas de chargement. Un nouveau modèle matériau numérique a été développé sous la forme d'un UMAT pour le logiciel éléments finis LS-DYNA. Ce modèle est capable de prédire un comportement linéaire ou non-linéaire, la déformation plastique, la rupture et l'endommagement progressif du matériau après rupture, tout en respectant les différentes restrictions imposées par le domaine industriel. Une formulation non-locale a été implémentée afin d'éviter les phénomènes indésirables de localisation de l'endommagement et d'offrir l'indépendance du modèle à la finesse du maillage. La première étape du processus de validation est la simulation des essais quasi-statiques sur des échantillons, ce qui permet de définir les paramètres numériques fondamentaux correspondant aux propriétés mécaniques et au comportement macroscopique du matériau étudié. Ensuite, les simulations d'impact assurent l'étude des modèles de dégradation. Une bonne corrélation a été obtenue pour les différents cas de chargement. / This study presents the development of an experimental and numerical material characterization method relying on the investigation of the material behavior of woven composites, from the microscopic scale up to the macroscopic level by increasing progressively the complexity of the geometries and load cases. Taking into account the different limitations imposed by the industrial and engineering fields, a new numerical mesh-insensitive material model has been developed as a UMAT for the explicit finite elements solver LS-DYNA. This model accounts for nonlinearity, plasticity, failure and post-failure damage in order to predict the macroscopic behavior and the damage tolerance of woven laminated composites. It includes also a smeared formulation to avoid undesirable localization phenomena. The first step of the validation process relies on the simulation of quasi-static coupon tests, which allows to set the fundamental numerical parameters corresponding to the effective material properties and to the macroscopic behavior observed experimentally. Then, the post-failure material degradation models have been investigated using the simulations of impact tests on elementary plates and on representative parts. A satisfying agreement has been obtained for the different load cases.

La demande en termes de miniaturisation des objets manufacturés est sans cesse croissante. La fabrication en grande série des pièces submillimétriques est rendue difficile par un certain nombre de phénomènes liés à l'échelle. L'objectif de cette thèse est de développer une approche complète permettant de prédire la réponse mécanique d'une tôle d'épaisseur submillimétrique lors d'une opération de mise en forme. Dans un premier temps, des essais expérimentaux en traction uniaxiale sur des tôles minces de cuivre, d'aluminium et d'acier inoxydable sont réalisés. Ils permettent, entre autres, de mettre en évidence l'influence de la surface libre et de la distribution de taille de grains sur la réponse mécanique des tôles. Dans un deuxième temps, une loi de comportement élasto-viscoplastique monocristalline, formulée en grandes déformations et prenant en compte les mécanismes physiques de déformation plastique, est développée et implantée dans le code éléments fins ABAQUS®. Cette loi de comportement utilise les densités de dislocations comme variables internes. Pour être fonctionnelle, cette approche nécessite la description de la microstructure du matériau simulé. Un traitement statistique de mesures réelles est développé afin de générer un modèle éléments finis représentant de façon précise la morphologie des grains, leurs tailles et leurs orientations cristallographiques. Dans un troisième temps, la mise en oeuvre de l'approche complète développée permet de reproduire les phénomènes que nous observons expérimentalement, à savoir : · un adoucissement du comportement global du matériau, · l'apparition d'un état de surface appelé peau d'orange, · une forte localisation de la déformation dans l'épaisseur de l'éprouvette, orientée perpendiculairement à la direction de traction. De plus, grâce à ce nouvel outil numérique, nous vérifions l'influence de la surface libre et la prédominance des gros grains sur la réponse globale des structures calculées. Finalement, une machine d'essais in situ est conçue et réalisée dans le but de valider l'approche numérique, dans des conditions proches de celles observées en micro-formage.

Dans cette thèse, deux méthodes de couplage sont proposées pour la simulation multi-échelle des procédés de mise en forme. Dans la première partie, une procédure simplifiée (couplage indirect) est adoptée pour coupler les codes éléments finis (Abaqus et LAM3) au modèle polycristallin avec un schéma de transition autocohérente basée sur le comportement élastoplastique du monocristal écrit dans le formalisme des grandes déformations. Cette procédure simplifiée consiste à lier le modèle polycristallin avec l'analyse EF par l'extraction de l'histoire de l'incrément de déformation et de contrainte macroscopique, obtenue à partir d'une simulation EF préliminaire avec une loi phénoménologique, et à l'utiliser comme trajet de chargement dans le modèle polycristallin. Cette méthode est appliquée pour la simulation multi-échelle du procédé de skin-pass. Le suivi du trajet de chargement extrait dans la demi-épaisseur de la tôle a permis de prédire l'évolution des grandeurs physiques associées au modèle de plasticité en particulier la texture cristallographique, la texture morphologique et l'écrouissage. Dans la seconde partie de cette thèse, un modèle polycristallin élastoplastique du type autocohérent en petites déformations est couplé au code EF Abaqus via la routine utilisateur UMAT. Ce couplage (dit couplage direct) consiste à utiliser la théorie de la plasticité cristalline comme loi de comportement à chaque point d'intégration du maillage EF. Le polycristal est représenté par un ensemble de N monocristaux. Chaque fois que le code EF a besoin d'information sur le comportement mécanique aux points d'intégration de chaque EF, le modèle polycristallin est appelé. Pour valider ce couplage développé, nous avons effectué des cas tests de simulation de trajets rhéologiques. Les résultats issus de ce couplage ont été validés avec des modèles de référence. À la différence des modèles phénoménologiques, ce couplage permet non seulement d'avoir des informations sur le comportement macroscopique de la structure mais aussi d'obtenir des informations sur l'état de la microstructure du matériau. / In this thesis, two coupling methods are proposed for the multiscale simulation of forming processes. In the first part, a simplified procedure (indirect coupling) is adopted to couple the finite element codes (Abaqus and LAM3) with a polycrystalline selfconsistent model based on the large strain elastoplastic behavior of single crystals. This simplified procedure consists in linking the polycrystalline model with the FE analysis by extracting the history of the increment of macroscopic strain and stress, obtained from a preliminary FE simulation with a phenomenological law, and then using it as loading path prescribed to the polycrystalline model. This method is applied to multiscale simulation of skin-pass processes. By following on the loading path extracted at the halfthickness of the sheet, we can predict the evolution of some physical parameters associated with the plasticity model, in particular the crystallographic texture, the morphological texture and hardening. In the second part on this thesis, a small strain version of the elastoplastic polycristalline self-consistent model is coupled to the Abaqus FE code via the user material subroutine UMAT. This coupling (called direct coupling) consists in using crystal plasticity theory as constitutive law at each integration point of the FE mesh. The polycristal is represented by a set of N single crystals. Each time the FE code needs information on the mechanical behavior at the integration points considered, the full polycrystalline constitutive model is called. In order to validate this coupling, simulations of simple mechanical tests have been conducted. The results of this coupling have been validated through comparison with reference models. Unlike phenomenological models, this coupling provides not only information on the overall macroscopic response of the structure, but also important information related to its microstructure

La technologie flip chip est de plus en plus répandue dans l'industrie électronique [trois dimensions (3D) System in Package] et est principalement utilisée pour la fabrication de réseaux détecteurs de grand format (mégapixels) et faible pas. Pour étudier la fiabilité de ces assemblages, des simulations numériques basées sur des méthodes d'éléments finis semblent être l'approche la moins chère. Cependant, de très grands assemblages contiennent plus d'un million de billes de brasure, et le processus d'optimisation de ces structures par des simulations numériques se révèle être une tâche très fastidieuse. Dans de nombreuses applications, la couche d'interconnexion de tels assemblages flip chip se compose de microbilles de brasure noyées dans de l'époxy. Pour ces configurations, nous proposons une approche alternative, qui consiste à remplacer cette couche d'interconnexion hétérogène par un matériau homogène équivalent (MHE). Un modèle micromécanique pour l'estimation de ses propriétés thermoélastiques équivalentes a été mis au point. La loi de comportement obtenue pour le MHE a ensuite été implémentée dans le logiciel par éléments finis (Abaqus®). Les propriétés élastiques des matériaux de l'assemblage sont définies par la littérature et également déterminées expérimentalement par une méthode de caractérisation mécanique : la nano-indentation. Les réponses thermomécaniques des assemblages testés soumis à des chargements correspondant aux conditions de fabrication ont été analysées. La technique d'homogénéisation-localisation a permis d'estimer les valeurs moyennes des contraintes et des déformations dans chaque phase de la couche d'interconnexion. Pour accéder plus précisément aux champs de contraintes et déformations dans ces phases, deux modèles de zoom structurel (couplage de modèles et submodeling), en tenant compte de la géométrie réelle de la bille de brasure, ont été testés. Les champs de contrainte et de déformation locaux obtenus corroborent avec les initiations de dommage observées expérimentalement sur les billes de brasure / Flip chip technology is increasingly prevalent in electronics assembly [threedimensional (3D) system in package] and is mainly used at fine pitch for manufacture of megapixel large focal-plane detector arrays. To estimate the reliability of these assemblies, numerical simulations based on finite-element methods appear to be the cheapest approach. However, very large assemblies contain more than one million solder bumps, and the optimization process of such structures through numerical simulations turns out to be a very time-consuming task. In many applications, the interconnection layer of such flip-chip assemblies consists of solder bumps embedded in epoxy filler. For such configurations, we propose an alternative approach, which consists in replacing this heterogeneous interconnection layer by a homogeneous equivalent material (HEM). A micromechanical model for the estimation of its equivalent thermoelastic properties has been developed. The constitutive law of the HEM obtained was then implemented in finite-element software (Abaqus®). Elastic properties of materials that compose the assembly were found in literature and by using mechanical characterization method especially nano-indentation. Thermomechanical responses of tested assemblies submitted to loads corresponding to manufacturing conditions have been analyzed. The homogenization-localization process allowed estimation of the mean values of stresses and strains in each phase of the interconnection layer. To access more precisely to the stress and strain fields in these phases, two models of structural zoom (model coupling and submodeling), taking into account the real solder bump geometry, have been tested. The local stress and strain fields obtained corroborate the experimentally damage initiation of the solder bumps observed

Despite the phenomena of fracture and fatigue having been the focus of academic research for more than 150 years, it remains in effect an empirical science lacking a complete and comprehensive set of predictive solutions. In this regard, the focus of the research in this thesis is on the development of new cohesive-zone models for fracture and fatigue that are afforded an ability to capture strain-rate effects. For the case of monotonic fracture in ductile material, different combinations of material response are examined with rate effects appearing either in the bulk material or localised to the cohesive-zone or in both. The development of a new rate-dependent CZM required first an analysis of two existing methods for incorporating rate dependency, i.e.either via a temporal critical stress or a temporal critical separation. The analysis revealed unrealistic crack behaviour at high loading rates. The new rate-dependent cohesive model introduced in the thesis couples the temporal responses of critical stress and critical separation and is shown to provide a stable and realistic solution to dynamic fracture. For the case of fatigue, a new frequency-dependent cohesive-zone model (FDCZM) has been developed for the simulation of both high and low-cycle fatigue-crack growth in elasto-plastic material. The developed model provides an alternative approach that delivers the accuracy of the loading-unloading hysteresis damage model along with the computational efficiency of the equally well-established envelope load-damage model by incorporating a fast-track feature. With the fast-track procedure, a particular damage state for one loading cycle is 'frozen in' over a predefined number of cycles. Stress and strain states are subsequently updated followed by an update on the damage state in the representative loading cycle which again is 'frozen in' and applied over the same number of cycles. The process is repeated up to failure. The technique is shown to be highly efficient in terms of time and cost and is particularly effective when a large number of frozen cycles can be applied without significant loss of accuracy. To demonstrate the practical worth of the approach, the effect that the frequency has on fatigue crack growth in austenitic stainless-steel 304 is analysed. It is found that the crack growth rate (da/dN) decreases with increasing frequency up to a frequency of 5 Hz after which it levels off. The behaviour, which can be linked to martensitic phase transformation, is shown to be accurately captured by the new FDCZM.

This study provides an overview of the implementation of a nonlinear microstructural constitutive model in ABAQUS employing a user subroutine at the level of the biomedical engineer. Two different element formulations are employed: a continuum incompressible and a plane stress incompressible. All examples are validated by performing a number of deformations on 2D and 3D square elements and comparing the analytical formulation in a programming language and the user subroutine in ABAQUS. Application models will be presented that provide a deeper look into the impacts of soft tissue deformation, damage, and fracture. Additionally, we investigate the mechanical behavior of skin layers in terms of the nominal stress-strain curve using uniaxial and cyclic loading tests on porcine skin specimens in two forms: dermis integrating epidermis and hypodermis. Experiments were performed on specimens from the belly and breast of the pigs and under both orthogonal orientations with respect to the spine direction. All tests were carried out at room temperature with cyclic loading at a constant strain rate and increasing stretch increments. Finally, data is fitted using microstructural constitutive model.

Every moment of every day, at least hundreds of thousands of tires roll across a surface throughout the world. Tires are indisputably important in our daily life. The tire's primary component is rubber, which consumes energy when it rotates on a substrate due to the viscoelastic material's internal friction: a phenomenon referred to as rolling resistance. The interaction between the tire and the road surface is one of the most intricate and crucial phenomena in an automobile, because it is responsible for creating forces, moments, and deformation in the tire. Additionally, the road's roughness interacts with the tire and contributes significantly to its performance. This dissertation aims to develop a comprehensive physics-based model for predicting the rolling resistance of a viscoelastic material due to dynamic deformations caused by tire rotation using an analytical approach. The model was developed by proposing a Gaussian wave function propagating across a tire circumference's viscoelastic medium. The wave function was selected to describe the displacement field produced by tire-road interaction. Additionally, by adopting a multi-scale modeling technique, the model was upgraded to estimate rolling resistance while taking into account surface roughness at all length scales, from macroscopic to microscopic. Additionally, another mathematical model was developed using the Fourier series approach to evaluate the steady-state stress response and energy dissipation for any harmonic and non-harmonic periodic strain signals. Additionally, the dissertation strove to build a continuum damage mathematical model using a combined testing/modeling methodology to predict the aging of Styrene-Butadiene Rubber (SBR) after continuous exposure to the atmosphere. The obtained model was developed through the implementation of optimization techniques while formulating a mathematical model, which was then combined with a physics-based model to predict rolling resistance while taking into account rubber aging. Calibration of hyperelastic and viscoelastic material models with testing data was performed using an optimization technique that yielded sufficient results. The results of all mathematical models obtained in this dissertation are reported subsequently. The stress response of a viscoelastic material under harmonic and non-harmonic strain input yielded good agreement with the FEA model obtained using ABAQUS. The rolling resistance behavior under various operating conditions, including texture and aging effects, was reported, and the results aligned with the experimental results found in the literature. / Doctor of Philosophy / Every moment of every day, hundreds of thousands of automobile tires roll across a surface somewhere in the world. A tire is an undeniably important part of everyday life. Rubber is the tire's main component, and when it rotates on a surface, it loses energy, resulting in a force that resists motion, known as rolling resistance force. The contact between the tire and the road is one of the most complicated and important phenomena that happens in an automobile because it is responsible for the vehicle's dynamic performance in areas such as acceleration, stopping distance, and stability. Another factor that affects tire and car performance and should be taken into account is the road's roughness. This dissertation used an analytical method to come up with an accurate physics-based model for predicting the rolling resistance force of a viscoelastic material caused by tire rotation. The model was developed by assuming a Gaussian wave function would move across the tire circumference. Additionally, using a multi-scale modeling technique, the model was improved so that it could calculate the value of rolling resistance force considering surface roughness in all lengths of scale. This project also developed an additional mathematical model using the Fourier series method to determine how the stress response and energy dissipation would behave for any harmonic and nonharmonic periodic strain signals. Additionally, the dissertation presents the developing of a continuum damage mathematical model that could predict the material property of styrene-butadiene rubber (SBR) after being exposed to the air for a long time (i.e., aged). The model was developed based on experimental data and optimization techniques. This model was then combined with a physics-based model to predict rolling resistance force while taking aging into account. The material models were defined using an optimization method that yielded good results. The stress response of a viscoelastic material when it was subjected to harmonic and non-harmonic strain was in good agreement with the Finite Element Analysis (FEA) model made with ABAQUS. Rolling resistance behavior was observed, and the results were consistent with those found in the literature.

Indonesia is a developing country that suffers from earthquakes and windstorms and where at least 60% of houses are non-engineered structures, built by unskilled workers using masonry and timber. The non-engineered housing units developed in urban region are also vulnerable to seismic hazard due to the use of low quality of material and constructions method. Those structures are not resistant to extreme lateral loads or ground movement and their failure during an earthquake or storm can lead to significant loss of life. This thesis is concerned with the structural performance of Indonesian low-rise buildings made of masonry and timber under lateral seismic load. The research presented includes a survey of forms of building structure and experimental, analytical and numerical work to predict the behaviour of masonry wall and traditional timber frame buildings. Experimental testing of both masonry and timber have been carried out in Indonesia to establish the quality of materials and to provide material properties for numerical simulations. The experimental study found that the strength of Indonesia-Bali clay brick masonry are below the minimum standard required for masonry structures built in seismic regions, being at least 50% lower than the requirement specified in British Standard and Eurocode-6 (BS EN 1996-1-1:2005). In contrast, Indonesian timber materials meet the strength classes specified in British Standard/Eurocode- 5 (BS EN 338:2009) in the range of strength grade D35-40 and C35).Structural tests under monotonic and cyclic loading have been conducted on building components in Indonesia, to determine the load-displacement capacity of local hand-made masonry wall panels and timber frames in order to: (1) evaluate the performance of masonry and timber frame structure, (2) investigate the dynamic behaviour of both structures, (3) observe the effect of in-plane stiffness and ductility level, and (4) examine the anchoring joint at the base of timber frame that resists the overturning moment. From these tests, the structural ductility was found to be less than two which is below the requirement of the relevant guidelines from the Federal Emergency Management Agency, USA (FEMA-306). It was also observed that the lateral stiffness of masonry wall is much higher than the equivalent timber frame of the same height and length. The experimental value of stiffness of the masonry wall panel was found to be one-twelfth of the recommended values given in FEMA-356 and the Canadian Building code. The masonry wall provides relatively low displacement compared to the large displacement of the timber frame at the full capacity level of lateral load, with structural framing members of the latter remaining intact. The weak point of the timber frame is the mechanical joint and the capacity of slip joint governs the lateral load capacity of the whole frame. Detailed numerical models of the experimental specimens were setup in Abaqus using three-dimensional solid elements. Cohesive elements were used to simulate the mortar behaviour, exhibiting cracking and the associated physical separation of the elements. Appropriate contact definitions were used where relevant, especially for the timber frame joints. A range of available material plasticity models were reviewed: Drucker-Prager, Crystalline Plasticity, and Cohesive Damage model. It was found that the combination of Crystalline Plasticity model for the brick unit and timber, and the Cohesive Damage model for the mortar is capable of simulating the experimental load-displacement behaviour fairly accurately. The validated numerical models have been used to (1) predict the lateral load capacity, (2) determine the cracking load and patterns, (3) carry out a detailed parametric study by changing the geometric and material properties different to the experimental specimens. The numerical models were used to assess different strengthening measures such as using bamboo as reinforcement in the masonry walls for a complete single storey, and a two-storey houses including openings for doors and windows. The traditional footing of the timber structures was analysed using Abaqus and was found to be an excellent base isolation system which partly explains the survival of those structures in the past earthquakes. The experimental and numerical results have finally been used to develop a design guideline for new construction as well as recommendations for retrofitting of existing structures for improved performance under seismic lateral load.