CN111339699A - 基于宏观应力与残余应力拓度的箱型结构重变形量预测方法 - Google Patents

基于宏观应力与残余应力拓度的箱型结构重变形量预测方法 Download PDF

Info

Publication number
CN111339699A
CN111339699A CN202010086804.3A CN202010086804A CN111339699A CN 111339699 A CN111339699 A CN 111339699A CN 202010086804 A CN202010086804 A CN 202010086804A CN 111339699 A CN111339699 A CN 111339699A
Authority
CN
China
Prior art keywords
stress
residual stress
box
macroscopic
point
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
CN202010086804.3A
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.)
Yangzhou University
Original Assignee
Yangzhou 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 Yangzhou University filed Critical Yangzhou University
Priority to CN202010086804.3A priority Critical patent/CN111339699A/zh
Publication of CN111339699A publication Critical patent/CN111339699A/zh
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L1/00Measuring force or stress, in general
    • G01L1/25Measuring force or stress, in general using wave or particle radiation, e.g. X-rays, microwaves, neutrons
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L5/00Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
    • G01L5/0047Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes measuring forces due to residual stresses

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Toxicology (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了基于宏观应力与残余应力拓度的箱型结构重变形量预测方法,该方法包括如下步骤:S1.箱型结构焊缝区域残余应力检测;S2.工况条件下结构宏观应力畸变区域系数的计算;S3.基于畸变区域系数的焊缝区域残余应力修正;S4.考虑工况应力峰谷效应的残余应力拓度计算;S5.基于宏观应力与残余应力拓度的箱型结构重变形量计算。该方法检测精度高,对于实现箱型结构重变形分布的实时评估具有重要的现实意义。

Description

基于宏观应力与残余应力拓度的箱型结构重变形量预测方法
技术领域
本发明涉及箱型结构重变形量检测领域,特别涉及基于宏观应力与残余应力拓度的箱型结构重变形量预测方法。
背景技术
随着社会的进步,箱型结构被越来越多的应用到大型工程机械结构中,其在使用过程中的二次变形问题越来越受到社会的关注。在这样的背景前提下,对箱型结构的重变形量进行实时预测是必要且紧迫的。大多箱型结构的检测方法往往单纯基于宏观应力或残余应力展开,存在计算准确性不高且检测效率低等缺点。如何将找出代表宏观应力与微观残余应力的关键因子,并在此基础上将二者结合来对箱型结构变形重分布进行预测是目前这一研究领域存在的共性问题。
发明内容
发明目的:本发明目的是提供基于宏观应力与残余应力拓度的箱型结构重变形量预测方法,该方法检测精度高。
技术方案:本发明提供基于宏观应力与残余应力拓度的箱型结构重变形量预测方法,包括如下步骤:
S1.箱型结构焊缝区域残余应力检测:
采用残余应力测试装置对箱型结构焊缝处每隔预设间隔取一个点i进行残余应力的测试,i=1,2,3...N,i为点的标号;获得某一点i对应XYZ三维坐标系中于X方向上的残余应力σxi、对应于Y方向上的残余应力σyi、对应于Z方向上的残余应力σzi
箱型结构焊缝处每个测量点i的残余应力总值σri表示为:
Figure BDA0002382348050000011
其中,σri为箱型结构焊缝处每个测量点的残余应力总值;βxi为箱型结构所对应的X向残余应力混合修正系数;βyi为箱型结构所对应的Y向残余应力混合修正系数;βzi为箱型结构所对应的Z向残余应力混合修正系数;Cmaxxiyizi)为βxi,βyi,βzi这三个参数中的最大值;i为取点的标号,i=1,2,3...N;
Figure BDA0002382348050000012
Figure BDA0002382348050000013
Figure BDA0002382348050000014
其中Cmaxxiyizi)为σxi,σyi,σzi这三个参数中的最大值;L为残余应力测量的点距;
得到箱型结构焊缝处每个测量点的残余应力总值σri后,对焊缝上的残余应力分布函数进行拟合获得F(p)=β·(A1P3+A2P2+A3P+A4)
其中,A1,A2,A3,A4为分布函数的拟合系数;F(p)为焊缝上残余应力分布函数;p为焊缝上任意一点距离焊缝初始点位置的距离;β为箱型结构所对应的名义残余应力混合修正系数;
Figure BDA0002382348050000021
其中
Figure BDA0002382348050000022
为i=1,2,3...N上βxi的叠加和;
Figure BDA0002382348050000023
为i=1,2,3...N上βyi的叠加和;
Figure BDA0002382348050000024
为i=1,2,3...N上βzi的叠加和;
S2.工况条件下结构宏观应力畸变区域系数的计算:
建立箱型结构的有限元分析模型,在有限元软件中对其进行网格划分,并按照实行工况条件输入载荷与边界约束,最后进行求解以确定箱型结构焊缝处的应力分布及最容易出现损伤的区域位置,并在有限元软件中自动生成该区域面积为A;
在实际工况条件下工作的箱型结构中,对上述通过有限元分析手段确定的最容易出现损伤的区域上选择n个均匀分布的点,并对这n个点在一段时间0到tm时间内所对应的宏观应力σqt进行测量,q为测量点的标记,q=1,2...n,t为测量的对应时间;
计算最容易出现损伤的区域中每个测量点宏观应力的最大间隔值:
Figure BDA0002382348050000025
其中,σJq为最容易出现损伤的区域中每个测量点宏观应力的最大间隔值;Cmaxqt)为第q个测量点在0到tm时刻内所有监测宏观应力中的最大值;Cminqt)为第q个测量点在0到tm时刻内所有监测宏观应力中的最小值;tmax为Cmaxqt)所对应的时间点;tmin为Cminqt)所对应的时间点;σm为数据监测的总时间;
Figure BDA0002382348050000026
为第q个测量点的宏观应力平均值;
Figure BDA0002382348050000027
然后计算工况应力条件下结构宏观应力的畸变区域系数B
Figure BDA0002382348050000031
其中,A为有限元软件中箱型结构焊缝处的应力分布及最容易出现损伤的区域的面积;C(σJq)max为n个测量点宏观应力最大间隔值中的最大值;C(σJq)min为n个测量点宏观应力最大间隔值中的最小值;y为最容易出现损伤的区域的面积变量;
S3.基于畸变区域系数的焊缝区域残余应力修正:
在上述过程的基础上,基于畸变系数对焊缝区域残余应力进行修正;
计算修正后的残余应力值σc
Figure BDA0002382348050000032
其中,B为工况应力条件下结构宏观应力的畸变区域系数;i为取点的标号,i=1,2,3...N;N为残余应力测量时取点的最大值;σri为箱型结构焊缝处每测量点的残余应力总值;β为箱型结构所对应的名义残余应力混合修正系数;
S4.考虑工况应力峰谷效应的残余应力拓度计算:
在上述过程的基础上,计算考虑工况应力峰谷效应的残余应力拓度D
Figure BDA0002382348050000033
其中,D为虑工况应力峰谷效应的残余应力拓度;
Figure BDA0002382348050000036
为为第q个测量点的宏观应力平均值;Cmaxqt)为第q个测量点在0到tm时刻内所有监测宏观应力中的最大值;Cminqt)为第q个测量点在0到tm时刻内所有监测宏观应力中的最小值;σc修正后的残余应力值;β为箱型结构所对应的名义残余应力混合修正系数;
S5.基于宏观应力与残余应力拓度的箱型结构重变形量计算:
在上述过程的基础上,基于宏观应力与残余应力拓度计算最容易出现损伤的区域中每个测量点的结构重变形值εq,q为测量点的标记,q=1,2,...n;
Figure BDA0002382348050000034
其中,
Figure BDA0002382348050000035
为为第q个测量点的宏观应力平均值;D为考虑工况应力峰谷效应的残余应力拓度;σc为修正后的残余应力值;E为弹性模量。
有益效果:本发明可以实现箱型结构重变形量的检测,并从基于宏观应力与残余应力拓度的角度出发来对箱型结构重变形量进行检测,从而避免了单纯依靠宏观应力或残余应力来检测整体箱型结构重变形量的局限性,更有效地提高箱型结构重变形量检测的准确性与效率。
附图说明
图1为本发明提供的基于宏观应力与残余应力拓度的箱型结构重变形量预测方法流程框图。
具体实施方式
如图1所示,本实施例基于宏观应力与残余应力拓度的箱型结构重变形量预测方法,该方法包括以下步骤:
S1.箱型结构焊缝区域残余应力检测:
采用X射线残余应力测试装置对箱型结构焊缝处每隔5mm取一个点i(i=1,2,3...N,i为点的标号)进行残余应力的测试,获得某一点i对应于X方向上的残余应力σx,对应于Y方向上的残余应力σy,对应于Z方向上的残余应力σz
箱型结构焊缝处每个测量点的残余应力总值σri可表示为:
Figure BDA0002382348050000041
其中,σri为箱型结构焊缝处每个测量点的残余应力总值;σxi为箱型结构焊缝处每个测量点在X方向上的残余应力;σyi为箱型结构焊缝处每个测量点在Y方向上的残余应力;σzi为箱型结构焊缝处每个测量点在Z方向上的残余应力;βxi为箱型结构所对应的X向残余应力混合修正系数;βyi为箱型结构所对应的Y向残余应力混合修正系数;βzi为箱型结构所对应的Z向残余应力混合修正系数;Cmaxxiyizi)为βxi,βyi,βzi这三个参数中的最大值;i为取点的标号(i=1,2,3...N)。
Figure BDA0002382348050000042
Figure BDA0002382348050000043
Figure BDA0002382348050000044
其中Cmaxxiyizi)为σxi,σyi,σzi这三个参数中的最大值;L为X射线残余应力测量的点距。
计算得到箱型结构焊缝处每个测量点的残余应力总值σri后,对焊缝上的残余应力分布函数进行拟合获得F(p)=β·(A1P3+A2P2+A3P+A4) (5)
其中,A1,A2,A3,A4为分布函数的拟合系数;F(p)为焊缝上残余应力分布函数;p为焊缝上任意一点距离焊缝初始点位置的距离;β为箱型结构所对应的名义残余应力混合修正系数。
Figure BDA0002382348050000045
其中
Figure BDA0002382348050000051
为i=1,2,3...N上βxi的叠加和;
Figure BDA0002382348050000052
为i=1,2,3...N上βyi的叠加和;
Figure BDA0002382348050000053
为i=1,2,3...N上βzi的叠加和。
S2.工况条件下结构宏观应力畸变区域系数的计算:
建立箱型结构的有限元分析模型,在有限元软件中对其进行网格划分,并按照实行工况条件输入载荷与边界约束,最后进行求解以确定箱型结构焊缝处的应力分布及最容易出现损伤的区域位置,并在有限元软件中自动生成该区域面积为A。
在实际工况条件下工作的箱型结构中,对上述通过有限元分析手段确定的最容易出现损伤的区域上选择7个均匀分布的点,并对这7个点在一段时间0到tm时间内所对应的宏观应力σqt进行测量(q为测量点的标记,q=1,2...7,t为测量的对应时间)。
计算最容易出现损伤的区域中每个测量点宏观应力的最大间隔值:
Figure BDA0002382348050000054
其中,σJq为最容易出现损伤的区域中每个测量点宏观应力的最大间隔值;Cmaxqt)为第q个测量点在0到tm时刻内所有监测宏观应力中的最大值;Cminqt)为第q个测量点在0到tm时刻内所有监测宏观应力中的最小值;tmax为Cmaxqt)所对应的时间点;tmin为Cminqt)所对应的时间点;σm为数据监测的总时间;
Figure BDA0002382348050000055
为第q个测量点的宏观应力平均值。
Figure BDA0002382348050000056
然后计算工况应力条件下结构宏观应力的畸变区域系数B
Figure BDA0002382348050000057
其中,A为有限元软件中箱型结构焊缝处的应力分布及最容易出现损伤的区域的面积;C(σJq)max为7个测量点宏观应力最大间隔值中的最大值;C(σJq)min为7个测量点宏观应力最大间隔值中的最小值;y为面积变量。
S3.基于畸变区域系数的焊缝区域残余应力修正:
在上述过程的基础上,基于畸变系数对焊缝区域残余应力进行修正。
计算修正后的残余应力值σc
Figure BDA0002382348050000061
其中,B为工况应力条件下结构宏观应力的畸变区域系数;i为取点的标号(i=1,2,3...N);N为残余应力测量时取点的最大值;σri为箱型结构焊缝处每测量点的残余应力总值;β为箱型结构所对应的名义残余应力混合修正系数。
S4.考虑工况应力峰谷效应的残余应力拓度计算:
在上述过程的基础上,计算考虑工况应力峰谷效应的残余应力拓度D
Figure BDA0002382348050000062
其中,D为虑工况应力峰谷效应的残余应力拓度;
Figure BDA0002382348050000063
为为第q个测量点的宏观应力平均值;Cmaxqt)为第q个测量点在0到tm时刻内所有监测宏观应力中的最大值;Cminqt)为第q个测量点在0到tm时刻内所有监测宏观应力中的最小值;σc修正后的残余应力值;β为箱型结构所对应的名义残余应力混合修正系数。
S5.基于宏观应力与残余应力拓度的箱型结构重变形量计算:
在上述过程的基础上,基于宏观应力与残余应力拓度计算最容易出现损伤的区域中每个测量点的结构重变形值εq(q为测量点的标记,q=1,2,...7)
Figure BDA0002382348050000064
其中,
Figure BDA0002382348050000065
为为第q个测量点的宏观应力平均值;D为考虑工况应力峰谷效应的残余应力拓度;σc为修正后的残余应力值;E为弹性模量,取值为2.06×1011Pa。

