Dans le cadre de la conception des structures automobiles et de la modélisation de leur tenue au crash, la prédiction du comportement mécanique et de la rupture des assemblages soudés par point par résistance, s’avère encore délicate. En introduisant une nouvelle génération d’assemblages multi-tôles en aciers multigrades soudés par point, la prédiction de leur tenue mécanique est devenue encore plus complexe. En effet, malgré une utilisation croissante de ces assemblages, l’étude expérimentale de leur comportement mécanique est encore très limitée, tandis que leur modélisation macroscopique n’a pas encore été proposée. Dans ce contexte, un dispositif expérimental approprié de type Arcan est proposé pour caractériser la tenue mécanique et la rupture des assemblages trois tôles soudés par points, en chargements purs et combinés et dans les différentes configurations possibles. Une approche de modélisation macroscopique est également proposée pour simuler le comportement des assemblages trois tôles soudés, pour des applications de type structure. Un modèle d’élément équivalent ’connecteur’ est tout d’abord proposé, identifié et validé par rapport un assemblage deux tôles soudé par point. Ensuite, une approche de modélisation utilisant trois connecteurs est proposée pour simuler le comportement mécanique et la rupture d’un assemblage trois tôles soudé par point par résistance. Une stratégie d’identification est aussi détaillée pour simplifier la procédure d’identification des paramètres du modèle. Le modèle est identifié et validé par rapport aux résultats expérimentaux des essais de caractérisation de type Arcan réalisés au cours de cette thèse. / Resistance spot welding (RSW) is widely used in the automotive industry. However, the prediction of their non-linear mechanical behavior and their failure in structural crashworthiness computations remains a challenge. The introduction of a new generation ofmulti-sheet multi-steel grade spot welded assemblies makes the prediction of their mechanical strength more complex.Despite the increasing use of these assemblies, the experimental study of their mechanical behavior is still very limited, while their macroscopic modeling has not yet been studied. An experimental device based on Arcan principle is thus proposed to characterize the mechanical strength and the failure of three-sheet spot welded assemblies, in pure and combined loads for different loading configurations. A macroscopic modeling approach is also proposed to simulate their mechanical behavior for structural applications. As a first step, a ’connector’ element is proposed, identified and validated for a two-sheet spot welded assembly. Then, a modeling approach using three connectors is proposed to simulate the mechanical behavior and the rupture of a three-sheet multi- steel grade spot welded assembly. An identification methodology is also presented to simplify the identification procedure of the model parameters. The model is identified and validated based on the experimental results of the Arcan tests carried out during this thesis.

This work includes the development of computational fluid dynamic (CFD) and finite element analysis (FEA) models to investigate fluid flow , solidification, and stress in the shell within the mold during continuous casting. The flow and solidification simulation is validated using breakout shell measurements provided by an industrial collaborator. The shell can be obtained by the solidification model and used in a FEA stress model. The stress model was validated by former research related to stress within a solidifying body presented by Koric and Thomas. The work also includes the application of these two models with a transient solidification model and a carbon percentage investigation on both solidification and deformation.

