基于模态缩聚法的脉冲推力器模态分析.pdfVIP

中国航天第三专业信息网第27居年会航天动力技术发展与应用学术会议论文集 基于模态缩聚法的脉冲推力器模态分析 隋欣方蜀州张平 (北京理j二大学机电工程学院，北京100081) 摘要：用MSC／Nastran软件对脉冲推力器壳体进行有限元模态分析，给出壳体的前5阶固有频 率和振型，同时用广义动态缩聚法对模型进行模态分析，与精确模型计算结果进行比较，并对模态分 析结果进行分析，找出容易产生共振和有害振型的频率，为进一步进行脉冲推力控制器的动力学分析 奠定基础。 主题词：脉冲推力器；有限元法；模态分析；动态缩聚法 1引言 微型圃体脉冲推力器是高速飞行器飞行姿态控制和弹道修正的一种执行装置。其尺寸微小、响应 快、工作时间短，推力冲量比大，是一种具有明显瞬态特性的特种固体发动机¨】。对脉冲推力器进行 有限元模态分析，确定其固有频率及振型，以避免系统在工作中或者运输过程中发生共振和出现有害 振型，是设计过程中的一个重要环节。本文对脉冲推力器进行了有限元模态分析，指出其薄弱环节， ● 为进一步进行脉冲控制器舱段及至飞行器总体设计的动力学分析奠定基础。 2建立脉冲推力器壳体结构有限元模型 脉冲推力器壳体由点火器壳体、燃烧室壳体、喷管座和喉衬组成，结构简图见图1。利用CATIA口J 软件建立推力器壳体的三维模型，三维模型图见图2。在建立三维模型时，为了简化模型，在不影响 构件强度、刚度和计算精度的情况下．去除小的倒角。 有限元计算的前处理软件采用MSC／Patran[”，用十节点四面体单元对模型进行网格划分，单元边 长定为O．5mm，局部边界区域由软件自动进行调整。模型生成14558个单元。 由于脉冲推力器通过喷管座固定到飞行器舱段上，因此，计算中将喷管座的圆形端面施加约束， 固定其位移为零，其余零件之间实际为螺纹连接，螺纹间隙用密封胶充填。本次计算中，由前处理软 件将其自动胶结。 隋欣等：基于模态缩聚法的脉冲推力器模杰分析 喷臂盔舀定面喉村 燃烧室壳体 点火具壳体 ’ 图1 脉冲推力器壳体结构简图 圈2脉冲推力器三维模型 脉冲推力器壳体材料为Lc4高强度铝合金，喉衬材料为石墨。材料特性见表l。其有限元模型见 圈3。 表1材料特性表 编号 材料 弹性横量(GPa)泊松比 密度(g／cm4)对应部件 1 嵋 70 O．3 2．85 其余结构 2 石墨 7．845 0．1 1．9 喉村 图3脉冲推力器有限元模型 3模态计算“堋 无阻尼自由振动时，弹性体的基本振动方程为： {0}={Ⅱ}陋]+{n}[M] (1) 式中：fM】——质量矩阵；陋】——刚度矩阵；每}——加速度向量；每卜一位移向量 弹性体的自由振动可以分解成一系列自由振动的叠加，对于线性结构，阻】，k】均为实对称矩 阵，方程具有下列简谐运动形式的解 扛G，Y，z，t)／--侈G，Y，z蜘她‘ (2) 式中，移k弘z*为位移矢量的幅值，噱为结构振动的固有频率。 将(2)式代入(1)式，经化简可得噼卜砖【扫】形}={0) (3) 移}有非零解的条件是其系数矩阵奇异，即det妊】一锑拟】)=0 (4) 通过求解方程(4)，可以得到结构系统振动的固有频率。 ～般的计算过程，就是利用结构真实的刚度矩阵和质量矩阵进行求解，这样对于小型的计算是可 行的。但是，当进行较大结构的动态分析时，用精确矩阵计算需要较高的硬件配置，并且用时较长。 436