基于径向基函数的载荷重分配算法

提出一种基于径向基函数(Radial Basis Functions,简记为RBF)的载荷与质量重新分配算法.通过局域支撑的径向基函数,构建函数空间中某子空间的一组基函数,选择合适的RBF影响半径,从而完成基于RBF的载荷与质量重分配流程,并开展数据插值算法在载荷质量分配过程的应用研究.以工程实际需要为出发点,以细网格...     

非结构嵌套网格中的一种改进型径向基函数插值方法

非结构混合网格拓扑关系相对复杂,在嵌套边界区域进行流场信息插值传递容易有精度的损失。利用已知的解析函数作为插值样本,对线性插值(LINE)、距离倒数权重插值(WA)、逆二次径向基函数插值(IQ)、Wendland暞sC2径向基函数插值四种方法的精度进行对比,筛选出精度较高的C2径向基函数;提出一种改进的C2径向基函数插值方法,通过调整径向基函数的作用半径来控制插值矩阵条件数,进而消除使用径向基函数插值在嵌套边界区域出现的数值奇性问题;选用MD30P/30N 多段翼与三维AEDC典型外挂物分离模型对该插值方法进行验证。结果表明:改进后的径向基函数插值方法能够有效消除嵌套边界区域数值奇性,同时计算收敛速率更快,插值精度更高。     

基于径向基函数的发动机瞬时推力计算方法

为了快速准确插值计算发动机的瞬时推力,提出了使用基于径向基函数(RBF)的推力计算方法。通过算例证明,采用径向基方法的插值结果精度高,算法成熟,使用计算机编程容易实现,可以快速求解发动机瞬时推力插值问题。该方法在飞机设计、特别是性能计算及参数优化过程中可以广泛使用。     

基于响应面插值的非线性气动弹性计算

飞翼、大展弦比低雷诺数气动布局容易在小迎角的条件下出现气流分离,会带来明显的非线性气动力问题,同时气动弹性带来的影响亦不可忽略。针对此类布局提出一种建立基于径向基函数插值的非线性压力系数分布模型的方法,利用径向基函数插值建立面元上压力系数对迎角导数的响应面,将压力系数积分并通过无限板样条(IPS)方法进行气动结构多次迭代插值实现非线性气弹分析。结果验证了该方法对于静气动弹性分析的有效性,同时能准确地反映弹性带来气动效率的降低和变形对升力及阻力的影响。     

一种基于径向基函数和峰值选择法的高效网格变形技术

基于径向基函数的网格变形方法因其具有诸多优点,而被广泛应用于气动外形优化设计等领域。对于大规模网格或复杂构型,该方法所需计算量是难以承受的。为了提高网格变形效率,可以通过减少建立插值模型所需支撑点数目来实现。为此,提出一种高效的选点算法--峰值选择法。该算法在选点过程的每个迭代步中对边界节点处的误差进行分析,从物面节点中选取多个峰值点来更新支撑点集,减少迭代步数,提高选点效率。在该算法的基础上,实现了网格的高效变形。三段翼型的网格变形算例证明：该方法可以在保证网格质量的同时实现复杂网格的变形。以DLR-F6复杂模型（约1 000万网格）的刚性运动和弹性大变形为算例对该方法的变形效率和变形后网格质量做了进一步评估：当相对误差设置为5.0×10^-7时,在保证变形后网格质量的前提下,该方法变形效率最快比传统贪婪算法提高了13倍,其中在选点效率方面最快提高了31倍。     

一种高效的CFD/CSD耦合飞行器多学科优化设计方法

建立了一种基于CFD/CSD (computational fluid dynamics/computational structure dynamics)耦合的多学科优化设计方法.该方法采用结构有限元方程进行结构特性分析,同时采用数值求解N-S方程方法获得气动载荷,保证了设计结果的可信性,并采用基于代理模型的优化设计方法解决了采用高精度模型后计算量太大的问题.采用RBF(径向基函数)插值方法进行流固耦合界面数据插值,保证了传递过程中的能量守恒.利用建立的优化设计方法对一典型的翼身组合体进行了气动和结构优化设计,算例表明在满足约束条件的前提下优化后阻力和质量分别减小了5.0%和4.2%,证明了该方法的有效性.      

基于径向基函数的网格变形及非线性气动弹性时域仿真研究

