Parameter estimation for chaotic systems using a geometric approach: theory and experiment

A method for estimating model parameters based on chaotic system response data is described. This estimation problem is made challenging by sensitive dependence to initial conditions. The standard maximum likelihood estimation method is practically infeasible due to the non-smooth nature of the likelihood function. We bypass the problem by introducing an alternative, smoother function that admits a better-defined maximum and show that the parameters that maximize this new function are asymptotically equivalent to maximum likelihood estimates. We use simulations to explore the influence of noise and available data on model Duffing and Lorenz oscillators. We then apply the approach to experimental data from a chaotic Duffing system. Our method does not require estimation of initial conditions and parameter estimates may be obtained even when system dynamics have been estimated from a delay embedding.     

Maximum likelihood Newton recursive and the Newton iterative estimation algorithms for Hammerstein CARAR systems

This paper discusses the identification problems of Hammerstein controlled autoregressive autoregressive (CARAR) systems using the maximum likelihood principle and Newton optimization method. A Newton recursive algorithm and a Newton iterative algorithm using the maximum likelihood principle are presented. The simulation results show that the proposed algorithms can effectively estimate the parameters of the Hammerstein CARAR systems.     

基于S变换的非平稳随机过程演变功率谱密度估计

提出了一种基于S变换的估计Priestley非平稳随机过程演变功率谱密度的方法。此方法的根本在于,相对于S变换的"变换核",Priestley非平稳随机过程的调制函数为慢变函数。因此,非平稳随机过程的S变换可视为相位修正后的另一非平稳随机过程。推导出了对应于特定频率点的S变换瞬时均方值和非平稳随机过程演变功率谱密度之间的关系式。将功率谱密度函数表达为有限个频率点的级数展开,通过求解一组代数方程,就能得到级数展开中每个频率点的时变系数,由此,可给出非平稳随机过程的演变功率谱密度。由于级数展开中的高斯形状函数不依赖于时间,因此,本文所提算法具有较高的计算效率。最后,给出了均匀调制和非均匀调制非平稳随机过程演变功率谱估计的算例。     

稀薄流到连续流的气体运动论模型方程算法研究

通过引入碰撞松弛参数和当地平衡态分布函数对BGK模型方程进行修正，确定含流态控制参数可描述不同流域气体流动特性的气体分子速度分布函数的简化控制方程。发展和应用离散速度坐标法于气体分子速度空间，利用一套在物理空间和时间上连续而速度空间离散的分布函数来代替原分布函数对速度空间的连续依赖性。基于非定常时间分裂数值计算方法和无波动、无自由参数的NND耗散差分格式，建立直接求解气体分子速度分布函数的气体运动论有限差分数值方法。推广应用改进的Gauss-Hermite无穷积分法和华罗庚－王元提出的以单和逼近重积分的黄金分割数论积分方法等，对离散速度空间进行宏观取矩获取物理空间各点的气体流动参数，由此发展一套从稀薄流到连续流各流域统一的气体运动论数值算法。通过对不同Knudsen数下一维激波管问题、二维圆柱绕流和三维球体绕流的初步数值实验表明文中发展的数值算法是可行的。     

基于广义卡尔曼滤波的桥梁结构物理参数识别

基于广义卡尔曼滤波提出了随机荷载作用下桥梁结构物理参数的识别方法。首先,以荷载为观测对象,推导出基于有限元模型的桥梁结构系统的观测方程,以结构待识别的物理参数为状态向量,建立系统状态方程;然后,对该状态方程和观测方程构成的非线性参数系统应用广义卡尔曼滤波,从而识别出结构的物理参数。对一座简支梁桥和一座三跨连续梁桥在不同工况下的物理参数识别进行了数值仿真,结果表明本文方法能够准确地识别桥梁结构全部刚度参数、质量参数和阻尼参数,且具有很强的抗噪性能,从而验证了本文方法的有效性和鲁棒性,可应用于识别大型桥梁结构的物理参数。     

The Kinematic Flow Structure for the Gvirtzman–Gorelick In Situ VOC Remediation System

