基于非确定性的翼型鲁棒反设计方法

实际的翼型是在一个非确定性的状态范围内工作的,为此应用鲁棒控制理论建立了翼型的鲁棒气动反设计方法及数值求解过程,基本的优化器为基于伴随方程的位流优化方法。对守恒型位流方程应用高效AF 2格式求解,对非守恒型的伴随方程采用了旋转差分和高效AF 3格式,提高了梯度计算的精度和效率。优化算例表明:本文方法所得到的翼型在一个较宽的工作状态范围内都具有良好的性能。     

Development of a Synchronization Method for Fluid-Thermal Study of Hypersonic Flow

A synchronization method is developed for the fluid-thermal study of hypersonic flow.Different from conventional loosely/tightly coupled methods which separately deal with the flow field and the structure temperature field,the presented method expresses the governing equations in a unified framework so that the two fields can be calculated simultaneously.For efficiently solving the unified equations,the finite volume method together with the dual-time stepping approach is employed.Like in the flow field,the local time step is also used in the temperature field,which is determined from thermal conductivity spectral radii.In order to treat the fluid-structure interface more conveniently,an expanded virtual boundary is introduced.For validation,several fluid-thermal hypersonic flow problems are simulated.The computed results are compared with those obtained from the coupled methods and the experiment.In the continuous heating problems,the stagnation temperatures predicted by both the coupled and synchronization methods are in good agreements with the experimental data.In the unsteady flowthermal hypersonic flows,the stagnation heat fluxes predicted by the presented method and tightly coupled method are basically the same,which agree better with the experimental data than those predicted by the loosely coupled method.In terms of prediction of the stagnation temperature,the synchronization method shows better accuracy than the tightly coupled method.     

耦合传导/辐射情况下的传热系数和表面热流辨识

耦合传导/辐射情况下的传热参数和表面热流辨识问题是一类特殊的传热逆问题,在工程上应用较广。本文首先建立了耦合传导/辐射的传热问题的数值方法,进行了算例验证。然后分别建立了同时辨识材料热物性参数和辐射吸收系数的算法、辨识非稳态表面热流密度的算法;分析了测点位置、测量误差、辐射吸收系数对辨识结果的影响。结果表明:同时辨识材料热物性参数和辐射吸收系数时,利用右端测点的温度辨识出热传导系数值的精度较高,利用中间测点温度辨识出的辐射吸收特性系数值精度较高。对于表面热流的辨识,测点越靠近加热面,辨识结果的精度越高;吸收系数对中间测点的辨识结果影响大于其他位置测点的辨识结果。     

多层壁瞬态温度分布计算

风挡玻璃电防冰和除冰装置、飞机机翼及直升机旋翼电除冰装量都是有内加热的多层壁。这些多层壁中温度分布计算不仅在这些装置的设计计算中是必需的，而且在考虑多层壁的热变形及热应力时也是必要的。本文用数值法对有内加热的多层壁温度分布进行了研究，根据导热微分方程及各种边界条件，用克兰克－尼可尔辛法推出了典型结构的有限差表示式，然后用高斯－赛德迭代法求解线性方程组，得到了可用于多层壁瞬态温度分布的计算软件。本文用该软件对实际的多层金属壁及多层玻璃作了计算。计算表明，计算结果与文［３—５］所得结果相符。本文所述的方法可推广到二维及三维瞬态多层壁温度分布计算，也可用于无内加热层的多层壁温度计算。     

相变材料热控系统内部接触热阻的辨识方法研究

在相变材料热控系统的吸热融化过程中,由于材料与容器壁面之间的空穴效应形成的接触热阻,将对熔化传热形成阻尼效应.采用反演热传导反问题的方法,将参数辨识中的灵敏度法和共轭梯度法应用到实验,建立了一种可用于材料熔化过程界面稳态和瞬态接触热阻的辨识方法.模拟辨识与初步应用表明:辨识方法计算精度高、稳定性好,为精细化分析相变材料热控系统中的熔化吸热过程奠定基础.     

连杆位置综合的一种解析方法

本文根据半角转动的原理,提出了一种连杆位置综合的解析方法。此法概念清晰,计算比较方便。而且易于与图解法相沟通。列出了连杆五个位置综合问题的四个非线性方程。在解法上把上述问题化归为用优化方法求解。     

