PSO节点寻优的样条逼近微分

样条节点分布对逼近精度的影响很大,样条节点寻优模型的雅可比矩阵求解复杂。利用经典参数的PSO（粒子群算法）解算样条最优节点分布模型,以逼近残差的平方和为目标函数,每步对节点序列排序后再计算目标函数,可以提高样条逼近精度,为样条节点寻优提供了一种较好的实现方法。仿真计算表明,在一定的节点个数情况下,PSO节点寻优的逼近效果比Carl de Boor的NEWNOT程序中的方法要好。     

双丝热线探头的一种标定新方法

对X型双丝热线探头标定工作进行了研究,提出使用一种新的基于B样条递推最小二乘的标定方法,该方法选取具有低阶光滑特性的B样条函数并采用递推最小二乘的方法估计控制参数,提高了校准精度。通过分析实际的风洞试验数据,对提出的标定方法的有效性进行了验证。     

弹道跟踪数据融合处理的快速算法

针对弹道跟踪数据融合处理中的大计算量环节研究了快速算法。用样条函数表示弹道参数,建立了多测元的联合观测模型和弹道参数的非线性融合计算模型,给出了弹道参数的求解算法,分析了弹道参数融合计算中的大型矩阵运算问题,利用基础线性代数函数库提高了大型矩阵的运算速度。建立了样条模型计算的非线性约束优化模型,给出了确定样条节点位置的优化算法,通过分析样条模型的计算原理设计了并行算法,实现了样条模型的并行化计算。仿真结果表明,弹道参数融合计算和样条模型计算的效率都得到了显著提高,计算时间减少了65.47%,对缩短数据处理周期有重要意义。     

逆向工程中曲面重构算法研究与实现

给出了基于坐标测量机测得的海量数据重新构造样板物体的数字化表面模型的算法流程。该算法首先利用部分测量数据构造张量积的B样条网格曲面,通过确定适当的检查点,并计算检查点到所构造的B样条网格曲面间的距离,从而控制重构曲面的构造精度。讨论了与文中算法相关的空间点到NURBS曲面间的距离计算等问题。给出了1个工程应用实例,表明该算法稳定可靠、效率较高。     

不同插值方法对残缺环阵声源定位的影响

研究了利用不完整麦克风阵列进行声源定位的算法，通过修正傅里叶插值、B样条插值和三次样条插值对残缺麦克风阵列的交叉互谱矩阵的插值补偿，获得了声源定位算法。通过数值模拟对声源定位算法进行了验证，发现在波束形成的声源强度上，傅里叶插值模拟结果最大偏差达到5.21 dB,B样条插值为1.17 dB，三次样条插值为0.80 dB；在声源位置偏差上，傅里叶插值为0.04 m,B样条插值为0.01 m，三次样条插值为0.01 m。结果表明傅里叶插值计算得到的声源定位的强度、位置精度和动态性能最差，三次样条插值最优，B样条插值表现一般。实验验证也得出类似的结论，因此利用三次样条插值计算不完整麦克风阵列的交叉互谱矩阵最佳。     

B-样条函数的一种快速递推法

 本文提出了常用的等距结点B-样条函数的一种快速递推法及递推计算格式和程序设计。阐明了此法具有计算简便、递推快速和更为直观的特点。它为B-样条函数的数学分析和推广应用提供了方便。此法将易于为一般工程技术人员掌握和运用。     

基于队列的参数域T网格Bézier分割方法

针对T样条曲面Bézier提取中的参数域T网格的Bézier分片问题,提出了一种基于队列的参数域T网格的Bézier分片算法。该算法首先通过延伸T型参数点对参数域网格进行切分,之后通过计算每个参数点所耦合的T样条混合函数非零区域对参数域网格进一步切分。上述切分过程只需借助一个队列数组即可实现。通过实例测试,验证了本算法的有效性。     

基于单纯B样条的航空发动机机载稳态模型研究

为了减小航空发动机稳态建模的模型误差、降低复杂度及提升其实时性,提出了一种基于单纯B样条函数的航空发动机稳态模型建模方法。该函数是局部多项式基函数的线性组合,因此求解该函数为线性回归问题,通过运用广义最小二乘方法来求解B系数,从而提高计算效率和提高模型精度。最后建立了基于该算法的二维和四维涡扇发动机稳态模型,并分别与相同建模样本条件下的多输入多输出约简迭代最小二乘支持向量机稳态模型进行了比较,表明了单纯B样条建模方法不仅继承了B样条的算法复杂度低、存储数据量小和实时性好等优点,同时避免了最小二乘支持向量回归机不能拟合大样本数据的缺点,且拟合效果优于最小二乘支持向量机。     

