CN108534753A - 一种基于车轮垂直动载荷的路面功率谱测量方法 - Google Patents

一种基于车轮垂直动载荷的路面功率谱测量方法 Download PDF

Info

Publication number
CN108534753A
CN108534753A CN201810286506.1A CN201810286506A CN108534753A CN 108534753 A CN108534753 A CN 108534753A CN 201810286506 A CN201810286506 A CN 201810286506A CN 108534753 A CN108534753 A CN 108534753A
Authority
CN
China
Prior art keywords
vertical dynamic
wheel
power spectrum
dynamic load
vehicle
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
Application number
CN201810286506.1A
Other languages
English (en)
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.)
Yanshan University
Original Assignee
Yanshan University
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 Yanshan University filed Critical Yanshan University
Priority to CN201810286506.1A priority Critical patent/CN108534753A/zh
Publication of CN108534753A publication Critical patent/CN108534753A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C7/00Tracing profiles
    • G01C7/02Tracing profiles of land surfaces
    • G01C7/04Tracing profiles of land surfaces involving a vehicle which moves along the profile to be traced

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)
  • Vehicle Body Suspensions (AREA)

Abstract

本发明公开了一种基于车轮垂直动载荷的路面功率谱测量方法,该方法是基于所设计装有六分力仪的实验车实现,其内容包括:首先利用车轮六分力仪采集实验车轮胎垂直动载荷信号,然后建立路面不平度与车轮垂直动载荷间关系,进一步地通过校准实验车的放大因子得到路面功率谱,最后根据垂直动载荷功率谱进行路面不平度特征识别。本发明基于车轮垂直动载荷的路面功率谱测量方法减小了道路表面清洁度对路面不平度真实值的影响,同时具有良好的测量动态稳定性;设计的实验车避免了由于发动机振动及传动系统回转引起的垂直动载荷,提高了车轮垂直动载荷谱测量的准确性;提出的路面不平度特征识别方法,简便直观,便于工作人员对路面进行评价。

Description

