Classical Chinese characters,presented through calligraphy,seal engraving,or painting,can exhibit different aesthetics and essences of Chinese characters,making them the most important asset of the Chinese people.Call...Classical Chinese characters,presented through calligraphy,seal engraving,or painting,can exhibit different aesthetics and essences of Chinese characters,making them the most important asset of the Chinese people.Calligraphy and seal engraving,as two closely related systems in traditional Chinese art,have developed through the ages.Due to changes in lifestyle and advancements in modern technology,their original functions of daily writing and verification have gradually diminished.Instead,they have increasingly played a significant role in commercial art.This study utilizes the Evaluation Grid Method(EGM)and the Analytic Hierarchy Process(AHP)to research the key preference factors in the application of calligraphy and seal engraving imagery.Different from the traditional 5-point equal interval semantic questionnaire,this study employs a non-equal interval semantic questionnaire with a golden ratio scale,distinguishing the importance ratio of adjacent semantic meanings and highlighting the weighted emphasis on visual aesthetics.Additionally,the study uses Importance-Performance Analysis(IPA)and Technique for Order Preference by Similarity to Ideal Solution(TOPSIS)to obtain the key preference sequence of calligraphy and seal engraving culture.Plus,the Choquet integral comprehensive evaluation is used as a reference for IPA comparison.It is hoped that this study can provide cultural imagery references and research methods,injecting further creativity into industrial design.展开更多
Due to the difficulties in obtaining large deformation mining subsidence using differential Interferometric Synthetic Aperture Radar (D-InSAR) alone, a new algorithm was proposed to extract large deformation mining ...Due to the difficulties in obtaining large deformation mining subsidence using differential Interferometric Synthetic Aperture Radar (D-InSAR) alone, a new algorithm was proposed to extract large deformation mining subsidence using D-InSAR technique and probability integral method. The details of the algorithm are as follows:the control points set, containing correct phase unwrapping points on the subsidence basin edge generated by D-InSAR and several observation points (near the maximum subsidence and inflection points), was established at first; genetic algorithm (GA) was then used to optimize the parameters of probability integral method; at last, the surface subsidence was deduced according to the optimum parameters. The results of the experiment in Huaibei mining area, China, show that the presented method can generate the correct mining subsidence basin with a few surface observations, and the relative error of maximum subsidence point is about 8.3%, which is much better than that of conventional D-InSAR (relative error is 68.0%).展开更多
In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a gene...In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.展开更多
The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the bas...The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the basic function and of the weight function,and is mainly determined by that of the weight function.Therefore,the weight function greatly affects the accuracy of results obtained.Different kinds of weight functions,such as the spline function, the Gauss function and so on,are proposed recently by many researchers.In the present work,the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method.The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed.Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and α in Gauss and exponential weight functions are in the range of reasonable values,respectively,and the higher the smoothness of the weight function,the better the features of the solutions.展开更多
On the basis of three geological models and several orebody boundaries, a method of grid subdivision and integral has been proposed to calculate and evaluate the resources of cobalt-rich crusts on the seamounts in the...On the basis of three geological models and several orebody boundaries, a method of grid subdivision and integral has been proposed to calculate and evaluate the resources of cobalt-rich crusts on the seamounts in the central Pacific Ocean. The formulas of this method are deduced and the interface of program module is designed. The method is carried out in the software "Auto mapping system of submarine topography and geomorphology MBChart". This method and program will possibly become a potential tool to calculate the resources of seamounts and determine the target diggings for China' s next Five-year Plan.展开更多
This paper determines the exact error order on optimization of adaptive direct methods of approximate solution of the class of Fredholm integral equations of the second kind with kernel belonging to the anisotropic So...This paper determines the exact error order on optimization of adaptive direct methods of approximate solution of the class of Fredholm integral equations of the second kind with kernel belonging to the anisotropic Sobolev classes, and also gives an optimal algorithm.展开更多
For deep tunnel projects,selecting an appropriate initial support distance is critical to improving the self-supporting capacity of surrounding rock.In this work,an intuitive method for determining the tunnel’s initi...For deep tunnel projects,selecting an appropriate initial support distance is critical to improving the self-supporting capacity of surrounding rock.In this work,an intuitive method for determining the tunnel’s initial support distance was proposed.First,based on the convergence-confinement method,a three-dimensional analytical model was constructed by combining an analytical solution of a non-circular tunnel with the Tecplot software.Then,according to the integral failure criteria of rock,the failure tendency coefficients of hard surrounding rock were computed and the spatial distribution plots of that were constructed.On this basis,the tunnel’s key failure positions were identified,and the relationship between the failure tendency coefficient at key failure positions and their distances from the working face was established.Finally,the distance from the working face that corresponds to the critical failure tendency coefficient was taken as the optimal support distance.A practical project was used as an example,and a reasonable initial support distance was successfully determined by applying the developed method.Moreover,it is found that the stability of hard surrounding rock decreases rapidly within the range of 1.0D(D is the tunnel diameter)from the working face,and tends to be stable outside the range of 1.0D.展开更多
For the longitudinal seismic response analysis of a tunnel structure under asynchronous earthquake excitations,a longitudinal integral response deformation method classified as a practical approach is proposed in this...For the longitudinal seismic response analysis of a tunnel structure under asynchronous earthquake excitations,a longitudinal integral response deformation method classified as a practical approach is proposed in this paper.The determinations of the structural critical moments when maximal deformations and internal forces in the longitudinal direction occur are deduced as well.When applying the proposed method,the static analysis of the free-field computation model subjected to the least favorable free-field deformation at the tunnel buried depth is performed first to calculate the equivalent input seismic loads.Then,the equivalent input seismic loads are imposed on the integral tunnel-foundation computation model to conduct the static calculation.Afterwards,the critical longitudinal seismic responses of the tunnel are obtained.The applicability of the new method is verified by comparing the seismic responses of a shield tunnel structure in Beijing,determined by the proposed procedure and by a dynamic time-history analysis under a series of obliquely incident out-of-plane and in-plane waves.The results show that the proposed method has a clear concept with high accuracy and simple progress.Meanwhile,this method provides a feasible way to determine the critical moments of the longitudinal seismic responses of a tunnel structure.Therefore,the proposed method can be effectively applied to analyze the seismic response of a long-line underground structure subjected to non-uniform excitations.展开更多
In this paper, we suggest a method for solving Fredholm integral equation of the first kind based on wavelet basis. The continuous Legendre and Chebyshev wavelets of the first, second, third and fourth kind on [0,1] a...In this paper, we suggest a method for solving Fredholm integral equation of the first kind based on wavelet basis. The continuous Legendre and Chebyshev wavelets of the first, second, third and fourth kind on [0,1] are used and are utilized as a basis in Galerkin method to approximate the solution of integral equations. Then, in some examples the mentioned wavelets are compared with each other.展开更多
A direct method to find the first integral for two-dimensional autonomous system in polar coordinates is suggested. It is shown that if the equation of motion expressed by differential 1-forms for a given autonomous H...A direct method to find the first integral for two-dimensional autonomous system in polar coordinates is suggested. It is shown that if the equation of motion expressed by differential 1-forms for a given autonomous Hamiltonian system is multiplied by a set of multiplicative functions, then the general expression of the first integral can be obtained, An example is given to illustrate the application of the results.展开更多
Using the Picard iteration method and treating the involved integration by numerical quadrature formulas, we propose a numerical scheme for the second kind nonlinear Volterra integral equations. For enlarging the conv...Using the Picard iteration method and treating the involved integration by numerical quadrature formulas, we propose a numerical scheme for the second kind nonlinear Volterra integral equations. For enlarging the convergence region of the Picard iteration method, multistage algorithm is devised. We also introduce an algorithm for problems with some singularities at the limits of integration including fractional integral equations. Numerical tests verify the validity of the proposed schemes.展开更多
In this article, a general formula of the first integral method has been extended to celebrate the exact solution of nonlinear time-space differential equations of fractional orders. The proposed method is easy, direc...In this article, a general formula of the first integral method has been extended to celebrate the exact solution of nonlinear time-space differential equations of fractional orders. The proposed method is easy, direct and concise as compared with other existent methods.展开更多
The Laguerre polynomial method has been successfully used to investigate the dynamic responses of a half-space.However,it fails to obtain the correct stress at the interfaces in a layered half-space,especially when th...The Laguerre polynomial method has been successfully used to investigate the dynamic responses of a half-space.However,it fails to obtain the correct stress at the interfaces in a layered half-space,especially when there are significant differences in material properties.Therefore,a coupled Legendre-Laguerre polynomial method with analytical integration is proposed.The Rayleigh waves in a one-dimensional(1D)hexagonal quasicrystal(QC)layered half-space with an imperfect interface are investigated.The correctness is validated by comparison with available results.Its computation efficiency is analyzed.The dispersion curves of the phase velocity,displacement distributions,and stress distributions are illustrated.The effects of the phonon-phason coupling and imperfect interface coefficients on the wave characteristics are investigated.Some novel findings reveal that the proposed method is highly efficient for addressing the Rayleigh waves in a QC layered half-space.It can save over 99%of the computation time.This method can be expanded to investigate waves in various layered half-spaces,including earth-layered media and surface acoustic wave(SAW)devices.展开更多
The meshless local boundary integral equation method is a currently developed numerical method, which combines the advantageous features of Galerkin finite element method(GFEM), boundary element method(BEM) and elemen...The meshless local boundary integral equation method is a currently developed numerical method, which combines the advantageous features of Galerkin finite element method(GFEM), boundary element method(BEM) and element free Galerkin method(EFGM), and is a truly meshless method possessing wide prospects in engineering applications. The companion solution and all the other formulas required in the meshless local boundary integral equation for a thin plate were presented, in order to make this method apply to solve the thin plate problem.展开更多
The paper applied the isogeometric boundary element method(IGABEM)to thermoelastic problems.The Non-Uniform Rational B-splines(NURBS)used to construct geometric models are employed to discretize the boundary integral ...The paper applied the isogeometric boundary element method(IGABEM)to thermoelastic problems.The Non-Uniform Rational B-splines(NURBS)used to construct geometric models are employed to discretize the boundary integral formulation of the governing equation.Due to the existence of thermal stress,the domain integral term appears in the boundary integral equation.We resolve this problem by incorporating radial integration method into IGABEM which converts the domain integral to the boundary integral.In this way,IGABEM can maintain its advantages in dimensionality reduction and more importantly,seamless integration of CAD and numerical analysis based on boundary representation.The algorithm is verified by numerical examples.展开更多
This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solve...This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.展开更多
In this paper, we present a new method, a mixture of homotopy perturbation method and a new integral transform to solve some nonlinear partial differential equations. The proposed method introduces also He’s polynomi...In this paper, we present a new method, a mixture of homotopy perturbation method and a new integral transform to solve some nonlinear partial differential equations. The proposed method introduces also He’s polynomials [1]. The analytical results of examples are calculated in terms of convergent series with easily computed components [2].展开更多
In this paper,the classical composite middle rectangle rule for the computation of Cauchy principal value integral(the singular kernel 1=(x-s))is discussed.With the density function approximated only while the singula...In this paper,the classical composite middle rectangle rule for the computation of Cauchy principal value integral(the singular kernel 1=(x-s))is discussed.With the density function approximated only while the singular kernel is calculated analysis,then the error functional of asymptotic expansion is obtained.We construct a series to approach the singular point.An extrapolation algorithm is presented and the convergence rate of extrapolation algorithm is proved.At last,some numerical results are presented to confirm the theoretical results and show the efficiency of the algorithms.展开更多
In recent years, many methods have been used to find the exact solutions of nonlinear partial differential equations. One of them is called the first integral method, which is based on the ring theory of commutative a...In recent years, many methods have been used to find the exact solutions of nonlinear partial differential equations. One of them is called the first integral method, which is based on the ring theory of commutative algebra. In this paper, exact travelling wave solutions of the Non-Boussinesq wavepacket model and the (2 + 1)-dimensional Zoomeron equation are studied by using the first integral method. From the solving process and results, the first integral method has the characteristics of simplicity, directness and effectiveness about solving the exact travelling wave solutions of nonlinear partial differential equations. In other words, tedious calculations can be avoided by Maple software;the solutions of more accurate and richer travelling wave solutions are obtained. Therefore, this method is an effective method for solving exact solutions of nonlinear partial differential equations.展开更多
The main aim of the paper is to examine the concentration of the longitudinal dispersion phenomenon arising in fluid flow through porous media. These phenomenon yields a partial differential equation namely Burger’s ...The main aim of the paper is to examine the concentration of the longitudinal dispersion phenomenon arising in fluid flow through porous media. These phenomenon yields a partial differential equation namely Burger’s equation, which is solved by mixture of the new integral transform and the homotopy perturbation method under suitable conditions and the standard assumption. This method provides an analytical approximation in a rapidly convergent sequence with in exclusive manner computed terms. Its rapid convergence shows that the method is trustworthy and introduces a significant improvement in solving nonlinear partial differential equations over existing methods. It is concluded that the behaviour of concentration in longitudinal dispersion phenomenon is decreases as distance x is increasing with fixed time t > 0 and slightly increases with time t.展开更多
文摘Classical Chinese characters,presented through calligraphy,seal engraving,or painting,can exhibit different aesthetics and essences of Chinese characters,making them the most important asset of the Chinese people.Calligraphy and seal engraving,as two closely related systems in traditional Chinese art,have developed through the ages.Due to changes in lifestyle and advancements in modern technology,their original functions of daily writing and verification have gradually diminished.Instead,they have increasingly played a significant role in commercial art.This study utilizes the Evaluation Grid Method(EGM)and the Analytic Hierarchy Process(AHP)to research the key preference factors in the application of calligraphy and seal engraving imagery.Different from the traditional 5-point equal interval semantic questionnaire,this study employs a non-equal interval semantic questionnaire with a golden ratio scale,distinguishing the importance ratio of adjacent semantic meanings and highlighting the weighted emphasis on visual aesthetics.Additionally,the study uses Importance-Performance Analysis(IPA)and Technique for Order Preference by Similarity to Ideal Solution(TOPSIS)to obtain the key preference sequence of calligraphy and seal engraving culture.Plus,the Choquet integral comprehensive evaluation is used as a reference for IPA comparison.It is hoped that this study can provide cultural imagery references and research methods,injecting further creativity into industrial design.
基金Project (BK20130174) supported by the Basic Research Project of Jiangsu Province (Natural Science Foundation) Project (1101109C) supported by Jiangsu Planned Projects for Postdoctoral Research Funds,China+1 种基金Project (201325) supported by the Key Laboratory of Geo-informatics of State Bureau of Surveying and Mapping,ChinaProject (SZBF2011-6-B35) supported by the Priority Academic Program Development of Jiangsu Higher Education Institutions,China
文摘Due to the difficulties in obtaining large deformation mining subsidence using differential Interferometric Synthetic Aperture Radar (D-InSAR) alone, a new algorithm was proposed to extract large deformation mining subsidence using D-InSAR technique and probability integral method. The details of the algorithm are as follows:the control points set, containing correct phase unwrapping points on the subsidence basin edge generated by D-InSAR and several observation points (near the maximum subsidence and inflection points), was established at first; genetic algorithm (GA) was then used to optimize the parameters of probability integral method; at last, the surface subsidence was deduced according to the optimum parameters. The results of the experiment in Huaibei mining area, China, show that the presented method can generate the correct mining subsidence basin with a few surface observations, and the relative error of maximum subsidence point is about 8.3%, which is much better than that of conventional D-InSAR (relative error is 68.0%).
基金supported by the Swiss National Science Foundation(Grant No.189882)the National Natural Science Foundation of China(Grant No.41961134032)support provided by the New Investigator Award grant from the UK Engineering and Physical Sciences Research Council(Grant No.EP/V012169/1).
文摘In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.
文摘The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the basic function and of the weight function,and is mainly determined by that of the weight function.Therefore,the weight function greatly affects the accuracy of results obtained.Different kinds of weight functions,such as the spline function, the Gauss function and so on,are proposed recently by many researchers.In the present work,the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method.The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed.Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and α in Gauss and exponential weight functions are in the range of reasonable values,respectively,and the higher the smoothness of the weight function,the better the features of the solutions.
基金This study was supported by Projects under contract Nos DY105 China's 0cean-03-01-01 and DY105-03-01-07the National Natural Science Foundation of China under contract No.40506017the Youth Foundation of Marine High-tech Project of China under contract No.2002AA616010.
文摘On the basis of three geological models and several orebody boundaries, a method of grid subdivision and integral has been proposed to calculate and evaluate the resources of cobalt-rich crusts on the seamounts in the central Pacific Ocean. The formulas of this method are deduced and the interface of program module is designed. The method is carried out in the software "Auto mapping system of submarine topography and geomorphology MBChart". This method and program will possibly become a potential tool to calculate the resources of seamounts and determine the target diggings for China' s next Five-year Plan.
基金Project supported by the Natural Science Foundation of China(10371009)Research Fund for the Doctoral Program Higher Education
文摘This paper determines the exact error order on optimization of adaptive direct methods of approximate solution of the class of Fredholm integral equations of the second kind with kernel belonging to the anisotropic Sobolev classes, and also gives an optimal algorithm.
基金Project(2021JLM-49) supported by Natural Science Basic Research Program of Shaanxi-Joint Fund of Hanjiang to Weihe River Valley Water Diversion Project,ChinaProject(42077248) supported by the National Natural Science Foundation of China
文摘For deep tunnel projects,selecting an appropriate initial support distance is critical to improving the self-supporting capacity of surrounding rock.In this work,an intuitive method for determining the tunnel’s initial support distance was proposed.First,based on the convergence-confinement method,a three-dimensional analytical model was constructed by combining an analytical solution of a non-circular tunnel with the Tecplot software.Then,according to the integral failure criteria of rock,the failure tendency coefficients of hard surrounding rock were computed and the spatial distribution plots of that were constructed.On this basis,the tunnel’s key failure positions were identified,and the relationship between the failure tendency coefficient at key failure positions and their distances from the working face was established.Finally,the distance from the working face that corresponds to the critical failure tendency coefficient was taken as the optimal support distance.A practical project was used as an example,and a reasonable initial support distance was successfully determined by applying the developed method.Moreover,it is found that the stability of hard surrounding rock decreases rapidly within the range of 1.0D(D is the tunnel diameter)from the working face,and tends to be stable outside the range of 1.0D.
基金National Natural Science Foundation of China under Grant No.51478247。
文摘For the longitudinal seismic response analysis of a tunnel structure under asynchronous earthquake excitations,a longitudinal integral response deformation method classified as a practical approach is proposed in this paper.The determinations of the structural critical moments when maximal deformations and internal forces in the longitudinal direction occur are deduced as well.When applying the proposed method,the static analysis of the free-field computation model subjected to the least favorable free-field deformation at the tunnel buried depth is performed first to calculate the equivalent input seismic loads.Then,the equivalent input seismic loads are imposed on the integral tunnel-foundation computation model to conduct the static calculation.Afterwards,the critical longitudinal seismic responses of the tunnel are obtained.The applicability of the new method is verified by comparing the seismic responses of a shield tunnel structure in Beijing,determined by the proposed procedure and by a dynamic time-history analysis under a series of obliquely incident out-of-plane and in-plane waves.The results show that the proposed method has a clear concept with high accuracy and simple progress.Meanwhile,this method provides a feasible way to determine the critical moments of the longitudinal seismic responses of a tunnel structure.Therefore,the proposed method can be effectively applied to analyze the seismic response of a long-line underground structure subjected to non-uniform excitations.
文摘In this paper, we suggest a method for solving Fredholm integral equation of the first kind based on wavelet basis. The continuous Legendre and Chebyshev wavelets of the first, second, third and fourth kind on [0,1] are used and are utilized as a basis in Galerkin method to approximate the solution of integral equations. Then, in some examples the mentioned wavelets are compared with each other.
文摘A direct method to find the first integral for two-dimensional autonomous system in polar coordinates is suggested. It is shown that if the equation of motion expressed by differential 1-forms for a given autonomous Hamiltonian system is multiplied by a set of multiplicative functions, then the general expression of the first integral can be obtained, An example is given to illustrate the application of the results.
文摘Using the Picard iteration method and treating the involved integration by numerical quadrature formulas, we propose a numerical scheme for the second kind nonlinear Volterra integral equations. For enlarging the convergence region of the Picard iteration method, multistage algorithm is devised. We also introduce an algorithm for problems with some singularities at the limits of integration including fractional integral equations. Numerical tests verify the validity of the proposed schemes.
文摘In this article, a general formula of the first integral method has been extended to celebrate the exact solution of nonlinear time-space differential equations of fractional orders. The proposed method is easy, direct and concise as compared with other existent methods.
基金Project supported by the National Natural Science Foundation of China(No.12102131)the Natural Science Foundation of Henan Province of China(No.242300420248)the International Science and Technology Cooperation Project of Henan Province of China(No.242102521010)。
文摘The Laguerre polynomial method has been successfully used to investigate the dynamic responses of a half-space.However,it fails to obtain the correct stress at the interfaces in a layered half-space,especially when there are significant differences in material properties.Therefore,a coupled Legendre-Laguerre polynomial method with analytical integration is proposed.The Rayleigh waves in a one-dimensional(1D)hexagonal quasicrystal(QC)layered half-space with an imperfect interface are investigated.The correctness is validated by comparison with available results.Its computation efficiency is analyzed.The dispersion curves of the phase velocity,displacement distributions,and stress distributions are illustrated.The effects of the phonon-phason coupling and imperfect interface coefficients on the wave characteristics are investigated.Some novel findings reveal that the proposed method is highly efficient for addressing the Rayleigh waves in a QC layered half-space.It can save over 99%of the computation time.This method can be expanded to investigate waves in various layered half-spaces,including earth-layered media and surface acoustic wave(SAW)devices.
文摘The meshless local boundary integral equation method is a currently developed numerical method, which combines the advantageous features of Galerkin finite element method(GFEM), boundary element method(BEM) and element free Galerkin method(EFGM), and is a truly meshless method possessing wide prospects in engineering applications. The companion solution and all the other formulas required in the meshless local boundary integral equation for a thin plate were presented, in order to make this method apply to solve the thin plate problem.
基金This study was funded by the National Natural Science Foundation of China(NSFC)(Grant Nos.11702238,51904202 and 11902212)and Nanhu Scholars Program for Young Scholars of XYNU.
文摘The paper applied the isogeometric boundary element method(IGABEM)to thermoelastic problems.The Non-Uniform Rational B-splines(NURBS)used to construct geometric models are employed to discretize the boundary integral formulation of the governing equation.Due to the existence of thermal stress,the domain integral term appears in the boundary integral equation.We resolve this problem by incorporating radial integration method into IGABEM which converts the domain integral to the boundary integral.In this way,IGABEM can maintain its advantages in dimensionality reduction and more importantly,seamless integration of CAD and numerical analysis based on boundary representation.The algorithm is verified by numerical examples.
基金the National Science and Tech-nology Council,Taiwan for their financial support(Grant Number NSTC 111-2221-E-019-048).
文摘This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.
文摘In this paper, we present a new method, a mixture of homotopy perturbation method and a new integral transform to solve some nonlinear partial differential equations. The proposed method introduces also He’s polynomials [1]. The analytical results of examples are calculated in terms of convergent series with easily computed components [2].
基金The work of Jin Li was supported by National Natural Science Foundation of China(Grant No.11471195)China Postdoctoral Science Foundation(Grant No.2015T80703)+1 种基金Shan-dong Provincial Natural Science Foundation of China(Grant No.ZR2016JL006)Na-tional Natural Science Foundation of China(Grant No.11771398).
文摘In this paper,the classical composite middle rectangle rule for the computation of Cauchy principal value integral(the singular kernel 1=(x-s))is discussed.With the density function approximated only while the singular kernel is calculated analysis,then the error functional of asymptotic expansion is obtained.We construct a series to approach the singular point.An extrapolation algorithm is presented and the convergence rate of extrapolation algorithm is proved.At last,some numerical results are presented to confirm the theoretical results and show the efficiency of the algorithms.
文摘In recent years, many methods have been used to find the exact solutions of nonlinear partial differential equations. One of them is called the first integral method, which is based on the ring theory of commutative algebra. In this paper, exact travelling wave solutions of the Non-Boussinesq wavepacket model and the (2 + 1)-dimensional Zoomeron equation are studied by using the first integral method. From the solving process and results, the first integral method has the characteristics of simplicity, directness and effectiveness about solving the exact travelling wave solutions of nonlinear partial differential equations. In other words, tedious calculations can be avoided by Maple software;the solutions of more accurate and richer travelling wave solutions are obtained. Therefore, this method is an effective method for solving exact solutions of nonlinear partial differential equations.
文摘The main aim of the paper is to examine the concentration of the longitudinal dispersion phenomenon arising in fluid flow through porous media. These phenomenon yields a partial differential equation namely Burger’s equation, which is solved by mixture of the new integral transform and the homotopy perturbation method under suitable conditions and the standard assumption. This method provides an analytical approximation in a rapidly convergent sequence with in exclusive manner computed terms. Its rapid convergence shows that the method is trustworthy and introduces a significant improvement in solving nonlinear partial differential equations over existing methods. It is concluded that the behaviour of concentration in longitudinal dispersion phenomenon is decreases as distance x is increasing with fixed time t > 0 and slightly increases with time t.