CN110929386A - 利用卫星陀螺数据辨识整星挠性振动模态参数的方法 - Google Patents
利用卫星陀螺数据辨识整星挠性振动模态参数的方法 Download PDFInfo
- Publication number
- CN110929386A CN110929386A CN201911060913.1A CN201911060913A CN110929386A CN 110929386 A CN110929386 A CN 110929386A CN 201911060913 A CN201911060913 A CN 201911060913A CN 110929386 A CN110929386 A CN 110929386A
- Authority
- CN
- China
- Prior art keywords
- satellite
- vibration
- whole
- modal
- flexible
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
- G01C25/005—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass initial alignment, calibration or starting-up of inertial devices
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Navigation (AREA)
Abstract
本发明提供了一种利用卫星陀螺数据辨识整星挠性振动模态参数的方法，利用卫星在轨喷气控制后欠阻尼自由振动时期陀螺测量的姿态角速度数据辨识整星挠性振动模态参数。整星挠性振动模态参数包括在轨挠性振动模态频率和模态阻尼比。辨识过程仅利用卫星平台装载的陀螺测量数据进行分析处理，不使用挠性附件上安装的位移、速度、加速度振动传感器。本方法仅利用卫星平台现有惯性姿态敏感器的测量数据进行分析处理。避免了因安装其他振动传感器而增加的卫星设计、制造难度和在轨运行的风险。辨识出的整星挠性振动模态参数可以为空间挠性部件的结构设计、结构健康监测、结构故障诊断、结构振动控制等方面的应用提供必要的支持。
Description
技术领域
本发明涉及卫星数据辨识领域，具体地，涉及一种利用卫星陀螺数据辨识整星挠性振动模态参数的方法。
背景技术
随着航天技术的不断发展，大型空间结构已经是航天领域中的一个重要发展方向，也将是空间开发必备的基础设施【于登云,夏人伟,孙国江.在轨航天器动力学参数辨识技术研究.中国空间科学技术，2008年2月，第1期】。为获取挠性卫星结构动力学特性，常采用有限元建模、地面试验和在轨辨识等方法计算振动模态参数。有限元建模过程中的模型简化、条件假设等影响了模型精度，特别是难以对铰链等接触式机构精确建模。由于重力和大气阻力等因素的影响，大型挠性结构很难在地面进行组装并实施全尺寸的动力学参数辨识试验。因此，对于带有大型挠性附件的航天器，难以通过有限元建模或地面试验得到精确的结构动力学参数。基于上述因素，对太阳阵、空间展开天线等大挠性结构进行在轨模态参数辨识研究是十分迫切和必要的，同时也具有较高的理论意义和实际应用价值。挠性卫星结构动力学参数包括模态频率、模态阻尼、阵型、耦合系数等，这些参数具有重要的物理意义，可以为空间挠性部件的结构设计、结构健康监测、结构故障诊断、结构振动控制等方面的应用提供必要的支持【CN102982196A】。如【CN103926840A】介绍了一种利用ZVD成型器主动抑制太阳帆板挠性振动的方法，以对挠性模态阻尼比和模态频率作为模型输入。
目前在结构动力学领域的模态参数辨识方法主要包括频域方法、时域方法和近年兴起的时频域方法三种。大多数方法需要利用在轨激励和在挠性附件上布置传感器获取信息，如【李逍然.大型太阳能帆板模态参数在轨辨识研究.哈尔滨工业大学硕士学位论文，2013年6月】研究了挠性附件上传感器的最小配置数量和最优布局方案；【CN105486474A】介绍了一种卫星挠性部件的在轨模态辨识的实现系统及方法，需要对挠性附件进行脉冲激励，接收和监测各个测点的脉冲响应信号，并采集卫星在轨稳态运行过程中卫星挠性部件所产生的加速度响应信号；【CN106557633A】和【CN107609296A】介绍了卫星太阳阵传感器布局的两种方法。【CN106408570A】基于双目视觉测量的太阳帆展开结构动特性在轨辨识方法，从图像中直接提取结构振动位移信息，进而通过工作模态分析技术实时获取结构的动态特性，实现结构动特性的在轨辨识。也有两项专利提出了利用陀螺测量的姿态角速度数据对挠性卫星模态参数的在轨辨识方法：【CN103970964A】介绍了一种挠性卫星模态参数在轨辨识方法，需要采集到执行机构施加到挠性卫星体上的力矩和挠性卫星体相对惯性坐标系的角速度信息，利用子空间辨识算法获得模态参数与力矩到角速度的传递函数。【CN105157728A】提出了一种可抑制陀螺噪声影响的挠性卫星模态参数辨识方法，也是利用卫星在轨飞行时卫星本体角速度的测量数据辨识整星的模态频率和模态阻尼比参数，对陀螺数据进行差分处理抑制这两部分噪声引起的辨识误差。
以上挠性卫星在轨动力学参数辨识方法在工程应用方面有两个约束：一是在轨航天器难于施加动力学参数辨识所需要的已知激励，只能利用在轨航天结构的展开和收拢，结构的对接和分离、发动机的点火等产生激励源进行激励，而这些激励源信号都难于测量【于登云,夏人伟,孙国江.在轨航天器动力学参数辨识技术研究.中国空间科学技术，2008年2月，第1期】。二是在挠性附件上安装振动传感器的数量收到工程实施的限制。一般的挠性附件均需要在轨展开，各类速度、加速度、位移传感器的安装固定、电缆铺设等均对展开机构产生不利影响，增加了设计难度。此外，对于太阳阵等旋转挠性部件还需要增加驱动机构滑环信号通道数量。
为克服以上方法的不足，本专利提出了一种利用自由响应进行辨识、不依赖挠性附件振动传感器的模态参数在轨辨识方法。
发明内容
针对现有技术中的缺陷，本发明的目的是提供一种利用卫星陀螺数据辨识整星挠性振动模态参数的方法。
根据本发明提供的一种利用卫星陀螺数据辨识整星挠性振动模态参数的方法，包括如下步骤：
步骤1：选取卫星在轨喷气闭环控制后的欠阻尼自由振动区的i轴(i＝X,Y,Z)姿态角速度测量数据ωi(t)进行分析；
步骤2：根据地面挠性卫星动力学分析计算结果，设置i轴方向的关注振动频带为[f1,f2]Hz，利用带通滤波器对欠阻尼自由振动区姿态角速度测量数据ωi(t)进行滤波处理，剥离i轴姿态角速度信号中的长周期运动和高频振动分量，仅保留关注频带内的振动信息，得到滤波后的i轴姿态角速度ωi′(t)；
步骤3：对滤波后i轴姿态角速度ωi′(t)进行功率谱密度分析，功率谱密度在关注振动频带[f1,f2]Hz区间内最大值的对应频率即为关注模态的在轨振动频率fi；
步骤4：根据挠性振动模态变量与整星姿态动力学关系式将滤波后的i轴姿态角速度ωi′(t)转化为模态变量的一阶导数
再对模态变量的一阶导数
进行积分，获得模态变量的时间序列ηi(t)；
步骤5：提取模态变量ηi(t)的振动外包络点yi以及其对应的时间ti，用指数函数
对其振动外包络进行拟合，得到拟合系数ai、bi，根据振动包络与阻尼比关系式可计算i轴振动模态的阻尼比ξi。
优选地，挠性振动模态变量与整星姿态动力学关系式的表达式为：
其中，Ji为卫星在i轴方向的转动惯量，Bi为挠性附件振动模态的转动耦合系数。
优选地，振动包络与阻尼比关系式关系式的表达式为：
其中，bi为模态变量振动外包络的指数拟合系数。
优选地，整星挠性振动模态参数通过利用卫星在轨喷气控制后欠阻尼自由振动时期陀螺测量的姿态角速度数据进行辨识。
优选地，所述的整星挠性振动模态参数，包括在轨挠性振动模态频率和模态阻尼比。
优选地，利用卫星陀螺数据辨识整星挠性振动模态参数的辨识过程仅利用卫星平台装载的陀螺测量数据进行分析处理。
优选地，利用带通滤波器对陀螺测量的姿态角速度数据进行滤波处理，剥离卫星本体的长周期运动和其他高频振动误差分量。
与现有技术相比，本发明具有如下的有益效果：
本发明满足挠性卫星在轨结构动力学特性辨识的需求，提取在轨挠性振动模态频率和模态阻尼比等整星挠性振动模态参数。与其他计算卫星在轨结构动力学参数的方法相比，本方法不需要在挠性附件上其他安装传感器，仅利用卫星平台现有惯性姿态敏感器的测量数据进行分析处理。避免了因安装其他振动传感器而增加的卫星设计、制造难度和在轨运行的风险。辨识出的整星挠性振动模态参数可以为空间挠性部件的结构设计、结构健康监测、结构故障诊断、结构振动控制等方面的应用提供必要的支持。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
图1为整星挠性振动模态参数的在轨辨识步骤流程图；
图2为某单翼太阳阵卫星构型示意图；
图3为卫星在轨喷气闭环控制时期的姿态角速度变化曲线图；
图4为卫星姿态角速度滤波前后的时域对比图；
图5为卫星姿态角速度滤波前后的频域对比图；
图6为挠性模态变量的变化曲线和外包络拟合图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变化和改进。这些都属于本发明的保护范围。
本发明提供一种利用卫星陀螺数据辨识整星挠性振动模态参数的方法，利用卫星在轨喷气控制后欠阻尼自由振动时期陀螺测量的姿态角速度数据辨识整星挠性振动模态参数。整星挠性振动模态参数包括在轨挠性振动模态频率和模态阻尼比。辨识过程仅利用卫星平台装载的陀螺测量数据进行分析处理，不使用挠性附件上安装的位移、速度、加速度振动传感器。具体实施例如下：
某遥感卫星安装有单翼太阳电池阵，构型见图1。发射入轨后，在喷气控制期间的滚动X轴姿态角速度测量数据如图2所示。
步骤1：选取卫星在轨喷气闭环控制后的欠阻尼自由振动区的滚动X轴姿态角速度测量数据ωX(t)进行分析(图2)。
步骤2：根据地面挠性卫星动力学分析计算结果，X方向主振动模态的频率约为0.37Hz，设置关注的振动频带为[0.259,0.481]Hz，利用5阶巴特沃兹带通滤波器对欠阻尼自由振动区姿态角速度测量数据ωX(t)进行滤波处理，剥离滚动X轴姿态角速度信号中的长周期运动和高频振动分量，仅保留关注频带内的振动信息，得到滤波后的滚动X轴姿态角速度ωX′(t)。
姿态角速度测量数据滤波前后的时域对比和功率谱密度对比分别见图3和图4。从时域对比和频域对比均可清晰看出，[0.259,0.481]Hz频带以内的信号被保留，其余频段的信号被大幅削弱。
步骤3：对滤波后滚动X轴姿态角速度ωX′(t)进行功率谱密度分析，功率谱密度在[0.259,0.481]Hz区间内最大值的对应频率即为关注模态的在轨振动频率fX。从图4可知，fX＝0.36621Hz。
步骤4：根据式，将滤波后的单轴姿态角速度ωX′(t)转化为模态变量的一阶导数。其中，JX＝6872.13kg·m2为卫星在滚动X轴方向的转动惯量，BX＝33.48762m·kg1/2为挠性附件振动模态的转动耦合系数。再对模态变量的一阶导数进行积分，获得模态变量的时间序列，见图5。
步骤5：提取模态变量的振动外包络点yX以及其对应的时间tX，用指数函数对其振动外包络进行拟合(图5)，得到拟合系数aX、bX。根据式可计算滚动X轴振动模态的阻尼比。本实例中，拟合系数结果为：aX＝3.375×1021、bX＝-0.009565，计算出的阻尼比结果为：ξX＝0.001522。
与其他计算卫星在轨结构动力学参数的方法相比，本方法不需要在挠性附件上其他安装传感器，仅利用卫星平台现有惯性姿态敏感器的测量数据进行分析处理。避免了因安装其他振动传感器而增加的卫星设计、制造难度和在轨运行的风险。辨识出的整星挠性振动模态参数可以为空间挠性部件的结构设计、结构健康监测、结构故障诊断、结构振动控制等方面的应用提供必要的支持。
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本申请的实施例和实施例中的特征可以任意相互组合。
Claims (7)
1.一种利用卫星陀螺数据辨识整星挠性振动模态参数的方法，其特征在于，包括如下步骤：
步骤1：选取卫星在轨喷气闭环控制后的欠阻尼自由振动区的i轴(i＝X,Y,Z)姿态角速度测量数据ωi(t)进行分析；
步骤2：根据地面挠性卫星动力学分析计算结果，设置i轴方向的关注振动频带为[f1,f2]Hz，利用带通滤波器对欠阻尼自由振动区姿态角速度测量数据ωi(t)进行滤波处理，剥离i轴姿态角速度信号中的长周期运动和高频振动分量，仅保留关注频带内的振动信息，得到滤波后的i轴姿态角速度ωi′(t)；
步骤3：对滤波后i轴姿态角速度ωi′(t)进行功率谱密度分析，功率谱密度在关注振动频带[f1,f2]Hz区间内最大值的对应频率即为关注模态的在轨振动频率fi；
步骤4：根据挠性振动模态变量与整星姿态动力学关系式将滤波后的i轴姿态角速度ωi′(t)转化为模态变量的一阶导数
再对模态变量的一阶导数
进行积分，获得模态变量的时间序列ηi(t)；
步骤5：提取模态变量ηi(t)的振动外包络点yi以及其对应的时间ti，用指数函数
对其振动外包络进行拟合，得到拟合系数ai、bi，根据振动包络与阻尼比关系式可计算i轴振动模态的阻尼比ξi。
2.根据权利要求1所述的利用卫星陀螺数据辨识整星挠性振动模态参数的方法，其特征在于，挠性振动模态变量与整星姿态动力学关系式的表达式为：
其中，Ji为卫星在i轴方向的转动惯量，Bi为挠性附件振动模态的转动耦合系数。
3.根据权利要求1所述的利用卫星陀螺数据辨识整星挠性振动模态参数的方法，其特征在于，振动包络与阻尼比关系式关系式的表达式为：
其中，bi为模态变量振动外包络的指数拟合系数。
4.根据权利要求1所述的利用卫星陀螺数据辨识整星挠性振动模态参数的方法，其特征在于，整星挠性振动模态参数通过利用卫星在轨喷气控制后欠阻尼自由振动时期陀螺测量的姿态角速度数据进行辨识。
5.根据权利要求1所述的利用卫星陀螺数据辨识整星挠性振动模态参数的方法，其特征在于，所述的整星挠性振动模态参数，包括在轨挠性振动模态频率和模态阻尼比。
6.根据权利要求1所述的利用卫星陀螺数据辨识整星挠性振动模态参数的方法，其特征在于，利用卫星陀螺数据辨识整星挠性振动模态参数的辨识过程仅利用卫星平台装载的陀螺测量数据进行分析处理。
7.根据权利要求1所述的利用卫星陀螺数据辨识整星挠性振动模态参数的方法，其特征在于，利用带通滤波器对陀螺测量的姿态角速度数据进行滤波处理，剥离卫星本体的长周期运动和其他高频振动误差分量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911060913.1A CN110929386A (zh) | 2019-11-01 | 2019-11-01 | 利用卫星陀螺数据辨识整星挠性振动模态参数的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911060913.1A CN110929386A (zh) | 2019-11-01 | 2019-11-01 | 利用卫星陀螺数据辨识整星挠性振动模态参数的方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110929386A true CN110929386A (zh) | 2020-03-27 |
Family
ID=69850155
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911060913.1A Pending CN110929386A (zh) | 2019-11-01 | 2019-11-01 | 利用卫星陀螺数据辨识整星挠性振动模态参数的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110929386A (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112506245A (zh) * | 2020-11-26 | 2021-03-16 | 西北工业大学 | 一种利用柔性部件根部转动的振动抑制方法 |
| CN113359431A (zh) * | 2021-06-17 | 2021-09-07 | 北京控制工程研究所 | 一种针对航天器挠性振动的在线辨识与抑制方法 |
| CN114184192A (zh) * | 2021-12-27 | 2022-03-15 | 北京计算机技术及应用研究所 | 一种惯性测量装置角速度测量通道传递函数的获取方法 |
| CN117968670A (zh) * | 2024-04-02 | 2024-05-03 | 北京航空航天大学 | 一种微小卫星大角度姿态机动中的高精度指向方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109612665A (zh) * | 2019-01-08 | 2019-04-12 | 上海卫星工程研究所 | 辨识整星挠性振动模态参数的方法与系统 |
| CN109655218A (zh) * | 2019-01-08 | 2019-04-19 | 上海卫星工程研究所 | 用卫星陀螺数据辨识整星挠性振动模态频率的方法及系统 |
| CN109917797A (zh) * | 2019-01-09 | 2019-06-21 | 上海卫星工程研究所 | 利用卫星陀螺数据辨识整星挠性振动模态阻尼方法与系统 |
-
2019
- 2019-11-01 CN CN201911060913.1A patent/CN110929386A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109612665A (zh) * | 2019-01-08 | 2019-04-12 | 上海卫星工程研究所 | 辨识整星挠性振动模态参数的方法与系统 |
| CN109655218A (zh) * | 2019-01-08 | 2019-04-19 | 上海卫星工程研究所 | 用卫星陀螺数据辨识整星挠性振动模态频率的方法及系统 |
| CN109917797A (zh) * | 2019-01-09 | 2019-06-21 | 上海卫星工程研究所 | 利用卫星陀螺数据辨识整星挠性振动模态阻尼方法与系统 |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112506245A (zh) * | 2020-11-26 | 2021-03-16 | 西北工业大学 | 一种利用柔性部件根部转动的振动抑制方法 |
| CN112506245B (zh) * | 2020-11-26 | 2021-10-26 | 西北工业大学 | 一种利用柔性部件根部转动的振动抑制方法 |
| CN113359431A (zh) * | 2021-06-17 | 2021-09-07 | 北京控制工程研究所 | 一种针对航天器挠性振动的在线辨识与抑制方法 |
| CN114184192A (zh) * | 2021-12-27 | 2022-03-15 | 北京计算机技术及应用研究所 | 一种惯性测量装置角速度测量通道传递函数的获取方法 |
| CN114184192B (zh) * | 2021-12-27 | 2023-09-26 | 北京计算机技术及应用研究所 | 一种惯性测量装置角速度测量通道传递函数的获取方法 |
| CN117968670A (zh) * | 2024-04-02 | 2024-05-03 | 北京航空航天大学 | 一种微小卫星大角度姿态机动中的高精度指向方法 |
| CN117968670B (zh) * | 2024-04-02 | 2024-05-28 | 北京航空航天大学 | 一种微小卫星大角度姿态机动中的高精度指向方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110929386A (zh) | 利用卫星陀螺数据辨识整星挠性振动模态参数的方法 | |
| CN108827299B (zh) | 一种基于改进四元数二阶互补滤波的飞行器姿态解算方法 | |
| CN102016605B (zh) | 加速计和用于误差补偿的方法 | |
| Liu et al. | Reaction wheel disturbance modeling, jitter analysis, and validation tests for solar dynamics observatory | |
| CN109612665B (zh) | 辨识整星挠性振动模态参数的方法与系统 | |
| CN109655218B (zh) | 用卫星陀螺数据辨识整星挠性振动模态频率的方法及系统 | |
| CN106989761B (zh) | 一种基于自适应滤波的空间飞行器制导工具在轨标定方法 | |
| CN104215244B (zh) | 基于发射惯性坐标系的空天飞行器组合导航鲁棒滤波方法 | |
| Morelli et al. | Aircraft dynamic modeling in turbulence | |
| Prem et al. | Pseudomeasurement-aided estimation of angle of attack in mini unmanned aerial vehicle | |
| CN109612664B (zh) | 利用陀螺数据辨识卫星挠性附件在轨振动状态方法与系统 | |
| Koehl et al. | Wind-disturbance and aerodynamic parameter estimation of an experimental launched micro air vehicle using an EKF-like observer | |
| CN112345199B (zh) | 一种暂冲式高速风洞迎角传感器振动影响修正方法 | |
| CN109612666B (zh) | 利用陀螺数据辨识卫星挠性附件在轨相对位移方法与系统 | |
| CN109612667B (zh) | 利用陀螺数据辨识卫星挠性附件在轨端部位移方法与系统 | |
| CN114370930A (zh) | 利用mhd微角振动传感器测量卫星结构微角振动的方法 | |
| CN107839903B (zh) | 单翼卫星转移轨道段帆板放气时间估计方法 | |
| Song et al. | Modelling a small-size unmanned helicopter using optimal estimation in the frequency domain | |
| Gülhan | Main achievements of the rocket technology flight experiment ROTEX-T | |
| McCrink | Development of Flight-Test Performance Estimation Techniques for Small Unmanned Aerial Systems | |
| CN109917797A (zh) | 利用卫星陀螺数据辨识整星挠性振动模态阻尼方法与系统 | |
| Banazadeh et al. | Development, instrumentation, and dynamics identification of a coanda air vehicle | |
| Fedosov et al. | Measurement of microaccelerations on board of the LEO spacecraft | |
| Yang et al. | Aviation Sensors and Their Calibration | |
| CN112858960B (zh) | 一种飞机尾杆振动引起的磁干扰的补偿方法及装置 |
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 |