一种快速获取复材翼面跨音速颤振速度的方法

专利名称：一种快速获取复材翼面跨音速颤振速度的方法
技术领域：
本发明是一种快速获取复材翼面跨音速颤振速度的方法，属于飞机复合材料翼面气动弹性剪裁领域，用于飞机复合材料翼面设计中的跨音速颤振速度优化。
背景技术：
翼面跨音速颤振是翼面设计的一种重要情况，在跨音速颤振分析过程中，求解得到跨音速颤振临界速度是重要步骤。一般的跨音速颤振计算采用非定常气动力方法，需要很长的计算时间，给翼面优化带来巨大的困难。当前国内外飞行器复合材料翼面跨音速颤振优化与分析尚未形成成熟的技术方案。过去大都采用实验和工程估算的方法得到翼面的跨音速颤振临界速度；近年来，随着计算流体动力学(CFD)、计算结构动力学(CSD)的发展和计算机性能的提高，将CFD和CSD结合起来进行复合材料翼面颤振的数值模拟有了相当大的进展，成为当前气动弹性研究的重要方向之一。和以往的方法相比，采用CFD/CSD耦合进行跨音速颤振求解的特点是精度较高，但其不足之处在于数值模拟过程的计算量很大， 耗时过长。在缩短CFD/CSD耦合计算的时间方面，除了依赖计算机硬件性能的提高、利用大型机参与计算外，还采用了计算机集群并行计算的方法。为设计出满足跨音速颤振要求的飞行器，需要对飞行器复合材料翼面结构进行跨音速颤振优化。优化中若直接采用求解非定常气动力的方法计算跨音速颤振临界速度，计算时间是极其惊人的。造成优化迭代的单步时间长达数月，这在工程上无法接受。利用响应面进行翼面优化，必须解决复合材料板参数和设计变量可行域描述。这样可以找到响应面的边界，实现高精度插值拟合。而常规复合材料板的描述方法无法满足上述要求。利用三向刚度方法描述板的力学特性，可以简化复合材料板的多种设计变量参数(包括各种角度和相应的铺层比例，铺层厚度等)，便于进行翼面优化设计。三向刚度特性是复合材料板元件固有的力学特性，对于一块给定厚度(层数)的复合材料平板，0°、士45°、90°铺层的比例直接影响到平板的轴向刚度特性EI1、EI2和剪切刚度特性GJ，如附图1所示。铺层角度、比例各异的复合材料板均可转化为相同厚度的、仅有0°、士45°、90°铺层的复合材料板，且保证板的三向刚度特性相同。由于复合材料板的三向刚度特性参数Ell、EI2、GJ之间存在约束关系，可用其中的两个参数表示另外一个参数，由此得到复合材料板元件刚度可行域的数学表述。具体的公式推导如下
2C(1 4- Cl丨(y.- — U,) — Gfl2El1 二-[(IZ1 ▲ U;) ’ 8 - G/] - Γ/ Β , CI )
1 I^C LV 1 ，.2[(IZ1 ^ υΒ) ·δ - β、J其中C = ，δ为板厚，U1, U4、U5为与材料特性相关的不变数，U1, U4, U5的计算
Ds i 2
方法如下首先计算纤维主向刚度系数Qh = ..,Q22 = Sas，Qfifi = G19Q12= μ 12Q22 = u21Qn，|l2S = 5lgil
其中Εη、Ε22、μ12、μ 21、G12为与材料相关的工程弹性常数则U, = ^(SQii + 3Q,2 + 2Q + 4Q.,)Ui = ^(.Qii 十 Q25 十 6Q” - 4Q.S!3)Us = ^Qii 十 Qs-将C值带入(I)，整理得到EIp EIjPGJ的约束关系，如(II)式所示。(EI1-EI2) 2GJ2+(C4-C6-C2) GJ+(C5+CrC3) = 0 (II)其中参数C1 C6的表达式如下，C1 = 2EI1EI2 (EI^EI2) (UfU5) δC2 = 2EI1EI2 (EI^EI2)C3 = 4EI1EI2 (UfU5)2 δ 2C4 = SEI1EI2(U^U5) δC5 = (EI^EI2)2 (U5-U4)2 δ 2C6 = 2 (EI^EI2)2 (U5-U4) δ复合材料平板三向刚度等效转换与刚度特性参数间约束关系式的推导为本方法中建立刚度可行域提供了理论基础，保证了刚度可行域的正确性与有效性。

发明内容
本发明正是针对上述现有技术存在的问题而设计提供了一种快速获取复材翼面跨音速颤振速度的方法，其目的是解决复合材料翼面跨音速颤振求解、分析与优化中复合材料板参数与设计变量可行域的描述问题，以及计算时间过长的问题。鉴于跨音速颤振优化能力的限制，本方法在亚音速颤振优化的基础上，选定有限个翼面敏感区称为敏感元件，利用三向刚度来描述敏感元件的力学特性，进而建立敏感元件的三向刚度可行域，在刚度可行域内选择特征点进行非定常气动力颤振计算，利用计算结果分别建立每个敏感元件与整个复合材料翼面的跨音速颤振临界速度响应面。在其后跨音速颤振求解、分析与优化的过程中，利用查询响应面的方式快速获取翼面的跨音速颤振速度，从而代替了长时间的非定常气动力计算，达到减少计算量，缩短计算时间的目的。本发明的目的是通过以下技术方案来实现的该种快速获取复材翼面跨音速颤振速度的方法，其特征在于该方法的步骤是(1)对飞机翼面蒙皮进行亚音速颤振优化，优化后上、下翼面蒙皮的有限元模型上，将铺层厚度δ相差5层以内的几个相邻设计单元组成一个敏感元件，单侧上、下翼面蒙皮上共选取2 10个敏感元件，左、右翼面蒙皮上选取的敏感元件对称，敏感元件包含的区域不重叠；(2)对于每个敏感元件，计算其所包含的所有设计单元的平均厚度，并将该敏感元件所包含的所有设计单元规整化。规整化的具体步骤为找出厚度与平均厚度最接近的设计单元，将同一敏感元件中其余设计单元的铺层去除，替换为该设计单元的铺层，并将该设计单元的厚度作为该敏感元件规整后的平均厚度d ；(3)对每一个敏感元件，建立刚度可行域并选取特征点
影响复合材料翼面跨音速颤振特性的变量是翼面蒙皮参数，通常包括铺层厚度 (层数)、铺层角度和各种铺层角之间的比例关系。在翼面跨音速颤振分析与优化的过程中，各种铺层之间比例的取值是不连续的，且铺层厚度(层数)、角度、比例的组合情况很多，参数变量的不同组合与跨音速颤振速度V。p之间很难找到特定的解析关系。鉴于上述原因，难以直接找到参数变量及变量组合的可行域，并建立参数变量与跨音速颤振速度V。p之间较为精确的响应面(全部为内插值的响应面)。本方法采用三向刚度力学参数代替复杂的铺层组合，建立敏感元件的三向刚度可行域，以可行域中的特征点与对应的颤振速度为基础建立该敏感元件与翼面跨音速颤振速度之间的响应面。根据各敏感元件自身特有的约束条件，建立不同的三向刚度可行域，进一步建立对应于每个敏感元件的跨音速颤振速度响应面。对每一个敏感元件，建立刚度可行域并选取特征点的过程如下3. 1确定在建立刚度可行域时，每层铺层可选用的铺层角θ在优化过程中，每层的铺层角可在0° 180°范围内任意选取。由于其他所有的铺层组合均可在厚度不变的情况下转换为0°、士45°和90°的铺层组合，且转换前后的刚度特性不变。因此在建立刚度可行域时，每层铺层可选用的铺层角θ为0°、士45°和 90°之间的一种，以此三种铺层角的组合代表该敏感元件所有可能的铺层状态，这同时也保证了工程上实施的可行性。3. 2确定每一个敏感元件的最小许用厚度dmin和最大许用厚度dmax本方法选用规整后的设计单元的平均厚度Cli为最小许用厚度dmin，根据敏感元件许用的结构空间和许用的结构应变水平，确定该敏感元件的最大许用厚度dmax。在建立刚度可行域时，利用最小许用厚度dmin和最大许用厚度dmax，可以保证可行域是有边界的，进而保证建立的跨音速颤振速度响应面是有边界的。3. 3根据每层铺层可选用的铺层角θ、敏感元件的最大许用厚度dmax和最小许用厚度dmin，确定每个敏感元件所有可能的铺层组合在实际应用中，由于工程设计上的要求，铺层板应为复合材料对称均衡板。即铺层角θ为士45°时，对应的铺层为4层的整数倍，铺层角θ为0°或90°时，对应的铺层为偶数层。在这种约束条件下，敏感元件的刚度可行域并非一个连续的曲面，而是由许多离散的点组成，且这些点所对应的三向刚度特性满足(II)式的约束关系，即这些离散点都在 (II)式表示的隐式曲面上。3. 4分别计算每个敏感元件在每种铺层角与铺层厚度组合下的三向刚度特性值 EIi、EI2、GJ，计算步骤如下首先计算纤维主向刚度系数(^.，二 li = li Q66 = G12Q12 = μ 12Q22 = υ 21Q12，P2.- =其中En、E22, μ 12、μ 21、G12为与材料相关的工程弹性常数然后计算层合板拉伸刚度矩阵，All =SLi^ ( ) k (i=l，2,6;j=l，2，6)其中N为层合板的层数，tk为单层厚度
对于每一层复合材料铺层，铺层角为θ时的刚度系数^计算公式，Qis = ra4Q:il + ZmsE2CQi2 ^ 2Q ) + H4Q22Q12 = Ri2R2CQil, + Q22 -Q2, = B4Q:,,+二 + 2QSi5)十 m^Q.,Qis = Ki'iitQ-a — Qi2) + SinHqt2 - Q::) — 2mn(m2 - n:)Q贴Qsg = KIII3(Qii-Q12) + m n(Q12 - Q22) + 2mn(m2 - n2)Q66Q6f5 = Krn^Qil+ Q,, -2Q12 - 2Qes>+(m4 + n4)QM其中m = cos θ，n = sin θ，θ为当前层的铺层角EI”EI2、GJ计算式如下EI1 = A11-A12^A12A22 ；EI2 = A22-A12^A12A11 ；GJ = A66 ；3. 5绘制刚度点，分别创建每一个敏感元件的刚度可行域步骤3. 4中，每一组EIp EI2, GJ数据作为敏感元件的一个刚度点，将敏感元件的所有刚度点绘制到EIpEI2、GJ构成的三维空间中，这些离散的刚度点即构成该敏感元件的刚度可行域，每个敏感元件的可行域所包络的空间均呈现类似三棱台的形状，在同一铺层厚度下，各刚度点在一帆状曲面上，从不同角度观察此帆状曲面如附图2、附图3、附图4 ；不同厚度的铺层对应的帆状曲面形状相似，彼此为比例缩放关系，缩放比为厚度比的四次方， 如附图5。3. 6对每一个敏感元件，分别在其刚度可行域内选取特征点在类似三棱台的刚度可行域的顶面、底面和内部等分剖面上选取特征点，选取位置为上述每面的三个顶点和三条边的中点，其中内部等分剖面的数量根据刚度可行域的大小选取，一般选取1 4个，特征点数目介于18 36个之间。如计算量可以承受，应考虑在上述各面上适量增选特征点，以提高跨音速颤振速度响应面的插值拟合精度。(4)对每个敏感元件，分别建立该敏感元件的跨音速颤振速度响应面一般情况下，翼面上会有多个敏感元件，对应每个敏感元件，有不同的刚度可行域和相应的特征点，因此需要分别针对每个敏感元件建立跨音速颤振速度响应面，其建立过程分为如下步骤4. 1创建刚度特征点对应的翼面有限元计算模型文件对于每个敏感元件，去除敏感元件上的原有铺层，用敏感元件的每一个特征点所代表的铺层铺覆在敏感元件上，翼面其他部位的铺层不变，由此建立敏感元件的每一个特征点对应的翼面有限元计算模型文件，每个敏感元件创建与特征点一一对应的18 36个有限元计算模型文件。如增选特征点，则对应增加相应的有限元计算模型文件；4. 2分别计算4. 1中每个翼面有限元计算模型的非定常气动力，得到相应的跨音速颤振临界速度V。p，如增选特征点，则v。p个数会对应增加；4. 3对每个敏感元件，分别建立跨音速颤振速度响应面将敏感元件的每个特征点的轴向刚度特性EI1、剪切刚度特性GJ和厚度δ作为自变量(18 36组或更多)，对应的翼面跨音速颤振速度Vcp (18 36个或更多)作为目标值，进行高维曲面拟合与插值，得到敏感元件的跨音速颤振速度响应面。由于另一方向的轴向刚度EI2与EI1和GJ存在(II)式所示的固定的约束关系，则EI2F需要作为自变量参与跨音速颤振速度响应面的建立。在响应面建立后，在敏感元件的刚度可行域内的任意一点，通过查询对应的响应面来代替复杂的非定常气动力计算，即可快速获得该状态下翼面的跨音速颤振速度V。p。本发明技术方案的优点和积极效果是1.利用三向刚度来描述敏感元件的铺层和质量特性，解决了如何建立覆盖复合材料板各种铺层角度与铺层比例可行域的技术难点；2.通过采用有限个特征点建立响应面的方法，保证了可行域里所有点进行跨音速颤振速度插值时，都落在所建立的响应面边界内，保证了较高的插值精度；3.利用查询响应面的方式代替复杂的非定常气动力计算用以获取翼面的跨音速颤振速度，大幅提高了计算速度，缩短了计算时间，解决了复合材料翼面跨音速颤振优化在工程上应用的难点；4.经过对某飞机复合材料翼面跨音速颤振优化计算的验证，建立的基于三向刚度的跨音速颤振速度响应面误差范围为士 lm/s，响应面比较光滑，表明了其鲁棒性好的特点， 在颤振分析与优化、以及实际工程上均具实用性。


