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.

Dynamic Phase Field Description on Ductile Fracture Process

In comparison to discrete descriptions of fracture process,the recently proposed phase field methodology averts the numerical tracking strategy of discontinuities in solids,which enables the numerical implement simplification.An implicit finite element formulation based on the diffuse phase field is extended for stable and efficient analysis of complex dynamic fracture process in ductile solids.This exhibited formulation is shown to capture entire range of the characteristics of ductile material presenting J2-plasticity,embracing plasticization,cracks initiation,propagation,branching and merging while fulfilling the basic principle of thermodynamics.Herein,we implement a staggered time integration scheme of the dynamic elasto-plastic phase field method into the commercial finite element code.The numerical performance of the present advanced phase field model has been examined through several classic dynamic fracture benchmarks,and in all cases simulation results are in good agreement with the associated experimental data and other numerical results in previous literature.

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 MIXED MODE FRACTURE CRITERION BASED ON THE MAXIMUM TANGENTIAL STRESS IN BRITTLE INCLUSION

A closed-form solution for predicting the tangential stress of an inclusion located in mixed mode Ⅰ and Ⅱ crack tip field was developed based on the Eshelby equivalent inclusion theory. Then a mixed mode fracture criterion, including the fracture direction and the critical load, was established based on the maximum tangential stress in the inclusion for brittle inclusioninduced fracture materials. The proposed fracture criterion is a function of the inclusion fracture stress, its size and volume fraction, as well as the elastic constants of the inclusion and the matrix material. The present criterion will reduce to the conventional one as the inclusion having the same elastic behavior as the matrix material. The proposed solutions are in good agreement with detailed finite element analysis and measurement.

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.

Evaluation of parameters affecting reflection cracking in geogrid-reinforced overlay

The influence of the most important parameters on the service life of reinforced asphalt overlay with geogrid materials in bending mode was examined by employing the Taguchi method and analysis of variance techniques. The objectives of this experiment was to investigate the effects of grid stiffness, tensile strength, coating type, amount of tack coat, overlay thickness, crack width and stiffnesses of asphalt overlay and existing asphalt concrete on propagation of the reflection cracking. Results indicate that the stiffnesses of cracked layer and overlay are the main significant factors that can directly improve the service life of an overlay against the reflection cracking. Generally, glass grid is more effective in reinforced overlay than polyester grid. Effect of crack width of the existing layer is significant when its magnitude increases from 6 to 9 mm.

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.

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.

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.

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.

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.

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.

THE ANALYSES OF THE THREE-DIMENSIONAL STRESS STRUCTURE NEAR THE CRACK TIP OF MODE I CT SPECIMENS IN ELASTICPLASTIC STATE(Ⅰ)--THE ANALYSES OF CONSTRAINT PARAMETERS AND FRACTURE PARAMETERS

In the present paper,three dimensional analyses of some general constraint parameters and fracture parameters near the crack tip of Mode I CT specimens in two different thicknesses are carried out by employing ADINA program.The results reveal that the constraints along the thickness direction are obviously separated into two parts:the keeping similar high constraint field(Z_(1))and rapid reducing constraints one(Z_(2)).The two fields are experimentally confiremed to correspond to the smooth region and the shear lip on the fracture face respectively.So the three dimensional stress structure of Mode I specimens can be derived through discussing the two fields respectively.The distribution of the Crack Tip Opening Displacement(CTOD)along the thickness direction and the three dimensional distribution of the void growth ratio(V_(g))near the crack tip are also obtained.The two fracture parameters are in similar trends along the thickness direction,and both of them can reflect the effect of thickness and that of the loading level to a certain degree.

The crack tip fields of Mode Ⅱ stationary growth crack on bimaterial interface

A mechanical model is established for mode II interfacial crack static growing along an elastic-elastic power law creeping bimaterial interface. For frictional contact of boundary conditions on crack faces, asymptotic solutions of the stresses and strains of near tip-crack are got. It was shown that in stable creep growing phase, elastic deformation and viscous deformation are equally dominant in the near-tip field, the stress and strain have the same singularity and there is not the oscillatory singularity the field. Through numerical calculation , it is shown that the frictional coefficient η notably influence the crack-tip field.