The classical continuum mechanics, which studies the mechanical behavior of structures based on partial differential equations, shows its deficiencies when it encounters a discontinuity. Peridynamics based on integral equations can simulate fracture but suffers from high computational costs. A hybrid local/nonlocal model combining the advantages of peridynamics with those of classical continuum mechanics can simulate fracture and reduce the computational cost. Under the framework of the hybrid local/nonlocal model, this research developed an approach and computational codes for fracture simulations. First, we developed the computational codes based on the hybrid model with a priori partition of the domain between local and nonlocal models to tackle engineering problems with relevant level of difficulty. Second, we developed a strength-induced approach to enhance the capability of the computational codes because the strength-induced approach can limit the peridynamic model to necessary computational steps at the time level and a relatively small zone at the space level during a simulation. The strength-induced approach also improved the hybrid models by enabling an automatic partition of the domain without manual involvement. At last, a strength-induced computational code was developed based on this research. This dissertation complemented and illustrated numerically some previous work of Cohmas laboratory, in which a new route was introduced toward simulating the whole process of material behaviors including elastic deformation, crack nucleation and propagation until structural failure.

No description available.

Ce travail concerne l'étude de l'intégrité des matériaux ou des structures assemblées en s'intéressant à la modélisation du mécanisme de rupture à l'aide des modèles de zones cohésives (MZC). Cette approche présente l'avantage d'incorporer une longueur caractéristique dans la description de la rupture, ce qui permet notamment d'évaluer des effets de taille. Trois paramètres caractérisent ces MZC : la contrainte de traction Tmax à laquelle l'interface ou le matériau peut résister avant d'amorcer sa décohésion, l'ouverture critique à partir de laquelle une fissure est créée localement et finalement une loi traction-ouverture qui décrit la répartition des efforts cohésifs selon le mécanisme opérant.L'objectif principal de cette thèse est d'identifier les paramètres cohésifs caractérisant la rupture interfaciale dans un joint de colle ou cohésif dans un matériau. Pour cela, une première étape était d'élaborer un modèle analytique, décrivant correctement la cinématique d'un essai DCB ou Wedge, pour caractériser la rupture mode I. Bien que le mode d'ouverture soit opérant, le champ de déplacement au voisinage de l'entaille n'est pas K-dominant pour ces essais, du moins pas toujours. Plusieurs lois de traction-ouverture ont été considérées afin d'étudier leur influence sur la réponse locale et globale de l'essai. Une méthodologie d'identification inverse a été proposée à partir d'un modèle analytique, qui consiste à extraire les paramètres cohésifs en minimisant l'erreur au sens des moindres carrés entre les déflections analytique et numérique. Une fois validée, elle a été par la suite appliquée à un cas réel, qui est l'essai d'insertion de lame. La mesure du champ de déplacement expérimental est possible grâce à une mesure du champ de déplacement par corrélation d'image.Une analyse approfondie a été aussi consacrée à l'étude de la "Process Zone" pour un essai DCB ou Wedge. Cette étude met en évidence la variation de Lcz en fonction de la géométrie des éprouvettes, des propriétés de la zone cohésive, des propriétés mécaniques du matériau ou encore la forme de la loi traction-ouverture utilisée. Une nouvelle expression pour estimer Lcz est établie pour les zones cohésives rectangulaires et triangulaires.Une deuxième approche d’identification locale, basée sur le travail de Réthoré et Estevez (2013), a été aussi proposée et discutée. Elle a été mise en oeuvre pour un essai d'insertion de lame, avant d'être appliquée à un essai de flexion 4 points entaillé. Un aller-retour entre simulation numérique et résultat expérimental permet d'identifier les propriétés cohésives du matériau ou de l'interface / This work concerns the study of materials and assemblies structures integrity using cohesive zone model (CZM) to analyze fracture. These models have the advantage to incorporate a characteristic length in the description of fracture initiation and propagation, which can lead to size effects studies. Three parameters characterize the CZM : the maximal cohesive traction Tmax to which the interface or the material can resist before the onset of debonding, the critical crack opening from which a crack is created locally and finally the traction-separation law which describe the cohesive traction distribution depending on fracture process.The main purpose of this thesis is to identify the cohesive zone parameters describing fracture at the interface or in the material. The first step was to elaborate an analytical model which can describe properly the DCB or Wedge Test kinematic, to characterize mode I fracture. Despite the fact of mode I fracture, the displacement field near the crack-tip is not K-dominant for these tests, at least not always. Various traction-separation laws were considered in order to study their influence on the local and global response of the test. An inverse identification methodology has been proposed from the analytical model, which can extract cohesive parameters through a least square error minimization between numerical and analytical deflection. Once validated, it was subsequently applied to a real Wedge test. The experimental displacement field measurement was done due to digital image correlation measurement.A deep analysis to evaluate the fracture process zone length has been also dedicated in the case of Wedge or DCB Tests. This analysis shows that Lcz is not an intrinsic interface or materials property and it can vary depending on the sample's geometry, the cohesive zone properties or the traction-separation law. A new expression to determine Lcz was established for rectangular and triangular cohesive zone.A second local identification approach, based on the work of Réhoré and Estevez (2013), has been also proposed. It was implemented to analyze the Wedge test, before applying it to a notched four points bending test. A round trip between numerical simulations and experimental results allow identifying the cohesive properties in the materiel or at the interface.

