柔性扑翼全机气动结构耦合特性及稳定性研究

由于微型扑翼飞行器具有体积小、重量轻等特点,在飞行过程中会产生明显的柔性变形,因此在对扑翼气动特性进行数值模拟时,有必要考虑柔性变形的影响。基于结构动力学理论,发展了一种适用于扑翼全机气动结构耦合特性的数值模拟方法与一种适用于扑翼的结构动力学求解方法。方法利用哈密顿原理,对动能和应变能进行变分,进而得到扑翼动力学方程,采用结构有限元方法对运动方程进行离散并求解;采用基于嵌套网格的CFD求解器对微型扑翼全机非定常绕流进行数值模拟。采用CFD/CSD耦合求解器对微型扑翼飞行器全机气动结构耦合特性进行数值模拟,分析了结构变形对柔性扑翼气动特性的影响,并分别对扑动过程中刚性扑翼和柔性扑翼的压心变化范围进行研究,分析了结构变形对柔性扑翼稳定性的影响。     

微型扑翼低雷诺数绕流气动特性研究

针对微型扑翼所处典型飞行状态,数值模拟研究了微型扑翼绕流的低雷诺数气动特性.基于结构化嵌套网格求解了预处理后的三维非定常雷诺平均Navier-Stokes方程,空间离散采用了中心格式有限体积法,非定常时间推进为双时间法.通过与文献低雷诺数扑翼的试验值对比验证了本文算法的有效性,接下来通过大量计算研究了关键气动参数和展向折叠扑动模式的影响规律.研究结果表明迎角、扭转角和折叠扑动对升力影响较大,扭转角和减缩频率对推力影响较大.本文的研究得到了微型扑翼气动力特性的机理性结论,有助于理解微型扑翼飞行原理,为设计微型扑翼提供了参考依据.     

飞行器气动/结构多学科延迟耦合伴随系统数值研究

基于自主研发的大规模并行结构化网格雷诺平均Navier-Stokes（RANS）求解器PMB3D以及流固耦合代码FSC3D建立了飞行器静气动弹性数值模拟技术,进一步基于并行化伴随方程求解器PADJ3D,开展了Navier-Stokes方程与结构动力学方程耦合离散伴随的推导、构造。对各个学科伴随变量进行延迟处理,进行耦合伴随系统的解耦,学科之间的影响通过右端强迫项的形式在方程中体现,通过松耦合形式进行各个学科方程右端项数据传递,各个学科的伴随方程一定程度上能够相对独立,进一步实现了基于LDLT方法的结构伴随方程的高效求解;对典型客机柔性机翼进行梯度信息求解,并与考虑气动弹性影响的差分结果进行对比分析。数值模拟结果表明,延迟耦合伴随形式更有利于保持原有程序结构形式以及程序模块化,梯度计算精度完全满足优化设计需要,能够为柔性机翼飞行器气动/结构多学科优化设计提供研究基础与技术平台。     

基于双向流固耦合的柔性扑翼变形数值研究

通过研究柔性扑翼变形、惯性力及气动力的相互关系和作用机理,为扑翼飞行器的设计提供参考。在考虑扑翼运动特点的情况下,发展了一种求解雷诺平均N-S方程和结构动力学方程的双向流固耦合数值方法。对柔性翼惯性力方程进行改进,使其更方便与气动升力进行数值对比。针对柔性翼变形研究的困难性,提出了一种通过全翼面应变去分析柔性翼变形的方法。研究结果表明:柔性翼的扑动频率对惯性力和气动升力的影响程度较为一致;柔性翼应变主要集中在前缘和翼根周围,随着远离前缘和翼根,其应变值逐渐减小,应变一定程度上可以反映柔性翼的弦向变形;柔性翼变形导致推力在整个周期上的主要方向为正,能够显著提升扑翼飞行能力。     

考虑弹性转轴的扑翼三维气动数值计算

通过求解雷诺平均Navier-Stokes方程,研究了同时具有扑动和俯仰运动的三维扑翼气动特性.为了模拟扑翼的弹性特性,在扑翼俯仰轴处设计了扭簧.研究了扑翼平面形状、俯仰弹性轴与质心的相对位置、俯仰刚度变化等参数的影响.结果表明:俯仰角相对扑动角有一个大的相位超前量,且刚度系数越大该超前量越大;考虑了俯仰弹性特征后,计...     

微型飞行器低雷诺数矩形扑翼非定常气动特性的数值模拟

采用双时间步方法求解三维可压缩非定常N-S方程,数值模拟了微型飞行器低雷诺数矩形扑翼的非定常绕流,首先将得到的结果与文献进行了对比,数据间具有较好的一致性.然后针对不同的展弦比、减缩频率及初始攻角,计算了矩形扑翼的非定常气动特性及表面流态和动态压力分布,并分析了翼尖涡对扑翼非定常气动特性的影响.     

翼伞平面形状对翼伞气动性能的影响

