CN113408072B - 风力机柔塔系统固有振动特性快速建模及仿真方法 - Google Patents
风力机柔塔系统固有振动特性快速建模及仿真方法 Download PDFInfo
- Publication number
- CN113408072B CN113408072B CN202110710893.9A CN202110710893A CN113408072B CN 113408072 B CN113408072 B CN 113408072B CN 202110710893 A CN202110710893 A CN 202110710893A CN 113408072 B CN113408072 B CN 113408072B
- Authority
- CN
- China
- Prior art keywords
- tower
- wind turbine
- equation
- flexible
- coordinate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000004088 simulation Methods 0.000 title claims abstract description 8
- 239000011159 matrix material Substances 0.000 claims abstract description 50
- 238000012546 transfer Methods 0.000 claims abstract description 50
- 230000005540 biological transmission Effects 0.000 claims abstract description 11
- 230000005484 gravity Effects 0.000 claims abstract description 7
- 230000009471 action Effects 0.000 claims abstract description 6
- 238000006073 displacement reaction Methods 0.000 claims description 27
- 230000033001 locomotion Effects 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 6
- 230000003993 interaction Effects 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 238000013016 damping Methods 0.000 claims description 4
- 238000005452 bending Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000010008 shearing Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 10
- 238000004458 analytical method Methods 0.000 abstract description 7
- 230000008569 process Effects 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- RZVAJINKPMORJF-UHFFFAOYSA-N Acetaminophen Chemical compound CC(=O)NC1=CC=C(O)C=C1 RZVAJINKPMORJF-UHFFFAOYSA-N 0.000 description 1
- 235000017166 Bambusa arundinacea Nutrition 0.000 description 1
- 235000017491 Bambusa tulda Nutrition 0.000 description 1
- 241001330002 Bambuseae Species 0.000 description 1
- 235000015334 Phyllostachys viridis Nutrition 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 239000011425 bamboo Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/06—Wind turbines or wind farms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了一种风力机柔塔系统固有振动特性快速建模及仿真方法,包括以下过程:对风力机柔塔系统进行合理的简化,将系统离散成各多体动力学元件;基于多体系统传递矩阵法推导考虑重力作用的塔筒梁传递矩阵;确定其它多体动力学元件的传递矩阵;确定边界条件,建立柔塔系统总传递方程;求解总传递方程,得到柔塔系统的振动特性,本发明充分考虑风力机柔塔系统的结构细节,计算效率高并且仿真精度高,能实现对风力机柔塔系统的快速建模与固有振动特性分析,同时,为工程上类似于柔塔结构的细长柔性体系统动力学建模与振动特性分析提供参考。
Description
技术领域
本发明涉及多体系统动力学技术领域,特别涉及一种风力机柔塔系统。
背景技术
随着风能产业的快速发展,为了更有效地捕获风能,大型化的兆瓦级风力机正逐步成为主流。近年来,兆瓦级风力机主要应用高柔性塔筒支撑。柔性塔可支撑风力机轮毂到110~150m的高度,成本较低。塔架一般为薄壁钢空心圆筒,具有柔性大、阻尼小、重量轻的特点。但是风力机柔塔引起的流体-结构相互作用(FSI)振动更加突出。风力机的生命周期主要受其动力学响应影响,为了避免风力机在流速范围内产生共振,必须知道其固有频率。而风力机的低阶固有频率主要受塔筒影响,因此预测风力机整机振动特性对研究风力机塔筒结构稳定性具有重要意义,为进一步研究塔筒减振提供技术支持。另外,风力机塔筒承受着复合荷载,主要包括风荷载、自身重量及其引起的偏心矩。这些载荷可导致机电设备的运行失效,也会显著增加塔架上的变形和应力。柔塔作为风力机的重要承载部件,其变形和振动会降低结构强度和稳定。在严重情况下,可能会导致疲劳失效,甚至塔架倒塌。因此,风力机柔塔的振动特性预估对于确定结构可靠性具有重要意义,同时也需要一种风力柔塔系统建模和振动特性快速计算方法。
风力机可归结为由多个刚体(机舱、轮毂等)和柔性体(叶片和塔架)组成的多刚柔体系统。有限元法(FEM)是描述柔性体变形的常用方法。已有文献对风力机塔筒进行了有限元建模,但在其模型中,刚体往往被简化为集中质量。有限元法虽然具有较强的通用性,但不能单独完成多刚柔体系统的动力学建模。此外,构造高阶连续场函数难度大,计算效率低。多体系统传递矩阵法(MSTMM)是芮筱亭院士提出的一种新的多体动力学方法,已应用于多种多刚柔体系统的动力学建模中。该方法可用于计算多体系统的振动特性和动力响应。MSTMM具有矩阵阶低、公式简单、计算成本低等优点。MSTMM在风力机柔塔中的应用,既能考虑柔性塔的几何非线性,又能考虑柔体变形与刚体运动的耦合。因此,该方法适用于风力机柔塔系统的建模和振动特性分析。
发明内容
针对现有技术中存在的不足,本发明提供了一种风力机柔塔系统固有振动特性快速建模及仿真方法,实现了风力机柔塔振动特性的快速建模和仿真,方便了风力机柔塔系统的振动特性分析。
本发明的目的是这样实现的:一种风力机柔塔系统固有振动特性快速建模方法,包括以下过程:
步骤一:对风力机柔塔系统进行简化,将系统离散成各多体动力学元件;
步骤二:基于多体系统传递矩阵法推导考虑重力作用的塔筒梁传递矩阵;
步骤三:确定其他多体动力学元件的传递矩阵;
步骤四:确定边界条件,建立柔塔系统总传递方程;
作为本发明的进一步限定,步骤一具体包括:忽略影响塔筒振动特性的次要部分,对风力机柔塔结构模型做出如下假设和简化:
(1)忽略塔筒模型扭转运动;
(2)假设风力机叶片、轮毂和机舱为一个存在质量偏心的空间刚体,忽略其中复杂部件的影响;
(3)忽略塔筒自身的结构阻尼;
(4)根据塔架的高质量比特性,忽略塔筒内部空气质量;
(5)假设塔筒内部法兰为集中质量,忽略螺栓、螺丝和螺孔造成的影响;
(6)假设塔筒截面始终保持圆形,不考虑壳体截面变形;
(7)假设塔筒表面平滑,不存在焊缝凹陷;
(8)假设塔筒为梁结构,假设不存在剪切效应和局部屈曲效应;
(9)假设塔筒与机舱为固定连接,忽略机舱相对于塔筒的扭转;
(10)忽略风力机内部爬梯、电缆架、平台等对塔架整体强度影响较小的复杂部件;
根据风力机柔塔系统的结构,将塔筒离散成多段考虑顶部构件压力作用、不同结构参数的Euler-Bernoulli梁;将风力机塔筒内部的法兰、螺栓和支撑其的平台处理为集中质量;将作用在塔筒顶部的风力机叶片、轮毂和机舱处理为一个存在质量偏心的空间刚体;
根据模型中所有元件的运动和相互作用,状态矢量可统一为Z=[Y,Z,Θy,Θz,My,Mz,Qy,Qz]T,其中,Z为沿坐标轴z位移对应的模态坐标列阵,Y为沿坐标轴y位移对应的模态坐标列阵,Θz为该点处相对于平衡位置相对于在z轴的角位移对应的模态坐标列阵,Θy为该点处相对于平衡位置相对于y轴的角位移对应的模态坐标列阵,Mz为沿坐标轴z内力矩对应的模态坐标列阵,Mx为沿坐标轴x内力矩对应的模态坐标列阵,Qz为沿坐标轴z内力对应的模态坐标列阵,Qy为沿坐标轴y内力对应的模态坐标列阵。
作为本发明的进一步限定,步骤二中所述推导塔筒的梁传递矩阵包括:
建立考虑重力的双自由度Euler-Bernoulli梁的振动微分方程:
式中,EI为梁的弯曲刚度,为单位长度质量,t为时间,z为塔筒流向位移,y为塔筒横向位移,x表示在x轴上的位置,G为梁轴向所受的压力;
由模态坐标转换,令
y(x,t)=Y(x)eiωt,z(x,t)=Z(x)eiωt (3)
式中,ω为角速度,Y(x)、Z(x)为y(x,t)、z(x,t)通过模态坐标转换得到,i代表虚数;
将式(3)其代入式(1)、(2),使偏微分方程变为4阶常微分方程得:
式(4-5)可计算得通解:
式中A1~A8为常数,并且
根据Euler-Bernoulli梁特性,可得到模态坐标系下的角位移Θy、Θz,力矩My、Mz和内力Qy、Qz的表达式:
将式(6)和式(9)写成矩阵的形式可得:
即Z(x)=B(x)·a,a=[A1,A2,A3,A4,A5,A6,A7,A8]T,其中Z(x)为状态矢量,即方程(10)等式左边项,B(x)代表方程(10)等式右边的第一个矩阵;令x=0代入方程(10)可得到Z(0)=B(0)·a;由于B-1(0)·Z(0)=a,因此有Z(x)=B(x)·B-1(0)·Z(0)=U(x)·Z(0);
所以该段Euler-Bernoulli梁从输入点到x长位置处的传递矩阵为:
U(x)=B(x)·B-1(0) (11)
此处,
根据式(11),可写出长度为x的双自由度Euler-Bernoulli梁传递矩阵形式为:
其中,
作为本发明的进一步限定,步骤三中确定其他多体动力学元件的传递矩阵包括:
法兰、螺栓和支撑其的平台被简化为集中质量,作用在塔筒内部;此类常规多体动力学元件,可以直接写出其柔塔系统坐标系中的传递矩阵形式:
式中,mi为该元件的质量,ω同样为角频率;
塔筒顶部支撑的风力机部件被简化为一空间刚体;作为常规多体动力学元件,根据状态矢量和坐标系,同样可直接写出质量为m的空间刚体传递矩阵形式:
式中,
其中J为刚体相对于输入点的惯量矩阵,下标代表转动轴。
作为本发明的进一步限定,步骤四中,所述建立柔塔系统总传递方程包括:
确认系统两端边界条件,塔筒底端为固定边界,位移和角位移都为0,则状态矢量为Z1,0=[0,0,0,0,My,Mz,Qy,Qz]T;顶部构件为自由振动,力和力矩都为0,则状态矢量为Zn+1,n=[Y,Z,Θy,Θz,0,0,0,0]T;根据MSTMM对链式多刚柔体系统的定义,按传递方向拼接各元件传递矩阵,写出系统总传递方程:
Zn+1,n=UnUn-1...Ui...U2U1Z1,0=UallZ1,0 (18)
式中,Ui下标为元件编号为i的元件传递矩阵,Uall=UnUn-1...Ui...U2U1为总传递矩阵。
一种风力机柔塔系统固有振动特性快速仿真方法,对所述建模方法得到的柔塔系统总传递方程求解,得到柔塔系统的振动特性。
作为本发明的进一步限定,所述求解总传递方程,得到柔塔系统的振动特性包括:假设Uall为:
将边界条件带入式(18),可得到:
根据式(19)中的线性齐次项,柔塔系统的特征方程可写成:
对方程(20)求解,即可求出柔塔系统固有频率;求解对应于固有频率的系统边界点状态矢量Z1,0、Zn+1,n,进而通过元件传递方程得到系统全部联接点的状态矢量,即为系统的振型。
与现有技术相比,本发明的有益效果在于:
(1)本发明提出一种风力机柔塔系统固有振动特性快速建模及仿真方法,推导过程简单且无需划分大量网格单元,省去了每改变一次参数就要对塔筒进行重新建模的过程;同时也能给工程上类似的细长柔性结构的动力学建模与快速仿真提供参考;
(2)本发明采用多体系统传递矩阵法,无需划分大量网格单元,具有矩阵阶次低、计算量小、计算效率高的优点,更能满足工程上快速计算的需求;
(3)本发明在风力机柔塔模型的建模中,考虑了顶部质量,考虑柔性塔截面非线性变化的特性,使模型更接近实际;同时,还考虑了柔塔与机舱等其它风力机部件之间的刚柔耦合效应。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1是本发明方法流程示意图。
图2是风力机柔塔系统简化物理模型示意图。
图3是柔塔的梁模型示意图。
图4是空间刚体模型。
图5是基于MSTMM计算出的系统圆频率计算结果图。
图6是计算出的柔塔第一阶振型(a)是z方向;(b)是y方向。
图7是计算出的柔塔第二阶振型(a)是z方向;(b)是y方向。
图8是计算出的柔塔第三阶振型(a)是z方向;(b)是y方向。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
现有技术中,多体系统传递矩阵法(MSTMM)由芮筱亭院士及其团队建立,用于多体动力学分析。该方法能实现一般多刚体、多刚柔体系统动力学研究,其高计算效率方便运用于工程中。运用于本发明中对风力机柔塔系统的多体动力学建模,可以实现振动特性的快速计算,本发明可为工程上类似的细长柔性结构的动力学建模与快速仿真提供参考。
本发明的一种风力机柔塔系统固有振动特性快速建模及仿真方法,参见图1所示。具体包括以下步骤:
步骤一:对风力机柔塔系统进行合理的简化,将系统离散成各多体动力学元件;
现有的三叶片水平轴风力机可归结为由多个刚体(机舱、轮毂等)和柔体(叶片和塔架)组成的多刚柔体系统。风力机塔筒底部靠桩土固定,塔筒顶部连接机舱,而机舱通过传动轴与外部的轮毂和叶片相连。为方便生产和运输,塔筒通常需要分成多段,在安装时通过法兰和螺栓分段连接。
考虑到风力机系统一个复杂的多刚柔体系统,忽略影响塔筒振动特性的次要部分,对风力机柔塔结构模型做出如下假设和简化:(1)忽略塔筒模型扭转运动;(2)假设风力机叶片、轮毂和机舱为一个存在质量偏心的空间刚体,忽略其中复杂部件的影响;(3)忽略塔筒自身的结构阻尼;(4)根据塔架的高质量比特性,忽略塔筒内部空气质量;(5)假设塔筒内部法兰为集中质量,忽略螺栓、螺丝和螺孔造成的影响;(6)假设塔筒截面始终保持圆形,不考虑壳体截面变形;(7)假设塔筒表面平滑,不存在焊缝凹陷;(8)假设塔筒为梁结构,假设不存在剪切效应和局部屈曲效应;(9)假设塔筒与机舱为固定连接,忽略机舱相对于塔筒的扭转;(10)忽略风力机内部爬梯、电缆架、平台等对塔架整体强度影响较小的复杂部件。
考虑柔塔系统各部件之间的相互作用,基于MSTMM将整个系统离散成各多体动力学元件,如图2所示。根据风力机柔塔系统的结构,将塔筒离散成多段考虑顶部构件压力作用、不同结构参数的Euler-Bernoulli梁;将风力机塔筒内部的法兰、螺栓和支撑其的平台处理为集中质量;将作用在塔筒顶部的风力机叶片、轮毂和机舱处理为一个存在质量偏心的空间刚体。建立坐标系在塔筒底部,以来流方向为z轴,以垂直来流方向、平行于地面的横向为y轴,以塔筒中心轴线为x轴。
根据MSTMM,定义传递方向为塔筒底部向顶端刚体方向,如图2中箭头指向,按传递方向对系统中的元件进行编号。根据模型中所有元件的运动和相互作用,状态矢量可统一为Z=[Y,Z,Θy,Θz,My,Mz,Qy,Qz]T。Z为沿坐标轴z位移对应的模态坐标列阵,Y为沿坐标轴y位移对应的模态坐标列阵,Θz为该点处相对于平衡位置相对于在z轴的角位移对应的模态坐标列阵,Θy为该点处相对于平衡位置相对于y轴的角位移对应的模态坐标列阵,Mz为沿坐标轴z内力矩对应的模态坐标列阵,Mx为沿坐标轴x内力矩对应的模态坐标列阵,Qz为沿坐标轴z内力对应的模态坐标列阵,Qy为沿坐标轴y内力对应的模态坐标列阵。
步骤二:基于MSTMM推导考虑重力作用的塔筒梁传递矩阵;
如图3所示为简化为Euler-Bernoulli梁元件的某段塔筒,该段梁长为l。图中G为受到的梁上方元件的压力,ZI为输入端状态矢量,ZO为输出端状态矢量,U(x)代表该段梁的传递矩阵。根据以上条件,建立考虑重力的双自由度Euler-Bernoulli梁的振动微分方程:
式中,EI为梁的弯曲刚度,为单位长度质量,t为时间,z为塔筒流向位移,y为塔筒横向位移,x表示在x轴上的位置。
由模态坐标转换,令
y(x,t)=Y(x)eiωt,z(x,t)=Z(x)eiωt (3)
式中,ω为角速度,Y(x)、Z(x)为y(x,t)、z(x,t)通过模态坐标转换得到。
将式(3)其代入式(1-2),使偏微分方程变为4阶常微分方程得:
式(4-5)可计算得通解:
式中A1~A8为常数,并且
根据Euler-Bernoulli梁特性,可得到模态坐标系下的角位移Θy、Θz,力矩My、Mz和内力Qy、Qz的表达式:
将式(6)和式(9)写成矩阵的形式可得:
即Z(x)=B(x)·a,a=[A1,A2,A3,A4,A5,A6,A7,A8]T,其中Z(x)为状态矢量,即方程(10)等式左边项,B(x)代表方程(10)等式右边的第一个矩阵。令x=0代入方程(10)可得到Z(0)=B(0)·a。由于B-1(0)·Z(0)=a,因此有Z(x)=B(x)·B-1(0)·Z(0)=U(x)·Z(0)。
所以该段Euler-Bernoulli梁从输入点到x长位置处的传递矩阵为:
U(x)=B(x)·B-1(0) (11)
此处,
根据式(11),可写出长度为x的双自由度Euler-Bernoulli梁传递矩阵形式为:
其中,
本发明中,由塔筒离散成的梁元件传递矩阵形式都为式(12)中的Ubeam,但各个梁元件结构参数不同,需根据各梁元件的参数填入矩阵中。
步骤三:确定其他多体动力学元件的传递矩阵;
法兰、螺栓和支撑其的平台被简化为集中质量,作用在塔筒内部。此类常规多体动力学元件,可以直接写出其柔塔系统坐标系中的传递矩阵形式:
式中,mi为该元件的质量,ω同样为角频率。
塔筒顶部支撑的风力机部件(机舱、轮毂、叶片等)被简化为一空间刚体,模型如图4所示。I(0,0,0)代表输入点坐标,O(b1,b2,b3)代表输出点坐标,C(c1,c2,c3)代表质心坐标。作为常规多体动力学元件,根据状态矢量和坐标系,同样可直接写出质量为m的空间刚体传递矩阵形式:
式中,
其中J为刚体相对于输入点的惯量矩阵,下标代表转动轴。
步骤四:确定边界条件,建立柔塔系统总传递方程;
确认系统两端边界条件,塔筒底端为固定边界,位移和角位移都为0,则状态矢量为Z1,0=[0,0,0,0,My,Mz,Qy,Qz]T;顶部构件为自由振动,力和力矩都为0,则状态矢量为Zn+1,n=[Y,Z,Θy,Θz,0,0,0,0]T。根据MSTMM对链式多刚柔体系统的定义,按传递方向拼接各元件传递矩阵,写出系统总传递方程:
Zn+1,n=UnUn-1...Ui...U2U1Z1,0=UallZ1,0 (18)
式中,Ui下标为元件编号为i的元件传递矩阵,Uall=UnUn-1...Ui...U2U1为总传递矩阵。
步骤五:求解总传递方程,得到柔塔系统的振动特性;
假设Uall为:
将边界条件带入式(18),可得到:
根据式(19)中的线性齐次项,柔塔系统的特征方程可写成:
方程(20)只与系统的结构参数和固有频率有关,对于实际的线性振动系统,式(20)必有非零解,即可求出柔塔系统固有频率。求解对应于固有频率的系统边界点状态矢量Z1,0、Zn+1,n,进而通过元件传递方程得到系统全部联接点的状态矢量,即为系统的振型。
通过该方法可以快速计算系统每个部件的参数对风力机柔塔振动特性,即振型和频率的影响。
实施例
针对型号为NORDEX S70/1500的三叶片水平轴风力机柔塔进行多体动力学建模并振动特性分析。得到的柔塔系统的固有圆频率计算结果如图5所示。图中横坐标为圆频率,纵坐标表示Δ值的大小,当Δ值接近于0时即可以求出圆频率。图中显示的竖线即搜根过程,每条竖线下对应的就是该阶模态的圆频率。图6为计算出的风力机柔塔的第一阶振型,(a)是z方向;(b)是y方向;图7为风力机柔塔的第二阶振型,(a)是z方向;(b)是y方向;图8为风力机柔塔的第三阶振型,(a)是z方向;(b)是y方向。
同样,本方法对其他风力机柔塔建立多体动力学模型并实现振动特性分析,只需修改简化多体动力学模型和其中具体参数即可。
上实施例的说明只是用于帮助理解本发明的方法及其核心思想。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以对本发明进行若干改进和修饰,这些改进和修饰也落入本发明权利要求的保护范围内。
Claims (2)
1.一种风力机柔塔系统固有振动特性快速建模方法,其特征在于,包括以下过程:
步骤一:对风力机柔塔系统进行简化,将系统离散成各多体动力学元件,具体包括:忽略影响塔筒振动特性的次要部分,对风力机柔塔结构模型做出如下假设和简化:
(1)忽略塔筒模型扭转运动;
(2)假设风力机叶片、轮毂和机舱为一个存在质量偏心的空间刚体,忽略其中复杂部件的影响;
(3)忽略塔筒自身的结构阻尼;
(4)根据塔架的高质量比特性,忽略塔筒内部空气质量;
(5)假设塔筒内部法兰为集中质量,忽略螺栓、螺丝和螺孔造成的影响;
(6)假设塔筒截面始终保持圆形,不考虑壳体截面变形;
(7)假设塔筒表面平滑,不存在焊缝凹陷;
(8)假设塔筒为梁结构,假设不存在剪切效应和局部屈曲效应;
(9)假设塔筒与机舱为固定连接,忽略机舱相对于塔筒的扭转;
(10)忽略风力机内部爬梯、电缆架、平台对塔架整体强度影响小的复杂部件;
根据风力机柔塔系统的结构,将塔筒离散成多段考虑顶部构件压力作用、不同结构参数的Euler-Bernoulli梁;将风力机塔筒内部的法兰、螺栓和支撑其的平台处理为集中质量;将作用在塔筒顶部的风力机叶片、轮毂和机舱处理为一个存在质量偏心的空间刚体;
根据模型中所有元件的运动和相互作用,状态矢量可统一为X=[Y,Z,Θy,Θz,My,Mz,Qy,Qz]T,其中,Z为沿坐标轴z位移对应的模态坐标列阵,Y为沿坐标轴y位移对应的模态坐标列阵,Θz为该点处相对于平衡位置相对于在z轴的角位移对应的模态坐标列阵,Θy为该点处相对于平衡位置相对于y轴的角位移对应的模态坐标列阵,Mz为沿坐标轴z内力矩对应的模态坐标列阵,Qz为沿坐标轴z内力对应的模态坐标列阵,Qy为沿坐标轴y内力对应的模态坐标列阵;
步骤二:基于多体系统传递矩阵法推导考虑重力作用的塔筒梁传递矩阵,包括:
建立考虑重力的双自由度Euler-Bernoulli梁的振动微分方程:
式中,EI为梁的弯曲刚度,为单位长度质量,t为时间,z为塔筒流向位移,y为塔筒横向位移,x表示在x轴上的位置,G为梁轴向所受的压力;
由模态坐标转换,令
y(x,t)=Y(x)eiωt,z(x,t)=Z(x)eiωt (3)
式中,ω为角速度,Y(x)、Z(x)为y(x,t)、z(x,t)通过模态坐标转换得到,i代表虚数;
将式(3)其代入式(1)、(2),使偏微分方程变为4阶常微分方程得:
式(4-5)可计算得通解:
式中A1~A8为常数,并且
根据Euler-Bernoulli梁特性,可得到模态坐标系下的角位移Θy、Θz,力矩My、Mz和内力Qy、Qz的表达式:
将式(6)和式(9)写成矩阵的形式可得:
即Z(x)=B(x)·a,a=[A1,A2,A3,A4,A5,A6,A7,A8]T,其中Z(x)为状态矢量,即方程(10)等式左边项,B(x)代表方程(10)等式右边的第一个矩阵;令x=0代入方程(10)可得到Z(0)=B(0)·a;由于B-1(0)·Z(0)=a,因此有Z(x)=B(x)·B-1(0)·Z(0)=U(x)·Z(0);
所以该段Euler-Bernoulli梁从输入点到x长位置处的传递矩阵为:
U(x)=B(x)·B-1(0) (11)
此处,
根据式(11),可写出长度为x的双自由度Euler-Bernoulli梁传递矩阵形式为:
其中,
步骤三:确定其他多体动力学元件的传递矩阵,包括:
法兰、螺栓和支撑其的平台被简化为集中质量,作用在塔筒内部;其柔塔系统坐标系中的传递矩阵形式:
式中,mi为该元件的质量,ω同样为角频率;
塔筒顶部支撑的风力机部件被简化为一空间刚体;根据状态矢量和坐标系,同样可直接写出质量为m的空间刚体传递矩阵形式:
式中,
其中J为刚体相对于输入点的惯量矩阵,下标代表转动轴;
步骤四:确定边界条件,建立柔塔系统总传递方程,所述建立柔塔系统总传递方程包括:
确认系统两端边界条件,塔筒底端为固定边界,位移和角位移都为0,则状态矢量为Z1,0=[0,0,0,0,My,Mz,Qy,Qz]T;顶部构件为自由振动,力和力矩都为0,则状态矢量为Zn+1,n=[Y,Z,Θy,Θz,0,0,0,0]T;根据MSTMM对链式多刚柔体系统的定义,按传递方向拼接各元件传递矩阵,写出系统总传递方程:
Zn+1,n=UnUn-1...Ui...U2U1Z1,0=UallZ1,0 (18)
式中,Ui下标为元件编号为i的元件传递矩阵,Uall=UnUn-1...Ui...U2U1为总传递矩阵。
2.一种风力机柔塔系统固有振动特性快速仿真方法,其特征在于,对权利要求1所述建模方法得到的柔塔系统总传递方程求解,得到柔塔系统的振动特性:假设Uall为:
将边界条件带入式(18),可得到:
根据式(19)中的线性齐次项,柔塔系统的特征方程可写成:
对方程(20)求解,即可求出柔塔系统固有频率;求解对应于固有频率的系统边界点状态矢量Z1,0、Zn+1,n,进而通过元件传递方程得到系统全部联接点的状态矢量,即为系统的振型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110710893.9A CN113408072B (zh) | 2021-06-25 | 2021-06-25 | 风力机柔塔系统固有振动特性快速建模及仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110710893.9A CN113408072B (zh) | 2021-06-25 | 2021-06-25 | 风力机柔塔系统固有振动特性快速建模及仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113408072A CN113408072A (zh) | 2021-09-17 |
CN113408072B true CN113408072B (zh) | 2023-10-13 |
Family
ID=77679423
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110710893.9A Active CN113408072B (zh) | 2021-06-25 | 2021-06-25 | 风力机柔塔系统固有振动特性快速建模及仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113408072B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117874943B (zh) * | 2023-12-21 | 2024-06-14 | 南京航空航天大学 | 基于柔顺工装优化布局的航空大型薄壁构件刚性强化方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104573255A (zh) * | 2015-01-22 | 2015-04-29 | 天津理工大学 | 一种基于改进多体系统传递矩阵的机械臂建模与求解方法 |
CN105156270A (zh) * | 2015-10-09 | 2015-12-16 | 上海电机学院 | 风力发电机独立变桨控制系统及方法 |
CN108733917A (zh) * | 2018-05-17 | 2018-11-02 | 桂林电子科技大学 | 一种基于传递矩阵的载货汽车动力学建模方法 |
CN109583107A (zh) * | 2018-12-06 | 2019-04-05 | 国电联合动力技术有限公司 | 柔性塔筒的低风速智能优化方法、装置及存储介质 |
CN110837676A (zh) * | 2019-11-05 | 2020-02-25 | 扬州大学 | 一种基于多体系统传递矩阵法的舵系统振动特性预测方法 |
CN110837678A (zh) * | 2019-11-08 | 2020-02-25 | 扬州大学 | 基于多体系统传递矩阵法的二元翼型频域颤振模型建模方法 |
CN110889169A (zh) * | 2019-11-22 | 2020-03-17 | 扬州大学 | 基于多体系统传递矩阵法的舵面系统非线性颤振模型建模方法 |
-
2021
- 2021-06-25 CN CN202110710893.9A patent/CN113408072B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104573255A (zh) * | 2015-01-22 | 2015-04-29 | 天津理工大学 | 一种基于改进多体系统传递矩阵的机械臂建模与求解方法 |
CN105156270A (zh) * | 2015-10-09 | 2015-12-16 | 上海电机学院 | 风力发电机独立变桨控制系统及方法 |
CN108733917A (zh) * | 2018-05-17 | 2018-11-02 | 桂林电子科技大学 | 一种基于传递矩阵的载货汽车动力学建模方法 |
CN109583107A (zh) * | 2018-12-06 | 2019-04-05 | 国电联合动力技术有限公司 | 柔性塔筒的低风速智能优化方法、装置及存储介质 |
CN110837676A (zh) * | 2019-11-05 | 2020-02-25 | 扬州大学 | 一种基于多体系统传递矩阵法的舵系统振动特性预测方法 |
CN110837678A (zh) * | 2019-11-08 | 2020-02-25 | 扬州大学 | 基于多体系统传递矩阵法的二元翼型频域颤振模型建模方法 |
CN110889169A (zh) * | 2019-11-22 | 2020-03-17 | 扬州大学 | 基于多体系统传递矩阵法的舵面系统非线性颤振模型建模方法 |
Non-Patent Citations (3)
Title |
---|
A computationally efficient modeling method for the vibration analyses of two-dimensional system structures using re duce d transfer matrix method for multibody system;Hanjing Lu et al.;《Journal of Sound and Vibration》;全文 * |
传递矩阵法在风力机塔架动力学分析中的应用;黄巍 等;《科学技术与工程》;第1548-1553页 * |
黄巍 等.传递矩阵法在风力机塔架动力学分析中的应用.《科学技术与工程》.2012,第1548-1553页. * |
Also Published As
Publication number | Publication date |
---|---|
CN113408072A (zh) | 2021-09-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | A Fourier series solution for the transverse vibration of rotating beams with elastic boundary supports | |
CN102759927B (zh) | 一种使用多级隔振平台提高光学载荷成像质量的方法 | |
Li et al. | Dynamic modeling and analysis of wind turbine drivetrain considering the effects of non-torque loads | |
Jokar et al. | Dynamic modeling and free vibration analysis of horizontal axis wind turbine blades in the flap-wise direction | |
CN113654756B (zh) | 一种海上浮式风机主动式实时混合模型试验方法 | |
CN113408072B (zh) | 风力机柔塔系统固有振动特性快速建模及仿真方法 | |
CN103678762A (zh) | 优化的复合材料机翼气动弹性风洞模型的缩比建模方法 | |
CN110889169A (zh) | 基于多体系统传递矩阵法的舵面系统非线性颤振模型建模方法 | |
CN110837678A (zh) | 基于多体系统传递矩阵法的二元翼型频域颤振模型建模方法 | |
Qian et al. | Active flutter suppression of a multiple-actuated-wing wind tunnel model | |
CN110516384B (zh) | 一种多输入多输出齿轮传动系统动力学建模方法 | |
Dai et al. | Wind-tunnel evaluation for an active sting damper using multimodal neural networks | |
Patil et al. | A scalable time-parallel solution of periodic rotor dynamics in X3D | |
Bernhammer et al. | Geometrically nonlinear structural modal analysis using fictitious masses | |
CN114154381A (zh) | 一种包含伺服传动链的机械机构谐振频率预测方法 | |
Tian et al. | Dynamical modeling and experimental validation for squeeze film damper in bevel gears | |
Su et al. | Analysis of dynamic characteristic for misalignment-spline gear shaft based on whole transfer matrix method | |
CN117236142B (zh) | 解算索杆系的有限质点分析方法、系统及介质 | |
Luo et al. | Static equilibrium form-finding analysis of cable-strut system based on nonlinear dynamic finite element method | |
CN104773305A (zh) | 一种非对称式飞机燃油系统模拟试验台结构及其设计方法 | |
CN109141872B (zh) | 基于复模态分析法的双排行星齿轮传动均载系数试验方法 | |
CN110990910A (zh) | 时程激励下线性耗能结构响应的快速迭代法 | |
Stabile et al. | Numerical and experimental validation of the breadboard model of a novel hexapod platform for high-performance micro-vibration mitigation | |
CN110704894B (zh) | 斜拉桥桥塔地震响应的计算方法 | |
CN104699913A (zh) | 一种风力发电机组主机架强度计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |