ANCF索梁单元应变耦合问题与模型解耦

索粱结构在土木工程、航空航天等领域有着广泛的应用.在各类索梁动力学建模方法中,由于绝对节点坐标方法(absolute nodal coordinate formulation,ANCF)能够描述柔性体的大变形和大转动问题,因此非常适合大变形索梁结构的动力学建模.对绝对节点坐标索梁单元的应变进行分析可知,弯曲变形会引起单元内部轴向应变的不均匀分布,即单元轴向应变与弯曲应变相互耦合.这种应变耦合效应使单元产生伪应变能,导致单元刚度增大,造成单元失真.分析不同弯曲角下的单元应变及应变能可知,弯曲变形越大,单元失真越严重.通过构造等效一维杆单元重新描述轴向应变,实现了轴向应变与弯曲应变解耦.在此基础上推导广义弹性力,得到了绝对节点坐标索梁单元的应变解耦模型.对解耦前后的两种梁模型进行静力学和动力学仿真,结果表明;解耦模型消除了单元伪应变,相比原模型表现出更好的收敛性和曲率连续性,在相同单元数目下具有更高的精度.同时由于解耦模型降低了单元刚度,因此相比原模型,速度曲线中不再有高频振动.     

Stretch and Rotation of a Suspended Cable Subject to Transverse Fluid Excitation

ＩｎｔｒｏｄｕｃｔｉｏｎＣａｂｌｅｓａｒｅｅｘｔｅｎｓｉｖｅｌｙｕｓｅｄｉｎｖａｒｉｏｕｓｆｉｅｌｄｓｏｆｅｎｇｉｎｅｅｒｉｎｇ .Ｔｈｅｏｂｓｅｒｖｅｄｄｙｎａｍｉｃｓｏｆｃａｂｌｅｓｉｎｐｒａｃｔｉｃｅｏｆｔｅｎｅｘｈｉｂｉｔｓｖｅｒｙｃｏｍｐｌｅｘｂｅｈａｖｉｏｒ.Ｈｅｒｅａｒｅａｆｅｗ ,ｂｕｔｔｙｐｉｃａｌｅｘａｍｐｌｅｓ:(ａ)Ｔｈｅｓｔｅｅｌｃａｂｌｅｏｆｓｈｉｐｐｉｎｇｃｒａｎｅｍａｙｈａｖｅｔｈｅａｐｅｒｉｏｄｉｃｓｗｉｎｇｔｈａｔｐｒｏｄｕｃｅｓｐｏｏｒｐｏｓｉｔｉｏｎａｃｃｕｒａｃ…     

Fluid Flow-Induced Nonlinear Vibration of Suspended Cables

The nonlinear interaction of the first two in-plane modes of a suspended cable with a moving fluid along the plane of the cable is studied. The governing equations of motion for two-mode interaction are derived on the basis of a general continuum model. The interaction causes the modal differential equations of the cable to be non-self-adjoint. As the flow speed increases above a certain critical value, the cable experiences oscillatory motion similar to the flutter of aeroelastic structures. A co-ordinate transformation in terms of the transverse and stretching motions of the cable is introduced to reduce the two nonlinearly coupled differential equations into a linear ordinary differential equation governing the stretching motion, and a strongly nonlinear differential equation for the transverse motion. For small values of the gravity-to-stiffness ratio the dynamics of the cable is examined using a two-time-scale approach. Numerical integration of the modal equations shows that the cable experiences stretching oscillations only when the flow speed exceeds a certain level. Above this level both stretching and transverse motions take place. The influences of system parameters such as gravity-to-stiffness ratio and density ratio on the response characteristics are also reported.     

旋转运动柔性梁的假设模态方法研究

采用假设模态法对旋转运动柔性梁的动力特性进行研究，给出简化的控制模型. 首先采用Hamilton原理和假设模态离散化方法，在计入柔性梁由于横向变形而引起的轴向变形的二阶耦合量的条件下，推导出基于柔性梁变形位移场一阶完备的一次近似耦合模型，然后对该模型进行简化，忽略柔性梁纵向变形的影响，给出一次近似简化模型，最后将采用假设模态离散化方法的结果与采用有限元离散化方法的结果进行了对比研究. 研究中考虑了两种情况：非惯性系下的动力特性研究和系统大范围运动为未知的动力特性研究. 研究结果显示，当系统大范围运动为高速时，在假设模态离散化方法中应增加模态数目，较少的模态数目将导致较大误差. 一次近似简化模型能够较好地反映出系统的动力学行为，可用于主动控制设计的研究.     

行进绳索在横向流体激励下的运动