对带气室的展弦比为3的不同平面形状翼伞模型的流场进行了三维定常数值模拟,详细考察了平面形状对翼伞气动性能的影响.运用有限体积法对三维坐标系下不可压雷诺时均Navier-Stokes(RANS)方程进行了直接求解,采用剪切应力输运(SST)k-ω二方程湍流模型进行湍流模拟.数值模拟得出的原始翼伞的气动性能参数与试验数据在...     

涡轮机械叶片的流固耦合数值计算方法

发展了一种流固弱耦合数值方法来判断涡轮机械叶片气弹稳定性.通过一种数据交换方法将计算结构动力学(CSD)的节点振动位移施加到计算流体力学(CFD)的叶片表面网格点上,CFD分析采用多层动网格技术实时更新可动域的网格点坐标,并通过有限体积法求解了用k-ε湍流模型封闭的雷诺平均Navier-Stokes(RANS)方程.以NASA 67转子叶片为例,在设计状态下,对比了叶片前4阶模态下一个振动周期内的气动功及模态气动阻尼比,初步分析了叶片的气弹稳定性.计算结果表明:所发展的流固弱耦合数值方法用于判断涡轮机械叶片气弹稳定性是可行的.      

跨声速风扇叶片变形对气动性能的影响

为了研究跨声速风扇叶片变形对气动性能的影响,发展了流固耦合方法计算风扇叶片在气动力作用下变形的数值程序,耦合求解了三维可压缩雷诺平均的Navier-Stokes (N-S)方程和描述叶片振动的结构动力学方程,并模拟了跨声速风扇Rotor 67转子叶片在气动载荷下的变形,重点分析了由于叶片变形带来的气动性能的变化.      

DLR-F6翼身组合体的高阶精度数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS)和剪切应力输运(SST)两方程湍流模型,开展了DLR-F6翼身组合体的高阶精度数值模拟研究。主要目的是确认WCNS模拟跨声速典型运输机构型的能力。采用粗、中、细3套网格开展了网格收敛性研究,从气动特性、压力分布、表面流态等方面研究了网格密度对DLR-F6翼身组合体气动特性的影响;采用中等网格开展了来流迎角对气动特性的影响研究。通过与试验数据、CFL3D软件和TRIP软件计算结果的对比,表明网格密度主要影响激波位置和压差阻力系数,同时对翼身结合部分离区大小有一定影响;采用高阶精度计算方法显著提高了气动力系数的模拟精度,力矩系数数值模拟结果与试验的差异有待进一步分析。     

ν-gap度量及其在飞行控制律评估中的应用

 传统控制律评估方法主要用于单输入单输出（SISO）系统，且对模型参数摄动考虑不够全面,针对这些不足，研究了ν-gap度量方法。在介绍系统广义稳定裕度相关概念的基础上，给出了ν-gap度量的定义、特点和性质以及近似摄动模型的计算，提出ν-gap度量评估控制律的步骤。实例结果表明，该方法不仅克服了上述传统评估方法的缺陷，而且还有根据所求的各摄动影响情况忽略影响小的元素，以减少计算量及可以找到最坏情况下的参数摄动组合等优点。     

Integrated Entry Guidance for Reusable Launch Vehicle

A method for the implementation of integrated three-degree-of-freedom constrained entry guidance for reusable launch vehicle is presented. Given any feasible entry conditions, terminal area energy management interface conditions, and the reference trajectory generated onboard then, the method can generate a longitudinal guidance profile rapidly, featuring linear quadratic regular method and a proportional-integral-derivative tracking law with time-varying gains, which satisfies all the entry corridor constraints and meets the requirements with high precision. Afterwards, by utilizing special features of crossrange parameter, establishing bank-reversal corridor, and determining bank-reversals according to dynamically adjusted method, the algorithm enables the lateral entry guidance system to fly a wide range of missions and provides reliable and good performance in the presence of significant aerodynamic modeling uncertainty. Fast trajectory guidance profiles and simulations with a reusable launch vehicle model for various missions and aerodynamic uncertain-ties are presented to demonstrate the capacity and reliability of this method.     

基于飞机油箱模型形状特征油量测量切片步长选择方法研究

 在分析飞机数字式油量测量过程中目前广泛使用的切片法油量测量原理的基础上,针对现有的定步长切片法无法得到准确、可靠的燃油质量特性数据库的缺陷,结合对飞机油箱模型形状特征的分析,提出了基于飞机油箱模型形状特征的油量测量切片步长选择方法。此方法包括切片步长整体和局部选择两个过程,整体选择以实现相邻两切片平面所夹油箱模型体积近似相等为目的来确定切片步长,以体现油箱模型截面整体变化规律;局部选择以设计切片平面与截面突变平面重合或尽可能接近的方式,突出油箱截面的局部变化特征。实验结果表明：该切片步长选择方法较定步长方法能够建立更为合理、可靠的燃油质量特性数据库,从而提高了油量测量精度。     

The effect of inlet conditions on the flow and heat transfer in multiple rotating cavity with axial throughflow