图1为复合材料板三向刚度特性示意图
图2为某固定厚度下的刚度点构成的帆状曲面
图3为图2的另一视角
图4为图2的另一视角
图5为不同厚度下的刚度点构成的帆状曲面组
图6为某个敏感元件对应的刚度可行域与其中的特征点
图7为本发明实施例中上下翼面敏感元件的分布图
图8为最大许用厚度面上特征点的选取示意图
图9为响应面的误差分析曲线
图10为某固定厚度下响应面的剖面图
具体实施例方式以下将结合附图和实施例对本发明技术方案作进一步地详述以下将以某型飞机复合材料尾翼为例，介绍采用本发明所述快速获取复材翼面跨音速颤振速度的方法的实施过程，该过程的步骤是1.首先对飞机翼面蒙皮进行亚音速颤振优化，优化方法和过程采用在先申请的国防专利技术方案，该专利申请的名称是飞机翼面蒙皮亚音速颤振优化方法，专利申请号是201110171379. 9，该技术方案的主要内容是用翼面每块蒙皮元素作为优化对象，同时将蒙皮的铺层角度，铺层厚度增量、铺层位置选作设计变量，保证遗传算法的搜索是在包括了所有可能解的空间域内进行，采用双循环优化流程，解决了大设计变量(> 1000)的解耦问题，同时采用“小增量填谷法”有效的避免了算法的震荡，确保算法的高重复性。
在经过亚音速颤振优化后的水平尾翼上、下翼面蒙皮的有限元模型上，将铺层厚度δ相差5层以内的几个相邻设计单元1组成一个敏感元件2，对于本例来说，共选取8个敏感元件2，敏感元件2位于上下翼面前缘中央略偏翼尖处，以及后缘升降舵作动器附近， 其分布如附图7所示。平尾左、右翼面蒙皮上选取的敏感元件2对称，敏感元件2之间的区
域不重叠。 方法中，需针对每一个敏感元件2分别建立跨音速颤振速度的响应面。以下内容以其中一块敏感元件2为例，说明建立刚度可行域，选取特征点以及建立相应的跨音速颤振速度响应面的过程。2.分别计算该敏感元件2所包含的所有设计单元1的平均厚度，找出厚度与平均厚度最接近的设计单元1，将其余设计单元1的铺层去除并替换为该设计单元1的铺层，同时将该设计单元1的厚度作为该敏感元件2规整后的平均厚度d ；3.建立该敏感元件2的刚度可行域并选取特征点，其过程如下3. 1按照工程常用的情况，确定在建立刚度可行域3时，敏感元件每层铺层可选用的铺层角θ为0°、士45°和90°之间的一种；3. 2确定该敏感元件的最小许用厚度dmin和最大许用厚度dmax3. 2. 1确定最小许用厚度dmin该敏感元件2的最小许用厚度dmin为8层(单层0. 12mm)；3. 2. 2确定最大许用厚度dmax该敏感元件2的最大许用厚度dmax为56层(单层0. 12mm)；3. 3根据铺层角θ、敏感元件的最小许用厚度dmin和最大许用厚度dmax，确定所有可能的铺层组合；3. 4分别计算敏感元件2在每种铺层组合下的三向刚度EIp EI2, GJ，计算步骤如下首先计算纤维主向刚度系数(^.，Qii =,Qii = ’ Q66 = G12
Mil^i^"KDD丄二
__么，-·Ui2Q12= P12Q22= U21Qn^ii = ^T"其中En，E22, μ 12、μ 21、G12为与材料相关的工程弹性常数然后计算层合板拉伸刚度矩阵，Aii=E^tk ( ) (i=l,2’6;j=，2,6)
‘ 式其中N为层合板的层数，tk为单层厚度对于每一层复合材料铺层，铺层角为θ时的系数兔3的计算公式，Q1, = Jir4Q15 + ZmzJiHQ12 + 2Q } +Ii4Q22Qli = IB2 B1CQli +Q22 - 4Q:sa)+(mi + n4)Qi:Q,2 = Ji4Q11+ 2rasn2(Qls+ 2QS6)+ IB4Q12Qif3 =— Qi:)十Imnims-G2JquQ2s = !ImHQ1, - Qi2) + m3n(Q;i2- Q22) + 2mn(m2 -Qi5g = IB2IiHQiI + Q22 - 2Qi2 - 2q6S)+(m^ + n4)QM
其中m = cos θ，n = sine，θ为当前层的铺层角EI1, EI2、GJ计算式如下EI1 = A11-A12^A12A22 ；EI2 = A22-A12^A12A11 ；GJ = A66 ；对于本实例来说，该尾翼选用的复合材料的工程弹性常数为=E11 = 1. 29E+llPa, E22 = 9. 82E+9Pa，μ 12 = 0. 3056，G12 = 5. 29Ε+9Ι^，铺层的单层厚度 tk = 0. 12mm，敏感元件 2的最小许用dmin厚度为8层，最大许用厚度dmax为56层。3. 5步骤3. 4中，每一组EIpEI2、GJ数据做为敏感元件2的一个刚度点，将敏感元件2的所有刚度点绘制到EIpEI2、GJ构成的三维空间中，即构成敏感元件2的刚度可行域 3，刚度可行域3所包络的空间呈现类似三棱台的形状，参见附图6。3. 6在刚度可行域3内选取特征点在类似三棱台的刚度可行域3的顶面4(最小许用厚度面)、底面5(最大许用厚度面)和内部等分剖面6上选取特征点7，本实例中选择两个三等分剖面6，厚度分别为M层和40层，特征点7的选取位置为每个面的三个顶点和三条边的中点，参见附图6所示。在类似三棱台的刚度可行域3的顶面4、底面5和内部等分剖面6上，根据计算能力，可再增选特征点7。如本例中在底面5 (最大许用厚度面)的每条边上分别增选两个边等分点71，在对称轴线上增选三个轴等分点72，以及在内部1/4高处，左右对称增选两个内部点73，具体的增选情况参见附图8所示。特征点7总计达到35个。4.建立跨音速颤振速度响应面4. 1创建刚度特征点7对应的翼面有限元计算模型文件去除敏感元件2上规整后的铺层，用敏感元件2的每一个刚度特征点7所代表的铺层逐一铺覆在敏感元件2上，建立敏感元件2的每一个刚度特征点7对应的翼面有限元计算模型文件；4. 2分别计算4. 1中每一个翼面有限元计算模型的非定常气动力，得到翼面跨音速颤振速度V。p，特征点的轴向刚度EI1、剪切GJ，厚度δ及相应的翼面跨音速颤振速度V。p 如下表表1.特征点的刚度与跨音速颤振速度
权利要求
1. 一种快速获取复材翼面跨音速颤振速度的方法，其特征在于该方法的步骤是(1)对飞机翼面蒙皮进行亚音速颤振优化，在优化后上、下翼面蒙皮的有限元模型上， 将铺层厚度δ相差5层以内的几个相邻设计单元(1)组成一个敏感元件O)，单侧上、下翼面蒙皮上共选取2 10个敏感元件O)，左、右翼面蒙皮上选取的敏感元件( 对称，敏感元件( 包含的区域不重叠；(2)对于每个敏感元件O)，首先计算其所包含的所有设计单元(1)的平均厚度，并将该敏感元件( 所包含的所有设计单元(1)规整化。规整化的具体步骤为找出厚度与平均厚度最接近的设计单元(1)，将同一敏感元件中其余设计单元(1)的铺层去除，替换为该设计单元(1)的铺层，并将该设计单元(1)的厚度作为该敏感元件( 规整后的平均厚度 d；(3)对每个敏感元件O)，分别建立各自的刚度可行域C3)并选取特征点，其过程如下·3. 1对每个敏感元件O)，按照工程常用的情况，确定在建立刚度可行域C3)时，每层铺层可选用的铺层角θ为0°、士45°和90°之间的一种；·3. 2确定每个敏感元件O)的最小许用厚度dmin和最大许用厚度dmax 3. 2. 1确定最小许用厚度dmin对每个敏感元件O)，选用该敏感元件( 规整后的平均厚度d为最小许用厚度dmin ； 3. 2. 2确定最大许用厚度dmax对每个敏感元件O)，根据该敏感元件(2)许用的结构空间和许用的结构应变水平，确定该敏感元件O)的最大许用厚度dmax ；·3. 3根据每个敏感元件( 每层铺层可选用的铺层角θ、最大许用厚度(1_和最小许用厚度dmin，确定每个敏感元件( 各自的所有可能的铺层组合；·3. 4分别计算每个敏感元件( 在每种铺层角与铺层厚度组合下的三向刚度特性值 EI1、EI2、GJ，计算步骤如下，首先计算纤维主向刚度系数Aj,
2.根据权利要求1所述的快速获取复材翼面跨音速颤振速度的方法，其特征在于该铺层板应符合工程设计要求，为复合材料对称均衡板。即铺层角θ为士45°时，对应的铺层为4层的整数倍，铺层角θ为0°或90°时，对应的铺层为偶数层。
3.根据权利要求1所述的快速获取复材翼面跨音速颤振速度的方法，其特征在于对于每个敏感元件O)，在三棱台的刚度可行域C3)的顶面、底面( 和内部等分剖面(6) 上，根据计算能力，可再增选特征点(7)。
全文摘要
本发明是一种快速获取复材翼面跨音速颤振速度的方法，本方法在亚音速颤振优化的基础上，选定有限个翼面敏感区称为敏感元件，利用三向刚度来描述敏感元件的力学特性，进而建立敏感元件的三向刚度可行域，在刚度可行域内选择特征点进行非定常气动力颤振计算，利用计算结果分别建立每个敏感元件与整个复合材料翼面的跨音速颤振临界速度响应面。在其后跨音速颤振求解、分析与优化的过程中，利用查询响应面的方式快速获取翼面的跨音速颤振速度，从而代替了长时间的非定常气动力计算，达到减少计算量，缩短计算时间的目的。
文档编号G06F17/50GK102385651SQ201110182009
公开日2012年3月21日 申请日期2011年7月1日 优先权日2011年7月1日
发明者崔德刚, 张睿, 杜海, 熊青岳, 钱卫 申请人:中国航空工业集团公司科学技术委员会, 中航复合材料有限责任公司