The asymptotic field of mode Ⅰ quasi-static crack growth on the interface between a rigid and a pressure-sensitive material

A mechanical model of the quasi-static interface of a mode I crack between a rigid and a pressure-sensitive viscoelastic material was established to investigate the mechanical characteristic of ship-building engineering hi-materials. In the stable growth stage, stress and strain have the same singularity, ie （σ, ε） ∝ r^-1/（n-1）. The variable-separable asymptotic solutions of stress and strain at the crack tip were obtained by adopting Airy＇s stress function and the numerical results of stress and strain in the crack-tip field were obtained by the shooting method. The results showed that the near-tip fields are mainly governed by the power-hardening exponent n and the Poisson ratio v of the pressure-sensitive material. The fracture criterion of mode I quasi-static crack growth in pressure-sensitive materials, according to the asymptotic analyses of the crack-tip field, can be viewed from the perspective of strain.

Interface crack problems for mode Ⅱ of double dissimilar orthotropic composite materials

The fracture problems near the similar orthotropic composite materials are interface crack tip for mode Ⅱ of double disstudied. The mechanical models of interface crack for mode Ⅱ are given. By translating the governing equations into the generalized hi-harmonic equations, the stress functions containing two stress singularity exponents are derived with the help of a complex function method. Based on the boundary conditions, a system of non-homogeneous linear equations is found. Two real stress singularity exponents are determined be solving this system under appropriate conditions about bimaterial engineering parameters. According to the uniqueness theorem of limit, both the formulae of stress intensity factors and theoretical solutions of stress field near the interface crack tip are derived. When the two orthotropic materials are the same, the stress singularity exponents, stress intensity factors and stresses for mode II crack of the orthotropic single material are obtained.

Asymmetrical dynamic propagation problems on mode III interface crack

By the application of the theory of complex functions, asymmetrical dynamic propagation problems on mode Ⅲ interface crack are studied. The universal representations of analytical solutions are obtained by the approaches of self-similar function. The problems researched can be facilely transformed into Riemann-Hilbert problems and analytical solution to an asymmetrical propagation crack under the condition of point loads and unit-step loads, respectively, is acquired. After those solutions were used by superposition theorem, the solutions of arbitrarily complex problems could be attained.

Crack-tip field on mode Ⅱ interface crack of double dissimilar orthotropic composite materials

Two systems of non-homogeneous linear equations with 8 unknowns are obtained.This is done by introducing two stress functions containing 16 undetermined coefficients and two real stress singularity exponents with the help of boundary conditions.By solving the above systems of non-homogeneous linear equations,the two real stress singularity exponents can be determined when the double material parameters meet certain conditions.The expression of the stress function and all coefficients are obtained based on the uniqueness theorem of limit.By substituting these parameters into the corresponding mechanics equations,theoretical solutions to the stress intensity factor,the stress field and the displacement field near the crack tip of each material can be obtained when both discriminants of the characteristic equations are less than zero.Stress and displacement near the crack tip show mixed crack characteristics without stress oscillation and crack surface overlapping.As an example,when the two orthotropic materials are the same,the stress singularity exponent,the stress intensity factor,and expressions for the stress and the displacement fields of the orthotropic single materials can be derived.

Elastic-viscoplastic field at mixed-mode interface crack-tip under compression and shear

For a compression-shear mixed mode interface crack, it is difficult to solve the stress and strain fields considering the material viscosity, the crack-tip singularity, the frictional effect, and the mixed loading level. In this paper, a mechanical model of the dynamic propagation interface crack for the compression-shear mixed mode is proposed using an elastic-viscoplastic constitutive model. The governing equations of propagation crack interface at the crack-tip are given. The numerical analysis is performed for the interface crack of the compression-shear mixed mode by introducing a displacement function and some boundary conditions. The distributed regularities of stress field of the interface crack-tip are discussed with several special parameters. The final results show that the viscosity effect and the frictional contact effect on the crack surface and the mixed-load parameter are important factors in studying the mixed mode interface crack- tip fields. These fields are controlled by the viscosity coefficient, the Mach number, and the singularity exponent.

