基于浸入边界法的低雷诺数流固耦合数值模拟(英文)

提出一种基于SIMPLE算法的非定常流固耦合计算方法。流体Navier-Stokes方程空间采用非结构化网格有限体积法离散,时间项采用了欧拉隐式方法。利用浸入式边界方法模拟静止或者运动固体区域,流固界面作用力通过流体体积(VOF)方法进行处理。从而可以用固定网格求解任意复杂区域中的流固耦合作用。本文模拟了低雷诺数静止及振荡圆柱绕流,所得结果与文献中贴体网格计算结果吻合,从而验证了本文方法的合理性和正确性。     

超声速平头圆柱绕流DSMC/EPSM混合算法研究

本文将EPSM算法与DSMC方法结合，构造了可模拟含近连续流区及过渡流区的DSMC／EPSM混合算法。运用混合算法模拟马赫数等于5．37时超声速平头圆柱绕流，并与DSMC结果进行比较，证明了DSMC／EPSM混合算法的有效性，并对两种算法的计算效率进行了比较。     

谐波平衡法在非定常流场中的应用

谐波平衡法（HBM）适用于模拟周期性非定常流场。求解这类问题时,只需要一个周期中几个等距时刻的流场解,即可重建整个周期流动的时间历程,计算效率较高;同时还采用谱算子精确求解流场变量对物理时间的导数项,保障了算法的计算精度。鉴于此,详细介绍了谐波平衡法的原理和实现过程,并给出了如何在现有计算程序下实现谐波平衡法的流程。取跨声速翼型和超声速钝锥两个周期性强迫俯仰振荡绕流为例,以双时间步方法为参考,对谐波平衡法的效率和精度作了详细的考核和分析。计算结果表明,谐波平衡法的内存消耗较大,谐波数N_H=2时约为双时间步方法的4倍,但谐波平衡法取得与双时间步方法相近结果的计算时间仅为后者的1/5。因此,谐波平衡法是一种有应用前景的工程实用方法。对于长周期问题,谐波平衡法的优势更加明显。     

Static aeroelastic analysis of very flexible wings based on non-planar vortex lattice method

A rapid and efficient method for static aeroelastic analysis of a flexible slender wing when considering the structural geometric nonlinearity has been developed in this paper. A non-planar vortex lattice method herein is used to compute the non-planar aerodynamics of flexible wings with large deformation. The finite element method is introduced for structural nonlinear statics analysis. The surface spline method is used for structure/aerodynamics coupling. The static aeroelastic characteristics of the wind tunnel model of a flexible wing are studied by the nonlinear method presented, and the nonlinear method is also evaluated by comparing the results with those obtained from two other methods and the wind tunnel test. The results indicate that the traditional linear method of static aeroelastic analysis is not applicable for cases with large deformation because it produces results that are not realistic. However, the nonlinear methodology, which involves combining the structure finite element method with the non-planar vortex lattice method, could be used to solve the aeroelastic deformation with considerable accuracy, which is in fair agreement with the test results. Moreover, the nonlinear finite element method could consider complex structures. The non-planar vortex lattice method has advantages in both the computational accuracy and efficiency. Consequently, the nonlinear method presented is suitable for the rapid and efficient analysis requirements of engineering practice. It could be used in the preliminary stage and also in the detailed stage of aircraft design.     

转子强度整体求解的有限元柔度计算法

本文对既有径向过盈又有轴向预紧、结构复杂的转子的整体强度提出一种计算方法,即先用有限元素法求转子各子结构联接节点柔度,再求联接节点内力,最后求出各子结构变形与应力。其优点是使转子与各子结构间的力学关系清晰,并可在小型计算机或微机上求解大题目。本文提供了该方法的力学模型与基本方程、计算方法和算例。      

面向混合网格高精度阻力预测的熵修正方法

