关于非线性动力系统的数值追踪

一个依赖于参数的动力系统，当参数改变时，如何追踪它的行为，特别是如何得到它与初条件无关的稳定解集，是人们长期以来很感兴趣的问题，它不仅对离散动力系统有意义，对于连续动力系统，如结构、流体流动等也都可以通过离散化这类问题来处理，近廿年来，笔地得及其他作者从事这类问题的研究。并取得了一些初步结果。这里我们作一个简要的介绍。     

解析与数值结合法及其应用

评述了解析与数值结合法发展的基础、研究的内容、分类、特点及应用和发展前景。     

三维对流问题的拟协调六面体单元解法

寻找一种高精度的空间单元插值模式是数值求解三维对流问题的关键。在前人研究的基础上,探讨了一种任意空间六面体的拟协凋单元,保证节点上的物理量函数及其一阶导数连续。算例表明,该方法具有良好的计算稳定性和低数值阻尼的优点,且计算工作量大大小于协调单元法,有利于推广应用于对流扩散方程的数值求解。     

一种建筑材料细观力学数值模拟的新方法

数值流形方法通过数学和物理双重网格,分析连续和非连续问题,已应用于模拟节理岩体裂隙的开裂与闭合问题.但对于裂纹尖端的局部化问题,数值流形方法需要像有限元那样在裂纹尖端设置细密单元.本文利用裂纹尖端解析解将数值流形方法的基函数进行扩展,推导了相应的试函数.从最小势能原理出发提出了断裂力学的数值流形方法,推导了相应的求解方程,将其应用于建筑材料细观力学数值模拟.最后给出两个数值算例,将计算结果与解析解对比,说明该方法的正确性和可行性.     

透射边界条件在波动谱元模拟中的实现:一维波动

多次透射公式(multi-transmitting formula,MTF)是一种具有普适性的局部人工边界条件,但其在近场波动数值模拟中一般与有限元法结合.由于波动谱元模拟的数值格式与有限元格式有极大的不同,传统的MTF在谱元离散格式中无法直接实现.为了使物理概念清楚、精度可控的多次透射人工边界条件能够适应波动谱元模拟的需求,首先指出多次透射边界与谱元离散格式结合的基本问题,并分析了空间内插和时间内插两种方案的可行性.然后从空间内插角度出发,提出基于拉格朗日多项式插值模式的MTF谱元格式,并采用一种简单内插方法实现高阶MTF.最后通过一维波动数值试验检验这些MTF谱元格式的精度,并讨论其数值稳定性.结果表明:对于一、二阶MTF,几种格式的精度相当;对于三、四阶MTF,基于谱单元位移模式插值的格式精度最高.相反,随着插值多项式阶次的升高,不同MTF格式的稳定临界值逐步降低,但是所有格式均在人工波速大大超过物理波速时才可能发生失稳.     

Reconstruction of high order derivatives by new mollification methods

In this paper, the problem of reconstructing numerical derivatives from noisy data is considered. A new framework of mollification methods based on the L generalized solution regularization methods is proposed. A specific algorithm for the first three derivatives is presented in the paper, in which a modification of TSVD, termed cTSVD is chosen as the regularization technique. Numerical examples given in the paper verify the theoretical results and show efficiency of the new method.     

数值流形方法研究及应用进展

基于有限覆盖技术的数值流形方法是一种新的广义的数值方法.该方法的场函数近似原理和有限元、无网格、单位分解等方法相似,但在网格划分、覆盖形式、近似函数等方面有其自身的特点和优势.对该方法近年来在理论研究和应用方面取得的重要进展进行了综述.在理论研究方面,日前已对不同形式物理覆盖流形单元的性能进行了研究,结果表明流形单元的精度较有限单元高,且提高覆盖函数的阶次能提高单元的精度;同时理论研究已由二维低阶流形方法推广到三维高阶流形方法,由线性流形方法推广到非线性流形方法,由基于能量原理的流形方法推广到基于加权余量的流形方法,非协调流形方法、无网格流形方法等也已开展了研究;此外,覆盖系统的自动生成、覆盖函数的形式以及边界条件的处理方法等流形方法相关理论的研究也取得了进展.在应用方面,开展了有关岩石破坏和裂纹扩展等非连续变形分析更深入的研究,并已逐步推广到金属塑性变形分析、多孔介质变形分析以及温度场的数值分析等多个领域.针对日前流形方法的研究和应用现状,该文展望了流形方法理论及实现方法的研究方向、及其在计算流体力学、金属成形等大变形问题、多物理场分析等领域的应用前景.     