Potential theory and Stokes' stream function techniques are used to investigate the flow structure around the recirculation system developed by Gvirtzman and Gorelick (1992, 1993), which consists of an extraction well and a gallery (trench) for the recharge of treated water to the aquifer. Analytical formulas are derived for the drawdown, velocity, and stream function for a model in which the extraction well is modeled as a uniformly distributed line sink and the gallery is modeled as a uniformly distributed ring source. Travel times are reported for water particles traveling along the streamlines containing 50 and 90% of the flow for various degrees of well penetration and various radii of the ring source. The travel times along the streamline resulting in the shortest travel time (not necessarily the shortest path) are also reported for various degrees of well penetration and various radii of the ring source. The method completely eliminates the use of numerical finite-difference or finite-element methods and can be used for optimization of technological parameters of this remediation system.     

Time series prediction of nonlinear ship structural responses in irregular seaways using a third-order Volterra model

To predict the nonlinear structural responses of a ship traveling through irregular waves, a third-order Volterra model was applied based on the given irregular data. A nonlinear wave–body interaction system was identified using the nonlinear autoregressive with exogenous input (NARX) technique, which is one of the most commonly used nonlinear system identification schemes. The harmonic probing method was applied to extract the first-, second- and third-order frequency response functions of the system. To achieve this, a given set of time history data of both the irregular wave excitation and the corresponding midship vertical bending moment for a certain sea state was fed into the three-layer perceptron neural network. The network parameters are determined based on the supervised training. Next, the harmonic probing method was applied to the identified system to extract the frequency response function of each order. While applying the harmonic probing method, the nonlinear activation function (i.e., the hyperbolic tangent function) was expanded into a Taylor series for harmonic component matching. After the frequency response functions were obtained, the structural responses of the ship under an arbitrary random wave excitation were easily calculated with rapidity using a third-order Volterra series. Additionally, the methodology was validated through the in-depth analysis of a nonlinear oscillator model for a weak quadratic and cubic stiffness term, whose analytic solutions are known. It was confirmed that the current method effectively predicts the nonlinear structural response of a large container carrier under arbitrary random wave excitation.     

A high‐order element‐based Galerkin method for the barotropic vorticity equation

A high‐order element‐based Galerkin method is developed to solve the non‐divergent barotropic vorticity equation (BVE). The solution process involves solving a conservative transport equation for the vorticity fields and a Poisson equation for the stream function fields. The discontinuous Galerkin method is employed for solving the transport equation and a spectral element method (continuous Galerkin) is used for the Poisson equation. A third‐order strong stability preserving explicit Runge–Kutta scheme is used for time integration. A series of tests have been performed to validate the model, which include the evolution of an idealized tropical cyclone and interaction of dual vortices in close proximity. The numerical convergence study is performed by solving the BVE on the sphere where the analytic solution is known. The test results are consistent with physical observations, and the model exhibits exponential convergence. Copyright © 2008 John Wiley & Sons, Ltd.     

Reduced order state‐space models from the pulse responses of a linearized CFD scheme

This paper describes a method for obtaining a time continuous reduced order model (ROM) from a system of time continuous linear differential equations. These equations are first put into a time discrete form using a finite difference approximation. The unit sample responses of the discrete system are calculated for each system input and these provide the Markov parameters of the system. An eigenvalue realization algorithm (ERA) is used to construct a discrete ROM. This ROM is then used to obtain a continuous ROM of the original continuous system. The focus of this paper is on the application of this method to the calculation of unsteady flows using the linearized Euler equations on moving meshes for aerofoils undergoing heave or linearized pitch motions. Applying a standard cell‐centre spatial discretization and taking account of mesh movement a continuous system of differential equations is obtained which are continuous in time. These are put into discrete time form using an implicit finite difference approximation. Results are presented demonstrating the efficiency of the system reduction method for this system. Copyright © 2003 John Wiley & Sons, Ltd.     

PARAMETER IDENTIFICATION IN OFFSHORE PLATFORM USING ARMA MODEL AND TECHNOLOGY OF EXTRACTING FREE VIBRATION SIGNAL

IntroductionOffshoreplatformisoneofthemostimportantresearchedstructuresofoceanengineeringin2 1stcentury .Forthemal_conditionandsufferedfromimpactingofwind ,wave,flowandiceandenvironmentalcorrosionfrequently ,thecontrolofvibrationandon_linemonitoringofstructurearebecomingmoreandmoreimportant.Parametersidentificationisoneofkerneltechnologieswhicharemadeofseriesmoderndesigntechniquesinoffshoreplatform[1].Thedynamicactionbetweenwind ,wave,flow ,iceandoffshoreplatformaretypicalnonstationarystochast…     

Time domain parameters identification of foundation-structure interaction system

The time domain parameter laenuncauon memoa oi me iounuauon-structure interaction system is presented. On the basis of building the computation mode and the motion equation of the foundation-structure interaction system, the system parameter identification method was established by using the extended Kalman filter （EKF） technique and taking the unknown parameters in the system as the augment state variables. And the time parameter identification process of the foundation-structure interaction system was implemented by using the data of the layer foundation-storehouse interaction system model test on the large vibration platform. The computation result shows that the established parameter identification method can induce good parameter estimation.     

半主动TLCD对偏心结构的离散模糊变结构控制

采用在结构水平双向设置TLCD半主动控制装置的方法,对偏心结构在多维地震作用下的离散模糊变结构控制问题进行了研究.首先建立了结构-TLCD扭转耦联控制系统微分方程及其离散化状态空间表达,然后阐述了基于模糊逻辑的变结构控制的基本策略,最后应用该模糊变结构控制策略对一个五层偏心结构的TLCD半主动控制问题进行了仿真分析.仿真分析的结果表明,该方法在对结构的平动反应和扭转反应都能起到较好的减震效果的同时,还具有对结构固有参数不确定性的较强的鲁棒性.     

Robust adaptive intelligent sliding model control for a class of uncertain chaotic systems with unknown time-delay

In this paper, a robust adaptive intelligent sliding model control (RAISMC) scheme for a class of uncertain chaotic systems with unknown time-delay is proposed. A sliding surface dynamic is appropriately constructed to guarantee the reachability of the specified sliding surface. Within this scheme, neuro-fuzzy network (NFN) is utilized to approximate the unknown continuous function. The robust controller is an adaptive controller used to dispel the unknown uncertainty and approximation errors. The adaptive parameters of the control system are tuned on-line by the derived adaptive laws based on a Lyapunov stability analysis. Using appropriate Lyapunov–Krasovskii (L–K) functional in the Lyapunov function candidate, the uncertainty caused by unknown time delay is compensated and the global asymptotic stability of the error dynamics system in the specified switching surface is accomplished. Finally, the proposed RAISMC system is applied to control a Hopfield neural network, Cellular neural networks, Rössler system, and to achieve synchronization between the Chen system with two time delays with Rössler system without time delay. The results are representative of outperformance of the proposed method in all cases.     

基于格子Boltzmann方法的粘弹性流体圆柱绕流分析

将光滑界面法引入到格子Boltzmann方法中分析粘弹性流体绕流问题,分别采用单松弛模型和对流扩散模型求解运动方程和Oldroyd-B本构方程,针对圆形和椭圆内部边界条件,给出连续界面插值函数,在此基础上,运用光滑界面法将内部边界转换为作用力项施加到演化方程中。首先分析圆柱绕流问题,给出不同材料参数情况下的流场分布和阻力系数计算结果,比较发现与宏观数值模拟结果相吻合。将模型拓展到绕椭圆流动中,分析椭圆形状和材料参数对粘弹性流体绕柱流的影响,发现随着椭圆长轴与短轴比值的增加和维森伯格数的增加,阻力系数逐渐下降,并且长短轴比对迭代收敛有较大影响。     

外共振耦合Duffing-van der Pol系统的首次穿越

