Dans le cadre de la lutte contre les nuisances sonores et vibratoires, cette thèse porte sur la modélisation numérique des structures amorties par dispositifs piézoélectriques shuntés. La première partie du travail concerne la modélisation par éléments finis de structures en vibrations avec des pastilles piézoélectriques shuntées. Dans un premier temps, une formulation éléments finis originale, qui utilise des variables électriques globales (différence de potentiel et charge dans chaque pastille piézoélectrique), est analysée et validée. Dans un second temps, différentes stratégies de réduction de modèle basées sur la méthode de projection modale sont proposées pour résoudre le problème électromécanique discrétisé par éléments finis à moindre coût. La convergence de ces modèles d’ordre réduits est ensuite analysée pour les cas de shunts résistif et résonant. La deuxième partie du travail est consacrée à l’optimisation du système électromécanique, dans le but de maximiser l’amortissement apporté par les dispositifs piézoélectriques shuntés. Pour cela, une procédure d’optimisation topologique, basée sur la méthode SIMP (Solid Isotropic Material with Penalization method), est développée pour déterminer les géométries et les emplacements optimaux des pastilles piézoélectriques. Cette procédure permet de maximiser le coefficient de couplage électromécanique modal entre les éléments piézoélectriques et la structure hôte, ceci de façon indépendante du choix des composants du circuit électrique. Les avantages de l’approche proposée sont mis en avant à travers un exemple de validation et un cas d'application industrielle. Enfin, la dernière partie du travail propose une approche numérique pour modéliser et optimiser la réduction du rayonnement acoustique de plaques minces dans le domaine des basses fréquences avec des éléments piézoélectriques shuntés. Cette approche est valable pour n’importe quelle plaque mince bafflée et non trouée, indépendamment des conditions aux limites. Un exemple d’application concernant l’atténuation du rayonnent acoustique d’une plaque avec renforts est présenté et analysé. / Passive structural vibration and noise reduction by means of shunted piezoelectric patches is addressed in this thesis. The first part of the work concerns the finite element modeling of shunted piezoelectric systems. Firstly, an original finite element formulation, with only a couple of electric variables per piezoelectric patch (the global charge/ voltage), is analyzed and validated. Secondly, several reduced order models based on a normal mode expansion are proposed to solve the electromechanical problem. The convergence of these reduced order models is then analyzed for a resistive and a resonant shunt circuits. In the second part of the work, the concept of topology optimization, based on the Solid Isotropic Material with Penalization method (SIMP), is employed to optimize, in terms of damping efficiency, the geometry of piezoelectric patches as well as their placement on the host elastic structure. The proposed optimization procedure consists of distributing the piezoelectric material in such a way as to maximize the modal electromechanical coupling factor of the mechanical vibration mode to which the shunt is tuned, independently of the choice of electric circuit components. Numerical examples validate and demonstrate the potential of the proposed approach for the design of piezoelectric shunt devices. Finally, the last part of the work concerns the numerical modeling of noise and vibration reduction of thin structures in the low frequency range by using shunted piezoelectric elements. An efficient approach that can be applied to any thin continuous plates in an infinite baffle, independently of the boundary conditions, is proposed. An application example of a thin plate with reinforcements is presented and analyzed.

Nous avons conçu un modèle pour étudier le rôle de la membrane dans les déformations cellulaires. Le cytoplasme est modélisé par un fluide selon les équations de Navier-Stokes supplémentées d'un nouveau terme pour rendre compte de l'activité de l'actine. La membrane est considérée comme une coque mince élastique. Les outils développés (utilisant la méthode des éléments finis, FEMLAB^{\copyright} et MATLAB^{\copyright}) permettent des simulations 2D. La simulation de déformations à tension constante ou en supposant la tension fonction décroissante de la densité de membrane montre que seul le second cas permet des déformations de grande amplitude avec une forte courbure. La définition de la tension de membrane comme une fonction décroissante donne une interprétation mécanique de l'ajout de membrane qui permet de rendre compte de déformations plus importantes encore. Les dynamiques proposées pour la courbure de repos permettent d'étudier le maintien des déformations en l'absence de poussée d'actine. Nous suggérons que les paramètres mécaniques de la membrane varient selon les types et les compartiments cellulaires. Résumé Notre travail montre que la dépendance de la tension de membrane vis à vis de la densité de membrane est cruciale pour obtenir de grandes déformations. De même, l'ajout de membrane est optimal s'il a lieu de part et d'autre de la poussée d'actine. Enfin, nous avons exploré l'hypothèse de l'adaptation de la courbure de repos comme mécanisme permettant le maintien de la déformation.

