CN106644779B - 一种星上火工品爆炸冲击谱的等效方法 - Google Patents
一种星上火工品爆炸冲击谱的等效方法 Download PDFInfo
- Publication number
- CN106644779B CN106644779B CN201611074392.1A CN201611074392A CN106644779B CN 106644779 B CN106644779 B CN 106644779B CN 201611074392 A CN201611074392 A CN 201611074392A CN 106644779 B CN106644779 B CN 106644779B
- Authority
- CN
- China
- Prior art keywords
- impact
- acceleration
- satellite
- initiating explosive
- explosive device
- 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
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/30—Investigating strength properties of solid materials by application of mechanical stress by applying a single impulsive force, e.g. by falling weight
- G01N3/313—Investigating strength properties of solid materials by application of mechanical stress by applying a single impulsive force, e.g. by falling weight generated by explosives
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/003—Generation of the force
- G01N2203/0032—Generation of the force using mechanical means
- G01N2203/0039—Hammer or pendulum
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
一种星上火工品爆炸冲击谱的等效方法如下:(1)获取星上设备火工品爆炸冲击时域加速度谱值,转换成冲击响应谱;(2)利用冲击响应谱求得获得不同频率点的振动位移;(3)对每一频率点的振动位移累计求和;(4)根据卫星通用地面冲击试验对星上施加的冲击加速谱值,按照步骤(2)和(3),获得地面通用冲击响应谱的位移累计和;(5)根据步骤(3)和(4)确定星上设备的冲击响应谱与地面通用冲击响应谱条件的位移比值。(6)利用获得的比值与星上设备地面冲击试验的通用冲击响应谱,求得星上设备火工品爆炸冲击等效的冲击响应谱。
Description
技术领域
本发明涉及一种星上火工品爆炸冲击谱的等效方法,适用于地面冲击台真实模拟星上火工品爆炸冲击环境,考核评估冲击环境对星上产品的影响。
背景技术
星箭包带、星上太阳翼、天线和伸杆等机构类产品为了实现特定的功能,需要采用火攻品装置解锁,利用装药燃烧爆炸产生的载荷来实现其功能,但是航天器火工品装置爆炸时会产生高量级、宽频带、短时间的复杂振荡性爆炸冲击载荷,在结构中通常以高频应力波的形式传播,会造成航天器电子产品内微型电子芯片、元器件焊点、晶振、电磁继电器等破坏,进而影响产品在轨性能。
由于火工品爆炸冲击环境作用下的产品应力和变形尚缺乏精确有效的工程分析手段,在产品的设计中,对于爆炸分离冲击环境作用下的适应性问题仅能定性地考虑。目前,星上火工品爆炸冲击试验中冲击响应较大的电子设备,需要进行冲击环境的考核和评估,但是使用火工装置模拟爆炸冲击环境十分昂贵,根据标准MIL-HDBK-340A、NASA-STD-7003与GJB 1027A,对产品进行鉴定试验,试验方法一般采用经典的冲击响应谱方法。但是火工品附近冲击较大的产品,冲击谱表现为低频量级较低,某些高频频率峰值较高,地面上的摆锤、电动振动台和气泡冲击台难以模拟该频谱的冲击响应。如按照实际测量的峰值进行冲击试验,由于试验设备的限制,难以达到,也会造成冲击过试验,对冲击敏感的器件造成危害,所以研究火工品爆炸冲击响应与一般试验规范冲击谱等效关系具有重要意义。
发明内容
本发明解决的技术问题为:针对星上火工品爆炸产生的某些频率较大量级的冲击谱,可利用本文提出的冲击等效量级方法,等价为较低量级的试验规范条件,在地面冲击台可真实模拟星上工品爆炸解锁的环境。
本发明解决的技术方案为:一种星上火工品爆炸冲击谱等效的方法,步骤如下:
星上火工品爆炸冲击测量得到的时域加速度是离散的,用表示离散后的第i个时域加速度代入公式(2),0≤i≤t,可得到不同频率单自由度系统受火工品爆炸冲击加速度激励相对基座的位移响应,即δi,根据δi,计算不同频率的单自由度系统加速度公式如下::
(4)把式(4)变换如下:
(5)在设定的频率范围内,对步骤(4)所得星上设备的不同频点受火工品爆炸冲击相对火工品的位移δi求和,得到星上设备的冲击累计位移Df,公式如下:
式中,δi通用为星上设备的不同频点受摆锤冲击相对摆锤的位移;
(7)根据步骤(5)的Df和步骤(6)的Df通用,确定星上设备的冲击响应谱与地面通用冲击响应谱条件的累计位移比值σ。
本发明与现有技术相比的优点在于:
(1)针对火工品爆炸冲击较大的产品,按照频率峰值冲击谱条件进行试验会造成冲击过试验,损伤产品,本发明可以将高量级的冲击等效成低量级的冲击试验,且可以在地面冲击台进行试验,试验手段简单、成熟,比火工品装置模拟爆炸冲击环境更加廉价,并已经在工业领域广泛应用。
附图说明
图1为本发明一种星上火工品爆炸冲击谱的等效方法流程图;
图2为本发明获得的星上火工品爆炸冲击时域加速度曲线;
图3为本发明计算得到的星上火工品爆炸星上设备冲击响应谱曲线;
图4为本发明计算得到的星上设备近似地面冲击试验的通用冲击谱曲线;
图5为本发明计算得到的冲击响应谱不同频率转化振动位移曲线。
具体实施方式
本发明的基本思路为:提出一种星上火工品爆炸冲击谱的等效方法如下:(1)获取星上设备火工品爆炸冲击时域加速度谱值,转换成冲击响应谱;(2)利用冲击响应谱求得获得不同频率点的振动位移,如图5所示;(3)对每一频率点的振动位移累计求和;(4)根据卫星通用地面冲击试验对星上施加的冲击加速谱值,按照步骤(2)和(3),获得地面通用冲击响应谱的位移累计和;(5)根据步骤(3)和(4)确定星上设备的冲击响应谱与地面通用冲击响应谱条件的位移比值。(6)利用获得的比值与星上设备地面冲击试验的通用冲击响应谱,求得星上设备火工品爆炸冲击等效的冲击响应谱。
下面结合附图和具体实施例对本发明进一步说明。
星上火工品爆炸冲击测量得到的时域加速度是离散的,用表示离散后的第i个时域加速度,0≤i≤t,代入公式(2),得到不同频率单自由度系统受火工品爆炸冲击加速度激励相对基座的位移响应,即δi,根据δi,计算不同频率的单自由度系统加速度公式如下:
将卫星火工品爆炸中测量得到的星上设备冲击时域加速度进行冲击响应谱(SRS)转换,如公式(4),得到星上设备受到火工品冲击在频域范围内冲击响应谱,如图2所示,图中可以看处,在星上火工品爆炸冲击中,测量星上6个设备受到冲击响应的冲击响应谱曲线,仅在某个频率点冲击响应较大。
把式(4)变换如下:
根据已经获得了多个星上设备的冲击响应谱值,按照公式(5)将冲击响应谱每一频率点的加速度转化成振动位移,得到的曲线如下图4所示,图中曲线的值相当于每个频率点处的位移大小,与应变正相关联。可以得出2500Hz以后,火工品爆炸冲击对星上设备产生的位移或者应变较小,设备的结构应力较小,说明高频冲击加速度对设备结构影响较小。
经过上述处理后,利用公式(6)对图4中所得的位移进行频域累计,得到星上火工品冲击响应累积位移Df,6个星上设备冲击累计位移值如表1所示。
表1等效冲击条件
计算结果见表1第四行所示,看出等效冲击量级条件比星上设备受到火工品爆炸的冲击谱最大值降低很多,如星上设备E在6703Hz的量级为4493g,经过等效计算后其等效冲击条件为1190g。
为了验证结果正确性,采用弹性力学中应力与应变正相关的理论,在星上设备安装根部粘贴应变片,测量星上火工品爆炸冲击时设备结构应变。将图2中设备A冲击响应谱曲线放大,如图3所示,可以看出,设备A受到的火工品冲击响应曲线与设备地面受摆锤冲击的试验通用冲击响应谱型相近,可以认为相当于480g的通用冲击谱条件。弹性力学中,应变与应力是正相关的,即火工品爆炸冲击测量的设备根部变形与冲击力是正相关的,与测量的冲击加速度值是正相关的,根据火工品爆炸冲击测量其他5个设备安装根部应变值与设备A的应变比例,可以计算出其他设备等效冲击加速度值,并与本发明冲击等效法计算值比较,见下表2所示。可见本发明计算冲击等效法与基于应变等效冲击法比较吻合,证明提出的冲击等效方法是正确合理的,可用于工程应用。
表2等效冲击条件计算与测量比较
星上设备 | 设备A | 设备B | 设备C | 设备D | 设备E | 设备F |
应变比例 | 1 | 1.2423 | 1.7089 | 1.3023 | 1.8938 | 2.1055 |
基于应变等效法(g) | 488 | 606 | 834 | 634 | 927 | 1030 |
冲击等效法(g) | 425 | 665 | 1100 | 650 | 1190 | 1245 |
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (2)
1.一种星上火工品爆炸冲击谱等效的方法,其特征在于步骤如下:
式中,δ(t)为单自由度系统受火工品冲击加速度激励相对其基座的位移,为单自由度系统受到火工品冲击加速度激励相对其基座的加速度;为爆炸冲击的时域加速度,作为基座输入的激励加速度,ωn为单自由度系统的固有圆周率;
式中,δ(0)为t=0时刻单自由度系统受火工品冲击加速度激励相对其基座的位移,为t=0时刻单自由度系统受到火工品冲击加速度激励相对其基座的速度,ωn为单自由度系统的固有圆周率,τ为时间,0~t,为τ时刻的星上设备受到火工品爆炸冲击的时域加速度;
星上火工品爆炸冲击测量得到的时域加速度是离散的,用表示离散后的第i个时域加速度,代入公式(2),0≤i≤t,得到不同频率单自由度系统受火工品爆炸冲击加速度激励相对基座的位移响应,即δi,根据δi,计算不同频率的单自由度系统加速度公式如下:
(4)把式(4)变换如下:
(5)在设定的频率范围内,对步骤(4)所得星上设备的不同频率点受火工品爆炸冲击相对火工品的位移δi求和,得到星上设备的冲击累计位移Df,公式如下:
式中,δi通用为星上设备的不同频率点受摆锤冲击相对摆锤的位移;
(7)根据步骤(5)的Df和步骤(6)的Df通用,确定星上设备的冲击响应谱与地面通用冲击响应谱条件的累计位移比值σ;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611074392.1A CN106644779B (zh) | 2016-11-29 | 2016-11-29 | 一种星上火工品爆炸冲击谱的等效方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611074392.1A CN106644779B (zh) | 2016-11-29 | 2016-11-29 | 一种星上火工品爆炸冲击谱的等效方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106644779A CN106644779A (zh) | 2017-05-10 |
CN106644779B true CN106644779B (zh) | 2020-06-09 |
Family
ID=58813174
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611074392.1A Active CN106644779B (zh) | 2016-11-29 | 2016-11-29 | 一种星上火工品爆炸冲击谱的等效方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106644779B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108507743B (zh) * | 2018-02-02 | 2022-08-12 | 北京航天斯达科技有限公司 | 一种基于lmd的冲击响应谱时域信号合成的方法及装置 |
CN110160894B (zh) * | 2018-04-16 | 2021-08-10 | 北京机电工程研究所 | 冲击载荷试验系统 |
CN110031171B (zh) * | 2019-04-30 | 2020-08-18 | 中国工程物理研究院总体工程研究所 | 一种电动振动台冲击响应谱试验的加权控制方法 |
CN111474384B (zh) * | 2020-04-20 | 2022-03-08 | 西安现代控制技术研究所 | 一种薄壁位置冲击信号数据修正方法 |
CN111964860B (zh) * | 2020-06-30 | 2022-04-12 | 中国空间技术研究院 | 一种新型冲击响应数据判读方法 |
CN112067224B (zh) * | 2020-09-18 | 2022-05-17 | 中煤科工开采研究院有限公司 | 落锤冲击模式下锚杆全程动载响应特征的检测方法 |
CN114492142B (zh) * | 2022-02-25 | 2023-03-28 | 清华大学 | 用于测试航天器元器件抗火工冲击能力的装置及方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103808481A (zh) * | 2012-11-07 | 2014-05-21 | 西北机器有限公司 | 冲击响应谱试验机用冲击回转摆锤机构 |
CN104458453A (zh) * | 2014-11-28 | 2015-03-25 | 中国人民解放军78618部队 | 一种加载于火工品上能量的测试系统及其测试方法 |
CN104568613A (zh) * | 2014-12-19 | 2015-04-29 | 西安近代化学研究所 | 一种基于等效缩比模型的舰船舱室内爆炸结构毁伤评价方法 |
CN105158090A (zh) * | 2015-07-23 | 2015-12-16 | 北京理工大学 | 一种测量材料抗爆炸冲击性能的系统及方法 |
-
2016
- 2016-11-29 CN CN201611074392.1A patent/CN106644779B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103808481A (zh) * | 2012-11-07 | 2014-05-21 | 西北机器有限公司 | 冲击响应谱试验机用冲击回转摆锤机构 |
CN104458453A (zh) * | 2014-11-28 | 2015-03-25 | 中国人民解放军78618部队 | 一种加载于火工品上能量的测试系统及其测试方法 |
CN104568613A (zh) * | 2014-12-19 | 2015-04-29 | 西安近代化学研究所 | 一种基于等效缩比模型的舰船舱室内爆炸结构毁伤评价方法 |
CN105158090A (zh) * | 2015-07-23 | 2015-12-16 | 北京理工大学 | 一种测量材料抗爆炸冲击性能的系统及方法 |
Non-Patent Citations (4)
Title |
---|
Review of pyroshock wave measurement and simulation for space systems;Jung-Ryul Lee et al.;《Measurement》;20120202;第45卷;第631-642页 * |
冲击响应谱试验规范述评;卢来洁等;《振动与冲击》;20021231;第21卷(第2期);第18-20页、第31页 * |
小卫星随机振动特性分析与试验验证方法探讨;杨新峰等;《航天器环境工程》;20140831;第31卷(第4期);第357-362页 * |
爆炸分离冲击环境的模拟试验技术进展;毛勇建等;《导弹与航天运载技术》;20071231(第4期);第37-44页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106644779A (zh) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106644779B (zh) | 一种星上火工品爆炸冲击谱的等效方法 | |
Aykan et al. | Vibration fatigue analysis and multi-axial effect in testing of aerospace structures | |
CN109506874B (zh) | 基于弹性应力波加载的冲击响应谱试验装置及试验方法 | |
CN106525226B (zh) | 一种基于现场振动载荷识别的评估方法及系统 | |
Habtour et al. | Review of response and damage of linear and nonlinear systems under multiaxial vibration | |
CN103983412A (zh) | 用于振动有限元模型修正航空电子设备工作模态测量方法 | |
RU2558688C1 (ru) | Способ испытаний многомассовых систем виброизоляции | |
Bateman et al. | Validation of pyroshock data | |
RU2607361C1 (ru) | Способ испытаний многомассовых систем виброизоляции | |
CN104792450A (zh) | 一种设备对安装基座激励力的间接测量方法 | |
CN111964860B (zh) | 一种新型冲击响应数据判读方法 | |
RU2603826C1 (ru) | Способ исследования двухмассовых систем виброизоляции | |
RU2596232C1 (ru) | Стенд для испытаний многомассовых систем виброизоляции | |
Draper et al. | Development and Implementation of a hybrid dynamic force measurement system at AEDC tunnel 9 | |
Daga et al. | Frequency domain convolutional model of a pyroshock plate for qualification of space equipment | |
Goyder et al. | Measurement of damping in a chain of bolted joints | |
Dang et al. | A novel shock life model for BGA solder joints based on evolutionary power spectral density | |
Suryam et al. | Coherence measurement for early contact detection between two components | |
RU2637718C1 (ru) | Способ исследования двухмассовых систем виброизоляции | |
George et al. | Shock response spectrum of launch vehicle using LabVIEW | |
Cho et al. | Fatigue Life Prediction of Sensor Pod for Aircraft Considering Aircraft Loads | |
Sellami et al. | 3D Finite Volume model for free and forced vibrations computation in on-shore wind turbines | |
Igolkin et al. | Modal analysis of a dynamic mock-up for spacecraft | |
Alexander et al. | Modal analysis of the small spacecraft dynamic model | |
RU2659306C1 (ru) | Способ исследования двухмассовых систем виброизоляции |
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 |