铣槽结构液体火箭发动机推力室壳体热应力分析

介绍了液体火箭发动机推力室铣槽结构热应力的数值分析方法,通过建立液体火箭发动机推力室的流场燃烧和导热理论模型,运用有限体积法考虑液膜冷却计算出发动机工作时的燃气、燃烧室壳体和冷却工质的温度场,将得出的结果作为壳体热应力计算模型的边界条件进行热应力场有限元分析。内、外壁温度的计算数据与实验结果基本相符。     

基于切削热-力耦合效应的表面强化技术及其工艺试验研究

经过对超高强度钢工件高速切削加工表面质量的系统研究,发现在某种特定的切削条件下,工件表面质量将会发生显著改变.该变化不仅表现为表面特征的大幅改善,如表面粗糙度可以控制在0.4μm以下,无宏观裂纹发生等,还表现为表层及次表层显微硬度的明显提高、最大残余压应力深度的大幅增加(达到0.35 mm)以及表面耐磨性能的显著增强等.研究结果表明,合理的切削工艺将极有可能成为一种能够在获得工件结构形状与精度指标的同时,对加工表面表层、次表层的组织特性与工件使用性能产生积极影响的表面强化手段.     

胶粘剂弹性模量对间隙接头应力和强度的影响

 运用弹塑性有限元法,研究了胶粘剂弹性模量对含间隙的铝合金单搭接接头的影响。通过分析胶层和被粘物的应力分布,结果表明：间隙连接对接头端部的应力分布无明显影响,但对间隙附近的应力峰值影响较大;随着胶粘剂弹性模量的增大,间隙附近的应力峰值反而减小;对于高弹性模量胶粘剂接头,间隙对接头的应力峰值无明显影响,在不降低接头名义强度的基础上极大地提高了接头的实际强度;而对于低弹性模量胶粘剂接头,间隙使得接头的应力峰值明显上升,显著地降低了接头的名义强度和实际强度。数值模拟结果与已发表的试验数据大致吻合。     

陶瓷纤维梯度增强活塞的梯度方程研究

 功能梯度材料零件具有单质材料零件无法比拟的理化性能优势,然而由于材料分布复杂以及对功能梯度材料本身性能研究不充分,使性能分析存在很多困难。论文应用复合材料热性能理论,采用有限元分析软件ADINA,分析陶瓷纤维梯度增强活塞（材料梯度方程的参数不同）的温度分布和应力分布,结果表明陶瓷纤维梯度层可以明显改变活塞温度分布,缓和由于热膨胀系数不匹配,在陶瓷纤维增强层与活塞本体交界处产生的应力。根据计算结果拟合出温度峰值、整体应力峰值和层间应力峰值与方程系数之间的曲线,并加以验证。     

一类非线性系统的混沌控制

 描述航天器、陀螺和气浮台等刚体姿态运动的欧拉动力学方程,是一个具有广泛代表意义的三阶非线性方程。当该方程中的参数取不同值时,可得到著名的Lorenz系统、Rssler系统、Newton-Leipnik系统、Chen系统及Lü系统。在不同的外力矩作用下,该动力学系统会呈现出相当复杂的动力学行为。从该系统中,发现了一大类新的混沌吸引子。本文分析了这一类混沌吸引子具有的共同特征,并采用基于输出反馈的PI型控制器将一种新的混沌运动稳定于指定平衡点。仿真结果表明,该控制器能够有效地抑制混沌,能将系统稳定于任意指定的不稳定平衡点。     

A Nonlinear Sub-grid Scale Model for Compressible Turbulent Flow

本文采用直接过滤的Navier-Stokes方程组作为可压缩湍流大涡模拟控制方程组。方程组中因过滤产生的高阶相关项用Taylor级数展开近似,但仅保留级数的一阶导数项。这样产生的误差相当于丢失了模型的部分耗散作用,本文用一种动力学非线性亚格子模型来补偿丢失的耗散影响。本文根据Caylay-Hamilton定理导出了一个非线性亚格子模型,模型中包含的系数由动力学模式确定。与传统的动力学Smagorinsky模式相比,这种动力学非线性模型(DNM)稳定性更好,且不需要在统计均匀方向进行统计平均来计算模型系数,因此减少了计算开销。本文用这种非线性亚格子模型对绕双椭球高超声速湍流进行了模拟,并将所得结果与实验值、计算值及理论值进行了对比。结果表明,本文模型可以有效地模拟可压缩湍流流场。     