给出横向流体对行进绳索的作用力描述,建立了绳索的动力学方程。由于该方程具有零刚度特征,引入Pilipchuk变换,以自平衡状态起度量的径向振动和回转运动来描述绳的运动,获得非零刚度系统。然后,用两变量参数摄动法求得绳索关于扩张振动和回转运动的约化型,使得绳索运动可近似由一个二维动力系统来描述。最后,用数值方法讨论了行进速度和重力对绳索运动形态的影响。     

Dynamic response of a deepwater riser subjected to combined axial and transverse excitation by the nonlinear coupled model

In offshore engineering long slender risers are simultaneously subjected to both axial and transverse excitations. The axial load is the fluctuating top tension which is induced by the floater’s heave motion, while the transverse excitation comes from environmental loads such as waves. As the time-varying axial load may trigger classical parametric resonance, dynamic analysis of a deepwater riser with combined axial and transverse excitations becomes more complex. In this study, to fully capture the coupling effect between the planar axial and transverse vibrations, the nonlinear coupled equations of a riser’s dynamic motion are formulated and then solved by the central difference method in the time domain. For comparison, numerical simulations are carried out for both linear and nonlinear models. The results show that the transverse displacements predicted by both models are similar to each other when only the random transverse excitation is applied. However, when the combined axial dynamic tension and transverse wave forces are both considered, the linear model underestimates the response because it ignores the coupling effect. Thus the coupled model is more appropriate for deep water. It is also found that the axial excitation can significantly increase the riser’s transverse response and hence the bending stress, especially for cases when the time-varying tension is located at the classical parametric resonance region. Such time-varying effects should be taken into account in fatigue safety assessment.     

风作用下覆冰悬索的非平面非线性动力学

主要研究侧向风载荷作用下小垂度覆冰悬索的非线性非平面运动的复杂动力学.根据分析力学、弹性力学和空气动力学理论,建立覆冰悬索3个自由度非线性振动的偏微分运动方程,并对其进行无量纲化,运用Galerkin方法对偏微分运动方程进行离散得到3个自由度的常微分方程,再利用多尺度法得到面内主共振2:1内共振的平均方程.利用数值方法研究悬索的非线性运动,结果表明系统呈现周期、多倍周期、概周期和混沌运动的规律.     

悬索主共振响应的多尺度解与同伦分析解

利用哈密顿变分原理,引入拟静态假设,建立了悬索面内非线性运动方程,并采用Galerkin方法对其进行离散。接着运用多尺度法和同伦分析法得到了悬索前两阶模态主共振响应的近似解。为验证这两种分析方法的适用性,同时采用龙格-库塔法对方程直接进行了数值积分。数值计算结果表明,随着悬索垂跨比以及振幅的增加,由多尺度法与同伦分析法得到的幅频响应曲线存在明显的定性与定量的差别,而同伦分析法结果与数值法的结果更加接近。最后比较了两种分析方法得到的位移场与索力时程响应曲线。     

Three-dimensional static and dynamic stiffness analyses of the cable driven parallel robot with non-negligible cable mass and elasticity

This study investigates the three-dimensional static and dynamic stiffness analyses of the cable driven parallel robot by considering cable mass, elasticity, and mass of end-effector. According to these models, optimization of cable tensions and cable lengths using fminimax solver is performed to determine the static stiffness. Static and dynamic stiffness of the cables are obtained with simulations. Results show that analyses in three-dimension are very important to measure the actual performance of the cable driven parallel robot. This demonstrates potential for general applicability and motion of the cable driven parallel robot.     

Global dynamics of the cable under combined parametrical and external excitations

The chaotic dynamics and global bifurcations of the suspended elastic cable under combined parametric and external excitations are investigated. The non-linear equations of motion of the elastic cable to small vibration of one support are derived. The averaged equations are obtained by using the method of multiple scales. Based on the averaged equations, the theory of normal form and Maple program are used to obtain the explicit expressions of normal form associated with a double zero and a pair of pure imaginary eigenvalues. On the basis of the normal form, global bifurcation analysis of the parametrically and externally excited suspended elastic cable is given by a global perturbation method developed by Kovacic and Wiggins. The chaotic motion of the elastic cable is also found by numerical simulation.     

Stochastic optimal control of cable vibration in plane by using axial support motion

A stochastic optimal control strategy for a slightly sagged cable using support motion in the cable axial direction is proposed.The nonlinear equation of cable motion in plane is derived and reduced to the equations for the first two modes of cable vibration by using the Galerkin method.The partially averaged Ito equation for controlled system energy is further derived by applying the stochastic averaging method for quasi-non-integrable Hamiltonian systems.The dynamical programming equation for the controlled system energy with a performance index is established by applying the stochastic dynamical programming principle and a stochastic optimal control law is obtained through solving the dynamical programming equation.A bilinear controller by using the direct method of Lyapunov is introduced.The comparison between the two controllers shows that the proposed stochastic optimal control strategy is superior to the bilinear control strategy in terms of higher control effectiveness and efficiency.     

Dynamic analysis of a flexible hub-beam system with tip mass

