CN107885911B - 一种卫星随机振动试验推力评估方法 - Google Patents
一种卫星随机振动试验推力评估方法 Download PDFInfo
- Publication number
- CN107885911B CN107885911B CN201711003130.0A CN201711003130A CN107885911B CN 107885911 B CN107885911 B CN 107885911B CN 201711003130 A CN201711003130 A CN 201711003130A CN 107885911 B CN107885911 B CN 107885911B
- Authority
- CN
- China
- Prior art keywords
- system model
- random vibration
- satellite
- transverse
- thrust
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M7/00—Vibration-testing of structures; Shock-testing of structures
- G01M7/02—Vibration-testing by means of a shake table
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种卫星随机振动试验推力评估方法。该方法沿传递路径方向构建包括卫星在内的系统模态,确定所有主模态及剩余模态,从控制点导出系统实际输入,利用模态与随机振动总均方根力的关系,最终给出卫星随机振动试验推力。本发明解决了目前卫星随机振动推力计算在考虑振动台台体结构影响、试验控制点影响以及高频模态影响等方面的欠缺与不足,实现了提高卫星随机振动试验推力的计算精度的目的。
Description
技术领域
本发明涉及卫星随机振动评估技术,尤其涉及一种卫星随机振动试验推力评估方法。
背景技术
卫星一般在振动台进行随机振动试验,以检验卫星承受随机振动环境的能力。一个振动台有额定的推力,因此,卫星随机振动时需要评估所需推力在振动台额定的推力内。
随机振动的特性是所有频率同时振动,随机振动产生的力是每个频率的力按物体模态特性而综合形成的。当前,卫星随机振动推力一般用有限元方法进行计算,即利用常见的有限元分析软件对卫星按固定边界进行建模,计算出每个频率的传递函数,然后在边界处施加随机振动加速度输入谱进行计算分析,得到卫星边界的总均方根力。然而,卫星安装在振动台上时,振动台台体结构以及夹具会参与到随机振动中的力传递,会影响卫星边界的随机振动推力,如何正确考虑振动台台体结构与夹具的影响还没有很好解决;其次、振动控制点也会影响卫星边界的随机振动推力,也有待探讨;再有,由于计算量的巨大以及高频的较大误差,卫星的基于有限元方法的随机振动分析很难对卫星高频段分析,很多时候没有覆盖全频段实际模态。
上述三个方面的不足,导致当前对卫星的随机振动有限元建模计算的推力结果可能误差很大,无法正确判断一个振动台能否推动卫星进行随机振动试验,出现了试验中振动台推力不足试验中断或者大质量卫星能通过试验而轻质量卫星而无法通过试验等现象。
发明内容
本发明的技术解决问题是:针对目前卫星随机振动推力计算在考虑振动台台体结构影响、试验控制点影响以及高频模态影响等方面的欠缺与不足,提供了一种全路径、全模态、反映控制与响应规律的卫星随机振动试验推力计算方法。
本发明的技术解决方案是:沿传递路径方向构建包括卫星在内的系统模态,确定所有主模态及剩余模态,从控制点导出系统实际输入,利用模态与随机振动总均方根力的关系,最终给出卫星随机振动试验推力,包括如下步骤:
建立卫星有限元模型、振动夹具有限元模型和扩展台面有限元模型;
利用所述卫星有限元模型、振动夹具有限元模型和扩展台面有限元模型,组建系统模型,所述系统模型包括横向系统模型和纵向系统模型;
分别对横向系统模型和纵向系统模型进行模态分析,得到系统模型的各阶主模态和剩余模态;
给定卫星随机振动试验条件,并对所述卫星随机振动试验条件在主频率处的共振放大进行下凹修正处理,以将所述卫星随机振动试验条件转换为卫星随机振动试验下凹条件;
对横向系统模型和纵向系统模型的频响传递函数进行分析,利用卫星随机振动试验下凹条件,确定横向系统模型输入条件和纵向系统模型输入条件;
针对横向系统模型,计算各阶主模态的随机振动等效加速度Aih和剩余模态的随机振动等效加速度Arh;针对纵向系统模型,计算各阶主模态的随机振动等效加速度Aiz和剩余模态的随机振动等效加速度Arz;
在横向振动试验中,计算所述横向系统模型输入条件的总均方根,得到横向刚体部分等效加速度ag,获取振动台滑台、加强台面、振动台牛头和动框的质量总和Mht,根据ag和Mht,计算出横向刚体部分的随机振动所需推力Fh2;在纵向振动试验中,计算所述纵向系统模型输入条件的总均方根,得到振动等效加速度ak,获取动框的质量Mk,根据ak和Mk,计算出纵向刚体部分的随机振动所需推力Fz2;
根据Aih和Arh,计算出横向系统模型的随机振动等效加速度Ah,然后获取横向系统模型的质量Mwj,根据Ah和Mwj,计算出所述横向系统模型所需推力Fh1;根据Fh1和Fh2,计算出卫星横向随机振动试验所需总推力Fh;
根据Aiz和Arz,计算出纵向系统模型的随机振动等效加速度Az,然后获取纵向系统模型的质量Mwjt,根据Az和Mwjt,计算出所述纵向系统模型所需推力Fz1;根据Fz1和Fz2,计算出卫星纵向随机振动试验所需总推力Fz。
进一步地,横向系统模型的各阶主模态的随机振动等效加速度Aih和剩余模态的随机振动等效加速度Arh的计算公式分别为:
Aih=(∑Sh×[(2ξffhi)2+fhi 4]/[(fhi 2-f2)2+(2ξffhi)2])1/2M′hi
Arh=(Ah5% 2/(1-qh 2))1/2
式中,Sh为横向系统模型输入条件中的随机振动加速度谱,f为频率,即随机振动试验中的频段,取值为20Hz-2000Hz,∑代表从20Hz到2000Hz计算求和,fhi为横向系统模型第i阶共振频率,ξ为阻尼,取值为0.05或由测试确定,M′hi为横向系统第i阶模态质量比;qh=1-M′h5%/M'hr;M′h5%和Ah5%分别为横向系统模型模态质量比5%以上的主模态中最后一阶模态的模态质量比和等效加速度,M'hr为横向系统模型剩余模态的模态质量比之和。
进一步地,纵向系统模型的各阶主模态的随机振动等效加速度Aiz和剩余模态的随机振动等效加速度Arz的计算公式分别为:
Aiz=(∑Sz×[(2ξffzi)2+fzi 4]/[(fzi 2-f2)2+(2ξffzi)2])1/2M′zi
Arz=(Az5% 2/(1-qz 2))1/2
式中,Sz为纵向系统模型输入条件中的随机振动加速度谱,f为频率,即随机振动试验中的频段,取值为20Hz-2000Hz,∑代表从20Hz到2000Hz计算求和,fzi为纵向系统模型第i阶共振频率,ξ为阻尼,取值为0.05或由测试确定,M′zi为纵向系统第i阶模态质量比;qz=1-M′z5%/M'zr;M′z5%和Az5%分别为纵向系统模型模态质量比5%以上的主模态中最后一阶模态的模态质量比和等效加速度,M'zr为纵向系统模型剩余模态的模态质量比之和。
进一步地,横向刚体部分的随机振动所需推力Fh2的计算公式为:Fh2=agMht;纵向刚体部分的随机振动所需推力Fz2的计算公式为:Fz2=akMk。
进一步地,横向系统模型的随机振动等效加速度Ah的计算公式为:Ah=(ΣAih 2+Arh 2)1/2;横向系统模型所需推力Fh1的计算公式为:Fh1=AhMwj;卫星横向随机振动试验所需总推力Fh的计算公式为:Fh=Fh1+Fh2。
进一步地,纵向系统模型的随机振动等效加速度Az的计算公式为:Az=(ΣAiz 2+Arz 2)1/2;纵向系统模型所需推力Fz1的计算公式为:Fz1=AzMwjt;卫星纵向随机振动试验所需总推力Fz的计算公式为:Fz=Fz1+Fz2。
本发明与现有技术相比具有如下有益效果:
(a)本发明提出的卫星随机振动试验推力评估方法,由原来单独卫星建模改为沿振动路径的全系统建模,从而可以考虑振动台台体及夹具引起的模态变化(组成与分布)及影响。随机振动中,模态参与不同,模态质量比的高低以及模态频率的位置都大大影响随机振动推力。本方法考虑了随机振动试验中整个振动路径中的参与模态及模态分布,计算更接近真实状态。
(b)本发明考虑解决了原来以卫星试验控制点作为输入而引起的误差。控制卫星的振动输入是运载火箭给定的试验条件所规定的,但卫星试验控制点的输入一般并不是系统的输入,尤其振动夹具刚性较差时。以卫星试验控制点作为输入进行计算一般高出实际随机振动推力。
(c)本发明能覆盖较宽频段和具有较高精度。卫星模态分析结果在卫星动力学分析中是较为准确可靠的,它一般是振动频响分析的基础,因此本发明直接基于模态结果来开展随机振动推力计算。原来的卫星有限元随机振动分析与推力计算,由于卫星的复杂结构,一般很难分析很宽的频段,而且高频段的结果欠准确,本发明利用主模态特性,统计足够的模态质量比,并利用剩余模态弥补原高阶频率计算的不足。
附图说明
图1为本发明提出的一种卫星随机振动试验推力评估方法的流程图;
图2为实施例中的卫星有限元模型;
图3为实施例中的夹具有限元模型;
图4为实施例中的系统有限元模型;
图5为实施例中的卫星随机振动试验条件;
图6为实施例中的卫星随机振动下凹试验条件。
具体实施方式
下面结合附图和实施例对本发明作进一步详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部结构。
图1为本发明提出的一种卫星随机振动试验推力评估方法的流程图。参考图1,本实施例提供的卫星随机振动试验推力评估方法具体可以包括如下步骤:
(1)建立各部分计算模型。分别建立卫星有限元模型,振动夹具有限元模型,振动台的扩展台面有限元模型;其它振动台所属结构如滑台、加强台面、牛头、动框,可不建立有限元模型,只统计质量;
(2)组建系统模型。卫星随机振动试验一般按空间三个方向进行,根据各部分模型组建三个方向的系统有限元模型,即分别建立两个横向振动(X,Y)的系统有限元模型和纵向振动(Z)的系统有限元模型;
(3)主模态与剩余模态计算。分别对横向(X,Y)和纵向(Z)系统模型进行模态分析,得到每个方向各阶模态质量比和对应频率;根据模态质量比确定每个方向的系统主模态,5%以下的模态归纳为剩余模态;
(4)确定卫星随机振动试验下凹条件。对运载火箭给定的卫星随机振动试验条件(三个方向)在主频率处的共振放大进行下凹修正,以消除或抑制卫星在真实发射环境中这些并不存在的共振放大效应。三个方向修正后,卫星随机振动试验条件转换为卫星随机振动试验下凹条件,此下凹条件是卫星随机振动试验时三个方向实际控制输入条件;
(5)确定系统模型的输入条件。系统模型的输入在振动路径上位于卫星振动控制点的前面,它经过一定的路径传递和放大后才是卫星随机振动试验下凹条件量级,根据系统输入至控制点的响应放大关系,由控制点的下凹条件反向计算出系统的输入条件(三个方向);
(6)各阶主模态的随机振动等效加速度计算。按横向和纵向分别计算每个方向各阶主模态对应系统的每个方向各阶主频率。横向系统模型和纵向系统模型中各阶主模态等效加速度计算公式如下:
Aih=(∑Sh×[(2ξffhi)2+fhi 4]/[(fhi 2-f2)2+(2ξffhi)2])1/2M′hi (1)
Aiz=(∑Sz×[(2ξffzi)2+fzi 4]/[(fzi 2-f2)2+(2ξffzi)2])1/2M′zi (2)
式中,Aih为横向系统模型各阶主模态等效加速度,Aiz为纵向系统模型各阶主模态等效加速度;Sh为横向系统模型输入条件中的随机振动加速度谱,Sz为纵向系统模型输入条件中的随机振动加速度谱,Sh和Sz可由步骤(5)得到;f为频率,即随机振动试验中的频段,一般20Hz-2000Hz,∑代表从20Hz到2000Hz计算求和;fhi为横向系统模型第i阶共振频率,fzi为纵向系统模型第i阶共振频率;ξ为阻尼,一般取为0.05或由测试确定;M′hi为横向系统第i阶模态质量比,M′zi为纵向系统第i阶模态质量比,M′hi和M′zi可由步骤(3)得到。
(7)剩余模态的随机振动等效加速度计算。一般模态质量比在5%以下的高阶模态归纳为剩余模态,横向系统模型和纵向系统模型的剩余模态的等效加速度分别按如下公式计算
Arh=(Ah5% 2/(1-qh 2))1/2 (3)
Arz=(Az5% 2/(1-qz 2))1/2 (4)
式中,Arh为横向系统模型剩余模态等效加速度;qh=1-M′h5%/M'hr;M′h5%和Ah5%分别为横向系统模型模态质量比5%以上的主模态中最后一阶模态的模态质量比和等效加速度,M'hr为横向系统模型剩余模态的模态质量比之和。Aiz为纵向系统模型剩余模态等效加速度;qz=1-M′z5%/M'zr;M′z5%和Az5%分别为纵向系统模型模态质量比5%以上的主模态中最后一阶模态的模态质量比和等效加速度,M'zr为纵向系统模型剩余模态的模态质量比之和。
(8)刚体部分等效加速度及推力计算。横向振动试验中,振动台滑台、加强台面、振动台牛头和动框作为刚体考虑,刚体总质量记为Mht,其振动等效加速度即为横向系统模型输入条件Sh的总均方根,可由步骤(5)得到,记为ag,则横向刚体部分(振动台滑台、加强台面、振动台牛头和动框)的随机振动所需推力为
Fh2=agMht (5)
纵向振动中没有振动台滑台、加强台面、振动台牛头,只有动框作为刚体考虑,动框质量记为Mk,其振动等效加速度即为纵向系统模型输入条件Sz的总均方根,可由步骤(5)得到,记为ak,则纵向刚体部分(动框)的随机振动所需推力为
Fz2=akMk (6)
(9)卫星横向随机振动试验推力计算。横向系统模型(卫星与夹具)的随机振动等效加速度Ah为横向各阶主模态等效加速度(由步骤(6)得到)和剩余模态等效加速度(由步骤(7)得到)之和,即
Ah=(ΣAih 2+Arh 2)1/2 (7)
横向系统模型(卫星与夹具)所需推力为
Fh1=AhMwj (8)
其中Mwj为横向系统模型(卫星与夹具)质量。
卫星横向随机振动试验所需总推力为
Fh=Fh1+Fh2=AhMwj+agMht (9)
(10)卫星纵向随机振动试验推力计算。纵向系统模型(卫星、夹具与扩展台面)的随机振动等效加速度Az为纵向各阶主模态等效加速度(由步骤(6)得到)和剩余模态等效加速度(由步骤(7)得到)之和,即
Az=(ΣAiz 2+Arz 2)1/2 (10)
纵向系统模型(卫星、夹具与扩展台面)的纵向推力为
Fz1=AzMwjt (11)
其中Mwjt为纵向系统模型(卫星、夹具与扩展台面)质量。
卫星纵向随机振动试验所需总推力为
Fz=Fz1+Fz2=AzMwjt+akMk (12)
下面进行详细说明:
建立各部分有限元模型
卫星结构较为复杂,有板结构,杆结构,圆壳结构等,材料有金属材料、复合材料等,卫星结构尺寸和材料信息可由卫星设计报告或图纸得到,或由设计三维图直接导入,材料参数可以查阅材料手册。卫星具体建模方法和要求参见中国空间技术研究院院标《航天器模态计算方法(Q/W 678a-2008)》中有关有限元建模的要求以及中国空间技术研究院院标《航天器有限元模型检验方法(Q/W 1403-2015)》;有限元软件一般可选择常用商业分析软件,如Nastran/Patran,Ansys,Abaqus等。
振动夹具用于把卫星固定于振动台上,其结构相对简单,一般为一体化铸造加工出来,材料多为铝合金,建模时有限元单元一般选为体单元,自动生成全部单元。振动台的扩展台面尺寸和材料可由振动台出厂资料查到,其有限元模型与振动夹具建模类似。
组建系统模型
根据卫星在振动态上安装状态来组建系统模型。横向随机振动(X,Y)试验时卫星通过夹具安装在振动台滑台与加强台面上,振动台本身处于横向,通过牛头与滑台和加强台面连接,振动路径为动框→牛头→滑台→加强台面→夹具→卫星,其中动框、牛头、滑台、加强台面在横向表现出较高刚度,可简化为刚体,因此横向振动的系统模型指把卫星、夹具的有限元模型连接一起的分析模型,连接方法一般采用有限元中的MPC连接。横向振动(X,Y)试验中,动框、牛头、滑台、加强台面保持安装状态不变,卫星与夹具在振动台上通过旋转90°形成X和Y两种横向安装。
纵向随机振动(Z)试验时,振动台从横向安装状态调整为纵向安装状态(参照振动台操作程序与要求),此时振动台朝上,卫星通过夹具安装在振动台扩展台面上,振动路径为动框→扩展台面→夹具→卫星,其中动框在纵向简化为刚体,因此纵向振动的系统模型指把卫星、夹具与扩展台面的有限元模型连接一起的分析模型,连接方法一般采用有限元中的MPC连接。
系统模态分析
在系统有限元模型根部施加约束边界,有限元分析软件中设定有效模态质量的输出,然后分别对横向(X,Y)和纵向(Z)系统模型进行模态分析。从模态分析结果中,按频率从低到高提取每个方向主要模态的有效模态质量比和对应频率。一般卫星结构一阶模态质量比最大,然后依次降低,如多数卫星的横向一阶模态质量比在50%-80%之间,这样大约前几阶之和则代表了卫星模态有效质量。因此,从模态结果中,提取模态质量比大于5%的模态即可得到主模态,5%以下的模态归纳为剩余模态。
试验输入条件修正
随机振动试验条件是运载火箭给定的卫星随机振动试验时的输入条件(一般给出三个方向的条件),此条件规定了20-2000Hz频段的加速度功率输入谱,由于卫星做振动试验时在卫星主频率处将产生共振放大效应,此共振效应并不是卫星发射时真实环境,因此,一般卫星随机振动试验时要对一些频段的条件进行下凹修正,以消除或减少非真实的共振放大。
下凹的量级一般参考类似卫星平台的以前随机振动试验下凹量级,带宽约共振频率前后的20Hz。对于没有参考的卫星平台,建议下凹量级为原量级的一半,并在实际试验中根据响应情况进行修正。
系统输入条件计算
随机振动试验下凹条件是卫星随机振动试验的输入条件,即是卫星根部控制点的输入条件,一般来说它不是系统分析模型的输入条件,系统分析模型的输入条件经过振动夹具的传递放大之后才是控制点的随机振动试验下凹条件,因此系统输入条件一般要低于随机振动试验下凹条件。系统输入条件的计算方法为在系统模型根部进行输入为1个单位的频响传递函数分析,由于响应放大主要发生在共振频率处,因此,根据卫星根部控制点的响应放大倍数,在共振频率处对随机振动试验下凹条件除以放大倍数的平方,得到初步的系统输入条件,把此初步的系统输入条件代入系统模型进行随机振动响应分析,根据控制点的响应对系统输入条件进行进一步修正,使控制点的响应与随机振动试验下凹条件符合。
分类计算随机振动推力
a)按式(1)和式(2)计算主模态的随机振动推力。一般计算系统模型前2~3阶主模态的随机振动力则足够,其余模态的模态质量比较小则归于剩余模态。由于随机振动试验输入谱一般为非均匀加速度谱,因此,计算每个主模态的等效加速度时,在20-2000Hz频段内依次对每个频率进行计算求和。所有主模态的等效加速度的平方根乘以系统质量即为系统主模态的随机振动推力。计算时,按三个方向分别计算。
b)按式(3)和式(4)计算剩余模态的随机振动推力。对于一些卫星,主模态的模态质量比不高,比较分散,这时剩余模态起很大作用。剩余模态的等效加速度的平方根乘以系统质量即为系统剩余模态的随机振动推力。计算时,同样按三个方向分别计算。
c)计算振动中的刚体部分随机振动推力。刚体部分的推力计算比较简单,即为刚体部分的随机振动加速度总均方根乘以刚体部分质量。刚体部分不需要建立有限元分析模型,系统分析模型只针对非刚性部分建模,振动计算中系统模型与刚体部分连接,刚体部分的加速度也即系统模型的输入。因此,用系统模型输入条件的总均方根乘以刚体部分质量得到刚体部分随机振动推力。不同方向使用不同方向的刚体质量。
d)总推力计算。卫星随机振动试验中产生的推力不仅仅卫星本身的推力,它是各部分推力之和。对于每一个方向,把此方向系统模型的主模态、剩余模态以及刚体部分的随机振动推力加在一起即是此方向卫星随机振动试验所需要的总推力,把它与振动台额定推力能力比较,则可以评估卫星在此振动台进行试验的安全性和可行性。
实施例:
某XX卫星整星重量约为1050kg,在额定推力9000kgf的振动台做随机振动试验,本例对其横向随机振动试验时的推力进行计算。根据卫星布局设计报告和三维构型图,使用Patran/Nastran有限元软件,建立卫星的有限元模型见图2。试验夹具重约200kg,材料为铝合金,根据夹具的设计报告,使用Patran/Nastran软件建立试验夹具的有限元模型见图3。把夹具的有限元模型导入卫星的有限元模型中,通过位置和方向变换的操作,使夹具与卫星的相对位置与试验时安装位置符合(即夹具与卫星对接环对接),然后使用MPC单元把夹具与卫星对接环连接,形成系统有限元分析模型,见图4。
对系统模型进行模态分析,分析得到横向一阶频率fh1=20.2Hz,一阶模态质量比M'h1=57.3%,二阶频率fh2=85Hz,二阶模态质量比M'h2=19.5%,其它模态质量较小,因此,主模态的计算只取前2阶,其余归为剩余模态。剩余模态的模态质量比之和为M'hr=1-57.3%-19.5%=23.2%。
运载火箭给出的XX卫星横向随机振动试验条件见图5,其中150Hz-600Hz频段的加速度功率谱为0.054g2/Hz,两端的加速度功率谱按规定的斜率变化,全频段(20-2000Hz)的加速度总均方根为6.56g。运载火箭给定的随机振动试验条件在试验中一般要进行下凹,试验中在一阶主频20.2Hz处下凹到0.004g2/Hz,在二阶主频85Hz下凹到0.014g2/Hz,其它频段由于星上设备的考虑也进行了下凹,最终的随机振动下凹条件见图6,随机振动下凹条件的加速度总均方根为6.51g。
随机振动下凹条件用于控制卫星的输入,控制点布置在卫星根部,即夹具上端面。而系统模型输入在夹具底部,夹具底部的输入比夹具上端低,主要表现在系统模型的主频率处降低。使用系统模型进行频响计算,根据夹具底部到夹具上端面的放大关系,确定系统模型在一阶频率附近20Hz-35Hz量级降为0.002g2/Hz,二阶频率附近80Hz-90Hz量级降为0.007g2/Hz,其它频段的量级与随机振动下凹条件相同,这样系统模型随机振动输入的总均方根降为5.6g,系统模型的输入条件各频段加速度功率谱Sh见表1。
表1系统随机振动试验输入条件
系统模型的主模态及模态质量、剩余模态及剩余模态质量、系统模型的输入条件都已确定之后,可以利用公式(1)、公式(3)、公式(7)计算得到系统模型随机振动等效加速度Ah=1.3946g。系统模型(卫星与夹具)的质量为Mwj=1250kg,根据公式(8),计算得到系统模型的随机振动推力为Fh1=1743kgf。
随机振动中刚体部分包括振动台的水平滑台235kg、加强台面300kg、牛头23kg、动框56kg,因此刚体部分质量Mht=614kg,刚体部分振动加速度ag=5.6g,则根据公式(5)刚体部分的随机振动推力为Fh2=3438kgf。
因此,XX卫星在此振动台进行横向振动试验所需总推力为Fh=Fh1+Fh2=5181kgf。本振动台额定最大推力9000kgf,额定最大电流1000安,根据计算得到的5181kgf推力换算振动台的驱动电流约为575安。试验时记录振动台驱动电流530安,说明推力计算结果基本正确。
注意,上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解,本发明不限于这里所述的特定实施例,对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此,虽然通过以上实施例对本发明进行了较为详细的说明,但是本发明不仅仅限于以上实施例,在不脱离本发明构思的情况下,还可以包括更多其他等效实施例,而本发明的范围由所附的权利要求范围决定。
Claims (1)
1.一种卫星随机振动试验推力评估方法,其特征在于,包括:
建立卫星有限元模型、振动夹具有限元模型和扩展台面有限元模型;
利用所述卫星有限元模型、振动夹具有限元模型和扩展台面有限元模型,组建系统模型,所述系统模型包括横向系统模型和纵向系统模型;
分别对横向系统模型和纵向系统模型进行模态分析,得到系统模型的各阶主模态和剩余模态;
给定卫星随机振动试验条件,并对所述卫星随机振动试验条件在主频率处的共振放大进行下凹修正处理,以将所述卫星随机振动试验条件转换为卫星随机振动试验下凹条件;
对横向系统模型和纵向系统模型的频响传递函数进行分析,利用卫星随机振动试验下凹条件,确定横向系统模型输入条件和纵向系统模型输入条件;
针对横向系统模型,计算各阶主模态的随机振动等效加速度Aih和剩余模态的随机振动等效加速度Arh;针对纵向系统模型,计算各阶主模态的随机振动等效加速度Aiz和剩余模态的随机振动等效加速度Arz;
在横向振动试验中,计算所述横向系统模型输入条件的总均方根,得到横向刚体部分等效加速度ag,获取振动台滑台、加强台面、振动台牛头和动框的质量总和Mht,根据ag和Mht,计算出横向刚体部分的随机振动所需推力Fh2;在纵向振动试验中,计算所述纵向系统模型输入条件的总均方根,得到振动等效加速度ak,获取动框的质量Mk,根据ak和Mk,计算出纵向刚体部分的随机振动所需推力Fz2;
根据Aih和Arh,计算出横向系统模型的随机振动等效加速度Ah,然后获取横向系统模型的质量Mwj,根据Ah和Mwj,计算出所述横向系统模型所需推力Fh1;根据Fh1和Fh2,计算出卫星横向随机振动试验所需总推力Fh;
根据Aiz和Arz,计算出纵向系统模型的随机振动等效加速度Az,然后获取纵向系统模型的质量Mwjt,根据Az和Mwjt,计算出所述纵向系统模型所需推力Fz1;根据Fz1和Fz2,计算出卫星纵向随机振动试验所需总推力Fz;
其中,横向系统模型的各阶主模态的随机振动等效加速度Aih和剩余模态的随机振动等效加速度Arh的计算公式分别为:
Aih=(∑Sh×[(2ξffhi)2+fhi 4]/[(fhi 2-f2)2+(2ξffhi)2])1/2M′hi
Arh=(Ah5% 2/(1-qh 2))1/2
式中,Sh为横向系统模型输入条件中的随机振动加速度谱,f为频率,即随机振动试验中的频段,取值为20Hz-2000Hz,∑代表从20Hz到2000Hz计算求和,fhi为横向系统模型第i阶共振频率,ξ为阻尼,取值为0.05或由测试确定,M′hi为横向系统第i阶模态质量比;qh=1-M′h5%/M'hr;M′h5%和Ah5%分别为横向系统模型模态质量比5%以上的主模态中最后一阶模态的模态质量比和等效加速度,M'hr为横向系统模型剩余模态的模态质量比之和;
纵向系统模型的各阶主模态的随机振动等效加速度Aiz和剩余模态的随机振动等效加速度Arz的计算公式分别为:
Aiz=(∑Sz×[(2ξffzi)2+fzi 4]/[(fzi 2-f2)2+(2ξffzi)2])1/2M′zi
Arz=(Az5% 2/(1-qz 2))1/2
式中,Sz为纵向系统模型输入条件中的随机振动加速度谱,f为频率,即随机振动试验中的频段,取值为20Hz-2000Hz,∑代表从20Hz到2000Hz计算求和,fzi为纵向系统模型第i阶共振频率,ξ为阻尼,取值为0.05或由测试确定,M′zi为纵向系统第i阶模态质量比;qz=1-M′z5%/M'zr;M′z5%和Az5%分别为纵向系统模型模态质量比5%以上的主模态中最后一阶模态的模态质量比和等效加速度,M'zr为纵向系统模型剩余模态的模态质量比之和;
所述横向刚体部分的随机振动所需推力Fh2的计算公式为:Fh2=agMht;纵向刚体部分的随机振动所需推力Fz2的计算公式为:Fz2=akMk;
所述横向系统模型的随机振动等效加速度Ah的计算公式为:Ah=(ΣAih 2+Arh 2)1/2;横向系统模型所需推力Fh1的计算公式为:Fh1=AhMwj;卫星横向随机振动试验所需总推力Fh的计算公式为:Fh=Fh1+Fh2;
所述纵向系统模型的随机振动等效加速度Az的计算公式为:Az=(ΣAiz 2+Arz 2)1/2;纵向系统模型所需推力Fz1的计算公式为:Fz1=AzMwjt;卫星纵向随机振动试验所需总推力Fz的计算公式为:Fz=Fz1+Fz2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711003130.0A CN107885911B (zh) | 2017-10-24 | 2017-10-24 | 一种卫星随机振动试验推力评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711003130.0A CN107885911B (zh) | 2017-10-24 | 2017-10-24 | 一种卫星随机振动试验推力评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107885911A CN107885911A (zh) | 2018-04-06 |
CN107885911B true CN107885911B (zh) | 2021-07-13 |
Family
ID=61782371
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711003130.0A Active CN107885911B (zh) | 2017-10-24 | 2017-10-24 | 一种卫星随机振动试验推力评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107885911B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109063237A (zh) * | 2018-06-19 | 2018-12-21 | 上海卫星工程研究所 | 一种适合多附件柔性航天器的系统模态计算方法 |
CN109818668A (zh) * | 2018-12-28 | 2019-05-28 | 航天东方红卫星有限公司 | 一种航天器在非均匀随机振动条件下的载荷预估方法 |
CN110826132B (zh) * | 2019-11-04 | 2021-06-01 | 重庆大学 | 一种结构分散振动控制系统设计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102129490A (zh) * | 2011-02-24 | 2011-07-20 | 航天东方红卫星有限公司 | 一种模块化卫星压紧参数的确定方法 |
CN103744090A (zh) * | 2013-12-31 | 2014-04-23 | 北京邮电大学 | 一种评估卫星导航系统之间兼容性的方法及装置 |
CN104077456A (zh) * | 2014-07-06 | 2014-10-01 | 哈尔滨工业大学 | 空间飞行器姿态控制地面仿真系统效能评估方法 |
-
2017
- 2017-10-24 CN CN201711003130.0A patent/CN107885911B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102129490A (zh) * | 2011-02-24 | 2011-07-20 | 航天东方红卫星有限公司 | 一种模块化卫星压紧参数的确定方法 |
CN103744090A (zh) * | 2013-12-31 | 2014-04-23 | 北京邮电大学 | 一种评估卫星导航系统之间兼容性的方法及装置 |
CN104077456A (zh) * | 2014-07-06 | 2014-10-01 | 哈尔滨工业大学 | 空间飞行器姿态控制地面仿真系统效能评估方法 |
Non-Patent Citations (3)
Title |
---|
Vibration Testing of Small Satellites;Tom Sarafin等;《2017 Instar Engineering and Consulting》;20170719;第1-21页 * |
卫星振动试验中推力的估算方法;李艳辉;《航天器环境工程》;20130430;第30卷(第2期);第160-164页 * |
小卫星随机振动特性分析与试验验证方法探讨;杨新峰等;《航天器环境工程》;20140831;第31卷(第4期);第357-362页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107885911A (zh) | 2018-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107885911B (zh) | 一种卫星随机振动试验推力评估方法 | |
CN107220407B (zh) | 一种变流器柜体振动仿真模型构建方法 | |
CN105912772B (zh) | 一种集频率、振型和频响函数为一体的模型修正方法 | |
CN109556765B (zh) | 一种叶片非接触振动应力测量值换算方法 | |
CN102867097B (zh) | 计及静弹性变形影响的光固化快速成型风洞模型设计方法 | |
CN115630558B (zh) | 一种复合材料构件装配变形预测方法 | |
CN110188417A (zh) | 一种基于多级超单元的复杂结构有限元模型修正方法 | |
CN109117954A (zh) | 基于混合径向基神经网络的黑箱系统设计优化方法 | |
CN105740547B (zh) | 一种基于弹簧-质量块串子模型的设备和开口加筋板耦合振动基频预报方法 | |
CN109408887A (zh) | 一种考虑接触弹性变形的计算转子装配轴线偏斜的有限元分析方法 | |
CN102810127A (zh) | 航天器虚拟振动试验系统 | |
CN110516365A (zh) | 一种测试螺栓连接刚度的方法 | |
CN109580193A (zh) | 一种在爆炸荷载作用下等截面简支梁荷载系数的计算方法 | |
CN104834795B (zh) | 包带连接结构接触摩擦非线性特性模拟方法及系统 | |
CN106354944A (zh) | 一种基于振动模态的磨机基础边界支承刚度的识别方法 | |
CN116698471B (zh) | 一种飞行器舵面静强度试验方法 | |
CN113051787A (zh) | 基于短时实测动态应力的空空导弹吊挂疲劳寿命估算方法和系统 | |
CN112347686B (zh) | 一种仪器安装结构的螺栓校核方法 | |
CN115422649A (zh) | 吊挂结构的腐蚀疲劳寿命预测方法及相关设备 | |
CN107784143B (zh) | 风力发电机组叶片气弹稳定性的确定方法及装置 | |
CN111488684A (zh) | 一种载荷平衡计算方法 | |
CN116011124B (zh) | 一种振动试验系统的仿真建模方法 | |
CN109740194A (zh) | 不等边角钢节间交叉斜材的分析方法 | |
CN114386201B (zh) | 一种轴承的振动特性分析方法 | |
Demenko et al. | Evaluation of spacecraft composite structure load-carrying capability by shock test |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |