CN109367765B - 基于路况识别的飞机防滑刹车自适应控制方法 - Google Patents

基于路况识别的飞机防滑刹车自适应控制方法 Download PDF

Info

Publication number
CN109367765B
CN109367765B CN201811098223.0A CN201811098223A CN109367765B CN 109367765 B CN109367765 B CN 109367765B CN 201811098223 A CN201811098223 A CN 201811098223A CN 109367765 B CN109367765 B CN 109367765B
Authority
CN
China
Prior art keywords
wheel
rate
brake
wheel speed
deceleration rate
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
CN201811098223.0A
Other languages
English (en)
Other versions
CN109367765A (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 hangchen Airborne Intelligent System Technology Co.,Ltd.
Original Assignee
Beijing Hangchen Airborne Intelligent System 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 Beijing Hangchen Airborne Intelligent System Technology Co ltd filed Critical Beijing Hangchen Airborne Intelligent System Technology Co ltd
Priority to CN201811098223.0A priority Critical patent/CN109367765B/zh
Publication of CN109367765A publication Critical patent/CN109367765A/zh
Application granted granted Critical
Publication of CN109367765B publication Critical patent/CN109367765B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B64—AIRCRAFT; AVIATION; COSMONAUTICS
    • B64C—AEROPLANES; HELICOPTERS
    • B64C25/00—Alighting gear
    • B64C25/32—Alighting gear characterised by elements which contact the ground or similar surface 
    • B64C25/42—Arrangement or adaptation of brakes
    • B64C25/44—Actuating mechanisms
    • B64C25/46—Brake regulators for preventing skidding or aircraft somersaulting

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Regulating Braking Force (AREA)

Abstract

本公开提供一种基于路况识别的飞机防滑刹车控制方法，包括：(1)根据实时采集的刹车压力和实时采集的机轮轮速获得初始的刹车系统关键参数；(2)根据初始的刹车系统关键参数和刹车压力，获得刹车结合力的变化率和滑移率的变化率，根据初始的刹车系统关键参数获得估计的机轮轮速；(3)获得实际机轮轮速和估计的机轮轮速的差值；根据实际轮速和估计轮速的差值以及刹车压力，获得刹车系统关键参数的调节量；(4)根据刹车结合力的变化率和滑移率的变化率生成目标机轮减速率；(5)根据目标机轮减速率生成目标机轮轮速；(6)根据目标机轮轮速和实际机轮轮速，生成防滑控制量并输出给飞机刹车系统。

Description

基于路况识别的飞机防滑刹车自适应控制方法
技术领域
本公开属于飞机液压刹车系统控制领域，具体涉及一种基于路况识别的飞机防滑刹车自适应控制方法。
背景技术
飞机刹车系统是保证飞机安全起降的关键子系统之一，飞机能否安全可靠起降，直接影响到飞机及机载人员的安全。随着飞机体积重量的不断增大，着陆环境日趋复杂，使得刹车系统受到诸多外部因素和内部不稳定因素的影响，令飞机起降的可靠性和安全性问题更加突出。
目前大多数飞机采用自动防滑液压刹车系统，它以调节作用在刹车机轮上的压力来改变刹车力矩以保证飞机机轮无抱死刹车(拖胎)，而刹车压力的调节应当与刹车时飞机滑行速度和路面情况相适应，从理论上要求，刹车控制装置应能自动地对轮胎和跑道摩擦系数的最大值进行寻优，使飞机的刹车滑跑距离最短，轮胎和刹车片磨损量最小，实现安全，高效，可靠，迅速地刹停飞机，然而在实际刹车控制器的设计中，实现起来面临诸多难以克服的困难。
在影响刹车效率的众多因素中，除了飞机自身在着陆前飞行速度受着陆环境的影响，以及着路场上风力，路况等自然因素影响，另一关键因素就是防滑刹车系统控制器的控制规律。飞机防滑刹车控制律的问题一直是国内外广泛关注的课题，国内如北航等知名高校，美国兰利中心，波音公司等均投入大量人力物力予以研究。
我国的防滑刹车系统从机械防滑发展到模拟电子防滑，目前已进入数字电子防滑的阶段，但国内常选用多门限PBM控制设计防滑刹车系统的控制律，这种方法通常存在低速打滑现象，所以应用新的控制理论进行防滑刹车系统的设计，也将成为今后防滑刹车的重要研究方向。
发明内容
为了解决至少一个上述技术问题，一方面，本公开提供一种基于路况识别的飞机防滑刹车控制方法，包括：
(1)根据实时采集的刹车压力和实时采集的机轮轮速，获得初始的刹车系统关键参数；
(2)根据初始的刹车系统关键参数和所述刹车压力，获得刹车结合力的变化率和滑移率的变化率，根据刹车系统关键参数获得估计的机轮轮速；
(3)获得所述机轮轮速和所述估计的机轮轮速的差值，即实际轮速和估计轮速的差值；根据所述实际轮速和估计轮速的差值以及刹车压力，获得刹车系统关键参数的调节量；刹车系统关键参数等于所述初始的刹车系统关键参数与所述调节量之和；
(4)根据所述刹车结合力的变化率和所述滑移率的变化率生成目标机轮减速率；
(5)根据所述目标机轮减速率生成目标机轮轮速；
(6)根据所述目标机轮轮速和所述机轮轮速，生成防滑控制量并输出给飞机刹车系统。
根据本公开的至少一个实施方式，用所述调节量对所述初始的刹车系统关键参数进行调节后，生成的刹车系统关键参数替代所述步骤(2)中的初始的刹车系统关键参数。
根据本公开的至少一个实施方式，基于所述调节量对由于刹车系统关键参数改变而导致的刹车压力变化进行补偿。
根据本公开的至少一个实施方式，根据所述刹车结合力的变化率和所述滑移率的变化率生成目标机轮减速率的具体步骤为：
根据设定的机轮减速率上限值、初始的指令-减速率比例系数以及飞行员刹车压力指令，判断目标机轮减速率是否超出所述设定的机轮减速率上限值；
如果目标机轮减速率超出所述设定的机轮减速率上限值，则判断所述刹车结合力的变化率和所述滑移率的变化率的比值是否大于δ，如果大于，则输出目标机轮减速率，并令所述设定的机轮减速率上限值等于所述目标机轮减速率，如果不大于，则输出刹车压力与刹车系统关键参数的乘积作为目标机轮减速率；
如果目标机轮减速率未超出所述设定的机轮减速率上限值，则判断所述刹车结合力的变化率和所述滑移率的变化率的比值是否大于δ，如果大于，则输出目标机轮减速率，如果不大于，如果不大于，则输出刹车压力与刹车系统关键参数的乘积作为目标机轮减速率；
其中，δ为刹车过程中用于判定是否找到路面最大摩擦力点的阈值。
根据本公开的至少一个实施方式，将所述防滑控制量连同刹车指令压力共同作为总输出，输出给所述飞机刹车系统。
另一方面，本公开提供一种基于路况识别的飞机防滑刹车控制系统，包括：
系统初始参数估计模块，根据实时采集的刹车压力和实时采集的机轮轮速，获得初始的刹车系统关键参数；
系统状态观测模块，根据初始的刹车系统关键参数和所述刹车压力，获得刹车结合力的变化率和滑移率的变化率，根据初始的刹车系统关键参数获得估计的机轮轮速；
系统参数自适应模块，获得所述机轮轮速和所述估计的机轮轮速的差值，即实际轮速和估计轮速的差值；根据所述实际轮速和估计轮速的差值以及刹车压力，获得刹车系统关键参数的调节量；刹车系统关键参数等于所述初始的刹车系统关键参数与所述调节量之和；
目标机轮减速率生成模块，根据所述刹车结合力的变化率和所述滑移率的变化率生成目标机轮减速率；
目标机轮轮速生成模块，根据所述目标机轮减速率生成目标机轮轮速；以及
机轮轮速控制模块，根据所述目标机轮轮速和所述机轮轮速，生成防滑控制量并输出给飞机刹车系统。
根据本公开的至少一个实施方式，用所述调节量对所述初始的刹车系统关键参数进行调节后，生成的刹车系统关键参数替代输入到所述系统状态观测模块中的初始的刹车系统关键参数。
根据本公开的至少一个实施方式，所述机轮轮速控制模块还基于所述调节量对由于刹车系统关键参数改变而导致的刹车压力变化进行补偿。
根据本公开的至少一个实施方式，所述目标机轮减速率生成模块根据所述刹车结合力的变化率和所述滑移率的变化率生成目标机轮减速率的具体步骤为：
根据设定的机轮减速率上限值、初始的指令-减速率比例系数以及飞行员刹车压力指令，判断目标机轮减速率是否超出所述设定的机轮减速率上限值；
如果目标机轮减速率超出所述设定的机轮减速率上限值，则判断所述刹车结合力的变化率和所述滑移率的变化率的比值是否大于δ，如果大于，则输出目标机轮减速率，并令所述设定的机轮减速率上限值等于所述目标机轮减速率，如果不大于，则输出刹车压力与刹车系统关键参数的乘积作为目标机轮减速率；
如果目标机轮减速率未超出所述设定的机轮减速率上限值，则判断所述刹车结合力的变化率和所述滑移率的变化率的比值是否大于δ，如果大于，则输出目标机轮减速率，如果不大于，如果不大于，则输出刹车压力与刹车系统关键参数的乘积作为目标机轮减速率；
其中，δ为刹车过程中用于判定是否找到路面最大摩擦力点的阈值。
根据本公开的至少一个实施方式，将所述防滑控制量连同刹车指令压力共同作为总输出，输出给所述飞机刹车系统。
本公开的上述方法和系统仅用有限实时采集的刹车压力信号及轮速信号，在系统参数估计和状态观测模块中通过合理的简化条件(忽略空气阻力与升力，初始刹车阶段刹车压力较小致使机轮减速率较小)，实现了对刹车过程关键参数K的初始化(K＝Kb/M,其中Kb为刹车盘力矩系数，M为飞机质量)。在前述条件下可得到机速减速率与轮速减速率近似相等，进而通过积分轮速减速率得到机速，利用观测得到的机速与采集到的轮速做计算得到滑移率，直接数值微分得到dλ；再利用机轮动力学方程结合初始阶段的Kb并进行数值微分得到刹车结合力的变化率dF，通过计算dλ与dF比的绝对值作为识别路况，进而控制机轮得到跑道所能提供最大结合力刹车的依据。
另外由于在初始阶段刹车压力小致使滑移率的变化也小，所以将滑移率中机速与轮速分别由机速的变化率和轮速的变化率代替，并结合机轮动力学方程解得机轮轮速变化率，再对其进行数值积分可得到模型观测出的轮速值。根据观测的轮速与实际轮速做差通过类似M.I.T类型的调参方法生成系统参数K的调节量ΔK即可实现K的自适应调节；系统初始化后结合上述计算参数通过设定判断初始化最大机轮计算率，初始指令-减速率系数，减速率安全阈量三个关键参数生成目标减速率；最后通过观测轮速与实际轮速的轮速差经过控制器结合ΔK产生防滑控制量，实现不同路况下的自适应防滑控制。
本公开的优点在于：
1)通过刹车压力信号及采集的轮速信号即可对刹车系统关键参数K完成初始化，并设计控制器实现自适应校正。
2)可通过计算所得的系统关键参数K结合刹车压力实现对系统状态量dF,dλ，以及轮速的观测。
3)通过系统观测参数实现刹车系统在不同路况最佳滑移率下最优减速率的轮速防滑控制。
附图说明
附图示出了本公开的示例性实施方式，并与其说明一起用于解释本公开的原理，其中包括了这些附图以提供对本公开的进一步理解，并且附图包括在本说明书中并构成本说明书的一部分。
图1是本公开至少一个实施方式的单一跑道状态下机轮与跑道结合力与滑移率的关系曲线(图1中纵坐标F为结合力，横坐标λ为滑移率)。
图2是本公开至少一个实施方式的基于路况识别的飞机防滑刹车控制方法和控制系统的控制框图。
图3是本公开至少一个实施方式的目标机轮减速率生成模块流程图。
图4是本公开至少一个实施方式的轮速控制模块流程图。
具体实施方式
下面结合附图和实施方式对本公开作进一步的详细说明。可以理解的是，此处所描述的具体实施方式仅用于解释相关内容，而非对本公开的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本公开相关的部分。
需要说明的是，在不冲突的情况下，本公开中的实施方式及实施方式中的特征可以相互组合。下面将参考附图并结合实施方式来详细说明本公开。
如图3-4所示，本实施方式的基于路况识别的飞机防滑刹车控制方法，包括：
(1)根据实时采集的刹车压力bc和实时采集的机轮轮速ω，获得初始的刹车系统关键参数initK；
(2)根据初始的刹车系统关键参数initK和刹车压力bc，获得刹车结合力的变化率dF和滑移率的变化率dλ，根据初始的刹车系统关键参数initK获得估计的机轮轮速
Figure BDA0001806041760000061
(3)获得机轮轮速ω和估计的机轮轮速
Figure BDA0001806041760000062
的差值Δω，即实际轮速和估计轮速的差值；根据实际轮速和估计轮速的差值以及刹车压力bc，获得刹车系统关键参数K的调节量ΔK；刹车系统关键参数K等于初始的刹车系统关键参数initK与调节量ΔK之和；
(4)根据刹车结合力的变化率dF和滑移率的变化率dλ生成目标机轮减速率
Figure BDA0001806041760000063
(5)根据目标机轮减速率
Figure BDA0001806041760000064
生成目标机轮轮速ωd；
(6)根据目标机轮轮速ωd和机轮轮速ω，生成防滑控制量Abc并输出给飞机刹车系统。
在本实施方式中，用调节量ΔK对初始的刹车系统关键参数initK进行调节后，生成的刹车系统关键参数K替代步骤(2)中的初始的刹车系统关键参数initK。
在本实施方式中，基于调节量ΔK对由于刹车系统关键参数K改变而导致的刹车压力bc变化进行补偿。
在本实施方式中，根据刹车结合力的变化率dF和滑移率的变化率dλ生成目标机轮减速率
Figure BDA0001806041760000065
的具体步骤为：
根据设定的机轮减速率上限值、初始的指令-减速率比例系数以及飞行员刹车压力指令，判断目标机轮减速率
Figure BDA0001806041760000071
是否超出设定的机轮减速率上限值；
如果目标机轮减速率
Figure BDA0001806041760000072
超出设定的机轮减速率上限值，则判断刹车结合力的变化率dF和滑移率的变化率dλ的比值是否大于δ，如果大于，则输出目标机轮减速率
Figure BDA0001806041760000073
并令设定的机轮减速率上限值等于目标机轮减速率
Figure BDA0001806041760000074
如果不大于，则输出刹车压力bc与刹车系统关键参数K的乘积作为目标机轮减速率
Figure BDA0001806041760000075
如果目标机轮减速率
Figure BDA0001806041760000076
未超出设定的机轮减速率上限值，则判断刹车结合力的变化率dF和滑移率的变化率dλ的比值是否大于δ，如果大于，则输出目标机轮减速率
Figure BDA0001806041760000077
如果不大于，如果不大于，则输出刹车压力bc与刹车系统关键参数K的乘积作为目标机轮减速率
Figure BDA0001806041760000078
其中，δ为刹车过程中用于判定是否找到路面最大摩擦力点的阈值。
在本实施方式中，将防滑控制量Abc连同刹车指令压力Pbc共同作为总输出bc，输出给飞机刹车系统。
本实施方式的基于路况识别的飞机防滑刹车控制系统，包括：
系统初始参数估计模块，根据实时采集的刹车压力bc和实时采集的机轮轮速ω，获得初始的刹车系统关键参数initK；
系统状态观测模块，根据初始的刹车系统关键参数initK和刹车压力bc，获得刹车结合力的变化率dF和滑移率的变化率dλ，根据刹车系统关键参数initK获得估计的机轮轮速
Figure BDA0001806041760000079
系统参数自适应模块，获得机轮轮速ω和估计的机轮轮速
Figure BDA00018060417600000710
的差值Δω，即实际轮速和估计轮速的差值；根据实际轮速和估计轮速的差值以及刹车压力bc，获得刹车系统关键参数K的调节量ΔK；刹车系统关键参数K等于初始的刹车系统关键参数initK与调节量ΔK之和；
目标机轮减速率生成模块，根据刹车结合力的变化率dF和滑移率的变化率dλ生成目标机轮减速率
Figure BDA00018060417600000711
目标机轮轮速生成模块，根据目标机轮减速率
Figure BDA00018060417600000712
生成目标机轮轮速ωd；以及
机轮轮速控制模块，根据目标机轮轮速ωd和机轮轮速ω，生成防滑控制量Abc并输出给飞机刹车系统。
在本实施方式中，用调节量ΔK对初始的刹车系统关键参数initK进行调节后，生成的刹车系统关键参数K替代输入到系统状态观测模块中的初始的刹车系统关键参数initK。
在本实施方式中，机轮轮速控制模块还基于调节量ΔK对由于刹车系统关键参数K改变而导致的刹车压力bc变化进行补偿。
在本实施方式中，目标机轮减速率生成模块根据刹车结合力的变化率dF和滑移率的变化率dλ生成目标机轮减速率
Figure BDA0001806041760000081
的具体步骤为：
根据设定的机轮减速率上限值、初始的指令-减速率比例系数以及飞行员刹车压力指令，判断目标机轮减速率
Figure BDA0001806041760000082
是否超出设定的机轮减速率上限值；
如果目标机轮减速率
Figure BDA0001806041760000083
超出设定的机轮减速率上限值，则判断刹车结合力的变化率dF和滑移率的变化率dλ的比值是否大于δ，如果大于，则输出目标机轮减速率
Figure BDA0001806041760000084
并令设定的机轮减速率上限值等于目标机轮减速率
Figure BDA0001806041760000085
如果不大于，则输出刹车压力bc与刹车系统关键参数K的乘积作为目标机轮减速率
Figure BDA0001806041760000086
如果目标机轮减速率
Figure BDA0001806041760000087
未超出设定的机轮减速率上限值，则判断刹车结合力的变化率dF和滑移率的变化率dλ的比值是否大于δ，如果大于，则输出目标机轮减速率
Figure BDA0001806041760000088
如果不大于，则输出刹车压力bc与刹车系统关键参数K的乘积作为目标机轮减速率
Figure BDA0001806041760000089
其中，δ为刹车过程中用于判定是否找到路面最大摩擦力点的阈值。
在本实施方式中，将防滑控制量Abc连同刹车指令压力Pbc共同作为总输出bc，输出给飞机刹车系统。
更详细的，结合图1对刹车问题进行描述：
由跑道与机轮结合系数影响因素，在单一跑道状态下结合力与滑移率的趋势关系如图1所示，在单一跑道状态下利用跑道所提供的最大结合力进行刹车，也即找曲线最高点所对应的滑移率，并控制刹车压力使得滑移率达到这一最大滑移率。路况发生改变或轮载发生变化，则曲线会发生跳变，需动态的平稳的追踪当前结合状态下的最大结合力，在这一过程中，机轮应尽可能不出现打滑状态。
另外，在此状态下，观察曲线形状可以得到下面的结论。当跑道状态不发生改变时，机轮工作点在曲线上来回运动。
在最高点的左侧运动时，无论工作点向上运动还是向下运动(1或4)，dF/dλ>0恒成立。反之，工作点在最高点右侧运动时，不论是向上运动还是向下运动(2或3)，dF/dλ<0恒成立。并且越接近最高点，dF/dλ的绝对值越小。根据这一特征，可以进行控制律设计，使得机轮与地面的结合力接近跑道所能提供的最大结合力，从而提高刹车系统效率。
若要进行滑移状态的判定，需要知道dF/dλ的值的大小。对于整个刹车系统来说，输入参数只有刹车指令和机轮转速两个参数，要得到结合力的变化率dF和滑移率的变化率dλ必须对输入的数据进行处理，利用算法估计出所需要的量进行检测和控制。
本公开的基于路况识别的飞机防滑刹车自适应控制方法和控制系统总体控制框图2所示，整个防滑刹车自适应控制系统有一个根据路况判定目标机轮减速率的生成器，目标机轮减速率生成器在机轮处于防滑状态时，输出目标的减速率
Figure BDA0001806041760000091
该减速率输入到下一单元目标轮速生成模块中去，该模块的作用是将输入的目标机轮减速率结合初始速度生成目标轮速，将目标轮速与实际轮速作差得到的误差值输入下一级的轮速控制器中去，由轮速控制器产生防滑量Abc，连同飞行员给出的刹车指令压力Pbc一起，共同作为总输出bc输出给刹车系统。再将总输出bc结合机轮轮速进行系统初始参数估计、系统参数自适应、系统状态观测，用于初始系统参数initK与参数模型观测轮速
Figure BDA0001806041760000092
以及dF与dλ的估计。
总体而言，整个控制系统分为3个部分，系统参数估计模块和系统状态观测模块，目标轮速生成部分，目标机轮减速率生成模块与轮速控制模块，下面分别对3个部分进行详细介绍。
一、系统参数估计模块和系统状态观测模块，功能是通过输入轮速和刹车指令，观测并自适应得到刹车系统关键参数K以及估计dF与dλ两个量的大小。主要包括如下几个部分：
(1)系统初始状态估计：
利用在初始刹车状态下系统参数如刹车盘结合系数等变化不大的情况下，进行系统参数的初步近似估计，用于路况检测。
由于机轮方程：
Figure BDA0001806041760000101
式中J为机轮转动惯量，ω为机轮转动角速度，Ff为结合力，r为机轮半径，Tb为刹车力矩，μ为轮胎与跑道表面结合系数，FN为作用在机轮上的垂直载荷。便于分析简单，现将空气阻力与升力忽略(加上空气阻力和升力计算方式不变)：
Figure BDA0001806041760000102
式中Vp为飞机速度，现将式1，2联立得：
Figure BDA0001806041760000103
式中bc·Kb为1式中刹车力矩Tb，Kb为刹车盘力矩系数，bc为刹车压力。在初始开始刹车时，由于刹车压力很小，即bc很小，所以导致
Figure BDA0001806041760000104
很小，又由于机轮转动惯量很小，
Figure BDA0001806041760000105
和bcKb两项差别可达到上百倍，故忽略
Figure BDA0001806041760000106
此时有：
Figure BDA0001806041760000107
即：
Figure BDA0001806041760000108
令系统参数
Figure BDA0001806041760000109
则：
Figure BDA00018060417600001010
此时，由于刹车压力小，滑移率也较小，此时可得
Figure BDA00018060417600001011
则代入反求系统参数K，在初始状态下，可得：
Figure BDA00018060417600001012
到此可得到初始参数initK。
(2)系统状态观测：
该部分主要目的是利用框图中输入的系统参数K，计算估计得到dF与dλ，以及系统根据系统参数K反算出的轮速
Figure BDA00018060417600001013
根据上面描述，在初始阶段，在打滑未发生时，由于系统压力较小，满足两个条件，即：
Figure BDA0001806041760000111
和bc·Kb两项差别可达到上百倍，且
Figure BDA0001806041760000112
知:
Figure BDA0001806041760000113
则可直接利用积分积出飞机速度Vp，在计算得到Vp后，由于轮速可以采集，可以由公式：
Figure BDA0001806041760000114
可计算滑移率λ，直接微分即可得到滑移率微分dλ。又由于：
Figure BDA0001806041760000115
根据实验结果，可观察出在初始阶段的Kb，则给定一个初始值Kb，即可进行dF的计算，在此可以指出，Kb的初值对dF的变化影响不大。到此即求出dF的值，此处及上处微分的处理均采用微分跟踪器(TD)，在此不再赘述。
由于：
Figure BDA0001806041760000116
在初始阶段，压力小，所以滑移率的变化也很小，此时有：
Figure BDA0001806041760000117
将式13带入式3进行整理有：
Figure BDA0001806041760000118
则，上下同时除以M得：
Figure BDA0001806041760000119
由于J/M非常小，所以将之忽略，得到
Figure BDA00018060417600001110
此处的滑移率λ可由上面求dλ时产生的λ代替，将上述值进行积分，即可得到由模型观测出的轮速
Figure BDA00018060417600001111
(3)系统参数自适应：
该部分的输入为
Figure BDA0001806041760000121
即估计轮速和实际轮速的差值，此处导致估计轮速和实际轮速有差的原因在很大程度上可归结为系统参数K的变化(在没有打滑的情况下)。故可设计自适应律，根据轮速的差值Δω，生成系统参数调节量ΔK，并同初始系统参数initK一起生成系统参数。
在这里，系统参数调节量ΔK的调节方式采样类似M.I.T类型的调参方法，将轮速的差值直接作为系统参数的输入，该误差越大，则说明K参数应当向误差变化相反的方向发展，误差越大，则K的调整量应当越多。为了防止参数最终出现抖动的现象，在此采用PID控制器，将轮速误差Δω作为控制器的输入，PID的控制器输出作为K参数的调整量ΔK。有：
ΔK＝Kp[Δω+1/KI∫Δωdt+KddΔω/dt] (17)
此时：
K＝initK+ΔK (18)
如此即达到了可以利用轮速自适应系统参数的目的。
二、目标机轮减速率生成模块，功能：通过输入的由输入结合力变化率和滑移率变化率，生成目标机轮减速率。
(1)在进行刹车控制前，需要设定3个参数，dωmax0，K0，dωs，含义分别为初始最大机轮减速率，初始指令-减速率比例系数，减速率安全阈量。
最大机轮减速率：设定的机轮最大减速率。
指令-减速率比例系数：刹车指令与减速率之间的线性对应系数，K0＝dωmax/bcmax。
减速率安全阈量：即将出现打滑时机轮减速率指令的减量。
在开始刹车时，即启动判定dF/dλ>δ若飞行员指令给到最大时，此判定仍然为真，那么防滑控制不再输出控制，输出使能信号为0，直到该判定条件为假。
δ为大于0但接近0的数，当条件dF/dλ<δ满足时，判定达到了路面最大摩擦力点，并进行目标减速率设定等一系列步骤，当条件不满足时，判定还未找到最大摩擦力点，压力仍需增加。
在开始刹车后，如果满足了该条件为假，那么记录此时得刹车压力指令为bcmax，即压力最大值。在整个轮速控制期间，bc必须都保证小于bcmax以保证机轮不出现打滑现象。此时，由于未发生打滑，则可以根据：
Figure BDA0001806041760000131
Figure BDA0001806041760000132
两式联立得到的：
Figure BDA0001806041760000133
并带入bcmax得到dωmax:
Figure BDA0001806041760000134
作为dωmax0。并记录此时的dc。此时K＝dωmax0/dc。dωs＝k*dωmax0。式中k为一安全系数，取值小于1一般取0.1。
进行完上述判定及赋值之后，开始如图3流程。
根据输入的指令，根据初始比例分配目标机轮减速率，判定目标机轮减速率是否超过最大值，若判定为真：
则继续判定dF/dλ>δ是否为真。按照图3中所述更新变量值，并输出机轮减速率作为模块输出。
若判定为假，则按照图3更新变量，输出机轮减速率同上。二者的最大区别在于在判定dF/dλ>δ为真的情况下是否要进行最大机轮减速率更新。
三、轮速控制模块，功能：输入机轮减速率，控制机轮转速，输出防滑量。结构框图如图4。
该模块输入为上个模块产生的目标机轮减速率与轮速采样，输出为防滑控制量，同时防滑控制量也是该系统的总输出。
该模块的作用是控制轮速ω到预期轮速ωd。当前轮速减去机轮减速率乘以控制周期可以得到所需的目标轮速，将这一轮速作为控制目标进行控制即可。
Figure BDA0001806041760000141
e(k)＝ωd-ω (24)
此处采用PID控制器：
u(t)＝kp[e(t)+1/TI∫e(t)dt+TDde(t)/dt] (25)
此处未避免系统发散不设置D参数。
打滑检测的部分也在此部分进行。在此通过上述步骤产生的滑移率和机轮减速率进行处理，调节参数kp以调节机轮打滑状态下松压程度的大小，在检测到λ＞Kλ·λmax且dω＜Kω·dωmax时，输出Δdc＝kp(dF/dλ)与dc做和运算即：dc＝dc+Δdc，Δdc＝kp(dF/dλ)。
除此之外，此模块输入还有ΔK，用于做由于系统参数改变导致的压力补偿，由于：
Figure BDA0001806041760000142
在刹车过程中K的改变其实也反应出了Kb的变化，因为整个刹车过程中，飞机质量基本是不发生改变的，改变的只是刹车盘结合系数Kb。根据K的改变的大小，乘以相应的系数Kp，即可进行压力补偿。
本领域的技术人员应当理解，上述实施方式仅仅是为了清楚地说明本公开，而并非是对本公开的范围进行限定。对于所属领域的技术人员而言，在上述公开的基础上还可以做出其它变化或变型，并且这些变化或变型仍处于本公开的范围内。

