CN111086520A - 一种适用于四轮驱动车辆多轮高滑移率时的速度估计算法 - Google Patents
一种适用于四轮驱动车辆多轮高滑移率时的速度估计算法 Download PDFInfo
- Publication number
- CN111086520A CN111086520A CN202010052297.1A CN202010052297A CN111086520A CN 111086520 A CN111086520 A CN 111086520A CN 202010052297 A CN202010052297 A CN 202010052297A CN 111086520 A CN111086520 A CN 111086520A
- Authority
- CN
- China
- Prior art keywords
- vehicle speed
- wheel
- speed
- slip
- estimating
- 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.)
- Pending
Links
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W40/00—Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
- B60W40/10—Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to vehicle motion
- B60W40/105—Speed
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W2520/00—Input parameters relating to overall vehicle dynamics
- B60W2520/26—Wheel slip
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Mathematical Physics (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Regulating Braking Force (AREA)
- Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
Abstract
本发明公开了一种适用于四轮驱动汽车多轮高滑移率时的速度估计算法,所述方法分两层,第一层包括:利用基于轮速的车速估算模块,将车轮角速度转换为车轮纵向速度并修正轮胎旋转半径和轮胎标记半径的偏差,根据模糊逻辑估计车速;利用基于加速度的车速估算模块,根据积分加速度变量来估计车速;多轮大滑移率时第一层算法无法准确估算车速,第二层算法中车速估算模块引入模糊逻辑以精确判断车轮打滑状况,并估算精确车速;按照第二层车速估算模块的加权系数,平衡车轮和加速度传感器的两种信号,估算最终车速。本方法引入了临时估算车速,以解决对车速和轮胎滑移率估算不准确的问题,本方法提高了精度,适用于四轮驱动车辆多轮高滑移率时估算速度。
Description
技术领域
发明涉及四轮驱动车辆状态估计技术领域,具体为一种适用于四轮驱动车辆多轮高滑移率时的速度估计算法。
背景技术
随着智能和先进的驾驶员辅助系统的发展,车辆状态估计对于汽车主动系统的意义越来越大,而车速是车辆状态估计的关键变量,多个车轮出现滑移时的速度估算是重点也是难点。
现有的技术通常针对低滑移率,在普通路面上进行车辆速度估计。然而,四驱车辆在雪地路面上出现高滑移率是更危险的情况,特别是四轮同时打滑。传统上,基于轮速传感器的纵向车速估算和基于加速度传感器的纵向车速估算都有缺点,基于轮速传感器的纵向车速估算在没有车轮滑移的情况下估算值是比较准确的;基于加速度传感器的车速估算在加速度比较小或是频繁加减速的情况下估算值精度差。另外,由于轮胎模型的误差而引起的动力学算法性能下降,也会影响车辆速度估计精度。
利用四轮编码器和加速度传感器,引入轮胎转动系数修正轮胎旋转半径与轮胎标记半径之间的偏差,再引入临时估算车速处理车速和轮胎滑移率计算不准确的问题,能在多轮高滑移率情况下准确估算车速。本发明在此基础上,采用二次模糊控制的方法提高最终估算车速的精度,提出了一种四驱车辆多轮高滑移率的速度估算方法。
发明内容
本发明的目的在于克服现有的技术缺陷,提供一种四轮驱动车辆多轮高滑移率时的速度估计算法。
为实现上述目的,可以通过一下技术方案来实现:一种四轮驱动车辆多轮高滑移率时的速度估计算法,包括轮速计算模块、基于轮速的车速估算模块、信号滤波模块、基于加速度的车速估算模块、滑移率置信度计算模块、车速估算模块和最终模糊控制模块。
具体包括以下步骤:1)根据第一层算法进行预估算,在大滑移率的情况下分别基于四个轮速传感器和纵向加速度传感器估算相对准确的纵向车速,作为参考值提供给第二层算法;2) 根据第二层算法进行最终的纵向车速估算,调节在不同车轮滑移情况下的置信程度,使得在大滑移率、多轮打滑的情况下得出精确的纵向车速估算值。
其中,所述的轮速计算模块接收传感器信号,信号为车轮角速度ωi,j,其数值会随实际车辆轮胎气压等变化而有偏差,因此引入轮胎转动系数η进行标定,再输出估算车轮线速度Vi,j给基于轮速的车速估算模块。
其中,所述的基于轮速的车速估算模块从传感器得到的四个轮速中选取最小速度,并计算四个车轮的滑移率,作为模糊控制的输入量,输出相对准确的车速估算值vE,w给基于加速度的车速估算模块和第二层算法。
其中,所述的信号滤波模块用于处理加速度传感器的信号,对车辆纵向加速度asensor进行滤波。
其中,所述的基于加速度的车速估算模块接收滤波后的加速度,采用小步长时间迭代积分,然后用线性补偿方式进行标定输出车速估算值VE a给第二层算法。
其中,所述的车速估算模块根据第一层算法得出的相对准确的估算车速,在考虑滑移率和滑移车轮个数基础上,输出估算值到滑移率置信度计算模块和最终模糊控制模块。
其中,所述的滑移率置信度计算模块利用置信参数平衡基于轮速的车速估算和基于加速度的车速估算中估算出的纵向车速,根据不同滑移率得出相对准确的纵向车速值VE temp,并重新计算出四个轮胎的滑移率输出到最终模糊控制模块。
与现有的技术相比,本发明具有以下优点:
本发明通过引入轮胎旋转系数,抵消轮胎旋转半径和标记轮胎半径之间的偏差,提高了基于轮速的车速估算模块的计算精度;并且通过噪音滤波处理加速度传感器的信号,再采用线性补偿方式消除估计纵向车速和实际车速之间的偏差,提高了车速估计精度。
本发明通过计算临时估算车速的方法来处理轮胎滑移率估算不准确的问题,临时车速估算值并不作为最终估算值,但作为参考估算滑移率非常必要,出现大滑移且有两个以上车轮滑移时,通过加速度积分和标定校正的车速估算值在车速较低且加速度变化较大时不精确,利用相对准确的临时估算车速再次计算轮胎滑移率,提高了滑移率估算精度。
本发明利用置信度参数平衡基于轮速的车速估算和基于加速度的车速估算,可仅仅使用低成本的速度编码器和加速度传感器,通过两次模糊控制得出最终车速估算值,大大减小了车速和滑移率误差,提高了控制精度。
附图说明
图1为本发明的算法总体流程图
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
如图1,本发明提供一种一种四轮驱动车辆多轮高滑移率时的速度估计算法,包括以下步骤:
1)在第一层算法中进行预估算,在大滑移率的情况下分别基于四个轮速传感器和纵向加速度传感器估算出相对准确的临时车速,作为第二层算法的参考值;
11)轮速计算模块根据传感器信号车轮角速度ωi,j,转化为车轮线速度Vi,j,引入轮胎转动系数η进行标定校正;
12)根据步骤11)得到的车轮速度Vi,j,基于轮速的车速估算模块从传感器得到的四个轮速中选取最小速度,并计算四个车轮的滑移率,作为模糊控制的输入量,输出量为车速估算值vE,w;
13)根据步骤12)得到的车速估算值vE,w,基于加速度的车速估算模块综合滤波后的加速度信号,采用小步长时间迭代积分,然后用线性补偿方式进行标定计算车速估算值VE a;
2)根据步骤1)得到的车速估算值vE,w和VE a在第二层算法中进行最终的纵向车速估算,调节在不同车轮滑移情况下的置信程度,使得在大滑移率、多轮打滑的情况下得出精确的纵向车速估算值;
21)首先通过车速估算模块根据相对准确的临时估算车速,考虑滑移率和滑移车轮个数,分别输出变量到滑移率置信度计算模块和最终模糊控制模块;
实施例:
本实施例中,根据上述方法实现的速度估算流程具体为:
首先,建立了基于轮速的车速估计器,通过模糊控制逻辑来估计车速,估计器包括将车轮角速度转换为车轮纵向速度,并且修正轮胎转动半径和轮胎标记半径的偏差。其次,设计了基于加速度的车速估计器,基于积分加速度变量来估计车速。第三,考虑了在四轮高滑移率的情况下,即使用模糊控制逻辑,在四轮高滑移率的情况下,也无法在第一次估算中获得准确的估计车速,因此在车速估计器的第二次模糊逻辑引入了相对准确的临时车速以更准确地判断车轮打滑状况,以获得准确的估计车速。第四,建立了车速估计器,通过置信度模块的加权系数,平衡从车轮编码器和加速度传感器获得的信号,然后得到最终估计车速。
所述的轮速编码器输出的信号为角速度信号,通过每个取样周期内的计数来计算车轮角速度,轮胎滚动半径计算公式为:
为了消除轮胎转动半径与轮胎标号半径的差异,需要引入车轮轮胎转动因数η,此系数由简单的车辆直线加速实验确定,车轮速度计算公式为:
Vi,j=3.6·ωij·Rij·η
其中Ri,j为轮胎半径,i代表前轮或后轮,j代表左轮或右轮;ωi,j为车轮角速度,i代表前轮或后轮,j代表左轮或右轮;η为轮胎转动系数;Vi,j为车轮速度,i代表前轮或后轮,j代表左轮或右轮。
所述的基于轮速的车速估算模块,选取四个轮速的最小值作为参考值,考虑到这个最小的轮速只作为车速的参考值,这个值并不是一个确定的轮胎给出的,而是由行驶过程中四个轮胎中的最小值确定的:
Vmin=min{Vi,j},1≤i,j≤2
分析看出实际分为一下几种情况:
1)所有轮胎都不打滑,那么在直线行驶的状态下,四个车轮的轮速应该是相等的,计算公式为:
vi,j=真实纵向车速
2)四个轮胎中至少有一个轮胎不打滑,那么假设根据这个车轮计算出的车速为Vm,m,计算公式为:
Vi,j≥Vm,m=真实纵向车速
其中m≠i,j,1≤m,i,j≤2;
3)四个轮胎都打滑,计算公式为:
Vi,j≥Vm,m≥真实纵向车速
其中m≠i,j,1≤m,i,j≤2。
这种情况很难准确计算出的车速准确地表征当前的纵向车速,因此需要比较接近的纵向车速数值作为参考,此数值与基于加速的的车速估算值比较后,根据置信度参数进行校准。
除四个车轮最小轮速外,滑移率也是最终模糊控制的输入变量,计算公式为:
其中m≠i,j,1≤m,i,j≤2。
所述的最终模糊控制模块,考虑了所有输入变量即最小轮速和四个轮胎的滑移率及滑移车轮的数量,判断出基于轮速的车速估算和基于加速度的车速估算模块计算出的车速的置信度,根据置信度赋予权重值,估算车速计算公式为:
VE_w=mon(Vi,j)/fuzzy logic weight
其中VE_w为基于轮速的估算车速。
所述的基于加速度的车速估算模块,采用小步长时间迭代积分,计算公式为:
其中asensor为车辆纵向加速度;Vtime为加速积分器的初始车速。
尽管经过噪声滤波并且积分过程由对应的处理来尽量使结果准确,但是由于传感器自身的精度问题,得到的估计车速仍然和实际车速有一定的偏差,采用线性补偿方式进行修正,计算公式为:
VE_a=V'E_a·标定系数
其中VE_a为基于加速度的估算车速。
所述的车速估算模块为了满足大滑移率且多轮打滑的情况,采用权重系数矩阵的方式平衡两种估算方法的纵向车速估算值,为最终的模糊算法提供精确的输入变量,车速估算公式为:
VE_temp=WVE_w+(1-W)VE_a
其中VE temp为相对精确的估算车速;W为fuzzy logicweight的矩阵。
所述的滑移率置信度计算模块考虑了滑移率和打滑车轮个数,如果轮胎打滑车轮数量小于2,置信参数偏向于基于轮速的估计车速;如果轮胎打滑数量大于2,置信参数偏向于基于加速度的估计车速,再进行逻辑估算最终车速vE。
Claims (7)
1.一种适用于四轮驱动车辆多轮高滑移率时的速度估计算法,其特征在于,包括轮速计算模块(1)、基于轮速的车速估算模块(2)、信号滤波模块(3)、基于加速度的车速估算模块(4)、滑移率置信度计算模块(5)、车速估算模块(6)和最终模糊控制模块(7),还包括以下步骤:
1)在第一层算法中进行预估算,在大滑移率的情况下分别基于四个轮速传感器和纵向加速度传感器估算出相对准确的临时车速,作为第二层算法的参考值;
2)根据步骤1)得到的车速估算值vE,w和VE a在第二层算法中进行最终的纵向车速估算,调节在不同车轮滑移情况下的置信程度,使得在大滑移率、多轮打滑的情况下得出精确的纵向车速估算值。
2.根据权利要求1所述的一种四轮驱动车辆多轮高滑移率时的速度估计算法,其特征在于,所述的步骤1)具体包括以下步骤:
11)轮速计算模块(1)根据传感器信号车轮角速度ωi,j,转化为车轮线速度Vi,j,引入轮胎转动系数η进行标定校正;
12)根据步骤11)得到的车轮速度Vi,j,基于轮速的车速估算模块(2)从传感器得到的四个轮速中选取最小速度,并计算四个车轮的滑移率,作为模糊控制的输入量,输出量为车速估算值vE,w;
13)根据步骤12)得到的车速估算值vE,w,基于加速度的车速估算模块(4)接收信号滤波模块(3)后的加速度信号,采用小步长时间迭代积分,然后用线性补偿方式进行标定计算车速估算值VE a。
4.根据权利要求2所述的一种四轮驱动车辆多轮高滑移率时的速度估计算法,其特征在于,所述的步骤12)中,基于轮速的车速估算模块的最小轮速的方式分别为:
第一种情况:所有轮胎都不打滑,那么在直线行驶的状态下,四个车轮的轮速应该是相等的,计算公式为:
Vi,j=真实纵向车速
第二种情况:四个轮胎中至少有一个轮胎不打滑,那么假设根据这个车轮计算出的车速为Vm,m,计算公式为:
Vi,j≥Vm,m=真实纵向车速
第三种情况:四个轮胎都打滑,计算公式为:
Vi,j≥Vm,m≥真实纵向车速
其中m≠i,j,1≤m,i,j≤2。
6.根据权利要求1所述的一种四轮驱动车辆多轮高滑移率时的速度估计算法,其特征在于,所述的步骤2)具体包括以下步骤:
21)首先通过车速估算模块(6)根据相对准确的临时估算车速,考虑滑移率和滑移车轮个数,分别输出变量到滑移率置信度计算模块(5)和最终模糊控制模块(7);
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010052297.1A CN111086520A (zh) | 2020-01-17 | 2020-01-17 | 一种适用于四轮驱动车辆多轮高滑移率时的速度估计算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010052297.1A CN111086520A (zh) | 2020-01-17 | 2020-01-17 | 一种适用于四轮驱动车辆多轮高滑移率时的速度估计算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111086520A true CN111086520A (zh) | 2020-05-01 |
Family
ID=70399288
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010052297.1A Pending CN111086520A (zh) | 2020-01-17 | 2020-01-17 | 一种适用于四轮驱动车辆多轮高滑移率时的速度估计算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111086520A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111619548A (zh) * | 2020-05-28 | 2020-09-04 | 清华大学 | 车辆驱动防滑控制方法、装置、计算机设备和存储介质 |
CN111775950A (zh) * | 2020-07-07 | 2020-10-16 | 清华大学苏州汽车研究院(吴江) | 车辆参考速度测算方法、装置、设备、存储介质和系统 |
CN112579966A (zh) * | 2021-03-01 | 2021-03-30 | 天津所托瑞安汽车科技有限公司 | Abs参考车速的计算方法、装置、电子设备和介质 |
CN113044041A (zh) * | 2021-03-09 | 2021-06-29 | 中国第一汽车股份有限公司 | 基于轮速及整车加速度的汽车车速及轮胎滑移率估计方法 |
CN113715836A (zh) * | 2021-11-04 | 2021-11-30 | 上海拿森汽车电子有限公司 | 车速估算方法、系统及计算机存储介质 |
CN114684157A (zh) * | 2022-04-15 | 2022-07-01 | 北京新能源汽车股份有限公司 | 一种车速估计方法、装置和电动汽车 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102009654A (zh) * | 2010-11-12 | 2011-04-13 | 清华大学 | 一种全轮电驱动车辆的纵向车速估计方法 |
CN104021310A (zh) * | 2014-06-23 | 2014-09-03 | 清华大学 | 基于UKF与修正Dugoff轮胎模型路面峰值附着系数估算方法 |
CN106394561A (zh) * | 2015-11-10 | 2017-02-15 | 北京中科易电信息科技股份有限公司 | 一种车辆的纵向车速的估计方法和装置 |
CN110562263A (zh) * | 2019-09-19 | 2019-12-13 | 中国人民解放军陆军装甲兵学院 | 一种基于多模型融合的轮毂电机驱动车辆车速估计方法 |
-
2020
- 2020-01-17 CN CN202010052297.1A patent/CN111086520A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102009654A (zh) * | 2010-11-12 | 2011-04-13 | 清华大学 | 一种全轮电驱动车辆的纵向车速估计方法 |
CN104021310A (zh) * | 2014-06-23 | 2014-09-03 | 清华大学 | 基于UKF与修正Dugoff轮胎模型路面峰值附着系数估算方法 |
CN106394561A (zh) * | 2015-11-10 | 2017-02-15 | 北京中科易电信息科技股份有限公司 | 一种车辆的纵向车速的估计方法和装置 |
CN110562263A (zh) * | 2019-09-19 | 2019-12-13 | 中国人民解放军陆军装甲兵学院 | 一种基于多模型融合的轮毂电机驱动车辆车速估计方法 |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111619548A (zh) * | 2020-05-28 | 2020-09-04 | 清华大学 | 车辆驱动防滑控制方法、装置、计算机设备和存储介质 |
CN111619548B (zh) * | 2020-05-28 | 2021-01-12 | 清华大学 | 车辆驱动防滑控制方法、装置、计算机设备和存储介质 |
CN111775950A (zh) * | 2020-07-07 | 2020-10-16 | 清华大学苏州汽车研究院(吴江) | 车辆参考速度测算方法、装置、设备、存储介质和系统 |
CN111775950B (zh) * | 2020-07-07 | 2022-02-11 | 清华大学苏州汽车研究院(吴江) | 车辆参考速度测算方法、装置、设备、存储介质和系统 |
CN112579966A (zh) * | 2021-03-01 | 2021-03-30 | 天津所托瑞安汽车科技有限公司 | Abs参考车速的计算方法、装置、电子设备和介质 |
CN113044041A (zh) * | 2021-03-09 | 2021-06-29 | 中国第一汽车股份有限公司 | 基于轮速及整车加速度的汽车车速及轮胎滑移率估计方法 |
CN113715836A (zh) * | 2021-11-04 | 2021-11-30 | 上海拿森汽车电子有限公司 | 车速估算方法、系统及计算机存储介质 |
CN114684157A (zh) * | 2022-04-15 | 2022-07-01 | 北京新能源汽车股份有限公司 | 一种车速估计方法、装置和电动汽车 |
CN114684157B (zh) * | 2022-04-15 | 2024-09-17 | 北京新能源汽车股份有限公司 | 一种车速估计方法、装置和电动汽车 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111086520A (zh) | 一种适用于四轮驱动车辆多轮高滑移率时的速度估计算法 | |
CN110562263B (zh) | 一种基于多模型融合的轮毂电机驱动车辆车速估计方法 | |
CN110987470B (zh) | 一种基于模型迭代的汽车质量在线估计方法 | |
CN110095635B (zh) | 一种全轮驱动车辆的纵向车速估计方法 | |
CN113002549B (zh) | 一种车辆状态估算方法、装置、设备及存储介质 | |
US6789017B2 (en) | Vehicle steering angle position determination method | |
US20030195689A1 (en) | Method of estimating quantities that represent state of vehicle | |
CN109606378A (zh) | 面向非高斯噪声环境的车辆行驶状态估计方法 | |
CN108819950B (zh) | 汽车稳定性控制系统的车速估计方法及系统 | |
CN109515445A (zh) | 一种用于全轮独立驱动车辆的纵向车速估计方法及其装置 | |
CN110979026A (zh) | 一种基于实时路况的分布式驱动公交车转矩分配方法 | |
CN113276862B (zh) | 车辆行驶状态估计方法 | |
CN114572224B (zh) | 一种路面最大附着系数的估算方法及终端 | |
CN113771857A (zh) | 一种用于车辆控制的纵向车速估计方法和系统 | |
CN113044041A (zh) | 基于轮速及整车加速度的汽车车速及轮胎滑移率估计方法 | |
CN113715836A (zh) | 车速估算方法、系统及计算机存储介质 | |
CN109282832A (zh) | 适用于特殊路面的惯性辅助里程计自适应标定方法及系统 | |
CN111942399A (zh) | 一种基于无迹卡尔曼滤波的车速估算方法及系统 | |
CN111806430A (zh) | 一种用于自动泊车的车速计算方法 | |
CN111062088B (zh) | 一种四轮独立驱动电动汽车参考车速容错自适应估计方法 | |
CN111231976A (zh) | 一种基于变步长的车辆状态估计方法 | |
CN113978476B (zh) | 一种考虑传感器数据丢失的线控汽车轮胎侧向力估计方法 | |
CN114043986B (zh) | 一种考虑质量失配的轮胎路面附着系数多模型融合估计方法 | |
CN113911130B (zh) | 一种车辆横摆角速度估计方法及系统 | |
CN114590264A (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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20200501 |
|
WD01 | Invention patent application deemed withdrawn after publication |