Finite crack propagation in a micropolar elastic solid

The dynamic propagation of a finite crack under mode-1 loading in a micropolar elastic solid is investigated. By using an integral transform method, a pair of two-dimensional singular integral equations governing stress and couple stress is formulated in terms of displacement transverse to the crack, macro and micro rotations, and microinertia. These equations are solved numerically, and solutions for dynamic stress intensity and couple stress intensity factors are obtained by utilizing the values of the strengths of the square root singularities in macrorotation and the gradient of microrotation at the crack tips. The motion of the crack tips and the load on the crack surface are not prescribed in the formulation of the problem. Therefore, the method of solution is applicable to nonuniform rates of propagation of a crack under an arbitrary time-dependent load on the crack surface. As an example, the diffraction of a micropolar dilatational wave by a stationary crack is considered. The behavior of the microrotation field and the dynamic couple stress intensity factor, influenced by microinertia, in addition to the dynamic stress intensity factor, are examined. The classical elasticity solution for the corresponding problem arises as a special case when the micropolar moduli are dropped from the present solution.     

Analysis of a moving crack in a functionally graded strip between two homogeneous layers

In this paper a finite crack with constant length (Yoffe-type crack) propagating in a functionally graded strip with spatially varying elastic properties between two dissimilar homogeneous layers under in-plane loading was studied. By utilizing the Fourier transformation technique, the mixed boundary problem is reduced to a system of singular integral equations that are solved numerically. The influences of the geometric parameters, the graded parameter, the crack length and speed on the stress intensity factors are investigated. The numerical results show that the graded parameters, the thicknesses of the functionally graded strip and the two homogeneous layers, the crack size and speed have significant effects on the dynamic fracture behavior.     

研究冷挤压孔附近残余应力强度因子分布的激光散斑—权函数混合法

根据确定裂纹应力强度因子权函数法基本理论的特征，提出实验力学－权函数混合法的构想，为解决更复杂和三维问题等难提供一种有效的新途径。用激光散斑－权函数混合法研究了冷挤压孔附近残余应力强度因子的变化规律，得到的结果对冷挤压孔加工工艺中控制裂纹长度，提高构件的疲劳寿命和可靠性设计可起到指导性的作用。     

Anti-plane shear behavior of an arbitrarily oriented crack in bonded materials with a nonhomogeneous interfacial zone

The anti-plane shear problem of bonded elastic materials containing a crack at an arbitrary angle to the graded interfacial zone is investigated in this paper. The interfacial zone is modeled as a nonhomogeneous interlayer of finite thickness with the continuously varying shear modulus between the two dissimilar, homogeneous half-planes. Formulation of the crack problem is based upon the use of the Fourier integral transform method and the coordinate transformations of basic field variables. The resulting Cauchy-type singular integral equation is solved numerically to provide the values of modeIII stress intensity factors. A comprehensive parametric study is then presented of the influence of crack obliquity on the stress intensity factors for different crack size and locations and for different material combinations, in conjunction with the material nonhomogeneity within the graded interfacial zone.     

Analysis of a crack approaching a circular hole in cross-ply laminates under biaxial loading

The problem of a crack approaching a circular hole in cross-ply laminates under uniaxial and biaxial loading is investigated in this paper. The effects of material orthotropy, geometry [R/d and a/d], and loading conditions on crack tip singularity are investigated. The stress intensity factors are obtained by the modified mapping collocation method. The present results for an isotropic infinite plate show good agreement with existing solutions. The results for cross-ply laminates show that the stress intensity factors strongly depend on material orthotropy, geometry, and loading condition. The stress intensity factors for cross-ply laminates exist between those for θ=0° and those for θ=90° in the whole range of crack length and decrease as the percentage of 0° plies increases. In the range of small crack length the stress intensity factors for biaxial tension are higher than those for uniaxial tension. In the range of large crack length the stress intensity factors for uniaxial tension are higher than those for biaxial tension.     

Moving interfacial crack between two dissimilar soft ferromagnetic materials in uniform magnetic field

