一种新的覆冰导线舞动非线性有限元分析方法

本文基于Hamilton原理,建立了耦合三个平动自由度和一个扭转自由度以及偏心覆冰作用下导线的非线性舞动方程,提出了一种新的覆冰导线舞动分析非线性有限元模型,将相邻跨导线和绝缘子等效为线性弹簧单元,同时考虑覆冰导线非线性气动力和大幅舞动的几何非线性,采用Newmark-β时间积分结合修正的Newton-Raphson非线性迭代求解舞动有限元方程。本文数值解与D形覆冰导线舞动实测值无论在振幅还是频率方面均吻合得非常好,证明了本文新提出的导线舞动有限元模型的准确性。本文研究还表明：舞动是一种以上下运动为主的低频振动,通常发生在一阶上下自振频率附近,其振幅和频率由输电线路的物理参数和风载荷唯一确定,与初始运动状态无关。     

GALLOPING BEHAVIORS OF SECTOR-SHAPE ICED EIGHT BUNDLE CONDUCTORS1)

考虑到重冰区的特高压输电线路更容易形成接近扇形的覆冰,结合风洞实验及数值模拟方法研究了扇形覆冰八分裂导线的舞动特征。通过扇形覆冰八分裂输电线的节段模型风洞实验,获得各覆冰子导线随风攻角变化的气动参数。在ABAQUS软件中建立单档扇形覆冰八分裂输电线路的有限元模型,通过程序UEL输入各覆冰子导线的气动参数,根据数值模拟结果可获得舞动轨迹及振幅等舞动特征,最后讨论了风速、档距和初始风攻角对八分裂导线舞动特征的影响。     

考虑温度效应的覆冰导线动力学建模及舞动特征研究

针对以往研究忽略了温度效应对覆冰导线舞动特性的影响,本文推导了考虑温度效应影响的覆冰导线舞动控制方程.基于悬链法、热应力理论推导了覆冰导线的偏微分舞动方程,接着通过Galerkin法将该偏微分方程转化为常微分方程.建立气动载荷模型,将气动力引入到舞动方程中,随后采用多尺度求得了覆冰导线的位移响应,最后进行了参数分析、算...     

AUTOMATED WIND TUNNEL MEASUREMENTS

It is common practice in power transmission to use bundled overhead electrical conductors because the lower rate of corona discharge, compared with single conductors, makes the bundle attractive for high voltages. However, secondary problems may arise from the vibrations induced by the aerodynamic forces which act on the leeward conductors as a result of their immersion in the wake of a windward conductor. Thus, in addition to possible aeolian and galloping oscillations observed for a single conductor, bundle conductors may also be subjected to wake-induced oscillations. ^1,2 Although extensive research has been done on the aeroelastic behavior of single and bundled conductors which are bare, little is known about the galloping or wake oscillation characteristics of given bundle conductors which, depending upon the environmental conditions during the formation of ice, may have five or so representative shaped cross-sections. One important task involves the aerodynamic quasi-static force and moment measurements in a wind tunnel using full scale conductor cross-sections each of which must be rotated through many angles at different spacings. Such measurements are repetitive as well as tediously time consuming and, hence, are beneficially automated. A reliable, accurate and inexpensive automated apparatus is described in this paper which permits the aerodynamic measurements to be performed efficiently and without intervention. Most importantly, it relieves the operator of the usual responsibility for manually controlling the separation and orientation of the sections as well as the need to manually compensate fluctuations in the free stream air speed. For simplicity, the apparatus is detailed in the context of a twin bundle conductor.     

连续体覆冰导线模型的舞动实验研究

建立了适合连续体覆冰导线舞动实验的专用风洞,以连续体单跨覆冰单导线模型为实验对象,采用激光传感器测量了导线不同位置处的位移响应,得到了导线在不同风速下的舞动振型,利用力传感器得到了导线舞动时的动张力。结果表明:覆冰导线在来流风场作用下进行舞动,在随风速增大的过程中先后经历了两个大幅舞动阶段,一阶模态和二阶模态在两个阶段中分别被激发;动张力幅值与舞动幅值基本成正比,动张力频率包含导线模型舞动频率的一倍频和二倍频。在大风速时,导线模型的舞动为多阶模态的耦合振动,在舞动过程中无固定的波峰和结点。     

连续档覆冰导线舞动数值模拟及参数分析简