Cette thèse est consacrée à la simulation du comportement à la rupture de bois sous des chargements à long terme et sous des conditions d'Humidité Relatives (HR) de l'air variables. Il est connu que le bois est un matériau fortement hygroscopique, ses propriétés mécaniques et de rupture sont en effet très dépendantes de sa teneur en eau. En outre, la stabilité d'une fissure existante dans un élément structural peut être fortement influencée parles variations, en particulier brusques, d'humidité relative qui peut conduire à la rupture inattendue de l'élément.L'approche thermodynamique proposée intègre l'effet de mécanosorption dans l'expression analytique de la déformation, en découplant les déformations mécaniques et celles dues au comportement mécanosorptif du matériau. En outre, la rupture quasi-fragile du matériau boisest traduite par un modèle de zone cohésive dont les paramètres de cohésion sont fonctions de la teneur en eau afin de simuler l’effet de l'humidité sur les propriétés de rupture. Sur cette base, une formulation incrémentale permet l'intégration de l'effet des variations soudaines d’humidité relative (autrement dit, le choc hydrique) sur la zone d’élaboration(zone cohésive) en introduisant un champ de contraintes supplémentaires le long de cette zone. Fonction de la variation de HR, ce champ de contraintes supplémentaires dépend de l'état de contrainte et de l'ouverture de la fissure le long de la zone cohésive, mais également de l'humidité en pointe de fissure (matériau non endommagé). Dans l'analyse par éléments finis, un opérateur tangent algorithmique est utilisé pour résoudre le problème non linéaire en combinant le modèle de mécanosorption et le modèle de zone cohésive et en intégrant l'effet du choc hydrique.La simulation du comportement d'une éprouvette entaillée soumise à un chargement constant et à des variations cycliques de HR montre un fort couplage entre le comportement mécanosorptif et l'effet du choc hydrique HR sur la zone d’élaboration. Ce couplage entraîne une augmentation de la propagation des fissures et conduit à une fissuration plus précoce par rapport à celle obtenue à partir du modèle de mécanosorption seul ou à partir du modèle de zone cohésive en intégrant l'effet des variations soudaines de HR. En outre, le couplage entre le modèle mécanosorptif et le modèle de zone cohésive en intégrant l'effet du chochydrique montre l'intérêt d'une telle approche numérique pour décrire le comportement complexe des éléments de charpente en bois soumis à des conditions climatiques variables,comportement qui ne peut être prédit par une simple superposition des deux modélisations. / This thesis is dedicated to the simulation of the fracture behavior of wood under long-termloading and variable relative humidity conditions. Indeed, wood is well-known to be a highlyhygroscopic material in so far as its mechanical and fracture properties are very dependenton moisture. Moreover, the stability of an existent crack in a structural element can bestrongly affected by the sudden variations of relative humidity (RH) and can lead tounexpected failure of the element.The thermodynamic approach proposed in this thesis includes the mechano-sorptive effect inthe analytical expression of the deformation, by operating a decoupling of the strain in amechanical part and a mechano-sorptive part in material. Moreover, the quasi-brittle fractureof wood is here simulated from a cohesive zone model whose cohesive parameters arefunctions of the moisture in order to mimic the moisture-dependent character of the fractureproperties. On this basis, an increment formulation allows the integration of the effect ofsudden RH variations on the fracture process zone (cohesive zone) by introducing anadditional stress field along this zone. As a function of the RH variation, this additional stressfield depends on not only the stress state and the crack opening along the cohesive zone butalso the material moisture ahead of the zone (undamaged material). In the finite elementanalysis, an algorithmic tangent operator is used to solve the non-linear problem combiningmechano-sorptive model and cohesive zone model including the effect of sudden RHvariations.The simulation of a notched structural element submitted to a constant load and cyclic RHvariations exhibits a strong coupling between the mechano-sorptive behavior and the effectof the RH variations on the fracture process zone (FPZ). This coupling results in an increaseof the crack propagation kinetic and leads to a precocious failure compared to those obtainedfrom the mechano-sorptive model or from the effect of sudden RH variations on the FPZ.Moreover, the coupling between the mechano-sorptive model and the effect of sudden RHvariations on the FPZ which cannot be predicted by a simple superposition of both effects,showing the interest of such a numerical approach in order to describe the complex behaviorof wood structural elements submitted to variable climatic conditions.

