Semi-permeable crack analysis in a 2D magnetoelectroelastic medium based on the EMPS model

For a crack in a magnetoelectroelastic plane under the electrically and magnetically semi-permeable boundary condition, we derive the non-linear analytical solution of the strip electric–magnetic polarization saturation (EMPS) model. Using the extended dislocation theory and integral equation method, we obtain the electric and magnetic yielding zones, as well as the field intensity factor and local J-integral. Adapting an iterative method, numerical examples were performed to analyze the effect of different boundary conditions and the electric–magnetic saturated properties on the electric displacement and magnetic induction in the crack cavity, electric and magnetic yielding zones, stress intensity factor and local J-integral.     

Anti-plane analysis of semi-infinite crack in piezoelectric strip

Using the complex variable function method and the technique of the conformal mapping, the fracture problem of a semi-infinite crack in a piezoelectric strip is studied under the anti-plane shear stress and the in-plane electric load. The analytic solutions of the field intensity factors and the mechanical strain energy release rate are presented under the assumption that the surface of the crack is electrically impermeable. When the height of the strip tends to infinity, the analytic solutions of an infinitely large piezoelectric solid with a semi-infinite crack are obtained. Moreover, the present results can be reduced to the well-known solutions for a purely elastic material in the absence of the electric loading. In addition, numerical examples are given to show the influences of the loaded crack length, the height of the strip, and the applied mechanical/electric loads on the mechanical strain energy release rate.     

Effect of Maxwell stress on a moving crack with polarization saturation region in ferroelectric solid

The focus of this work is on a generalized two-dimensional problem of a crack moving in a piezoelectric solid subjected to uniform electrical load at infinity. The novel point includes that the electric field inside the crack is taken into account when polarization saturation region exists. Based on the extended Stroh formalism and complex function method, explicit expressions of both the stress fields in the solid and electric fields inside the crack are derived by using semi-permeable crack model, respectively. Effect of Maxwell stress along the crack surface is investigated and the results are illustrated graphically. It is shown that the moving speed of the crack cannot exceed the lowest bulk wave speed. It is also found that the medium properties inside the crack and surrounding the ferroelectric solid at infinity directly affect the Maxwell stress, and as a result the Maxwell stresses are remarkable and cannot be ignored under different electric load.     

Anti-plane analysis on a finite crack in a one-dimensional hexagonal quasicrystal strip

The anti-plane fracture problem for a finite crack in a one-dimensional hexagonal quasicrystal strip is analyzed. By using Fourier transforms, the mixed boundary value problems are reduced to the dual integral equations. The solution of the dual integral equations is then expressed by the complete elliptic integrals of the first and the third kinds. The expressions for stress, strains, displacements and field intensity factors of the phonon and phason fields near the crack tip are obtained exactly. The path-independent integral derived by a conservation law equals the energy release rate, which can be used as a fracture criterion for a mode III fracture problem.     

Geometrical and interfacial non-linearities in the analysis of delamination in composites

The subject of this paper is the study of interactions between delamination and geometrical nonlinearities. This problem generally addressed by means of Fracture Mechanics, is here treated by modelling interfacial degradation through irreversible, softening interface constitutive laws. Geometrical nonlinearities are associated to large displacements of thin layers which are assumed to be hyperelastic. A general formulation is presented and then specialised to the case of rectilinear laminated beams. Numerical aspects concerning finite element implementation are discussed together with a local control algorithm for the treatment of unstable paths. Numerical examples showing the effect of interaction between geometrical and material nonlinearities are presented.     

Evaluation of the Lazarus-Leblond constants in the asymptotic model of the interfacial wavy crack

