非光滑优化的双层规划模型及内点算法

对于包含接触约束的非光滑结构优化问题,其非光滑性体现在状态函数并不是处处可微的,针对含有应力约束及接触约束的非光滑结构优化问题,建立了一种双层规划模型,避免了求解时非光滑性所带来的问题,同时提出了一种迭代算法,用对偶内点二次规划进行分析,线性规划进行优化,算例表明这种方法十分有效。     

NON-INTERIOR SMOOTHING ALGORITHM FOR FRICTIONAL CONTACT PROBLEMS

A new algorithm for solving the three-dimensional elastic contact problem with friction is presented. The algorithm is a non-interior smoothing algorithm based on an NCP-function. The parametric variational principle and parametric quadratic programming method were applied to the analysis of three-dimensional frictional contact problem. The solution of the contact problem was finally reduced to a linear complementarity problem, which was reformulated as a system of nonsmooth equations via an NCP-function. A smoothing approximation to the nonsmooth equations was given by the aggregate function. A Newton method was used to solve the resulting smoothing nonlinear equations. The algorithm presented is easy to understand and implement. The reliability and efficiency of this algorithm are demonstrated both by the numerical experiments of LCP in mathematical way and the examples of contact problems in mechanics.     

Chaos synchronization of discrete-time dynamic systems with a limited capacity communication channel

This paper offers a new control strategy for discrete-time chaos synchronization where the drive system and the response system are coupled via a limited capacity communication channel (LCCC for short). One simple condition is presented to ensure synchronization between the two chaotic systems coupled by a LCCC. Based on this condition, an explicit coder–decoder pair for the coding algorithm is designed and the synchronization error between the two chaotic systems decays to zero exponentially based on this coding algorithm. Finally, the proposed control strategy is applied to the well-known H\′{e}non system, and numerical simulations illustrate the validity of the obtained result.     

粘弹性阻尼修改后结构系统动特性的预测

无阻尼或比例阻尼结构系统经粘弹阻尼修改后,可变为粘弹性阻尼系统。要获得其动特性,需求解复特征值问题。但是,随之带来了计算量大、费用高等问题。尤其是粘弹性材料特性随频率而变化,需求解高阶非线性复特征值问题,这对于一个自由度较大的结构,计算量太大,通常较难实现。本文在特征值修改方法的基础上,提出粘弹性阻尼局部动力修改方法,即仅需已知原结构系统的实模态参数,就可求出粘弹阻尼修改后系统的复模态参数。还发展了特征值和修改量同时迭代方法,有效地解决了粘弹材料复模量随频率变化引起的非线性复特征值问题。     

振动系统动力学设计迭代算法及解的存在性讨论

振动系统动力学设计被抽象为高维广义非线性特征值反问题。若系统构成以可变参数表示，则可构造一个非线性多元函数。基于多元函数极小原理，提出了一套求解这一特征值反问题的迭代算法。该算法不受系统给定阶数和给定方向的限制，也适用于具有重特征值的退化情况，系统或结构的构成材料可以是任意的。文中同时讨论了解的存在条件，且以显式表达，可方便地应用于工程实际。结合某直升机旋翼桨叶的动力学设计，给出了应用的数值算例。大量数值仿真结果及应用实践表明，本文算法具很好的收敛性，并有较高的计算精度。     

Using Fractional Order Adjustment Rules and Fractional Order Reference Models in Model-Reference Adaptive Control

This paper investigates the use of Fractional Order Calculus (FOC) inconventional Model Reference Adaptive Control (MRAC) systems. Twomodifications to the conventional MRAC are presented, i.e., the use offractional order parameter adjustment rule and the employment offractional order reference model. Through examples, benefits from theuse of FOC are illustrated together with some remarks for furtherresearch.     

Observer-based model reference adaptive control for unknown time-delay chaotic systems with input nonlinearity

Existence of unknown time-delay in the systems is a drastic restriction that it can menace the stability criteria and even deteriorate the performance system. This undesired case would be more intensified if that the uncertain input nonlinearity effects are also considered. To handle the input nonlinearities effects (results in dead-zone and/or hysteresis phenomena) and also unknown time-delay in the chaotic systems, this paper presents an observer-based Model Reference Adaptive Control (MRAC) scheme for a class of unknown time-delay chaotic systems with disturbances. This new method is a delay-independent variable-structure control method which is integrated with an observer system. The main task of the proposed approach is to accomplish a perfect tracking procedure such that unknown parameters are adapted via output estimation error. Furthermore, stability of the closed-loop system is achieved by means of the Lyapunov stability theory. Finally, the proposed methods are applied to some famous chaotic systems to verify the effectiveness of the proposed methods.     

Formulation and Solution of Dynamic Problems of Elastic Rod Systems Subjected to Boundary Conditions Described by Multivalued Relations

