冲击响应谱绝对校准研究中冲击响应谱的计算方法

介绍了应用于冲击响应谱绝对校准中的冲击响应谱的计算方法,该方法的实现基于MATLAB编程语言,即通过MATLAB中的仿真直接求解单自由度二阶微分方程的数值解以获得冲击响应谱的绝对复现.     

Foreign Object Damage to Fan Rotor Blades of Aeroengine Part Ⅱ： Numerical Simulation of Bird Impact

Bird impact is one of the most dangerous threats to flight safety. The consequences of bird impact can be severe and, therefore, the aircraft components have to be certified for a proven level of bird impact resistance before being put into service. The fan rotor blades of aeroengine are the components being easily impacted by birds. It is necessary to ensure that the fan rotor blades should have adequate resistance against the bird impact, to reduce the flying accidents caused by bird impacts. Using the contacting-impacting algorithm, the numerical simulation is carded out to simulate bird impact. A three-blade computational model is set up for the fan rotor blade having shrouds. The transient response curves of the points corresponding to measured points in experiments, displacements and equivalent stresses on the blades are obtained during the simulation. From the comparison of the transient response curves obtained from numerical simulation with that obtained from experiments, it can be found that the variations in measured points and the corresponding points of simulation are basically the same. The deforming process, the maximum displacements and the maximum equivalent stresses on blades are analyzed. The numerical simulation verifies and complements the experiment results.     

面向未知区域深度测量的序列图像稠密点特征生成算法

对未知着降区平坦度测量是无人机在复杂地形下安全着陆的关键问题。首先,根据小孔成像原理推导出基于单目序列图像的未知区域深度计算方程;其次,针对稀疏匹配存在深度信息重构误差大而稠密匹配在平滑区域误匹配率高的问题,提出一种基于Delaunay三角剖分的稠密点特征生成算法;然后,分别对序列图像中的2帧图像提取亚像素级Harris角点和尺度不变特征变换(SIFT)特征点,并分别进行特征点匹配;再以2种特征点间的欧氏距离作为约束条件将2种特征点进行融合,生成准稠密特征点;最后,将准稠密特征点进行Delaunay三角剖分,并根据每个剖分三角形上3个顶点像素偏差的方差值制定稠密特征点的生成策略,并结合所提出的深度计算方程计算整个未知区域各点的深度信息。通过Vega Prime(VP)搭建仿真演示验证系统,实验结果表明在机载相机距地面400m处计算高度分别为90m和55m的物体深度信息时,其深度测量相对误差不超过0.89%,具有较高的精度。     

无网格法耦合RNG k-ε湍流模型在亚、跨声速翼型黏性绕流中的数值模拟

基于移动最小二乘无网格方法,耦合RNG(Re-Normalisation Group)k-ε湍流模型求解雷诺平均Navier-Stokes方程。采用AUSM(Advection Upstream Splitting Method)+-up迎风格式求解数值通量,应用在高度各向异性点云结构中取得良好结果的点云重构技术结合移动最小二乘法拟合空间导数,并用三阶SSP(Strong Stability Preserving)型Runge-Kutta显式时间推进格式求解离散后的控制方程。在此基础之上,实现了对NACA0012、RAE2822翼型亚、跨声速黏性绕流的数值模拟,给出了翼型表面压力系数分布曲线、不同位置处的平均速度剖面、马赫数等值线等计算结果,并与实验值及相关文献数值模拟结果进行比较,结果吻合较好。表明所发展的结合点云重构技术的无网格方法耦合RNGk-ε湍流模型能够成功模拟翼型亚、跨声速黏性绕流,验证了所提算法的有效性,并拓展了无网格方法求解湍流流动的途径。     

遥测数据时间序列滑动窗口动态分割技术

经由传感器采集的航天器遥测数据是产品验证、算法检验、数据挖掘等学科和方向重要的数据来源之一。由于原始测量数据采集密度高、信息量大,通常采用滑动窗口技术进行逐段处理和算法应用。但当滑动窗口的宽度固定时,数据的局部信息不能被充分提取出来;为此,提出遥测时序滑动窗口的动态分割流程。仿真结果表明:提出的时序分割和窗口宽度计算方法,可以根据数据自身的特点得到更加准确合理的非等长窗口集合;在保持曲线线段曲率的基础上,滑动窗口动态计算方法可以极大地提取拐点信息,合理地给出窗口宽度,对数据进行了简化与提纯,降低了建模难度,为下一步聚类分析、曲线匹配等工作提供了技术支持,提高了遥测数据的诊断与分析效率。     

