结冰风洞中过冷大水滴云雾演化特性数值研究

为明晰结冰风洞中过冷大水滴（SLD）云雾演化特性,发展了基于欧拉法的SLD液滴运动、传热和传质耦合计算方法,针对3 m×2 m结冰风洞主试验段水平收缩构型,分析了SLD云雾沉降收缩特性、动量平衡特性和热平衡特性,探索了液滴变形破碎的影响,评估了构型出口处SLD液滴动量平衡和热平衡状态。研究结果表明：直径超过250 μm的SLD液滴在构型内会发生显著形变,液滴尺寸越大则变形程度越强,尤其在160 m/s工况下,当液滴直径超过750 μm后,SLD液滴甚至会发生破碎;液滴变形破碎效应会增大液滴加速度和液滴温度下降率,促使SLD液滴趋近动量平衡和热平衡状态;SLD云雾（最大液滴直径小于1 000 μm）在构型出口处会出现显著的粒径浓度分层、动量分层和热分层现象,其中直径小于100 μm的小尺寸液滴速度快、温度低且不断凝结,趋于平衡态,但直径超过500 μm的大尺寸SLD液滴速度慢、温度高且不断蒸发,则显著偏离平衡态;增大试验段气流速度尽管会减弱SLD云雾粒径浓度分层程度,但会增强动量分层和热分层程度,尤其在160 m/s工况下,SLD云雾会均匀分布在构型出口中心区域内（-0.75 m < Y < 0.75 m且-0.5 m < Z < 0.5 m）,与其平衡态间的最大速度差和温度差将分别超过18 m/s和20℃。     

基于代理模型全局优化的自适应参数化方法

对于翼型气动隐身设计问题,设计变量的配置对设计结果影响很大,而简单地增加设计变量不能保证得到理想的结果。提出一种适用于代理模型全局优化的自适应参数化方法：利用全局敏感性分析方法——基本效应法,得到设计空间关于目标函数的敏感区域信息,并以此为根据增加设计变量;利用节点插入算法将低维样本在高维空间内进行重构,避免了重新取样的工作量。相对于传统固定设计空间维度方法,自适应参数化方法在设计空间的敏感区域扩展维度,能够更加精准地描述外形并反映目标的变化趋势。通过飞翼布局翼型的气动隐身优化算例,证实自适应参数化方法可以大幅提高优化设计质量和效率。     

变循环发动机过渡态性能直接模拟方法

为了实现多变量可调的变循环发动机（VCE）过渡态性能模拟和控制规律设计,在稳态逆算法和过渡态性能隐式格式计算方法的基础上开发了变循环发动机过渡态直接模拟方法。在稳态逆算法模型中加入容积效应和转子动力学方程,实现了同时给定加速率、涡轮前温度和压缩部件工作点的条件下直接模拟计算过渡态过程中变几何参数和燃油流量的调节规律,验证了方法的精度和可行性,并将该方法用于变循环发动机转模态性能模拟和控制规律设计。计算结果表明,过渡态直接模拟方法的误差在0.58%以内,超声速巡航状态下由单外涵模态转换到双外涵模态的时间约为1 s,海平面静止状态下双外涵转单外涵的时间约为3 s,且推力、转速、涡轮前温度、喘振裕度等参数过渡平稳。该方法可简化转模态控制规律设计流程,并提高设计精度。     

高温铝合金电磁超声检测回波特性及因素分析

针对高温铝合金在线检测条件下,温度对铝合金电磁超声检测回波特性的影响规律尚不明确、高温检测时缺陷定量/定位补偿困难这一难题,以螺旋线圈电磁超声换能器（EMAT）为例,建立了高温铝合金EMAT检测过程的场路耦合有限元模型;研究了温度对EMAT激励/接收换能效率、EMAT激励/接收电路的功率分配特性、超声传播过程中的扩散/介质衰减特性、回波幅值和超声声速等因素的影响规律;研制了耐高温EMAT探头,对20~500℃高温铝合金试样进行了检测实验,并测定了高温铝合金的超声介质衰减系数和超声声速。在仿真和实验相结合的基础上,分析了高温检测时超声回波幅值变化特性及其影响因素。结果表明：对于铝合金这类非铁磁性金属材料,导致高温时超声回波幅值下降的主要原因是超声介质衰减系数随着温度的升高而增大,其次为高温时EMAT激励/接收电路的功率分配特性的改变。在激励EMAT在试样表面形成的洛伦兹力不变的条件下,其所激励的超声波回波幅值具有随着温度的增加而增加的特点,可以有效减缓超声回波幅值下降的趋势。     

Multiscale topology optimization using feature-driven method

