Detonation-driven fracture in thin shell structures: Numerical studies

Detonation-driven fracture of thin structures is studied numerically by a 3D discrete crack meshfree method. These types of failure mechanisms play an important role in pipes and vessels. I therefore proposed a three-dimensional meshfree method and an efficient discrete crack model to describe crack propagation. The method is based on separation of particles similar to the visibility method but its implementation is more efficient. I assume here through-thickness cracks though the method can be extended to crack growth in arbitrary directions. The load is applied as travelling pressure wave obtained from pure fluid simulation in accordance with experimental measurements. Numerical results to experimental data show good agreement.     

Out-of-plane displacement measurement near the crack tip during a crack propagation: Validation of a 3D formulation for a specimen in PMMA loaded in mode I.

The fundamental aim of this study is the determination zone of the 3D effects and the transient one at the vicinity of the crack tip during a crack propagation in brittle materials ( PMMA ) using an optical method (Michelson interferometer). With the obtained interferograms, we can extract the phase (thus the relief) by using a new numerical approach based on the principle of images correlation between real fringes and virtual fringes. Different dynamic tests are realized by a plate loaded in mode I under a constant loading. We compare the obtained data with the two-dimensional theory of Westergaard (plane stress hypothesis) [1]. With the divergence is established, we propose a new 3D formulation, based on a formulation employed for static crack, which takes into account 3D and transient effects. For the static cracks, the 3D effects relate to a presence of the state of three-dimensional stresses. However in dynamics, the transient effects appear and are related to the crack propagation velocity. The 3D effects and transient effects lead to results equivalent to experimental ones in terms of displacement but are completely different to results given by the two-dimensional theory near the crack tip. It is possible to quantify the zone when the plane stress hypothesis is not valid according to the crack propagation speed V. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Shape sensitivity of a plane crack front

The 3D‐elasticity model of a solid with a plane crack under the stress‐free boundary conditions at the crack is considered. We investigate variations of a solution and of energy functionals with respect to perturbations of the crack front in the plane. The corresponding expansions at least up to the second‐order terms are obtained. The strong derivatives of the solution are constructed as an iterative solution of the same elasticity problem with specified right‐hand sides. Using the expansion of the potential and surface energy, we consider an approximate quadratic form for local shape optimization of the crack front defined by the Griffith criterion. To specify its properties, a procedure of discrete optimization is proposed, which reduces to a matrix variational inequality. At least for a small load we prove its solvability and find a quasi‐static model of the crack growth depending on the loading parameter. Copyright © 2003 John Wiley & Sons, Ltd.     

Analysis solution method for 3D planar crack problems of two-dimensional hexagonal quasicrystals with thermal effects

An analysis solution method (ASM) is proposed for analyzing arbitrarily shaped planar cracks in two-dimensional (2D) hexagonal quasicrystal (QC) media. The extended displacement discontinuity (EDD) boundary integral equations governing three-dimensional (3D) crack problems are transferred to simplified integral-differential forms by introducing some complex quantities. The proposed ASM is based on the analogy between these EDD boundary equations for 3D planar cracks problems of 2D hexagonal QCs and those in isotropic thermoelastic materials. Mixed model crack problems under combined normal, tangential and thermal loadings are considered in 2D hexagonal QC media. By virtue of ASM, the solutions to 3D planar crack problems under various types of loadings for 2D hexagonal QCs are formulated through comparison to the corresponding solutions of isotropic thermoelastic materials which have been studied intensively and extensively. As an application, analytical solutions of a penny-shaped crack subjected uniform distributed combined loadings are obtained. Especially, the analytical solutions to a penny-shaped crack subjected to the anti-symmetric uniform thermal loading are first derived for 2D hexagonal QCs. Numerical solutions obtained by EDD boundary element method provide a way to verify the validity of the presented formulation. The influences of phonon-phason coupling effect on fracture parameters of 2D hexagonal QCs are assessed.     

3D modeling of crack growth in electro-magneto-thermo-elastic coupled viscoplastic multiphase composites