[pt] A importância do sal em diversas atividades humanas tem acompanhado as fases principais da civilização. Aplicações estratégicas e inovadoras envolvendo rochas salinas são planejadas até os dias atuais, tais como barreiras geológicas para abandono de poços e cavernas de estocagem para novas fontes de energia ou descarbonização. Esta tese foca no comportamento geomecânico do sal com atenção especial à fluência, dilatância, efeitos térmicos e seus impactos no comportamento hidráulico. Modelos deram suporte ao desenvolvimento de metodologias para avaliação de condições de integridade e estanqueidade de projetos estratégicos no contexto de energia. Eles incluem o abandono de poços do pré-sal, cavernas de sal e sal como rocha capeadora de reservatórios carbonáticos. Simulações numéricas avançadas utilizando o framework GeMA investigaram cenários representativos de cada projeto. Para tanto, modelos constitutivos de fluência, condições de contorno customizadas e variáveis de saída específicas foram implementados. Dentre os modelos de fluência, o EDMT foi desenvolvido no contexto desta pesquisa de doutorado, com foco no sal brasileiro. Os estudos de caso conduziram a observações importantes. Um poço do pré-sal pode fechar completamente por fluência após a remoção do revestimento. A redução da pressão e o aquecimento do poço aceleraram seu fechamento de modo significativo. No cenário mais crítico, a permeabilidade final foi ligeiramente superior a 100 vezes o valor inicial. Uma caverna de sal pode suportar mecanicamente ciclos de pressão causados por um esquema contingencial de fornecimento de hidrocarbonetos, embora uma permeabilidade final de 6,5 vezes o valor inicial tenha sido observada. Na geomecânica de reservatórios, a contribuição da fluência da rocha capeadora na subsidência é pequena e irreversível mesmo que a pressão do reservatório seja recuperada. As mudanças na permeabilidade foram sutis apesar da contribuição da fluência do sal e da variação de pressão do reservatório. A revisão da literatura, os modelos de fluência, as metodologias de análise, as funcionalidades incorporadas ao framework e as discussões motivadas pelos estudos de caso representam o valor agregado da tese para a comunidade científica e para a indústria no que tange ao aproveitamento estratégico de rochas salinas. / [en] The importance of salt in many human activities has accompanied the main phases of civilization. Strategic and innovative applications in salt rocks are planned until the present day, such as geological barriers for well abandonment and storage caverns for new energy sources or decarbonization. This thesis focuses on the geomechanical behavior of salt with special regards to creep, dilation, thermal effects and their impact on the hydraulic behavior. Models have supported the development of methodologies to assess integrity and tightness conditions of strategic projects in the context of energy. These include Pre-salt well abandonment, salt caverns and salt as caprock of carbonate reservoirs. Advanced numerical simulations using the framework GeMA have investigated representative scenarios of each project. To this end, creep constitutive models, customized boundary conditions and specific output variables have been implemented. Among the creep models, the EDMT has been developed in the context of this doctoral research, focusing on Brazilian salt. The case studies have led to important findings. A Pre-salt well can close completely by creep after casing removal. Pressure reduction and well heating have accelerated closure significantly. In the most critical scenario, the final permeability was slightly higher than 100 times the initial value. A salt cavern can mechanically support pressure cycles caused by a contingency hydrocarbon supply scheme; however, a final permeability of 6.5 times the initial value has been noticed. In reservoir geomechanics, the caprock creep contribution in the subsidence is small and irreversible even if the reservoir pressure is recovered. Permeability changes have been subtle despite the creep contribution and the reservoir pressure variation. The literature review, the creep models, the analysis methodologies, the capabilities incorporated into the framework and the discussions motivated by the case studies represent the added value of the thesis to the scientific community and industry regarding the strategic use of salt rocks.

[pt] O presente trabalho investigou o comportamento mecânico do concreto projetado com fibras aplicado nas escavações da mina Cuiabá e o seu controle tecnológico através do ensaio de duplo puncionamento. Para essa pesquisa utilizouse o traço desenvolvido pelos engenheiros da mina Cuiabá com dois tipos diferentes de fibra sintética e uma fibra de aço. Inicialmente, caracterizou-se o material compósito através de ensaios de flexão em prismas (EN 14651) e de arrancamento. Em seguida, os resultados obtidos foram utilizados em conjunto com os ensaios de flexão em painéis circulares (ASTM C1550), quadrados (EN 14488-5) e de larga escala para o dimensionamento do sistema de suporte utilizado na mina Cuiabá. As frações mais baixas de fibras sintéticas são mais indicadas nos locais onde a qualidade da rocha é razoável ou boa, enquanto as fibras de aço podem ser aplicadas onde a qualidade do maciço rochoso é classificada como muito fraca. Em relação ao controle tecnológico do concreto com fibras, os ensaios de duplo puncionamento realizados em máquinas mais rígidas são mais eficientes em diminuir a extensão da instabilidade pós-fissuração, limitar os erros nas correlações com os ensaios de flexão (EN 14651) e em reduzir a variação das propriedades mecânicas do material. A aplicação do ensaio de duplo puncionamento com controle fechado pela corrente extensométrica foi o método mais efetivo em limitar a instabilidade pós-fissuração independente da máquina de ensaios selecionada. / [en] The present work investigated the mechanical behavior of the fiber reinforced shotcrete applied at the Cuiabá mine excavation and its quality control through the double punch test. For this research, the mix composition was developed by Cuiabá mine Engineers with the addition of two different types of synthetic fibers and one steel fiber. Firstly, the cement based composite was studied at the material level through bending (EN 14651) and pullout tests. Thereafter, the support system of the Cuiabá mine excavation was designed based on round (ASTM C1550), square (EN 14488-5) and large scale panel tests. The application of lower synthetic fiber volume fractions was more suited for the areas with good quality rock indexes, while steel fibers were better suited at field with lower rock classes. Regarding the quality control of FRC, the double punch tests that were conducted in stiffer testing machines were more effective in reducing the post-peak instability, in limiting the errors in the correlations with bending tests (EN 14651) and in reducing the mechanical properties variations of the material. Finally, the application of the double punch test with closed-loop control of the circumferential extensometer turned to be the most effective method in limiting the post-peak instability regardless of the selected machine.

