基于空间算子代数的航天器多体动力学递推实时仿真算法

进一步发展了空间算子代数符号体系,给出了一组坐标系无关空间矢量,张量和矢阵符号表示和相应的坐标投影规则.基于空间算子代数符号体系,推导了适用于航天器实时仿真的无根、树状拓扑、开环多体动力学递推算法,并根据投影规则给出了递推算法在惯性坐标系的投影方程.新的空间算子代数符号体系有符号定义严格,推导过程力学概念清晰,推导结果易于实时仿真程序实现的优点.推导的递推算法非常适合于正逆混合问题和变拓扑问题的求解,可以满足多体航天器动力学实时仿真的要求.     

地下电缆传热分析及其内部温度场数值计算

高压电力输送在城市一般都采用地下电缆的形式.为确保电缆能够安全正常的工作,采用精确方法预计地下电缆的载流量和热特性是很有意义的.采用有限差分法,以单芯电缆为例,进行传热分析,最终可以确定电缆允许的载流量.计算与实验结果表明,完全可以满足工程需要.     

随机集的概率假设密度粒子滤波

多目标跟踪问题中,当目标数已知时,可以用概率数据互联(PDA)或联合概率数据互联(JPDA)算法。而当目标数未知或随时间变化时,需要对不同目标数的跟踪进行比较。可以把目标集看作随机集进行讨论,目标数N是随机变量。随机集的跟踪通过有限集统计(FISST)理论来完成。文中讨论了用粒子滤波实现跟踪随机集的方法。实验表明,在杂波环境下,粒子滤波可以稳健跟踪目标状态和目标数。     

直升机桨叶刚柔耦合特性及计算方法分析

为了研究桨叶刚柔耦合特性,采用中等变形梁模型处理桨叶弹性变形,分别采用小转角和有限转角两种方法处理桨叶绕铰的刚性运动。以桨叶扬起下坠碰撞问题为研究对象。计算分析表明:(1)弹性桨叶存在较明显的刚性运动和弹性变形耦合;(2)挥舞角较小时,两种处理方法计算结果一致,均与试验数据吻合较好;(3)挥舞角较大时,刚柔耦合增强,两种处理方法的计算值差别较为明显;(4)桨叶弯曲刚度对桨叶刚柔耦合特性有较明显影响,刚度越小,耦合越强;(5)桨叶如果比较柔软,小转角处理方法计算值相对有限转角处理方法,相位滞后、幅值增加均较为明显,计算桨叶动响应时需引起注意。      

基于Euler／N—S方程的跨音速非线性静气动弹性问题研究

在C-H网格的基础上,采用Jameson的中心差分有限体积法求解Euler/N-S方程,采用结构影响系数法计算结构的弹性变形,用三角元面积加权法和常体积转换法(CVT)实现流固耦合,计算了弹性后掠机翼在跨音速状态的静气动弹性问题     

应用幂变换法构造低周疲劳寿命预测的幂指函数模型

以GH4049合金低周疲劳数据分析为例,以数据残差检验为基础,说明应用Manson-Coffin方程时普遍存在的问题.从塑性应变幅与疲劳失效反向数在双对数坐标系下的二次曲线特征及残差稳定化角度考虑,基于幂变换方法构造了低周疲劳寿命预测的幂指函数模型,来改善残差图性态及模型预测精度.给出了13种材料的幂变换指数p值及模型寿命预测能力对比结果.分析表明:幂指函数模型对低周疲劳寿命预测结果分散带均在2倍因子以内,具有良好的精度.Manson-Coffin方程实际是幂指函数模型在双对数坐标系下的一阶线性近似,这是幂指函数模型寿命预测精度较高的原因.     

求解Euler方程的无网格与有限体积混合算法研究

研究了无网格方法与有限体积法相结合的、求解Euler方程的混合算法.与单一的无网格方法相比较,由于在大部分计算区域采用了有限体积法,使得发展的算法在运算效率上能与有限体积法相当;同时在物体附近嵌入了无网格区,使得在处理几何外形上更具灵活性.有关算法涉及的无网格与网格区域搭接处理,通过在交界面处引入辅助卫星点和网格单元构成边界信息,实现了区域间的流动信息传递.Euler方程的空间导数分别在两区域用有限体积法和无网格法离散近似,时间方向都采用四步显式Runge-Kutta格式推进求解,数值模拟了喷管内流和绕翼型外流,并分别与整体有限体积法和整体无网格方法进行了比较.算例展示出,提出的混合算法能有效捕捉激波间断,且两区域等值线过渡光滑,算法效率如预期与有限体积法相当,表明该方法是可行的.     

基于库伦摩擦效应的箔片轴承极限承载力理论与试验

基于有限单元法建立了考虑库伦摩擦的波箔型径向气体箔片轴承的箔片结构模型,采用有限差分法和有限单元法耦合求解Reynolds方程和气膜厚度方程,通过求解轴颈达到极限偏心率时的轴承极限承载力,研究了箔片结构库伦摩擦效应对轴承极限承载力的影响规律,并搭建了轴承极限承载力测试试验台,利用温度法测量了两个具有不同轴承壳内表面粗糙度的波箔型径向气体箔片轴承的轴承极限承载力.通过对比分析仿真结果与试验结果表明:轴承壳圆柱孔内表面粗糙度为0.4μm的轴承在10000r/min和20000r/min下,轴承极限承载力分别为15.5N和42.3N;而表面粗糙度为1.6μm的轴承极限承载力为10.9N和29.6N,这是由于波纹箔片和轴承壳体之间的库伦摩擦力增大了波纹箔片的刚度,因此增大箔片结构摩擦因数使得轴承极限承载力降低,并且仿真结果变化趋势和试验结果变化趋势吻合.      

2.5维机织复合材料疲劳寿命预测方法

针对疲劳载荷作用下的2.5维机织复合材料,建立了疲劳寿命预测方法.该方法主要包括单胞模型、疲劳失效判定准则和材料性能退化方法3部分.选取单胞模型为研究对象,利用三维有限元技术进行应力分析;引入改进的三维Hashin疲劳失效准则和Mises准则作为纤维束和树脂基体的疲劳失效判据;采用刚度性能突降准则描述疲劳失效后的材料性能,采用考虑纤维体积分数影响的剩余刚度和剩余强度退化模型描述失效前材料的性能.通过疲劳寿命预测值与试验值的对比,验证了疲劳寿命预测方法的有效性.研究表明:经向拉-拉疲劳寿命随经纱纤维体积分数增大而增加,纬向拉-拉疲劳寿命受纬纱纤维体积分数影响较小.      

基于MIMU的适应多步态的步行导航算法

针对目前国内外步态导航算法中引入的零速检测大都不能很好的对多种步态的行走进行辨识的问题,提出了基于低成本MIMU（微惯性测量单元）、且能兼容多种步态的步行导航算法。算法通过采集安置在行人脚部MIMU输出的测量信息,用捷联惯性积分算法进行导航解算。其间,提出一种新的零速检测方法对行走时变化的步速和步型进行准确辨识,进而找到脚步的零速时刻点,并通过设计的扩展卡尔曼滤波（EKF）对导航解算结果进行零速修正（ZUPT）,实现系统的反馈。最后进行两组实验对算法验证。结果表明,该步态导航算法能对行走时的多步态问题有很好的兼容性,零速修正时刻辨识准确度高,两组实验的导航解算误差均达到0.6%以内,进一步提高了步态导航算法的精度和实用性。