Study of differential control method for solving chaotic solutions of nonlinear dynamic system

In this paper, we focus on the need to solve chaotic solutions of high-dimensional nonlinear dynamic systems of which the analytic solution is difficult to obtain. For this purpose, a Differential Control Method (DCM) is proposed based on the Mechanized Mathematics-Wu Elimination Method (WEM). By sampling, the computer time of the differential operator of the nonlinear differential equation can be substituted by the differential quotient of solving the variable time of the sample. The main emphasis of DCM is placed on substituting the differential quotient of a small neighborhood of the sample time of the computer system for the differential operator of the equations studied. The approximate analytical chaotic solutions of the nonlinear differential equations governing the high-dimensional dynamic system can be obtained by the method proposed. In order to increase the computational efficiency of the method proposed, a thermodynamics modeling method is used to decouple the variable and reduce the dimension of the system studied. The validity of the method proposed for obtaining approximate analytical chaotic solutions of the nonlinear differential equations is illustrated by the example of a turbo-generator system. This work is applied to solving a type of nonlinear system of which the dynamic behaviors can be described by nonlinear differential equations.     

Nonlinear implicit iterative method for solving nonlinear ill-posed problems

In the paper, we extend the implicit iterative method for linear ill-posed operator equations to solve nonlinear ill-posed problems. We show that under some conditions the error sequence of solutions of the nonlinear implicit iterative method is monotonically decreasing and, with this monotonicity, prove convergence of the new method for both the exact and perturbed equations.     

轮对非线性动力学系统蛇行运动的解析解

在对铁路车辆系统的极限环幅值和非线性临界速度进行分析时通常采用数值方法, 不便于研究其随系统参数的变化规律. 轮对系统保留了影响车辆系统动力学性能的几个关键要素: 如轮轨几何非线性约束、轮轨接触蠕滑关系和悬挂系统等, 可以反映铁路车辆系统蛇行运动的本质特性. 轮对系统自由度少、参数少, 可以采用解析方法进行分析. 本文选取合适的特征量把轮对非线性动力学方程无量纲化, 得到了带有小参数的两自由度微分方程; 采用多尺度方法对该方程进行了解析求解; 给出了轮对系统极限环幅值的解析表达式并对其稳定性进行了判定; 给出了轮对系统的分岔速度解析表达式, 并进而获得系统的非线性临界速度的解析表达式. 在对得到的解析解用数值结果进行验证后, 用得到的解析解进行了系统参数影响分析. 传统的分岔图计算方法(如降速法、路径跟踪法等)需对微分方程进行大量数值积分计算方可求解系统的非线性临界速度值, 而通过本文获得的解析表达式可直接给出系统的非线性临界速度值和极限环幅值, 便于研究轮对系统动力学特性随参数的变化规律,进行快速方案比对和筛选, 为转向架结构优化设计提供参考.      

Eigenfunction expansion method and its application to two-dimensional elasticity problems based on stress formulation

This paper proposes an eigenfunction expansion method to solve twodimensional （2D） elasticity problems based on stress formulation. By introducing appropriate state functions, the fundamental system of partial differential equations of the above 2D problems is rewritten as an upper triangular differential system. For the associated operator matrix, the existence and the completeness of two normed orthogonal eigenfunction systems in some space are obtained, which belong to the two block operators arising in the operator matrix. Moreover, the general solution to the above 2D problem is given by the eigenfunction expansion method.     

一类偏微分方程的Hamilton正则表示

主要给出一系列关于力学中的偏微分方程的无穷维Ｈａｍｉｌｔｏｎ正则表示．其中包括变系数线性偏微分方程,ＫｄＶ方程,ＭＫｄＶ方程,ＫＰ方程,Ｂｏｕｓｉｎｅｓｑ方程等的无穷维Ｈａｍｉｌｔｏｎ正则表示．     

HIGH-ORDER DISCONTINUOUS GALERKIN SOLUTION OF N-S EQUATIONS ON HYBRID MESH

针对层流NS方程发展了混合网格上的高阶间断有限元方法,给出了物面边界高阶近似的具体步骤以及近物面弯曲单元的处理方法。对数值离散产生的非线性方程组采用牛顿迭代进行求解,每个牛顿循环采用预处理广义最小余量法求解产生的大型稀疏线性系统。使用该方法得到了典型算例的数值结果,并跟前人的计算结果进行了比较。计算结果表明,混合网格上应用高阶间断有限元方法求解黏性流动具有很好的应用前景。     

Eigenfunction expansion method of upper triangular operator matrix and application to two-dimensional elasticity problems based on stress formulation

