学科间载荷参数空间插值传递方法

针对多学科优化分析中面面耦合学科网格间载荷传递的问题, 提出了在参数空间中的三次函数插值传递方法.该方法将三维空间的源学科载荷点投影到二维参数空间平面上, 用三次插值函数对目标点载荷进行计算.此方法解决了叶片表面弯曲对载荷传递精度的影响.某涡轮叶片温度传递的结果表明本文的方法具有很高的精度.      

基于网格变形技术的涡轮叶片变形传递

为实现学科耦合,需要在涡轮叶片多学科设计优化的气动模型和结构模型之间传递温度、气压和变形等载荷.研究结构变形向气动网格传递的实现方法.参数空间插值方法解决耦合面网格不一致时的变形插值问题;网格变形技术用来调整气动网格内部节点的位置,保证变形传递后的网格质量.某涡轮叶片变形传递实例的表明:该方法能够用于多学科设计优化问题中的变形传递.      

弹性载荷设计中插值方法的研究

通过对考虑静气动弹性影响的载荷设计中多种插值方法的研究与比较,对适用于三维情况下的气动和结构模型之间耦合计算时的载荷传递和位移插值方法进行原理分析和优劣比较,最终确定了曲面样条函数和形函数面积坐标加权两种数据传递方法,分别对其进行了公式推导与编程实现,并利用已有的算例模型计算得出的各项气动力数据调用程序进行插值计算,最后将插值前后的结果进行对比来验证程序的合理性。     

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

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

两种高精度高效率静气弹数据传递函数

遵循静气弹问题中气动/结构网格之间的数据传递要求,引入了基于径向基函数的流固耦合界面数据传递方法,确定了紧支半径的选取原则,并编制相应的计算程序。以M6机翼为例,将不同径向基函数的插值精度及计算效率进行比较,结果表明体样条（Volume Spline）函数和C2函数具有良好的插值精度和计算效率。这两种流固耦合界面数据传递方法非常适用于工程应用。     

基于MDO体系的涡轮叶片热-结构耦合分析

在涡轮叶片参数化建模和气动分析基础上,开发了三维坐标插值程序,实现了学科间载荷信息传递;基于经验准则公式,开发了换热系数计算程序;在研究ANSYS软件参数化设计语言的基础上,利用插值及换热程序输出APDL参数化加载宏文件,实现了边界条件的精确加载;基于ANSYS软件APDL命令流,设计了涡轮叶片热分析模块及热-结构耦合分析模块,为建立涡轮叶片MDO体系奠定了基础。     

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

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

基于多重网格方法的跨声速颤振数值模拟研究

使用多重网格的方法对二维翼型和三维机翼的跨声速颤振进行了数值求解.流场控制方程为N-S方程,湍流模型采用SA模型网格.计算网格采用结构网格,空间离散采用中心格式,使用双时间法进行时间推进.流场与结构之间的数据通过径向基函数(RBF)插值方法进行交换.通过耦合求解流场控制方程和结构运动方程得到Isogai(case A)机翼模型和Agard445.6机翼的跨声速颤振边界.通过与使用单重网格时得到的结果相比较得出,多重网格方法能够节省计算时间,提高颤振分析计算效率.     

基于迁移学习的SAR目标超分辨重建

针对合成孔径雷达(SAR)目标超分辨重建问题,提出了一种基于迁移学习的超分辨方法。在光学图像梯度域中联合训练超完备字典与稀疏编码映射,利用半耦合字典联系SAR图像与光学图像,寻找SAR图像在半耦合字典下的稀疏编码,并在高分辨率字典下完成重建。结合SAR图像的先验信息,使用正则化方法对SAR目标进行特征增强。所提方法在TerraSAR-X数据和MSTAR数据上进行了仿真实验,重建结果表明,相比目前的插值方法和稀疏表示方法,所提方法空间分辨率可提高0.5~1.5个像素。正则化增强结果表明,引入稀疏先验的正则化增强能够进一步提高空间分辨率并抑制杂波比,最后分析了正则化参数的选取对图像质量的影响。     

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

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