This thesis has discussed both properties and geometry of concrete slabs used in bridges. It gave understanding on behavior of concrete in both tension and compression zones and how crack propagates in specimens by presenting both theory of fracture and performing concrete tests like tension splitting, uniaxial compression and uniaxial tension tests. Furthermore, it supported experimental tests with finite elements modelling for each test, and illustrated both boundary conditions and loads. The thesis has used ARAMIS cameras to observe crack propagations in all experimental tests, and its first study at LNU that emphasized on Brazilian test, because of importance of this test to describe both crushing and cracking behavior of concrete under loading. It’s an excellent opportunity to understand how concrete and steel behave individually and in combination with each other, and to understand fracture process zone, and this has been discussed in theory chapter. The geometry change that could affect stresses distributions has also described in literature and modelled to give good idea on how to model slabs in different angles in the methodology chapter. Thus, thesis will use finite elements program (Abaqus) to model both experimental specimens and concrete slabs without reinforcement to emphasize on concrete behavior and skewness effect. This means studying both properties of concrete and geometry of concrete slabs. This thesis has expanded experimental tests and chose bridges as an application.

Notre nouvelle approche expérimentale consiste à étudier la fissuration de matériaux mous, principalement des gels polymériques et colloidaux, qui ont des tailles microstructurales micrométriques. Cette augmentation de la taille microscopique va avoir pour conséquence d’augmenter la taille de la zone de process et va rendre son observation plus facile avec des moyens standard de microscopie (à transmission et confocale).Pour se faire, nous avons mis au point un nouveau dispositif expérimental pour étudier la propagation de fissures dans des matériaux mous. Cette expérience permet de faire croître une fissure de manière contrôlée dans un échantillon mou et d’inspecter la pointe de fissure à haute résolution pour des fissures se propageant entre 1 µm/s and 1cm/s. En travaillant avec des gels de polymère physiques, nous avons analyse la forme de fissure ainsi que les champs de déplacement proches pointe (en utilisant des techniques de corrélation d’image) à petites et grandes échelles et à différentes vitesses. Nous avons montré qu’il existait une séparation d’échelles spatiales entre les échelles où l’élasticité linéaire s’applique, les échelles auxquelles les non linéarités émergent et les échelles auxquelles la dissipation se produit. Cette dernière échelle n’a pas pu être investigué dans le cas de gels polymériques. De récentes expériences sur des gels colloïdaux, ayant une longueur micro-structurale plus grande que celle des gels polymers, montre que nous sommes capables de sonder en temps réel les échelles d’endommagement lors de la fissuration. / Our novel experimental approach consists in studying fracture mechanics of soft materials, mainly polymer and colloidal gels, which have microstructures with large typical length scales. This increase in the microscopic length scale will consequently increase the typical size of the process zone and make its observation easier with standard microscopy techniques (optical or confocal).To do so, we designed a novel experimental device to study crack propagation in such soft materials. This experiment enables us to grow a unique crack in a controlled way in a soft specimen and to look at the crack tip at high magnification for crack velocities between 1 µm/s and 1cm/s. Working on physical polymer gels, we analysed the crack shape and crack displacement fields (using Digital Image Correlation) at large and intermediate scales for various velocities. We realized there was a separation of scales between the scale at which LEFM applies, the scale at which elastic nonlinearities emerge and the scale at which dissipation occurs. This last scale could not be investigated with the polymer gel. Recent experiments on colloidal gels, which have a microscopic length scale bigger than the one of polymer gels, show that we are able to probe damage at the microstructural scale.