This work presents a time-domain hypersingular integral equation (TD-HIE) method for modeling 3D crack growth in electro-magneto-thermo-elastic coupled viscoplastic multiphase composites (EMTE-CVP-MCs) under extended incremental loads rate through intricate theoretical analysis and numerical simulations. Using Green’s functions, the extended general incremental displacement rate solutions are obtained by time-domain boundary element method. Three-dimensional arbitrary crack growth problem in EMTE-CVP-MCs is reduced to solving a set of TD-HIEs coupled with boundary integral equations, in which the unknown functions are the extended incremental displacement discontinuities gradient. Then, the behavior of the extended incremental displacement discontinuities gradient around the crack front terminating at the interface is analyzed by the time-domain main-part analysis method of TD-HIE. Also, analytical solutions of the extended singular incremental stresses gradient and extended incremental integral near the crack fronts in EMTE-CVP-MCs are provided. In addition, a numerical method of the TD-HIE for a 3D crack subjected to extended incremental loads rate is put forward with the extended incremental displacement discontinuities gradient approximated by the product of time-domain basic density functions and polynomials. Finally, examples are presented to demonstrate the application of the proposed method.     

Efficient crack growth analyzes by combining fast methods for BEM

The combination of fast methods for the boundary element method (BEM) for efficient crack growth analyzes is presented. Due to the nonlinearity of fatigue crack growth an incremental procedure has to be applied. Within each increment a stress analysis is needed. Based on the asymptotic stress field the stress intensity factors (SIFs) are calculated by an extrapolation method. Then, a new crack front is determined by a reliable 3D crack growth criterion. Finally, the numerical model has to be updated for the next increment. The time dominant factor in each increment is the computation of the stress field. Due to the stress concentration problem the BEM is utilized. To speed-up the calculation several independent fast methods are exploited. An algebraic technique is the adaptive cross approximation (ACA) method which is acting on the system matrix itself. The application of the substructure technique leads to a blockwise band matrix and therefore to reduced memory requirements. Further savings in memory and computation time are reached by modelling cracks with the dual discontinuity method (DDM) and using the ACA method in each substructure. The efficiency of the combined methods is shown by a complex industrial example. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Influence of Multiple Micro Cracks on the Damage Behavior of a Macro-Crack Tip北大核心CSCD

The solution of an infinite plane containing a macro crack and a cluster of micro cracks under uniaxial tensile load was presented based on Muskhelishvili’s complex function method and the stepwise recursive method. The stress field and stress intensity factor K were obtained. Combined with the damage mechanics, damage parameter D of the macro-crack tip and the micro-crack tip under uniaxial tension was redefined, and the influence of different damage zone forms on the damage of the crack tip was analyzed. The results show that, both the chain-distribution and the reverse-chain-distribution micro cracks have an amplifying effect on the macro crack growth, and the damage parameter increases with the decrease of the inclination angle of the micro crack and the reduction of the distance between the macro crack and the micro cracks. For a relatively small inclination angle of the micro crack, the damage parameters of the macro crack and the micro crack heightens, and the damage parameter of the macro crack increases with the micro-crack length. For evenly distributed micro cracks in the continuous damage zone, the micro cracks have an amplifying effect on the macro-crack growth, and the damage parameter of the macro crack increases with the micro-crack number. © 2022 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Direct FEM Simulation and Replacement Crack Approach for Mixed Mode Crack Growth in a Unit Cell

A whole class of continuum damage models uses microcracks as the main source of reduction of stiffness. For the growth of these cracks mostly only mode I is considered. We want to present a method to describe mixed mode crack growth inside a unit cell with a crack, without the need of a direct FEM simulation of crack growth per integration point. We replace the infinitesimal grown and kinked crack with the help of a replacement crack model. This replacement method is mainly based on the equivalence of the dissipation of the original kinking and the replacement crack. The resulting evolution of the stiffness of the unit cell is compared to a direct FEM simulation of mixed mode crack growth. The crack growth criterion used is the principle of maximum energy release rate, which has shown to be a direct consequence of a variational principle of a body with a crack [1]. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Destruction of epoxy resin by 10−3-sec laser pulses

High-speed photography is used to investigate the destruction of epoxy resin by laser pulses. The results suggest the following mechanism is at work. Initially, structural changes (gas-bubble formation) take place near the focus of the beam. These lead to the formation of an initial crack of fixed minimum size. At first the crack grows rapidly, then more slowly, as it is screened by other cracks. The discontinuous crack growth predicted by the theory is observed, the jump velocities being an order greater than the mean rate of crack growth on the same interval. The total area of destruction increases linearly with time; at the beginning of the process this increase is primarily associated with an increase in the number of cracks, in the final stages with the further growth of the cracks already formed.Scientific-Research Institute of Mechanics, Moscow Lomonosov State University. Translated from Mekhanika Polimerov, No. 3, pp. 460–464, May–June, 1969.     

Crack nucleation sensitivity analysis

