CN111141520A - 一种基于改进经验小波变换的滚动轴承故障诊断方法 - Google Patents
一种基于改进经验小波变换的滚动轴承故障诊断方法 Download PDFInfo
- Publication number
- CN111141520A CN111141520A CN202010113154.7A CN202010113154A CN111141520A CN 111141520 A CN111141520 A CN 111141520A CN 202010113154 A CN202010113154 A CN 202010113154A CN 111141520 A CN111141520 A CN 111141520A
- Authority
- CN
- China
- Prior art keywords
- sub
- frequency
- empirical
- band
- rolling bearing
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/04—Bearings
- G01M13/045—Acoustic or vibration analysis
Abstract
本发明公开了一种基于改进经验小波变换的滚动轴承故障诊断方法,涉及故障诊断、信号处理技术领域,该方法包括:采集样本振动信号,对振动信号进行傅里叶变换从时域信号变换到频域信号,得到频域信号的主频带;然后利用尺度空间法对主频带进行自适应划分,划分为一系列子频带;计算主频带和各子频带的均值、方差,将均值和方差均小于主频带的子频带定义为冗余分量,将冗余分量与其下一个子频带合并;然后,在优化后的每个频带上建立自适应经验小波滤波器,获得对应的经验模式;最后,筛选出裕度因子最大的经验模式进行包络谱分析;该方法能够有效从强背景噪声下提取到滚动轴承的故障特征频率,很大程度上提高了滚动轴承故障诊断的准确度。
Description
技术领域
本发明涉及故障诊断、信号处理技术领域,尤其是一种基于改进经验小波变换的滚动轴承故障诊断方法。
背景技术
在长期重压和非平稳运行条件下,滚动轴承的早期微弱故障往往难以得到有效监控。早期微弱故障的产生和不能得到有效监控和诊断会使旋转机械不能有效运行甚至会导致机械部件的损坏。振动分析法被认为是滚动轴承故障检测和诊断的有效方法之一,通过对采集的振动信号进行分析和处理来判断滚动轴承是否存在故障。通常,故障特征频率是滚动轴承出现故障的关键性指标,但滚动轴承工作环境往往存在大量背景噪声和干扰信号,采用常规时域或频域分析方法往往会难以提取到滚动轴承的故障特征。
经验小波变换(Experience Wavelet Transform,EWT)被认为是滚动轴承故障特征提取的有效时频分析方法之一。但应用经典频谱分割方法来分割频带时,容易过分割导致冗余分量的产生不利于后续完成滚动轴承故障特征频率的提取,难以从强背景噪声和干扰信号下提取到滚动轴承的故障特征频率。
发明内容
本发明人针对上述问题及技术需求,提出了一种基于改进经验小波变换的滚动轴承故障诊断方法,该方法能自适应的对频谱进行划分,将被定义为冗余分量的子频带进行合并,在优化后的子频带上建立自适应小波滤波器得到对应的经验模式,筛选出裕度因子最大的经验模式进行包络谱分析,所得结果能很好的从强背景噪声和干扰信号下提取到滚动轴承的故障特征频率。
本发明的技术方案如下:
一种基于改进经验小波变换的滚动轴承故障诊断方法,包括如下步骤:
获取滚动轴承故障样本的振动信号,对振动信号进行FFT变换,得到频域信号的主频带;利用尺度空间法对主频带的频谱进行自适应划分得到一系列子频带;将一系列子频带的频域特征参数与主频带的频域特征参数进行比较,确定一系列子频带中的冗余分量,将每一个冗余分量与下一个子频带进行合并,以优化尺度空间法划分的子频带的频谱;在优化后的子频带上建立经验小波滤波器组,获得一系列经验模式;计算每个经验模式的裕度因子,筛选出裕度因子最大的经验模式进行包络谱分析;从包络谱中提取出滚动轴承的故障特征频率进行分析,从而完成滚动轴承故障样本的故障诊断。
其进一步的技术方案为,利用尺度空间法对主频带的频谱进行自适应划分得到一系列子频带,包括,
尺度空间法的定义如下:
离散函数f(m)的尺度空间L(m,t)定义为:
首先在离散化频谱中寻找所有的局部极小值点,每一初始极小值点都会对应一条曲线;在尺度空间L(ω,t)中,随着尺度参数t的增长,局部极小值点的数量减少,随着移动步长的变化,每个局部极小值点对应的曲线进行叠加生成一条直线,直线用Ci(i∈[1,N0])表示,直线长度为Li,采用Otsu方法自适应确定一个阈值T,若直线长度Li大于阈值T,则将对应的局部极小值点ωi作为频带划分点。
其进一步的技术方案为,将一系列子频带的频域特征参数与主频带的频域特征参数进行比较,确定一系列子频带中的冗余分量,包括,
频域特征参数集包括主频带的第一均值和第一方差、以及每个子频带的第二均值和第二方差,将第二均值和第二方差均小于第一均值和第一方差的子频带定义为冗余分量;
第二均值V和第二方差M的计算公式分别如下:
其中,X(k)为频域幅值谱,X(k)的定义如下:
其进一步的技术方案为,在优化后的子频带上建立经验小波滤波器组,获得一系列经验模式,包括,依照经验小波的定义,在优化后的每个子频带上构造带通经验小波滤波器,形成滤波器组,由Meyer小波理论来构造经验小波,经验小波的尺度函数φn(ω)和小波函数ψn(ω)定义为:
其中,x(ω)、ψn(ω)和φ1(ω)分别代表振动信号的离散化频谱、离散化频谱的小波函数和经验小波的尺度函数,和分别表示x(ω)、ψn(ω)和φ1(ω)的傅里叶变换,和则分别代表各自的共轭,F-1(·)表示傅里叶逆变换;
振动信号x(t)经过经验小波变换可以得到如下经验模式:
基于上述细节系数,振动信号x(t)被重构为:
其进一步的技术方案为,计算每个经验模式的裕度因子,包括,裕度因子CL为对冲击脉冲信号较为敏感的参数,裕度因子CL的定义如下:
其中,XP=max{|xi|}表示峰值,故障振动信号用xi表示。
本发明的有益技术效果是:
根据本申请公开的一种基于改进经验小波变换的滚动轴承故障诊断方法,在尺度空间法划分完各个子频带后,通过计算主频带和各子频带的均值、方差,将均值和方差均小于主频带的子频带定义为冗余分量,将冗余分量和下一个子频带进行合并,以优化尺度空间法划分的子频带的频谱,通过减少冗余分量的产生以避免发生主频带过分割的现象,优化后的子频带通过经验小波滤波器组进行经验小波变换后产生一系列经验模式,通过计算每个经验模式的裕度因子,筛选出裕度因子最大的经验模式并进行包络谱分析,使用本申请的方法能够有效从强背景噪声和干扰信号中提取到滚动轴承的故障特征频率,很大程度上提高了滚动轴承故障诊断的准确度,对于滚动轴承的故障诊断技术有着重要的意义,并且能够推动故障诊断技术进一步发展。
附图说明
图1是本申请提供的基于改进经验小波变换的滚动轴承故障诊断方法的流程图。
具体实施方式
下面结合附图对本发明的具体实施方式做进一步说明。
请参考图1,本申请公开了一种基于改进经验小波变换的滚动轴承故障诊断方法,该方法针对故障轴承旋转工作下提取的振动信号进行分析,即针对强背景噪声和干扰信号下滚动轴承故障特征提取和降噪处理,该方法包括:
S1:获取滚动轴承故障样本的振动信号x(t),对振动信号x(t)进行FFT变换,使其由时域信号x(t)变换到频域信号X(f),并且将频域信号X(f)的频率范围归一化到[0,π]范围内,得到振动频域信号的主频带。
其中,定义频域信号X(f)为:
S2:利用尺度空间法对主频带的频谱进行自适应划分得到一系列子频带。
尺度空间法的相关步骤如下:
S21:定义一个离散函数f(m),离散函数f(m)的尺度空间L(m,t)定义为:
式(2)中,代表高斯核函数,t为尺度参数,m为离散函数的自变量,M为离散函数移动步长最大值,且其中3≤C≤6,n为离散函数的移动步长。其中,参数C影响着计算误差,当C=6时,可以使误差最小化,达到最大精度。
S23:确定自适应寻找频带划分点的方法如下:
首先在离散化频谱X(ω)中寻找所有的局部极小值点,每一初始极小值点都会对应一条曲线。在尺度空间L(ω,t)中,随着尺度参数t的增长,局部极小值点的数量减少。随着尺度步长的变化,每个局部极小值点对应的曲线进行叠加生成一条直线,此直线用Ci(i∈[1,N0])表示,直线长度为Li。采用Otsu方法自适应确定一个阈值T,若直线长度Li大于阈值T,则将对应的局部极小值点ωi作为频带划分点。
S3:将一系列子频带的频域特征参数与主频带的频域特征参数进行比较,确定一系列子频带中的冗余分量,将每一个冗余分量与下一个子频带进行合并,以优化尺度空间法划分的子频带的频谱。
频域特征参数集包括主频带的第一均值和第一方差、以及每个子频带的第二均值和第二方差,将第二均值和第二方差均小于第一均值和第一方差的子频带定义为所述冗余分量。
第二均值V和第二方差M的计算公式分别如下:
其中,X(k)为振动信号的傅里叶频域谱幅值,X(k)的定义如下:
在本申请中,主频带的第一均值和第一方差也通过公式(3)-公式(5)计算得到。
在故障振动信号频谱上,从左至右依次划分出一系列子频带,下一个子频带表示为与冗余分量相连的右边的子频带。若定义从左到右所有子频带为x1,x2,…xn,冗余分量为xm(1≤m<n),则下一个子频带为xm+1。
S4:在优化后的子频带上建立经验小波滤波器组,获得一系列经验模式。
依照经验小波的定义,在优化后的每个子频带上构造带通经验小波滤波器,形成滤波器组,以此得到每个子频带对应的经验模式。由Meyer小波理论来构造经验小波,经验小波的尺度函数φn(ω)和小波函数ψn(ω)定义为:
其中,x(ω)、ψn(ω)和φ1(ω)分别代表振动信号的离散化频谱、离散化频谱的小波函数和经验小波的尺度函数,和分别表示x(ω)、ψn(ω)和φ1(ω)的傅里叶变换,和则分别代表各自的共轭,F-1(·)表示傅里叶逆变换。
振动信号x(t)经过经验小波变换可以得到如下经验模式:
基于上述细节系数,振动信号x(t)被重构为:
S5:计算每个所述经验模式的裕度因子CL,筛选出所述裕度因子CL最大的经验模式进行包络谱分析。在本申请中,采用Hilbert包络谱对筛选出的经验模式进行分析。
裕度因子CL为对冲击脉冲信号较为敏感的参数,裕度因子CL的定义如下:
其中,XP=max{|xi|}表示峰值,故障振动信号用xi表示。
S6:从包络谱中提取出滚动轴承的故障特征频率进行分析,从而完成滚动轴承故障样本的故障诊断。
使用本申请提供的方法将原始的振动信号分解为一系列包含多种频率的经验模式,经过筛选出有效裕度因子从而提取出最佳经验模式,通过对最佳经验模式进行Hilbert包络谱分析,能够有效从强背景噪声和干扰信号之下提取到滚动轴承的故障特征频率,很大程度上提高了滚动轴承故障诊断的准确度,对于滚动轴承的故障诊断技术有着重要的意义,并且能够推动故障诊断技术进一步发展。
以上所述的仅是本申请的优选实施方式,本发明不限于以上实施例。可以理解,本领域技术人员在不脱离本发明的精神和构思的前提下直接导出或联想到的其他改进和变化,均应认为包含在本发明的保护范围之内。
Claims (5)
1.一种基于改进经验小波变换的滚动轴承故障诊断方法,其特征在于,所述滚动轴承故障诊断方法包括:
获取滚动轴承故障样本的振动信号,对所述振动信号进行FFT变换,得到频域信号的主频带;利用尺度空间法对所述主频带的频谱进行自适应划分得到一系列子频带;将所述一系列子频带的频域特征参数与所述主频带的频域特征参数进行比较,确定所述一系列子频带中的冗余分量,将每一个所述冗余分量与下一个子频带进行合并,以优化所述尺度空间法划分的子频带的频谱;在优化后的子频带上建立经验小波滤波器组,获得一系列经验模式;计算每个所述经验模式的裕度因子,筛选出所述裕度因子最大的经验模式进行包络谱分析;从包络谱中提取出滚动轴承的故障特征频率进行分析,从而完成滚动轴承故障样本的故障诊断。
2.根据权利要求1所述的滚动轴承故障特征提取方法,其特征在于,所述利用尺度空间法对所述主频带的频谱进行自适应划分得到一系列子频带,包括,
所述尺度空间法的定义如下:
离散函数f(m)的尺度空间L(m,t)定义为:
确定自适应寻找所述频带划分点的方法如下:
首先在所述离散化频谱中寻找所有的局部极小值点,每一初始极小值点都会对应一条曲线;在所述尺度空间L(ω,t)中,随着尺度参数t的增长,局部极小值点的数量减少,随着移动步长的变化,每个所述局部极小值点对应的曲线进行叠加生成一条直线,所述直线用Ci(i∈[1,N0])表示,直线长度为Li,采用Otsu方法自适应确定一个阈值T,若所述直线长度Li大于所述阈值T,则将对应的局部极小值点ωi作为所述频带划分点。
4.根据权利要求1所述的滚动轴承故障特征提取方法,其特征在于,所述在优化后的子频带上建立经验小波滤波器组,获得一系列经验模式,包括,依照经验小波的定义,在优化后的每个所述子频带上构造带通经验小波滤波器,形成滤波器组,由Meyer小波理论来构造所述经验小波,所述经验小波的尺度函数φn(ω)和小波函数ψn(ω)定义为:
其中,x(ω)、ψn(ω)和φ1(ω)分别代表振动信号的离散化频谱、离散化频谱的小波函数和经验小波的尺度函数,和分别表示x(ω)、ψn(ω)和φ1(ω)的傅里叶变换,和则分别代表各自的共轭,F-1(·)表示傅里叶逆变换;
所述振动信号x(t)经过经验小波变换可以得到如下经验模式:
基于上述细节系数,所述振动信号x(t)被重构为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010113154.7A CN111141520A (zh) | 2020-02-24 | 2020-02-24 | 一种基于改进经验小波变换的滚动轴承故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010113154.7A CN111141520A (zh) | 2020-02-24 | 2020-02-24 | 一种基于改进经验小波变换的滚动轴承故障诊断方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111141520A true CN111141520A (zh) | 2020-05-12 |
Family
ID=70527854
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010113154.7A Pending CN111141520A (zh) | 2020-02-24 | 2020-02-24 | 一种基于改进经验小波变换的滚动轴承故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111141520A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111623982A (zh) * | 2020-06-15 | 2020-09-04 | 大连理工大学 | 一种基于apewt和imomeda的行星齿轮箱早期故障诊断方法 |
CN111738155A (zh) * | 2020-06-23 | 2020-10-02 | 东北石油大学 | 往复压缩机轴承间隙故障诊断方法 |
CN112417966A (zh) * | 2020-10-22 | 2021-02-26 | 南京邮电大学 | 一种基于音频信号特征和ai算法的多维度机器故障检测方法 |
CN113111618A (zh) * | 2021-03-09 | 2021-07-13 | 电子科技大学 | 一种基于改进的经验小波变换的模拟电路故障诊断方法 |
CN113176460A (zh) * | 2021-04-28 | 2021-07-27 | 河海大学 | 一种基于改进经验小波变换的电能质量扰动信号检测方法 |
CN113589208A (zh) * | 2021-07-23 | 2021-11-02 | 深圳市联影高端医疗装备创新研究院 | 射频系统的频率确定方法、装置、磁共振设备及存储介质 |
CN114757060A (zh) * | 2022-06-15 | 2022-07-15 | 石家庄铁道大学 | 一种基于形态学函数构建尺度空间的轴承故障诊断方法 |
CN116358873A (zh) * | 2023-05-31 | 2023-06-30 | 厚普智慧物联科技有限公司 | 滚动轴承故障诊断方法、装置、设备及介质 |
CN117686226A (zh) * | 2024-02-04 | 2024-03-12 | 南京凯奥思数据技术有限公司 | 一种基于能量比和能量和的轴承自动故障诊断方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108375472A (zh) * | 2018-02-12 | 2018-08-07 | 武汉科技大学 | 基于改进经验小波变换的轴承故障诊断方法及系统装置 |
CN109525215A (zh) * | 2018-09-29 | 2019-03-26 | 长安大学 | 一种采用峭度谱确定子频带边界的经验小波变换方法 |
CN110006652A (zh) * | 2019-04-12 | 2019-07-12 | 燕山大学 | 一种滚动轴承故障诊断方法及系统 |
EP3550278A1 (en) * | 2018-04-06 | 2019-10-09 | United Technologies Corporation | Integrated vibe/odm fusion and trending analysis for prognostic health management of engine bearing |
-
2020
- 2020-02-24 CN CN202010113154.7A patent/CN111141520A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108375472A (zh) * | 2018-02-12 | 2018-08-07 | 武汉科技大学 | 基于改进经验小波变换的轴承故障诊断方法及系统装置 |
EP3550278A1 (en) * | 2018-04-06 | 2019-10-09 | United Technologies Corporation | Integrated vibe/odm fusion and trending analysis for prognostic health management of engine bearing |
CN109525215A (zh) * | 2018-09-29 | 2019-03-26 | 长安大学 | 一种采用峭度谱确定子频带边界的经验小波变换方法 |
CN110006652A (zh) * | 2019-04-12 | 2019-07-12 | 燕山大学 | 一种滚动轴承故障诊断方法及系统 |
Non-Patent Citations (4)
Title |
---|
张雨: "《时间序列的混沌和符号分析及实践》", 31 March 2007, 《国防科技大学》 * |
李兴林: "滚动轴承故障诊断技术现状及发展", 《2009全国青年摩擦学学术会议论文集》 * |
栗蕴琦: "经验小波变换和支持向量机在滚动轴承故障诊断中的应用研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
郭辉等: "基于峭度与互相关的IEWT轴承故障诊断方法研究", 《铁道科学与工程学报》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111623982A (zh) * | 2020-06-15 | 2020-09-04 | 大连理工大学 | 一种基于apewt和imomeda的行星齿轮箱早期故障诊断方法 |
CN111738155A (zh) * | 2020-06-23 | 2020-10-02 | 东北石油大学 | 往复压缩机轴承间隙故障诊断方法 |
CN112417966A (zh) * | 2020-10-22 | 2021-02-26 | 南京邮电大学 | 一种基于音频信号特征和ai算法的多维度机器故障检测方法 |
CN112417966B (zh) * | 2020-10-22 | 2022-08-26 | 南京邮电大学 | 一种基于音频信号特征和ai算法的多维度机器故障检测方法 |
CN113111618B (zh) * | 2021-03-09 | 2022-10-18 | 电子科技大学 | 一种基于改进的经验小波变换的模拟电路故障诊断方法 |
CN113111618A (zh) * | 2021-03-09 | 2021-07-13 | 电子科技大学 | 一种基于改进的经验小波变换的模拟电路故障诊断方法 |
CN113176460A (zh) * | 2021-04-28 | 2021-07-27 | 河海大学 | 一种基于改进经验小波变换的电能质量扰动信号检测方法 |
CN113589208A (zh) * | 2021-07-23 | 2021-11-02 | 深圳市联影高端医疗装备创新研究院 | 射频系统的频率确定方法、装置、磁共振设备及存储介质 |
CN114757060A (zh) * | 2022-06-15 | 2022-07-15 | 石家庄铁道大学 | 一种基于形态学函数构建尺度空间的轴承故障诊断方法 |
CN116358873A (zh) * | 2023-05-31 | 2023-06-30 | 厚普智慧物联科技有限公司 | 滚动轴承故障诊断方法、装置、设备及介质 |
CN116358873B (zh) * | 2023-05-31 | 2023-09-15 | 厚普清洁能源(集团)股份有限公司 | 滚动轴承故障诊断方法、装置、设备及介质 |
CN117686226A (zh) * | 2024-02-04 | 2024-03-12 | 南京凯奥思数据技术有限公司 | 一种基于能量比和能量和的轴承自动故障诊断方法及系统 |
CN117686226B (zh) * | 2024-02-04 | 2024-04-16 | 南京凯奥思数据技术有限公司 | 一种基于能量比和能量和的轴承自动故障诊断方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111141520A (zh) | 一种基于改进经验小波变换的滚动轴承故障诊断方法 | |
Qin et al. | M-band flexible wavelet transform and its application to the fault diagnosis of planetary gear transmission systems | |
Wang et al. | Complete ensemble local mean decomposition with adaptive noise and its application to fault diagnosis for rolling bearings | |
Singh et al. | Pseudo-fault signal assisted EMD for fault detection and isolation in rotating machines | |
Wang et al. | Fault diagnosis of diesel engine based on adaptive wavelet packets and EEMD-fractal dimension | |
Lei et al. | A review on empirical mode decomposition in fault diagnosis of rotating machinery | |
Yang et al. | Vibration feature extraction techniques for fault diagnosis of rotating machinery: a literature survey | |
CN113375939B (zh) | 基于svd和vmd的机械件故障诊断方法 | |
Popescu et al. | Basic tools for vibration analysis with applications to predictive maintenance of rotating machines: an overview | |
Barbini et al. | Phase editing as a signal pre-processing step for automated bearing fault detection | |
CN111582128B (zh) | 一种基于狼群参数化联合字典的机械故障稀疏表示方法 | |
CN116750290B (zh) | 一种袋式包装机的运行状态智能监测方法 | |
Zhang et al. | Improved local cepstrum and its applications for gearbox and rolling bearing fault detection | |
CN112098102A (zh) | 一种基于ewt-scwt的内燃机异响识别与诊断方法 | |
Liu et al. | Asymmetric penalty sparse model based cepstrum analysis for bearing fault detections | |
CN110987431B (zh) | 一种基于tqwt辅助spc的轴承状态监测及故障诊断的方法 | |
Antoni et al. | Cyclostationarity in condition monitoring: 10 years after | |
CN109630908B (zh) | 一种多次降噪的管道泄漏定位方法 | |
Lv et al. | Longitudinal synchroextracting transform: A useful tool for characterizing signals with strong frequency modulation and application to machine fault diagnosis | |
Gan et al. | Fault feature enhancement for rotating machinery based on quality factor analysis and manifold learning | |
Han et al. | Second-order multisynchrosqueezing wavelet transform for bearing fault detection | |
CN110147637A (zh) | 基于小波及谐波成分贪婪稀疏识别的碰摩故障诊断方法 | |
Li et al. | A new method of nonlinear feature extraction for multi-fault diagnosis of rotor systems | |
CN115655719A (zh) | 一种轴承振动信号分阶段降噪方法及轴承故障识别方法 | |
CN114894478A (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200512 |