带颗粒型金属橡胶夹层阻尼结构的振动响应研究

针对冲压发动机中高温板壳结构的振动抑制需求,提出新型颗粒金属橡胶夹层阻尼结构,基于模态应变能法建立了其动力学理论模型和数值求解方法。与试验结果对比表明,共振幅值的相对误差不大于20%。进一步采用理论分析与试验相结合的手段,论述了颗粒型金属橡胶夹层阻尼的有效性,可有效降低振动响应3~8倍;随着填充密度的增加,减振效果增加;对于一弯和二弯振动,下部填充减振效果最好,而在下部填充,振动响应会被放大;依据数值计算获得的应变能分布结果,进行夹层设计可达到良好的减振效果。     

基于阶跃响应的飞行器单独静、动导数数值求解方法

基于阶跃响应方法,结合刚性动网格技术,对飞行器的单独静、动导数的精细化数值计算进行了研究.以纵向为例,通过给物面施加恒定附加攻角,求解得到阶跃响应运动过程的非定常气动力,求导得到静导数.同样给物面施加恒定的俯仰角速度,并同时强迫物面平动以抵消俯仰转动产生的附加攻角影响,可由非定常气动力求导得到动导数值.分别利用NACA0012翼型和三维SACCON飞翼无人机进行了计算验证,各攻角下的静、动导数值与文献、试验结果吻合得很好,最大误差不超过5%.结论表明:基于阶跃响应的单独静、动导数直接模拟方法计算耗时仅为传统强迫振动方法的21%,效率相对较高,且可推广到横航向的动导数计算,为飞行器的稳定性研究提供参考.      

基于ANSYS Workbench鸟撞飞机风挡有限元分析

利用ANSYS Workbench三维有限元软件,建立了鸟撞飞机风挡的力学分析模型,形成了一套完整合理的鸟撞飞机动态响应分析方法。其中,对模型建立和约束施加等提出了简便的处理方法,并且进行了比较详细的描述。通过有限元分析模拟计算鸟撞飞机前风挡的动态响应,得出了风挡的应力、位移及其分布规律,并与试验结果进行比较,两者吻合较好,可为新型风挡结构设计提供参考。     

柔性接头有效摆心漂移特性

设计了带缠裹式防热套的柔性接头,弹性件和防热套材料均采用硅橡胶.研究了柔性接头有效摆心和几何球心的关系,采用2阶四项式Mooney-Rivlin模型模拟弹性件的本构关系,利用ANSYS软件分别仿真计算柔性接头在压强为0,9MPa下不同摆角条件的有效摆心和几何球心漂移,得到了两者的漂移包络线,把仿真的有效摆心漂移结果与试验测量值进行比较,并利用间隙碰撞试验验证了几何球心漂移仿真结果的准确性,在此基础上研究了有效摆心漂移的原因.结果表明:摆心包络线的仿真结果与试验结果吻合较好,有效摆心和几何球心是两个不同的概念,柔性接头的有效摆心漂移在半径为2.5mm,高为1.4mm的圆柱包络面内,而几何球心的漂移在半径为0.013mm,高为0.6mm的圆柱包络面内,有效摆心漂移是由接头各层增强件和弹性件的变形不同引起的.      

基于三维裂纹尖端应力场的应力强度因子计算方法

提出一种基于无限大体裂纹尖端弹性应力场理论解的前几项多项式函数,对实际裂纹体弹性应力场有限元解进行拟合来计算应力强度因子的方法.该方法在计算应力强度因子时不需要预先假设裂纹尖端的应力应变状态,应力强度因子计算结果更符合三维裂纹体裂纹尖端实际的应力应变状态.首先基于二维无限大板中心穿透裂纹应力场理论解验证了方法的有效性,探讨了拟合确定应力强度因子需要的多项式应力函数的项数.然后分别以二维大板中心穿透裂纹、三维大体内埋圆裂纹和三维有限厚板中心穿透裂纹的应力强度因子计算为例,通过与无限大板和无限大体应力强度因子理论解以及基于位移外推法和1/4节点张开位移法的应力强度因子有限元解的对比分析,验证了该方法的有效性和合理性.研究表明该方法能够合理反映三维裂纹体裂纹尖端的实际应力应变状态,计算得到的应力强度因子数值更合理.