CN112834018A - 一种助航灯工作状态的检测方法、存储介质和电子设备 - Google Patents

一种助航灯工作状态的检测方法、存储介质和电子设备 Download PDF

Info

Publication number
CN112834018A
CN112834018A CN202011511303.1A CN202011511303A CN112834018A CN 112834018 A CN112834018 A CN 112834018A CN 202011511303 A CN202011511303 A CN 202011511303A CN 112834018 A CN112834018 A CN 112834018A
Authority
CN
China
Prior art keywords
vibration
conduction coefficient
acceleration
threshold value
navigation aid
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
CN202011511303.1A
Other languages
English (en)
Other versions
CN112834018B (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.)
Harbin Zhengyuan Technology Information Technology Co ltd
Original Assignee
Harbin Zhengyuan Technology Information Technology Co ltd
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 Harbin Zhengyuan Technology Information Technology Co ltd filed Critical Harbin Zhengyuan Technology Information Technology Co ltd
Priority to CN202011511303.1A priority Critical patent/CN112834018B/zh
Publication of CN112834018A publication Critical patent/CN112834018A/zh
Application granted granted Critical
Publication of CN112834018B publication Critical patent/CN112834018B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明实施例公开了一种助航灯工作状态检测方法、存储介质和电子设备。本发明实施例的技术方案通过助航灯的第一方向振动加速度、助航灯的第二方向振动加速度和第三方向振动加速度,结合在飞机起飞降落过程中飞机的质量、飞机第一方向振动加速度、飞机第二方向振动加速度以及飞机第三方法振动加速,获得第一振动传导系数、第二振动传导系数和第三振动传导系数。根据所述第一振动传导系数以及第二振动传导系数的函数与第一预定值的关系,以及将所述第三振动传导系数与多个门限值比较,确定所述助航灯的工作状态。由此,通过本发明实施例的技术方案检测助航灯工作状态的效率高、可靠性高。

Description