The paper addresses the problem of a semi-infinite plane crack along the interface between two isotropic half-spaces. Two methods of solution have been considered in the past: Lazarus and Leblond [1998a. Three-dimensional crack-face weight functions for the semi-infinite interface crack-I: variation of the stress intensity factors due to some small perturbation of the crack front. J. Mech. Phys. Solids 46, 489-511, 1998b. Three-dimensional crack-face weight functions for the semi-infinite interface crack-II: integrodifferential equations on the weight functions and resolution J. Mech. Phys. Solids 46, 513-536] applied the “special” method by Bueckner [1987. Weight functions and fundamental fields for the penny-shaped and the half-plane crack in three space. Int. J. Solids Struct. 23, 57-93] and found the expression of the variation of the stress intensity factors for a wavy crack without solving the complete elasticity problem; their solution is expressed in terms of the physical variables, and it involves five constants whose analytical representation was unknown; on the other hand, the “general” solution to the problem has been recently addressed by Bercial-Velez et al. [2005. High-order asymptotics and perturbation problems for 3D interfacial cracks. J. Mech. Phys. Solids 53, 1128-1162], using a Wiener-Hopf analysis and singular asymptotics near the crack front.The main goal of the present paper is to complete the solution to the problem by providing the connection between the two methods. This is done by constructing an integral representation for Lazarus-Leblond's weight functions and by deriving the closed form representations of Lazarus-Leblond's constants.     

Wave diffraction-induced crack/cohesive strip growth based on rudimentary non-linear models

The reaction of a Mode I 2D semi-infinite crack to the diffraction of a plane tension wave in an unbounded solid is treated. A vanishingly thin cohesive zone, or strip, represents rudimentary crack edge yield, and two generic response laws that incorporate rudimentary non-linear effects are studied. Simultaneous crack/zone growth at constant rates is treated, with pure zone growth recovered as a limit case.Both models involve multiple solution sets that can be obtained analytically and studied. Attention focuses on those associated with uniform strip fields.Results based on parameters for 2024-T4 aluminum and 4130-HT steel suggest that the two models involve tension loads critical for crack/zone growth lower than standard model predictions, especially for high growth rates. The results also suggest that dynamic crack/zone growth may be more sensitive to multiaxiality near the crack edge, and proceed at lower loads. All of these results are general in the sense that they are obtained without specifying a complete set of growth criteria.     

Fracture analysis on the arc-shaped interfacial crack between a homogeneous cylinder and its coating

The mechanical model is established for a mode I arc-shaped crack along the interface between a functionally graded layer and a homogeneous cylinder. The displacement and stress fields are expressed in the form of infinite series by the method of variable separation. By introducing two dislocation density functions, the mixed boundary conditions of the crack are recast into a system of Cauchy singular integral equations, which are then solved numerically. Parametric studies on the stress intensity factors (SIFs) yield two aspects of guidelines for engineering designs. (a) The SIFs depend on the ratio between the outer and inner radii rather than their respective values, and 1.5 seems to be an optimal value for this ratio. (b) The outer graded layer should be stiffer than the inner homogeneous cylinder and large non-homogeneity parameter is beneficial to SIF reduction.     

Numerical analysis of crack development in structurally nonuniform composites

Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, No. 6, pp. 115–121, November–December, 1992.     

Fatigue crack opening stress based on the strip-yield model

The modified strip-yield model based on the Dugdale model and two-dimensional approximate weight function method were utilized to evaluate the effect of in-plane constraint, transverse stress, on the fatigue crack closure. The plastic zone sizes and the crack opening stresses considering transverse stress were calculated for four specimens: single edge-notched tension (SENT) specimen, single edge-notched bend (SENB) specimen, center-cracked tension (CCT) specimen, double edge-notched tension (DENT) specimen under uniaxial loading. And the crack opening behavior of the center-cracked specimen under biaxial loading was also evaluated. Normalized crack opening stresses σ_op/σ_max for four specimens were successfully described by the normalized plastic zone parameter Δω^′_rev/ω^′ considering transverse stress, where Δω^′_rev and ω^′ are the size of the reversed plastic zone at the moment of first crack tip closure and the size of the forward plastic zone for maximum stress, respectively. The normalized plastic zone parameter with transverse stress also was satisfactorily correlated with the behavior of crack closure for CCT specimen under biaxial loading.     

Experimental and numerical analysis of fretting crack formation based on 3D X-FEM frictional contact fatigue crack model

