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.     

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.     

Nonlinear transient response of stay cable with viscoelasticity damper in cable-stayed bridge

Taking the bending stiffness, static sag, and geometric non-linearity into consideration, the space nonlinear vibration partial differential equations were derived. The partical differential equations were discretized in space by finite center difference approximation, then the nonlinear ordinal differential equations were obtained. A hybrid method involving the combination of the Newmark method and the pseudo-force strategy was proposed to analyze the nonlinear transient response of the inclined cable-dampers system subjected to arbitrary dynamic loading. As an example, two typical stay cables were calculated by the present method. The results reveal both the validity and the deficiency of the viscoelasticity damper for vibration control of stay cables. The efficiency and accuracy of the proposed method is also verified by comparing the results with those obtained by using Runge-Kutta direct integration technique. A new time history analysis method is provided for the research on the stay cable vibration control.     

Analysis of in-plane 1:1:1 internal resonance of a double cable-stayed shallow arch model with cables' external excitations

The nonlinear dynamic behaviors of a double cable-stayed shallow arch model are investigated under the one-to-one-to-one internal resonance among the lowest modes of cables and the shallow arch and external primary resonance of cables. The in-plane governing equations of the system are obtained when the harmonic excitation is applied to cables. The excitation mechanism due to the angle-variation of cable tension during motion is newly introduced. Galerkin's method and the multi-scale method are used to obtain ordinary differential equations(ODEs) of the system and their modulation equations, respectively. Frequency-and force-response curves are used to explore dynamic behaviors of the system when harmonic excitations are symmetrically and asymmetrically applied to cables. More importantly, comparisons of frequency-response curves of the system obtained by two types of trial functions, namely, a common sine function and an exact piecewise function, of the shallow arch in Galerkin's integration are conducted.The analysis shows that the two results have a slight difference; however, they both have sufficient accuracy to solve the proposed dynamic system.     

NONLINEAR VIBRATION OF THE COUPLED STRUCTURE OF SUSPENDED-CABLE-STAYED BEAM-1:2 INTERNAL RESONANCE

Through the Galerkin method the nonlinear ordinary differential equations （ODEs） in time are obtained from the nonlinear partial differential equations （PDEs） to describe the mo- tion of the coupled structure of a suspended-cable-stayed beam. In the PDEs, the curvature of main cables and the deformation of cable stays are taken into account. The dynamics of the struc- ture is investigated based on the ODEs when the structure is subjected to a harmonic excitation in the presence of both high-frequency principle resonance and 1：2 internal resonance. It is found that there are typical jumps and saturation phenomena of the vibration amplitude in the struc- ture. And the structure may present quasi-periodic vibration or chaos, if the stiffness of the cable stays membrane and frequency of external excitation are disturbed.     

斜拉桥索力优化的强次可行序列二次规划法

提出了一种斜拉桥索力优化的实用方法-强次可行序列二次规划法. 该方法通过建立斜拉 桥索力优化的非线性规划模型，以主梁和索塔的弯曲应变能为目标函数，斜拉索的索力为设 计变量，结构应力和索力为约束条件，计入大跨度斜拉桥各种几何非线性因素的影响，采用 强次可行序列二次规划算法进行优化求解，确定斜拉桥成桥合理状态的索力. 运用该方法对 某斜拉桥进行索力优化，结果表明该方法简单、有效.     

STATIC ANALYSIS OF CABLE STRUCTURE

Based on the nonlinear geometric relation between strain and displacement for flexible cable, the equilibrium equation under self-weight and influence of temperature was established and an analytical solution of displacement and tension distribution defined in Eulerian coordinate system was accurately obtained. The nonlinear algebraic equations caused by cable structure were solved directly using the modified Powell hybrid algorithm with high precision routine DNEQNE of Fortran. For example, a cable structure consisting of three cables jointly supported by a vertical spring and all the other ends fixed was calculated and compared with various methods by other scholars.     

