Boundary Element Analysis forModeⅢCrack Problems of Thin-Walled Structures from Micro-to Nano-Scales

This paper develops a new numerical framework for modeⅢcrack problems of thin-walled structures by integrating multiple advanced techniques in the boundary element literature.The details of special crack-tip elements for displacement and stress are derived.An exponential transformation technique is introduced to accurately calculate the nearly singular integral,which is the key task of the boundary element simulation of thin-walled structures.Three numerical experiments with different types of cracks are provided to verify the performance of the present numerical framework.Numerical results demonstrate that the present scheme is valid for modeⅢcrack problems of thin-walled structures with the thickness-to-length ratio in the microscale,even nanoscale,regime.

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.

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.

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.

MESHLESS METHOD FOR 2D MIXED-MODE CRACK PROPAGATION BASED ON VORONOI CELL

A meshless method integrated with linear elastic fracture mechanics(LEFM)is presented for 2D mixed-mode crack propagation analysis.The domain is divided automatically into sub-domains based on Voronoi cells,which are used for quadrature for the potential energy. The continuous crack propagation is simulated with an incremental crack-extension method which assumes a piecewise linear discretization of the unknown crack path.For each increment of the crack extension,the meshless method is applied to carry out a stress analysis of the cracked structure.The J-integral,which can be decomposed into mode Ⅰ and mode Ⅱ for mixed-mode crack,is used for the evaluation of the stress intensity factors(SIFs).The crack-propagation direction,predicted on an incremental basis, is computed by a criterion defined in terms of the SIFs. The flowchart of the proposed procedure is presented and two numerical problems are analyzed with this method.The meshless results agree well with the experimental ones,which validates the accuracy and efficiency of the method.

EXACT SOLUTIONS OF NEAR CRACK LINE FIELDS FOR MODE I CRACK UNDER PLANE STRESS CONDITION IN AN ELASTIC－PERFECTLY PLASTIC SOLID

The near crack line analysis method has been used in the present paper,The classical small scale yielding conditions have been completely abandoned in the analyses and one inappropriate matching condition used to be used at the elasticplastic boundary has been corrected.The reasonable solution of the plastic stresses near the crack line region has been established.By matching the plastic stresses with the exact elastic stresses at the elastic-plastic boundary,the plastic stresses the length of the plastic zone and the unit normal vector of the elastic-plastic boundary near the crock line region have been obtained for a mode I crack under uniaxial tension,as well as a mode I crack under biaxial tension,which shows that for both conditions the plastic stress componentsσy, and σsy.he length of the plastic zone and the unit normal vector of the elastic-plastic boundary are quite the same while the plastic stress σs is different.

A HALF PLANE CRACK UNDER THREE-DIMENSIONAL COMBINED MODE IMPACT LOADING

The dynamic stress intensity factor history for a half plane crack in an otherwise unbounded elastic body, with the crack faces subjected to a traction distribution consisting of two pairs of suddenly-applied shear line loads is consid- ered. The analytic expression for the combined mode stress intensity factors as a function of time is obtained. The method of solution is based on the application of integral transforms and the Wiener-Hopf technique. Some features of the solutions are discussed and graphical numerical results are presented.

The effect of fractional thermoelasticity on a two-dimensional problem of a mode I crack in a rotating fiber-reinforced thermoelastic medium

This article is concerned with the effect of rotation on the general model of the equations of the generalized thermoe- lasticity for a homogeneous isotropic elastic half-space solid, whose surface is subjected to a Mode-I crack problem. The fractional order theory of thermoelasticity is used to obtain the analytical solutions for displacement components, force stresses, and temperature. The boundary of the crack is subjected to a prescribed stress distribution and temperature. The normal mode analysis technique is used to solve the resulting non-dimensional coupled governing equations of the problem. The variations of the considered variables with the horizontal distance are illustrated graphically. Some particular cases are also discussed in the context of the problem. Effects of the fractional parameter, reinforcement, and rotation on the varia- tions of different field quantities inside the elastic medium are analyzed graphically. Comparisons are made between the results in the presence and those in the absence of fiber-reinforcing, rotating and fractional parameters.

Simultaneously measuring initiation toughness of pure mode Ⅰ and mode Ⅱ cracks subjected to impact loads

