直升机旋翼／机身耦合系统的气弹响应分析（二）方程的求解

在文（一）得到的直升机旋翼系统运动方程的基础上,运用动态子结构的方法对旋翼系统和机身作为相对独立的部分进行分析,通过界布的力平衡和几何协调将子系统耦合系统整体系统,用分割－迭代法求解直至旋翼系统和机身的响应同时收敛到精度要求,并研制相应的计算程序,给出了工程算例。     

转子—非线性支承系统振动响应的优化计算

本文用一种新的优化方案计算装有非线性弹性支承－挤压油膜阻尼器的转子系统的振动响应。首先根据转子系统的结构特点，建立其无量纲形式的非线性运动微分方程；然后由微分方程构造－控制目标函数，最后对此目标函数进行优化计算，求得转子系统的振动响应。     

Combined resonance of low pressure cylinder-generator rotor system with bending-torsion coupling

A nonlinear model of a low pressure cylinder-generator rotor system is presented to study sub-synchronous resonance and combined resonance. Analytical results are obtained by an averaging method. Transition sets and bifurcation diagrams are obtained based on the singularity theory for the two-state variable system. The bifurcation characteristics are analyzed to provide a basis for the optimal design and fault diagnosis of the rotor system. Finally, the theoretical results are verified with the numerical results.     

Stabilization control of a bumblebee in hovering and forward flight

Our previous study shows that the hovering and forward flight of a bumblebee do not have inherent stability （passive stability）. But the bumblebees are observed to fly stably. Stabilization control must have been applied. In this study, we investigate the longitudinal stabilization control of the bumblebee. The method of computational fluid dynamics is used to compute the control derivatives and the techniques of eigenvalue and eigenvector analysis and modal decomposition are used for solving the equations of motion. Controllability analysis shows that at all flight speeds considered, although inherently unstable, the flight is controllable. By feedbacking the state variables, i.e. vertical and horizontal velocities, pitching rate and pitch angle （which can be measured by the sensory system of the insect）, to produce changes in stroke angle and angle of attack of the wings, the flight can be stabilized, explaining why the bumblebees can fly stably even if they are passively unstable.     

Nonlinear aeroelastic coupled trim and stability analysis of rotor-fuselage

Based on the Hamilton principle and the moderate deflection beam theory, discretizing the helicopter blade into a number of beam elements with 15 degrees of freedora, and using a quasi-steady aero-model, a nonlinear coupled rotor/fuselage equation is established. A periodic solution of blades and fuselage is obtained through aeroelastic coupled trim using the temporal finite element method （TEM）. The Peters dynamic inflow model is used for vehicle stability. A program for computation is developed, which produces the blade responses, hub loads, and rotor pitch controls. The correlation between the analytical results and related literature is good. The converged solution simultaneously satisfies the blade and the vehicle equilibrium equations.     

Bifurcation and chaotic response of a cracked rotor system with viscoelastic supports

Cracks appearing in the shaft of a rotary system are one of the main causes of accidents for large rotary machine systems. This research focuses on investigating the bifurcation and chaotic behavior of a rotating system with considerations of various crack depth and rotating speed of the system’s shaft. An equivalent linear-spring model is utilized to describe the cracks on the shaft. The breathing of the cracks due to the rotation of the shaft is represented with a series truncated time-varying cosine series. The geometric nonlinearity of the shaft, the masses of the shaft and a disc mounted on the shaft, and the viscoelasticity of the supports are taken into account in modeling the nonlinear dynamic rotor system. Numerical simulations are performed to study the bifurcation and chaos of the system. Effects of the shaft’s rotational speed, various crack depths and viscosity coefficients on the nonlinear dynamic properties of the system are investigated in detail. The system shows the existence of rich bifurcation and chaos characteristics with various system parameters. The results of this research may provide guidance for rotary machine design, machining on rotary machines, and monitoring or diagnosing of rotor system cracks.     

Modeling and analysis of a coupled rigid-flexible system