空心涡轮叶片复杂陶芯弯扭变形分析方法比较

陶芯作为熔模铸造空心涡轮叶片的内腔转接件,其弯扭变形程度直接关系到空心涡轮叶片的壁厚尺寸精度。为了研究陶芯制造过程中出现的弯曲和扭曲变形,首先,借鉴叶片给出了陶芯弯曲变形和扭曲变形的定义,并根据陶芯结构特点研究了陶芯弯扭变形的计算方法,包括基于轮廓线的弯扭变形计算方法:主动轮廓模型法和样条拟合轮廓线法,与基于数据点的弯扭变形计算方法:凸包算法、二维配准算法和距离权值算法。然后,分别通过仿真与实测数据比较了5种算法,结果表明:样条拟合轮廓线法具有更高的稳定性和精度。     

基于B样条模型的参数化卡尔曼滤波器设计

以标准B样条函数为基础,建立了以位置、样条系数为状态变量的参数化卡尔曼滤波器,用于解决外测数据的实时滤波问题。按照时间更新跨节点与否,分2种情况给出了状态转移方程。在时间更新跨样条节点时,使用样条函数的一阶连续导数条件,估计新增样条节点系数,由此实现滤波器在跨节点处的平滑过渡。通过仿真数据对新方法进行验证,并与已有的2类典型滤波方法进行比较,结果表明,本方法的滤波精度与另一类直接基于弹道信号表示的样条递推滤波方法精度相当,且可表现出更优的收敛性。新方法具有样条参数化模型的相同优点,可对时域信号全时段建模,可利用先验信息设计弹道优选节点而实现滤波性能优化,缺点在于状态更新的策略较为复杂。     

超声速翼型气动优化设计

首先分析了几何外形和相对厚度对超声速翼型气动特性的影响。基于遗传算法(GA)和气动力快速工程算法,对于相对厚度为3.5%的多边形翼型进行优化设计,多边形翼型的优化外形趋于四边形,最大厚度点后移到翼型弦线的60%左右,随着迎角或者马赫数增大下翼面会变薄,上翼面变厚,最大厚度点相应稍有后移。对于相对厚度为4%的双圆弧翼型,采用两步优化设计方法,第1步优化结合基于B样条的类别形状函数变换(CST)参数化方法与小波分解方法,实现几何外形的局部控制与光顺处理,并且采用本征正交分解(POD)代理模型降低优化过程中流场计算的工作量;第2步优化采用基于Navier-Stokes方程的最速下降法(SDA),修正第1步优化中代理模型和小波光顺引入的误差;优化设计得到的翼型近似为四边形,其相对厚度最大点后移到翼型弦线的60%~65%处,升阻比可以提高7%。     

基于解析法的用户差分距离误差解算方法(英文)

为了进一步增强星基导航系统的完好性性能,提出了一种基于解析法的用户差分距离误差计算方法。文中首先介绍了卫星时钟和星历改正数的计算方法,认为计算用户差分距离误差的关键是查找卫星服务区域内的最差用户位置。通过变换卫星时钟和星历误差协方差矩阵,将查找最差用户位置变换成一个解析几何问题,通过数学推导得出最差用户位置的解析式,进而求解用户差分距离误差。通过计算机仿真对比了解析法和遍历法之间的性能差异,结果表明:解析法的正确性由遍历法得到了验证,解析法可以减少90%的运行时间并且具有更低的计算复杂度,便于工程实现。     

基于闭路制导的弹道导弹发射诸元的仿真算法

提出了一种基于闭路制导的弹道导弹的发射诸元算法.给定发射点与瞄准点的地理位置,以压低弹道为目标,通过一种迭代方法求解关机点弹道倾角,并基于地球引力级数展开模型提出了一种虚拟目标的快速计算方法,以主动段弹道平面和被动段弹道平面二面角为指标完成了发射方位角的迭代计算,最后通过仿真算例论证了该算法的可行性.     

Polynomial Representation of NURBS and Its Application to High Frequency Scattering Prediction