The dynamic behavior of rod systems under the action of external force factors described by multivalued (subdifferential) relations is studied. The mathematical formulation of the problem is given in the form of a dynamic quasivariational inequality. With the use of the Newmark difference scheme, successive approximations, and finite-element discretization, the problem is reduced to minimization of a convex nonsmooth finite-dimensional functional with respect to velocities at each time step. Introduction of auxiliary variables by the method of a modified Lagrangian reduces the problem of minimization of this functional to a sequence of smooth problems of nonlinear programming. The algorithm is verified using the numerical solution for a problem with one degree of freedom. The algorithm proposed is used to calculate the rods of deep-well pumps.     

Symplectic solution system for reissner plate bending

Based on the Hellinger-Reissner variatonal principle for Reissner plate bendingand introducing dual variables, Hamiltonian dual equations for Reissner plate bending werepresented. Therefore Hamiltonian solution system can also be applied to Reissner platebending problem, and the transformation from Euclidian space to symplectic space and fromLagrangian system to Hamiltonian system was realized. So in the symplectic space whichconsists of the original variables and their dual variables, the problem can be solved viaeffective mathematical physics methods such as the method of separation of variables andeigenfunction-vector expansion. All the eigensolutions and Jordan canonical formeigensolutions for zero eigenvalue of the Hamiltonian operator matrix are solved in detail, and their physical meanings are showed clearly. The adjoint symplectic orthonormal relation of the eigenfunction vectors for zero eigenvalue are formed. It is showed that the alleigensolutions for zero eigenvalue are basic solutions of the Saint-Venant problem and theyform a perfect symplectic subspace for zero eigenvalue. And the eigensolutions for nonzeroeigenvalue are covered by the Saint-Venant theorem. The symplectic solution method is notthe same as the classical semi-inverse method and breaks through the limit of the traditional semi-inverse solution. The symplectic solution method will have vast application.     

Double Eigenvalue Problems for Hemivariational Inequalities

The aim of the present paper is to study a new type of eigenvalue problem, called a double eigenvalue problem, which arises in hemivariational inequalities related to nonconvex nonsmooth energy functionals. The paper provides existence results as well as some qualitative properties for the solutions to double eigenvalue problems for hemivariational inequalities under the presence of given nonlinear compact operators which are not necessarily of a variational structure. It presents three different approaches to such problems: minimization, minimax methods and (sub) critical point theory on a sphere. Applications illustrate the theory. (Accepted April 2, 1996)     

Numerical Continuation of Homoclinic Orbits to Non-Hyperbolic Equilibria in Planar Systems

In this paper we develop numerical algorithms for thecontinuation of degenerate homoclinic orbits to non-hyperbolicequilibria in planar systems. The first situation corresponds to asaddle-node equilibrium (a zero eigenvalue) and the second one is theso-called cuspidal loop (double-zero eigenvalue). The methods proposedmay deal with codimension-two and -three homoclinic connections.Application of the algorithms to several examples supports its validityand demonstrates its usefulness.     

重分析边界处理与广义主从控制关系

本文讨论了在结构重分析中局部网格细分时，修改区域与未修改区域交界处的位移连续性问题．文中提出了广义主从控制的概念，把边界上的新增节点作为从节点，同时受多点控制，保证了结构的连续性，并对有限元结构中的常见单元给出了边界的位移控制关系，而且此关系已应用于DDJ程序系统，实现了算法及程序设计．     

多体系统摩擦接触问题的DAE-LCP方法

当多体系统的约束全部是摩擦接触时, 其动力学问题可归结为一个常微分方程(ordinarydifferential equation, ODE)与线性互补问题(linear complementarityproblem, LCP)的混合动力学问题. 如果除了摩擦接触之外还增加了光滑的双边约束, 则需要将ODE-LCP混合动力学模型推广为微分代数方程(differential algebra equation, DAE)与LCP的混合动力学模型. 该文采用DAE与LCP混合动力学方法求解不考虑碰撞但同时含有持续摩擦接触及光滑等式约束的多体系统动力学问题. 在建立系统动力学模型时,首先将含摩擦的约束从系统中移去得到基本动力学系统. 由于基本系统中带有等式约束, 所以基本系统的动力学方程为一组DAE. 结合基本系统的DAE与约束的互补条件便可以得到DAE-LCP混合动力学模型. 数值计算采用基于DAE与LCP的步进(time-stepping)算法, 将系统动力学方程及其约束离散化并转化为一个混合LCP进行求解. 该算法无需进行滞-滑状态检测, 避免了事件检测导致的繁复计算. 利用所提方法对典型机构的非光滑非线性特征进行了数值分析,验证了该文方法的正确有效性.      

受迫简谐响应分析中的特征值问题

无阻尼结构的受迫振动的共振频率与自由振动的特征值直接相关。在频域响应谱中,共振频率对应于响应峰值位置。指出频谱中的低谷(相对最小值)对应的频率也可用特征值问题求解。当最小值为0时,对应的频率是著名的反共振频率。另一种可能是,处于两个共振频率之间存在非零的最小响应,对应的频率称为最小响应频率。基于特征值问题的列式,反共振频率或最小响应频率的灵敏度分析可以直接通过已有的特征值灵敏度分析方法求解。给出了详细的数学推导并通过数值算例验证。     

Event-triggered robust adaptive control for uncertain nonlinear systems preceded by actuator dead-zone

