Quantification of non-viscous damping in discrete linear systems

The damping forces in a multiple-degree-of-freedom engineering dynamic system may not be accurately described by the familiar ‘viscous damping model’. The purpose of this paper is to develop indices to quantify the extent of any departures from this model, in other words the amount of ‘non-viscosity’ of damping in discrete linear systems. Four indices are proposed. Two of these indices are based on the non-viscous damping matrix of the system. A third index is based on the residue matrices of the system transfer functions and the fourth is based on the (measured) complex modes of the system. The performance of the proposed indices is examined by considering numerical examples.     

Optimum distribution of additive damping for vibrating beams

Cost and weight effectiveness of concentrated and distributed additive damping is studied for linear systems (discrete and continuous) under prescribed harmonic loads and/or displacements. Stiffness and mass changes due to additive damping are included. From a numerical example studied it can be concluded that optimal damping distributions can reduce resonant responses by about 40% as compared to uniformly distributed damping of the same cost or weight. The optimization technique as well as an exact displacement method for analysis of harmonically vibrating beams and frames are presented.     

Use of the Hilbert transform in modal analysis of linear and non-linear structures

An initial study into the application of the Hilbert transform in modal analysis procedures is presented. It is shown that typical structural non-linearities such as non-linear damping and stiffness can be detected and identified directly without the need to generate explicit models. No assumptions regarding the degree of non-linearity are made, which is a restriction in the classical methods for dealing with non-linearities. The properties of the Hilbert transform are discussed with respect to linear and non-linear dynamical systems, and a discrete transform, developed from the continuous functions, is derived in the frequency domain and adapted to modal analysis data in the form of mobility transfer functions. Truncation effects arising from limited frequency ranges of the mobility transfer functions are accounted for by employing correction terms in the frequency domain. Several examples are studied of single and multi-mode systems with non-linearities such as friction, clearance and non-linear stiffness. These examples indicate that the Hilbert transform offers a new method for extending modal analysis to the domains of non-linear systems.     

Influences of Structure Disorder and Temperature on Properties of Proton Conductivity in Hydrogen-Bond Molecular Systems

The dynamic properties of proton conductivity along hydrogen-bonded molecular systems,for example,ice crystal,with structure disorder or damping and finite temperatures exposed in an externally applied electric-field have been numerically studied by Runge-Kutta way in our Soliton model.The results obtained show that the proton-soliton is very robust against the structure disorder including the fluctuation of the force constant and disorder in the sequence of masses and thermal perturbation and damping of medium,the velocity of its conductivity increases with increasing of the externally applied electric-field and decreasing of the damping coefficient of medium,but the proton-soliton disperses for quite great fluctuation of the force constant and damping coefficient.In the numerical simulation we find that the proton-soliton in our model is thermally stable in a large region of temperature of T ≤ 273 K under influences of damping and externally applied electric-field in ice crystal.This shows that our model is available and appropriate to ice.     

A model for vibration transmission of light-heavy structures

The vibration transmission of light-heavy structures is investigated in this paper. The light-heavy structure consists of a thin beam and a mass block. Based on numerical simulations with the finite element method and experiments, the block's effect on the thin beam is defined. A theoretical model for this beam-block structure is successfully developed, which is validated and agrees very well with the numerical and experimental models. Two kinds of transfer functions of velocities between any two points on the beam-block structure are studied experimentally and theoretically. The theoretical transfer functions agree well with the experimental results. There are peaks and valleys in the transfer functions, where the peaks occur at the anti-resonant frequencies of the second point and the valleys at the anti-resonant frequencies of the first point. Away from these peaks and valleys the magnitude of the transfer functions are about 0 dB for two points on the beam, and about 20 dB in our experiments for a point on the beam and another point on the block (close to the theoretical prediction of 18 dB determined by the mass ratio of the beam and the block). With these transfer functions, new techniques might be developed for indirect measurement of the vibration of the thin beam by measuring the vibration of the block.     

Novel schemes for inserting seismic dampers in shear-type systems based upon the mass proportional component of the Rayleigh damping matrix

