CN106197807B - 一种用于动态力的测量方法 - Google Patents
一种用于动态力的测量方法 Download PDFInfo
- Publication number
- CN106197807B CN106197807B CN201610669413.8A CN201610669413A CN106197807B CN 106197807 B CN106197807 B CN 106197807B CN 201610669413 A CN201610669413 A CN 201610669413A CN 106197807 B CN106197807 B CN 106197807B
- Authority
- CN
- China
- Prior art keywords
- signal
- dynamic
- frequency
- function
- dynamic force
- 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.)
- Expired - Fee Related
Links
- 238000000691 measurement method Methods 0.000 title claims abstract description 22
- 238000005259 measurement Methods 0.000 claims abstract description 39
- 230000005540 biological transmission Effects 0.000 claims abstract description 27
- 238000012545 processing Methods 0.000 claims abstract description 18
- 238000005516 engineering process Methods 0.000 claims abstract description 7
- 238000012546 transfer Methods 0.000 claims description 11
- 230000004044 response Effects 0.000 claims description 10
- 230000003068 static effect Effects 0.000 claims description 8
- 238000010183 spectrum analysis Methods 0.000 claims description 6
- 238000005314 correlation function Methods 0.000 claims description 5
- 238000005316 response function Methods 0.000 claims description 5
- 238000011068 loading method Methods 0.000 claims description 4
- 230000000737 periodic effect Effects 0.000 claims description 4
- 108010076504 Protein Sorting Signals Proteins 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000004880 explosion Methods 0.000 claims description 3
- 239000003721 gunpowder Substances 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 2
- 230000009191 jumping Effects 0.000 claims description 2
- 238000003672 processing method Methods 0.000 claims description 2
- 238000000034 method Methods 0.000 description 17
- 230000008859 change Effects 0.000 description 13
- 238000004458 analytical method Methods 0.000 description 8
- 238000013016 damping Methods 0.000 description 8
- 238000006073 displacement reaction Methods 0.000 description 8
- 238000012360 testing method Methods 0.000 description 8
- 238000002474 experimental method Methods 0.000 description 6
- 230000001133 acceleration Effects 0.000 description 4
- 238000005311 autocorrelation function Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000009499 grossing Methods 0.000 description 4
- 230000033001 locomotion Effects 0.000 description 4
- 241001269238 Data Species 0.000 description 3
- 230000003321 amplification Effects 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 238000003199 nucleic acid amplification method Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000010998 test method Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 241000184339 Nemophila maculata Species 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000002737 fuel gas Substances 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000035935 pregnancy Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 239000007921 spray Substances 0.000 description 1
- 238000005507 spraying Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L23/00—Devices or apparatus for measuring or indicating or recording rapid changes, such as oscillations, in the pressure of steam, gas, or liquid; Indicators for determining work or energy of steam, internal-combustion, or other fluid-pressure engines from the condition of the working fluid
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L23/00—Devices or apparatus for measuring or indicating or recording rapid changes, such as oscillations, in the pressure of steam, gas, or liquid; Indicators for determining work or energy of steam, internal-combustion, or other fluid-pressure engines from the condition of the working fluid
- G01L23/22—Devices or apparatus for measuring or indicating or recording rapid changes, such as oscillations, in the pressure of steam, gas, or liquid; Indicators for determining work or energy of steam, internal-combustion, or other fluid-pressure engines from the condition of the working fluid for detecting or indicating knocks in internal-combustion engines; Units comprising pressure-sensitive members combined with ignitors for firing internal-combustion engines
Landscapes
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Combustion & Propulsion (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种用于动态力的测量方法,属于力学和信号处理技术领域。所述测量方法通过动态力测量***收集输入信号、输出信号以及其他信号,利用绘图软件绘制出曲线,初步判断各信号的,据此获取信号特征值和特征函数;推导输入信号和输出信号的传递函数关系;应用传递函数,推算出输出信号或输入信号。本发明适用范围广、易于操作、结果可靠;本发明适合于不能直接测量动态力场合,所得结果的准确性高。
Description
技术领域
本发明属于力学和信号处理技术领域,具体涉及一种用于动态力的测量方法。
背景技术
动态力指零件或部件等物体在进行加速或减速运动(直线或旋转运动)、单自由度物体受外力作强迫振动时产生的力,以及物体作加速或减速旋转运动时产生的力或旋转力矩。力的测量分为静态和动态两种形式。其中,静态测量指测量期间其值可认为是恒定的量的测量;而动态测量指为确定量的瞬时值及其随时间变化所进行的测量,即被测量是随时间而变化的。动态力具有各种各样的形式表现,例如,踏跳力(指运动员起跳瞬间对地面的作用力)、脉冲爆震发动机推力(这种发动机以周期的、脉动的方式产生推力)、姿态控制火箭发动机改变姿态时产生的动态推力、轨道控制火箭发动机在轨道变换时产生的动态推力(后文简写为,姿/轨控火箭发动机)、减震器上受的载荷等等。
一个完整的测试***应包括以下几部分:被测参数的检测部分、信息转换部分和参数的显示部分。测试***的检测部分直接与被测对象***,直接感受被测参数的变化,并把被测参数变换为易于传输、处理、测量的信号,通常称完成这一任务的装置(或器件)为传感器(或一次变换器)。信号的放大、变换、运算部分是将传感器的输出信号进行再变换、放大和衰减、调制和解调、阻抗变换、分析和运算等项处理,使信号变成适于需要而又便于记录和显示的信号,它是联系传感器和参数显示部分的中间环节,称为中间变换器(二次变换器或测量电路)。测试***的最后部分是参数显示装置,即记录器、显示器,又称为输出变换器(胡寿松.自动控制原理(胡寿松+第五版)[M].)。因此,一个测试***可由三个基本环节组成,即传感器、中间变换器、记录显示仪器(李永新,吴健.信号与动态测量***[M].人民邮电出版社,2014.)。
现有技术中常用的测量方法为弹簧-质量-阻尼***法(王济.matlab在振动信号处理中的应用[M].中国水利水电出版社,知识产权出版社,2005.),如图1所示原理图,弹簧2将动架和力产生装置3连接在推力测试台静架,同时设置阻尼1,使其形成一个弹簧-质量-阻尼***,动架和力产生装置3工作时,通过测量该振动响应来推算出动态力4,并进一步求得平均力。动架和力产生装置3的底部和两侧安装有滚动轴承,以保证力产生装置工作时整个动架和力产生装置3只能在轴向产生振动。这个***可看为一个单自由度的弹簧-质量-阻尼***。依据简化的原理图,由牛顿运动定律列写,推力求解表达式,
令
上式中,ω为***的角频率;ξ为阻尼比。这两个特性参数可以通过***识别的方法获得。
F(t)为发动机推力;
md+me为动架、发动机的质量和;
C为阻尼系数;
K为弹簧刚度系数;
为质量块的振动响应。
在已知质量时,力产生装置的瞬时推力和平均推力可以通过测量质量块的振动响应求得。
注意事项:一、选用的传感器频响范围要宽;二、对某一时刻各个振动响应的测量应同步进行,这样才能保证计算时的正确性和所测推力信号的不失真。
现有动态力测量方法的不足:
1.脉动式的产生动态力,导致力产生装置和动架部分,突然在轴向产生振动,因间断产生力时间极短,又接着导致力产生装置和动架连续的振动,动态力作用于弹簧上,产生轴向位移,因此用位移一次求导得到速度,二次求导得到加速度,再求出动态力,动态力因克服自身惯性力而小于实际动态力值。
2.据现有技术方案中所述,对振动***在模拟多循环脉冲力作用下的强迫振动响应进行动态仿真。选取不同的弹簧刚度系数、阻尼系数,获得了振动***的位移、速度、加速度响应曲线。理论分析和数值仿真都发现,根据振动***位移响应平均值求出的平均力,即为力产生装置的平均力。发现弹簧刚度系数过大,会导致振动***平均位移的减小,不利于准确测量;反之,弹簧刚度系数过小则会引起振动***位移量峰值的增大,很难选择到合适量程范围的位移传感器。现有技术方案中选取合适刚度的弹簧,没有明确可靠的方法或者经验。文章通过多次仿真试验与个人经验来选择弹簧。因此,试验准备阶段工作复杂,可操作性差。
3.现有技术方案中,选取适合刚度弹簧的做法是,实验前预估力产生装置的力值范围并多次试验不同刚度的弹簧,以求试验得到使位移传感器产生位移大小适中的弹簧。如果换一台不同力值范围的力产生装置意味着,需要重新选取合适刚度的弹簧,相对于其他推力测量方法,这会导致试验程序过于繁琐。
发明内容
本发明为了解决现有技术中存在的问题,提供一种用于动态力的测量方法,对现有的动态力测量方法进行改进,可以高精度、快响应的测量动态力。
本发明提供的用于动态力的测量方法基于如下两点假设:
假设1、动态力测量***经过校准或者动、静态标定且符合技术要求。
假设2、动态力测量***的数据采集***数据量足够且精度符合要求。
所述的用于动态力的测量方法,具体包括如下步骤:
步骤1:通过动态力测量***收集输入信号和输出信号。输出信号即为所述的需要测量的动态力。
步骤2:将收集到的输入信号、输出信号,利用绘图软件绘制出曲线,初步判断各信号的类型。
步骤3:根据信号的类型,获取信号特征值和特征函数,如峰值、均值、功率、能量、平均功率、方均根值等特征值以及频谱、概率密度函数等。
步骤4:信号处理,推导输入信号和输出信号的传递函数关系。
所述信号处理具体实施流程如下:
首先,对信号进行预处理,例如,去趋势项、数据平滑、滤波等等。
其次,信号的类型须是稳态确定性的,针对周期信号和非周期信号的处理方法不同。周期信号频谱分析采用离散傅里叶级数分析,非周期信号频谱分析采用离散傅里叶变换。先辨识求解***(或理论分析)的传递函数,接着求解传递函数的零极点并绘制测量***的零极点分布图,判断动态力测量***稳定性。求解出动态力测量***的频率响应函数的幅频特性(线性增益或分贝增益)、相频特性,估计动态力测量***的工作频带,定量评价动态力测量***的动态测量能力。
步骤5:应用传递函数,推算出输出信号或输入信号。
使用动态力测量***测量输入信号,应用步骤4中的传递函数,得出输出信号即需要测量的动态力,具体为:
首先,根据步骤4中传递函数,判断动态力测量***的传递函数是否可以分解为典型***(理想放大***、典型一阶***、典型二阶***、理想微分***、理想积分***、一阶微分***和二阶微分***等)的串联组合,以简化应用工作;
其次,对于不同组试验,输入信号作为输入量利用传递函数求解到输出信号并与实测输出信号对比。
本发明测量方法中所述的动态力,包括运动员的踏跳力、动力机械(分为,风力机械、水力机械和热力发动机三大类)承受的动载荷或产生的动态力(比如,风力机的叶片承受的动载荷、水轮机发电时叶片承受的动载荷、以热力发动机为动力源的机械装置产生的动态力)和火药***时产生的冲击力。
本发明的优点在于:
(1)分解传递函数,将复杂函数简单化,成熟的典型***的特性理论易于理解,加快数据处理进程。软件仿真能在一定程度上验证***传递函数的准确性。
(2)本发明适用范围广、易于操作、结果可靠。
(3)本发明方法适合于不能直接测量动态力场合,例如,运动员的踏跳力,可运用本发明方法,测量运动员作用接触面上的压力反推踏跳力;以热力发动机为动力源的机械装置产生的动态力(例如,脉冲爆震发动机的推力测量,可以通过测量发动机的燃气压力来近似求解推力;姿态/轨道控制火箭发动机在调整姿态或者变换轨道时产生的动态推力,可以通过测量燃气压力来推算瞬态推力);火药***时产生的冲击,可通过测量瞬时气压压力来测算冲击力。
(4)本发明方法的应用所得结果能作为参考标准。本发明方法对数据进行了有效化处理,所得结果的准确性高。
附图说明
图1为现有技术中的弹簧-质量-阻尼***法原理图。
图2为本发明方法总体步骤流程图。
图3为姿态/轨道控制火箭发动机动态推力测量***原理图。
图4为***的传递函数形象图。
图5为本发明的测量方法中信号处理和推导关系流程图。
图中:
1.阻尼; 2.弹簧; 3.动架和力产生装置; 4.动态力; 5.静架;
6.动态力传感器; 7.力产生装置; 8.动架; 9.测量输入信号的传感器。
具体实施方式
下面以姿态/轨道控制火箭发动机动态推力测量(后文简写为,姿轨控火箭发动机)为实例说明应用本发明动态力的测量方法,不同领域的动态力测量需灵活运用本发明方法。
在本实例中,姿态/轨道控制火箭发动机的动态推力测量***经静态和动态标定过且在误差范围内,静态标定采用砝码增减法或机械螺旋加载方式,动态标定采用力锤锤击法。动态推力测量***经过校准或者动、静态标定且符合技术要求。动态推力测量***的数据采集***数据量足够且精度符合要求。
姿轨控火箭发动机动态推力测量***如图3所示,所述动态推力测量***包括静架5、动态力传感器6、力产生装置7、动架8、测量输入信号的传感器9。所述的动态力传感器6设置在静架5和力产生装置7之间,动架8设置在静架5上,力产生装置7位于动架8上,在所述的力产生装置7上设置有测量输入信号的传感器9。
结合图2所示流程,应用本发明提供的用于动态力的测量方法,对姿轨控火箭发动机动态推力进行测量,具体包括如下步骤:
步骤1:通过动态推力测量***收集输入、输出信号。
姿轨控火箭发动机动态推力测量***测量得到两个物理参数:推力(输出信号)和压力(输入信号)。
在保证测量精度与采样数据足够的前提下,测量动态推力采用压电式动态力传感器,测量压力采用压电式动态压力传感器。
所述的压电式动态力传感器位于发动机封闭端,测量轴向主矢推力,压电式动态压力传感器位于喷管出口,测量发动机内部压力。
发动机工作时,燃烧室的高温、高压燃气是产生推力的主要因素。因此,为了获得发动机的推力(输出信号),可将发动机内部压力作为输入信号。
姿轨控火箭发动机做地面试验时,可以近似看成单自由度的振动***。这时发动机只有轴向一个自由度,发动机推力主要沿着轴向产生,因此压电式动态力传感器安置于发动机封闭端,这样测量简单高效。
发动机工作时,发动机内部压力急剧变化,高温燃气喷出喷管时产生推力,测量喷管出口部位的压力变化,有利于后续探索压力-推力的关系。
步骤2:将测量收集到的输入信号和输出信号,利用绘图软件绘制出曲线,初步判断信号的类型。
初步判断信号的类型,具体实施流程如下:
步骤2.1:判定动态推力测量***的类型;
从不同角度来判断,如果信号是离散时间信号,则动态推力测量***就是离散时间***;如果是连续时间信号,则动态推力测量***就是连续时间***。如果所测物理量之间的约束关系全部是线性(不仅仅是简单的比例关系)的,则动态推力测量***是线性***,反之,所测物理量之间不全是线性关系,则动态推力测量***是非线性***。如果动态推力测量***的结构参数保持常量,不随时间变化,则动态推力测量***为定常***(时不变***);反之,则动态推力测量***为时变***。动态推力测量***类型主要从以上三方面来判断。
一般而言,如果信号为数值形式存储,则动态推力测量***为离散时间***;动态推力测量***的结构/性能参数比如,各阶固有频率、相应的阻尼比、振型、响应时间、采样周期及采样数据精度等等几乎不随时间变化的,故认为动态推力测量***为定常***;测量输入物理量与输出物理量之间存在某种确定的函数关系(一般精确的通用的数学表达式无法理论推导得出),认为动态推力测量***是线性的。总之,当前的推力测量***(动态推力/静态推力)绝大部分是离散时间线性定常***。
步骤2.2:根据信号曲线,先判断信号稳态性与否,接着判断信号确定性与否;
步骤2.3:判断确定性信号的周期性与否;
将不同试验条件下或同试验条件下的多组压力、推力数据分别在两个坐标系(时间为自变量)绘制出;将同试验条件下同一组压力和推力数据同坐标系显示。观察同一物理量的不同组实验数据的峰值、周期、噪声、上升沿、衰减等特征,得出信号类型。
在实例中,观察同一物理量(压力、推力)的不同组实验数据的峰值、周期、噪声、上升沿、衰减等特征。结果显示,同一试验条件下,同一物理量的不同组试验数据,峰值趋同、周期相近、上升沿变化趋势一致、衰减趋势相同、噪声情况相仿。不同试验条件下,不同物理量的变化趋势基本一致,压力、推力的变化曲线呈周期性、脉冲式、负指数衰减的变化特征。结合信号的分类,初步判定两种物理量信号是稳态确定性周期性信号。
上述处理过程中,将数值形式的信号以曲线的形式显示,便于发现规律和判别信号类型。判定信号的类型,便于后续采取相应的信号分析方法获取信号的幅值相频特性。
步骤3:根据信号的类型,利用数据分析软件显示并分析输入信号和输出信号。
具体实施流程:
步骤3.1:使用数据分析软件(Matlab、Tecplot、Excel、Python等)显示所有信号,每种信号单独显示,输入和输出信号同坐标显示,信号曲线均以时间为自变量。再次详细获取信号数据的变化情况,包括时域的上升时间、过冲和异常值等衡量信号数据的动态特性参数。
步骤3.2:使用数据分析软件求出输入、输出等信号的特征统计量和进行频谱分析,初步研判频域的有效频带、截止频率等,又例如,峰值、有效值、平均绝对值和高、低频特性等。
获取时域的信号动态特性,可以初步判断动态推力测量***的动态测量性能(响应快慢、测量精度两方面)以及了解信号的特征参数的量级。
获取上升时间、衰减趋势与稳态状况以了解信号的瞬态变化所处的量级、预判衰减函数刻画与得出极值的量级。判断有无过冲、趋势项、异常值等,以了解动态推力测量***固有误差或者随机动态误差,以便于下一步骤采取相应的信号处理措施(消除趋势项、剔除异常值等措施)。
输入信号与输出信号同坐标显示,以获得幅值、时延的情况,判断动态推力测量***响应质量、输入和输出信号的相关程度,以研判在下一步骤中是否需要采取信号相关分析定量的考察关系密切程度。
通过傅里叶变换对,将信号从时域变换到频域分析,以了解信号的频域范围、不同频率的幅值大小与相位变化情况、有效频带的范围。
实例中,压力和推力均以时间为自变量在两个坐标系分开绘制变化曲线。另外将压力和推力曲线在同一个坐标系即第三个坐标系中显示;
求出压力、推力的特征统计量并分析频谱特性。
步骤4:处理输入信号和输出信号,推导输入信号和输出信号的关系。
输入信号和输出信号的关系,即传递函数对动态推力测量***的形象化见图4,输入信号x(t)经过动态推力测量***得到输出信号y(t),用h(t-k)来表征动态推力测量***因输入信号的激励所产生的响应。
实例中,已知压力、推力信号数据(数字信号、离散时间、离散量化数值),利用数据处理软件分别求出输入和输出信号序列的z变换,然后用输出信号z函数比输入信号z函数得到动态推力测量***的z传递函数。将z传递函数变换为零极点形式,求出零极点。做出零极点图,判断动态推力测量***的稳定性。稳定动态力推力测量***的z传递函数,可令z=ejω,ω表角频率,求出动态推力测量***的频率响应函数和幅相特性曲线。
用公式表达上述的处理流程为:
已知:压力p[n1]、推力f[n2],两种物理信号均是离散信号序列,n1、n2分别是每种信号数据的序列号;
未知(需求解量):传递函数H(z)(z可为任意的复数,z域是含有复数的频域)。
结合图5所示流程,所述的处理输入信号、输出信号,具体方法步骤为:
步骤4.1,信号的预处理:
对输入信号、输出信号进行去趋势项、数据平滑、滤波等处理。
对于步骤3中的输入信号和输出信号,如果信号因偏离基线而存在波形畸变,则需要进行去趋势项处理(常用多项式最小二乘法);如果信号存在的噪声,则需要对信号进行数据平滑处理(平均法、五点三次平滑法)。通过滤波过滤信号次要甚至无用频段(脉冲响应信号一般选用指数窗,脉冲信号推荐用矩形窗)。
步骤4.2,对离散信号做单边z变换,将离散时间信号分解为离散复指数的线性组合(复变函数的罗朗级数展开发现:离散时间信号一般都可以分解为xa[n]=zn的线性组合),求解线性时不变离散时间***的传递函数H(z)。
将压力p[n1]表示为zn线性组合:
其中:压力信号的单边z变换为
推力f[n2]的单边z变换同上,用FS(z)表示,n2∈[0,+∞)
可得,和
bk、ak指不随时间变化的***参数,M、N分别是分子项和分母项z变量多项式的阶数。
仿真一个可实现的、稳定连续时间***(其数学特征包括:t<0时,h(t)=0,limt→+∞h(t)=0)的实用离散时间***的传递函数H(z)的定义域一般都为:包含|z|=1的圆外开区间:|z|>r1,r1<1,r1是小于1的数值。
动态推力测量***内部的结构、原理不清楚时(看成黑箱),根据输入和输出信号,建立一种能衡量动态推力测量***动态特性的函数关系,并依据这种关系来推算任意输入情况下的输出量。这种动态力的测量方法,简便高效。
求解动态推力测量***的传递函数,据此判断动态推力测量***的稳定性。根据***频率响应函数得出的幅相特性,来定量评价***的动态响应性能。
步骤4.3,求解动态推力测量***传递函数的零极点,传递函数变换成零极点形式,作出零极点分布图;
变换离散时间***的传递函数表达式为零极点(零点:使H(z)的值为零的值,以zk表示,极点:使H(z)的值为∞的值,以pk表示。)形式:
其中,a0、b0表示不随时间变化的z的零阶***参数;M、N分别表示分子多项式、分母多项式的阶数,ck、dk分别指分子项、分母项提取z变量后的系数。
绘制零极点图,判断动态推力测量***的稳定性(一个线性时不变***,当且仅当它的传递函数H(z)的收敛域包括单位圆|z|=1时,该***就是稳定的)。
步骤4.4,求离散频率响应函数,***的幅频、相频特性函数或图像:
对于定义域包含|z|=1的H(z),取其中TS是动态推力测量***(离散时间***)的取样时间,ω表示角频率,得到离散频率响应函数:
ωS是采样角频率。定义,Ω=ωTS为无量纲化(角)频率,取值范围是[0,π]。
由于取样时间TS量级上一般足够小,故离散频率响应函数认为是连续时不变***的频响函数。
幅频函数:
相频函数:
绘制出相应的幅相频曲线。
步骤4.5,求出输入、输出信号的自相关和互相关函数:
根据步骤3判断的信号类型及信号关联性的预判,采用不同的函数定量的衡量信号的相关性。
自相关函数:
对于能量信号x(t)(时限信号)的自相关函数Rxx(τ),
对于功率信号(大部分的周期信号)的自相关函数,
互相关函数Ryx(τ):
对于能量信号x(t)与y(t),
对于功率信号,
τ指时间延迟,T指信号。
该步骤中,对周期或非周期的信号进行离散时间傅里叶级数或离散时间傅里叶变换进行频谱分析,是为了了解信号的频率成分、判断信号的高低频属性和噪声频率成分。
求解自相关函数,可以知道信号x(t)或x[n]与信号时延后x(t-τ)或x[n-k]的结果呈线性比例关系的程度,能有效的描述信号波形前后关联性。互相关函数能描述一种信号与另一种信号的相似程度,从而客观的评价测量***的动态特性。
在实例中,利用数据处理软件对压力信号、推力信号(均为功率信号)分别求取自相关和互相关函数。
步骤5:应用传递函数,推算出推力或者压力。
在实例中,测量姿轨控火箭发动机的喷管出口的压力,压力作为输入信号,应用传递函数,得出推力。
将传递函数各项分解为经典测量***的串联组合,使处理数据步骤简化。进一步的使用仿真模块,模拟输入信号,观测输出信号的波形变化。
Claims (5)
1.一种用于动态力的测量方法,其特征在于:
假设1、动态力测量***经过校准或者动、静态标定且符合技术要求;
假设2、动态力测量***的数据采集***数据量足够且精度符合要求;
所述的用于动态力的测量方法,具体包括如下步骤:
步骤1:通过动态力测量***收集输入信号和输出信号;
步骤2:将收集到的输入信号、输出信号,利用绘图软件绘制出曲线,初步判断各信号的类型;
步骤3:根据信号的类型,获取信号特征值和特征函数;
步骤4:信号处理,推导输入信号和输出信号的传递函数关系;
所述信号处理具体实施流程如下:
首先,对信号进行预处理;
其次,信号的类型须是稳态确定性的,针对周期信号和非周期信号的处理方法不同;周期信号频谱分析采用离散傅里叶级数分析,非周期信号频谱分析采用离散傅里叶变换;先辨识求解传递函数,接着求解传递函数的零极点并绘制测量***的零极点分布图,判断动态力测量***稳定性;求解出动态力测量***的频率响应函数的幅频特性、相频特性,估计动态力测量***的工作频带,定量评价动态力测量***的动态测量能力;
具体地,所述的信号处理的步骤如下:
已知:压力p[n1]、推力f[n2],两种物理信号均是离散信号序列,n1、n2分别是每种信号数据的序列号;
未知:传递函数H(z),z为任意的复数,z域是含有复数的频域;
步骤4.1,信号的预处理:
步骤4.2,对离散信号做单边z变换,将离散时间信号分解为离散复指数的线性组合,求解线性时不变离散时间***的传递函数H(z);
将压力p[n1]表示为zn线性组合:
其中:压力信号的单边z变换为
推力信号f[n2]的单边z变换同上,用FS(z)表示,
得,和所以有传递函数如下:
bk、ak指不随时间变化的***参数,M、N分别是分子项和分母项z变量多项式的阶数;
步骤4.3,求解动态推力测量***传递函数的零极点,传递函数变换成零极点形式,作出零极点分布图;
变换离散时间***的传递函数表达式为零极点形式:
其中,a0、b0表示不随时间变化的z的零阶***参数;零点是指使H(z)的值为零的值,以zk表示,极点是指使H(z)的值为∞的值,以pk表示;M、N分别表示分子多项式、分母多项式的阶数,ck、dk分别指分子项、分母项提取z变量后的系数;
步骤4.4,求离散频率响应函数,***的幅频、相频特性函数或图像:
对于定义域包含|z|=1的H(z),取其中TS是动态推力测量***的取样时间,ω表示角频率,得到离散频率响应函数:
ωs是采样角频率;定义,Ω=ωTS为无量纲化频率,取值范围是[0,π];
由于取样时间TS量级上足够小,故离散频率响应函数认为是连续时不变***的频响函数;
幅频函数:
相频函数:
绘制出相应的幅相频曲线;
步骤4.5,求出输入、输出信号的自相关和互相关函数;步骤5:应用传递函数,推算出输出信号或输入信号。
2.根据权利要求1所述的一种用于动态力的测量方法,其特征在于:所述的动态力,包括运动员的踏跳力、动力机械承受的动载荷或产生的动态力和火药***时产生的冲击力。
3.根据权利要求1所述的一种用于动态力的测量方法,其特征在于:所述的判断信号的类型包括判断信号的离散连续性、稳定性、确定性和周期性。
4.根据权利要求1所述的一种用于动态力的测量方法,其特征在于:所述的步骤3具体为,
步骤3.1:使用数据分析软件显示所有信号,每种信号单独显示,输入和输出信号同坐标显示,信号曲线均以时间为自变量;再次详细获取信号数据的时域的上升时间、过冲和异常值;
步骤3.2:使用数据分析软件求出输入、输出信号的特征统计量和进行频谱分析,初步研判频域的有效频带和截止频率。
5.根据权利要求1所述的一种用于动态力的测量方法,其特征在于:步骤5中,先将动态力测量***的传递函数分解为经典测量***的串联组合,再进行推算、仿真和验证。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610669413.8A CN106197807B (zh) | 2016-08-15 | 2016-08-15 | 一种用于动态力的测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610669413.8A CN106197807B (zh) | 2016-08-15 | 2016-08-15 | 一种用于动态力的测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106197807A CN106197807A (zh) | 2016-12-07 |
CN106197807B true CN106197807B (zh) | 2018-10-16 |
Family
ID=57521493
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610669413.8A Expired - Fee Related CN106197807B (zh) | 2016-08-15 | 2016-08-15 | 一种用于动态力的测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106197807B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103622015A (zh) * | 2013-11-28 | 2014-03-12 | 江西仙客来生物科技有限公司 | 灵芝孢子油ve软胶囊 |
CN109060239A (zh) * | 2018-06-15 | 2018-12-21 | 西北工业大学 | 一种脉冲爆震发动机试车台推力测试动态标定方法和装置 |
CN110823346B (zh) * | 2018-08-10 | 2021-05-18 | 河南工业大学 | 基于底面单圈标准差指数模型的粮仓检测方法及*** |
CN110595683B (zh) * | 2019-08-12 | 2021-03-09 | 中北大学 | 一种冲击波测试***的低频补偿方法 |
CN111776247A (zh) * | 2020-06-04 | 2020-10-16 | 核工业北京地质研究院 | 一种滑翔无人机发射***推力测量***及方法 |
CN113091795B (zh) | 2021-03-29 | 2023-02-28 | 上海橙科微电子科技有限公司 | 光电器件与信道的测量方法及***、装置、介质 |
CN114323364B (zh) * | 2021-11-11 | 2024-01-30 | 浙江中控技术股份有限公司 | 一种带测振功能的高精度压力传感器及修正方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1079308A (zh) * | 1992-05-23 | 1993-12-08 | 贺守正 | 用电子计算机的电力谐波监测方法 |
CN101460901A (zh) * | 2006-05-15 | 2009-06-17 | 易格斯创新技术***有限公司 | 监测风力发电装置的转子叶片的负荷的方法 |
CN101660957A (zh) * | 2009-09-22 | 2010-03-03 | 南京理工大学 | 滚动直线导轨副摩擦力动态检测装置及其检测方法 |
CN104458116A (zh) * | 2014-12-03 | 2015-03-25 | 沈阳工业大学 | 一种三角波力发生器及其力检测***特性测试方法 |
WO2015187027A1 (en) * | 2014-06-05 | 2015-12-10 | National Oilwell Varco Norway As | Method and device for estimating downhole string variables |
-
2016
- 2016-08-15 CN CN201610669413.8A patent/CN106197807B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1079308A (zh) * | 1992-05-23 | 1993-12-08 | 贺守正 | 用电子计算机的电力谐波监测方法 |
CN101460901A (zh) * | 2006-05-15 | 2009-06-17 | 易格斯创新技术***有限公司 | 监测风力发电装置的转子叶片的负荷的方法 |
CN101660957A (zh) * | 2009-09-22 | 2010-03-03 | 南京理工大学 | 滚动直线导轨副摩擦力动态检测装置及其检测方法 |
WO2015187027A1 (en) * | 2014-06-05 | 2015-12-10 | National Oilwell Varco Norway As | Method and device for estimating downhole string variables |
CN104458116A (zh) * | 2014-12-03 | 2015-03-25 | 沈阳工业大学 | 一种三角波力发生器及其力检测***特性测试方法 |
Non-Patent Citations (1)
Title |
---|
动态压力测量***非线性模型辨识;中国博士学位论文全文数据库 工程科技Ⅱ辑;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20050615(第2期);第3、10-11、55、104-105页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106197807A (zh) | 2016-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106197807B (zh) | 一种用于动态力的测量方法 | |
CN106960068B (zh) | 一种基于脉冲激励响应频谱的模态阻尼比快速计算方法 | |
Peeters et al. | Stochastic system identification for operational modal analysis: a review | |
Kegerise et al. | Mode-switching and nonlinear effects in compressible flow over a cavity | |
Li et al. | Substructural response reconstruction in wavelet domain | |
Kaiser | Advancements in the split Hopkinson bar test | |
Hosder et al. | A non-intrusive polynomial chaos method for uncertainty propagation in CFD simulations | |
Spyrou et al. | Detection of surf-riding behavior of ships in irregular seas | |
He et al. | Approximate entropy analysis of the acoustic emission from defects in rolling element bearings | |
CN105865793B (zh) | 一种提高多转子航空发动机振动监测精度的方法 | |
CN110596247B (zh) | 一种温度变化环境下的超声结构健康监测方法 | |
Martinez-Ríos et al. | Applications of the generalized Morse wavelets: a review | |
CN105138770A (zh) | 基于间接可靠性特征量的航天产品可靠性仿真评估方法 | |
Chaiko et al. | Models for analysis of water hammer in piping with entrapped air | |
Bosma et al. | Wave energy converter modeling in the time domain: A design guide | |
CN105044702B (zh) | 脉冲波形的拟合方法 | |
Thalmann et al. | Trotter derivation of algorithms for Brownian and dissipative particle dynamics | |
Fan et al. | Research on Dynamic Load Identification Based on Explicit Wilson‐θ and Improved Regularization Algorithm | |
CN106446320A (zh) | 基于imf瞬时能量曲率差值的井架钢结构损伤识别方法 | |
JP6977970B2 (ja) | 時系列データ解析装置及び時系列データ解析用プログラム | |
Borminskii et al. | Integrated acoustic monitoring system for jet fuel tanks | |
Stadlmair et al. | Thermoacoustic damping rate determination from combustion noise using Bayesian statistics | |
Waisman et al. | Non-Gaussian model for ringing phenomena in offshore structures | |
JP6869514B2 (ja) | 時系列データ解析装置及び時系列データ解析用プログラム | |
Zhu et al. | A stochastic model for the random impact series method in modal testing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181016 |