It is both theoretically and practically important to investigate the problem of event-triggered adaptive tracking control for a class of uncertain nonlinear systems subject to actuator dead-zone, which aims at reducing communication rate and compensating actuator nonlinearity simultaneously. In this paper, to handle such a problem, an event-trigger based adaptive compensation scheme is proposed for the system preceded by actuator dead-zone. The challenges of this work can be roughly classified into two categories: how to compensate the nonsmooth dead-zone nonlinearity and how to eliminate the quantization signal effects caused by event-triggered strategy. To resolve the first challenge, a new decomposition of dead-zone mathematical model is employed so that dead-zone nonlinearity can be successively compensated by using robust approach. In addition, an adaptive controller and its triggering event are co-designed based on the relative threshold strategy, such that an asymptotic tracking performance can be ensured. The proposed scheme is proved to guarantee the globally bounded of all closed-loop signals and the asymptotic convergence performance of tracking error toward zero. The simulation results illustrate the effectiveness of our proposed control scheme.     

Order reduction of structural dynamic systems with static piecewise linear nonlinearities

A technique for order reduction of dynamic systems in structural form with static piecewise linear nonlinearities is presented. By utilizing two methods which approximate the nonlinear normal mode (NNM) frequencies and mode shapes, reduced-order models are constructed which more accurately represent the dynamics of the full model than do reduced models obtained via standard linear transformations. One method builds a reduced-order model which is dependent on the amplitude (initial conditions) while the other method results in an amplitude-independent reduced model. The two techniques are first applied to reduce two-degree-of-freedom undamped systems with clearance, deadzone, bang-bang, and saturation stiffness nonlinearities to single-mode reduced models which are compared by direct numerical simulation with the full models. It is then shown via a damped four-degree-of-freedom system with two deadzone nonlinearities that one of the proposed techniques allows for reduction to multi-mode reduced models and can accommodate multiple nonsmooth static nonlinearities with several surfaces of discontinuity. The advantages of the proposed methods include obtaining a reduced-order model which is signal-independent (doesn’t require direct integration of the full model), uses a subset of the original physical coordinates, retains the form of the nonsmooth nonlinearities, and closely tracks the actual NNMs of the full model.     

分层介质中Love波的扩展W-W算法

应用精细积分法(PIM)和扩展Wittrick-Williams(W-W)算法求解横观各向同性分层半空间中的Love波问题.Love波对应于波数-频率域线性常微分方程的本征值问题.精细积分法是求解线性常微分方程两端边值问题和初值问题的高精度算法.利用本征值计数技术,扩展W-W算法可以不遗漏地找到所有本征值.因此,文中使用的方法可以得到计算机精度意义下的精确解.     

Methods for obtaining sufficient conditions for the stability of autonomous conservative systems

A computational method for obtaining sufficient conditions for the stability of the stationary solution of autonomous conservative systems is proposed in the paper. This method is adapted to linear autonomous gyroscopic systems with three degrees of freedom. It is based on the positive definiteness of a parametric quadratic form composed of the gyroscopic force matrices and the potential function. The control parameters for the stability of the zero solution of the gyroscopic system are the entries of the gyroscopic force matrix. The algorithm of the computational method includes estimating one gyroscopic force parameter in the equation constructed from a necessary stability condition.A special example is used to demonstrate the application of this algorithm. Comparison is performed with some well-known methods for obtaining sufficient conditions on the basis of an incomplete set of first integrals of motion. It is shown that the positive definiteness of the modified potential energy may result in stable as well as unstable motions.     

Analytical and numerical methods of symplectic system for Stokes flow in two-dimensional rectangular domain

In this paper,a new analytical method of symplectic system.Hamiltonian system,is introduced for solving the problem of the Stokes flow in a two-dimensional rectangular domain.In the system,the fundamental problem is reduced to all eigenvalue and eigensolution problem.The solution and boundary conditions call be expanded by eigensolutions using ad.ioint relationships of the symplectic ortho-normalization between the eigensolutions.A closed method of the symplectic eigensolution is presented based on completeness of the symplectic eigensolution space.The results show that fundamental flows can be described by zero eigenvalue eigensolutions,and local effects by nonzero eigenvalue eigensolutions.Numerical examples give various flows in a rectangular domain and show effectivenees of the method for solving a variety of problems.Meanwhile.the method can be used in solving other problems.     

Analytical and numerical methods of symplectic system for Stokes flow in two-dimensional rectangular domain

In this paper, a new analytical method of symplectic system, Hamiltonian system, is introduced for solving the problem of the Stokes flow in a two-dimensional rectangular domain. In the system, the fundamental problem is reduced to an eigenvalue and eigensolution problem. The solution and boundary conditions can be expanded by eigensolutions using adjoint relationships of the symplectic ortho-normalization between the eigensolutions. A closed method of the symplectic eigensolution is presented based on completeness of the symplectic eigensolution space. The results show that fundamental flows can be described by zero eigenvalue eigensolutions, and local effects by nonzero eigenvalue eigensolutions. Numerical examples give various flows in a rectangular domain and show effectiveness of the method for solving a variety of problems. Meanwhile, the method can be used in solving other problems.