Highly Accurate Golden Section Search Algorithms and Fictitious Time Integration Method for Solving Nonlinear Eigenvalue Problems

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.

Exponential Time Differencing Method for a Reaction-Diffusion System with Free Boundary

For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geometries may lead to difficulties in the accuracy when discretizing the high-order derivatives on grid points near the boundary.It is very challenging to design numerical methods that can efficiently and accurately handle both difficulties.Applying an implicit scheme may be able to remove the stability constraints on the time step,however,it usually requires solving a large global system of nonlinear equations for each time step,and the computational cost could be significant.Integration factor(IF)or exponential time differencing(ETD)methods are one of the popular methods for temporal partial differential equations(PDEs)among many other methods.In our paper,we couple ETD methods with an embedded boundary method to solve a system of reaction-diffusion equations with complex geometries.In particular,we rewrite all ETD schemes into a linear combination of specificФ-functions and apply one state-of-the-art algorithm to compute the matrix-vector multiplications,which offers significant computational advantages with adaptive Krylov subspaces.In addition,we extend this method by incorporating the level set method to solve the free boundary problem.The accuracy,stability,and efficiency of the developed method are demonstrated by numerical examples.

Integrating Geochemical and Geophysical Method in Coexistence of Oil and Potassium to Identify K-rich Brine:Research and Application in Southwestern Sichuan

Lithology of Triassic in southwestern Sichuan is consistent with the whole basin,and there is no discussion about stratum division,the difference is stratum denudation which is made by the uplifting of Luzhou uplift,especially

Multiresolution method for bending of plates with complex shapes

A high-accuracy multiresolution method is proposed to solve mechanics problems subject to complex shapes or irregular domains.To realize this method,we design a new wavelet basis function,by which we construct a fifth-order numerical scheme for the approximation of multi-dimensional functions and their multiple integrals defined in complex domains.In the solution of differential equations,various derivatives of the unknown function are denoted as new functions.Then,the integral relations between these functions are applied in terms of wavelet approximation of multiple integrals.Therefore,the original equation with derivatives of various orders can be converted to a system of algebraic equations with discrete nodal values of the highest-order derivative.During the application of the proposed method,boundary conditions can be automatically included in the integration operations,and relevant matrices can be assured to exhibit perfect sparse patterns.As examples,we consider several second-order mathematics problems defined on regular and irregular domains and the fourth-order bending problems of plates with various shapes.By comparing the solutions obtained by the proposed method with the exact solutions,the new multiresolution method is found to have a convergence rate of fifth order.The solution accuracy of this method with only a few hundreds of nodes can be much higher than that of the finite element method(FEM)with tens of thousands of elements.In addition,because the accuracy order for direct approximation of a function using the proposed basis function is also fifth order,we may conclude that the accuracy of the proposed method is almost independent of the equation order and domain complexity.

A SOLVING METHOD FOR A SYSTEM OF PARTIAL DIFFERENTIAL EQUATIONS WITH AN APPLICATION TO THE BENDING PROBLEM OF A THICK PLATE

A theorem of solving a system of linear non-homogeneous differential equations through integrating and adding its basic solutions is put forward and proved, the mathematical role and physical nature of the theorem is interpreted briefly. As an example, the theorem is applied to solve the problem of thermo-force bending of a thick plate.

Fuzzy Inventory Model under Selling Price Dependent Demand and Variable Deterioration with Fully Backlogged Shortages

The objective is to develop a model considering demand dependent on selling price and deterioration occurs after a certain period of time, which follows two-parameter Weibull distribution. Shortages are allowed and fully backlogged. Fuzzy optimal solution is obtained by considering hexagonal fuzzy numbers and for defuzzification Graded Mean Integration Representation Method. A numerical example is provided for the illustration of crisp and fuzzy, both models. To observe the effect of changes in parameters, sensitivity analysis is carried out.

Study on Characteristics of 3-D Translating-Pulsating Source Green Function of Deep-Water Havelock Form and Its Fast Integration Method

The singularities, oscillatory performances and the contributing factors to the 3-＇D translating-pulsating source Green function of deep-water Havelock form which consists of a local disturbance part and a far-field wave-like part, are analyzed systematically. Relative numerical integral methods about the two parts are presented in this paper. An improved method based on LOBATTO rule is used to eliminate singularities caused respectively by infinite discontinuity and jump discontinuous node from the local disturbance part function, which makes the improvement of calculation efficiency and accuracy possible. And variable substitution is applied to remove the singularity existing at the end of the integral interval of the far-field wave-like part function. Two auxiliary techniques such as valid interval calculation and local refinement of integral steps technique in narrow zones near false singularities are applied so as to avoid unnecessary integration of invalid interval and improve integral accordance. Numerical test results have proved the efficiency and accuracy in these integral methods that thus can be applied to calculate hydrodynamic performance of floating structures moving in waves.

PARAMETRIC EQUATIONS OF NONHOLONOMIC NONCONSERVATIVE SYSTEMS IN THE EVENT SPACE AND THE METHOD OF THEIR INTEGRATION

In this paper,the parametric equations with multipliers of nonholonomic nonconservative sys- tems in the event space are established,their properties are studied,and their explicit formulation is obtained. And then the field method for integrating these equations is given.Finally,an example illustrating the appli- cation of the integration method is given.

Fast precise integration method for hyperbolic heat conduction problems

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.

Evaluation of integration methods for hybrid simulation of complex structural systems through collapse

