基于直升机/涡轴发动机综合仿真平台的发动机非线性模型预测控制

 基于具有可靠置信度的直升机/涡轴发动机综合仿真平台,研究了涡轴发动机带约束优化的非线性模型预测控制(NMPC)技术。首先通过设计多输出迭代约简最小二乘支持向量回归机(RRLSSVR),训练具有较好实时性、精度及泛化能力的内嵌式预测模型,在高度0~5 km、前飞速度0~75 m/s范围内模型精度达5‰。其次,考虑了扭矩、燃油流量、动力涡轮转速、燃气涡轮转速等综合信息及相关约束对控制效果的影响,利用在线序列二次规划(SQP)算法实现滚动优化控制,而后加入目标转速偏差的积分项以消除静差,保证输出恒定。最后,通过对直升机进行机动飞行大扰动仿真验证了该预测控制器对扰动的抑制能力,相比传统串级PID控制,能够显著降低动力涡轮转速下垂/超调量,达到更好的控制品质。     

阻尼不相似动力学模型的时域响应修正方法

 动力学相似的缩比模型试验中,阻尼往往难以保证相似,使得缩比模型的动响应测试数据难以直接换算到结构原型上,因此提出一种阻尼不相似动力学模型的时域响应修正方法来解决这一问题。该修正方法假定缩比模型的阻尼可用比例阻尼模型近似表征,并直接从实际缩比模型与理想缩比模型的一般强迫振动响应计算方程出发,利用线性系统的叠加原理和模态叠加法,将动响应修正量的求解转换为理想模型在一个等效附加激励力作用下的时域响应求解。同时,针对实际工程中响应测点数目有限的问题,利用模态缩聚法进行了未测量点的响应反演。该修正方法仅需已知结构原型和实际缩比模型的模态阻尼比,以及准确建立的实际缩比模型的质量矩阵和刚度矩阵,即可实现实际缩比模型在任意激励工况下的测试动响应修正。以某型飞机的翼梁缩比模型为研究对象,对所提出的响应修正方法进行了验证。试验和计算结果对比分析表明,修正后的响应功率谱密度(PSD)和响应的均方根(RMS)值与理想模型基本一致,表明了本文方法的可行性和有效性。     

安装角对压气机叶栅气动噪声特性的影响

 基于大涡模拟(LES)和边界元方法对轴流压气机叶栅湍流流场以及流场诱导的噪声进行计算,在不同叶栅安装角下研究来流攻角和来流雷诺数对叶栅气动噪声产生、辐射的影响。研究表明:来流雷诺数不变时,同一安装角下,随着来流攻角从-5°~20°变化,叶栅监测曲线上的声压级先减小后增大,在0°来流攻角下声压级达到最小。安装角为45°时,外场总声压级随来流攻角的分布与30°安装角变化趋势相近。但安装角为60°时,总声压级的变化则明显变缓。在0°来流攻角下,总声压级比安装角为30°和45°时增加了近6 dB,但在其他正来流攻角下,变化并不明显。叶栅的最小声压值出现在弦线方向附近,安装角改变时,最小声压级出现的位置也不同。安装角不变,随着来流雷诺数的增大,叶栅表面的分离减小,损失降低。但叶栅表面的压力脉动随着来流雷诺数的增大而增大,使外场辐射噪声增加。     

IIR数字滤波器的序列约束Chebyshev设计方法

 进行无限冲击响应(IIR)数字滤波器的无约束Chebyshev设计时,得到的滤波器幅值响应通常在过渡带存在明显过冲现象,且通带边缘的群延迟误差较大。本文提出一种序列约束Chebyshev(SCC)方法,在通带和阻带最大频率响应(FR)误差小于某给定值ρ的条件下,使过渡带频率响应误差(相对某个具有单调幅值特性的过渡响应函数)的Chebyshev模最小,并通过收缩和对分找出ρ的最小值。该方法不但可以使通带和阻带最大频率响应误差尽可能小,还能减小过渡带幅值响应的过冲以及带边群延迟误差。仿真实例表明了该方法的有效性。     

悬停状态直升机桨叶扭转分布的优化数值计算

建立了一套基于高精度计算流体力学(CFD)技术和代理模型优化算法的旋翼气动外形设计方法。在该方法中,旋翼流场气动性能的计算采用了基于Navier-Stokes/Euler方程的CFD方法,并根据流场特点、精度和效率的要求采用Baldwin-Lomax(B-L)湍流模型,通量计算采用Roe-MUSCL格式进行。为提高网格生成质量和便于流场控制方程的求解,将流场分成两个区域,即围绕旋翼的黏性区和无黏的背景网格区。其中,桨叶网格使用了基于二维翼型网格的参数化方法生成,数值计算结果表明该方法有效地提高了网格生成质量及效率。在参考旋翼流场及桨叶细节流动分析的基础上给出设计变量及范围,有效减小了优化问题的规模;为满足优化和机理分析的需要,将基于置换遗传算法优化的拉丁超立方(PermGA LHS)方法和径向基函数(RBF)的代理模型优化方法引入到桨叶外形的优化设计中。首先以Helishape 7A旋翼为算例,检验了数值模拟方法的准确性。然后,应用所建立的优化方法针对旋翼负扭转分布进行了优化计算,结果表明优化后的旋翼悬停气动性能比优化前有了明显提高。     

