A New Finite Element Method for Strain Gradient Theories and Applications to Fracture Analyses

A new compatible finite element method for strain gradient theories is presented. In the new finite element method, pure displacement derivatives are taken as the fundamental variables. The new numerical method is successfully used to analyze the simple strain gradient problems – the fundamental fracture problems. Through comparing the numerical solutions with the existed exact solutions, the effectiveness of the new finite element method is tested and confirmed. Additionally, an application of the Zienkiewicz–Taylor C1 finite element method to the strain gradient problem is discussed. By using the new finite element method, plane-strain mode I and mode II crack tip fields are calculated based on a constitutive law which is a simple generalization of the conventional J2 deformation plasticity theory to include strain gradient effects. Three new constitutive parameters enter to characterize the scale over which strain gradient effects become important. During the analysis the general compressible version of Fleck–Hutchinson strain gradient plasticity is adopted. Crack tip solutions, the traction distributions along the plane ahead of the crack tip are calculated. The solutions display the considerable elevation of traction within the zone near the crack tip.

Fractals and Scaling in Fracture Induced by Microcrack Coalescence

A numerical model is proposed to simulate fracture induced by the coalescence of numerous microcracks, in which the condition for coalescence between two randomly nucleated microcracks is determined in terms of a load-sharing principle. The results of the simulation show that, as the number density of nucleated microcracks increases, stochastic coalescence first occurs followed by a small fluctuation, and finally a newly nucleated microcrack triggers a cascade coalescence of microcracks resulting in catastrophic failure. The fracture profiles exhibit self-affine fractal characteristics with a universal roughness exponent, but the critical damage threshold is sensitive to details of the model. The spatiotemporal distribution of nucleated microcracks in the vicinity of critical failure follows a power-law behaviour, which implies that the microcrack system may evolve to a critical state.

Effects of Nanostructures on the Fracture Strength of the Interfaces in Nacre

A strengthening mechanism arising from a type of inorganic nanostructure in the organic matrix layers is presented by studying the structural and mechanical properties of the interfaces in nacre. This nanostructural mechanism not only averagely increases the fracture strength of the organic matrix interfaces by about 5 times, but also effectively arrests the cracks in the organic matrix layers and causes the crack deflection in this biomaterial. The present investigation shows that the main mechanism governing the strength of the organic matrix interfaces relies on the inorganic nanostructures rather than the organic matrix. This study provides a guide to the interfacial design of synthetic materials.

Essential and Non-Essential Work of Fracture of PI/SiO2 Hybrid Thin Films

Essential work of fracture (EWF) analysis is used to study the effect of the silica doping level on fracture toughness of polyimide/silica (PI/SiO2) hybrid films. By using double-edge-notched-tension (DENT) specimens with different ligament lengths, it seems that the introduction of silica additive can improve the specific essential work of fracture (w (e) ) of PI thin films, but the specific non-essential work of fracture (beta w (p) ) will decease significantly as the silica doping level increasing from 1 to 5 wt.%, and even lower than that of neat PI. The failure process of the fracture is investigated with online scanning electron microscope (SEM) observation and the parameters of non-essential work of fracture, beta and w (p) , are calculated based on finite element (FE) method.

The relationship of SIF between plate and plane fracture problems and the effect of the plate thickness on SIF

In the present paper, it is shown that the zero series eigenfunctions of Reissner plate cracks/notches fracture problems are analogous to the eigenfunctions of anti-plane and in-plane. The singularity in the double series expression of plate problems only arises in zero series parts. In view of the relationship with eigen-values of anti-plane and in-plane problem, the solution of eigen-values for Reissner plates consists of two parts: anti-plane problem and in-plane problem. As a result the corresponding eigen-values or the corresponding eigen-value solving programs with respect to the anti-plane and in-plane problems can be employed and many aggressive SIF computed methods of plane problems can be employed in the plate. Based on those, the approximate relationship of SIFs between the plate and the plane fracture problems is figured out, and the effect relationship of the plate thickness on SIF is given.

Micro- and macroscopic characterizations of the viscoelastic fracture of resin-based fibre composites

Micro- and macroscopic characterizations of the viscoelastic fracture of a unidirectional carbon-fibre-reinforced epoxy composite are presented. First, the micro-cracking behavior of the material is studied by the use of scanning electron microscopy; the in situ creep cracking process is observed and the crack propagation is measured. In order to obtain insight into the mechanisms of the observed creep cracking, macroscopic investigations were also carried out. Finite-element method simulations were carried out to calculate the stress distribution and the variation of stresses with time. A theoretical analysis of the orthotropy of viscoelastic fracture behavior of the material is also conducted.

Modeling Nonlinear Peeling of Ductile Thin Films-Critical Assessment of Analytical Bending Models Using FE Simulations

Three analytical double-parameter criteria based on a bending model and a two-dimensional finite element analysis model are presented for the modeling of ductile thin film undergoing a nonlinear peeling process. The bending model is based on different governing parameters: (1) the interfacial fracture toughness and the separation strength, (2) the interfacial fracture toughness and the crack tip slope angle, and (3) the interfacial fracture toughness and the critical Mises effective strain of the delaminated thin film at the crack tip. Thin film nonlinear peeling under steady-state condition is solved with the different governing parameters. In addition, the peeling test problem is simulated by using the elastic-plastic finite element analysis model. A critical assessment of the three analytical bending models is made by comparison of the bending model solutions with the finite element analysis model solutions. Furthermore, through analyses and comparisons for solutions based on both the bending model and the finite element analysis model, some connections between the bending model and the finite element analysis model are developed. Moreover, in the present research, the effect of different selections for cohesive zone shape on the ductile film peeling solutions is discussed.

