The compact implicit integration factor （cIIF） method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF me...The compact implicit integration factor （cIIF） method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF method to some complex-valued nonlinear evolutionary equations such as the nonlinear SchrSdinger （NLS） equation and the complex Ginzburg-Landau （GL） equation. Detailed algorithm formulation and practical implementation of cIIF method are performed. The numerical results indicate that this method is very accurate and efficient.展开更多
The stress intensity factors(SIFs)for two-dimensional cracks are extracted using the p-version finite element method(P-FEM)and the contour integral method.Several numerical experiments,e.g.,crack initiating from the e...The stress intensity factors(SIFs)for two-dimensional cracks are extracted using the p-version finite element method(P-FEM)and the contour integral method.Several numerical experiments,e.g.,crack initiating from the edge of a circular hole under an unidirectional uniform tension and two equal-length,unequal-length hole-edge cracks,respectively,at a rectangular plate,an inclined centered crack under uniaxial tension at a square plate and a pipeline crack model,are used to demonstrate the accuracy and effectiveness of the approaches.SIFs are presented for the effects of various crack lengths and length-width ratio.Numerical results are analyzed and compared with reference solutions and results obtained by the Voronoi cell finite element method,boundary element method,high-order extended finite element method(high-order XFEM)and commercial finite element software ABAQUS in the available literature.Numerical results are in good agreement with the benchmark problems and show faster convergence rate,higher accuracy and better numerical stability.展开更多
Using the single crack solution and the regular solution of plane harmonic function, the problem of Saint_Venant bending of a cracked cylinder by a transverse force was reduced to solving two sets of integral equation...Using the single crack solution and the regular solution of plane harmonic function, the problem of Saint_Venant bending of a cracked cylinder by a transverse force was reduced to solving two sets of integral equations and its general solution was then obtained. Based on the obtained solution, a method to calculate the bending center and the stress intensity factors of the cracked cylinger whose cross_section is not thin_walled, but of small torsion rigidity is proposed. Some numerical examples are given.展开更多
In this study,to develop a benefit-allocation model,in-depth analysis of a distributed photovoltaic-powergeneration carport and energy-storage charging-pile project was performed;the model was developed using Shapley ...In this study,to develop a benefit-allocation model,in-depth analysis of a distributed photovoltaic-powergeneration carport and energy-storage charging-pile project was performed;the model was developed using Shapley integrated-empowerment benefit-distribution method.First,through literature survey and expert interview to identify the risk factors at various stages of the project,a dynamic risk-factor indicator system is developed.Second,to obtain a more meaningful risk-calculation result,the subjective and objective weights are combined,the weights of the risk factors at each stage are determined by the expert scoring method and entropy weight method,and the interest distribution model based on multi-dimensional risk factors is established.Finally,an example is used to verify the rationality of the method for the benefit distribution of the charging-pile project.The results of the example indicate that the limitations of the Shapley method can be reasonably avoided,and the applicability of the model for the benefit distribution of the charging-pile project is verified.展开更多
Extended finite element method （XFEM） implementation of the interaction integral methodology for evaluating the stress intensity factors （SIF） of the mixed-mode crack problem is presented. A discontinuous function...Extended finite element method （XFEM） implementation of the interaction integral methodology for evaluating the stress intensity factors （SIF） of the mixed-mode crack problem is presented. A discontinuous function and the near-tip asymptotic function are added to the classic finite element approximation to model the crack behavior. Two-state integral by the superposition of actual and auxiliary fields is derived to calculate the SIFs. Applications of the proposed technique to the inclined centre crack plate with inclined angle from 0° to 90° and slant edge crack plate with slant angle 45°, 67.5° and 90° are presented, and comparisons are made with closed form solutions. The results show that the proposed method is convenient, accurate and computationallv efficient.展开更多
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.展开更多
In order to improve the precision of mining subsidence prediction, a mathematical model using Support Vector Machine （SVM） was established to calculate the displacement factor. The study is based on a comprehensive ...In order to improve the precision of mining subsidence prediction, a mathematical model using Support Vector Machine （SVM） was established to calculate the displacement factor. The study is based on a comprehensive analysis of factors affecting the displacement factor, such as mechanical properties of the cover rock, the ratio of mining depth to seam thickness, dip angle of the coal seam and the thickness of loose layer. Data of 63 typical observation stations were used as a training and testing sample set. A SVM regression model of the displacement factor and the factors affecting it was established with a kernel function, an insensitive loss factor and a properly selected penalty factor. Given an accurate calculation algorithm for testing and analysis, the results show that an SVM regression model can calcu- late displacement factor precisely and reliable precision can be obtained which meets engineering requirements. The experimental results show that the method to calculation of the displacement factor, based on the SVM method, is feasible. The many factors affecting the displacement factor can be consid- ered with this method. The research provides an efficient and accurate approach for the calculation of displacement in mining subsidence orediction.展开更多
Using the single crack solution and the regular solution elf plane harmonic function, the problem of Saint-Venant bending of a cracked cylinder by a transverse force was reduced to solving two sets of integral equatio...Using the single crack solution and the regular solution elf plane harmonic function, the problem of Saint-Venant bending of a cracked cylinder by a transverse force was reduced to solving two sets of integral equations and its general solution was then obtained. Based on the obtained solution, a method to calculate the bending center and the stress intensity factors of the cracked cylinger whose cross-section is not thin-walled, but of small torsion rigidity is proposed. Some numerical examples are given.展开更多
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.展开更多
This paper presents a direct traction boundary integral equation method(DTBIEM)for two-dimensional crack problems of materials.The traction boundary integral equation was collocated on both the external boundary and e...This paper presents a direct traction boundary integral equation method(DTBIEM)for two-dimensional crack problems of materials.The traction boundary integral equation was collocated on both the external boundary and either side of the crack surfaces.The displacements and tractions were used as unknowns on the external boundary,while the relative crack opening displacement(RCOD)was chosen as unknowns on either side of crack surfaces to keep the single-domain merit.Only one side of the crack surfaces was concerned and needed to be discretized,thus the proposed method resulted in a smaller system of algebraic equations compared with the dual boundary element method(DBEM).A new set of crack-tip shape functions was constructed to represent the strain field singularity exactly,and the SIFs were evaluated by the extrapolation of the RCOD.Numerical examples for both straight and curved cracks are given to validate the accuracy and efficiency of the presented method.展开更多
A singular integral equation method is proposed to analyze the two-dimensional(2D)multiple cracks in anisotropic piezoelectric bimaterial.Using the Somigliana formula,a set of singular integral equations for the multi...A singular integral equation method is proposed to analyze the two-dimensional(2D)multiple cracks in anisotropic piezoelectric bimaterial.Using the Somigliana formula,a set of singular integral equations for the multiple crack problems are derived,in which the unknown functions are the derivatives of the generalized displacement discontinuities of the crack surfaces.Then,the exact analytical solution of the extended singular stresses and extended stress intensity factors near the crack tip is obtained.Singular integrals of the singular integral equations are computed by the Gauss-Chebyshev collocation method.Finally,numerical solutions of the extended stress intensity factors of some examples are presented and discussed.展开更多
[目的/意义]学术话语权是国家文化软实力的重要体现,研究学者学术话语权,丰富和完善学者学术话语权评价理论和指标体系,有助于国家学术话语权提升。[方法/过程]本文融合传统文献指标与Altmetrics指标,构建学者学术话语权综合评价体系。...[目的/意义]学术话语权是国家文化软实力的重要体现,研究学者学术话语权,丰富和完善学者学术话语权评价理论和指标体系,有助于国家学术话语权提升。[方法/过程]本文融合传统文献指标与Altmetrics指标,构建学者学术话语权综合评价体系。以Web of Science基因编辑领域数据为例,从学者学术话语影响力和社交话语引导力两维度出发,采用集成因子分析、熵权法、TOPSIS及二维评价法等方法进行实证分析。[结果/结论]研究表明,融合Altmetrics指标的综合评价体系具有一定的可信性,其中“篇均被引”和“Patent Mentions”是影响学者学术话语权的关键指标,且社交媒体和百度百科等传播平台能有效提高学者学术话语权。展开更多
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.展开更多
The behavior of two collinear anti-plane shear cracks in a piezoelectric layer bonded to two half spaces is investigated by the Schmidt method. The cracks are vertically to the interfaces of the piezoelectric layer. B...The behavior of two collinear anti-plane shear cracks in a piezoelectric layer bonded to two half spaces is investigated by the Schmidt method. The cracks are vertically to the interfaces of the piezoelectric layer. By using the Fourier transform, the problem can be solved with two pairs of triple integral equations. These equations are solved using the Schmidt method. This process is quite different from that adopted previously. Numerical examples are provided to show the effect of the geometry of the interacting cracks and the piezoelectric constants of the material upon the stress intensity factor of the cracks.展开更多
Stress intensity factors for a three dimensional rectangular interfacial crack were considered using the body force method. In the numerical calculations, unknown body force densities were approximated by the products...Stress intensity factors for a three dimensional rectangular interfacial crack were considered using the body force method. In the numerical calculations, unknown body force densities were approximated by the products of the fundamental densities and power series; here the fundamental densities are chosen to express singular stress fields due to an interface crack exactly. The calculation shows that the numerical results are satisfied. The stress intensity factors for a rectangular interface crack were indicated accurately with the varying aspect ratio, and bimaterial parameter.展开更多
For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geomet...For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geometries may lead to difficulties in the accuracy when discretizing the high-order derivatives on grid points near the boundary.It is very challenging to design numerical methods that can efficiently and accurately handle both difficulties.Applying an implicit scheme may be able to remove the stability constraints on the time step,however,it usually requires solving a large global system of nonlinear equations for each time step,and the computational cost could be significant.Integration factor(IF)or exponential time differencing(ETD)methods are one of the popular methods for temporal partial differential equations(PDEs)among many other methods.In our paper,we couple ETD methods with an embedded boundary method to solve a system of reaction-diffusion equations with complex geometries.In particular,we rewrite all ETD schemes into a linear combination of specificФ-functions and apply one state-of-the-art algorithm to compute the matrix-vector multiplications,which offers significant computational advantages with adaptive Krylov subspaces.In addition,we extend this method by incorporating the level set method to solve the free boundary problem.The accuracy,stability,and efficiency of the developed method are demonstrated by numerical examples.展开更多
文摘The compact implicit integration factor （cIIF） method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF method to some complex-valued nonlinear evolutionary equations such as the nonlinear SchrSdinger （NLS） equation and the complex Ginzburg-Landau （GL） equation. Detailed algorithm formulation and practical implementation of cIIF method are performed. The numerical results indicate that this method is very accurate and efficient.
基金the firancinal support of the National Natural Science Foundation of China(Grant No:51769011)for this work,and the authors are also deeply grateful to the editors and revewerse for tbeir rigorous work and valuable comments.
文摘The stress intensity factors(SIFs)for two-dimensional cracks are extracted using the p-version finite element method(P-FEM)and the contour integral method.Several numerical experiments,e.g.,crack initiating from the edge of a circular hole under an unidirectional uniform tension and two equal-length,unequal-length hole-edge cracks,respectively,at a rectangular plate,an inclined centered crack under uniaxial tension at a square plate and a pipeline crack model,are used to demonstrate the accuracy and effectiveness of the approaches.SIFs are presented for the effects of various crack lengths and length-width ratio.Numerical results are analyzed and compared with reference solutions and results obtained by the Voronoi cell finite element method,boundary element method,high-order extended finite element method(high-order XFEM)and commercial finite element software ABAQUS in the available literature.Numerical results are in good agreement with the benchmark problems and show faster convergence rate,higher accuracy and better numerical stability.
文摘Using the single crack solution and the regular solution of plane harmonic function, the problem of Saint_Venant bending of a cracked cylinder by a transverse force was reduced to solving two sets of integral equations and its general solution was then obtained. Based on the obtained solution, a method to calculate the bending center and the stress intensity factors of the cracked cylinger whose cross_section is not thin_walled, but of small torsion rigidity is proposed. Some numerical examples are given.
基金Supported by Science and Technology Foundation of SGCC Research and development of key models for decision support of energy internet companies(NO.SGSDJY00GPJS1900057).
文摘In this study,to develop a benefit-allocation model,in-depth analysis of a distributed photovoltaic-powergeneration carport and energy-storage charging-pile project was performed;the model was developed using Shapley integrated-empowerment benefit-distribution method.First,through literature survey and expert interview to identify the risk factors at various stages of the project,a dynamic risk-factor indicator system is developed.Second,to obtain a more meaningful risk-calculation result,the subjective and objective weights are combined,the weights of the risk factors at each stage are determined by the expert scoring method and entropy weight method,and the interest distribution model based on multi-dimensional risk factors is established.Finally,an example is used to verify the rationality of the method for the benefit distribution of the charging-pile project.The results of the example indicate that the limitations of the Shapley method can be reasonably avoided,and the applicability of the model for the benefit distribution of the charging-pile project is verified.
基金Projects(41172244,41072224) supported by the National Natural Science Foundation of ChinaProject(2009GGJS-037) supported by the Foundation of Youths Key Teacher by the Henan Educational Committee,China
文摘Extended finite element method （XFEM） implementation of the interaction integral methodology for evaluating the stress intensity factors （SIF） of the mixed-mode crack problem is presented. A discontinuous function and the near-tip asymptotic function are added to the classic finite element approximation to model the crack behavior. Two-state integral by the superposition of actual and auxiliary fields is derived to calculate the SIFs. Applications of the proposed technique to the inclined centre crack plate with inclined angle from 0° to 90° and slant edge crack plate with slant angle 45°, 67.5° and 90° are presented, and comparisons are made with closed form solutions. The results show that the proposed method is convenient, accurate and computationallv efficient.
基金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.
基金the Research and Innovation Program for College and University Graduate Students in Jiangsu Province (No.CX10B_141Z)the National Natural Science Foundation of China (No.41071273) for support of this project
文摘In order to improve the precision of mining subsidence prediction, a mathematical model using Support Vector Machine （SVM） was established to calculate the displacement factor. The study is based on a comprehensive analysis of factors affecting the displacement factor, such as mechanical properties of the cover rock, the ratio of mining depth to seam thickness, dip angle of the coal seam and the thickness of loose layer. Data of 63 typical observation stations were used as a training and testing sample set. A SVM regression model of the displacement factor and the factors affecting it was established with a kernel function, an insensitive loss factor and a properly selected penalty factor. Given an accurate calculation algorithm for testing and analysis, the results show that an SVM regression model can calcu- late displacement factor precisely and reliable precision can be obtained which meets engineering requirements. The experimental results show that the method to calculation of the displacement factor, based on the SVM method, is feasible. The many factors affecting the displacement factor can be consid- ered with this method. The research provides an efficient and accurate approach for the calculation of displacement in mining subsidence orediction.
文摘Using the single crack solution and the regular solution elf plane harmonic function, the problem of Saint-Venant bending of a cracked cylinder by a transverse force was reduced to solving two sets of integral equations and its general solution was then obtained. Based on the obtained solution, a method to calculate the bending center and the stress intensity factors of the cracked cylinger whose cross-section is not thin-walled, but of small torsion rigidity is proposed. Some numerical examples are given.
文摘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.
基金This work was supported by The National Key R&D Program of China(Grant No.2017YFC0804601)the National Natural Science Foundation of China(No.51741410)Open Research Fund of State Key Laboratory of Geomechanics and Geotechnical Engineering,Institute of Rock and Soil Mechanics,Chinese Academy of Sciences(Grant No.Z017017).
文摘This paper presents a direct traction boundary integral equation method(DTBIEM)for two-dimensional crack problems of materials.The traction boundary integral equation was collocated on both the external boundary and either side of the crack surfaces.The displacements and tractions were used as unknowns on the external boundary,while the relative crack opening displacement(RCOD)was chosen as unknowns on either side of crack surfaces to keep the single-domain merit.Only one side of the crack surfaces was concerned and needed to be discretized,thus the proposed method resulted in a smaller system of algebraic equations compared with the dual boundary element method(DBEM).A new set of crack-tip shape functions was constructed to represent the strain field singularity exactly,and the SIFs were evaluated by the extrapolation of the RCOD.Numerical examples for both straight and curved cracks are given to validate the accuracy and efficiency of the presented method.
基金The authors would like to express their special thanks to the National Natural Science Foundation of China(Project No.11172320).
文摘A singular integral equation method is proposed to analyze the two-dimensional(2D)multiple cracks in anisotropic piezoelectric bimaterial.Using the Somigliana formula,a set of singular integral equations for the multiple crack problems are derived,in which the unknown functions are the derivatives of the generalized displacement discontinuities of the crack surfaces.Then,the exact analytical solution of the extended singular stresses and extended stress intensity factors near the crack tip is obtained.Singular integrals of the singular integral equations are computed by the Gauss-Chebyshev collocation method.Finally,numerical solutions of the extended stress intensity factors of some examples are presented and discussed.
文摘[目的/意义]学术话语权是国家文化软实力的重要体现,研究学者学术话语权,丰富和完善学者学术话语权评价理论和指标体系,有助于国家学术话语权提升。[方法/过程]本文融合传统文献指标与Altmetrics指标,构建学者学术话语权综合评价体系。以Web of Science基因编辑领域数据为例,从学者学术话语影响力和社交话语引导力两维度出发,采用集成因子分析、熵权法、TOPSIS及二维评价法等方法进行实证分析。[结果/结论]研究表明,融合Altmetrics指标的综合评价体系具有一定的可信性,其中“篇均被引”和“Patent Mentions”是影响学者学术话语权的关键指标,且社交媒体和百度百科等传播平台能有效提高学者学术话语权。
基金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 Post Doctoral Science Foundation of Heilongjiang Province, the Natural Science Foundation of Heilongjiang Province，哈尔滨工业大学校科研和教改项目
文摘The behavior of two collinear anti-plane shear cracks in a piezoelectric layer bonded to two half spaces is investigated by the Schmidt method. The cracks are vertically to the interfaces of the piezoelectric layer. By using the Fourier transform, the problem can be solved with two pairs of triple integral equations. These equations are solved using the Schmidt method. This process is quite different from that adopted previously. Numerical examples are provided to show the effect of the geometry of the interacting cracks and the piezoelectric constants of the material upon the stress intensity factor of the cracks.
文摘Stress intensity factors for a three dimensional rectangular interfacial crack were considered using the body force method. In the numerical calculations, unknown body force densities were approximated by the products of the fundamental densities and power series; here the fundamental densities are chosen to express singular stress fields due to an interface crack exactly. The calculation shows that the numerical results are satisfied. The stress intensity factors for a rectangular interface crack were indicated accurately with the varying aspect ratio, and bimaterial parameter.
文摘For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geometries may lead to difficulties in the accuracy when discretizing the high-order derivatives on grid points near the boundary.It is very challenging to design numerical methods that can efficiently and accurately handle both difficulties.Applying an implicit scheme may be able to remove the stability constraints on the time step,however,it usually requires solving a large global system of nonlinear equations for each time step,and the computational cost could be significant.Integration factor(IF)or exponential time differencing(ETD)methods are one of the popular methods for temporal partial differential equations(PDEs)among many other methods.In our paper,we couple ETD methods with an embedded boundary method to solve a system of reaction-diffusion equations with complex geometries.In particular,we rewrite all ETD schemes into a linear combination of specificФ-functions and apply one state-of-the-art algorithm to compute the matrix-vector multiplications,which offers significant computational advantages with adaptive Krylov subspaces.In addition,we extend this method by incorporating the level set method to solve the free boundary problem.The accuracy,stability,and efficiency of the developed method are demonstrated by numerical examples.