一种基于车轮垂直动载荷的路面功率谱测量方法
技术领域
本发明属于路面特征识别领域,具体涉及一种基于车轮垂直动载荷的路面功率谱测量方法。
背景技术
路面功率谱密度函数常用来表示路面不平度,并作为车辆动态响应、悬架的优化和控制等计算的输入单位,其描述精度和准确性直接影响后续研究。
目前,针对路面功率谱的采集与测量多采用加速度传感器和激光传感器的间接式测量方法,导致测量误差较大、设备价格昂贵以及适用范围小等问题。其中常用的基于惯性基准的激光断面仪具有很多缺点:(1)在路面检测中,道路的表面清洁度会影响到数据的真实值;(2)静态路面检测精度较高,但未考虑由车体振动而引起加速度传感器产生的检测偏差,故其难以检测复杂路面;(3)在雨天或路面有水的情况下,不能进行路面检测,故其受天气影响加大。这些缺点大幅影响路面功率谱的描述精度和准确性,也极大地限制了应用各类传感器的间接式测量方法。
应用车轮六分力仪的直接式测量路面功率谱方法,可根据车轮垂直动载荷的大小和变化快慢反应路面不平度状况,其测量误差较小,具有良好的动态特性,可满足复杂路面功率谱测量要求且不受路面清洁度及天气的影响。
发明内容
为了解决现有路面功率谱测量技术中存在的测量误差较大、动态稳定性差以及适用范围小等问题,本发明提供一种基于车轮垂直动载荷的路面功率谱测量方法。当车辆在不平路面匀速行驶时,路面不平度将直接影响车轮所受垂直动载荷。首先利用车轮六分力仪采集实验车轮胎垂直动载荷数据,然后建立路面不平度与车轮垂直动载荷间关系,进一步地通过校准实验车的放大因子得到路面功率谱,最后根据垂直动载荷功率谱进行路面不平度特征识别。
本发明为达到上述目的,采取如下技术方案:
一种基于车轮垂直动载荷的路面功率谱测量方法,该方法是基于所设计装有六分力仪的实验车实现的,该方法内容包括以下步骤:
步骤1:测量实验车轮胎垂直动载荷信号,并进行信号预处理,得到垂直动载荷功率谱;
步骤2:基于四分之一车辆模型建立路面不平度与车轮垂直动载荷间关系表达式;
步骤3:校准实验车的放大因子,并计算路面不平度功率谱;
步骤4:根据垂直动载荷功率谱进行路面不平度特征识别。
进一步的,在步骤1中,所述实验车包括车轮和该车轮对应的车身,车轮为非簧载质量,车身为簧载质量;该车轮指实验车中的一个车轮;车身是指车轮所对应的实验车的相应车身部分;在车身与车轮之间安装有减振器,其阻尼为cs;车身与车轮之间的连接等效于悬架弹簧,其刚度为ks;牵引车通过拖车钩拉动实验车行进;其中,车轮中装有车轮六分力仪;
所述测量实验车轮胎垂直动载荷信号,并进行信号预处理,得到垂直动载荷功率谱;其具体实现过程如下:
所述牵引车拉动实验车在待测路面上行驶,测量实验车轮胎振动的垂直动载荷信号;
对测量的垂直动载荷信号进行预处理,去除信号中的噪声因子及趋势项得到预处理垂直动载荷信号F(t);
根据所得预处理垂直动载荷信号F(t),可以得到垂直动载荷功率谱密度函数,即下式:
其中,GF(w)表示垂直动载荷功率谱密度函数,T表示预处理垂直动载荷信号F(t)的时间区间。
进一步的,在步骤2中,所述基于四分之一车辆模型建立路面不平度与车轮垂直动载荷间关系表达式,是基于四分之一车辆模型实现的,就是对于四分之一车辆振动系统,由路面不平度函数推导车轮垂直动载荷表达式,所得表达式表明车轮垂直动载荷与路面不平度激励具有相同的频率,只是幅值不同;
设车轮垂直动载荷与路面不平度幅值间存在放大因子P,
其中,kt为轮胎刚度,A和ω分别为路面不平度函数的振幅和频率,|Zw|和分别为复振幅Zw的模和相位角;
(2)式表明放大因子P与车辆结构有关。
进一步的,在步骤3中,所述校准实验车的放大因子就是经多次实验,通过数学统计理论去除测量过程中的不准确数据,得到校准后的实验车的放大因子P;
牵引车拉动实验车在已知实验路面上行驶,测量实验车轮胎垂直动载荷信号,测量进行多次,得到一组垂直动载荷信号Fi(t);
已知实验路面的随机不平度信号y(t);
根据车辆放大因子定义,即设车轮垂直动载荷与路面不平度幅值间存在放大因子P,可得到实验车多次测量的一组放大因子:
通过数学统计理论去除测量的不准确数据,得到校准后的实验车放大因子P。
进一步的,在步骤4中,所述根据垂直动载荷功率谱进行路面不平度特征识别,就是依据垂直动载荷功率谱构成的不同,定义路面识别特征值ξ,定量的进行路面不平度特征识别;
不同路面不平度引起的垂直动载荷的差异主要体现在功率谱结构的区别,故可依据垂直动载荷功率谱构成的不同,作为对路面不平度识别的标准;
定义垂直动载荷功率谱在低频部分有特征量H1,在高频部分有特征量H2,则有:
其中,GF(ω)为垂直动载荷功率谱密度函数;
根据垂直动载荷功率谱在低频部分和高频部分的构成,定义路面识别特征值ξ:
根据车辆在不同路面行驶的识别特征值ξ对路面不平度做出特征识别。
本发明与现有技术相比具有如下有益效果:
(1)本发明通过车轮垂直动载荷直接测量的路面功率谱,具有测量误差较小、动态特性良好以及适用广泛等优点。垂直动载荷的大小和变化的快慢能够准确地反应路面的不平度状况,减小了道路表面清洁度对路面不平度真实值的影响,同时具有良好的测量动态稳定性。
(2)本发明所设计的实验车仅由车轮和简易减震装置组成,避免了由于发动机振动及传动系统回转所引起的垂直动载荷的影响,减少了产生垂直动载荷的激励源,提高了垂直载荷谱测量的准确性。并且所设计的实验车结构简单,自由度少,可近似于线性结构,便于车辆结构参数计算,提高了路面功率谱测量的准确性。实验车可由不同牵引车拉动,适用于各种复杂路面测量且不受特殊天气情况限制,对提高复杂路面检测准确性,促进高质量路面建设具有重要意义。
(3)本发明所提出的根据垂直动载荷功率谱构成不同进行路面不平度特征识别的方法,直接从垂直动载荷功率谱中寻找能反映路面不平度特征的信息,以此来实现对路面不平度的定量识别。这种特征识别方法具有简便直观的优点,便于工作人员对路面进行评价。
附图说明
图1为本发明方法的流程示意图;
图2为测量小车的整体示意图;
图3为四分之一车辆被动悬架系统的示意图;
图4为车轮垂直动载荷功率谱的示意图。
具体实施方式
下面结合附图与具体实施方式对本发明作进一步详细描述:
当车辆在不平路面匀速行驶时,路面不平度状况将直接表现为车轮动载荷的大小及变化快慢。首先利用车轮六分力仪采集实验车轮胎垂直动载荷数据,然后建立路面不平度与车轮垂直动载荷间关系,进一步地通过校准实验车的放大因子得到路面功率谱,最后根据垂直动载荷功率谱构成不同进行路面不平度特征识别。据此提出本发明的一种基于车轮垂直动载荷的路面功率谱测量方法,如图1所示,该方法步骤内容详细说明如下:
步骤1:测量实验车轮胎垂直动载荷信号,并进行信号预处理,得到垂直动载荷功率谱;
所述实验车整体示意图如图2所示,包括车轮1和该车轮1对应的车身2,车轮1为非簧载质量,车身2为簧载质量。该车轮1指实验车中的一个车轮。车身2是指车轮1所对应的实验车的相应车身部分。在车身2与车轮1之间安装有减振器3,其阻尼为cs。车身2与车轮1之间的连接等效于悬架弹簧4,其刚度为ks。牵引车通过拖车钩5拉动实验车行进。其中,车轮1中装有车轮六分力仪。
牵引车拉动实验车在待测路面上行驶,测量实验车轮胎振动的垂直动载荷信号。
对测量的垂直动载荷信号进行预处理,去除信号中的噪声因子及趋势项得到预处理垂直动载荷信号F(t)。
根据所得预处理垂直动载荷信号F(t),可以得到垂直动载荷功率谱密度函数,即下式:
其中,GF(w)表示垂直动载荷功率谱密度函数,T表示预处理垂直动载荷信号F(t)的时间区间。
步骤2:基于四分之一车辆模型建立路面不平度与车轮垂直动载荷间关系表达式;
所述基于四分之一车辆模型建立路面不平度与车轮垂直动载荷间关系表达式,是基于四分之一车辆模型实现的,就是对于四分之一车辆振动系统,由路面不平度函数推导车轮垂直动载荷表达式,所得表达式表明车轮垂直动载荷与路面不平度激励具有相同的频率,只是幅值不同。
建立四分之一车辆被动悬架系统动力学方程:
其中,mb为簧载质量,mw为非簧载质量,cs为减振器阻尼系数,ks为弹簧刚度,kt为轮胎刚度,zb为车身垂直位移,zw为轮胎垂直位移,y为路面的不平度函数。
车轮垂直力:
Fz=kt(y-zw) (4)
假设路面不平度为单位简谐函数:
y=Asin(ωt) (5)
将(4)式代入悬架系统动力学方程,并改写为矩阵形式:
利用复指数法求解,令Aktsin(ωt)=Aktejωt,设方程的解为:
zb(t)=Z1ejωt,zw(t)=Z2ejωt (7)
将(7)式代入矩阵方程(6)中,化简得:
求得:
式中|Zw|和分别为复振幅Zw的模和相位角。取实部,路面不平度激励下的轮胎垂直位移Zw
代入(4)式得车轮垂直动载荷:
上式表明车轮垂直动载荷与路面不平度激励具有相同的频率,只是幅值不同。
设车轮垂直动载荷与路面不平度幅值间存在放大因子P,则有:
上式表明放大因子P与车辆结构有关。
步骤3:校准实验车的放大因子,并计算路面不平度功率谱;
所述校准实验车的放大因子就是经多次实验,通过数学统计理论去除测量过程中的不准确数据,得到校准后的实验车放大因子P。
牵引车拉动实验车在已知实验路面上行驶,测量实验车轮胎垂直动载荷信号,测量进行多次,得到一组垂直动载荷信号Fi(t)。
已知实验路面的随机不平度信号y(t)。
根据所述车辆放大因子定义,可得到实验车多次测量的放大因子
通过数学统计理论去除测量的不准确数据,得到校准后的实验车放大因子P。根据步骤1中预处理垂直动载荷信号F(t)和校准后的实验车放大因子P,得到路面不平度信号Y(t),公式如下:
进一步地,可以得到悬架行程功率谱密度函数,即下式:
步骤4:根据垂直动载荷功率谱进行路面不平度特征识别。
所述根据垂直动载荷功率谱进行路面不平度特征识别就是依据垂直动载荷功率谱构成的不同,定义路面识别特征值ξ,定量的进行路面不平度特征识别。
如图4所示为车轮垂直动载荷功率谱的示意图。不同路面不平度引起的垂直动载荷的差异主要体现在功率谱结构的区别,故可依据垂直动载荷功率谱构成的不同,作为对路面不平度识别的标准。
由路面不平度引起的垂直动载荷主要集中在低频,少数部分作用于高频。
在频率区间[0,β]内,将垂直动载荷功率谱分为低频部分和高频部分,通过实验数据定义以频率点α为分段节点。如在图4所示车轮垂直动载荷功率谱中,以100Hz为分段节点,左侧为低频部分,右侧为高频部分。
定义垂直动载荷功率谱在低频部分有特征量H1,在高频部分有特征量H2
其中,GF(ω)为垂直动载荷功率谱密度函数。
根据垂直动载荷功率谱在低频部分和高频部分的构成,定义路面识别特征值ξ
根据车辆在不同路面行驶的识别特征值ξ对路面不平度做出特征识别。
本发明利用车轮垂直动载荷来进行路面功率谱测量,具有测量误差较小,适用广泛,减小了道路表面清洁度对路面不平度真实值的影响,良好的测量动态稳定性等优点;同时装有车轮六分力仪的实验车结构简单,避免了由于发动机振动及传动系统回转引起的垂直动载荷,提高了车轮垂直动载荷谱测量的准确性;另外所提出的路面不平度特征识别方法,简便直观,便于工作人员对路面进行评价。