La Tomographie Optique Diffuse Résolue en Temps (TOD-RT) est une technique d'imagerie clinique et préclinique en pleine croissance. Elle fournit les cartes d'absorption et de diffusion optique des organes explorés, et les paramètres physiologiques associés. La Tomographie Optique Diffuse de Fluorescence Résolue en Temps (TODF-RT) est basée sur la détection de photons de fluorescence. Elle permet de déterminer les cartes de la concentration et du temps de vie de sondes fluorescentes et d'accéder à une imagerie métabolique et moléculaire très importante pour le diagnostic et le suivi thérapeutique, en particulier en cancérologie. L'objectif de cette thèse était de réaliser des images 3D de TOD/TODF-RT sur des rongeurs en utilisant une technologie optique résolue en temps, les acquisitions étant réalisées à l'aide de fibres optiques disposées autour de l'animal et sans contact avec sa surface. Le travail a été mené en quatre étapes : 1- mise en place d'un dispositif d'imagerie de la surface de l'animal et reconstruction de son contour 3D, 2- modélisation de l'approche sans contact pour la résolution du problème direct, 3- traitement des mesures prenant en compte la réponse impulsionnelle de l'appareil, 4- établissement d'une méthode de reconstruction des images basée sur une sélection de points judicieusement choisis sur les profils temporels. Ces travaux ont permis d'obtenir des images optiques 3D de bonne qualité en réduisant la diaphonie entre l'absorption et la diffusion. Ces améliorations ont été obtenues tout en diminuant le temps de calcul, par comparaison avec les méthodes utilisant la totalité des profils temporels.

Ce travail de thèse est consacré à l'étude des renforts textiles techniques 2D à l'échelle mésoscopique. La méthode des éléments finis est utilisée pour résoudre un problème aux limites, fortement non linéaire, dans le domaine du renfort fibreux. Deux nouveaux modes de déformations pour le comportement des mèches de renforts ont été développés et caractérisés. Le premier mode est la compression longitudinale, qui permet de traduire la faible rigidité des mèches lorsqu'elles subissent une dilatation négative dans le sens des fibres. La relation conflictuelle sur le plan numérique entre la rigidité en tension, très forte, et la rigidité en compression longitudinale, très faible, peut être résolue via trois méthodes : réduction du pas de temps critique, addition de la contribution en tension ou avec une nouvelle stratégie pour l'actualisation du champ de contrainte. Le second mode de déformation est la dilatation transversale des mèches considérée comme conséquence directe de la compression longitudinale. Ce phénomène d'expansion de matière dans les directions transverses peut être observé avec un essai de compression longitudinale in-situ sous tomographie X et est communément ignorer. Un protocole expérimental a été mise en place pour mesurer cette dilatation transversale des mèches et un coefficient de Poisson a pu être caractérisé par identification inverse. Une campagne expérimentale a permis d'identifier les paramètres matériau du modèle mésoscopique et les résultats de simulations sont comparés aux images issues d'essai mécanique in situ sous tomographe. / This thesis is devoted to the mesoscopic study on the performance of textile reinforcements. F.E. simulation is carried out on a mesoscale model for the fibrous material, based on which two kinds of new deformation modes are developed. The first one is a longitudinal compression mode, which is used to reflect the small stiffness when the yarn is compressed longitudinally. The incompatibility problem between the small longitudinal compression stiffness and the large tension stiffness are solved by three different strategies: constraining the critical step time, adding the nonlinear tension part, or using a new strategy to update the stress. The second one is transverse expansion mode that could reflect the influence from longitudinal deformation to transverse deformation. This deformation could be found in tomography view but was ignored by the former researches. An experiment is designed to measure the expansion magnitude, and the geometrical inverse fitting process is applied to measure the value of the longitudinal-transverse Poisson ratio. The parameters of the mesoscale model are measured by a series of mechanical experiments and the simulation results are verified by the tomography methodology.

