多自由度强非线性耦合参激系统随机响应计算方法

本文采用非高斯矩闭合方法计算了一个含参数激振的二自由度强非线性耦合系统受随机激励的响应方差。用逐次消除法消除方程中的惯性耦合后推导出直到6阶的矩方程,采用中心累积量截断技术计算了白噪音激励下的响应方差,并与时域直接积分和数字模拟方法的结果进行比较,取得了一些有意义的结论。     

基于等效线性化法的弱非线性体系下的地震响应分析

对弱非线性单自由度非自由振动系统结构在白噪声地震激励下的随机响应问题进行了研究。首先,建立了结构的弱非线性响应方程;然后根据能量平衡原理对运动方程进行了等效线性化;最后获得了金井清谱(Kanai-Tajimi谱)下等效线性化方程的位移响应及速度响应的方差,将弱非线性随机响应问题简化为求解线性随机响应问题,并给出了算例验证。结果表明,该方法可以准确得到弱非线性随机响应问题的解。同时,建立了单自由度结构弱非线性随机响应频域分析的统一解析解法。     

一种考虑结构非线性的连续突风气弹响应计算方法

连续突风气弹响应在工程中一般采用随机理论进行分析求解,通过将全机频响函数,与激励功率谱密度积分,即可得到各自由度均方根值。该方法基于频域线性求解,无法直接进行非线性响应分析。为此本研究基于匹配滤波器理论,提出了一套可考虑结构非线性的连续突风响应工程简化计算方法。通过引入突风滤波器,串联至全机气弹动力学系统,构成了连续突风的时域分析模型。针对所选取自由度的均方根值,通过匹配滤波器理论,计算得到一套各自由度“时间匹配”的响应,在此基础上,针对非线性结构模型,通过调整脉冲激励幅值,寻找最大的响应作为非线性结果,最终建立了一套适合于工程应用的、考虑结构非线性的连续突风响应计算方法。算例表明,本研究所建立的方法针对连续突风,可获得满足工程精度要求的各自由度“时间匹配”的响应,同时可用于结构非线性响应计算。     

加性非平稳激励下结构速度响应的FPK方程降维

结构在随机激励下的非线性响应分析是具有高度挑战性的困难问题. 对于白噪声或过滤白噪声激励,求解FPK方程将获得结构响应 的精确解. 遗憾的是,对于非线性多自由度系统,FPK方程难以直接求解. 事实上,其数值解法严重受限于方程维度,而解析求解 则仅适用于少数特定的系统,且多是稳态解. 因此,将FPK方程进行降维,是求解高维随机动力响应分析问题的重要途径. 本文针 对幅值调制的加性白噪声激励下多自由度非线性结构的非平稳随机响应分析问题,将联合概率密度函数满足的高维FPK方程进行降 维. 针对结构速度响应概率密度函数求解,通过引入等价漂移系数,原FPK方程可转化为一维FPK型方程. 建议了构造等价漂移系数 的条件均值函数方法. 进而,采用路径积分方法求解降维FPK型方程,得到速度概率密度函数的数值解答. 结合单自由度Rayleigh 振子、十层线性剪切型框架和非线性剪切型框架结构在幅值调制的加性白噪声激励下的非平稳速度响应求解,讨论了本文方法的精 度和效率,验证了其有效性.      

加性非平稳激励下结构速度响应的FPK方程降维

结构在随机激励下的非线性响应分析是具有高度挑战性的困难问题.对于白噪声或过滤白噪声激励,求解FPK方程将获得结构响应的精确解.遗憾的是,对于非线性多自由度系统,FPK方程难以直接求解.事实上,其数值解法严重受限于方程维度,而解析求解则仅适用于少数特定的系统,且多是稳态解.因此,将FPK方程进行降维,是求解高维随机动力响应分析问题的重要途径.本文针对幅值调制的加性白噪声激励下多自由度非线性结构的非平稳随机响应分析问题,将联合概率密度函数满足的高维FPK方程进行降维.针对结构速度响应概率密度函数求解,通过引入等价漂移系数,原FPK方程可转化为一维FPK型方程.建议了构造等价漂移系数的条件均值函数方法.进而,采用路径积分方法求解降维FPK型方程,得到速度概率密度函数的数值解答.结合单自由度Rayleigh振子、十层线性剪切型框架和非线性剪切型框架结构在幅值调制的加性白噪声激励下的非平稳速度响应求解,讨论了本文方法的精度和效率,验证了其有效性.     

具有组合非线性阻尼的非线性能量阱振动抑制响应分析

非线性能量阱是一种振动能量吸收装置,其在结构振动抑制中具有十分重要的作用.文章对具有组合非线性阻尼非线性能量阱的系统进行振动抑制相关的分析.首先对具有组合非线性阻尼非线性能量阱的系统进行理论模型的描述,对系统模型的运动方程利用复变量平均法进行推导,得到系统的慢变方程.其次对系统的慢变方程运用多尺度法进行强调制响应的分析,通过对系统进行慢不变流形和相轨迹的研究,描述系统强调制响应发生的条件基础.此外,还利用一维映射对系统进行分析,揭示外激励幅值对强调制响应存在时频率失谐系数取值区间的影响规律.最后利用能量谱、时间响应和庞加莱映射对耦合组合非线性阻尼非线性能量阱系统进行了振动抑制的相关研究,揭示组合非线性阻尼的非线性能量阱不同阻尼比、阻尼和刚度对其振动抑制效果的影响规律,得出组合非线性阻尼非线性能量阱和主结构响应存在一致性的现象,并验证所提出的组合非线性阻尼非线性能量阱模型具有较好的振动抑制能力.     

车辆-轨道系统垂向随机振动的辛方法分析

将轨道视为无限长的周期结构,建立车辆轨道垂向耦合模型. 使用虚拟激励法将随机的轨道不平顺激励转化为确定性的简谐激励,再用辛数学方法求解轨道结构的频率响应特性和耦合系统的响应功率谱. 整个计算模型只有26个自由度,求解过程快速而精确.数值算例中,将该方法与常规有限元方法进行了比较,验证了方法的高效性和正确性,讨论了车辆速度对系统随机响应的影响.      

参数激励下非线性受控系统的动力学和稳定性

研究两自由度参数激励系统的非线性动力学与控制问题．利用Lagrange方程建立含反馈控制的参激捅及其驱动机构组成的系统动力学方程，以多尺度方法获得一阶近似控制方程．然后，对系统受一阶摸态参激主共振与一、二阶模态间3：1内共振联合作用下的幅额响应及其稳定性，以及反馈参数对系统稳态行为的影响作了详细分析．结果表明，响应的稳定域位置和大小取决于位移反馈，位移立方反馈改变了系统的非线性程度，速度反馈类似于阻尼，可使系统呈现自激振动特性．     

考虑固体域几何非线性的固液耦合系统瞬态有限元分析

本文基于理想可压缩流体压力传输方程及增量型ＴｏｔａＩＬａｇｒａｎｇｉａｎ列式，建立了考虑固体域几何非线性的固液耦合系统有限元方程，并开展了非线性分析工作。本文的方法可用于流体域远大于固体域，从而可仅仅研究固体域非线性响应的情况。由于流体域每个节点仅有一个未知量（压力），故流体域的计算规模得以大大减缩。同时，因为本文方法仅需对固体域进行非线性分析，故显著地简化了固液耦合非线性问题的求解工作。本文的方法亦可用于三维分析。     

非线性振动中的非高斯矩方法

本文在Fokker-Planck方程的基础上,把非高斯矩方法推广到更具有一般性的非线性系统,用来研究非线性系统对高斯白噪声激励的响应。此法是一种有实用价值的方法。它对系统的非线性没有限制,而且,适用于平稳和非平稳情况。对本质非线性问题,仍能达到较高的精度,数值计算也是方便的。     

基于非线性谐振电路的双稳态俘能器的俘能与动力学特性研究

双稳态俘能器可实现宽频和高效的俘能效果.目前的研究主要在双稳态结构中接入单一电阻电路进行俘能.本文将非线性RLC (电阻-电感-电容)谐振电路引入到三弹簧式双稳态结构中,构建两自由度非线性系统,以实现俘能特性的提升.设计永磁体与线圈的构型,获得了非线性机电耦合系数.推导并得到了两自由度非线性俘能器的控制方程.利用谐波平衡法推导得到了系统的电流与位移的频率响应关系.基于雅可比矩阵对解的稳定性进行了判别.将解析解与数值解进行了对比验证.结果表明,在双稳态俘能器中引入非线性二阶谐振电路不仅有利于低频俘能,还可进一步提升俘能响应,拓宽俘能带宽.相同的电路参数下,与线性电路相比非线性电路可通过电流的倍频现象实现结构更低频率的能量俘获.减小谐振电路与双稳态结构共振频率之比,增加基础激励幅值,减小静平衡点之间的距离均可提升俘能器的俘能效果.通过调控谐振电路与双稳态共振频率之比和基础激励幅值等参数,可实现系统单倍周期响应、多倍周期响应及混沌响应之间的切换.     