La défaillance des matériaux et structures d'ingénierie peut être considéré comme le résultat d'une interaction complexe entre différents phénomènes physiques tels que la nucléation des cavités, les microfissures, les microvides et d'autres processus irréversibles. Ces micro-défauts se fondent éventuellement en une ou plusieurs macro-fissures conduisant à une diminution de la capacité portante et finalement à une défaillance de la structure considérée. La prévention des défaillances des structures et des composants structurels a toujours été un sujet important et une préoccupation majeure en ingénierie. Cette thèse vise à représenter une défaillance localisée dans des matériaux non linéaires sans dépendance de maillage. Un intérêt particulier sera le cas de l’adoucissement dynamique des déformations. Les phénomènes localisés sont pris en compte en utilisant l'approche des discontinuités embarquées fortes dans laquelle le champ de déplacement est amélioré pour capturer la discontinuité. Sur la base de cette approche, on a d'abord développé un modèle unidimensionnel de barres élasto-plastique capable de représenter une défaillance pour des matériaux ductiles avec un durcissement combiné dans une zone de processus de fracture FPZ et un adoucissement avec des discontinuités fortes encastrées. Les résultats comparant le modèle unidimensionnel proposé aux travaux (semi-) analytiques sont présentés. Il a été démontré que la stratégie proposée offre des solutions indépendantes de maillage. La déformation augmente dans le domaine de l’adoucissement avec une diminution simultanée de la contrainte. Le problème se décharge élastiquement à l'extérieur de la zone d’adoucissement de déformation. L'énergie dissipée se trouve à disparaître. Le modèle a également été comparé à un modèle de dommage unidimensionnel capable de représenter la fracture dynamique de la barre d'endommagementélasto-endommagée dans la zone de traitement de fracture - FPZ et de adoucissement avec de discontinuités fortes encastrées pour trouver un bon accord entre deux modèles. Un modèle d'éléments finis bidimensionnel a été développé, capable de décrire à la fois le mécanisme de dommage diffus accompagné d'un durcissement initial et d'une réponse d’adoucissement ultérieure de la structure. On a analysé les résultats de plusieurs simulations numériques effectuées sur des essais mécaniques classiques sous des charges progressivement croissantes telles que le test Brésilien ou le test de flexion en trois points. Le cadre de dynamique proposé est montré pour augmenter la robustesse de calcul. On a constaté que la direction finale des macro-fissures est assez bien prédite et que l'influence des effets d'inertie sur les solutions obtenues est assez modeste notamment en comparaison entre différentes mailles. Ce modèle bidimensionnel a été étendu plus loin dans le modèle bidimensionnel de discontinuité intégrée en viscodamage pour aider à explorer brièvement la mise en œuvre du schéma de point intermédiaire de second ordre qui peut fournir des résultats améliorés sous limitation de la régularisation visqueuse du modèle de dégâts localisés. / Failure of engineering materials and structures can be considered as a result of a complex interplay between different physical phenomena such as nucleation of cavities, microcracks,microvoids and other irreversible processes. These micro-defects eventually coalesce into one or more macro-cracks leading to a decrease in the load-bearing capability and finally, to failure of the structure under consideration. Prevention of failure of structures and structural parts has always been a critical subject and a major concern in engineering. This thesis aims to represent localized failure in non linear materials without mesh dependency. Of special interest will be the case of dynamic strain-softening. Localized phenomena are taken into account by using the embedded strong discontinuities approach in which the displacement field is enhanced to capture the discontinuity. Based upon this approach, a one-dimensional model for elasto-plastic bar capable of representing failure for ductile materials with combined hardening in FPZ-fracture process zone and softening with embedded strong discontinuities was first developed. Results comparing the proposed one-dimensional model to (semi-) analytical works are presented. It was shown that the proposed strategy provides mesh independent solutions. Strain increases in the softening domain with a simultaneous decrease of stress. The problem unloads elastically outside the strain softening region. The strain energy is found to vanish. The model was also compared with a one dimensional damage model capable of representing the dynamic fracture for elasto-damage bar with combined hardening in fracture process zone - FPZ and softening with strong embedded discontinuities to find a good agreement between two models. A two-dimensional finite element model was developed, capable of describing both the diffuse damage mechanism accompanied by initial strain hardening and subsequent softening response of the structure. The results of several numerical simulations, performed on classical mechanical tests under slowly increasing loads such as Brazilian test or three-point bending test were analyzed. The proposed dynamics framework is shown to increase computational robustness. It was found that the final direction of macro-cracks is predicted quite well and that influence of inertia effects on the obtained solutions is fairly modest especially in comparison among different meshes. This two-dimensional model was expanded further into the two dimensional continuum viscodamage-embedded discontinuity model to help briefly explore the implementation of the second order mid-point scheme that can provide improved results under limitation of viscous regularization of localized failure damage model.

