L'os est un matériau vivant avec une structure hiérarchique complexe qui lui confère des propriétés mécaniques remarquables. L'os subit perpétuellement des contraintes mécaniques et physiologiques, ainsi sa qualité et sa résistance à la fracture évoluent constamment au cours du temps à travers le processus de remodelage osseux. La qualité osseuse est non seulement définie par la densité minérale osseuse mais également par les propriétés mécaniques ainsi que la microarchitecture. Dans le cadre de la présente thèse, on a développé une modélisation multiéchelle unifiée couplant à la fois les activités cellulaires au comportement mécanique de l'os tenant compte des différents niveaux hiérarchiques de l'os: de l'ultrastructure au remodelage osseux. Ce modèle permet d'étudier le comportement mécano-bibliologique de l'os et de prédire ses propriétés mécaniques apparentes à différentes échelles allant du nanoscopique au macroscopique en fonction des constituants élémentaires de l'os. Pour atteindre cet objectif, une démarche en quatre phases a été adoptée. La première phase consiste à décrire les constituants élémentaires de l'os. La deuxième phase avait pour objectif la modélisation multiéchelle de l'ultrastructure osseuse constituée de trois échelles nanoscopiques (microfibrille, fibrille et fibre) par la méthode des éléments finis et des réseaux de neurones. La troisième phase correspond à la modélisation des échelles micro-macroscopiques de l'os cortical (lamelle, ostéon, os cortical) en utilisant comme paramètres d'entrée les propriétés de la fibre déterminées dans la deuxième phase. Enfin, dans la dernière phase, on a développé un modèle mécano-biologique du remodelage osseux permettant de simuler le processus d'adaptation osseuse tenant compte explicitement des activités biologiques des cellules osseuses. Les propriétés mécaniques prédites par nos algorithmes multiéchelles ont servi pour alimenter le modèle de remodelage. Ce modèle a été implémenté au code de calcul d'éléments finis ABAQUS/Standard à travers sa routine utilisateur UMAT. Finalement, le modèle EF mécano-biologique multiéchelle du remodelage osseux a été appliqué pour simuler différents scénarii de remodelage sur des fémurs humains (2D et 3D). Différents facteurs ont été ainsi analysés tels que l'âge, le genre, l'amplitude des activités physiques, etc. Les résultats obtenus sont conformes (qualitativement) avec les observations cliniques et cohérents avec les différentes études expérimentales. En conclusion: (i) Les modèles unifiés ainsi développés (modèle multiéchelle, modèle mécano-biologique de remodelage osseux) contribuent à l'analyse fine du comportement de l'os humain. (ii) L'application des algorithmes a permis d'effectuer des essais virtuels pour analyser les effets combinés de nombreux facteurs caractérisant la qualité osseuse.