The need for energy efficient, clean and quiet, energy conversion devices for mobile and stationary applications has presented proton exchange membrane (PEM) fuel cells as a potential energy source. The use of PEM fuel cells for automotive and other transient applications, where there are rapid changes in load, presents a need for better understanding of transient behavior. In particular at low humidity operations; one of the factors critical to the performance and durability of fuel cell systems is water transport in various fuel cell layers, including water absorption in membrane. An essential aspect to optimization of transient behavior of fuel cells is a fundamental understanding of response of fuel cell system to dynamic changes in load and operating parameters. This forms the first objective of the dissertation. An insight in to the time scales associated with various transport phenomena will be discussed in detail. In the second component on the study, the effects of membrane properties on the dynamic behavior of the fuel cells are analyzed with focus on membrane dry-out for low humidity operations. The mechanical behavior of the membrane is directly related to the changes in humidity levels in membrane and is explored as a part third objective of the dissertation. Numerical studies addressing this objective will be presented. Finally, porous media undergoing physical deposition (or erosion) are common in many applications, including electrochemical systems such as fuel cells (for example, electrodes, catalyst layer s, etc.) and batteries. The transport properties of these porous media are a function of the deposition and the change in the porous structures with time. A dynamic fractal model is introduced to describe such structures undergoing deposition and, in turn, to evaluate the changes in their physical properties as a function of the deposition. / Ph. D.