A simple analytical expression for crack nucleation sensitivity analysis is proposed relying on the concept of topological derivative and applied within a two‐dimensional linear elastic fracture mechanics theory (LEFM). In particular, the topological asymptotic expansion of the total potential energy together with a Griffith‐type energy of an elastic cracked body is calculated. As a main result, we derive a crack nucleation criterion based on the topological derivative and a criterion for determining the direction of crack growth based on the topological gradient. The proposed methodology leads to an axiomatic approach of crack nucleation sensitivity analysis. Copyright © 2010 John Wiley & Sons, Ltd.     

Bayesian prediction of crack growth based on a hierarchical diffusion model

A general Bayesian approach for stochastic versions of deterministic growth models is presented to provide predictions for crack propagation in an early stage of the growth process. To improve the prediction, the information of other crack growth processes is used in a hierarchical (mixed‐effects) model. Two stochastic versions of a deterministic growth model are compared. One is a nonlinear regression setup where the trajectory is assumed to be the solution of an ordinary differential equation with additive errors. The other is a diffusion model defined by a stochastic differential equation where increments have additive errors. While Bayesian prediction is known for hierarchical models based on nonlinear regression, we propose a new Bayesian prediction method for hierarchical diffusion models. Six growth models for each of the two approaches are compared with respect to their ability to predict the crack propagation in a large data example. Surprisingly, the stochastic differential equation approach has no advantage concerning the prediction compared with the nonlinear regression setup, although the diffusion model seems more appropriate for crack growth. Copyright © 2016 John Wiley & Sons, Ltd.     

A BE based finite element for crack propagation processes

A boundary element based finite macro element for the simulation of 3D crack propagation in the framework of linear elastic fracture mechanics is presented. While the major part of the numerical model is discretized with finite elements, a small domain containing the crack is meshed with boundary elements. By means of the Symmetric Galerkin BEM a stiffness formulation for the cracked BE domain is obtained which enables a direct FEM/BEM coupling. All necessary operations for the crack propagation are carried out within this boundary element based finite macro element and exploit the potential of the boundary integral formulation. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Analysis of arbitrarily shaped planar cracks in two-dimensional hexagonal quasicrystals with thermal effects. Part II: Numerical solutions

The extended displacement discontinuity (EDD) boundary element method is developed to analyze an arbitrarily shaped planar crack in two-dimensional (2D) hexagonal quasicrystals (QCs) with thermal effects. The EDDs include the phonon and phason displacement discontinuities and the temperature discontinuity on the crack face. Green's functions for uniformly distributed EDDs over triangular and rectangular elements for 2D hexagonal QCs are derived. Employing the proposed EDD boundary element method, a rectangular crack is analyzed to verify the Green's functions by discretizing the crack with rectangular and triangular elements. Furthermore, the elliptical crack problem for 2D hexagonal QCs is investigated. Normal, tangential, and thermal loads are applied on the crack face, and the numerical results are presented graphically.     

Three-Dimensional Variational Analysis of Composite Structures Using Bernstein Polynomial Approximations. Report 2

The application of the previously developed 3D varionational analysis approach to the investigation of crack propagation in composite bonded joints is presented. In this application, the propagation of three different types of a 2D planar crack (adhesive, cohesive, and interfacial) is modeled by relaxing the respective continuity conditions for displacements between adjacent bricks in the mosaic structure. The crack propagation process is then characterized by the release rate of the total potential energy between two consecutive states of the mosaic body with different crack lengths. Numerical examples illustrate the 3D analysis of double-lap adhesively bonded joints with unidirectional and cross-ply laminated composite adherends. The numerical results provide an illustration of various characteristics of the crack propagation process. The values of the ultimate failure load predicted by analyzing the initial stage of crack propagation are found to be in a good agreement with experimental data.     

Multiple crack fatigue growth modeling by displacement discontinuity method with crack-tip elements

This paper presents a numerical approach for modeling multiple crack fatigue growth in a plane elastic infinite plate. It involves a generation of Bueckner’s principle, a displacement discontinuity method with crack-tip elements (a boundary element method) proposed recently by the author and an extension of Paris’ law to a multiple crack problem under mixed-mode loading. Because of an intrinsic feature of the boundary element method, a general multiple crack growth problem can be solved in a single-region formulation. In the numerical simulation, for each increment of crack extension, remeshing of existing boundaries is not necessary. Crack extension is conveniently modeled by adding new boundary elements on the incremental crack extension to the previous crack boundaries. Fatigue growth modeling of an inclined crack in an infinite plate under biaxial cyclic loads is taken into account to illustrate the effectiveness of the present numerical approach. As an example, the present numerical approach is used to study the fatigue growth of three parallel cracks with same length under uniaxial cyclic load. Many numerical results are given.     

