运输类飞机货物空投时货舱地板载荷的确定方法

标题：运输类飞机货物空投时货舱地板载荷的确定方法

摘要：本发明属于飞机飞行载荷设计技术领域，涉及对运输类飞机货物空投时货舱地板载荷确定方法的改进。其特征在于，确定货舱地板载荷的步骤如下：飞机空投货物参数初始化；飞机有限元模态分析；飞机非定常气动力计算；飞机空投货物动响应计算；计算货舱地板结构载荷动响应；确定空投时货舱地板的设计载荷。本发明提供了一种改进的运输类飞机货物空投时货舱地板载荷确定方法，提高了货舱地板载荷的确定精度，保证了货舱地板的结构安全性和飞机的飞行安全。

申请号：CN201310564115.9

申请日：2013/11/13

申请人：中国航空工业集团公司西安飞机设计研究所

首项权利要求：1.运输类飞机货物空投时货舱地板载荷的确定方法，其特征在于，确定货舱地板载荷的步骤如下：
1.1、飞机空投货物参数初始化：依据国军标GJB67.2A-2008中的空投阵风准则进行飞机空投货物参数初始化，需要初始化的空投货物参数为：输入飞机的有限元模型FEM、飞机的质量ma、飞机飞行高度H、飞机飞行速度V、平均气动弦长c、货物的质量mb、货物在飞机货舱内的初始装载位置x0、货物的转动惯量Iby、货物在货舱内运动的加速度ab、模态数目k、离散阵风强度uw、离散阵风尺度Lw；
1.2、飞机有限元模态分析 : 以上述飞机的有限元模型FEM、模态数目k为控制参数，采用Nastran软件对飞机的有限元模型FEM进行模态分析，得到飞机的模态矩阵Φ、每阶模态的振动圆频率、广义质量矩阵Mqq和广义刚度矩阵Kqq，模态矩阵Φ包含飞机的刚体运动模态和弹性振动模态；
1.3、飞机非定常气动力系数矩阵计算：以上述模态矩阵Φ、飞机飞行高度H、飞机飞行速度V、平均气动弦长c、离散阵风强度uw和离散阵风尺度Lw为输入变量进行时域非定常气动力系数矩阵计算，采用亚音速偶极子格网法进行计算得到结构振动引起的非定常气动力系数矩阵Qa，采用阵风混合建模方法进行计算得到阵风引起的激励力系数矩阵Qw；
1.4、飞机空投货物动响应计算：
1.4.1、重构弹性飞机货物空投动响应微分矩阵方程 : 在步骤1.1的12个参数、步骤1.2的模态矩阵Φ、每阶模态的振动圆频率、广义质量矩阵Mqq和广义刚度矩阵Kqq和步骤1.3的结构振动引起的非定常气动力系数矩阵Qa及阵风引起的激励力系数矩阵Qw的基础上，重构弹性飞机货物空投动响应微分矩阵方程：
其中，ξ为对应于模态矩阵Φ的广义坐标位移，包含对应于飞机的刚体运动模态的广义坐标位移和弹性振动模态的广义坐标位移；
为对应于模态矩阵Φ的广义坐标速度，包含对应于飞机的刚体运动模态的广义坐标速度和弹性振动模态的广义坐标速度；
为对应于模态矩阵Φ的广义坐标加速度，包含对应于飞机的刚体运动模态的广义坐标加速度和弹性振动模态的广义坐标加速度；Φbz为货物所处位置处的飞机模态矩阵分量；Vb2a为货物在货舱内相对飞机的运动速度，通过对货物在货舱内运动的加速度ab的一次积分可以得到；x为货物在货舱内的移动距离，通过对货物在货舱内运动的加速度ab的两次积分可以得到；Qaξ为飞机的结构振动引起的非定常气动力；Qwuw为阵风引起的激励力；g为重力加速度；
为一个长度为k的列向量；
为利用Runge-Kutta法求解微分矩阵方程，将式(1)整理成一阶微分矩阵方程：
其中，Aae为弹性飞机本体与非定常气动力组成系统的状态矩阵；Baw为离散阵风对弹性飞机动响应的扰动矩阵；Bag为由于货物移动对飞机产生作用力的扰动矩阵；uw为离散阵风强度；xae为状态矢量：
其中，xa为非定常气动力的滞后根项目，模拟从飞机翼面上脱落的涡;
1.4.2、求解飞机1g平飞时的广义坐标位移ξTrim : 飞机1g平飞时满足以下条件：
1.4.2.1、广义坐标加速度
和广义坐标速度
均为零向量；
1.4.2.2、货物相对飞机运动速度和位移均为0；
1.4.2.3、无离散阵风激励力；
根据上述条件，飞机1g平飞时的广义坐标位移ξTrim的计算如公式(4)所示：
其中，Φbz0为货物在货舱内初始装载位置处的飞机的模态矩阵分量；
1.4.3、求解飞机货物空投动响应微分矩阵方程 : 以飞机1g平飞时的广义坐标位移ξTrim为公式(2)的初始求解条件，采用变步长Runge-Kutta法求解飞机货物空投动响应微分矩阵方程，得到飞机的广义坐标位移ξ的动响应、广义坐标速度
的动响应和广义坐标加速度
的动响应 :
其中，Cae为弹性飞机本体与非定常气动力组成系统的输出矩阵；Daw为离散阵风对弹性飞机动响应输出的传递矩阵；Dag为由于货物移动对飞机动响应输出的传递矩阵；
1.5、计算货舱地板载荷动响应 : 根据步骤1.4的飞机的广义坐标位移ξ的动响应、广义坐标速度
的动响应和广义坐标加速度
的动响应，计算货物对飞机货舱地板的作用载荷
动响应，如式(6)所示：
其中，
为飞机质心处的加速度；
为飞机的俯仰角加速度；
为飞机的俯仰角速度；
为货物距离飞机质心的距离；
1.6、确定空投时货舱地板的设计载荷：根据步骤1.5中得到的货舱地板载荷动响应
选择
的极大值和极小值作为货舱地板的空投设计载荷。

专利类型：发明申请