Le dimensionnement des joints collés est une préoccupation majeure du CNES pour lesapplications spatiales des futurs lanceurs. Pour dimensionner une structure collée, il est nécessaire depouvoir apprécier les caractéristiques mécaniques du joint collé.Dans cette étude, trois adhésifs structuraux ont été sélectionnés (Hysol®EA 9321, Hysol®EA9394 et Hysol EA® 9395). Après leur caractérisation massique, une étude statistique pour mettre enévidence les effets des différents paramètres (vitesse d’essai, géométrie éprouvette, le degré depolymérisation…) a été entreprise.La deuxième étape a pour objectif de fiabiliser l’analyse des essais de fissuration etd’améliorer la compréhension des mécanismes d’endommagement et de propagation de fissure dansles liaisons collées. Trois types d’essai ont été utilisés, à savoir, l’essai Double Cantilever Beam(DCB), pour l’étude du mode I, l’essai End Notched Flexure (ENF), pour le mode II, et l’essai MixedMode Bending (MMB), pour les chargements en mode mixte I/II. Nous avons développé de nouvellesinstrumentations et méthodologies d’analyse. Pour affiner le protocole de test standard, la techniquedite de « backface strain monitoring » a été utilisée. Elle consiste à positionnées des jauges dedéformation sur les surfaces de l’éprouvette de façon à enregistrer l’évolution du signalextensométrique durant la propagation de la fissure. Cette méthode permet une meilleure estimation dela position front de fissure ainsi que l'étude de la répartition des contraintes le long du joint de colle.La corrélation d'images numériques (DIC) a également été utilisée afin de proposer un nouveauprotocole de calibrage de la longueur de fissure et pour comparer un modèle analytique (poutre deTimoshenko sur fondation élastique) avec les résultats expérimentaux. / Adhesive bonding is being strongly considered in space applications CNES as anadvantageous assembly technique for future launchers. Correct design of adhesive joints is of majorconcern. Aerospace adhesives are tough viscoelastic matrices (special epoxy resins) reinforced withnano-, or microparticles. Extended use of adhesive joints in structural applications is limited due to thedifficulties in predicting in-service performance, frequently leading to over-conservative design.Three structural adhesives (Hysol®EA 9321, Hysol®EA 9394 and Hysol®EA 9395) wereselected. After their bulk characterization, statistical studies to highlight effects of different parameterse.g. speed, test piece geometry, degree of polymerization were undertaken.In the second stage, fracture mechanics tests were effected employing: the double cantileverbeam (DCB) configuration (mode I characterisation), the three point bending end-notched flexure(ENF) (mode II) and the mixed-mode bending (MMB) (combined mode I/II loading). Crack growth inbonded joints was investigated in a novel way. To refine standard test protocol, the backface strainmonitoring technique was used. Strain gauges were used to measure the strain on the exposed skin ofthe adherends during crack onset and propagation. This method allows better estimation of the crackfront position as well as fine investigation of the stress distribution along the bondline and in the crackfront vicinity. Digital image correlation (DIC) was also used to compare analytical models, e.g.Timoshenko beam on elastic foundation model with experimental results.

