基于AR法的输电塔线体系风速时程模拟
基于AR模型的风速时程模拟
关键词 ：Ａ 模型； Ｒ 风速时程； 空问相关性；ＯＴＡ ＦＲ ＲＮ
ＯＴＡ 随着 我 国经 济快 速发展 和 施工 工艺 的不 断进 步 ， 各 Ｆ Ｒ Ｒ Ｎ语 言 编 制 了模 拟 空 间 点脉 动 风 速 时 程 的计 算 种 大型 、复杂 的高层 结构和 大跨 空 间结 构不 断兴 建 ， 风 程 序 ，对 计 算 所 得 的风 速 时 程 进 行 快 速 傅 立 叶 变 换 ＦＴ， 以 荷载 已成 为结 构抗风 计算 、 防灾 减灾 分析 的重要 设计荷 （Ｆ ） 并把 获得 的功 率谱谱 与 目标功 率谱 进行 比较 ， 载， 结构 风振 响应分析 研 究 日益 受 到学术 界和 工程 界 的 验证基 于 Ａ Ｒ模 型风速 时程 模拟 的可行 性 。
＾
研究与探讨
式中，
个数 ） Ｌ １ ＋ Ｋ 一２ （ ３） ２
（ １ 式 ）
＾
１ × ｌ … １】 ｌ Ｍ Ｉ ， ， ｒ １， ｌＩ ｌ Ｐ ｒ Ｍ ｒ Ｐ Ｔ’
用。
５结论
【 考文 献】 参
［］ １ 王星 华 ． 土 固化 浆 液 在 地 下 工 程 中 的应 用 ［］ 北 京 ： 国铁 粘 Ｍ． 中
（） １ 粘土 固化 浆液 主要材 料是 粘土 ，可 以就 地取 材， 道 出版 社 ，９ ８ １９ ． 水泥 、 水玻 璃用量 少， 因此材 料成 本低 。 ［］ ２ 阮文 军 ． 液基 本 性 能 与 岩 体 裂 隙注 浆 扩 散 研 究 ［］２ ０ ． 浆 Ｄ．０３
后塑性 强度 增大， 既保证 了可 注性 ， 又不会 扩散 太远 。 而
［］ 吉林 工程技术师范学院学报， ０ ８ ２ （０ ：２ ６ Ｊ． ２ ０ ，４ １） ６ — ４
基于AR法的输电塔风振响应时程分析
动风速时程 ， 不需要针对特定结构 ， 因此 更 具 一 般 性 ． 借 助 于 计 算 机 的迅 速 普 及 和数 值 分 析 方 法 的深 入 研
究， 风 荷 载 的数 值 模 拟 理 论 研 究 成 果 丰 硕 ］ ， Ｉ ａ ｎ ｎ ｕ ｚ ｚ ｉ 等Ｅ ７ ］ 针 对 高耸 塔 桅 结 构 用不 同方 法 获 得 的 风 速 时 程 作 用 下 的风 振 相 应结 果 作 出 比较后 ，建 议采 用 线 性 滤 波 法 的 自回 归 （ Ａｕ ｔ ｏ Ｒ ｅ ｇ ｒ ｅ ｓ ｓ ｉ ｖ ｅ ， Ａ Ｒ） 方法＿ ８ ］ 模 拟 互 相关 的 多重 脉 动 风 速 时程 ．
摘 要 ： 采用 Ｍ Ａ ｒ Ｌ Ａ Ｂ编 制 了基 于 Ａ Ｒ 法 的风 荷 载 模 拟 程 序 ， 以广 西 某 工 程 实 际 ２ ２ ０ ｋ Ｖ 输 电塔 为 对 象 ． 用Ａ Ｎ Ｓ Ｙ Ｓ分 析 了 输 电塔 架 结 构 在 实 际 风 场 中不 同风 向角 的 振 动 响应 以及 不 同跨 度 下 时 程 分 析 和 规 范 计 算 两 者 １ ０ （ ２ ０ １ ３ ） ０ ２ — ０ ０ ７ ４ ． ０ ４
基 于 ＡＲ 法 的 输 电塔 风 振 响应 时程分 析
陈建 芳 ａ ＇ ｂ ， 陈家 豪 ， 何 志
（ 广 西 大学 ａ ． 土木 建 筑 工 程 学 院 ； ｂ ． 工 程 防灾 与 结 构 安 全 教 育 部 重 点 实 验 室 ， 广西 南 宁 ５ ３ ０ ０ ０ ４ ）
塔 架 结 构 的抗 风 设 计 研 究提 供 了必要 的技 术 参 数 ．
１ 塔 架计 算模 型
采用 ２ ２ ０ ｋ Ｖ广西主干电网中工程实例输 电塔架模型 ， 塔高 ５ ６ ． ２ｍ， 呼
基于AR法的沿海地区强台风风场模拟研究
等提 出的基于上海地 区实测台风风速 的风谱模型． 该模型实质上是根据台风特性， 对 Ｄ ａ ｖ ｅ ｎ ｐ ｏ ｔ风谱 ｒ
模型进行了修正．
３ 脉 动风场 模 拟方 法
响较大的台风： （ １ ） ２ ０ ０ ９ 年台风 “ 莫拉 克”，中心最
大风 力 １ ２级 ；（ ２ ） ２ ０ １ ０年 台风 “ 圆规 ”，中心最 大
“ ） ＝ ∑∑ （ ） ｌ ・
ｋ＝ｌ／ ＝１
式 中，
ｃ ｏ ｓ （ ｔ ｏ ｒ ｔ ＋九 ＋ （ ） ） ， Ｊ＝１ ， ２ ， …， ｍ， （ ３ ） ， ） 为互功率谱密度矩阵 ∞ ， ） 的Ｃ ｈ ｏ ｌ ｅ ｓ ｋ ｙ
分解矩阵 ／ （ ｃ ｏ ， ） 中的元素， 肼 为两个不 同作用作用 点之间的相位角， 为介于 ０ ～ ２ ７ ｃ 之间均匀分布的
风力 １ ３ 级； （ ３ ） ２ ０ １ １ 年台风 “ 梅花” ， 中心最大风力 １ ６ 级； （ ４ ） ２ ０ １ ２年台风 “ 海葵”， 中心最大风力 ｌ １ 级； （ ５ ） ２ ０ １ ３ 年台风 “ 菲特 ” ， 中心最大风力 １ ４ 级． 本 次研究选取 台风级别为强台风， 取 １ ０ ｍ 高度处平
２ ＿ ３ 石 沅谱
据 中国气象局规 定， 热带气旋按 中心附近 地 面最大风速 划分为六个等级 ， 其 中易产生大规模 破坏 的有 ： （ １ ） 台风：底层 中心附近最 大平均风速
３ ２ ． ７ ￣ ４ １ ． ４ ｎ － ｉ ・ ｓ ～ ， 即１ ２ ～ ｌ ３ 级； （ ２ ） 强台风： 底层中心
附 近最 大平 均风 速 ４ １ ． ５ ～ ５ ０ ． ９ １ Ｔ Ｉ ・ Ｓ ‘ 。 ，即 １ ４ ～ １ ５ 级； （ ３ ） 超 强 台 风 ：底 层 中 心 附 近 最 大 平 均 风 速 ≥５ １ ． ０
基于AR模型模拟超高层建筑的脉动风速时程
基于AR模型模拟超高层建筑的脉动风速时程李春祥;都敏;韩兵康【期刊名称】《地震工程与工程振动》【年(卷),期】2008(28)3【摘要】风荷载是超高层建筑设计的主要荷载之一,而且通过风振时域分析可以更全面地了解超高层建筑风振响应特性,更直观地反映超高层建筑风致振动控制的有效性。
因此,本文使用线性滤波法即白噪声滤波法(white noise filtration method,WNFM)中的自回归(auto-regressive,AR)模型模拟超高层建筑的风速时程。
首先,考虑超高层建筑风速时程的时间和空间相关性,导出了四阶AR模型的参数表达式。
接着,基于AR模型模拟了一幢高度为200 m超高层建筑的风速时程。
最后,通过比较模拟风速功率谱、模拟自相关函数和互相关函数与目标风速功率谱、目标自相关函数和互相关函数的吻合程度,验证基于AR模型模拟超高层建筑风速时程的可行性。
【总页数】8页(P87-94)【关键词】超高层建筑;风荷载;风速时程;自回归模型;随机过程;模拟【作者】李春祥;都敏;韩兵康【作者单位】上海大学土木工程系,上海200072;同济大学土木工程学院,上海200092【正文语种】中文【中图分类】TU311【相关文献】1.基于AR模型的空间脉动风速时程模拟方法研究 [J], 赵海霞2.基于ARMA模型模拟高架桥的脉动风速时程 [J], 李春祥;谈雅雅;李锦华3.基于AR模型模拟的转体桥梁脉动风数时程 [J], 宋浩4.基于AR模型的大跨悬索桥脉动风速时程模拟 [J], 白泉;徐樊;杨少波5.基于线性滤波法的超高层建筑脉动风速时程模拟 [J], 常乐;郭小飞因版权原因，仅展示原文概要，查看原文内容请购买。
基于AR模型方法的多维脉动风荷载时程的模拟
通过设 计好 的过滤 器 ，使 其输 出为具 有 给定谱 的 随机过 程 。该 法 占用 内存 少 ，计算 快捷 。近年来 ，线 性
滤波器 法 中 自 回归滑 动 平 均 模 型 （ ｕｏｅｒｓｉ ｖ ｖ ｒ ｅＭｏｅｓ Ａ ｔ ｇｅｓ ｅＭｏｅＡ ｅａ ｄ ｌ，简 记 为 Ａ Ｍ ｒ ｖ ｇ Ｒ Ａ） 和 自回归 模 型 （ ｕｏｅｒｓｉ ｏ ｅ ，简记 为 Ａ Ａ ｔｒｇ ｓ ｅＭ ｄｌ ｅ ｖ ｓ Ｒ） 被广 泛用 于描 述平稳 随机 过程 ，取得 了 良好 结果 。
Ⅳ（） ＝Ｃ ‘ （ ） ｔ ｎｔ
（） ２
式 中，（）＝ ［ ｔ ，：ｔ ， ，ｎｆ ］ ｎ（）是 均 值为 ０ 方差 为 １且彼 此相 互 独立 的正 态 随机过 ｎｔ ｎ （） ｎ（） … ａ（） ，ｉｆ 、 程 ， １２ … ， ； ｉ＝ ，， ｎＣ为 ｎ阶下 三角矩 阵 ， 通过 ｎｘｎ阶协方差 矩 阵 Ｒ 的乔里斯 基 （ ｈｌｓｙ 解确 定 ： Ｃ ｏｅｋ ）分
收稿 日期 ： ０ ７—０ ２０ ９—０ ９ 作者 简 介 ： 建 平 （ ９ ２一） 王 １６ ，男 ，副 教 授 ，研 究 方 向 ：结 构非 线 性 分 析 。
维普资讯
第 ５期
王建 平 ： 于 Ａ 基 Ｒ模 型方法 的多 维脉动风荷载时程的模拟
ｐ
（ ，）＝ ｚｔ
ｕ ｔ Ａ ）＋Ｎ（） （ —ｋ ｔ ｔ
（） １
式 中 ， 间第 点坐标 ， 空 ＝１２ … ，； Ａ ，， ｎｐ为 Ｒ模 型阶数 ；ｔ Ａ 是模 拟 风速 时程 的时 间步长 ； 为 Ａ Ｒ模 型
风电机组塔架的脉动风速时程模拟
风电机组塔架的脉动风速时程模拟叶赟;宫兆宇【摘要】In this paper, an autoregressive model(AR model) is used to simulate wind speed time series. e spectrum of simulated wind speed time series is found in agreement with the target spectrum, Davenport wind speed spectrum. Samples of the uctuating wind load on the nodes of a structure are obtained. Using the WAWS, the article builds an AR model to calculate the model order and edit a simulation program. rough the analysis on some wind turbines tower, the feasibility and e ciency of this simulation model is veri ed.% 本文简述了谐波合成法中的自回归模型(AR)模拟出给定风速功率谱的风速时程序列，并验证其与目标谱(Davenport谱)的一致性，从而得到作用在各节点的脉动风荷载时程样本的方法。
本文采用谐波合成法，建立了脉动风速时程的 AR 模型，编辑出脉动风速时程模拟程序，并对某风电机组塔架进行脉动时程分析，验证了该脉动风速时程模拟的可行性与有效性。
【期刊名称】《风能》【年(卷),期】2013(000)001【总页数】6页(P72-77)【关键词】脉动风;数值模拟;自回归模型;风电机组塔架【作者】叶赟;宫兆宇【作者单位】内蒙古科技大学建筑与土木工程学院，包头 014010;内蒙古科技大学建筑与土木工程学院，包头 014010【正文语种】中文【中图分类】TM614高耸结构风荷载是结构设计时必须要考虑的一类重要的随机荷载，风振响应成为控制结构设计的重要因素。
基于AR线性滤波法脉动风模拟
基于AR线性滤波法脉动风模拟随着我国大跨度桥梁的出现，风荷载必然是其主要考虑对象。
风荷载分为平均风和脉动风，然而在研究分析桥梁一系列的震动问题上，脉动风又起主要控制作用。
有些桥址因气象资料不足，给分析研究带来了很多困难，以重庆朝天门长江大桥桥面为研究对象采用线性滤波法模拟脉动风。
标签：风场特性；大跨度桥梁；脉动风模拟1、引言我国多数地区以山地为主的国家，尤其是在西南片区。
中国拥有960万平方公里的国土面积，其中有近65%的山地。
我们要想在这种山地形修建公路，必然得修建桥梁。
要想研究向引起对桥梁的风致震动，必须研究清楚该桥址的风场特性。
我国中西部地区尤其是西南地区的大跨度桥梁通常位于高山峡谷的山区地带，然而这些地区通常缺乏历史气象记录，除此之外山区峡谷的阵风强烈，风速变化较为频繁，湍流强度较大，非平稳特性较为突出，桥址区风速场空间分布不但具有三维特征还很复杂，同跨海和跨江大桥桥址区的风场特性相比有着明显差异。
如果采用基于相关规范和规定的风特性进行分析，而不考慮桥址区地形的具体特征，将导致相当大的误差，甚至使计算结果发生质的变化。
由于风在山区峡谷地形条件下较为复杂，不同环境条件下桥址区的风场特性相差较大，没有类似风场可借鉴。
目前业界相关研究对复杂地形地貌区强风分布还少之又少，而风参数的不确定性是影响桥梁结构抗风研究精度最重要的因素[1]。
重庆朝天门长江大桥连接重庆三大中央商务区：弹子石、江北城、解放碑，大桥位于溉澜溪青草坪。
该桥为钢桁架拱桥形式，其主跨达552m，具有世界第一钢拱桥的称号。
大桥西接江北区五里店立交，东接南岸区渝黔高速公路黄桷湾立交，全长1741米，是主城一条东西向快速干道。
其重庆朝天门长江大桥周围环境如图1所示，由重庆交通科研设计院提供的主桁立面结构图，假设水面距桥墩顶面为20m，则主梁的离地高度为60米。
2、AR线性滤波法谐波合成法，主要是基于三角级数求和，以散谱逼近目标随机过程的随机模拟，其最大不足就是在运算时准备就绪比较低，对时间相关性的模拟是没有考虑在其中的，并且我们在进行模拟时每个频率都要做大量的运算[2]。
风速时程AR模型及其快速实现
R ( j ∃ t) = -
∑7
k= 1
k
R [ (j -
k ) ∃ t ], j = 1, …, p
( 7)
p
R ( 0) = -
∑7
k= 1
k
R ( k ∃ t) + R N
( 8)
写成 A R 模型的正则方程, 即:
R ・7 = RN Op
( 9)
V (X , Y , Z , t) = -
R 12 ( ∃ t) R 22 ( 0) R 32 ( ∃ t) R
( p + 1) 2
R 13 ( 2∃ t) R 23 ( ∃ t ) R 33 ( 0) R
…
R 1 (p + 1) ( p ∃ t)
… R 2 (p + 1) [ ( p - 1) ∃ t ] … R 3 (p + 1) [ ( p - 2) ∃ t ] ω …
R 11 ( 0) R 21 ( ∃ t) R= R R 31 ( 2∃ t)
( p + 1) 1
式中, 7 = [ I , 7 1 , …, 7 p ]T , 为 ( p + 1 ) M ×M 阶矩 阵, I 是 M 阶单位阵; R N 意义同前; O p 为 pM ×M 阶矩阵, 其元素全部为 0; R 为 ( p + 1)M × ( p + 1)M 阶自相关 Toep litz 矩阵, 写成分块矩阵的形式, 即:
p
( 4) j ∃ t) + N ( t) V T ( t - j ∃ t) ( 式中, j = 0, …, p. 作 数 学 期 望 运 算 Exp ecta t ion [ ] ) 并结合自相关函数的如下性质: T ( 5) R ( - j ∃ t) = E [V ( t) V ( t- j ∃ t) ] ( 6) R ( - j ∃ t ) = R ( j ∃ t) 得两个方程 ( 组) , 即:
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
i
3 . 2 风速时程分析 ( 1) 点 1 的脉动风速时程曲线见图 2、 风速模 拟谱与目标谱拟合曲线见图 3 。由图可看出, 采 用 AR 法编 制 程 序 模 拟 的脉 动 风 速 谱 与采 用 Kaim al 谱计算获得的目标谱拟合效果好。
( 6)
T
为 M ! M 阶方阵;
RN 为 M ! M 阶方阵; Op 为 p - 1 M ! M 阶矩 阵; 其元素全部为 0。 AR 模 型阶 数 根据 最 小 AIC 准则 确 定[ 6] 。 AIC 函数为: AIC ( p ) = N lg 其中
收稿日期 : 2010 09 16, 修回日期 : 2010 10 12
2 基于线性滤波器法的脉动风模拟
空间中 M 个点相关脉动风速时程 V( X, Y, Z, t) 列向量的 AR 模型可表示为:
p
V X, Y, Z, t = k= 1
k
V X, Y, Z, t - k t + N t ( 4)
基金项目 : 东北电力大学研究生创新基 金资助项目 ( 东电研字 200921) 作者简介 : 秦力 ( 1970 ) , 男 , 教授 , 研究方向为输电塔结构可靠性分析 , E mail: jilinql@ 163. co m
表 1 提取 风速点坐标
式中 , ( x i , y i , z i ) 为空间第 i 点坐标, i = 1, 2, 长; 为 AR 模 型自回归系数 矩阵, k = 1 , 2,
, ,
1 6 14 26 30 78
Ta b. 1 Co o rdina t e o f e xt rac t e d point s
& 170 & 其中 X = [ x 1, x 2, Y = [ y1 , y 2 , Z = [ z 1, z 2, , xM] T , yM] , zM]
T T
水
电
能
源
科
学
2011 年
10 H z; % 模型参数。节点设置总数为 78 个 , 计算 模型阶数 p = 4, 表 1 为模型部分节点坐标值。
2 ! 2 !
图2
点 1 脉动风速时程曲线
+ 2( p + 1)
( 7)
Fig . 2 T ime his t o ry curv e o f f luc t ua t ing w ind spe ed o f po int 1
= 2R 0 - R N
式中, N 为样本容量。从一阶模型开始求 AIC ( p ) 的函数值, 直至找到使其最小的 p 为止, 一般取 4~ 5 阶即可满足要求。
1 脉动风的基本特性
脉动风荷载为随机荷载[ 4] , 是风荷载中的动 力部分, 其振动速度和方向具有随时间和空间随 机变化的特点, 可采用脉动自动率谱、 互动率谱描 述脉动风速时程。 1 . 1 脉动风自功率谱 设计中 较常 用的 脉 动风 速谱 为 Davenport 谱, 以 10 m 高度处的风速为基准风速 , 且风速谱 不随高度变化。输电塔线体系作为高耸结构 , 用
Fig . 3
Fit t e d c urv e of s imula t e d spe ct rum a nd t arg e t s pe c t rum o f po int 1
( 2) 点 1 、 6、 14 脉动风速时程曲线比较。为便 于比较 , 将 点 6、 14 的 脉动 风速值 分别加 20、 40 m/ s, 比较结果见图 4 。由图可看出: # 不同高度 处脉动风速变化趋势相同 , 但各时刻的速度不同, 表明脉动风速具有随机性 ; ∃ 随高度增大, 平均风 速变大, 但脉动风的波动区间变小 , 表明输电塔线
3 算例
图 3 点 1 风速模拟谱与目标 谱拟合曲线
3 . 1 风速时程模型 500 kV 栖霞 ∀ 文登 ( 昆嵛 ) 送电工程直线塔 为 5D SZ 1 双回路直线塔, 塔高 66. 4 m, 档距 500 m, 建立三塔四线模型见图 1。基于 Mat lab 软件 编制脉动风速时程模拟程序, 各参数分别为: # 基 本参数。根据文献[ 7] 求得标准高度 ( 10 m ) 处平 均风速为 v 10 = 29. 665 m/ s, 地面粗糙度系数 k = 0. 005; ∃ 时间 和 频 率参 数。时 间 步长 0. 1 s, 时程总长 t = 300 s, 初始频 率 0. 01H z, 截止频率
m
M; p 为 AR 模型阶 数; t 为模拟风 速的时间步
k
p ; N t 为独立随机过程向量。 根据 风 速 时 程 假 定 , 式 ( 4 ) 两 边 同 时 乘 VT X, Y, Z, t- j t , 并求数学期望有 :
p
R j t =k= 1
k
R
j- k
t
( j = 1, 2 , , p ) ( 5) 式中 , R 为 p M ! p M 阶自相关 T oeplitz 矩阵。 则 AR 模型的正则方程为 : RN R = Op 其中 = 1, 2, , p 式中 , 为 p M ! M 阶矩阵 ,
b. 脉动风具有空间相关性, 两点之间距离越 近相关性越强 , 距离越远相关性越弱, 且水平距离 变化时相位亦有差别。 c. 采用 AR 法模拟输电塔线体系的脉动风速 时程 , 结果表明该方法可行, 在选择适当的时间步 长方面值得进一步研究。 参考文献 :
[ 1] 国家电力公司华东电 力设计 院 , 国 家电力 公司电 力 规划设计 总院 . 110~ 500 kV 架空 送电线 路设计 技 术规程 ( D L/ T 5092 1999P ) [ S] . 北 京 : 中 国电 力 出 版社 , 1999. [ 2] [ 3] [ 4] [ 5] [ 6] 王 之 宏 . 风 荷 载 的模 拟 研 究 [ J] . 建 筑 结 构 学 报 , 1994, 15( 1) : 44 52. 刘锡良 , 周颖 . 风 荷载 的几 种 模拟 方法 [ J] . 工 业 建 筑 , 2005, 35( 5) : 81 84. Dy rbye C, H ansen S O. 结构风荷载作用 [ M ] . 薛 素 铎 , 李雄彦 , 译 . 北京 : 中国建筑工业出版社 , 2006. 黄本 才 , 汪 从 军 . 结 构 抗风 分 析原 理 及应 用 ( 第 二 版 ) [ M ] . 上海 : 同 济大学出版社 , 2008. 马文平 , 李 兵兵 , 田 红心 , 等 . 随机 信号 分析 与应 用 [ M ] . 北 京 : 科学出版社 , 2006. [ 7] 中国建筑 科学 研究 院 , 同 济大 学 , 建设 部建 筑设 计 院 , 等 . 建筑 结构 荷载 规 范 ( G B50009 2001) [ S] . 北 京 : 中国建 筑工业出版社 , 2002.
模拟点高度; n 为脉动风频率; v z 为 z 高度处的平 均风速; k 为地面粗糙度系数 ; v 10 为标准高度 10 m 处的平均风速。 1 . 2 脉动风互功率谱 观测表明 , 脉动风并非完全同步, 与风速及各 点的相对位置有关 , 因此需考虑互功率谱相干函 数。互功率谱数学表达式为: S ij f = S ii f S j j f r ij f ( 2) 式中, S ii f 、 S j j f 分别为 i 、 j 点的自谱密度函 数; r ij f 为 i 、 j 点的相干函数, 根据 Shiotani 的建 议取值为 : r ij x i , x j , y i , y j , z i , z j = exp{- [ ( x j - x i ) 2 / L x + ( y j - y i ) 2 / L y + ( z j - z i ) 2 / L z ] } 1/ 2 ( 3) 其中 L x = L y = 50; L z = 60 式中, x 、 y、 z 为模拟风速 点的空间坐标 ; L x 、 L y、 L z 分别为考虑风速空间各方向相关性的系数。
目前 , 在输电线路设计过程中 , 风荷载一般按 静力风考虑 。但输电塔线体系中的铁塔、 导线、 绝缘子之间的动力耦合作用强烈, 脉动风对塔线 耦合体系的作用不能忽略。通常采用风洞试验和 现场实测的方法确定风荷载, 但耗资大、 周期长、 适用性有限。因此, 脉动风速时程的计算机模拟 便具有重要意义。目前 , 常用脉动风速时程获得 风荷载时程, 模拟方法主要有谐波叠加法、 线性滤 [ 2, 3] 波法两种 , 均基于蒙特卡洛法, 将脉动风速谱 模拟成脉动风速时程 , 再在准定常假设的基础上 将风速时程换为风荷载过程。它们从单一脉动风 速时程模拟发展到多个相关风速时程 , 各有其优 缺点。其中线性滤波法中的自回归法 ( AR) 模型 因速度快、 计算量小 , 已广泛应用于随机振动的时 域模拟。由于输电塔线体系模型大、 自由度多 , 风 速时程模拟时计算量极大。鉴此, 本文采用线性 滤波法模拟输电塔线体系风速时程, 并以 500 kV 栖霞 文登( 昆嵛) 送电工程为例进行了模拟分析。
第 29 卷第 2 期
秦
力等 : 基于 A R 法的输电塔线体系风速 时程模拟
& 171 &
体系脉动风振作用随高度增加而减弱。 ( 3) 点 26、 30 、 78 脉动风速时程曲线比较。为 便于比较, 将点 30、 78 的脉动风速值分别加 20、 40 m/ s, 比较结果见图 5 。由图可看出, 在高度不 变时 , 各点的平均风速相同 , 但脉动风速不同 , 相 位有差异 , 表明脉动风具有空间相关性。