为开展非线性气动弹性研究,基于非线性结构有限元软件NASTRAN和自主研制的多块结构化计算流体力学(CFD)求解器,开发了一套基于计算流体力学/计算结构动力学(CFD/CSD)耦合求解方法的气动弹性时域仿真程序.该程序采用径向基函数(RBF)交换两套求解器之间的数据并进行网格变形.为提高RBF方法的效率,构造了基于多次插值的空间待插值点精简算法.在多次插值过程中,每次插值的对象为上次插值的误差,并同时限制插值区域,以此实现了空间待插值网格数的精简.数个网格变形的算例表明该方法可支持大变形运动,并且具有较高的计算效率.采用此程序开展了AGARD 445.6机翼颤振计算、大展弦比机翼的静气动弹性计算与切尖三角翼极限环振荡(LCO)现象的动气动弹性仿真,结果揭示了当机翼展弦比较大或者响应幅值较大时,结构非线性对于气动弹性有显著影响.     

基于物理量梯度修正的RBF数据传递方法

径向基函数是复杂外形飞行器多场耦合计算中重要的数据传递方法。复杂飞行器物理量分布通常与局部流动特征密切相关且在不同区域变化剧烈,如何在具备各向同性特征的径向基函数中考虑实际物理量各向异性分布因素的影响成为制约其精度进一步提升的关键。针对径向基函数的问题,提出了一种基于局部物理量梯度对径向基三方向权重进行自适应修正的数据传递新方法。采用高超声速舵面、翼身组合体算例对该方法在单次数据传递中的可行性及效果进行了验证,同时基于该方法对压缩拐角外形开展了气动力/热/结构多场耦合数值模拟,并与风洞试验结果进行了对比分析,以验证该方法在实际耦合问题计算长时间高频次数据传递中的适用性及可靠性。研究结果表明该方法能够很好地提升如薄板样条插值（TPS）、多重二次函数插值（MQ）等全域基函数的单次插值效果,并且有效地改进了MQ方法对形态参数的鲁棒性,而对于梯度修正后的紧支基函数,其改进效果更佳,以较少的点达到了TPS、MQ等全域基函数的插值效果,实现了数据传递效率和精度的兼顾。采用该方法获得的耦合计算结果与试验对比一致,表明了该方法的可用性以及在复杂外形多场耦合问题中的良好工程应用前景。     

基于RBF神经网络的压气机特性仿真

为克服传统仿真方法误差较大的问题，提出了一种基于RBF（多变量插值的径向基函数）神经网络的压气机特性仿真方法。利用RBF神经网络能逼近任意非线性系统的特点，对压气机特性进行了拟合。试验结果表明，此方法具有精度高，收敛速度快等优点，可广泛运用于发动机数值仿真及控制模拟等领域。     

基于本征正交分解的流场快速预测方法研究

对于飞行器气动设计,通过实验或CFD计算获得所有设计状态处流场信息成本高、耗时长,难以满足工程应用要求。利用基于本征正交分解(POD)降阶模型结合相应的插值方法,只需要通过实验或CFD计算获取有限个数的采样流场,可以快速预测出具有满足精度要求的设计状态处流场信息。首先,证明对POD基系数插值获取预测流场与对采样流场直接插值获取预测流场的等价性,并分析采样流场分布的要求。然后,以三段翼型流场作为研究对象,POD降阶模型分别与三次样条插值和径向基插值结合,实现设计状态处流场信息快速预测。最后,通过将流场预测结果与相同条件下的流场CFD计算结果对比,表明了POD降阶模型结合两种插值方法的各自特性。     

基于径向基神经网络的气动热预测代理模型

为快速获取高超声速飞行器表面热流数据并缩短飞行器气动热设计周期,提出了一种基于径向基神经网络的气动热快速预测代理模型方法。首先,在飞行器表面每一个离散化的网格节点单独构造一种正则化的径向基神经网络。随后,通过训练集对所有网络同时进行训练,获得各自网络的连接权值。最后,所有网格节点的神经网络协同预测飞行器表面不同位置的热流。对NASA火星实验室的椭圆钝化高超声速飞行器的应用表明,所提出的代理模型方法在模型训练完成后能够快速进行飞行器表面热流预测,并且模型具有良好的泛化能力,在驻点及迎风大面积区域热流预测结果与数值模拟的偏差在10%以内。     

典型尖楔结构燃气流热试验工况确定方法

在传统冷壁热流模拟方法的基础上,进一步提出以热壁温度及热流密度的时序变化曲线为控制目标的燃气流热试验工况确定方法,即利用壁温控制目标与实测值的偏差对热壁热流控制目标做一定修正,以尽可能消除和弥补前期试验误差,同时利用300K冷壁边界热流密度数据库插值迭代方法,快速确定一定气动热模拟所需燃气流温度,解决了沿飞行轨迹瞬态热试验技术难题之一。利用CFD数值模拟方法,建立了典型尖楔结构高/中温双路燃气流组合热试验300K冷壁边界热流密度数据库,并针对典型尖楔结构沿某飞行轨迹9个典型状态气动热模拟需求,确定相应双路燃气流热模拟参数。相关数值计算结果显示,驻点区域热流密度平均模拟偏差为4.5%,平板区热流密度平均模拟偏差为4.6%,两者最大模拟偏差均不大于8%,满足工程试验精度要求。同时,瞬态热分析结果显示第45s时,距驻点1mm处最大温度梯度达到21K/mm,距驻点10.1mm处最大温度梯度达到18K/mm,满足气动热大温度梯度效应需求。      

