Stress waves affect the stress field at the crack tip and dominate the dynamic crack propagation.Therefore,evaluating the influence of blasting stress waves on the crack propagation behavior and the mechanical charact...Stress waves affect the stress field at the crack tip and dominate the dynamic crack propagation.Therefore,evaluating the influence of blasting stress waves on the crack propagation behavior and the mechanical characteristics of crack propagation is of great significance for engineering blasting.In this study,ANSYS/LS-DYNA was used for blasting numerical simulation,in which the propagation characteristics of blasting stress waves and stress field distribution at the crack tip were closely observed.Moreover,ABAQUS was applied for simulating the crack propagation path and calculating dynamic stress intensity factors(DSIFs).The universal function was calculated by the fractalmethod.The results show that:the compressive wave causes the crack to close and the reflected tensile wave drives the crack to initiate and propagate,and failure mode is mainly tensile failure.The crack propagation velocity varies with time,which increases at first and then decreases,and the crack arrest occurs due to the attenuation of stress waves and dissipation of the blasting energy.In addition,crack arrest toughness is smaller than the crack initiation toughness,applied pressure waveforms(such as the peak pressure,duration,waveforms,wavelengths and loading rates)have a great influence on DSIFs.It is conducive to our deep understanding or the study of blasting stress waves dominated fracture,suggesting a broad reference for the further development of rock blasting in engineering practice.展开更多
This paper studies the dynamic stress intensity factor (DSIF) at the interface in an adhesive joint under shear loading. Material damage is considered. By introducing the dislocation density function and using the i...This paper studies the dynamic stress intensity factor (DSIF) at the interface in an adhesive joint under shear loading. Material damage is considered. By introducing the dislocation density function and using the integral transform, the problem is reduced to algebraic equations and can be solved with the collocation dots method in the Laplace domain. Time response of DSIF is calculated with the inverse Laplace integral transform. The results show that the mode Ⅱ DSIF increases with the shear relaxation parameter, shear module and Poisson ratio, while decreases with the swell relaxation parameter. Damage shielding only occurs at the initial stage of crack propagation. The singular index of crack tip is -0.5 and independent on the material parameters, damage conditions of materials, and time. The oscillatory index is controlled by viscoelastic material parameters.展开更多
A dynamic weight function method is presented for dynamic stress intensity factors of circular disk with a radial edge crack under external impulsive pressure. The dynamic stresses in a circular disk are solved under ...A dynamic weight function method is presented for dynamic stress intensity factors of circular disk with a radial edge crack under external impulsive pressure. The dynamic stresses in a circular disk are solved under abrupt step external pressure using the eigenfunction method. The solution consists of a quasi-static solution satisfying inhomogeneous boundary conditions and a dynamic solution satisfying homogeneous boundary conditions. By making use of Fourier- Bessel series expansion, the history and distribution of dynamic stresses in the circular disk are derived. Furthermore, the equation for stress intensity factors under uniform pressure is used as the reference case, the weight function equation for the circular disk containing an edge crack is worked out, and the dynamic stress intensity factor equation for the circular disk containing a radial edge crack can be given. The results indicate that the stress intensity factors under sudden step external pressure vary periodically with time, and the ratio of the maximum value of dynamic stress intensity factors to the corresponding static value is about 2.0.展开更多
The time-domain BEM was developed to analyze the dynamic stress intensity factor ( DSIF) of 3-D elastodynamic crack problems. To simulate the stress singularity along the front of a crack, eight-node isoparametric sin...The time-domain BEM was developed to analyze the dynamic stress intensity factor ( DSIF) of 3-D elastodynamic crack problems. To simulate the stress singularity along the front of a crack, eight-node isoparametric singular elements were used, and the DSIF for a semi-circular surface crack was firstly calculated based on displacement equation using the time-domain BEM formulation. The new scheme to determine the time step was brought forward. By the dynamic analysis program of time-domain BEM compiled by its, several numerical examples are presented, which demonstrate the unconditional stability and high accuracy of time-domain BEM applied to 3-D elastodynamic crack problems.展开更多
Dynamic stress intensity factors are evaluated for thick-walled cylinder with a radial edge crack under internal impulsive pressure. Firstly, the equation for stress intensity factors under static uniform pressure is ...Dynamic stress intensity factors are evaluated for thick-walled cylinder with a radial edge crack under internal impulsive pressure. Firstly, the equation for stress intensity factors under static uniform pressure is used as the reference case, and then the weight function for a thick-walled cylinder containing a radial edge crack can be worked out. Secondly, the dynamic stresses in uncracked thick-walled cylinders are solved under internal impulsive pressure by using mode shape function method. The solution consists of a quasi-static solution satisfying inhomogeneous boundary conditions and a dynamic solution satisfying homogeneous boundary condi- tions, and the history and distribution of dynamic stresses in thick-walled cylinders are derived in terms of Fourier-Bessel series. Finally, the dynamic stress intensity factor equations for thick-walled cylinder containing a radial edge crack sub- jected to internal impulsive pressure are given by dynamic weight function method. The finite element method is utilized to verify the results of numerical examples, showing the validity and feasibility of the proposed method.展开更多
A new formula is obtained to calculate dynamic stress intensity factors of the three-point bending specimen containing a single edge crack in this study. Firstly, the weight function for three-point bending specimen c...A new formula is obtained to calculate dynamic stress intensity factors of the three-point bending specimen containing a single edge crack in this study. Firstly, the weight function for three-point bending specimen containing a single edge crack is derived from a general weight function form and two reference stress intensity factors, the coefficients of the weight function are given. Secondly, the history and distribution of dynamic stresses in uncracked three-point bending specimen are derived based on the vibration theory. Finally~ the dynamic stress intensity factors equations for three-pointing specimen with a single edge crack subjected to impact loadings are obtained by the weight function method. The obtained formula is verified by the comparison with the numerical results of the finite element method (FEM). Good agreements have been achieved. The law of dynamic stress intensity factors of the three-point bending specimen under impact loadings varing with crack depths and loading rates is studied.展开更多
This paper presents a formulation for three-dimensional elasto-dynamics with an elliptic crack based on the Laplace and Fourier transforms and the convolution theorem. The dynamic stress intensity factor for the crack...This paper presents a formulation for three-dimensional elasto-dynamics with an elliptic crack based on the Laplace and Fourier transforms and the convolution theorem. The dynamic stress intensity factor for the crack is determined by solving a Fredholm integral equation of the first kind. The results of this paper are very close to those given by the two-dimensional dual integral equation method.展开更多
The problem of a Griffith crack in an unbounded orthotropic functionally graded material subjected to antipole shear impact was studied. The shear moduli in two directions of the functionally graded material were assu...The problem of a Griffith crack in an unbounded orthotropic functionally graded material subjected to antipole shear impact was studied. The shear moduli in two directions of the functionally graded material were assumed to vary proportionately as definite gradient. By using integral transforms and dual integral equations, the local dynamic stress field was obtained. The results of dynamic stress intensity factor show that increasing shear moduli's gradient of FGM or increasing the shear modulus in direction perpendicular to crack surface can restrain the magnitude of dynamic stress intensity factor.展开更多
The method of wave-function expansion in elliptical coordinates,elliptical cosine half-range expansion and Mathieu function were applied to obtain an exact analytical solution of the dynamic stress concentration facto...The method of wave-function expansion in elliptical coordinates,elliptical cosine half-range expansion and Mathieu function were applied to obtain an exact analytical solution of the dynamic stress concentration factor(DSCF)around an elliptical cavity in a shallow,semi-elliptical hill.An infinite system of simultaneous linear equations for solving this problem was established by substituting the wave expression obtained by the Mathieu function including the standing wave expression of elliptical lining given herein into the boundary condition obtained by the region-matching method.The finite equations system with unknown coefficients obtained by truncation were solved numerically,and the results in the case of an ellipse degenerating into a circle were compared with previous results to verify the accuracy of the method.The effects of different aspect ratios,incident wave angles and aperture ratios on the dynamic stress concentration around the elliptical cavity were described.Some numerical results,when the elliptical hill was changed into a circular one,were analyzed and compared in detail.In engineering,this model can be regarded as a semi-cylindrical hill with an elliptical cylindrical unlined tunnel under the action of SH waves,and the results are significant in aseismic design.展开更多
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 combined mode point loads ...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 combined mode point loads that move in a direction perpendicular to the crack edge is considered.The analytic expression for the combined mode stress intensity factors as a function of time for any point along the crack edge is obtained.The method of solution is based on the application of integral transform together with the Wiener-Hopf technique and the Cagniard-de Hoop method. Some features of the solution are discussed and graphical results for various point load speeds are presented.展开更多
Using digital laser dynamic caustics experimental system and conducting simulation experiment researched the influence rule of blasting excavation of a new roadway on neighboring existed different cross-section roadwa...Using digital laser dynamic caustics experimental system and conducting simulation experiment researched the influence rule of blasting excavation of a new roadway on neighboring existed different cross-section roadways. The experimental results show that the influence of blast load on adjacent roadway has a good relationship with the cross-section of roadway. The expansion distance of precrack existed in circular, arch-wall, rectangular roadway is respectively 1.76, 1.61 and 0 cm under blast load.At the same time, the direct-blast side of rectangular roadway has more obvious damage compared with circular and arch-wall roadway. It explains that plane reflects more stress wave than arc, so that it exerts more tensile failure in the direct-blast side, which leads to less stress wave diffracting to the precrack in the back-blast side. When the precrack extends, higher value dynamic stress intensity factor in circular roadway works longer than that of arch-wall roadway. Indirectly, it explains that plane's weakening function on stress wave is significantly stronger than arc. Stress wave brings about self-evident influence on the upper and bottom endpoints of the rectangular roadway, and it respectively extends 1.03, 2.06 cm along the line link direction of the center of the blasthole and the upper and bottom endpoints on the right wall.展开更多
The present paper is exposed theoretically to the influence on the dynamic stress intensity factor (DSIF) in the piezoelectric bi-materials model with two symmet- rically permeable interracial cracks near the edges ...The present paper is exposed theoretically to the influence on the dynamic stress intensity factor (DSIF) in the piezoelectric bi-materials model with two symmet- rically permeable interracial cracks near the edges of a circular cavity, subjected to the dynamic incident anti-plane shearing wave (SH-wave). An available theoretical method to dynamic analysis in the related research field is provided. The formulations are based on Green's function method. The DSIFs at the inner and outer tips of the left crack are obtained by solving the boundary value problems with the conjunction and crack- simulation technique. The numerical results are obtained by the FORTRAN language program and plotted to show the influence of the variations of the physical parameters, the structural geometry, and the wave frequencies of incident wave on the dimensionless DSIFs. Comparisons with previous work and between the inner and outer tips are con- cluded.展开更多
The scaled boundary finite element method (SBFEM) is a recently developed numerical method combining advantages of both finite element methods (FEM) and boundary element methods (BEM) and with its own special fe...The scaled boundary finite element method (SBFEM) is a recently developed numerical method combining advantages of both finite element methods (FEM) and boundary element methods (BEM) and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors (SIFs) directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.展开更多
In this paper,dynamic caustic method is applied to analyze the blast-induced crack propagation and distribution of the dynamic stress field around an empty circular hole in polymethyl methacrylate(PMMA)material under ...In this paper,dynamic caustic method is applied to analyze the blast-induced crack propagation and distribution of the dynamic stress field around an empty circular hole in polymethyl methacrylate(PMMA)material under delayed ignition blasting loads.The following experimental results are obtained.(1)In directional-fracture-controlled blasting,the dynamic stress intensity factors(DSIFs)and the propagation paths of the blast-induced cracks are obviously influenced by the delayed ignition.(2) The circular hole situated between the two boreholes poses a strong guiding effect on the coelesence of the cracks,causing them to propagate towards each other when cracks are reaching the circular hole area.(3)Blast-induced cracks are not initiated preferentially because of the superimposed effect from the explosive stress waves on the cracking area.(4) By using the scanning electron microscopy(SEM)method,it is verified that the roughness of crack surfaces changes along the crack propagation paths.展开更多
The scattering of shear waves (SH waves) by nano-scale arbitrary shape inclusion in infinite plane is studied by complex variable function theory. Firstly, the governing equation and the relationships between stress a...The scattering of shear waves (SH waves) by nano-scale arbitrary shape inclusion in infinite plane is studied by complex variable function theory. Firstly, the governing equation and the relationships between stress and displacement are given by classical elastic theory. Secondly, the arbitrary shape inclusion in the two-dimensional plane is transformed into a unit circle domain by conformal mapping, the incident wave field and the scattered wave field are presented. Next, the stress and displacement boundary conditions are established by considering surface elasticity theory, The infinite algebraic equations for solving the unknown coefficients of the scattered and standing waves are obtained. Finally, the influence of surface effect, non-dimensional wave number, Shear modulus and hole curvature on the dynamic stress concentration factor are analyzed by some examples, the numerical results show that the surface effect weakens the dynamic stress concentration. With the increase of wave number, the dynamic stress concentration factor (DSCF) decreases. Shear modulus and hole curvature have significant effects on DSCF.展开更多
Three dimensional dynamic stress intensity factors are analyzed for a curved crack with a second order perturbation method. The method is extended to obtain an approximate representation of a three dimensional dynamic...Three dimensional dynamic stress intensity factors are analyzed for a curved crack with a second order perturbation method. The method is extended to obtain an approximate representation of a three dimensional dynamic stress intensity factors at the tip of a curved crack. Due to three dimensional curved crack growth the dynamic energy release rate can be calculated by using the Irwin's formula. A three dimensional curved crack in materials with inhomogeneous fracture toughness are considered. Paths of a brittle three dimensional curved crack propagating along a welded joint are predicted via the present method, where the effects of dynamic applied stresses, residual stresses, and material deterioration due to welding are taken into considerations.展开更多
The transient fracture behavior of a functionally graded layered structure subjected to an in-plane impact load is investigated. The studied structure is composed of two homogeneous layers and a functionally gradedint...The transient fracture behavior of a functionally graded layered structure subjected to an in-plane impact load is investigated. The studied structure is composed of two homogeneous layers and a functionally gradedinterlayer with a crack perpendicular to the boundaries. The impact load is applied on the face of the crack. Fourier transform and Laplace transform methods are used to formulate the present problem in terms of a singular integral equation in Laplace transform domain. Considering variations of parameters such as the nonhomogeneity constant, the thickness ratio and the crack length, the dynamic stress intensity factors (DSIFs) in time domain are studied and some meaningful conclusions are obtained.展开更多
The scattering of SH wave by a cylindrical piezoelectric inclusion partially debonded from its surrounding piezoelectric material is investigated using the wave function expansion method and singular integral ...The scattering of SH wave by a cylindrical piezoelectric inclusion partially debonded from its surrounding piezoelectric material is investigated using the wave function expansion method and singular integral equation technique. The debonding regions are modeled as mul- tiple arc-shaped interface cracks with non-contacting faces. By expressing the scattered ?elds as wave function expansions with unknown coe?cients, the mixed boundary value problem is ?rstly reduced to a set of simultaneous dual series equations. Then dislocation density functions are introduced as unknowns to transform these dual series equations into Cauchy singular integral equations of the ?rst type, which can be numerically solved easily. The solution is valid for arbi- trary number and size of the debonds. Finally, numerical results of the dynamic stress intensity factors are presented for the cases of one debond and two debonds. The e?ects of incidence direc- tion, crack con?guration and various material parameters on the dynamic stress intensity factors are respectively discussed. The solution of this problem is expected to ?nd applications in the investigation of dynamic fracture properties of piezoelectric materials with cracks.展开更多
Due to the incompatibility of the interlaminar deformations,the interface debonding or cracking usually happens in a layered magnetoelectric(ME)structure under an applied load.In this paper,the transient responses of ...Due to the incompatibility of the interlaminar deformations,the interface debonding or cracking usually happens in a layered magnetoelectric(ME)structure under an applied load.In this paper,the transient responses of the anti-plane interface cracks in piezoelectric(PE)-piezomagnetic(PM)sandwich structures are studied by the standard methods of the integral transform and singular integral equation.Discussion on the numerical examples indicates that the PE-PM-PE structure under electric impact is more likely to fracture than the PM-PE-PM structure under a magnetic impact.The dynamic stress intensity factors(DSIFs)are more sensitive to the variation of the active layer thickness.The effects of the material constants on the DSIFs are dependent on the roles played by PE and PM media during the deformation process.展开更多
A mathematical formulation is presented for the dynamic stress intensity factor (mode I) of a finite permeable crack subjected to a time-harmonic propagating longitudinal wave in an infinite poroelastic solid. In part...A mathematical formulation is presented for the dynamic stress intensity factor (mode I) of a finite permeable crack subjected to a time-harmonic propagating longitudinal wave in an infinite poroelastic solid. In particular, the effect of the wave-induced fluid flow due to the presence of a liquid-saturated crack on the dynamic stress intensity factor is analyzed. Fourier sine and cosine integral transforms in conjunction with Helmholtz potential theory are used to formulate the mixed boundary-value problem as dual integral equations in the frequency domain. The dual integral equations are reduced to a Fredholm integral equation of the second kind. It is found that the stress intensity factor monotonically decreases with increasing frequency, decreasing the fastest when the crack width and the slow wave wavelength are of the same order. The characteristic frequency at which the stress intensity factor decays the fastest shifts to higher frequency values when the crack width decreases.展开更多
基金This researchwas supported by the National Natural Science Foundation of China(No.52227805)the Fundamental Research Funds for Central Universities(No.2022JCCXLJ01).Awards were granted to the author Liyun Yang.
文摘Stress waves affect the stress field at the crack tip and dominate the dynamic crack propagation.Therefore,evaluating the influence of blasting stress waves on the crack propagation behavior and the mechanical characteristics of crack propagation is of great significance for engineering blasting.In this study,ANSYS/LS-DYNA was used for blasting numerical simulation,in which the propagation characteristics of blasting stress waves and stress field distribution at the crack tip were closely observed.Moreover,ABAQUS was applied for simulating the crack propagation path and calculating dynamic stress intensity factors(DSIFs).The universal function was calculated by the fractalmethod.The results show that:the compressive wave causes the crack to close and the reflected tensile wave drives the crack to initiate and propagate,and failure mode is mainly tensile failure.The crack propagation velocity varies with time,which increases at first and then decreases,and the crack arrest occurs due to the attenuation of stress waves and dissipation of the blasting energy.In addition,crack arrest toughness is smaller than the crack initiation toughness,applied pressure waveforms(such as the peak pressure,duration,waveforms,wavelengths and loading rates)have a great influence on DSIFs.It is conducive to our deep understanding or the study of blasting stress waves dominated fracture,suggesting a broad reference for the further development of rock blasting in engineering practice.
基金the National Natural Science Foundation of China(No.10672027)the National Basic Research Program of China(No.2006CB601205)the National Science Fund for Distin-guished Young Scholars of China(No.50625414)
文摘This paper studies the dynamic stress intensity factor (DSIF) at the interface in an adhesive joint under shear loading. Material damage is considered. By introducing the dislocation density function and using the integral transform, the problem is reduced to algebraic equations and can be solved with the collocation dots method in the Laplace domain. Time response of DSIF is calculated with the inverse Laplace integral transform. The results show that the mode Ⅱ DSIF increases with the shear relaxation parameter, shear module and Poisson ratio, while decreases with the swell relaxation parameter. Damage shielding only occurs at the initial stage of crack propagation. The singular index of crack tip is -0.5 and independent on the material parameters, damage conditions of materials, and time. The oscillatory index is controlled by viscoelastic material parameters.
文摘A dynamic weight function method is presented for dynamic stress intensity factors of circular disk with a radial edge crack under external impulsive pressure. The dynamic stresses in a circular disk are solved under abrupt step external pressure using the eigenfunction method. The solution consists of a quasi-static solution satisfying inhomogeneous boundary conditions and a dynamic solution satisfying homogeneous boundary conditions. By making use of Fourier- Bessel series expansion, the history and distribution of dynamic stresses in the circular disk are derived. Furthermore, the equation for stress intensity factors under uniform pressure is used as the reference case, the weight function equation for the circular disk containing an edge crack is worked out, and the dynamic stress intensity factor equation for the circular disk containing a radial edge crack can be given. The results indicate that the stress intensity factors under sudden step external pressure vary periodically with time, and the ratio of the maximum value of dynamic stress intensity factors to the corresponding static value is about 2.0.
文摘The time-domain BEM was developed to analyze the dynamic stress intensity factor ( DSIF) of 3-D elastodynamic crack problems. To simulate the stress singularity along the front of a crack, eight-node isoparametric singular elements were used, and the DSIF for a semi-circular surface crack was firstly calculated based on displacement equation using the time-domain BEM formulation. The new scheme to determine the time step was brought forward. By the dynamic analysis program of time-domain BEM compiled by its, several numerical examples are presented, which demonstrate the unconditional stability and high accuracy of time-domain BEM applied to 3-D elastodynamic crack problems.
基金supported by the China Aviation Industry Corporation I Program (ATPD-1104-02).
文摘Dynamic stress intensity factors are evaluated for thick-walled cylinder with a radial edge crack under internal impulsive pressure. Firstly, the equation for stress intensity factors under static uniform pressure is used as the reference case, and then the weight function for a thick-walled cylinder containing a radial edge crack can be worked out. Secondly, the dynamic stresses in uncracked thick-walled cylinders are solved under internal impulsive pressure by using mode shape function method. The solution consists of a quasi-static solution satisfying inhomogeneous boundary conditions and a dynamic solution satisfying homogeneous boundary condi- tions, and the history and distribution of dynamic stresses in thick-walled cylinders are derived in terms of Fourier-Bessel series. Finally, the dynamic stress intensity factor equations for thick-walled cylinder containing a radial edge crack sub- jected to internal impulsive pressure are given by dynamic weight function method. The finite element method is utilized to verify the results of numerical examples, showing the validity and feasibility of the proposed method.
基金supported by the China Aviation Industry Corporation I Program (No.ATPD-1104-02)the Science Foundation of Nanjing University of Science and Technology (No.2010GJPY026)
文摘A new formula is obtained to calculate dynamic stress intensity factors of the three-point bending specimen containing a single edge crack in this study. Firstly, the weight function for three-point bending specimen containing a single edge crack is derived from a general weight function form and two reference stress intensity factors, the coefficients of the weight function are given. Secondly, the history and distribution of dynamic stresses in uncracked three-point bending specimen are derived based on the vibration theory. Finally~ the dynamic stress intensity factors equations for three-pointing specimen with a single edge crack subjected to impact loadings are obtained by the weight function method. The obtained formula is verified by the comparison with the numerical results of the finite element method (FEM). Good agreements have been achieved. The law of dynamic stress intensity factors of the three-point bending specimen under impact loadings varing with crack depths and loading rates is studied.
基金The project supported by the National Natural Science Foundation of China (K19672007)
文摘This paper presents a formulation for three-dimensional elasto-dynamics with an elliptic crack based on the Laplace and Fourier transforms and the convolution theorem. The dynamic stress intensity factor for the crack is determined by solving a Fredholm integral equation of the first kind. The results of this paper are very close to those given by the two-dimensional dual integral equation method.
文摘The problem of a Griffith crack in an unbounded orthotropic functionally graded material subjected to antipole shear impact was studied. The shear moduli in two directions of the functionally graded material were assumed to vary proportionately as definite gradient. By using integral transforms and dual integral equations, the local dynamic stress field was obtained. The results of dynamic stress intensity factor show that increasing shear moduli's gradient of FGM or increasing the shear modulus in direction perpendicular to crack surface can restrain the magnitude of dynamic stress intensity factor.
基金Fundamental Research Funds for Central Universities under Grant No.3072019CF0205。
文摘The method of wave-function expansion in elliptical coordinates,elliptical cosine half-range expansion and Mathieu function were applied to obtain an exact analytical solution of the dynamic stress concentration factor(DSCF)around an elliptical cavity in a shallow,semi-elliptical hill.An infinite system of simultaneous linear equations for solving this problem was established by substituting the wave expression obtained by the Mathieu function including the standing wave expression of elliptical lining given herein into the boundary condition obtained by the region-matching method.The finite equations system with unknown coefficients obtained by truncation were solved numerically,and the results in the case of an ellipse degenerating into a circle were compared with previous results to verify the accuracy of the method.The effects of different aspect ratios,incident wave angles and aperture ratios on the dynamic stress concentration around the elliptical cavity were described.Some numerical results,when the elliptical hill was changed into a circular one,were analyzed and compared in detail.In engineering,this model can be regarded as a semi-cylindrical hill with an elliptical cylindrical unlined tunnel under the action of SH waves,and the results are significant in aseismic design.
基金the National Natural Science Foundation of China
文摘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 combined mode point loads that move in a direction perpendicular to the crack edge is considered.The analytic expression for the combined mode stress intensity factors as a function of time for any point along the crack edge is obtained.The method of solution is based on the application of integral transform together with the Wiener-Hopf technique and the Cagniard-de Hoop method. Some features of the solution are discussed and graphical results for various point load speeds are presented.
基金provided by the National Natural Science Foundation of China (Nos. 51274204 and 51134025)National Key Basic Research Program (No. 2010CB732002)The Ministry of Education Program for New Century Excellent Talents to Support Project of China (No. NCET-12-0965)
文摘Using digital laser dynamic caustics experimental system and conducting simulation experiment researched the influence rule of blasting excavation of a new roadway on neighboring existed different cross-section roadways. The experimental results show that the influence of blast load on adjacent roadway has a good relationship with the cross-section of roadway. The expansion distance of precrack existed in circular, arch-wall, rectangular roadway is respectively 1.76, 1.61 and 0 cm under blast load.At the same time, the direct-blast side of rectangular roadway has more obvious damage compared with circular and arch-wall roadway. It explains that plane reflects more stress wave than arc, so that it exerts more tensile failure in the direct-blast side, which leads to less stress wave diffracting to the precrack in the back-blast side. When the precrack extends, higher value dynamic stress intensity factor in circular roadway works longer than that of arch-wall roadway. Indirectly, it explains that plane's weakening function on stress wave is significantly stronger than arc. Stress wave brings about self-evident influence on the upper and bottom endpoints of the rectangular roadway, and it respectively extends 1.03, 2.06 cm along the line link direction of the center of the blasthole and the upper and bottom endpoints on the right wall.
基金supported by the National Natural Science Foundation of China(No.51108113)
文摘The present paper is exposed theoretically to the influence on the dynamic stress intensity factor (DSIF) in the piezoelectric bi-materials model with two symmet- rically permeable interracial cracks near the edges of a circular cavity, subjected to the dynamic incident anti-plane shearing wave (SH-wave). An available theoretical method to dynamic analysis in the related research field is provided. The formulations are based on Green's function method. The DSIFs at the inner and outer tips of the left crack are obtained by solving the boundary value problems with the conjunction and crack- simulation technique. The numerical results are obtained by the FORTRAN language program and plotted to show the influence of the variations of the physical parameters, the structural geometry, and the wave frequencies of incident wave on the dimensionless DSIFs. Comparisons with previous work and between the inner and outer tips are con- cluded.
基金The project supported by the National Natural Science Foundation of China (50579081)the Australian Research Council (DP0452681)The English text was polished by Keren Wang
文摘The scaled boundary finite element method (SBFEM) is a recently developed numerical method combining advantages of both finite element methods (FEM) and boundary element methods (BEM) and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors (SIFs) directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.
基金Supported by the Opening Fund of the State Key Laboratory of Explosion Science and Technology of Beijing Institute of Technology(KFJJ17-01M)
文摘In this paper,dynamic caustic method is applied to analyze the blast-induced crack propagation and distribution of the dynamic stress field around an empty circular hole in polymethyl methacrylate(PMMA)material under delayed ignition blasting loads.The following experimental results are obtained.(1)In directional-fracture-controlled blasting,the dynamic stress intensity factors(DSIFs)and the propagation paths of the blast-induced cracks are obviously influenced by the delayed ignition.(2) The circular hole situated between the two boreholes poses a strong guiding effect on the coelesence of the cracks,causing them to propagate towards each other when cracks are reaching the circular hole area.(3)Blast-induced cracks are not initiated preferentially because of the superimposed effect from the explosive stress waves on the cracking area.(4) By using the scanning electron microscopy(SEM)method,it is verified that the roughness of crack surfaces changes along the crack propagation paths.
文摘The scattering of shear waves (SH waves) by nano-scale arbitrary shape inclusion in infinite plane is studied by complex variable function theory. Firstly, the governing equation and the relationships between stress and displacement are given by classical elastic theory. Secondly, the arbitrary shape inclusion in the two-dimensional plane is transformed into a unit circle domain by conformal mapping, the incident wave field and the scattered wave field are presented. Next, the stress and displacement boundary conditions are established by considering surface elasticity theory, The infinite algebraic equations for solving the unknown coefficients of the scattered and standing waves are obtained. Finally, the influence of surface effect, non-dimensional wave number, Shear modulus and hole curvature on the dynamic stress concentration factor are analyzed by some examples, the numerical results show that the surface effect weakens the dynamic stress concentration. With the increase of wave number, the dynamic stress concentration factor (DSCF) decreases. Shear modulus and hole curvature have significant effects on DSCF.
基金supported by National Natural Science Foundation of China(No.91016026)Henan Province Natural Science Foundation Subsidy Project(No.152300410003)
文摘Three dimensional dynamic stress intensity factors are analyzed for a curved crack with a second order perturbation method. The method is extended to obtain an approximate representation of a three dimensional dynamic stress intensity factors at the tip of a curved crack. Due to three dimensional curved crack growth the dynamic energy release rate can be calculated by using the Irwin's formula. A three dimensional curved crack in materials with inhomogeneous fracture toughness are considered. Paths of a brittle three dimensional curved crack propagating along a welded joint are predicted via the present method, where the effects of dynamic applied stresses, residual stresses, and material deterioration due to welding are taken into considerations.
基金the National Science Foundation for Excellent Young Investigators(10325208)the National Natural Science Foundation of China(10432030)the China Postdoctoral Science Foundation(2004036018)
文摘The transient fracture behavior of a functionally graded layered structure subjected to an in-plane impact load is investigated. The studied structure is composed of two homogeneous layers and a functionally gradedinterlayer with a crack perpendicular to the boundaries. The impact load is applied on the face of the crack. Fourier transform and Laplace transform methods are used to formulate the present problem in terms of a singular integral equation in Laplace transform domain. Considering variations of parameters such as the nonhomogeneity constant, the thickness ratio and the crack length, the dynamic stress intensity factors (DSIFs) in time domain are studied and some meaningful conclusions are obtained.
基金Project supported by the Research Fund for Doctors of Hebei Province China (No. B2001213).
文摘The scattering of SH wave by a cylindrical piezoelectric inclusion partially debonded from its surrounding piezoelectric material is investigated using the wave function expansion method and singular integral equation technique. The debonding regions are modeled as mul- tiple arc-shaped interface cracks with non-contacting faces. By expressing the scattered ?elds as wave function expansions with unknown coe?cients, the mixed boundary value problem is ?rstly reduced to a set of simultaneous dual series equations. Then dislocation density functions are introduced as unknowns to transform these dual series equations into Cauchy singular integral equations of the ?rst type, which can be numerically solved easily. The solution is valid for arbi- trary number and size of the debonds. Finally, numerical results of the dynamic stress intensity factors are presented for the cases of one debond and two debonds. The e?ects of incidence direc- tion, crack con?guration and various material parameters on the dynamic stress intensity factors are respectively discussed. The solution of this problem is expected to ?nd applications in the investigation of dynamic fracture properties of piezoelectric materials with cracks.
基金Project supported by the National Natural Science Foundation of China(Nos.11272222,11502108,and 11611530686)the Natural Science Foundation for Distinguished Young Scholars of Jiangsu Province of China(No.BK20140037)
文摘Due to the incompatibility of the interlaminar deformations,the interface debonding or cracking usually happens in a layered magnetoelectric(ME)structure under an applied load.In this paper,the transient responses of the anti-plane interface cracks in piezoelectric(PE)-piezomagnetic(PM)sandwich structures are studied by the standard methods of the integral transform and singular integral equation.Discussion on the numerical examples indicates that the PE-PM-PE structure under electric impact is more likely to fracture than the PM-PE-PM structure under a magnetic impact.The dynamic stress intensity factors(DSIFs)are more sensitive to the variation of the active layer thickness.The effects of the material constants on the DSIFs are dependent on the roles played by PE and PM media during the deformation process.
基金supported by the National Natural Science Foundation of China (Grant 11372091)China Scholarship Council (Grant 201406120086)
文摘A mathematical formulation is presented for the dynamic stress intensity factor (mode I) of a finite permeable crack subjected to a time-harmonic propagating longitudinal wave in an infinite poroelastic solid. In particular, the effect of the wave-induced fluid flow due to the presence of a liquid-saturated crack on the dynamic stress intensity factor is analyzed. Fourier sine and cosine integral transforms in conjunction with Helmholtz potential theory are used to formulate the mixed boundary-value problem as dual integral equations in the frequency domain. The dual integral equations are reduced to a Fredholm integral equation of the second kind. It is found that the stress intensity factor monotonically decreases with increasing frequency, decreasing the fastest when the crack width and the slow wave wavelength are of the same order. The characteristic frequency at which the stress intensity factor decays the fastest shifts to higher frequency values when the crack width decreases.