Identification of Weakly Nonlinear Systems Using Describing Function Inversion

In this paper describing functions inversion is used and the restoring force of a nonlinear element in a MDOF system is characterized. The describing functions can be obtained using linearized frequency response functions (FRFs). The response of the system to harmonic excitation forces at distinct frequencies close to the resonant frequency results in linearized FRFs. The nonlinear system can be approximated at each excitation frequency by an equivalent linear system. This approximation leads to calculation of the first-order describing functions. By having the experimental describing functions calculated and the system’s responses corresponding to the nonlinear element (measured or interpolated), nonlinear parameter identification can be performed. Two numerical and experimental case studies are provided to show the applicability of this method.     

Deterministic and Stochastic Response of Nonlinear Coupled Bending-Torsion Modes in a Cantilever Beam

The purpose of this study is to understand the main differences between the deterministic and random response characteristics of an inextensible cantilever beam (with a tip mass) in the neighborhood of combination parametric resonance. The excitation is applied in the plane of largest rigidity such that the bending and torsion modes are cross-coupled through the excitation. In the absence of excitation, the two modes are also coupled due to inertia nonlinearities. For sinusoidal parametric excitation, the beam experiences instability in the neighborhood of the combination parametric resonance of the summed type, i.e., when the excitation frequency is in the neighborhood of the sum of the first bending and torsion natural frequencies. The dependence of the response amplitude on the excitation level reveals three distinct regions: nearly linear behavior, jump phenomena, and energy transfer. In the absence of nonlinear coupling, the stochastic stability boundaries are obtained in terms of sample Lyapunov exponent. The response statistics are estimated using Monte Carlo simulation, and measured experimentally. The excitation center frequency is selected to be close to the sum of the bending and torsion mode frequencies. The beam is found to experience a single response, two possible responses, or non-stationary responses, depending on excitation level. Experimentally, it is possible to obtain two different responses for the same excitation level by providing a small perturbation to the beam during the test.     

Nonlinear dynamic response of a functionally graded plate with a through-width surface crack

A nonlinear vibration analysis of a simply supported functionally graded rectangular plate with a through-width surface crack is presented in this paper. The plate is subjected to a transverse excitation force. Material properties are graded in the thickness direction according to exponential distributions. The cracked plate is treated as an assembly of two sub-plates connected by a rotational spring at the cracked section whose stiffness is calculated through stress intensity factor. Based on Reddy’s third-order shear deformation plate theory, the nonlinear governing equations of motion for the FGM plate are derived by using the Hamilton’s principle. The deflection of each sub-plate is assumed to be a combination of the first two mode shape functions with unknown constants to be determined from boundary and compatibility conditions. The Galerkin’s method is then utilized to convert the governing equations to a two-degree-of-freedom nonlinear system including quadratic and cubic nonlinear terms under the external excitation, which is numerically solved to obtain the nonlinear responses of cracked FGM rectangular plates. The influences of material property gradient, crack depth, crack location and plate thickness ratio on the vibration frequencies and transient response of the surface-racked FGM plate are discussed in detail through a parametric study.     

Reduced-order models for nonlinear vibrations of fluid-filled circular cylindrical shells: Comparison of POD and asymptotic nonlinear normal modes methods

The aim of the present paper is to compare two different methods available for reducing the complicated dynamics exhibited by large amplitude, geometrically nonlinear vibrations of a thin shell. The two methods are: the proper orthogonal decomposition (POD), and an asymptotic approximation of the nonlinear normal modes (NNMs) of the system. The structure used to perform comparisons is a water-filled, simply supported circular cylindrical shell subjected to harmonic excitation in the spectral neighbourhood of the fundamental natural frequency. A reference solution is obtained by discretizing the partial differential equations (PDEs) of motion with a Galerkin expansion containing 16 eigenmodes. The POD model is built by using responses computed with the Galerkin model; the NNM model is built by using the discretized equations of motion obtained with the Galerkin method, and taking into account also the transformation of damping terms. Both the POD and NNMs allow to reduce significantly the dimension of the original Galerkin model. The computed nonlinear responses are compared in order to verify the accuracy and the limits of these two methods. For vibration amplitudes equal to 1.5 times the shell thickness, the two methods give very close results to the original Galerkin model. By increasing the excitation and vibration amplitude, significant differences are observed and discussed. The response is investigated also for a fixed excitation frequency by using the excitation amplitude as bifurcation parameter for a wide range of variation. Bifurcation diagrams of Poincaré maps obtained from direct time integration and calculation of the maximum Lyapunov exponent have been used to characterize the system.     

