Multi-scaled analysis of the damped dynamics of an elastic rod with an essentially nonlinear end attachment

We study multi-frequency transitions in the transient dynamics of a viscously damped dispersive finite rod with an essentially nonlinear end attachment. The attachment consists of a small mass connected to the rod by means of an essentially nonlinear stiffness in parallel to a viscous damper. First, the periodic orbits of the underlying hamiltonian system with no damping are computed, and depicted in a frequency–energy plot (FEP). This representation enables one to clearly distinguish between the different types of periodic motions, forming back bone curves and subharmonic tongues. Then the damped dynamics of the system is computed; the rod and attachment responses are initially analyzed by the numerical Morlet wavelet transform (WT), and then by the empirical mode decomposition (EMD) or Hilbert–Huang transform (HTT), whereby, the time series are decomposed in terms of intrinsic mode functions (IMFs) at different characteristic time scales (or, equivalently, frequency scales). Comparisons of the evolutions of the instantaneous frequencies of the IMFs to the WT spectra of the time series enables one to identify the dominant IMFs of the signals, as well as, the time scales at which the dominant dynamics evolve at different time windows of the responses; hence, it is possible to reconstruct complex transient responses as superposition of the dominant IMFs involving different time scales of the dynamical response. Moreover, by superimposing the WT spectra and the instantaneous frequencies of the IMFs to the FEPs of the underlying hamiltonian system, one is able to clearly identify the multi-scaled transitions that occur in the transient damped dynamics, and to interpret them as ‘jumps’ between different branches of periodic orbits of the underlying hamiltonian system. As a result, this work develops a physics-based, multi-scaled framework and provides the necessary computational tools for multi-scaled analysis of complex multi-frequency transitions of essentially nonlinear dynamical systems.     

Broadband energy exchanges between a dissipative elastic rod and a multi-degree-of-freedom dissipative essentially non-linear attachment

We analyze complex, multi-frequency, non-linear modal interactions in the damped dynamics of a viscously damped dispersive finite rod coupled to a multi-degree-of-freedom essentially non-linear attachment. We perform a parametric study to show that the attachment can be an effective broadband energy absorber and dissipater of shock energy from the rod. It is shown that strong targeted energy transfer from the rod to the attachment occurs when there is strong stiffness asymmetry in the attachment. For weak viscous dissipation, a clear understanding of dynamical transitions in the integrated rod-non-linear attachment system can be gained by wavelet transforming the time series and superimposing the resulting wavelet spectra in the frequency-energy plot (FEP) of the periodic orbits of the underlying Hamiltonian system. Two distinct NES configurations are analyzed in detail, and their damped responses are analyzed by the Hilbert-Huang transform (HHT). We show that the HHT is capable of analyzing even complex non-linear damped transitions, by providing the dominant frequency components (or equivalently, time scales) at which the non-linear phenomena take place, and clarifying the series of non-linear resonance captures between the rod and attachment dynamics that are responsible for the broadband energy exchanges in this system.     

含摩擦与碰撞平面多刚体系统动力学线性互补算法

基于接触力学理论和线性互补问题的算法, 给出了一种含接触、碰撞以及库伦干摩擦, 同时具有理想定常约束(铰链约束) 和非定常约束(驱动约束) 的平面多刚体系统动力学的建模与数值计算方法. 将系统中的每个物体视为刚体, 但考虑物体接触点的局部变形, 将物体间的法向接触力表示成嵌入量与嵌入速度的非线性函数,其切向摩擦力采用库伦干摩擦模型. 利用摩擦余量和接触点的切向加速度等概念, 给出了摩擦定律的互补关系式; 并利用事件驱动法, 将接触点的黏滞-滑移状态切换的判断及黏滞状态下摩擦力的计算问题转化成线性互补问题的求解. 利用第一类拉格朗日方程和鲍姆加藤约束稳定化方法建立了系统的动力学方程, 由此可降低约束的漂移, 并可求解该系统的运动、法向接触力和切向摩擦力, 还可以求解理想铰链约束力和驱动约束力. 最后以一个类似夯机的平面多刚体系统为例, 分析了其动力学特性, 并说明了相关算法的有效性.      

Influence of friction and asymmetric freeplay on the limit cycle oscillation in aeroelastic system: An extended Hénon’s technique to temporal integration

Nonlinear effects such as friction and freeplay on the control surfaces can affect aeroelastic dynamics during flight. In particular, these nonlinearities can induce limit cycle oscillations (LCO), changing the system stability, and because of this it is essential to employ computational methods to predict this type of motion during the aircraft development cycle. In this context, the present article presents a matrix notation for describing the Hénon’s method used to reduce errors when considering piecewise linear nonlinearities in the numerical integration process. In addition, a new coordinate system is used to write the aeroelastic system of equations. The proposal defines a displacement vector with generalized and physical variables to simplify the computational implementation of the Hénon’s technique. Additionally, the article discusses the influence of asymmetric freeplay and friction on the LCO of an airfoil with control surface. The results show that the extended Hénon’s technique provides more accurate LCO predictions, that friction can change the frequency and amplitude of these motions, and the asymmetry of freeplay is important to determine the LCO behavior.     