Claims (10)

1.一种基于路况识别的飞机防滑刹车控制方法，其特征在于，包括：
（1）根据实时采集的刹车压力
Figure 880277DEST_PATH_IMAGE001
和实时采集的机轮轮速
Figure 141363DEST_PATH_IMAGE002
，获得初始的刹车系统关键参数
Figure 699383DEST_PATH_IMAGE003
;
（2）根据初始的刹车系统关键参数
Figure 87639DEST_PATH_IMAGE003
和所述刹车压力
Figure 47505DEST_PATH_IMAGE001
，获得刹车结合力的变化率
Figure 230224DEST_PATH_IMAGE004
和滑移率的变化率
Figure 541120DEST_PATH_IMAGE005
，根据初始的刹车系统关键参数
Figure 218220DEST_PATH_IMAGE003
获得估计的机轮轮速
Figure 32592DEST_PATH_IMAGE006
；
（3）获得所述机轮轮速
Figure 651792DEST_PATH_IMAGE002
和所述估计的机轮轮速
Figure 449984DEST_PATH_IMAGE007
的差值Δω，即实际轮速和估计轮速的差值；根据所述实际轮速和估计轮速的差值以及刹车压力
Figure 648884DEST_PATH_IMAGE001
，获得刹车系统关键参数
Figure 98189DEST_PATH_IMAGE008
的调节量Δ
Figure 622711DEST_PATH_IMAGE008
；刹车系统关键参数
Figure 642620DEST_PATH_IMAGE008
等于所述初始的刹车系统关键参数
Figure 910790DEST_PATH_IMAGE003
与所述调节量Δ
Figure 965334DEST_PATH_IMAGE008
之和；
（4）根据所述刹车结合力的变化率
Figure 145911DEST_PATH_IMAGE009
和所述滑移率的变化率
Figure 918695DEST_PATH_IMAGE005
生成目标机轮减速率
Figure 724977DEST_PATH_IMAGE010
；
（5）根据所述目标机轮减速率
Figure 634027DEST_PATH_IMAGE011
生成目标机轮轮速
Figure 500352DEST_PATH_IMAGE012
；
（6）根据所述目标机轮轮速
Figure 744120DEST_PATH_IMAGE012
和所述机轮轮速
Figure 619672DEST_PATH_IMAGE002
，生成防滑控制量
Figure 383229DEST_PATH_IMAGE013
并输出给飞机刹车系统；
其中，刹车系统关键参数
Figure 420455DEST_PATH_IMAGE008
为：
K=Kb/M,
式中，Kb为刹车盘力矩系数，且Kb=Tb/bc，Tb为刹车力矩，bc为刹车压力；M为飞机质量；
初始的刹车系统关键参数
Figure 902252DEST_PATH_IMAGE003
为：
Figure 66648DEST_PATH_IMAGE014
式中，
Figure 950290DEST_PATH_IMAGE016
为机轮转动角速度，r为机轮半径，
Figure 892839DEST_PATH_IMAGE001
为刹车压力；
所述调节量Δ
Figure 861932DEST_PATH_IMAGE008
的获得过程为，采用PID控制器，将轮速误差Δω作为控制器的输入，PID的控制器输出作为K参数的调整量Δ
Figure 813707DEST_PATH_IMAGE008
，有：
Figure 801124DEST_PATH_IMAGE017
。
2.根据权利要求1所述的基于路况识别的飞机防滑刹车控制方法，其特征在于，用所述调节量Δ
Figure 445732DEST_PATH_IMAGE008
对所述初始的刹车系统关键参数
Figure 636541DEST_PATH_IMAGE003
进行调节后，生成的刹车系统关键参数
Figure 392008DEST_PATH_IMAGE008
替代所述步骤（2）中的初始的刹车系统关键参数
Figure 984663DEST_PATH_IMAGE003
。
3.根据权利要求1或2所述的基于路况识别的飞机防滑刹车控制方法，其特征在于，基于所述调节量Δ
Figure 550905DEST_PATH_IMAGE008
对由于刹车系统关键参数
Figure 963432DEST_PATH_IMAGE008
改变而导致的刹车压力
Figure 522589DEST_PATH_IMAGE001
变化进行补偿。
4.根据权利要求1或2所述的基于路况识别的飞机防滑刹车控制方法，其特征在于，根据所述刹车结合力的变化率
Figure 969751DEST_PATH_IMAGE004
和所述滑移率的变化率
Figure 956161DEST_PATH_IMAGE005
生成目标机轮减速率
Figure 121563DEST_PATH_IMAGE011
的具体步骤为：
根据设定的机轮减速率上限值、初始的指令-减速率比例系数以及飞行员刹车压力指令，判断目标机轮减速率
Figure 468100DEST_PATH_IMAGE011
是否超出所述设定的机轮减速率上限值；
如果目标机轮减速率
Figure 35348DEST_PATH_IMAGE011
超出所述设定的机轮减速率上限值，则判断所述刹车结合力的变化率
Figure 927080DEST_PATH_IMAGE004
和所述滑移率的变化率
Figure 579778DEST_PATH_IMAGE005
的比值是否大于
Figure 215159DEST_PATH_IMAGE018
，如果大于，则输出目标机轮减速率
Figure 653225DEST_PATH_IMAGE011
，并令所述设定的机轮减速率上限值等于所述目标机轮减速率
Figure 715859DEST_PATH_IMAGE011
，如果不大于，则输出刹车压力
Figure 855853DEST_PATH_IMAGE001
与刹车系统关键参数
Figure 29346DEST_PATH_IMAGE008
的乘积作为目标机轮减速率
Figure 571185DEST_PATH_IMAGE011
；
如果目标机轮减速率
Figure 53988DEST_PATH_IMAGE011
未超出所述设定的机轮减速率上限值，则判断所述刹车结合力的变化率
Figure 681279DEST_PATH_IMAGE004
和所述滑移率的变化率
Figure 658462DEST_PATH_IMAGE005
的比值是否大于
Figure 320387DEST_PATH_IMAGE018
，如果大于，则输出目标机轮减速率
Figure 724824DEST_PATH_IMAGE011
，如果不大于，则输出刹车压力
Figure 324564DEST_PATH_IMAGE001
与刹车系统关键参数
Figure 105438DEST_PATH_IMAGE008
的乘积作为目标机轮减速率
Figure 621870DEST_PATH_IMAGE011
；
其中，
Figure 197208DEST_PATH_IMAGE018
为刹车过程中用于判定是否找到路面最大摩擦力点的阈值。
5.根据权利要求1所述的基于路况识别的飞机防滑刹车控制方法，其特征在于，将所述防滑控制量
Figure 533511DEST_PATH_IMAGE019
连同刹车指令压力
Figure 852497DEST_PATH_IMAGE020
共同作为总输出
Figure 472703DEST_PATH_IMAGE001
，输出给所述飞机刹车系统。
6.一种基于路况识别的飞机防滑刹车控制系统，其特征在于，包括：
系统初始参数估计模块，根据实时采集的刹车压力
Figure 484521DEST_PATH_IMAGE001
和实时采集的机轮轮速
Figure 308121DEST_PATH_IMAGE021
，获得初始的刹车系统关键参数
Figure 430798DEST_PATH_IMAGE022
；
系统状态观测模块，根据初始的刹车系统关键参数
Figure 672554DEST_PATH_IMAGE022
和所述刹车压力
Figure 589695DEST_PATH_IMAGE001
，获得刹车结合力的变化率
Figure 635011DEST_PATH_IMAGE004
和滑移率的变化率
Figure 561379DEST_PATH_IMAGE005
，根据初始的刹车系统关键参数
Figure 906909DEST_PATH_IMAGE023
获得估计的机轮轮速
Figure 994951DEST_PATH_IMAGE007
；
系统参数自适应模块，获得所述机轮轮速
Figure 42410DEST_PATH_IMAGE021
和所述估计的机轮轮速
Figure 506890DEST_PATH_IMAGE007
的差值Δω，即实际轮速和估计轮速的差值；根据所述实际轮速和估计轮速的差值以及刹车压力
Figure 706927DEST_PATH_IMAGE001
，获得刹车系统关键参数
Figure 965870DEST_PATH_IMAGE008
的调节量Δ
Figure 251358DEST_PATH_IMAGE008
；刹车系统关键参数
Figure 270261DEST_PATH_IMAGE008
等于所述初始的刹车系统关键参数
Figure 324804DEST_PATH_IMAGE022
与所述调节量Δ
Figure 754649DEST_PATH_IMAGE008
之和；
目标机轮减速率生成模块，根据所述刹车结合力的变化率
Figure 527433DEST_PATH_IMAGE004
和所述滑移率的变化率
Figure 333714DEST_PATH_IMAGE005
生成目标机轮减速率
Figure 492032DEST_PATH_IMAGE011
；
目标机轮轮速生成模块，根据所述目标机轮减速率
Figure 358357DEST_PATH_IMAGE011
生成目标机轮轮速
Figure 352858DEST_PATH_IMAGE012
；以及
机轮轮速控制模块，根据所述目标机轮轮速
Figure 962831DEST_PATH_IMAGE012
和所述机轮轮速
Figure 726387DEST_PATH_IMAGE021
，生成防滑控制量
Figure 779925DEST_PATH_IMAGE019
并输出给飞机刹车系统；
其中，刹车系统关键参数
Figure 996143DEST_PATH_IMAGE008
为：
K=Kb/M
式中，Kb为刹车盘力矩系数，且Kb=Tb/bc，Tb为刹车力矩，bc为刹车压力；M为飞机质量；
初始的刹车系统关键参数
Figure 409807DEST_PATH_IMAGE022
为：
Figure 27870DEST_PATH_IMAGE014
式中，
Figure 501577DEST_PATH_IMAGE016
为机轮转动角速度，r为机轮半径，
Figure 205090DEST_PATH_IMAGE001
为刹车压力；
所述调节量Δ
Figure 507650DEST_PATH_IMAGE008
的获得过程为，采用PID控制器，将轮速误差Δω作为控制器的输入，PID的控制器输出作为K参数的调整量Δ
Figure 980220DEST_PATH_IMAGE008
，有：
Figure 624828DEST_PATH_IMAGE017
。
7.根据权利要求6所述的基于路况识别的飞机防滑刹车控制系统，其特征在于，用所述调节量Δ
Figure 815638DEST_PATH_IMAGE008
对所述初始的刹车系统关键参数
Figure 571104DEST_PATH_IMAGE022
进行调节后，生成的刹车系统关键参数
Figure 180071DEST_PATH_IMAGE008
替代输入到所述系统状态观测模块中的初始的刹车系统关键参数
Figure 730001DEST_PATH_IMAGE023
。
8.根据权利要求6或7所述的基于路况识别的飞机防滑刹车控制系统，其特征在于，所述机轮轮速控制模块还基于所述调节量Δ
Figure 408107DEST_PATH_IMAGE008
对由于刹车系统关键参数
Figure 701685DEST_PATH_IMAGE008
改变而导致的刹车压力
Figure 414426DEST_PATH_IMAGE001
变化进行补偿。
9.根据权利要求6或7所述的基于路况识别的飞机防滑刹车控制系统，其特征在于，所述目标机轮减速率生成模块根据所述刹车结合力的变化率
Figure 650104DEST_PATH_IMAGE009
和所述滑移率的变化率
Figure 549927DEST_PATH_IMAGE005
生成目标机轮减速率
Figure 647196DEST_PATH_IMAGE011
的具体步骤为：
根据设定的机轮减速率上限值、初始的指令-减速率比例系数以及飞行员刹车压力指令，判断目标机轮减速率
Figure 214444DEST_PATH_IMAGE011
是否超出所述设定的机轮减速率上限值；
如果目标机轮减速率
Figure 371756DEST_PATH_IMAGE011
超出所述设定的机轮减速率上限值，则判断所述刹车结合力的变化率
Figure 775186DEST_PATH_IMAGE004
和所述滑移率的变化率
Figure 410567DEST_PATH_IMAGE005
的比值是否大于
Figure 832321DEST_PATH_IMAGE018
，如果大于，则输出目标机轮减速率
Figure 426113DEST_PATH_IMAGE011
，并令所述设定的机轮减速率上限值等于所述目标机轮减速率
Figure 300529DEST_PATH_IMAGE011
，如果不大于，则输出刹车压力
Figure 988868DEST_PATH_IMAGE001
与刹车系统关键参数
Figure 265128DEST_PATH_IMAGE008
的乘积作为目标机轮减速率
Figure 764243DEST_PATH_IMAGE011
；
如果目标机轮减速率
Figure 125954DEST_PATH_IMAGE011
未超出所述设定的机轮减速率上限值，则判断所述刹车结合力的变化率
Figure 368716DEST_PATH_IMAGE004
和所述滑移率的变化率
Figure 250216DEST_PATH_IMAGE005
的比值是否大于
Figure 920232DEST_PATH_IMAGE018
，如果大于，则输出目标机轮减速率
Figure 769239DEST_PATH_IMAGE011
，如果不大于，则输出刹车压力
Figure 815692DEST_PATH_IMAGE001
与刹车系统关键参数
Figure 66545DEST_PATH_IMAGE008
的乘积作为目标机轮减速率
Figure 891151DEST_PATH_IMAGE011
；
其中，
Figure 227454DEST_PATH_IMAGE018
为刹车过程中用于判定是否找到路面最大摩擦力点的阈值。
10.根据权利要求6所述的基于路况识别的飞机防滑刹车控制系统，其特征在于，将所述防滑控制量
Figure 812019DEST_PATH_IMAGE019
连同刹车指令压力
Figure 182958DEST_PATH_IMAGE020
共同作为总输出
Figure 929197DEST_PATH_IMAGE024
，输出给所述飞机刹车系统。
CN201811098223.0A 2018-09-20 2018-09-20 基于路况识别的飞机防滑刹车自适应控制方法 Active CN109367765B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811098223.0A CN109367765B (zh) 2018-09-20 2018-09-20 基于路况识别的飞机防滑刹车自适应控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811098223.0A CN109367765B (zh) 2018-09-20 2018-09-20 基于路况识别的飞机防滑刹车自适应控制方法