Claims (4)

1.一种基于车轮垂直动载荷的路面功率谱测量方法,其特征在于:该方法是基于所设计装有六分力仪的实验车实现的,该方法内容包括以下步骤:
步骤1:测量实验车轮胎垂直动载荷信号,并进行信号预处理,得到垂直动载荷功率谱;
步骤2:基于四分之一车辆模型建立路面不平度与车轮垂直动载荷间关系表达式;
步骤3:校准实验车的放大因子,并计算路面不平度功率谱;
步骤4:根据垂直动载荷功率谱进行路面不平度特征识别。
2.根据权利要求1所述的一种基于车轮垂直动载荷的路面功率谱测量方法,其特征在于:在步骤1中,所述实验车包括车轮和该车轮对应的车身,车轮为非簧载质量,车身为簧载质量;该车轮指实验车中的一个车轮;车身是指车轮所对应的实验车的相应车身部分;在车身与车轮之间安装有减振器,其阻尼为cs;车身与车轮之间的连接等效于悬架弹簧,其刚度为ks;牵引车通过拖车钩拉动实验车行进;其中,车轮中装有车轮六分力仪;
所述测量实验车轮胎垂直动载荷信号,并进行信号预处理,得到垂直动载荷功率谱;其具体实现过程如下:
所述牵引车拉动实验车在待测路面上行驶,测量实验车轮胎振动的垂直动载荷信号;
对测量的垂直动载荷信号进行预处理,去除信号中的噪声因子及趋势项得到预处理垂直动载荷信号F(t);
根据所得预处理垂直动载荷信号F(t),可以得到垂直动载荷功率谱密度函数,即下式:
其中,GF(w)表示垂直动载荷功率谱密度函数,T表示预处理垂直动载荷信号F(t)的时间区间。
3.根据权利要求1所述的一种基于车轮垂直动载荷的路面功率谱测量方法,其特征在于:在步骤3中,所述校准实验车的放大因子就是经多次实验,通过数学统计理论去除测量过程中的不准确数据,得到校准后的实验车的放大因子P;
牵引车拉动实验车在已知实验路面上行驶,测量实验车轮胎垂直动载荷信号,测量进行多次,得到一组垂直动载荷信号Fi(t);
已知实验路面的随机不平度信号y(t);
根据车辆放大因子定义,即设车轮垂直动载荷与路面不平度幅值间存在放大因子P,可得到实验车多次测量的一组放大因子:
通过数学统计理论去除测量的不准确数据,得到校准后的实验车放大因子P。
4.根据权利要求1所述的一种基于车轮垂直动载荷的路面功率谱测量方法,其特征在于:在步骤4中,所述根据垂直动载荷功率谱进行路面不平度特征识别,就是依据垂直动载荷功率谱构成的不同,定义路面识别特征值ξ,定量的进行路面不平度特征识别;
不同路面不平度引起的垂直动载荷的差异主要体现在功率谱结构的区别,故可依据垂直动载荷功率谱构成的不同,作为对路面不平度识别的标准;
定义垂直动载荷功率谱在低频部分有特征量H1,在高频部分有特征量H2,则有:
其中,GF(ω)为垂直动载荷功率谱密度函数;
根据垂直动载荷功率谱在低频部分和高频部分的构成,定义路面识别特征值ξ:
根据车辆在不同路面行驶的识别特征值ξ对路面不平度做出特征识别。
CN201810286506.1A 2018-03-30 2018-03-30 一种基于车轮垂直动载荷的路面功率谱测量方法 Pending CN108534753A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810286506.1A CN108534753A (zh) 2018-03-30 2018-03-30 一种基于车轮垂直动载荷的路面功率谱测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810286506.1A CN108534753A (zh) 2018-03-30 2018-03-30 一种基于车轮垂直动载荷的路面功率谱测量方法

