A multi-scale based cohesive zone model for the analysis of thickness scaling effect in fiber bridging

John Botsis, Luis Pablo Canal Casado
Elsevier, 2017
Article

Résumé

A cohesive zone model (CZM) is proposed to assess the thickness scaling effect associated with fiber bridging during fracture. The CZM was developed through a multi-scale simulation approach and utilizes an embedded cell model of the Double Cantilever Beam (DCB) that explicitly accounts for the bridging bundles on the fracture plane. In particular, micromechanical simulations of failure were carried out, for varying arms thickness, in order to determine the homogenized fracture behavior. To model the observed scaling effect, the conventional cohesive law, formulated as an opening-stress relation, is enriched with information on the crack opening angle. Continuum finite element simulations indicated that the proposed CZM was able to mimic very well the essential features observed in the experiments, e.g. raising R-curve behavior and thickness scaling effect on the energy dissipated at steady state. (C) 2016 Elsevier Ltd. All rights reserved.

Official source

https://infoscience.epfl.ch/record/226903?ln=fr

À propos de ce résultat

Cette page est générée automatiquement et peut contenir des informations qui ne sont pas correctes, complètes, à jour ou pertinentes par rapport à votre recherche. Il en va de même pour toutes les autres pages de ce site. Veillez à vérifier les informations auprès des sources officielles de l'EPFL.

Concepts associés

Chargement

Publications associées

Chargement

Mixed-Mode Static and Fatigue Failure Criteria for Adhesively-Bonded FRP Joints

Moslem Shahverdi

Fiber-reinforced polymer (FRP) composites are increasingly used in many load-bearing structures. Joints are the most critical structural elements as they often become the weak links in engineering structures. Bonding and bolting are two of the main joining techniques used for composite elements but, although bolted and bonded joints have their advantages and disadvantages, adhesively-bonded joints are often preferred. This is thanks to their better performance in terms of cost-effective structures with uniformly distributed stresses without the need for cutting fastener holes that result in high stress concentrations in the FRP adherends at the edges of the holes. In addition, due to the absence of stress concentrations, adhesively-bonded joints exhibit better fatigue behavior. However, the efficient and reliable use of adhesively-bonded joints requires a thorough understanding of their failure mechanisms under mixed-mode quasi-static and fatigue loading conditions. In real structural applications, failure of adhesively-bonded joints occurs due to crack initiation and propagation under varying mixed-mode ratios. Since it is difficult to describe this progressive fracture in a structural joint, fracture joints with constant mode-mixity ratios are used to determine the strain energy release rates for crack initiation and propagation. Subsequently, mixed-mode failure criteria based on the determined strain energy release rates can be established and used to model/predict the fatigue and fracture behaviors of structural joints. The main objectives of this thesis were to understand the fracture and fatigue behaviors of adhesively-bonded FRP fracture joints under Mode I, Mode II, and mixed-Mode I/II loading conditions and, consequently, develop mixed-mode static and fatigue failure criteria applicable to structural joints. To fulfill these objectives a combined experimental-analytical-numerical study was undertaken. Double cantilever beam (DCB-Mode I), end-load split (ELS-Mode II), and mixed-mode bending (MMB-Mode I/II) specimens were examined under quasi-static and constant amplitude fatigue loading. Analytical approaches were used to determine the strain energy release rates under quasi-static and fatigue loading, whereas phenomenological models were developed to characterize the fatigue and fracture behaviors of the joints. A new approach designated the “extended global method” was established and applied for the analysis of the quasi-static and fatigue experimental data and the fracture mode partitioning in the mixed-mode experiments. In parallel, non-linear finite element models were developed to simulate the quasi-static fracture behavior and investigate the effects of both the existing asymmetry and the observed fiber bridging on the fracture behavior of the examined joints. The bridging zone was modeled by zero-thickness cohesive elements. An exponential traction-separation description of the cohesive zone model was shown to appropriately model the fiber bridging and calculate its contribution to the fracture energy under quasi-static loading. The FE models were successfully employed for the separation of the fracture parameters, i.e. strain energy released at the crack tip (Gtip) and the amount of energy released along the crack-bridging zone (Gbr), in the quasi-static mixed-mode failure criterion for crack propagation. In addition to the studies on the quasi-static fracture behavior of the joints, a new phenomenological fatigue crack growth formulation for the modeling and prediction of Rratio effects on the Mode I fatigue behavior of adhesively-bonded joints was also derived. The model was used for prediction of the fatigue behavior of the examined joints under different R-ratios. The results proved that the model could accurately simulate the fatigue behavior of the examined joints and was also capable of predicting behavior exhibited under unseen loading conditions, i.e. the R-ratios used for estimating the model parameters were different from those used for validating its prediction capability. The experimental results and numerical analyses combined with the phenomenological models were used to establish mixed-mode static and fatigue failure criteria for crack initiation and crack propagation for the adhesively-bonded joints. The derived mixed-mode failure criteria can be used for simulating/predicting the crack initiation and progressive crack propagation in structural joints comprising the same adhesive and adherends.

EPFL2013

Chargement

Modèle numérique du comportement non-linéaire d'ouvrages massifs en béton non armé

Patrice Droz