含多处损伤搭接结构应力强度因子有限元分析

 采用铆钉的两种表现形式（详细杆元和简化的弹簧元）的组合方式建立了含多处损伤(MSD)多铆钉搭接结构的有限元模型,通过分析得到了不同损伤模式下的应力强度因子(SIF)随裂纹扩展历程变化规律的曲线。结果表明,相对于单裂纹和无腐蚀情况,MSD和腐蚀损伤都会增加SIF值,当裂纹间距大于孔间距一半时,可不考虑裂纹间干涉效应;在给定的循环应力幅下,应力强度因子幅随干涉量的增加而降低,但干涉量大到一定值时,SIF对其不再敏感;铆钉材料对SIF的影响很小,而铆钉直径增大20%,SIF约增加6%。反平面翘曲对SIF有很大的影响,相对于无反平面翘曲,SIF下降一个数量级。     

舵机滚珠丝杠在内部热源作用下响应特性

舵机滚珠丝杠在运行过程中由于摩擦作用和电机运行产生的内部热源会导致其结构温度上升和变形,是影响滚珠丝杠副运行精度的因素之一。本文基于有限元理论建立舵机滚珠丝杠副模型,并通过在ANSYS Workbench中插入命令流的方式实现了在滚珠丝杠上施加移动摩擦热载荷,获得了滚珠丝杠副在内部热源作用下的温度、变形和应力响应特性。结果表明:滚珠丝杠轴向热流密度呈双峰趋势,且其高温区域主要集中在螺母行程范围内的中间区域;随着转速的增加,滚珠丝杠各处温度及其波动幅度加大;丝杠最大热应力位于与轴承配合的轴肩最外侧处。     

An improved numerical method for constructing Halo/Lissajous orbits in a full solar system model

An improved numerical method that can construct Halo/Lissajous orbits in the vicinity of collinear libration points in a full solar system model is investigated. A full solar system gravitational model in the geocentric rotating coordinate system with a clear presentation of the angular velocity relative to the inertial coordinate system is proposed. An alternative way to determine patch points in the multiple shooting method is provided based on a dynamical analysis with Poincare′sections. By employing the new patch points and sequential quadratic programming, Halo orbits for L_1, L_2, and L_3 points as well as Lissajous orbits for L_1 and L_2 points in the EarthMoon system are generated with the proposed full solar system gravitational model to verify the effectiveness of the proposed method.     

单推力航天器交会对接轨迹规划及跟踪控制

针对单推力航天器交会对接问题,提出一种轨迹规划及跟踪算法。首先,考虑到追踪航天器只沿本体X轴安装推力器,且推力方向固定,为了实现从起始位置转移至期望位置并满足姿态要求,基于三维螺旋线设计两阶段转移轨迹,根据初末位置以及末端速度方向要求,求解螺旋线参数。该螺旋线可以保证在初末速度方向固定情况下,曲率积分最小。其次,为了降低轨迹跟踪难度并减小初始时刻的位置跟踪控制力,需要将转移轨迹初始速度与追踪星X轴重合。传统螺旋线无法满足该约束条件。本文对传统螺旋线进行改进,提出一种旋转螺旋线轨迹设计方法。通过引入姿态旋转矩阵,将螺旋线在三维空间旋转,在不改变曲线形状的前提下满足初末位置及速度方向要求。然后,为了跟踪转移轨迹以及跟踪期望推力方向,提出基于CLF（Control Lyapunov Function）的滑模控制策略,当追踪星X轴与期望推力方向夹角较大时,采用CLF,保证最优性;当姿态误差收敛至滑模面附近时,切换为滑模控制,以提升系统鲁棒性。最后,通过仿真验证旋转螺旋线相比于传统螺旋线的优势。