本文对仅有一档满足舞动条件的连续档输电线路的振动进行了详细地研究. 利用ABAQUS 通用软件建立考虑绝缘子串连续档覆冰输电线的精细有限元模型,通过用户自定义单元UEL 施加空气动力载荷,获得考虑几何非线性的连续档覆冰导线大幅振动的时域解. 根据数值结果研究了连续档输电导线档距及档数变化对连续档内不同档输电线舞动的影响规律,所得结果对连续档覆冰导线舞动和防舞技术的研究具有一定的参考价值.     

构造杂交应力单元的柔度矩阵H对角化方法

证明了杂交元柔度矩阵 H非奇异的充分必要条件是假设应力模式线性无关 ;以及等价应力模式形成相同的杂交元。在此基础上建立了假设应力模式的 Hilbert子空间 ,从而可以利用 Schmidt方法简单地得到等价的正交应力模式 ,实现了柔度矩阵 H对角化 ,使得杂交元形成过程中完全避免了繁杂的矩阵求逆运算 ,提高了杂交元分析的计算效率 ,特别在柔度矩阵不容易显式求逆的材料非线性分析中更具有实际意义     

Calculation of magnetic fields and currents in axisymmetric systems of inductively coupled moving conductors

An algorithm for calculating magnetic fields and currents in axisymmetric systems of inductively coupled moving and stationary conductors is developed using a hybrid method which combines the finite and boundary element methods. The finite element method is used to approximate the unsteady diffusion equation for the θ-component of the magnetic vector potential in the conductors, and the boundary-element method is employed to eliminate the space around the conductors. The proposed method takes into account the connections of the conductors with each other or with external energy sources by means of ideal electrical circuits with lumped parameters R, L, and C. An effective method is developed to take into account the external circuits by an appropriate modification of the mass matrix and the source vector of the obtained system of ordinary differential equations. Examples of using the method to calculate the fields of single- and multi-turn solenoids, magnetic flux concentrators, and induction accelerators with various methods of delivering external electromagnetic energy are considered. The high computational efficiency of the method is shown, in particular, for the case of constant electrothermal properties and sizes of the conductors. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 1, pp. 22–29, January–February, 2009.     

杂交元本征应力模式和应力子空间的性质研究

详细讨论了有限元本征应力模式和应力子空间的性质，并着重讨论和进一步完善了与杂交应力有限元应力子空间有关的一些定理，为提出新方法提供了理论基础，主要包括：（1）证明了杂交元特征值不大于对应位移元的特征值；（2）证明了矩阵H非奇异的充分必要条件是假设应力模式线性无关；（3）证明了杂交元所对应位移元的本征应力模式形成的杂交元与该位移元相同；（4）证明了等价假设应力模式形成相同的杂交元；（5）证明了确定杂交元本征应力模式的充分必要条件是其范数平方等于所形成杂交元的变形模态特征值；（6）证明了杂交元假设应力模式与变形模态的能量一一对应的充分必要条件是假设应力模式彼此正交且与所对应位移元的本征应力模式除了一一对应者之外都正交。     

Control of cross-flow-induced vibrations of square cylinders using linear and nonlinear delayed feedbacks

We investigate the effectiveness of linear and nonlinear time-delay feedback controls to suppress high amplitude oscillations of an elastically mounted square cylinder undergoing galloping oscillations. A representative model that couples the transverse displacement and the aerodynamic force is used. The quasi-steady approximation is used to model the galloping force. A linear analysis is performed to investigate the effect of linear time-delay controls on the onset speed of galloping and natural frequencies. It is demonstrated that a linear time-delay control can be used to delay the onset speed of galloping. The normal form of the Hopf bifurcation is then derived to characterize the type of the instability (supercritical or subcritical) and to determine the effects of the linear and nonlinear time-delay parameters on their outputs near the bifurcation. The results show that the nonlinear time-delay control can be efficiently implemented to significantly reduce the galloping amplitude and suppress any dangerous behavior by converting any subcritical Hopf bifurcation into a supercritical one.     

基于3自由度的新月形覆冰输电线舞动稳定性研究