多刚体系统动力学的正则方程与辛算法

针对多刚体系统动力学数值计算精度的关键问题，首先采用Ｋａｎｅ方程导出了二阶形式的动力学方程，并研究了其中系数矩阵的恒等关系，进而建立了多刚体系统动力学的哈密顿体系并获得正则方程。在此基础上，分析了多刚体系统动力学数值积分的能量耗散现象，揭示其根源在于，数值算法用前一时刻值代替当前值。理论分析和计算结果表明，采用多刚体系统动力学的正则方程与辛算法的结合，可使多体系统动力学的数值计算性态得以改善，能有效消除“能量耗散”现象，从而使动力学计算的误差累积与违约得以控制，保证了计算精度。     

Euler方程无网格算法及布点技术

研究了Bitina的显式无网格算法。在该算法的基础上,给出了Euler方程无网格离散形式,运用RungeKutta显式时间推进格式推进求解。文中的耗散项采用了与非结构网格上类似的守恒型耗散算子,边界条件的处理借鉴了结构化网格处理技术。此外,本文还描述了一种区域离散布点方法,研究了点云生成的选点准则,并成功地数值模拟了二维翼型的典型绕流。     

强非线性问题的新解法

用改进了的LP法求解一类强非线性自由振动的初值问题及一类非线性二阶常微分方程的初值问题（未知函数取值较小时）。算例表明,本文方法对于求解强非线性问题,具有很好的精度,而通用的奇异摄动法（如多尺度法）或M．E．Shvez法则差多了。     

高精度微分求积曲梁单元的建立与应用

首先由能量原理导出曲梁弯曲问题的控制微分方程，在此基础上应用微分求积法原理分别给出了曲梁内点和端点的微分求积方程，由此形成曲梁单元的刚度方程，从而建立了微分求积曲梁单元，并给出了曲梁结构刚度方程的边界条件。通过算例分析，得到了微分求积单元法结构离散时应使单元数量少的原则和求解精度与单元长度比基本无关的性质。与有限元方法的结果比较表明，本文导出的曲梁单元是一种具有很高求解精度的单元。》     

High-Order Discontinuous Galerkin Method for Hovering Rotor Simulations Based on a Rotating Reference Frame

An implicit higher ? order discontinuous Galerkin(DG) spatial discretization for the compressible Euler equations in a rotating frame of reference is presented and applied to a rotor in hover using hexahedral grids. Instead of auxiliary methods like grid adaptation,higher ? order simulations(fourth ? and fifth ? order accuracy) are adopted.Rigorous numerical experiments are carefully designed,conducted and analyzed. The results show generally excellent consistence with references and vigorously demonstrate the higher?order DG method's better performance in loading distribution computations and tip vortex capturing, with much fewer degrees of freedom(DoF). Detailed investigations on the outer boundary conditions for hovering rotors are presented as well. A simple but effective speed smooth procedure is developed specially for the DG method. Further results reveal that the rarely used pressure restriction for outlet speed has a considerable advantage over the extensively adopted vertical speed restriction.     

COMPUTATION OF UNSTEADY AERODYNAMIC FORCES ON WINGS WITH NAVIER-STOCKS EQUATIONS

INTRODUCTIONFor the requirement of instantaneity of un-steady flows,the minimal time step is needed,which leads to a requirement for a large numberof time steps.To an explicit method,it marchesvery slowly and many of the acceleration tech-niques used for steady flows cannot be used be-cause they destroy the time accuracy.The“dual-time” method allows an implicit discretisation tobe used in real time,but at each real time stepmarches the solution in a pseudo time to steadystate,through an …     

大口径高精度复合材料反射器制造技术

为满足星载碳纤维蜂窝夹层结构天线反射器的大口径、高精度要求,提出了一种基于设计、工艺和测量的高精度反射器制造方法。该方法通过高热稳定性、高精度模具来提高反射面形面精度,通过蜂窝拼接来降低反射面固化时残余应力,通过高低温时效来降低固化后反射器部件残余应力,通过优化测量方法来提高反射器测量精度。实测结果表明,该技术能有效提高反射器制造精度和改善反射器热变形。     

飞机液压系统柱塞泵热分析

