基于Broyden法的旋翼多体系统气动弹性分析 Aeroelasticity Analysis of Rotor Multibody System Based on Broyden Method.pdfVIP

航 空 学报 Dec．252012V01．33No．122171．2182 ActaAeronautica etAstronauticaSinica ISSN1000—6893 CN¨-1929／V http：Hhkxb．buaa．edu．Cflhkxb@buaa．edu．cn 71—12 文章编号：1000—6893(2012)12—21 基于Broyden法的旋翼多体系统气动弹性分析 虞志浩，杨卫东*，张呈林 南京航空航天大学直升机旋翼动力学国家级重点实验室，江苏南京 210016 摘 要：建立了旋翼多体系统气动弹性模型并给出了一种适合于该模型响应计算的数值计算方法。采用柔性多体系统 动力学方法建立旋翼气动弹性模型，利用驱动约束显著简化约束方程形式，集成大变形桨叶模型，准确考虑变形的非线 性，适合于对采用柔性结构的先进旋翼进行气动弹性分析。基于Broyden法改进隐式积分法积分一步中非线性方程的 求解，避免求取切线矩阵和矩阵求逆运算，保持隐式积分法具有较好稳定性的同时提高计算效率，解决了旋翼多体系统 气动弹性力学方程隐式表达且具有较强非线性和较高刚性比造成的响应计算困难。通过模型旋翼桨叶响应计算验证了 结构模型与气动弹性响应求解方法。采用建立的气动弹性模型计算悬停和前飞状态旋翼气动弹性稳定性，与试验结果 对比验证了模型的正确性。研究了不同的稳定性计算方法、桨叶结构模型和入流模型等对悬停和前飞稳定性计算的影 响，结果表明本文所采用的结构、气动模型及气动弹性稳定性计算方法提高了气动弹性稳定性分析精度。 关键词：旋翼；气动弹性；Broyden法；多体系统动力学；稳定性 中图分类号：V211．47文献标识码：A 随着对直升机性能要求的提高，旋翼作为影 法[101”]，通过直接积分计算旋翼响应，由响应信 响直升机性能的主要部件，结构不断改进与创新， 号分析旋翼气动弹性特性。相比采用模态坐标变 球柔旋翼、无铰式旋翼和无轴承旋翼等先进构型 换并截取低阶模态的方法n3。¨]，瞬态响应分析法 旋翼的技术不断完善，应用日益广泛。这些先进 保留了结构模型中高频和非线性对气动弹性分析 构型旋翼桨叶具有显著的弹性变形，几何非线性 的影响，提高了分析精度。 与结构耦合对旋翼气动弹性分析有显著影响。为 由于旋翼气动弹性方程刚性比较高[15‘17]，响 准确预估旋翼的气动弹性稳定性、减小振动载荷， 应求解一般采用稳定性较好的隐式积分法，如平 近年来在气动弹性模型中开始采用柔性多体系统 动力学方法u‘4J，集成大变形桨叶模型H_9]、非定常分中非线性方程采用Newton法求解。这些方法 翼型气动力模型和非定常人流模型口“¨3等精度 适于求解二阶常微分方程，且需要在积分一步中 更高的模型。这些更为精细的结构、气动模型的 反复计算非线性方程求解所需的Jacobi矩阵，因 方程一般以隐式表达，在提高气动弹性分析精度 而求解含动力人流方程等一阶微分方程在内的旋 的同时，也使得旋翼气动弹性方程具有较强的非 翼气动弹性方程时有困难[1{]。此外与可以写出 线性和较高的刚性比。为准确考虑非线性和高频 解析切线矩阵的中等变形梁模型、准定常气动力模 成分的影响，气动弹性分析采用瞬态响应分析 型不同，高精度的结构、气动等模型采用隐式表达， 13：02 收稿日期：2012—03—02；退修日期：2012．04—27；录用日期：2012一06—11；网络出版时间：2012—08．30 网络出版地址：WWW，cnki．fiet／kcms／detail／11．1929．V13