This paper presents a multiscale design method for simultaneous topology optimization of both macrostructures and microstructures. Geometric features are extended as design primitives at both macro and micro scales and represented by Level Set Functions (LSFs). Parameters related to the locations, sizes, and orientations of macro and micro features are considered as design variables and optimized simultaneously. In the overlapping areas of different macro features, embedded microstructures are optimally figured out as the solution of the corresponding sub-optimization problem. In this study, the eXtended Finite Element Method (XFEM) is implemented for structural and sensitivity analyses with respect to design variables. This method has the advantage of using a fixed grid independent of the topology optimization process. The homogenization procedure is applied to calculate the effective properties of considered microstructures in each macro feature. Numerical examples are presented to illustrate the effectiveness of the proposed method. Results depict that the multiscale design cannot obviously improve structural stiffness compared with a solid-material design under the linear elastic condition.     

Effect of bristle pack position on the rotordynamic characteristics of brush-labyrinth seals at various operating conditions

Over the last few decades, the research on the effect of bristle pack position on the rotordynamic characteristics of the brush-labyrinth seals is not sufficient. To this end, two kinds of brush-labyrinth seals for the bristle pack element installed upstream of the labyrinth teeth named BSU and installed downstream of the labyrinth teeth called BSD were used to investigate the effect of bristle pack position on the rotordynamic characteristics of the brush-labyrinth seals. Using the numerical model combining the porous medium model and the whirling rotor method, the rotordynamic characteristics of the BSU and BSD at various operating conditions including four kinds of pressure ratios, five kinds of inlet preswirl speeds and four kinds of rotor spinning speeds were conducted. The obtained results show that the effects of operating conditions on rotordynamic coefficients for the different seal configurations are different. The direct stiffness, cross-coupled stiffness and direct damping of the BSU are lower than those of the BSD. The rotordynamic coefficients of the BSU are more insensitive to the operating conditions variation. From the perspective of the seal stability, the BSU is a better brush-labyrinth seal configuration at high pressure ratio, high positive preswirl or high rotor spinning speed conditions. While in the case of low pressure ratio, low positive preswirl or low rotor spinning speed conditions, the BSD is a better choice.     

民用飞机颤振试飞技术研究

为了能够顺利完成ARJ21-700飞机颤振试飞,研究民用运输类飞机颤振试飞的依据与要求,提出民用运输类飞机颤振试飞的测试方法、试飞方法、激励方法和数据处理方法;结合ARJ21-700飞机合格审定试飞进行试飞验证,给出该型飞机颤振特性试飞的符合性验证结果。结果表明:该套方法适合民用运输类飞机合格审定颤振试飞。     

基于LS-DYNA滚动轴承使役可靠度计算方法

利用软件LS-DYNA对滚动轴承使役过程进行了动态有限元计算,基于响应面法（RSM）建立6408滚动轴承输入参数与输出响应的函数关系方程,设定各类随机误差分布,利用矩阵试验法确定抽样点在输入变量抽样空间的位置,进行一系列确定性拟合试验.每次试验借助LS-DYNA进行精确求解,运用求解数据对方程进行最小二乘法回归分析确定方程的组成项及系数.最后利用函数关系方程计算轴承任意时刻的使役可靠度.结果表明:计算得出该轴承在0.020008s时的使役可靠度为97.062%.      

基于库伦摩擦效应的箔片轴承极限承载力理论与试验

基于有限单元法建立了考虑库伦摩擦的波箔型径向气体箔片轴承的箔片结构模型,采用有限差分法和有限单元法耦合求解Reynolds方程和气膜厚度方程,通过求解轴颈达到极限偏心率时的轴承极限承载力,研究了箔片结构库伦摩擦效应对轴承极限承载力的影响规律,并搭建了轴承极限承载力测试试验台,利用温度法测量了两个具有不同轴承壳内表面粗糙度的波箔型径向气体箔片轴承的轴承极限承载力.通过对比分析仿真结果与试验结果表明:轴承壳圆柱孔内表面粗糙度为0.4μm的轴承在10000r/min和20000r/min下,轴承极限承载力分别为15.5N和42.3N;而表面粗糙度为1.6μm的轴承极限承载力为10.9N和29.6N,这是由于波纹箔片和轴承壳体之间的库伦摩擦力增大了波纹箔片的刚度,因此增大箔片结构摩擦因数使得轴承极限承载力降低,并且仿真结果变化趋势和试验结果变化趋势吻合.      

考虑旋翼调制影响的直升机RCS特性分析及评估

考虑旋翼调制的影响,建立了适合于直升机全机雷达散射截面(RCS)特性计算的“面元-边缘法”.兼顾计算效率和精度,在外形变化剧烈地区域,网格进行加密处理,并保证电磁网格密度和尺度满足雷达波波长大小的比例关系.以桨盘倾倒方式计入桨叶的挥舞和变距运动,采用准静态法模拟旋翼对全机RCS特性的影响.分析了某直升机RCS极化、姿态、频率响应特性,并根据直升机RCS和雷达探测距离关系,提出了4级预警机制和角域范围.研究表明:旋翼转动时全机RCS动态响应具有连续性和对称性,振荡区散射水平强,RCS幅值为-5~12dB ·m2;奇数片桨叶比偶数片桨叶的RCS减缩2~5dB ·m2,且有利于控制直升机RCS包络线和散射峰值的时域响应,增强雷达隐身性能.