This article presents a method that uses physical optics (PO) techniques to compute the monostatic radar cross section (RCS) of electrically large conducting objects modeled by non-uniform rational B-spline (NURBS) surfaces. At the beginning, a new algorithm to convert recursive B-spline basis function into piecewise polynomials in power form is presented. Then, algorithm computes the polynomial representation of B-spline basis functions and NURBS surface geometric parameters are obtained. The PO integral over NURBS surfaces of an electrically large conducting object is used to predict the object's RCS. The NURBS surface is divided into small piecewise polynomial parametric patches by isoparametric curves, and the PO integral expression over the parametric domain of each polynomial parametric patch is reduced to an analytical expression which permits an accurate and effective computation of the PO integral by using a modified Ludwig's algorithm. The RCS of the object can be obtained by adding up the PO integral contribution of each polynomial parametric patch. The effectiveness of this method is verified by numerical examples.     

k/N（G）结构系统初始备件配置方法

针对现有备件冗余度建模中无限呼叫总体假设不合理、供应可用度计算困难的问题,提出一种k/N（G）结构系统初始备件配置方法。首先,采用有限源排队论获取备件短缺数概率分布。然后,给出计算k/N（G）结构系统供应可用度的方案,并证明了方案的正确性。该方案在计算随机挑选的一套系统未因备件短缺而停机的概率基础上,利用k（或更多）套系统良好的二项分布概率估计系统供应可用度。最后,应用边际分析算法确定各项备件的数量。实例分析表明：该方法能优化k/N（G）结构系统各组成部件的库存数量,且与单项法相比具有明显的优势,能以较少的备件费用达到相同的供应可用度目标。     

星点图像的目标提取算法研究

对于星点灰度重心法等亚像素提取算法来说,需要提取星点区域内的每一个像素。基于扫描的边界搜索法只能获取目标的边界像素,文中对其进行了扩展,通过判断边界像素中最大最小行列数来确定目标区域的具体范围,从而提取目标区域内的每一个像素。针对大视场星点图像目标范围小、相互间距离远的特点,为减小计算量,提出了基于视框搜索的目标提取算法。通过实验证明,基于视框的算法在计算量上比基于边界搜索的算法减少3/4。     

涡轮叶片流固耦合协作优化方法

针对结构和流场相互影响的流固耦合多学科设计问题,提出一种耦合量曲面拟合的协作优化方法。其特点是用曲面函数拟合结构-流场界面耦合量(温度、气压和变形)的空间分布,并用曲面拟合参数的一致性代替原有耦合量界面分布的学科一致性。涡轮叶片的算例表明:曲面拟合函数可以表征结构和流场界面的耦合量;用拟合参数的一致性代替界面直接耦合是可行的;该协作优化方法是可以用来解决流固耦合多学科优化设计问题的。     

A combined method for calculating the trajectories of suspended particles

A combined method for calculating the trajectories of suspended particles in the hydrodynamic velocity fields is proposed. The method utilizes capabilities of existing CFD codes and parallel computing on a computer with a graphics card that supports CUDA technology. The software package that implements the method is described.     

基于偏导数的全局灵敏度指标的高效求解方法

在全局灵敏度分析领域,基于偏导数的全局灵敏度指标由于其优良的特性得到了广泛的关注。针对目前求解基于偏导数的全局灵敏度指标计算效率低的问题,提出了一种高效的求解方法。该方法利用乘法降维近似地将响应函数展开为连乘积的形式,从而将求解偏导数全局灵敏度的高维积分问题转化为一维积分的连乘积,在保证计算精度的前提下大大降低了求解基于偏导数的全局灵敏度指标的计算量。在求解特定点的偏导数时,采用了复数步长方法,提高了计算精度。最后,通过算例验证了所提方法的准确性和高效性。     

Dynamic modeling and analysis of vortex filament motion using a novel curve-fitting method

Applications of a novel curve-fitting technique are presented to efficiently predict the motion of the vortex filament, which is trailed from a rigid body such as wings and rotors. The governing equations of the motion, when a Lagrangian approach with the present curve-fitting method is applied, can be transformed into an easily solvable form of the system of nonlinear ordinary differential equations. The applicability of Be′zier curves, B-spline, and Lagrange interpolating polynomials is investigated. Local Lagrange interpolating polynomials with a shift operator are proposed as the best selection for applications, since it provides superior system characteristics with minimum computing time, compared to other methods. In addition, the Gauss quadrature formula with local refinement strategy has been developed for an accurate prediction of the induced velocity computed with the line integration of the Biot–Savart law. Rotary-wing problems including a vortex ring problem are analyzed to show the efficiency, accuracy, and flexibility in the applications of the proposed method.