With the linear interpolation method, an improved absorbing boundary condition（ABC）is introduced and derived, which is suitable for the alternating-direction-implicit finite- difference time-domain （ADI-FDTD） meth...With the linear interpolation method, an improved absorbing boundary condition（ABC）is introduced and derived, which is suitable for the alternating-direction-implicit finite- difference time-domain （ADI-FDTD） method. The reflection of the ABC caused by both the truncated error and the phase velocity error is analyzed. Based on the phase velocity estimation and the nonuniform cell, two methods are studied and then adopted to improve the performance of the ABC. A calculation case of a rectangular waveguide which is a typical dispersive transmission line is carried out using the ADI-FDTD method with the improved ABC for evaluation. According to the calculated case, the comparison is given between the reflection coefficients of the ABC with and without the velocity estimation and also the comparison between the reflection coefficients of the ABC with and without the nonuniform processing. The reflection variation of the ABC under different time steps is also analyzed and the acceptable worsening will not obscure the improvement on the absorption. Numerical results obviously show that efficient improvement on the absorbing performance of the ABC is achieved based on these methods for the ADI-FDTD.展开更多
Euler-Bernoulli beam equation is very important that can be applied in the field of mechanics, science and technology. Some authors have put forward many different numerical methods, but the precision is not enough hi...Euler-Bernoulli beam equation is very important that can be applied in the field of mechanics, science and technology. Some authors have put forward many different numerical methods, but the precision is not enough high. In this paper, we will illustrate the high-precision numerical method to solve Euler-Bernoulli beam equation. Three numerical examples are studied to demonstrate the accuracy of the present method. Results obtained by our method indicate new algorithm has the following advantages: small computational work, fast convergence speed and high precision.展开更多
A new seismic ray-tracing method is put forward based on parabolic travel-time interpolation(PTI) method, which is more accurate than the linear travel-time interpolation (LTI) method. Both PTI method and LTI method a...A new seismic ray-tracing method is put forward based on parabolic travel-time interpolation(PTI) method, which is more accurate than the linear travel-time interpolation (LTI) method. Both PTI method and LTI method are used to compute seismic travel-time and ray-path in a 2-D grid cell model. Firstly, some basic concepts are introduced. The calculations of travel-time and ray-path are carried out only at cell boundaries. So, the ray-path is always straight in the same cells with uniform velocity. Two steps are applied in PTI and LTI method, step 1 computes travel-time and step 2 traces ray-path. Then, the derivation of LTI formulas is described. Because of the presence of refraction wave in shot cell, the formula aiming at shot cell is also derived. Finally, PTI method is presented. The calculation of PTI method is more complex than that of LTI method, but the error is limited. The results of numerical model show that PTI method can trace ray-path more accurately and efficiently than LTI method does.展开更多
C^1 natural element method （C^1 NEM） is applied to strain gradient linear elasticity, and size effects on mi crostructures are analyzed. The shape functions in C^1 NEM are built upon the natural neighbor interpolati...C^1 natural element method （C^1 NEM） is applied to strain gradient linear elasticity, and size effects on mi crostructures are analyzed. The shape functions in C^1 NEM are built upon the natural neighbor interpolation （NNI）, with interpolation realized to nodal function and nodal gradient values, so that the essential boundary conditions （EBCs） can be imposed directly in a Galerkin scheme for partial differential equations （PDEs）. In the present paper, C^1 NEM for strain gradient linear elasticity is constructed, and sev- eral typical examples which have analytical solutions are presented to illustrate the effectiveness of the constructed method. In its application to microstructures, the size effects of bending stiffness and stress concentration factor （SCF） are studied for microspeciem and microgripper, respectively. It is observed that the size effects become rather strong when the width of spring for microgripper, the radius of circular perforation and the long axis of elliptical perforation for microspeciem come close to the material characteristic length scales. For the U-shaped notch, the size effects decline obviously with increasing notch radius, and decline mildly with increasing length of notch.展开更多
Multistep integration methods are being extensively used in the simulations of high dimensional systems due to their lower computational cost.The block methods were developed with the intent of obtaining numerical res...Multistep integration methods are being extensively used in the simulations of high dimensional systems due to their lower computational cost.The block methods were developed with the intent of obtaining numerical results on numerous points at a time and improving computational efficiency.Hybrid block methods for instance are specifically used in numerical integration of initial value problems.In this paper,an optimized hybrid block Adams block method is designed for the solutions of linear and nonlinear first-order initial value problems in ordinary differential equations(ODEs).In deriving themethod,the Lagrange interpolation polynomial was employed based on some data points to replace the differential equation function and it was integrated over a specified interval.Furthermore,the convergence properties along with the region of stability of the method were examined.It was concluded that the newly derived method is convergent,consistent,and zero-stable.The method was also found to be A-stable implying that it covers the whole of the left/negative half plane.From the numerical computations of absolute errors carried out using the newly derived method,it was found that the method performed better than the ones with which we compared our results with.Themethod also showed its superiority over the existing methods in terms of stability and convergence.展开更多
By combining the classical appropriate functions “1, x, x 2” with the method of multiplier enlargement, this paper establishes a theorem to approximate any unbounded continuous functions with modified positive...By combining the classical appropriate functions “1, x, x 2” with the method of multiplier enlargement, this paper establishes a theorem to approximate any unbounded continuous functions with modified positive linear operators. As an example, Hermite Fejér interpolation polynomial operators are analysed and studied, and a general conclusion is obtained.展开更多
In this paper,an optimality condition for nonlinear programming problems with box constraints is given by using linear transformation and Lagrange interpolating polynomials.Based on this condition,two new local optim...In this paper,an optimality condition for nonlinear programming problems with box constraints is given by using linear transformation and Lagrange interpolating polynomials.Based on this condition,two new local optimization methods are developed.The solution points obtained by the new local optimization methods can improve the Karush–Kuhn–Tucker(KKT)points in general.Two global optimization methods then are proposed by combining the two new local optimization methods with a filled function method.Some numerical examples are reported to show the effectiveness of the proposed methods.展开更多
基金The National Natural Science Foundation of China(No.60702027)the Free Research Fund of the National Mobile Communications Research Laboratory of Southeast University (No.2008B07)the National Basic Research Program of China(973 Program)(No.2007CB310603)
文摘With the linear interpolation method, an improved absorbing boundary condition（ABC）is introduced and derived, which is suitable for the alternating-direction-implicit finite- difference time-domain （ADI-FDTD） method. The reflection of the ABC caused by both the truncated error and the phase velocity error is analyzed. Based on the phase velocity estimation and the nonuniform cell, two methods are studied and then adopted to improve the performance of the ABC. A calculation case of a rectangular waveguide which is a typical dispersive transmission line is carried out using the ADI-FDTD method with the improved ABC for evaluation. According to the calculated case, the comparison is given between the reflection coefficients of the ABC with and without the velocity estimation and also the comparison between the reflection coefficients of the ABC with and without the nonuniform processing. The reflection variation of the ABC under different time steps is also analyzed and the acceptable worsening will not obscure the improvement on the absorption. Numerical results obviously show that efficient improvement on the absorbing performance of the ABC is achieved based on these methods for the ADI-FDTD.
文摘Euler-Bernoulli beam equation is very important that can be applied in the field of mechanics, science and technology. Some authors have put forward many different numerical methods, but the precision is not enough high. In this paper, we will illustrate the high-precision numerical method to solve Euler-Bernoulli beam equation. Three numerical examples are studied to demonstrate the accuracy of the present method. Results obtained by our method indicate new algorithm has the following advantages: small computational work, fast convergence speed and high precision.
文摘A new seismic ray-tracing method is put forward based on parabolic travel-time interpolation(PTI) method, which is more accurate than the linear travel-time interpolation (LTI) method. Both PTI method and LTI method are used to compute seismic travel-time and ray-path in a 2-D grid cell model. Firstly, some basic concepts are introduced. The calculations of travel-time and ray-path are carried out only at cell boundaries. So, the ray-path is always straight in the same cells with uniform velocity. Two steps are applied in PTI and LTI method, step 1 computes travel-time and step 2 traces ray-path. Then, the derivation of LTI formulas is described. Because of the presence of refraction wave in shot cell, the formula aiming at shot cell is also derived. Finally, PTI method is presented. The calculation of PTI method is more complex than that of LTI method, but the error is limited. The results of numerical model show that PTI method can trace ray-path more accurately and efficiently than LTI method does.
基金supported by the SDUST Spring Bud (2009AZZ021)Taian Science and Technology Development (20112001)
文摘C^1 natural element method （C^1 NEM） is applied to strain gradient linear elasticity, and size effects on mi crostructures are analyzed. The shape functions in C^1 NEM are built upon the natural neighbor interpolation （NNI）, with interpolation realized to nodal function and nodal gradient values, so that the essential boundary conditions （EBCs） can be imposed directly in a Galerkin scheme for partial differential equations （PDEs）. In the present paper, C^1 NEM for strain gradient linear elasticity is constructed, and sev- eral typical examples which have analytical solutions are presented to illustrate the effectiveness of the constructed method. In its application to microstructures, the size effects of bending stiffness and stress concentration factor （SCF） are studied for microspeciem and microgripper, respectively. It is observed that the size effects become rather strong when the width of spring for microgripper, the radius of circular perforation and the long axis of elliptical perforation for microspeciem come close to the material characteristic length scales. For the U-shaped notch, the size effects decline obviously with increasing notch radius, and decline mildly with increasing length of notch.
基金This research was funded by Fundamental Research Grant Scheme(FRGS)under the Ministry of Higher Education Malaysia,grant number with project ref:FRGS/1/2019/STG06/UTP/03/2.
文摘Multistep integration methods are being extensively used in the simulations of high dimensional systems due to their lower computational cost.The block methods were developed with the intent of obtaining numerical results on numerous points at a time and improving computational efficiency.Hybrid block methods for instance are specifically used in numerical integration of initial value problems.In this paper,an optimized hybrid block Adams block method is designed for the solutions of linear and nonlinear first-order initial value problems in ordinary differential equations(ODEs).In deriving themethod,the Lagrange interpolation polynomial was employed based on some data points to replace the differential equation function and it was integrated over a specified interval.Furthermore,the convergence properties along with the region of stability of the method were examined.It was concluded that the newly derived method is convergent,consistent,and zero-stable.The method was also found to be A-stable implying that it covers the whole of the left/negative half plane.From the numerical computations of absolute errors carried out using the newly derived method,it was found that the method performed better than the ones with which we compared our results with.Themethod also showed its superiority over the existing methods in terms of stability and convergence.
文摘By combining the classical appropriate functions “1, x, x 2” with the method of multiplier enlargement, this paper establishes a theorem to approximate any unbounded continuous functions with modified positive linear operators. As an example, Hermite Fejér interpolation polynomial operators are analysed and studied, and a general conclusion is obtained.
基金the National Natural Science Foundation of China(No.11471062).
文摘In this paper,an optimality condition for nonlinear programming problems with box constraints is given by using linear transformation and Lagrange interpolating polynomials.Based on this condition,two new local optimization methods are developed.The solution points obtained by the new local optimization methods can improve the Karush–Kuhn–Tucker(KKT)points in general.Two global optimization methods then are proposed by combining the two new local optimization methods with a filled function method.Some numerical examples are reported to show the effectiveness of the proposed methods.
