风力机桨叶模态分析的das实验研究.docxVIP

风力机桨叶模态分析的das实验研究 0 风力机桨叶模态分析的研究 由于风动机叶片上负荷的交叉性和随机性，其自身的弹性结构是不可避免的，叶片的振动是不可避免的。当叶片运行的旋转频率接近其固有频率时就会产生较大的动应力,发生共振。振动带来的疲劳会降低材料强度,缩短整机的使用寿命。因而对于桨叶的结构动力学特性分析是风力机研究工作的一项重要环节。它直接影响到风力机的性能和稳定可靠性。 叶片结构动力学研究主要应用的分析方法——模态分析是近年来进行结构动力学研究的主要手段,分为计算模态分析和试验模态分析。国内风力机领域在这方面的研究多是进行有限元分析及通过降阶进行数值求解的计算模态分析,实验模态分析的研究还比较缺乏。 实验模态分析是通过对输入和响应信号的参数识别获得模态参数——频率、阻尼比及振型—的实验方法。对风力机桨叶进行实验模态分析能较精确的确定风力机叶片的固有频率和振型,有助于对计算模态分析的计算模型和边界条件进行校核和修正。同时避免低阶固有频率与转速频率重合,为风力机桨叶的设计计算提供必要的实验参数依据。 1 单模态参数的识别 桨叶的内侧一端与轮毂的连接视为固定连接,另一端为自由端。在实验中水平固定安装桨叶内侧端,在自由尖端施加激振源。将桨叶简化为一根弹性简支梁,并处于小变形范围内,因而对桨叶仅考虑在挥舞方向(z方向)上的挥舞振动。 桨叶上某一点施加激振力Fj,桨叶各点的振动响应xi,以及任意两点之间复频响应函数hij之间的关系如下: [x1x2?xn]=[h11h12?h1nh21h22?h2n???hn1hn2?hnn]?[0?Fj?0](1)??????x1x2?xn??????=??????h11h21?hn1h12h22?hn2???h1nh2n?hnn???????????????0?Fj?0????????(1) 在实验中复频响应函数hij(ω)=xi(ω)Fj(ω)hij(ω)=xi(ω)Fj(ω)可以通过对所测得的激振力Fj和振动响应xi的值,进行参数识别而得到,其基本原理如下: 系统的传递函数矩阵定义为G(s)=(s2M+sC+K)-1,所以可以表示成: G(s)=ΦΦ-1(s2M+sC+K)-1(ΦT)-1ΦT =Φ[Φ-1(s2M+sC+K)Φ-1]ΦT =Φ(s2Mp+sCp+Kp)-1ΦT(2) =n∑i=1?i?ΤiΚpis2+Cpis+Κpi=∑i=1n?i?TiKpis2+Cpis+Kpi 其中Mp、Cp、Kp表示经过解偶变换之后的主质量矩阵、主阻尼矩阵和主刚度矩阵,Φ、?i表示振型矩阵和第i个振型向量。由(2)式得到复频响应函数。 Η(ω)=n∑i=1?i?ΤiΚpi-ω2Μpi+iωCpiH(ω)=∑i=1n?i?TiKpi?ω2Mpi+iωCpi(3) 复频响应函数的虚部和实部都独立地包含了系统的模态参数,由于虚部受模态的影响较实部要小而且在固有频率附近呈现明显的峰值,所以选择虚部进行参数识别。设实测的复频响应函数虚部是H(jω),要求第i阶模态参数,即Hi(jω)作为主导模态,设其它模态对它的影响为H0(jω),就可以将第i阶单模态复频响应函数的虚部简单地表示为 Ηi(jω)=Η(jω)-Η0(jω)=Cω-Dω3A-Bω2+ω4(4)Hi(jω)=H(jω)?H0(jω)=Cω?Dω3A?Bω2+ω4(4) 其中A、B、C、D为表征复模态参数的常数。但(4)式是位移响应的传递函数表达式,实验中所采用的是加速度传感器测响应,所以由两者之间的关系 H*i(jω)=Hi(jω)2·(jω)2=-Hi(jω)ω2 得到加速度响应传递函数的虚部表达式。 Η*i(jω)=Dω5-Cω3A-Bω2+ω4(5)H?i(jω)=Dω5?Cω3A?Bω2+ω4(5) 设对于采样频率ωk第i阶模态的估计值与真实值之间的误差是εk=Hi-Hik,对其加权得到ek,再对所有频率采样点(k=1,2,……,m)取均方误差为 E2=m∑k=1e2k=m∑k=1[εk(A-Bω2k+ω4k)]2=m∑k=1[Dω5k-Cω3k-Ηik(A-Bω2k+ω4k)]2(6)E2=∑k=1me2k=∑k=1m[εk(A?Bω2k+ω4k)]2=∑k=1m[Dω5k?Cω3k?Hik(A?Bω2k+ω4k)]2(6) 均方误差分别对A、B、C、D求偏导数并令其等于零就可得到这4个参数的值,再利用(3)式计算出模态参数阻尼比、振型向量等。 2 测试系统和设备 2.1 结构几何参数 本实验中选用玻璃钢复合纤维叶片,桨叶翼型为NACA632,长0.85m,由保定螺旋桨厂制造。风力机FD2-300,整机功率300W,额定转速400r/min。各种几何参数见表1,截面序号从根部开始,等间距展开。 考虑到桨叶从根部开始共有11个等距变化的不