For a dynamic system of a rigid hub and a flexible cantilever beam, the traditional hybrid coordinate model assumes the small deformation in structural dynamics where axial and transverse displacements at any point in the beam are uncoupled. This traditional hybrid coordinate model is referred as the zeroth-order approximation coupling model in this paper, which may result in divergence to the dynamic problem of some rigid–flexible coupling systems with high rotational speed. In this paper, characteristics of a flexible hub-beam system with a tip mass is studied. Based on the Hamilton theory and the finite element discretization method, and in consideration of the second-order coupling quantity of the axial displacement caused by the transverse displacement of the beam, the rigid–flexible coupling dynamic model (referred as the first-order approximation coupling (FOAC) model in this paper) and the corresponding model in non-inertial system for the flexible hub-beam system with a tip mass are presented firstly, then the dynamic characteristics of the system are studied through numerical simulations under twos cases: the large motion of the system is known and is unknown. Simulation and comparison studies using both the traditional zeroth-order model and the proposed first-order model show that even small tip mass may affect dynamic characteristics of the system significantly, which may result in the largening of vibrating amplitude and the descending of vibrating frequency of the beam, and may affect end position of the hub-beam system as well. The effect of the tip mass becomes large along with the increasing of rotating speed of large motion of the system. When the large motion of the system is at low speed, the traditional ZOAC model may lead to a large error, whereas the proposed FOAC model is valid. When the large motion is at high speed, the ZOAC model may result in divergence to the dynamic problem of the flexible hub-beam system, while the proposed second model can still accurately describe the dynamic hub-beam system.     

Two-to-one resonant multi-modal dynamics of horizontal/inclined cables. Part I: Theoretical formulation and model validation

This paper is first of the two papers dealing with analytical investigation of resonant multi-modal dynamics due to 2:1 internal resonances in the finite-amplitude free vibrations of horizontal/inclined cables. Part I deals with theoretical formulation and validation of the general cable model. Approximate nonlinear partial differential equations of 3-D coupled motion of small sagged cables – which account for both spatio-temporal variation of nonlinear dynamic tension and system asymmetry due to inclined sagged configurations – are presented. A multi-dimensional Galerkin expansion of the solution of nonplanar/planar motion is performed, yielding a complete set of system quadratic/cubic coefficients. With the aim of parametrically studying the behavior of horizontal/inclined cables in Part II [25], a second-order asymptotic analysis under planar 2:1 resonance is accomplished by the method of multiple scales. On accounting for higher-order effects of quadratic/cubic nonlinearities, approximate closed-form solutions of nonlinear amplitudes, frequencies and dynamic configurations of resonant nonlinear normal modes reveal the dependence of cable response on resonant/nonresonant modal contributions. Depending on simplifying kinematic modeling and assigned system parameters, approximate horizontal/inclined cable models are thoroughly validated by numerically evaluating statics and non-planar/planar linear/non-linear dynamics against those of the exact model. Moreover, the modal coupling role and contribution of system longitudinal dynamics are discussed for horizontal cables, showing some meaningful effects due to kinematic condensation.     

Large Amplitude Three-Dimensional Free Vibrations of Inclined Sagged Elastic Cables

The nonlinear characteristics in the large amplitude three-dimensionalfree vibrations of inclined sagged elastic cables are investigated. Amodel formulation which is not limited to cables having smallsag-to-span ratios and takes into account the axial deformation effectis considered. Based on a multi-degree-of-freedom cable model, a finitedifference discretization is employed within a numerical solution of thegoverning equations of three-dimensional coupled motion. Variousnumerical examples of arbitrarily inclined sagged cables with initialout-of-plane or in-plane motions are carried out for the case of aspecified end tension. The major findings consist of highlighting theextent of two-and three-dimensional nonlinear couplings, the occurrenceof nonlinear dynamic tensions, and the meaningfulness of modaltransition phenomena ensuing from the activation of various internalresonance conditions. The influence of cable inclination on thenonlinear dynamic behavior is also evaluated. Comprehensive discussionand comparison of large amplitude free vibrations of horizontal andinclined sagged cables are presented.     

变化外形的水线对斜拉索风雨激振的激发作用

讨论变化外形的水线对斜拉桥拉索风雨激振的激发作用.首先,建立考虑水线外形变化的斜拉桥拉索风雨激振计算模型,并将由采用变化水线外形的模型得到的拉索最大振幅与固定水线外形模型的计算结果及试验实测结果对比,对模型进行验证;在模型得到验证的基础上,进一步讨论拉索和变化水线外形的水线振动稳定性及水线对拉索风雨激振贡献.计算结果表明,在四个典型工况下,考虑水线外形变化的模型得到的拉索计算振幅固定水线外形模型能更好地与试验实测振幅吻合;理论分析发现拉索速度和水线速度的耦合作用对水线运动的激发起到主导作用;水线的加速发散发生在拉索大幅激振之前.没有水线加速发散的时程中,拉索均不发生大幅激振.     