Nonlinear heave-roll coupling and ship rolling

A nonlinear model for simulating the heave-roll motions of ships in following waves is presented. The parametric excitation is modeled by a Hill's type equation, instead of the conventional Mathieu's equation. The model includes not only the linear but also the quadratic coupling term. Instability conditions for parametrically excited rolling motions are derived using the harmonic balance method. The results are verified by numerical analyses. The effects of including the quadratic coupling term on the instability conditions and nonlinear responses are studied. The complex dynamic behaviour of the coupled system in the various instability regions is also investigated. Bifurcations of the flip, fold and pitchfork types are observed in the Poincaré mapping of the numerically simulated responses. Chaotic motions leading to boundary crises and inevitable capsize are also reported.     

Nonlinear dynamic analysis of a quarter vehicle system with external periodic excitation

In this paper, a nonlinear dynamic model of a quarter vehicle with nonlinear spring and damping is established. The dynamic characteristics of the vehicle system with external periodic excitation are theoretically investigated by the incremental harmonic balance method and Newmark method, and the accuracy of the incremental harmonic balance method is verified by comparing with the result of Newmark method. The influences of the damping coefficient, excitation amplitude and excitation frequency on the dynamic responses are analyzed. The results show that the vibration behaviors of the vehicle system can be control by adjusting appropriately system parameters with the damping coefficient, excitation amplitude and excitation frequency. The multi-valued properties, spur-harmonic response and hardening type nonlinear behavior are revealed in the presented amplitude-frequency curves. With the changing parameters, the transformation of chaotic motion, quasi-periodic motion and periodic motion is also observed. The conclusions can provide some available evidences for the design and improvement of the vehicle 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.     

On energy pumping, synchronization and beat phenomenon in a nonideal structure coupled to an essentially nonlinear oscillator

In this paper, we investigated the effectiveness of the linear electromechanical vibration absorber (LEVA) and a nonlinear electromechanical vibration absorber (NEVA) in the vibration attenuation for nonideal structures (NIS). This electromechanical damping device consists of an electrical system coupled magnetically to a mechanical structure under a nonideal excitation. An analysis of the effects of the parameters of coupling and of nonlinear coefficients with increasing of constant torque of the DC motor is presented.     

Static and dynamic behaviors of belt-drive dynamic systems with a one-way clutch

This paper focuses on the nontrivial equilibrium and the steady-state periodic response of belt-drive system with a one-way clutch and belt flexural rigidity. A nonlinear piecewise discrete–continuous dynamic model is established by modeling the motions of the translating belt spans as transverse vibrations of axially moving viscoelastic beams. The rotations of the pulleys and the accessory are also considered. Furthermore, the transverse vibrations and the rotation motions are coupled by nonlinear dynamic tension. The nontrivial equilibriums of the belt-drive system are obtained by an iterative scheme via the differential and integral quadrature methods (DQM and IQM). Moreover, the periodic fluctuation of the driving pulley is modeled as the excitation of the belt-drive system. The steady-state periodic responses of the dynamic system are, respectively, studied via the high-order Galerkin truncation as well as the DQM and IQM. The time histories of the system are numerically calculated based on the 4th Runge–Kutta time discretization method. Furthermore, the frequency–response curves are presented from the numerical solutions. Based on the steady-state periodic response, the resonance areas of the dynamic system are obtained by using the frequency sweep. Moreover, the influences of the truncation terms of the Galerkin method, such as 6-term, 8-term, 10-term, 12-term, and 16-term, are investigated by comparing with the DQM and IQM. Numerical results demonstrate that the one-way clutch reduces the resonance responses of the belt-drive system via the torque-transmitting directional function. Furthermore, the comparisons in numerical examples show that the investigation on steady-state responses of the belt-drive system with a one-way clutch and belt flexural rigidity needs 16-term truncation