适用于独立变桨控制方式的风机气动模型和风模型的研究.pdf

水 力 发 电 第 3 7 卷第 4 期 2 0 1 1年 4月 适 用 于 独 立 变 桨 控 制 方 式 的 风 机 气 动 模 型 和 风 模 型 的 研 究 张 迪 ，王 维庆 ，王 海云 ，李 强 1 ． 新疆大学 电气工程 学院 ，新疆乌鲁木齐8 3 0 0 4 7 ； 2 ． 北京金风科创 风电设备有 限公 司，北京1 0 0 1 7 6 摘 要 建立 了一种新 的基于叶素理论的适用 于独立变桨方式 的风机气动模型 ，并与传统基于风能利用系数经验公 式的风机气动模型进行 了比较分析。同时基 于风剪塔影理论建立 了风模型 ，同新气动模型联 合仿 真 ，通过新气动模 型的输出仿真结果验证 了风剪塔影造成单只叶片的拍击力 1 P波动 ，3只叶片合成主轴气动转矩 的 3 P波动 。最后通 过独立设定桨距角 ，联合对风模 型与新气动模型仿真 ，验证 了新模 型对 于独立变桨方式 的适用性。 关 键 词 风 机 气 动模 型 ；风 剪 切 ；塔 影 效 应 ；独 立变 桨 S t u d y o n W i n d M o d e l a n d W i n d Tu r b i n e Ae r o d y n a m i c M o d e l S u i t a b l e f o r I n d i v i d u a l Pi t c h Co n t r o l Z h a n g D i ， Wa n g We i q i n g ， Wa n g H a i y u n ， L i Q i a n g 1 ． E l e c t ri c a l E n g i n e e ri n g C o l l e g e ， X i n j i a n g U n i v e r s i t y ， U r u mq i 8 3 0 0 4 7 ， X i n j i a n g ， C h i n a ； 2 ． G o l d w i n d S c i e n c e ＆ C r e a t i o n Wi n d p o w e r E q u i p me n t C o ． ， L t d ． ， B e ij i n g 1 0 0 1 7 6 ， C h i n a Ab s t r a c t A w i n d t u r b i n e a e r o d y n a mi c mo d e l 。 wh i c h b a s e d o n Bl a d e E l e me n t Mo me n t u m a n d s u i t a b l e f o r i n d i v i d u a l p i t c h c o n t r o 1 ． i S d e v e l o p e d a n d c o mp a r e d w i t h t h e t r a d i t i o n a 1 a e r o d y n a mi c mo d e l b a s e d o n r o t o r p o w e r c o e f f i c i e n t ． T h e wi n d mo d e l b a s e d o n t h e t h e o r y o f wi n d s h e a r a n d t o we r s h a d o w i s a l S O d e v e l o p e d ． By t h e c o s i mu l a t i n g ．山e c o n c l u s i o n s o f t h e wi n d s h e a r a n d t o we r s h a d o w l e a d i n g t o t h e 1 P fl u c t u a t i o n o f fl a p f o r c e i n o n e b l a d e a n d t h e 3 P f l u c t u a t i o n o f ma i n s h a f t t o r q u e b y 3 b l a d e s a r e v e r i fi e d t h r o u g h t h e o u t p u t s o f n e w a e r o d y n a mi c mo d e 1 ． T a k i n g i n d i v i d u a l p i t c h a n g l e a s a n i n p u t t o t h e mo d e l s ， t h e C O s i mu l a t i o n a l s o v e r i fie s t h e a p pl i ca bi l i t y o f ne w mo de l s i n i n d i v i du a l pi t c h c o n t r o 1 ． Ke v W 0 r d s w i n d t u r b i n e a e r o d y n a mi c mo d e l ； wi n d s h e a r ； t o we r e ff e c t ； i n d i v i d u a l p i t c h c o n t r o l 中 图分 类 号 T M3 1 5 文 献 标 识码 A 文章 编号 0 5 5 9 9 3 4 2 2 01 1 0 3 0 0 9 0 0 5 0 引言 风 机 变桨 控 制分 为 集体 变 桨 和独 立 变桨 ，随 着 风机 的容 量 。叶轮 和塔 架 的尺 寸 的增 加 ，如何 降低 风 机 载 荷 的 问 题 日益 凸 显 出 来 ， 独 立 变 桨 同 集 体 变 桨 相 比 的优势 在 于 降低 风机 载 荷 同时 保证 风 力发 电 机 的输 出 额 定 功 率 。 独 立 变 桨 是 变 桨 控 制 发 展 的 趋 势 。少 数 国外 先 进 的风 机 制造 商 已经 将独 立 变 桨应 用 于 兆 瓦 级 风 机 上 ， 国 内先 进 的 风 机 制 造 商 也 加 大 了对 独 立 变 桨 控 制 技 术 的研 究 。 随 着 风 机 变 桨 控 制 方 式 从 集 体 变 桨 过 渡 到 独 立 变 桨 ． 由 于 在 传 统 的 应 用 于 3叶 片 集 体 变 桨 控 制 策 略 的 风 机 气 动 模 型 中 ，3只 桨 叶 的 桨 距 角 对 风 机 气 动转矩 的影 响是 耦合 在一 起 的 ， 无法 利用 模 型来 分 析单 只 叶 片 的桨距 角 对 单 只 叶片 的拍 击 力 和升 力 的 影 响 ，所 以 需 要 回 到 叶 素 理 论 中来 研 究 适 用 于 独 立 变 桨 风 机 气 动 模 型 的 建 立 方 法 ， 因 此 需 要 一 种 新 的 建 模方 式 1 风机气动模型的建立 1 ． 1 应 用于集体 变 桨的风 机气 动模 型 参 考文献 【 1 】 中的风机 风 能利 用系 数 的经验 计 算 公 式 为 收 稿 日期 2 0 1 0 1 2 0 4 基金项 目国家 8 6 3计划支持项 目 2 0 0 9 A A 0 5 Z 4 4 5 作者简介 张迪 1 9 8 5 一 ，男 ，内蒙古包头人 ，硕士研究生 ，研 究方向为风机并网与控制技术． 第 3 7卷第 4期 张迪， 等 适用于独立变桨控制方式的风机气动模型和风模型的研究 G a ， ／ 3 C 一 C 一 。 鲁 c A 1 出增大桨距角对降低气动转矩的效果是非常明显的。 A 1 1 A l_一 2 1 A O ．08 一 式 中 ，A 为 叶 尖 速 比 ；A 为 中 间 变 量 ； 为 桨 距 角 ； CI o51 76， C2 I 1 6， C3 -- 0 ．4， C4 5， C 1 ， C6 0 ．0 06 8。 风 轮 吸 收 的 功 率 表 达 式 为 P P 竹尺 c A， V 3 式 中 ，P为 空 气 密 度 ；R 为 叶 片 半 径 ； V为 风 速 。 风 轮 产 生 的 转 矩 为 4 ∞ 式 中 ， 为 风 轮 转 速 。 叶 尖 速 比 A 的 计 算 公 式 A 5 V 将 式 4 和 5 代 人 式 3 中 得 出 风 机 气 动 转 矩 r厂的 表 达 式 为 盯 R G a， ／ A 6 基 于 式 1 、 2 、 6在 s i mu l i n k中 搭 建 模 型 结 构 ，仿 真 条件 为桨 距 角 取 0 。 ，空 气 密 度 为 1 ． 2 9 k g ／ m’ ， 叶 片 径 4 0 m 风 轮 转 速 从 0 r ／ mi n变 化 到 3 5 r ／ mi n，风 速 从 3 I n ／ S到 1 4 m／ s 。 仿 真 结 果 ① 风 机 的输 出转 矩 变化 情 况 如 图 】中 曲线 所 示 ；② 根 据文 献 [ 1 ] ， 最 大 值 对应 的 减 速 比 A 为 8 ． 1 ， 带 入 式 5 可 以 计 算 出 最 优 转 速 和 与 其 对 应 的 最 优 转 矩 ⋯用 光 滑 的 曲 线 连 接 形 成 图 1巾虚 线 。 虚 线 中 的 低 风 速 下 段 可 以 作 为 低 风 速 下 最 大 功 率 跟 踪 的 参 考 曲 线 。 O 5 l O 1 5 2O 2 5 30 35 转速／ r ． mi n 图 1 桨 距 角 为 0 。．不 同 风 速 下 叶 轮 转 矩 随 转 速 变 化 曲线 当 风 速 高 于 额 定 风 速 时 ．通 过 变 桨 电 机 驱 动 桨 叶 增 赶桨 距 角 ， 降 低 了 C ⋯ 从 而 降 低 风 机 气 动 转 矩 ，维 持 风 轮 转 速 额 定 。 罔 2绘 制 了 风 速 为 1 2 m／ s 、 桨 距 角 为 0 。 ，风 速 为 1 3 r r d s 、 桨 距 角 为 l 。 ，风 速 1 4 m／ s 、桨 距 角 5 。 转 矩 同 转 速 变 化 的 对 比 曲 线 ，呵 以 看 磐4 -K 5 霎4 2 1 图 2变 桨 对 气 动 转 矩 的影 响 这 种 建 模 的 方 式 可 以用 于 对 集 体 变 桨 策 略 的 仿 真 验 证 。但 是 对 于 独 立 变 桨 ， 由 于 模 型 没 有 独 立 的 桨 距 角 输 入 ， 而 只 能 输 入 3个 桨 叶 集 体 的 桨 距 角 ， 所 以 无 法 应 用 于 独 立 变 桨 策 略 的 验 证 。 1 _ 2应用 于独 立 变桨 的风 机气 动模 型 根 据 文 献 『 2 ] 的 叶 素 理 论 如 图 3所 示 。 若 以 叶 素 为 参 考 系 ， 风 将 会 叠 加 一 个 同 叶 素 转 动 方 向 相 反 ， 大小 相 等 的线 速 度 一 “，风 速 和 一 u矢 量 合 为 叶素 的 相 对 风 速 。 为 叶 素 真 正 感 受 到 的 风 速 。 由 于 相 对 风 速 的 作 用 ，在 风 机 叶 素 产 生 了 与 W 同 一 方 向 的 阻 力 ‰与 垂 直 于 阻 力 的 升 力 。 扫风平面 、 ／ d 、 浆 距 ii g ，＼ 角 相对转婆5 u 一 一 一 i n ＼＼相 风速 V d l ,o 风 速 一 W 图 3叶 素 受 力 分 析 1 d p C d w a L d r 7 1 d 去 p C L w 2 L d r 8 二 ， 9 式 中 ， 为 阻 力 系 数 ， 为 升 力 系 数 ， 同 一 样 没 有 量 纲 ， 二 者 的 取 值 变 化 决 定 于 攻 角 i ， 影 响 、 的 因 素 有 弯 曲 度 ， 叶 片 厚 度 ，表 面 粗 糙 度 和 雷 诺 数 ， 、 数 据 参 考 翼 型 为 NACA0 01 5；L为 弦 长 ；d r可 以理 解 为 一 段 很 短 的 叶 展 。攻 角 变 化 范 同 为 0 2 叮 r 。 为 丁 计 算 出 作 用 在 风 机 主 轴 上 的 气 动 转 矩 和 作 用 在 叶 片 上 的 拍 击 力 首 先 需 要 将 叶 素 上 的 阻 力 和 推 力 经 行 轴 向 和 切 向 分 解 与 合 成 如 图 4 。 每 一 小 段 叶 素 上 产 生 的 升 力 d 和 拍 击 力 d 5 O 5 O I ． 。 【 ／ 簿辞 _ lu n* s i n ／ d ， c o s ／“ 图 4 和 F 的分解与合成 由公 式 1 0 和 公 式 1 1 计 算 ⋯ 来 dFaaF t c ‘ sldF *s i n| d Fud F s i n／一 d h s I 式 中 ，，为人 流 角 。 1 0 1 1 为 _ _r实 现 仿 真 ，将 微 分 放 大 为 变 化 。 将 每 只 桨 叶 被 分 为 均 匀 的 1 0段 ，每 段 的 弦 长 为 此 段 端 弦 长 的平 均 值 。 分 别 建 立 每 段 叶 片 的 彳 『 l 力 AF a l不 1 f 升 力 △几 的 模 型 。 △ h p c ， L △ r 1 2 二 △ p C w 2 L A r 1 3 AF a AFIi ft C O S I 十 △F s i ni 0 、 4 AFu A s i n， △ O O S ， 1 5 将 一 只 叶 片 』 二 每 段 的 手 f l 击 力 △ 对 轮 毅 处 的 转 矩 △ t求 和 ， 最 终 汁算 f f } 每 只 川‘ 片 』 二 的 拍 山 力 与 3只 叶 片 总 体 轴 I J 的气 动 } 0 7 I j F f h i t △F i 1 6 1 3 1 0 ∑ ∑ △ 1 7 ，J J 式 巾 ，i 为 每 段 叶 片 的 标 ； 第 i 段 l ‘ 片 【 l I 心 距离轮 缎 的距 离 ； 为叶片数 、 基 于 式 1 2 一 1 7 s i mu l i n k r I I 建 ● 一 仿 真 模 ． 仿 真 条 件 为 空 气 密 度 1 ． 2 5 k g ／ m ， 【 】 t ‘ ” K 4 0 13 1 ， 风 速 从 3 m／ s至 1 3 m／ s ， 风 轮 转 速 从 0 钏 3 5 r ／ rai n． 桨 距 角 保 持 一 1 。 仿 真 结 果 如 I 5， t 轴 动 转 随 菥 风 速 ，风 轮 转 速 ， 变化 曲 线 。 『 f 1 于 桨 距 们 为一 l 。 ， 所 以 黎 体 的 主 轴 气 动 转 矩 大 于 l中 的 主 轴 气 动 转 矩 ， 罔 5巾横 坐标 为主轴 转速 ，纵 l 标 为气动 转 矩 ， 曲 线 下 方 的 面 积 为 叶 轮 吸 收 的 风 能 ，根 掂 最 人 面 积 的 方 法 同样 可 以 得 到 最 佳 转 速 值 以 及 对 应 的 最 佳 转 矩 值 作 为 低 风 速 下 的 最 大 功 牢 跟 踪 山 参 曲 线 1 ． 3两种建模 方 式的 比较 1 第 一 种 风 机 气 动 模 型 足 否 可 以 反 映 风 机 的 气 动 转 矩 关 键 在 于 对 ， A， l i t I 线 的 拟 合 准 确 度 ， 对 于 不 同 类 型 的 叶 片 需 要 大 馈 的 数 据 进 行 拟 合 ， l 第 二 种 气 动 模 型 的 建 立 只 需 要 知 道 片 的 升 力 系 数 团W a t e r P o w e r V o 1 ． 3 7 N o ． 4 图 5新 模型 不 同风 速 下 叶 轮 转 矩 随 转 速变 化 曲线 C ， 和 阻 力 系 数 c 。 ，这 两个 参 数 可 以 在 对 应 的 叶 片 参 数 叶 1 找 到 ， 能 够 直 接 放 入 模 型 中 ，而 且 由 于 单 只 叶 片 模 型 南一 段 段 叶 片 拼 接 而 成 ． 南 此 就 能 实 现 对 单 独 一 只 叶 片 的 桨 距 角 和 风 速 进 行 独 立 修 改 ， 可 应 用 于 独 立 变 桨 中 、 2 第 种 风 机 气 动 模 型 无 法 计 算 单 只 叶 片 的 拍 击 力 而 第 二 种 模 型 町 以计 算 。 ． 3第 二 利1 模 型 的 分 段 数 可 以 增 加 ，将 会 更 加 准 确 模 拟 叶 片 ，具 有 发 腰 潜 力 2 风剪塔影对风机气动模型输出的影响 2 ． 1 风 模 型 同风 机气 动模 型 的影响 与配 合 风 塔 影 最 终 刈 于 风 机 气 动 模 型 模 型 的 影 响 包 括 ① 单 叶片 L会产 生 1 P的拍 击 力波动 ，波 动 频 率 为 叶轮 旋 转周 期 的倒 数 。② 在 风 机 主轴 转 矩 中 产 生 3 P 的 波 动 ， 波 动 频 率 为 1 P 波 动 的 频 率 的 倍 。 塔 影 引 起 的 波 动 要 比风 剪 更 大 。 j 风 机 的 每 只 叶 片 被 分 为 l 0段 ，风 速 模 型 输 ⋯ 的 风 速 也 I』 这 1 0段 堪 一 段 巾 心 的 位 置 对 应 。 即 风 速 最 终 都 腹 化 为 叶 片 位 置 角 度 0 j 叶 片 位 置 点 到 轮 毅 I t 心 距 离 的 厕 数 从 而 风 速 模 型 输 f f I 3 0个 风 速 变 I |j 3只 叶 片 r f 1 的 3 0个 叶 片 段 位 对 应 。 2 ． 2风 剪切 以及 与模型 配合 风 是 指 当 随 着 高 度 的 不 同 ． 风 速 随 之 变 化 从 f 造 成 转 矩 和 拍 出 力 的 周 期 性 波 动 ， 参 考 文 献 [ 3 】 t的 风 切 公 式 L V h ； “ 1 8 凡 【 】 式 巾 ，v h 为 某 一 高 度 处 J x L 速 ； V h 。 为 轮 毅 高 度 处 风 速 ； 为 风 剪 切 经 验 指 数 ，取 1 ／ 7 。 将 式 1 8 转 化 为每 段 叶 片 位 置 角 度 的 数 如 式 1 9 V ， V 0 1 “ 1 9 凡 f 1 0t o t 2 0 式 叶l ，r 为 第 i 段 叶 片 距 离 轮 毅 的 径 长 J 蔓； 为 风 第 3 7卷第 4期 张迪 ， 等 适用于独立 变桨控制方式 的风机气 动模 型和风模 型的研究 轮 转 速 ；根 据 式 2 0 和 式 2 1 ， 作 用 存 单 只 叶 片 某 段 上 的 风 速 ， 会 叠 加 一 个 频 率 为 t o ／ 6 0 Hz的 分 量 。 仿 真 设 定 条 件 为 叶 片 半 径 4 0 m， 轮 毂 高 度 7 0 I n。转 速 为 1 9 r ／ mi n，桨 距 角 为 4 。 ，轮 毂 处 的 风 速 为 l 3 n t ／ s ． 观 察 与 一 段 叶 片 对 应 的 风 剪 模 型 的 输 出 风 速 如 降i 6。 图 6风 剪 切 效 应 下 的风 速 波 动 对 风 速 模 型 输 出 的 数 据 经 行 快 速 傅 里 叶 变 化 得 到 频 谱 图 如 图 7 。 图 7风 速 频 谱 图 m 7可 以 看 出 ． 风 速 中 除 了 直 流 分 量 外 ， 还 有 明 显 叠 加 了 0 ． 2 Hz 谐 波 含 量 。 风 剪 切 造 成 叶 片 拍 击 力 的 l P 波 动 如 图 8， 频 率 为 0 ． 2 Hz ， 即 等 于 叶 轮 旋 转 周 期 5 S的 倒 数 ， 而 仿 真 设 定 转 速 为 l 2 r ／ rai n，也 就 是 0 ． 2 r ／ s ， 即 5 s ／ r ， 即 叶 轮 旋 转 周 期 周 期 为 5 S仿 真 结 果 验 证 了 风 剪 对 单 叶 片 扣 击 力 和 气 动 转 矩 影 响 的 结 论 。 、- 一 ’＼一一 ， 。＼ ／ ， 、 ／ 、／＼ ／＼ 4 } 姆 3 ． 一 ⋯一一一⋯一 0 5 1 0 l 5 2 0 25 3 0 3 5 叫 问 ／ s 图 8风 剪切 造 成 叶 片 拍 击 力 的 1 P 波 动 风 切 造 成 主 轴 转 矩 3 P波 动 ． 频 率 为 0 ． 6 Hz如 罔 9 2 ． 3 塔 影效 应及 与模 型 配合 塔 影 效 应 指 塔 架 正 面 的 风 速 由 丁塔 架 的 作 川 而 g 1 1 55 主 1 1 5 0 釜 1 1 4 5 ＼1 1 4O 1 1 3 5 时 ／ s 图 9风剪切造成主轴转矩 3 P波动 减 小 ．所 以 当 叶 片 经 过 塔 架 的 竖 直 位 置 时 ． 叶 片 上 的 拍 击 力 会 减 小 ，单 只 叶 片 旋 转 一 周 后 由 风 速 引 起 的 拍 击 力 会 产 生 l P波 动 。 而 对 于 3Ⅱ 十 ‘ 片 风 机 ， 叶 轮 旋 转 一 周 后 ， 每 只 叶 片 郁 经 过 一 次 塔 架 位 置 ，所 以 在 主 轴 转 矩 中 会 产 生 3 P波 动 。 结 合 文 献 [ 3 】 的 塔 影 效 应 计 算 公 式 得 出 每 段 叶 片 对 应 位 置 的 风 速 表 达 式 V h o ， fiwt t s e 2 1 t ts e m “ 2 赫 2 2 ml 2 3 式 中 ，a为 塔 架 的 半 径；R 为 叶 片 半 径 ；H 为 塔 架 高 度 为 叶 片 与 塔 架 的 距 离 ； 塔 影 效 应 只 在 位 置 角 度 0为 叮 T到 耵 内 起 作 用 。 仿 真 的 其 他 条 件 不 变 考 虑 塔 影 效 应 的 情 况 下 观 察 到 单 只 叶 片 的 拍 击 中 产 生 了 1 P 波 动 ． 主 轴 转 矩 中 产 生 3 P 波 动 ， 塔 影 效 应 造 成 叶 片 拍 击 力 的 1 P波 动 ， 频 率 为 0 ． 2 Hz ， 如 冈 1 0 5． 0 g 4 5 主4 0 鱼 3 ． 5 童 3 ． 0 磐2 ． 5 2 0 1 5 图 1 0塔影效应造成叶片拍击力的 1 P波动 塔 影 效 应 造 成 主 轴 转 矩 3 I “波 动 ， 频 率 为 0 ． 6 Hz 如 图 l 1 一 1 2 l 1 孑 l o 宝9 耋 8 7 0 5 1 0 j 5 20 州 ／ s 2 5 3 0 3 5 图 11 塔 影 效 应 造 成 主 轴 转 矩 3 P 波 动 2． 4风剪 塔 影对风 机 的影 响 风 塔 影 都 会 造 成 只 n f 片 的 拍 击 力 中 的 l P波 动 ． 主 轴 转 矩 【 f 1 的 3 P波 动 ，但 是 ⋯ 丁塔 影 效 心 对 风 速 的 影 响 比 风 剪 切 更 加 严 乖 。 所 以 塔 影 效 应 产 生 的 波 动 更 加 严 重 。 3 独立变桨对叶片拍击力的影响 为 了减 小 风对 叶 片拍 击 力 的影 响 ，需 要对 每 只 叶 片 的 桨 距 角 进 行 独 立 的 控 制 ， 风 速 为 1 3 m ／ s ，转 速 为 1 2 r ／ mi n的 条 件 下 ，3只 桨 叶 的 桨 距 角 分 别 为 2 。 ，4 。 ，6 。 ，从 仿 真 结 果 可 以 看 出 桨 距 角 的 增 大 对 叶 片 的 拍 击 力 的 影 响 如 冈 1 2 Z ＼ 伯 图 1 2独立 变 桨 对 叶 片 拍 击 力 的影 响 4 结论 1 综合 比较 了基 于 C 经 验公 式与 基于 叶 素理 论 的 风 机 气 动 模 型 ， 通 过 分 析 两 种 模 型 构 架 思 想 和 仿 真结果 ．说 明 了基 于 叶素理 论 的风机 气 动模 型 应 用 于独 立 变 桨 的 正 确 性 。 2 在 基 于 叶 素 理 论 的 风 机 气 动 模 型 与 基 于 风 剪 切 塔 影 效 应 理 论 的 风 模 型 基 础 上 ．通 过 分 析 二 者 联 合 仿 真 的 结 果 ， 验 证 了 风 剪 塔 影 导 致 单 只 叶 片 上 接 第 4 5页 拍 击 力 的 1 P波 动 和 3只 叶 片 产 生 主 轴 转 矩 中 的 3 P 波 动 。 参考文献 [ 1 ] He i e r S ． G r i d I n t e g r a t i o n o f wi n d E n e r g y C o n v e r s i o n S y s t e m s [ M] ． 1 9 9 8 ． [ 2 ] 叶杭冶． 风 力发 电机组 的控 制技术 [ M] ． 机 械工业 出版社 2 0 0 2 5 3 3 3 ． [ 3 ] 李少林 ，张兴 ，谢震 ，等． 双馈风力发 电系统 3次功率脉动 的 研究 [ J ] ． 电网技术 2 0 1 0 ，3 4 4 3 8 3 9 ． [ 4 ] S o n g S H ， J e o n g I ； C ， O h J H， e t a 1 ． S i mu l a t i o n o f o u t p u t c h a r a c t e r - i s t i e s o f r o t o r h l a d e s u s i n g a h a r d wa r e i n- l o o p wind t u r b i n e s i mu l a -- t o r [ J ] ． I E E E T r a n s 0 11 E n e r g y C o n v e r s i o n ， 2 0 0 5 ， 3 6 1 7 9 1 - 1 7 9 6 ． [ 5 ] A n d e r s o n C R i c h o n J ， C a m p b e l l T J ． A n A e r o d y n a m i c Mo m e n t Co n t r o l l e d S u rfa c e f o r Gu s t L o a d Al l e v i a t i o n o n W i n d Tu r b i n e Ro t o r s [ J ] ． 1 E E E T R AN S A C q 、I O N S O N C O N T R O L S Y S T E MS T E C HN O L OGY 1 9 9 8 6 5 5 7 8 5 8 0 ． [ 6 ] D o l a n D S L ， L e h n P W． S i mu l a t i o n m o d e l o f w i n d t u r b i n e 3 p t o r q u e o s c i l l a t i o n s d u e t o w i n d s h e a r a n d t o w e r s h a d o w[ J ] ． I E E E T r a n s o n En e r g y C o n v e r s i o n，2 0 0 6，2 1 3 71 7 - 7 2 4 ． [ 7 ] G e y l e r M， C a s e l i t z P I n d i v id u a l B l a d e P i t c h C o n t r o l D e s i g n for L o a d R e d u c t i o n o n L a r g e Wi n d T u r b i n e s [ c y ／ P r o c e e d i n g s o f E u r o p e a n W i n d E n e r gy Co n f e r e n c e ． 2 0 0 7 1 4 ． [ 8 ] 董萍 ，吴捷 ，杨金明 ，等 ． 风力发 电机组 建模研究 现状[ J ] ． 太 阳 能 学 报 2 0 0 4 ，2 5 5 6 1 3 6 1 5 ． [ 9 』 V a n B ld a r s G E ， B o n g e r s P M M． F l e x i b l e w i n d t u r b i n e mo d e l v a l i d a t i o n [ J ] ． Wi n d E n g i n e e r i n g ，1 9 9 2 ， l 6 4 2 4 7 2 5 6 ． [ 1 0 ] Mu lj a d i E ， B u t t e r f i e l d C P ． P i t c h C o n t r o l l e d V a r i a b l e - S p e e d Wi n d T u r b i n e G e n e r a t i o n [ J ] ． I E E E T r a n s o n I n d u s t r y A p p l i c a t i o n s ， 2 0 0 1 ， 3 7 1 2 4 0 2 4 6 ． 责任 编辑高 瑜 表 5库水 位 升 降 循 环 下 的坝 体 应 力 变 形 特 征 值 于 距 坝 顶 约 2 ／ 3坝 高处 ； 同 时 ，坝 基 深 厚 覆 盖 层 受 坝 体 的 影 响 ，也 会 向 上 、下 游 发 生 水 平 位 移 。 与 粘 土 心 墙 土 石 坝 一 样 。 沥 青 混 凝 土 心 墙 坝 也 存 在 明 湿 的 应 力 拱 现 象 ， 心 墙 的应 力 比 坝 壳 堆 石 的 席 力 低 ， 不 过 蓄 水 后 这 种 拱 效 应 减 弱 ； 蓄 水 后 ， 由 于 库 水 压 力 作 用 在 心 墙 上 ，从 而 导 致 上 游 坝 壳 及 上 游 围 堰 的 应 力 水 平 增 大 ，沥 青 混 凝 土 心 墙 及 下 游 堆 石 的 应 力 水 平 明 降 低 。 另 外 ，坝 体 的 应 力 变 形 分 布 特 征 受 水 田W a t e r P o w e r V o 1 ． 3 7 N o ． 4 位 升 降 的影 响 不 大 。 参考文献 [ 1 ] 徐 晗 ，汀 明元 ，等． 深厚覆盖 层 3 0 0 i n级超 高土质心墙 坝应力 变形 特征 [ J ] ． 岩土力学，2 0 0 8 ，2 9 t 6 4 6 8 ． [ 2 ] 温 州，邵磊． 深厚覆盖层 3 5 0 m 级心墙堆石 坝动强度计算 分析 [ J ] ．i 峡大学学报 自然科学版 ，2 0 0 9 ，3 1 2 3 6 3 9 ． 责任 编辑杨健