CN117817417A - 一种基于数字孪生的大型筒体最优加工位姿调整方法 - Google Patents

一种基于数字孪生的大型筒体最优加工位姿调整方法 Download PDF

Info

Publication number
CN117817417A
CN117817417A CN202410249588.8A CN202410249588A CN117817417A CN 117817417 A CN117817417 A CN 117817417A CN 202410249588 A CN202410249588 A CN 202410249588A CN 117817417 A CN117817417 A CN 117817417A
Authority
CN
China
Prior art keywords
cylinder
processed
point
point cloud
cloud 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.)
Pending
Application number
CN202410249588.8A
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.)
North University of China
Original Assignee
North University of China
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 North University of China filed Critical North University of China
Priority to CN202410249588.8A priority Critical patent/CN117817417A/zh
Publication of CN117817417A publication Critical patent/CN117817417A/zh
Pending legal-status Critical Current

Links

Landscapes

  • Numerical Control (AREA)

Abstract

本发明属于筒体位姿调整技术领域,解决了大型筒体辅助找正时间长、加工精度低的问题。提供了一种基于数字孪生的大型筒体最优加工位姿调整方法,包括以下步骤:构建待加工筒体的重构点云模型和理想点云模型;通过基于双向KD‑tree的ICP算法,获取待加工筒体的重构点云模型和理想点云模型间的位置差值矩阵;得到待加工筒体在坐标轴X、Y、Z方向的三轴调整参数,控制调姿伺服电机进行待加工筒体的位姿调整;待数字孪生系统中待加工筒体的重构点云模型和理想点云模型完全重合后;再次采集待加工筒体的点云数据进行校核,保证待加工筒体的位姿在误差要求范围内。本发明够解决大型筒体加工过程中的精度低、调姿时间长的问题。

Description

一种基于数字孪生的大型筒体最优加工位姿调整方法
技术领域
本发明属于筒体位姿调整技术领域,具体涉及一种基于数字孪生的大型筒体最优加工位姿调整方法。
背景技术
在大型筒体内孔加工过程中,针对筒体内径φ2110mm,端框外径φ2300mm,筒体长度2000mm~4800mm,最大质量3000kg以上的筒体,加工后要求筒壁厚度大于14mm,内圆圆度为0.3mm,端框端面与轴线的垂直度为0.08mm。
以上关键参数均由立车加工工序完成,操作过程涉及到筒体的找正;目前主要采用钳工划线预估筒体偏差,立车参照加工线进行找正,整个过程均在机床内完成;该找正方法中,机床辅助找正时间长、加工精度低且存在人为误差。
发明内容
本发明为了解决现有技术中存在的上述至少一个技术问题,提供了一种基于数字孪生的大型筒体最优加工位姿调整方法。
本发明采用如下的技术方案实现:一种基于数字孪生的大型筒体最优加工位姿调整方法,包括以下步骤:
S1:采集待加工工位上筒体的点云数据,构建待加工筒体的重构点云模型;
S2:根据待加工筒体理想加工状态下的位姿数据,构建待加工筒体的理想点云模型;
S3:通过基于双向KD-tree的ICP算法,获取待加工筒体的重构点云模型和理想点云模型间的位置差值矩阵;
S4:由位置差值矩阵得到待加工筒体在坐标轴X、Y、Z方向的三轴调整参数,控制调姿伺服电机进行待加工筒体的位姿调整;
S5:构建待加工筒体调姿的数字孪生系统,将待加工筒体的物理实体的实时位姿变化映射给待加工筒体的重构点云模型;
S6:待数字孪生系统中待加工筒体的重构点云模型和理想点云模型完全重合后;再次采集待加工筒体的点云数据进行校核,保证待加工筒体的位姿在误差要求范围内。
优选地,基于双向KD-tree的ICP算法包括下述步骤:
S31:分别构建待加工筒体的重构点云模型和理想点云模型上的点集合的KD-tree;
S32:选取重构点云模型上的任意一点P作为目标点,通过基于双向KD-tree的ICP算法搜索待加工筒体的理想点云模型上距离点P的最近点Q;
S33:通过基于双向KD-tree的ICP算法搜索待加工筒体的重构点云模型上距离点Q的最近点,若点Q的最近点为点P,则点Q作为点P的配准点;若点Q的最近点为除点P以外的其他点,则将重构点云模型上的点P作为噪声点剔除;
S34:选取重构点云模型上除P以外的任意一点作为目标点,根据步骤S32、S33中的方法,选取该点在理想点云模型上的配准点;
S35:遍历重构点云模型上的所有点,获取在理想点云模型上对应的配准点,实现待加工筒体的重构点云模型和理想点云模型的配准,获取两者的位置差值矩阵。
优选地,设置目标点和与其对应的配准点之间的距离阈值,若目标点与配准点之间的欧式距离小于等于距离阈值时,保留该组配准后的点;若目标点与配准点之间的欧式距离大于距离阈值时,则将该目标点作为噪声点剔除。
优选地,待加工筒体的物理实体的移动通过PLC控制调姿伺服电机的进给量实现;在Unity软件中通过借助Nmodbus库函数读写控制待加工筒体的物理实体移动的调姿伺服电机的寄存器地址,进而实现待加工筒体的物理实体和重构点云模型的虚实交互。
优选地,调姿伺服电机设置在待加工筒体的底部,包括控制待加工筒体沿X、Y轴方向移动的伺服电机和控制待加工筒体沿Z轴方向移动的伺服电机。
优选地,通过设置在待加工筒体的物理实体内的点云扫描仪在水平面的旋转和竖直面的平移完成筒体内表面点云数据的采集,将采集的数据进行点云重构,得到待加工筒体的重构点云模型。
与现有技术相比,本发明的有益效果是:
本发明通过基于双向KD-tree的ICP算法,获取待加工筒体的重构点云模型和理想点云模型间的位置差值矩阵,反应出两者间的位姿差值;并将位姿差值转换为调姿伺服电机的进给量,实现待加工筒体的位姿调整,该方法能够解决大型筒体加工过程中的精度低、调姿时间长的问题。并通过构建待加工筒体的物理实体和重构点云模型间的数字孪生系统,实现待加工筒体的物理实体调姿过程的可视化。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单的介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明中待加工筒体调姿的三维模型图;
图2是传统的ICP点云配准算法的配准结果图;
图3是本发明的基于双向KD-tree的ICP算法的配准结果图;
图4是本发明中点云扫描的俯视图;
图5是本发明中点云扫描的正视图。
图中:1-待加工筒体的物理实体;2-点云扫描仪;3-调姿伺服电机。
具体实施方式
结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚,完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部实施例。基于本发明的实施例,本领域的普通技术人员在没有做出创造性劳动的前提下所得到的所有其他实施方式,都属于本发明所保护的范围。
须知,本说明书附图所绘示的结构、比例、大小等,均仅用以配合说明书所揭示的内容,以供熟悉此技术的人士了解与阅读,并非用以限定本发明可实施的限定条件,故不具技术上的实质意义,任何结构的修饰、比例关系的改变或大小的调整,在不影响本发明所能产生的功效及所能达成的目的下,均应落在本发明所揭示的技术内容能涵盖的范围内,需要说明的是,在本说明书中,诸如第一和第二之类的关系术语仅仅用来将一个实体与另外几个实体区分开来,而不一定要求或者暗示这些实体之间存在任何实际的关系或者顺序。
本发明提供了一种实施例:
如图1所示,一种基于数字孪生的大型筒体最优加工位姿调整方法,包括以下步骤:
S1:采集待加工工位上筒体的点云数据,构建待加工筒体的重构点云模型;
S2:建立待加工筒体调姿的三维模型,包括底座工装夹具和待加工筒体的理想点云模型;根据待加工筒体理想加工状态下的位姿数据,构建待加工筒体的理想点云模型;
S3:通过基于双向KD-tree的ICP算法,获取待加工筒体的重构点云模型和理想点云模型间的位置差值矩阵;
S4:由位置差值矩阵得到待加工筒体在坐标轴X、Y、Z方向的三轴调整参数,控制调姿伺服电机3进行待加工筒体的位姿调整;
S5:构建待加工筒体调姿的数字孪生系统,将待加工筒体的物理实体1的实时位姿变化映射给待加工筒体的重构点云模型;
S6:待数字孪生系统中待加工筒体的重构点云模型和理想点云模型完全重合后;再次采集待加工筒体的点云数据进行校核,保证待加工筒体的位姿在误差要求范围内。
如图4、图5所示,构建待加工筒体的重构点云模型:点云扫描仪2采集到自身到筒体内壁的距离和自身到机床主轴的距离/>;通过设置在待加工筒体的物理实体1内的点云扫描仪2在水平面的旋转和竖直面的平移完成筒体内表面点云数据的采集,每次旋转角度为/>,旋转次数记为/>(/>),相机距机床加工平面为/>;获取的点云数据的中点坐标计算方式为:
将采集的数据进行点云重构,得到待加工筒体的重构点云模型。
本实施例中,基于双向KD-tree的ICP算法包括下述步骤:
S31:分别构建待加工筒体的重构点云模型和理想点云模型上的点集合的KD-tree;
S32:选取重构点云模型上的任意一点P作为目标点,通过基于双向KD-tree的ICP算法搜索待加工筒体的理想点云模型上距离点P的最近点Q;
S33:通过基于双向KD-tree的ICP算法搜索待加工筒体的重构点云模型上距离点Q的最近点,若点Q的最近点为点P,则点Q作为点P的配准点;若点Q的最近点为除点P以外的其他点,则将重构点云模型上的点P作为噪声点剔除;
S34:选取重构点云模型上除P以外的任意一点作为目标点,根据步骤S32、S33中的方法,选取该点在理想点云模型上的配准点;
S35:遍历重构点云模型上的所有点,获取在理想点云模型上对应的配准点,实现待加工筒体的重构点云模型和理想点云模型的配准,获取两者的位置差值矩阵。
目前应用最广泛的点云精配准算法是迭代最近点算法(Iterative ClosestPoint, ICP)。ICP算法的基本逻辑是先通过对应关系估计、关键点检测等方法找到源点云和目标点云的待匹配点对,然后计算旋转和平移矩阵,对source点云进行旋转平移到target点云上,根据设置的阈值进行判断,如果不满足阈值要求就重复以上过程继续计算,最终达到最大迭代次数或者满足设定的均方差阈值才能停止迭代。
采用传统的ICP点云配准算法得出的大型筒体配准结果如图2所示(图中竖直的筒体为理想点云模型,倾斜角度大的筒体为重构点云模型,倾斜角度小的筒体为采用传统ICP点云配准算法配准后的点云模型),配准精度为:;对应的位置差值矩阵为:
本发明采用的基于双向KD-tree的ICP算法,适用于配准数据量大的点云模型,且配准精度相比传统的IPC点云配准算法大幅提升。
采用基于双向KD-tree的ICP算法得出的大型筒体配准结果如图3所示(图中竖直的筒体为理想点云模型,倾斜角度大的筒体为重构点云模型,与竖直的筒体基本重合的筒体为采用双向KD-tree的点云配准算法配准后的点云模型,在图中基本重合,难以分辨),配准精度为:;对应的位置差值矩阵为:
位置差值矩阵是一个由一个的旋转矩阵和一个/>的平移向量组成的矩阵,可表示为:
其中,是/>的旋转矩阵的元素,分别表示物体绕X、Y、Z轴旋转的角度。/>是平移向量的元素,表示物体在X、Y、Z方向上的位移。
故位置差值矩阵中即为待加工筒体在坐标轴X、Y、Z方向的三轴调整参数;将坐标轴X、Y、Z方向的三轴调整参数转换为电机进给量的方法为现有技术,再此不做赘述。
由此可知,传统的ICP点云配准算法精度为,基于双向KD-tree的ICP算法精度为/>,精度提升1倍;基于双向KD-tree的ICP算法适用于配准数据量大的点云模型,且配准精度大幅提升。
本实施例中,还对应设置有目标点和与其对应的配准点之间的距离阈值,若目标点与配准点之间的欧式距离小于等于距离阈值时,保留该组配准后的点;若目标点与配准点之间的欧式距离大于距离阈值时,则将该目标点作为噪声点剔除。
根据基于双向KD-tree的ICP算法和距离阈值的设置,可快速剔除配准精度低的点,保留配准精度高的点。
本实施例中,通过构建待加工筒体调姿的数字孪生系统,实现大型筒体调姿的可视化,达到实时监测大型筒体的位姿,实时读写PLC寄存器地址得到调姿伺服电机3的进给量,将物理实体的运动位移映射给重构点云模型,实现物理实体与重构点云模型的虚实映射,在可视化界面精确反应大型筒体实时位姿,监测筒体的位姿是否满足加工需求。
待加工筒体的物理实体1的移动通过PLC控制调姿伺服电机3的进给量实现,调姿伺服电机3设置在待加工筒体的底部,包括控制待加工筒体沿X、Y轴方向移动的伺服电机和控制待加工筒体沿Z轴方向移动的伺服电机。在Unity客户端中通过读写PLC寄存器地址就可以做到控制待加工筒体的物理实体1的运动,在Unity中如果直接引入读写串口的命名空间,会涉及到复杂的进制转换,因此选择在Unity中使用Nmodbus库函数来利用串口对寄存器进行读写。NModbus是用C#编写的Modbus通信协议库,它支持的Modbus协议包括ModbusRTU、ASCII、TCP和UDP,可用于编程读写Modbus设备的寄存器。在Unity软件中通过借助Nmodbus库函数读写控制待加工筒体的物理实体1移动的调姿伺服电机3的寄存器地址,进而实现待加工筒体的物理实体1和重构点云模型的虚实交互。读写寄存器时所产生的数据,同样也可用于控制Unity中的重构点云模型。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应该涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。

