基于ALE法的输电塔-风流固耦合振动特性研究.pdf

东 北 电 力 大 学 学 报 第29卷第 6期 JournalOfNortheastDianliUniversity Vo1．29．No．6 2009年 l2月 NaturaIScienceEl‘ition Dec．，2009 文章编号：1005—2992(2009)06—0001—05 基于 ALE法的输 电塔 一风流 固耦合振动特性研究 裴延学，祝 贺 (东北电力大学 建筑T程学院，吉林 吉林 132012) 摘 要：基于ALE(拉格朗 日一欧拉)法推导了结构一风场联合作用的流固耦合方程，并对结构的 运动方程与流体控制方程进行了空间离散及迭代求解。应用此方法，分析输电塔 一风耦联 自由振动特 性 ，得到在耦合场作用下输电塔结构的振动特性 ，计算结果表明了本文采用方法求解非线性耦合振动问 题 的有效性 。 关 键 词：ALE法；输 电塔；风；流固耦合；振动 中图分类号：TM753 文献标识码 ：A 风对输电塔结构的作用属于流固耦合力学的研究范畴，其特点是风(流体)与输电塔(固体)之间的 相互作用，输电塔在风的作用下发生振动，其动力特性为位移和结构变形，而输 电塔又会反过来影响风 流场，改变风荷的风向及大小，严重时将发生风致灾变。 国内外的研究多进行风洞试验及数值模拟，鉴于数值模拟方法的实用性，在本文中应用 ALE法进 行流固耦合计算 ，推导了流体、固体运动控制方程及流固耦合方程，并介绍ALE法通过约束加速度及速 度实现其算法，最后通过算例分析证明本方法的有效性。 1 流体控制方程 非静止的不可压缩 Navier—Stocks流体控制方程可描述为 卜 ： +u·2v ·Vu一2v ·re(Ⅱ)+ P=b V ·u = 0 式中，“为网格速度 ，为物质速度 边界条件和初始条件为： =一pl+2 () (u)=÷( +(Vu)) 2 结构运动方程 结构运动方程可统一表达为如下形式，并采用中心差分法进行求解。 Mb +cf+kX = F (2) 收稿 日期 ：2009—11—12 基金项 目：中国电机工程学会电力青年科技创新项 目资助2008—16，吉林省教育厅”十一五”重点科研项 目资助 2007—51 作者简介：裴延学(1961一)，吉林市人，从事工程管理研究工作． 2 东北电力大学学报 第29誊 中心差分法基于用有限差分代替位移对时间的导数，若采用等时间步长，At=At，则加速度和速度 的中心差分近似为： +1一 一1 ． +l一2v 十／)i — — — 广 - 1 (3) N (3)式代A (2)式，则有： +c +kX (4) 3 ALE法基本原理 ALE法引入拉格朗13和欧拉坐标之外的第三个任意参考坐标，如下式 ： Xi，t)／Ot= ，t)／Ot+ of( ，t)／Ox (5) 式中，五为拉格朗13坐标， 为欧拉坐标，i为相对速度。 ‘ 连续方程：