The energy norm convergence rate of the finite element solution of the heat equation is reduced by the time-regularity of the exact solution. This paper presents an adaptive finite element treatment of time-dependent ...The energy norm convergence rate of the finite element solution of the heat equation is reduced by the time-regularity of the exact solution. This paper presents an adaptive finite element treatment of time-dependent singularities on the one-dimensional heat equation. The method is based on a Fourier decomposition of the solution and an extraction formula of the coefficients of the singularities coupled with a predictor-corrector algorithm. The method recovers the optimal convergence rate of the finite element method on a quasi-uniform mesh refinement. Numerical results are carried out to show the efficiency of the method.展开更多
The multi-variable finite element algorithm based on the generalized Gulerkin's method is more flexible to establish a finite element model in the continuum mechanics. By using this algorithm and numerical tests a...The multi-variable finite element algorithm based on the generalized Gulerkin's method is more flexible to establish a finite element model in the continuum mechanics. By using this algorithm and numerical tests a new singular finite element for elasto-plastic fracture analysis has been formulated. The results of numerical tests show that the new element possesses high accuracy and good performance. Some rules for formulating a multi-variable singular finite element are also discussed in this paper.展开更多
Using the method of the boundary integral equation, a set of singular integral equations of the hear transfer problems and the thermo-elastic problems of a crack embedded in a two-dimensional finite body is derived, a...Using the method of the boundary integral equation, a set of singular integral equations of the hear transfer problems and the thermo-elastic problems of a crack embedded in a two-dimensional finite body is derived, and then,its numerical method is proposed by the numerical method of the singular integral equations combined with boundary element method. Moreover, the singular nature of temperature gradient field near the crack front is proved by the main-part analysis method of the singular integral equation, and the singular temperature gradients are exactly obtained. Finally, several typical examples calculated.展开更多
Quadrature rules for evaluating singular integrals that typically occur in the boundary element method (BEM) for two-dimensional and axisymmetric three-dimensional problems are considered. This paper focuses on the nu...Quadrature rules for evaluating singular integrals that typically occur in the boundary element method (BEM) for two-dimensional and axisymmetric three-dimensional problems are considered. This paper focuses on the numerical integration of the functions on the standard domain [-1, 1], with a logarithmic singularity at the centre. The substitution x = tp, where p (≥ 3) is an odd integer is given particular attention, as this returns a regular integral and the domain unchanged. Gauss-Legendre quadrature rules are applied to the transformed integrals for a number of values of p. It is shown that a high value for p typically gives more accurate results.展开更多
A general algorithm is applied to the regularization of nearly singular integrals in the boundary element method of planar potential problems. For linear elements, the strongly singular and hypersingular integrals of ...A general algorithm is applied to the regularization of nearly singular integrals in the boundary element method of planar potential problems. For linear elements, the strongly singular and hypersingular integrals of the interior points very close to boundary were categorized into two forms. The factor leading to the singularity was transformed out of the integral representations with integration by parts, so non-singular regularized formulas were presented for the two forms of integrals. Furthermore, quadratic elements are used in addition to linear ones. The quadratic element very close to the internal point can be divided into two linear ones, so that the algorithm is still valid. Numerical examples demonstrate the effectiveness and accuracy of this algorithm. Especially for problems with curved boundaries, the combination of quadratic elements and linear elements can give more accurate results.展开更多
A singularly perturbed advection-diffusion two-point Robin boundary value problem whose solution has a single boundary layer is considered. Based on the piecewise linear polynomial approximation, the finite element me...A singularly perturbed advection-diffusion two-point Robin boundary value problem whose solution has a single boundary layer is considered. Based on the piecewise linear polynomial approximation, the finite element method is applied to the problem. Estimation of the error between solution and the finite element approximation is given in energy norm on shishkin-type mesh.展开更多
The basic principle and numerical technique for simulating two three-dimensional bubbles near a free surface are studied in detail by using boundary element method. The singularities of influence coefficient matrix ar...The basic principle and numerical technique for simulating two three-dimensional bubbles near a free surface are studied in detail by using boundary element method. The singularities of influence coefficient matrix are eliminated using coordinate transformation and so-called 4π rule. The solid angle for the open surface is treated in direct method based on its definition. Several kinds of configurations for the bubbles and free surface have been investigated. The pressure contours during the evolution of bubbles are obtained in our model and can better illuminate the mechanism underlying the motions of bubbles and free surface. The bubble dynamics and their interactions have close relation with the standoff distances, buoyancy parameters and initial sizes of bubbles. Completely different bubble shapes, free surface motions, jetting patterns and pressure distributions under different parameters can be observed in our model, as demon- strated in our calculation results.展开更多
In this paper, we shall give a group of new isoparametric elements suitable .for St.Wnant's torsion of a bar with vertical cracks. These elements are eight-pointisoparametric element, quarter eight-Point isoparame...In this paper, we shall give a group of new isoparametric elements suitable .for St.Wnant's torsion of a bar with vertical cracks. These elements are eight-pointisoparametric element, quarter eight-Point isoparametric element and eight-pointisoparametric transition element. Among these elements, the second and the thirdelements possess the singularity of order r-1/2 at crack tip. Using these elements, wehave completed the calculations of St. Venant's torsion for a cylinder with a radialvertical crack. The calculated results show that the isoparametric elements given bythis paper have ideal accuracy. good convergence. high speed of convergence, lowfreedom and little computational time, and so they can be widely applied to practice.展开更多
In this article, a direct stress approach based on finite element analysis to determine the stress intensity factor is improved. Firstly, by comparing the rigorous solution against the asymptotic solution for a proble...In this article, a direct stress approach based on finite element analysis to determine the stress intensity factor is improved. Firstly, by comparing the rigorous solution against the asymptotic solution for a problem of an infinite plate embedded a central crack, we found that the stresses in a restrictive interval near the crack tip given by the rigorous solution can be used to determine the stress intensity factor, which is nearly equal to the stress intensity factor given by the asymptotic solution. Secondly, the crack problem is solved numerically by the finite element method. Depending on the modeling capability of the software, we designed an adaptive mesh model to simulate the stress singularity. Thus, the stress result in an appropriate interval near the crack tip is fairly approximated to the rigorous solution of the corresponding crack problem. Therefore, the stress intensity factor may be calculated from the stress distribution in the appropriate interval, with a high accuracy.展开更多
Although boundary displacement and traction are independent field variables in boundary conditions of an elasticity problem at a non-singular boundary point, there exist definite relations of singularity intensities b...Although boundary displacement and traction are independent field variables in boundary conditions of an elasticity problem at a non-singular boundary point, there exist definite relations of singularity intensities between boundary displacement derivatives and tractions at a singular boundary point. The analytical forms of the relations at a singular smooth point for 2D isotropic elastic problems have been established in this work. By using the relations, positions of the singular boundary points and the corresponding singularity intensities of the unknown boundary field variables can be determined a priori. Therefore, more appropriate shape functions of the unknown boundary field variables in singular elements can be constructed. A numerical example shows that the accuracy of the BEM analysis using the developed theory is greatly increased.展开更多
The paper discusses how to reduce higher singularity order of a boundary integral equation. The approach will be discussed in some detail for plane elasticity.Numerical results for the meshes of unequal length boundar...The paper discusses how to reduce higher singularity order of a boundary integral equation. The approach will be discussed in some detail for plane elasticity.Numerical results for the meshes of unequal length boundary elements are reported.Higher precision for both deflection and force is obtained than that obtained with a general boundary element method.展开更多
In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body...In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body-fitted meshes are used.For homogeneous jump conditions,both non-conforming and conforming basis functions are constructed in such a way that they satisfy the natural jump conditions. For non-homogeneous jump conditions,a pair of functions that satisfy the same non-homogeneous jump conditions are constructed using a level-set representation of the interface.With such a pair of functions,the discontinuities across the interface in the solution and flux are removed;and an equivalent elasticity interface problem with homogeneous jump conditions is formulated.Numerical examples are presented to demonstrate that such methods have second order convergence.展开更多
By using the concept of finite-part integral, a set of hypersingular integro-differential equations for multiple interracial cracks in a three-dimensional infinite bimaterial subjected to arbitrary loads is derived. I...By using the concept of finite-part integral, a set of hypersingular integro-differential equations for multiple interracial cracks in a three-dimensional infinite bimaterial subjected to arbitrary loads is derived. In the numerical analysis, unknown displacement discontinuities are approximated with the products of the fundamental density functions and power series. The fundamental functions are chosen to express a two-dimensional interface crack rigorously. As illustrative examples, the stress intensity factors for two rectangular interface cracks are calculated for various spacing, crack shape and elastic constants. It is shown that the stress intensity factors decrease with the crack spacing.展开更多
WT5”BZ]In this paper,we have introudced the notion of fuzzy semi-(inf-)singular,and dicussed their algebraic structure and properties.Mainly,we use these concepts and properties to research cancellative properties of...WT5”BZ]In this paper,we have introudced the notion of fuzzy semi-(inf-)singular,and dicussed their algebraic structure and properties.Mainly,we use these concepts and properties to research cancellative properties of fuzzy languages. [WT5”HZ]展开更多
The uniqueness of solution of field point, inside a convex region due to singular source(s) with kernel function decreasing with distance increasing, outside-region-distribution(s) such that the boundary condition exp...The uniqueness of solution of field point, inside a convex region due to singular source(s) with kernel function decreasing with distance increasing, outside-region-distribution(s) such that the boundary condition expressed by the response of the source(s) is satisfied, is proved by using the condition of kernel function decreasing with distance increasing anal an integral inequality. Examples of part of these singular sources such as Kelvin's point force, Point-Ring-Couple (PRC) etc. are given. The proof of uniqueness of solution of field point in a twisted shaft of revolution due to PRC distribution is given as an example of application.展开更多
文摘The energy norm convergence rate of the finite element solution of the heat equation is reduced by the time-regularity of the exact solution. This paper presents an adaptive finite element treatment of time-dependent singularities on the one-dimensional heat equation. The method is based on a Fourier decomposition of the solution and an extraction formula of the coefficients of the singularities coupled with a predictor-corrector algorithm. The method recovers the optimal convergence rate of the finite element method on a quasi-uniform mesh refinement. Numerical results are carried out to show the efficiency of the method.
文摘The multi-variable finite element algorithm based on the generalized Gulerkin's method is more flexible to establish a finite element model in the continuum mechanics. By using this algorithm and numerical tests a new singular finite element for elasto-plastic fracture analysis has been formulated. The results of numerical tests show that the new element possesses high accuracy and good performance. Some rules for formulating a multi-variable singular finite element are also discussed in this paper.
文摘Using the method of the boundary integral equation, a set of singular integral equations of the hear transfer problems and the thermo-elastic problems of a crack embedded in a two-dimensional finite body is derived, and then,its numerical method is proposed by the numerical method of the singular integral equations combined with boundary element method. Moreover, the singular nature of temperature gradient field near the crack front is proved by the main-part analysis method of the singular integral equation, and the singular temperature gradients are exactly obtained. Finally, several typical examples calculated.
文摘Quadrature rules for evaluating singular integrals that typically occur in the boundary element method (BEM) for two-dimensional and axisymmetric three-dimensional problems are considered. This paper focuses on the numerical integration of the functions on the standard domain [-1, 1], with a logarithmic singularity at the centre. The substitution x = tp, where p (≥ 3) is an odd integer is given particular attention, as this returns a regular integral and the domain unchanged. Gauss-Legendre quadrature rules are applied to the transformed integrals for a number of values of p. It is shown that a high value for p typically gives more accurate results.
文摘A general algorithm is applied to the regularization of nearly singular integrals in the boundary element method of planar potential problems. For linear elements, the strongly singular and hypersingular integrals of the interior points very close to boundary were categorized into two forms. The factor leading to the singularity was transformed out of the integral representations with integration by parts, so non-singular regularized formulas were presented for the two forms of integrals. Furthermore, quadratic elements are used in addition to linear ones. The quadratic element very close to the internal point can be divided into two linear ones, so that the algorithm is still valid. Numerical examples demonstrate the effectiveness and accuracy of this algorithm. Especially for problems with curved boundaries, the combination of quadratic elements and linear elements can give more accurate results.
文摘A singularly perturbed advection-diffusion two-point Robin boundary value problem whose solution has a single boundary layer is considered. Based on the piecewise linear polynomial approximation, the finite element method is applied to the problem. Estimation of the error between solution and the finite element approximation is given in energy norm on shishkin-type mesh.
基金supported by the Funds for Creative Research Groups of China (50921001)the State Key Development Program for Basic Research of China (2010CB832704)
文摘The basic principle and numerical technique for simulating two three-dimensional bubbles near a free surface are studied in detail by using boundary element method. The singularities of influence coefficient matrix are eliminated using coordinate transformation and so-called 4π rule. The solid angle for the open surface is treated in direct method based on its definition. Several kinds of configurations for the bubbles and free surface have been investigated. The pressure contours during the evolution of bubbles are obtained in our model and can better illuminate the mechanism underlying the motions of bubbles and free surface. The bubble dynamics and their interactions have close relation with the standoff distances, buoyancy parameters and initial sizes of bubbles. Completely different bubble shapes, free surface motions, jetting patterns and pressure distributions under different parameters can be observed in our model, as demon- strated in our calculation results.
文摘In this paper, we shall give a group of new isoparametric elements suitable .for St.Wnant's torsion of a bar with vertical cracks. These elements are eight-pointisoparametric element, quarter eight-Point isoparametric element and eight-pointisoparametric transition element. Among these elements, the second and the thirdelements possess the singularity of order r-1/2 at crack tip. Using these elements, wehave completed the calculations of St. Venant's torsion for a cylinder with a radialvertical crack. The calculated results show that the isoparametric elements given bythis paper have ideal accuracy. good convergence. high speed of convergence, lowfreedom and little computational time, and so they can be widely applied to practice.
基金financial support of the National Natural Science Foundation of China (Grant 11572226)
文摘In this article, a direct stress approach based on finite element analysis to determine the stress intensity factor is improved. Firstly, by comparing the rigorous solution against the asymptotic solution for a problem of an infinite plate embedded a central crack, we found that the stresses in a restrictive interval near the crack tip given by the rigorous solution can be used to determine the stress intensity factor, which is nearly equal to the stress intensity factor given by the asymptotic solution. Secondly, the crack problem is solved numerically by the finite element method. Depending on the modeling capability of the software, we designed an adaptive mesh model to simulate the stress singularity. Thus, the stress result in an appropriate interval near the crack tip is fairly approximated to the rigorous solution of the corresponding crack problem. Therefore, the stress intensity factor may be calculated from the stress distribution in the appropriate interval, with a high accuracy.
文摘Although boundary displacement and traction are independent field variables in boundary conditions of an elasticity problem at a non-singular boundary point, there exist definite relations of singularity intensities between boundary displacement derivatives and tractions at a singular boundary point. The analytical forms of the relations at a singular smooth point for 2D isotropic elastic problems have been established in this work. By using the relations, positions of the singular boundary points and the corresponding singularity intensities of the unknown boundary field variables can be determined a priori. Therefore, more appropriate shape functions of the unknown boundary field variables in singular elements can be constructed. A numerical example shows that the accuracy of the BEM analysis using the developed theory is greatly increased.
基金The project supported by National Natural Science Foundation of China.
文摘The paper discusses how to reduce higher singularity order of a boundary integral equation. The approach will be discussed in some detail for plane elasticity.Numerical results for the meshes of unequal length boundary elements are reported.Higher precision for both deflection and force is obtained than that obtained with a general boundary element method.
基金supported by the US ARO grants 49308-MA and 56349-MAthe US AFSOR grant FA9550-06-1-024+1 种基金he US NSF grant DMS-0911434the State Key Laboratory of Scientific and Engineering Computing of Chinese Academy of Sciences during a visit by Z.Li between July-August,2008.
文摘In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body-fitted meshes are used.For homogeneous jump conditions,both non-conforming and conforming basis functions are constructed in such a way that they satisfy the natural jump conditions. For non-homogeneous jump conditions,a pair of functions that satisfy the same non-homogeneous jump conditions are constructed using a level-set representation of the interface.With such a pair of functions,the discontinuities across the interface in the solution and flux are removed;and an equivalent elasticity interface problem with homogeneous jump conditions is formulated.Numerical examples are presented to demonstrate that such methods have second order convergence.
基金supported by the National Natural Science Foundation of China (No. 10872213)
文摘By using the concept of finite-part integral, a set of hypersingular integro-differential equations for multiple interracial cracks in a three-dimensional infinite bimaterial subjected to arbitrary loads is derived. In the numerical analysis, unknown displacement discontinuities are approximated with the products of the fundamental density functions and power series. The fundamental functions are chosen to express a two-dimensional interface crack rigorously. As illustrative examples, the stress intensity factors for two rectangular interface cracks are calculated for various spacing, crack shape and elastic constants. It is shown that the stress intensity factors decrease with the crack spacing.
文摘WT5”BZ]In this paper,we have introudced the notion of fuzzy semi-(inf-)singular,and dicussed their algebraic structure and properties.Mainly,we use these concepts and properties to research cancellative properties of fuzzy languages. [WT5”HZ]
文摘The uniqueness of solution of field point, inside a convex region due to singular source(s) with kernel function decreasing with distance increasing, outside-region-distribution(s) such that the boundary condition expressed by the response of the source(s) is satisfied, is proved by using the condition of kernel function decreasing with distance increasing anal an integral inequality. Examples of part of these singular sources such as Kelvin's point force, Point-Ring-Couple (PRC) etc. are given. The proof of uniqueness of solution of field point in a twisted shaft of revolution due to PRC distribution is given as an example of application.