This paper presents the results of an extensive study carried out to investigate the applicability of a novel scheme for inserting added viscous dampers in shear-type systems. The findings, even though developed with specific reference to civil building structures, provide useful insight also for the effective addition of viscous dampers in mechanical dynamic systems (of similar characteristics) when excited at the base.The novel scheme proposed (referred to as the MPD system) is based upon the mass proportional component of the Rayleigh damping matrix (MPD matrix) and is characterised by a peculiar damper placement which sees the dampers placed so that they connect each mass to a fixed point.Firstly, the paper briefly recalls (a) the physical principles and (b) selected results of numerical investigations which show that the MPD system is characterised by superior dissipative properties.Secondly, the paper investigates the implementation of the MPD system in civil building structures. Two solutions are envisaged herein: direct implementation (through the use of long buckling-resistant dampers which connect each storey to the ground) and indirect implementation (by placing common dampers between the structure and a very stiff lateral-resisting element adjacent or internal to the structure). The first solution leads to the implementation in the structure of an exact MPD matrix, if damper sizing is chosen appropriately. The second solution (simpler than the first one to implement in building structures) leads to an exact MPD matrix, if, in addition to appropriate damper sizing, the lateral-resisting element is infinitely stiff. As far as the direct implementation is concerned, this paper shows how long buckling-resistant braces are available for structural systems up to three storey high. As far as the indirect implementation is concerned, this paper shows (through extensive numerical parametric investigations) how this solution is capable of providing damping effects which are similar to those offered by the direct implementation, even for lateral-resisting elements characterised by finite lateral stiffness. The results obtained also provide insight for the optimal insertion of viscous dampers in coupled mechanical dynamic systems.     

核磁共振中的量子控制

近年来, 随着量子信息科学的发展, 对由量子力学原理描述的微观世界的主动调控已成为重要的前沿研究领域. 为构造实际的量子信息处理器, 一个关键的挑战是: 如何对处于噪声环境下的量子体系实现一系列高精度的任意操作, 以完成目标量子信息处理任务. 为此, 人们将经典系统控制论的思想方法延伸到量子体系的领域, 提出了大量的量子控制方法以及相关的数值技术(如量子优化控制、量子反馈控制等), 并取得了丰富的研究成果. 核磁共振自旋体系具备成熟的系统理论和操控技术, 为量子控制方法的实用性研究提供了优秀的实验测试平台. 因此, 基于核磁共振的量子控制成为量子控制领域的重要方向. 本文简要介绍了量子控制的基本概念和方法; 从系统控制论的角度对核磁共振自旋体系的基本原理和重要控制任务做了阐述; 介绍了近些年来在该领域发展的相关控制方法及其应用; 对基于核磁共振体系的量子控制的进一步的研究做了几点展望.     

On eigenproblem solution of damped vibrations associated with gyroscopic moments

A new efficient approach is presented for solving the quadratic eigenvalue problem of weakly, nonproportionally damped vibration systems. In the analysis of these systems, gyroscopic moments and external damping are both considered. Traditional restriction of symmetry of inertia, damping and stiffness matrices is slightly relaxed. A second-order perturbation theory is developed such that the perturbed solution is based on the eigensolution of an unperturbed subproblem. This subproblem considers the unperturbed system in two different forms: (i) a conservative, gyroscopic part of an original problem, or (ii) a nonconservative, gyroscopic part of an original problem that is proportionally damped. To cope with asymmetry of the system matrices, a Duncan's like state formulation is used to bring these matrices into a suitable form for perturbations. Two numerical examples are introduced for explaining the detailed implementation of the presented approach. Additionally, a practical problem of rotor supported by two tilting pad-bearings is investigated. The eigensolutions obtained by the current approach match, to a great extent, other solutions obtained by time-consuming exact methods. The investigation procedure given here gives a framework to handle vibration problems of weakly nonproportional damping and/or weakly asymmetric inertia, damping and stiffness matrices.     

Transverse emittance blow-up from beam injection errors in synchrotrons with nonlinear feedback systems

The problem of transverse emittance blow-up from beam injection errors in synchrotrons with nonlinear feedback systems is considered. The relative emittance growth is calculated for linear and nonlinear feedback transfer functions. Effects of an increase of the damping decrement of the beam coherent oscillations and of a decrease of the coherent transverse amplitude spread of different bunches in case of the damper with a positive cubic term in the feedback transfer function are discussed. The text was submitted by the author in English.     

