矩的显式积分算法研究与应用

矩的求解通常被用于求解有限元、体积、惯性矩等问题中.基于矩的叠加性,首先给出了在三维空间中计算域的离散方式,并推导了矩的显式积分公式,随后将其推广到n维空间中,该表达式易于在计算机上实现;设计了矩的并行计算算法,并通过Fortran和Python混编的方式,实现了矩的并行计算;对多重精度下的样例数据给出了一个算例,实现了零阶矩和二阶矩的计算,并和串行算法、逐次降维算法作出比较,进行了效率分析和误差分析.结果显示,矩的显式积分并行计算算法易于程序实现,并且在效率上高于串行算法,能够很容易推广到高维空间,该算法具有高度可并行性,误差主要来自计算域离散.      

接收模式下加罩天线系统特性的快速电磁仿真

针对接收模式下电大尺寸的天线-天线罩系统电磁特性分析问题,提出了一种全波法和高频法结合的混合算法:积分方程/修正型表面积分+多层快速多极子(IE/MSI+MLFMA)方法.该方法将由高频方法MSI得到的天线罩内场分布作为天线阵列的激励场,再利用基于体-面混合积分方程(VSIE)和MLFMA的快速全波分析方法精确计算天线阵;同时应用球谐函数展开、预处理技术、混合并行等技术进一步改进算法的计算效率.同传统的全波数值法相比,该方法在保证计算精度、满足工程需要的前提下,解决了求解计算时间长、计算效率低的问题,实现了对电大尺寸天线-天线罩系统的快速、高效仿真.      

瞬态传热问题的微分求积和精细积分求解方法

给出了瞬态传热问题的高效高精度求解方法,该方法分别用微分求积法(DQM)和精细积分法(PIM)离散空间域和时间域.微分求积方法除了精度高、效率高之外,处理复杂边界条件的灵活性也优于有限元法(FEM).用精细积分法处理一阶瞬态传热微分控制方程,不需要增加额外自由度,还可以达到计算机精度.给出了隔热结构4种边界条件下的数值结果.并就上表面恒温、其他面绝热边界条件计算结果与有限元分析结果进行了对比,算例分析表明,采用微分求积和精细积分法布置少量的节点就可以达到较高的精度.      

一种新的矩独立重要性测度分析方法及高效算法

为了更加合理地分析输入随机变量对结构系统失效概率的影响,提出了一种新的矩独立重要性测度分析方法。传统的重要性测度指标只能估计输入随机变量在固定点时对结构系统输出响应的影响,而所提新指标能够充分反映输入随机变量在其分布区域的所有缩减区间上变化时对结构系统输出响应的平均影响程度,更加符合工程实际。为了求解新指标,给出了2种算法:传统的双层重复抽样蒙特卡罗（DLRS MC）方法和自适应超球重要抽样（ARBIS）方法。双层重复抽样蒙特卡罗方法计算结果可以作为对比参照解,但求解效率较低,计算量很大;自适应超球重要抽样方法在满足新指标求解精度的前提下,计算效率得到很大提高。应用数值算例和工程算例证明了所提新指标的意义和所提新算法的高效性。      

像素对图形电磁计算精度的影响和解决方法

在图形电磁计算中,通过控制图像绘制分辨率提高物理光学积分的计算精度.使用像素级的面元近似,把物理光学积分转化为具有数值积分形式的基于像素的求和.使用Nyquist 空间采样理论,提出了能够满足高频RCS(Radar Cross Section)计算精度需要的像素对应的电尺寸值;由此给出了图形电磁计算中的图像分辨率的确定准则,并提出一种改变图像分辨率的方法.通过分析典型目标和复杂目标的计算实例,这种确定图像分辨率的准则可以成立.      

一种计算碰撞概率的分析解

碰撞概率的大小表征航天器相遇过程的危险程度,碰撞概率是航天器规避机动策略制定的参考依据.如果两个飞行器的相对速度很大且相遇过程很短,相对运动可近似为线性运动,这样就可以将三维概率求解问题转化为二维概率问题.这类的碰撞概率求解通常应用数值计算方法.基于概率积分推导一种碰撞概率计算的分析解,能有效降低计算量并不损失太多精度.最终给出模拟算例验证分析解的准确性.     