Publications (1)

Publication Number Publication Date
CN108534753A true CN108534753A (zh) 2018-09-14

Family

ID=63482897

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810286506.1A Pending CN108534753A (zh) 2018-03-30 2018-03-30 一种基于车轮垂直动载荷的路面功率谱测量方法

Country Status (1)

Country Link
CN (1) CN108534753A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108303066A (zh) * 2018-01-04 2018-07-20 燕山大学 一种基于悬架行程测量路面功率谱的方法
CN112696436A (zh) * 2020-12-28 2021-04-23 吉林大学 一种基于路面载荷谱的轮毂电机轴承热管理系统及其控制方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5065618A (en) * 1990-11-14 1991-11-19 Hodges Transportation Inc. Method and apparatus for determining terrain surface profiles
CN103981795A (zh) * 2014-05-28 2014-08-13 江苏科技大学 一种利用车辆悬架传感器实现路面谱软测量的方法
CN104233935A (zh) * 2014-08-28 2014-12-24 吉林大学 一种基于道路纵断面信息的路面质量等级辨识方法
CN104864846A (zh) * 2015-03-24 2015-08-26 江苏科技大学 基于加速度传感器的车载式路面不平度采集系统及工作方法
CN105427375A (zh) * 2015-12-16 2016-03-23 北京汽车研究总院有限公司 一种三维路面轮廓生成的方法及系统
CN105973189A (zh) * 2016-05-05 2016-09-28 广西科技大学 软路面路面不平度加载谱的获取方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5065618A (en) * 1990-11-14 1991-11-19 Hodges Transportation Inc. Method and apparatus for determining terrain surface profiles
CN103981795A (zh) * 2014-05-28 2014-08-13 江苏科技大学 一种利用车辆悬架传感器实现路面谱软测量的方法
CN104233935A (zh) * 2014-08-28 2014-12-24 吉林大学 一种基于道路纵断面信息的路面质量等级辨识方法
CN104864846A (zh) * 2015-03-24 2015-08-26 江苏科技大学 基于加速度传感器的车载式路面不平度采集系统及工作方法
CN105427375A (zh) * 2015-12-16 2016-03-23 北京汽车研究总院有限公司 一种三维路面轮廓生成的方法及系统
CN105973189A (zh) * 2016-05-05 2016-09-28 广西科技大学 软路面路面不平度加载谱的获取方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108303066A (zh) * 2018-01-04 2018-07-20 燕山大学 一种基于悬架行程测量路面功率谱的方法
CN112696436A (zh) * 2020-12-28 2021-04-23 吉林大学 一种基于路面载荷谱的轮毂电机轴承热管理系统及其控制方法