This paper presents the dynamic magnetoelastic stress intensity factors of a Yoffe-type moving crack at the interface between two dissimilar soft ferromagnetic elastic half-planes. The solids are subjected to a uniform in-plane magnetic field and the crack is opened by internal normal and shear tractions. The problem is considered within the framework of linear magnetoelasticity. By application of the Fourier integral transform, the mixed boundary problem is reduced to a pair of integral equations of the second kind with Cauchy-type singularities. The singular integral equations are solved by means of a Jacobi polynomial expansion method. For a particular case, closed-form solutions are obtained. It is shown that the magnetoelastic stress intensity factors depend on the moving velocity of the crack, the magnetic field and the magnetoelastic properties of the materials.     

Analysis of a flat annular crack under shear loading

An annular crack in an infinite isotropic elastic solid under shear loading is analyzed. General solution to the Navier's equilibrium equation is expressed in terms of three harmonic functions. Employing the Hankel transform the harmonic functions are represented by the solution of a pair of triple integral equations. The triple integral equations are reduced to a pair of mixed Volterra-Fredholm integral equations, which are numerically solved. The stress intensity factors of the annular crack under various shear loadings such as uniform radial shear, linearly varying radial shear, uniform shear and linearly varying shear are calculated as the Poisson's ratio ν anda/b (a; inner radius,b; outer radius) vary.     

φ30×108接链环的疲劳断裂分析

阐述利用有限元法对机械结构零部件进行疲劳断裂分析的相关理论和方法。并在此基础上对φ30×108接链环应用有限元分析软件ANSYS进行分析,确定裂纹易产生的危险截面,求解危险截面处不同深度和形态的表面椭圆裂纹的应力强度因子。     

风电齿轮齿根弯曲应力及裂纹特性分析

以风电齿轮为例,利用有限元分析软件ABAQUS,建立了太阳轮简化力学模型,并进行了齿根弯曲应力的分析计算。分析结果与传统计算方法得出的结果基本一致,从而验证了简化模型的正确性。在此基础上,研究了齿轮齿根裂纹特性,分析了初始裂纹长度和外加载荷对应力强度因子(SIF)的影响。结果表明,随着初始裂纹长度的增加,应力强度因子也随之增加,并且应力强度因子与载荷等比例增加。在初始裂纹长度和载荷相同的情况下,应力强度因子KⅠ远大于KⅡ和KⅢ,即在弯曲应力作用下张开型裂纹为风电齿轮轮齿折断失效的主要原因。     

Optimisation Method for Determination of Crack Tip Position Based on Gauss-Newton Iterative Technique

In the digital image correlation research of fatigue crack growth rate,the accuracy of the crack tip position determines the accuracy of the calculation of the stress intensity factor,thereby affecting the life prediction.This paper proposes a Gauss-Newton iteration method for solving the crack tip position.The conventional linear fitting method provides an iterative initial solution for this method,and the preconditioned conjugate gradient method is used to solve the ill-conditioned matrix.A noise-added artificial displacement field is used to verify the feasibility of the method,which shows that all parameters can be solved with satisfactory results.The actual stress intensity factor solution case shows that the stress intensity factor value obtained by the method in this paper is very close to the finite element result,and the relative error between the two is only-0.621％;The Williams coefficient obtained by this method can also better define the contour of the plastic zone at the crack tip,and the maximum relative error with the test plastic zone area is-11.29％.The relative error between the contour of the plastic zone defined by the conventional method and the area of the experimental plastic zone reached a maximum of 26.05％.The crack tip coordinates,stress intensity factors,and plastic zone contour changes in the loading and unloading phases are explored.The results show that the crack tip change during the loading process is faster than the change during the unloading process;the stress intensity factor during the unloading process under the same load condition is larger than that during the loading process;under the same load,the theoretical plastic zone during the unloading process is higher than that during the loading process.     

Numerical studies on crack growth in a steel tubular T-joint

The results of a numerical investigation on fatigue crack growth of an offshore tubular T-joint under the action of axial, in-plane and out-of-plane bending loads are presented in this paper. Extensive stress analysis has been carried out to determine the location of the hot spot stress along the brace—chord intersection for each load case. Semi-elliptical cracks with varying crack lengths and crack depths were introduced at the hot spot location by means of line spring elements for stress intensity factor evaluation. The line springs were properly constrained to prevent the problem of crack surface penetration. The stress intensity factors obtained are then used in a crack growth law for life estimation.     