随着飞机液压系统向高压化、大功率化方向发展,液压系统的发热及温升问题受到了广泛关注。柱塞泵作为液压系统中重要的能量转换装置,对液压系统的功率损失和温度有重大的影响,良好的柱塞泵的热分析对液压系统的热设计和温度控制有重要的意义。然而传统的柱塞泵热分析方法,如平均油温法、软件仿真法和热节点网络法,因不够准确、高效和简洁而不能满足现代飞机液压系统热分析和热设计的相关需求。本文针对飞机液压系统柱塞泵发热的相关问题,采用基于集总参数的热网络分析方法,建立了柱塞泵热网络模型,通过对柱塞泵不同工况下的仿真,验证了模型的有效性和准确性,为柱塞泵及液压系统的热分析提供了参考。     

求解代数Riccati方程的一种新算法

本文讨论了Hamiltonian矩阵在辛相似变换下的标准形,由此给出代数Riccati方程存在非负定解的一个充分条件,提出了求解代数Riccati方程的一种新的算法。该方法节省运算量,且精度较高,尤其适合于解阶数不太高、系数矩阵为满阵的Riccati方程。最后给出了一个数值例子,并将该方法与其他方法作了比较.     

正交各向异性圆柱层壳振动和屈曲分析

本文从正交各向异性圆柱层壳的动力和相容方程出发，导出了层壳振动和屈曲的控制微分方程，再用加权残值法进行分析，得支了层壳的振动和屈曲的一般表达式，由此可以求得层壳的自振频率和临界载荷，文中引用算例说明应用，并通过３种不同圆管自振频率的比较，验证了本文所介绍方法的有效性。     

CALCULATION OF THREE DIMENSIONAL BOUNDARY LAYER EQUATIONS USING DIFFERENTIAL QUADRATURE METHOD

ＣＡＬＣＵＬＡＴＩＯＮＯＦＴＨＲＥＥＤＩＭＥＮＳＩＯＮＡＬＢＯＵＮＤＡＲＹＬＡＹＥＲＥＱＵＡＴＩＯＮＳＵＳＩＮＧＤＩＦＦＥＲＥＮＴＩＡＬＱＵＡＤＲＡＴＵＲＥ　ＭＥＴＨＯＤＴａｎｇＤｅｎｇｂｉｎ；ＣｈｅｎｇＧｕｏｗｅｉ（ＤｅｐａｒｔｍｅｎｔｏｆＡｅｒｏ...     

新型无铰式旋翼气弹综合分析的全本征方程方法

针对新型无铰式旋翼构型存在的大变形几何非线性及根部多路传力的特点,并考虑非定常气动载荷作用,建立了无铰式旋翼气弹综合分析的全本征方程。发展了高精度的伽辽金法变阶有限单元,利用Newmark平均速度法和牛顿松弛迭代,实现了旋翼动力学与气动力的紧耦合求解。经过算例的对比验证表明,本文方法能够准确模拟新型无铰式旋翼根部多路传力条件,并在稳态气弹响应、振动载荷以及气弹稳定性计算等方面具有较好的精度。     

微分求积单元法分析变曲率井内受径向约束管柱的非线性屈曲

基于Love的空间弯扭杆平衡方程,通过引入井壁约束条件导出了平面变曲率井内受径向约束管柱的平衡方程。采用微分求积单元和增量迭代法求解了曲率线性变化的井内管柱的非线性屈曲问题,通过与有限元计算结果的对比验证了所构建方法和编写程序的正确性,而且还表明微分求积单元法有方法简单、易于实施,计算量少、精度较高等优点。计算结果表明,等曲率井内管柱屈曲的临界载荷明显大于曲率线性变化井内管柱屈曲的临界载荷,另外,变曲率井眼的曲率变化对管柱弯矩、井壁约束力有显著的影响。     

求解织物热湿耦合方程的控制体-时域递归展开算法

为求解织物热湿传递耦合方程,给出了一种基于控制体积法和时域递归展开的求解织物热湿传递耦合方程的算法。首先,在时域对方程变量和参数进行级数展开,然后使用控制体积法对方程空间域进行离散,从而化连续性的非线性微分方程组为一系列的递归形式的线性代数方程组。文中给出了问题的求解步骤和算例。结论表明该方法的预测结果具有不依赖时间步长的特点。