A FORMAL REPRESENTATION FOR FEATURE-BASED DESIGN

从三个方面研究了基于特征设计的形式化描述方法：基于特征设计的基本定义、设计过程模型及设计算法。文中给出了特征分类、特征操作、特征模型验证和特征映射策略及其形式化描述方式，建立了兼顾并行工程和顺序工程设计方式的基于特征设计过程形式化描述模型，提出了基于特征的设计算法，为解决基于特征设计系统开发和运行缺乏统一方法学这一课题作了有益的尝试。     

可重复使用运载器级间分离过程研究

级间分离过程中,可重复使用运载器受到非定常扰动,可能会造成分离失败。为研究级间分离过程中可重复使用运载器运动特性及周围流场分布随分离高度的变化情况,建立航天器和运载器的运动方程,使用动网格技术和有限体积法求解N-S方程。研究结果表明：随着级间分离高度的升高,气动力不断减小,非定常扰动对分离产生的影响越来越小;分离过程中,航天器和运载器的平动速度差较小,角速度差越来越小,最大差值从30 km高空的0.427 rad/s减小到50 km 高空的0.1 rad/s,二者之间的最小距离变大,更有利于安全分离。     

开关磁阻电机非线性磁参数建模方法

建立了一种基于MATLAB/SIMULINK的非线性磁参数法的开关磁阻电机数学模型,在电磁场有限元(FEM)数值计算结果的基础上,解算出电机磁特性函数矩阵ψ(i,θ),把ψ(i,θ)函数矩阵三次样条插值变换成电机电流特性矩阵i(ψ,θ),再利用MATLAB/SIMULINK中的查表模块建立了开关磁阻电机系统的非线性模型。该模型中的电机本体模块、功率变换器模块、控制模块等都是通过SIMULINK中PSB通用功能模块搭建起来的,修改容易且直观性强,试验结果验证了该模型的准确性。     

低空风切变系统建模及其对直升机飞行安全威胁定性分析

系统地建立了含湍流的三维低空风切变模型,并综合分析了其对直升机飞行安全的威胁。建立了微下击暴流风场及大气湍流场组合的风切变模型,在不增加计算量的前提下,选取特征点发展直升机飞行动力学模型,有效捕捉了风切变的切变特性,提高了在风场中的动态响应计算精度。模型配备姿态保持功能的控制增稳系统已符合一般直升机的飞行状态,并改善机体响应。根据风速分布的特点,选取不同飞行速度、不同风场位置进行飞行仿真,定性地对比分析状态量变化与风场对应关系,并且以垂向通道为例,从动力学角度分析验证了响应的理论计算表达式。结果表明：湍流主要导致高频姿态角响应,风切变对飞行状态量变化占主导作用,且垂向风是引起威胁的主要因素,据此提出危险风场规避建议。     

多变量线性时变系统的特征模型及自适应模糊控制方法

对于阶数及参数未知的高阶系统，如何设计一个工程上易于实现的低阶控制器，以实现有效的高性能控制，这是目前航天等领域亟待解决的问题。该文采用将对象特征和控制性能要求相结合的方法，对多输入一多输出高阶线性时变系统和一类非线性时变系统从理论上详细推导出了其特征模型，在此基础上设计了基于特征模型的稳定的自适应模糊广义预测控制方案。由于控制结构中使用了分层模糊系统，从而减少了模糊规则数目。通过对一个挠性航天器控制的仿真实验验证了所给方法的有效性。     

A transient bulk flow model with circular whirl motion for rotordynamic coefficients of annular seals

Bulk flow model with perturbation simplification has been used to calculate rotordynamic coefficients in annular seals which have significant influences on the dynamic behavior of rotors in turbomachinery. In this work, a transient bulk flow model with arbitrary rotor motion is developed, and the boundary conditions and friction factor in the model are calibrated with steady Computational Fluid Dynamics (CFD) analysis. The numerical solution scheme is developed based on the finite element method to obtain the transient reaction force in the seal clearance. With a periodic circular rotor orbit, the transient forces at multiple whirling frequencies are used to evaluate the rotordynamic coefficients. The leakage flowrate of CFD analysis has good agreement with experimental results and the calibrated parameters in bulk flow model are dependent on operating conditions. Although CFD calibration improves the accuracy of the perturbed bulk flow model, the direct damping is overestimated and the cross-coupled damping is underestimated. Compared with the perturbed model, the predictions of the transient bulk flow model are more agreeable with the experiment.     

