预测-五阶容积卡尔曼滤波方法

为适用于强非线性、非高斯过程噪声系统，结合预测滤波（PF）与高阶容积卡尔曼滤波（HCKF），提出一种预测-五阶容积卡尔曼滤波（P5thCKF）方法。通过预测滤波方法对系统模型中的过程噪声及其方差阵进行实时调整，进而将新模型代入到五阶容积卡尔曼滤波框架中进行实时递推状态估计。推导了五阶球面单形-径向积分准则，采用五阶球面单形积分准则处理球面积分，广义高斯-拉盖尔积分准则处理径向积分；描述了预测滤波方法并对模型误差调整量进行了推导。通过2个仿真实验验证了本文方法在强非线性、非高斯过程噪声系统中的可行性以及应用于工程实践的可能性。     

PVM平台下通量守恒的嵌入自适应网格解激波的底部流

运用嵌入式自适应网格技术来提高结构网格求解复杂流场的分辨率。该方法根据流场的形成,将密集网格嵌入到物理量剧烈变化的区域,使求解精度得以提高的同时,避免较大的计算开销,从而达到求解效率和求解精度之间的平衡,ＰＶＭ平台使得被嵌入的网格可以在网络中不同的机器上完成,大幅度降低了墙上时间。原始网格民被嵌网格之间数值通量衬恒的实现增强了该方法在捕捉间断时的有效性和可靠性     

复合材料层合梁分层问题中应力强度因子与守恒积分的解析变分解法

利用虚功原理导出了涉及两种不同介质的守恒积分表达式以及守恒积分J_Ⅰ与应力强度因子K_Ⅱ的关系式。对于对称正交铺层复合材料层合梁的分层问题,在利用分区广义变分原理确定了应力强度因子K_Ⅱ以及应力与位移本征展开式中的所有待定系数的基础上,进一步求出守恒积分值。最后,作为验证,由J_Ⅰ求出K_Ⅱ,并与由解析变分方法求出的K_Ⅱ相比较。计算表明,本方法前期准备工作少,计算节省机时,结果收敛迅速;沿不同积分路径求出的守恒积分保持恒定,并且用两种方法求出的K_Ⅱ具有很好的一致性。     

动态裂纹扩展韧性 K_(ID) 的估算及方法的评估

通过对30CrMnSi2A钢制成的DCB（带裂尖小孔）试件的试验结果和利用动态有限元法（DFEM）与远场路径无关的J′积分相结合的方法计算了动态裂纹扩展韧性。较全面地对各种估算KID的工程方法进行了分析比较和评估,讨论了使用范围。此外,改进了Kanninen的半经验公式,使它能对于不同泊松比的材料能方便地估算极限速度VL.同时,指出材料常数m应与裂纹的初始长度和止裂长度aa及裂尖半径有关,给出了计算公式,使得估算KID更为简便和合理。     

计算地球流体力学若干新进展

本文简要介绍了有关计算地球流体力学问题研究的若干新进展,其中包括显式平方守恒格式,经济算法,辛算子法,泥沙冲积和三角洲发育的数值模拟等,从中也可以看到计算地球流体力学的发展是具有重要理论意义和实用价值的。     

守恒型可压缩一维湍流方法及其在超声速标量混合层中的应用

一维湍流（ODT）方法是一种能在一维计算域上遵循湍流基本物理规律的湍流建模方法。通过结合确定性和随机性求解方法,能够在一维计算域上准确捕捉到湍流统计规律,且降维建模可显著减小计算量。ODT方法主要被广泛用于不可压湍流和湍流燃烧研究,若要将其拓展用于模拟高速可压缩湍流,需对建模方法进行深度改进。相比于不可压ODT方法,本文基于欧拉参考框架,针对可压缩湍流的特性,将因变量由原始变量改为有利于减小可压缩湍流模拟误差的守恒通量,并加入了组分求解模块。对确定性和随机性求解模块均进行了相应的深度改进,开发出具有标量混合模拟功能的守恒型可压缩ODT方法。在确定性模块中改为求解以守恒通量为变量的一维截断控制方程,在随机性模块中构造一维涡时,将三联映射的作用对象也相应地由原始变量改为守恒通量,并选用了可保证变密度情况下动量守恒的双核变换。通过模拟空间发展超声速平面湍流混合层并将自相似阶段结果与实验结果比对,验证该方法对可压缩剪切湍流场中标量混合的捕捉精度。守恒型可压缩ODT方法模拟得到的速度场和组分场的平均剖面和脉动强度分布与实验结果准确吻合,精度明显优于传统的耦合梯度扩散亚格子模型的大涡模拟方法（LES-GRAD.DIFF.）以及耦合线性涡（LEM）亚格子模型的大涡模拟方法（LES-LEM）,且该方法的降维处理使其在降低计算成本方面具有显著优势。     

基于超网格的重叠网格守恒插值方法

保证重叠网格边界数据插值的守恒性是计算流体力学面临的一大难题。基于超网格方法与格心格式有限体积法,发展了一种新型混合重叠网格守恒插值方法。详细研究了用网格切割技术在重叠网格边界构造局部超网格,用网格求交算法合理地扩大贡献单元模板,建立了一种适用于任意网格类型的隐式并行重叠网格守恒插值方法,以超网格为媒介可实现重叠网格二阶精度的守恒插值。数值结果表明,本文方法对二阶分布的流场变量具有严格的守恒性,相比三线性插值方法和逆向距离权插值方法,本文方法减小了数值误差,提高了重叠网格边界的插值精度,加快了计算收敛速度,改善了重叠区域网格尺度相差较大时流场的光滑性和连续性。     

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

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

控制力矩陀螺磁悬浮转子的指数趋近积分滑模控制研究

针对微重力、动框架等多源扰动下磁悬浮控制力矩陀螺磁轴承转子系统的强非线性、参数摄动及未建模动态等问题,开展控制力矩陀螺的磁轴承转子鲁棒控制研究。建立了磁轴承转子轴向通道动力学模型与电磁力工作点线性化模型,分析了系统参数摄动及未建模动态,并在此基础上提出了一种非线性系统指数趋近积分滑模控制方法,设计了含有电流及位移的积分滑模平面,采用指数趋近律与饱和函数抑制抖振实现到达阶段滑模控制律。实验结果表明：采用该方法可有效抑制滑模的抖振现象,参数摄动时系统具有较佳动态性能,同时在转子转动、框架工作时保证系统良好鲁棒性。     

求解结构可靠度的直接积分方法

针对一类极限状态函数为随机变量的和差积商的有理函数形式,提出了求解结构可靠度的直接积分方法。给出了一些基本函数的概率密度函数。通过定义中间变量将这类特殊的极限状态函数转化为基本函数的函数,进而求出中间变量的概率密度函数,最终求出极限状态函数的概率密度函数.用积分的方法求出结构可靠度。该方法不用考虑随机变量服从何种概率分布,且计算结果为精确解。