Roe通量差分分裂格式具有耗散小、分辨率高等特点,在亚跨超声速流场模拟中得到广泛应用,但空间离散时需要对Roe平均矩阵的特征值进行熵修正,使用熵修正会增大耗散,影响阻力的求解精度。本文针对非结构混合网格中Roe格式熵修正的特点,通过改进传统的Roe格式Harten-Yee熵修正方法,提出了一种可提高非结构混合网格黏性计算精度的Harten-Yee熵修正改进方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析,改进后的熵修正方法残差收敛特性与原始Harten-Yee熵修正一致,计算精度提高,计算结果和无熵修正时基本一致,说明改进后的熵修正方法既保留了使用熵修正带来的程序鲁棒性等优点,同时把熵修正对阻力预测精度的影响降到了最低,验证了方法的有效性。采用改进后方法对第3届AIAA阻力会议的计算模型DLR-F6翼身组合体进行了详细的模拟,分析了网格收敛性和雷诺数的影响。结果表明：改进的Harten-Yee熵修正,更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。     

旋转离心叶轮内湍流的非结构化网格数值计算

给出了同位非结构化网格数值求解旋转离心叶轮内湍流流场的SIMPLE算法。网格划分采用改进的阵面推进法。控制方程采用通用的对流扩散方程并利用有限控制容积法离散。湍流模型采用标准的k-ε模型。形成的代数方程组用带有预条件器的共轭梯度平方法（CGS）求解。最后利用所给出的算法对一实验用的旋转离心叶轮进行计算。计算结果和实验数据基本相符,说明该算法合理,编程正确,可作为流场预测的理论工具。     

高超声速平头圆柱绕流化学反应流场的DSMC/EPSM混合算法研究

本文发展了平衡粒子模拟方法(EPSM)，建立了与高温气体化学反应动力学理论相匹配的EPSM耦合模型，并通过混合参数进行流区的自动识别，将EPSM方法与蒙特卡罗直接模拟方法(DSMC)结合，构造了可模拟化学反应流动的DSMC／EPSM混合算法。应用该算法对汲及化学反应的轴对称情况下高超声速平头圆柱绕流流场进行模拟，将结果与DSMC方法的结果进行比较，验证了新算法对求解化学反应流动的可行性。将混合算法的计算效率与DSMC方法的计算效率进行比较，发现混合算法能够大大提高计算效率。     

基于弹簧质量模型的复合材料螺接修理载荷传递计算方法

复合材料螺接修理具有操作简便,对连接件表面处理的要求不高,施工快速、性能可靠等优点,在复合材料的临时修理、尤其是战伤修理中应用较为广泛。然而其修理设计过程较为复杂,难以满足工程快速定参的需要。在螺接接头弹簧质量模型的基础上,针对蒙皮和补片等宽的情况,提出了一种复合材料螺接修理结构载荷传递比例计算方法。通过引入载荷按刚度分配的原则对模型进行修正,将模型进一步推广到蒙皮宽度大于补片宽度的情况。然后讨论了模型中不同弹簧刚度的计算方法。给出了采用细观力学方法以及均匀化方法等计算含损伤蒙皮等效刚度的计算方法,推导得到了规则排列多列螺栓连接中各排螺栓等效刚度的计算方法,并证明该刚度满足叠加原理。最后,以含圆形损伤孔的复合材料蒙皮板的螺接修理问题为例,对模型进行了考核,并与有限元法（FEM）预测结果进行了对比。结果表明：模型预测结果与有限元结果吻合较好,预测误差最大为7.7%。采用该模型可以高效、准确地实现复合材料螺接修理的分析设计。     

面向非结构混合网格高精度阻力预测的梯度求解方法

