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

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

Info

Publication number
CN113820067A
CN113820067A CN202111382834.XA CN202111382834A CN113820067A CN 113820067 A CN113820067 A CN 113820067A CN 202111382834 A CN202111382834 A CN 202111382834A CN 113820067 A CN113820067 A CN 113820067A
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.)
Granted
Application number
CN202111382834.XA
Other languages
English (en)
Other versions
CN113820067B (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 (8)

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

* 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 過渡応答波形特徴量算出方法
JP2003083738A (ja) * 2001-08-20 2003-03-19 Tektronix Inc ステップ関数応答導出方法
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 高通股份有限公司 基于高频带信号特征的时间增益调整

Patent Citations (6)

* 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 過渡応答波形特徴量算出方法
JP2003083738A (ja) * 2001-08-20 2003-03-19 Tektronix Inc ステップ関数応答導出方法
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
杨凡等: "基于激波管校准的冲击波压力传感器动态特性研究", 《南京理工大学学报》 *
林然等: "高g值加速度计高冲击校准技术综述", 《探测与控制学报》 *

Also Published As

Publication number Publication date
CN113820067B (zh) 2022-02-18

Similar Documents

Publication Publication Date Title
Joseph et al. Near field zones of quiet
JP3250801B2 (ja) ディジタルフィルタに使用するための係数を生成するためのシステム
CN104935293A (zh) 大功率变压器自适应有源降噪控制方法及控制系统
CN106199185B (zh) 一种基于连续对数扫频的线性脉冲响应测量方法及系统
CN113820067B (zh) 强冲击传感器下阶跃响应动态特性计算方法及发生装置
CN109186693B (zh) 一种自适应超声波的回波信号检测电路
Kostek et al. Hybrid noise control in ducts
WO2020158068A1 (ja) 振動制御装置
Tell et al. Acoustic waves in granular packings at low confinement pressure
CN111211726A (zh) 用于生成马达驱动信号的系统
CN111697952B (zh) 一种基于数字pzc系统调节脉冲宽度的方法及系统
CN106706190B (zh) 一种确定冲击载荷等效静载的方法
KR100727488B1 (ko) 충격 응답 스펙트럼이 나타나는 시간들의 차이가 짧은충격파형의 합성 방법 및 충격 응답 내역을 구하는 디지털필터와 그 응용
CN107084222B (zh) 一种航炮缓冲控制方法
JP7476150B2 (ja) 診断方法、診断装置、及び診断プログラム
Kim et al. Robust broadband vibration control of a flexible structure using an electrical dynamic absorber
CN111931288B (zh) 一种火工冲击响应时频域分布相似性测度方法
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) 加速度时域截断信号波形再现的处理方法
Li et al. Dual-driver standing wave tube: acoustic impedance matching with robust repetitive control
JP4411034B2 (ja) 多重飛翔体を用いた衝撃加速度発生による加速度センサの校正評価装置とその方法
JP3371477B2 (ja) 車両用能動型振動制御装置及び車両用能動型騒音制御装置
Gabrielson Measurement of Blast Waveforms with Condenser Microphones: Measurement Interpretation and Correction
Wisniewiski Dynamic calibration of pressure transducers with an improved shock tube system
US20240017653A1 (en) Semi active suspension control methods using skyhook and endstop control

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