CN113820067B - 强冲击传感器下阶跃响应动态特性计算方法及发生装置 - Google Patents

强冲击传感器下阶跃响应动态特性计算方法及发生装置 Download PDF

Info

Publication number
CN113820067B
CN113820067B CN202111382834.XA CN202111382834A CN113820067B CN 113820067 B CN113820067 B CN 113820067B CN 202111382834 A CN202111382834 A CN 202111382834A CN 113820067 B CN113820067 B CN 113820067B
Authority
CN
China
Prior art keywords
impact sensor
frequency characteristic
order system
expression
amplitude
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
Application number
CN202111382834.XA
Other languages
English (en)
Other versions
CN113820067A (zh
Inventor
张振海
张文一
张振山
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Haitai Weina Technology Development Co ltd
Beijing Institute of Technology BIT
Original Assignee
Beijing Haitai Weina Technology Development Co ltd
Beijing Institute of Technology BIT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Haitai Weina Technology Development Co ltd, Beijing Institute of Technology BIT filed Critical Beijing Haitai Weina Technology Development Co ltd
Priority to CN202111382834.XA priority Critical patent/CN113820067B/zh
Publication of CN113820067A publication Critical patent/CN113820067A/zh
Application granted granted Critical
Publication of CN113820067B publication Critical patent/CN113820067B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L25/00—Testing or calibrating of apparatus for measuring force, torque, work, mechanical power, or mechanical efficiency
    • 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/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开的强冲击传感器下阶跃响应动态特性计算方法及发生装置，涉及冲击传感器动态特性计算领域。本发明的方法，对高量程加载条件下产生的下阶跃激励响应信号进行滤波、截取、归一化、翻转处理，得到近似为一阶或二阶系统的强冲击传感器阶跃响应，再通过计算获得传感器阶跃响应动态特性参数，实现强冲击传感器动态特性计算。本发明的装置，使用受冲击发生塑性变形无回弹的脉宽调节缓冲垫作为高量程冲击加载的惯性力传递介质；碰撞过程中，弹丸通过脉宽调节缓冲垫推动砧体加速运动，缓冲结构不断塌陷压缩，压缩过程为塑性变形无回弹；当缓冲结构压缩至极限，弹丸被止挡盘结构阻挡，砧体与弹丸瞬间分离，使得与砧体连接的冲击传感器所受冲击力突然消失为零，产生下阶跃激励及响应。

Description

强冲击传感器下阶跃响应动态特性计算方法及发生装置
技术领域
本发明涉及一种强冲击传感器下阶跃响应动态特性计算方法及发生装置，涉及冲击传感器动态特性计算领域。
背景技术
阶跃函数是一种特殊的连续时间函数，是一个从0跳变到1的过程，属于奇异函数。在传感器动态特性测试中，阶跃激励是理想的激励方法之一。
冲击传感器的静态特性和动态特性在不同量程的激励条件下呈非线性关系，不同量程激励之间无法相互替代。在强冲击传感器测试中，只有高量程的激励才能保证结果的准确度，所以高量程加载条件下强冲击传感器的动态特性计算尤为重要。目前使用激波管可以获得低量程的阶跃激励信号，但高量程的阶跃激励信号很难得到。
发明内容
针对冲击传感器高量程加载条件下，尚无下阶跃响应动态特性计算方法的问题，本发明目的之一是提供一种强冲击传感器下阶跃响应动态特性计算方法，对高量程加载条件下产生的下阶跃激励响应信号进行处理，得到强冲击传感器的阶跃响应，计算获得传感器动态特性参数，进而实现强冲击传感器下阶跃响应动态特性计算。
本发明的另一个目的是提供一种强冲击传感器下阶跃激励发生装置，能够在高量程加载条件下产生下阶跃激励及响应。
本发明是通过以下技术方案实现的。
本发明公开的一种强冲击传感器下阶跃响应动态特性计算方法，对冲击传感器下阶跃响应信号进行滤波处理，截取所得曲线峰值点之后的部分进行数据归一化处理后再上下翻转曲线，产生近似一阶或二阶系统时域曲线，将时域曲线各参数代入一阶或二阶系统幅频特性及相频特性表达式得到冲击传感器幅频特性和相频特性，进而获取工作频带得到冲击传感器的动态特性参数，实现强冲击传感器下阶跃响应动态特性计算。
本发明公开的一种强冲击传感器下阶跃响应动态特性计算方法，包括如下步骤：
步骤一，对冲击传感器下阶跃响应信号进行滤波处理，截取所得曲线峰值点之后的部分进行数据归一化处理后再上下翻转曲线得到冲击传感器下阶跃响应的归一化镜像数据曲线D。
步骤1.1：通过滤波器对冲击传感器下阶跃响应信号进行滤波处理，得到冲击传感器下阶跃响应的滤波曲线A。
步骤1.2：截取所得滤波曲线A峰值点之后的时域信号数据，得到冲击传感器下阶跃响应的截断曲线B。
步骤1.3：以峰值点幅值为标准，将所得截断曲线B所有数据点幅值除以峰值点幅值进行归一化处理，得到冲击传感器下阶跃响应的归一化数据曲线C；
步骤1.4：对所得归一化数据曲线C以峰值点幅值的一半为对称轴上下翻转得到冲击传感器下阶跃响应的归一化镜像数据曲线D。
步骤二，步骤一所得归一化镜像数据曲线D近似为一阶或二阶系统时域曲线，将归一化镜像数据曲线D各参数代入一阶或二阶系统幅频特性及相频特性表达式得到冲击传感器幅频特性和相频特性，进而获取工作频带得到冲击传感器的动态特性参数，实现强冲击传感器下阶跃响应动态特性计算。
步骤2.1：步骤一所得归一化镜像数据曲线D近似为一阶或二阶系统时域曲线。
步骤2.2：当归一化镜像数据曲线D近似为一阶系统时域曲线时，将归一化镜像数据曲线D各参数代入一阶系统幅频特性及相频特性表达式得到冲击传感器幅频特性和相频特性，进而获取工作频带得到冲击传感器的动态特性参数，实现强冲击传感器下阶跃响应动态特性计算。
步骤2.2.1：由归一化镜像数据曲线D得到时间常数τ。
步骤2.2.2：将时间常数τ代入一阶系统幅频特性表达式（1），得到冲击传感器幅频特性。
Figure 454179DEST_PATH_IMAGE002
其中：ω为角频率。
步骤2.2.3：将时间常数τ代入一阶系统相频特性表达式（2），得到冲击传感器相频特性。
Figure 949883DEST_PATH_IMAGE004
步骤2.2.4：将一阶系统幅频特性表达式（1）变换为一阶系统对数幅频特性表达式（3），由一阶系统对数幅频特性表达式（3）得到工作频带ωc表达式（4）。
Figure 10242DEST_PATH_IMAGE006
Figure 173371DEST_PATH_IMAGE008
步骤2.3：当归一化镜像数据曲线D近似为二阶系统时域曲线时，将归一化镜像数据曲线D各参数代入二阶系统幅频特性及相频特性表达式得到冲击传感器幅频特性和相频特性，进而获取工作频带得到冲击传感器的动态特性参数，实现强冲击传感器下阶跃响应动态特性计算。
步骤2.3.1：由归一化镜像数据曲线D得到最大超调量σ和过渡时间ts，将表达式（5）变换为表达式（6），将最大超调量σ和过渡时间ts代入表达式（6）得到阻尼比ζ，并将阻尼比ζ代入表达式（7）得到固有频率ωn。
Figure 90511DEST_PATH_IMAGE010
Figure 339090DEST_PATH_IMAGE012
Figure 999878DEST_PATH_IMAGE014
其中：系数b根据过渡时间误差带范围取值确定。
步骤2.3.2：将阻尼比ζ和固有频率ωn代入二阶系统归一化传递函数表达式（8），并将二阶系统归一化传递函数表达式（8）进行变换得到二阶系统频率响应函数表达式（9），根据二阶系统频率响应函数表达式（9）得到二阶系统归一化幅频特性表达式（10），利用二阶系统归一化幅频特性表达式（10）得到冲击传感器幅频特性。
Figure 17513DEST_PATH_IMAGE016
Figure 308817DEST_PATH_IMAGE018
其中：ω为角频率。
Figure 841429DEST_PATH_IMAGE020
步骤2.3.3：由二阶系统频率响应函数表达式（9）得到二阶系统相频特性表达式（11），利用二阶系统相频特性表达式（11）得到冲击传感器相频特性。
Figure 243592DEST_PATH_IMAGE022
步骤2.3.4：将二阶系统归一化幅频特性表达式（10）变换为二阶系统对数幅频特性表达式（12），由二阶系统对数幅频特性表达式（12）得到二阶系统工作频带ωc表达式（13），二阶系统工作频带ωc表达式（13）有解条件为式（14）所示，则可解得工作频带ωc如式（15）所示。
Figure 912471DEST_PATH_IMAGE024
Figure 374676DEST_PATH_IMAGE026
Figure 394584DEST_PATH_IMAGE028
Figure 334859DEST_PATH_IMAGE030
本发明公开的一种强冲击传感器下阶跃激励发生装置，使用受冲击发生塑性变形无回弹的脉宽调节缓冲垫作为高量程冲击加载的惯性力传递介质。弹丸和脉宽调节缓冲垫碰撞过程中，弹丸通过脉宽调节缓冲垫推动砧体一起加速运动，脉宽调节缓冲垫结构不断塌陷压缩，该塌陷压缩过程为塑性变形无回弹。当脉宽调节缓冲垫结构被压缩至极限，弹丸被止挡盘结构阻挡停止向前运动时，砧体与弹丸瞬间分离，使得与砧体连接的冲击传感器所受冲击力突然消失为零，产生近似理想的下阶跃激励及响应。
作为优选，所述脉宽调节缓冲垫材料选用泡沫铝。
本发明公开的一种强冲击传感器下阶跃激励发生装置，包括弹丸、脉宽调节缓冲垫、止挡盘、砧体和冲击传感器。脉宽调节缓冲垫放置在砧体凹槽中。脉宽调节缓冲垫和砧体放置在止挡盘筒内。止挡盘与炮管固定连接。弹丸、脉宽调节缓冲垫、止挡盘、砧体和冲击传感器沿空气炮炮管轴线方向同轴安装布置。弹丸被压缩空气推动由弹舱位置沿着空气炮炮管加速向前运动，弹丸和脉宽调节缓冲垫碰撞过程中，弹丸通过脉宽调节缓冲垫推动砧体一起加速运动，脉宽调节缓冲垫结构不断塌陷压缩，该塌陷压缩过程为塑性变形无回弹。当脉宽调节缓冲垫被压缩至极限，弹丸被止挡盘结构阻挡停止向前运动时，砧体与弹丸瞬间分离，使得与砧体连接的冲击传感器所受冲击力突然消失为零，产生近似理想的下阶跃激励及响应。
有益效果：
1、针对现有技术中尚无专门针对强冲击传感器下阶跃响应信号的处理方法，本发明公开的一种强冲击传感器下阶跃响应动态特性计算方法，对下阶跃激励响应信号进行滤波处理，截取所得曲线峰值点之后的部分进行数据归一化处理后再上下翻转曲线，产生近似一阶或二阶系统时域曲线，将时域曲线各参数代入一阶或二阶系统幅频特性及相频特性表达式得到冲击传感器幅频特性和相频特性，进而获取工作频带得到冲击传感器的动态特性参数，实现强冲击传感器下阶跃响应动态特性计算。
2、本发明公开的一种强冲击传感器下阶跃响应动态特性计算方法，通过滤波器对冲击传感器下阶跃响应信号进行滤波处理，得到冲击传感器下阶跃响应的滤波曲线A；截取所得曲线峰值点之后的时域信号数据，得到冲击传感器下阶跃响应的截断曲线B；以峰值点幅值为标准，将所得截断曲线B所有数据点幅值除以峰值点幅值进行归一化处理，得到冲击传感器下阶跃响应的归一化数据曲线C；对所得归一化数据曲线C以峰值点幅值的一半为轴上下翻转得到冲击传感器下阶跃响应的归一化镜像数据曲线D，使得归一化镜像数据曲线D能够采用现有阶跃信号处理方法获取下阶跃响应动态特性。
3、针对冲击传感器高量程加载条件下，阶跃激励不易实现的问题，本发明公开的一种强冲击传感器下阶跃激励发生装置，使用受冲击发生塑性变形无回弹的脉宽调节缓冲垫作为高量程冲击加载的惯性力传递介质。弹丸和脉宽调节缓冲垫碰撞过程中，弹丸通过脉宽调节缓冲垫推动砧体一起加速运动，脉宽调节缓冲垫结构不断塌陷压缩，该塌陷压缩过程为塑性变形无回弹。当脉宽调节缓冲垫结构被压缩至极限，弹丸被止挡盘结构阻挡停止向前运动时，砧体与弹丸瞬间分离，使得与砧体连接的冲击传感器所受冲击力突然消失为零，产生近似理想的下阶跃激励及响应。
4、本发明公开的一种强冲击传感器下阶跃激励发生装置，通过大量试验和数据分析，选出能够满足受冲击发生塑性变形无回弹的泡沫铝作为脉宽调节缓冲垫材料，经实验验证选用泡沫铝作为脉宽调节缓冲垫材料，能够产生近似理想的下阶跃激励及响应。
附图说明
图1为本发明所述计算方法对下阶跃响应信号处理过程的流程图。
图2为通过本发明所述计算方法步骤1.1得到的冲击传感器下阶跃响应的滤波曲线A。
图3为通过本发明所述计算方法步骤1.2得到的冲击传感器下阶跃响应的截断曲线B。
图4为通过本发明所述计算方法步骤1.3得到的冲击传感器下阶跃响应的归一化数据曲线C。
图5为通过本发明所述计算方法步骤1.4得到的冲击传感器下阶跃响应的归一化镜像数据曲线D。
图6为通过本发明所述计算方法处理后得到的近似一阶系统时域曲线。
图7为通过本发明所述计算方法计算后得到的近似一阶系统幅频响应曲线。
图8为通过本发明所述计算方法计算后得到的近似一阶系统相频响应曲线。
图9为通过本发明所述计算方法处理后得到的近似二阶系统时域曲线。
图10为通过本发明所述计算方法计算后得到的近似二阶系统幅频响应曲线。
图11为通过本发明所述计算方法计算后得到的近似二阶系统相频响应曲线。
图12为本发明公开的一种强冲击传感器下阶跃激励发生装置的结构示意图。
其中：1—弹丸、2—脉宽调节缓冲垫、3—止挡盘、4—砧体、5—冲击传感器。
图13为本发明所述下阶跃激励发生过程。
具体实施方式
为了表明本发明解决的技术问题及有益效果，下面结合附图和实施例对发明内容做进一步说明。
图1为本发明所述计算方法对下阶跃响应信号处理过程的流程图。
参照图1，本发明所述的利用下阶跃激励响应信号对高量程加载条件下强冲击传感器进行动态特性计算的方法，其具体步骤为：
步骤一，对冲击传感器下阶跃响应信号进行滤波处理，截取所得曲线峰值点之后的部分进行数据归一化处理后再上下翻转曲线得到冲击传感器下阶跃响应的归一化镜像数据曲线D。
步骤1.1：记下阶跃激励条件下冲击传感器加速时间即下阶跃响应信号波形持续时间为T，则选择FIR低通滤波器通带截止频率为10/T，阻带频率为15/T，通带波动为3 dB，阻带衰减为50 dB。使用该参数滤波器对冲击传感器下阶跃响应信号进行滤波处理，得到冲击传感器下阶跃响应的滤波曲线A如图2所示。
步骤1.2：截取所得曲线A峰值点之后的时域信号数据，得到冲击传感器下阶跃响应的截断曲线B如图3所示。
步骤1.3：以峰值点幅值为标准，将所得曲线B所有数据点幅值除以峰值点幅值进行归一化处理，得到冲击传感器下阶跃响应的归一化数据曲线C如图4所示。
步骤1.4：对所得归一化数据曲线C以峰值点幅值的一半即0.5为对称轴上下翻转得到冲击传感器下阶跃响应的归一化镜像数据曲线D如图5所示。
步骤二，步骤一所得归一化镜像数据曲线D近似为一阶或二阶系统时域曲线，将归一化镜像数据曲线D各参数代入一阶或二阶系统幅频特性及相频特性表达式得到冲击传感器幅频特性和相频特性，进而获取工作频带得到冲击传感器的动态特性参数，实现强冲击传感器下阶跃响应动态特性计算。
步骤2.1：步骤一所得归一化镜像数据曲线D近似为一阶或二阶系统时域曲线。
步骤2.2：当归一化镜像数据曲线D近似为一阶系统时域曲线时，将归一化镜像数据曲线D各参数代入一阶系统幅频特性及相频特性表达式得到冲击传感器幅频特性和相频特性，进而获取工作频带得到冲击传感器的动态特性参数，实现强冲击传感器下阶跃响应动态特性计算。
步骤2.2.1：取图6所示的归一化镜像数据曲线D到达稳定幅值的63.2%，即幅值归一化后的0.632处对应的时间值τ，τ即为时间常数。
步骤2.2.2：将时间常数τ代入一阶系统幅频特性表达式（1），得到冲击传感器幅频特性。
Figure 327085DEST_PATH_IMAGE032
其中：ω为角频率。
步骤2.2.3：将时间常数τ代入一阶系统相频特性表达式（2），得到冲击传感器相频特性。
Figure 756930DEST_PATH_IMAGE034
步骤2.2.4：将一阶系统幅频特性表达式（1）变换为一阶系统对数幅频特性表达式（3），将式（3）和式（2）角频率取对数后得到近似一阶系统幅频响应曲线如图7所示，近似一阶系统相频响应曲线如图8所示。由一阶系统对数幅频特性表达式（3）得到工作频带ωc表达式（4），常用于评估强冲击传感器动态特性的10 dB工作频带由工作频带ωc表达式（4）解得结果如式（5）所示。
Figure 732976DEST_PATH_IMAGE036
Figure 273679DEST_PATH_IMAGE038
Figure 854833DEST_PATH_IMAGE040
步骤2.3：当归一化镜像数据曲线D近似为二阶系统时域曲线时，将归一化镜像数据曲线D各参数代入二阶系统幅频特性及相频特性表达式得到冲击传感器幅频特性和相频特性，进而获取工作频带得到冲击传感器的动态特性参数，实现强冲击传感器下阶跃响应动态特性计算。
步骤2.3.1：图9所示的归一化镜像数据曲线D幅值的瞬时最大偏差值与稳态值之比即为最大超调量σ，如式（6）所示。由阻尼比ζ与超调量σ的关系式（7）得到阻尼比如式（8）所示。过渡时间即调节时间，指系统受到扰动后回到稳定状态所经历的最短时间。由过渡时间ts、阻尼比ζ和固有频率ωn的关系式（9）得到固有频率如式（10）所示。
Figure 189999DEST_PATH_IMAGE042
其中：Ymax为瞬时最大偏差值；Y(∞)为稳态值。
Figure 653342DEST_PATH_IMAGE044
Figure 732156DEST_PATH_IMAGE046
Figure 230133DEST_PATH_IMAGE048
其中：系数b在误差带范围取±5%时为3.5，在误差带范围取±2%时为4.4。
Figure 736201DEST_PATH_IMAGE050
步骤2.3.2：将阻尼比ζ和固有频率ωn代入二阶系统归一化传递函数表达式（11），并将二阶系统归一化传递函数表达式（11）进行变换得到二阶系统频率响应函数表达式（12），根据二阶系统频率响应函数表达式（12）得到二阶系统归一化幅频特性表达式（13），利用二阶系统归一化幅频特性表达式（13）得到冲击传感器幅频特性。
Figure 90435DEST_PATH_IMAGE052
Figure 238519DEST_PATH_IMAGE054
其中：ω为角频率。
Figure 794265DEST_PATH_IMAGE056
步骤2.3.3：由二阶系统频率响应函数表达式（12）得到二阶系统相频特性表达式（14），利用二阶系统相频特性表达式（14）得到冲击传感器相频特性。
Figure 2393DEST_PATH_IMAGE058
步骤2.3.4：将二阶系统归一化幅频特性表达式（13）变换为二阶系统对数幅频特性表达式（15），将式（15）和式（4）角频率取对数后得到近似二阶系统幅频响应曲线如图10所示，近似二阶系统相频响应曲线如图11所示。由二阶系统对数幅频特性表达式（15）得到二阶系统工作频带ωc表达式（16），二阶系统工作频带ωc表达式（16）有解条件为式（17）所示，则可解得工作频带ωc如式（18）所示。常用于评估强冲击传感器动态特性的10 dB工作频带由式（17）和式（18）解得结果如式（19）和式（20）所示。
Figure 909169DEST_PATH_IMAGE060
Figure 533048DEST_PATH_IMAGE062
Figure 271197DEST_PATH_IMAGE064
Figure 322330DEST_PATH_IMAGE066
若
Figure 513140DEST_PATH_IMAGE068
，则工作频带：
Figure 206289DEST_PATH_IMAGE070
若
Figure 736628DEST_PATH_IMAGE072
，则工作频带：
Figure 286558DEST_PATH_IMAGE074
如图12所示，本实施例公开的一种强冲击传感器下阶跃激励发生装置，包括弹丸1、脉宽调节缓冲垫2、止挡盘3、砧体4和冲击传感器5。所述脉宽调节缓冲垫2材料选用泡沫铝。脉宽调节缓冲垫2放置在砧体4凹槽中。脉宽调节缓冲垫2和砧体4放置在止挡盘3筒内。止挡盘3与炮管固定连接。弹丸1、脉宽调节缓冲垫2、止挡盘3、砧体4和冲击传感器5沿空气炮炮管轴线方向同轴安装布置。如图13所示，弹丸1被压缩空气推动由弹舱位置沿着空气炮炮管加速向前运动，弹丸1和脉宽调节缓冲垫2碰撞过程中，弹丸1通过脉宽调节缓冲垫2推动砧体4一起加速运动，脉宽调节缓冲垫2结构不断塌陷压缩，该塌陷压缩过程为塑性变形无回弹。当脉宽调节缓冲垫2被压缩至极限，弹丸1被止挡盘3结构阻挡停止向前运动时，砧体4与弹丸1瞬间分离，使得与砧体4连接的冲击传感器5受冲击力突然消失为零，产生近似理想的下阶跃激励及响应。
以上所述的具体描述，对发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (3)

1.一种强冲击传感器下阶跃响应动态特性计算方法，其特征在于：对冲击传感器下阶跃响应信号进行滤波处理，截取所得曲线峰值点之后的部分进行数据归一化处理后再上下翻转曲线，产生近似一阶或二阶系统时域曲线，将时域曲线各参数代入一阶或二阶系统幅频特性及相频特性表达式得到冲击传感器幅频特性和相频特性，进而获取工作频带得到冲击传感器的动态特性参数，实现强冲击传感器下阶跃响应动态特性计算；
所述的一种强冲击传感器下阶跃响应动态特性计算方法，包括如下步骤，
步骤一，对冲击传感器下阶跃响应信号进行滤波处理，截取所得曲线峰值点之后的部分进行数据归一化处理后再上下翻转曲线得到冲击传感器下阶跃响应的归一化镜像数据曲线D；
步骤一实现方法为，
步骤1.1：通过滤波器对冲击传感器下阶跃响应信号进行滤波处理，得到冲击传感器下阶跃响应的滤波曲线A；
步骤1.2：截取所得滤波曲线A峰值点之后的时域信号数据，得到冲击传感器下阶跃响应的截断曲线B；
步骤1.3：以峰值点幅值为标准，将所得截断曲线B所有数据点幅值除以峰值点幅值进行归一化处理，得到冲击传感器下阶跃响应的归一化数据曲线C；
步骤1.4：对所得归一化数据曲线C以峰值点幅值的一半为对称轴上下翻转得到冲击传感器下阶跃响应的归一化镜像数据曲线D；
步骤二，步骤一所得归一化镜像数据曲线D近似为一阶或二阶系统时域曲线，将归一化镜像数据曲线D各参数代入一阶或二阶系统幅频特性及相频特性表达式得到冲击传感器幅频特性和相频特性，进而获取工作频带得到冲击传感器的动态特性参数，实现强冲击传感器下阶跃响应动态特性计算；
步骤二实现方法为，
步骤2.1：步骤一所得归一化镜像数据曲线D近似为一阶或二阶系统时域曲线；
步骤2.2：当归一化镜像数据曲线D近似为一阶系统时域曲线时，将归一化镜像数据曲线D各参数代入一阶系统幅频特性及相频特性表达式得到冲击传感器幅频特性和相频特性，进而获取工作频带得到冲击传感器的动态特性参数，实现强冲击传感器下阶跃响应动态特性计算；
步骤2.2.1：由归一化镜像数据曲线D得到时间常数τ；
步骤2.2.2：将时间常数τ代入一阶系统幅频特性表达式(1)，得到冲击传感器幅频特性；
Figure FDA0003462606790000011
其中：ω为角频率；
步骤2.2.3：将时间常数τ代入一阶系统相频特性表达式(2)，得到冲击传感器相频特性；
Figure FDA0003462606790000021
步骤2.2.4：将一阶系统幅频特性表达式(1)变换为一阶系统对数幅频特性表达式(3)，由一阶系统对数幅频特性表达式(3)得到工作频带ωc表达式(4)；
G(ω)＝20lg A(ω)＝-10lg(ω2τ2+1) (3)
Figure FDA0003462606790000022
步骤2.3：当归一化镜像数据曲线D近似为二阶系统时域曲线时，将归一化镜像数据曲线D各参数代入二阶系统幅频特性及相频特性表达式得到冲击传感器幅频特性和相频特性，进而获取工作频带得到冲击传感器的动态特性参数，实现强冲击传感器下阶跃响应动态特性计算；
步骤2.3.1：由归一化镜像数据曲线D得到最大超调量σ和过渡时间ts，将表达式(5)变换为表达式(6)，将最大超调量σ和过渡时间ts代入表达式(6)得到阻尼比ζ，并将阻尼比ζ代入表达式(7)得到固有频率ωn；
Figure FDA0003462606790000023
Figure FDA0003462606790000024
Figure FDA0003462606790000025
其中：系数b根据过渡时间误差带范围取值确定；
步骤2.3.2：将阻尼比ζ和固有频率ωn代入二阶系统归一化传递函数表达式(8)，并将二阶系统归一化传递函数表达式(8)进行变换得到二阶系统频率响应函数表达式(9)，根据二阶系统频率响应函数表达式(9)得到二阶系统归一化幅频特性表达式(10)，利用二阶系统归一化幅频特性表达式(10)得到冲击传感器幅频特性；
Figure FDA0003462606790000026
Figure FDA0003462606790000027
其中：ω为角频率；
Figure FDA0003462606790000031
步骤2.3.3：由二阶系统频率响应函数表达式(9)得到二阶系统相频特性表达式(11)，利用二阶系统相频特性表达式(11)得到冲击传感器相频特性；
Figure FDA0003462606790000032
步骤2.3.4：将二阶系统归一化幅频特性表达式(10)变换为二阶系统对数幅频特性表达式(12)，由二阶系统对数幅频特性表达式(12)得到二阶系统工作频带ωc表达式(13)，二阶系统工作频带ωc表达式(13)有解条件为式(14)所示，则可解得工作频带ωc如式(15)所示；
Figure FDA0003462606790000033
Figure FDA0003462606790000034
G(ω)≤-10lg(4ζ2-4ζ4) (14)
Figure FDA0003462606790000035
所述的一种强冲击传感器下阶跃响应动态特性计算方法所需要的下阶跃激励反应信号由强冲击传感器下阶跃激励发生装置产生，所述强冲击传感器下阶跃激励发生装置使用受冲击发生塑性变形无回弹的脉宽调节缓冲垫作为高量程冲击加载的惯性力传递介质；弹丸和脉宽调节缓冲垫碰撞过程中，弹丸通过脉宽调节缓冲垫推动砧体一起加速运动，脉宽调节缓冲垫结构不断塌陷压缩，该塌陷压缩过程为塑性变形无回弹；当脉宽调节缓冲垫结构被压缩至极限，弹丸被止挡盘结构阻挡停止向前运动时，砧体与弹丸瞬间分离，使得与砧体连接的冲击传感器所受冲击力突然消失为零，产生近似理想的下阶跃激励及响应。
2.如权利要求1所述的一种强冲击传感器下阶跃响应动态特性计算方法，其特征在于：所述脉宽调节缓冲垫材料选用泡沫铝。
3.如权利要求1所述的一种强冲击传感器下阶跃响应动态特性计算方法，其特征在于：强冲击传感器下阶跃激励发生装置包括弹丸、脉宽调节缓冲垫、止挡盘、砧体和冲击传感器；脉宽调节缓冲垫放置在砧体凹槽中；脉宽调节缓冲垫和砧体放置在止挡盘中心孔内；止挡盘与炮管固定连接；弹丸、脉宽调节缓冲垫、止挡盘、砧体和冲击传感器沿空气炮炮管轴线方向同轴安装布置；弹丸被压缩空气推动由弹舱位置沿着空气炮炮管加速向前运动，弹丸和脉宽调节缓冲垫碰撞过程中，弹丸通过脉宽调节缓冲垫推动砧体一起加速运动，脉宽调节缓冲垫结构不断塌陷压缩，该塌陷压缩过程为塑性变形无回弹；当脉宽调节缓冲垫被压缩至极限，弹丸被止挡盘结构阻挡停止向前运动时，砧体与弹丸瞬间分离，使得与砧体连接的冲击传感器所受冲击力突然消失为零，产生近似理想的下阶跃激励及响应。
CN202111382834.XA 2021-11-22 2021-11-22 强冲击传感器下阶跃响应动态特性计算方法及发生装置 Active CN113820067B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111382834.XA CN113820067B (zh) 2021-11-22 2021-11-22 强冲击传感器下阶跃响应动态特性计算方法及发生装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111382834.XA CN113820067B (zh) 2021-11-22 2021-11-22 强冲击传感器下阶跃响应动态特性计算方法及发生装置

Publications (2)

Publication Number Publication Date
CN113820067A CN113820067A (zh) 2021-12-21
CN113820067B true CN113820067B (zh) 2022-02-18

Family

ID=78917975

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111382834.XA Active CN113820067B (zh) 2021-11-22 2021-11-22 强冲击传感器下阶跃响应动态特性计算方法及发生装置

Country Status (1)

Country Link
CN (1) CN113820067B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101458263A (zh) * 2009-01-08 2009-06-17 中北大学 高冲击加速度传感器动态特性的溯源校准方法
US20090319251A1 (en) * 2008-06-23 2009-12-24 Choshu Ito Circuit Simulation Using Step Response Analysis in the Frequency Domain
CN102788613A (zh) * 2012-07-30 2012-11-21 中国船舶重工集团公司第七○二研究所 实船水下爆炸强冲击环境的分频段测量方法与装置
CN106463136A (zh) * 2014-06-26 2017-02-22 高通股份有限公司 基于高频带信号特征的时间增益调整

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07239711A (ja) * 1994-03-01 1995-09-12 Fuji Electric Co Ltd 過渡応答波形特徴量算出方法
US20030035376A1 (en) * 2001-08-20 2003-02-20 Xiaofen Chen Derivation of composite step-function response

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090319251A1 (en) * 2008-06-23 2009-12-24 Choshu Ito Circuit Simulation Using Step Response Analysis in the Frequency Domain
CN101458263A (zh) * 2009-01-08 2009-06-17 中北大学 高冲击加速度传感器动态特性的溯源校准方法
CN102788613A (zh) * 2012-07-30 2012-11-21 中国船舶重工集团公司第七○二研究所 实船水下爆炸强冲击环境的分频段测量方法与装置
CN106463136A (zh) * 2014-06-26 2017-02-22 高通股份有限公司 基于高频带信号特征的时间增益调整

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于激波管校准的冲击波压力传感器动态特性研究;杨凡等;《南京理工大学学报》;20170630;第41卷(第3期);全文 *
高g值加速度计高冲击校准技术综述;林然等;《探测与控制学报》;20150831;第37卷(第04期);第110页 *

Also Published As

Publication number Publication date
CN113820067A (zh) 2021-12-21

Similar Documents

Publication Publication Date Title
Joseph et al. Near field zones of quiet
JP3250801B2 (ja) ディジタルフィルタに使用するための係数を生成するためのシステム
US10789932B2 (en) Noise control
Surace et al. On the non-linear characteristics of automotive shock absorbers
CN113820067B (zh) 强冲击传感器下阶跃响应动态特性计算方法及发生装置
Kostek et al. Hybrid noise control in ducts
McInerny et al. High-intensity rocket noise: Nonlinear propagation, atmospheric absorption, and characterization
CN109311380B (zh) 用于混合动力电动车辆不点火检测的响应幅度修正
CN111211726B (zh) 用于生成马达驱动信号的系统
JP2020122734A (ja) 振動制御装置
CN111628848B (zh) 通信符号的检测方法、设备及计算机可读存储介质
US3420098A (en) Transient synthesizing system
FR3064066A1 (fr) Banc d'essai pour generer des charges dynamiques sur un actionneur et methode pour realiser des essais sur un actionneur en utilisant ledit banc d'essai
CN106197917B (zh) 加速度时域截断信号波形再现的处理方法
JP4411034B2 (ja) 多重飛翔体を用いた衝撃加速度発生による加速度センサの校正評価装置とその方法
US20240077632A1 (en) Seismic acquisition system based on frequency domain expansion mems sensor
JP3371477B2 (ja) 車両用能動型振動制御装置及び車両用能動型騒音制御装置
US10533911B2 (en) Signal processing device, air pressure sensor assembly and electronics apparatus
KR102096002B1 (ko) 데이터 축약 통합 계측 장치, 이의 제어 방법, 그리고 이 방법을 저장한 컴퓨터 판독 가능 저장 매체
Wisniewiski Dynamic calibration of pressure transducers with an improved shock tube system
Gabrielson Measurement of Blast Waveforms with Condenser Microphones: Measurement Interpretation and Correction
şi Mecatronică Experimental study of an electromechanical system used to control the mechanical mobility
US11835426B2 (en) Vibration control system using kurtosis response spectrum
KR102121191B1 (ko) 2관성계 구동기 제어 장치
Mazzola et al. On the generation of transient acoustic pulses in water

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