This paper studies the eigenfunction expansion method to solve the two-dimensional(2D) elasticity problems based on the stress formulation.The fundamentalsystem of partial differential equations of the 2D problems is rewritten as an upper tri-angular differential system based on the known results,and then the associated uppertriangular operator matrix matrix is obtained.By further research,the two simpler com-plete orthogonal systems of eigenfunctions in some space are obtained,which belong tothe two block operators arising in the operator matrix.Then,a more simple and conve-nient general solution to the 2D problem is given by the eigenfunction expansion method.Furthermore,the boundary conditions for the 2D problem,which can be solved by thismethod,are indicated.Finally,the validity of the obtained results is verified by a specificexample.     

General solution for large deflection problems of nonhomogeneous circular plates on elastic foundation

In this paper, a new method is presented based on [1]. It can be used to solve the arbitrary nonlinear system of differential equations with variable coefficients. By this method, the general solution for large deformation of nonhomogeneous circular plates resting on an elastic foundation is derived. The convergence of the solution is proved. Finally, it is only necessary to solve a set of nonlinear algebraic equations with three unknowns. The solution obtained by the present method has large convergence range and the computation is simpler and more rapid than other numerical methods.Numerical examples given at the end of this paper indicate that satisfactory results of stress resullants and displacements can be obtained by the present method. The correctness of the theory in this paper is, confirmed.     

Direct Integration Methods with Integral Model for Dynamic Systems

ＩｎｔｒｏｄｕｃｔｉｏｎＴｒａｎｓｉｅｎｔｓｔａｔｅａｎａｌｙｓｉｓｈａｓｂｅｅｎａｎａｃｔｉｖｅｒｅｓｅａｒｃｈａｒｅａｉｎｍａｎｙｅｎｇｉｎｅｅｒｉｎｇｐｒｏｂｌｅｍｓ.Ｉｎｐｒａｃｔｉｃａｌｓｉｔｕａｔｉｏｎ ,ｌｉｋｅｓｔｒｕｃｔｕｒａｌｍｅｃｈａｎｉｃｓ,ｔｈｅｓｙｓｔｅｍｓｂｅｉｎｇｓｔｕｄｉｅｄａｒｅｕｓｕａｌｌｙｎｏｎｌｉｎｅａｒａｎｄｔｉｍｅ＿ｄｅｐｅｎｄｅｎｔ.Ｔｈｅａｎａｌｙｔｉｃａｌｍｅｔｈｏｄｓｔｏｓｏｌｖｅｔｈｅｓｅｐｒｏｂｌｅｍｓａｒｅｉｍｐｏｓｓｉｂｌｅ .Ｔｈｅｎｕｍ…     

Hilbert伴随算子逆定理在有限元法中的应用

固体力学或其它学科的大量问题均归结为求解偏微分方程组。本文把 Hilbert 伴随算子逆定理用于有限元法,求解非正定和正定偏微分方程组。它可以带有任意变系数及复杂的边界条件。文小给出了收敛性证明。并给出统一的计算公式。利用本文的方法,可以给出一个非协调有限元。单元之间的协调连续条件仅需在节点上满足,因此很容易处理。和一般的有限元法相比,有更高的精度。文末给出算例,表明利用本文方法获得的解可以收敛于精确解,并有较高的收敛精度。     

Numerical solution of stiff systems from nonlinear dynamics using single-term Haar wavelet series

The paper presents single-term Haar wavelet series (STHWS) approach to the solution of nonlinear stiff differential equations arising in nonlinear dynamics. The properties of STHWS are given. The method of implementation is discussed. Numerical solutions of some model equations are investigated for their stiffness and stability and solutions are obtained to demonstrate the suitability and applicability of the method. The results in the form of block-pulse and discrete solutions are given for typical nonlinear stiff systems. As compared with the TR BDF2 method of Shampine and Gill’s method, the STHWS turns out to be more effective in its ability to solve systems ranging from mildly to highly stiff equations and is free from stability constraints.     

Analysis of a two-grid method for semiconductor device problem

The mathematical model of a semiconductor device is governed by a system of quasi-linear partial differential equations.The electric potential equation is approximated by a mixed finite element method,and the concentration equations are approximated by a standard Galerkin method.We estimate the error of the numerical solutions in the sense of the Lqnorm.To linearize the full discrete scheme of the problem,we present an efficient two-grid method based on the idea of Newton iteration.The main procedures are to solve the small scaled nonlinear equations on the coarse grid and then deal with the linear equations on the fine grid.Error estimation for the two-grid solutions is analyzed in detail.It is shown that this method still achieves asymptotically optimal approximations as long as a mesh size satisfies H=O(h^1/2).Numerical experiments are given to illustrate the efficiency of the two-grid method.     

Projection-algebraic approximation of linear and nonlinear operator differential equations in Banach spaces