In order to simultaneously measure the initiation toughness of pure mode Ⅰ and mode Ⅱ cracks in one specimen,a large-size double-cracked concave-convex plate(DCCP)specimen configuration was proposed.Impacting tests were implemented in the drop plate impact device.Strain gauges were employed to measure impact loads and crack initiation time.The corresponding numerical model was established by using the dynamic finite difference program AUTODYN,and the experimental-numerical method and ABAQUS code were utilized to obtain the initial fracture toughness of the crack.Using experiments and numerical research,we concluded that the DCCP specimen is suitable for measuring the initial fracture toughness of pure mode Ⅰ and mode Ⅱ cracks at the same time;the dynamic initiation toughness increases with the increase of loading rate and the crack initiation time decreases with increasing loading rate;the initiation toughness of mode Ⅱ crack is 0.5 times that of mode Ⅰ crack when subjected to the same loading rate.For the pre-crack in the vicinity of the bottom of a sample,when its length increases from 20 to 100 mm,the dynamic initiation toughness of the pure mode Ⅰ crack gradually decreases,and the longer the lower crack length is,the easier the crack would initiate,but the dynamic initiation toughness of pure mode Ⅱ crack varies little.

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.

FATIGUE GROWTH MODELING OF MIXED-MODE CRACK IN PLANE ELASTIC MEDIA

This paper presents an extension of a displacement discontinuity method with cracktip elements （a boundary element method） proposed by the author for fatigue crack growth analysis in plane elastic media under mixed-mode conditions. The boundary element method consists of the non-singular displacement discontinuity elements presented by Crouch and Starfield and the crack-tip displacement discontinuity elements due to the author. In the boundary element implementation the left or right crack-tip element is placed locally at the corresponding left or right crack tip on top of the non-singular displacement discontinuity elements that cover the entire crack surface and the other boundaries. Crack growth is simulated with an incremental crack extension analysis based on the maximum circumferential stress criterion. In the numerical simulation, for each increment of crack extension, remeshing of existing boundaries is not required because of an intrinsic feature of the numerical approach. Crack growth is modeled by adding new boundary elements on the incremental crack extension to the previous crack boundaries. At the same time, the element characteristics of some related elements are adjusted according to the manner in which the boundary element method is implemented. As an example, the fatigue growth process of cracks emanating from a circular hole in a plane elastic plate is simulated using the numerical simulation approach.

A MODE Ⅱ CRACK IN A TWO-DIMENSIONAL OCTAGONAL QUASICRYSTALS

The present study develops the fracture theory for a two-dimensional octagonal quasicrystals. The exact analytic solution of a Mode Ⅱ Griffith crack in the material was obtained by using the Fourier transform and dual integral equations theory, then the displacement and stress fields, stress intensity factor and strain energy release rate were determined, the physical sense of the results relative to phason and the difference between mechanical behaviors of the crack problem in crystal and quasicrystal were figured out. These provide important information for studying the deformation and fracture of the new solid phase.

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.

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.

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.

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.

Microscopic study of mode I crack tip deformation

By simulating edge dislocation emissions from a mode I crack tip along multiple inclined slip planes, the plastic zone and dislocation-free zone around the crack tip are obtained. It is found that the shape of the mode I plastic zone consists of two leaning forward loops which is better agreement with experimental observations. Except at the crack tip there are also stress peaks in front of the crack tip. A formula of the maximum peak stress as a function of the applied stress intensity factor and the friction stress has been regressed.

ASYMPTOTIC ANALYSIS OF MODE Ⅱ STATIONARY GROWTH CRACK ON ELASTIC-ELASTIC POWER LAW CREEPING BIMATERIAL INTERFACE

A mechanical model was established for modeⅡinterfacial crack static growing along an elastic_elastic power law creeping bimaterial interface. For two kinds of boundary conditions on crack faces, traction free and frictional contact, asymptotic solutions of the stress and strain near tip_crack were given. Results derived indicate that the stress and strain have the same singularity, there is not the oscillatory singularity in the field; the creep power_hardening index n and the ratio of Young's module notably influence the crack_tip field in region of elastic power law creeping material and n only influences distribution of stresses and strains in region of elastic material. When n is bigger, the creeping deformation is dominant and stress fields become steady,which does not change with n. Poisson's ratio does not affect the distributing of the crack_tip field.

In-plane shear(ModeⅡ) crack sub-critical propagation of rock at high temperature

In-plane shear crack sub-critical propagation of rock at high temperature was studied by finite element method and shear-box(i.e.compression-shear) test with newly designed electrically conductive adhesive method.Numerical and experimental results show that the normalized shear(Mode Ⅱ) stress intensity factors,K ⅡT/KT0 is decreased as the temperature increases because high temperature can improve stress distribution at crack tip and reduce the Mode Ⅱ stress intensity factor.Microscopic features of fractured surface are of little pits and secondary micro-cracks in the vicinity(1.5-4.0 mm) of the crack tip.The chevron-shape secondary cracks gradually merge in the length of about 4-5 mm and disappear along the direction of crack propagation.Stable shear crack propagation time is increased with the increasing temperature while the stable shear crack propagation rate is decreased with the increasing temperature,since high temperature can increase the shear(Mode Ⅱ) fracture toughness and prevent the crack growth.It is necessary to ensure the ligament of specimen long enough to measure the maximum unstable crack propagation rate of rock.