针对覆冰输电线舞动问题提出了一种基于非对称空气动力系数矩阵的临界风速计算方法．基于拟静态理论得到覆冰输电线的气动载荷,该气动载荷考虑了横向运动以及扭转运动对相对风攻角的影响,最后建立等效的3自由度覆冰输电线舞动模型．在初始风攻角处对气动载荷进行泰勒展开,得到非对称的线性空气动力系数矩阵．结合3自由度振动方程以及非对称空气动力系数矩阵,采用Rourh-Hurwitz准则计算覆冰输电线舞动发生的临界风速．通过风洞实验测得新月型覆冰单导线的空气动力系数,根据本文提出的理论分析了竖向振动频率、面外振动频率以及扭转振动频率对临界风速的影响,最后与DenHartog理论得到的临界风速进行了对比．本文研究成果对于指导覆冰输电线路防舞设计具有理论意义．     

温度效应对覆冰导线舞动特征的影响研究

以某耐张段单跨导线作为研究对象,基于哈密顿原理和增量热场理论,建立了考虑温度效应的覆冰导线非线性舞动模型并推导出其舞动方程。分别使用近似模态和复杂模态得到离散后的面内-面外耦合的有限维度舞动方程,接着采用多尺度法得到该方程的近似解。数值算例研究表明,温度对频率的影响与Irvine参数有关,对张力的影响与初始张力的范围有关,而且降温会增大导线舞动幅值;采用精确模态计算得到的舞动幅值显著小于使用近似模态计算得到的幅值。     

THE STRESS SUBSPACE OF HYBRID STRESS ELEMENT AND THE DIAGONALIZATION METHOD FOR FLEXIBILITY MATRIX H

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｆｉｎｉｔｅｅｌｅｍｅｎｔａｎａｌｙｓｉｓ,ｔｈｅｄｉｓｐｌａｃｅｍｅｎｔｅｌｅｍｅｎｔｓａｒｅｖｅｒｙｐｏｐｕｌａｒｂｅｃａｕｓｅｏｆｔｈｅｉｒｅａｓｙｃｏｎｓｔｒｕｃｔｉｏｎ .Ｈｏｗｅｖｅｒ,ｉｔｉｓｖｅｒｙｄｉｆｆｉｃｕｌｔｔｏａｓｓｕｍｅｔｈｅｃｏｎｔｉｎｕｏｕｓｄｉｓｐｌａｃｅｍｅｎｔｆｉｅｌｄｆｏｒｓａｔｉｓｆａｃｔｏｒｙｓｈｅｌｌｅｌｅｍｅｎｔｓｔｈａｔｉｓｆｒｅｅｆｒｏｍｓｈｅａｒｌｏｃｋｉｎｇ .Ｉｎ 1 964 ,ｔｈｅｈｙｂｒｉｄｓｔｒｅｓｓｅｌｅｍｅｎｔｗ…     

Orthogonal basic deformation mode method for zero-energy mode suppression of hybrid stress elements

A set of basic deformation modes for hybrid stress finite elements are directly derived from the element displacement field. Subsequently, by employing the so-called united orthogonal conditions, a new orthogonalization method is proposed. The resulting orthogonal basic deformation modes exhibit simple and clear physical meanings. In addition, they do not involve any material parameters, and thus can be efficiently used to examine the element performance and serve as a unified tool to assess different hybrid elements. Thereafter, a convenient approach for the identification of spurious zero-energy modes is presented using the positive definiteness property of a flexibility matrix. Moreover, based on the orthogonality relationship between the given initial stress modes and the orthogonal basic deformation modes, an alternative method of assumed stress modes to formulate a hybrid element free of spurious modes is discussed. It is found that the orthogonality of the basic deformation modes is the sufficient and necessary condition for the suppression of spurious zero-energy modes. Numerical examples of 2D 4-node quadrilateral elements and 3D 8-node hexahedral elements are illustrated in detail to demonstrate the efficiency of the proposed orthogonal basic deformation mode method.     

开口薄壁杆件结构稳定分析的精确单元和两步求解算法

从控制微分方程的通解出发,构造受偏心压力作用开口薄壁杆件的精确形函数,建立用于开口薄壁杆件结构稳定性分析的精确有限元,得到了单元刚度矩阵和几何刚度矩阵的显式表达,提出了计算给定区间内各阶临界荷载以及相应失稳模态的两步计算方法。计算结果表明,与常规单元相比,采用精确单元无需进行网格细分就可以获得精确的数值结果,结合本文的两步求解算法,可以准确获得给定区间内全部临界荷载和失稳模态。     