Les matériaux architecturés attirent de plus en plus d’attentions de par leur capacité à combiner différentes propriétés ciblées. Dans ce contexte, les matériaux enchevêtrés, et plus particulièrement les matériaux monofilamentaires enchevêtrés, présentent des propriétés intéressantes en terme de légèreté, de ductilité, et de facteur de perte. En raison de l’architecture interne complexe de ces matériaux, leur caractérisation et la compréhension des mécanismes de déformation nécessitent une méthodologie adaptée. Dans cette étude, l’enchevêtrement est réalisé manuellement pour différents fils d’acier et soumis à une compression oedométrique. De manière à étudier le comportement sous charge de ce type de matériaux, un dispositif de compression uniaxiale guidée a été mis en place dans le tomographe. Il est ainsi possible de suivre, à l’aide de mesures quantitatives, la déformation de l’échantillon et l’évolution du nombre de contacts pour différentes fraction volumiques. L’utilisation de ces données microstructurales a permis un meilleure compréhension du comportement mécanique de tels enchevêtrements. Une rigidité pouvant varier de 20 à 200 MPa en fonction des paramètres de mise en forme (diamètre et forme du fil, fraction volumique, matériau constitutif) a été déterminé. Un matériau homogène de rigidité plus faible a pu être obtenu en pré-déformant le fil sous forme de ressort avant enchevêtrement. Le facteur de perte du matériau a ensuite été mesuré à la fois sous chargement statique et dynamique. L’analyse mécanique dynamique a mis en évidence la capacité de ce matériau à absorber de l’énergie avec une valeur de facteur de perte d’environ 0.25. Les propriétés mécaniques du matériau ont tout d’abord été modélisées analytiquement par un modèle de poutres et un bon accord avec les résultats expérimentaux a pu être obtenu en définissant un paramètre d’orientation equivalent, spécifique à la compression oedométrique de matériaux enchevêtrés. En parallèle, un modéle éléments discrets a été developé afin de simuler le comportement en compression de matériaux monofilamentaires enchevêtrés. Ce modèle s’appuie sur une discrétisation du fil en éléments sphériques, acquise à partir de données de tomographie. Bien que seul le comportement élastique du fil constitutif ait été pris en compte, une bonne adéquation entre résultats numériques et expérimentaux a été obtenu en ajustant les coefficients de frottement du modèle. / Playing with the architecture of a material is a clever way of tailoring its properties for multi-functional applications. A lot of research have been made, in the past few years, on what is now referred to as “architectured materials” (metal foams, entangled materials, steel wool, etc), mostly for their capacity to be engineered in order to present specific properties, inherent to their architecture. In this context, some studies have been carried out concerning entangled materials but only a few on monofilament entangled materials. Such a material, with no filament ends, could exhibit interesting properties for shock absorption, vibration damping and ductility. In this study, entanglements were manually produced, using different types of wire, and submitted to constrained (inside a PTFE die) in-situ compressive tests within the laboratory tomograph. This technique enabled a 3D, non destructive, microstructural characterization of the complex architecture of these materials, along with the analysis of their macroscopic mechanical properties. The stiffness of this material was found to be in a 20-200 MPa range and homogeneous samples could be obtained, while lowering their stiffness, by pre-deforming the initial wire as a spring. Damping measurements were performed using different types of entanglements (constitutive materials, volume fraction, wire diameter, wire shape) under both monotonic and dynamic loadings and directly linked to the measurements of the number of contacts. The Dynamic Mechanical Analysis underlined the great capacity of this material to absorb energy with a loss factor of about 0.25 and damping was found to decrease with the stiffness of the entanglement. The mechanical properties of this material were first modeled using an analytical “beam” model based on the experimental evolution of the mean distance between contacts and a good agreement was found with the experimental results. In parallel, a Discrete Element Method was used in order to model the compressive behaviour of Monofilament Entangled Materials. Although purely elastic properties were taken into account in the model, a very good agreement with the experimental results was obtained by adjusting the friction coefficients of the model. This tends to prove that the plasticity of these entangled materials is rather due to the structure (friction) than to the constitutive material itself.