This paper discusses experimental results from two different build configurations of a heated multiple rotating cavity test rig.Measurements of heat transfer from the discs and tangential velocities are presented.The test rig is a 70% full scale version of a high pressure compressor stack of an axial gas turbine engine.Of particular interest are the internal cylindrical cavities formed by adjacent discs and the interaction of these with a central axial throughflow of cooling air.Tests were carried out for a range of non-dimensional parameters representative of high pressure compressor internal air system flows(Re up to 5×106 and Rez up to 2×105).Two different builds have been tested.The most significant difference between these two build configurations is the size of the annular gap between the(non-rotating) drive shaft and the bores of the discs.The heat transfer data were obtained from thermocouple measurements of surface temperature and a conduction solution method.The velocity measurements were made using a two component,LDA system.The heat transfer results from the discs show differences between the two builds.This is attributed to the wider annular gap allowing more of the throughflow to penetrate into the cavity.There are also significant differences between the radial distributions of tangential velocity in the two builds of the test rig.For the narrow annular gap,there is an increase of non-dimensional tangential velocity V/Ωr with radial location to solid body rotation V/Ωr=1.For the wider annular gap,the non-dimensional velocities show a decrease with radial location to solid body rotation.      

二维翼段颤振的μ控制

 采用超声电机作为作动器来实现含控制面的翼段颤振鲁棒抑制。针对作者设计的二维翼段颤振主动抑制系统，通过理论与实验相结合的方法，建立了考虑沉浮方向阻尼和作动器模型参数不确定性的控制系统模型，设计了μ控制器，并对控制器做了降阶处理。数值仿真和风洞试验表明，μ控制器可有效地抑制颤振的发生，将颤振临界速度提高23.4%。相对于H_∞控制器，μ控制器的控制效果和鲁棒性更好。     

航空发动机管路测量数据分割方法

 为了提高航空发动机管路测量数据的反求建模效率,提出了一种区域增长分割算法。该方法主要是利用管路表面在曲面索引系数映射、主曲率映射、高斯映射上的特性,并基于均值漂移算法和遗传算法提出了区域分割算法中种子区域的选择策略。然后利用区域增长分割方法实现了对管路测量数据的分割。经过仿真和实测数据验证,所提管路分割算法具有较好的分割质量和效率。     

Abnormal Shape Mould Winding

为解决网格化芯模的缠绕问题,本文提出了复合材料面片缠绕机理;接着详细分析了面片缠绕过程中的芯模凹曲面上纤维滑线和架空现象,应用微分几何曲面理论和空间几何理论,提出判据及其解决方案;最后,针对飞机发动机进气道的缠绕成型,编制缠绕控制程序并进行相应的实验,验证了面片缠绕方法的可行性。     

Interfacial Properties Modification of Carbon Fiber/ Polyarylacetylene Composites

This work was dedicated to performing surface oxidation and coating treatments on carbon fibers (CF) and investigating the changes of fiber surface properties after these treatments, including surface composition, relative volume of functional groups, and sur- face topography with X-ray photoelectron spectroscopy (XPS) and atom force microscopy (AFM) technology. The results show that, after oxidation treatments, interfacial properties between CF and non-polar polyarylacetylene (PAA) resin are remarkably modified by removing weak surface layers and increasing fiber surface roughness. Coating treatment by high char phenolic resin solution after oxida- tion makes interface of CF/PAA composites to be upgraded and the interfacial properties further bettered.     

Thermal Stresses and Cracks During the Growth of Large-sized Sapphire with SAPMAC Method

本文利用有限元法对冷心放肩微量提拉法(SAPMAC)生长大尺寸蓝宝石单晶过程中的热应力进行了模拟计算,应用临界缺陷理论解释了裂纹的萌生与扩展机理。研究表明,晶体内的热应力主要与晶体生长速率,环境温度和结晶取向有关;较大热应力多出现在结晶界面,放肩、收尾以及直径急剧变化等位置;最大热应力总是出现在籽晶与新生晶体的界面附近。裂纹将在临界缺陷位置产生,并在应力作用下沿a或m面扩展。计算结果与实验结果基本吻合,通过晶体生长系统和生长工艺的改进,对晶体的开裂问题得到了有效的抑制。     

Optimization Design System for Composite Structures Based on Grid Technology

To solve the topology optimization of complicated multi-objective continuous/discrete design variables in aircmit structure design, a Parallel Pareto Genetic Algorithm （PPGA） is presented based on grid platform in this paper. In the algorithm, the commercial finite element analysis （FEA） software is integrated as the calculating tool for analyzing the objective functions and the filter of Pareto solution set based on weight information is introduced to deal with the relationships among all objectives. Grid technology is utilized in PPGA to realize the distributed computations and the user interface is developed to realize the job submission and job management locally/remotely. Taking the aero-elastic tailoring of a composite wing for optimization as an example, a set of Pareto solutions are obtained for the decision-maker. The numerical results show that the aileron reversal problem can be solved by adding the limited skin weight in this system. The algorithm can be used to solve complicated topology optimization for composite structures in engineering and the computation efficiency can be improved greatly by using the grid platform that aggregates numerous idle resources.