CN115048882A - 一种激光焊缝形貌预测方法 - Google Patents

一种激光焊缝形貌预测方法 Download PDF

Info

Publication number
CN115048882A
CN115048882A CN202210819398.6A CN202210819398A CN115048882A CN 115048882 A CN115048882 A CN 115048882A CN 202210819398 A CN202210819398 A CN 202210819398A CN 115048882 A CN115048882 A CN 115048882A
Authority
CN
China
Prior art keywords
weld
welding
next moment
moment
initial
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.)
Withdrawn
Application number
CN202210819398.6A
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.)
Changshu Institute of Technology
Original Assignee
Changshu Institute of Technology
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 Changshu Institute of Technology filed Critical Changshu Institute of Technology
Priority to CN202210819398.6A priority Critical patent/CN115048882A/zh
Publication of CN115048882A publication Critical patent/CN115048882A/zh
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23K—SOLDERING OR UNSOLDERING; WELDING; CLADDING OR PLATING BY SOLDERING OR WELDING; CUTTING BY APPLYING HEAT LOCALLY, e.g. FLAME CUTTING; WORKING BY LASER BEAM
    • B23K26/00—Working by laser beam, e.g. welding, cutting or boring
    • B23K26/20—Bonding
    • B23K26/21—Bonding by welding
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23K—SOLDERING OR UNSOLDERING; WELDING; CLADDING OR PLATING BY SOLDERING OR WELDING; CUTTING BY APPLYING HEAT LOCALLY, e.g. FLAME CUTTING; WORKING BY LASER BEAM
    • B23K26/00—Working by laser beam, e.g. welding, cutting or boring
    • B23K26/70—Auxiliary operations or equipment
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Mechanical Engineering (AREA)
  • Plasma & Fusion (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Fluid Mechanics (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了一种激光焊缝形貌预测方法，针对已有激光焊缝形貌预测方法中存在的数据和场景依赖性问题，本发明从经典物理学理论出发，提出了一种普适的激光焊缝形貌预测方法，从而实现对不同焊接设备、焊接材料及焊接结构下的焊缝形貌的准确预测。首先，输出初始时刻的焊缝形貌函数和焊接参数；其次，基于经典物理学理论求解下一时刻的焊缝形貌函数；进而，提取下一时刻的焊缝形貌并输出下一时刻的焊缝形貌特征。与现有焊缝形貌预测方法相比，本发明具有较好的普适性，算法简单容易实施，计算效率高，数值稳定性和鲁棒性较好。

Description

一种激光焊缝形貌预测方法
技术领域
本发明涉及一种普适的激光焊缝形貌预测方法，其属于焊缝形貌预测方法领域。
背景技术
焊缝形貌是评价焊接工艺质量的关键指标。为获得良好的焊缝形貌，科研人员和工程师需进行大量的焊接工艺实验。通过对比不同工艺条件下的焊缝形貌，筛选出最优的焊缝形貌。由于焊接工艺实验需要消耗较多的人力和物力成本，此方法存在一定的不足。数值模拟方法可以预测不同工艺条件下的的焊缝形貌，减少焊接工艺实验次数，降低焊接实验成本，对于提高激光焊接研发和生产效率有着重要的意义。
在激光焊接过程中，材料受高能量密度激光辐照作用会产生剧烈的蒸发，该物理过程会强烈改变激光焊接的焊缝形貌。针对激光焊接的蒸发过程，现有数值模拟方法主要采用反冲压力模型进行描述(Pang S,Chen L,Zhou J,et al.Journal of Physics D:Applied Physics,2011,44:025301；Ai Y,Jiang P,Shao X,et al.Applied ThermalEngineering,2016,113:980-993；Feng Y,Gao X,Zhang Y,et al.The InternationalJournal of Advanced Manufacturing Technology,2021,112:2301-2312.)。此类方法基于经验模型描述蒸发过程给熔池带来的动量和能量作用，存在一定的数据和场景依赖性。当焊接设备、焊接材料或焊接结构改变后，此方法很难准确预测焊缝形貌。上述难点制约了焊缝形貌预测方法在激光焊接中的应用，因此迫切需要发展一种普适的激光焊缝形貌预测方法。
发明内容
为解决已有激光焊缝形貌预测方法存在的数据和场景依赖性问题，本发明提出了一种普适的激光焊缝形貌预测方法。
本发明采用如下技术方案：一种普适的激光焊缝形貌预测方法，包括如下步骤：
S1,输入初始时刻的焊缝形貌函数Hn：
Figure BDA0003743566350000011
其中，公式(1)中，Δxmin为最小网格尺寸，φn为初始时刻工件截面内各点到焊缝的距离，tanh为双曲正切函数；
S2,确定初始时刻的焊接参数，各焊接参数表达如下：
Figure BDA0003743566350000012
Figure BDA0003743566350000021
S3,获取下一时刻的焊缝形貌函数Hn+1，获取过程如下：
S3.1,计算初始时刻的动量源Sn：
Figure BDA0003743566350000022
其中，
Figure BDA0003743566350000023
分别代表由惯性、表面力和体积力引起的动量源，
Figure BDA0003743566350000024
为梯度算子，
Figure BDA0003743566350000025
为散度算子，上标T代表矩阵转置，un、Tn、α、Tenv、σ为步骤S2中列举的初始时刻的焊接参数，g为重力加速度常数，ρ、μ分别为气液混合相的密度和粘度，n为气液界面法向量；
S3.2,计算对流后的速度u*：
Figure BDA0003743566350000026
其中，Δt为计算时间步长，un为步骤S2中列举的初始时刻的焊接参数；
S3.3,计算下一时刻的压力pn+1：
Figure BDA0003743566350000027
其中，
Figure BDA0003743566350000028
为材料蒸发时的质量流率；
S3.4,计算下一时刻的速度un+1：
Figure BDA0003743566350000031
S3.5,计算初始时刻的能量源Qn：
Figure BDA0003743566350000032
其中，e为自然对数，Tn、Lvap、Cliq、Cvap、Tsat、I、r0为步骤S2中列举的初始时刻的焊接参数，r为计算域内各点到激光光斑中心的距离，C和k分别为气液混合相的比热容和热导率；
S3.6,计算下一时刻的温度Tn+1：
Figure BDA0003743566350000033
S3.7,计算下一时刻的焊缝形貌函数Hn+1：
Figure BDA0003743566350000034
S4,绘制下一时刻焊缝形貌函数Hn+1的分布云图；
S5,提取下一时刻的焊缝形貌；
S6,输出下一时刻的焊缝形貌特征。
与现有方法相比，本发明的优点如下：
(1)本方法基于经典物理学理论求解激光焊缝形貌演化过程，有效解决已有经验性方法中的数据和场景依赖性问题，可适用于不同的焊接设备、材料及结构，具有较好的普适性；
(2)本方法基于不可压缩算法求解激光焊缝形貌演化过程，具有方法简单、容易实施、计算效率高等优点；
(3)本方法采用一套控制方程求解液相和气相两种流体，具有较好的数值稳定性和鲁棒性。
附图说明
图1是本发明提出的一种普适的激光焊缝形貌预测方法的流程过程图。
图2是初始时刻的焊缝形貌函数定义的示意图。
图3是下一时刻焊缝形貌函数分布云图的示意图。
图4是下一时刻的焊缝形貌示意图。
图5是下一时刻的焊缝形貌特征示意图。
具体实施方式
下面结合附图和实施例对本发明进行进一步阐述。
针对已有激光焊缝形貌预测方法存在的数据和场景依赖性问题，本发明从经典物理学理论出发，提出了一种普适的激光焊缝形貌预测方法。
结合图1，本发明所述的一种普适的激光焊缝形貌预测方法，包括如下步骤：
第一步：输入初始时刻的焊缝形貌函数Hn
本发明以双曲正切函数表示初始时刻的焊缝形貌Hn，该初始时刻的焊缝形貌定义为双曲正切函数的0.5等值面。因此，初始时刻的焊缝形貌函数Hn见公式(1)：
Figure BDA0003743566350000041
其中，Δxmin为最小网格尺寸，φn为初始时刻工件截面内各点到焊缝的距离，tanh为双曲正切函数，如图2所示。
第二步：确定初始时刻的焊接参数
初始时刻的焊接参数，见表1：
表1初始时刻的焊接参数
Figure BDA0003743566350000042
第三步：获取下一时刻的焊缝形貌函数Hn+1
本发明基于经典物理学理论求解下一时刻的焊缝形貌函数Hn+1，具体实现过程如下：
①计算初始时刻的动量源Sn，见公式(2)：
Figure BDA0003743566350000051
其中，
Figure BDA0003743566350000052
分别代表由惯性、表面力和体积力引起的动量源，
Figure BDA0003743566350000053
为梯度算子，
Figure BDA0003743566350000054
为散度算子，上标T代表矩阵转置，un、Tn、α、Tenv、σ为初始时刻的焊接参数，由表1给出，g为重力加速度常数，ρ、μ分别为气液混合相的密度和粘度，见公式(2-1)：
Figure BDA0003743566350000055
其中，ρliq、ρvap、μliq、μvap为初始时刻的焊接参数，均由表1给出。
n为气液界面法向量，见公式(2-2)：
Figure BDA0003743566350000056
其中，Hn为初始时刻的焊缝形貌，由公式(1)获得。
②计算对流后的速度u*，见公式(3)：
Figure BDA0003743566350000057
其中，
Figure BDA00037435663500000512
为梯度算子，Δt为计算时间步长，un为初始时刻的焊接参数，由表1给出，Sn由公式(2)获得。
③计算下一时刻的压力pn+1，见公式(4)：
Figure BDA0003743566350000058
其中，
Figure BDA0003743566350000059
为梯度算子，
Figure BDA00037435663500000510
为散度算子，Δt为计算时间步长，ρ为气液混合相的密度，由公式(2-1)获得，u*为对流后的速度，由公式(3)获得，Hn为初始时刻的焊缝形貌，由公式(1)获得，
Figure BDA00037435663500000513
为材料蒸发时的质量流率，见公式(4-1)：
Figure BDA00037435663500000511
其中，
Figure BDA00037435663500000514
为梯度算子，Tn、kvap、kliq、Lvap为初始时刻的焊接参数，由表1给出，n为气液界面法向量，由公式(2-2)获得。
④计算下一时刻的速度un+1，见公式(5)：
Figure BDA0003743566350000061
其中，
Figure BDA0003743566350000062
为梯度算子，Δt为计算时间步长，u*为对流后的速度，由公式(3)获得，pn +1为下一时刻的压力，由公式(4)获得，ρ由公式(2-1)获得。
⑤计算初始时刻的能量源Qn，见公式(6)：
Figure BDA0003743566350000063
其中，
Figure BDA0003743566350000064
为梯度算子，
Figure BDA0003743566350000065
为散度算子，e为自然对数，Tn、Lvap、Cliq、Cvap、Tsat、I、r0为初始时刻的焊接参数，由表1给出，ρ为气液混合相的密度，由公式(2-1)获得，
Figure BDA00037435663500000613
为材料蒸发时的质量流率，由公式(4-1)获得，Hn为初始时刻的焊缝形貌，由公式(1)获得，r为计算域内各点到激光光斑中心的距离，C和k分别为气液混合相的比热容和热导率：
Figure BDA0003743566350000066
其中，Cliq、Cvap、kliq、kvap为初始时刻的焊接参数，均由表1给出，Hn为初始时刻的焊缝形貌，由公式(1)获得。
⑥计算下一时刻的温度Tn+1，见公式(7)：
Figure BDA0003743566350000067
其中，
Figure BDA0003743566350000068
为梯度算子，Δt为计算时间步长，un、Tn为初始时刻的焊接参数，由表1给出，Qn为初始时刻的能量源，由公式(6)获得。
⑦计算下一时刻的焊缝形貌函数Hn+1，见公式(8)：
Figure BDA0003743566350000069
其中，
Figure BDA00037435663500000610
为梯度算子，Δt为计算时间步长，Hn为初始时刻的焊缝形貌，由公式(1)获得，un+1为下一时刻的速度，由公式(5)获得，ρ为气液混合相的密度，由公式(2-1)获得，
Figure BDA00037435663500000614
为材料蒸发时的质量流率，由公式(4-1)获得，n为气液界面法向量，由公式(2-2)获得。
其中，公式(2)至公式(8)中的梯度算子
Figure BDA00037435663500000611
和散度算子
Figure BDA00037435663500000612
采用有限体积法(FiniteVolume Method)离散，以保证严格的物理场守恒型。
公式(4)可采用超松弛法(Successive Over Relaxation Method)、共轭梯度法(Conjugate Gradient Method)或快速傅里叶变换(Fast Fourier Transform)等方法进行求解。
第四步：绘制下一时刻的焊缝形貌函数Hn+1的分布云图。
根据下一时刻的焊缝形貌函数Hn+1，绘制下一时刻的焊缝形貌函数Hn+1的分布云图，如图3所示。
第五步：提取下一时刻的焊缝形貌。
根据下一时刻的焊缝形貌函数Hn+1的分布云图，提取Hn+1的0.5等值面作为下一时刻的焊缝形貌，如图4所示。
第六步：输出下一时刻的焊缝形貌特征。
根据下一时刻的焊缝形貌，输出下一时刻的焊缝形貌特征，包括焊缝熔宽、熔深和余高等，如图5所示。
与现有焊缝形貌预测方法相比，本发明基于经典物理学理论求解激光焊缝形貌演化过程，有效解决已有经验性方法中存在的数据和场景依赖性问题，可适用于不同的焊接设备、材料及结构，具有较好的普适性。本方法基于不可压缩解法求解激光焊缝形貌演化过程，具有算法简单、容易实施、计算效率高等优点。本方法采用一套控制方程求解液相和气相两种流体，具有较好的数值稳定性和鲁棒性。

Claims (9)

1.一种激光焊缝形貌预测方法，其特征在于，包括如下步骤：
S1,输入初始时刻的焊缝形貌函数Hn：
Figure FDA0003743566340000011
其中，公式(1)中，Δxmin为最小网格尺寸，φn为初始时刻工件截面内各点到焊缝的距离，tanh为双曲正切函数；
S2,确定初始时刻的焊接参数，各焊接参数表达如下：
Figure FDA0003743566340000012
S3,获取下一时刻的焊缝形貌函数Hn+1，获取过程如下：
S3.1,计算初始时刻的动量源Sn：
Figure FDA0003743566340000013
其中，
Figure FDA0003743566340000014
分别代表由惯性、表面力和体积力引起的动量源，
Figure FDA0003743566340000015
为梯度算子，
Figure FDA0003743566340000016
为散度算子，上标T代表矩阵转置，un、Tn、α、Tenv、σ为步骤S2中列举的初始时刻的焊接参数，g为重力加速度常数，ρ、μ分别为气液混合相的密度和粘度，n为气液界面法向量；
S3.2,计算对流后的速度u*：
Figure FDA0003743566340000021
其中，Δt为计算时间步长，un为步骤S2中列举的初始时刻的焊接参数；
S3.3,计算下一时刻的压力pn+1：
Figure FDA0003743566340000022
其中，
Figure FDA0003743566340000027
为材料蒸发时的质量流率；
S3.4,计算下一时刻的速度un+1：
Figure FDA0003743566340000023
S3.5,计算初始时刻的能量源Qn：
Figure FDA0003743566340000024
其中，e为自然对数，Tn、Lvap、Cliq、Cvap、Tsat、I、r0为步骤S2中列举的初始时刻的焊接参数，r为计算域内各点到激光光斑中心的距离，C和k分别为气液混合相的比热容和热导率；
S3.6,计算下一时刻的温度Tn+1：
Figure FDA0003743566340000025
S3.7,计算下一时刻的焊缝形貌函数Hn+1：
Figure FDA0003743566340000026
S4,绘制下一时刻焊缝形貌函数Hn+1的分布云图；
S5,提取下一时刻的焊缝形貌；
S6,输出下一时刻的焊缝形貌特征。
2.如权利要求1所述的方法，其特征在于，步骤S3.1中，ρ、μ分别为气液混合相的密度和粘度，根据公式(2-1)获得：
Figure FDA0003743566340000031
其中，ρliq、ρvap、μliq、μvap为步骤S2中列举的初始时刻的焊接参数；
n为气液界面法向量，根据公式(2-2)获得：
Figure FDA0003743566340000032
3.如权利要求1所述的方法，其特征在于，步骤S3.3中，
Figure FDA0003743566340000036
为材料蒸发时的质量流率，根据公式(4-1)获得：
Figure FDA0003743566340000033
其中，Tn、kvap、kliq、Lvap为步骤S2中列举的初始时刻的焊接参数，n为气液界面法向量。
4.如权利要求1所述的方法，其特征在于，步骤S3.5中，C和k分别为气液混合相的比热容和热导率，根据公式(6-1)获得：
Figure FDA0003743566340000034
其中，Cliq、Cvap、kliq、kvap为步骤S2中列举的初始时刻的焊接参数。
5.如权利要求1所述的方法，其特征在于，公式(2)至公式(8)中的梯度算子
Figure FDA0003743566340000037
和散度算子
Figure FDA0003743566340000035
采用有限体积法离散。
6.如权利要求1所述的方法，其特征在于，公式(4)采用超松弛法、共轭梯度法或快速傅里叶变换中任一方法进行求解。
7.如权利要求1所述的方法，其特征在于，步骤S1中，以双曲正切函数表示初始时刻的焊缝形貌函数Hn，焊缝形貌定义为双曲正切函数的0.5等值面。
8.如权利要求1所述的方法，其特征在于，步骤S5中，根据下一时刻的焊缝形貌函数Hn +1，提取Hn+1的0.5等值面作为下一时刻的焊缝形貌。
9.如权利要求1所述的方法，其特征在于，步骤S6中，根据下一时刻的焊缝形貌，输出下一时刻的焊缝形貌特征，包括焊缝熔宽、熔深和余高。
CN202210819398.6A 2022-07-13 2022-07-13 一种激光焊缝形貌预测方法 Withdrawn CN115048882A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210819398.6A CN115048882A (zh) 2022-07-13 2022-07-13 一种激光焊缝形貌预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210819398.6A CN115048882A (zh) 2022-07-13 2022-07-13 一种激光焊缝形貌预测方法

Publications (1)

Publication Number Publication Date
CN115048882A true CN115048882A (zh) 2022-09-13

Family

ID=83164905

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210819398.6A Withdrawn CN115048882A (zh) 2022-07-13 2022-07-13 一种激光焊缝形貌预测方法

Country Status (1)

Country Link
CN (1) CN115048882A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114841002A (zh) * 2022-05-11 2022-08-02 常熟理工学院 一种适用于万瓦级激光焊接的高精度焊缝形貌预测方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114841002A (zh) * 2022-05-11 2022-08-02 常熟理工学院 一种适用于万瓦级激光焊接的高精度焊缝形貌预测方法
CN114841002B (zh) * 2022-05-11 2024-05-03 常熟理工学院 一种适用于万瓦级激光焊接的高精度焊缝形貌预测方法

Similar Documents

Publication Publication Date Title
Ji et al. Large eddy simulation of turbulent attached cavitating flow with special emphasis on large scale structures of the hydrofoil wake and turbulence-cavitation interactions
Di Ilio et al. Fluid flow around NACA 0012 airfoil at low-Reynolds numbers with hybrid lattice Boltzmann method
Ji et al. Numerical analysis of unsteady cavitating turbulent flow and shedding horse-shoe vortex structure around a twisted hydrofoil
Wen et al. Structure and stability of steady porous medium convection at large Rayleigh number
Shen et al. Direct numerical simulations of turbulence and hyporheic mixing near sediment–water interfaces
CN115048882A (zh) 一种激光焊缝形貌预测方法
Elsakka et al. Comparison of the computational fluid dynamics predictions of vertical axis wind turbine performance against detailed pressure measurements
Li et al. Very large eddy simulation of cavitation from inception to sheet/cloud regimes by a multiscale model
Anika et al. Hall current effects on magnetohydrodynamics fluid over an infinite rotating vertical porous plate embedded in unsteady laminar flow
Issakhov Numerical modelling of the thermal effects on the aquatic environment from the thermal power plant by using two water discharge pipes
Issakhov Numerical modelling of distribution the discharged heat water from thermal power plant on the aquatic environment
Yin et al. Large eddy simulation of cloud cavitation and wake vortex cavitation around a trailing-truncated hydrofoil
EP4150499A1 (en) Computational analysis of physical systems
Yang et al. Thermal and fluid transport in micro-open-cell metal foams: Effect of node size
Diao et al. Three-dimensional prediction of reservoir water temperature by the lattice Boltzmann method: Validation
Martin et al. Numerical investigation of slamming loads on floating offshore wind turbines
Chen et al. Construction and validation of an urban area flow and dispersion model on building scales
Athota et al. Simulations of thermal wind formation in idealised mountain-valley systems using OpenFOAM
Kanehira et al. The effects of smoothing length on the onset of wave breaking in smoothed particle hydrodynamics (SPH) simulations of highly directionally spread waves
Hidalgo et al. Cavitating flow simulation with mesh development using Salome open source software
Jiang et al. Numerical simulation of three-dimensional high-speed flows using a second-order nonlinear model
Lu et al. A robust scheme for numerical simulation of heat transfer in two-fluid flows with high volumetric heat capacity contrasts
CN114841002B (zh) 一种适用于万瓦级激光焊接的高精度焊缝形貌预测方法
Gordeev et al. Investigation on development of free-surface structures in turbulent liquid lithium flow
Zhang et al. Research Status on Inflow turbulence generation method with Large Eddy Simulation of CFD numerical wind tunnel

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
WW01 Invention patent application withdrawn after publication
WW01 Invention patent application withdrawn after publication

Application publication date: 20220913