Les matériaux cellulaires sont des échantillons à très forte porosité qui peuvent être décrits à deux échelles : la mésostructure et la microstructure. Le lien entre l'architecture des matériaux et les propriétés mécaniques a été largement étudié dans la littérature. Les caractéristiques microstructurales peuvent avoir une influence importante sur les propriétés macroscopiques. Le but de ce travail est de relier les caractéristiques architecturales et microstructurales des matériaux cellulaires à leurs propriétés mécaniques grâce notamment à la tomographie aux rayons X. Une nouvelle approche combinant l'imagerie 3D à plusieurs résolutions, le traitement d'images et la modélisation éléments finis a permis de prendre en compte la microstructure de la phase solide. Quatre matériaux cellulaires modèles ont ainsi été étudiés : des mousses d'aluminium, des structures cellulaires périodiques en alliage de cobalt-chrome, des échantillons de β-TCP et des composites hydroxyapatite/β-TCP. Les matériaux métalliques ont été fournis par des collègues d'autres laboratoires, tandis que les matériaux céramiques ont été fabriqués dans le cadre de cette étude. Pour chaque type de matériaux (métaux et céramiques), une structure régulière et une stochastique ont été comparées. Pour utiliser la méthode multi-échelle développée dans ce travail, les échantillons ont d'abord été scannés grâce à la tomographie locale dans laquelle l'échantillon est placé près de la source de rayons X. La tomographie locale permet de scanner la petite partie irradiée de l'échantillon et d'obtenir une image agrandie par rapport aux images à plus basse résolution. Ces images permettent d'observer certains détails de la phase solide non visibles à plus basse résolution. Différentes étapes de traitement d'images ont ensuite été mises en œuvre pour obtenir une image à basse résolution incluant les informations provenant des images à haute résolution. Ceci a été réalisé grâce à une série d'opération de seuillage et sous-résolution des images à haute résolution. Le résultat de ces différentes étapes de traitement d'images donne une image de l'échantillon initial à basse résolution mais qui inclut l'information supplémentaire décelée à haute résolution. Ensuite, des essais mécaniques in situ ont été réalisés dans le tomographe pour suivre à basse résolution l'évolution des échantillons pendant la déformation. Les images initiales citées plus haut ont été utilisées pour produire des maillages éléments finis. Des programmes Java ont été adaptés pour créer des fichiers d'entrée pour les modèles éléments finis à partir des images initiales et des maillages. Les images initiales contenant les informations à propos de la phase solide, les images des essais mécaniques et les modèles éléments finis ont permis d'expliquer le comportement mécanique des échantillons en reliant les sites d'endommagement expérimentaux et les lieux de concentrations de contraintes calculés. / Cellular materials are highly porous systems for which two scales are mainly important: the mesostructure and the microstructure. The mesostructure corresponds to the architecture of the materials: distribution of solid phase “walls” and macroporosity and can be characterized by X-ray tomographic low resolution images. The link between the architecture of the materials and the mechanical properties has been frequently studied. The microstructure refers to the characteristics of the solid phase. Its microstructural features (presence of a secondary phase or of defects due to the sintering) can have a strong influence on the macroscopic properties. The aim of this work is to link the morphological and microstructural features of metallic and ceramic based cellular materials and their mechanical properties thanks to X-ray tomography and finite element modelling. A new method combining X-ray tomography at different resolutions, image processing and creation of finite element modelling enabled to take into account some microstuctural features of the cellular samples. Four different cellular materials were studied as model materials: aluminium foam fabricated by a liquid state process, cobalt periodic structures made by additive manufacturing, β-TCP porous samples fabricated by conventional sacrificial template processing route and hydroxyapatite/β-TCP composites made by additive manufacturing (robocasting). The metal based materials were provided by colleagues while the ceramic based porous materials were fabricated in the frame of the current study. For each type (metals or ceramics), a stochastic and a regular structure have been compared. For implementing the multiscale method developed in this work, the samples were firstly scanned in a so called “local” tomography mode, in which the specimen is placed close to the X-ray source. This allowed to reconstruct only the small irradiated part of the sample and to obtain a magnified image of a subregion. These images enable to observe some details which are not visible in lower resolution. Different image processing steps were performed to generate low resolution images including microstructural features imaged at high resolution. This was done by a series of thresholding and scaling of the high resolution images. The result of these processing steps was an image of the initial sample. Then, in situ mechanical tests were performed in the tomograph to follow the deformation of the sample at low resolution. The above mentioned initial images were used to produce finite element meshes. Special Java programs were adapted to create finite element input files from initial images and meshes. The initial images containing information about the solid phase, the images from the mechanical tests and the finite element models were combined to explain the mechanical behaviour of the sample by linking the experimental damage locations in the sample and the simulated stress concentration sites.

This study is devoted to the thermo-mechanical constitutive modeling for saturated stiff clays and the development of a corresponding efficient stress integration algorithm. The mechanical behavior of natural Boom Clay in isothermal conditions was first characterized. The Modified Cam Clay model (MCC) was then applied to simulate the natural Boom Clay behavior. It has been found that the MCC gives poor-quality predictions of the natural Boom Clay behavior. Thereby, an adapted Cam Clay model (ACC-2) was developed by introducing a new yield surface and a new plastic potential as well as a Two-surface plastic mechanism. This model allows satisfactory prediction of the main features of the mechanical behavior of natural Boom Clay. Moreover, the constitutive equations of this model can be formulated mathematically as in a classic elasto-plastic model. Thus, the classic stress integration algorithm can be applied. The thermal effects were considered by assessing the performance of some advanced thermo-mechanical models (Cui et al., 2000; Abuel-Naga et al., 2007; Laloui and François, 2008; 2009). It appears that all the three models can capture the main features of the thermo-mechanical behavior of saturated clays. However, each constitutive model has its own limitations or unclear points from the theoretical point of view. The stress integration algorithm of the thermo-mechanical model proposed by Cui et al. (2000) at the stress point level was also developed using a specifically designed adaptive time-stepping scheme. The computation time required to achieve a given accuracy is largely reduced with the adaptive sub-stepping considered for both mechanical and thermal loadings. A Two-surface thermo-mechanical model (TEAM model) was developed based on the Two-surface plastic mechanism. The proposed model extends the model of Cui et al. (2000) to a Two-surface formulation, considering the plastic strain coupling between the thermal and the mechanical loading paths. The simulation of drained tests shows that this model can capture the main thermo-mechanical features of natural Boom Clay along different loading paths. The TEAM model was finally extended to undrained conditions. After setting up an appropriate effective stress principle and defining a volumetric strain condition, the undrained heating process was analyzed. The validity of the thermo-hydro-mechanical constitutive equations was examined based on the data from typical tests

