Spelling suggestions: "subject:"mécanique dess milieu continus"" "subject:"mécanique deus milieu continus""
1 |
Modélisation des effets de l'évolution microstructurale sur le comportement mécanique du superalliage monocristallin AM1Gaubert, Anaïs 30 November 2009 (has links) (PDF)
Ce travail s'inscrit dans le contexte de la modélisation et de la prévision de la durée de vie des aubes de turbines haute pression des moteurs aéronautiques. Ces pièces sont réalisées en superalliage monocristallin base nickel tel que l'AM1, matérieu de l'étude. Dans le cadre de ce travail, nous nous sommes intéressés aux évolutions microstructurales se produisant à haute température sous chargement mécanique connues sous le nom de coalescence orientée des précipités. Ce travail a consisté, dans un premier temps, a étudier le comportement du matérieu non-vieilli. Des essais à 950°C ont été réalisés afin d'enrichir la base d'essais existente. Ils ont permis l'identification d'un modèle de viscosité de type sinus hyperbolique sur une large gamme de vitesses de sollicitations. Le comportement initial du matériau a également été étudié à l'échelle mésoscopique (échelle des phases). Un modèle de comportement a été identifié pour chacune des phases afin de reproduire le comportement macroscopique du matériau. Cette étude, effectuée dans le cadre de la viscoplasticité classique a permis de mettre en évidence les limites de cette approche. Notamment, elle ne permet pas de modéliser les effets d'échelle observés en plasticité. Parallèlement, nous nous sommes intéressés au phénomène de coalescence orientée, tant d'un point de vue expérimental que numérique. Des essais mécaniques après vieillissement ont été réalisés. Différentes conditions de vieillissement ont été étudiées, en fluage à différentes température et contraintes et suivant différentes orientations cristallographiques mais également sous chargement cyclique. Les essais mécaniques ont montré un effet prépondérant de la coalescence orientée sur l'écrouissage du matériau, allant dans le sens d'un adoucissement. La coalescence orientée a également été modélisée par la méthode des champs de phases. Nous avons proposé un couplage du modèle champs de phases avec un modèle de comportement de viscoplasticité cristalline afin de prendre en compte l'influence de l'activité plastique dans les couloirs de matrice. Les simulations ont montré une accélération de la cinétique de mise en radeaux due à la plasticité ainsi que des précipités ayant une forme plus réaliste vis-à-vis de l'expérience, notamment en vieillissement sous chargement cyclique lent. La prise en compte des effets de la coalescence orientée passe par la modélisation des effets d'échelle en plasticité. C'est pourquoi nous nous sommes intéressés aux milieux continus généralisés. Nous avons étudié analytiquement les effets d'échelle produits par un modèle de Cosserat et le modèle microcurl dans le cas d'un matériau biphasé en cisaillement. Les deux modèles prévoient un écrouissage cinématique linéaire dépendant de la taille de la microstructure. Enfin, les développement précédents et les essais réalisés dans le cadre de ce travail ont permis la construction et l'identification d'un modèle macroscopique étendu qui rend compte des effets de la coalescence orientée sur le comportement mécanique de l'AM1.
|
2 |
Méthodes de contrôle de la qualité de solutions éléments finis (application à l'acoustique)Bouillard, Philippe 05 December 1997 (has links)
This work is dedicated to the control of the accuracy of computational simulations of sound propagation and scattering. Assuming time-harmonic behaviour, the mathematical models are given as boundary value problems for the Helmholtz equation <i>Delta u+k2u=0 </i> in <i>Oméga</i>. A distinction is made between interior, exterior and coupled problems and this work focuses mainly on interior uncoupled problems for which the Helmholtz equation becomes singular at eigenfrequencies.
As in other application fields, error control is an important issue in acoustic computations. It is clear that the numerical parameters (mesh size h and degree of approximation p) must be adapted to the physical parameter k. The well known ‘rule of the thumb’ for the h version with linear elements is to resolve the wavelength <i>lambda=2 pi k-1</i> by six elements characterising the approximability of the finite element mesh. If the numerical model is stable, the quality of the numerical solution is entirely controlled by the approximability of the finite element mesh. The situation is quite different in the presence of singularities. In that case, <i>stability</i> (or the lack thereof) is equally (sometimes more) important. In our application, the solutions are ‘rough’, i.e., highly oscillatory if the wavenumber is large. This is a singularity inherent to the differential operator rather than to the domain or the boundary conditions. This effect is called the <i>k-singularity</i>. Similarly, the discrete operator (“stiffness” matrix) becomes singular at eigenvalues of the discretised interior problem (or nearly singular at damped eigenvalues in solid-fluid interaction). This type of singularities is called the <i>lambda-singularities</i>. Both singularities are of global character. Without adaptive correction, their destabilizing effect generally leads to large error of the finite element results, even if the finite element mesh satisfies the ‘rule of the thumb’.
The k- and lambda-singularities are first extensively demonstrated by numerical examples. Then, two <i>a posteriori</i> error estimators are developed and the numerical tests show that, due to these specific phenomena of dynamo-acoustic computations, <i>error control cannot, in general, be accomplished by just ‘transplanting’ methods that worked well in static computations</i>. However, for low wavenumbers, it is necessary to also control the influence of the geometric (reentrants corners) or physical (discontinuities of the boundary conditions) singularities. An <i>h</i>-adaptive version with refinements has been implemented. These tools have been applied to two industrial examples : the GLT, a bi-mode bus from Bombardier Eurorail, and the Vertigo, a sport car from Gillet Automobiles.
As a conclusion, it is recommanded to replace the rule of the thumb by a criterion based on the control of the influence of the specific singularities of the Helmholtz operator. As this aim cannot be achieved by the <i>a posteriori</i> error estimators, it is suggested to minimize the influence of the singularities by modifying the formulation of the finite element method or by formulating a “meshless” method.
|
3 |
Rupture et délamination de films minces / Fracture and delamination in thin filmsMarthelot, Joel 18 June 2014 (has links)
Les films minces sont abondamment utilisés pour protéger ou fonctionnaliser les surfaces. Les contraintes mécaniques peuvent cependant entraîner la ruine de ces structures par la propagation de fissures ou leur décollement. Nous présentons deux systèmes expérimentaux qui amènent à s'interroger sur les limites de validité des critères de stabilité communément admis. Pour une couche de sol- gel présentant une énergie d'adhésion modérée avec son substrat de silicium, nous observons la propagation de fissures dans des films réputés stables. Les motifs de propagation sont étonnement réguliers (en forme de spirales ou d'allées de croissants). Nous caractérisons expérimentalement cette propagation. Nous mettons en évidence un nouveau mode de propagation de fissure où la délamination et la fissure transverse collaborent. Nous montrons l'existence d'une taille d'interaction robuste entre fissures indépendantes du chargement et des énergies de rupture. Nous proposons un modèle simple permettant de prédire les conditions de propagation et la sélection des motifs. Nous caractérisons enfin la dynamique de propagation de telles fissures. Dans un second système expérimental motivé par les applications en électronique flexible, nous nous intéressons à des films nanométriques d'oxyde de zinc déposés sur des substrats de polymère. La rigidité du substrat est alors faible devant celle de la couche. Nous décrivons la fragmentation du film et proposons une extension d'un modèle de couche élastique sur fondation élasto-plastique. Enfin, nous nous intéressons à la compression et à la fragmentation d'un multicouche comportant une fine épaisseur d'argent. / Nanometer to micrometer thin film coatings are extensively used in material science to protect and functionalize surfaces. Coating processes however generally result into residual strains and may induce the propagation of cracks with dramatic consequences on film reliability. We present two experimental models which challenges the usual stability criteria for thin films. For a thin sol-gel layer moderatly adherent to a silicon substrate, we observe fracture propagation in coatings considered stable. Moreover, intriguing cracks of regular geometry are found to self-replicate an initial triggering pattern independently from the local material properties. We present a novel mode of fracture which involves the simultaneous progression of both, a crack and a front of delamination. We describe the fracture mechanism with universal energy arguments and provide a simple model describing the crack path. We eventually characterize the dynamics of the propagation. In a second model experiment, we study the fragmentation of coating deposited on a compliant substrate (motivated by flexible electronics application). We describe the fragmentation and propose an extension to an elastic layer deposited on an elastic foundation to an elasto-plastic case. Eventually, we study the compression and fragmentation of a zinc oxyde and silver multilayer
|
4 |
Contribution to the improvement of meshless methods applied to continuum mechanics / Contribution à l’amélioration des méthodes sans maillage appliquées à la mécanique des milieux continusFougeron, Gabriel 02 October 2018 (has links)
Cette thèse présente un cadre général pour l’étude de schémas de discrétisation nodaux sans maillageformulé en termes d’opérateurs discrets définis sur un nuage de points : intégration volumique et de bord, gradientet opérateur de reconstruction. Ces définitions dotent le nuage de points d’une structure plus faible que celledéfinie par un maillage, mais partageant avec elle certain concepts fondamentaux. Le plus important d’entre euxest la condition de compatibilité intégro-différentielle. Avec la consistance linéaire du gradient discret, cet analoguediscret de la formule de Stokes constitue une condition nécessaire à la consistance linéaire des opérateurs elliptiquesen formulation faible. Sa vérification, au moins de manière approchée, permet d’écrire des discrétisations dont le tauxde convergence est optimal. La construction d’opérateurs discrets compatibles est si difficile que nous conjecturons– sans parvenir à le démontrer – qu’elle nécessite soit la résolution d’un système linéaire global, soit la constructiond’un maillage : c’est "la malédiction sans-maillage". Trois grandes approches pour la construction d’opérateursdiscrets compatibles sont étudiées. Premièrement, nous proposons une méthode de correction permettant de calculerl’opérateur gradient compatible le plus proche – au sens des moindres carrés – sans mettre à mal la consistancelinéaire. Dans le cas particulier des gradients DMLS, nous montrons que le gradient corrigé est en réalité globalementoptimal. Deuxièmement, nous adaptons l’approche SFEM au cadre opérateur et constatons qu’elle définit desopérateurs consistants à l’ordre un et compatibles. Nous proposons une méthode d’intégration discrète exploitantla relation topologique entre les cellules et les faces d’un maillage qui préserve ces caractéristiques. Troisièmement,nous montrons qu’il est possible de générer tous les opérateurs sans maillage rien qu’avec la donnée d’une formuled’intégration volumique nodale en exploitant la dépendance fonctionnelle des poids volumiques nodaux par rapportà la position des noeuds du nuage, l’espace continu sous-jacent et le nombre de noeuds. Les notions de consistance desdifférents opérateurs sont caractérisées en termes des poids volumiques initiaux, formant un jeu de recommandationpour la mise au point de bonnes formules d’intégration. Dans ce cadre, nous réinterprétons les méthodes classiquesde stabilisation de la communauté SPH comme cherchant à annuler l’erreur sur la formule de Stokes discrète.L’exemple des opérateurs SFEM trouve un équivalent en formulation volume, ainsi que la méthode d’intégrationdiscrète s’appuyant sur un maillage. Son écriture nécessite toutefois une description très précise de la géométriedes cellules du maillage, en particulier dans le cas où les faces ne sont pas planes. Nous menons donc à bienune caractérisation complète de la forme de telles cellules uniquement en fonction de la position des noeuds dumaillage et des relations topologiques entre les cellules, permettant une définition sans ambigüité de leur volume etcentre de gravité. Enfin, nous décrivons des schémas de discrétisation d’équations elliptiques utilisant les opérateurssans-maillage et proposons plusieurs possibilités pour traiter les conditions au bord tout en imposant le moinsde contraintes sur la position des noeuds du nuage de points. Nous donnons des résultats numériques confirmantl’importance capitale de vérifier les conditions de compatibilité, au moins de manière approchée. Cette simple recommandation permet dans tous les cas d’obtenir des discrétisations dont le taux de convergence est optimal. / This thesis introduces a general framework for the study of nodal meshless discretization schemes. Itsfundamental objects are the discrete operators defined on a point cloud : volume and boundary integration, discretegradient and reconstruction operator. These definitions endow the point cloud with a weaker structure than thatdefined by a mesh, but share several fundamental concepts with it, the most important of them being integrationdifferentiationcompatibility. Along with linear consistency of the discrete gradient, this discrete analogue of Stokes’sformula is a necessary condition to the linear consistency of weakly discretized elliptic operators. Its satisfaction, atleast in an approximate fashion, yields optimally convergent discretizations. However, building compatible discreteoperators is so difficult that we conjecture – without managing to prove it – that it either requires to solve a globallinear system, or to build a mesh. We dub this conjecture the "meshless curse". Three main approaches for theconstruction of discrete meshless operators are studied. Firstly, we propose a correction method seeking the closestcompatible gradient – in the least squares sense – that does not hurt linear consistency. In the special case ofMLS gradients, we show that the corrected gradient is globally optimal. Secondly, we adapt the SFEM approachto our meshless framework and notice that it defines first order consistent compatible operators. We propose adiscrete integration method exploiting the topological relation between cells and faces of a mesh preserving thesecharacteristics. Thirdly, we show that it is possible to generate each of the meshless operators from a nodal discretevolume integration formula. This is made possible with the exploitation of the functional dependency of nodal volumeweights with respect to node positions, the continuous underlying space and the total number of nodes. Consistencyof the operators is characterized in terms of the initial volume weights, effectively constituting guidelines for thedesign of proper integration formulae. In this framework, we re-interpret the classical stabilization methods of theSPH community as actually seeking to cancel the error on the discrete version of Stokes’s formula. The example ofSFEM operators has a volume-based equivalent, and so does its discrete mesh-based integration. Actually computingit requires a very precise description of the geometry of cells of the mesh, in particular in the case where its facesare not planar. We thus fully characterize the shape of such cells, only as a function of nodes of the mesh andtopological relations between cells, allowing unambiguous definition of their volumes and centroids. Finally, wedescribe meshless discretization schemes of elliptic partial differential equations. We propose several alternatives forthe treatment of boundary conditions with the concern of imposing as few constraints on nodes of the point cloudas possible. We give numerical results confirming the crucial importance of verifying the compatibility conditions,at least in an approximate fashion. This simple guideline systematically allows the recovery of optimal convergencerates of the studied discretizations.
|
5 |
Biomécanique de la croissance de la plaque d'athérosclérose : contribution à l'étude des contraintes résiduellesMesnier, Nicolas 22 June 2011 (has links) (PDF)
Ce travail de thèse a pour ambition d'apporter des éléments de compréhension du rôle des contraintes résiduelles au sein des plaques d'athérosclérose. Il s'inscrit dans un programme de recherche qui vise à développer de nouveaux outils favorisant l'aide au diagnostic de patients atteints de cette maladie cardiovasculaire. Suivant les évolutions majeures de ces dernières années, la première partie de ce travail présente une vision contemporaine de la maladie d'athérosclérose et des nombreux enjeux de modélisation mécanique auxquels la compréhension de la maladie et son diagnostic font appel. Nous faisons émerger un manque certain quant au rôle des contraintes résiduelles dans cette maladie. Les contraintes résiduelles sont un marqueur naturel du développement des tissus. La deuxième partie de ce travail aborde leur quantification expérimentale au sein des artères saines et pathologiques de souris. De façon intéressante nous montrons leur importance et la nécessité de les prendre en compte pour avoir une image raisonnable de l'état de contraintes des artères en conditions physiologiques. Ces résultats originaux n'ont toutefois pas permis d'établir de corrélation avec l'intensité du processus inflammatoire au sein des plaques. Enfin, la dernière partie de ce travail aborde la modélisation du développement des contraintes résiduelles. Un effort de description de la structure géométrique des corps avec des contraintes résiduelles a été fait afin de faire émerger de nouveaux outils de modélisation de leur développement. Une première application d'un modèle de croissance intrinsèque est présentée et illustrée. Cette approche semble prometteuse et apporte déjà quelques éléments de compréhension du développement des contraintes résiduelles. Toutefois, les travaux effectués ne constituent qu'une étape préliminaire et n'ont pas encore été exploités pleinement. Par contre ces travaux ouvrent assurément de nouvelles perspectives de recherche.
|
6 |
Biomécanique de la croissance de la plaque d'athérosclérose : contribution à l'étude des contraintes résiduelles / Biomechanics of atheroma plaque growth : contribution to the study of residual stressesMesnier, Nicolas 22 June 2011 (has links)
Ce travail de thèse a pour ambition d'apporter des éléments de compréhension du rôle des contraintes résiduelles au sein des plaques d'athérosclérose. Il s'inscrit dans un programme de recherche qui vise à développer de nouveaux outils favorisant l'aide au diagnostic de patients atteints de cette maladie cardiovasculaire. Suivant les évolutions majeures de ces dernières années, la première partie de ce travail présente une vision contemporaine de la maladie d'athérosclérose et des nombreux enjeux de modélisation mécanique auxquels la compréhension de la maladie et son diagnostic font appel. Nous faisons émerger un manque certain quant au rôle des contraintes résiduelles dans cette maladie. Les contraintes résiduelles sont un marqueur naturel du développement des tissus. La deuxième partie de ce travail aborde leur quantification expérimentale au sein des artères saines et pathologiques de souris. De façon intéressante nous montrons leur importance et la nécessité de les prendre en compte pour avoir une image raisonnable de l'état de contraintes des artères en conditions physiologiques. Ces résultats originaux n'ont toutefois pas permis d'établir de corrélation avec l'intensité du processus inflammatoire au sein des plaques. Enfin, la dernière partie de ce travail aborde la modélisation du développement des contraintes résiduelles. Un effort de description de la structure géométrique des corps avec des contraintes résiduelles a été fait afin de faire émerger de nouveaux outils de modélisation de leur développement. Une première application d'un modèle de croissance intrinsèque est présentée et illustrée. Cette approche semble prometteuse et apporte déjà quelques éléments de compréhension du développement des contraintes résiduelles. Toutefois, les travaux effectués ne constituent qu'une étape préliminaire et n'ont pas encore été exploités pleinement. Par contre ces travaux ouvrent assurément de nouvelles perspectives de recherche. / This graduate research has the ambition to bring some understandings of the role of residual stresses in atheroma plaques. It falls under a research program which aims at developing new tools in supporting the diagnosis of patients concerned with this cardiovascular disease. Considering the major changes in the past decade, the first part of this work presents a contemporary vision of the atherosclerosis disease and many mechanical modeling challenges brought by its understanding and diagnosis appeal. We emphasized a lack in the description of the role of residual stresses in this disease even if they are a natural marker of tissues development. The second part of this work approaches their experimental quantification within healthy and pathological mice arteries. Interestingly, we have shown their importance and the need to take them into account in order to have a reasonable image of the stresses state of arteries in physiological conditions. These original results however did not make it possible to establish correlation with the intensity of the inflammatory process within the plaques. Finally, the last part of this work approaches the modeling of the development of residual stresses. An effort in the description of the geometrical structure of material bodies with residual stresses allowed the emergence of new tools to model their development. A first application of an intrinsic growth model is presented and illustrated. This approach seems promising and already brings some key explaining elements in the development of residual stresses. However, carried out work constitutes only one preliminary stage and were not yet fully exploited. Hence, this work undoubtedly opens new prospects of research.
|
7 |
Modélisation du comportement des matériaux granulaires par des approches discrètes et continuesRahmoun, Jamila 07 December 2006 (has links) (PDF)
Les propriétés physiques des milieux granulaires trouvent leur origine à l'échelle locale des contacts entre les grains. Cette étude est consacrée à l'analyse de l'influence de ces contacts sur le comportement global du matériau à différentes échelles. Nous développons dans la première partie de ce mémoire, une approche continue qui s'affranchit des limitations de la théorie de Janssen et permet de calculer les contraintes dans un matériau granulaire ensilé. Cette approche permet également de représenter aussi bien qualitativement que quantitativement l'effet d'écrantage dans les silos. Dans une seconde partie, une modélisation du phénomène d'écrantage est effectuée à partir des simulations numériques discrètes réalisées avec le code MULTICOR. Nous calculons les contraintes moyennes s'exerçant au niveau des parois lors du remplissage d'un silo par un milieu granulaire polydisperse. Une bonne concordance est observée entre les résultats des simulations numériques discrètes et ceux de l'approche continue développée dans la première partie. Une autre conséquence de l'existence de contacts privilégiés entre les grains est l'anisotropie caractéristique des matériaux granulaires. Dans la dernière partie, nous développons une approche micromécanique qui permet de modéliser cette anisotropie par l'intermédiaire d'un tenseur de texture d'ordre quatre. Nous proposons une hypothèse cinématique générale qui est incorporée dans un schéma d'homogénéisation des milieux granulaires anisotropes et comparée, dans le cas isotrope, à des simulations numériques discrètes existant dans la littérature. Les tendances qualitatives et quantitatives des résultats sont tout à fait satisfaisantes.
|
8 |
Rupture et délamination de films mincesMarthelot, Joel 18 June 2014 (has links) (PDF)
Les films minces sont abondamment utilisés pour protéger ou fonctionnaliser les surfaces. Les contraintes mécaniques peuvent cependant entraîner la ruine de ces structures par la propagation de fissures ou leur décollement. Nous présentons deux systèmes expérimentaux qui amènent à s'interroger sur les limites de validité des critères de stabilité communément admis. Pour une couche de sol- gel présentant une énergie d'adhésion modérée avec son substrat de silicium, nous observons la propagation de fissures dans des films réputés stables. Les motifs de propagation sont étonnement réguliers (en forme de spirales ou d'allées de croissants). Nous caractérisons expérimentalement cette propagation. Nous mettons en évidence un nouveau mode de propagation de fissure où la délamination et la fissure transverse collaborent. Nous montrons l'existence d'une taille d'interaction robuste entre fissures indépendantes du chargement et des énergies de rupture. Nous proposons un modèle simple permettant de prédire les conditions de propagation et la sélection des motifs. Nous caractérisons enfin la dynamique de propagation de telles fissures. Dans un second système expérimental motivé par les applications en électronique flexible, nous nous intéressons à des films nanométriques d'oxyde de zinc déposés sur des substrats de polymère. La rigidité du substrat est alors faible devant celle de la couche. Nous décrivons la fragmentation du film et proposons une extension d'un modèle de couche élastique sur fondation élasto-plastique. Enfin, nous nous intéressons à la compression et à la fragmentation d'un multicouche comportant une fine épaisseur d'argent.
|
9 |
Application de la Méthode des Points Matériels aux phénomènes gravitaires / Application of the Material Point Method to gravitational phenomenaGracia Danies, Fabio 12 January 2018 (has links)
Dans les régions de montagne, la prévision des évènements gravitaires reste un défi pour la gestion des risques. Des méthodes de calcul telles que la méthode des éléments discrets (DEM), où les particules interagissent les unes avec les autres pour restituer un comportement global d’une masse granulaire, ont été utilisées pour aborder ce type de problématique. L’application de la DEM reste normalement limitée aux évènements de petits volumes impliquant un nombre de blocs plutôt faible, puisque les temps de calcul peuvent devenir rapidement prohibitifs avec l’augmentation du nombre de particules. Les méthodes de calcul continues sont donc une alternative intéressante car elles permettent de réduire les temps de calcul. Elles nécessitent cependant la définition d’une loi de comportement macroscopique capable de représenter correctement les principaux traits de comportement mécanique du matériau au sein de la masse. L'objectif principal du travail de thèse réside dans le développement d’un outil numérique permettant de modéliser certains aléas gravitaires tels que les écoulements en masse. Notre choix s’est porté sur une méthode Lagrangienne-Eulérienne (méthode des points matériels – MPM) capable de gérer de grandes déformations tout en bénéficiant des principaux avantages de la méthode des éléments finis (FEM). La méthode utilise une grille Eulérienne fixe sur laquelle se déplacent des points matériels pendant les simulations. Un outil numérique, nommé MPMbox (2D et 3D), a été développé entièrement durant la thèse en C++. Le code a été validé à l'aide d'une série de solutions analytiques en quasi-statique (tests géotechniques standards) ainsi que par des applications de la littérature incluant des déformations importantes et rapides (tests d'affaissement). Après validation, le code a été confronté aux prédictions d’un outil de calcul DEM (DEMbox) dans le cadre de simulations numériques impliquant l'écoulement (initiation, régimes transitoires, propagation et arrêt) d'un matériau granulaire (particules sphero-polyhédriques) sur un plan incliné. Les résultats ont été comparés en termes de distance de propagation, de forme du dépôt et d'énergies dissipées à l'interface et dans la masse pendant l'écoulement. Pour les applications qui ont suivies, des éléments discrets ont été couplés à la MPM afin qu'un bloc rigide (DEM) puisse interagir avec un sol déformable (MPM). Cette application a consisté en l'analyse (2D) de la collision entre un bloc rocheux rigide (rond ou carré) et un sol bicouche élastoplastique. Les investigations ont été largement basées sur la mesure de coefficients de restitution (rapport des énergies cinétiques avant et après impact) qui reste difficile à déterminer expérimentalement. / In mountainous regions, the prediction of gravitational phenomena remains a challenge for the management of risk. Computational methods such as the Discrete Element Method (DEM) have been used for the modeling of these types of phenomena, where particles interact with each other to give an overall behavior of the mass. Its application can be somewhat restricted to small and medium number of blocks, since the computational time can easily become too large. Continuum analyses are therefore an attractive approach, which can reduce the computational times, but that rely on a constitutive law to represent the behavior within the mass. The main objective of this PhD was to develop a numerical tool that allowed the modeling of some specific gravitational hazards, such as the flowing of mass. A Lagrangian-Eulerian method such as the Material Point Method (MPM) is able to handle large deformations, while preserving most of the capabilities of the Finite Element Method (FEM). The method uses an Eulerian grid which is only used as a numerical scratch-pad, and remains fixed during simulations. A numerical tool named MPMbox (2D and 3D) was then developed from the ground up using C++. The code was validated using a series of analytical solutions for quasi-static analysis (some standard geotechnical tests), as well as simulations including large and rather rapid deformations (slump tests). After validation, the code was first used to make a numerical comparison with the DEM. In the comparison, a parametric survey was carried out during which the flow of a granular material on a sloped surface was simulated. Results were compared in terms of run-out distance, spread of the deposit and energy dissipated at the interface and within the mass during the flow. For a second study, discrete elements were coupled with MPM so that a rigid block could interact with a deformable soil. This application consisted in the (2D) analysis of the collision between a discrete block (round and squared) and a bounded elasto-plastic double-layered soil (soft over hard layers). The investigations were largely based on the measurement of the restitution coefficient (ratio of kinetic energies before to after the impact), which cannot be easily determined experimentally.
|
10 |
Modélisation mécanique des tissus biologiques : application à la croissance des tumeurs solides et à la reconstruction multiéchelles des propriétés élastiques de la cuticule d'arthropode / Mechanical modeling of biological tissues : application to solid tumor growth and multiscale reconstruction of the elastic properties of arthropod cuticleLhadi, Safaa 21 September 2015 (has links)
De nos jours, l’enjeu de la mécanobiologie ne cesse de grandir. On s’intéresse à la description des problèmes biophysiques d’un point de vue mécanique avec des approches multiéchelles. Dans ce travail, nous proposons d’étudier deux exemples mettant en évidence le rôle important de la mécanique sur des processus purement biologiques. 1) La croissance tumorale dans son stade avasculaire : nous proposons un modèle continu où le tissu tumoral est considéré capable de croître et de se déformer tout en obéissant aux lois de conservation. Nous proposons ensuite pour étudier l’effet des propriétés mécaniques du microenvironnement -où réside la tumeur- sur le développement tumoral d’intégrer certaines conditions aux interfaces tumeur/microenvironnement. 2) La reconstruction des propriétés élastiques de la cuticule d’arthropode : nous proposons un modèle multiéchelles de son comportement mécanique fondé sur la structure hiérarchique établie dans la littérature. Pour remédier à la sous-estimation du modèle des propriétés élastiques de la cuticule, nous proposons d’inclure les interfaces à certaines échelles qui pourraient améliorer la transmission des efforts aux constituants multiéchelles du composite (cuticule) et donc améliorer les propriétés élastiques macroscopiques de ce dernier. / Nowadays, the challenge of mechanobiology keeps growing. We are interested in the description of biophysical problems from a mechanical point of view with multiscale approaches.In the present study, we propose to study two examples highlighting the substantial role of mechanics on purely biological processes. 1) Tumor growth in the avascular stage: we propose a continuous model where tumor tissue is considered able to grow and to deform while obeying to conservation laws. Then, we propose to study the effect of the mechanical properties of the microenvironment- where lives the tumor- on the tumor development by integration of certain interfaces conditions tumor/microenvironment. 2) Reconstruction of the elastic properties of the arthropod cuticle: we propose a multiscale model of its mechanical behavior based on the hierarchical structure established in the literature. To remedy the under-estimation of the cuticle elastic properties of the model, we propose to include the interfaces to some scales that could improve the transmission of forces to the multiscale components of the composite (cuticle) and thus improve their macroscopic elastic properties.
|
Page generated in 0.1069 seconds