变截面Timoshenko梁的单元刚度矩阵

变截面构件在工程中应用广泛,在对变截面梁进行数值计算时,需要建立变截面梁单元的刚度矩阵。该文采用势能驻值原理,考虑了轴力引起的几何非线性和剪切变形的影响,将梁截面刚度的变化率作为小量,得到了近似到二阶的单元刚度矩阵。在构造位移模式时,从梁的微分平衡方程出发,得到同样近似到二阶、分别以三次和五次多项式表示的剪切和弯曲位移模式。该文还证明了单元刚度矩阵的奇异性,给出了轴压刚度的表达式,定量论证了与某些精确解的误差,表明在一定范围内,该文的结果具有足够的精度。最后以一个计算实例说明该文的单元刚度矩阵具有较快的收敛性。     

Periodic and quasiperiodic galloping of a wind-excited tower under external excitation

The galloping of tall structures excited by steady and unsteady wind may be periodic or quasiperiodic (QP) with amplitudes having the same order of magnitude. While the onset of periodic and QP galloping was studied, their control on the other hand has received less attention. In this paper, we conduct analytical study on the effect of a fast harmonic excitation on the onset of periodic and QP galloping in the presence of steady and unsteady wind. We consider the cases where the unsteady wind activates either external excitation, parametric one or both. A perturbation analysis is performed to obtain close expressions of QP solution and the corresponding modulation envelopes. We show that at various loading situations, the periodic and QP galloping onset is significantly influenced by the amplitude of the fast external excitation. In the case where the unsteady wind activates parametric excitation, the QP galloping occurs with higher frequency modulation compared to the case where the unsteady wind activates external excitation. In the case where external and parametric excitations are activated simultaneously, fast harmonic excitation eliminates bistability in the amplitude response and gives rise to a new small QP modulation envelope.     

ON THE AEROELASTIC BEHAVIOUR OF RECTANGULAR CYLINDERS IN CROSS-FLOW

An experimental and numerical study of the aeroelastic behaviour of elongated rectangular and square cylinders is presented. The main results are for a rectangular section with an aspect ratio of 2. The experiments were performed with a flexible cylinder clamped at both ends. This configuration leads to unusual lock-in of the vortex shedding with different bending modes, although the final steady oscillations occur in the fundamental mode. The galloping regime is also investigated, and the effect of free-stream turbulence intensity. Critical velocities are detected which do not correspond to calculations using the quasi-steady theory. A simple modelling of galloping is proposed to better fit the experiments, but it is shown that some of the configurations, in turbulent flow, are probably interacting with the vortex shedding and make the modelling inefficient. Numerical simulations on a 2-D rectangular section are presented and the resulting wall pressure distributions are analysed using the proper orthogonal decomposition technique. Indicators are proposed in order to link the proper functions with their contribution to the aerodynamic force components, and then a classification of the proper shapes of the decomposition is done. It is shown by comparison between the static case and forced oscillations, in the galloping range, that secondary vortices inside the shear layer become symmetrical and their effect on the forces is cancelled.     

A numerical study of rotational and transverse galloping rectangular bodies

In this paper we investigate rotational and translational galloping instabilities due to fluid/structure interaction using a previously developed algorithm. This numerical technique utilizes a two-dimensional spectral/hp element method and a frame of reference transformation to ensure efficient computations. Both transverse and rotational motion of rectangular sections of varying aspect ratio are simulated for a Reynolds number of 250 and at reduced velocities which promote a galloping response. Qualitative comparisons with quasi-steady theory and experimental data are found to be favourable.     

大跨屋盖结构风振响应主要参振模态确定方法研究

合理选取主要参振模态是保证大跨屋盖结构风振响应计算精度的关键问题。本文以模态空间分布与风荷载空间分布相关性强弱作为选取主要参振模态的依据,考虑到屋盖结构表面风荷载分布的不确定性,提出了通过低阶主要贡献模态间接寻找风荷载强相关高阶模态的思想,首先实现对高阶参振模态数量的缩减。在此基础上,结合低阶参振模态,通过对模态响应方差矩阵的简化处理,构造了其等效矩阵,根据等效矩阵中对角线元素定义了模态参与系数,通过该系数能够方便的实现对主要参振模态的选取。最后,通过国家体育场屋盖主结构的风振响应分析对所提出方法的有效性进行了验证。