Damping uncertainty due to noise and exponential windowing

Estimation of damping levels in structures is required in many applications and the theoretical damping predictions in many cases are either difficult or not very reliable. Therefore, determination of damping levels based on experimental data is employed quite often. However, in some cases, experimental approach may not be able to provide definite answers either in the sense that the identified damping level may exhibit high level of uncertainty. In experimental approach, the most widely used methods for damping determination require vibration spectrums or the Frequency Response Function(s) which are obtained by Fourier transformation of the time-domain data. During this process, it is often necessary, especially for lightly damped structures, to modify the time-domain data by using exponential windowing so as to minimise the leakage effect in the spectrum. The so-called numerical damping artificially added by this process can be subtracted later in order to obtain the correct damping level. However, for lightly damped structures, the artificially introduced numerical damping can be significantly greater than the actual damping level. This inevitably brings the accuracy and the reliability issues, especially when the data are contaminated by noise. This paper addresses damping uncertainty in frequency-domain estimation when (i) the data are contaminated by noise and (ii) numerical damping via exponential windowing is introduced during the signal processing phase of the spectrum estimation. Some numerical simulations are performed first in order to assess the adverse effects of noise on damping estimations and resulting damping uncertainty is examined as a function of noise level in the data. Then, damping uncertainty due to the use of exponential windowing is investigated using experimental data. A relationship between damping uncertainty and the level of added numerical damping is presented when the so-called Line-Fit method is used for damping estimation from measured Frequency Response Functions.     

Practical improvements to real and imaginary spectral based modal parameter measurements of SDOF systems

For a single degree of freedom system, especially with non-light damping, the use of the real spectral part of either the displacement or velocity responses (or the transfer functions based on them) has the advantages of determining the natural frequency (f_n) directly, independent of the response parameter, and providing an accurate measurement for damping (ζ). However, this method is sensitive to spectral phase errors due to temporal misalignment of the signal or due to net inter-channel delay differences caused by signal filtering. Two techniques are presented to correct for these errors; one based on the correct temporal alignment of the impulsive part, and the other on the infimum of the imaginary spectral part. These are first demonstrated using a numerical model, and are shown to facilitate the correct measurement of f_n and bring ζ within the expected error limits due to quantisation. Secondly, they are applied to an experimental system and are seen to greatly improve the consistency between measurements using different methods.     

Damping lateral vibrations in rotary machinery using motor speed modulation

This paper outlines a new principle for damping lateral vibrations of rotary systems. According to this principle, no changes in the visco-elastic properties of the system to be damped are required. The method is based on the generation of a harmonic additive to the constant speed of rotation that provides significant damping of lateral vibrations at critical speeds of rotation. This concept is validated analytically using the method of averaging and additionally with the help of direct numerical integration. The solution is shown to represent Fourier series containing Bessel functions. Consequently, proper choice of the parameters of the additional harmonic component ensuring that the Bessel functions have minimum values results from a minimization of the solution itself. Thus, the analytical solution and numerical results prove this concept by showing an essential decrease of the amplitudes of lateral vibrations of the damped system compared with those of the undamped system. The physical explanation of this effect is presented.     

Modal analysis by free vibration response only for discrete and continuous systems

This work aims to develop the algorithm for modal analysis by free vibration response only (MAFVRO), in particular for the general or non-proportional viscous damping system model. If the structural displacement or acceleration response due to free vibration can be measured, the system response matrices, including the displacement, velocity and acceleration, can be obtained through numerical differential or integration methods. These response matrices can then be applied to the developed MAFVRO method to determine the structural modal parameters. The numerical differential and integration methods are introduced and adopted to establish the modal parameter prediction program for the non-proportional damping model of MAFVRO. This work also shows the applications of MAFVRO to the multiple degree-of-freedom (mdof) systems and the cantilever beam, respectively. Both the discrete and continuous systems are demonstrated for the feasibility of the MAFVRO algorithm. The developed method uses the free vibration output response only and can obtain the structural modal parameters successfully.     

Response and stability of a SDOF strongly nonlinear stochastic system with light damping modeled by a fractional derivative