Investigation of surface texture influence on hydrodynamic performance of parallel slider bearing under transient condition

Although surface texture has been becoming an attractive technique for improving hydrodynamic performance in lubrication system, the study of surface texture mainly focus on the steady-state load conditions. This investigation evaluates the effect of surface texture on hydrodynamic performance under transient lubrication in parallel bearing. In the case, considering cavitation phenomena in liquid lubricating film, a mass conserving formulation based on Elrod–Adams model with JFO conditions is then employed. The implementation of numerical simulation is discretized in spatial and time domain for calculating the pressure distribution. The effect of different geometrical configurations of lubricated contact on hydrodynamic performance is evaluated. The results show that surface texture in parallel bearing yields a great improvement in terms of friction force and friction coefficient.     

Piercing a liquid surface with an elastic rod: Buckling under capillary forces

When a thin elastic structure comes in contact with a liquid interface, capillary forces can be large enough to induce elastic deformations. This effect becomes particularly relevant at small scales where capillary forces are predominant, for example in microsystems (micro-electro-mechanical systems or microfluidic devices) under humid environments. In order to explore the interaction between capillarity and elasticity, we have developed a macroscopic model system in which an initially immersed vertical elastic rod is raised through a horizontal liquid surface. We follow a combined approach of experiments, theory and numerical simulations to study this system. In spite of its apparent simplicity, our experiment reveals a complex phase diagram, involving large hysteretic behaviour. We employ Kirchhoff equations for thin elastic rods and use path-following methods from which we obtain a variety of equilibrium states and associated transitions that are in excellent qualitative and quantitative agreement with those observed experimentally.     

The follower load problem and bifurcation to periodic motions

The problem of a follower load applied to a thin cantilevered elastic rod is presented. The stability of the straight static solution is considered, and a post-critical analysis is performed to describe non-trivial periodic motions. As the magnitude of the follower load is increased past a critical value, the real valued frequencies of the linearized system coalesce and then split apart with complex-conjugate values, rendering the straight state unstable. A Liapunov-Schmidt reduction is used to analyze the non-linear governing equations and to investigate the possibility of non-trivial periodic solutions. The rod is assumed to be incompressible and unshearable while only general requirements are made on the constitutive rule for bending. It is found that in the post-critical regime the higher order terms in the constitutive rule play an important role in determining whether non-trivial periodic solutions are possible.     

ANALYTICAL SOLUTIONS FOR RESPONSE OF COLLISION OF PARTICLE WITH CONICAL ROD CAUSED BY LONGITUDINAL VIBRATION

The objective is to present exact analytical solutions of longitudinal impact analysis for slender conical rods struck by a particle and a new method is proposed for conical rod-particle impact analysis, in which the superposition method is used and the response of the rod is presented. These analytical results are exact and can be used to validate the numerical methods or other analytical results. The numerical example shows that one of the advantages of the present method is that the analytical form is very simple. The result is that mass ratio and some variables describing the geometrical shape of rods such as taper, length and radius play an important role in impact dynamic system.     

Sparse matrix implicit numerical integration of the Stiff differential/algebraic equations: Implementation

The finite element absolute nodal coordinate formulation (ANCF) is often used in modeling very flexible bodies in multibody system (MBS) applications. This formulation leads to a constant mass matrix, allowing for an efficient sparse matrix implementation. Nonetheless, the use of the ANCF finite elements to model stiff structures can lead to high frequencies associated with ANCF coupled deformation modes, as discussed in the literature. Implicit numerical integration methods can be effectively used to develop efficient procedures for the solution of MBS differential/algebraic equations. Most existing implicit integration algorithms, however, require numerical differentiation of the equations of motion, and some of these integration methods do not ensure that the kinematic algebraic constraint equations are satisfied at all levels (position, velocity, and acceleration). Because of these limitations, existing implicit integration methods can be less accurate and less efficient when used to solve large scale MBS applications. In order to circumvent this problem, the two-loop implicit sparse matrix numerical integration (TLISMNI) method was proposed for the solution of MBS differential/algebraic equations. The TLISMNI method does not require numerical differentiation of the forces and allows for an efficient sparse matrix implementation. This paper discusses TLISMNI implementation issues including the step size selection, the error control, and the effect of the numerical damping. The relation between the step size selection and the structure stiffness is also discussed. The use of the computer implementation described in this paper is demonstrated by solving very stiff structure problems using the Hilber?CHughes?CTaylor (HHT) method, which includes numerical damping. An eigenvalue analysis and Fast Fourier Transform (FFT) are performed in order to identify the fundamental modes of deformation and demonstrate that the contributions of these fundamental modes can be erroneously damped out when some other implicit integration methods are used. The TLISMNI method, on the other hand, captures the contributions of these fundamental modes. The results, obtained using the TLISMNI method, are compared with the results obtained using other methods including the implicit HHT-I3 and the explicit Adams integration methods. The results obtained show that the TLISMNI method can be five times faster than the other two methods when no numerical damping is considered.     

