The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data...The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data. However, because of its discretization, the FD method is only stable under certain conditions. The Arbitrary Difference Precise Integration （ADPI） method is based on the FD method and adopts an integration scheme in the time domain and an arbitrary difference scheme in the space domain. Therefore, the ADPI method is a semi-analytical method. In this paper, we deduce the formula for the ADPI method based on the 3D elastic equation and improve its stability. In forward modeling cases, the ADPI method was implemented in 2D and 3D elastic wave equation forward modeling. Results show that the travel time of the reflected seismic wave is accurate. Compared with the acoustic wave field, the elastic wave field contains more wave types, including PS- and PP- reflected waves, transmitted waves, and diffracted waves, which is important to interpretation of seismic data. The method can be easily applied to elastic wave equation numerical simulations for eoloical models.展开更多
A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to t...A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to the fact that the matrix exponential is sparse. The presented method employs the sparsity of the matrix exponential to improve the original precise integration method. The merits are that the proposed method is suitable for large hyperbolic heat equations and inherits the accuracy of the original version and the good computational efficiency, which are verified by two numerical examples.展开更多
The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcom...The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcome the result mesh-sensitivity problem in the dynamic strain softening or strain localization analysis. The equations for the dynamic elastic-plastic problems are derived in terms of the parametric variational principle, which is valid for associated, non-associated and strain softening plastic constitutive models in the finite element analysis. The precise integration method, which has been widely used for discretization in time domain of the linear problems, is introduced for the solution of dynamic nonlinear equations. The new algorithm proposed is based on the combination of the parametric quadratic programming method and the precise integration method and has all the advantages in both of the algorithms. Results of numerical examples demonstrate not only the validity, but also the advantages of the algorithm proposed for the numerical solution of nonlinear dynamic problems.展开更多
This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using prec...This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability.展开更多
The quasi-Shannon interval wavelet is constructed based on the interpolation wavelet theory, and an adaptive precise integration method, which is based on extrapolation method is presented for nonlinear ordinary diffe...The quasi-Shannon interval wavelet is constructed based on the interpolation wavelet theory, and an adaptive precise integration method, which is based on extrapolation method is presented for nonlinear ordinary differential equations ( ODEs). And then, an adaptive interval wavelet precise integration method (AIWPIM) for nonlinear partial differential equations(PDEs) is proposed. The numerical results show that the computational precision of AIWPIM is higher than that of the method constructed by combining the wavelet and the 4th Runge-Kutta method, and the computational amounts of these two methods are almost equal. For convenience, the Burgers equation is taken as an example in introducing this method, which is also valid for more general cases.展开更多
An improved precise integration method （IPIM） for solving the differential Riccati equation （DRE） is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise in...An improved precise integration method （IPIM） for solving the differential Riccati equation （DRE） is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise integration method （PIM） for solving the DRE is connected with the scaling and squaring method for computing the exponential of a matrix. The error analysis of the scaling and squaring method for the exponential of a matrix is applied to the PIM of the DRE. Based ,on the error analysis, the criterion for choosing two parameters of the PIM is given. Three kinds of IPIMs for solving the DRE are proposed. The numerical examples machine accuracy solutions. show that the IPIM is stable and gives the展开更多
General purpose graphic processing unit （GPU） calculation technology is gradually widely used in various fields. Its mode of single instruction, multiple threads is capable of seismic numerical simulation which has ...General purpose graphic processing unit （GPU） calculation technology is gradually widely used in various fields. Its mode of single instruction, multiple threads is capable of seismic numerical simulation which has a huge quantity of data and calculation steps. In this study, we introduce a GPU-based parallel calculation method of a precise integration method （PIM） for seismic forward modeling. Compared with CPU single-core calculation, GPU parallel calculating perfectly keeps the features of PIM, which has small bandwidth, high accuracy and capability of modeling complex substructures, and GPU calculation brings high computational efficiency, which means that high-performing GPU parallel calculation can make seismic forward modeling closer to real seismic records.展开更多
The difficulty in solving stochastic dynamics problems lies in the need for a large number of repeated computations of deterministic dynamic equations,which has been a challenge in stochastic dynamics analysis and was...The difficulty in solving stochastic dynamics problems lies in the need for a large number of repeated computations of deterministic dynamic equations,which has been a challenge in stochastic dynamics analysis and was discussed in this study.To efficiently and accurately compute the exponential of the dynamics state matrix and the matrix functions due to external loads,an adaptively filtered precise integration method was proposed,which inherits the high precision of the precise integrationmethod,improves the computational efficiency and saves the memory required.Moreover,the perturbation method was introduced to avoid repeated computations of matrix exponential and terms due to external loads.Based on the filtering and perturbation techniques,an adaptively filtered precise integration method considering perturbation for stochastic dynamics problems was developed.Two numerical experiments,including a model of phononic crystal and a bridge model considering random parameters,were performed to test the performance of the proposed method in terms of accuracy and efficiency.Numerical results show that the accuracy and efficiency of the proposed method are better than those of the existing precise integration method,the Newmark-βmethod and the Wilson-θmethod.展开更多
In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the met...In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the method of variable coefficient dimensional expanding,the non-homogeneous ordinary differential equations（ODEs） are transformed into homogeneous ODEs.Then the interval is divided evenly,and the transfer matrix in every subinterval is worked out using the high order multiple perturbation method,and a set of algebraic equations is given in the form of matrix by the precise integration relation for each segment,which is worked out by the reduction method.Finally numerical examples are elaboratedd to validate the present method.展开更多
This paper presents a precise method for solving singularly perturbed boundary-value problems with the boundary layer at one end. The method divides the interval evenly and gives a set of algebraic equations in a matr...This paper presents a precise method for solving singularly perturbed boundary-value problems with the boundary layer at one end. The method divides the interval evenly and gives a set of algebraic equations in a matrix form by the precise integration relationship of each segment. Substituting the boundary conditions into the algebraic equations, the coefficient matrix can be transformed to the block tridiagonal matrix. Considering the nature of the problem, an efficient reduction method is given for solving singular perturbation problems. Since the precise integration relationship introduces no discrete error in the discrete process, the present method has high precision. Numerical examples show the validity of the present method.展开更多
The Non-uniform rational B-spline （NURBS） enhanced scaled boundary finite element method in combination with the modified precise integration method is proposed for the transient heat conduction problems in this pap...The Non-uniform rational B-spline （NURBS） enhanced scaled boundary finite element method in combination with the modified precise integration method is proposed for the transient heat conduction problems in this paper. The scaled boundary finite element method is a semi-analytical technique, which weakens the governing differential equations along the circumferential direction and solves those analytically in the radial direction. In this method, only the boundary is discretized in the finite element sense leading to a re- duction of the spatial dimension by one with no fundamental solution required. Neverthe- less, in case of the complex geometry, a huge number of elements are generally required to properly approximate the exact shape of the domain and distorted meshes are often un- avoidable in the conventional finite element approach, which leads to huge computational efforts and loss of accuracy. NURBS are the most popular mathematical tool in CAD industry due to its flexibility to fit any free-form shape. In the proposed methodology, the arbitrary curved boundary of problem domain is exactly represented with NURBS basis functions, while the straight part of the boundary is discretized by the conventional Lagrange shape functions. Both the concepts of isogeometric analysis and scaled boundary finite element method are combined to form the governing equations of transient heat conduction analy- sis and the solution is obtained using the modified precise integration method. The stiffness matrix is obtained from a standard quadratic eigenvalue problem and the mass matrix is determined from the low-frequency expansion. Finally the governing equations become a system of first-order ordinary differential equations and the time domain response is solved numerically by the modified precise integration method. The accuracy and stability of the proposed method to deal with the transient heat conduction problems are demonstrated by numerical examples.展开更多
In view of the limitation of the difference method,the adjustment model of CPⅢprecise trigonometric leveling control network based on the parameter method was proposed in the present paper.The experiment results show...In view of the limitation of the difference method,the adjustment model of CPⅢprecise trigonometric leveling control network based on the parameter method was proposed in the present paper.The experiment results show that this model has a simple algorithm and high data utilization,avoids the negative influences caused by the correlation among the data acquired from the difference method and its accuracy is improved compared with the difference method.In addition,the strict weight of CPⅢprecise trigonometric leveling control network was also discussed in this paper.The results demonstrate that the ranging error of trigonometric leveling can be neglected when the vertical angle is less than 3 degrees.The accuracy of CPⅢprecise trigonometric leveling control network has not changed significantly before and after strict weight.展开更多
With a three-dimensional semiclassical ensemble method, we theoretically investigated the nonsequential double ionization of Ar driven by the spatially inhomogeneous few-cycle negatively chirped laser pulses. Our resu...With a three-dimensional semiclassical ensemble method, we theoretically investigated the nonsequential double ionization of Ar driven by the spatially inhomogeneous few-cycle negatively chirped laser pulses. Our results show that the recollision time window can be precisely controlled within an isolated time interval of several hundred attoseconds, which is useful for understanding the subcycle correlated electron dynamics. More interestingly, the correlated electron momentum distribution (CEMD) exhibits a strong dependence on laser intensity. That is, at lower laser intensity, CEMD is located in the first quadrant. As the laser intensity increases,CEMD shifts almost completely to the second and fourth quadrants, and then gradually to the third quadrant.The underlying physics governing the CEMD's dependence on laser intensity is explained.展开更多
The HY-2 satellite carrying a satellite-borne GPS receiver is the first Chinese radar altimeter satellite, whose radial orbit determination precision must reach the centimeter level. Now HY-2 is in the test phase so t...The HY-2 satellite carrying a satellite-borne GPS receiver is the first Chinese radar altimeter satellite, whose radial orbit determination precision must reach the centimeter level. Now HY-2 is in the test phase so that the observations are not openly released. In order to study the precise orbit determination precision and procedure for HY-2 based on the satellite- borne GPS technique, the satellite-borne GPS data are simulated in this paper. The HY-2 satellite-borne GPS antenna can receive at least seven GPS satellites each epoch, which can validate the GPS receiver and antenna design. What＇s more, the precise orbit determination processing flow is given and precise orbit determination experiments are conducted using the HY-2-borne GPS data with both the reduced-dynamic method and the kinematic geometry method. With the 1 and 3 mm phase data random errors, the radial orbit determination precision can achieve the centimeter level using these two methods and the kinematic orbit accuracy is slightly lower than that of the reduced-dynamic orbit. The earth gravity field model is an important factor which seriously affects the precise orbit determination of altimeter satellites. The reduced-dynamic orbit determination experiments are made with different earth gravity field models, such as EIGEN2, EGM96, TEG4, and GEMT3. Using a large number of high precision satellite-bome GPS data, the HY-2 precise orbit determination can reach the centimeter level with commonly used earth gravity field models up to above 50 degrees and orders.展开更多
After the trajectory simulation model of rudder control rocket with six degrees of freedom is established by Matlab/ Simulink, the simulated targeting of rudder control rocket with rudder angle error and starting cont...After the trajectory simulation model of rudder control rocket with six degrees of freedom is established by Matlab/ Simulink, the simulated targeting of rudder control rocket with rudder angle error and starting control moment error is carried out respectively by means of Monte Carlo method and the distribution of impact points of rudder control rocket is counted from all the successful subsamples. In the case of adding interference errors associated with rudder angle error and starting time error, the simulation analysis of impact point dispersion is done and its lateral and longitudinal correction abilities at different targeting angles are simulated to identify the effects of these factors on characteristics and control precision of the rudder control rocket, which provides the relevant reference for high-precision design of rudder control system.展开更多
This paper presents an improved precise integration algorithm fortransient analysis of heat transfer and some other problems. Theoriginal precise integration method is improved by means of the inve-rse accuracy analys...This paper presents an improved precise integration algorithm fortransient analysis of heat transfer and some other problems. Theoriginal precise integration method is improved by means of the inve-rse accuracy analysis so that the parameter N, which has been takenas a constant and an independent pa- rameter without consideration ofthe problems in the original method, can be generated automaticallyby the algorithm itself.展开更多
Nonlinear dynamic equations can be solved accurately using a precise integration method. Some algorithms exist, but the inversion of a matrix must be calculated for these al- gorithms. If the inversion of the matrix d...Nonlinear dynamic equations can be solved accurately using a precise integration method. Some algorithms exist, but the inversion of a matrix must be calculated for these al- gorithms. If the inversion of the matrix doesn’t exist or isn’t stable, the precision and stability of the algorithms will be afected. An explicit series solution of the state equation has been pre- sented. The solution avoids calculating the inversion of a matrix and its precision can be easily controlled. In this paper, an implicit series solution of nonlinear dynamic equations is presented. The algorithm is more precise and stable than the explicit series solution and isn’t sensitive to the time-step. Finally, a numerical example is presented to demonstrate the efectiveness of the algorithm.展开更多
In this paper, the spectral element method(SEM)is improved to solve the moving load problem. In this method, a structure with uniform geometry and material properties is considered as a spectral element, which means t...In this paper, the spectral element method(SEM)is improved to solve the moving load problem. In this method, a structure with uniform geometry and material properties is considered as a spectral element, which means that the element number and the degree of freedom can be reduced significantly. Based on the variational method and the Laplace transform theory, the spectral stiffness matrix and the equivalent nodal force of the beam-column element are established. The static Green function is employed to deduce the improved function. The proposed method is applied to two typical engineering practices—the one-span bridge and the horizontal jib of the tower crane. The results have revealed the following. First, the new method can yield extremely high-precision results of the dynamic deflection, the bending moment and the shear force in the moving load problem.In most cases, the relative errors are smaller than 1%. Second, by comparing with the finite element method, one can obtain the highly accurate results using the improved SEM with smaller element numbers. Moreover, the method can be widely used for statically determinate as well as statically indeterminate structures. Third, the dynamic deflection of the twin-lift jib decreases with the increase in the moving load speed, whereas the curvature of the deflection increases.Finally, the dynamic deflection, the bending moment and the shear force of the jib will all increase as the magnitude of the moving load increases.展开更多
The investigation was carried out on the technical problems of finishing the inner surface of elbow parts and the action mechanism of particles in elbow precision machining by abrasive flow.This work was analyzed and ...The investigation was carried out on the technical problems of finishing the inner surface of elbow parts and the action mechanism of particles in elbow precision machining by abrasive flow.This work was analyzed and researched by combining theory,numerical and experimental methods.The direct simulation Monte Carlo(DSMC)method and the finite element analysis method were combined to reveal the random collision of particles during the precision machining of abrasive flow.Under different inlet velocity,volume fraction and abrasive particle size,the dynamic pressure and turbulence flow energy of abrasive flow in elbow were analyzed,and the machining mechanism of particles on the wall and the influence of different machining parameters on the precision machining quality of abrasive flow were obtained.The test results show the order of the influence of different parameters on the quality of abrasive flow precision machining and establish the optimal process parameters.The results of the surface morphology before and after the precision machining of the inner surface of the elbow are discussed,and the surface roughness Ra value is reduced from 1.125μm to 0.295μm after the precision machining of the abrasive flow.The application of DSMC method provides special insights for the development of abrasive flow technology.展开更多
High-speed and precision positioning are fundamental requirements for high-acceleration low-load mechanisms in integrated circuit （IC） packaging equipment. In this paper, we derive the transient nonlinear dynamicres...High-speed and precision positioning are fundamental requirements for high-acceleration low-load mechanisms in integrated circuit （IC） packaging equipment. In this paper, we derive the transient nonlinear dynamicresponse equations of high-acceleration mechanisms, which reveal that stiffness, frequency, damping, and driving frequency are the primary factors. Therefore, we propose a new structural optimization and velocity-planning method for the precision positioning of a high-acceleration mechanism based on optimal spatial and temporal distribution of inertial energy. For structural optimization, we first reviewed the commonly flexible multibody dynamic optimization using equivalent static loads method （ESLM）, and then we selected the modified ESLM for optimal spatial distribution of inertial energy; hence, not only the stiffness but also the inertia and frequency of the real modal shapes are considered. For velocity planning, we developed a new velocity-planning method based on nonlinear dynamic-response optimization with varying motion conditions. Our method was verified on a high-acceleration die bonder. The amplitude of residual vibration could be decreased by more than 20% via structural optimization and the positioning time could be reduced by more than 40% via asymmetric variable velocity planning. This method provides an effective theoretical support for the precision positioning of high-acceleration low-load mechanisms.展开更多
基金supported by the National Science and Technology Major Project of China(Grant No. 2011ZX05004-003,2011ZX05014-006-006)the National Key Basic Research Program of China(Grant No. 2013CB228602)the Natural Science Foundation of China(Grant No. 40974066)
文摘The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data. However, because of its discretization, the FD method is only stable under certain conditions. The Arbitrary Difference Precise Integration （ADPI） method is based on the FD method and adopts an integration scheme in the time domain and an arbitrary difference scheme in the space domain. Therefore, the ADPI method is a semi-analytical method. In this paper, we deduce the formula for the ADPI method based on the 3D elastic equation and improve its stability. In forward modeling cases, the ADPI method was implemented in 2D and 3D elastic wave equation forward modeling. Results show that the travel time of the reflected seismic wave is accurate. Compared with the acoustic wave field, the elastic wave field contains more wave types, including PS- and PP- reflected waves, transmitted waves, and diffracted waves, which is important to interpretation of seismic data. The method can be easily applied to elastic wave equation numerical simulations for eoloical models.
基金supported by the National Natural Science Foundation of China (Nos. 10902020 and 10721062)
文摘A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to the fact that the matrix exponential is sparse. The presented method employs the sparsity of the matrix exponential to improve the original precise integration method. The merits are that the proposed method is suitable for large hyperbolic heat equations and inherits the accuracy of the original version and the good computational efficiency, which are verified by two numerical examples.
文摘The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcome the result mesh-sensitivity problem in the dynamic strain softening or strain localization analysis. The equations for the dynamic elastic-plastic problems are derived in terms of the parametric variational principle, which is valid for associated, non-associated and strain softening plastic constitutive models in the finite element analysis. The precise integration method, which has been widely used for discretization in time domain of the linear problems, is introduced for the solution of dynamic nonlinear equations. The new algorithm proposed is based on the combination of the parametric quadratic programming method and the precise integration method and has all the advantages in both of the algorithms. Results of numerical examples demonstrate not only the validity, but also the advantages of the algorithm proposed for the numerical solution of nonlinear dynamic problems.
文摘This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability.
文摘The quasi-Shannon interval wavelet is constructed based on the interpolation wavelet theory, and an adaptive precise integration method, which is based on extrapolation method is presented for nonlinear ordinary differential equations ( ODEs). And then, an adaptive interval wavelet precise integration method (AIWPIM) for nonlinear partial differential equations(PDEs) is proposed. The numerical results show that the computational precision of AIWPIM is higher than that of the method constructed by combining the wavelet and the 4th Runge-Kutta method, and the computational amounts of these two methods are almost equal. For convenience, the Burgers equation is taken as an example in introducing this method, which is also valid for more general cases.
基金Project supported by the National Natural Science Foundation of China(Nos.10902020 and 10721062)
文摘An improved precise integration method （IPIM） for solving the differential Riccati equation （DRE） is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise integration method （PIM） for solving the DRE is connected with the scaling and squaring method for computing the exponential of a matrix. The error analysis of the scaling and squaring method for the exponential of a matrix is applied to the PIM of the DRE. Based ,on the error analysis, the criterion for choosing two parameters of the PIM is given. Three kinds of IPIMs for solving the DRE are proposed. The numerical examples machine accuracy solutions. show that the IPIM is stable and gives the
基金supported by the National Natural Science Foundation of China (Nos 40974066 and 40821062)National Basic Research Program of China (No 2007CB209602)
文摘General purpose graphic processing unit （GPU） calculation technology is gradually widely used in various fields. Its mode of single instruction, multiple threads is capable of seismic numerical simulation which has a huge quantity of data and calculation steps. In this study, we introduce a GPU-based parallel calculation method of a precise integration method （PIM） for seismic forward modeling. Compared with CPU single-core calculation, GPU parallel calculating perfectly keeps the features of PIM, which has small bandwidth, high accuracy and capability of modeling complex substructures, and GPU calculation brings high computational efficiency, which means that high-performing GPU parallel calculation can make seismic forward modeling closer to real seismic records.
基金the support of the National Natural Science Foundation of China(Grant Nos.11472067 and 51609034)the Science Foundation of Liaoning Province of China(No.2021-MS-119)+1 种基金the Dalian Youth Science and Technology Star Project(No.2018RQ06)the Fundamental Research Funds for the Central Universities(Grant No.DUT20GJ216).
文摘The difficulty in solving stochastic dynamics problems lies in the need for a large number of repeated computations of deterministic dynamic equations,which has been a challenge in stochastic dynamics analysis and was discussed in this study.To efficiently and accurately compute the exponential of the dynamics state matrix and the matrix functions due to external loads,an adaptively filtered precise integration method was proposed,which inherits the high precision of the precise integrationmethod,improves the computational efficiency and saves the memory required.Moreover,the perturbation method was introduced to avoid repeated computations of matrix exponential and terms due to external loads.Based on the filtering and perturbation techniques,an adaptively filtered precise integration method considering perturbation for stochastic dynamics problems was developed.Two numerical experiments,including a model of phononic crystal and a bridge model considering random parameters,were performed to test the performance of the proposed method in terms of accuracy and efficiency.Numerical results show that the accuracy and efficiency of the proposed method are better than those of the existing precise integration method,the Newmark-βmethod and the Wilson-θmethod.
基金supported by the National Natural Science Foundation of China (11132004 and 51078145)the Natural Science Foundation of Guangdong Province (9251064101000016)
文摘In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the method of variable coefficient dimensional expanding,the non-homogeneous ordinary differential equations（ODEs） are transformed into homogeneous ODEs.Then the interval is divided evenly,and the transfer matrix in every subinterval is worked out using the high order multiple perturbation method,and a set of algebraic equations is given in the form of matrix by the precise integration relation for each segment,which is worked out by the reduction method.Finally numerical examples are elaboratedd to validate the present method.
基金Project supported by the National Natural Science Foundation of China(No.10672194)the China-Russia Cooperative Project(the National Natural Science Foundation of China and the Russian Foundation for Basic Research)(No.10811120012)
文摘This paper presents a precise method for solving singularly perturbed boundary-value problems with the boundary layer at one end. The method divides the interval evenly and gives a set of algebraic equations in a matrix form by the precise integration relationship of each segment. Substituting the boundary conditions into the algebraic equations, the coefficient matrix can be transformed to the block tridiagonal matrix. Considering the nature of the problem, an efficient reduction method is given for solving singular perturbation problems. Since the precise integration relationship introduces no discrete error in the discrete process, the present method has high precision. Numerical examples show the validity of the present method.
基金support by the National Natural Science Foundation of China(grant No.51779033,51409038)the National Key Research and Development Plan(grant No.2016YFB0201001)the National Natural Science Foundation of China(grant No.51421064)
文摘The Non-uniform rational B-spline （NURBS） enhanced scaled boundary finite element method in combination with the modified precise integration method is proposed for the transient heat conduction problems in this paper. The scaled boundary finite element method is a semi-analytical technique, which weakens the governing differential equations along the circumferential direction and solves those analytically in the radial direction. In this method, only the boundary is discretized in the finite element sense leading to a re- duction of the spatial dimension by one with no fundamental solution required. Neverthe- less, in case of the complex geometry, a huge number of elements are generally required to properly approximate the exact shape of the domain and distorted meshes are often un- avoidable in the conventional finite element approach, which leads to huge computational efforts and loss of accuracy. NURBS are the most popular mathematical tool in CAD industry due to its flexibility to fit any free-form shape. In the proposed methodology, the arbitrary curved boundary of problem domain is exactly represented with NURBS basis functions, while the straight part of the boundary is discretized by the conventional Lagrange shape functions. Both the concepts of isogeometric analysis and scaled boundary finite element method are combined to form the governing equations of transient heat conduction analy- sis and the solution is obtained using the modified precise integration method. The stiffness matrix is obtained from a standard quadratic eigenvalue problem and the mass matrix is determined from the low-frequency expansion. Finally the governing equations become a system of first-order ordinary differential equations and the time domain response is solved numerically by the modified precise integration method. The accuracy and stability of the proposed method to deal with the transient heat conduction problems are demonstrated by numerical examples.
基金National Natural Science Foundation of China(No.41661091)。
文摘In view of the limitation of the difference method,the adjustment model of CPⅢprecise trigonometric leveling control network based on the parameter method was proposed in the present paper.The experiment results show that this model has a simple algorithm and high data utilization,avoids the negative influences caused by the correlation among the data acquired from the difference method and its accuracy is improved compared with the difference method.In addition,the strict weight of CPⅢprecise trigonometric leveling control network was also discussed in this paper.The results demonstrate that the ranging error of trigonometric leveling can be neglected when the vertical angle is less than 3 degrees.The accuracy of CPⅢprecise trigonometric leveling control network has not changed significantly before and after strict weight.
基金supported by the National Natural Science Foundation of China (Grant No. 12074329)Nanhu Scholars Program for Young Scholars of Xinyang Normal University。
文摘With a three-dimensional semiclassical ensemble method, we theoretically investigated the nonsequential double ionization of Ar driven by the spatially inhomogeneous few-cycle negatively chirped laser pulses. Our results show that the recollision time window can be precisely controlled within an isolated time interval of several hundred attoseconds, which is useful for understanding the subcycle correlated electron dynamics. More interestingly, the correlated electron momentum distribution (CEMD) exhibits a strong dependence on laser intensity. That is, at lower laser intensity, CEMD is located in the first quadrant. As the laser intensity increases,CEMD shifts almost completely to the second and fourth quadrants, and then gradually to the third quadrant.The underlying physics governing the CEMD's dependence on laser intensity is explained.
基金supported partially by the National Natural Science Foundation of China (Nos. 40974004 and 40974016)Key Laboratory of Dynamic Geodesy of CAS, China (No. L09-01) R&I Team Support Program and the Graduate Science and Technology Foundation of SDUST, China (No. YCA110403)
文摘The HY-2 satellite carrying a satellite-borne GPS receiver is the first Chinese radar altimeter satellite, whose radial orbit determination precision must reach the centimeter level. Now HY-2 is in the test phase so that the observations are not openly released. In order to study the precise orbit determination precision and procedure for HY-2 based on the satellite- borne GPS technique, the satellite-borne GPS data are simulated in this paper. The HY-2 satellite-borne GPS antenna can receive at least seven GPS satellites each epoch, which can validate the GPS receiver and antenna design. What＇s more, the precise orbit determination processing flow is given and precise orbit determination experiments are conducted using the HY-2-borne GPS data with both the reduced-dynamic method and the kinematic geometry method. With the 1 and 3 mm phase data random errors, the radial orbit determination precision can achieve the centimeter level using these two methods and the kinematic orbit accuracy is slightly lower than that of the reduced-dynamic orbit. The earth gravity field model is an important factor which seriously affects the precise orbit determination of altimeter satellites. The reduced-dynamic orbit determination experiments are made with different earth gravity field models, such as EIGEN2, EGM96, TEG4, and GEMT3. Using a large number of high precision satellite-bome GPS data, the HY-2 precise orbit determination can reach the centimeter level with commonly used earth gravity field models up to above 50 degrees and orders.
文摘After the trajectory simulation model of rudder control rocket with six degrees of freedom is established by Matlab/ Simulink, the simulated targeting of rudder control rocket with rudder angle error and starting control moment error is carried out respectively by means of Monte Carlo method and the distribution of impact points of rudder control rocket is counted from all the successful subsamples. In the case of adding interference errors associated with rudder angle error and starting time error, the simulation analysis of impact point dispersion is done and its lateral and longitudinal correction abilities at different targeting angles are simulated to identify the effects of these factors on characteristics and control precision of the rudder control rocket, which provides the relevant reference for high-precision design of rudder control system.
基金the National Natural Science Foundation of China (No.19872016,19872017)the National Key Basic Research Special Foundation (G1999032805)the Foundation for University Key Teachers by the Ministry of Education of China
文摘This paper presents an improved precise integration algorithm fortransient analysis of heat transfer and some other problems. Theoriginal precise integration method is improved by means of the inve-rse accuracy analysis so that the parameter N, which has been takenas a constant and an independent pa- rameter without consideration ofthe problems in the original method, can be generated automaticallyby the algorithm itself.
基金Project supported by the National Natural Science Foundation of China(Nos.60273048and60174023).
文摘Nonlinear dynamic equations can be solved accurately using a precise integration method. Some algorithms exist, but the inversion of a matrix must be calculated for these al- gorithms. If the inversion of the matrix doesn’t exist or isn’t stable, the precision and stability of the algorithms will be afected. An explicit series solution of the state equation has been pre- sented. The solution avoids calculating the inversion of a matrix and its precision can be easily controlled. In this paper, an implicit series solution of nonlinear dynamic equations is presented. The algorithm is more precise and stable than the explicit series solution and isn’t sensitive to the time-step. Finally, a numerical example is presented to demonstrate the efectiveness of the algorithm.
基金supported by the National Key Technology R&D Program (Grant 2011BAJ02B01-02)the National Natural Science Foundation of China (Grant 11602065)
文摘In this paper, the spectral element method(SEM)is improved to solve the moving load problem. In this method, a structure with uniform geometry and material properties is considered as a spectral element, which means that the element number and the degree of freedom can be reduced significantly. Based on the variational method and the Laplace transform theory, the spectral stiffness matrix and the equivalent nodal force of the beam-column element are established. The static Green function is employed to deduce the improved function. The proposed method is applied to two typical engineering practices—the one-span bridge and the horizontal jib of the tower crane. The results have revealed the following. First, the new method can yield extremely high-precision results of the dynamic deflection, the bending moment and the shear force in the moving load problem.In most cases, the relative errors are smaller than 1%. Second, by comparing with the finite element method, one can obtain the highly accurate results using the improved SEM with smaller element numbers. Moreover, the method can be widely used for statically determinate as well as statically indeterminate structures. Third, the dynamic deflection of the twin-lift jib decreases with the increase in the moving load speed, whereas the curvature of the deflection increases.Finally, the dynamic deflection, the bending moment and the shear force of the jib will all increase as the magnitude of the moving load increases.
基金Projects(51206011,U1937201)supported by the National Natural Science Foundation of ChinaProject(20200301040RQ)supported by the Science and Technology Development Program of Jilin Province,China+1 种基金Project(JJKH20190541KJ)supported by the Education Department of Jilin Province,ChinaProject(18DY017)supported by Changchun Science and Technology Program of Changchun City,China。
文摘The investigation was carried out on the technical problems of finishing the inner surface of elbow parts and the action mechanism of particles in elbow precision machining by abrasive flow.This work was analyzed and researched by combining theory,numerical and experimental methods.The direct simulation Monte Carlo(DSMC)method and the finite element analysis method were combined to reveal the random collision of particles during the precision machining of abrasive flow.Under different inlet velocity,volume fraction and abrasive particle size,the dynamic pressure and turbulence flow energy of abrasive flow in elbow were analyzed,and the machining mechanism of particles on the wall and the influence of different machining parameters on the precision machining quality of abrasive flow were obtained.The test results show the order of the influence of different parameters on the quality of abrasive flow precision machining and establish the optimal process parameters.The results of the surface morphology before and after the precision machining of the inner surface of the elbow are discussed,and the surface roughness Ra value is reduced from 1.125μm to 0.295μm after the precision machining of the abrasive flow.The application of DSMC method provides special insights for the development of abrasive flow technology.
基金supported by the National Key Basic Research Program of China (2011CB013104)National Natural Science Foundation of China (U1134004)+2 种基金Guangdong Provincial Natural Science Foundation (2015A030312008)Science and Technology Program of Guangzhou (201510010281)Guangdong Provincial Science and Technology Plan (2013B010402014)
文摘High-speed and precision positioning are fundamental requirements for high-acceleration low-load mechanisms in integrated circuit （IC） packaging equipment. In this paper, we derive the transient nonlinear dynamicresponse equations of high-acceleration mechanisms, which reveal that stiffness, frequency, damping, and driving frequency are the primary factors. Therefore, we propose a new structural optimization and velocity-planning method for the precision positioning of a high-acceleration mechanism based on optimal spatial and temporal distribution of inertial energy. For structural optimization, we first reviewed the commonly flexible multibody dynamic optimization using equivalent static loads method （ESLM）, and then we selected the modified ESLM for optimal spatial distribution of inertial energy; hence, not only the stiffness but also the inertia and frequency of the real modal shapes are considered. For velocity planning, we developed a new velocity-planning method based on nonlinear dynamic-response optimization with varying motion conditions. Our method was verified on a high-acceleration die bonder. The amplitude of residual vibration could be decreased by more than 20% via structural optimization and the positioning time could be reduced by more than 40% via asymmetric variable velocity planning. This method provides an effective theoretical support for the precision positioning of high-acceleration low-load mechanisms.