The singular stress field and stress intensity factors of a crack terminating at a bimaterial interface

For the fracture evaluation of inclined cracks terminating at the dissimilar material interface, not only the singularities, but also the detailed stress field and its stress intensity factors are necessary. However, though there are many researches reported on the singularity analysis, the stress field and its stress intensity factors are still not clear. This paper has deduced theoretically the singular stress and displacement fields near the tip of a crack terminating at the interface between bonded dissimilar materials, for both cases of real and oscillatory singularities. From the deduced singular stress field, the stress intensity factors are defined for such a crack, and the corresponding numerical extrapolation methods are also proposed. Through the numerical examinations, it is found that the theoretical stress distributions agree well with the numerical results obtained by the finite element method. Moreover, the proposed extrapolation method shows a good linearity, thus it can be used as an efficient way to determine the characteristics of the stress and displacement fields near the tip of a crack terminating at interface.     

Analysis of an unsteadily propagating crack under mode I and II loading

Stress and displacement fields for an unsteadily propagating crack under mode I and II loading are developed through an asymptotic analysis. Dynamic equilibrium equations for the unsteady state are developed and the solution to the displacement fields and the stress fields for a crack propagating with high crack tip acceleration, deceleration and rapidly varying stress intensity factor. The influence of transients on the higher order terms of the stress fields are explicitly revealed. Using these stress components, isochromatic fringes around the propagating crack are generated for different crack speeds, crack tip accelerations and the time rate of change of stress intensity factor, and the effects of the transients on these fringes are discussed. The effects of the transients on the dynamic stress intensity factor are discussed when a crack propagates with high acceleration and deceleration. The effect of transient on the time rate of change of dynamic stress intensity factor below Rayleigh wave speed in an infinite body is also studied.     

Discussion on calculation of the local flexibility due to the crack in a pipe

The stress intensity factor plays an important role in the calculation of the local flexibility due to the crack in a cracked structure. Many researches on the stress intensity factors in a cracked beam or rotor have been made, but there are some difficulties in calculating the stress intensity factors in a cracked pipe. Maiti et al. obtained the local flexibility due to the crack in a pipe experimentally by deflection and natural frequency methods without calculating the stress intensity factor. In this paper, the stress intensity factor in a cracked pipe can be calculated by dividing a cracked pipe into a series of thin annuli, and a method to calculate the local flexibility due to the crack in a pipe is presented. The experiment results are given to verify the effectiveness of such a method.     

Transient analysis of several cracks in an orthotropic material under impact load

In this paper, the dynamic response of an orthotropic plane with multiple cracks is analyzed analytically. The anti-plane impact loading condition is considered. The dislocation solution is utilized to derive integral equations for multiple interacting cracks with any location and orientation in an orthotropic plane. These equations are solved numerically thereby obtaining the dislocation density function on the crack faces and stress intensity factors of cracks. Numerical results are presented to illustrate the influences of parameters such as cracks geometry and orientation on the dynamic stress intensity factors.     

Enriched finite element method for three-dimensional viscoelastic interface crack problems

The enriched finite element method is developed for three-dimensional problems of an interface crack between elastic and viscoelastic (including dissimilar viscoelastic) materials. According to the displacement fields of elastic interface crack, the displacement fields of viscoelastic interface crack are derived through the correspondence principle. By incorporating the displacement expressions into the displacement model of regular element, the incremental formulations of enriched element are derived. The stress intensity factors and strain energy release rates can be solved based on the enriched degree of freedoms. A 3-D through interface crack at the center of jointed dissimilar viscoelastic plate subjected to remote tension and a quarter-circular viscoelastic interface corner crack subjected to uniform thermal loading was investigated using the enriched finite element method. It is shown that the present solutions are consistent with the analytical solutions, which indicates the present method is correct and efficient.

     

Analysis of a conducting crack in an electrostrictive ceramic under combined electric and mechanical loading