Ce travail présente quelques résultats concernant le comportement en temps fini et en temps long de méthodes numériques pour des équations stochastiques. On s'intéresse d'abord aux équations différentielles stochastiques de Langevin et de Langevin amorti. On montre un résultat concernant l'analyse d'erreur faible rétrograde de ses équations par des schémas numériques implicites. En particulier, on montre que l'erreur entre le générateur associé au schéma numérique et la solution d'une équation de Kolmogorov modifiée est d'ordre élevé par rapport au pas de discrétisation. On montre aussi que la dynamique associée au schéma numérique est exponentiellement mélangeante. Dans un deuxième temps, on étudie le comportement en temps long d'une discrétisation en temps et en espace d'une EDPS semi-linéaire avec un bruit blanc additif, qui possède une unique mesure invariante . On considère une discrétisation en temps par un schéma d'Euler et en espace par une méthode des éléments finis. On montre que la moyenne, par rapport aux lois invariantes (qui n'est pas forcément unique) associées à l'approximation, par des fonctions tests suffisamment régulières est proche de la quantité correspondante pour . Plus précisément, on étudie la vitesse de convergence par rapport aux différents paramètres de discrétisation. Enfin, on s'intéresse à une EDPS semi-linéaire avec un bruit blanc additif dont le terme non-linéaire est un polynôme. On étudie la convergence au sens faible d'une approximation en temps par un schéma de splitting implicite.

Le développement de capillaires sanguins à partir d'un réseau pré-existant, l'angiogenèse, joue un rôle fondamental dans de nombreux contextes physiopathologiques, tels la cicatrisation des tissus ou le développement d'une tumeur solide. Dans cette thèse, nous nous sommes intéressés à la régulation de ce phénomène par les facteurs mécaniques (rigidité, viscosité, traction cellulaire). Dans un dialogue permanent entre l'expérimentation et la modélisation, nous avons développé un modèle théorique biomécanique minimal des premières étapes de l'angiogenèse in vitro, où l'angiogenèse est supposée issue d'une instabilité mécanique entre les forces actives de traction cellulaire et la résistance passive viscoélastique de la matrice extracellulaire. Notre modèle consiste en un système d'équations aux dérivées partielles non-linéaires couplées, résolu par la méthode des éléments finis. Nous avons mené des analyses de stabilité linéaire et non-linéaire de l'état d'équilibre homogène pour déterminer les points de bifurcation du système correspondant à une instabilité de Turing. Nous avons ensuite effectué une étude approfondie de l'influence des différents paramètres sur la formation du réseau. Les résultats des simulations numériques sont comparés avec succès aux résultats expérimentaux, obtenus par notre équipe ou extraits de la littérature. Dans une seconde partie de nos travaux, nous avons étudié des voies de régulation possibles, par les effets mécaniques, de la dégradation de la matrice extracellulaire. Nous avons alors montré que la régulation mécanique de la dégradation pouvait être un processus clé de l'angiogenèse in vitro.

Cette thèse analyse le spectre d'opérateurs de Schrödinger avec champ magnétique constant dans des ouverts de type diédraux. Pour comprendre l'influence d'une arête courbe sur la première valeur propre de l'opérateur dans la limite semi-classique, il faut connaître le bas du spectre de l'opérateur de Schrödinger magnétique avec champ constant sur un dièdre infini. Par transformation de Fourier ce problème se ramène à l'étude d'une famille d'opérateurs à paramètre sur un secteur infini. On calcule le spectre essentiel de ces opérateurs sur le secteur et on montre que dans certains cas il y a des valeurs propres discrètes sous le spectre essentiel. Par comparaison avec des opérateurs de Sturm-Liouville singuliers sur le demi-axe on obtient des majorations du bas du spectre de l'opérateur sur le dièdre : pour un angle d'ouverture assez petit et certaines orientations du champ magnétique, celui-ci est strictement inférieur aux quantités spectrales issues du cas régulier. Finalement on applique ces résultats à l'opérateur de Schrödinger avec champ magnétique constant et petit paramètre dans des domaines bornés de l'espace possédant des arêtes courbes. Pour déterminer une asymptotique de la première valeur propre dans la limite semi-classique, on construit des quasi-modes près de l'arête à l'aide des fonctions propres du problème à paramètre sur le secteur. En utilisant une partition du domaine selon que l'on soit près de l'arête ou du bord régulier, on obtient le premier terme de l'asymptotique pour diverses orientations du champ magnétique et on montre dans certains cas que la première valeur propre est inférieure aux valeurs propres associées à des ouverts réguliers.

L'objectif de ce mémoire de recherche est de déterminer les propriétés acoustiques des milieux poroélastiques à partir d'une démarche multi-échelle et multi-physique. Il traite d'échantillons réels de mousses, à cellules ouvertes ou partiellement fermées, dont les propriétés microstructurales sont caractérisées par des techniques d'imagerie. Cette information est utilisée afin d'identifier une cellule périodique idéalisée tridimensionnelle, qui soit représentative du comportement acoustique du milieu poreux réel. Les paramètres gouvernant les propriétés acoustiques du milieu sont obtenus en appliquant la méthode d'homogénéisation des structures périodiques. Dans une première étape, la structure des mousses est supposée indéformable. Il a été montré que pour le cas d'une distribution étroite de tailles caractéristiques de la géométrie locale, le comportement macroscopique d'une mousse à cellule ouverte peut être calculé à partir des propriétés géométriques locales de manière directe. Dans le cas d'une distribution étendue, le comportement acoustique du milieu est gouverné par des tailles critiques qui sont déterminées à partir de la porosité et de la perméabilité statique pour une mousse à cellules ouvertes ; pour une mousse à cellules partiellement fermées il est nécessaire d'identifier en plus une dimension connue de la géométrie locale. Nos résultats sont comparés avec succès à des données expérimentales obtenues par des mesures au tube d'impédance. Dans une seconde étape, les propriétés élastiques effectives du milieu poreux sont déterminées. Une modélisation par éléments finis de la cellule représentative a été mise en œuvre. Les paramètres élastiques calculés sont finalement comparés avec les données de la littérature, ainsi qu'à des essais mécaniques

Les phénomènes de couplage électromagnétique entre les lignes aériennes de transmission d'énergie et des structures voisines sont inévitables, surtout dans les zones densément peuplées. Les effets indésirables découlants de cette proximité sont variés, allant de l'établissement des tensions dangereuses à l'apparition de phénomènes de corrosion dus au courant alternatif. L'étude de cette classe de problèmes est nécessaire pour assurer la sécurité dans les alentours de la zone d'interaction et aussi pour préserver l'intégrité des équipements et des dispositifs présents. Cependant, la modélisation complète de ce type d'application implique la représentation tridimensionnelle de la région d'intérêt et nécessite des méthodes numériques de calcul de champs spécifiques. Dans ces travaux, des problèmes liés à la circulation de courants électriques dans le sol (ou de couplage dit conductif) seront abordés avec la méthode des éléments finis. D'autres problèmes résultants de la variation temporelle des champs électromagnétiques (ou de couplage dit inductif) seront aussi considérés et traités avec la méthode PEEC (Partial Element Equivalent Circuit) généralisée. Plus précisément, une condition limite particulière sur le potentiel électrique est proposée pour tronquer le domaine de calcul dans l'analyse par éléments finis des problèmes de couplage conductif et une formulation PEEC complète pour la modélisation de problèmes de couplage inductif est présentée. Des configurations tests de complexités croissantes sont considérées pour valider les approches précédentes. Ces travaux visent ainsi à apporter une contribution à la modélisation de cette classe de problèmes, qui tendent à devenir communs avec l'expansion des réseaux électriques. / Electromagnetic coupling phenomena between overhead power transmission lines and other nearby structures are inevitable, especially in densely populated areas. The undesired effects resulting from this proximity are manifold and range from the establishment of hazardous potentials to the outbreak of alternate current corrosion phenomena. The study of this class of problems is necessary for ensuring security in the vicinities of the interaction zone and also to preserve the integrity of equipment and devices there present. However, the complete modeling of this type of application requires the three-dimensional representation of the region of interest and needs specific numerical methods for field computation. In this work, the modeling of problems arising from the flow of electrical currents in the ground (the so-called conductive coupling) will be addressed with the finite element method. Those resulting from the time variation of the electromagnetic fields (the so-called inductive coupling) will be considered as well, and they will be treated with the generalized PEEC (Partial Element Equivalent Circuit) method. More specifically, a special boundary condition on the electric potential is proposed for truncating the computational domain in the finite element analysis of conductive coupling problems, and a complete PEEC formulation for modeling inductive coupling problems is presented. Test configurations of increasing complexities are considered for validating the foregoing approaches. These works aim to provide a contribution to the modeling of this class of problems, which tend to become common with the expansion of power grids. / Fenômenos de acoplamento eletromagnético entre linhas aéreas de transmissão de energia e outras estruturas vizinhas são inevitáveis, sobretudo emáreas densamente povoadas. Os efeitos indesejados decorrentes desta proximidadesão variados, indo desde o estabelecimento de potenciais perigosos até o surgimento de processos de corrosão por corrente alternada. O estudo desta classe de problemas é necessária para a garantia da segurança nas imediações da zona de interação e também para se preservar a integridade de equipamentos e dispositivos ali presentes. Entretanto, a modelagem completa deste tipo de aplicação requer a representação tridimensional da região de interesse e necessita de métodos numéricos de cálculo de campos específicos. Neste trabalho, serão abordadas as modelagens de problemas decorrentes da circulação de correntes elétricas no solo (ditos de acoplamentocondutivo) com o método dos elementos finitos. Também serão considerados problemas produzidos pela variação temporal dos campos eletromagnéticos (ditos de acoplamento indutivo), que serão tratados com o método PEEC(Partial Element Equivalent Circuit) generalizado. Mais especificamente, uma condição de contorno particular sobre o potencial elétrico é proposta para o truncamento do domínio de cálculo na análise de problemas de acoplamento condutivo com o método dos elementos finitos, e uma formulação completa tipo PEEC para a modelagem de problemas de acoplamento indutivo é apresentada. Problemas teste de complexidades crescentes são considerados para a validação das abordagens precedentes. Estes trabalhos visam fornecer desta forma uma contribuição à modelagem desta classe de problemas, que tendem a se tornar comuns com a expansão das redes elétricas.

Le but de ce travail est d'étudier l'influence de la microstructure de roches poreuses hétérogènes sur le comportement à l'échelle macroscopique. Ainsi, nous avons caractérisé la microstructure et les propriétés micromécaniques (grâce à des tests de nano-indentation) de deux roches oolithiques poreuses (calcaire de Lavoux et minerai de fer) pour calculer leurs propriétés mécaniques et thermiques effectives. Les roches oolithiques sont constituées d'un assemblage de grains poreux (oolithes), de pores et de cristaux intergranulaires. La microscopie électronique à balayage et la tomographie 3D aux rayons X ont été utilisées pour identifier les différents composants de ces roches. Une attention particulière a été accordée à la tomographie aux rayons X car cette méthode analytique permet de caractériser le réseau poreux (taille, distribution spatiale et fraction volumique), ainsi que la forme des oolithes et des cristaux inter-oolithiques. La nouveauté de ce travail réside dans la prise en compte de la forme 3D réelle des pores. Par conséquent, nous avons approximé les oolites poreuses par des sphères et les pores de forme irrégulière par des ellipsoïdes. Cette approximation a été réalisée grâce à l'analyse en composantes principales (ACP), qui fournit les propriétés géométriques telles que la longueur des demi-axes et l'orientation des ellipsoïdes résultants. La sphéricité des oolites approximées a été calculée et les valeurs proches de 1 nous ont permis de considérer les oolithes comme des sphères. Pour vérifier l'approximation dans le cas des pores, nous avons évalué la contribution de ces pores tridimensionnels de forme irrégulière aux propriétés élastiques et thermiques effectives. Ainsi, les tenseurs de contribution de souplesse pour les pores irréguliers 3D et leurs approximations ellipsoïdales ont été calculés en utilisant la méthode des éléments finis (FEM). Ces tenseurs ont été comparés et une erreur relative a été estimée pour évaluer la précision de l'approximation. Cette erreur produit une distance maximale de 4,5% entre les deux solutions pour les pores et les ellipsoïdes, ce qui vérifie la procédure d'approximation proposée basée sur ACP. La méthode numérique FEM a été vérifiée en comparant la solution numérique des tenseurs de contribution des ellipsoïdes à la solution analytique basée sur la théorie d'Eshelby. La différence entre ces deux solutions ne dépasse pas 3%. La même méthode numérique a été utilisée pour calculer les tenseurs de contribution de résistivité thermique. Les tenseurs de souplesse et de résistivité calculés ont été utilisés pour évaluer les propriétés élastiques effectives (module élastique et coefficient de cisaillement) et la conductivité thermique effective en considérant le schéma d'homogénéisation de Maxwell en deux étapes. Les résultats ont montré une influence importante de la porosité sur les propriétés effectives. Enfin, les résultats obtenus pour les pores irréguliers ont été comparés à ceux des ellipsoïdes et ils ont montré un bon accord avec un écart maximal de 4% ce qui vérifie l'approximation des pores de forme irrégulière par des ellipsoïdes triaxiaux / The aim of this work is to study the influence of the microstructure of heterogeneous porous rocks on the behavior at the macroscopic scale. Thus, we characterized the microstructure and micromechanical properties (thanks to nano-indentation tests) of two porous oolitic rocks (Lavoux limestone and iron ore) to calculate their effective mechanical and thermal properties. Oolitic rocks are constituted by an assemblage of porous grains (oolites), pores and inter-granular crystals. Scanning electron microscopy and X-ray 3D Computed Tomography were used to identify the different components of these rocks. Particular attention was given to X-Ray computed tomography since this analytical method allows the characterization of the porous network (size, spatial distribution, and volume fraction), and the shapes of oolites and inter-oolitic crystals. The novelty of this work lies in taking into account the 3D real shape of pores. Hence, we approximated porous oolites by spheres and irregularly shaped pores by ellipsoids. This approximation was performed thanks to the principal component analysis (PCA), which provides the geometrical properties such as length of semi-axes and orientation of resulting ellipsoids. The sphericity of the approximated oolites was calculated and the values close to 1 allowed us to consider oolites as spheres. To verify the approximation in the case of pores, we evaluated the contribution of these irregularly shaped three-dimensional pores to the overall elastic properties. Thus, compliance contribution tensors for 3D irregular pores and their ellipsoidal approximations were calculated using the finite element method (FEM). These tensors were compared and a relative error was estimated to evaluate the accuracy of the approximation. This error produces a maximum discrepancy of 4.5% between the two solutions for pores and ellipsoids which verifies the proposed approximation procedure based on PCA. The FEM numerical method was verified by comparing the numerical solution for compliance contribution tensors of ellipsoids to the analytical solution based on Eshelby’s theory. The difference between these two solutions does not exceed 3%. The same numerical method was used to calculate thermal resistivity contribution tensors. Calculated compliance and resistivity contribution tensors were used to evaluate effective elastic properties (bulk modulus and shear coefficient) and effective thermal conductivity by considering the two-step Maxwell homogenization scheme. The results showed an important influence of the porosity on effective properties. Finally, the results obtained for irregular pores were compared to those for ellipsoidal ones and they showed a good agreement with a maximum deviation of 4% which verifies once again the approximation of irregularly shaped pores by tri-axial ellipsoids

Beaucoup de méthodes numériques ont été développées pour résoudre l'équation de Vlasov, car obtenir des simulations numériques précises en un temps raisonnable pour cette équation est un véritable défi. Cette équation décrit en effet l'évolution de la fonction de distribution de particules (électrons/ions) qui dépend de 3 variables d'espace, 3 variables de vitesse et du temps. L'idée principale de cette thèse est de réécrire l'équation de Vlasov sous forme d'un système hyperbolique par semi-discrétisation en vitesse. Cette semi-discrétisation est effectuée par méthode d'éléments finis. Le modèle ainsi obtenu est appelé équation de Vlasov réduite. Nous proposons différentes méthodes numériques pour résoudre efficacement ce modèle: méthodes des volumes finis, méthodes semi-Lagrangiennes et méthodes Galerkin discontinus. / Many numerical methods have been developed in order to selve the Vlasov equation, because computing precise simulations in a reasonable time is a real challenge. This equation describes the time evolution of the distribution function of charged particles (electrons/ions), which depends on 3 variables in space, 3 in velocity and time. The main idea of this thesis is to rewrite the Vlasov equation in the form of a hyperbolic system using a semi-discretization of the velocity. This semi-discretization is achieved using the finite element method. The resulting model is called the reduced Vlasov equation. We propose different numerical methods to salve this new model efficiently: finite volume methods, semi-Lagrangian methods and discontinuous Galerkin methods.

Cette étude porte sur la détermination d’un critère de rupture des gaines de combustible de Zircaloy-4 détendu hydruré contenant un blister d’hydrures, en conditions accidentelles représentatives d’un accident d’injection de réactivité. Deux plages de comportement différentes en fonction de la température ont clairement été mises en évidence grâce à l’étude bibliographique, aux différentes campagnes d’essais mécaniques et aux analyses des faciès de rupture des éprouvettes rompues : une rupture de type fragile pour la gaine à 25°C et une rupture ductile à 350°C.A 25°C, la rupture fragile a été traitée par une analyse globale en mécanique élasto-plastique de la rupture. A partir des essais mécaniques effectués à 25°C sur les gaines contenant des blisters, des simulations numériques par éléments finis ont été réalisées avec le code CAST3M. Des calculs d’intégrales-J en pointe de fissure ont alors permis d’identifier un critère de rupture en ténacité moyenne de 13,8 +/- 3,1 MPa.m1/2.A 350°C, une campagne d’essais biaxés de type pression interne couplée à la traction axiale a été réalisée sur des tronçons de Zircaloy-4 contenant des blisters, à des biaxialités des contraintes représentatives du RIA. Il a été montré que la rupture de la gaine, avec et sans blister, avait lieu de façon ductile, que la déformation diamétrale à rupture de la gaine diminuait lorsque la profondeur de blister augmentait, et que la biaxialité des contraintes n’avait pas d’effet sur la rupture des gaines contenant un blister suffisamment profond.Un modèle d’endommagement ductile couplé à la plasticité, basé sur un formalisme de type GTN, a été utilisé. Afin d’améliorer la description de l’endommagement des gaines de Zircaloy-4, une nouvelle source de germination de porosités liée au paramètre de Lode a été intégrée dans le modèle. L’évaluation de la triaxialité des contraintes et du paramètre de Lode dans les simulations numériques de la rupture ductile des gaines à 350°C a notamment permis de comprendre certaines tendances expérimentales. / This study deals with the determination of a fracture criterion for hydrided, cold worked and stress relieved Zircaloy-4 fuel cladding tubes with hydride blister, during a reactivity initiated accident. Two types of fracture profiles were identified, depending on the temperature, thanks to a bibliographical study, mechanical tests and fracture profiles analysis : brittle fracture at 25°C, and ductile fracture at 350°C.At 25°C, brittle fracture was studied by a global analysis in elasto-plastic fracture mechanic. Numerical simulations were performed by a finite element method with the CAST3M code, based on mechanical tests on fuel cladding tubes with blisters. Crack tip J-integral calculations were carried out to identify a mean fracture toughness of 13,8 +/- 3,1 MPa.m1/2.At 350°C, internal pressure combined to axial tensile tests were performed on Zircaloy-4 fuel cladding tubes with hydride blisters, at stress biaxialities corresponding to those of a RIA. It was observed a ductile fracture for tubes with and without blister. It was shown that hoop strain at failure decreases when blister thickness increases, and that stress biaxiality has no effect on cladding tubes bearing a thick blister. A ductile fracture model based on the GTN model was employed and a nucleation of voids due to shear stress was introduced, based on the Lode parameter. Stress triaxiality and Lode parameter were assessed in numerical simulations to understand some experimental observations.

La simulation des procédés de mise en forme des composites à renforts tissés de type RTM est un enjeu majeur pour les industries de pointe mettant en œuvre ce type de matériaux. Au cours de ces procédés, la préforme tissée est souvent soumise à des déformations importantes. La connaissance et la simulation du comportement mécanique de la préforme à l’échelle macroscopique et à l’échelle mésoscopique s’avère souvent nécessaire pour optimiser la phase de conception de pièces composites formées par de tels procédés. Une analyse du comportement mésoscopique des préformes tissées de composites est d’abord proposée. Une loi de comportement hyperélastique isotrope transverse est développée, permettant de décrire le comportement mécanique de chacun des modes de déformation de la mèche : élongation dans la direction des fibres, compaction et distorsion dans le plan d’isotropie de la mèche, cisaillement le long des fibres. Une méthodologie est proposée pour identifier les paramètres de cette loi de comportement à l’aide d’essais sur la mèche et sur le tissu, et une validation par comparaison avec des essais expérimentaux est présentée. Une analyse du comportement macroscopique des renforts interlocks est ensuite proposée : une loi de comportement hyperélastique orthotrope est développée et implémentée. Cette loi, extension de la loi de comportement pour la mèche, est également basée sur une description phénoménologique des modes de déformation de la préforme. Une méthode d’identification des paramètres de cette loi de comportement est mise en œuvre, utilisant des essais expérimentaux classiques dans le contexte des renforts tissés (tension uniaxiale, compression, bias extension test, flexion). Cette seconde loi de comportement est validée par comparaison avec des essais de flexion et d’emboutissage hémisphérique. / Simulating the preforming stage of RTM-like fabric-reinforced composites manufactoring processes is a major stake for industries which use such materials. During such processes, the woven preform often undergoes finite deformations. Simulation methods are then required to optimize the conception of composites parts formed by RTM. An analysis of the mechanical behaviour of woven preforms at mesoscale is first presented. A transversely isotropic hyperelastic behaviour law is developped in order to describe the mechanical behaviour of each deformation mode of the yarn : elongation in the direction of fibres, compaction and distorsion in the transverse plane and along-fibres shear. An identification method is set up for this behaviour law which allows to compute its parameters by use of simple experimental tests on the yarn and on the fabric. The behaviour law is then validated by comparizon between simulations ans experimental tests. An analysis of the mechanical behaviour of interlock woven preforms at macroscale is the presented. An orthotropic hyperelastic behaviour law is developped and implemented as an extension of the behaviour law for the yarn. A phenomenological approach is also used to describe the mechanical behaviour of each deformation mode of the preform. An identification method is set up and put into place, based on tests well known in the field of fabric reinforcements : tensile test, crushing test, bias extension test, flexure test. A hemispherical stamping simulation is set up and compared to experiment for validation purpose.