基于弹簧-悬臂梁模型最小变形能的气动载荷分配方法

在进行气动载荷到结构载荷转换时使用弹簧-悬臂梁模型来代替悬臂梁模型,解决了如果气动载荷点与部分结构节点重合则会导致难以求解的问题,载荷分配结果符合“近大远小”的规律,且分配前后完全符合等效性。研究了模型参数对载荷分配结果的影响,讨论了利用悬臂梁与弹簧刚度比,刚度与目标点距离的比来控制分配结果的方法。另外,讨论了使用构造函数代替物理模型的方法。     

连分式扩充的粒子群神经网络压气机特性重构方法

对航空发动机压气机原始二维等转速线的数据点进行连分式扩充,通过两次网络训练,增加转速特性数据,在三维空间中进行BP(back propagation)网络模型重构.根据压气机特性数据空间分布的特点,引入压力比函数,调整计算区域,定义网络的输入输出数据,利用试探法确定隐含层维数.采用基于趋利避害原则的粒子群算法对网络的初始权值和阈值进行优化,建立了压气机压比和效率特性的整体代理模型.最后以某型发动机的低压压气机为例进行了压气机特性模型的重构.通过模型的校核与验证表明:采用这种方法建立的模型精度较高,优于普遍采用的传统二维插值方法和普通BP神经网络模型.最终建立的重构模型对于采用选配法、坐标法和部件法等以压气机通用特性曲线为基础的发动机模型的求解,可提高计算精度和迭代速度,具有一定的工程应用价值.      

An adaptive sequential experiment design method for model validation

Efficient experiment design is of great significance for the validation of simulation model with high nonlinearity and large input space. Excessive validation experiment raises the cost while insufficient test increases the risks of accepting an invalid model. In this paper, an adaptive sequential experiment design method combining global exploration criterion and local exploitation criterion is proposed. The exploration criterion utilizes discrepancy metric to improve the space-filling property of the design points while the exploitation criterion employs the leave one out error to discover informative points. To avoid the clustering of samples in the local region, an adaptive weight updating approach is provided to maintain the balance between exploration and exploitation. Besides, the credibility distribution function characterizing the relationship between the input and result credibility is introduced to support the model validation experiment design. Finally, six benchmark problems and an engineering case are applied to examine the performance of the proposed method. The experiments indicate that the proposed method achieves satisfactory performance for function approximation in accuracy and convergence.     

基于Hankel阵的荧光油膜灰度与厚度预测模型改进

针对荧光油膜灰度与厚度关系数据采集方法较为繁琐、耗时耗力这一问题,研究了基于Hankel阵的系统辨识算法,并在此基础提出了Hankel阵误差修正模型和Hankel阵高阶迭代误差修正模型等两种改进方法,利用了极少数据量建立模型,实现了对其余未知荧光油膜厚度值的预测,且保持了较高的精度。试验结果表明:基于极少数据量建模以预测该数据量外的数据点这一特殊背景下,插值法的外插能力显得并不适用。而传统Hankel阵预测模型的预测精度为76.69%,Hankel阵误差修正模型和Hankel阵高阶迭代误差修正模型的预测精度分别为85.69%和89.25%,较之传统方法预测精度分别提高了9%和12.56%,为荧光油流技术领域针对荧光油膜灰度与厚度建模问题提供了一种可行技术路线,具有一定的实际工程应用意义。      

基于径向基函数改进的Delaunay图映射动网格方法

为了提高大变形下动网格生成的效率和质量,基于径向基函数插值技术对原始Delaunay图映射动网格方法进行了改进。首先通过带多项式的径向基函数插值方法计算出背景网格远场边界点的位移;然后将背景网格位移插值到计算网格;最后利用衰减函数将计算网格远场位移衰减为零。基于矩形旋转网格变形实例比较了改进方法与原始Delaunay图映射动网格方法之间的差别,并研究了径向基基点数目和衰减函数类型对变形网格质量的影响。矩形旋转网格变形实例说明该方法能够精确恢复出Delaunay背景网格的旋转特性。NACA 0012翼型、NLR 7301两段翼和M6机翼网格变形算例进一步证明,通过添加控制点,该方法能够不重构背景网格实现大变形下高质量动网格的生成。