基于格子Boltzmann方法表面形貌对微通道对流换热的影响

为推动微通道冷却技术在航空发动机中的应用,进一步提高航空发动机性能,本文通过扫描电子显微镜、粗糙度扫描仪等多方式扫描得到实际加工物体表面形貌,结合分形理论进行处理,得到真实、光滑、分形插值、康托尔集四类圆形截面管模型,管径100μm~400μm;采用格子Boltzmann方法进行数值模拟,首先通过与文献中实验结果的对比,验证了此方法的正确性,然后对空气在这四类圆形截面微通道中的对流换热特性进行了数值模拟,计算雷诺数Re范围80~640;结果发现：表面形貌是微通道中对流换热不可忽略的因素,相对粗糙度越大越有利于换热,真实扫描管的换热性能要比相同条件下光滑管高2.22%;分形理论可用于微通道表面形貌的构建,在相同迭代次数下,分形插值模型与真实扫描形貌换热性能相差0.03%,而康托尔集模型与真实扫描形貌换热性能相差1.88%,分形插值模型较康托尔集模型更能真实地反映实际形貌;微通道中,随着雷诺数的增加Nu将不再是常数,而是有增加的趋势。     

基于压缩量偏差约束的整体叶盘砂布轮数控抛光路径规划

数控抛光是实现整体叶盘叶片表面光整加工自动化、高效化的关键技术,而抛光路径规划对表面质量有着显著的影响。在抛光曲面时,砂布轮沿接触曲线上的径向压缩量偏差导致材料去除量一致性较差。针对该问题,提出一种基于压缩量偏差约束的路径规划方法。通过对砂布轮受曲面挤压变形几何分析,建立了自由曲面抛光接触线上砂布轮压缩变形量分布模型。在给定刀触点处,以接触曲线上最大压缩量偏差最小化为评价指标构建了刀轴矢量优化模型,并利用改进的粒子群优化（PSO）算法进行求解。基于最优刀轴矢量,完成自由曲面上等参数相邻刀触轨迹曲线的计算。最后,通过数据对比和试验,验证了所提方法的正确性和有效性。结果表明,抛光区域的轮廓度、粗糙度均得到明显的改善。     

非正交弧线齿面齿轮齿面设计及根切研究

非正交弧线齿面齿轮是以端面渐开线弧线齿圆柱齿轮为假想刀具包络展成的面齿轮。根据面齿轮啮合原理推导了非正交弧线齿面齿轮齿面方程,结合Matlab数学模型和根切理论研究了该种面齿轮的根切现象,计算出面齿轮不产生根切的最大内径,并通过Catia仿真滚齿来验证计算的准确性。该研究旨在提出一种新型非正交面齿轮,通过理论推导及仿真模拟来研究该种面齿轮的根切现象,在此基础上探究非正交弧线齿面齿轮的内径影响因素。研究结果表明:在同等参数下,当位置参数增大时,刀具渐开线截面高度和内径逐渐减小,当刀具圆弧半径增大时,内径逐渐减小,刀具渐开线截面高度逐渐增大。     

带非线性支撑的转子有限元模型求解方法

用数值方法研究了非线性支撑的柔性转子系统的动学行为,提出了一种将有限元与非线性支撑结合的模型和求解方法。利用有限元法(FEM)构建转轴和转盘部分的模型,通过矩阵进行组合;利用离散元方法对包含滚动轴承和挤压油膜阻尼器(SFD)的支撑部分进行建模,此部分包含4个单元,分别为轴承内圈、外圈、SFD内圈和支撑鼠笼。有限元部分和离散元部分通过轴端节点相连,仿真过程中轴端位移传递给非线性支撑部分,支撑部分通过位移计算得到的非线性力反过来作用于有限元转子轴端部分。为了耦合求解有限元转子和非线性支撑组成的数学模型,提出了一种综合的迭代求解方法,克服传统的有限元求解方法对轴端隐性非线性支撑的求解局限性。由于转轴部分采用了Timoshenko梁单元建模,对比与简单转子模型,可以考虑陀螺力矩和轴的柔性特征,更能体现非线性支撑对振动真实影响。在建立的20个轴单元的有限元转子模型中,非线性响应更多体现在靠近非线性支撑的节点1和节点21处,响应频谱中靠近轴端的节点能体现出滚动轴承的2倍和3倍变柔振动频率。     

发动机唇口电热防冰系统性能仿真

采用三维耦合传热方法对发动机唇口电热防冰系统进行了多状态下的性能仿真,获得了唇口表面平衡温度及动态温升响应时间,以此评估防冰系统性能,并通过与冰风洞试验结果对比,对仿真过程进行了修正。研究表明:唇口表面加热区平衡温度均高于273.15 K,温升响应时间小于60 s,满足防冰需要,较为严酷的防冰状态点主要为飞行速度较大、环境温度较低的状态点;修正唇口材料导热系数后,仿真与试验结果吻合良好。研究结果能有效指导工程研制,对提高防除冰仿真分析结果的准确性提供了很好的借鉴。     