A conducting crack in an electrostrictive ceramic under combined electric and mechanical loading is investigated. Analysis based on linear dielectric model predicts that the surfaces of the crack are not open completely but they are contact near the crack tip. The complete solution for the crack with a contact zone in a linear electrostrictive ceramic under combined electric and mechanical loading is obtained by using the complex variable formula. The asymptotic problems for a semi-infinite crack with a partial opening zone as well as for a fully open semi-infinite crack in a nonlinear electrostrictive ceramic are analyzed in order to investigate the effect of the electrical nonlinearity on the stress intensity factor under small scale nonlinear conditions. Particular attention is devoted to a finite crack in the nonlinear electrostrictive ceramic subjected to combined electric and mechanical loading. The stress intensity factor for the finite crack under small scale nonlinear conditions is obtained from the asymptotic analysis.     

Asymptotic analysis of a wedge crack in an anisotropic dielectric material

The asymptotic problem of a conducting crack emanating from the vertex of a wedge in an anisotropic dielectric material under purely electric loading is investigated. The wedge crack in an isotropic material is solved using a conformal mapping technique. The solution of the wedge crack in an anisotropic dielectric material is obtained from that for the transformed isotropic problem after applying a linear transformation method. The electric field intensity factor for the anisotropic wedge crack is obtained in the closed form. The effects of crack and wedge angles as well as anisotropic parameters on the electric field intensity factor are illustrated. The electric field intensity factor is numerically calculated by using the J integral and finite element analysis to validate the exact solution of the electric field intensity factor.     

Stress intensity factors of a crack lying in a viscoelastic layer embedded in two different media under an anti-plane concentrated load

A general solution to the problem of a crack lying in a viscoelastic layer sandwiched between two different media under an anti-plane concentrated load is presented. The fundamental elastic solution is obtained as a rapidly convergent series in terms of complex potentials via successive iterations of Möbius’ transformation in order to satisfy continuity conditions on two different interfaces. Based on the correspondence principle, the Laplace-transformed viscoelastic solution is then directly determined from the corresponding elastic solution. Therefore, the time-dependent mode-III stress intensity factor (SIF) of the crack tip can be solved numerically in a straightforward manner. Furthermore, the viscoelastic moduli of polyurethane at various temperatures are regressed from a uni-axial test with a constant load while a standard solid model is applied to formulate the viscoelastic constitutive equation. Finally, the effects of working temperature, boundary conditions and the geometric relation on the SIF are also discussed.     

Characteristics of a transiently propagating crack in functionally graded materials

When a crack propagates with acceleration, deceleration and time rates of change of stress intensity factors, it is very important for us to understand the effects of acceleration, deceleration and time rates of change of stress intensity factors on the individual stresses and displacements at the crack tip. Therefore, the crack tip stress and displacement fields for a transiently propagating crack along gradient in functionally graded materials (FGMs) with an exponential variation of shear modulus and density are developed and the characteristics of a transiently propagating crack from the fields are analyzed. The effects of the rate of change of the stress intensity factor and the crack tip acceleration on the individual stresses at the crack tip are opposite each other. Specially, the isochromatics (constant maximum shear stress) of Mode I tilt backward around the crack tip with an increase of crack tip acceleration, and tilt forward around the crack tip with an increase of the rate of change of the dynamic mode I stress intensity factor. This paper was recommended for publication in revised form by Associate Editor Chongdu Cho Kwang-Ho Lee received a Ph.D. degree in Yeungnam University in 1993. Dr. Lee is currently a professor at the School of Mechanical and Automotive Engineering at Kyungpook National University in Korea. He also had worked in KOMSCO as an engineer and researcher (1982.3–1996.2). He is interested in the fields of fracture and stress analysis on the composite, interface, nano and functionally graded materials by theoretical and experimental mechanics. Specially, his major interest is analysis of dynamic crack tip fields. Young-Jae Lee received his B.S degree in Agricultural Civil Engineering from Gyeongsang National University (GNU) in 1982. He then received his M.S. and Ph.D. degrees from GNU in 1984 and 1995, respectively. Dr. Lee is currently a professor at the department of Civil Engineering at Kyungpook National University in Korea. From 2005 to 2006, he had served as an editor of Korea Institute for Structure Maintenance and Inspection. His research interests are in the area of evaluation, diagnosis and optimum design of structure. Sang-Bong Cho received a Ph. D. degree from Tokyo University in 1989. Dr. Cho is currently a professor at the division of Mechanical and Automation Engineering at Kyungnam University in Korea. His research interests are in the area of fracture mechanics, FEM stress analysis and fretting fatigue.