(1.云南电网公司,云南 昆明 650011; 2.云南省电力设计院,云南 昆明 650051; 3.国网智能电网研究院,北京 100124; 4.北京市劳动保护科学研究所 国家环境保护城市噪声与振动控制工程技术中心,北京 100054)
(1. Yunnan Power Grid Corporation,Kunming 650011,Yunnan,China)(2. Yunnan Electric Power Design Institute,Kunming 650051,Yunnan,China)(3. Smart Grid Research Institute,Beijing 100192,China)(4. State Environmental Protection Engineering Center for City Nois
electrical equipment; disconnector; modal analysis; under the earthquake; dynamic time history analysis
备注
利用试验和数值仿真的方法,对高压隔离开关的模态动力学特征及其在水平与竖向地震波作用下的动力时程响应特征进行了研究。结果 表明,隔离开关的低阶共振频率与地震波较为接近,容易受到地震激励力的影响。在Ⅸ度罕遇地震烈度条件下,其瓷柱根部区域可能会因为最大弯矩超出设计要求而发生断裂失效,从而导致隔离开关设备的整体失效。经计算分析,采用安装隔振器的方法,可以降低瓷柱根部区域的最大弯矩值,从而保证结构的抗震需求。
The modal dynamic characteristic and dynamical time-history response of high-voltage disconnector under horizontal and vertical seismic waves are studied by experiment and numerical simulation method. The result shows that the low-order resonance frequency of the disconnector is similar with that of the seismic wave,so it is easily to be effect by the seismic excitation force. Under the strong earthquake in Ⅸ degree,the porcelain column insulator of the high-voltage disconnector breaks because of its maximum bending moment exceeds the design requirements. Through calculation and analysis,it could reduce the maximum bending moment of porcelain column insulator to ensure the seismic demand of the structure by installing the isolator.
引言
变电站高压设备大多具有较低的固有频率,部分低阶固有频率低于10 Hz,与地震波频率接近。因此,在地震发生时,高压设备通常会在地震波作用下发生共振,从而引发设备的破坏和失效(Jean-Bernard,2007; 韩军科等,2008; 卢智成等,2009)。
近年来,地震灾害已在国内外许多地区造成电力设备损坏事故,并引发了大面积停电。1994年美国Northridge地震中,230 kV与550 kV变电站的电气设备严重破坏(ASCE-TCLEE,1997)。1995年日本阪神地震造成48座变电站设备损毁,100万户停电(日本建筑学会等,1998)。1999年土耳其Kocaeli地震中,多个变电站变压器、开关设备受到不同程度破坏(AIJ et al.,2001)。2000年日本新泻西部地震造成500 kV日野变电站瓷柱断裂和11座变电站受损①。2008年汶川地震造成110 kV及以上变电站停运90座,线路停运181条,电力损失负荷685万 kW,高压设备损毁严重(谢强,2008)。
为了有效避免地震灾害对高压设备的不良影响,需要对地震波作用下设备的动力特性进行深入研究。然而,当前变电站高压设备模态与动力时程特性的研究工作仍然不够充分,这使得高压设备的抗震改进工作存在一定困难。鉴于此,本文主要针对地震波作用下隔离开关的动力特性进行研究,通过模态试验和数值模拟分析,掌握此类设备的模态动力特性和地震激励下的动力时程响应特点,从而为改进此类设备抗震性能和抗震措施提供有效依据和指导。
1 模态动力学特性分析
1.1 基于特征系统理论的模态分析方法特征系统理论是一种基于多输入多输出的时域模态识别理论方法,最早由美国航天局Langley研究中心提出(张枢文等,2007; Juang,Pappa,1985)。根据这一理论和实际分析要求,本次模态计算分析模型如下:
隔离开关结构的系统振动方程式:
M(¨overq)+C(·overq)+Kq=f(q,t)=Bfu(t).(1)
其中,M、C、K分别为系统的质量阵,阻尼阵和刚度阵,Bf为输入分配矩阵,将式(1)写成状态方程的形式:
{(·overx)=Ax+Bu,
y=Cx.(2)
其中,A、B、C分别为系统矩阵、控制矩阵和测量矩阵,其表达式如下:
x=[q
(·overq)], A=[0 I
-M-1K -M-1C], B=[ 0
M-1Bf].(3)
若结构没有发生变化,则系统矩阵和控制矩阵为常数矩阵,由此可得一组一阶常微分方程组,根据初始条件积分可得:
x(t)=eA(t-t0)+∫t∫t0eA(t-τ)Bu(τ)dτ, t≥t0.(4)
如果使用等间隔采样,再加之零阶保持在一个采样间隔内采样值不变的特性,可得:
x[(k+1)T]=eATx(kT)+∫t∫0eAsdsBu(kT).(5)
这样就可以得到系统离散时刻的脉冲响应矩阵式:
h(0)=CB,h(T)=CAB,…,h(k)=CAkB.(6)
系统的Hankel矩阵为
H(k)=[h(k)h(k+1)… h(k+s)
h(k+1)h(k+2)… h(k+s+1)
h(k+r)h(k+r+1)… h(k+r+s)].(7)
将H(0)做奇异值分解,选取前r个最大奇异值对系统进行缩减,最终可得系统的脉冲响应函数为
h(k)=ITpPrD-1/2r[D-1/2rPTrH(1)QrD-1/2r]kD-1/2rQTrIm.(8)
其中,Pr、Qr为正交矩阵中对应前r个最大奇异值的SVD分解的列,Dr为r个最大奇异值组成的对角矩阵,Ip为p阶单位矩阵,p为输出向量的维数。
由此可得,系统的最小实现可采用:
{A=D-1/2rPTrH(1)QrD-1/2r,
B=D-1/2rQTrIm,
C=ITpPrD-1/2r.(9)
通过上述模态分析方法进行系统实现,就可以求解系统的模态参数,从而直接对系统矩阵A进行特征分解,得到动力学系统的固有频率和阻尼比。
1.2 模态试验分析为了准确掌握隔离开关结构的模态动力特性,本文采用力锤激励、多点输入输出的方法进行隔离开关结构模态试验。某隔离开关模态试验在实际设备现场进行,其结构图纸和分析模型如图1所示。
对隔离开关结构进行模态试验进行数据分析,得到前6阶模态的主要频率特征,如表1所示,其结构振型如图2所示。
six order experimental modal of the disconnector1.3 结构模态计算与验证利用有限元模型计算方法对隔离开关设备的模态特征进行分析,可以得到前6阶固有频率的计算结果,如表2和图3所示。
experimental modal of the disconnector计算结果显示,模态试验与仿真结果具有较好的一致性,试验达到了预期效果。模态分析也揭示了隔离开关易受地震影响的重要原因:一方面,隔离开关的前6阶模态处于10 Hz以下的低频区域,当该设备处于地震II类区域或相近类别的场地时,
会因为场地传递的地震波频率与设备频率相近(Jean-Bernard,2007; 建筑抗震设计规范,2010),而在地震时出现设备结构共振并引发失效。另一方面,隔离开关的低阶模态以弯扭变形为主,在地震时瓷柱根部区域和连接位置容易出现应力集中,当地震烈度较高时,这些位置会逐渐发生塑性或脆性变形,从而引发失效。
2 地震波作用下隔离开关的时程动力响应分析
从近年来地震现场户外高压设备的破坏情况来看,高压隔离开关的主要破坏区域集中在瓷柱根部区域(谢强,2008; AIJ et al.,2001; 日本建筑学会等,1998; ASCE-TCLEE,1997)佐藤诰章,佐藤清隆,当麻纯一.2001.2000年鸟取西部地震变压器被害要因的究明.日本电力中央研究所报告NO.U01023.。因此,
图3 隔离开关的前6阶模态振型计算结果
Fig.3 Calculation results of the first six order numerical modal shapes of the disconnector因此,在研究地震波作用对设备的影响时,需要重点关注瓷柱根部最大应力,即地震作用下瓷柱根部最大弯矩。
在地震作用下,高压隔离开关的瓷柱顶部会出现较大的位移,这可能引发瓷柱柱状结构的碰撞和断裂破坏,因而也需关注高压隔离开关结构顶部的最大位移。
为了有效计算隔离开关在地震波作用下的时程动力响应,本文在上述模态试验与模拟验证基础上,建立了地震时程计算的有限元模型(图4)。模型的本体、支座和支架利用实体单元构建,主刀和地刀利用质量单元构建,隔离开关之间的连接法兰根据《电力设施抗震设计规范》(2013)构建,其结构等效为弯曲刚度与胶装部位法兰弯曲刚度相同的圆形截面进行模拟,计算公式如下:
Kc=(6.54dch2c)/(te)×107.(10)
其中,dc为瓷柱胶装部位外径,hc为瓷柱与法兰胶装高度,te为法兰与瓷柱之间的胶装厚度。在隔离开关有限元模型中,混凝土立柱底部采用全约束固定边界条件,即三方向位移(UX,UY,UZ)和转角(ROTX,ROTY,ROTZ)为0。
为了准确掌握隔离开关在强地震中随时间变化的动力响应情况,笔者采用时程分析法进行地震响应分析,分别以阪神地震波和El Centro地震波作为数值模拟计算的输入载荷,如图5所示。隔离开关按抗震设防烈度Ⅸ度进行时程分析,地震波的计算输入参数如表3所示。
图4 隔离开关的动力时程计算有限元模型
Fig.4 Finite element model of dynamic time history calculation of the disconnector图5 输入地震波的水平向和竖向时程曲线
(a)阪神地震波(水平向);(b)El Centro地震波(水平向);(c)阪神地震波(竖向);(d)El Centro地震波(竖向)
Fig.5 Horizontal and vertical time history curves of the imported earthquake wave (a)Kobe earthquake wave(horizontal);(b)El Centro earthquake wave(horizontal);(c)Kobe earthquake根据上述隔离开关模态分析可知,隔离开关低阶模态的振型弯曲特征明显,瓷柱根部区域应力集中较其他区域更为严重,因而本文主要关注隔离开关瓷柱根部区域的地震响应情况。对结构输入水平与竖向地震波并计算,可以得到隔离开关在地震波作用下的时程响应特征。
图6 不同强度水平和竖向地震波作用下的瓷柱根部弯矩时程响应曲线
(a-1)阪神地震波(140 gal)水平作用;(b-1)El Centro地震波(140 gal)水平作用;(a-2)阪神
地震波(140 gal)竖向作用;(b-2)El Centro地震波(140 gal)竖向作用;(c-1)阪神地震波
(620 gal)水平作用;(d-1)El Centro地震波(620 gal)水平作用;(c-2)阪神地震波
(620 gal)竖向作用;(d-2)El Centro地震波(620 gal)竖向作用
Fig.6 Time history curve of bending moment of porcelain column insulator under the different horizontal and vertical earthquake wave (a-1)under the horizontal Kobe earthquake(140 gal);(b-1)under the horizontal EI Centro earthquake(140 gal);(a-2) under the vertical Kobe earthquake(140 gal);(b-2)under the vertical EI Centro earthquake(140 gal);(c-1)under the horizontal Kobe earthquake(620 gal);(d-1)under the horizontal EI Centro earthquake(620 gal);(c-2)under the vertical Kobe earthquake(620 gal);(d-2)under the vertical EI Centro earthquake(620 gal)如图6所示,从地震响应的计算结果可以看出,隔离开关在阪神地震波和El Centro地震波作用下表现出相应的时程响应特征。在Ⅸ度多遇地震条件下,隔离开关处于弹性变形阶段,其结构动力响应曲线的幅值较小; 而在Ⅸ度罕遇地震作用下,其结构响应曲线的幅值明显增加。
表4显示了隔离开关的地震响应时程分析结果。根据本设备支柱绝缘子的设计要求,瓷柱的破坏弯矩为20 kN·m。
从时程分析结果来看,在Ⅸ度多遇地震(加速度峰值140 gal)条件下,隔离开关各部分仍具有足够的刚度,可以保证设备的使用要求; 但在Ⅸ度罕遇地震(加速度峰值620 gal)条件下,根
据《电力设施抗震设计规范》(2013),瓷柱承受的最大弯矩不应大于20 kN·m /1.67,而隔离开关瓷柱的根部最大弯矩则已接近或超出相关设计要求,进入塑性或脆性变形的危险状态,随时可能发生破坏和失效。因此,在此种情况下应采取必要措施,阻止地震能量向瓷柱区域传递,保证设备的安全平稳运行。
为此,考虑在隔离开关横梁位置安装金属橡胶隔振器。如图7a、b所示,在隔离开关4个横梁位置,每个位置分别安装2个水平向隔振器和2个竖向隔振器,并采用combin14单元模拟阻尼器减震效果。隔离开关的混凝土立柱底部仍采用全约束固定边界条件,即三方向位移(UX,UY,UZ)和转角(ROTX,ROTY,ROTZ)为0。隔振器的滞回曲线如图8所示。
完成金属橡胶隔振器的建模安装后,在Ⅸ度罕遇地震条件下进行隔离开关的地震响应时程计算,得到隔离开关瓷柱根部的弯矩时程响应曲线,如图9所示。
表5显示了安装隔振器后的隔离开关地震响应时程分析结果。从时程分析结果来看,在Ⅸ度罕遇地震(加速度峰值620 gal)条件下,隔离开关瓷柱最大弯矩比安装隔振器前明显下降,满足最
图8 隔振器滞回曲线
(a)静态加载;(b)1 Hz加载;(c)2 Hz加载;(d)3 Hz加载
Fig.8 Hysteresis curve of vibration isolator(a)static loading;(b)1 Hz loading;(c)2 Hz loading;(d)3 Hz loading图9 在Ⅸ度罕遇地震下安装隔振器的隔离开关瓷柱根部弯矩时程响应曲线
(a-1)阪神地震波水平向作用;(a-2)阪神地震波竖向作用;(b-1)El Centro地震波
水平向作用;(b-2)El Centro地震波竖向作用
Fig.9 Bending moment-time curve of Porcelain column insulator under the Ⅸ degree earthquake after installing the vibration isolators (a-1)under horizontal Kobe earthquake wave;(a-2)under vertical Kobe earthquake wave;(b-1)under horizontal El Centro earthquake wave; (b-2)under vertical El Centro earthquake wave大弯矩不大于20 kN·m /1.67的取值设计要求。另外,瓷柱顶部的最大位移虽然有所增加,但与瓷柱间的距离相比仍然处于可接受范围内。因而,可以认为隔离开关在此条件下能够安全稳定的运行。
3 结论
本文对高压隔离开关结构开展了模态试验分析、模态仿真分析以及Ⅸ度多遇/罕遇地震条件下的动力响应分析。结果表明,隔离开关的前六阶模态处于10 Hz以下的低频区域,当该设备在地震II类及附近类别的场地运行时,会因为其频率与地震波频率相近,而引发结构共振和失效。同时,隔离开关的低阶模态以弯扭变形为主,当地震烈度较高时,这些位置会逐渐发生塑性或脆性变形,引发失效。在水平和竖向地震波作用下进行结构响应分析表明,隔离开关在Ⅸ度多遇地震条件下,具有足够的刚度来保证设备安全运行,但在Ⅸ度罕遇地震条件下,其瓷柱区域可能会首先处于塑性或脆性变形状态,随时发生失效。为此,在Ⅸ度罕遇地震条件下,本文选用金属橡胶隔振器进行了隔离开关的隔震计算分析,取得了较好的隔震效果,将为进一步开展高压隔离开关的抗震研究提供重要参考依据。
本文在撰写过程中得到康钟绪博士的帮助,在此表示衷心感谢!
- 韩军科,朱全军,杨风利.2008.基底隔震高压隔离开关抗震性能分析[J].电力建设,(6):39-42.
- 卢智成,程永锋,代泽兵.2009.基底隔震电气隔离开关的动力反应分析[J].武汉大学学报(工学版),10(S1):262-266.
- 日本建筑学会,日本土木工程学会,日本地盘工学会.1998.阪神淡路大震灾调查报告共通编3:都市防灾报告[R].日本丸善株式会社.
- 谢强.2008.“5·12”汶川大地震电力系统震害调查初步报告[R].上海:同济大学.
- 张枢文,嵇春艳,吴立人.2007.ERA法识别大型结构损伤与定位[J].江苏科技大学学报(自然科学版),21(3):5-9.
- AIJ,JSCE,JGS.2001.Report on the Damage Investigation of the Kocaeli Earthquake in Turkey[R].Maruse Corp.,Tokyo,Japan.
- ASCE-TCLEE.1997.Northridge Earthquake:Lifeline Performance and Post-earthquake Response[J].American Society of Civil Engineering Technical Council on Lifeline Earthquake Engineering.
- Jean-Bernard Dastous.2007.Guidelines for Seismic Design of Flexible Buswork between Substation Equipment[J].Earthquake Eng. Struct.Dyn.,36:191-208.
- Juang J.,Pappa R.1985.An Eigensystem Realization Algorithm for Modal Parameter Identification and Model Reduction[J].Journal of Guidance,8(5):620-627.
- GB50011-2010,建筑抗震设计规范[S].
- GB50260-2013,电力设施抗震设计规范[S].