一种助航灯工作状态的检测方法、存储介质和电子设备
技术领域
本发明涉及振动检测领域,具体涉及一种助航灯工作状态的检测方法、存储介质和电子设备。
背景技术
随着航空运输量的持续增长,各大机场的跑道也一直保持着接近满负荷的运行状态,助航地灯作为机场的重要组成部分,其质量好坏及稳定性直接影响着机场的运行安全。受飞机频繁的起飞、降落和滑行作用,助航地灯长时间承受着无规律的随机荷载作用。助航灯受到不同频带、不同大小的冲击作用力。其中,飞机降落到跑到时,由于道面不平,飞机产生振动,飞机的振动作用到助航灯上。这是助航灯主要振动来源。也就是说,当飞机振动系统以一定的速度滑过具有不同平整度的水泥混凝土道面时,便会对助航灯产生较为强烈的冲击荷载,由此导致助航灯的松动。
现有技术中,通常通过人工方式检测助航灯的工作状态。人工检查的具体方式为操作人员每天敲击助航灯具,从而判断助航灯螺丝松动情况。这种方式容易造成助航灯具的人为损坏。通过人工检查的方式费时费力且容易出现人为失误。尤其在运输繁忙的季节,由于飞机停航时间较短,需要全体人员加班加点才能完成每天的作业任务。
发明内容
有鉴于此,为了解决人工检查助航灯工作状态效率低、可靠性差的问题,本发明实施例的目的是提供一种助航灯工作状态的检测方法、存储介质和电子设备。
第一方面,本发明实施例提供一种助航灯工作状态的检测方法,包括:
获取待检测助航灯的振动加速度;
获取降落于所述待检测助航灯所在跑道的飞机加速度和飞机质量;
根据所述助航灯振动加速度和所述飞机加速度以及飞机质量计算第一振动传导系数、第二振动传导系数和第三振动传导系数,其中,所述第一振动传导系数对应于水平面第一坐标轴,所述第二振动传导系数对应于水平面第二坐标轴,所述第三振动传导系数对应于垂直坐标轴;
根据历史时间段所述第三振动传导系数确定多个门限值;
根据所述第一振动传导系数以及第二振动传导系数的函数与第一预定值的关系,以及将所述第三振动传导系数与多个门限值比较,确定所述助航灯的工作状态。
优选地,所述门限值包括第一门限值、第二门限值和第三门限值,其中,第一门限值小于第二门限值,第二门限值小于第三门限值;
根据上述助航灯工作状态的检测方法,所述方法还包括:
确定第一门限值;
确定第二门限值;以及
确定第三门限值;
所述根据历史时间段所述第三振动传导系数确定多个门限值包括:
获取历史时间段多个第三振动传导系数组成第三振动传导系数序列;
计算第三振动传导系数序列的期望值和标准差;
响应于所述期望值小于第二预定值,根据所述期望值和所述标准差确定第一概率密度函数;
根据第一区间在所述第一概率密度函数的边界值确定第一门限值。
优选地,所述方法还包括:
获取历史时间段多个第三振动传导系数组成第三振动传导系数序列;
计算第三振动传导系数序列的期望值和标准差;
响应于所述期望值大于等于第二预定值且小于第三预定值,根据所述期望值和所述标准差确定第二概率密度函数;
根据第二区间在所述第二概率密度函数的边界值确定第二门限值。
优选地,所述方法还包括:
获取历史时间段多个第三振动传导系数组成第三振动传导系数序列;
计算第三振动传导系数序列的期望值和标准差;
响应于所述期望值小于等于第三预定值,根据所述期望值和所述标准差确定第三概率密度函数;
根据第二区间在所述第三概率密度函数的边界值确定第三门限值。
优选地,所述第一振动传导系数以及第二振动传导系数的函数为:
第一振动传导系数与第二振动传导系数的平方和的开方。
优选地,所述确定所述助航灯的工作状态包括:
响应于第一振动传导系数与第二振动传导系数的平方和的开方小于第一预定值,且第三振动传导系数大于第一门限值小于等于第二门限值,将所述助航灯的工作状态确定为正常工作状态;
响应于第一振动传导系数与第二振动传导系数的平方和的开方小于第一预定值,且第三振动传导系数大于第一门限值小于等于第二门限值,将所述助航灯的工作状态确定为无碾压松动状态;
响应于第一振动传导系数与第二振动传导系数的平方和的开方大于第一预定值且小于第四预定值,且第一振动传导系数大于第二门限值小于等于第三门限值,将所述助航灯的工作状态确定为碾压状态;
响应于第一振动传导系数与第二振动传导系数的平方和的开方大于第一预定值且小于第四预定值,且第一振动传导系数大于第三门限值,将所述助航灯的工作状态确定为碾压松动状态。
优选地,所述获取待检测助航灯的振动加速度度包括:
获取助航灯第一方向振动加速度;
获取助航灯第二方向振动加速度;以及
获取助航灯第三方向振动加速度;
所述获取降落于所述待检测助航灯所在跑道的飞机加速度和飞机质量包括:
获取飞机质量;
获取飞机第一方向振动加速度;
获取飞机第二方向振动加速度;以及
获取飞机第三方向振动加速度。
优选地,所述根据所述助航灯振动加速度和所述飞机加速度以及飞机质量计算第一振动传导系数、第二振动传导系数和第三振动传导系数包括:
获取助航灯第一方向振动加速度,获取飞机第一方向振动加速度,根据所述飞机第一方向振动加速度与所述助航灯第一方向振动加速度的比值计算所述第一振动传导系数;
获取助航灯第二方向振动加速度,获取飞机第二方向振动加速度,根据所述飞机第二方向振动加速度与所述助航灯第二方向振动加速度的比值计算所述第二振动传导系数;
获取助航灯第三方向振动加速度,获取飞机第三方向振动加速度,获取所述飞机质量,以及获取所述助航灯质量,根据所述助航灯第三方向振动加速度与所述助航灯质量的乘积与所述飞机质量和所述飞机第三方向振动加速度的比值计算所述第三振动传导系数。
第二方面,本发明实施例提供一种存储介质,包括:用于存储计算机程序指令,所述计算机程序指令在被处理器执行时实现上述任一项所述的方法。
第三方面,本发明实施例提供一种电子设备,包括存储器和处理器,其特征在于,所述存储器用于存储一条或多条计算机程序指令,其中,所述一条或多条计算机程序指令被所述处理器执行以上述任一项所述的方法。
本发明实施例的技术方案通过助航灯的第一方向振动加速度、助航灯的第二方向振动加速度和第三方向振动加速度,结合在飞机起飞降落过程中飞机的质量、飞机第一方向振动加速度、飞机第二方向振动加速度以及飞机第三方法振动加速,获得第一振动传导系数、第二振动传导系数和第三振动传导系数。根据历史时间段所述第三振动传导系数确定第一门限值、第二门限值和第三门限值。根据所述第一振动传导系数以及第二振动传导系数的函数与第一预定值的关系,以及将所述第三振动传导系数与多个门限值比较,确定所述助航灯的工作状态。由此,通过本发明实施例的技术方案检测助航灯工作状态的效率高、可靠性高。
附图说明
通过以下参照附图对本发明实施例的描述,本发明的上述以及其它目的、特征和优点将更为清楚,在附图中:
图1是本发明实施例助航灯工作状态检测系统的示意图;
图2是本发明实施例助航灯工作状态检测流程图;
图3是本发明实施例计算振动传导系数的流程图;
图4是本发明实施例计算第一门限值、第二门限值的和第三门限值的流程图;
图5是本发明实施例确定第一门限值的流程图;
图6是本发明实施例根据第一区间在所述第一概率密度函数的边界值的示意图;
图7是本发明实施例确定第二门限值的流程图;
图8是本发明实施例根据第二区间在所述第二概率密度函数的边界值的示意图;
图9是本发明实施例确定第三门限值的流程图;
图10是本发明实施例根据第二区间在所述第三概率密度函数的边界值的示意图;
图11是本发明实施例的电子设备的示意图。
具体实施方式
以下基于实施例对本发明进行描述,但是本发明并不仅仅限于这些实施例。在下文对本发明的细节描述中,详尽描述了一些特定的细节部分。对本领域技术人员来说没有这些细节部分的描述也可以完全理解本发明。为了避免混淆本发明的实质,公知的方法、过程、流程、元件和电路并没有详细叙述。
此外,本领域普通技术人员应当理解,在此提供的附图都是为了说明的目的,并且附图不一定是按比例绘制的。
除非上下文明确要求,否则在说明书的“包括”、“包含”等类似词语应当解释为包含的含义而不是排他或穷举的含义;也就是说,是“包括但不限于”的含义。
在本发明的描述中,需要理解的是,术语“第一”、“第二”等仅用于描述目的,而不能理解为指示或暗示相对重要性。此外,在本发明的描述中,除非另有说明,“多个”的含义是两个或两个以上。
飞机以一定速度通过机场水泥混凝土道面时,由于道面不平将引起飞机轮胎、起落架及车身与装载货物所构成的悬挂质量的振动,并通过具有一定弹性和阻尼作用的起落架传递,将时变激励荷载作用于道面结构,从而引起道面结构的随机振动。此振动通过道面传递到助航地灯上,引起助航地灯跟随振动,通过获取助航地灯的振动信号,并对振动信号处理,可以判定助航地灯的螺丝紧固性。
通过对助行灯振动数据的历史数据进行分析,并结合飞机降落于所述待检测助航灯所在跑道的飞机加速度和飞机质量能够检测助行灯的工作状态。
图1是本发明实施例助航灯工作状态检测系统的示意图。
助航灯是引导飞机在夜间和复杂天气条件下起飞、着陆和滑行的目视助航设备。助航灯工作状态异常有可能导致飞机事故。为了检测助航灯的工作状态,本发明实施例提供一种助行灯工作状态检测系统。
参照图1,本发明实施例的助航灯工作状态检测系统包括:助航灯1、传感器2、振动检测转换电路3、服务器4和飞机监视系统5。
传感器2用于采集助航灯的振动信号,振动检测转换电路3用于将振动信号转换为服务器4能够识别的振动加速度。
传感器2包括:第一传感器20、第二传感器21和第三传感器22。
具体地,传感器20采集第一方向的振动信号,振动检测转换电路3将其转换为助航灯第一振动加速度adx
传感器21采集第二方向的振动信号,振动检测转换电路3将其转换为助航灯第二方向的振动加速度ady
传感器22采集第三方向的振动信号,振动检测转换电路3将其转换为助航灯第三方向的振动加速度adz
振动检测转换电路3与服务器4通信连接。振动检测转换电路3将助航灯第一振动加速度adx、助航灯第二振动加速度ady和助航灯第三振动加速度adz发送到服务器4。在一种可选的实现方式中,振动检测转换电路3通过有线电缆的方式与服务器4通信连接。从而将助航灯第一振动加速度adx、助航灯第二振动加速度ady和助航灯第三振动加速度adz发送到服务器4。
在另一种可选的实现方式中,振动检测转换电路3通过无线的方式与服务器4通信连接。从而将助航灯第一振动加速度adx、助航灯第二振动加速度ady和助航灯第三振动加速度adz发送到服务器4。
服务器4可以是集群也可以是云服务器。
服务器4与飞机监视系统5(飞机监视系统,ADS-B,Automatic DependentSurveillance-Broadcast)通信连接。飞机监视系统5能够实时获得飞机起飞、滑行和降落过程中的各项数据。飞机监视系统5还能够将实时获得飞机起飞、滑行和降落过程中的各项数据存储在数据库中保存形成飞机起飞、滑行和降落过程中的历史数据。
由于飞机起飞、滑行和降落过程对助航灯影响最大,通过服务器4与飞机监视系统5通信连接,服务器4能够获取飞机起飞、滑行和降落过程的历史数据。
具体地,通服务器4与飞机监视系统5通信连接,服务器4获取飞机的质量和飞机的加速度。
服务器4通过振动检测转换电路3获取的助航灯第一振动加速度adx、助航灯第二振动加速度ady和助航灯第三振动加速度adz,以及通过飞机监视系统5获取飞机的质量和飞机的加速度,经本发明实施例助航灯工作状态检测后,服务器4能够检测助行灯的工作状态。
图2是本发明实施例助航灯工作状态检测流程图。
本发明实施例的助航灯工作状态检测在服务器4上运行,包括步骤100-步骤500。
步骤100,获取待检测助航灯的振动加速度。
助航灯的振动加速度包括助航灯第一方向振动加速度adx、助航灯第二方向振动加速度ady以及助航灯第三方向振动加速度adz
结合图1,在一种可选的实现方式中,传感器2用于获取助航灯的振动信号。传感器与振动检测转换电路3通信连接。振动检测转换电路3将振动信号转换为振动加速度,并发送给服务器4。从而,服务器4能够获取助航灯的振动加速度。在另一种可选的实现方式中,传感器2为多轴传感器,通过多轴传感器能够获得助航灯的振动信号。也即传感器2能够直接获取X方向的振动信号、Y方向的振动信号和Z方向的振动信号。传感器2与振动检测转换电路3电连接。振动检测转换电路3将X方向的振动信号处理后转换为助航灯第一方向振动加速度adx。振动检测转换电路3将Y方向的振动信号处理后转换为助航灯第二方向振动加速度ady。振动检测转换电路3将Z方向的振动信号处理后转换为助航灯第三方向振动加速度adz
在步骤100中,获取助航灯第一方向振动加速度adx
具体地,第一传感器20获取助航灯X方向的振动信号。
振动检测转换电路3与第一传感器20电连接。第一传感器20将助航灯X方向的振动信号发送到振动检测转换电路3。振动检测转换电路3将助航灯X方向的振动信号处理后转换为助航灯第一方向振动加速度adx。服务器4与振动检测转换电路3通信连接。从而服务器4能够获取助航灯第一方向振动加速度adx
在步骤100中获取助航灯第二方向振动加速度ady
具体地,第二传感器21获取助航灯Y方向的振动信号。
振动检测转换电路3与第二传感器20电连接。第二传感器21将助航灯Y方向的振动信号发送到振动检测转换电路3。振动检测转换电路3将助航灯Y方向的振动信号处理后转换为助航灯第二方向振动加速度ady。服务器4与振动检测转换电路3通信连接。从而服务器4能够获取助航灯第二方向振动加速度ady
在步骤100中,获取助航灯第三方向振动加速度adz
第三传感器22获取助航灯Z方向的振动信号。
振动检测转换电路3与第三传感器20电连接。第三传感器21将助航灯Z方向的振动信号发送到振动检测转换电路3。振动检测转换电路3将助航灯Z方向的振动信号处理后转换为助航灯第三方向振动加速度adz。服务器4与振动检测转换电路3通信连接。从而服务器4能够获取助航灯第三方向振动加速度adz
步骤200,获取降落于所述待检测助航灯所在跑道的飞机加速度和飞机质量。其中,飞机加速度包括:飞机第一方向振动加速度apx、飞机第二方向振动加速度apy、获取飞机第三方向振动加速度apz
在步骤200中获取飞机质量mp、获取飞机第一方向振动加速度apx、获取飞机第二方向振动加速度apy、以及获取飞机第三方向振动加速度apz
结合图1,在一种可选的实现方式中,服务器4与服务器监视系统5通信连接。从而,服务器4能够获取监视系统5中的飞机质量mp、飞机第一方向振动加速度apx、飞机第二方向振动加速度apy、以及飞机第三方向振动加速度apz
为了检测助航灯工作状态,需要确定第一振动传导系数kx、第二振动传导系数ky和第三振动传导系数kz
步骤300,根据所述助航灯振动加速度和所述飞机加速度以及飞机质量计算第一振动传导系数、第二振动传导系数和第三振动传导系数,其中,所述第一振动传导系数对应于水平面第一坐标轴,所述第二振动传导系数对应于水平面第二坐标轴,所述第三振动传导系数对应于垂直坐标轴。
图3是本发明实施例计算振动传导系数的流程图。
参照图3,步骤300包括3个子步骤,分别是步骤310-步骤330。
步骤310,计算第一振动传导系数。
在步骤310中,获取助航灯第一方向振动加速度adx,获取飞机第一方向振动加速度apx,根据所述飞机第一方向振动加速度apx与所述助航灯第一方向振动加速度adx的比值计算所述第一振动传导系数kx
其中,第一振动传导系数kx=adx/apx
步骤320计算第二振动传导系数。
在步骤320中,获取助航灯第二方向振动加速度ady,获取飞机第二方向振动加速度apy,根据所述飞机第二方向振动加速度apy与所述助航灯第二方向振动加速度ady的比值计算所述第二振动传导系数ky
其中,第二振动传导系数ky=ady/apy
步骤330计算第三振动传导系数。
在步骤330中,获取助航灯第三方向振动加速度adz,获取飞机第三方向振动加速度apz,获取所述飞机质量mp,以及获取所述助航灯质量md,根据所述助航灯第三方向振动加速度adz与所述助航灯质量md的乘积与所述飞机质量mp和所述飞机第三方向振动加速度apz的比值计算所述第三振动传导系数kz
第三振动传导系数kz=md*adz/(mp*apz)
为了检测助航灯工作状态,还需要预先确定多个门限值。
步骤400,根据历史时间段所述第三振动传导系数确定多个门限值。
所述门限值包括第一门限值、第二门限值和第三门限值,其中,第一门限值小于第二门限值,第二门限值小于第三门限值。
图4是本发明实施例计算第一门限值、第二门限值的和第三门限值的流程图。
参照图4,步骤400包括3个子步骤,分别是步骤410-步骤430。
步骤410,确定第一门限值。
图5是确定第一门限值的流程图。
参照图5,步骤410包括4个子步骤,分别是步骤411-步骤414。
步骤411,获取历史时间段多个第三振动传导系数组成第三振动传导系数序列。
在一种可选的实现方式中,服务器4预先通过飞机监视系统5获取历史时段的多个第三振动传导系数kz组成第三振动传导系数序列P1(kz)。
历史时段可以是以天为单位的历史时段,或者以星期为单位的历史时段,或者是以月为单位的历史时段。
例如,服务器4通过飞机监视系统5获取前一天的多个第三振动传导系数kz1,kz2,kz3,……kzi,组成第三振动传导系数序列P1(kz):
P1(kz)={kz1,kz2,kz3,……kzi},在一种可选的实现方式中,将历史时段中i个第三振动传导系数组成第三振动传导系数系列P1(kz)。
步骤412,计算第三振动传导系数序列的期望值和标准差。
第三振动传导系数序列期望值的计算方法为:
μ1=E(kzi)
由于kz是离散值,kz 1,kz 2,kz 3,……kz i对应的概率为
Figure BDA0002846490780000111
p1,p2,p3,……pi,第三振动传导系数序列期望值E(kzi)=∑i kzi*pi。也即第三振动传导系数序列期望值是各第三振动传导系数与对应概率的乘积的和。第三振动传导系数序列标准差的计算方法为:
其中sqrt()是算数平方根。
步骤413,响应于所述期望值小于第二预定值,根据所述期望值和所述标准差确定第一概率密度函数。
优选地,根据经验值,第二预定值为0.5。
优选地,第一概率密度函数为高斯分布的概率密度函数。
当μ1<0.5,将期望值μ1和标准差σ1带入高斯分布的确定的概率密度函数,确定第一概率密度函数。第一概率密度函数为:
Figure BDA0002846490780000112
其中,exp()是指数函数。
步骤414,根据第一区间在所述第一概率密度函数的边界值确定第一门限值。
图6是本发明实施例根据第一区间在所述第一概率密度函数的边界值的示意图。
参照图6优选地,根据经验值,第一区间为[μ1-3σ1,μ1+3σ1]。
由此,计算出第一区间在所述第一概率密度函数的边界值A1。将边界值A1确定为第一门限值。
接下来,确定第二门限值。
步骤420,确定第二门限值。
图7是本发明实施例确定第二门限值的流程图。
参照图7,步骤420包括4个子步骤,分别是步骤421-步骤424。
步骤421,获取历史时间段多个第三振动传导系数组成第三振动传导系数序列P2(kz):
P2(kz)={(kz1,kz2,kz3,……kzm}。在一种可选的实现方式中,将历史时段中m个第三振动传导系数组成第三振动传导系数系列。
步骤422,计算第三振动传导系数序列的期望值μ2和标准差σ2
第三振动传导系数序列期望值为:
μ2=E(kzm)
由于kz是离散值,kz 1,kz 2,kz 3,……kz m对应的概率为p1,p2,p3,……pm,第三振动传导系数序列期望值E(kzi)=∑i kzm*pm。也即第三振动传导系数序列期望值是各第三振动传导系数与对应概率的乘积的和。
第三振动传导系数序列标准差为:
Figure BDA0002846490780000121
其中,sqrt()是算数平方根。
步骤423,响应于所述期望值大于等于第二预定值且小于第三预定值,根据所述期望值和所述标准差确定第二概率密度函数。
优选地,根据经验值,第二预定值为0.5,第三预定值为1。
优选地,第二概率密度函数为高斯分布的概率密度函数。
当0.5≤μ3<1,确定第三概率密度函数。第三概率密度函数为:
将期望值μ2和标准差σ2带入高斯分布的确定的概率密度函数,确定第二概率密度函数。第二概率密度函数为:
Figure BDA0002846490780000122
其中,exp()是指数函数。
步骤424,根据第二区间在所述第二概率密度函数的边界值确定第二门限值。
图8是本发明实施例根据第二区间在所述第二概率密度函数的边界值的示意图。
参照图8,优选地,根据经验值,第二区间为[μ21,μ1+3σ1]。
根据第二区间在所述第二概率密度函数的边界值A2确定为第二门限值。
由此,计算出第二区间在所述第二概率密度函数的边界值A2。将边界值A2确定为第二门限值。
接下来,确定第三门限值。
步骤430,确定第三门限值。
图9是确定第三门限值的流程图。
参照图9,步骤430包括4个子步骤,分别是步骤431-步骤434。
步骤431,获取历史时间段多个第三振动传导系数组成第三振动传导系数序列P3(kz):
P3(kz)={kz1,kz2,kz3,……kzj}。例如,将历史时段中j个第三振动传导系数组成第三振动传导系数系列。
步骤432,计算第三振动传导系数序列的期望值和标准差。
第三振动传导系数序列期望值μ3为:
μ3=E(kzj)
由于kz是离散值,kz 1,kz 2,kz 3,……kzj对应的概率为p1,p2,p3,……pj,第三振动传导系数序列期望值E(kzi)=∑ikzj*pj。也即第三振动传导系数序列期望值是各第三振动传导系数与对应概率的乘积的和。
第三振动传导系数序列标准差σ3为:
Figure BDA0002846490780000131
其中,sqrt()是算数平方根。
步骤433,响应于所述期望值μ3大于等于第三预定值,根据所述期望值μ3和所述标准差σ3确定第三概率密度函数。
优选地,根据经验值,第三预定值为1。
优选地,第三概率密度函数为高斯分布的概率密度函数。
当μ3≥1,确定第三概率密度函数。确定第三概率密度函数为:
Figure BDA0002846490780000141
其中,exp()是指数函数。
步骤434,根据第二区间在所述第三概率密度函数的边界值确定第三门限值。
图10是本发明实施例根据第二区间在所述第三概率密度函数的边界值的示意图。
参照图10,优选地,根据经验值,第二区间为
Figure BDA0002846490780000142
根据第二区间在所述第三概率密度函数的边界值A3确定第三门限值。
由此,计算出第二区间在所述第三概率密度函数的边界值A3。将边界值A3确定为第三门限值。
接下来,确定助航灯的工作状态。
步骤500,根据所述第一振动传导系数以及第二振动传导系数的函数与第一预定值的关系,以及将所述第三振动传导系数与多个门限值比较,确定所述助航灯的工作状态。
步骤500包括6个子步骤,分别是步骤510-步骤560。
步骤510,计算第一振动传导系数以及第二振动传导系数的函数。
第一振动传导系数以及第二振动传导系数的函数为:
第一振动传导系数与第二振动传导系数的平方和的开方
Figure BDA0002846490780000143
优选地,根据经验值,将第一预定值确定为0.707。
步骤520,判断第一振动传导系数与第二振动传导系数的平方和与第一预定值或第二预定值的关系,以及判断第三振动传导系数与第一门限值、第二门限值、第三门限值的关系。
步骤530,响应于第一振动传导系数与第二振动传导系数的平方和的开方小于第一预定值,且第三振动传导系数大于第一门限值小于等于第二门限值,将所述助航灯的工作状态确定为正常工作状态。
也就是说,当
Figure BDA0002846490780000151
kz≤第一门限值时,将所述助航灯的工作状态确定为正常工作状态。
步骤540,响应于第一振动传导系数与第二振动传导系数的平方和的开方大于第一预定值,且第三振动传导系数大于第二门限值小于等于第三门限值,将所述助航灯的工作状态确定为碾压状态。
也就是说,当
Figure BDA0002846490780000152
第一门限值<kz≤第二门限值时,将所述助航灯的工作状态确定为无松动碾压状态。
步骤550,响应于第一振动传导系数与第二振动传导系数的平方和的开方大于第一预定值且小于第四预定值,且第三振动传导系数大于第三门限值,将所述助航灯的工作状态确定为碾压状态。
优选地,根据经验值,将第四预定值确定为1.414。
也就是说,当
Figure BDA0002846490780000153
时,第二门限值<kz≤第三门限值时,将所述助航灯的工作状态确定为碾压松动状态。
步骤560,响应于第一振动传导系数与第二振动传导系数的平方和的开方大于第一预定值且小于第四预定值时,且第三振动传导系数大于第三门限值,将所述助航灯的工作状态确定为碾压松动状态。
也就是说,当
Figure BDA0002846490780000154
时,kz>第三门限值时,将所述助航灯的工作状态确定为碾压松动状态。
本发明实施例的技术方案通过助航灯的第一方向振动加速度、助航灯的第二方向振动加速度和第三方向振动加速度,结合在飞机起飞降落过程中飞机的质量、飞机第一方向振动加速度、飞机第二方向振动加速度以及飞机第三方法振动加速,获得第一振动传导系数、第二振动传导系数和第三振动传导系数。根据历史时间段所述第三振动传导系数确定第一门限值、第二门限值和第三门限值。根据所述第一振动传导系数以及第二振动传导系数的函数与第一预定值的关系,以及将所述第三振动传导系数与多个门限值比较,确定所述助航灯的工作状态。由此,通过本发明实施例的技术方案检测助航灯工作状态效率高、可靠性高。
图11是本发明实施例的电子设备的示意图。
如图11所示的电子设备11包括通用的硬件结构,其至少包括处理器111和存储器112。处理器111和存储器112通过总线113连接。存储器112适于存储处理器111可执行的指令或程序。处理器111可以是独立的微处理器,也可以是一个或者多个微处理器集合。由此,处理器111通过执行存储器112所存储的指令,从而执行如上所述的本发明实施例的方法流程实现对于数据的处理和对于其它装置的控制。总线113将上述多个组件连接在一起,同时将上述组件连接到显示控制器114和显示装置以及输入/输出(I/O)装置115。输入/输出(I/O)装置115可以是鼠标、键盘、调制解调器、网络接口、触控输入装置、体感输入装置、打印机以及本领域公知的其他装置。典型地,输入/输出装置115通过输入/输出(I/O)控制器116与系统相连。
本领域的技术人员应明白,本申请的实施例可提供方法、装置(设备)或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可读存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品。
本申请是参照根据本申请实施例的方法、装置(设备)和计算机程序产品的流程图来描述的。应理解可由计算机程序指令实现流程图中的每一流程。
这些计算机程序指令可以存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现流程图一个流程或多个流程中指定的功能。
也可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程中指定的功能的装置。
本发明的另一实施例涉及一种非易失性可读存储介质,用于存储计算机可读程序,所述计算机可读程序用于供计算机执行上述部分或全部方法的实施例。
即,本领域技术人员可以理解,实现上述实施例的方法中的全部或部分步骤是可以通过程序来指定相关的硬件来完成,该程序存储在一个可读存储介质中,包括若干指令用以使得一个设备(可以是单片机,芯片等)或处理器(processor)执行本申请各实施例所述方法的全部或部分步骤。而前述的可读存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述仅为本发明的优选实施例,并不用于限制本发明,对于本领域技术人员而言,本发明可以有各种改动和变化。凡在本发明的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种助航灯工作状态的检测方法,其特征在于,所述方法包括:
获取待检测助航灯的振动加速度;
获取降落于所述待检测助航灯所在跑道的飞机加速度和飞机质量;
根据所述助航灯振动加速度和所述飞机加速度以及飞机质量计算第一振动传导系数、第二振动传导系数和第三振动传导系数,其中,所述第一振动传导系数对应于水平面第一坐标轴,所述第二振动传导系数对应于水平面第二坐标轴,所述第三振动传导系数对应于垂直坐标轴;
根据历史时间段所述第三振动传导系数确定多个门限值;
根据所述第一振动传导系数以及第二振动传导系数的函数与第一预定值的关系,以及将所述第三振动传导系数与多个门限值比较,确定所述助航灯的工作状态。
2.根据权利要求1所述的助航灯工作状态的检测方法,其特征在于,所述门限值包括第一门限值、第二门限值和第三门限值,其中,第一门限值小于第二门限值,第二门限值小于第三门限值;
所述方法还包括:
确定第一门限值;
确定第二门限值;以及
确定第三门限值;
所述根据历史时间段所述第三振动传导系数确定多个门限值包括:
获取历史时间段多个第三振动传导系数组成第三振动传导系数序列;
计算第三振动传导系数序列的期望值和标准差;
响应于所述期望值小于第二预定值,根据所述期望值和所述标准差确定第一概率密度函数;
根据第一区间在所述第一概率密度函数的边界值确定第一门限值。
3.根据权利要求2所述的助航灯工作状态的检测方法,其特征在于,所述方法还包括:
获取历史时间段多个第三振动传导系数组成第三振动传导系数序列;
计算第三振动传导系数序列的期望值和标准差;
响应于所述期望值大于等于第二预定值且小于第三预定值,根据所述期望值和所述标准差确定第二概率密度函数;
根据第二区间在所述第二概率密度函数的边界值确定第二门限值。
4.根据权利要求3所述的助航灯工作状态的检测方法,其特征在于,所述方法还包括:
获取历史时间段多个第三振动传导系数组成第三振动传导系数序列;
计算第三振动传导系数序列的期望值和标准差;
响应于所述期望值小于等于第三预定值,根据所述期望值和所述标准差确定第三概率密度函数;
根据第二区间在所述第三概率密度函数的边界值确定第三门限值。
5.根据权利要求4所述的助航灯工作状态的检测方法,其特征在于,所述第一振动传导系数以及第二振动传导系数的函数为:
第一振动传导系数与第二振动传导系数的平方和的开方。
6.根据权利要求5所述的助航灯工作状态的检测方法,其特征在于,所述确定所述助航灯的工作状态包括:
响应于第一振动传导系数与第二振动传导系数的平方和的开方小于第一预定值,且第三振动传导系数大于第一门限值小于等于第二门限值,将所述助航灯的工作状态确定为正常工作状态;
响应于第一振动传导系数与第二振动传导系数的平方和的开方小于第一预定值,且第三振动传导系数大于第一门限值小于等于第二门限值,将所述助航灯的工作状态确定为无碾压松动状态;
响应于第一振动传导系数与第二振动传导系数的平方和的开方大于第一预定值且小于第四预定值,且第一振动传导系数大于第二门限值小于等于第三门限值,将所述助航灯的工作状态确定为碾压状态;
响应于第一振动传导系数与第二振动传导系数的平方和的开方大于第一预定值且小于第四预定值,且第一振动传导系数大于第三门限值,将所述助航灯的工作状态确定为碾压松动状态。
7.根据权利要求1所述的助航灯工作状态的检测方法,其特征在于,所述获取待检测助航灯的振动加速度度包括:
获取助航灯第一方向振动加速度;
获取助航灯第二方向振动加速度;以及
获取助航灯第三方向振动加速度;
所述获取降落于所述待检测助航灯所在跑道的飞机加速度和飞机质量包括:
获取飞机质量;
获取飞机第一方向振动加速度;
获取飞机第二方向振动加速度;以及
获取飞机第三方向振动加速度。
8.根据权利要求7所述的助航灯工作状态的检测方法,其特征在于,所述根据所述助航灯振动加速度和所述飞机加速度以及飞机质量计算第一振动传导系数、第二振动传导系数和第三振动传导系数包括:
获取助航灯第一方向振动加速度,获取飞机第一方向振动加速度,根据所述飞机第一方向振动加速度与所述助航灯第一方向振动加速度的比值计算所述第一振动传导系数;
获取助航灯第二方向振动加速度,获取飞机第二方向振动加速度,根据所述飞机第二方向振动加速度与所述助航灯第二方向振动加速度的比值计算所述第二振动传导系数;
获取助航灯第三方向振动加速度,获取飞机第三方向振动加速度,获取所述飞机质量,以及获取所述助航灯质量,根据所述助航灯第三方向振动加速度与所述助航灯质量的乘积与所述飞机质量和所述飞机第三方向振动加速度的比值计算所述第三振动传导系数。
9.一种计算机可读存储介质,用于存储计算机程序指令,其特征在于,所述计算机程序指令在被处理器执行时实现如权利要求1-8中任一项所述的方法。
10.一种电子设备,包括存储器和处理器,其特征在于,所述存储器用于存储一条或多条计算机程序指令,其中,所述一条或多条计算机程序指令被所述处理器执行以实现如权利要求1-8中任一项所述的方法。
CN202011511303.1A 2020-12-18 2020-12-18 一种助航灯工作状态的检测方法、存储介质和电子设备 Active CN112834018B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011511303.1A CN112834018B (zh) 2020-12-18 2020-12-18 一种助航灯工作状态的检测方法、存储介质和电子设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011511303.1A CN112834018B (zh) 2020-12-18 2020-12-18 一种助航灯工作状态的检测方法、存储介质和电子设备

Publications (2)

Publication Number Publication Date
CN112834018A true CN112834018A (zh) 2021-05-25
CN112834018B CN112834018B (zh) 2022-11-04

Family

ID=75923633

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011511303.1A Active CN112834018B (zh) 2020-12-18 2020-12-18 一种助航灯工作状态的检测方法、存储介质和电子设备

Country Status (1)

Country Link
CN (1) CN112834018B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024045412A1 (zh) * 2022-09-01 2024-03-07 中国港湾工程有限责任公司 用于飞机场目视助航设备的管理系统及方法

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5530764A (en) * 1993-03-19 1996-06-25 Mazda Motor Corporation Vibration control system for an automotive vehicle
RU2098764C1 (ru) * 1996-05-29 1997-12-10 Русланов Александр Семенович Способ определения местоположения подвижных объектов и устройство для его реализации
CN102680078A (zh) * 2012-04-18 2012-09-19 上海现代建筑设计(集团)有限公司 地铁交通引起的邻近建筑振动三维频域分析方法
CN106786673A (zh) * 2017-02-07 2017-05-31 华北电力科学研究院有限责任公司 双馈风机串补输电系统次同步谐振的抑制方法及装置
RU2672529C1 (ru) * 2017-12-18 2018-11-15 Российская Федерация, от имени которой выступает Государственная корпорация по атомной энергии "Росатом" (Госкорпорация "Росатом") Способ формирования спектров случайной вибрации
US20190126959A1 (en) * 2017-10-31 2019-05-02 Mitsubishi Heavy Industries, Ltd. Abnormality detection device, abnormality detection method, and program
CN110132406A (zh) * 2019-06-04 2019-08-16 中车工业研究院有限公司 一种永磁直驱电力机车联轴器振动检测方法及系统
CN110411757A (zh) * 2019-07-30 2019-11-05 安徽江淮汽车集团股份有限公司 轴头动态载荷计算方法、装置、设备及存储介质
CN110488042A (zh) * 2019-07-22 2019-11-22 中车青岛四方机车车辆股份有限公司 一种列车加速度检测方法、系统、电子设备及存储介质
RU2715345C1 (ru) * 2019-08-07 2020-02-26 Общество с ограниченной ответственностью научно-производственное предприятие "ТИК" Пьезоэлектрический измерительный преобразователь
CN111538311A (zh) * 2020-04-22 2020-08-14 北京航空航天大学 一种基于数据挖掘的机械设备柔性多状态自适应预警方法及装置
CN111551352A (zh) * 2020-04-01 2020-08-18 国网宁夏电力有限公司电力科学研究院 一种gis设备的断路器的状态的检测方法及系统
CN111829742A (zh) * 2019-03-29 2020-10-27 Oppo(重庆)智能科技有限公司 电子设备振动测试方法、装置、电子设备及存储介质
CN111954143A (zh) * 2020-08-29 2020-11-17 深圳市韶音科技有限公司 一种获取振动传递函数的方法和系统

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5530764A (en) * 1993-03-19 1996-06-25 Mazda Motor Corporation Vibration control system for an automotive vehicle
RU2098764C1 (ru) * 1996-05-29 1997-12-10 Русланов Александр Семенович Способ определения местоположения подвижных объектов и устройство для его реализации
CN102680078A (zh) * 2012-04-18 2012-09-19 上海现代建筑设计(集团)有限公司 地铁交通引起的邻近建筑振动三维频域分析方法
CN106786673A (zh) * 2017-02-07 2017-05-31 华北电力科学研究院有限责任公司 双馈风机串补输电系统次同步谐振的抑制方法及装置
US20190126959A1 (en) * 2017-10-31 2019-05-02 Mitsubishi Heavy Industries, Ltd. Abnormality detection device, abnormality detection method, and program
RU2672529C1 (ru) * 2017-12-18 2018-11-15 Российская Федерация, от имени которой выступает Государственная корпорация по атомной энергии "Росатом" (Госкорпорация "Росатом") Способ формирования спектров случайной вибрации
CN111829742A (zh) * 2019-03-29 2020-10-27 Oppo(重庆)智能科技有限公司 电子设备振动测试方法、装置、电子设备及存储介质
CN110132406A (zh) * 2019-06-04 2019-08-16 中车工业研究院有限公司 一种永磁直驱电力机车联轴器振动检测方法及系统
CN110488042A (zh) * 2019-07-22 2019-11-22 中车青岛四方机车车辆股份有限公司 一种列车加速度检测方法、系统、电子设备及存储介质
CN110411757A (zh) * 2019-07-30 2019-11-05 安徽江淮汽车集团股份有限公司 轴头动态载荷计算方法、装置、设备及存储介质
RU2715345C1 (ru) * 2019-08-07 2020-02-26 Общество с ограниченной ответственностью научно-производственное предприятие "ТИК" Пьезоэлектрический измерительный преобразователь
CN111551352A (zh) * 2020-04-01 2020-08-18 国网宁夏电力有限公司电力科学研究院 一种gis设备的断路器的状态的检测方法及系统
CN111538311A (zh) * 2020-04-22 2020-08-14 北京航空航天大学 一种基于数据挖掘的机械设备柔性多状态自适应预警方法及装置
CN111954143A (zh) * 2020-08-29 2020-11-17 深圳市韶音科技有限公司 一种获取振动传递函数的方法和系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李玉学 等: "大跨屋盖结构风致背景响应和共振响应实用组合方法", 《振动与冲击》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024045412A1 (zh) * 2022-09-01 2024-03-07 中国港湾工程有限责任公司 用于飞机场目视助航设备的管理系统及方法

Also Published As

Publication number Publication date
CN112834018B (zh) 2022-11-04

Similar Documents

Publication Publication Date Title
CN102812502B (zh) 监控跑道状况的装置及方法
US9058743B2 (en) Automatic landing method and device for an aircraft on a strong slope runway
US8355850B2 (en) Method and device of calculating aircraft braking friction and other relating landing performance parameters based on the data received from aircraft&#39;s on board flight data management system
EP2092401B1 (en) Hard landing detection
JP2017515716A (ja) ブレーキ制御システムを用いて滑走路状態を報告する方法
US11231521B2 (en) Method and server for supporting aviation safety
CN112834018B (zh) 一种助航灯工作状态的检测方法、存储介质和电子设备
CN103377565A (zh) 利用起飞和着陆性能数据来改进跑道感知的系统和方法
DK201300589A1 (en) Dynamic alarm zones for bird detection systems
JP2013249052A (ja) 翼端センサ情報を選別するためのシステムおよび方法
US10930163B2 (en) Systems and methods for validating real-time condition of a landing field using aircraft data
EP1811887A2 (en) Takeoff and landing performance indicator for fixed wing aircraft
LU102309B1 (en) A Method for Managing the Clearance around an Airport, an Intelligent Terminal and a Storage Medium
CN112579464A (zh) 自动驾驶算法的校验方法、装置、设备以及存储介质
CN110645899B (zh) 一种结构裂缝分布式光纤监测方法及装置
CN112329101A (zh) 机场跑道的净空评定方法
CN113538660A (zh) 一种基于bim和点云测量的输变电工程验收方法
CN111524324B (zh) Fod报警级别判定方法、电子设备及存储介质
CN109754193A (zh) 基于航空器性能的ads-b航迹去噪方法
CN115079293A (zh) 一种用于机场跑道异物检测和道面检测的设备及方法
CN111477036B (zh) 一种空管自动化系统航空器高度异常检测方法
CN107522052B (zh) 一种电梯异常预警系统和方法
CN115293225A (zh) 飞行员平飘顶杆成因分析方法和装置
Balicki et al. Safety performance indicators assessment for small aircraft airframe systems
CN114399498A (zh) 一种机场跑道表面状况评估方法及系统

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