Theermal Stress and Fracture in Transparent Optical Materials due to Laser Energy-Absorption by a Defect

利用热弹性理论分析了在光学材料中由于缺陷吸收激光能量引起的温度和热应力分布,并且针对一个简单的裂纹模型分析了热应力产生的应力强度因子,给出了一些主要参数对于应力强度因子的影响的规律。

Latin hypercube sampling used in the calculation of the fracture probability

The paper proposes Latin hypercube sampling combined with the stratified sampling of variance reduction technique to calculate accurate fracture probability. In the compound sampling, the number of simulations is relatively small and the calculation error is satisfactory.

The fracture energy of metal fibre reinforced ceramic composites (MFCs)

A model is presented for prediction of the fracture energy of ceramic-matrix composites containing dispersed metallic fibres. It is assumed that the work of fracture comes entirely from pull-out and/or plastic deformation of fibres bridging the crack plane. Comparisons are presented between these predictions and experimental measurements made on a commercially-available composite material of this type, containing stainless steel (304) fibres in a matrix predominantly comprising alumina and alumino-silicate phases. Good agreement is observed, and it's noted that there is scope for the fracture energy levels to be high (~20kJm-2). Higher toughness levels are both predicted and observed for coarser fibres, up to a practical limit for the fibre diameter of the order of 0.5mm. Other deductions are also made concerning strategies for optimisation of the toughness of this type of material. © 2010 Elsevier Ltd.

The investigation on dynamic fracture behavour of materials under compressive-shear combined stress waves

In this paper, an improved plate impact experimental technique is presented for studying dynamic fracture mechanism of materials, under the conditions that the impacting loading is provided by a single pulse and the loading time is in the sub-microsecond range. The impacting tests are carried out on the pressure-shear gas gun. The loading rate achieved is dK/dt similar to 10(8) MPa m(1/2) s(-1). With the elimination of influence of the specimen boundary, the plane strain state of a semi-infinite crack in an infinite elastic plate is used to simulate the deformation fields of crack tip. The single pulses are obtained by using the "momentum trap" technique. Therefore, the one-time actions of the single pulse are achieved by eradicating the stress waves reflected from the specimen boundary or diffracted from the crack surfaces. In the current study, some important phenomena have been observed. The special loading of the single pulse can bring about material damage around crack tip, and affect the material behavior, such as kinking and branching of the crack propagation. Failure mode transitions from mode I to mode II crack are observed under asymmetrical impact conditions. The mechanisms of the dynamic crack propagation are consistent with the damage failure model.

Characterization of the Fracture Work for Ductile Film Undergoing the Micro-Scratch

The interface adhesion strength (or interface toughness) of a thin film/substrate system is often assessed by the micro-scratch test. For a brittle film material, the interface adhesion strength is easily obtained through measuring the scratch driving forces. However, to measure the interface adhesion strength (or interface toughness) for a metal thin film material (the ductile material) by the microscratch test is very difficult, because intense plastic deformation is involved and the problem is a three-dimensional elastic-plastic one. In the present research, using a double-cohesive zone model, the failure characteristics of the thin film/substrate system can be described and further simulated. For a steady-state scratching process, a three-dimensional elastic-plastic finite element method based on the double cohesive zone model is developed and adopted, and the steady-state fracture work of the total system is calculated. The parameter relations between the horizontal driving forces (or energy release rate of the scratching process) and the separation strength of thin film/substrate interface, and the material shear strength, as well as the material parameters are developed. Furthermore, a scratch experiment for the Al/Si film/substrate system is carried out and the failure mechanisms are explored. Finally, the prediction results are applied to a scratch experiment for the Pt/NiO material system given in the literature.

Steady-state crack growth and fracture work based on the theory of mechanism-based strain gradient plasticity

Mode I steady-state crack growth is analyzed under plane strain conditions in small scale yielding. The elastic-plastic solid is characterized by the mechanism-based strain gradient (MSG) plasticity theory [J. Mech. Phys. Solids 47 (1999) 1239, J. Mech. Phys. Solids 48 (2000) 99]. The distributions of the normal separation stress and the effective stress along the plane ahead of the crack tip are computed using a special finite element method based on the steady-state fundamental relations and the MSG flow theory. The results show that during the steady-state crack growth, the normal separation stress on the plane ahead of the crack tip can achieve considerably high value within the MSG strain gradient sensitive zone. The results also show that the crack tip fields are insensitive to the cell size parameter in the MSG theory. Moreover, in the present research, the steady-state fracture toughness is computed by adopting the embedded process zone (EPZ) model. The results display that the steady-state fracture toughness strongly depends on the separation strength parameter of the EPZ model and the length scale parameter in the MSG theory. Furthermore, in order for the results of steady crack growth to be comparable, an approximate relation between the length scale parameters in the MSG theory and in the Fleck-Hutchinson strain gradient plasticity theory is obtained.

Thermal fracture characteristics induced by laser beam

The damage mechanism of a cracked material due to laser beam thermal shock is an important problem when the interactions of high power laser beam with materials are studied. The transient thermal stress intensity factors (TSIFs) for a center crack in an infinite plate subjected to laser beam thermal shock are investigated. When the crack is in the heat affected zone, the compressive thermal stress causes the crack surface to come into contact with each other over a certain contact length, but the high tensile stresses around the crack tip tend to open the crack. In this case, the problem may be treated as a pair of embedded cracks problem with a smooth closure condition of the center crack. The TSIFs and the crack contact lengths are calculated with different laser beam spatial shapes. The TSIFs induced by thermal shock are in marked different from those induced by general mechanical loading.