针对非结构混合网格的特点,通过改进传统的Green-Gauss梯度求解方法,提出了一种可提高非结构混合网格黏性计算精度的节点型Green-Gauss梯度求解方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析。改进后的梯度求解方法残差收敛更好,下降量级更多,阻力系数和试验吻合更好,激波区域压力分布和分离区域流场细节的模拟更精确,说明改进后的梯度求解方法有效提高了程序的鲁棒性和阻力预测精度,验证了方法的有效性。采用改进后的方法对第5届AIAA阻力预测研讨会的通用研究模型（CRM）进行了详细的模拟分析,结果表明：改进的梯度求解方法更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。     

高超声速风洞变雷诺数试验技术研究

为满足高超声速飞行器气动力雷诺数效应研究需求,在CARDC的Φ1米高超声速风洞中开展了变雷诺数试验技术研究.该项试验技术是利用Φ1米高超声速风洞采用高压下吹-真空抽吸驱动运行方式、风洞运行参数范围宽的特点,通过宽范围内调节风洞运行总压而大幅改变模拟雷诺数.研究采用了单点变雷诺数试验技术和连续变雷诺数试验技术两种手段来开展高超声速飞行器气动力雷诺数效应模拟.单点变雷诺数试验是通过一系列不同雷诺数条件、不同试验车次的试验结果,获得气动特性随雷诺数的变化规律;连续变雷诺数试验时,控制风洞总压从高到低连续变化,测量获取模型处于某一姿态角条件时气动力随雷诺数的变化规律.本文介绍了变雷诺数试验的风洞开车方式、试验及数据处理方法等,并开展了某升力体飞行器和某弹头模型雷诺数效应试验研究.研究结果表明:采用单点和连续变雷诺数试验技术相结合的方式,能较为完整、准确地获得飞行器模型气动力随雷诺数的变化规律.     

多个空间高斯源信号情况下成组三维图像特征提取方法

针对功能磁共振成像数据中含有多个高斯信号源的盲源分离问题,介绍了一种成组典型相关分析方法(Group BSS-CCA)。这个方法的组分析框架参照了GIFT工具箱中的Group ICA算法,具体的典型相关分析方法应用的是Friman等人提出的BSS-CCA算法。以验证该方法的有效性为目的,设计了仿真实验;结果表明,该方法能较好识别出混合在人脑功能磁共振成像数据的2个空间高斯信号。Group BSS-CCA算法对研究人脑的功能磁共振成像数据具有较高的实用价值。     

Real time drag minimization using redundant control surfaces

A method for minimizing the drag of a wind tunnel model with multiple control surfaces is presented. The minimization is performed in the wind tunnel and measurements are performed in real time. Real time measurements introduce difficulties such as noise in the signals, hysteresis and problems with repeatability of the function evaluations. The lack of a numerical function to minimize therefore puts certain demands on the optimization method, and hence a derivative free method, often referred to as a generating set search method, is used. The proposed method is generalized to take both linear equality constraints as well as linear inequality constraints into account. The generating set search method is implemented in the wind tunnel and tests show that the drag can be decreased while satisfying the constraints.     

高升阻比翼型的设计

使用解析形状函数法和参数化翼型表示方法,将求解绕翼型流场的N-S方程解与优化方法相结合,可设计出新的高升阻比翼型.计算结果表明,设计的新翼型满足设计要求,且获得比原始翼型高得多的升阻比,翼型设计方法是实用的.参数化翼型表示法比解析形状函数法对给定的初始翼型离散数据的准确度影响更小,更适合用于新翼型的设计.     

离心叶轮全三维反问题方法设计误差稳定性研究

用流线曲率全三维反问题方法和引入滑移因子模型的准三维流线曲率反问题方法设计了叶片数不同的和增压比不同的小型高速离心叶轮.用商用计算流体力学(CFD)求解器估算了这些叶轮的性能.对比了这两种反问题方法的设计误差稳定性,考察了分流叶片对设计误差的影响.结果表明,流线曲率全三维反问题方法能够根据叶片负荷对叶片设计做出恰当调整.较准三维反问题方法,流线曲率全三维反问题方法所设计叶轮的增压比普遍偏高.随叶片负荷增加,较准三维反问题方法,流线曲率全三维反问题方法的设计误差变化较小,即其设计误差稳定性明显较好.分流叶片对准三维反问题方法设计结果的影响明显较小.      

基于Kirchhoff方法的跨音速螺旋桨的气动声学计算

本文采用Kirchhof方法计算双桨叶Hartzel跨音速螺旋桨的远场气动噪声。Kirchhoff面为包围整个螺旋桨的固定圆柱面,首先采用FLOWer软件包直接求解Euler方程得到圆柱内和圆柱面上的近场气动流场解,然后根据Kirchhof方法求解远场声场。通过计算详细研究了圆柱的直径和高度对远场声场的影响,以及圆柱各部分对远场声压的贡献,文中同时给出了Farassat线性方法的计算结果。      

跨音速叶栅粘流计算的多级广义有限元法

将 Taylor-Galerkin广义有限元法和多级有限元的思想结合起来 ,并在人工粘性的处理上作了改进 ,构成了收敛速度和稳定性均较好的多级广义有限元法。利用该方法 ,并基于 N-S方程研究了二维跨音速叶栅流动。计算结果与实验结果符合较好。通过计算表明 ,该方法计算稳定、收敛较快 ,是叶轮机械内部跨音速流动强有力的计算手段     

空间柔性网—膜结构分析的数学规划算法

 提出了以参数变分原理为基础的数学规划算法, 以计算分析太空用品中网- 膜结构起皱现象。针对网状柔性结构分析, 算法的原理是将结构化成多个具有一定非线性本构关系的杆单元, 并在此基础上进行杆系结构非线性分析的规划算法构造, 可以较好地模拟网状柔性结构的起皱性质。而对于一般的膜结构, 则构造了用于直接分析的膜单元相对应的本构模型, 并在此基础上构造问题的求解算法。不同于常规的数值计算方法, 提出的算法不依赖于应力的迭代过程。参数变分原理和数学规划算法的应用, 使本算法可以很有效地分析预测网- 膜起皱的区域、模式和应力分布, 从而为解决网- 膜结构起皱问题提供依据。给出的数值算例说明了方法正确与有效性。     

Analytical solution methods for eigenbuckling of symmetric cross-ply composite laminates

Based on the first-order shear deformation theory, this paper explores the analytical methods for eigenbuckling of symmetric cross-ply rectangular composite laminates with a pair of parallel edges simply supported and the remaining two edges arbitrarily constrained. The main con-tribution of present paper lies in two aspects: one is to present a simple and effective analytical method, namely, the separation-of-variables method, which can generate the closed-form buckling solutions without any computational difficulty;the other is to incorporate the accurate computation method of exponential matrix into the state space technique to avoid the inevitable numerically ill-conditioned problems reported in several literatures. The results obtained via both analytical meth-ods are identical, and a good agreement with their counterparts in literature is observed. The separation-of-variables method can generate exact solutions within 1 s, which is impossible if the state space method is employed. Besides, the combination of the accurate computation method of exponential matrix and the state space method greatly improves the computational efficiency and gives correct results compared with the straightforward use of state space method.     

结构系统可靠性及可靠性灵敏度分析的改进子集模拟法

在再生自适应子集模拟(RASS)法的基础上,提出了一种改进的再生自适应子集模拟(MRASS)法以用于结构系统的可靠性及可靠性灵敏度分析。MRASS法继承了RASS法中马尔可夫链再生、延迟拒绝、自适应马尔可夫过程及分量各自采样等优点,并改进了自适应采样过程。MRASS法通过对产生的样本点的接受率与最佳接受率进行比较,有针对性地寻找合适的建议分布方差,提高了抽样的效率。工程算例的数值结果表明:MRASS法相比于RASS法及传统的子集模拟(SS)法,在处理具有高维随机变量、小失效概率及高度非线性特点的结构系统时,有更好的适应性、稳健性及精度。