基于CUDA的超声二维声场EFIT仿真

随着图形处理器（GPU）的快速发展,基于计算设备统一构架（CUDA）可以方便地将并行计算技术应用于超声声场数值仿真计算,极大地提升计算效率。阐述了弹性动力学有限积分算法（EFIT）的原理,在采用CPU实现带吸收边界的钢材料二维点源激励声场仿真的基础上,基于GPU实现了仿真模型的并行计算,介绍了GPU程序的设计流程和参数优化方法,包括纹理内存使用、吸收边界优化和数据传输优化。对比了相同条件下CPU和GPU仿真计算的耗时和平均计算效率,定量分析了GPU对于EFIT模型效率的提升。比对结果表明,EFIT具有良好的并行计算条件,采用并行计算方法能够有效提升模型计算速度,对于复杂声场仿真应用具有广阔的应用前景。      

基于连续小波变换的飞行器结构模态参数辨识

给出了一种基于连续小波变换的多输入多输出MIMO(Multiple-Input Multiple-Output)飞行器结构模态参数的辨识方法.对结构离散运动方程进行连续小波变换建立了小波域内的系统AR(Auto Regression)模型,AR模型的系数矩阵决定着系统的动力学特性,可以通过最小二乘法求得.模态参数可以通过求解由AR系数矩阵构成的特征矩阵的特征值来获得.与已有基于小波变换的模态参数辨识方法相比,该方法应用了连续小波的时移共变性和小波变换的滤波能力来确保辨识的效率.在辨识过程中,采用优化算法提高了辨识的精度和稳定性.算例仿真结果表明该方法具有较高的计算精度和稳定性,能用于飞行器结构模态参数的辨识.      

多星对合作目标的分布式协同导航滤波算法

针对多颗在轨卫星对空间合作目标的协同导航问题,提出了一种适用于协同导航的分布式球面单形-径向容积求积分卡尔曼滤波（DSSRCQKF）算法。为了计算非线性滤波中的高斯加权积分,分别使用球面单形准则和二阶高斯-拉盖尔求积分准则计算球面积分和径向积分,提出了一种新的球面单形-径向容积求积分准则。将该准则嵌入分布式卡尔曼滤波框架中,结合协同导航的非线性数学模型,给出适用于协同导航的DSSRCQKF算法,该算法要求每颗导航星仅与其邻居星进行通信,通过数据的分布式融合实现对目标星轨道状态的一致估计,从而避免了传统集中式处理中较高的通信和计算压力。仿真实验结果表明,与分布式卡尔曼滤波相比,本文算法将对合作目标的实时定位精度提高了11 m,定速精度提高了0.02 m/s,从而验证了本文算法的有效性。      

气动力计算的积分技术讨论

讨论了"尾迹面积分法"的气动力计算方法及其改进技术.基于欧拉方程求解流场,通过对飞行器以外的一个控制体采用动量定理而得,是传统远场积分法的改进.在飞行器所在流场的下游一垂直于来流速度的截面上进行积分计算,得到升力、诱导阻力、激波阻力和总阻力.这种方法的优点是有利于外形复杂物体的气动力积分计算,并可将总阻力按产生的物理机理进行分解,以使设计师对飞行器的气动特点有更为明确的了解.详细讨论了影响这种方法计算精度和效率的多种因素及解决途径.各种数值模拟结果证明了该方法和改进技术的正确性和实用性.      

精细积分方法的稳定性和精度分析

分析了结构动力分析的精细积分方法的稳定性、精度和计算工作量,讨论了离散时间间隔、指数矩阵幂级数展开式的截断阶数L以及2N类算法的阶数N的优化问题.说明了精细积分方法是条件稳定的.综合考虑稳定性、精度和计算工作量,判定截断阶数L取4时精细积分方法的总体效果最好,并给出了N的参数优化公式.最后给出2个例题验证了稳定性和精度分析的正确性.     