分数算子描述的黏弹性体力学问题数值方法

讨论由黎曼-刘维尔 (Riemann-Liouville)分数导数描述的黏弹性体力学问题的数 值方法. 该方法利用黎曼-刘维尔分数 导数定义中核函数的特性，并结合被积函数在单步中的逼近以及Newmark型数值法，建立 了分数导数计算公式. 算例表明，该算法具有收敛快、精度高、稳定性好和易于 应用和改进的优点. 在对动态系统的瞬态响应分析和有限元分析格式中，算法都获得了 满意的结果.     

数值流形方法的粘性边界问题初探

在实际工程数值流形方法分析中,采用固定约束边界的方法处理无限域或者半无限域的情况,边界处应力波的反射造成模拟结果与实际情况不符.本文基于Lysmer等人提出的粘性边界理论,在边界上设置阻尼器,推导相应粘性边界条件下流形单元刚度矩阵的数值计算格式,经岩石长条中弹性波传播算例,并与有限元结果对比,验证了该粘性边界的有效性,有利于数值流形方法的工程中推广应用.     

Study on the Prediction Method of Low-Dimension Time Series that Arise from the Intrinsic Nonlinear Dynamics

The prediction methods and its applications of the nonlinear dynamic systems determined from chaotic time series of low-dimension are discussed mainly. Based on the work of the foreign researchers, the chaotic time series in the phase space adopting one kind of nonlinear chaotic model were reconstructed. At first, the model parameters were estimated by using the improved least square method. Then as the precision was satisfied, the optimization method was used to estimate these parameters. At the end by using the obtained chaotic model, the future data of the chaotic time series in the phase space was predicted. Some representative experimental examples were analyzed to testify the models and the algorithms developed in this paper. The results show that if the algorithms developed here are adopted, the parameters of the corresponding chaotic model will be easily calculated well and true. Predictions of chaotic series in phase space make the traditional methods change from outer iteration to interpolations. And if the optimal model rank is chosen, the prediction precision will increase notably. Long term superior predictability of nonlinear chaotic models is proved to be irrational and unreasonable. Paper from Chen Yu-shu, Member of Editorial of Committee, AMM Foundation item: the National Natural Science Foundation of China (19990510); the National Key Basic Research Special Fund(G1998020316) Biography: Ma Jun-hai(1965-), Professor, Doctor     

Optimal nonlinear observers for chaotic synchronization with message embedded

This paper presents an optimal nonlinear observer for synchronizing the transmitter-receiver pair with guaranteed optimal performance. In the proposed scheme, a generalized nonlinear state-space observer via uniform matrix transformations is constructed to estimate the transmitter state and the information signal, simultaneously. A nonlinear optimal design approach is used to synchronize chaotic systems. Solving the Hamilton–Jacobi–Bellman (H–J–B) equations we can obtain a linear optimal feedback scheme for piecewise-linear chaotic systems. Moreover, a robust scheme derived from the H _∞ optimization theory improves the synchronization performance of general nonlinear chaotic systems by suppressing the influence of their high order residual terms. Finally, two numerical simulation examples are illustrated by the chaotic Chua’s circuit system and the Lorenz chaotic system to demonstrate the effectiveness of our scheme.     

MIXED FINITE ELEMENT METHODS FOR THE SHALLOW WATER EQUATIONS INCLUDING CURRENT AND SILT SEDIMENTATION （Ⅰ）——THE CONTINUOUS-TIME CASE

An initial-boundary value problem for shallow equation system consisting of water dynamics equations, silt transport equation, the equation of bottom topography change, and of some boundary and initial conditions is studied, the existence of its generalized solution and semidiscrete mixed finite element (MFE) solution was discussed, and the error estimates of the semidiscrete MFE solution was derived. The error estimates are optimal.     

MODAL SYNTHESIS METHOD FOR NORM COMPUTATION OF H∞ DECENTRALIZED CONTROL SYSTEMS （II）

