Nous présentons deux approches différentes pour décrire le couplage hydromécanique des géomatériaux. Dans une approche de type phénoménologique nous traitons le milieu poreux comme un milieu continu équivalent dont les interactions entre la phase fluide et le squelette solide constituent le couplage du mélange à l'échelle macroscopique. En caractérisant le comportement de chaque phase nous arrivons à décrire le comportement couplé du milieu couplé saturé.Nous utilisons cette approche pour modéliser des essais expérimentaux faits sur un cylindre creux pour une roche argileuse (argile de Boom). Les résultats expérimentaux montrent de façon claire que le comportement de cette roche est fortement anisotrope. Nous avons choisi de modéliser ces essais en utilisant une lois de comportement élasto-plastique pour laquelle la partie élastique est transversalement isotrope.Le problème aux conditions aux limites étudié met en évidence des déformations localisées autour du forage intérieur. Afin de décrire de façon objective le développement de ces bandes de cisaillement nous avons considéré un milieu continu local de type second gradient qui permet d'introduire une longueur interne. De ce fait nous avons pu étudier le problème d'unicité en montrant qu'un changement de la discrétisation temporelle du problème aux limites peut conduire à des solutions différentes.Dans la deuxième approche étudiée nous caractérisons la microstructure du matériau avec des grains et un réseau de canaux pour la phase fluide. À l'aide d'un processus numérique d'homogénéisation nous arrivons à calculer numériquement la contrainte du mélange et le flux massique. Cette méthode d'homogénéisation numérique a été implémentée dans un code aux éléments finis afin d'obtenir des résultats macro. Une validation de l'implentation est proposée pour des calculs en mecanique pure et en hydromécanique. / We present two different approaches to describe the hydromechanical behaviour of geomaterials. In the first approach the porous media is studied through an equivalent continuum media where the interaction between the fluide and solid phases caracterize the coupling behaviour at the macroscale.We take into account this approach to model experimental tests performed over a hollow cylinder sample of clay rock (Boom Clay), considered for nuclear waste storage. The experimental results clearly show that the mechanical behaviour of the material is strongly anisotropic. For this reason we chose an elasto-plastic model based on Drucker-Prager criterion where the elastic part is characterized by cross anisotropy.The numerical results of boundary value problem clearly show localised strains around the inner hollow section. In order to regularize the numerical problem we consider a second gradient local continuum media with an enriched kinematic where an internal lenght can be introduced making the results mesh independent. The uniqueness study is carried out showing that changing the temporal discretization of the problem leads to different solutions.In the second approach we study the hydromechanical behaviour of a porous media that it is characterised by the microstructure of the material. The microstructure taken into account is composed by elastic grains, cohesives interfaces and a network of fluid channels. Using a periodic media a numerical homogenization (square finite element method) is considered to compute mass flux, stress and density of the mixture. In this way a pure numerical constitutive law is built from the microstructure of the media. This method has been implemented into a finite element code (Lagamine, Université de Liège) to obtain results at the macroscale. A validation of this implementation is performed for a pure mechanical boundary value problem and for a hydromechanical one.