模糊变分原理在求解结构固有频率中的应用

在处理含模糊参数的结构固有频率问题时,现有的方法局限于先将模糊参数转化为区间数求解,再由区间解构造出相应的模糊解.由于现有方法计算量较大,所以提出了基于模糊变分原理的模糊里兹法和模糊有限元法.通过在瑞利商变分中直接引入模糊参数,得到了结构固有频率问题的模糊变分原理.作为模糊变分原理的应用,推导了模糊里兹法和模糊有限元法.这两种方法都可以直接得到问题的模糊解,与传统的区间解法相比,提高了计算效率.算例表明,这两种方法不但具有较高的计算精度,而且可以很大程度上减少计算量.     

基于改进的动态Kriging模型的结构可靠度算法

对于复杂航空航天机械产品,极限状态方程往往表现出隐式、高度非线性的特点,而且通常需要调用有限元分析,从而耗费大量时间。将混合粒子群-模拟退火（PSOSA）算法应用到Kriging模型中相关参数的寻优过程,提高了预测精度。同时结合动态更新机制,逐渐加入样本点,尽可能减少函数的调用次数,从而提高了计算效率,并将该算法应用到结构可靠性分析中。通过案例分析,和传统蒙特卡罗模拟方法、响应面等经典方法进行对比,所提算法与蒙特卡罗模拟方法计算结果更加接近,计算时间大大缩短,效率和精度都明显改进。      

泵阀协调控制电动静液作动器方案分析

针对典型的EHA (Electro-Hydrostatic Actuator)系统存在的频响较低的问题,为了兼顾作动系统的效率和频响,将控制阀引入了EHA系统,提出了3种泵阀协调控制的EHA方案,分别是:采用EHSV(Electro-Hydraulic Servo Valve)的EHA系统,采用DDV(Direct Drive Valve)的EHA系统以及采用TPCV(Total Pressure Control Valve)的EHA系统.阐述了这3种方案的系统组成及工作原理,采用AMESim对这3种方案及典型的EHA进行了仿真对比分析.从仿真结果可以看出:泵阀协调控制的EHA系统可以大大提高系统的频响,同时还具有较高的效率.作为3种过渡方案,将对目前机载电动静液作动系统的研制具有实际指导意义.      

多级导流诱导轮与叶轮一体型线优化

通过分析型线参数化系数与结构参数的相关性联系,实现基于型线的参数化改型,从而达到多级导流一体式诱导轮与叶轮性能优化的目的。采用最小二乘法实现叶轮型线拟合,将拟合型线模型的计算流体力学(CFD)仿真结果与原模型的试验结果进行对比,结果表明:仿真结果与原模型试验结果扬程值和效率值误差均小于2%,使用的拟合方法能够实现型线优化。通过分析拟合型线参数化系数与叶轮进出口安装角等的相关性联系,对不同结构参数模型进行型线优化,并对其性能进行仿真对比。结果表明:3组型线模型的内流场整体上正常稳定,其叶轮的各个流道出口和蜗壳隔舌处存在一定范围的低压区,且由于隔舌处漩涡流动的影响使得扩散管出口存在一定区域的高速团。型线3为最优的型线模型,其扬程值与效率值均达到了优化目的,且在设计工况下,其性能提高最为明显。      

计算故障先验概率的最大熵方法

为了解决故障先验概率估算不准的问题,提出了基于最大熵的故障先验概率的计算模型.该模型以相关的先验信息作为最大概率估计的约束条件,并通过拉格朗日函数,将故障先验概率估算问题转化成无约束优化问题.为了实现对无约束优化问题的快速求解,提出了一种基于最速下降法和牛顿法的混合梯度算法;并且,针对大规模系统中故障变量过多的情况,依据系统分解的原则,将高维故障空间分解为多个低维故障空间,给出了低维故障空间求解的快速计算方法.通过最大熵方法和故障平均间隔(MTTF,Mean Time To Failure)方法的结果比较,证明最大熵方法更具准确性.