CN110208708A - 一种采用二分法仿真锂离子电池极限充放电电流的方法 - Google Patents
一种采用二分法仿真锂离子电池极限充放电电流的方法 Download PDFInfo
- Publication number
- CN110208708A CN110208708A CN201910517159.3A CN201910517159A CN110208708A CN 110208708 A CN110208708 A CN 110208708A CN 201910517159 A CN201910517159 A CN 201910517159A CN 110208708 A CN110208708 A CN 110208708A
- Authority
- CN
- China
- Prior art keywords
- emulation
- max
- charging
- voltage
- time
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/367—Software therefor, e.g. for battery testing using modelling or look-up tables
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/385—Arrangements for measuring battery or accumulator variables
- G01R31/387—Determining ampere-hour charge capacity or SoC
- G01R31/388—Determining ampere-hour charge capacity or SoC involving voltage measurements
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Secondary Cells (AREA)
Abstract
本发明公开了一种采用二分法仿真锂离子电池极限充放电电流的方法,获取被测电池常温倍率放电的电压‑时间实验数据,根据被测电池的参数建立被测电池的三维电化学‑热耦合模型,并进行常温倍率放电仿真获取常温倍率放电的电压‑时间仿真数据,比较被测电池的电压‑时间实验数据与三维电化学‑热耦合模型的电压‑时间仿真数据对三维电化学‑热耦合模型进行修正得到修正后的三维电化学‑热耦合模型,利用修正后的三维电化学‑热耦合模型,设置被测电池的初始电压、充放电停止条件与计算时间,采用二分法进行极限充放电仿真;通过仿真测试极限充放电电流,不仅缩短测试时间,而且节约实验成本。
Description
技术领域
本发明属于锂电池管理技术领域,具体涉及一种采用二分法仿真锂离子电池极限充放电电流的方法。
背景技术
随着科技的不断发展以及环境和能源危机的不断加重,新能源汽车的应用也越来越广泛。动力电池是新能源汽车中常用的动力来源,是汽车中最重要的部件之一,动力电池性能的好坏直接影响汽车的续航能力和安全性。电动汽车在急加速或急刹车时,需要电池在短时间内放出或吸收很大的功率,将导致流过电池的电流瞬间增大或减小,有可能引起电池的过充、过放等滥用行为。为了实现电池冲放电功率的优化管理,需要准确估计电池的极限冲放电电流,以及时对电池的冲放电电流进行限制,防止电池发生过充或过放。电池的极限充放电电流一般用一段时间内(如10s)电池能充入或放出的最大电流来表示。
目前电池极限冲放电电流的确定方法主要有两种,一种是实验测试方法,对电池进行不同温度和不同荷电状态下的充放电测试,得到电池在不同状态下的极限充放电电流,这种方法需要的实验量很大,实施起来需要的测试资源很多,此外还受电池一致性的影响。另外一种方法是基于模型的估计方法,如中国发明专利(申请号201611201357.1)提出了一种基于扩展等效电路模型的电池极限充放电电流估计方法,此方法的主要思路是比较等效电路模型计算得到的电池电压与截止电压,放电若小于截止电压就减小电流,反之就增大电流,充电则与之相反。此种方法,充电与放电各需要一个模型;此外,电流增大与减小的程度没有具体的控制方法,会影响估计效率。
发明内容
本发明的目的就是针对上述技术的不足,提供一种既能缩短测试时间又能节约实验成本的采用二分法仿真锂离子电池极限充放电电流的方法,既可用于充电极限电流仿真,也可用于放电极限电流仿真。
为实现上述目的,本发明所设计的采用二分法仿真锂离子电池极限充放电电流的方法,包括如下步骤:
S1)获取被测电池常温倍率放电的电压-时间实验数据,并在恒温箱中对被测电池进行不同温度下的SOC-OCV实验,获取被测电池在不同温度下的SOC-OCV实验数据;
S2)根据被测电池的参数建立被测电池的三维电化学-热耦合模型,并进行常温倍率放电仿真获取常温倍率放电的电压-时间仿真数据;
S3)比较步骤S1)中被测电池的电压-时间实验数据与步骤S2)中三维电化学-热耦合模型的电压-时间仿真数据对三维电化学-热耦合模型进行修正得到修正后的三维电化学-热耦合模型;
S4)利用修正后的三维电化学-热耦合模型,设置被测电池的初始电压、充放电停止条件与计算时间,采用二分法进行极限充放电仿真;
其中,被测电池的初始电压为被测电池的开路电压,为步骤S1)SOC-OCV实验数据中测的开路电压OCV值;充电的停止条件为被测电池的电压大于被测电池上限截止电压,放电的停止条件为被测电池的电压小于被测电池下限截止电压;计算时间为实验要求极限充电持续时间t1和实验要求极限放电持续时间t1’,记充电仿真计算停止时实际计算时间为t2和放电仿真计算停止时实际计算时间为t2’,充电时满足t2=t1的电流I为Imin、满足t2<t1的电流I为Imax,预设的停止标准Imax-Imin<a,充电电流I在[Imin,Imax]之间取值,放电时满足t2’=t1’的电流I’为Imin’、满足t2’<t1’的电流I’为Imax’,预设的停止标准Imax’-Imin’<a,放电电流I’在[Imin’,Imax’]之间取值。
进一步地,所述步骤S4)中,利用修正后的三维电化学-热耦合模型采用二分法进行极限充电仿真具体过程为:首先预设一个I0进行充电仿真计算,根据计算停止时实际计算时间t2与t1的大小关系,来减小或增大I0,若t2=t1,则增大I0,直到t2<t1时停止;若t2<t1,则减小I0直到t2=t1时停止,这些I0值中,最后两个中的较小者为Imin、较大者为Imax;然后将[Imin,Imax]均分为n份获得n+1个电流值,对这n+1个电流值依次进行充电仿真,根据每个电流值充电仿真持续时间t2与t1的大小关系,判断得出新的Imax-1与Imin-1,n+1个I值所对应的t2中最后一个满足t2=t1的I为新的Imin即Imin-1、第一个满足t2<t1的I为新的Imax即Imax-1;依次重复进行,然后继续将[Imin-1,Imax-1]均分为n份获得n+1个电流值,对这n+1个电流值依次进行充电仿真,根据每个电流值充电仿真持续时间t2与t1的大小关系,判断得出新的Imax-2与Imin-2,直至新的Imax-i与Imin-i满足Imax-i-Imin-i<a停止计算,最后一次得到的Imin-i即为所需状态下充电的极限电流值,n取值为大于等于2的整数;
放电的极限电流值的仿真方法同充电的极限电流值的仿真方法。
进一步地,所述步骤S2)中,参数包括设计参数、电极动力学参数、热物性参数等,设计参数包括被测电池的长、宽、高等。
与现有技术相比,本发明具有以下优点:本发明采用二分法仿真锂离子电池极限充放电电流的方法,通过仿真测试极限充放电电流,不仅缩短测试时间,而且节约实验成本;本发明提供的方法只需一个模型,即可用于充电极限电流仿真,也可用于放电极限电流仿真,且采用二分法的思想,极大提高测试效率。
附图说明
图1为本发明采用二分法仿真锂离子电池极限充放电电流的方法流程图;
图2为本发明实施例23度10s极限电流放电仿真结果与实验结果的比较图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步的详细说明。
如图1所示一种采用二分法仿真锂离子电池极限充放电电流的方法,包括如下步骤:
S1)获取被测电池常温倍率放电的电压-时间实验数据,并在恒温箱中对被测电池进行不同温度下的SOC-OCV实验,获取被测电池在不同温度下的SOC-OCV实验数据;
S2)根据被测电池的参数建立被测电池的三维电化学-热耦合模型,并进行常温倍率放电仿真获取常温倍率放电的电压-时间仿真数据;其中,参数包括设计参数、电极动力学参数、热物性参数等,设计参数包括被测电池的长、宽、高等;
S3)比较步骤S1)中被测电池的电压-时间实验数据与步骤S2)中三维电化学-热耦合模型的电压-时间仿真数据对三维电化学-热耦合模型进行修正得到修正后的三维电化学-热耦合模型;
S4)利用修正后的三维电化学-热耦合模型,设置被测电池的初始电压、充放电停止条件与计算时间,采用二分法进行极限充放电仿真;其中,被测电池的初始电压为被测电池的开路电压,为步骤S1)SOC-OCV实验数据中测的开路电压OCV值;充电的停止条件为被测电池的电压大于被测电池上限截止电压,放电的停止条件为被测电池的电压小于被测电池下限截止电压;计算时间为实验要求极限充电持续时间t1和实验要求极限放电持续时间t1’,记充电仿真计算停止时实际计算时间为t2和放电仿真计算停止时实际计算时间为t2’,且充电时满足t2=t1的电流I为Imin、满足t2<t1的电流I为Imax,预设的停止标准Imax-Imin<a,充电电流I在[Imin,Imax]之间取值,放电时满足t2’=t1’的电流I’为Imin’、满足t2’<t1’的电流I’为Imax’,预设的停止标准Imax’-Imin’<a,放电电流I’在[Imin’,Imax’]之间取值
利用修正后的三维电化学-热耦合模型采用二分法进行极限充电仿真具体过程为:利用修正后的三维电化学-热耦合模型采用二分法进行极限充电仿真具体过程为:首先预设一个I0进行充电仿真计算,根据计算停止时实际计算时间t2与t1的大小关系,来减小或增大I0,若t2=t1,则增大I0,直到t2<t1时停止;若t2<t1,则减小I0直到t2=t1时停止,这些I0值中,最后两个中的较小者为Imin、较大者为Imax;然后将[Imin,Imax]均分为n份获得n+1个电流值,对这n+1个电流值依次进行充电仿真,根据每个电流值充电仿真持续时间t2与t1的大小关系,判断得出新的Imax-1与Imin-1,n+1个I值所对应的t2中最后一个满足t2=t1的I为新的Imin即Imin-1、第一个满足t2<t1的I为新的Imax即Imax-1;依次重复进行,然后继续将[Imin-1,Imax-1]均分为n份获得n+1个电流值,对这n+1个电流值依次进行充电仿真,根据每个电流值充电仿真持续时间t2与t1的大小关系,判断得出新的Imax-2与Imin-2,直至新的Imax-i与Imin-i满足Imax-i-Imin-i<a停止计算,最后一次得到的Imin-i即为所需状态(温度与SOC)下充电的极限电流值,本实施例中,n取值为大于等于2的整数;
放电的极限电流值的仿真方法同充电的极限电流值的仿真方法,在此不再赘述。
下面以48Ah三元锂离子电池的23度10s极限电流放电(10%、20%、35%、50%、65%、80%、95%SOC)为例,进行说明。
S1)获取被测电池常温倍率放电的电压-时间实验数据,并在恒温箱中对被测电池进行不同温度下的SOC-OCV实验,获取被测电池在不同温度下的SOC-OCV实验数据;
S2)根据被测电池的参数建立被测电池的三维电化学-热耦合模型,并进行常温倍率放电仿真获取常温倍率放电的电压-时间仿真数据;其中,参数包括设计参数、电极动力学参数、热物性参数等,设计参数包括被测电池的长、宽、高等;
S3)比较步骤S1)中被测电池的电压-时间实验数据与步骤S2)中三维电化学-热耦合模型的电压-时间仿真数据对三维电化学-热耦合模型进行修正得到修正后的三维电化学-热耦合模型;
S4)利用修正后的三维电化学-热耦合模型,设置被测电池的初始电压、充放电停止条件与计算时间,采用二分法进行极限放电仿真;其中,被测电池的初始电压为23度SOC-OCV实验测试得到的OCV值;放电的停止条件为被测电池的电压小于被测电池下限截止电压3v;实验要求极限放电持续时间t1’为10s,记计算停止时实际计算时间为t2’,放电时满足t2’=t1’的电流I’为Imin’、满足t2’<t1’的电流I’为Imax’,预设的停止标准Imax’-Imin’≤0.5。
以25度50%SOC 10s极限放电电流为例进行具体说明,设置被测电池的初始电压为23度SOC-OCV实验测试得到的50%SOC状态的OCV值3.635V,首先通过尝试找到一个Imin与Imax,分别为470A与502A;然后I在[470、502]之间取值,把[470、502]均分为4份得到5个I值:470、478、486、494、502,依次对这5个I值进行仿真计算,得到停止计算时实际计算时间t2’分别为10s、10s、10s、9s、8s,判断得出新的Imin与Imax分别为486A、494A;然后I在[486、494]之间取值,把[486、494]均分为4份得到5个I值:486、488、490、492、494,依次对这5个I值进行仿真计算,得到停止计算时实际计算时间t2’分别为10s、10s、9s、9s、9s,判断得出新的Imin与Imax分别为488A、490A;然后I在[488、490]之间取值,把[488、490]均分为4份得到5个I值:488、488.5、489、489.5、490,依次对这5个I值进行仿真计算,得到停止计算时实际计算时间t2’分别为10s、10s、10s、9s、9s,判断得出新的Imin与Imax分别为489A、489.5A,此时Imax-Imin=0.5,满足预设的停止标准,停止计算。得到待测电池在25度50%SOC状态下的10s极限放电电流为489A。按照这种方法,依次计算得到待测电池在23度其它SOC状态下的10s极限放电电流。
图2为实施例23度10s放电的极限电流仿真结果与实验结果的比较图,两者结果基本一致,因此,通过仿真测试极限充放电电流,不仅缩短测试时间,而且节约实验成本;另外,本发明提供的方法既可用于充电极限电流仿真,也可用于放电极限电流仿真。
Claims (3)
1.一种采用二分法仿真锂离子电池极限充放电电流的方法,其特征在于:包括如下步骤:
S1)获取被测电池常温倍率放电的电压-时间实验数据,并在恒温箱中对被测电池进行不同温度下的SOC-OCV实验,获取被测电池在不同温度下的SOC-OCV实验数据;
S2)根据被测电池的参数建立被测电池的三维电化学-热耦合模型,并进行常温倍率放电仿真获取常温倍率放电的电压-时间仿真数据;
S3)比较步骤S1)中被测电池的电压-时间实验数据与步骤S2)中三维电化学-热耦合模型的电压-时间仿真数据对三维电化学-热耦合模型进行修正得到修正后的三维电化学-热耦合模型;
S4)利用修正后的三维电化学-热耦合模型,设置被测电池的初始电压、充放电停止条件与计算时间,采用二分法进行极限充放电仿真;
其中,被测电池的初始电压为被测电池的开路电压,为步骤S1)SOC-OCV实验数据中测的开路电压OCV值;充电的停止条件为被测电池的电压大于被测电池上限截止电压,放电的停止条件为被测电池的电压小于被测电池下限截止电压;计算时间为实验要求极限充电持续时间t1和实验要求极限放电持续时间t1’,记充电仿真计算停止时实际计算时间为t2和放电仿真计算停止时实际计算时间为t2’,充电时满足t2=t1的电流I为Imin、满足t2<t1的电流I为Imax,预设的停止标准Imax-Imin<a,充电电流I在[Imin,Imax]之间取值,放电时满足t2’=t1’的电流I’为Imin’、满足t2’<t1’的电流I’为Imax’,预设的停止标准Imax’-Imin’<a,放电电流I’在[Imin’,Imax’]之间取值。
2.根据权利要求1所述采用二分法仿真锂离子电池极限充放电电流的方法,其特征在于:所述步骤S4)中,利用修正后的三维电化学-热耦合模型采用二分法进行极限充电仿真具体过程为:首先预设一个I0进行充电仿真计算,根据计算停止时实际计算时间t2与t1的大小关系,来减小或增大I0,若t2=t1,则增大I0,直到t2<t1时停止;若t2<t1,则减小I0直到t2=t1时停止,这些I0值中,最后两个中的较小者为Imin、较大者为Imax;然后将[Imin,Imax]均分为n份获得n+1个电流值,对这n+1个电流值依次进行充电仿真,根据每个电流值充电仿真持续时间t2与t1的大小关系,判断得出新的Imax-1与Imin-1,n+1个I值所对应的t2中最后一个满足t2=t1的I为新的Imin即Imin-1、第一个满足t2<t1的I为新的Imax即Imax-1;依次重复进行,然后继续将[Imin-1,Imax-1]均分为n份获得n+1个电流值,对这n+1个电流值依次进行充电仿真,根据每个电流值充电仿真持续时间t2与t1的大小关系,判断得出新的Imax-2与Imin-2,直至新的Imax-i与Imin-i满足Imax-i-Imin-i<a停止计算,最后一次得到的Imin-i即为所需状态下充电的极限电流值,n取值为大于等于2的整数;
放电的极限电流值的仿真方法与充电的极限电流值的仿真方法相同。
3.根据权利要求1所述采用二分法仿真锂离子电池极限充放电电流的方法,其特征在于:所述步骤S2)中,参数包括设计参数、电极动力学参数、热物性参数等,设计参数包括被测电池的长、宽、高等。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910517159.3A CN110208708A (zh) | 2019-06-14 | 2019-06-14 | 一种采用二分法仿真锂离子电池极限充放电电流的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910517159.3A CN110208708A (zh) | 2019-06-14 | 2019-06-14 | 一种采用二分法仿真锂离子电池极限充放电电流的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110208708A true CN110208708A (zh) | 2019-09-06 |
Family
ID=67792832
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910517159.3A Pending CN110208708A (zh) | 2019-06-14 | 2019-06-14 | 一种采用二分法仿真锂离子电池极限充放电电流的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110208708A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110728056A (zh) * | 2019-10-14 | 2020-01-24 | 合肥国轩高科动力能源有限公司 | 一种锂离子电池充放电的最大电流仿真测试方法 |
CN114062946A (zh) * | 2021-10-21 | 2022-02-18 | 合肥国轩高科动力能源有限公司 | 一种锂离子电池操作极限电流的测试方法及系统 |
CN116070466A (zh) * | 2023-03-08 | 2023-05-05 | 上海泰矽微电子有限公司 | 一种电池截止电压下的最优soc仿真寻值方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103935260A (zh) * | 2014-05-08 | 2014-07-23 | 山东大学 | 一种基于电池安全保护的电池管理方法 |
CN106772081A (zh) * | 2016-12-22 | 2017-05-31 | 清华大学 | 基于扩展等效电路模型的电池极限充放电电流估计方法 |
US20190123574A1 (en) * | 2017-10-24 | 2019-04-25 | Samsung Electronics Co., Ltd. | Method and apparatus for charging battery |
CN109713762A (zh) * | 2019-02-28 | 2019-05-03 | 深圳猛犸电动科技有限公司 | 一种充电控制方法、装置及终端设备 |
-
2019
- 2019-06-14 CN CN201910517159.3A patent/CN110208708A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103935260A (zh) * | 2014-05-08 | 2014-07-23 | 山东大学 | 一种基于电池安全保护的电池管理方法 |
CN106772081A (zh) * | 2016-12-22 | 2017-05-31 | 清华大学 | 基于扩展等效电路模型的电池极限充放电电流估计方法 |
US20190123574A1 (en) * | 2017-10-24 | 2019-04-25 | Samsung Electronics Co., Ltd. | Method and apparatus for charging battery |
CN109713762A (zh) * | 2019-02-28 | 2019-05-03 | 深圳猛犸电动科技有限公司 | 一种充电控制方法、装置及终端设备 |
Non-Patent Citations (1)
Title |
---|
张彩萍 等: "动力电池组峰值功率估计算法研究", 《系统仿真学报》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110728056A (zh) * | 2019-10-14 | 2020-01-24 | 合肥国轩高科动力能源有限公司 | 一种锂离子电池充放电的最大电流仿真测试方法 |
CN114062946A (zh) * | 2021-10-21 | 2022-02-18 | 合肥国轩高科动力能源有限公司 | 一种锂离子电池操作极限电流的测试方法及系统 |
CN116070466A (zh) * | 2023-03-08 | 2023-05-05 | 上海泰矽微电子有限公司 | 一种电池截止电压下的最优soc仿真寻值方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ouyang et al. | A dynamic capacity degradation model and its applications considering varying load for a large format Li-ion battery | |
CN105633472B (zh) | 一种锂离子电池自放电率一致性配组筛选方法 | |
Lam et al. | Practical capacity fading model for Li-ion battery cells in electric vehicles | |
Stroe et al. | Accelerated aging of Lithium-ion batteries based on electric vehicle mission profile | |
CN107861075B (zh) | 一种确定动力电池sop的方法 | |
CN108767909A (zh) | 一种标准的充电曲线及充电方法 | |
CN110208708A (zh) | 一种采用二分法仿真锂离子电池极限充放电电流的方法 | |
CN110244234A (zh) | 一种电池加速寿命测试方法 | |
CN105929336B (zh) | 一种动力锂离子电池健康状态估算方法 | |
Thirugnanam et al. | Mathematical modeling of Li-ion battery for charge/discharge rate and capacity fading characteristics using genetic algorithm approach | |
CN103698716B (zh) | 一种基于衰减系数的串联电池组可放出电量衰减评价方法 | |
CN112464571B (zh) | 基于多约束条件粒子群优化算法的锂电池组参数辨识方法 | |
CN105911471B (zh) | 一种计算电动汽车电池的健康状态参数的装置和方法 | |
CN105912799A (zh) | 一种液态或半液态金属电池的建模方法 | |
CN109581234B (zh) | 一种锂离子电池一致性筛选方法 | |
CN101458310A (zh) | 一种电池一致性的评价方法 | |
Aizpuru et al. | Passive balancing design for Li-ion battery packs based on single cell experimental tests for a CCCV charging mode | |
CN104051810A (zh) | 一种锂离子储能电池系统soc估算快速修正方法 | |
CN104600381B (zh) | 一种锂离子电池组单体布置结构的优化方法 | |
CN106443483A (zh) | 一种电动汽车用镍氢电池的脉冲循环寿命测试方法 | |
CN111438077A (zh) | 一种退役三元软包电池的梯次利用性能快速筛选及检测的方法 | |
CN105703022A (zh) | 一种基于温度控制电池衰减的锂离子动力电池充电方法 | |
CN103018683A (zh) | 一种电池循环性能加速评估方法 | |
CN110729520A (zh) | 一种电池快速充电方法 | |
CN110308399A (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 | ||
AD01 | Patent right deemed abandoned | ||
AD01 | Patent right deemed abandoned |
Effective date of abandoning: 20211217 |