CN103971019B - 一种基于几何特征的工件加工能耗预测方法 - Google Patents
一种基于几何特征的工件加工能耗预测方法 Download PDFInfo
- Publication number
- CN103971019B CN103971019B CN201410221458.XA CN201410221458A CN103971019B CN 103971019 B CN103971019 B CN 103971019B CN 201410221458 A CN201410221458 A CN 201410221458A CN 103971019 B CN103971019 B CN 103971019B
- Authority
- CN
- China
- Prior art keywords
- energy consumption
- workpiece
- geometric properties
- feature
- model
- 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.)
- Expired - Fee Related
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Numerical Control (AREA)
Abstract
本发明公开了一种基于几何特征的工件加工能耗预测方法,包括:根据CAD模型中集成的几何特征信息对工件CAD模型拆分建立特征树,用特征类树表示各几何特征类及其层次关系;利用物元法建立工件几何特征物元模型,再从物元模型中提取工件几何特征和加工方式,建立映射关系,结合加工过程中机床各能耗部件的能耗特性,建立基于几何特征的工件加工能耗模型,并通过解析几何特征的NC代码获取对应数控机床各能耗部件的运行参数;将各参数信息带入能耗模型,计算出加工此几何特征的能耗值,再通过对每个几何特征的能耗求和获取工件加工过程的总能耗。本发明围绕工件几何特征的能耗属性展开,只需根据拆分出的几何特征,就可以对工件加工过程能耗进行预测。
Description
技术领域
本发明涉及机械制造加工领域,尤其涉及一种基于几何特征的工件加工能耗预测方法。
背景技术
机械加工过程主要是对金属或其他材料的坯料或工件进行加工,使之获得所要求的几何形状、尺寸精度和表面质量。随着能效制造及低碳制造等一系列降低制造过程能耗、减少环境排放的制造模式的推进,实施节能优化已成为机械制造业发展的趋势。机械加工过程中工件能耗预测是其降低能耗、提高能效的首要步骤之一,近年来正在国内外迅速兴起并逐渐成为研究热点问题。
工件加工过程的能量消耗包括错综复杂的能源流、物料流和废物流,并贯穿于工件的各种加工工艺、加工设备和辅助设施中,具有复杂的动态特性。CN103235554A公开《一种基于NC代码的数控车床加工工件能耗获取方法》,其公开的方案是通过对数控机床NC代码的解析获取各能耗部件运行参数,并通过实验获取基础功率数据,最后将上述运行参数和基础功率输入到建立的数控机床各能耗部件的能耗模型中,通过计算获取工件的能耗。CN102621932A公开的《一种数控机床服役过程的能量消耗预测方法》,其公开的方案是根据工件加工工艺,将数控机床服役过程分为启动、空载和加工三个子过程,对每个子过程进行能耗预测,该方法只适用于工件加工工艺信息详细,且需要人工划分服役过程。
现有预测方法中,工件能耗预测主要通过实验测量获取机床加工过程中的主要功率和工艺参数,再利用模型公式对工件能耗进行估算,目前还缺乏利用工件几何特征来预测其能耗的方法。
发明内容
本发明要解决的技术问题在于针对现有技术中的缺陷,提供一种基于几何特征的工件加工能耗预测方法,使用本方法能方便的获取工件加工能耗。
本发明解决其技术问题所采用的技术方案是:一种基于几何特征的工件加工能耗预测方法,包括以下步骤:
1)根据cad模型中集成的几何特征信息对工件CAD模型拆分建立特征类树,所述特征类树用于表示各几何特征类及其层次关系和结构;
2)利用物元法建立工件几何特征物元模型,从物元模型中识别并提取工件几何特征和加工方式,建立几何特征和加工方式间的映射关系;
3)根据加工过程中机床各能耗部件的能耗特性,通过物元法建立的工件几何特征物元模型对工件各几何特征能量信息进行提取;
4)根据几何特征和加工方式间的映射关系和工件各几何特征能量信息,建立基于CAD模型几何特征的工件加工能耗模型;
5)通过解析几何特征的NC代码获取对应数控机床各能耗部件的运行参数;
6)根据各能耗部件的运行参数和能耗模型,计算各加工工件几何特征的能耗值,再通过对每个几何特征的能耗求和获取工件加工过程的总能耗。
按上述方案,所述步骤1)中利用CAD模型中集成的几何特征信息,建立工件几何特征类树的具体方法为:
首先按照工件特征的复杂程度将工件几何特征分为简单特征和复杂特征,简单特征是特征的最小单元,简单特征是独立的形状结构;复杂特征为简单特征的组合结构,然后再将复杂特征按加工要求分为基本特征和组合特征,将简单特征分为主特征和辅特征两种,并以几何特征类树的结构形式将所有特征递进排列,每个结点代表一个几何特征类。
按上述方案,所述步骤2)中建立几何特征与加工行为之间的映射关系为:
利用物元法建立工件几何特征的能耗物元模型R=(P,F,E),物元表示如下:
其中P代表工件,F代表工件几何特征,E代表特征能耗;所述特征能耗信息包括各几何特征对应的加工方法,还包括加工方法中选用的设备及相应的加工参数。
按上述方案,所述步骤2)中工件几何特征物元模型为:
按上述方案,步骤4)中基于几何特征的数控机床加工过程能耗模型为:
设工件由n个几何特征MF构成:MF={MF1,MF2,…,MFn},其某一几何特征和工件能耗的数学模型可以表述为:
其中EMF可表示为:
于是可得工件加工过程总能耗:
其中代表机床加工几何特征i时工序j的能耗,m表示加工某一几何特征时所需的加工工序数。
本发明产生的有益效果是:
1.本发明围绕工件几何特征的能耗属性展开,只需根据拆分出的几何特征,结合数控机床NC代码与加工几何特征有关的参数信息,分别分析工件主辅特征的能量消耗状况及与加工过程能耗间的映射关系,进而建立基于CAD模型几何特征的工件加工能耗预测模型,就可以对工件加工过程能耗进行预测;
2.本发明还可以提供分解到工件各加工特征的能耗数据,这些数据信息可为机床加工过程节能优化分析提供支持。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1是本发明实施例的流程示意图;
图2是工件几何特征分类图;
图3是本发明实施例的机床刀具运行轨迹图;
图4是本发明实施例的机床进给轴直线插补图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
如图1所示,本发明公开了一种基于几何特征的工件加工能耗预测方法,首先,利用CAD模型中集成的几何特征信息对工件CAD模型拆分建立特征树,用几何特征类树表示各几何特征类及其层次关系和结构;然后,先利用物元法建立工件几何特征物元模型,再从物元模型中识别和提取工件几何特征和加工方式,建立映射关系,结合加工过程中机床各能耗部件的能耗特性,建立基于CAD模型几何特征的工件加工能耗模型,并通过解析几何特征的NC代码获取对应数控机床各能耗部件的运行参数;最后,将各参数信息带入能耗模型,计算出加工此几何特征的能耗值,再通过对每个几何特征的能耗求和获取工件加工过程的总能耗。
1、工件CAD模型拆分建立特征树:
工件几何特征分为简单特征和复杂特征,复杂特征可分为基本特征和组合特征,而简单特征又可分为主特征和辅特征两种,基本特征由一个主特征和多个辅特征组成,如轴段和孔段等,组合特征由多个基本特征组成,如同轴孔和均布孔等,主特征通常表示用以构成工件的基本形体结构,如孔、槽等;辅特征则依附于主特征,构成零件模型的完整结构,如倒角、退刀槽等。
设工件P由M个复杂特征组成,每一个复杂特征又是由N个组合特征构成,每个组合特征由K个基本特征组成,每个基本特征由一个主特征和U个辅特征组成,经过层层分解,一个复杂工件可分割为由若干个主特征和辅特征组合排列而构成,如图2所示。
2、识别和提取工件几何特征和加工方式,建立工件几何特征的能耗模型:
2.1、建立几何特征与加工行为之间的映射关系为:
一个工件(记为P),可以通过各种各样的几何特征(记为F)来体现。描述事物实体的特征称为实义特征.如长、宽、质量等。事物关于某一特征的数量、程度或范围等称为该事物关于这一特征的量值。特征F的取值范围.称为它的量域,记作V。给定事物的名称P,它关于特征F的量值为V,以有序三元组作为描述事物的基本元,简称为物元.记作R(P,F,V)。
对步骤1中拆分出的主特征和辅特征利用上述介绍的物元法建立工件加工特征的能耗物元模型R=(P,F,E),然后对工件主特征和辅特征进行提取,物元表示如下:
其中P代表工件,F代表工件加工特征,E代表特征能耗;
进一步将上述几何特征信息和特征能耗信息利用物元法进行分解。几何特征信息主要包括几何信息和加工信息,特征能耗信息包括各几何特征对应的加工方法,还包括加工方法中选用的设备及相应的加工参数,为了能够建立几何特性到能耗的映射关系,还需要对物元模型进行进一步分解。利用物元的逐层迭代,建立三层的统一物元模型,将工件几何特征与各具体的加工参数进行关联,统一单元物元模型如下:
工件几何特征不但具有几何信息,还必须具有加工信息,具有可加工性,因此工件的几何特征与加工行为之间存在着对应关系。
2.2、数控机床几何特征NC代码的解析:
数控机床一般采用NC代码(又称数控程序)编程,NC代码编程语言属过程控制语言,其语句是控制数控机床运作的指令集。在数控机床用ISO编码表中,英文字母大都赋给了确切的含义。一个英文字母后边跟上一定范围的数码就构成了一个N C代码指令,主要包括了G代码、M代码、S代码、F代码、X代码、Y代码和Z代码,数控加工工件NC代码与机床运行的关联性如表1所示。
表1
2.3、建立基于几何特征的数控机床加工过程能耗模型为:
设工件由n个几何特征MF构成:MF={MF1,MF2,…,MFn},其中MF可表示为:
MF={X(t),Y(t),M(t),S(t)}
其中,X(t)为切削用量特性;Y(t)为加工设备特性;M(t)为加工工件的材料特性;S(t)为工件形状特性。
某一几何特征与加工能耗的数学模型可以表述为:
其中EMF可进一步表示为:
表示几何特征MFi需进行m工序加工,对每步工序的能耗求和就可以得出该几何特征的能耗。
于是可得工件加工过程总能耗:
其中代表机床加工特征i时工序j的能耗。
2.4、基于CAD模型几何特征的工件加工能耗的获取:
工件加工过程中能耗主要由数控机床主轴电机切削能耗EC和进给电机能耗EF构成,主轴电机能耗主要与工件材料去除过程相关,进给电机消耗的能量用于工作台和切削刀具在给定的进给速度下的进给。
1)主轴电机切削能耗模型为:
其中切削功率Pc和切削时间参数需要获取,切削时间t可以通过NC代码文件中获取的刀具路径和进给速度VC来计算。切削功率Pc可以表示式(5),
PC=FC×V (5)
其中Fc表示切削力,切削力主要与切削厚度和长度有关,在实际切削过程中切削力的影响因素很多,主要有工件、刀具材料、切削用量、刀具几何参数、刀具磨损状态以及切削液等,如标准切削实验,车削时的Fc可以表示为式(6),
其中ap、f、vc分别为背吃刀量、进给速度、切削速度,KFc详见《切削用量手册》表1.28,CFc为修正系数,xFc、yFc、nFc为指数,具体常用数值可见表2所示。
表2
2)进给电机能耗模型为:
数控机床进给电机一般包括X轴进给电机和Y轴进给电机以及Z轴进给电机,进给电机能耗可以表达为(式7),
其中Pf,tFs、tFe分别代表第i个进给轴的功率、开始时间和结束时间。
进给系统能耗一般包括两部分,一部分为工件加工之前,三轴进行快进和进给运动消耗的能量主要是为了进行对刀和将夹具移动到指定位置等准备工作,另一部分为工件加工过程中给定进给轴速度,平面移动加工工件所消耗的能量
①的获取:
主要与进给速度、进给时间和刀具运行轨迹有关,现结合图3讲述具体获取方法,假设刀具的路径为A→C→D→B,进给系统快给速度为Vr,首先三个方向的进给电机同时运转,刀具从点A移动到点C;然后Z向进给电机停止运转,X和Y向进给电机继续运转,刀具从点C移动到点D;最后只有X向电机运转,刀具从点D移动到点B。
假设三个进给电机的功率分别为则可由(式8)求得,
由于三个进给电机都是相同的速度Vr,因此(式8)可以改写为(式9),
其中:
②的获取:
和相似,也主要与进给速度、进给时间和刀具运行轨迹有关,但在加工过程中刀具不能严格地按照要求加工的曲线运动,只能用折线轨迹逼近所要加工的曲线,常用计算规则有直线插补计算规则、圆弧插补计算规则以及特殊循环计算规则。
本发明以直线插补计算规则为例,现结合图4来进一步解释的获取方法,在点A(0,0)、点B(xb,yb)和点C(xc,yc)之间进行直线插补,进给系统X轴电机和Y轴电机的能耗如(式11)所示:
其中分别为进给轴X和进给轴Y在电机速度vx和vy下的功率,vx和vy可以通过进给速度vf求得,如(式12)所示:
在忽略掉进给轴加速和减速的影响后,进给矢量速度vf为刀具从点A运动到点B时的速度,由于各分量上的运动时间是相同的,则(式11)可进一步改写为(式13):
其中时间可表示为(式14)所示:
公式(10)和(14)中的参数可以通过前述步骤中对NC代码的解析获取。
3、机械加工过程工件能耗的获取:
将式(5)和式(6)带入式(4)可获得在工件加工特征i下主轴电机切削能耗模型:
将式(8)~(14)经过整理后可得在工件加工特征i进给轴电机能耗模型:
最后将式(15)和式(16)带入式(3)可得基于加工特征的工件加工过程总能耗模型:
其中ECij和EFij分别代表机床加工特征i时工序j的主轴能耗和进给轴能耗。
通过上述基于几何特征的工件加工能耗获取方法可以看出,本发明可以应用于数控加工工件能耗评估与预测,可为工件加工过程能耗的节能优化提供一定的支持。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (3)
1.一种基于几何特征的工件加工能耗预测方法,其特征在于,包括以下步骤:
1)根据CAD模型中集成的几何特征信息对工件CAD模型拆分建立特征类树,所述特征类树用于表示各几何特征类及其层次关系和结构;
2)利用物元法建立工件几何特征物元模型,从物元模型中识别并提取工件几何特征和加工方式,建立几何特征和加工方式间的映射关系;
3)根据加工过程中机床各能耗部件的能耗特性,通过物元法建立的工件几何特征物元模型对工件各几何特征能量信息进行提取;
4)根据几何特征和加工方式间的映射关系和工件各几何特征能量信息,建立基于CAD模型几何特征的工件加工能耗模型;
5)通过解析几何特征的NC代码获取对应数控机床各能耗部件的运行参数;
6)根据各能耗部件的运行参数和能耗模型,计算各加工工件几何特征的能耗值,再通过对每个几何特征的能耗求和获取工件加工过程的总能耗。
2.根据权利要求1所述的加工能耗预测方法,其特征在于,所述步骤1)中利用CAD模型中集成的几何特征信息,建立工件几何特征类树的具体方法为:
首先按照工件特征的复杂程度将工件几何特征分为简单特征和复杂特征,其中简单特征是特征的最小单元;然后再将复杂特征按加工要求分为基本特征和组合特征,将简单特征分为主特征和辅特征两种,并以几何特征类树的结构形式将所有特征递进排列,每个结点代表一个几何特征类。
3.根据权利要求1所述的加工能耗预测方法,其特征在于,所述步骤2)中建立几何特征与加工方式之间的映射关系为:
利用物元法建立工件几何特征的能耗物元模型R=(P,F,E),物元表示如下:
其中P代表工件,F代表工件几何特征,E代表特征能耗;所述特征能耗信息包括各几何特征对应的加工方法,还包括加工方法中选用的设备及相应的加工参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410221458.XA CN103971019B (zh) | 2014-05-23 | 2014-05-23 | 一种基于几何特征的工件加工能耗预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410221458.XA CN103971019B (zh) | 2014-05-23 | 2014-05-23 | 一种基于几何特征的工件加工能耗预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103971019A CN103971019A (zh) | 2014-08-06 |
CN103971019B true CN103971019B (zh) | 2017-02-08 |
Family
ID=51240507
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410221458.XA Expired - Fee Related CN103971019B (zh) | 2014-05-23 | 2014-05-23 | 一种基于几何特征的工件加工能耗预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103971019B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI571820B (zh) * | 2014-11-06 | 2017-02-21 | 財團法人資訊工業策進會 | 工具機耗電量預測系統與方法 |
CN105278455B (zh) * | 2015-11-20 | 2017-11-24 | 华中农业大学 | 基于Gcode的可持续加工操作执行顺序及刀具路径优化方法 |
CN106475908A (zh) * | 2016-11-08 | 2017-03-08 | 上海大学 | 基于标准g代码的随动磨削过程机床运行能耗预测方法 |
CN107193258B (zh) * | 2017-06-22 | 2019-10-18 | 重庆大学 | 面向能耗的数控加工工艺路线与切削参数集成优化方法 |
EP3582044B1 (en) * | 2018-06-14 | 2022-04-13 | Sandvik Machining Solutions AB | Machining based on strategies selected from a database |
CN109034452A (zh) * | 2018-06-20 | 2018-12-18 | 上海安悦节能技术有限公司 | 用于汽车零配件生产企业的能耗预测方法 |
CN110560922A (zh) * | 2019-08-22 | 2019-12-13 | 浙江科技学院 | 一种基于几何特征的激光切割的加工能耗预测方法 |
CN113344244B (zh) * | 2021-04-29 | 2023-04-07 | 武汉科技大学 | 一种数字孪生驱动的工件加工能耗预测装置 |
CN116842740B (zh) * | 2023-07-08 | 2024-05-07 | 四川大学 | 基于物元模型的加工特征智能识别方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1102138A1 (de) * | 1999-11-22 | 2001-05-23 | Paul Schindele | Verfahren zum Herstellen von 3-dimensionalen Formteilen |
US6298279B1 (en) * | 1997-04-28 | 2001-10-02 | Seiko Instruments Inc. | CAD/CAM apparatus and machining apparatus for processing position/force information |
CN102621932A (zh) * | 2012-05-02 | 2012-08-01 | 重庆大学 | 一种数控机床服役过程的能量消耗预测方法 |
CN103235554A (zh) * | 2013-03-29 | 2013-08-07 | 重庆大学 | 一种基于nc代码的数控车床加工工件能耗获取方法 |
-
2014
- 2014-05-23 CN CN201410221458.XA patent/CN103971019B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6298279B1 (en) * | 1997-04-28 | 2001-10-02 | Seiko Instruments Inc. | CAD/CAM apparatus and machining apparatus for processing position/force information |
EP1102138A1 (de) * | 1999-11-22 | 2001-05-23 | Paul Schindele | Verfahren zum Herstellen von 3-dimensionalen Formteilen |
CN102621932A (zh) * | 2012-05-02 | 2012-08-01 | 重庆大学 | 一种数控机床服役过程的能量消耗预测方法 |
CN103235554A (zh) * | 2013-03-29 | 2013-08-07 | 重庆大学 | 一种基于nc代码的数控车床加工工件能耗获取方法 |
Non-Patent Citations (2)
Title |
---|
《基于加工特征的零件设计方法研究》;姜斌;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑 》;20120515(第5期);第C029-146页 * |
《机械加工系统节能降噪型综合任务分配模型及其应用》;曹华军;《机械工程学报》;20060531;第42卷(第5期);第97-102页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103971019A (zh) | 2014-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103971019B (zh) | 一种基于几何特征的工件加工能耗预测方法 | |
Liu et al. | A method for predicting the energy consumption of the main driving system of a machine tool in a machining process | |
CN102621932B (zh) | 一种数控机床服役过程的能量消耗预测方法 | |
Aramcharoen et al. | Critical factors in energy demand modelling for CNC milling and impact of toolpath strategy | |
CN103941644B (zh) | 一种基于时间参数的数控铣床能耗预测方法 | |
CN103699056B (zh) | 高速高精度数控加工的小线段实时平滑过渡插补方法 | |
Ma et al. | An energy consumption optimization strategy for CNC milling | |
WO2023272836A1 (zh) | 一种基于加工知识的数控工艺设计及优化方法 | |
CN106094729B (zh) | 一种基于制造特征的零件机械加工过程碳排放量化方法 | |
CN105488282A (zh) | 一种基于动态加工特征的切削参数分段和变切深优化方法 | |
CN104133417A (zh) | 叶片式流体机械数控加工切削力的快速预测方法 | |
CN102179727A (zh) | 机床主传动系统加工过程能耗信息在线检测方法 | |
CN104460515A (zh) | 一种基于后置处理五轴刀具长度补偿方法 | |
Shao et al. | A virtual machining model for sustainability analysis | |
CN106777877A (zh) | 一种基于隐含碳能的机床装备产品碳排放量化方法 | |
Zheng et al. | KBE-based stamping process paths generated for automobile panels | |
CN104698964B (zh) | 一种基于映射的复杂曲面数控加工运动分析方法 | |
Wang et al. | Intelligent manufacturing system of impeller for computer numerical control (CNC) programming based on KBE | |
He et al. | A framework of scheduling models in machining workshop for green manufacturing | |
Xu et al. | Modeling of process parameter selection with mathematical logic for process planning | |
Jasthi et al. | Studies on process plan representation in CAPP systems | |
Shao et al. | Tool path generation method for five-axis flank milling of corner by considering dynamic characteristics of machine tool | |
Rachieru et al. | An improved method for risk evaluation in failure modes and effects analysis of CNC lathe | |
CN106903555A (zh) | 一种刀具-工件切触区域的判定方法 | |
Li et al. | Development of NC Program Simulation Software Based on AutoCAD |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170208 Termination date: 20180523 |