ＩｎｔｒｏｄｕｃｔｉｏｎＭａｎｙｍｅｃｈａｎｉｃａｌｓｙｓｔｅｍｓｓｕｃｈａｓｒｏｂｏｔｍａｎｉｐｕｌａｔｏｒ,ｔｕｒｂｉｎｅｂｌａｄｅｓ,ｈｅｌｉｃｏｐｔｅｒｐｒｏｐｅｌｌｅｒａｎｄｍａｎ＿ｍａｄｅｓａｔｅｌｌｉｔｅｗｉｔｈｆｌｅｘｉｂｌｅａｔｔａｃｈｍｅｎｔｓａｒｅｃｏｍｐｏｓｅｄｏｆｒｉｇｉｄａｎｄｆｌｅｘｉｂｌｅｂｏｄｉｅｓ．Ｔｈｅｙａｒｅｃａｌｌｅｄｃｏｕｐｌｅｄｒｉｇｉｄ＿ｆｌｅｘｉｂｌｅｓｙｓｔｅｍｓ．Ｃｏｒｒｅｃｔｐｒｅｄｉｃｔｉｏｎｓｏｆｔｈｅｂｅｈａｖｉｏｒｏｆｆｌｅｘｉｂｌ…     

Active control of coupled flexural-torsional vibration in a flexible rotor-bearing system using electromagnetic actuator

Rotor-shaft systems are subject to non-uniform spin speed during start-up, coast-down or any non-stationary situation changing the spin speed suddenly, e.g., load fluctuation or sudden mass-loss like loss of a blade or a part thereof. For a flexurally and torsionally compliant rotor-shaft, the dynamics under non-uniform spin-speed shows inertial coupling among transverse and torsional coordinates through mass-unbalance and gyroscopic effect. This results into coupled transverse-torsional vibration, where torsional response consists of significant harmonic components at bisynchronous spin frequency, torsional natural frequency of the shaft, and at combination frequencies corresponding to sum and difference of spin and transverse natural frequencies and twice the transverse natural frequency of the rotor-shaft. As a result of the coupling, transverse rotor motion also influences the torsional motion. The Method of Multiple Scales (MMS) is used in this work to carry out an analysis of a simplified system to get an idea about the dominant frequencies of excitation. Results of numerical simulation are presented next to show the effectiveness and influence of actively controlling the transverse rotor motion on its torsional motion, at the dominant frequencies, with the help of non-contact electromagnetic force from an actuator. Transverse vibration control is also observed to control the torsional oscillations due to coupled nature of the problem. The Stability Limit Speed (SLS) of the system is also increased as a result of application of the active control action. Constant axial torque is observed to diminish the influence of coupling, and protect the system against torsional instability, but control action is a must to stabilize the transverse vibration of the system above SLS.     

Numerical investigation of structural geometric nonlinearity effect in high-aspect-ratio wing using CFD/CSD coupled approach

An efficient and robust fluid–structure coupled methodology has been developed to investigate the linear and non-linear static aeroelastic behavior of flexible high-aspect-ratio wing. A three-dimensional open source finite element solver has been loosely coupled with an in-house Reynolds-averaged Navier–Stokes solver, designed for hybrid-unstructured meshes, to perform aero-structural coupled simulations. For volume mesh deformation and two-way data interpolation over non-matching grids interface, a radial basis function methodology combined with a data reduction algorithm has been used. This technique is efficient in handling large deflections and provides high-quality deformed meshes. Structural geometric nonlinearity has been considered to predict the deformations in the vertical and torsional directions caused by gravitational and aerodynamic loading. A multi-material finite element model has been generated to match the experimental configuration. Computational aeroelastic simulations were performed on an experimental high-aspect-ratio aeroelastic wing model with a slender body at the tip to get non-linear static deflections, twist and structure natural frequencies. The effect of the geometric nonlinearity is significant for large deformation analysis and has been highlighted in the predicted maximum tip deflection and twist. Good qualitative and quantitative agreement has been achieved between the predicted results and the available experimental data.     