Control of disturbances of a hypersonic viscous shock layer on a flat plate

Receptivity of a viscous shock layer on a flat plate aligned at an angle of attack to external multiwave acoustic perturbations is studied. It is shown that external acoustic waves and periodic controlled perturbations introduced from the surface of the plate mounted at an angle of attack smaller than 20° generate entropy-vortex disturbances with a similar spatial distribution in the viscous shock layer. This result allows numerical implementation of the interference method of controlling disturbances generated in the viscous shock layer on the plate by external acoustic waves at one frequency and at a spectrum of frequencies by introducing blowing-suction perturbations on the plate surface with appropriate amplitudes and phases.     

Wave-based defect detection and interwire friction modeling for overhead transmission lines

In this study, the feasibility of continuous, online monitoring of power lines using ultrasonic waves is considered. Local and global wave-based approaches for wire break detection in overhead transmission lines are presented. Both methods use a sending/receiving transducer to generate an ultrasonic, longitudinal, elastic wave in the cable. Defects in the cable cause a portion of the incident ultrasonic wave to be reflected back to the transducer, which when received, can be used to identify the presence of the defect. Although the transducers can only be attached to the surface of the cable, subsurface wires can also be interrogated since elastic energy spreads to these wires through friction contact. This study also explores how the elastic energy of a propagating wave becomes distributed among contacting rods via friction contact. This work focuses specifically on a two-rod system in which the wave energy from an excited “active” rod is transmitted to a neighboring “passive” rod through friction contact. An energy-based model is used to approximate the time average elastic wave power in the two rods as a function of propagation distance. Power predictions from the energy-based model compare well with experimental measurements and finite element simulations.     

Modeling and analysis of rigid multibody systems with driving constraints and frictional translation joints

An approach is proposed for modeling and anal- yses of rigid multibody systems with frictional translation joints and driving constraints. The geometric constraints of translational joints with small clearance are treated as bilat- eral constraints by neglecting the impact between sliders and guides. Firstly, the normal forces acting on sliders, the driv- ing constraint forces （or moments） and the constraint forces of smooth revolute joints are all described by complementary conditions. The frictional contacts are characterized by a set- valued force law of Coulomb＇s dry friction. Combined with the theory of the horizontal linear complementarity problem （HLCP）, an event-driven scheme is used to detect the transi- tions of the contact situation between sliders and guides, and the stick-slip transitions of sliders, respectively. And then, all constraint forces in the system can be computed easily. Secondly, the dynamic equations of multibody systems are written at the acceleration-force level by the Lagrange multiplier technique, and the Baumgarte stabilization method is used to reduce the constraint drift. Finally, a numerical example is given to show some non-smooth dynamical behaviors of the studied system. The obtained results validate the feasibility of algorithm and the effect of constraint stabilization.     

Subharmonic and grazing bifurcations for a simple bilinear oscillator

In this paper, subharmonic and grazing bifurcations for a simple bilinear oscillator, namely the limit discontinuous case of the smooth and discontinuous (SD) oscillator are studied. This system is an important model that can be used to investigate the transition from smooth to discontinuous dynamics. A combination of analytical and numerical methods is used to investigate the existence, stability and bifurcations of symmetric and asymmetric subharmonic orbits. Grazing bifurcations for a particular periodic orbit are also discussed and numerical results suggest that the bifurcations are discontinuous. We show via concrete numerical experiments that the dynamics of the system for the case of large dissipation is quite different from that for the case of small dissipation.     

Transient responses of repeated impact of a beam against a stop

The transient behavior of a cantilever beam, driven by periodic force and repeated impacting against a rod-like stop, is the subject of this investigation. As impact and separation phase take place alternately, the transient waves induced either by impacts or by separations will travel in more complicated ways. Thus the transient responses of both the beam and the rod during repeated impact become an important issue. In both impact phase and separation phase, the transient wave propagations are solved by the expansion of transient wave functions in a series of Eigenfunctions (wave modes). From the solutions, the answer of impact force is derived directly, so that the divergence problem, encountered in solving impact force numerically by a strongly non-linear equation coupled the unknown impact force with motions, has been avoided. Numerical results show the convergence of the time-step size and truncation number of wave modes in the calculations of impact force by the present method. As the transient wave effect is considered, the numerical results can show several transient phenomena involving the propagation of transient impact-induced waves, sub-impact phases, long-term impact motion, chatter, sticking motion, synchronous impact, non-synchronous impact (including asynchronous impact) and impact loss.     