Ce travail concerne des expériences mécaniques, des analyses numériques et des modélisations analytiques de la rupture cohésives (Mode I), vis-à-vis de l’étude du comportement mis en évidence par la courbe de Résistance (Courbe-R) et l’effet d’échelle de structures entaillées en bois massif. Des expériences de fissuration sont combinées à des analyses numériques pour déterminer les propriétés de rupture au moyen d’une procédure appelée Théorie de la Mécanique de la Rupture Linéaire Élastique équivalente (TMRLE), basée sur la complaisance de la structure. La courbe-R, obtenue à partir des expériences, selon une méthode de correction du poids propre, montre l’existence d’un domaine endommagé (Zone de Processus de Rupture) de taille non négligeable se développant en fond de fissure. Dans des conditions de fissuration stationnaire, ce domaine atteint une taille critique, et l’énergie nécessaire pour faire propager la fissure avec ce domaine endommagé (par unité de surface de rupture), reste constante. Le taux de libération de l’énergie de fissuration ainsi attendu, joue un rôle important en Mécanique de la Rupture, car il est possible simuler le comportement quasi-fragile du matériau en combinaison avec les autres propriétés de cohésion. La loi d’effet d’échelle de Bažant, utilisée pour prévoir l’influence de la taille sur la contrainte nominale, est estimée à partir de la réunion de deux comportements asymptotiques réalisés sur de petites tailles (Analyse limite ou RdM) et des grandes tailles. Une procédure analytique est présentée pour évaluer le comportement asymptotique additionnel exhibé par la contrainte nominale dans le régime intermédiaire, de façon plus exacte. Une validation numérique est présentée, et l’information expérimentale vient confirmer ce comportement asymptotique. / This work concerns the mechanical testing, numerical analysis and modelling of cohesive fracture (Mode I) on the purpose to study the Resistance-curve behaviour and the size effect in wooden notched structures. The mechanical testing is combined with the numerical analysis to evaluate fracture properties by means of an equivalent LEFM approach based on the structure compliance. The Resistance-curve being revealed from the experiments, by means of a self-weight compensation method, correction puts into evidence that a non-negligible damaged domain (Fracture Process Zone) is under development in the crack front during the loading process. This being the case, among other fracture parameters issued from the Resistance-curve, the critical (asymptotic) energy release rate is determined, turning possible to use it in combination with other cohesive crack properties in the crack modelling (in Mode I). Thus, for a given geometry it is possible to monitor the critical dimension being revealed by the Fracture Process Zone (FPZ) during the crack propagation. The well known Bažant’s size effect law provides the scaling of the nominal strength through the asymptotic matching performed both on the small (Strength Theory) and on the large (LEFM) structure sizes. An analytical procedure is proposed to determine an additional asymptotic regime in the intermediate size range through a more accurate manner. Numerical validations of the proposed procedure are made and experimental data is presented revealing the scaling of the nominal strength through an envelop of values.

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.