旋翼桨叶固有特性计算.docxVIP

旋翼桨叶固有特性计算 摘要：本文基于弹性桨叶动力学知识，将分别对挥舞，摆振，扭转进行模型建立；然后采用变量分离方法结合模态叠加法，计算额定转速状态，桨叶的挥舞前3阶、摆振前2阶和扭转1阶的固有特性；采用MATLAB编程带入数据进行求解分析。 关键词：弹性桨叶分，离变量法，模态叠加 引言 对于桨叶剖面重心无偏心的这种常规铰接式旋翼，桨叶不存在各自由度之间的耦合，或者各自由度间的耦合很小，故在动力特性分析时，可以将挥舞、摆振、扭转单独分析，分析结果的精度度误差不大。 模型建立及求解 弹性桨叶挥舞固有特性计算 挥舞运动微分方程 e e F F F r r r 图1.1 桨叶挥舞运动受力图 采用积分形式的牛顿法求解：半径为r处外段上各力及力矩平衡，如图1.1取外段处一微段dρ，建立对r处的力矩平衡方程： = 1 \* GB3 ① 惯性力Fi=mZρdρ = 2 \* GB3 ② 离心力FN=mΩ2ρdρ = 3 \* GB3 ③ 气动力Fz 力臂ρ-r Mr=rRFz 由工程梁理论，结构弯矩可以表示为： Mr=EIyZ 其中Z=?2 将(1.2)式带入(1.1)式，整理得， E （1.3） 方程(1.3)两边对r求两次导数得到 m （1.4） 即桨叶的挥舞运动微分方程。 模型求解 这是一个变微分方程，通过分离变量法求解。Z(r,t)的解可以表示为 Z （1.5） Znr,t表示第n阶固有振动，Wnr表示第n阶固有振型函数，Tnt可以 Wnr 满足 （1.6） （1.7） 那么满足边界条件（1.6）、（1.7）连续系统的通解，可以采用模态叠加法，表示为 全部模态函数的线性组合 Z （1.8） 将（1.8）式带入（1.4）式中得 n=1 （1.9） 令Fz=0 E （1.10） 则由（1.9）式可得到自由振动方程 n=1 （1.11） 将（1.9）、（1.10）式等号两边作用算子eR?Wi e （1.12） 得到 M （1.13） K （1.14） 式中Mn=eRmWn2dr为第n阶主质量； M （1.15） 固有频率 ω （1.16） 对于（1.8）式取前N个有限项作为近似解，则有 Z （1.17） 由于实际中无法得到真实振型函数Wnr，故在计算中代以满足分边界条件的N个假设模态函数w1r、w2r、w3r?w W （1.18） 其中c1n, w （1.19） （1.20） （1.21） 设Tnt=sinωβnt+?, K （1.22） 其中M=M1nM2nM3n n=1,2,3为特征向量。 M （1.23） K （1.24） 桨叶挥舞固有振动归结为求解（1.22）式广义特征值问题，可求得3个特征值ωβ2，即固有频率和3个特征向量C，则前三 W （1.25） 计算结果 将已知参数带入上述模型，利用Matlab很容易求得桨叶挥舞前三阶固有频率，和振型函数。 固有频率为：ωβ1=1.0152Ω， 振型函数为： W1 W2r=1.12416*(0.204082*r -0.0204082)^3 - 0.190552*r + 3.85872*(0.204082*r - 0.0204082)^4 W3r=0.528103*r- 42.1677*(0.204082*r - 0.0204082)^3 + 94.1201*(0.204082*r - 图1.2 弹性桨叶挥舞振动前三阶振型 摆振固有特性结算 摆振运动微分方程 FNF F F F r 图2.1 桨叶摆振运动受力图 同理采用积分形式牛顿法，建立摆振运动微分方程 = 1 \* GB3 ①惯性力Fi=mxρdρ = 2 \* GB3 ②离心力FN=mΩ2ρdρ 力臂 = 3 \* GB3 ③气动力FX，力臂ρ-r M （2.1） 整理得 m （2.2） 方程求解 X(r,t)的解可以表示为 X （2.3） Xnr,t表示第n阶固有振动，Wnr表示第n阶固有振型函数，Tnt Wnr 满足 （2.4） （2.5） 那么满足边界条件（2.4）、（2.5）连续系统的通解，可以采用模态叠加法，表示为 全部模态函数的线性组合 X （2.6） 采用相同的处理方法,得到摆振自由振动方程为 M （2.7） 其中 M （2.8） K （2.9） Mn为第n阶主质量；Kn为第n 同样可以假设 W （2.10） 其中c1n, w （2.11） （2.12） 设Tnt=sinωβnt+? K