Finite Element Modeling of Crack Tip Blunting for Estimation of Energy Release Rate Component of Mode I Crack near a Strength Mismatched Interface

The paper presents finite element modeling of crack tip blunting for numerical estimation of fracture parameter of a Mode I crack, in weak alloy steel, which is near and normal to the interface of elastically identical but stronger maraging steel. The bimetallic body is subjected to monotonic load in K dominated regime. Crack tip yield zone across the interface, treated as Dugdale’s cohesive zone, is isolated from the bimetallic domain and is modeled alone under the action of respective cohesive stress over yield zones for obtaining the contribution of mismatch between yield strength of the steels in crack energy release rate component, Jinterface. Effect of far field load on Jinterface is found separately from a theoretical model. Numerical and theoretical results of Jinterface are in good agreement.

Fracture analysis of mode-II crack perpendicular to imperfect bimaterial interface

The problem of a mode-II crack interface of two bonded dissimilar materials close to and perpendicular to an imperfect is investigated. The imperfect interface is modelled by a linear spring with the vanishing thickness. The Fourier transform is used to solve the boundary-value problem and to derive a singular integral equation with the Cauchy kernel. The stress intensity factors near the left and right crack tips are evaluated by numerically solving the resulting equation. SeverM special cases of the mode-II crack problem with an imperfect interface are studied in detail. The effects of the interfacial imperfection on the stress intensity factors for a bimaterial system of aluminum and steel are shown graphically. The obtained observation reveals that the stress intensity factors are dependent on the interface parameters and vary between those with a fully debonded interface and those with a perfect interface.

Anti-Plane Elasticity Problem and Mode Ⅲ Crack Problem of Cubic Quasicrystal

The fracture theory of cubic quasicrystal was developed. The exact analytic solution of a Mode Ⅲ Griffith crack in the material was obtained by using the Fourier transform and dual integral equations theory, and so the displacement and stress fields, the stress intensity factor and strain energy release rate were determined. The results show that the stress intensity factor is independent of material constants, and the strain energy release rate is dependent on all material constants. These provide important information for studying the deformation and fracture of the new solid material.

Evaluation of stress intensity factors for bi-material interface cracks using displacement jump methods

The aim of the present work is to investigate the numerical modeling of interfacial cracks that may appear at the interface between two isotropic elastic materials. The extended finite element method is employed to analyze brittle and bi-material interfacial fatigue crack growth by computing the mixed mode stress intensity factors(SIF). Three different approaches are introduced to compute the SIFs. In the first one, mixed mode SIF is deduced from the computation of the contour integral as per the classical J-integral method,whereas a displacement method is used to evaluate the SIF by using either one or two displacement jumps located along the crack path in the second and third approaches. The displacement jump method is rather classical for mono-materials,but has to our knowledge not been used up to now for a bimaterial. Hence, use of displacement jump for characterizing bi-material cracks constitutes the main contribution of the present study. Several benchmark tests including parametric studies are performed to show the effectiveness of these computational methodologies for SIF considering static and fatigue problems of bi-material structures. It is found that results based on the displacement jump methods are in a very good agreement with those of exact solutions, such as for the J-integral method, but with a larger domain of applicability and a better numerical efficiency(less time consuming and less spurious boundary effect).

A PENNY-SHAPED CRACK IN A TRANSVERSELY ISOTROPIC PIEZOELECTRIC SOLID:MODES Ⅱ AND Ⅲ PROBLEMS

An exact analysis of the modes Ⅱ and Ⅲ problems of a penny- shaped crack in a transversely isotropic piezoelectric medium is performed in this paper.The potential theory method is employed based on the general solution of three-dimensional piezoelasticity and the four harmonics involved are represented by one complex potential.Previous results in potential theory are then utilized to obtain the exact solution that is expressed in terms of elementary functions.Comparison is made between the current results with those published and good agreement is obtained.

STRESS INTENSITY FACTOR OF AN ANTI-PLANE CRACK PARALLEL TO THE WEAK/MICRO-DISCONTINUOUS INTERFACE IN A BI-FGM COMPOSITE

The problem considered is a mode Ⅲ crack lying parallel to the interface of an exponential-type functional graded material （FGM） strip bonded to a linear-type FGM substrate with infinite thickness. By applying the Fourier integral transform, the problem was reduced as a Cauchy singular integral equation with an unknown dislocation density function. The collocation method based on Chebyshev polynomials proposed by Erdogan and Gupta was used to solve the singular integral equation numerically. With the numerical solution, the effects of the geometrical and physical parameters on the stress intensity factor （SIF） were analyzed and the following conclusions were drawn： （a） The region affected by the interface or free surface varies with the material rigidity, and higher material rigidity will lead to bigger affected region. （b） The SIF of the crack in the affected region and parallel to the micro-discontinuous interface is lower than those of the weak discontinuous cases. Reducing the weak-discontinuity of the interface will be beneficial to decrease the SIF of the interface-parallel crack in the region affected by the interface. （c） The effect of the free surface on SIF is more remarkable than that of the interface, and the latter is still more notable than that of the material rigidity. When the effects of the interface and free surface are fixed, increase of the material rigidity will enhance the value of SIF.

General forms of elastic-plastic matching equations for mode-Ⅲ cracks near crack line

Crack line analysis is an effective way to solve elastic-plastic crack problems. Application of the method does not need the traditional small-scale yielding conditions and can obtain sufficiently accurate solutions near the crack line. To address mode- Ⅲ crack problems under the perfect elastic-plastic condition, matching procedures of the crack line analysis method axe summarized and refined to give general forms and formulation steps of plastic field, elastic-plastic boundary, and elastic-plastic matching equations near the crack line. The research unifies mode-III crack problems under different conditions into a problem of determining four integral constants with four matching equations. An example is given to verify correctness, conciseness, and generality of the procedure.

Analysis of composite material interface crack face contact and friction effects using a new node-pairs contact algorithm

A new node-pairs contact algorithm is proposed to deal with a composite material or bi-material interface crack face contact and friction problem （e.g., resistant coating and thermal barrier coatings） subjected to complicated load conditions. To decrease the calculation scale and calculation errors, the local Lagrange multipliers are solved only on a pair of contact nodes using the Jacobi iteration method, and the constraint modification of the tangential multipliers are required. After the calculation of the present node-pairs Lagrange multiplier, it is turned to next contact node-pairs until all node-pairs have finished. Compared with an ordinary contact algorithm, the new local node-pairs contact algorithm is allowed a more precise element on the contact face without the stiffness matrix singularity. The stress intensity factors （SIFs） and the contact region of an infinite plate central crack are calculated and show good agreement with those in the literature. The contact zone near the crack tip as well as its influence on singularity of stress fields are studied. Furthermore, the frictional contacts are also considered and found to have a significant influence on the SIFs. The normalized mode-II stress intensity factors KII for the friction coefficient decrease by 16% when f changes from 1 to 0.

Dislocation distribution model of mode Ⅲ propagation crack

Dislocation distribution functions of the edges of mode Ⅲ propagation crack subjected to three types of loads were studied by the methods of the theory of complex variable functions,by which,the problems researched were facilely transformed into Riemann-Hilbert problems and Keldysh-Sedov mixed boundary value problems. Analytical solutions of stresses,displacement and dynamic stress intensity factor were obtained by the measures of the theory of self-similar functions and corresponding differential and integral operation. In terms of the relationship between dislocation distribution functions and displacements,analytical solutions of dislocation distribution functions were obttained,and variational rules of dislocation distribution functions were depicted.

Field structure at mode Ⅲ dynamically propagating crack tip in elastic-viscoplastic materials

An elastic-viscoplastic mechanics model is used to investigate asymptotically the mode Ⅲ dynamically propagating crack tip field in elastic-viscoplastic materials. The stress and strain fields at the crack tip possess the same power-law singularity under a linear-hardening condition. The singularity exponent is uniquely determined by the viscosity coefficient of the material. Numerical results indicate that the motion parameter of the crack propagating speed has little effect on the zone structure at the crack tip. The hardening coefficient dominates the structure of the crack-tip field. However, the secondary plastic zone has little influence on the field. The viscosity of the material dominates the strength of stress and strain fields at the crack tip while it does have certain influence on the crack-tip field structure. The dynamic crack-tip field degenerates into the relevant quasi-static solution when the crack moving speed is zero. The corresponding perfectly-plastic solution is recovered from the linear-hardening solution when the hardening coefficient becomes zero.