Linear Elastic Fracture Mechanics (LEFM) are used to study crack propagation in unreinforced concrete of massive structures. This approach is well suited to the nonlinear analysis of structures of large dimensions, taking into account the size effect, i.e., the influence of the dimension of the structure on its behaviour at rupture. The Finite Element method is used to calculate the structure and the hypothesis of plane strain is made. The use of a criterion of crack propagation, based on LEFM requires the calculation of stress intensity factors. The evaluation of the latter is performed by the mean of a surface integral defined around the tip of the crack studied. It has been shown in this work that this integral is derived from the path integral J. The use of the surface integral has also been extended to the cases where body forces (gravity, inertie) act, or when the edges of the crack are subjected to pressure. This method is precise and numerically efficient. A smeared crack model is used in order to avoid continuous remeshing during crack propagation. But, as it has been shown in this research, classical smeared crack elements do not give satisfactory results when the crack is submitted to shear loading, a new finite element, using smeared cracking but with discontinuous shape functions has been developed. The model which combines LEFM and the smeared crack approach has been applied to different classical problems of fracture mechanics. It leads to good results under static or dynamic loading. As an example, the nonlinear behaviour of the vertical cross section of a gravity dam during an earthquake has been calculated and the crack pattern identified. Many further developments could be done starting from this original approach, in order to simulate more closely the real behaviour of structures, taking into account friction in the cracks and the three-dimensional development of then.

EPFL1987

Chargement

In this study, the tensile and fracture properties and the microstructure of the reduced-activation tempered martensitic steel Eurofer97 have been investigated. This technical alloy is a 9%Cr steel developed within the European fusion material research program. To a lesser extend, the plastic flow properties of a equiaxed ferritic Fe9%Cr model alloy were also studied for comparison with those of the tempered martensitic structure. The main objectives of this work are described hereafter: to correlate the microstructural features with the plastic flow properties measured by the tensile tests for both the Eurofer97 steel and the model alloy. The correlation established should be reflected in a physically-based model of plastic flow. to study the fracture properties of the Eurofer97 steel in details in the lower ductile-to-brittle transition region. to calculate with finite element modeling the stress fields at the crack tip. This information is further used in conjunction with a local approach model for quasi-cleavage to reconstruct the fracture toughness-temperature curve. Tensile tests were carried out at different imposed nominal strain rate at several temperatures from 77K up to 473K, on both the Eurofer97 steel and the Fe9%Cr model alloy. The temperature dependence of the yield stress was precisely determined. As expected for body-centered cubic (BCC) materials, a strong increase of the yield stress by decreasing temperature, below 200 K, was observed. At higher temperatures, the temperature dependence of the yield stress was found much weaker, being associated with the temperature dependence of the shear modulus. Efforts were made to analyze in details the post-yield behavior (strain-hardening) as a function of temperature. The post-yield behavior was modeled using the Kocks phenomenological model based on the competition of storage and annihilation of dislocations. While this model was originally developed for face-centered cubic (FCC) metals where the rate-controlling mechanism of dislocation motion is the dislocation-dislocation interaction, we used is model in the high temperature domain (T>200K) of BCC materials, to model the strain-hardening evolution of the Eurofer97 steel and the model alloy. The values obtained for the key parameters of the model, namely the dislocation mean free path and annihilation coefficient, were found consistent with the microstructural features. The parameters temperature-dependence observed was also consistent with the physics of two basic mechanisms of dislocation storage and annihilation (dynamical recovery). For the low temperature domain, the strain-hardening model was modified to account for the strong Peierls lattice friction, based on an original idea of Rauch. Transmission electron microscopy observations were done to characterize the undeformed microstructures and their evolution with strain. A clear correlation was established between the stress-dependence of the strain-hardening and the microstructures. Fracture toughness tests on the Eurofer97 steel were performed with 0.2T C(T) and 0.4T C(T) specimens in the lower transition. Five temperatures were selected at which many tests were repeated to determine the amplitude of the inherent scatter of quasi-cleavage. The temperature were chosen in such a way that the measure fracture toughness remain below 150 MPa m1/2, which correspond for the 0.4T C(T) at a M value larger than 70 at the highest temperature. Such a M value for 0.4T C(T) specimens is known to ensure enough constraint. An attempt to analyze the experimental data in the framework of the master-curve approach by following the ASTM E-1921 standard was done. It was clearly demonstrated that the assumed shape of the toughness-temperature curve as described in the ASTM E-1921 standard for the reactor pressure vessel steels is not adequate for the tempered martensitic Eurofer97 steel, which present a particularly steep transition. Our Eurofer97 fracture database was then compared to the existing one on another similar steel, the F82H. Differences in the amplitude of the scatter of both steels was found while the lower bound of the toughness-temperature curve, describing the 1% failure probability was shown to be the same. 2D finite element simulations of the compact tension specimens were performed at various temperatures using the constitutive laws determined previously. The stress field around the crack tip were calculated and used to determine a local criterion of quasi-cleavage. The criterion is defined by the attainment of a critical stress encompassing a critical area. The lower bound of the toughness-temperature curve was then successfully reconstructed by using this local criterion. Finally, the relationship between the critical area and the applied stress intensity factor for the C(T) specimen was shown to follow a power law whose coefficients are dependent on the real dimensions of the specimen. Such a relationship allows scaling the C(T) toughness data from one size to another in case of in-plane constraint loss.

EPFL2006

Chargement

Mixed-Mode Static and Fatigue Failure Criteria for Adhesively-Bonded FRP Joints

Moslem Shahverdi

EPFL2013

EPFL2006

Modèle numérique du comportement non-linéaire d'ouvrages massifs en béton non armé

Patrice Droz

EPFL1987