Coupled global dynamics of an axially moving viscoelastic beam

The nonlinear global forced dynamics of an axially moving viscoelastic beam, while both longitudinal and transverse displacements are taken into account, is examined employing a numerical technique. The equations of motion are derived using Newton′s second law of motion, resulting in two partial differential equations for the longitudinal and transverse motions. A two-parameter rheological Kelvin–Voigt energy dissipation mechanism is employed for the viscoelastic structural model, in which the material, not partial, time derivative is used in the viscoelastic constitutive relations; this gives additional terms due to the simultaneous presence of the material damping and the axial speed. The equations of motion for both longitudinal and transverse motions are then discretized via Galerkin’s method, in which the eigenfunctions for the transverse motion of a hinged-hinged linear stationary beam are chosen as the basis functions. The subsequent set of nonlinear ordinary equations is solved numerically by means of the direct time integration via modified Rosenbrock method, resulting in the bifurcation diagrams of Poincaré maps. The results are also presented in the form of time histories, phase-plane portraits, and fast Fourier transform (FFTs) for specific sets of parameters.     

Multi-mode interactions in vortex-induced vibrations of flexible curved/straight structures with geometric nonlinearities

A general low-order fluid–structure interaction model capable of evaluating the multi-mode interactions in vortex-induced vibrations of flexible curved/straight structures is presented. Cross-flow motions due to unsteady lift forces of inclined sagged cables and tensioned beams in uniform currents are investigated. In contrast to a linear equation governing the transverse motion of straight beams or cables typically considered in the literature, coupled horizontal/vertical (axial/transverse) displacements and geometric nonlinearities of curved cable (straight beam) are accounted for. A distributed nonlinear wake oscillator is considered in the approximation of space–time varying hydrodynamics. This semi-empirical fluid force model in general depends on the mass-damping parameter and has further been modified to capture both the effects of varying initial curvatures of the inclined cylinder and the Reynolds number. Numerical simulations are performed in the case of varying flow velocities and parametric results highlight several meaningful aspects of vortex-induced vibrations of long flexible cylinders. These comprise multi-mode lock-in, sharing, switching and interaction features in the space and time domains, the estimated maximum modal and total amplitudes, the resonant nonlinear modes of flexible cylinders and their space–time modifications, and the influence of fluid/structure parameters. A shortcoming of single-mode or linear structural model is underlined. Some quantitative and qualitative comparisons of numerical/experimental results are discussed to demonstrate the validity and required improvement of the proposed modelling and analysis predictions.     

Numerical simulations of chaotic dynamics in a model of an elastic cable

The finite motions of a suspended elastic cable subjected to a planar harmonic excitation can be studied accurately enough through a single ordinary-differential equation with quadratic and cubic nonlinearities.The possible onset of chaotic motion for the cable in the region between the one-half subharmonic resonance condition and the primary one is analysed via numerical simulations. Chaotic charts in the parameter space of the excitation are obtained and the transition from periodic to chaotic regimes is analysed in detail by using phase-plane portraits, Poincaré maps, frequency-power spectra, Lyapunov exponents and fractal dimensions as chaotic measures. Period-doubling, sudden changes and intermittency bifurcations are observed.Part of this work was presented at the XVIIth Int. Congr. of Theor. and Appl. Mech., Grenoble, August 1988.     

Nonlinear dynamic response and active vibration control of the viscoelastic cable with small sag

IntroductionCablesareveryefficientstructuralmembersandhencehavebeenwidelyusedinmanylong_spanstructures,includingcable_supportbridges,guyedtowersandcable_supportroofs.Sincecablesarelight,veryflexibleandlightlydamped ,structuresutilizingcables,i.e .,cable_structuresystems,usuallyhavevariousdynamicproblems.Theirmodelsarethereforeverimportantinpredictingandcontrollingtheirresponses.Inthelastdecade,thenonlineardynamicvibrationandstabilitybehaviorofcablesandcable_structureshavedrawntheattentionofman…     

Nonlinear vibration analysis of viscoelastic cable with small sag

Both the inplane and out-of-plane transverse vibrations of a viscoelastic cable subjected to an initial stress distributing uniform on the cross section are studied. The constitution of the cable material is assumed to be of the hereditary integral type. The partial differential-integral equations of motion are derived first. Then by applying Galerkin's method, the governing equations are reduced to a set of second-order nonlinear differential-integral equations which are solved by finite difference numerical integration procedures. Finally, the effects of the viscosity parameter and the elastic parameter on the transient amplitudes of the first mode are investigated by numerical simulation. Project supported by the National Natural Science Foundation of China (No. 59635140) and the National Postdoctoral Foundation of China.