细长体分布式载荷的缩聚及其在气动伺服弹性地面模拟中的应用.pdf

引用格式:吴志刚 尹烈鹏 余长坤.细长体分布式载荷的缩聚及其在气动伺服弹性地面模拟中的应用.航空兵器 ():.():.()细长体分布式载荷的缩聚及其在气动伺服弹性地面模拟中的应用吴志刚 尹烈鹏 余长坤(北京航空航天大学 航空科学与工程学院 北京)摘 要:气动伺服弹性地面模拟试验是近年发展出的一种新型气动弹性地面试验技术 其中的一个关键问题就是非定常气动力的缩聚 即将飞行器所受的分布式气动载荷缩聚为少数几个集中载荷传统缩聚方法以简单的力、力矩平衡准则对分布力进行处理 或者是通过求解静不定问题 把支反力作为缩聚力 均存在一些缺点 本文提出一种基于一维样条插值的细长体分布式载荷缩聚方法 采用关键模态相似准则来优化缩聚点位置 以达到细长体在缩聚点集中式载荷作用下的动响应与分布式载荷作用下的动响应最为接近的效果 利用简支梁、悬臂梁两个算例验证了所提出的分布式载荷缩聚方法的精度 并将该方法应用于一个细长体导弹的气动伺服弹性地面模拟中 数值仿真表明 所提方法能够快速准确将导弹的非定常气动力缩聚成实时的集中力 满足气动伺服弹性地面模拟试验的精度和快速性要求关键词:气动伺服弹性 地面试验 载荷缩聚 细长体 样条插值 中图分类号:.文献标识码:文章编号:():./.引 言气动弹性地面模拟试验是近年来发展的新型气动弹性地面试验技术 该方法以真实飞行器为研究对象 根据传感器采集到的结构振动信息 实时计算缩聚点的集中非定常气动力 使用力加载装置进行集中力的实时加载 从而在地面无风洞情况下模拟飞行器在空中飞行时的气动弹性特性 世纪 年代 美国的 最早提出了颤振地面仿真的概念 初步探索了模拟导弹飞行时机翼气动力的地面试验方法 国内早期研究代表为潘树祥和齐丕骞 主要进行了结构热颤振的地面试验研究 受限于当时的理论方法及硬件水平 试验误差较大近年来 诸多学者开展了进一步的研究 美国 公司的 等 于 年开展了颤振地面模拟试验并命名为干风洞技术()试验使用少数激振器对结构施加集中力来模拟分布式气动载荷并进行了颤振预测 年 许云涛等 对颤振地面模拟试验中的非定常气动力模拟进行了细致的研究 随后 和张仁嘉等 开展了舵系统颤振地面模拟的实验研究 设计了只具有沉浮和旋转自由度的舵系统实验模型 年 等开展了矩形平板机翼颤振地面模拟试验 其他国内学者如胡巍和 等 在气动力缩聚理论、力控制方法、试验误差干扰影响等方面都对颤振地面模拟试验方法做出了贡献 对于气动伺服弹性地面模拟试验 和楚龙飞等 针对带有控制回路的细长体导弹进行了试验研究 受力加载设备特性和气动力计算理论限制 未考虑分布式非定常气动力的缩聚 仅将弹体划分为前后两个气动段 一般认为细长体气动段划分数量越多 气动力计算越准确缩聚()在力学上是指为满足某种等效条件 将分布力转化成作用在若干结点上集中力的处理方法 在气动伺服弹性地面模拟试验中 研究对象可能涉及长细比大于 的细长体 如导弹、火箭 细长体所受分布式载荷一般包括分布式气动载荷、惯性载荷等收稿日期:基金项目:航空科学基金项目()作者简介:吴志刚()男 江西萍乡人 教授 年第 卷第 期 年 月 航空兵器 .其通常是非定常的 且可以简化为沿细长体中心轴呈一维分布的分布力 该试验需要通过一定的加载设备来模拟这种分布式载荷 由于加载设备数量及加载空间的限制 需要将细长体分布式载荷等效缩聚为少数个作用点的集中载荷 因此 分布式载荷缩聚方法的优劣决定了试验结果的准确性工程上常用的一种缩聚方法是静力学等效方法 即将试验结构视为刚体 依据刚体静力学平衡的原则进行分布式载荷的等效缩聚 该方法在飞行器静力试验中应用广泛 但对动力学试验和气动弹性地面模拟试验并不适用 在这类试验中 分布式非定常气动力与试验结构的弹性变形相关 传统方法无法保证试验结构在缩聚后集中力作用下的动响应与分布式载荷作用下动响应等效 目前已有的文献 主要针对平板机翼等开展气动力缩聚研究 文献研究的也是一种小展弦比翼面的非定常气动力缩聚方法 而对于细长体的非定常气动力缩聚方法还鲜见有文献报道 大多数气动弹性地面模拟试验研究关注的是颤振问题 关于导弹气动伺服弹性的地面模拟研究也较少本文旨在提出一种细长体分布力的等效缩聚方法既可替代传统方法应用于细长体静力试验 也可应用于细长体动力学试验 与传统方法相比 该等效缩聚方法只与外载荷分布情况以及缩聚点的位置有关 不需要研究对象的结构信息 为验证方法的精度 通过两个算例分别验证缩聚方法在静力学、动力学中的效果 最后将本文方法应用于细长体导弹的气动伺服弹性地面模拟中 通过数值仿真来评估应用效果 载荷缩聚的问题描述细长体分布载荷的缩聚问题为:如图 所示的一维梁式结构中 作用于梁上有集度为()的分布力 在梁上寻找 个缩聚点(横坐标记为 )各缩聚点的集中作用力为 使得集中载荷的作用与原分布载荷等效 求缩聚点的位置及其集中力图 分布载荷缩聚问题示意图.这里需要说明一下“等效”的含义“等效”在不同的情况下含义是不同的 在传统方法中 静力学等效认为结构是一个刚体 从力的平衡角度来处理 即若缩聚前后 作用在结构上的合力、合力矩均相同 则视为等效若满足合力、合力矩相同的同时 要求缩聚点处位移、转角连续 则是另一种标准下的等效 本文方法的等效是基于动力学响应特性 传统方法的回顾传统方法主要有两种:静力学等效、静力学求解 下面分别简述其基本原理和处理方法.静力学等效静力学等效是假设梁为刚体 缩聚点的集中力与原分布作用力是静力等效的(在二维平面内即合力相等对参考点的合力矩相等)在实际处理中 首先将梁分成 段 分别对各段上的分布力进行集中 再按静力等效分配到各段梁的 个结点上 如图 所示图 静力学等效示意图.图中 与()的关系满足:()()()对每一段进行相同的求解 得到各缩聚点的集中力大小这种方法的问题是不能保证结构变形等效 需要有较多的缩聚点 才可满足工程中的等效加载要求.静力学求解静力学求解是在梁上选取 个缩聚结点 将这些缩聚点的平动位移自由度约束住 求解一个静不定问题如图 所示 求解得出的各缩聚点的约束反力即为缩聚后的集中力图 静力学求解方法示意图.相比.节方法 该方法是在合力等效、合力矩等效两个方程的基础上 增加了结构变形协调关系 该方法获得的缩聚力对原分布载荷的模拟精度较高 但如果缩聚点的位置不合适 可能导致求出的缩聚集中力出现不合理的值 此外 该方法与结构的刚度分布有关 事先需要知道结构的刚度分布信息 再通过有限元法完成静不定结构的静力求解 对于非定常分布载荷的实时缩聚并不合适航空兵器 年第 卷第 期基于样条插值的分布力缩聚方法.分布载荷的离散化为了数值处理的方便 不失一般性 首先应用第.节所述方法 将细长体上的连续分布式载荷()离散化到 个密集结点上 设这些密集结点的坐标形成的向量为 作用在其上的作用力形成的向量记为 即有()()当 较大(即结点密集)时 可以认为这种离散化所带来的误差是可以忽略的.分布载荷的缩聚在细长体上设置 个()缩聚点 其坐标形成的向量记为 将细长体(梁)在载荷 作用下密集点和缩聚点的位移分别记为 和 设二者有如下关系:()式中:为维数为 的转换矩阵 其元素与密集点和缩聚点的坐标 和 有关设经过载荷缩聚后 作用于缩聚点上的集中力组成的向量为 为了使得 与 等效 采用以下等效原则:这两种载荷在对应结点虚位移上做的虚功相等 即 ()将式()代入式()中 有 ()比较式()的两边 得到 ()因此 载荷缩聚问题转化为如何建立起密集点与缩聚点之间的位移转换关系 一旦获得 则可利用式()求得等效缩聚后的载荷.一维梁式结构的样条插值建立密集点与缩聚点之间位移转换关系的问题实际上是一个一维插值问题 细长体结构常由沿中轴线的梁作为有限元进行建模 因而常用的一维插值方法是梁样条插值法()但是梁样条插值法需要获得已知点的位移及斜率信息 这里采用一种更加简便的方法 即采用曲面样条插值法的一维情况 该方法假定拟合插值函数为()()()式中:()为引入的修正因子 一般与 的大小相当 这里取 和 为待定的系数 对于给定的 个结点 则有 ()()式中:()在该公式中 有个方程 个待定的未知数 为了确定未知数 还需补充两个方程 这里考虑力和力矩的平衡条件 补充方程为 ()()为此 构成的方程写成矩阵形式为 ()式中:()()式()可以简写为矩阵和向量的形式:()由式()求得向量 和 后 再利用式()可以求得密集点的位移及斜率:()()式中:()()()()()()式()()记为 ()()将式()代入式()()中 可以得到密集点插值位移 和插值斜率 与缩聚点位移 之间的关系:()()式中:为矩阵 删除前 列得到的矩阵 为矩阵吴志刚 等:细长体分布式载荷的缩聚及其在气动伺服弹性地面模拟中的应用删除前 列得到的矩阵.缩聚点的位置优化为了评价缩聚的效果 这里给出缩聚等效的评价准则 即根据缩聚点位移 经过插值得到的密集点位移和斜率 使得如下定义的误差 最小 即 ()()()式中:表示向量范数 该误差函数同时考虑了变形及转角 并分别做了归一化处理在实际处理中 为了避免事先求解结构变形 和 可以利用试验测得的结构固有模态振型或假设模态振型()来代替 即 ()()()式中:()()()()()()若要考虑前 阶模态 则 可以写为 ()()()这里忽略各阶模态的影响权重的差异 即均取为 以上误差实际上与缩聚点坐标 有关 因此缩聚点位置优化可以转化为设计变量为、目标函数为的优化问题 对于以上的优化问题 可以采用非梯度的智能优化算法来寻优 例如遗传算法 作为一种通用的问题求解方法 遗传算法采用简单的编码技术来表示各种复杂的结构 并通过对一组编码表示进行简单的遗传操作和优胜劣汰的自然选择来指导学习和确定搜索的方向载荷缩聚方法的验证本节包含算例(均匀简支梁模型)和算例(变截面悬臂梁模型)梁的模型长 沿轴向被划分为 个杆单元 包含 个节点 高度方向为 向 如图 所示 图中 均匀梁的截面为宽.、高.的矩形变截面梁的根部截面宽.、高.自由端截面宽.、高.截面尺寸线性变化 结构的阻尼比为.材料属性设置如表 所示.算例 及结果分析算例:一根 长的两端简支均匀梁 受到集度为()/的均匀分布载荷作用 缩聚点数量为 均匀分布在梁上 分别采用静力等效法、静力求解法和本文提出的缩聚方法得到缩聚载荷 并将各种缩聚载荷作用在梁上产生的位移变形与原均布载荷产生的变形做比较图 梁及其所处坐标系.表 算例使用的材料属性 算例泊松比杨氏模量/密度/(/).结果分析:()原均布载荷作用将分布载荷()向密集点离散化 当密集点足够多时 可以认为离散前后的效果是相同的 密集点选择梁上的 个有限元节点 离散结果为第 和 个点的集中力为.其余点的集中力为.在等效密集点载荷作用下 通过有限元软件 得到简支梁的静变形 并以此为参照()静力等效法由.节介绍的静力等效法 计算出作用在缩聚点上的集中力 将此集中力作为载荷 得出简支梁变形 经计算与密集点载荷的误差为 .()静力求解法采用静力求解法 在缩聚点的位置添加垂直于梁的约束 将静定问题转化为静不定问题 使用 求解此静不定问题 得到个约束力 缩聚力即取约束力的负方向 将得到的缩聚力加在梁上 进行静力分析 经计算与密集点载荷的误差为 .()等效缩聚法采用.节所述的静力学等效法 将分布载荷()向密集点 离散化 为分析方便 密集点即取梁上的有限元节点 即横坐标从 至、间隔.的均匀分布的 个点根据密集点与缩聚点的坐标 得到从密集点载荷 插值得到缩聚点载荷 所需的插值矩阵 根据 得出缩聚力 加载后 经计算与密集点载荷的误差为 .将 种缩聚方法与原密集点载荷产生的静变形分别进行对比 结果如图 所示 误差大小如表 所示 在本算例的静力学问题下 静力等效法的误差相对较大 而静力求解法和等效缩聚法的误差均较小在本算例中可以发现 静力求解法的精度略高于本文提出的等效缩聚法 但静力求解法的使用需要事先建航空兵器 年第 卷第 期立有限元模型 而对于较为复杂的结构并不方便 除此之外 静力求解法受限于有限元软件的求解速度 在针对多数场景下的非定常动载荷时并不实用 而等效缩聚法在动载荷的实时计算上体现出了更大的优势图 三种缩聚法集中载荷下简支梁的变形.表 采用不同方法得到的缩聚载荷及误差 缩聚方法缩聚点载荷静变形误差/静力等效法.静力求解法.等效缩聚法.算例 及结果分析算例:一根 长的悬臂变截面梁 受到集度为()/的均匀分布载荷作用 初始缩聚点数量为 均匀分布 首先采用本文提出的缩聚方法进行载荷缩聚 比较缩聚载荷与原分布载荷产生的位移和转角变形更进一步 采用悬臂变截面梁的前两阶固有模态振型作为假设模态振型 进行缩聚点的位置优化 并与之前的结果比较 最后将载荷设置成谐振力 比较载荷缩聚前后结构的时域响应和频率响应曲线结果分析:对于初始均匀分布的个缩聚点 通过本文方法 得到缩聚力 .将得到的缩聚力加在梁上 进行静力分析校验 梁位移变形、转角变形的结果如图 所示 二者的相对误差 .误差控制很好现对梁在 平面内的前两阶模态(一阶弯曲和二阶弯曲)进行插值 初始缩聚点数量为 位置均匀分布 采用遗传算法进行缩聚点位置优化 优化变量为缩聚点位置 即一个包含个坐标的位置向量 约束条件为每个坐标仅能取 以内的正整数 优化目标为使式()最小化本算例使用的遗传算法采用二进制编码 初始种群数量设置为 变异概率.迭代代 优化前后的结果如表 所示 模态的插值效果如图 所示 图中曲线上的标记代表插值点的位置图 等效缩聚法载荷下悬臂梁的位移转角变形.表 优化前后结果对比 缩聚点位置模态插值误差/静变形误差/优化前.优化后.图 一阶弯曲模态及插值结果.图 二阶弯曲模态及插值结果.可以发现 以模态插值误差为目标函数的遗传算法吴志刚 等:细长体分布式载荷的缩聚及其在气动伺服弹性地面模拟中的应用优化降低了模态插值误差 但用此结果重新进行悬臂变截面矩形梁的静力学实验时 静变形误差由原来的.增大到.原因是作为目标函数 模态插值误差反映的是结构在非定常气动力下的动响应精度 而静变形误差反映的是结构在静力作用下的变形 二者不能混为一谈 然而 尽管静变形误差在优化之后有所增加 但误差仍然在可接受范围之内接下来进行结构在动载荷下的响应比较 模拟矩形机翼上的典型椭圆气动力分布 在悬臂梁上施加一系列呈椭圆分布的谐振气动集中力 如图 所示图 模拟机翼表面气动力的椭圆分布.气动力表达式为().()式中:为离散的气动点坐标 所有气动点的气动力载荷频率相同且没有相位差本算例中的悬臂变截面矩形梁在 平面内的前三阶模态频率分别为.和.载荷的频率分别取为 和 缩聚点位置均匀分布 则由 得出缩聚点载荷为 .()分别对密集气动点载荷(式()和缩聚点载荷(式()进行瞬态响应分析 取梁自由端端点作位移和转角的响应曲线 如图 所示 激振力频率及对应误差(按式()计算)如表 所示 缩聚动载荷下的结构响应与原椭圆分布载荷下的结构响应几乎完全吻合对结构进行频率响应分析 取梁自由端端点 得到端点位移随载荷频率变化的频响曲线 载荷分别采用密集力载荷和缩聚力载荷 即式()()结果如图所示 结果表明在三阶模态频率内 由个缩聚点插值得到的缩聚力能很好地模拟悬臂梁上分布力的作用效果图 梁自由端时域响应曲线对比.表 激振力频率及对应误差 参数数值激振力频率/误差/.由此可见 本文提出的缩聚法能精确地模拟分布载荷下结构的动响应 同时需要指出的是 如果所关心的模态包含了高阶模态 即振型较为复杂 则将导致从结构点(缩聚点)到气动点(密集点)的振型插值 即式()的误差增加 因此需要适当增加缩聚点的个数载荷缩聚方法的工程应用.气动伺服弹性地面模拟试验系统建模针对细长体导弹气动伺服弹性问题 使用所提出的分布载荷缩聚方法 可将弹体和舵面上的密集非定常气动力等效到载荷缩聚点上的集中非定常气动力 从而采用少量力加载设备(如激振器)进行非定常气动力的实时加载 实现气动伺服弹性系统的地面模拟 建立的气动伺服弹性地面模拟试验仿真系统 可以分为缩聚非定常气动力计算系统、弹性导弹系统、飞行控制系统及舵机系统图 悬臂梁自由端 向位移的幅频响应曲线.航空兵器 年第 卷第 期.缩聚非定常气动力计算系统采用气动导数法计算细长体导弹的气动力 将弹体划分为若干密集气动段 舵面视为一个气动段 第 个气动段上作用的密集非定常气动力记为 通过气动导数法可将其表示为()()()式中:为第 段的升力线斜率 为气动力计算参考面积 式()写成矩阵为 ()式中:()()对于这一密集非定常气动力的缩聚 可选取个测量缩聚点和 个载荷缩聚点 其中测量缩聚点需要考虑位移和斜率插值 而载荷缩聚点只考虑位移插值 因此一般有 密集点的位移和斜率可由测量缩聚点的位移插值得到 即 ()根据虚功原理 载荷缩聚点的集中非定常气动力可以表示为 ()将式()()带入式()可得 ().弹性导弹系统弹性导弹上作用的载荷包括载荷缩聚点的集中非定常气动和舵面偏转惯性力 采用广义坐标的气动弹性方程可以写为 ()式中:为广义质量、广义阻尼和广义刚度矩阵为载荷缩聚点的模态振型矩阵 为舵面耦合质量矩阵 为舵面偏转角为了计算弹性导弹上作用的缩聚非定常气动力 弹性导弹系统需要输出物理坐标下的测量缩聚点位移和速度 采用广义坐标表示为 ()式中:为测量缩聚点的结构振动模态振型矩阵 为舵面偏转 构造的测量缩聚点的模态振型矩阵不失一般性 飞行控制系统以俯仰角速度()反馈为例 状态变量取 输入变量取 输出变量取 则气动弹性方程写成状态空间表达式为 ().飞行控制系统和舵机系统引入飞行控制系统的传递函数 则控制律解算指令为()舵机作为气动伺服弹性系统的执行机构 其传递函数一般可用三阶有理式表示:()写成状态空间表达式为 ()式中:状态变量 输入变量为控制解算指令 输出变量 由此可以建立细长体导弹的气动伺服弹性地面模拟试验系统 包括式()的缩聚非定常气动力计算系统、式()的弹性导弹系统、式()的飞行控制系统和式()的舵机系统.仿真模型仿真模型是一个具有一定结构刚度分布的细长体导弹 其示意图如图 所示 弹体划分为 个气动段 舵面视为一个气动段 气动伺服弹性分析时考虑沉浮和俯仰刚体模态、以及一阶弯曲(.)和二阶弯曲(.)弹性模态 弹性模态阻尼比取.在气动伺服弹性分析的关心频率范围内 舵面视为刚体 各气动段的气动导数采用商用 求解器计算得到 马赫数为.海平面高度图 仿真模型示意图.以气动网格节点的模态振型和振型斜率作为目标振型和振型斜率 根据两阶弹性模态插值误差最小进行缩吴志刚 等:细长体分布式载荷的缩聚及其在气动伺服弹性地面模拟中的应用聚点的位置优化 由于这一仿真模型的气动伺服弹性特性主要受刚体模态及一阶弯曲模态影响 而二阶弯曲模态的影响极小 因此两阶弹性模态权重系数分别取为.和.采用遗传算法优化得到 个载荷缩聚点和 个测量缩聚点(个位于弹体 个位于舵面)缩聚点位置均已标注在图 中 采用缩聚点模态振型插值得到的气动网格节点模态振型和振型斜率如图 所示 图中也绘制了缩聚点的模态振型(振型斜率)进一步分析插值误差 其中载荷缩聚点模态插值误差为.测量缩聚点模态插值误差为.具有相当高的插值精度 需要说明的是 位于弹体尾部的载荷缩聚点既是弹体非定常气动力的作用点 也是舵面非定常气动力的作用点弹体缩聚点采用本文所提梁的样条插值方法 而舵面缩聚点则采用简单的刚体插值方法图 气动网格节点模态振型及振型斜率插值结果(一阶弯曲模态和二阶弯曲模态).()舵机模型采用如下形式:.()飞行控制系统采用俯仰角速率反馈 控制系统传递函数为 .()将弹性导弹系统、缩聚非定常气动力计算系统、舵机系统和飞行控制系统进行组装 得到气动伺服弹性地面模拟时域仿真系统 如图 所示 给定闭环系统一个初始扰动 弹性导弹在扰动作用下产生变形 变形信号经过缩聚非定常气动力计算系统可得到载荷缩聚点的非定常气动力 而陀螺仪感知到的俯仰角速度信号传输给飞行控制系统生成控制指令信号 驱动舵机进行偏转非定常气动力和舵面偏转再作用于弹性导弹系统 如此往复形成闭环系统 随着风速的增大 系统在扰动作用下的振动幅值由收敛变为发散.仿真结果采用建立的气动伺服弹性地面模拟时域仿真系统进行时域仿真 仿真步长为.给弹性导弹系统一个微小的初始扰动 选取三个不同飞行速度的俯仰角速度响应信号 如图 所示 当飞行速度为./时导弹处于收敛状态 当飞行速度为./时 导弹处于发散状态 当飞行速度为./时 导弹处于等幅振荡的临界稳定状态 即为该细长体导弹在 .时的气动伺服弹性临界稳定速度 相应的失稳频率为.采用密集点非定常气动力进行理论气动伺服弹性稳定性分析 对比结果如表 所示 与理论气动伺服弹性临界速度相比 缩聚结果的临界速度误差仅为.仿真结果表明 采用所提分布载荷缩聚方法 结合遗传算法进行缩聚点的位置优化 将密集点的非定常气动力缩聚到有限数量缩聚点的非定常气动力 可以获得极高的缩聚精度 满足细长体导弹气动伺服弹性地面模拟试验需要图 气动伺服弹性地面模拟时域仿真系统.航空兵器 年第 卷第 期图 不同飞行速度下俯仰角速度时域响应.表 气动伺服弹性临界状态对比 临界速度/(/)频率/缩聚结果.密集点结果.结 论()本文提出了一种基于一维样条插值的细长体分布载荷缩聚方法 该方法精度高、速度快 且无需获知结构信息 适用于动力学试验和气动弹性地面试验()通过两个算例证明了本文缩聚方法在静、动载荷缩聚方面的准确性 配合缩聚点的位置优化 气动力模拟的误差能够进一步降低()将缩聚方法应用于细长体导弹的气动伺服弹性地面模拟中 数值仿真表明缩聚后系统的临界稳定速度和失稳频率与分布力情况一致 能够满足气动伺服弹性地面模拟试验的要求参考文献:.潘树祥 齐丕骞.地面模拟热颤振试验研究.强度与环境():.():.()./.许云涛 吴志刚 杨超.地面颤振模拟试验中的非定常气动力模拟.航空学报 ():.():.()许云涛.地面颤振模拟试验方法研究.北京:北京航空航天大学.:.().():.张仁嘉.飞行器气动伺服弹性若干关键问题研究.北京:北京航空航天大学.:.().():.胡巍 杨智春 谷迎松.带操纵面机翼气动弹性地面试验仿真系统中的气动力降阶方法.西北工业大学学报 ():.():.()./.():.():.楚龙飞.高超声速飞行器气动伺服弹性问题研究.北京:北京航空航天大学.:.().():.():.吴志刚 等:细长体分布式载荷的缩聚及其在气动伺服弹性地面模拟中的应用.():.宋巧治 李晓东.平板翼颤振地面模拟试验机理研究.结构强度研究():.():.().:.:.():.周明 孙树栋.遗传算法原理及应用.北京:国防工业出版社.:.:.()():.:航空兵器 年第 卷第 期