中国是地震多发国家之一，而且众多风电场分布在地震高烈度区。风电机组的叶轮叶片位于塔筒顶部，刚度较小，在地震力作用下易产生端部的“鞭梢效应”，不利于叶片安全运行，而塔架顶部支撑机舱和叶轮，结构上头重脚轻，也不利于抗震。因此在风电机组结构安全设计时，考虑地震作用的影响至关重要。

 

　　近年来，不少研究人员研究了风电机组塔筒的动力特性。棚邉隆进行了600kW大型风力发电设施的运转振动与地震作用的组合荷载的动力分析。祝磊用SAP2000计算了停机状态下7台风电机组的地震反应谱响应，提出了钢筒塔架的剪力和弯矩计算公式，并研究了叶片的弯矩分布规律。张斌用ANSYS建立了某1MW风电机组钢塔筒模型，研究了其抗风、抗震特性。余长海建立了某1.25MW风电机组的ANSYS有限元模型，研究了钢筒和钢筋混凝土筒的动力特性及风荷载下的静、动态反应。贺广零考虑钢筒和钢筋混凝土筒两种塔架形式，用ANSYS建立了某1.25MW风电机组的有限元模型，并分别计算了风、地震作用下的动力响应。毕继红设计了某2MW风电机组预应力钢筋混凝土塔筒，并用ANSYS分析其地震响应。

 

　　目前，大部分研究都集中在风电机组钢制塔筒，对混凝土塔筒的研究还很有限。本文建立风电机组混凝土塔筒有限元模型，进行模态分析，并按照中国抗震规范中的标准设计地震反应谱，选择某一特定的设计分组以及设防烈度，进行地震反应谱分析。

 

 

　　本文分别选取钢筋混凝土塔筒底部外径6m、7m、8m、9m、10m、11m，壁厚从200mm-600mm每100mm一个间隔的塔筒形式，总共分析了30种规格。钢筋混凝土塔高100m，外径从下到上按0.02的坡度逐渐减小，筒壁采用强度等级C30混凝土，弹性模量为3.0×104MPa，质量密度为2410kg/m3;钢筋采用HRB400级，弹性模量为200GPa。运用ANSYS有限元软件建立模型，参照某大型风电机组参数，把叶轮和机舱简化成一个550t的集中质量，采用mass21单元加在塔筒顶部;塔筒采用beam189单元模拟变截面，《高耸结构设计规范》6.2.2规定计算结构自振特性时，混凝土高耸结构的截面刚度取0.85EcI，因此ANSYS建模时混凝土的弹性模量取0.85Ec，材料弹性模量E=2.55×104MPa,泊松比为0.167，材料密度为2410kg/m3;塔筒与地面刚接。定义沿塔筒高度方向为Z轴正方向，塔筒横截面方向分别为X轴和Y轴。以塔筒外径7m，壁厚200mm为例，有限元模型如图1所示。

 

 

　　模态分析用来确定结构或构件的振动特性，即固有频率和振型。同时，模态分析也是其他动力学分析(如反应谱分析)前期必须完成的环节。在ANSYS中提取模态的方法有：BlockLanczos法、Subspace法、Reduced法、Unsymmetric法、Damped法等。本文采用BlockLanczos法进行模态分析，提取前10阶频率和对应的振型。部分频率及相应振型图如表1及图2所示。
 

 

 

　　由于把风电机组的机舱和叶轮看作集中质量加载在混凝土塔筒的顶端，因此本文所建立模型属对称型。模型在X、Y方向质量和刚度分布相同，分别导致第1阶频率和第2阶频率以及第3阶频率和第4阶频率一样。通过观察动态显示，其中第1阶振型特征和第3阶振型特征沿Y方向弯曲振动;第2阶振型特征和第4阶振型特征沿X方向振动;第7阶振型特征沿Z方向轴向振动，第10阶振型特征沿X、Y平面扭转振动。

 

　　选取塔顶外径7m、11m壁厚200mm的模型，对X、Y方向的1阶和2阶模态进行无量纲化处理，横坐标为振型位移与塔顶位移的比值，纵坐标为塔架Z坐标与塔架高度的比值，见图3。

 

 

　　因1阶频率和2阶频率相同，由图3可知，风电机组混凝土塔筒的1阶、2阶模态重合;且不同尺寸的混凝土塔筒1阶、2阶模态曲线走向相似。

 

　　经建模计算，本文所列其他尺寸混凝土塔筒振型特征与之相似，在此不一一列出。

 

 

　　依据中华人民共和国国家标准《建筑抗震设计规范》(GB50011-2010)所采用的相对于重力加速度的单质点绝对最大加速度，即Sa/g与体系自振周期T之间的关系作为设计反应谱。本文取抗震设防烈度8度，多遇地震，场地类别IV，设计地震分组三，特征周期Tg=0.90s。依照中华人民共和国国家标准《烟囱设计规范》(GB50051-2002)中混凝土烟囱规定的阻尼比，取0.05。经计算本文所对应的地震影响系数曲线图如图4所示。
 

 

　　模态分析后，继续用ANSYS有限元软件进行地震反应谱分析，地震加速度分别从X、Y两个方向分别输入。

 

　　地震加速度从X方向输入，进行有限元振型分解反应谱法分析得到位移场分布等值线图如图5所示。

 

　　图5显示，将风电机组的叶轮和机舱简化成集中质量加载在塔筒顶部，得出高100m，底部外径7m，壁厚200mm的混凝土塔筒顶部位移最大值为0.256m，为高度的0.256%。

 

　　将叶片和机舱简化成集中质量加载在塔筒顶部导致该模型属对称型，因此从X方向和Y方向输入地震加速度得出的结果是一样的，本文不再列出。

 

　　输入地震反应谱计算后，选取混凝土塔筒不同高度截面所对应的弯矩，并将其进行无量纲化处理，横坐标为塔筒各个截面弯矩与塔底弯矩的比值，纵坐标为塔架Z坐标与塔架高度的比值，得出混凝土塔筒各个截面弯矩沿塔筒高度的分布情况如图6所示。

 

　　图6显示，不同尺寸的风电机组钢筋混凝土塔筒在地震作用下，弯矩沿着塔筒高度的分布曲线大致相同且接近重合，呈近似线性关系。另外，最大弯矩发生在塔筒底部。

 

　　图7显示，不同尺寸的风电机组钢筋混凝土塔筒在地震作用下，剪力沿着塔筒高度的分布曲线走向趋势大致相同且接近重合。另外，最大剪力均发生在塔筒底部。
 

 

 

　　(1)通过将风电机组的叶片和机舱简化成集中质量加载在塔筒顶部，用有限元软件建立混凝土塔筒简化模型，模拟了混凝土塔筒的动力特性。

 

　　(2)进行模态分析计算，得出风电机组混凝土塔筒的各阶频率以及相应的振型图。由于简化后的混凝土塔筒属对称型，因此不同尺寸混凝土塔筒的1阶频率和2阶频率以及3阶频率和4阶频率相同，所对应的振型图也一致，只是方向不同。另外，不同尺寸的混凝土塔筒1阶、2阶模态曲线走向相似。

 

　　(3)通过输入地震加速度谱进行反应谱分析，得出多种规格混凝土塔筒的截面弯矩、剪力沿塔筒高度的分布情况，结果显示其走向大致相同且接近重合，弯矩、剪力沿塔筒分布呈近似直线关系。弯矩和剪力在风电机组钢筋混凝土塔筒底部最大。