This study examines the performance of integration methods for hybrid simulation of large and complex structural systems in the context of structural collapse due to seismic excitations. The target application is not necessarily for real-time testing, but rather for models that involve large-scale physical sub-structures and highly nonlinear numerical models. Four case studies are presented and discussed. In the first case study, the accuracy of integration schemes including two widely used methods, namely, modified version of the implicit Newmark with fixed-number of iteration （iterative） and the operator-splitting （non-iterative） is examined through pure numerical simulations. The second case study presents the results of 10 hybrid simulations repeated with the two aforementioned integration methods considering various time steps and fixed-number of iterations for the iterative integration method. The physical sub-structure in these tests consists of a single-degree-of-freedom （SDOF） cantilever column with replaceable steel coupons that provides repeatable highly- nonlinear behavior including fracture-type strength and stiffness degradations. In case study three, the implicit Newmark with fixed-number of iterations is applied for hybrid simulations of a 1：2 scale steel moment frame that includes a relatively complex nonlinear numerical substructure. Lastly, a more complex numerical substructure is considered by constructing a nonlinear computational model of a moment frame coupled to a hybrid model ofa 1：2 scale steel gravity frame. The last two case studies are conducted on the same porotype structure and the selection of time steps and fixed number of iterations are closely examined in pre-test simulations. The generated unbalance forces is used as an index to track the equilibrium error and predict the accuracy and stability of the simulations.

Modified precise time step integration method of structural dynamic analysis

The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method （e.g., an algorithm with an arbitrary order of accuracy） is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.

Response Sensitivity Analysis of the Dynamic Milling Process Based on the Numerical Integration Method

As one of the bases of gradient-based optimization algorithms, sensitivity analysis is usually required to calculate the derivatives of the system response with respect to the machining parameters. The most widely used approaches for sensitivity analysis are based on time-consuming numerical methods, such as finite difference methods. This paper presents a semi-analytical method for calculation of the sensitivity of the stability boundary in milling. After transforming the delay-differential equation with time-periodic coefficients governing the dynamic milling process into the integral form, the Floquet transition matrix is constructed by using the numerical integration method. Then, the analytical expressions of derivatives of the Floquet transition matrix with respect to the machining parameters are obtained. Thereafter, the classical analytical expression of the sensitivity of matrix eigenvalues is employed to calculate the sensitivity of the stability lobe diagram. The two-degree-of-freedom milling example illustrates the accuracy and efficiency of the proposed method. Compared with the existing methods, the unique merit of the proposed method is that it can be used for analytically computing the sensitivity of the stability boundary in milling, without employing any finite difference methods. Therefore, the high accuracy and high efficiency are both achieved. The proposed method can serve as an effective tool for machining parameter optimization and uncertainty analysis in high-speed milling.

A COMBINED PARAMETRIC QUADRATIC PROGRAMMING AND PRECISE INTEGRATION METHOD BASED DYNAMIC ANALYSIS OF ELASTIC-PLASTIC HARDENING/SOFTENING PROBLEMS

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.

A FIELD METHOD FOR SOLVING THE EQUATIONS OF MOTION OF NONHOLONOMIC SYSTEMS

In this paper,the field method for solving the equations of motion of holonomic nonconservative systems is extended to nonholonomic systems with constant mass and with variable mass.Two examples are given to illustrate its application.

Isogeometric Boundary Element Analysis for 2D Transient Heat Conduction Problem with Radial Integration Method

This paper presents an isogeometric boundary element method(IGABEM)for transient heat conduction analysis.The Non-Uniform Rational B-spline(NURBS)basis functions,which are used to construct the geometry of the structures,are employed to discretize the physical unknowns in the boundary integral formulations of the governing equations.Bezier extraction technique is employed to accelerate the evaluation of NURBS basis functions.We adopt a radial integration method to address the additional domain integrals.The numerical examples demonstrate the advantage of IGABEM in dimension reduction and the seamless connection between CAD and numerical analysis.

Improved precise integration method for differential Riccati equation

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

Poisson theory and integration method of Birkhoffian systems in the event space

This paper focuses on studying the Poisson theory and the integration method of a Birkhoffian system in the event space. The Birkhoff＇s equations in the event space are given. The Poisson theory of the Birkhoffian system in the event space is established. The definition of the Jacobi last multiplier of the system is given, and the relation between the Jacobi last multiplier and the first integrals of the system is discussed. The researches show that for a Birkhoffian system in the event space, whose configuration is determined by （2n ＋ 1） Birkhoff＇s variables, the solution of the system can be found by the Jacobi last multiplier if 2n first integrals are known. An example is given to illustrate the application of the results.

Analytical mechanics methods for solving Whittaker equations

The purpose of this paper is to study the solution of the celebrated Whittaker equations by using analytical mechanics methods, including the Lagrange-Noether method, Hamilton-Poisson method and potential integral method.

A STUDY ON THE WEIGHT FUNCTION OF THE MOVING LEAST SQUARE APPROXIMATION IN THE LOCAL BOUNDARY INTEGRAL EQUATION METHOD

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.

Selection of the best initial orbital elements of satellite based on fuzzy integration evaluation method

The fuzzy integration evaluation method （FIEM） is studied in order to select the best orbital elements from the multi-group initial orbits determined by a satellite TT＆C （Tracking, Telemetry and Control） center with all kinds of data sources. By employing FIEM together with the experience of TT＆C experts, the index system to evaluate the selection of the best initial orbits is established after the data sources and orbit determination theories are studied. Besides, the concrete steps in employing the method are presented. Moreover, by taking the objects to be evaluated as evaluation experts, the problem of how to generate evaluation matrices is solved. Through practical application, the method to select the best initial orbital elements has been proved to be flexible and effective The originality of the method is to find a new evaluation criterion （comparing the actually tracked orbits） replacing the traditional one （comparing the nominal orbits） for selecting the best orbital elements.