Finding vibrations of inclined cable structures by approximately solving governing equations for exterior matrix

In this paper, how to compute the eigenfrequencies of the structures composed of a series of inclined cables is shown. The physics of an inclined cable can be complicated, so solving the differential equations even approximately is difficult. However, rather than solving the system of 4 first-order equations governing the dynamics of each cable, the governing equations are instead converted to a set of equations that the exterior matrix satisfies. Therefore, the exterior matrix method (EMM) is used without solving the original governing equations. Even though this produces a system of 6 first-order equations, the simple asymptotic techniques to find the first three terms of the perturbative solution can be used. The solutions can then be assembled to produce a 6 × 6 exterior matrix for a cable section. The matrices for each cable in the structure are multiplied together, along with the exterior matrices for each joint. The roots of the product give us the eigenfrequencies of the system.     

LT50m缩比模型悬索舱体系统的风振分析

在新一代大型射电望远镜（LT）的悬索舱体系统中增加了抑制风激振动的下拉结构控制索系，并分析了50m缩比模型的风振响应。首先，由非线性静力分析确定结构初始静态参考位形和初应力；其次，针对结构特点模拟作用在悬索和馈源舱上的随机风荷。在上述研究的基础上，在时域中模拟分析下拉索系对结构风振的抑制效果，同时通过典型节点时域响应的FFT变换，给出响应自功率谱图，进一步描述了振动控制的效果。     

Nonlocal continuum model for large deformation analysis of SLGSs using the kp-Ritz element-free method

A geometrically nonlinear large deformation analysis of SLGSs is presented using the element-free kp-Ritz method. Classical plate theory (CLP) is applied to describe the geometrically nonlinear behavior of SLGSs. Nonlocal elasticity theory is incorporated into CLP to take the small-scale effect into consideration. The system nonlinear equations are derived from the Ritz procedure based on the total Lagrangian formulation. The modified Newton–Raphson method and arc-length continuation are employed to solve the nonlinear equations. The efficiency of the element-free kp-Ritz method is verified through comparison with results reported in previous research. Numerical cases are studied to examine the influence of boundary conditions, aspect ratio, side length and nonlocal parameters on the nonlinear large deformation behavior of SLGSs. An interesting phenomenon is observed in that the nonlocal parameter effect is related to the mathematical expression of the transverse load.     

THE PARAMETRIC VARATIONAL PRINCIPLE AND NON-LINEAR FINITE ELEMENT METHOD FOR ANALYSIS OF ASTROMESH ANTENNA STRUCTURES

针对大型周边桁架式索网天线由拉索拉压模量不同引起的本构非线性和结构大变形引起的几何非线性问题,给出了基于参变量变分原理的几何非线性有限元方法. 首先针对含预应力索单元拉压模量不同分段描述的本构关系,通过引入参变量,导出了基于参变量及其互补方程的统一描述形式,避免了传统算法需要根据当前变形对索单元张紧/松弛状态的预测,提高了算法收敛性. 然后利用拉格朗日应变描述索网天线结构大变形问题,结合几何非线性有限元法,建立了基于参变量的非线性平衡方程和线性互补方程;并给出了牛顿-拉斐逊迭代法与莱姆算法相结合的求解算法. 数值算例验证了本文提出的算法比传统算法具有更稳定的收敛性和更高的求解精度,特别适合于大型索网天线结构的高精度变形分析和预测.     

斜拉桥合理成桥状态确定的一阶分析法

本文建立了斜拉桥索力优化的非线性规划模型，其中以斜拉桥主梁和索塔的弯曲应变能为目标函数，以各斜拉索的索力为设计变量，结构的应力及索力为约束条件，采用一阶分析法进行求解，用以确定成桥合理状态的索力。在计算中，考虑大跨度斜拉桥各种几何非线性因素的影响，并列出了优化模型的具体表达式及优化过程中的关键求解策略。应用该法和空间非线性有限元分析程序分析了某千米级斜拉桥方案的合理成桥状态，计算结果表明：该方法简单、有效。     