Nowadays, numerical simulation of 3D fatigue crack growth is easily handled using the eXtended Finite Element Method coupled with level set techniques. The finite element mesh does not need to conform to the crack geometry. Most difficulties associated to complex mesh generation around the crack and the re-meshing steps during the possible propagation are hence avoided. A 3D two-scale frictional contact fatigue crack model developed within the X-FEM framework is presented in this article. It allows the use of a refined discretization of the crack interface independent from the underlying finite element mesh and adapted to the frictional contact crack scale. A stabilized three-field weak formulation is also proposed to avoid possible oscillations in the local solution linked to the LBB condition when tangential slip is occurring. Two basic three-dimensional numerical examples are presented. They aim at illustrating the capacities and the high level of accuracy of the proposed X-FEM model. Stress intensity factors are computed along the crack front. Finally an experimental 3D ball/plate fretting fatigue test with running conditions inducing crack nucleation and propagation is modeled. 3D crack shapes defined from actual experimental ones and fretting loading cycle are considered. This latter numerical simulation demonstrates the model ability to deal with challenging actual complex problems and the possibility to achieve tribological fatigue prediction at a design stage based on the fatigue crack modeling.     

Complex variable approach in studying modified polarization saturation model in two-dimensional semipermeable piezoelectric media

A modified polarization saturation model is proposed and addressed mathematically using a complex variable approach in two-dimensional(2 D) semipermeable piezoelectric media. In this model, an existing polarization saturation(PS) model in 2D piezoelectric media is modified by considering a linearly varying saturated normal electric displacement load in place of a constant normal electric displacement load, applied on a saturated electric zone. A centre cracked infinite 2D piezoelectric domain subject to an arbitrary poling direction and in-plane electromechanical loadings is considered for the analytical and numerical studies. Here, the problem is mathematically modeled as a non-homogeneous Riemann-Hilbert problem in terms of unknown complex potential functions representing electric displacement and stress components. Having solved the Hilbert problem, the solutions to the saturated zone length, the crack opening displacement(COD), the crack opening potential(COP), and the local stress intensity factors(SIFs) are obtained in explicit forms. A numerical study is also presented for the proposed modified model, showing the effects of the saturation condition on the applied electrical loading, the saturation zone length, and the COP. The results of fracture parameters obtained from the proposed model are compared with the existing PS model subject to electrical loading, crack face conditions, and polarization angles.     

Effects of specimen thickness, hardening and crack closure for the plastic strip model

The triaxial stress constraint T_z and the effective yield stress distribution in the plastic zone for strain hardening materials are considered for the strip where the thickness effect is investigated by introducing a plastic constraint factor α. This factor depends on the specimen thickness, crack length, load level and hardening exponent. A simple expression of the plastic zone length and an expression involving α are obtained. Application of the strip model to crack closure shows that the specimen thickness has an influence and the results are compared with those found finite element.     

一种基于SERR法的拟合裂纹呼吸模型

根据断裂力学理论,采用SERR法对含横向裂纹转轴的裂纹呼吸行为进行了研究.采用函数逼近的方法,将转轴旋转过程中裂纹的呼吸行为对刚度的影响规律拟合为一个周期函数,提出了一种描述裂纹呼吸行为的拟合模型,并讨论了裂纹深度变化对呼吸行为的影响.结果表明:拟合的裂纹呼吸模型考虑了裂纹呼吸行为的非线性特征与有限元以及实验结果更接近;在a/R<0.5情况下裂纹深度的变化对裂纹呼吸行为的变化规律影响较小,动力学仿真研究时可以忽略这种影响,只考虑刚度的变化.     

Numerical solution of polarization saturation/dielectric breakdown model in 2D finite piezoelectric media

