• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 20
    • 下载费用:30 金币  

    重庆时时彩群骗局: 一种快速获取复材翼面跨音速颤振速度的方法.pdf

    关 键 词:
    一种 快速 获取 复材翼面 跨音速 速度 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201110182009.5

    申请日:

    2011.07.01

    公开号:

    CN102385651A

    公开日:

    2012.03.21

    当前法律状态:

    授权

    有效性:

    有权

    法律详情: 授权|||实质审查的生效IPC(主分类):G06F 17/50申请日:20110701|||公开
    IPC分类号: G06F17/50 主分类号: G06F17/50
    申请人: 中国航空工业集团公司科学技术委员会; 中航复合材料有限责任公司
    发明人: 崔德刚; 杜海; 张睿; 熊青岳; 钱卫
    地址: 100012 北京市朝阳区安定门外北苑2号院
    优先权:
    专利代理机构: 中国航空专利中心 11008 代理人: 陈宏林
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201110182009.5

    授权公告号:

    102385651B||||||

    法律状态公告日:

    2013.05.15|||2012.06.06|||2012.03.21

    法律状态类型:

    授权|||实质审查的生效|||公开

    摘要

    本发明是一种快速获取复材翼面跨音速颤振速度的方法,本方法在亚音速颤振优化的基础上,选定有限个翼面敏感区称为敏感元件,利用三向刚度来描述敏感元件的力学特性,进而建立敏感元件的三向刚度可行域,在刚度可行域内选择特征点进行非定常气动力颤振计算,利用计算结果分别建立每个敏感元件与整个复合材料翼面的跨音速颤振临界速度响应面。在其后跨音速颤振求解、分析与优化的过程中,利用查询响应面的方式快速获取翼面的跨音速颤振速度,从而代替了长时间的非定常气动力计算,达到减少计算量,缩短计算时间的目的。

    权利要求书

    1.一种快速获取复材翼面跨音速颤振速度的方法,其特征在于:该
    方法的步骤是:
    (1)对飞机翼面蒙皮进行亚音速颤振优化,在优化后上、下翼面蒙皮
    的有限元模型上,将铺层厚度δ相差5层以内的几个相邻设计单元(1)组
    成一个敏感元件(2),单侧上、下翼面蒙皮上共选取2~10个敏感元件(2),
    左、右翼面蒙皮上选取的敏感元件(2)对称,敏感元件(2)包含的区域不
    重叠;
    (2)对于每个敏感元件(2),首先计算其所包含的所有设计单元(1)的
    平均厚度,并将该敏感元件(2)所包含的所有设计单元(1)规整化。规整
    化的具体步骤为:找出厚度与平均厚度最接近的设计单元(1),将同一敏
    感元件中其余设计单元(1)的铺层去除,替换为该设计单元(1)的铺层,
    并将该设计单元(1)的厚度作为该敏感元件(2)规整后的平均厚度d;
    (3)对每个敏感元件(2),分别建立各自的刚度可行域(3)并选取特征
    点,其过程如下:
    3.1对每个敏感元件(2),按照工程常用的情况,确定在建立刚度可
    行域(3)时,每层铺层可选用的铺层角θ为0°、±45°和90°之间的一
    种;
    3.2确定每个敏感元件(2)的最小许用厚度dmin和最大许用厚度dmax
    3.2.1确定最小许用厚度dmin
    对每个敏感元件(2),选用该敏感元件(2)规整后的平均厚度d为最
    小许用厚度dmin;
    3.2.2确定最大许用厚度dmax
    对每个敏感元件(2),根据该敏感元件(2)许用的结构空间和许用的
    结构应变水平,确定该敏感元件(2)的最大许用厚度dmax;
    3.3根据每个敏感元件(2)每层铺层可选用的铺层角θ、最大许用厚
    度dmax和最小许用厚度dmin,确定每个敏感元件(2)各自的所有可能的铺层
    组合;
    3.4分别计算每个敏感元件(2)在每种铺层角与铺层厚度组合下的
    三向刚度特性值EI1、EI2、GJ,计算步骤如下,
    首先计算纤维主向刚度系数Qij,
    Q66=G12
    Q12=μ12Q22=υ21Q11,
    其中E11、E22、μ12、μ21、G12为与材料相关的工程弹性常数
    然后计算层合板拉伸刚度矩阵,
    A ij = Σ k = 1 N t k ( Q ij ) k , ( i = 1,2,6 ; j = 1,2,6 ) ]]>
    其中N为层合板的层数,tk为单层厚度
    对于每层复合材料铺层,铺层角为θ时的刚度系数计算公式,
    Q 11 = m 4 Q 11 + 2 m 2 n 2 ( Q 12 + 2 Q 66 ) + n 4 Q 22 ]]>
    Q 12 = m 2 n 2 ( Q 11 + Q 22 - 4 Q 66 ) + ( m 4 + n 4 ) Q 12 ]]>
    Q 22 = n 4 Q 11 + 2 m 2 n 2 ( Q 12 + 2 Q 66 ) + m 4 Q 22 ]]>
    Q 16 = m 2 n ( Q 11 - Q 12 ) + m n 2 ( Q 12 - Q 22 ) - 2 mn ( m 2 - n 2 ) Q 66 ]]>
    Q 26 = mn 2 ( Q 11 - Q 12 ) + m 2 n ( Q 12 - Q 22 ) + 2 mn ( m 2 - n 2 ) Q 66 ]]>
    Q 66 = m 2 n 2 ( Q 11 + Q 22 - 2 Q 12 - 2 Q 66 ) + ( m 4 + n 4 ) Q 66 ]]>
    其中m=cosθ,n=sinθ,θ为当前层的铺层角
    板的三向刚度EI1、EI2、GJ值的计算式如下
    EI1=A11-A12*A12/A22;EI2=A22-A12*A12/A11;GJ=A66;
    3.5在三维空间中绘制刚度点,分别创建每个敏感元件(2)的刚度可
    行域,其操作过程如下:步骤3.4中,每个敏感元件(2)对应的每一组EI1、
    EI2、GJ数据做为该敏感元件(2)的一个刚度点,将该敏感元件(2)的所有
    刚度点绘制到EI1、EI2、GJ构成的三维空间中,这些离散的刚度点即构
    成该敏感元件(2)的刚度可行域(3),每个敏感元件(2)对应的刚度可行域
    所包络的空间均呈现类似三棱台的形状;
    3.6对每个敏感元件(2),分别在其刚度可行域(3)内选取特征点(7)
    在每个敏感元件(2)对应的类似三棱台的刚度可行域(3)的顶面(4)、
    底面(5)和内部等分剖面(6)上选取特征点(7),选取位置为上述每面的三
    个顶点和三条边的中点,其中内部等分剖面(6)的数量根据刚度可行域(3)
    的大小选取,一般选取1~4个,特征点(7)数目介于18~36个之间;
    (4)对于每个敏感元件(2),分别建立跨音速颤振速度响应面
    4.1创建刚度特征点(7)对应的翼面有限元计算模型文件
    对于每个敏感元件(2),去除敏感元件(2)上的原有铺层,用敏感元
    件(2)的每个特征点(7)所代表的铺层铺覆在敏感元件(2)上,翼面其他部
    位的铺层不变,由此建立该敏感元件(2)的每个特征点(7)对应的翼面有
    限元计算模型文件,每个敏感元件(2)创建与其特征点(7)一一对应的有
    限元计算模型文件;
    4.2分别计算4.1中每个有限元计算模型的非定常气动力,得到相
    应的翼面跨音速颤振临界速度Vcp;
    4.3对每个敏感元件(2),分别建立跨音速颤振速度的响应面
    将敏感元件(2)的每个特征点(7)的轴向刚度特性EI1、剪切刚度特性
    GJ和厚度δ作为自变量(18~36组),对应的翼面跨音速颤振速度Vcp
    (18~36个)作为目标值,进行高维曲面拟合与插值,得到敏感元件(2)
    的跨音速颤振速度响应面(8)。在响应面建立后,在敏感元件(2)刚度可
    行域(3)内的任意一点,通过查询对应的响应面(8)即可快速获得该状态
    下的翼面跨音速颤振速度Vcp。
    2.根据权利要求1所述的快速获取复材翼面跨音速颤振速度的方法,
    其特征在于:该铺层板应符合工程设计要求,为复合材料对称均衡板。
    即铺层角θ为±45°时,对应的铺层为4层的整数倍,铺层角θ为0°或
    90°时,对应的铺层为偶数层。
    3.根据权利要求1所述的快速获取复材翼面跨音速颤振速度的方法,
    其特征在于:对于每个敏感元件(2),在三棱台的刚度可行域(3)的顶面
    (4)、底面(5)和内部等分剖面(6)上,根据计算能力,可再增选特征点(7)。

    说明书

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

    技术领域

    本发明是一种快速获取复材翼面跨音速颤振速度的方法,属于飞机
    复合材料翼面气动弹性剪裁领域,用于飞机复合材料翼面设计中的跨音
    速颤振速度优化。

    背景技术

    翼面跨音速颤振是翼面设计的一种重要情况,在跨音速颤振分析过
    程中,求解得到跨音速颤振临界速度是重要步骤。一般的跨音速颤振计
    算采用非定常气动力方法,需要很长的计算时间,给翼面优化带来巨大
    的困难。当前国内外飞行器复合材料翼面跨音速颤振优化与分析尚未形
    成成熟的技术方案。过去大都采用实验和工程估算的方法得到翼面的跨
    音速颤振临界速度;近年来,随着计算流体动力学(CFD)、计算结构动力
    学(CSD)的发展和计算机性能的提高,将CFD和CSD结合起来进行复合材
    料翼面颤振的数值模拟有了相当大的进展,成为当前气动弹性研究的重
    要方向之一。和以往的方法相比,采用CFD/CSD耦合进行跨音速颤振求
    解的特点是精度较高,但其不足之处在于数值模拟过程的计算量很大,
    耗时过长。在缩短CFD/CSD耦合计算的时间方面,除了依赖计算机硬件
    性能的提高、利用大型机参与计算外,还采用了计算机集群并行计算的
    方法。

    为设计出满足跨音速颤振要求的飞行器,需要对飞行器复合材料翼
    面结构进行跨音速颤振优化。优化中若直接采用求解非定常气动力的方
    法计算跨音速颤振临界速度,计算时间是极其惊人的。造成优化迭代的
    单步时间长达数月,这在工程上无法接受。

    利用响应面进行翼面优化,必须解决复合材料板参数和设计变量可
    行域描述。这样可以找到响应面的边界,实现高精度插值拟合。而常规
    复合材料板的描述方法无法满足上述要求。利用三向刚度方法描述板的
    力学特性,可以简化复合材料板的多种设计变量参数(包括各种角度和
    相应的铺层比例,铺层厚度等),便于进行翼面优化设计。

    三向刚度特性是复合材料板元件固有的力学特性,对于一块给定厚
    度(层数)的复合材料平板,0°、±45°、90°铺层的比例直接影响到
    平板的轴向刚度特性EI1、EI2和剪切刚度特性GJ,如附图1所示。铺层
    角度、比例各异的复合材料板均可转化为相同厚度的、仅有0°、±45
    °、90°铺层的复合材料板,且保证板的三向刚度特性相同。由于复合
    材料板的三向刚度特性参数EI1、EI2、GJ之间存在约束关系,可用其中
    的两个参数表示另外一个参数,由此得到复合材料板元件刚度可行域的
    数学表述。具体的公式推导如下:

    EI 1 = 2 C 1 + C [ ( U 1 + U 5 ) · δ - GJ ] - ( 1 + C ) [ ( U 5 - U 4 ) - GJ ] 2 2 [ ( U 1 + U 5 ) · δ - GJ ] - - - ( 1 ) ]]>

    其中δ为板厚,U1、U4、U5为与材料特性相关的不变
    数,U1、U4、U5的计算方法如下:

    首先计算纤维主向刚度系数Qij

    Q66=G12

    Q12=μ12Q22=υ21Q11,

    其中E11、E22、μ12、μ21、G12为与材料相关的工程弹性常数

    U 1 = 1 8 ( 3 Q 11 + 3 Q 22 + 2 Q 12 + 4 Q 66 ) ]]>

    U 4 = 1 8 ( Q 11 + Q 22 + 6 Q 12 - 4 Q 66 ) ]]>

    U 5 = 1 8 ( Q 11 + Q 22 - 2 Q 12 + 4 Q 66 ) ]]>

    将C值带入(I),整理得到EI1、EI2和GJ的约束关系,如(II)式所示。

    (EI1-EI2)2GJ2+(C4-C6-C2)GJ+(C5+C1-C3)=0?????(II)

    其中参数C1~C6的表达式如下,

    C1=2EI1EI2(EI1+EI2)(U1+U5)δ

    C2=2EI1EI2(EI1+EI2)

    C3=4EI1EI2(U1+U5)2δ2

    C4=8EI1EI2(U1+U5)δ

    C5=(EI1+EI2)2(U5-U4)2δ2

    C6=2(EI1+EI2)2(U5-U4)δ

    复合材料平板三向刚度等效转换与刚度特性参数间约束关系式的推
    导为本方法中建立刚度可行域提供了理论基础,保证了刚度可行域的正
    确性与有效性。

    发明内容

    本发明正是针对上述现有技术存在的问题而设计提供了一种快速获
    取复材翼面跨音速颤振速度的方法,其目的是解决复合材料翼面跨音速
    颤振求解、分析与优化中复合材料板参数与设计变量可行域的描述问题,
    以及计算时间过长的问题。

    鉴于跨音速颤振优化能力的限制,本方法在亚音速颤振优化的基础
    上,选定有限个翼面敏感区称为敏感元件,利用三向刚度来描述敏感元
    件的力学特性,进而建立敏感元件的三向刚度可行域,在刚度可行域内
    选择特征点进行非定常气动力颤振计算,利用计算结果分别建立每个敏
    感元件与整个复合材料翼面的跨音速颤振临界速度响应面。在其后跨音
    速颤振求解、分析与优化的过程中,利用查询响应面的方式快速获取翼
    面的跨音速颤振速度,从而代替了长时间的非定常气动力计算,达到减
    少计算量,缩短计算时间的目的。

    本发明的目的是通过以下技术方案来实现的:

    该种快速获取复材翼面跨音速颤振速度的方法,其特征在于:该方
    法的步骤是:

    (1)对飞机翼面蒙皮进行亚音速颤振优化,优化后上、下翼面蒙皮的
    有限元模型上,将铺层厚度δ相差5层以内的几个相邻设计单元组成一个
    敏感元件,单侧上、下翼面蒙皮上共选取2~10个敏感元件,左、右翼
    面蒙皮上选取的敏感元件对称,敏感元件包含的区域不重叠;

    (2)对于每个敏感元件,计算其所包含的所有设计单元的平均厚度,
    并将该敏感元件所包含的所有设计单元规整化。规整化的具体步骤为:
    找出厚度与平均厚度最接近的设计单元,将同一敏感元件中其余设计单
    元的铺层去除,替换为该设计单元的铺层,并将该设计单元的厚度作为
    该敏感元件规整后的平均厚度d;

    (3)对每一个敏感元件,建立刚度可行域并选取特征点

    影响复合材料翼面跨音速颤振特性的变量是翼面蒙皮参数,通常包
    括铺层厚度(层数)、铺层角度和各种铺层角之间的比例关系。在翼面跨
    音速颤振分析与优化的过程中,各种铺层之间比例的取值是不连续的,
    且铺层厚度(层数)、角度、比例的组合情况很多,参数变量的不同组合
    与跨音速颤振速度Vcp之间很难找到特定的解析关系。鉴于上述原因,难
    以直接找到参数变量及变量组合的可行域,并建立参数变量与跨音速颤
    振速度Vcp之间较为精确的响应面(全部为内插值的响应面)。

    本方法采用三向刚度力学参数代替复杂的铺层组合,建立敏感元件
    的三向刚度可行域,以可行域中的特征点与对应的颤振速度为基础建立
    该敏感元件与翼面跨音速颤振速度之间的响应面。根据各敏感元件自身
    特有的约束条件,建立不同的三向刚度可行域,进一步建立对应于每个
    敏感元件的跨音速颤振速度响应面。

    对每一个敏感元件,建立刚度可行域并选取特征点的过程如下:

    3.1确定在建立刚度可行域时,每层铺层可选用的铺层角θ

    在优化过程中,每层的铺层角可在0°~180°范围内任意选取。由
    于其他所有的铺层组合均可在厚度不变的情况下转换为0°、±45°和
    90°的铺层组合,且转换前后的刚度特性不变。因此在建立刚度可行域
    时,每层铺层可选用的铺层角θ为0°、±45°和90°之间的一种,以
    此三种铺层角的组合代表该敏感元件所有可能的铺层状态,这同时也保
    证了工程上实施的可行性。

    3.2确定每一个敏感元件的最小许用厚度dmin和最大许用厚度dmax

    本方法选用规整后的设计单元的平均厚度di为最小许用厚度dmin,根
    据敏感元件许用的结构空间和许用的结构应变水平,确定该敏感元件的
    最大许用厚度dmax。在建立刚度可行域时,利用最小许用厚度dmin和最大
    许用厚度dmax,可以保证可行域是有边界的,进而保证建立的跨音速颤振
    速度响应面是有边界的。

    3.3根据每层铺层可选用的铺层角θ、敏感元件的最大许用厚度
    dmax和最小许用厚度dmin,确定每个敏感元件所有可能的铺层组合

    在实际应用中,由于工程设计上的要求,铺层板应为复合材料对称
    均衡板。即铺层角θ为±45°时,对应的铺层为4层的整数倍,铺层角
    θ为0°或90°时,对应的铺层为偶数层。在这种约束条件下,敏感元
    件的刚度可行域并非一个连续的曲面,而是由许多离散的点组成,且这
    些点所对应的三向刚度特性满足(II)式的约束关系,即这些离散点都在
    (II)式表示的隐式曲面上。

    3.4分别计算每个敏感元件在每种铺层角与铺层厚度组合下的三向
    刚度特性值EI1、EI2、GJ,计算步骤如下:

    首先计算纤维主向刚度系数Qij,

    Q66=G12

    Q12=μ12Q22=υ21Q12,

    其中E11、E22、μ12、μ21、G12为与材料相关的工程弹性常数
    然后计算层合板拉伸刚度矩阵,

    A ij = Σ k = 1 N t k ( Q ij ) k ( i = 1,2,6 ; j = 1,2,6 ) ]]>

    其中N为层合板的层数,tk为单层厚度

    对于每一层复合材料铺层,铺层角为θ时的刚度系数计算公式,

    Q 11 = m 4 Q 11 + 2 m 2 n 2 ( Q 12 + 2 Q 66 ) + n 4 Q 22 ]]>

    Q 12 = m 2 n 2 ( Q 11 + Q 22 - 4 Q 66 ) + ( m 4 + n 4 ) Q 12 ]]>

    Q 22 = n 4 Q 11 + 2 m 2 n 2 ( Q 12 + 2 Q 66 ) + m 4 Q 22 ]]>

    Q 16 = m 3 n ( Q 11 - Q 12 ) + mn 3 ( Q 12 - Q 22 ) - 2 mn ( m 2 - n 2 ) Q 66 ]]>

    Q 26 = mn 3 ( Q 11 - Q 12 ) + m 3 n ( Q 12 - Q 22 ) + 2 mn ( m 2 - n 2 ) Q 66 ]]>

    Q 66 = m 2 n 2 ( Q 11 + Q 22 - 2 Q 12 - 2 Q 66 ) + ( m 4 + n 4 ) Q 66 ]]>

    其中m=cosθ,n=sinθ,θ为当前层的铺层角

    EI1、EI2、GJ计算式如下

    EI1=A11-A12*A12/A22;

    EI2=A22-A12*Al2/A11;

    GJ=A66;

    3.5绘制刚度点,分别创建每一个敏感元件的刚度可行域

    步骤3.4中,每一组EI1、EI2、GJ数据作为敏感元件的一个刚度点,
    将敏感元件的所有刚度点绘制到EI1、EI2、GJ构成的三维空间中,这些
    离散的刚度点即构成该敏感元件的刚度可行域,每个敏感元件的可行域
    所包络的空间均呈现类似三棱台的形状,在同一铺层厚度下,各刚度点
    在一帆状曲面上,从不同角度观察此帆状曲面如附图2、附图3、附图4;
    不同厚度的铺层对应的帆状曲面形状相似,彼此为比例缩放关系,缩放
    比为厚度比的四次方,如附图5。

    3.6对每一个敏感元件,分别在其刚度可行域内选取特征点

    在类似三棱台的刚度可行域的顶面、底面和内部等分剖面上选取特
    征点,选取位置为上述每面的三个顶点和三条边的中点,其中内部等分
    剖面的数量根据刚度可行域的大小选取,一般选取1~4个,特征点数目
    介于18~36个之间。如计算量可以承受,应考虑在上述各面上适量增选
    特征点,以提高跨音速颤振速度响应面的插值拟合精度。

    (4)对每个敏感元件,分别建立该敏感元件的跨音速颤振速度响应面

    一般情况下,翼面上会有多个敏感元件,对应每个敏感元件,有不同
    的刚度可行域和相应的特征点,因此需要分别针对每个敏感元件建立跨
    音速颤振速度响应面,其建立过程分为如下步骤:

    4.1创建刚度特征点对应的翼面有限元计算模型文件

    对于每个敏感元件,去除敏感元件上的原有铺层,用敏感元件的每
    一个特征点所代表的铺层铺覆在敏感元件上,翼面其他部位的铺层不变,
    由此建立敏感元件的每一个特征点对应的翼面有限元计算模型文件,每
    个敏感元件创建与特征点一一对应的18~36个有限元计算模型文件。如
    增选特征点,则对应增加相应的有限元计算模型文件;

    4.2分别计算4.1中每个翼面有限元计算模型的非定常气动力,得
    到相应的跨音速颤振临界速度Vcp,如增选特征点,则Vcp个数会对应增加;

    4.3对每个敏感元件,分别建立跨音速颤振速度响应面

    将敏感元件的每个特征点的轴向刚度特性EI1、剪切刚度特性GJ和
    厚度δ作为自变量(18~36组或更多),对应的翼面跨音速颤振速度Vcp
    (18~36个或更多)作为目标值,进行高维曲面拟合与插值,得到敏感
    元件的跨音速颤振速度响应面。由于另一方向的轴向刚度EI2与EI1和GJ
    存在(II)式所示的固定的约束关系,则EI2不需要作为自变量参与跨音速
    颤振速度响应面的建立。

    在响应面建立后,在敏感元件的刚度可行域内的任意一点,通过查
    询对应的响应面来代替复杂的非定常气动力计算,即可快速获得该状态
    下翼面的跨音速颤振速度Vcp。

    本发明技术方案的优点和积极效果是:

    1.利用三向刚度来描述敏感元件的铺层和质量特性,解决了如何建
    立覆盖复合材料板各种铺层角度与铺层比例可行域的技术难点;

    2.通过采用有限个特征点建立响应面的方法,保证了可行域里所有
    点进行跨音速颤振速度插值时,都落在所建立的响应面边界内,保证了
    较高的插值精度;

    3.利用查询响应面的方式代替复杂的非定常气动力计算用以获取翼
    面的跨音速颤振速度,大幅提高了计算速度,缩短了计算时间,解决了
    复合材料翼面跨音速颤振优化在工程上应用的难点;

    4.经过对某飞机复合材料翼面跨音速颤振优化计算的验证,建立的
    基于三向刚度的跨音速颤振速度响应面误差范围为±1m/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和最大许用厚度dmax

    3.2.1确定最小许用厚度dmin

    该敏感元件2的最小许用厚度dmin为8层(单层0.12mm);

    3.2.2确定最大许用厚度dmax

    该敏感元件2的最大许用厚度dmax为56层(单层0.12mm);

    3.3根据铺层角θ、敏感元件的最小许用厚度dmin和最大许用厚度dmax,
    确定所有可能的铺层组合;

    3.4分别计算敏感元件2在每种铺层组合下的三向刚度EI1、EI2、GJ,
    计算步骤如下:

    首先计算纤维主向刚度系数Qij,

    Q66=G12

    Q12=μ12Q22=υ21Q11,

    其中E11,E22、μ12、μ21、G12为与材料相关的工程弹性常数
    然后计算层合板拉伸刚度矩阵,

    A ij = Σ k = 1 N t k ( Q ij ) k ( i = 1,2,6 ; j = 1,2,6 ) ]]>

    其中N为层合板的层数,tk为单层厚度

    对于每一层复合材料铺层,铺层角为θ时的系数的计算公式,

    Q 11 = m 4 Q 11 + 2 m 2 n 2 ( Q 12 + 2 Q 66 ) + n 4 Q 22 ]]>

    Q 12 = m 2 n 2 ( Q 11 + Q 22 - 4 Q 66 ) + ( m 4 + n 4 ) Q 12 ]]>

    Q 22 = n 4 Q 11 + Zm 2 n 2 ( Q 12 + 2 Q 66 ) + m 4 Q 22 ]]>

    Q 16 = m 3 n ( Q 11 - Q 12 ) + mn 3 ( Q 12 - Q 22 ) - 2 mn ( m 2 - n 2 ) Q 66 ]]>

    Q 26 = mn 3 ( Q 11 - Q 12 ) + m 3 n ( Q 12 - Q 22 ) + 2 mn ( m 2 - n 2 ) Q 66 ]]>

    Q 66 = m 2 n 2 ( Q 11 + Q 22 - 2 Q 12 - 2 Q 66 ) + ( m 4 + n 4 ) Q 66 ]]>

    其中m=cosθ,n=sinθ,θ为当前层的铺层角
    EI1、EI2、GJ计算式如下

    EI1=A11-A12*A12/A22;

    EI2=A22-A12*A12/A11;

    GJ=A66;

    对于本实例来说,该尾翼选用的复合材料的工程弹性常数为:
    E11=1.29E+11Pa,E22=9.82E+9Pa,μ12=0.3056,G12=5.29E+9Pa,铺层
    的单层厚度tk=0.12mm,敏感元件2的最小许用dmin厚度为8层,最大
    许用厚度dmax为56层。

    3.5步骤3.4中,每一组EI1、EI2、GJ数据做为敏感元件2的一个
    刚度点,将敏感元件2的所有刚度点绘制到EI1、EI2、GJ构成的三维空
    间中,即构成敏感元件2的刚度可行域3,刚度可行域3所包络的空间呈
    现类似三棱台的形状,参见附图6。

    3.6在刚度可行域3内选取特征点

    在类似三棱台的刚度可行域3的顶面4(最小许用厚度面)、底面5
    (最大许用厚度面)和内部等分剖面6上选取特征点7,本实例中选择两
    个三等分剖面6,厚度分别为24层和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中每一个翼面有限元计算模型的非定常气动力,
    得到翼面跨音速颤振速度Vcp,特征点的轴向刚度EI1、剪切GJ,厚度δ
    及相应的翼面跨音速颤振速度Vcp如下表:

    表1.特征点的刚度与跨音速颤振速度


    EI1:轴向刚度(MPa),GJ:剪切刚度(Mpa),δ:厚度(mm),Vcp:跨音速颤振速度(m/s)

    4.3建立敏感元件2跨音速颤振速度的响应面8

    将表1中每一个刚度特征点7的轴向刚度EI1、剪切刚度GJ和厚度
    δ作为自变量,对应的翼面跨音速颤振速度Vcp作为目标值,输入至响应
    面软件DesParO中进行高维曲面拟合与插值,得到敏感元件2的跨音速
    颤振速度响应面8。DesParO软件是由德国科学院数学与科学计算研究所
    开发的,由德国scapos?AG公司代销的软件,该软件是用于建立优化设
    计响应面的应用软件,是一种用于多目标参数优化设计的工具,软件可
    以在高维空间建立响应面,并有很高的拟合精度。软件还能够进行响应
    面的误差分析,保证计算结果具有3西格玛级置信度。

    在响应面8建立后,在该敏感元件2的刚度可行域3内的任意一点,
    通过查询对应的响应面8来代替复杂耗时的非定常气动力计算,即可快
    速获取该状态下的翼面跨音速颤振速度Vcp。

    4.4敏感元件2跨音速颤振速度响应面8的误差分析

    以表1中35个特征点的数据为基础,每次选用34个点建立跨音速
    颤振速度响应面,另一点作为验证点。采用对比分析的方法,以验证点
    的刚度EI1、GJ和厚度δ作为输入的自变量,查询响应面得到跨音速颤振
    速度Vcp1,Vcp1与该点经过非定常气动力求解得到的真实的跨音速颤振速度
    Vcp相比,即可得到该点的插值误差。此外,响应面软件DesParO有插值
    误差估算功能,也可大致估算出该点的插值误差范围。分析中获得的数
    据如下表2所示,误差曲线见附图9。

    表2.跨音速颤振速度响应面误差分析



    Vcp:实际跨音速颤振速度,Vcp1:插值跨音速颤振速度,参考误差均为±值
    Vcp、Vcp1、参考误差、实际误差的单位均为m/s

    在厚度δ确定时,某个敏感元件1对应的跨音速颤振速度响应面剖
    面形状如附图10所示。由跨音速颤振速度响应面建立的过程及相关的数
    据与图表,可以得出如下结论:

    1.软件估算的插值误差为±0.5m/s,实际验证的误差为±1m/s,实
    际误差大部分在估算的插值误差范围内,响应面插值拟合的精度较高;

    2.所建立的基于三向刚度可行域的跨音速颤振速度响应面表面光
    滑,没有突出的尖点,响应面的鲁棒性好;

    3.该跨音速颤振速度响应面的建立花费数天时间,计算量在工程上
    可以接受;

    4.响应面建立后,在其上查询一次跨音速颤振速度的时间介于
    0.1s~0.3s之间,远低于完成一次非定常气动力计算所需时间,用查询
    响应面的方式来代替复杂耗时的非定常气动力计算,可以实现复材翼面
    跨音速颤振速度的快速获取。

    关于本文
    本文标题:一种快速获取复材翼面跨音速颤振速度的方法.pdf
    链接地址://www.4mum.com.cn/p-5817370.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 今日重庆百变王牌开状 浙江十一选五最大遗漏数据 金牛娱乐 股票 腾讯股票行情 新疆11选5开奖查询 西安福彩中心网址 申城麻将 湖北11选5中奖新闻 2015年上证指数最高点 广西十一选五开奖结果走势图百度 江西十一选五多乐彩 山东11选5任二计划 股票分析师资格证 赛车北京pk10是什么 股票指数期货在到期日以成交股票进行交割