Ce travail, à l'interface entre les mathématiques appliquées et la physique, s'articule autour de la modélisation numérique des vésicules biologiques, un modéle pour les globules rouges du sang. Pour cela, le modéle de Canham et Helfrich est adopté pour décrire le comportement des vésicules. La modélisation numérique utilise la méthode des fonctions de niveau dans un cadre éléments finis. Un nouvel algorithme de résolution numérique combinant une technique de multiplicateurs de Lagrange avec une adaptation automatique de maillages garantit la conservation exacte des volumes et des surfaces. Cet algorithme permet donc de dépasser une limitation cruciale actuelle de la méthode des fonctions de niveau, à savoir les pertes de masse couramment observées dans ce type de problémes. De plus, les propriétés de convergence de la méthode des fonctions de niveau se trouvent ainsi grandement améliorées, comme l'indiquent de nombreux tests numériques. Ces tests comprennent notamment des problémes d'advection élémentaires, des mouvements par courbure moyenne ainsi que des mouvements par diffusion de surface. Concernant l'équilibre statique des vésicules, une condition générale d'équilibre d'Euler-Lagrange est obtenue à l'aide d'outils de dérivation de forme. En dynamique, le mouvement d'une vésicule sous l'action d'un écoulement de cisaillement est étudié dans le cadre des nombres de Reynolds élevés. L'effet du confinement est considéré, et les régimes classiques de chenille de char et de basculement sont retrouvés. Finalement, pour la premiére fois, l'effet des termes inertiels est étudié et on montre qu'au delà d'une valeur critique du nombre de Reynolds, la vésicule passe d'un mouvement de basculement à un mouvement de chenille de char. / This work, at the interface between the Applied Mathematics and Physics is connected about the numerical modelisation of biological vesicles, a pattern for the red blood cells. For this reason, the pattern of Canham and Helfrich is adopted to describe the behaviour of the vesicles. The numerical modelisation uses the Level Set method in finite element framework. A new algorithm of numerical resolution combining one technique of Lagrange multipliers with an automatic mesh adaptation ensures the accurate conservation of volumes and surfaces. Thus this algorithm enables to exceed an existing crucial restriction of the Level Set method, that's to say, the wastes of mass usually noticed in this kind of problems. Moreover, the proprieties of convergence of the Level Set method are thus much more improved, as shown in many numerical tests. Those tests chiefly include elementary problems of advection, motions by mean curvature just as motions by spread of surface. Concerning the static equilibrum of the vesicles, a mechanical equilibrum equation (Euler-Lagrange equation) of a vesicle membrane under a generalized elastic bending energy is obtained and the approach is based on shape optimization tools. In dynamics, the motion of a vesicle under the effect of a shear flow is elaborated in the frames of reference of high Reynolds numbers. The effect of confinement is respected, and the standard regimes of tank-treading and of tumbling motion are found again. Finally, for the first time, the effect of the inertia terms is elaborated and we show that beyond a critical value of the number of Reynolds the vesicle passes from a tumbling motion to a tank-treading motion.

Le but de cette thèse est de développer les outils numériques et semi-analytiques qui nous permettent d’étudier des guides non-linéaires complexes et réalistes qui pourraient être fabriqués et caractérisés expérimentalement.Nous présentons une étude complète d’une version améliorée du guide d’onde plasmonique non-linéaire à une dimension, en ajoutant deux couches tampons de diélectrique linéaire entre le cœur non-linéaire isotrope et les deux gaines métalliques. Ses couches réduisent les pertes et permettent leur diminution avec la puissance contrairement aux guides simples. De plus, les principaux modes plasmoniques non-linéaires peuvent présenter une transition spatiale vers des modes de types différents, qui peut être contrôlée par la puissance.Par la suite, nous étudions un nouveau guide plasmonique non-linéaire à fente utilisant un métamatériau, soit dans le cœur non-linéaire soit dans les gaines linéaires. Nous avons mis au point une méthode semi-analytique et une méthode numérique afin d’étudier les solutions stationnaires non-linéaires dans ce nouveau guide non-linéaire anisotrope. Nous avons montré analytiquement et numériquement que le cœur non-linéaire anisotrope peut être conçu afin d’atteindre de forts effets non-linéaires à faible puissance. Pour certains métamatériaux dans les gaines linéaires, la figure de mérite de ce guide d’onde augmente de plus de 50 fois par rapport aux guides isotropes.Pour conclure, nous présentons une nouvelle méthode basée sur la méthode des éléments finis de type vectoriel couplée à l’algorithme à puissance fixée pouvant calculer rigoureusement les effets non-linéaires dans des guides d’onde plasmoniques 2D. / The main goal of this PhD is to develop the semi-analytical and the numerical tools that allow us to study complicated and realistic nonlinear plasmonic waveguides which can be fabricated and characterized experimentally.First, we present a full study of an improved version of the one-dimensional nonlinear plasmonic slot waveguide, by adding two linear dielectric buffer layers between the isotropic nonlinear core and the two metal claddings. These additional layers reduce the overall losses and allow the losses to decrease with the power for some configurations unlike the usual slot. Furthermore, the main plasmonic modes can exhibit nonlinear spatial modal transitions towards new families of modes that can be controlled with the power.Second, we propose and study new one-dimensional nonlinear plasmonic slot waveguides with metamaterial regions either in the nonlinear core or in the linear claddings. For the metamaterial nonlinear core, we developed semi-analytical and fully numerical methods in order to study the nonlinear stationary solutions propagating in this anisotropic nonlinear waveguide. We have demonstrated both analytically and numerically that the anisotropic nonlinear core can be designed in order to achieve strong nonlinear effects at low power.For the structures with metamaterial linear claddings, the figure of merit can be extremely enhanced by more than 50 times compared with the simple one. Finally, we present the full derivation of a new nonlinear full vectorial finite element method based on the fixed power algorithm in order to quantify rigorously the nonlinear characteristics of realistic two-dimensional nonlinear plasmonic structures.