Claims (6)

1.一种基于数字孪生的大型筒体最优加工位姿调整方法,其特征在于,包括以下步骤:
S1:采集待加工工位上筒体的点云数据,构建待加工筒体的重构点云模型;
S2:根据待加工筒体理想加工状态下的位姿数据,构建待加工筒体的理想点云模型;
S3:通过基于双向KD-tree的ICP算法,获取待加工筒体的重构点云模型和理想点云模型间的位置差值矩阵;
S4:由位置差值矩阵得到待加工筒体在坐标轴X、Y、Z方向的三轴调整参数,控制调姿伺服电机(3)进行待加工筒体的位姿调整;
S5:构建待加工筒体调姿的数字孪生系统,将待加工筒体的物理实体(1)的实时位姿变化映射给待加工筒体的重构点云模型;
S6:待数字孪生系统中待加工筒体的重构点云模型和理想点云模型完全重合后;再次采集待加工筒体的点云数据进行校核,保证待加工筒体的位姿在误差要求范围内。
2.根据权利要求1所述的一种基于数字孪生的大型筒体最优加工位姿调整方法,其特征在于:基于双向KD-tree的ICP算法包括下述步骤:
S31:分别构建待加工筒体的重构点云模型和理想点云模型上的点集合的KD-tree;
S32:选取重构点云模型上的任意一点P作为目标点,通过基于双向KD-tree的ICP算法搜索待加工筒体的理想点云模型上距离点P的最近点Q;
S33:通过基于双向KD-tree的ICP算法搜索待加工筒体的重构点云模型上距离点Q的最近点,若距离点Q的最近点为点P,则点Q作为点P的配准点;若距离点Q的最近点为除点P以外的其他点,则将重构点云模型上的点P作为噪声点剔除;
S34:选取重构点云模型上除P以外的任意一点作为目标点,根据步骤S32、S33中的方法,选取该点在理想点云模型上的配准点;
S35:遍历重构点云模型上的所有点,获取在理想点云模型上对应的配准点,实现待加工筒体的重构点云模型和理想点云模型的配准,进而获取两者的位置差值矩阵。
3.根据权利要求2所述的一种基于数字孪生的大型筒体最优加工位姿调整方法,其特征在于:设置目标点和与其对应的配准点之间的距离阈值,若目标点与配准点之间的欧式距离小于等于距离阈值时,保留该组配准后的点;若目标点与配准点之间的欧式距离大于距离阈值时,则将该目标点作为噪声点剔除。
4.根据权利要求1所述的一种基于数字孪生的大型筒体最优加工位姿调整方法,其特征在于:待加工筒体的物理实体(1)的移动通过PLC控制调姿伺服电机(3)的进给量实现;在Unity软件中通过借助Nmodbus库函数读写控制待加工筒体的物理实体(1)移动的调姿伺服电机(3)的寄存器地址,进而实现待加工筒体的物理实体(1)和重构点云模型的虚实交互。
5.根据权利要求1所述的一种基于数字孪生的大型筒体最优加工位姿调整方法,其特征在于:调姿伺服电机(3)设置在待加工筒体的底部,包括控制待加工筒体沿X、Y轴方向移动的伺服电机和控制待加工筒体沿Z轴方向移动的伺服电机。
6.根据权利要求1所述的一种基于数字孪生的大型筒体最优加工位姿调整方法,其特征在于:通过设置在待加工筒体的物理实体(1)内的点云扫描仪(2)在水平面的旋转和竖直面的平移完成筒体内表面点云数据的采集,根据采集的数据进行点云重构,得到待加工筒体的重构点云模型。
CN202410249588.8A 2024-03-05 2024-03-05 一种基于数字孪生的大型筒体最优加工位姿调整方法 Pending CN117817417A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202410249588.8A CN117817417A (zh) 2024-03-05 2024-03-05 一种基于数字孪生的大型筒体最优加工位姿调整方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202410249588.8A CN117817417A (zh) 2024-03-05 2024-03-05 一种基于数字孪生的大型筒体最优加工位姿调整方法