研究了谐和力与宽带噪声激励下二自由度强非线性Duffing-van derPol系统的首次穿越问题. 在外共振情形, 应用随机平均法将系统动力学方程化为关于振幅与角变量的Itô随机微分方程. 然后建立了系统的可靠性函数满足的后向Kolmogorov方程以及平均首次穿越时间满足的Pontryagin方程. 在一定的边界条件和初始条件下, 用有限差分法求解了这两个高维偏微分方程, 得到系统的条件可靠性函数、平均首次穿越时间以及平均首次穿越时间的条件概率密度. 讨论了不同参数对系统可靠性以及平均首次穿越时间的影响. 用Monte Carlo数值模拟验证了理论方法的有效性.     

Experimental investigation and analytical description of a vibro-impact NES coupled to a single-degree-of-freedom linear oscillator harmonically forced

In this paper, the dynamics of a system composed of a harmonically forced single-degree-of-freedom linear oscillator coupled to a vibro-impact nonlinear energy sink (VI-NES) is experimentally investigated. The mass ratio between the VI-NES and the primary system is about \(1\%\). Depending on the external force’s amplitude and frequency, either a strongly modulated response (SMR) or a constant amplitude response (CAR) is observed. In both cases, an irreversible transfer of energy occurs from the linear oscillator toward the VI-NES: process known in the literature as passive targeted energy transfer. Furthermore, the problem is analytically studied by using the method of multiple scales. The obtained slow invariant manifold shows the existence of a stable and of an unstable branch of solutions, as well as of an energy threshold (a saddle-node bifurcation) for the solutions to appear. Subsequently, the fixed points of the problem are calculated. When a stable fixed point is reached, the system is naturally drawn to it and a CAR is established, whereas when no stable point is attained, the system exhibits a SMR regime. Finally, a good correlation between the experimental and the analytical results is presented.     

多点非均匀调制演变随机激励下结构地震的响应

针对大跨度结构在非均匀调制演变随机激励作用下,考虑行波效应时的非平衡随机地震响应问题,应用虚拟激励法进行了分析,由于虚拟激励法自动计及了参振振型的互相关项以及激励之间的互相关项,理论上是精确解,时变功率谱的计算采用精细逐步积分格式,使计算效率进一步得到提高。     

基于静力响应面的结构有限元模型修正方法

提出了基于静力响应面的结构有限元模型修正方法.运用响应面方法,将结构静力响应和结构参数之间复杂的隐式关系用显式函数近似表达出来;在此响应面模型(函数)基础上,通过优化计算对结构有限元模型参数进行修正.阐述了基于静力响应面的结构有限元模型修正方法的基本理论和一般实现过程.对两跨连续梁结构的静力模型修正数值算例分析结果表明:基于静力响应面的有限元模型修正方法可以减少结构有限元计算的次数、提高模型修正的优化效率,结构有限元模型修正结果具有可接受的精度.     

SELF-CORRECTED PID VIBRATION CONTROL BASED ON DUAL-FREEDOM DUAL-DRIVE INTELLIGENT CANTILEVER BEAM 1)

随着科技不断进步,智能结构的振动控制在航天航空、机械制造、车辆与船舶等领域得到了广泛应用。由于多输入多输出存在多样性和复杂性,严重威胁系统稳定性。为了解决这一问题,针对两输入单输出的双驱动智能悬臂梁系统提出一种自适应控制策略,首先基于压电线性本构方程,应用假设模态方法建立双驱动智能悬臂梁的力学模型,得到了基于闭环控制系统的状态方程,同时利用递推最小二乘法在线辨识系统参数设计比例积分微分(proportional--integral--derivative, PID)控制器实现自校正PID控制。通过数值仿真对比在有无PID 控制下两输入单输出双驱动智能悬臂梁系统的振动情况,分析自校正PID 控制的控制效果。通过实验验证自校正PID 控制对双输入单输出的双驱动智能悬臂梁系统的控制效果;再设置两组不同的单输入单输出自校正PID控制实验作对比。结果表明：自校正PID 控制方法可以较为有效地抑制智能悬臂梁的自由振动,相比单输入单输出的两组,两输入单输出自校正PID控制的效果更为明显和有效。