Analysis of arbitrarily shaped planar cracks in two-dimensional hexagonal quasicrystals with thermal effects. Part I: Theoretical solutions

An extended displacement discontinuity (EDD) boundary integral equation method is proposed for analysis of arbitrarily shaped planar cracks in two-dimensional (2D) hexagonal quasicrystals (QCs) with thermal effects. The EDDs include the phonon and phason displacement discontinuities and the temperature discontinuity on the crack surface. Green's functions for unit point EDDs in an infinite three-dimensional medium of 2D hexagonal QC are derived using the Hankel transform method. Based on the Green's functions and the superposition theorem, the EDD boundary integral equations for an arbitrarily shaped planar crack in an infinite 2D hexagonal QC body are established. Using the EDD boundary integral equation method, the asymptotic behavior along the crack front is studied and the classical singular index of 1/2 is obtained at the crack edge. The extended stress intensity factors are expressed in terms of the EDDs across crack surfaces. Finally, the energy release rate is obtained using the definitions of the stress intensity factors.     

A strip yield model solution for an internally cracked piezoelectric strip

An analysis of the crack closure and fatigue crack growth rate have been carried out for an infinitely long poled piezoelectric ceramic strip weakened by a straight hair line internal crack. The ceramic under consideration is assumed to be mechanically more brittle. The crack faces are perpendicular to the poled direction of the strip. The crack faces open in Mode-I deformation on account of in-plane tension applied to the edges of the strip together with either an in-plane electric displacement prescribed on edges of the strip or a uniform constant electric field prescribed on its edges. As a result, a yield zone is formed ahead of each tip of the crack. The yield zones developed are then arrested by applying a normal, cohesive, linearly varying yield point-stress to their rims. For each case, the Fourier transform method is used to find a solution. The resulting integral equations are solved numerically. Expressions are derived for the crack opening displacement and the crack growth rate. The variations in these quantities are plotted in relation to the affecting parameters, viz., the strip thickness, the yield zone length, the electric displacement, and material constants. A case study is presented graphically for PZT-4, PZT-5H, and BaTiO₃ ceramics. Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 44, No. 5, pp. 647–664, September–October, 2008.     

Crack growth modeling via 3D automatic adaptive mesh refinement based on modified-SPR technique

In this paper, the three-dimensional automatic adaptive mesh refinement is presented in modeling the crack propagation based on the modified superconvergent patch recovery technique. The technique is developed for the mixed mode fracture analysis of different fracture specimens. The stress intensity factors are calculated at the crack tip region and the crack propagation is determined by applying a proper crack growth criterion. An automatic adaptive mesh refinement is employed on the basis of modified superconvergent patch recovery (MSPR) technique to simulate the crack growth by applying the asymptotic crack tip solution and using the collapsed quarter-point singular tetrahedral elements at the crack tip region. A-posteriori error estimator is used based on the Zienkiewicz–Zhu method to estimate the error of fracture parameters and predict the crack path pattern. Finally, the efficiency and accuracy of proposed computational algorithm is demonstrated by several numerical examples.     

Transverse cracks in cross-ply laminates. 1. Stress analysis

During service loading of cross-ply laminates, transverse cracks occur in plies. The cracks parallel to the fiber direction are extended over the full thickness of transverse plies and often cross the entire test specimen width. It is widely recognized that the changes of laminate thermomechanical constants, caused by the transverse cracking of composite laminates, can be significant. Theoretical stress analysis in the cross-ply laminates in the vicinity of cracks is performed using numerical (FE) and analytical methods. The effect of transverse cracks on the degradation of elastic properties will be discussed in Part 2 [1]. Approximate analytical micromechanical models based on shear lag predictions, variational analysis, and numerical 2D finite element calculations were verified in their predictive abilities. The three variational models used are based on the principle of minimum complementary energy and have different degrees of accuracy with respect to the stress assumptions used (Hashin's, 2D 0° and 2D 0°/90° models). Using FEM, the plane stress and strain state were analyzed. The effect of material properties and layer thickness on the stress distribution in a 90° layer was evaluated by varying the crack spacing. The crack opening displacement (COD), normalized with respect to the far field strain, is proposed as a measure of reduction of the mechanical properties. Since the CODs are rather insensitive to the crack spacing (crack density) in a wide region, they will be used in modeling the stiffness reduction in these laminates [1].Translated from Mekhanika Kompozitnykh Materialov, Vol. 33, No. 6, pp. 796–820, November–December, 1997.