A parallel monolithic algorithm for the numerical simulation of large‐scale fluid structure interaction problems

A novel parallel monolithic algorithm has been developed for the numerical simulation of large‐scale fluid structure interaction problems. The governing incompressible Navier–Stokes equations for the fluid domain are discretized using the arbitrary Lagrangian–Eulerian formulation‐based side‐centered unstructured finite volume method. The deformation of the solid domain is governed by the constitutive laws for the nonlinear Saint Venant–Kirchhoff material, and the classical Galerkin finite element method is used to discretize the governing equations in a Lagrangian frame. A special attention is given to construct an algorithm with exact total fluid volume conservation while obeying both the global and the local discrete geometric conservation law. The resulting large‐scale algebraic nonlinear equations are multiplied with an upper triangular right preconditioner that results in a scaled discrete Laplacian instead of a zero block in the original system. Then, a one‐level restricted additive Schwarz preconditioner with a block‐incomplete factorization within each partitioned sub‐domains is utilized for the modified system. The accuracy and performance of the proposed algorithm are verified for the several benchmark problems including a pressure pulse in a flexible circular tube, a flag interacting with an incompressible viscous flow, and so on. John Wiley & Sons, Ltd.     

温度变化对端部激励斜拉索共振响应影响

基于增量热场理论,利用Hamilton变分原理,通过引入与张拉力和垂度相关的无量纲参数,建立了考虑温度变化影响下斜拉索非线性动力学模型,并推导其面内/外非线性运动微分方程。考虑斜拉索受端部激励,利用Galerkin法得到离散后的无穷维常微分方程组。面内和面外运动各取前两阶模态,向前和向后扫频,利用龙格-库塔法数值积分求解常微分方程组,得到共振区域的幅频响应曲线。算例分析表明,温度变化和斜拉索固有频率呈反比例关系;温度变化会导致斜拉索共振特性发生定性和定量的改变,如共振区间发生漂移、跳跃点位置发生移动、共振响应幅值发生改变;端部位移激励下,温度变化有可能导致斜拉索更多模态受到激发,从而影响各阶模态的能量以及模态间的能量传递。     

Accurate Higher-Order Analytical Approximate Solutions to Large-Amplitude Oscillating Systems with a General Non-Rational Restoring Force

A new approximate analytical approach for accurate higher-order nonlinear solutions of oscillations with large amplitude is presented in this paper. The oscillatory system is subjected to a non-rational restoring force. This approach is built upon linearization of the governing dynamic equation associated with the method of harmonic balance. Unlike the classical harmonic balance method, simple linear algebraic equations instead of nonlinear algebraic equations are obtained upon linearization prior to harmonic balancing. This approach also explores large parameter regions beyond the classical perturbation methods which in principle are confined to problems with small parameters. It has significant contribution as there exist many nonlinear problems without small parameters. Through some examples in this paper, we establish the general approximate analytical formulas for the exact period and periodic solution which are valid for small as well as large amplitudes of oscillation.     

Prediction of unwinding behaviors and problems of cables from inner-winding spool dispensers