Evolution of periodic orbits near the Lagrangian point L2

A study of the evolution of the periodic and the quasi-periodic orbits near the Lagrangian point L₂, which is located to the right of the smaller primary on the line joining the primaries and whose distance from the more massive primary is greater than the distance between the primaries, in the framework of restricted three-body problem for the Sun–Jupiter, Earth–Moon (relatively large mass ratio) and Saturn–Titan (relatively small mass ratio) systems is made. Two families of periodic orbits around the smaller primary are identified using the Poincaré surface of section method – family I (initially elliptical, gradually becomes egg-shaped with the increase in the Jacobi constant C and elongated towards the more massive primary) and family II (initially egg-shaped orbits elongated towards L₂ and gradually becomes elliptical with the increase in C). The family I in the Sun–Jupiter and Saturn–Titan systems contains two separatrix caused by third-order and fourth-order resonances, while the Earth–Moon system has only one separatrix which is caused by third-order resonances. Also in the Sun–Jupiter and the Saturn–Titan systems, family I merge with family II, around Jacobian constant 3.0393 and 3.0163, respectively, while in the Earth–Moon system, family II evolves separately from two different branches. The two branches merge at C = 3.184515. In the Earth–Moon system, the family II contains a separatrix due to third-order resonances which is absent in the other two systems.     

Numerical Simulation of Rotor Flow Field Based on Overset Grids and Several Spatial and Temporal Discretization Schemes

A numerical method based on solutions of Euler/Navier-Stokes (N-S) equations is developed for calculating the flow field over a rotor in hover. Jameson central scheme, van Leer flux-vector splitting scheme, advection upwind splitting method (AUSM) scheme, upwind AUSM/van Leer scheme, AUSM+ scheme and AUSMDV scheme are implemented for spatial discretization, and van Albada limiter is also applied. For temporal discretization, both explicit Runge-Kutta method and implicit lower-upper symmetric Gauss-Seidel (LU-SGS) method are attempted. Simultaneously, overset grid technique is adopted. In detail, hole-map method is utilized to identify intergrid boundary points (IGBPs). Furthermore, aimed at identification issue of donor elements, inverse-map method is implemented. Eventually, blade surface pressure distributions derived from numerical simulation are validated compared with experimental data, showing that all the schemes mentioned above have the capability to predict the rotor flow field accurately. At the same time, vorticity contours are illustrated for analysis, and other characteristics are also analyzed.     

Mesh-dependence of Material with Softening Behavior

The issue of mesh-dependence emerges when the conventional continuum damage model is applied to handling the softening behavior. In order to circumvent the mesh-dependence, the non-local theory is introduced into the conventional damage model and the finite element formulas are derived for two-dimensional gradient-enhanced damage model. A new element is proposed in which the basic unknown quantities are displacement, non-local equivalent strain and the gradient of non-local equivalent strain. The element and constitutive equation proposed in this article are added to the finite element software ABAQUS through user subroutine UEL. Numerical results show that the gradient-enhanced damage model can eliminate the mesh-dependence and is effective for dealing with the issue of softening behavior.     

Sliding mode control of reaction flywheel-based brushless DC motor with buck converter

Reaction flywheel is a significant actuator for satellites’ attitude control. To improve output torque and rotational speed accuracy for reaction flywheel, this paper reviews the modeling and control approaches of DC-DC converters and presents an application of the variable structure system theory with associated sliding regimes. Firstly, the topology of reaction flywheel is constructed. The small signal linearization process for a buck converter is illustrated. Then, based on the state averaging models and reaching qualification expressed by the Lee derivative, the general results of the sliding mode control (SMC) are analyzed. The analytical equivalent control laws for reaction flywheel are deduced detailedly by selecting various sliding surfaces at electromotion, energy consumption braking, reverse connection braking stages. Finally, numerical and experimental examples are presented for illustrative purposes. The results demonstrate that favorable agreement is established between the simulations and experiments. The proposed control strategy achieves preferable rotational speed regulation, strong rejection of modest disturbances, and high-precision output torque and rotational speed tracking abilities.     

提高飞行器升力的若干方法研究

本研究用数值模拟方法研究了四种增升方法。数据模拟的出发方程是可压流的Ｎ－Ｓ方程，格式为Ｂｅａｍ－Ｗａｒｍｉｎｇ格式的改进型。通过上述方法分别研究了：（１）动界效应对提高升力的作用；（２）涡的有利干扰对主升力面的影响；（３）吸气的驻涡增效应；（４）升力面下表面向下喷气的增升作用。研究结果表明，采用上述四种方法可有效地提高飞行器的升力，对改善飞行器的性能有着十分重要的意义。