热障涂层喷涂质量微焦点CT检测

采用微焦点X射线源对热障涂层样品进行了高分辨率X射线成像研究,采用多尺度对比度增强算法提高了热障涂层与基体之间的对比度.X射线图像显示热障涂层各处厚度差异较大,热障涂层内部可能存在大量孔隙.通过微焦点计算机层析成像(CT)对样品进行了高分辨率CT扫描与重建,得到了样品的CT重建图像.从热障涂层CT图像观察到热障涂层内部存在大量孔隙,且不同层孔隙尺寸、数量不同.结合3个坐标方向的切片图像分析了热障涂层深度方向的热障涂层厚度,得到了热障涂层厚度的二维分布,测量得到热障涂层平均厚度约为40μm.对沿热障涂层深度方向850~1131层切片图像,利用可视化软件VGstudio MAX得到了热障涂层内部孔隙三维分布.研究结果表明:微焦点CT可用于热障涂层制备工艺质量检测.      

复合材料与金属胶接结构的次弯曲效应研究

复合材料与金属胶接单搭接的次弯曲效应引起的剥离应力,会大大降低接头的强度。采用有限元方法分析复合材料铺向角、铺层顺序及斜削对次弯曲的影响,并进一步研究分析次弯曲和胶层剥离应力、板的弯曲刚度之间的关系。结果表明：次弯曲影响剥离应力的大小;弯曲刚度越大,次弯曲越小;斜削有助于减小次弯曲。     

基于时间谱方法的振荡翼型和机翼非定常黏性绕流数值模拟

 传统的双时间方法在非定常计算中长时间的过渡迭代推进求解导致其计算效率相对较低,针对周期性非定常流动问题的流动特征,发展了一种基于离散傅里叶变换的高效时间谱方法,用于求解振荡翼型和机翼的非定常黏性绕流。在时空耦合的雷诺平均Navier-Stokes(RANS)方程的求解中,对流项的离散应用了Roe的通量差分格式,物理时间项的离散方法为时间谱方法,伪时间推进采用了隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式。考虑到湍流的时空耦合效应,时空耦合的Spalart-Allmaras一方程湍流模型的物理时间项同样采用时间谱方法进行离散。为了进一步提高计算效率,当地时间步长和多重网格技术等加速收敛的措施均被采用。算例对俯仰振荡NACA0012翼型和Lann机翼的周期性非定常流场进行了数值计算。结果表明:对于周期性非定常流场的数值模拟,相比于传统的双时间方法,用时间谱方法近似物理时间项,不仅能够提高流场的计算精度,而且更能够大幅度提高计算效率。     

纤维桥连对复合材料Ⅰ型层间断裂韧性的影响

 采用ASTM D 5528标准对T700/9916复合材料单向板进行了试验,研究了双悬臂梁(DCB)试验中桥连纤维对Ⅰ型层间断裂韧性的影响。根据试验观察,提出桥连区存在单(s区域)、双纤维桥连区域(d区域)的分析模型,其中d区域内一点上桥连纤维的数量是s区域的两倍;通过对桥连纤维和悬臂梁相互作用的分析,推导出呈指数形式递减的单根纤维桥连力-转角函数。由此得到整个悬臂梁的桥连力函数,并对DCB试验进行了分析,得到的模拟曲线与试验曲线吻合较好。研究结果表明本文提出的纤维桥连模型能够真实地反映桥连纤维在裂纹扩展过程中的作用方式。     

水滴撞击特性的粒子统计法研究

基于水滴撞击特性研究的拉格朗日方法,提出了适合于考察非定常情况的粒子统计法.通过统计非定常流场中不同时刻撞击在防冰表面上不同区域的粒子数,获得防冰表面上粒子浓度分布随时间的变化情况,与基于染色法水滴撞击试验获得的色度分布具有直接对应关系;采用基于染色法的动态图像测量方法对某型发动机帽罩缩比模型开展了水滴撞击特性试验,对所提粒子统计法进行了验证.基于拉格朗日方法的粒子统计法获得的计算结果与试验结果符合较好,与采用拉格朗日轨迹法获得的结果具有较好的一致性.粒子统计法所模拟的物理过程与试验过程一致,在处理湍流随机性和来流液态水含量分布不均等方面具有优势.