Publications (1)

Publication Number Publication Date
CN117817417A true CN117817417A (zh) 2024-04-05

Family

ID=90509964

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202410249588.8A Pending CN117817417A (zh) 2024-03-05 2024-03-05 一种基于数字孪生的大型筒体最优加工位姿调整方法

Country Status (1)

Country Link
CN (1) CN117817417A (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110103071A (zh) * 2019-05-28 2019-08-09 内蒙古大学 一种变形复杂零件的数字化寻位加工方法
CN111145240A (zh) * 2019-11-18 2020-05-12 西宁市动物疫病预防控制中心(挂西宁市畜牧兽医站牌子) 一种基于3d相机的活体西门塔尔牛体尺在线测量方法
CN116108722A (zh) * 2023-02-28 2023-05-12 南京理工大学 基于数字孪生的大型结构件面形调控方法
CN116402866A (zh) * 2023-03-03 2023-07-07 西安交通大学 基于点云的零件数字孪生几何建模与误差评定方法及系统
US20230230312A1 (en) * 2022-01-20 2023-07-20 Qingdao university of technology Digital twin modeling method and system for assembling a robotic teleoperation environment
CN117020246A (zh) * 2023-10-09 2023-11-10 中北大学 一种大型筒体姿态调整方法
CN117274339A (zh) * 2023-10-24 2023-12-22 南京信息工程大学 一种基于改进的iss-3dsc特征结合icp的点云配准方法
CN117291959A (zh) * 2023-09-25 2023-12-26 太原理工大学 基于激光slam的工作面整体工作空间虚拟重构方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110103071A (zh) * 2019-05-28 2019-08-09 内蒙古大学 一种变形复杂零件的数字化寻位加工方法
CN111145240A (zh) * 2019-11-18 2020-05-12 西宁市动物疫病预防控制中心(挂西宁市畜牧兽医站牌子) 一种基于3d相机的活体西门塔尔牛体尺在线测量方法
US20230230312A1 (en) * 2022-01-20 2023-07-20 Qingdao university of technology Digital twin modeling method and system for assembling a robotic teleoperation environment
CN116108722A (zh) * 2023-02-28 2023-05-12 南京理工大学 基于数字孪生的大型结构件面形调控方法
CN116402866A (zh) * 2023-03-03 2023-07-07 西安交通大学 基于点云的零件数字孪生几何建模与误差评定方法及系统
CN117291959A (zh) * 2023-09-25 2023-12-26 太原理工大学 基于激光slam的工作面整体工作空间虚拟重构方法
CN117020246A (zh) * 2023-10-09 2023-11-10 中北大学 一种大型筒体姿态调整方法
CN117274339A (zh) * 2023-10-24 2023-12-22 南京信息工程大学 一种基于改进的iss-3dsc特征结合icp的点云配准方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
乔世权;张坤;: "基于改进ICP算法的曲面特征配准研究", 数学的实践与认识, no. 08, 23 April 2019 (2019-04-23), pages 137 - 145 *
何黎明等: "中国物流技术发展报告", vol. 1, 30 April 2022, 中国财富出版社, pages: 240 - 242 *
杜世立等: "多源数据融合的三维建模技术研究", 淮阴工学院学报, vol. 32, no. 3, 15 June 2023 (2023-06-15), pages 78 - 87 *
郑旻辉等: "基于改进ICP 拼接算法的柔性抛光气囊 工具磨损检测方法", 激光与光电子学进展, vol. 60, no. 16, 25 August 2023 (2023-08-25), pages 1612001 - 10 *

Similar Documents

Publication Publication Date Title
CN109598762B (zh) 一种高精度双目相机标定方法
CN109272574B (zh) 基于投影变换的线阵旋转扫描相机成像模型构建方法和标定方法
CN112669385B (zh) 基于三维点云特征的工业机器人工件识别与位姿估计方法
JP4785880B2 (ja) 三次元オブジェクト認識のためのシステムおよび方法
CN109584288B (zh) 一种五轴系统中三维模型的重构方法及系统
JP5606786B2 (ja) キャリブレーション方法及びキャリブレーション装置
CN110717943A (zh) 用于二维平面的眼在手上机械手手眼标定方法及系统
CN111028340B (zh) 精密装配中的三维重构方法、装置、设备及系统
CN112614098B (zh) 一种基于增强现实的毛坯定位与加工余量分析方法
CN112894815B (zh) 视觉伺服机械臂物品抓取最佳位姿检测方法
CN112781496A (zh) 一种非接触测量系统的测头位姿标定技术
WO2018201677A1 (zh) 基于光束平差的远心镜头三维成像系统的标定方法及装置
CN109015110A (zh) 一种机床回转轴位置误差建模与辨识方法
CN111400830B (zh) 一种三维毛坯工件的加工校准方法及装置
CN111027484A (zh) 基于三维成像的隧道钢拱识别方法
CN112361958B (zh) 一种线激光与机械臂标定方法
CN115077378A (zh) 一种应用于大型工件测量的线扫激光的手眼标定方法
CN110363801B (zh) 工件实物与工件三维cad模型的对应点匹配方法
Hesch et al. Mirror-based extrinsic camera calibration
CN110992416A (zh) 基于双目视觉与cad模型的高反光面金属零件位姿测量方法
CN117817417A (zh) 一种基于数字孪生的大型筒体最优加工位姿调整方法
CN107421476A (zh) 一种空间孔位测量基准误差补偿方法
CN109035238B (zh) 一种面向自由曲面零件的加工余量离线分析方法
CN113393507B (zh) 无人机点云与地面三维激光扫描仪点云配准方法
CN116399262A (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