Effects of interface adhesive and resin cohesive strength on tensile strength of resin sand based on numerical analysis

A plane-strain unit-cell finite element model was proposed to study the effects of resin/sand interface adhesive and resin cohesive strength on the overall tensile strength of resin sand,as well as the fracture modes.The main micro-scale characteristics of the numerical model were extracted from the micrograph of resin sand specimens by three-dimensional X-ray microscopy(3 D-XRM).The extended finite element method(XFEM)and cohesive behavior method were employed to explicitly describe the resin fracture and sand/resin interface debonding,separately.The corresponding experimental observation of micro-scale failure behavior based on the scanning electron microscopy(SEM)was presented for a comparison.The numerical results show that the initial failure of the model occurs at the sand/resin interface,followed by consequent resin failure.Dependent on the resin cohesive strength,the location of resin failure varies from the central zones to resin neck arc zones.A typical mixed mode fracture is observed,which is consistent with the corresponding micro-scale experimental observation.When the resin cohesive strength ranges between 8 and 12 MPa,the resin cracks occur at the central zone of resin bridges and propagate perpendicularly to the tensile direction until through cracks happen.At a higher range(between 12 and 16 MPa),interface cracks cross with resin cracks,bonding bridges of resin sand are broken.The interface adhesive strength has a more significant effect on the overall tensile strength of resin sand than the resin cohesive strength.The overall tensile strength of resin sand increases first then keeps stable with the increase of the resin cohesive strength.This work attempts to establish a numerical model which accurately describes the complicated mixed mode fracture of resin sand,which is beneficial to understand deeply the fracture mechanism of resin sand.

Dynamic propagation problems concerning surfaces of asymmetrical mode III crack subjected to moving loads

With the theory of complex functions, dynamic propagation problems concerning surfaces of asymmetrical mode Ⅲ crack subjected to moving loads are investigated. General representations of analytical solutions are obtained with self-similar functions. The problems can be easily converted into Riemann-Hilbert problems using this technique. Analytical solutions to stress, displacement and dynamic stress intensity factor under constant and unit-step moving loads on the surfaces of asymmetrical extension crack, respectively, are obtained. By applying these solutions, together with the superposition principle, solutions of discretionarily intricate problems can be found.

混凝土Ⅲ型裂缝动静力损伤断裂分析

本文把混凝土的动静力损伤本构模型通过应变等价性假设引入到线弹性断裂力学中，推出了混凝土Ⅲ型裂缝在动静力荷载作用下缝端附近的损伤场，给出了相应的判据，实现了损伤与断裂的耦合分析，为把损伤理论应用到混凝土断裂分析中提供了一条可行的途径。文末还给出了混凝土Ⅲ型裂缝允许损伤尺度的估计值。

Ⅲ型裂纹弹塑性场在裂纹线附近匹配方程的一般形式

针对理想弹塑性Ⅲ型裂纹问题,对线场分析方法的步骤和匹配过程进行了凝练和归纳,给出了裂纹线附近塑性场、弹塑性边界、弹塑性匹配方程的一般形式及其一般求解步骤,将不同条件下的Ⅲ型裂纹问题归结为由4个匹配方程确定4个待定常数,并通过一个具体问题,验证了这一方法的正确性、简明性和通用性.

Ⅲ型界面裂纹的非对称动态扩展问题

采用复变函数论的方法,对Ⅲ型界面裂纹的非对称动态扩展问题进行了研究.通过自相似函数的方法可以获得解析解的一般表达式.应用该法可以迅速地将所论问题转化为Riemann-Hilbert问题,并求得了非对称扩展裂纹分别在集中载荷、阶跃载荷作用下的解析解.利用这些解并采用叠加原理,就可以求得任意复杂问题的解.