Les structures maçonnées sont largement connues en génie civil comme constituant d’une partie des bâtiments, mais également en tant que garnissages réfractaires dans des structures utilisées à hautes températures, par exemple en sidérurgie. Malheureusement, les outils actuels ne sont pas suffisamment puissants pour prédire le comportement de ces structures avec l’apparition de fissures et pour tenir compte du comportement non linéaire d’un des deux constituants (le mortier par exemple). Ce travail de thèse contribue à la modélisation multi-échelles des maçonneries classiques et des garnissages réfractaires avec un faible coût numérique grâce à la technique d’homogénéisation périodique. Les techniques de modélisation et de simulation du comportement des maçonneries sont présentées et développées. L’influence des lois d’interface entre briques et mortier, des paramètres géométriques et matériels, ainsi que de la densité des fissures sur le comportement effectif des maçonneries est étudié. Trois approches (une extension analytique de Cecchi et Tralli, une approche numérique et un modèle micromécanique) sont proposées pour la détermination du comportement effectif d’une cellule périodique dans le cas de maçonneries avec mortier viscoélastique microfissuré et briques saines élastiques ou rigides. Les résultats des calculs sur deux exemples de maçonneries (1D et 2D) ont confirmé que l’approche multi-échelle est une solution appropriée avec une grande capacité à exprimer le comportement des maçonneries viscoélastiques microfissurées. Ce travail, limité au cas sans propagation de fissures, peut être étendu aux mortiers à comportement viscoplastique. / Masonry structures are widely used in civil engineering as part of buildings or in refractory linings of structures working at high temperatures, like in steel industry. Unfortunately, the present tools are not powerful enough to predict the behavior of these structures at their micro-cracked state and/or if one of their constituents behaves nonlinearly (e.g. the mortar). This research contributes to the multi-level modeling of classical masonries and refractory linings with a low numerical cost using basically the periodic homogenization technique. Modeling and simulation techniques of masonry behavior are presented and developped. The influence of interface law between bricks and mortar, of geometrical and material parameters, and of crack density on the effective masonry behavior is studied. Three approaches (analytical extension of Cecchi and Tralli, numerical approach and micromechanical modeling) were proposed to determine the effective behavior of a periodic masonry cell with micro-cracked viscoelastic mortar and safe elastic or rigid bricks. The results obtained on two examples of masonry (1D and 2D) confirmed that the multi-scale approach is a suitable solution with a great ability to model the effective behavior of microcracked viscoelastic masonry. This work, actually limited to the case without crack propagation, could be extended to mortars with viscoplastic behavior.

Dans le cadre de la lutte contre les nuisances sonores et vibratoires, cette thèse porte sur la modélisation numérique des structures amorties par dispositifs piézoélectriques shuntés. La première partie du travail concerne la modélisation par éléments finis de structures en vibrations avec des pastilles piézoélectriques shuntées. Dans un premier temps, une formulation éléments finis originale, qui utilise des variables électriques globales (différence de potentiel et charge dans chaque pastille piézoélectrique), est analysée et validée. Dans un second temps, différentes stratégies de réduction de modèle basées sur la méthode de projection modale sont proposées pour résoudre le problème électromécanique discrétisé par éléments finis à moindre coût. La convergence de ces modèles d’ordre réduits est ensuite analysée pour les cas de shunts résistif et résonant. La deuxième partie du travail est consacrée à l’optimisation du système électromécanique, dans le but de maximiser l’amortissement apporté par les dispositifs piézoélectriques shuntés. Pour cela, une procédure d’optimisation topologique, basée sur la méthode SIMP (Solid Isotropic Material with Penalization method), est développée pour déterminer les géométries et les emplacements optimaux des pastilles piézoélectriques. Cette procédure permet de maximiser le coefficient de couplage électromécanique modal entre les éléments piézoélectriques et la structure hôte, ceci de façon indépendante du choix des composants du circuit électrique. Les avantages de l’approche proposée sont mis en avant à travers un exemple de validation et un cas d'application industrielle. Enfin, la dernière partie du travail propose une approche numérique pour modéliser et optimiser la réduction du rayonnement acoustique de plaques minces dans le domaine des basses fréquences avec des éléments piézoélectriques shuntés. Cette approche est valable pour n’importe quelle plaque mince bafflée et non trouée, indépendamment des conditions aux limites. Un exemple d’application concernant l’atténuation du rayonnent acoustique d’une plaque avec renforts est présenté et analysé. / Passive structural vibration and noise reduction by means of shunted piezoelectric patches is addressed in this thesis. The first part of the work concerns the finite element modeling of shunted piezoelectric systems. Firstly, an original finite element formulation, with only a couple of electric variables per piezoelectric patch (the global charge/ voltage), is analyzed and validated. Secondly, several reduced order models based on a normal mode expansion are proposed to solve the electromechanical problem. The convergence of these reduced order models is then analyzed for a resistive and a resonant shunt circuits. In the second part of the work, the concept of topology optimization, based on the Solid Isotropic Material with Penalization method (SIMP), is employed to optimize, in terms of damping efficiency, the geometry of piezoelectric patches as well as their placement on the host elastic structure. The proposed optimization procedure consists of distributing the piezoelectric material in such a way as to maximize the modal electromechanical coupling factor of the mechanical vibration mode to which the shunt is tuned, independently of the choice of electric circuit components. Numerical examples validate and demonstrate the potential of the proposed approach for the design of piezoelectric shunt devices. Finally, the last part of the work concerns the numerical modeling of noise and vibration reduction of thin structures in the low frequency range by using shunted piezoelectric elements. An efficient approach that can be applied to any thin continuous plates in an infinite baffle, independently of the boundary conditions, is proposed. An application example of a thin plate with reinforcements is presented and analyzed.