We develop a Calogero-type projection-algebraic method of discrete approximations for linear differential equations in Banach spaces and analyze the convergence of finite-dimensional approximations based on the functional-analytic approach to discrete approximations and methods of operator theory in Banach spaces. Applications of the obtained results to the functional-interpolation scheme of the projection-algebraic method of discrete approximations are considered. Based on a generalized Leray–Schauder-type theorem, we consider the projection-algebraic scheme of discrete approximations and analyze its solvability and convergence for a special class of nonlinear operator equations.     

Construction of Basis Functions and Their Application to Boundary-Value Problems of Mechanics of Continuous Media

A unified method for constructing basis (eigen) functions is proposed to solve problems of mechanics of continuous media, problems of cubature and quadrature, and problems of approximation of hypersurfaces. Numerical-analytical methods are described, which allow obtaining approximate solutions of internal and external boundary-value problems of mechanics of continuous media of a certain class (both linear and nonlinear). The method is based on decomposition of the sought solutions of the considered partial differential equations into series in basis functions. An algorithm is presented for linearization of partial differential equations and reduction of nonlinear boundary-value problems, which are reduced to systems of linear algebraic equations with respect to unknown coefficients without using traditional methods of linearization.     

The Generalized Nonlinear Galerkin Method Applied to Hydrodynamics of the Open Sea

In this paper we present new numerical algorithms based on a generalized nonlinear Galerkin method in order to solve coastal and oceanic circulation problems. The equations system is based on the primitive equations of the ocean under Boussinesq and hydrostatic approximations. These equations are transformed using, at the same time, the classical σ transformation and an original homogenization of the boundary conditions. We use a well adapted special basis to apply the usual Galerkin method and the nonlinear Galerkin method. This basis is built on a modelization of the energetic transfers through the different scales of flow. Two approaches are proposed to solve the continuity equation: the (nonlinear) Galerkin method and the method of the characteristics. We present the advantages and drawbacks of both methods.     

METHOD OF GREEN＇S FUNCTION OF CORRUGATED SHELLS

By using the fundamental equations of axisymmetric shallow shells of revolution, the nonlinear bending of a shallow corrugated shell with taper under arbitrary load has been investigated. The nonlinear boundary value problem of the corrugated shell was reduced to the nonlinear integral equations by using the method of Green＇s function. To solve the integral equations, expansion method was used to obtain Green＇s function. Then the integral equations were reduced to the form with degenerate core by expanding Green＇s function as series of characteristic function. Therefore, the integral equations become nonlinear algebraic equations. Newton＇ s iterative method was utilized to solve the nonlinear algebraic equations. To guarantee the convergence of the iterative method, deflection at center was taken as control parameter. Corresponding loads were obtained by increasing deflection one by one. As a numerical example,elastic characteristic of shallow corrugated shells with spherical taper was studied.Calculation results show that characteristic of corrugated shells changes remarkably. The snapping instability which is analogous to shallow spherical shells occurs with increasing load if the taper is relatively large. The solution is close to the experimental results.     

Optimal control of quadratic functionals for affine nonlinear systems

In this paper we analyze the optimal control problem for a class of affine nonlinear systems under the assumption that the associated Lie algebra is nilpotent. The Lie brackets generated by the vector fields which define the nonlinear system represent a remarkable mathematical instrument for the control of affine systems. We determine the optimal control which corresponds to the nilpotent operator of the first order. In particular, we obtain the control that minimizes the energy of the given nonlinear system. Applications of this control to bilinear systems with first order nilpotent operator are considered.     

Sorting-free Hill-based stability analysis of periodic solutions through Koopman analysis

In this paper, we aim to study nonlinear time-periodic systems using the Koopman operator, which provides a way to approximate the dynamics of a nonlinear system by a linear time-invariant system of higher order. We propose for the considered system class a specific choice of Koopman basis functions combining the Taylor and Fourier bases. This basis allows to recover all equations necessary to perform the harmonic balance method as well as the Hill analysis directly from the linear lifted dynamics. The key idea of this paper is using this lifted dynamics to formulate a new method to obtain stability information from the Hill matrix. The error-prone and computationally intense task known by sorting, which means identifying the best subset of approximate Floquet exponents from all available candidates, is circumvented in the proposed method. The Mathieu equation and an n-DOF generalization are used to exemplify these findings.

     

Statics of nonlinear elastic cable-stayed systems with slipping flexible cables

A design model is proposed to describe spatial cable-stayed systems with slipping cables. The basic variables of the displacement method are selected. A system of nonlinear equilibrium equations is derived from the condition of minimum potential energy, taking into account the large displacements of joints and the large deformations of cables. A method and formulas needed to iteratively solve the system of nonlinear equations are presented. Individual consideration is given to a cable net with prescribed tension. The computed results make it possible to fabricate and mark out cable rods and to assemble cable nets based on their geometry alone, without the need for force measurement __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 5, pp. 79–87, May 2006.