Granular materials are characterized by its discrete nature which makes their behaviour very complex to understand when subjected to various loading situations. Comparing other materials, the understanding of granular materials is poor. This is because experimental analysis provides the macroscopic responses of the considered assembly whereas the discrete nature of the particles point to the fact of understanding the micro scale details and correlating it with the overall behaviour. Among the various modeling tools viz. analytical, physical or numerical, Discrete Element Method (DEM) a numerical technique, originally developed by Cundall (1971, 1974) and modified by Cundall and Strack (1979a, 1979b) is widely used for granular materials. Later a thorough description of DEM was given by Cundall (1988) and Hart et al (1988). Moreover Cundall & Hart (1992) reported discrete element code as one which allows finite displacements and rotations of discrete bodies along with recognition of new contacts as the calculation progresses which is followed in particle flow code and is used for this study. Generally the discrete particles are modeled as discs or spheres in 2-D and 3-D simulations respectively. The discs or spheres were considered as it is very easy to characterize the grain interactions and the contact detection. Even though the significance of particle shape has been reported in literature, a comprehensive 3-D study of the effect of particle shape on the various aspects of soil behaviour is lacking and is not reported. Particle shape is generally defined in terms of elongation, roundness and texture. Elongation is an indication of the particle aspect ratio whereas roundness measures the sharpness or angularity of particle’s edges and corners. Texture is related to the roughness of the surface. Particle gradation also plays a role in the mechanical behaviour. The influence of each of these factors on the mechanical behaviour of the assembly is important. Hence the major factors like elongation, texture etc which are used to define the particle shape are incorporated in this study. Evaluating the particle shape is another hurdle. In this study, the particle shape is analyzed using a 3D laser scanner which helps to identify the major and minor axis lengths of the particle. Additionally, the effective use of 3D DEM on large scale real life applications incorporating the particle shape effect is also not dealt with very extensively. Hence in this research, a comprehensive study on the calibration of DEM using glass beads, effect of particle shape on drained and undrained monotonic behaviour, liquefaction, post liquefaction and dynamic properties and the application of DEM to a grain polishing machine and an underground tunnel assembly is presented. In the present study, a set of drained triaxial tests were done on glass bead assembly using a laboratory triaxial set up. The glass beads used for the test were spherical and ellipsoidal in shape. The shape of glass beads was characterized through a sophisticated method of 3D laser scanning. In this scanning, the shape of the image of the glass bead is captured through an array of digitized points. These images obtained as unstructured 3D triangular meshes on processing will render the long and short axes of the particle which can be used for proper modeling of the particle shape. After obtaining the long and short axes for the particles, the same is used for the numerical modeling of the glass beads. The numerical simulation results have shown that the assembly modeled with clumped particles gives results qualitatively and quantitatively similar to the observed experimental macro responses. Hence this is used to demonstrate the power of DEM to realistically model the granular behaviour by incorporating the particle shape effect. In addition, undrained simulation of granular material has been numerically predicted from drained triaxial tests which are similar to the approach proposed by Norris et al (1997). An excellent correlation between undrained results predicted from drained triaxial test and undrained test (performed under constant volume conditions) has been observed. This further underlines the fact that the constant volume simulations are equivalent to undrained tests. Having validated the DEM results to the values obtained from the experiments on glass beads, a series of monotonic drained and undrained triaxial tests were performed on cylindrical assemblies of height to diameter ratio 2:1. Four different sets of assemblies were prepared which consists of particles of different aspect ratio to study the influence of particle shape. The behaviour of these assemblies under drained shearing indicates that the strength of the clumped assemblies is higher than that of the spherical assembly at all confining pressures. This has been explained from the magnitude of the anisotropic vi coefficients associated with the fabric and normal contact force tensors. It is also noted that eventhough both assemblies reach the peak strength at the same axial strain, the strain softening associated with the clumped assembly is very rapid which is due to the fact that clumps try to push each other apart as it offers more resistance to rotation resulting in dilation. Another significant observation is that a general increase in aspect ratio will not keep on increasing the strength. As the aspect ratio increases, the particles have a tendency to orient along their larger dimensions. This helps them to attain the lowest potential energy leading to a stable equilibrium and resulting in inherent fabric anisotropy. But when the particles try to align along the larger dimension, the formation of strong contact forces along the direction of loading is hindered. In addition, the lower strength associated with the higher aspect ratio particle assembly can also be attributed to the formation of unexpected void spaces when these longer particles bridge gaps over the underlying grains. The critical state studies indicate that the clumped assembly is having a higher residual strength compared to that of the spherical assembly. In the case of clumped assemblies also, irrespective of the initial loose or dense state of the assembly and the confining pressure applied, the samples reached the same critical state which underlines that the critical state is unique for a granular material. Moreover, the critical state line is non-linear for both the spherical and clumped assemblies. The studies conducted on the liquefaction behaviour indicates that at lower confining pressures the assemblies with particles consisting of lower aspect ratios loses its strength at less number of cycles which can be attributed to the interlocking of non-spherical particles resulting in higher number of contacts per particle. Moreover, during the initial cycles of loading, it is seen that the strong contacts are predominantly in the vertical direction or more precisely along the direction of maximum axial strain. But the plots extracted at higher cycles indicated that the strong contacts along the vertical direction have diminished considerably. This reduction in contact force magnitude and force chain indicates a drop in the number of contacts and is clearly visible in the gradual decrease of average coordination number. Another significant observation is that as the confining stresses increases, the rate of pore pressure generation of the assembly vii consisting of only spherical particles is less compared to the other two samples. Furthermore, at higher confining pressures, when the load direction reverses, the fabric of the clumped assemblies fails to change to a new orientation immediately. But to retain equilibrium the force anisotropy will quickly adjust itself. This mismatch results in losing the contacts and resulting in lower strength and less resistance to liquefaction at higher stresses for assemblies consisting of clumped particles. The post liquefaction study of the numerically liquefied samples shows that the assembly consisting of clumped shaped particles gained strength at a much faster rate compared to the assembly consisting of only spheres. This may be attributed to the ability of the clumps to rearrange themselves on a faster rate compared to that of the spherical particles. The rate of development of average coordination number is very significant as it explains the ability of the assembly to build up the deviatoric stress from a complete collapsed structure. As the contacts develop, the average coordination number as well as the deviatoric stress starts increasing with both the values higher for the assembly consisting of clumped particles. The evaluation of the dynamic properties viz. shear modulus and damping ratio showed a trend similar to the experimental observations on real granular materials. It is observed that the normalized shear modulus reduces with an increase in shear strain and the rate of reduction is very high at low strains for all the samples. It can be seen that as the confining pressure increases, the normalized shear modulus value also increases and the rate of increment is higher for the assemblies consisting of non-spherical particles. Furthermore, for all the samples the threshold shear strain is about 0.001 up to which the behaviour is elastic. Beyond the threshold shear strain, the variation of the normalized shear modulus ratio is non-linear. At small shear strains, the energy dissipation is low resulting in smaller values of damping. As the strains increase, the non-linearity in the constitutive behaviour results in higher material damping leading to high damping value. The simulations of the food polishing machine helps to understand the pattern of hitting of clumped grains on the wall with due importance to the velocity of hit, angle of hit, force of hit, and the number of grains hitting the wall. The modeling and subsequent extraction of the data helped to identify that the wear and tear of the machine was not uniform and was clustered to specific regions due to the non-uniform distribution of the considered parameters. This helped to design a more sophisticated system such that the parts which are subjected to more deterioration are provided with additional support. To bring out the effect of the particle shape, simulations are performed using spherical particles and the results show that the pattern of variation is same, but the magnitudes are different owing to the less surface area associated with the spherical particles. The 3-D simulations of an underground tunnel assembly in a weak weathered rock helped to understand the variation in the stability of the system with and without lining. It was observed that the introduction of lining resulted in a more stable configuration and the circumferential stresses were found to be distributed uniformly around the tunnel. FEM simulations also show a similar trend of stress and strain variations but were unable to capture the ground loosening around the tunnel and the formation of the ground arch whereas DEM could realistically capture all these phenomena. It was observed that as the shape changes from sphere to non-spherical particles, circumferential stresses around the tunnels increased. In addition, as the distance from the tunnel face increases, the strains are reduced. The maximum vertical strain near the crown of the tunnel is observed for the assembly consisting of spherical particles. In short, this research is focused on a comprehensive understanding of the particle shape effect on the mechanical behaviour of granular mass. Numerical simulations incorporating the shape effect has been done on drained and undrained monotonic shear tests, critical state, liquefaction, post liquefaction and dynamic properties. Besides, the granular dynamics simulation of the movement of long food grains in a food polishing machine and the behaviour of an underground tunnel in a granular assembly is also reported.