变截面弹性直杆纵振动分析的小波——DQ法

在经典微分求积(DQ)法基础上, 根据多分辨分析理论, 以尺度函数为基础构造插 值基函数, 形成了新的微分方程边值问题的求解方法------小波--DQ法, 并应用该方法分析了变截面弹性直杆的纵振动问题, 给出了其频率方程, 计算出了不同参数下固 支--固支, 自由--自由楔形直杆和锥形直杆的固有频率, 数值结果表明该方法是一个简单易行高精度的方法, 该方法可以推广应用于其他力学问题的数值分析.     

车辆纵向非光滑多体动力学建模与数值算法研究

在建立车辆纵向多体系统的动力学模型中, 将车身与车轮视为刚体, 两者通过减振器链接; 将传动系统视为一个圆盘通过扭簧和阻尼器与驱动轮连接; 将车轮与路面间的接触力简化为法向约束力、摩擦力和滚阻力偶,其中摩擦力的模型采用库仑干摩擦模型. 采用笛卡尔坐标作为该系统的广义坐标用于描述该系统的位形, 给出系统单双边的约束方程, 应用第一类拉格朗日方法建立了系统的动力学方程. 由于摩擦与滚阻的非光滑性, 使得该系统动力学方程不连续. 为便于计算, 建立了车轮与路面接触点的相对切向加速度与摩擦力余量的互补条件、车轮角加速度与滚阻力偶余量的互补条件, 以及车轮轮心法向加速度与路面法向约束力的互补条件. 将接触—分离、黏滞—滑移的判断问题转化成线性互补问题的求解, 并给出了具有约束稳定化的基于事件驱动法的数值计算方法. 最后, 应用该方法对车辆纵向多体系统进行了仿真, 分析了输出扭矩、摩擦及滚阻系数对其动力学行为的影响.      

考虑分布轴向力的细长杆横向振动与失稳分析

为量化梁、杆、柱的自重(下称分布轴向力)对静力失稳和动力横向振动的影响,在《材料力学》和《机械系统动力学》教材的基础上,建立了分布轴向力下的杆柱失稳和横向振动的力学、数学模型.采用有限差分法、伽辽金法和数值积分法获取计算结果.结果表明:考虑分布轴向力的杆柱横向振动固有频率随杆长增加而减小,杆柱失稳时一阶固有频率为0;分...     

表面微织构对球盘点接触润滑摩擦性能的影响

基于统一Reynolds方程系统模型开展了富油点接触工况下微织构表面润滑摩擦性能的数值模拟研究.在通过实验标定数值模拟中润滑剂流变参数的基础上,系统分析了微织构表面摩擦系数周期变化的全过程,初步揭示了微织构的减摩机理.结果表明：数值模拟结果与实验结果有较好的吻合;瞬时摩擦系数达到最小值时,微坑单元一般处于名义Hertz接触区域的前边界;当微坑运动到Hertz接触区域内时,微坑前沿局部膜厚减小,而微坑后边沿膜厚局部增大,形成局部膜厚增大区;局部膜厚增大区的大小对微织构的润滑摩擦性能有较大影响,其面积越大,减摩效果越好.     

一类双自由度碰振系统运动分析

基于Poincare映射方法对一类两自由度碰撞系统进行了分析。经过详细的理论演算得到单碰周期n的次谐运动的存在性判据和稳定性条件,给出计算Jacobi矩阵特征值的公式。数值模拟表明,该方法具有令人满意的结果。此外,还讨论了当不满足所提出的单碰周期n次谐运动的存在性条件时,可能会出现的运动形式。     

Non-linear dynamics of a mechanical system with a frictional unilateral constraint

We analyze the dynamics of a two-dimensional system constituted by two masses subjected to elastic, gravitational and viscous forces and constrained by a moving frictional mono-lateral surface. The model exhibits a time-varying dynamics capable of reproducing the hopping phenomenon, an unwanted phenomenon observed in many applications such as the motion of a robotic arm on a surface or that of a wiper on a windscreen. The system dynamics, besides being affected by geometrical non-linearities, has a non-smooth nature due to the impact and friction laws involved in the model. The complexity of the resulting equations and of the transition conditions require the problem to be solved numerically. Various periodic motions are found and the effect of varying the system parameters, in particular the friction coefficient, is investigated. Finally, simulations are used to gain some insight the behavior of the windscreen wiper.