Claims (5)

1.基于宏观应力与残余应力拓度的箱型结构重变形量预测方法,其特征在于:包括如下步骤:
S1.箱型结构焊缝区域残余应力检测:
采用残余应力测试装置对箱型结构焊缝处每隔预设间隔取一个点i进行残余应力的测试,i=1,2,3...N,i为点的标号;获得某一点i对应XYZ三维坐标系中于X方向上的残余应力σxi、对应于Y方向上的残余应力σyi、对应于Z方向上的残余应力σzi
箱型结构焊缝处每个测量点i的残余应力总值σri表示为:
Figure FDA0002382348040000011
其中,σri为箱型结构焊缝处每个测量点的残余应力总值;βxi为箱型结构所对应的X向残余应力混合修正系数;βyi为箱型结构所对应的Y向残余应力混合修正系数;βzi为箱型结构所对应的Z向残余应力混合修正系数;Cmaxxiyizi)为βxi,βyi,βzi这三个参数中的最大值;i为取点的标号,i=1,2,3...N;
Figure FDA0002382348040000012
Figure FDA0002382348040000013
Figure FDA0002382348040000014
其中Cmaxxiyizi)为σxi,σyi,σzi这三个参数中的最大值;L为残余应力测量的点距;
得到箱型结构焊缝处每个测量点的残余应力总值σri后,对焊缝上的残余应力分布函数进行拟合获得F(p)=β·(A1P3+A2P2+A3P+A4)
其中,A1,A2,A3,A4为分布函数的拟合系数;F(p)为焊缝上残余应力分布函数;p为焊缝上任意一点距离焊缝初始点位置的距离;β为箱型结构所对应的名义残余应力混合修正系数;
Figure FDA0002382348040000015
其中
Figure FDA0002382348040000016
为i=1,2,3...N上βxi的叠加和;
Figure FDA0002382348040000017
为i=1,2,3...N上βyi的叠加和;
Figure FDA0002382348040000018
为i=1,2,3...N上βzi的叠加和;
S2.工况条件下结构宏观应力畸变区域系数的计算:
建立箱型结构的有限元分析模型,在有限元软件中对其进行网格划分,并按照实行工况条件输入载荷与边界约束,最后进行求解以确定箱型结构焊缝处的应力分布及最容易出现损伤的区域位置,并在有限元软件中自动生成该区域面积为A;
在实际工况条件下工作的箱型结构中,对上述通过有限元分析手段确定的最容易出现损伤的区域上选择n个均匀分布的点,并对这n个点在一段时间0到tm时间内所对应的宏观应力σqt进行测量,q为测量点的标记,q=1,2...n,t为测量的对应时间;
计算最容易出现损伤的区域中每个测量点宏观应力的最大间隔值:
Figure FDA0002382348040000021
其中,σJq为最容易出现损伤的区域中每个测量点宏观应力的最大间隔值;Cmaxqt)为第q个测量点在0到tm时刻内所有监测宏观应力中的最大值;Cminqt)为第q个测量点在0到tm时刻内所有监测宏观应力中的最小值;tmax为Cmaxqt)所对应的时间点;tmin为Cminqt)所对应的时间点;σm为数据监测的总时间;
Figure FDA0002382348040000022
为第q个测量点的宏观应力平均值;
Figure FDA0002382348040000023
然后计算工况应力条件下结构宏观应力的畸变区域系数B
Figure FDA0002382348040000024
其中,A为有限元软件中箱型结构焊缝处的应力分布及最容易出现损伤的区域的面积;C(σJq)max为n个测量点宏观应力最大间隔值中的最大值;C(σJq)min为n个测量点宏观应力最大间隔值中的最小值;y为最容易出现损伤的区域的面积变量;
S3.基于畸变区域系数的焊缝区域残余应力修正:
在上述过程的基础上,基于畸变系数对焊缝区域残余应力进行修正;
计算修正后的残余应力值σc
Figure FDA0002382348040000025
其中,B为工况应力条件下结构宏观应力的畸变区域系数;i为取点的标号,i=1,2,3...N;N为残余应力测量时取点的最大值;σri为箱型结构焊缝处每测量点的残余应力总值;β为箱型结构所对应的名义残余应力混合修正系数;
S4.考虑工况应力峰谷效应的残余应力拓度计算:
在上述过程的基础上,计算考虑工况应力峰谷效应的残余应力拓度D
Figure FDA0002382348040000031
其中,D为虑工况应力峰谷效应的残余应力拓度;
Figure FDA0002382348040000032
为为第q个测量点的宏观应力平均值;Cmaxqt)为第q个测量点在0到tm时刻内所有监测宏观应力中的最大值;Cminqt)为第q个测量点在0到tm时刻内所有监测宏观应力中的最小值;σc修正后的残余应力值;β为箱型结构所对应的名义残余应力混合修正系数;
S5.基于宏观应力与残余应力拓度的箱型结构重变形量计算:
在上述过程的基础上,基于宏观应力与残余应力拓度计算最容易出现损伤的区域中每个测量点的结构重变形值εq,q为测量点的标记,q=1,2,...n;
Figure FDA0002382348040000033
其中,
Figure FDA0002382348040000034
为为第q个测量点的宏观应力平均值;D为考虑工况应力峰谷效应的残余应力拓度;σc为修正后的残余应力值;E为弹性模量。
2.根据权利要求1所述的预测方法,其特征在于:S1中,残余应力测试装置对箱型结构焊缝处每隔5mm取一个点i进行残余应力的测试。
3.根据权利要求1或2所述的预测方法,其特征在于:残余应力测试装置为X射线残余应力测试装置。
4.根据权利要求1所述的预测方法,其特征在于:S2至S5中,n取值为7。
5.根据权利要求1所述的预测方法,其特征在于:S5中,E取值为2.06×1011Pa。
CN202010086804.3A 2020-02-11 2020-02-11 基于宏观应力与残余应力拓度的箱型结构重变形量预测方法 Withdrawn CN111339699A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010086804.3A CN111339699A (zh) 2020-02-11 2020-02-11 基于宏观应力与残余应力拓度的箱型结构重变形量预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010086804.3A CN111339699A (zh) 2020-02-11 2020-02-11 基于宏观应力与残余应力拓度的箱型结构重变形量预测方法