When using H∞ techniques to design decentralized controllers for large systems, the whole system is divided into subsystems, which are analysed using H∞ control theorybefore being recombined. An analogy was established with substructural analysis instructural mechanics, in which H∞ decentralized control theory corresponds to substructuralmodal synthesis theory so that the optimal H∞ norm of the whole system corresponds to thefundamental vibration frequency of the whole structure. Hence, modal synthesismethodology and the extended Wittrick-Williams algorithm were transplanted from structuralmechanics to compute the optimal H∞ norm of the control system. The orthogonality and theexpansion theorem of eigenfunctions of the subsystems H~ control are presented in part (I)of the paper. The modal synthesis method for computation of the optimal H∞ norm ofdecentralized control systems and numerical examples are presented in part (Ⅱ).     

MODAL SYNTHESIS METHOD FOR NORM COMPUTATION OF H∞ DECENTRALIZED CONTROL SYSTEMS （ I ）

When using H∞ techniques to design decentralized controllers for large systems,the whole system is divided into subsystems, which are analysed using H∞ control theorybefore being recombined. An analogy was established with substructural analysis instructural mechanics, in which H∞ decentralized control theory corresponds to substructuralmodal synthesis theory so that the optimal H∞ norm of the whole system corresponds to thefundamental vibration frequency of the whole structure. Hence, modal synthesismethodology and the extended Wittrick-Williams algorithm were transplanted from structuralmechanics to compute the optimal H∞ norm of the control system. The orthogonality and theexpansion theorem of eigenfunctions of the subsystems H∞ control are presented in part(I) of the paper. The modal synthesis method for computation of the optimal H∞ norm ofdecentralized control systems and numerical examples are presented in part (Ⅱ).     

Modelling, Dynamics and Control of Tethered Satellite Systems

Tethered satellite systems (TSS) pose quite challenging problems concerning their modelling, derivation of the equations of motion, numerical simulation of their dynamics, deciding on stability of relative equilibria provided the system moves on a circular orbit around the Earth and the occurrence of chaotic dynamics. Moreover, for the processes of deployment or retrieval of one satellite from or to another satellite certain control strategies, for example time or energy optimal control, are necessary. All these problems are considered in this paper.     

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. The project supported by the National Key Basic Research Special Foundation (G1999032805), the National Natural Science Foundation of China (19872016, 50178016, 19832010) and the Foundation for University Key Teacher by the Ministry of Education of China     

Observer-based synchronization scheme for a class of chaotic systems using contraction theory

In this paper, an adaptive synchronization scheme is proposed for a class of nonlinear systems. The design utilizes an adaptive observer, which is quite useful in establishing a transmitter–receiver kind of synchronization scheme. The proposed approach is based on contraction theory and provides a very simple way of establishing exponential convergence of observer states to actual system states. The class of systems addressed here has uncertain parameters, associated with the part of system dynamics that is a function of measurable output only. The explicit conditions for the stability of the observer are derived in terms of gain selection of the observer. Initially, the case without uncertainty is considered and then the results are extended to the case with uncertainty in parameters of the system. An application of the proposed approach is presented to synchronize the family of N chaotic systems which are coupled through the output variable only. The numerical results are presented for designing an adaptive observer for the chaotic Chua system to verify the efficacy of the proposed approach. Explicit bounds on observer gains are derived by exploiting the properties of the chaotic attractor exhibited by Chua’s system. Convergence of uncertain parameters is also analyzed for this case and numerical simulations depict the convergence of parameter estimates to their true value.     

SDRE-based near optimal nonlinear controller design for unified chaotic systems

This paper proposes a near optimal controller design method for unified chaotic systems based on state-dependent Riccati equation (SDRE) approach. A parameterization of the optimal nonlinear control gain is given in terms of the solution matrix of an SDRE. A simple algorithm to compute the near optimal control gain is proposed. The proposed near optimal control design method is also extended to the synchronization problem for unified chaotic systems. Finally, the effectiveness of the proposed design method is verified via numerical simulations.     

Retrospective time integral scheme and its applications to the advection equation

To put more information into a difference scheme of a differential equation for making an accurate prediction, a new kind of time integration scheme, known as the retrospective (RT) scheme, is proposed on the basis of the memorial dynamics. Stability criteria of the scheme for an advection equation in certain conditions are derived mathematically. The computations for the advection equation have been conducted with its RT scheme. It is shown that the accuracy of the scheme is much higher than that of the leapfrog (LF) difference scheme. The project supported by the National Key Program for Developing Basic Sciences (G1999043408 and G1998040901-1) and the National Natural Sciences Foundation of China (40175024 and 40035010)