翼面颤振形状敏度分析与研究.pdfVIP

第十■届全国空气弹性学术交流会会议论文集 cARs一201l-2113 翼面颤振形状敏度分析与研究 杨佑绪，吴志刚，杨超 北京航空航天大学航空科学与工程学院，北京100191 摘要：差分敏度不能准确地反映翼面颤振响应随形状设计变量的非线性变化趋势且计算耗费巨大。针对这种 现象提出一种颤振形状敏度的解析计算方法，考虑的变量主要包括展弦比、梢根比、机翼面积和后掠角。利用 等效板模型，并结合一阶活塞理论求解质量、刚度和气动矩阵关于设计变量的解析导数。利用Lancaster伴随 法求得特征值敏度和颤振解析形状敏度。在颤振求解中采用状态空间形式运动方程，通过跟踪气动弹性极点在 复平面的轨迹确定颤振速度。通过和重新分析的比较，发现在基本设计点附近利用解析敏度的Taylor近似能够 准确预测颤振速度。 关键词：颤振；形状敏度；等效板模型；活塞理论；翼面设计； 对于现代飞行器，不论是从强度、控制还是性能的角度，柔性都成为设计人员关心的话题。为 了说明结构的变化对柔性的影响，可靠的气动弹性敏度求解成为关键。 根据变量类型的不同，一般将敏度分为大小、形状(构型)和拓扑三种。大小类型设计变量主 要有结构部件的横截面积、厚度，气动面弯度，控制系统传递函数的系数等；形状设计变量有展弦 比、梢根比、前缘后掠角、机翼面积，传递函数的阶次等；拓扑设计变量包括梁／肋数量，舵面的数 量等。 颤振形状敏度指出了颤振边界随形状参数变化的移动情况。在翼面早期设计阶段，借助于它设 计人员可以避免大量的重复分析。同时，敏度技术是基于近似概念的现代飞行器优化设计的基础。 颤振敏度是响应值敏度中较难的问题IlJ，颤振形状敏度更是如此。 敏度的求解方法主要有有限差分法，解析方法和介于两者之间的半解析方法。在翼面结构分析 中，结构的刚度和质量矩阵都是大小设计变量的线性函数，但却是形状设计变量的非线性函数【2，31， 差分法不能得到准确的形状敏度：同时，在差分敏度计算中，需要计算不同参数时的结构质量、刚 度矩阵以及气动矩阵，甚至需要重新划分计算网格、控制网格变形等。若采用差分法计算形状敏度， 计算耗费过于巨大；另外，差分敏度的准确性易受所选步长的影响。当步长较大时，容易产生截断 误差：步长过小又易产生舍去误差【4】。 国内，大量的研究者在从事气动弹性优化设计。高萍掣5】研究了翼面结构／颤振主动控制律一体 化设计，采用了颤振速度关于大小类型设计变量的解析敏度。林梦鹤等Il】推导了气动弹性剪裁中结 构响应敏度的解析表达式，设计变量考虑了复合材料结构的铺设角。刘湘宁等【6】在考虑结构几何非 线性和气动非线性的基础上，推导了大展弦比复合材料机翼颤振特性关于复合材料铺设角的敏度。 由前述可以看出，国内的敏度研究主要集中于大小类型敏度的求解，并未涉及形状敏度。 国外对翼面形状敏度进行了大量的研究。Kapania等人171计算了颤振速度关于形状设计变量的解 析敏度，但解析只限于气动部分，采用有限差分法计算结构矩阵的形状敏度。Livne【2】采用解析方法 求得机翼位移、应力和固有频率关于形状设计变量的解析敏度。Livne等人f8】采用完全解析的方法计 算了超音速壁板颤振动压关于形状和大小设计变量的敏度。但壁板结构简单，并没有将此方法推广 应用到复杂翼面结构。在上述研究中，虽然采用Lancaster的伴随法【9】求解特征值敏度。但实际上特 征值敏度的求解依赖于结构质量、刚度以及气动刚度和阻尼矩阵关于设计变量的敏度，而求解此类 系统矩阵敏度都是基于有限差分法进行。 面向设计的非定常气动力分析方法的发展明显滞后，这大大限制了气动弹性形状优化的发展【101。 49 杨佑绪，吴志刚等：翼面颤振形状敏度分析与研究 鉴于这种情况，本文采用简化的气动力计算理论．一阶活塞理论，并结合等效板结构计算模型¨--41， 求解真实翼面结构关于形状设计变量完全解析的颤振敏度，并和重新分析结果进行对比，验证解析 颤振形状敏度的准确性。 1翼面模型几何描述 为了求得气动和结构矩阵关于形状变量的解析敏度，对翼面构型形状进行参数化是必要的n51。 本文对翼面结构作以下简化：整个翼面采用一个翼型；梁的数量及弦向无量纲站位给定；肋的 数量及的展向