流固耦合数据交换的插值精度影响因素研究

提出了一种投影-形函数插值的新方法,用发展了三维线性插值方法。又包括使用优化插值和常体积转化等6种插值方法,从表征多种不同网格质量和待插值函数值情况的角度,进行了插值精度验证的应用研究,探索了网格质量和待插值函数值对插值精度的影响程度。研究结果表明,投影-形函数插值方法的效果很好;流/固的网格质量及其匹配状况在很多情况下对插值精度的影响要远大于插值方法本身。     

曲面样条插值法在飞机测压试验中的应用

飞机载荷设计中会涉及求解大量不同攻角和襟翼偏度的机翼气动压力分布,而风洞测压试验所测得的压力分布是有限的,如何选择一种合适的插值算法求出载荷设计中所需攻角和偏度的压力分布,是飞机气动载荷设计中的重要环节。采用样条曲面系数插值法,通过构造机翼压力分布的样条曲面函数,得到两种攻角和偏度情况下的压力曲面,进而插值出在两种攻角和偏度之间的任意一种攻角和偏度的压力曲面。将插值得到的压力数值和试验数据进行相似性对比分析发现:插值和试验值的相似系数较高与试验数据吻合度好。     

离心式压气机流-热-固耦合分析

利用松散耦合方法建立离心式压气机稳态流-热-固耦合模型,将离心式压气机解耦成气动、传热和结构三个子系统,子系统间存在热、力和变形的传递。由于不匹配的网格划分,利用反距离平均方法构造插值函数实现热与力向结构的传递,采用统一有限元模型实现热-固顺序耦合分析,利用网格重生成技术实现结构变形到气动分析模型的传递。经过3次迭代离心式压气机耦合系统达到收敛。结果表明,分析的离心式压气机内部流-热-固耦合强烈,耦合解法与非耦合解法得到压气机响应存在较大差异,对于多学科耦合强烈的构件,需要考虑学科间的耦合效应反应真实的工作状况。     

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

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

复杂载荷下偏心穿透裂纹应力强度因子的计算方法

提出了三维矩形平板内偏心穿透裂纹承受复杂非线性载荷下,仅含有3个系数的通用权函数。采用有限元法获得了三维矩形平板裂纹体模型下的3组参考应力强度因子,结合二元拉格朗日插值法获得了通用权函数系数。进行了自洽性验证以及裂纹面上承受简单3次、5次、7次、复杂幂函数应力分布和残余应力分布下通用权函数的验证。结果表明,通用权函数的自洽性误差在0.9%以内;裂纹面上承受复杂幂函数应力分布时,通用权函数法的误差在4.5%以内;裂纹面承受残余应力分布载荷时,通用权函数法的误差在8.5%以内。说明所提出的通用权函数具有较高的计算精度,可以满足工程中高效、准确地计算偏心穿透裂纹承受任意复杂非线性载荷下的应力强度因子需求。      

基于径向基函数与混合背景网格的动态网格变形方法

发展了一种基于径向基函数(RBF)插值和混合背景网格映射的动态网格变形方法。混合背景网格由运动物面附近的各向异性网格和偏离运动物面的各向同性网格组成。物面变形或运动后,包含在各向异性背景网格中的所有物面网格点的位移由指定的物面运动规律得到,而各向同性背景网格点的位移通过径向基函数插值获得。然后变形后的计算网格利用计算网格与混合背景网格之间固定的映射关系代数插值得到。各向同性网格距离物面一定的距离,消除了物面保形对径向基函数插值的要求,从而可以减少径向基函数基点的数量,以提高插值的效率和可靠性。最后,通过二维和三维算例对该方法的变形网格质量、参数影响、网格拓扑依赖性、复杂外形变形能力及网格变形效率进行了测试。测试结果表明,该方法能够高效生成大变形下高质量的变形网格。     

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

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