Wire-guided control technologies are widely used to increase the targeting accuracy of advanced military weapons through the use of unwinding dispensers to guarantee that unwinding occurs without any problems, such as tangling or cutting. In this study, the transient behaviors of cables unwinding from inner-winding cylindrical spool dispensers are investigated. The cable is withdrawn from the spool dispenser at a constant velocity through a fixed point located along the axis of the spool dispenser. And when the cable is flown out of the dispenser, because several dynamic forces such as inertial forces, Coriolis forces, centrifugal forces, tensile forces, and fluid-resistance forces act on the cables, the cables exhibit highly nonlinear and complex unwinding behaviors which are called as unwinding balloons. For predicting these complex unwinding motions, the governing equations of motion for cables unwinding from a cylindrical spool dispenser of the inner-winding type are derived using the extended Hamilton’s principles for an open system in which mass can be transported at each boundary. Modified finite difference methods are used to discretize the spatial variables of the derived nonlinear partial differential equations. The time responses of the unwinding cables are calculated using Newmark time integration methods. Finally, in order to present numerical examples, an inner-winding spool dispenser that can be unwound up to a length of 50 km is designed by using simple geometrical relationships. The behaviors during unwinding from the designed inner spool dispenser are presented and modifications of the spool dispenser for avoiding unwinding problems are proposed by using the suggested numerical methods.     

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

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

Dynamics of an Elastic Cable Carrying a Moving Mass Particle

The dynamic behavior of an elastic catenary cable due to a moving mass alongits length is investigated. The equations of motions are derived using theHamilton's principle for general supports that include the horizontal andinclined cables with small and large sags and for variable velocity of themoving mass. Those equations of motions are in general nonlinear partialdifferential equations due to the initial curvature of the cable. Theequations are also complex due to the presence of three different types ofaccelerations of the moving mass. Those are the normal, Coriolis andcentrifugal accelerations. Therefore, we used the Galerkin procedure withsine function (Fourier representation) and anti-derivative functions of thecompactly supported wavelets as trial basis and used direct integrationmethods to integrate the discretized equations of motions. Newton–Raphsonmethod is used for iterations. Several examples are studied and the resultsas obtained by Fourier and wavelet representations are compared. Because ofthe localization feature, wavelets are proven to minimize the spuriousoscillations specially those appearing in the cable tension.     

Elastic cables–rigid body coupled dynamics: asymptotic modeling and analysis

An elastic cables–rigid body coupled model is proposed for investigating dynamic interactions between cables’ nonlinear transversal vibrations and boundary tower’s torsional dynamics, arising in large transmission line–tower systems and suspended cable–bridge tower systems. By introducing a weak torsion assumption and a large moment of inertia for the tower, an asymptotic expansion of cables–tower coupled dynamics is conducted in a weakly nonlinear framework, and a cables–tower reduced coupled model is eventually established. After model’s validations using direct numerical simulations, two distinct kinds of coupled dynamics are fully investigated. The first is that an external torque is applied to the tower and the two cables would both be indirectly excited, asymmetrically, by the torsional/oscillating tower. The two cables’ responses are the same in this case. The second is that only one of the two cables, i.e., the leader cable, is directly excited, and the other cable, i.e., the follower one, is only indirectly excited through cables–tower dynamic interactions. In such kind of leader–follower dynamics, the leader cable is quite different from the follower one. Nonlinear coupled frequency response diagrams for both systems are constructed using numerical continuation algorithms, mainly focused on the coupled steady solutions’ stabilities and bifurcations. Furthermore, the dynamic effects of tower’s moment of inertia, wing span and damping are thoroughly investigated.

     

索网天线的参变量变分及非线性有限元方法

针对大型周边桁架式索网天线由拉索拉压模量不同引起的本构非线性和结构大变形引起的几何非线性问题,给出了基于参变量变分原理的几何非线性有限元方法. 首先针对含预应力索单元拉压模量不同分段描述的本构关系,通过引入参变量,导出了基于参变量及其互补方程的统一描述形式,避免了传统算法需要根据当前变形对索单元张紧/松弛状态的预测,提高了算法收敛性. 然后利用拉格朗日应变描述索网天线结构大变形问题,结合几何非线性有限元法,建立了基于参变量的非线性平衡方程和线性互补方程;并给出了牛顿-拉斐逊迭代法与莱姆算法相结合的求解算法. 数值算例验证了本文提出的算法比传统算法具有更稳定的收敛性和更高的求解精度,特别适合于大型索网天线结构的高精度变形分析和预测.