Publications (2)

Publication Number Publication Date
CN109367765A CN109367765A (zh) 2019-02-22
CN109367765B true CN109367765B (zh) 2021-03-05

Family

ID=65405321

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811098223.0A Active CN109367765B (zh) 2018-09-20 2018-09-20 基于路况识别的飞机防滑刹车自适应控制方法

Country Status (1)

Country Link
CN (1) CN109367765B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112623199B (zh) * 2020-12-09 2022-08-05 北京北摩高科摩擦材料股份有限公司 一种飞机自适应防滑刹车控制方法
CN113044207B (zh) * 2020-12-29 2022-09-20 中国航空工业集团公司西安飞机设计研究所 一种提高飞机防滑刹车系统安全性的方法
CN113246934B (zh) * 2021-05-11 2022-06-28 北京航空航天大学 一种基于自适应跑道的双轮控制方法、设备及存储器
CN113879518A (zh) * 2021-11-19 2022-01-04 中国商用飞机有限责任公司 防滑控制方法、装置、电子设备和可读存储介质
CN115783252B (zh) * 2023-02-08 2023-04-25 四川腾盾科技有限公司 一种低成本的无人机自主刹车压力控制计算方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102092373B (zh) * 2010-12-16 2013-09-04 西安航空制动科技有限公司 一种飞机自动刹车方法及装置
CN102991488B (zh) * 2012-11-26 2014-08-13 西安航空制动科技有限公司 具有自适应能力的刹车系统恒力矩控制方法
CN103612750A (zh) * 2013-11-20 2014-03-05 渭南高新区晨星专利技术咨询有限公司 一种飞机防滑刹车控制方法
US9975628B2 (en) * 2015-12-28 2018-05-22 Goodrich Corporation Anti-skid protection with undetected pressure sensor failure
CN105752325B (zh) * 2016-03-21 2017-07-14 北京航空航天大学 基于刹车力矩反馈的飞机防滑刹车控制方法
CN106515695B (zh) * 2016-12-26 2019-04-30 北京航空航天大学 基于摩擦特性的防滑刹车控制方法