A stochastic averaging procedure for a single-degree-of-freedom (SDOF) strongly nonlinear system with light damping modeled by a fractional derivative under Gaussian white noise excitations is developed by using the so-called generalized harmonic functions. The approximate stationary probability density and the largest Lyapunov exponent of the system are obtained from the averaged Itô stochastic differential equation of the system. It is shown that the approximate stationary solutions obtained by using the stochastic averaging procedure agree well with those from the numerical simulation of original systems. The effects of system parameters on the approxiamte stationary probability density and the largest Lyapunov exponent of the system are also discussed.     

Stochastic weighted particle methods for population balance equations

A class of coagulation weight transfer functions is constructed, each member of which leads to a stochastic particle algorithm for the numerical treatment of population balance equations. These algorithms are based on systems of weighted computational particles and the weight transfer functions are constructed such that the number of computational particles does not change during coagulation events. The algorithms also facilitate the simulation of physical processes that change single particles, such as growth, or other surface reactions.     

Time-dependent variance and covariance of responses of structures to non-stationary random excitations

In this paper techniques for the analysis of non-stationary random responses of linear structures, discretized by the finite element method so that they can be analyzed as multi-degree of freedom systems, subjected to non-stationary random excitation are developed. The non-stationary random excitation is represented as a product of (a) an exponentially decaying function and a white noise process, and (b) a modulating function in the form of an exponential envelope and a white noise process. Closed form expressions for the time-dependent variance and covariance of response of structures are presented. Application of these expressions is made for the analysis of non-stationary random responses of a physical model of a class of mast antenna structures subjected to base excitation. It is concluded that (a) the coupling terms do have a definite influence on the response; the magnitude of the influence is proportional to the amount of damping in the structure and proximity of the modes excited; (b) the non-stationary random excitations considered are general in that the modulating functions are not necessarily identical, and therefore the influence of various modulating functions of the excitations applied to different locations of the structure on responses can be examined quantivatively; and (c) for a given damping parameter the magnitudes of the modulating function parameters cannot be chosen arbitrarily though the shapes of normalized modulating functions can be selected to best fit the excitation realizations.     

Stationary response of combined linear dynamic systems to stationary excitation

The stationary response of a broad class of combined linear systems to stationary random excitation is determined by the normal mode method. The systems are characterized by a viscously damped simple beam (or string, membrane, thin plate or shell, etc.) connected at discrete points to a multiplicity of viscously damped linear oscillators and/or masses. The solution of the free vibration problem by way of Green functions and the deterministic forced vibration problem by modal analysis for both proportional and non-proportional damping is reviewed. The orthogonality relation for the natural modes of vibration is used to derive a unique relationship between the cross-spectral density functions of the applied forces and the cross-spectral density functions of the generalized forces. Finally, the response spectral density functions and the mean square responses of the beam and oscillators are derived in closed form, exact for the proportionally damped system and approximate for the non-proportionally damped system.     

A new analytical technique for vibration analysis of non-proportionally damped beams

Vibrating linear mechanical systems, in particular continuous systems, are often modelled considering proportional damping distributions only, although in many real situations this simplified approach does not describe the dynamics of the system with sufficient accuracy. In this paper an analytical method is given to take into account the effects of a more general viscous damping model, referred to as non-proportional damping, on a class of vibrating continuous systems. A state-form expansion applied in conjunction with a transfer matrix technique is adopted to extract the eigenvalues and to express the eigenfunctions in analytical form, i.e., complex modes corresponding to non-synchronous motions. Numerical examples are included in order to show the efficiency of the proposed method; non-proportional damping distributions of different type, such as internal and external lumped or distributed viscous damping, are tested on non-homogeneous Euler-Bernoulli beams in bending vibration with different boundary conditions. Finally, a discussion on root locus diagrams behaviour and on modal damping ratio significance for non-proportionally damped systems is presented.     

一类相对转动时滞非线性动力系统的稳定性分析

建立了具有时变刚度、非线性阻尼和谐波激励的一类相对转动时滞非线性动力系统的动力学方程.采用多尺度法推导出时滞动力系统的分岔响应方程,运用奇异性理论研究系统结构稳定性,得到主共振稳态响应方程的转迁集以及不同参数下分岔曲线的拓扑结构.应用Hopf分岔理论讨论了时滞动力系统动态稳定性,给出了系统产生极限环的条件,最后用数值模拟的方法研究了时滞参数对系统极限环幅值的影响。