Une méthode d’analyse des topographies résiduelles, obtenues par essais de nanoindentation Berkovich sur trois échantillons de nickel polycristallins cubiques à faces entrées (CFC), a été développée dans cette thèse. L’objectif de la méthode est d’évaluer la richesse de l’information contenue dans les empreintes pour l’identification de tout ou partie des paramètres d’une loi de plasticité monocristalline. La création d’une base de données de topographies résiduelles, mesurées par microscopie à force atomique (AFM), constitue la partie expérimentale du travail. Les distributions et dimensions de bourrelets montrent une grande sensibilité à l’orientation relative indenteur/grain et au taux d’écrouissage. Les topographies obtenues s’avèrent être de véritables « empreintes digitales » du mécanisme de plasticité à l’échelle du grain. L’élaboration sous le code ZéBuLon d’un modèle éléments finis (EF) 3D de l’essai de nanoindentation Berkovich, intégrant la loi de Méric-Cailletaud, permet de retrouver les observations expérimentales. Une étude numérique confirme la sensibilité de la topographie à l’orientation relative indenteur/grain et aux paramètres plastiques, notamment aux coefficients de la matrice d’interaction des dislocations présentes sur les systèmes de glissement. Afin d’évaluer la richesse du contenu informatif des empreintes, un indice d’identifiabilité est proposé. Son calcul est basé sur la multi-colinéarité des vecteurs de sensibilité des topographies résiduelles aux paramètres de la loi. Il permet de quantifier, en fonction des données topographiques prises en compte, le caractère mal posé du problème d’identification paramétrique. Les résultats obtenus montrent que l’identification de quatre à cinq paramètres de la loi de Méric-Cailletaud est envisageable en exploitant seulement deux empreintes. Ces travaux ouvrent la voie à l’identification du comportement à l’échelle du cristal, guidée par l’identifiabilité paramétrique. / In this thesis, a method is developed to analyse the residual topographies obtained by Berkovich nanoindentation tests on three face-centered cubic (FCC) polycrystalline nickel samples. The purpose is to measure the information richness of imprints for identifying all or part of parameters of a single crystal plasticity law. The experimental part consists in creating a residual topography database by atomic force microscopy (AFM) measurements on the samples. Pile-up distributions and sizes are largely sensitive to the indenter/grain relative orientation and the hardening rate. The topographies are true “fingerprints” of the plasticity mechanism at the grain scale. A 3D finite element (FE) modelling of nanoindentation test is developed using the code ZeBuLon and making use of the Méric-Cailletaud law. Numerical results show a good agreement with experimental observations and are largely sensitive to the indenter/grain relative orientation and the plastic parameters, including the interaction matrix coefficient specifying the interaction between dislocations on different slip systems. To measure the imprint information content, an identifiability index is proposed. Its calculation is based on the multicollinearity among the sensitivity vectors of topographies to the law parameters. According to the considered topographies, it measures if the numerical model updating problem is ill-posed. The results show that four to five parameters of the Méric-Cailletaud law can be identified by considering two topographies. This work paves the way for identifying the material behaviour at the grain scale using parametric identifiability

Le béton présente une rupture quasi-fragile du fait de la présence d’une zone d’élaboration de la rupture principalement microfissurée et de taille conséquence, se développant en avant du front de fissure. L’objectif de ces travaux consiste à mettre en évidence le processus d’évolution de la zone d’élaboration et l’étendue de cette dernière, et ce, sur la base d’une campagne d’essai Wedge Splitting suivie de simulations numériques aux éléments-finis utilisant le modèle d’endommagement isotrope de Fichant – La Borderie à l’échelle mésoscopique. Expérimentalement, les courbes de résistance obtenues dans le cadre de la Mécanique Linéaire Élastique de la Rupture équivalente attestent d’un phénomène de confinement précoce de la zone d’élaboration manifestement lié à la géométrie testée et aux propriétés du béton étudié.Les simulations numériques obtenues sur la base du modèle d’endommagement sont en accord avec les résultats expérimentaux et conduisent également à soupçonner ce phénomène de confinement. Par la suite, la simulation numérique est à nouveau utilisée sur cette même géométrie de spécimen mais en réduisant les propriétés de rupture de la matrice cimentaire afin de diminuer la taille de la zone d’élaboration. Les résultats de cette nouvelle simulation montrent un développement libre de la zone d’élaboration suivie d’une propagation auto-similaire de la fissure principale attestée par la présence d’un régime plateau de la courbe de résistance correspondante.Une analyse détaillée du champ d’endommagement puis du champ d’énergie restituée est réalisée et permet de mettre en évidence un critère de développement de la zone d’élaboration fondée sur une valeur maximale du taux local de restitution d’énergie d’endommagement. / Concrete exhibits a quasibrittle fracture due to the existence of a large fracture process zone (FPZ), mainly microcracked, which develops ahead of the crack front. The aim of the current work consists in highlighting the FPZ development and its extent. For that purpose,an experimental campaign using Wedge Splitting Test was carried out and followed by finite element simulation using Fichant – La Borderie isotropic damage model acting at the mesoscale of concrete. Experimental analysis exhibits that by the use from Resistance curves estimated in the framework of equivalent Linear Elastic Fracture Mechanics, the used geometry combined to the studied concrete properties are subjected to a confinement of FPZ. Numerical simulations achieved by the damage model are in agreement with experimental results and also seem to show FPZ confinement. There after, numerical simulations are used again on the same specimen geometry but by decreasing fracture properties of cementitious matrix, in order to minimize the FPZ size. Numerical results exhibit that the FPZ was now freely developed and followed by a self-similar propagation of the main crack as shown by the existence of a “plateau” value on the corresponding Resistance curve. A numerical analysis is performed and leads to the propositionof a FPZ development criteria based on a maximal value of the local damage energy release rate.