Also Published As

Publication number Publication date
CN109367765A (zh) 2019-02-22

Similar Documents

Publication Publication Date Title
CN109367765B (zh) 基于路况识别的飞机防滑刹车自适应控制方法
EP3453613B1 (en) System and method for aircraft electric taxi brake optimization
EP1732796B1 (en) Advanced braking system for aircraft
CN1086173C (zh) 用于自适配制动和起始滑动检测的系统和方法
US6890041B1 (en) Antilock brake systems employing a sliding mode observer based estimation of differential wheel torque
EP2920069B1 (en) Aircraft landing gear longitudinal force control
CA2656171C (en) Aircraft braking control
EP3363697B1 (en) Antiskid operation during degraded operation
CN109094772A (zh) 一种用于控制交通工具的制动的设备及飞行器
CN106828894B (zh) 基于轮速控制的飞机防滑刹车控制方法和系统
CN106828893B (zh) 基于结合力估计的防滑刹车控制方法和系统
JPH0114062B2 (zh)
CN106335634B (zh) 飞行器及其转向系统、转向方法和转向系统控制器
EP2998819B1 (en) Variable maximum commandable roll rate for directional control during engine-out rolling maneuver
EP3875323B1 (en) Systems and methods for aircraft antiskid braking
EP3653456B1 (en) System and method for rapid convergence antiskid initialization
GB2411934A (en) Braking system for aircraft
US11878789B2 (en) Aircraft brake control system
Bai et al. Neural network-based reference speed self-adjusting brake control algorithm
CN117799826A (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
TA01 Transfer of patent application right

Effective date of registration: 20191220

Address after: 100090 A018, floor 4, block a, Zhizhen building, No.7, Zhichun Road, Haidian District, Beijing

Applicant after: Beijing Hangchen System Technology Co., Ltd.

Address before: 100191 A018, 4th floor, Building A, Zhizhen Building, 7 Zhichun Road, Haidian District, Beijing

Applicant before: Beijing Hangchen System Technology Co., Ltd.

Applicant before: Shang Yaoxing

TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20210201

Address after: Room 101, unit 1, building 7, Fangshan science and Technology Park, Liandong u Valley, Yanxin South Road, Yanshan Petrochemical new material science and technology industrial base, Fangshan District, Beijing 102400

Applicant after: Beijing hangchen Airborne Intelligent System Technology Co.,Ltd.

Address before: 100090 A018, 4th floor, block a, Zhizhen building, No.7 Zhichun Road, Haidian District, Beijing

Applicant before: BEIJING HANGCHEN SYSTEM TECHNOLOGY Co.,Ltd.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant