发明名称 一种降低直升机旋翼桨毂振动载荷的方法
摘要 本发明公开了一种降低直升机旋翼浆毂振动载荷的方法,所提出的方法针对实际复合材料旋翼桨叶,桨叶截面具有C型翼梁、D型翼盒、前肋、后肋和蒙皮等结构,桨叶后端具有浆尖结构。桨叶优化设计以桨叶翼型典型剖面参数、复合材料铺层、翼尖后掠角和集中质量等作为设计变量,桨毂载荷作为目标函数,并以固有频率、气弹稳定性和自旋转动惯量作为约束条件建立优化模型。优化问题求解中使用了分级优化策略,并采用了复形法和改进的可行方向法相结合的算法。采用本发明提出的方法对旋翼浆毂进行减振优化设计能够显著降低桨毂载荷。
申请公布号 CN102722606B 申请公布日期 2014.05.14
申请号 CN201210164672.7 申请日期 2012.05.24
申请人 北京航空航天大学 发明人 向锦武;任毅如;罗漳平;张连鸿;郭俊贤;张亚军
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 姜荣丽
主权项 1.一种降低直升机旋翼桨毂振动载荷的方法,所述的方法用梁模型模拟实际的复合材料旋翼桨叶,该梁模型具有C型翼梁、D型翼盒、前肋、后肋和蒙皮结构,其特征在于:所述方法包括如下步骤:第一步:根据实际的旋翼桨叶结构建立相对应的桨叶优化模型;所述的桨叶优化模型中的目标函数为桨毂载荷f(D);旋翼桨叶设计变量包括剖面设计变量、铺层设计变量、集中质量设计变量和桨尖设计变量;剖面设计变量采用了如下所示V1设计变量确定,即桨叶在坐标系中的坐标值确定典型剖面,而典型剖面之间的剖面形状则通过线性连接得到;复合材料铺层设计变量为铺层角度和厚度;集中质量设计变量为各个质量块的质量及其展向位置,配重条的单位长度质量及其起始位置;浆尖设计变量包括后掠角、起始部位、尖削比和终止位置;V<sub>1</sub>=(x<sub>1</sub>,x<sub>2</sub>,y<sub>2</sub>,x<sub>3</sub>,X<sub>1</sub>,X<sub>2</sub>)其中x<sub>1</sub>,x<sub>2</sub>,y<sub>2</sub>和x<sub>3</sub>为C型翼梁上确定桨叶剖面形状的设计点坐标,X<sub>1</sub>和X<sub>2</sub>为前肋和后肋的位置,通过这几个参数能够确定桨叶剖面形状;约束函数为固有频率、桨叶惯量和气弹稳定性;具体如下:目标函数:Min(f(D))其中:f(D)=K<sub>F</sub>[(F<sub>x</sub>)<sup>2</sup>+(F<sub>y</sub>)<sup>2</sup>+(F<sub>z</sub>)<sup>2</sup>]<sup>1/2</sup>+K<sub>M</sub>[(M<sub>x</sub>)<sup>2</sup>+(M<sub>y</sub>)<sup>2</sup>+(M<sub>z</sub>)<sup>2</sup>]<sup>1/2</sup>或f(D)=K<sub>1</sub>F<sub>x</sub>+K<sub>2</sub>F<sub>y</sub>+K<sub>3</sub>F<sub>z</sub>固有频率约束条件:<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>g</mi><mrow><mi>i</mi><mrow><mo>(</mo><mi>D</mi><mo>)</mo></mrow></mrow><mi>U</mi></msubsup><mo>=</mo><mfrac><msub><mi>&omega;</mi><mi>i</mi></msub><msubsup><mi>&omega;</mi><mi>i</mi><mi>U</mi></msubsup></mfrac><mn>1</mn><mo>&le;</mo><mn>0</mn></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msubsup><mi>g</mi><mrow><mi>i</mi><mrow><mo>(</mo><mi>D</mi><mo>)</mo></mrow></mrow><mi>L</mi></msubsup><mo>=</mo><mn>1</mn><mo>-</mo><mfrac><msub><mi>&omega;</mi><mi>i</mi></msub><msubsup><mi>&omega;</mi><mi>i</mi><mi>L</mi></msubsup></mfrac><mo>&le;</mo><mn>0</mn></mrow></math>]]></maths>气弹稳定性约束条件:q(D)=ξ<sub>k</sub>+ε<sub>k</sub>≤0桨叶惯量约束:<img file="FDA0000459903850000013.GIF" wi="385" he="140" />其中F和M分别代表桨毂力和力矩,脚标x,y和z对应了三个方向上的分量,而K<sub>F</sub>,K<sub>M</sub>,K<sub>1</sub>,K<sub>2</sub>和K<sub>3</sub>分别为权重系数;i=1,2...m,ω<sub>i</sub>为结构固有频率,<img file="FDA0000459903850000014.GIF" wi="72" he="85" />为设计的频率上限,<img file="FDA0000459903850000015.GIF" wi="73" he="85" />为设计的频率下限,m为频率的阶数;k=1,2...m′,ξ<sub>k</sub>为悬停情况下第k阶模态特征值的实部,ε<sub>k</sub>为悬停情况下的可接受的k阶模态最小阻尼值,m′为模态数量;其中I<sub>0</sub>为直升机所必需的最小自旋目标值,I<sub>b</sub>为设计旋翼桨叶结构的自旋惯量值;第二步:求解桨叶优化模型,采用了分级优化策略,利用复形法和改进的可行方向法相结合的组合优化算法,具体如下:(1)给定原有旋翼桨叶设计方案,作为给定的当前设计方案;(2)计算当前设计方案的二维桨叶剖面特性,首先得到当前桨叶剖面的固有特性参数,从而得到桨叶惯量;然后通过求解由哈密尔顿原理建立一维梁的运动方程,得到桨叶固有频率及其振型;最后计算得到桨叶气弹响应和稳定性;(3)计算目标函数、约束条件:首先由桨叶气弹响应计算得到桨根部位的力,然后通过坐标变换得到桨毂载荷,约束条件函数有步骤(2)已得到;(4)采用线性近似方法对目标函数和约束条件进行近似,并求解近似问题,获得优化解;(5)若得到的优化解满足收敛条件,则结束,否则以得到的优化解作为新的设计方案重复(2)~(4)步,直到得到收敛的结果;第三步:输出桨叶设计参数,得到最优的旋翼桨叶结构模型。
地址 100191 北京市海淀区学院路37号