非自治系统用瞬态响应求周期解稳定度的方法

根据Floquet理论关于线性周期系数系统解的性质及稳定性条件 ,定义了衡量非线性非自治系统周期解受扰后的衰减速率指标—稳定度。从动力系统流的概念出发 ,给出了利用非线性非自治系统稳态周期解受扰后的瞬态响应信息计算周期解稳定度的方法。以不平衡滑动轴承 弹性转子系统为例 ,说明了该方法的有效性。将稳定度等于零作为临界判据 ,该方法不仅解决了工频周期解失稳边界的确定问题 ,而且解决了渐进稳定域的估计和抗冲击扰动裕度的计算问题     

UKF在深组合GPS/INS导航系统中的应用

采用联邦滤波的深组合GPS/INS导航系统预滤波器量测模型具有很强的非线性,导致扩展卡尔曼滤波（EKF）的预滤波器估计精度不高。Unscented卡尔曼滤波（UKF）方法是一种非线性分布近似方法,它使用有限数量的sigma点去逼近整个非线性动态系统的分布可能,从而避免了对非线性测量模型进行线性化,具有较高的精度和较好的鲁棒性。在分析深组合导航系统预滤波器模型和UKF原理的基础上,设计了基于UKF滤波算法的预滤波器,对码相位误差、载波相位误差、载波频率误差、载波频率变化率等参数进行估计,同时将UKF和EKF算法进行了仿真比较。结果表明,在深组合导航系统中使用UKF滤波比EKF有更高的导航定位精度。     

Mechanical modelling and stability investigation of a non-material system in rotor dynamics

For the first time the behaviour of a Timoshenko-rotor-model with a non-material constraint is investigated. The constraint is caused by an axially-moving disc guided by the flexible shaft. Both, the development of the equations of motion (including the additionally occuring jump conditions) and the analysis of stability are essentially influenced by the non-classical character of the system. As result some stability diagrams are shown. They are based on statistical methods of theory of stability. The results allow the conclusion that most of the non-material constraints lead to a system behaviour as well-known from parametric excitations.     

不同轮轨接触模型在车桥耦合振动中的比较

以工程实例为研究对象,建立了整车-整桥系统耦合振动数值分析模型。考虑车轮的跳轨和挤密情况,建立了单边弹簧-阻尼系统弹性轮轨接触模型。采用基于多体系统动力学和有限元法结合的联合仿真技术,计算了两种轮轨接触时动车组列车以不同车速通过大跨度连续桥梁的耦合振动响应。数值计算结果表明:两种轮轨接触模型的桥梁动力响应比较接近;列车的横向轮轨力、轮重减载率和脱轨系数相差较大,当速度为350km/h时,横向轮轨力增大了46.5%,轮重减载率增大了130.8%,脱轨系数增大了24.66%;用单边-弹簧阻尼系统弹性轮轨接触模型更符合实际。     

Non-Linear response of a long,discontinuous fiber/melt system in elongational flows

The authors investigated the transient elongational behavior of a highly-aligned 600% volume fraction long, discontinuous fiber filled poly-ether-ketone-ketone melt with a computer-controlled extensional rheometer at 370°C. Prior experiments at controlled strain rate and stress produced _E ⁺ (t, ) and ^– (t, _E) similar to a shear dominated flow of a non-linear viscoelastic fluid. Stress relaxation following steady extension showed nonlinear effects in the change in stress decay rate with increasing strain rate. Continuous relaxation spectra showed a shift in the spectral peak to smaller values of with increasing strain rate. The Giesekus nonlinear constitutive relation modeled the elongation and stress relaxation with shearing rate at the fiber surface set by a strain rate magnification factor. Suitable for elongation, the model produced insufficient shift in the stress relaxation spectrum to account for the large change in stress decay rate exhibited in the experiments.English alphabet a _r aspect ratio of the fibers or l/d - A ₀ initial uniform cross-section area of the specimen - d fiber diameter - f fiber volume fraction - H() relaxation spectrum found by the method of Ferry and William l length of the fiber - L(t) time function specimen length - L ₀ initial specimen length - r radial coordinate across the shear cell - R _i fiber radius and inner cell dimension - R _o outer cell radius - t time in s - t _max duration of the extension - T _g glass transition temperature of the polymer - v velocity of the moving end of the test specimen - x axial position where is calculated Greek alphabet nonlinearity parameter in the Giesekus relation - axial mass distribution along the specimen major axis - shear strain rate - strain tensor - ₍₁₎ first convected derivative of the strain tensor - ₍₂₎ second convected derivative of the strain tensor - average strain at the end of extension as determined from - extension strain rate - average extension strain rate determined from - ^– transient strain rate under controlled stress, creep, test - _E elongational viscosity - _Eapp apparent elongational viscosity determined from - _E ⁺ transient elongational viscosity - ₀ zero shear rate viscosity - relaxation parameter - ₁ relaxation parameter in either Jeffrey's or Giesekus fluid - ₂ retardation parameter in either Jeffrey's or Giesekus fluid - _max relaxation value at which 99.9% of the H spectrum had occurred - _p relaxation value at which H reaches a maximum - volumetric composite density - _E elongational stress - _E ⁺ transient elongational stress - _E ^– controlled elongational stress, creep stress - _E ^y peak elongational stress in controlled experiment - shear stress at surface of the fiber in a shear cell - _yx simple shear component of the strain rate tensor - stress tensor - ₁ first convected derivative of the stress tensor     

变质量、变阻尼、变刚度结构系统的动力响应

本文首次提出了变质量、变阻尼、变刚度结构系统的动力响应问题,并根据抓斗卸船机的工作过程导出了结构系统的动力方程及解法。最后给出了工程实际算例。     

Fractional order spring/spring-pot/actuator element in a multibody system: Application of an expansion formula

Fractional order models of a spring/spring-pot and spring/spring-pot/actuator element connected into a multibody system are proposed in order to represent smart materials and components in adaptronic systems by introducing new tuning parameter. The models are introduced into dynamic equations via generalized forces and using the Lagrange's equations of the second kind in covariant form. Generalized forces are derived by taking into account fractional order derivatives in force–displacement relations and by using the principle of virtual work. The numerical scheme for solving fractional order differential equations proposed in Atanacković and Stanković (2008) is used in order to approximate fractional order derivative of a composite function appearing in the presented fractional order model. Numerical example for the multibody system with three degrees of freedom is presented. The results obtained for generalized forces are compared for different values of parameters in the fractional order derivative model.     

Nonlinear normal modes and their superposition in a two degrees of freedom asymmetric system with cubic nonlinearities

This paper investigates nonlinear normal modes and their superposition in a two degrees of freedom asymmetric system with cubic nonlinearities for all nonsingular conditions, based on the invariant subspace in nonlinear normal modes for the nonlinear equations of motion. The focus of attention is to consider relation between the validity of superposition and the static bifurcation of modal dynamics. The numerical results show that the validity has something to do not only with its local restriction, but also with the static bifurcation of modal dynamics. Project Supported by the National Natural Science Foundation and PSF of China     

Stability of equilibrium states in a simple system with unilateral contact and Coulomb friction

The aim of this paper is to study the stability of equilibrium states in a mechanical system involving unilateral contact with Coulomb friction. Since the assumptions made in classical stability theorems are not satisfied with this class of systems, we return to the basic definitions of stability by studying the time evolution of the distance between a given equilibrium and the solution of a Cauchy problem where the initial conditions are in a neighborhood of the equilibrium. It was recently established that the dynamics is well posed in the case of analytical data. In the present study, we focus in particular on the stability of the equilibrium states under a constant force and deal only with a simple mass-spring system in .     

圆柱坐标系中弹性力学的哈密顿正则方程及其哈密顿单元方法

给出圆柱坐标系中弹性力学的哈密顿形式体系，用有限元法得到了哈密顿单元的半解析解，可用于求解复合材料叠层圆柱曲板问题。