Similar Documents

Publication Publication Date Title
CN104120644B (zh) 一种基于重力加速度传感器的路面平整度检测方法
JP6132304B2 (ja) 路面評価装置及び方法
CN101619968B (zh) 一种路面平整度检测方法和装置
CN103981795B (zh) 一种利用车辆悬架传感器实现路面谱软测量的方法
CN108839657B (zh) 一种基于汽车振动响应在线识别路面不平度信息的方法
CN110184885B (zh) 一种基于智能手机测试路面平整度的方法
CN106644519B (zh) 一种车辆垂向动力学模型参数辨识的方法及装置
CN109186528B (zh) 一种基于手机测得的加速度确定国际平整度指数的理论方法
CN102700551A (zh) 一种车辆行驶过程中路面坡度的实时估计方法
CN105000033B (zh) 一种轨道几何不平顺的检测评估系统
CN103625475B (zh) 一种基于递归最小二乘的车辆侧倾角与俯仰角估计方法
CN106840519A (zh) 一种两轴车辆质心高度的静态测量方法
CN108556850B (zh) 一种多轴重型车辆车轮垂向载荷通用计算方法
CN108534753A (zh) 一种基于车轮垂直动载荷的路面功率谱测量方法
CN105946866A (zh) 一种汽车行驶状态预警方法
D'Hooge et al. Application of real-world wind conditions for assessing aerodynamic drag for on-road range prediction
CN201530980U (zh) 一种路面平整度检测装置
CN104005324B (zh) 一种路面构造信息的检测系统
CN111006884B (zh) 基于傅立叶变换的车轮车轴侧偏角与侧偏刚度的测量方法
CN113389115B (zh) 车辆特性和路面平整度检测方法、装置、设备及存储介质
CN109849932B (zh) 一种路面自适应的车轮动载荷估计方法
CN108228994A (zh) 在越野道路随机激励下车辆及工作装置的应力计算方法
CN107144384A (zh) 一种轮胎经过错台时产生冲击力的测试计算方法
CN112526880B (zh) 一种车辆行驶过程中的路面高度实时估计方法
CN104949649A (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
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20180914