Le phénomène de radoucissement d'une roche est étudié à partir de l'observation des modes de localisation des déformations. A partir d'essais de compression triaxiale sur un grès de Fontainebleau, on caractérise à la fois de comportement pré-localisation et la formation des bandes de cisaillement. Les conditions d'essais sont optimisées par un choix adapté de la géométrie de l'échantillon (élancement, taille, forme) et des conditions à l'interface roche-plateau. On montre que le frettage à l'interface modifie l'allure des courbes contraintes-déformations et affecte le mode de rupture. Les caractéristiques de bandes de cisaillement (orientation et épaisseur) sont également modifiées par le frettage. Une technique anti-frettage par lubrification de l'interface est utilisée et validée. L'effet de la présence d'un fluide interstitiel est étudié. On montre que la loi de contrainte effective de Terzaghi est vérifiée dans le cas des essais réalisés en conditions drainées. Les phases de fissuration des échantillons correspondant au seuil de microfissuration, au seuil de dilatance et au seuil de localisation sont identifiées. Le problème du contrôle expérimental de la phase post-pic est également discuté. Les observations de bandes de cisaillement montrent que l'épaisseur et l'angle d'orientation de la bande par rapport à la contrainte principale mineure ont tendance à diminuer quand on augmente la pression de confinement. La microstructure des bandes est étudiée à partir de microphotographies obtenues au microscope électronique à balayage (MEB) et à la loupe binoculaire. L'évolution de la porosité et de la granulométrie dans la zone de localisation est quantifiée. Sous faible pression de confinement, on observe un phénomène de dilatance à l'intérieur et à proximité de la bande caractérisé par une augmentation de la porosité. A forte pression de confinement, on observe un phénomène de dilatance à proximité de la bande et de contractance au co eur de la bande caractérisé par une intense fragmentation des grains et une réduction importante de la porosité. Un modèle de comportement élastoplastique du second-gradient avec écrouissage et radoucissement pour le grès de Fontainebleau est développé. La phase d'écrouissage est étalonnée sur la phase pré-localisation des essais triaxiaux. Le comportement radoucissant est évalué de façon indirecte à partir d'une analyse de l'épaisseur des bandes de cisaillement. (Résumé de l'auteur)

Nous présentons deux approches différentes pour décrire le couplage hydromécanique des géomatériaux. Dans une approche de type phénoménologique nous traitons le milieu poreux comme un milieu continu équivalent dont les interactions entre la phase fluide et le squelette solide constituent le couplage du mélange à l'échelle macroscopique. En caractérisant le comportement de chaque phase nous arrivons à décrire le comportement couplé du milieu couplé saturé.Nous utilisons cette approche pour modéliser des essais expérimentaux faits sur un cylindre creux pour une roche argileuse (argile de Boom). Les résultats expérimentaux montrent de façon claire que le comportement de cette roche est fortement anisotrope. Nous avons choisi de modéliser ces essais en utilisant une lois de comportement élasto-plastique pour laquelle la partie élastique est transversalement isotrope.Le problème aux conditions aux limites étudié met en évidence des déformations localisées autour du forage intérieur. Afin de décrire de façon objective le développement de ces bandes de cisaillement nous avons considéré un milieu continu local de type second gradient qui permet d'introduire une longueur interne. De ce fait nous avons pu étudier le problème d'unicité en montrant qu'un changement de la discrétisation temporelle du problème aux limites peut conduire à des solutions différentes.Dans la deuxième approche étudiée nous caractérisons la microstructure du matériau avec des grains et un réseau de canaux pour la phase fluide. À l'aide d'un processus numérique d'homogénéisation nous arrivons à calculer numériquement la contrainte du mélange et le flux massique. Cette méthode d'homogénéisation numérique a été implémentée dans un code aux éléments finis afin d'obtenir des résultats macro. Une validation de l'implentation est proposée pour des calculs en mecanique pure et en hydromécanique.

Trois stratégies de modélisation de structures sont proposées dans ce document. La première consiste à l'utilisation des éléments poutres multifibres de cinématique Timoshenko et des lois de comportement basées sur la mécanique de l'endommagement pour le béton et la plasticité pour les aciers. Après une présentation succincte de plusieurs éléments poutres multifibres existants dans la littérature, un nouvel élément Timoshenko avec des fonctions d'interpolation d'ordre supérieure est proposé. La deuxième stratégie concerne la prise en compte de l'Interaction du Sol avec la Structure (ISS) et le développement d'un élément d'interface 3D modélisant une fondation rigide superficielle reposant sur un massif de sol. Basé sur le concept des macro-éléments, ce nouvel outil est formulé en variables globales (forces et déplacements) et selon la théorie de plasticité, permettant ainsi de diminuer le temps de calcul. Enfin, la troisième stratégie traite de problème de la localisation des déformations et de l'objectivité des résultats des calculs dans un milieu enrichi (milieu avec microstructure). L'évolution de la zone de localisation est étudiée avec des modèles second gradient locaux pour différents types de lois simples basées sur la théorie de plasticité ou de l'endommagement. Des comparaisons avec des résultats expérimentaux et/ou des exemples numériques prouvent la performance des stratégies de modélisation proposées.

Au cours des dernières décennies, le développement de nouveaux matériaux a progressé pour les applications liées à la mécanique. De nouvelles générations de composites ont été développées, qui peut offrir des avantages par rapport aux tapis unidirectionnels renforcés de fibres couramment utilisés les matériaux prennent alors le nom de woven fabrics. Le comportement de ce matériau est fortement influencé par la micro-structure du matériau. Dans la thèse, les modèles mécaniques et les schémas numériques capables de modéliser les comportement des tissus et des matériaux de réseau généraux ont été développés. Le modèle prend en compte la micro-structure au moyen d'une technique d'homogénéisation. Les fibres dans le réseau ont été traités comme des micro-poutres, ayant une rigidité à la fois en extension et en flexion, avec différents types de connexions. La procédure développée a été appliquée pour obtenir les modèles mécaniques homogénéisés pour certains types de réseaux de fibres biaxiaux et quadriaxiaux, simulant soit des réseaux de fibres (en ce cas a été supposé parmi les fibres) ou des tissus avec une interaction négligeable entre les faisceaux de fibres et en empêchant tout glissement relatif (dans ce cas, les connexions ont été simulés au moyen de pivots). Différentes géométries ont été analysées, y compris la cas dans lesquels les fibres ne sont pas orthogonales. On obtient généralement un premier milieu à gradient mais, dans certains cas, la procédure d'homogénéisation lui-même indique qu'un continuum d'ordre supérieur est mieux adapté pour représenter la déformation de la micro-structure. Des résultats spéciaux ont été obtenus dans le cas de fibres reliées par pivots. Dans ce cas, un matériau orthotrope à module de cisaillement nul a été obtenu. Un tel matériau a un tenseur constitutif elliptique, il peut donc conduire à des concentrations de contrainte. Cependant, il a été montré que certaines considérations sur le comportement physique de tels réseaux indiqué que les termes d'ordre supérieur inclus dans l'expansion des forces internes et des déformations, de sorte qu'un matériau de gradient de déformation a été obtenu. Les résultats obtenus peuvent être utilisés pour la conception de matériaux spécifiques nécessitant des propriétés. Bien que le modèle de référence soit un matériau de réseau, les résultats obtenus peuvent être appliqué à d'autres types similaires de microstructures, comme des matériaux pantographiques, des micro-dispositifs composé de micro-poutres, etc. Ils étaient limités à la gamme d'élasticité linéaire, qui est petite déformation et comportement élastique linéaire. Ensuite, les simulations numériques ont été axées sur les tests d'extension et les tests de biais. Le obtenu configurations déformées sont conformes aux tests expérimentaux de la littérature, tant pour tissus équilibrés et non équilibrés. De plus, une comparaison entre les premier et deuxième gradients des prédictions numériques ont été effectuées. Il a été observé que les prédictions de deuxième gradient mieux simuler les preuves expérimentales. / In the past decades there has been an impressive progress in the development of new materials for mechanical related applications. New generations of composites have been developed, that can offer advantages over the unidirectional fibre-reinforced mats commonly used then materials take the name of woven fabrics. The behaviour of this material is strongly influenced by the micro-structure of the material. In the thesis mechanical models and a numerical scheme able to model the mechanical behaviour of woven fabrics and general network materials have been developed. The model takes in to account the micro-structure by means of a homogenization technique. The fibres in the network have been treated like microbeams, having both extensional and bending stiffness, with different types of connection, according to the pattern and detail of the network. The developed procedure was applied for obtaining the homogenized mechanical models for some types of biaxial and quadriaxial networks of fibres, simulating either fibre nets (in this case rigid connection were assumed among the fibres) or tissues with negligible interaction between the fibre bundles, and with relative sliding prevented (in this case the connections were simulated by means of pivots). Different geometries were analysed, including the cases in which the fibres are not orthogonal. A first gradient medium is usually obtained but, in some cases, the homogenization procedure itself indicates that a higher order continuum is better fit to represent the deformation of the micro-structure. Special results were obtained for the case of fibres connected by pivots. In this cases an orthotropic material with zero shear modulus was obtained. Such a material has a not elliptic constitutive tensor, thus it can lead to strain concentrations. However, it was shown that some considerations about the physical behaviour of such networks indicated that higher order terms had to be included in the expansion of the internal forces and deformations, so that a strain gradient material was obtained. The results obtained can be used for the design of specific materials requiring ad-hoc properties. Although the reference model is a network material, the results obtained can be applied to other similar kinds of microstructures, like pantographic materials, micro devices composed by microbeams etc. They have been limited at the range of linear elasticity, that is small deformation and linear elastic behaviour. Then, numerical simulations were focused on extension tests and bias tests. The obtained deformed configurations are consistent with the literature experimental tests, both for balanced and unbalanced tissues. Moreover, a comparison between first and second gradient numerical predictions was performed. It was observed that second gradient predictions better simulate the experimental evidences.

L'objectif technique principal auquel répond cette thèse est la mise au point d'une méthode de régularisation, donnant des résultats objectifs par rapport au maillage, pour traiter les problèmes couplés hydromécaniques dans le cas des géomatériaux. La modélisation proposée s'inscrit dans le cadre des milieux à microstructure dilatants et s'inspire, du point de vue numérique, des formulations second gradient. Elle permet de prédire de façon robuste les comportements hydrauliques et mécaniques produits par la dégradation d'un milieu naturel. Le modèle ainsi obtenu, dit second gradient de dilatation, se distingue par le faible nombre de degrés de liberté ajouté dans la discrétisation éléments finis par rapport à celui des milieux continus. L'objectif numérique est de réduire les temps de calcul pour rendre les études à portée industrielle acceptables. On montre son efficacité en réalisant des simulations couplées hydromécaniques d'excavations souterraines. Enfin, on présente un algorithme de recherche de solutions multiples dans la direction des modes singuliers basé sur les principes de la théorie de la bifurcation pour traiter des non-linéarités dues à des comportements irréversibles de matériaux adoucissants. Le cadre de l'analyse de bifurcation proposée se limite au cas des opérateurs symétriques. On montre ainsi, sur des simulations d'essais biaxiaux homogènes et d'excavations souterraines en conditions drainées, que cet algorithme est un outil de calcul efficace et robuste pour détecter plusieurs solutions mais également pour franchir des instabilités numériques liées au mauvais conditionnement des matrices tangentes au voisinage des points singuliers ou à la présence de snap-back.

Modéliser la localisation de la déformation dans les géomatériaux de manière objective nécessite l'utilisation de méthodes capables de régulariser le problème aux limites en introduisant une longueur caractéristique. Dans le cadre de ce travail de thèse, nous avons choisi d'utiliser les milieux à microstructures de types second gradient. Une question, se pose alors quant à l'identification des paramètres constitutifs qui interviennent dans la formulation de ces milieux. L'objectif de cette thèse est de mettre en place une méthode d'identification d'une loi de comportement enrichie de type second gradient. Une partie de l'identification a été réalisée à partir d'essais homogènes, cependant l'identification complète nécessite de considérer des modes de déformation à forts gradients, comme cela est le cas en présence d'une localisation de la déformation. La procédure d'identification développée s'appuie sur des résultats expérimentaux d'essais mécaniques sur le grès de Vosges conduits à deux confinements différents (20 et 30 MPa), pour lesquels le mode de déformation des échantillons a été caractérisé à l'aide de mesures de champs cinématiques (corrélation d'image), y compris en régime localisé. Un certain nombre d'observables ont été extraits de ces essais, qui ont servi à établir la comparaison entre calculs numériques et observations expérimentales. Une étude de sensibilité exhaustive a été menée pour établir une matrice de sensibilité utile par la suite pour l'optimisation des observables. Afin de calculer cette matrice, deux études de sensibilité sont effectuées. Ces études consistent à évaluer l'influence de la variation de chaque paramètre constitutif sur les données sélectionnées. La première étude de sensibilité porte sur la partie homogène des essais, elle permet l'optimisation d'un certain nombre de paramètres qui jouent un rôle uniquement dans cette partie. La deuxième étude concerne le régime de déformation localisé. Celle-ci permet le calcul de la matrice de sensibilité. Grâce à cette matrice, nous avons réalisé une optimisation des observables. Ainsi, chaque observable pourra être optimisé indépendamment des autres. A l'issue de cette optimisation un jeu de paramètre est proposé. Il permet de reproduire de manière fiable, les essais expérimentaux.

Les études de faisabilité concernant le stockage géologique profond des déchets radioactifs ont conduit un intérêt accru concernant la modélisation géomécanique de la roche hte. En France, une roche hte potentielle est l'argilite du Callovo-Oxfordien du site de Meuse/Haute Marne. Etant donné que le principe de stockage géologique profond repose fortement sur la capacité de confinement de la formation hte, sa faible perméabilité est d'une importance clé. La perméabilité étant dépendante de la microstructure du matériau et de son évolution sous chargement, le comportement couplé hydro-mécanique de l'argilite est important. En effet, des modifications mécaniques sont induites par le creusement de la galerie d'entreposage, générant une zone endommagée (EDZ), pouvant conduire une modification de la perméabilité dans le voisinage de la galerie. Dans les matériaux microstructure complexe comme l'argilite du Callovo-Oxfordien, le comportement macroscopique trouve son origine dans l'interaction des constituants micro-mécaniques. En plus du couplage entre le comportement hydraulique et mécanique, un couplage entre les échelles micro et macro existe. Par le biais de l'élaboration d'un cadre d'homogénéisation du couplage hydro-mécanique, une approche de modélisation deuxéchelles est développée dans ce travail, dans laquelle la relation constitutive macroscopique découle directement du comportement à l'échelle microscopique. Un modèle existant du couplage hydro-mécanique, reposant sur l'identification de grains et d'espaces poreux intergranulaires à l'échelle micro est adopté comme point de départ. Ce modèle repose sur une homogénéisation numérique du comportement à la petite échelle afin d'obtenir à l'échelle macroscopique la réponse en contrainte et de transport du fluide interstitiel. Ce modèle est basé sur un VER périodique qui permet de déduire le comportement macroscopique local de l'argilite. En réponse, en un point d'intégration macro donné, à un incrément de la déformation et du gradient de pression, la réponse du VER permet d'exprimer l'incrément de contrainte et de flux associé, constituant de fait un équivalent numérique de la relation constitutive. Les problèmes aux conditions limites macro et micro sont traités simultanément par la méthode élément fini. Pour obtenir les opérateurs tangents consistants à l'échelle macro, la méthode d'homogénéisation par condensation statique des opérateurs tangeants micro est étendu au cas avec couplage hydro-mécanique. L'implémentation du modèle double échelle et la mise en uvre des développements théoriques d'homogénéisation ont été effectués dans le code élément fini Lagamine (Université de Liège). Pour la modélisation de la localisation de la déformation à l'échelle macro, qui, dans un formalisme de milieu continu classique, souffre de la dépendance au maillage, l'approche double-échelle a été utilisée dans un formalisme de milieu enrichi de type milieu de second gradient pour matériau poreux saturé. Les capacités du modèle homogénéisé numériquement, utilisé dans un cadre de milieu de second gradient, sont ensuite démontrées par des simulations d'essais dométriques et d'essais de compression biaxiaux. L'approche se confirme être un moyen puissant pour modéliser l'anisotropie initiale et induite du comportement mécanique et du comportement hydraulique. Pour la modélisation du comportement de l'argilite du Callovo-Oxfordien, des VER sont construits en tenant compte des travaux de caractérisation de la géométrie des inclusions microscopiques et des résultats expérimentaux d'essais macroscopiques.La loi de comportement homogénéisée numériquement ainsi calibrée est utilisée dans des simulations de creusement de galerie jusqu'à des niveaux d'endommagement générant une localisation de la déformation.Ces calculs montrent à la fois la pertinence et l'applicabilité du concept double échelle pour l'évaluation du comportement hydromécanique des EDZ dans un contexte du stockage des déchets radioactifs. / Feasibility studies for deep geological radioactive waste disposal facilities have led to an increased interest in the geomechanical modelling of its host rock. In France, a potential host rock is the Callovo-Oxfordian claystone. The low permeability of this material is of key importance, as the principal of deep geological disposal strongly relies on the sealing capacity of the host formation. The permeability being coupled to the mechanical material state, hydromechanical coupled behaviour of the claystone becomes important when mechanical alterations are induced by gallery excavation in the so-called excavation damaged zone (EDZ). In materials with microstructure such as the Callovo-Oxfordian claystone [Robinet et al., 2012], the macroscopic behaviour has its origin in the interaction of its mi- cromechanical constituents. In addition to the coupling between hydraulic and mech- anical behaviour, a coupling between the micro (material microstructure) and macro will be made. By means of the development of a framework of computational homo- genization for hydromechanical coupling, a doublescale modelling approach is formu- lated, for which the macroscale constitutive relations are derived from the microscale by homogenization. An existing model for the modelling of hydromechanical coupling based on the distinct definition of grains and intergranular pore space [Frey, 2010] is adopted and modified to enable the application of first order computational homogenization for obtaining macroscale stress and fluid transport responses. This model is used to constitute a periodic representative elementary volume (REV) that allows the rep- resentation of the local macroscopic behaviour of the claystone. As a response to deformation loading, the behaviour of the REV represents the numerical equivalent of a constitutive relation at the macroscale. For the required consistent tangent operators, the framework of computational homogenization by static condensation [Kouznetsova et al., 2001] is extended to hy- dromechanical coupling. The theoretical developments of this extension are imple- mented in the finite element code Lagamine (Li` ege) as an independent constitutive relation. For the modelling of localization of deformation, which in classical FE meth- ods suffers from the well-known mesh dependency, the doublescale approach of hy- dromechanical coupling is combined with a local second gradient model [Collin et al., 2006] to control the internal length scale of localized deformation. By accepting the periodic boundary conditions as a regularization of the microscale deformation, the use of the multiscale model in combination with the local second gradient model can be used for modelling localization phenomena in HM-coupled settings with material softening. The modelling capacities of the approach are demonstrated by means of simula- tions of oedometer tests and biaxial compression tests. The approach is demonstrated to be a powerful way to model anisotropy in the mechanical as well as the hydraulic behaviour of the material both in the initial material state and as an effect of hy- dromechanical alterations. For the application to the modelling of Callovo-Oxfordian claystone, microstructural REVs are calibrated to geometrical characteristics of the inclusion that form the microstructure under consideration and to macroscale ex- perimental results of the mechanical behaviour. The calibrated constitutive relation is used in the simulation of gallery excavation processes. These computations give a proof of concept of the doublescale assessment of the hydromechanical behaviour of the excavation damaged zones around galleries in the context of nuclear waste disposal.

Ce travail porte sur le développement de modèles micromécaniques pour le calcul de la réponse homogénéisée de matériaux architecturés, en particulier des matériaux se présentant sous forme de treillis répétitifs. Les matériaux architecturés et micro-architecturés couvrent un domaine très large de de propriétés mécaniques, selon la connectivité nodale, la disposition géométrique des éléments structuraux, leurs propriétés mécaniques, et l'existence d'une possible hiérarchie structurale. L'objectif principal de la thèse est la prise en compte des nonlinéarités géométriques résultant des évolutions importantes de la géométrie initiale du treillis, causée par une rigidité de flexion des éléments structuraux faible en regard de leur rigidité en extension. La méthode dite d'homogénéisation discrète est développée pour prendre en compte les non linéarités géométriques pour des treillis quais périodiques; des schémas incrémentaux sont construits qui reposent sur la résolution incrémentale et séquentielle des problèmes de localisation - homogénéisation posés sur une cellule de base identifiée, soumise à un chargement contrôlé en déformation. Le milieu continu effectif obtenu est en général un milieu micropolaire anisotrope, dont les propriétés effectives reflètent la disposition des éléments structuraux et leurs propriétés mécaniques. La réponse non affine des treillis conduit à des effets de taille qui sont pris en compte soit par un enrichissement de la cinématique par des variables de microrotation ou par la prise en compte des seconds gradients du déplacement. La construction de milieux effectifs du second gradient est faite dans un formalisme de petites perturbations. Il est montré que ces deux types de milieu effectif sont complémentaires en raison de l'analogie existant lors de la construction théorique des réponses homogénéisées, et par le fait qu'ils fournissent des longueurs internes en extension, flexion et torsion. Des applications à des structures tissées et des membranes biologiques décrites comme des réseaux de filaments quais-périodiques ont été faites. Les réponses homogénéisées obtenues sont validées par des comparaisons avec des simulations par éléments finis réalisées sur un volume élémentaire représentatif de la structure. Les schémas d'homogénéisation ont été implémentés dans un code de calcul dédié, alimenté par un fichier de données d'entrée de la géométrie du treillis et de ses propriétés mécaniques. Les modèles micromécaniques développés laissent envisager du fait de leur caractère prédictif la conception de nouveaux matériaux architecturés permettant d'élargir les frontières de l'espace 'matériaux-propriétés' / The present thesis deals with the development of micromechanical schemes for the computation of the homogenized response of architectured materials, focusing on periodical lattice materials. Architectured and micro-architectured materials cover a wide range of mechanical properties according to the nodal connectivity, geometrical arrangement of the structural elements, their moduli, and a possible structural hierarchy. The principal objective of the thesis is the consideration of geometrical nonlinearities accounting for the large changes of the initial lattice geometry, due to the small bending stiffness of the structural elements, in comparison to their tensile rigidity. The so-called discrete homogenization method is extended to the geometrically nonlinear setting for periodical lattices; incremental schemes are constructed based on a staggered localization-homogenization computation of the lattice response over a repetitive unit cell submitted to a controlled deformation loading. The obtained effective medium is a micropolar anisotropic continuum, the effective properties of which accounting for the geometrical arrangement of the structural elements within the lattice and their mechanical properties. The non affine response of the lattice leads to possible size effects which can be captured by an enrichment of the classical Cauchy continuum either by adding rotational degrees of freedom as for the micropolar effective continuum, or by considering second order gradients of the displacement field. Both strategies are followed in this work, the construction of second order grade continua by discrete homogenization being done in a small perturbations framework. We show that both strategies for the enrichment of the effective continuum are complementary due to the existing analogy in the construction of the micropolar and second order grade continua by homogenization. The combination of both schemes further delivers tension, bending and torsion internal lengths, which reflect the lattice topology and the mechanical properties of its structural elements. Applications to textiles and biological membranes described as quasi periodical networks of filaments are considered. The computed effective response is validated by comparison with FE simulations performed over a representative unit cell of the lattice. The homogenization schemes have been implemented in a dedicated code written in combined symbolic and numerical language, and using as an input the lattice geometry and microstructural mechanical properties. The developed predictive micromechanical schemes offer a design tool to conceive new architectured materials to expand the boundaries of the 'material-property' space