By using the analytic theory of hypersingular integral equations in three- dimensional fracture mechanics, the interactions between two parallel planar cracks under arbitrary loads are investigated. According to the c...By using the analytic theory of hypersingular integral equations in three- dimensional fracture mechanics, the interactions between two parallel planar cracks under arbitrary loads are investigated. According to the concepts and method of finite- part integrals, a set of hypersingular integral equations is derived, in which the unknown functions are the displacement discontinuities of the crack surfaces. Then its numerical method is proposed by combining the finite-part integral method with the boundary element method. Based on the above results, the method for calculating the stress intensity factors with the displacement discontinuities of the crack surfaces is presented. Finally, several typical examples are calculated and the numerical results are satisfactory.展开更多
Based on the fact that the singular boundary integrals in the sense of Cauchy principal value can be represented approximately by the mean values of two companion nearly singular boundary integrals, a vary general app...Based on the fact that the singular boundary integrals in the sense of Cauchy principal value can be represented approximately by the mean values of two companion nearly singular boundary integrals, a vary general approach was developed in the paper. In the approach, the approximate formulation before discretization was constructed to cope with the difficulties encountered in the corner treatment in the formulations of hypersingular boundary integral equations. This makes it possible to solve the hypersingular boundary integral equation numerically in a non regularized form and in a local manner by using conforming C 0 quadratic boundary elements and standard Gaussian quadratures similar to those employed in the conventional displacement BIE formulations. The approximate formulation is very convenient to use because the corner information is comprised naturally in the representations of those approximate integrals. Numerical examples in plane elasticity show that with the present approach, the compatible or better results can be achieved in comparison with those of the conventional BIE formulations.展开更多
The aim of this paper is to investigate the numerical solution of the hypersingular integral equation reduced by the harmonic equation. First, we transform the hypersingular integral equation into 2π-periodic hypersi...The aim of this paper is to investigate the numerical solution of the hypersingular integral equation reduced by the harmonic equation. First, we transform the hypersingular integral equation into 2π-periodic hypersingular integral equation with the map x=cot(θ/2). Second, we initiate the study of the multiscale Galerkin method for the 2π-periodic hypersingular integral equation. The trigonometric wavelets are used as trial functions. Consequently, the 2j+1 × 2j+1 stiffness matrix Kj can be partitioned j×j block matrices. Furthermore, these block matrices are zeros except main diagonal block matrices. These main diagonal block matrices are symmetrical and circulant matrices, and hence the solution of the associated linear algebraic system can be solved with the fast Fourier transform and the inverse fast Fourier transform instead of the inverse matrix. Finally, we provide several numerical examples to demonstrate our method has good accuracy even though the exact solutions are multi-peak and almost singular.展开更多
Using the Somigliana formula and the concepts of finite-part integral, a set of hypersingular integral equations to solve the arbitrary flat crack in three-dimensional elasticity is derived and its numerical method is...Using the Somigliana formula and the concepts of finite-part integral, a set of hypersingular integral equations to solve the arbitrary flat crack in three-dimensional elasticity is derived and its numerical method is then proposed by combining the finite-part integral method with boundary element method. In order to verify the method, several numerical examples are carried out. The results of the displacement discontinuities of the crack surface and the stress intensity factors at the crack front are in good agrement with the theoretical solutions.展开更多
A prestressed elastic medium containing a mode-Ⅲcrack is studied by means of the couple stress theory(CST).Based on the CST under initial stresses,a governing differential equation along with a mixed boundary value p...A prestressed elastic medium containing a mode-Ⅲcrack is studied by means of the couple stress theory(CST).Based on the CST under initial stresses,a governing differential equation along with a mixed boundary value problem is established.The singularities of the couple stress and force stress near the crack tips are analyzed through the asymptotic crack-tip fields resulting from the characteristic expansion method.To determine their intensity,a hypersingular integral equation is derived and numerically solved with the help of the Chebyshev polynomial.The obtained results show a strong size-dependence of the out-of-plane displacement on the crack and the couple stress intensity factor(CSIF)and the force stress intensity factor(FSIF)around the crack tips.The symmetric part of the shear stress has no singularity,and the skew-symmetric part related to the couple stress exhibits an r^(-3/2)singularity,in which r is the distance from the crack tip.The initial stresses also affect the crack tearing displacement and the CSIF and FSIF.展开更多
In this paper we have investigated the reflection and the transmission of a system of two symmetric circular-arc-shaped thin porous plates submerged in deep water within the context of linear theory. The hypersingular...In this paper we have investigated the reflection and the transmission of a system of two symmetric circular-arc-shaped thin porous plates submerged in deep water within the context of linear theory. The hypersingular integral equation technique has been used to analyze the problem mathematically. The integral equations are formulated by applying Green's integral theorem to the fundamental potential function and the scattered potential function into a suitable fluid region, and then using the boundary condition on the porous plate surface. These are solved approximately using an expansion-cure-collocation method where the behaviour of the potential functions at the tips of the plates have been used. This method ultimately produces a very good numerical approximation for the reflection and the transmission coefficients and hydrodynamic force components. The numerical results are depicted graphically against the wave number for a variety of layouts of the arc. Some results are compared with known results for similar configurations of dual rigid plate systems available in the literature with good agreement.展开更多
A new type of dual boundary integral equations(DBIE)is presented first,through which,a smaller system of equations needs to be solved in fracture analysis.Then a non-conforming crack tip element in two-dimensional pro...A new type of dual boundary integral equations(DBIE)is presented first,through which,a smaller system of equations needs to be solved in fracture analysis.Then a non-conforming crack tip element in two-dimensional problems is proposed.The exact formula for the hypersingular integral over the non-con- forming crack tip element is given next.By virtue of Green's-function-library strategy,a series of stress in- tensity factors(SIF)of different crack orientations,locations and/or sizes in a complicated structure can be obtained easily and efficiently.Finally,several examples of fracture analysis in two dimensions are given to demonstrate the accuracy and efficiency of the method proposed.展开更多
According to the constitutive relationship in linear piezoceramics, elliptical crack problems in the impermeable case are reconsidered with the hypersingular integral equation method. Unknown displacement and electric...According to the constitutive relationship in linear piezoceramics, elliptical crack problems in the impermeable case are reconsidered with the hypersingular integral equation method. Unknown displacement and electric potential jumps in the integral equations are approximated with a product of the fundamental density function and polynomials, in which the fundamental density function reflects the singular behavior of electroelastic fields near the crack front and the polynomials can be reduced to a real constant under uniform loading. Ellipsoidal coordinates are cleverly introduced to solve the unknown displacement and electric potential jumps in the integral equations under uniform loading. With the help of these solutions and definitions of electroelastic field intensity factors, exact expressions for mode Ⅰ, mode Ⅱ and mode Ⅲ stress intensity factors as well as the mode Ⅳ electric displacement intensity factor are obtained. The present results under uniform normal loading are the same as the available exact solutions, but those under uniform shear loading have not been found in the literature as yet.展开更多
The present paper is concerned with scattering of water waves from a vertical plate, modeled as an elastic plate, submerged in deep water covered with a thin uniform sheet of ice. The problem is formulated in terms of...The present paper is concerned with scattering of water waves from a vertical plate, modeled as an elastic plate, submerged in deep water covered with a thin uniform sheet of ice. The problem is formulated in terms of a hypersingular integral equation by a suitable application of Green's integral theorem in terms of difference of potential functicns across the barrier. This integral equation is solved by a collocation method using a finite series involving Chebyshev polynomials. Reflection and transmission coefficients are obtained numerically and presented graphically for various values of the wave number and ice-cover parameter.展开更多
Formulation in terms of hypersingular integral equations for the interaction between straight and curved cracks in plane elasticity is obtained using the complex variable functions method. The curved length coordinate...Formulation in terms of hypersingular integral equations for the interaction between straight and curved cracks in plane elasticity is obtained using the complex variable functions method. The curved length coordinate method and a suitable numerical scheme are used to solve such integrals numerically for the unknown function, which are later used to find the stress intensity factor, SIF.展开更多
In this paper, we study the behavior of the solution at the crack edges for a nearly circular crack with developing cusps subject to shear loading. The problem of finding the resulting force can be written in the form...In this paper, we study the behavior of the solution at the crack edges for a nearly circular crack with developing cusps subject to shear loading. The problem of finding the resulting force can be written in the form of a hypersingular integral equation. The equation is then trans-formed into a similar equation over a circular region using conformal mapping. The equation is solved numerically for the unknown coefficients, which will later be used in finding the stress intensity factors. The sliding and tearing mode stress intensity factors are evaluated for cracks and displayed graphically. Our results seem to agree with the existing asymptotic solution.展开更多
We present in this paper a numerical method for hypersingular boundary integral equations. This method was developed for planar crack problems: additional edge singularities are known to develop in that case. This pa...We present in this paper a numerical method for hypersingular boundary integral equations. This method was developed for planar crack problems: additional edge singularities are known to develop in that case. This paper includes a rigorous error analysis proving the convergence of our numerical scheme. Three types of examples are covered: the Laplace equation in free space, the linear elasticity equation in free space, and in half space.展开更多
文摘By using the analytic theory of hypersingular integral equations in three- dimensional fracture mechanics, the interactions between two parallel planar cracks under arbitrary loads are investigated. According to the concepts and method of finite- part integrals, a set of hypersingular integral equations is derived, in which the unknown functions are the displacement discontinuities of the crack surfaces. Then its numerical method is proposed by combining the finite-part integral method with the boundary element method. Based on the above results, the method for calculating the stress intensity factors with the displacement discontinuities of the crack surfaces is presented. Finally, several typical examples are calculated and the numerical results are satisfactory.
文摘Based on the fact that the singular boundary integrals in the sense of Cauchy principal value can be represented approximately by the mean values of two companion nearly singular boundary integrals, a vary general approach was developed in the paper. In the approach, the approximate formulation before discretization was constructed to cope with the difficulties encountered in the corner treatment in the formulations of hypersingular boundary integral equations. This makes it possible to solve the hypersingular boundary integral equation numerically in a non regularized form and in a local manner by using conforming C 0 quadratic boundary elements and standard Gaussian quadratures similar to those employed in the conventional displacement BIE formulations. The approximate formulation is very convenient to use because the corner information is comprised naturally in the representations of those approximate integrals. Numerical examples in plane elasticity show that with the present approach, the compatible or better results can be achieved in comparison with those of the conventional BIE formulations.
文摘The aim of this paper is to investigate the numerical solution of the hypersingular integral equation reduced by the harmonic equation. First, we transform the hypersingular integral equation into 2π-periodic hypersingular integral equation with the map x=cot(θ/2). Second, we initiate the study of the multiscale Galerkin method for the 2π-periodic hypersingular integral equation. The trigonometric wavelets are used as trial functions. Consequently, the 2j+1 × 2j+1 stiffness matrix Kj can be partitioned j×j block matrices. Furthermore, these block matrices are zeros except main diagonal block matrices. These main diagonal block matrices are symmetrical and circulant matrices, and hence the solution of the associated linear algebraic system can be solved with the fast Fourier transform and the inverse fast Fourier transform instead of the inverse matrix. Finally, we provide several numerical examples to demonstrate our method has good accuracy even though the exact solutions are multi-peak and almost singular.
基金The project supported by the Foundation of the Science Research of the State Education Commission of the People's Republic of China
文摘Using the Somigliana formula and the concepts of finite-part integral, a set of hypersingular integral equations to solve the arbitrary flat crack in three-dimensional elasticity is derived and its numerical method is then proposed by combining the finite-part integral method with boundary element method. In order to verify the method, several numerical examples are carried out. The results of the displacement discontinuities of the crack surface and the stress intensity factors at the crack front are in good agrement with the theoretical solutions.
基金Project supported by the National Natural Science Foundation of China(Nos.11672336,12072374)。
文摘A prestressed elastic medium containing a mode-Ⅲcrack is studied by means of the couple stress theory(CST).Based on the CST under initial stresses,a governing differential equation along with a mixed boundary value problem is established.The singularities of the couple stress and force stress near the crack tips are analyzed through the asymptotic crack-tip fields resulting from the characteristic expansion method.To determine their intensity,a hypersingular integral equation is derived and numerically solved with the help of the Chebyshev polynomial.The obtained results show a strong size-dependence of the out-of-plane displacement on the crack and the couple stress intensity factor(CSIF)and the force stress intensity factor(FSIF)around the crack tips.The symmetric part of the shear stress has no singularity,and the skew-symmetric part related to the couple stress exhibits an r^(-3/2)singularity,in which r is the distance from the crack tip.The initial stresses also affect the crack tearing displacement and the CSIF and FSIF.
基金Partially Supported by the Department of Science and Technology Through a Research Grant to RG(No.SR/FTP/MS-020/2010)
文摘In this paper we have investigated the reflection and the transmission of a system of two symmetric circular-arc-shaped thin porous plates submerged in deep water within the context of linear theory. The hypersingular integral equation technique has been used to analyze the problem mathematically. The integral equations are formulated by applying Green's integral theorem to the fundamental potential function and the scattered potential function into a suitable fluid region, and then using the boundary condition on the porous plate surface. These are solved approximately using an expansion-cure-collocation method where the behaviour of the potential functions at the tips of the plates have been used. This method ultimately produces a very good numerical approximation for the reflection and the transmission coefficients and hydrodynamic force components. The numerical results are depicted graphically against the wave number for a variety of layouts of the arc. Some results are compared with known results for similar configurations of dual rigid plate systems available in the literature with good agreement.
基金the Aeronautical Science Foundation of China (No.99C53026).
文摘A new type of dual boundary integral equations(DBIE)is presented first,through which,a smaller system of equations needs to be solved in fracture analysis.Then a non-conforming crack tip element in two-dimensional problems is proposed.The exact formula for the hypersingular integral over the non-con- forming crack tip element is given next.By virtue of Green's-function-library strategy,a series of stress in- tensity factors(SIF)of different crack orientations,locations and/or sizes in a complicated structure can be obtained easily and efficiently.Finally,several examples of fracture analysis in two dimensions are given to demonstrate the accuracy and efficiency of the method proposed.
基金Project supported by the Jiangxi Provincial Natural Science Foundation (No.0112001)the Japan Society for the Promotion of Science Postdoctoral Fellowship (No.P01205).
文摘According to the constitutive relationship in linear piezoceramics, elliptical crack problems in the impermeable case are reconsidered with the hypersingular integral equation method. Unknown displacement and electric potential jumps in the integral equations are approximated with a product of the fundamental density function and polynomials, in which the fundamental density function reflects the singular behavior of electroelastic fields near the crack front and the polynomials can be reduced to a real constant under uniform loading. Ellipsoidal coordinates are cleverly introduced to solve the unknown displacement and electric potential jumps in the integral equations under uniform loading. With the help of these solutions and definitions of electroelastic field intensity factors, exact expressions for mode Ⅰ, mode Ⅱ and mode Ⅲ stress intensity factors as well as the mode Ⅳ electric displacement intensity factor are obtained. The present results under uniform normal loading are the same as the available exact solutions, but those under uniform shear loading have not been found in the literature as yet.
基金supported by the Department of Science and Technology of New Delhi (No.SR/SY/MS:521/08)
文摘The present paper is concerned with scattering of water waves from a vertical plate, modeled as an elastic plate, submerged in deep water covered with a thin uniform sheet of ice. The problem is formulated in terms of a hypersingular integral equation by a suitable application of Green's integral theorem in terms of difference of potential functicns across the barrier. This integral equation is solved by a collocation method using a finite series involving Chebyshev polynomials. Reflection and transmission coefficients are obtained numerically and presented graphically for various values of the wave number and ice-cover parameter.
基金Ministry of Science,Technology and Innovation(MOSTI),Malaysia for the Science Fund,Vot No.5450657
文摘Formulation in terms of hypersingular integral equations for the interaction between straight and curved cracks in plane elasticity is obtained using the complex variable functions method. The curved length coordinate method and a suitable numerical scheme are used to solve such integrals numerically for the unknown function, which are later used to find the stress intensity factor, SIF.
基金supported by the Ministry Of Higher Education Malaysia for the Fundamental Research Grant scheme,project No. 01-04-10-897FRthe NSF scholarship
文摘In this paper, we study the behavior of the solution at the crack edges for a nearly circular crack with developing cusps subject to shear loading. The problem of finding the resulting force can be written in the form of a hypersingular integral equation. The equation is then trans-formed into a similar equation over a circular region using conformal mapping. The equation is solved numerically for the unknown coefficients, which will later be used in finding the stress intensity factors. The sliding and tearing mode stress intensity factors are evaluated for cracks and displayed graphically. Our results seem to agree with the existing asymptotic solution.
文摘We present in this paper a numerical method for hypersingular boundary integral equations. This method was developed for planar crack problems: additional edge singularities are known to develop in that case. This paper includes a rigorous error analysis proving the convergence of our numerical scheme. Three types of examples are covered: the Laplace equation in free space, the linear elasticity equation in free space, and in half space.