The polarization saturation (PS) model [Gao, H., Barnett, D.M., 1996. An invariance property of local energy release rates in a strip saturation model of piezoelectric fracture. Int. J. Fract. 79, R25–R29; Gao, H., Zhang, T.Y., Tong, P., 1997. Local and global energy release rates for an electrically yielded crack in a piezoelectric ceramic. J. Mech. Phys. Solids 45, 491–510], and the dielectric breakdown (DB) model [Zhang, T.Y., Zhao, M.H., Cao, C.F., 2005. The strip dielectric breakdown model. Int. J. Fract. 132, 311–327] explain very well some experimental observations of fracture of piezoelectric ceramics. In this paper, the nonlinear hybrid extended displacement discontinuity-fundamental solution method (NLHEDD-FSM) is presented for numerical analysis of both the PS and DB models of two-dimensional (2D) finite piezoelectric media under impermeable and semi-permeable electric boundary conditions. In this NLHEDD-FSM, the solution is expressed approximately by a linear combination of fundamental solutions of the governing equations, which includes the extended point force fundamental solutions with sources placed at chosen points outside the domain of the problem under consideration, and the extended Crouch fundamental solutions with extended displacement discontinuities placed on the crack and the electric yielding zone. The coefficients of the fundamental solutions are determined by letting the approximated solution satisfy certain conditions on the boundary of the domain, on the crack face and the electric yielding zone. The zero electric displacement intensity factor in the PS model or the zero electric field strength intensity factor in the DB model at the outer tips of the electric yielding zone is used as a supplementary condition to determine the size of the electric yielding zone. Iteration approaches are adopted in the NLHEDD-FSM. The electric yielding zone is determined, and the extended intensity factors and the local J-integral are calculated for center cracks in piezoelectric strips. The effects of finite domain size, saturation property and different electric boundary conditions, as well as different models on the electric yielding zone and the local J-integral, are studied.     

Boundary element weight function analysis of a strip yield crack in a rotating disk

The boundary element method combined with subtration of Bueckner singular fields are used to obtain weight functions for an internal edge crack in a rotating annular disk. A previously developed, general representation of the weight function is used which leads to integrals that can be evaluated analytically to obtain the stress intensity factor and surface displacements of the crack. The determination of crack tip opening displacements for the strip yield crack is reduced to a non-singular integral which can be evaluated in closed form. The strip yield zone length and crack tip opening displacement are determined for an internal radial crack in a rotating annular disk for a range of crack lengths and rotational speeds.     

A rate-dependent damage model for brittle materials based on the dominant crack

A rate-dependent, continuum damage model is developed for brittle materials under dynamic loading. This model improves on the approach (ISOSCM) of [Addessio, F.L., Johnson, J.N., 1990. A constitutive model for the dynamic response of brittle materials. Journal of Applied Physics 67, 3275–3286] in several respects. (1) A new damage surface is found by applying the generalized Griffith instability criterion to the dominant crack (having the most unstable orientation), rather than by averaging the instability condition over all crack orientations as done previously. The new surface removes a discontinuity in the damage surface in ISOSCM when the pressure changes sign. (2) The strain due to crack opening is more consistent with crack mechanics, with only the tensile principal stresses contributing to the crack opening strain. This is achieved by incorporating a projection operator in the equation for the crack opening strain. One consequence of incorporating the projection operator is a prediction of shear dilatancy, which is not accounted for in ISOSCM. (3) The evolution of damage, which is based on the energy-release rate for the dominant crack, has a physical basis, whereas in the previous approach the damage growth rate was assumed to be an exponential function of the distance from the stress state to the damage surface without specific physical justification.An implicit algorithm has been developed so that a larger time step can be used than with the explicit algorithm used in ISOSCM. The numerical results of a silicon carbide (SiC) ceramic under several loading paths (hydrostatic tension/compression, uniaxial strain, uniaxial stress, and shear) and strain rates are presented to illustrate the main features of the model.     

一种考虑界面相物性影响的界面单元法

界面是由复杂的界面相简化而成的,界面破坏实际是界面相材料的破坏。数值计算为了方便,如经典模型和内聚力模型等,都把很薄的界面相作无厚度化处理。导致只能考虑界面的面力,而无法考虑界面相内的应力(平行于界面方向的应力)。使界面失效准则先天性地排除了界面相内部应力的影响,从界面相材料失效机理的角度来看这是不够严谨的。本文将界面相材料等效为一种弹性连续体,由界面本构关系推导得到了一种新的界面单元。该单元具有界面参数易确定、对界面相物性可以进行等效描述等优点。通过商用有限元软件ABAQUS和用户子程序UEL实现了数值分析,并与直接物理模型的数值模拟结果进行对比,证明了本方法的简便及准确性。通过对不同界面相厚度结构的进一步分析,探讨了本文方法的可行范围。