CN111474384B - 一种薄壁位置冲击信号数据修正方法 - Google Patents
一种薄壁位置冲击信号数据修正方法 Download PDFInfo
- Publication number
- CN111474384B CN111474384B CN202010310606.0A CN202010310606A CN111474384B CN 111474384 B CN111474384 B CN 111474384B CN 202010310606 A CN202010310606 A CN 202010310606A CN 111474384 B CN111474384 B CN 111474384B
- Authority
- CN
- China
- Prior art keywords
- frequency
- data
- thin
- impact
- impact signal
- 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
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P15/00—Measuring acceleration; Measuring deceleration; Measuring shock, i.e. sudden change of acceleration
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L5/00—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
- G01L5/0052—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes measuring forces due to impact
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Discrete Mathematics (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于冲击信号数据修正技术领域,具体涉及一种导弹导引头薄壁位置冲击信号数据修正方法。在某些情况下加速度传感器不得不安装在薄壁件上,薄壁件的局部共振和传感器的附加质量使测到的冲击信号在某些频点被放大,测得的冲击信号并非作用到产品上的真实环境,本发明公开的方法适用于此类冲击数据的修正。本发明主要提高了实测数据的准确性,减弱薄壁件局部模态引起的信号失真,避免根据实测冲击数据制定的试验条件量值偏大。本发明将薄壁件简化为单自由度弹簧阻尼系统,计算系统的幅频响应将其转化为修正系数,通过修正系数还原实测冲击信号的频域幅值从而得到更准确的冲击信号。
Description
技术领域
本发明属于冲击信号数据修正技术领域,具体涉及一种薄壁位置冲击信号数据修正方法。
背景技术
冲击指一个结构受到的瞬态载荷,能激起瞬态扰动的力、位移、速度或加速度的突然变化,并且在很小的空间内释放出相当大的能量。由于冲击运动会给设备带来损伤和破坏,因此需要研究设备所处的冲击环境,并通过模拟环境检查其耐冲击能力。为获取设备所处的冲击环境,需要在设备上的一些关键位置安装加速度传感器,处理传感器获取到的信号得到设备典型的冲击环境。如何使测到的冲击环境更接近设备真实承受的冲击环境显得至关重要。
通常,传感器应该选择安装在结构刚性最好的位置,以免因结构局部模态造成信号失真从而测到的冲击环境不能很好的模拟设备真实经历的冲击环境。对于导弹来说,导引头是导弹发现、识别与跟踪目标的部件,是导弹的核心部件,关乎作战任务的成败,但同时也是导弹上抗冲击能力较差的部件,因此导引头位置的冲击环境备受关注。由于导弹要求轻量化,通常其内部结构非常紧凑,有时为测量导引头的冲击,传感器不得不安装在一些薄壁件上。薄壁件属于结构刚性很差的部位,冲击信号在低频的某一范围内能量被远远放大,以该信号归纳的导引头冲击环境比其实际承受的冲击环境更恶劣。用这种放大后的冲击环境考核导引头的耐冲击能力往往造成导引头的损坏,带来巨大的经济损失。
目前,国内外几乎没有对薄壁件上测得的冲击信号进行修正的有效方法,对此类问题,更多的是采取另行选择其他位置或者改变产品的结构的措施。
发明内容
(一)要解决的技术问题
本发明要解决的技术问题是:如何提供一种因局部模态导致失真的冲击信号修正方法,要求该方法能够在不改变现有结构布局的条件下,使修正后的冲击信号更接近设备真实承受的冲击环境。
(二)技术方案
为解决上述技术问题,本发明提供一种导弹导引头薄壁位置冲击信号数据修正方法,所述方法包括如下步骤:
步骤1:通过理论计算或者实测得到系统的幅频响应特性系数;
其中,所述幅频响应特性的理论计算过程如下:通过实测信号确定被测部位的固有频率fn和阻尼系数ξ;然后将被测部位的薄壁结构件当作单自由度弹簧阻尼系统,把固有频率fn和阻尼系数ξ代入公式(1),可求得被测部位的幅频响应特性系数,即冲击信号在各个频点的放大系数;
公式(1)中η为各个频点的幅频响应特性系数,λ=f/fn为各个频点f与固有频率fn的频率比;
步骤2:对传感器测到的冲击信号进行傅里叶变换,将冲击信号的数据转换为频域数据;
由于实测得到的冲击信号实际上是离散信号,因此采用离散傅里叶变换将其转换为频域数据;假设y={y0,y1,y2,…,yN-1}是冲击信号经过等时间间隔取样后的时间序列,其离散傅里叶变换的频域数据形式为:
其中,i为复数的基本单位,k为时域数据变换到频域数据后的点;
将变换后的频域数据除以对应频点的幅频响应特性系数,即可得到修正后的频域数据;修正后的频域数据表达式为:
步骤3:通过逆傅里叶变换将修正后的频域数据转换到时域,得到修正后的冲击信号;逆傅里叶变换的过程如下:
通过上述三个步骤,即可对受局部模态影响的原始的冲击信号进行修正。
其中,所述步骤1中,通过实测数据得到系统幅频响应系数的方法包括激振器法。
其中,所述步骤1中,通过实测数据得到系统幅频响应系数的方法包括锤击法。
其中,所述步骤1中,通过实测数据得到系统幅频响应系数的方法包括扫频法。
此外,本发明还提供一种导弹导引头薄壁位置冲击信号数据修正方法,其将薄壁结构件简化为单自由度弹簧阻尼系统,计算该单自由度弹簧阻尼系统的幅频响应特性系数,将幅频响应特性系数作为修正系数,将实测信号时域数据通过傅里叶变换转为频域数据,利用该幅频响应特性系数修正频域数据,对修正后的频域数据做逆傅里叶变换得到时域信号。
其中,该方法可推广至其他薄壁结构件的冲击信号修正。
其中,该方法的应用场合并不局限于导弹的导引头。
其中,所述幅频响应特性系数利用固有频率和阻尼比信息理论计算得到,通过有限元计算得到,也可通过试验方法测量得到。
(三)有益效果
与现有技术相比,本发明技术方案将冲击信号的修正问题转化为已知传递系统和系统输出求系统输入的问题,系统的幅频响应特性和系统的输出都可通过试验测得。通过幅频响应可知输入信号经过传递系统后在各个频点的放大系数,再利用系统的输出即传感器测到的信号依据各频点的放大系数将其还原。
最终,本发明可实现在不改变现有结构布局的条件下,测到更接近设备真实情况的冲击信号。
附图说明
图1是典型的幅频响应曲线示意图。
图2-1及图2-2是导引头位置实测信号与冲击响应谱示意图。
图3是导引头薄壁位置沿导弹轴向的幅频响应曲线示意图。
图4-1及图4-2是导引头实测数据和修正信号的对比结果示意图。
图5-1至图5-3是修正信号与仪器舱参考点位实测数据的对比结果示意图。
具体实施方式
为使本发明的目的、内容、和优点更加清楚,下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。
为解决现有技术问题,本发明提供一种导弹导引头薄壁位置冲击信号数据修正方法,所述方法包括如下步骤:
步骤1:通过理论计算或者实测得到系统的幅频响应特性系数;
其中,所述幅频响应特性的理论计算过程如下:通过实测信号确定被测部位的固有频率fn和阻尼系数ξ;然后只考虑主要因素,将被测部位的薄壁结构件当作单自由度弹簧阻尼系统,把固有频率fn和阻尼系数ξ代入公式(1),可求得被测部位的幅频响应特性系数,即冲击信号在各个频点的放大系数;
公式(1)中η为各个频点的幅频响应特性系数,λ=f/fn为各个频点f与固有频率fn的频率比;典型的幅频响应曲线如图1所示。
其中,一个冲击信号是由不同频率成分信号组成的,要求各个频率成分(频点)的真实值。此处,f指代各个频率成分,也就是各个频点,也可称为频段。
步骤2:对传感器测到的冲击信号进行傅里叶变换,将冲击信号的数据转换为频域数据;
由于实测得到的冲击信号实际上是离散信号,因此采用离散傅里叶变换将其转换为频域数据;假设y={y0,y1,y2,…,yN-1}是冲击信号经过等时间间隔取样后的时间序列,其离散傅里叶变换的频域数据形式为:
其中,i为复数的基本单位,k为时域数据变换到频域数据后的点;
将变换后的频域数据除以对应频点的幅频响应特性系数,即可得到修正后的频域数据;修正后的频域数据表达式为:
步骤3:通过逆傅里叶变换将修正后的频域数据转换到时域,得到修正后的冲击信号;逆傅里叶变换的过程如下:
通过上述三个步骤,即可对受局部模态影响的原始的冲击信号进行修正。
所述步骤1中,通过实测数据得到系统幅频响应系数的方法包括激振器法、锤击法、扫频法等,这些方法都比较成熟,在工程上得到了广泛的使用,因此本发明对此不展开介绍。
此外,本发明还提供一种导弹导引头薄壁位置冲击信号数据修正方法,其将薄壁结构件简化为单自由度弹簧阻尼系统,计算该单自由度弹簧阻尼系统的幅频响应特性系数,将幅频响应特性系数作为修正系数,将实测信号时域数据通过傅里叶变换转为频域数据,利用该幅频响应特性系数修正频域数据,对修正后的频域数据做逆傅里叶变换得到时域信号。
其中,该方法可推广至其他薄壁结构件的冲击信号修正。
其中,该方法的应用场合并不局限于导弹的导引头。
其中,所述幅频响应特性系数利用固有频率和阻尼比信息理论计算得到,通过有限元计算得到,也可通过试验方法测量得到。
实施例1
本实施例中,如图2-1及图2-2所示,为导引头上薄壁位置测到的一发射冲击信号,及其冲击响应谱。从冲击响应谱可看到在806Hz附近有一凸点,该频率恰好对应导引头薄壁位置沿导弹轴向方向固有频率,此信号为该位置的冲击响应,并不是导引头位置的承受的冲击载荷。根据该响应信号制定的试验条件考核导引头的耐冲击能力将会造成过试验,因此需要对其进行修正。修正过程如下:
步骤1:通过理论计算或者实测得到系统的幅频响应特性曲线。
已知导引头处薄壁位置沿导弹轴向方向的固有频率为805Hz,参考国军标150.18A-2009可将阻尼系数设置为0.05。按照公式(1)可求得系统的幅频响应即冲击信号在各个频点的放大系数。该薄壁结构的幅频响应曲线见图3。
步骤2:对传感器测到的冲击信号进行傅里叶变换将数据转换到频域,按照公式(3)进行修正。
步骤3:通过逆傅里叶变换将修正后的频域数据转换到时域,得到修正后的时域信号修正后的信号如图4-1及图4-2所示,从图4-1及图4-2中信号的冲击响应谱可看出在806Hz附近的频率凸点得到了修正。
通常,在导弹的仪器舱和导引头舱壁连接段是最佳的测点位置,测到的信号可认为是导引头实际承受的冲击载荷,但往往因为弹体的空间有限此处无法安装传感器。为验证本发明的准确性,此次飞行试验在仪器舱和导引头连接段位置安装了振动加速度传感器作为修正结果的对照,对比结果如图5-1至图5-3所示,从图中可看出修正后的时域信号与仪器舱壁连接段的实测信号几乎一致,因此证明本发明的正确性。该方法可进一步推广应用于导弹导引头薄壁位置发射冲击信号修正。
综上,本发明属于冲击信号数据修正技术领域,具体涉及一种导弹导引头薄壁位置冲击信号数据修正方法。在某些情况下加速度传感器不得不安装在薄壁件上,薄壁件的局部共振和传感器的附加质量使测到的冲击信号在某些频点被放大,测得的冲击信号并非作用到产品上的真实环境,本发明公开的方法适用于此类冲击数据的修正。本发明主要提高了实测数据的准确性,减弱薄壁件局部模态引起的信号失真,避免根据实测冲击数据制定的试验条件量值偏大。本发明将薄壁件简化为单自由度弹簧阻尼系统,计算系统的幅频响应将其转化为修正系数,通过修正系数还原实测冲击信号的频域幅值从而得到更准确的冲击信号。本发明的突出优点是理论依据充分、需要的输入参数少、计算时间短、修正后的数据准确,因此本方法具有推广应用空间。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (6)
1.一种薄壁位置冲击信号数据修正方法,其特征在于,所述方法包括如下步骤:
步骤1:通过理论计算或者实测得到系统的幅频响应特性系数;
其中,所述幅频响应特性的理论计算过程如下:通过实测信号确定被测部位的固有频率fn和阻尼系数ξ;然后将被测部位的薄壁结构件当作单自由度弹簧阻尼系统,把固有频率fn和阻尼系数ξ代入公式(1),可求得被测部位的幅频响应特性系数,即冲击信号在各个频点的放大系数;
公式(1)中η为各个频点的幅频响应特性系数,λ=f/fn为各个频点f与固有频率fn的频率比;
步骤2:对传感器测到的冲击信号进行傅里叶变换,将冲击信号的数据转换为频域数据;
由于实测得到的冲击信号实际上是离散信号,因此采用离散傅里叶变换将其转换为频域数据;假设y={y0,y1,y2,…,yN-1}是冲击信号经过等时间间隔取样后的时间序列,其离散傅里叶变换的频域数据形式为:
其中,i为复数的基本单位,k为时域数据变换到频域数据后的点;
将变换后的频域数据除以对应频点的幅频响应特性系数,即可得到修正后的频域数据;修正后的频域数据表达式为:
步骤3:通过逆傅里叶变换将修正后的频域数据转换到时域,得到修正后的冲击信号;逆傅里叶变换的过程如下:
通过上述三个步骤,即可对受局部模态影响的原始的冲击信号进行修正。
2.如权利要求1所述的薄壁位置冲击信号数据修正方法,其特征在于,所述步骤1中,通过实测数据得到系统幅频响应系数的方法包括激振器法。
3.如权利要求1所述的薄壁位置冲击信号数据修正方法,其特征在于,所述步骤1中,通过实测数据得到系统幅频响应系数的方法包括锤击法。
4.如权利要求1所述的薄壁位置冲击信号数据修正方法,其特征在于,所述步骤1中,通过实测数据得到系统幅频响应系数的方法包括扫频法。
5.如权利要求1所述的薄壁位置冲击信号数据修正方法,其特征在于,该方法可推广至其他薄壁结构件的冲击信号修正。
6.如权利要求1所述的薄壁位置冲击信号数据修正方法,其特征在于,该方法的应用场合并不局限于导弹的导引头。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010310606.0A CN111474384B (zh) | 2020-04-20 | 2020-04-20 | 一种薄壁位置冲击信号数据修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010310606.0A CN111474384B (zh) | 2020-04-20 | 2020-04-20 | 一种薄壁位置冲击信号数据修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111474384A CN111474384A (zh) | 2020-07-31 |
CN111474384B true CN111474384B (zh) | 2022-03-08 |
Family
ID=71755969
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010310606.0A Active CN111474384B (zh) | 2020-04-20 | 2020-04-20 | 一种薄壁位置冲击信号数据修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111474384B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112629637B (zh) * | 2020-11-27 | 2021-10-26 | 华南理工大学 | 一种高频底座力天平信号的时域校准方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100727488B1 (ko) * | 2004-11-05 | 2007-06-13 | 국방과학연구소 | 충격 응답 스펙트럼이 나타나는 시간들의 차이가 짧은충격파형의 합성 방법 및 충격 응답 내역을 구하는 디지털필터와 그 응용 |
CN102346812B (zh) * | 2011-08-11 | 2014-07-09 | 上海理工大学 | 汽车座椅的复杂冲击加速度能量评估方法 |
CN102788613A (zh) * | 2012-07-30 | 2012-11-21 | 中国船舶重工集团公司第七○二研究所 | 实船水下爆炸强冲击环境的分频段测量方法与装置 |
JP6471580B2 (ja) * | 2015-03-27 | 2019-02-20 | 日本電気株式会社 | 衝撃試験装置および衝撃試験装置の調整値算出方法 |
CN105547619A (zh) * | 2015-12-04 | 2016-05-04 | 东北大学 | 一种薄壁构件高阶模态频率及阻尼测试方法及系统 |
CN106644779B (zh) * | 2016-11-29 | 2020-06-09 | 航天东方红卫星有限公司 | 一种星上火工品爆炸冲击谱的等效方法 |
CN106708095B (zh) * | 2016-12-15 | 2019-05-10 | 西安现代控制技术研究所 | 一种导引头动态跟踪数据采集方法 |
CN106706190B (zh) * | 2016-12-15 | 2019-12-24 | 中国航空工业集团公司西安飞机设计研究所 | 一种确定冲击载荷等效静载的方法 |
CN110362928B (zh) * | 2019-07-17 | 2023-01-24 | 大连民族大学 | 确定超弹性薄壁圆柱壳的内共振特性的方法 |
-
2020
- 2020-04-20 CN CN202010310606.0A patent/CN111474384B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111474384A (zh) | 2020-07-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110794170B (zh) | 一种加速度计两自由度动态模型参数辨识的方法 | |
Walter | Air-blast and the science of dynamic pressure measurements | |
CN111474384B (zh) | 一种薄壁位置冲击信号数据修正方法 | |
Yang et al. | Accurate measurement of high-frequency blast waves through dynamic compensation of miniature piezoelectric pressure sensors | |
CN104537251A (zh) | 一种风机叶片的冲击载荷识别方法 | |
CN109902408A (zh) | 一种基于数值运算和改进的正则化算法的载荷识别方法 | |
Yang et al. | Dynamic characterization method of accelerometers based on the Hopkinson bar calibration system | |
Svete et al. | A method for correcting the high-frequency mechanical vibration effects in the dynamic calibration of pressure measurement systems using a shock tube | |
CN110619173B (zh) | 一种利用曲率指标辨识一维杆结构刚度薄弱环节的方法 | |
CN109901393B (zh) | 一种阀控缸电液加速度伺服系统随机波再现控制方法 | |
CN111931288B (zh) | 一种火工冲击响应时频域分布相似性测度方法 | |
Park et al. | Performance improvement of a pressure sensor using fast median algorithm | |
CN110569618B (zh) | 一种卡箍-管路系统频响函数的不确定性分析方法 | |
Chen et al. | The implementation of the data validation process in a gas turbine performance monitoring system | |
Farago et al. | Experimental study on free vibratory behavior of nonlinear structure | |
Coca et al. | Time domain analysis of plasma turbulence observed upstream of a quasi‐parallel shock | |
CN112781722B (zh) | 一种基于有限测点的管路振动全面检测预示方法 | |
Draper et al. | Development of hypersonic wind tunnel force measurement using deconvolution of acceleration data | |
CN117076973B (zh) | 基于对抗性网络的轴承包络阶谱生成方法、设备及介质 | |
Zhong et al. | The Liquid Rocket Engine Experiment Data Quality Improvement Based on 3σ-LMBP | |
Zhang et al. | Dynamic responses of a nonlinear time-varying multi-degree-of-freedom system | |
KR102096002B1 (ko) | 데이터 축약 통합 계측 장치, 이의 제어 방법, 그리고 이 방법을 저장한 컴퓨터 판독 가능 저장 매체 | |
Matsuoka | Research on the Dynamic Performance of Engine Indicators | |
Wisniewiski | Dynamic calibration of pressure transducers with an improved shock tube system | |
Liu et al. | Generation and denoising of vibration noise in shock wave testing |
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 |