Dans le domaine du renforcement et/ou de la réparation des structures en béton armé par des matériaux composites à l'aide de la méthode du collage extérieur au moyen d'un adhésif époxy, une des préoccupations de la communauté scientifique est l'intégrité structurelle de ce système dans le cas d'incendie dans lequel la haute température est une caractéristique essentielle et peut atteindre jusqu'à 1200°C. Ce travail de recherche est axé sur le comportement thermomécanique à très haute température des matériaux composites [un composite à base de polymère carbone/ époxy (Carbon Fiber Reinforced Polymer- CFRP), un composite textile/ mortier cimentaire (Textile Reinforced Concrete- TRC) et un adhésif à base d'époxy]. L'évolution des propriétés mécaniques et d'autres aspects mécaniques de ces matériaux composites avec la température a été caractérisée. Une nouvelle procédure expérimentale concernant la mesure de la déformation de l'éprouvette à l'aide du capteur laser est développée et validée. Une étude numérique et expérimentale a été réalisée dans le but de déterminer principalement la température à la rupture des joints « composite/ adhésif/ composite » sous les sollicitations mécaniques et thermiques. L'efficacité de la protection thermique de deux isolants [PROMASPRAY®T (produit commercial de la société PROMAT] et Isolant A (produit développé par le LGCIE site Tusset) a aussi été étudiée dans cette thèse. Enfin, une approche numérique, à l'aide du logiciel ANSYS, est utilisée afin de déterminer, de façon préliminaire et approximative, à l'échelle matériau, les propriétés thermiques des matériaux (composite textile/ mortier cimentaire -TRC et Isolant A) / In the area of the strengthening and/or the reparation of reinforced concrete structures with composites by means of the external bonding method using an epoxy adhesive, one of the preoccupation of the scientific community is the structural integrity of this system in the event of fire in which the high temperature is the essential feature et can reach up to 1200°C. This research focuses on the thermo-mechanical behavior of composite materials [carbon/epoxy adhesive composite (or carbon fiber reinforced polymer (CFRP), textile/cementitious mortar composite (or textile reinforced concrete (TRC)] and an epoxy-based adhesive. The evolution of mechanical properties and other mechanical aspects of these materials with the temperature has been characterized. A new experimental procedure concerning the measurement of sample strain by the laser sensor is developed and validated. An experimental and numerical study has been realized in order to mainly determine the temperature at the failure of "composite/adhesive/composite" joints under thermal and mechanical loadings. The effectiveness of the thermal protection of two insulators [PROMASPRAY®T (a commercial product of the PROMAT company and the insulator A (product developed by the LGCIE site Tuset)] has also been investigated in this PhD thesis. Finally, a numerical approach, using ANSYS software, is used to determine, in the preliminary and approximate way, at material scale, thermal properties of the materials [the textile reinforced concrete (TRC) and the insulator A]