Interface flux reconstruction method based on optimized weight essentially non-oscillatory scheme

Aimed at the computational aeroacoustics multi-scale problem of complex configurations discretized with multi-size mesh, the flux reconstruction method based on modified Weight Essentially Non-Oscillatory (WENO) scheme is proposed at the interfaces of multi-block grids. With the idea of Dispersion-Relation-Preserving (DRP) scheme, different weight coefficients are obtained by optimization, so that it is in WENO schemes with various characteristics of dispersion and dissipation. On the basis, hybrid flux vector splitting method is utilized to intelligently judge the amplitude of the gap between grid interfaces. After the simulation and analysis of 1D convection equation with different initial conditions, modified WENO scheme is proved to be able to independently distinguish the gap amplitude and generate corresponding dissipation according to the grid resolution. Using the idea of flux reconstruction at grid interfaces, modified WENO scheme with increasing dissipation is applied at grid points, while DRP scheme with low dispersion and dissipation is applied at the inner part of grids. Moreover, Gauss impulse spread and periodic point sound source flow among three cylinders with multi-scale grids are carried out. The results show that the flux reconstruction method at grid interfaces is capable of dealing with Computational AeroAcoustics (CAA) multi-scale problems.     

图形交互自动编程及对三维曲面的处理

作为一种新型的计算机辅助编程技术，图形交互自动编程正在走进更多的企业。本文介绍了这种编程方法的特点和基本步骤，重点阐述了利用这种办法处理传统编程方法不易解决的三维曲面编程问题。     

工序问题的数学模型

本文讨论了加工的次序问题，对于不同的、具有实际意义的问题，我们得到了相应的优化准则，这些准则在实际应用中便于操作，利用这些优化准则得到了加工次序问题的最优的结果。     

欧洲计算流体力学——个人观点

介绍了70年代初到20世纪末欧洲科学家对CFD发展的主要贡献,并重点关注了对国际流体力学数值方法会议(ICNMFD)和国际计算力学学术会议(ISCFD)的成功举办具有重大贡献的科学家.ICNMFD始于1969年(新西伯利亚),由参与了美国和苏联太空竞赛的两国科学家联合创办,而ISCFD则作为ICNMFD在亚太地区的同类会议初创于1995年(东京).由Oshima教授和本人提议,得益于上述两个会议全体科学委员会成员的共同努力,ICNMFD和ISCFD自2000年开始,合并成为每两年举办一次的国际计算流体力学会议(ICCFD).绝大部分本文讨论所涉及的学者通常都在受邀进行大会报告后,成为了前述三大国际会议中的一个或多个科学委员会的成员.作者对文中讨论所涉及的研究人员的圈定负有全权责任,上述人员的选择更多地基于应用而非理论计算流体力学,因此尚存在不完善之处.此外由于文中所讨论的学者大都是作者在40多年计算流体力学研究生涯中的私交好友,因此并没有将年轻一代的科学家考虑进来,但这并非由于他们的贡献不重要.另外值得注意的是欧洲的计算流体力学同美国的计算流体力学有非常紧密的联系,本文讨论所涉及的研究学者有些最早从欧洲开始涉足计算流体力学而现在则生活和工作于美国(例如Glowinski,Roe,Van Leer等人),而有些则是在美国取得了博士学位(例如Launder,Napolitano),他们都同美国的科学家进行过合作,因此很难清晰地对欧洲和美国的计算流体力学进行界定.对那些现已移居美国的欧洲科学家们,本文优先考虑了他们移居之前的贡献.最后但并非是最不重要的,作者通过电子邮件一一联系了文中所涉及的科学家,希望他们本人能够给出其自身贡献的评价草稿,从而帮助作者以最权威的视角在文中相应章节对其工作进行评述.本文向提供自述评价的诸位致以谢意,同时也希望其他同仁能够对本人不恰当的擅自评述给予谅解.