Publications (1)

Publication Number Publication Date
CN111339699A true CN111339699A (zh) 2020-06-26

Family

ID=71185139

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010086804.3A Withdrawn CN111339699A (zh) 2020-02-11 2020-02-11 基于宏观应力与残余应力拓度的箱型结构重变形量预测方法

Country Status (1)

Country Link
CN (1) CN111339699A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111241741A (zh) * 2020-03-25 2020-06-05 南京市特种设备安全监督检验研究院 一种基于残余应力影响区间修正的起重机轮轨磨损监测方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111241741A (zh) * 2020-03-25 2020-06-05 南京市特种设备安全监督检验研究院 一种基于残余应力影响区间修正的起重机轮轨磨损监测方法
CN111241741B (zh) * 2020-03-25 2023-04-28 南京市特种设备安全监督检验研究院 一种基于残余应力影响区间修正的起重机轮轨磨损监测方法

Similar Documents

Publication Publication Date Title
CN108895974B (zh) 一种结构变形光纤监测与重构方法及***
CN106897543B (zh) 模态柔度曲率矩阵范数的梁结构损伤识别方法
CN110319808B (zh) 倾角传感器测量大跨度拱桥拱肋变形精度预测及评估方法
CN102778403B (zh) 一种焊缝材料参数识别方法
CN107817395B (zh) 一种基于电网线路运行误差远程校准的窃电排查方法
CN108801141B (zh) 基于激光测量技术的超高压隔离开关智能调试校准***
CN105956216A (zh) 大跨钢桥基于均匀温度响应监测值的有限元模型修正方法
CN101349731B (zh) 一种电压稳定性实时评估方法
CN103838931A (zh) 一种工程机械臂架类结构再制造准入期评估方法
CN104535011B (zh) 一种三维线激光车辙检测设备的室内参数标定方法
CN110487576B (zh) 损伤状态倾角对称斜率的等截面梁损伤识别方法
Povarov et al. Data support system for controlling decentralised nuclear power industry facilities through uninterruptible condition monitoring
CN116577051A (zh) 一种虑及损伤等效的随机振动疲劳寿命分析方法
CN109992579A (zh) 一种公路基础设施多源异质数据的数据修复方法及***
CN111339699A (zh) 基于宏观应力与残余应力拓度的箱型结构重变形量预测方法
CN105068035A (zh) 一种电压互感器误差水平动态检测方法及***
CN1782672A (zh) 改善发电设备故障检测的方法和装置
CN102269581B (zh) 一种建筑构件加工精度测评方法
CN111597742B (zh) 一种基于焊接构件联动效应的农机结构残余变形预测方法
CN107679630A (zh) 一种基于比例维修模型的维修作业时间估计方法
CN117195662A (zh) 一种大型结构件焊接变形预测方法及***
CN109979622B (zh) 核电厂稳压器疲劳寿命在线监测评估***与方法
CN112966358B (zh) 一种基于数据修补的主动配电网状态感知方法
CN113868925A (zh) 一种折减系数确定方法和托架结构的设计方法
CN110954349B (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
WW01 Invention patent application withdrawn after publication
WW01 Invention patent application withdrawn after publication

Application publication date: 20200626