CN111230887B - 一种基于数字孪生技术的工业涂胶机器人运行状态监测方法 - Google Patents

一种基于数字孪生技术的工业涂胶机器人运行状态监测方法 Download PDF

Info

Publication number
CN111230887B
CN111230887B CN202010161152.5A CN202010161152A CN111230887B CN 111230887 B CN111230887 B CN 111230887B CN 202010161152 A CN202010161152 A CN 202010161152A CN 111230887 B CN111230887 B CN 111230887B
Authority
CN
China
Prior art keywords
state
digital twin
data
fatigue
gluing robot
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.)
Active
Application number
CN202010161152.5A
Other languages
English (en)
Other versions
CN111230887A (zh
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.)
Hefei University
Original Assignee
Hefei 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 Hefei University filed Critical Hefei University
Priority to CN202010161152.5A priority Critical patent/CN111230887B/zh
Publication of CN111230887A publication Critical patent/CN111230887A/zh
Application granted granted Critical
Publication of CN111230887B publication Critical patent/CN111230887B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J9/00Programme-controlled manipulators
    • B25J9/16Programme controls
    • B25J9/1628Programme controls characterised by the control loop
    • B25J9/1653Programme controls characterised by the control loop parameters identification, estimation, stiffness, accuracy, error analysis
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J11/00Manipulators not otherwise provided for
    • B25J11/0075Manipulators for painting or coating
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J19/00Accessories fitted to manipulators, e.g. for monitoring, for viewing; Safety devices combined with or specially adapted for use in connection with manipulators
    • B25J19/0095Means or methods for testing manipulators

Abstract

本发明涉及一种基于数字孪生技术的工业涂胶机器人运行状态监测方法,与现有技术相比解决了难以针对工业涂胶机器人进行运行状态预测的缺陷。本发明包括以下步骤:获取传感器所采集的工业涂胶机器人运行状态数据;工业涂胶机器人运行状态数据的预处理;构建数字孪生设备驱动模型;设备运行状态的初期评价;数字孪生设备驱动模型的修正;工业涂胶机器人运行状态的监测。本发明通过利用数字孪生驱动的设备运行模型,实现了工业涂胶机器人运行状态的实时评价,从而获得更为精准的操作建议、运维建议,实现设备运行的生命周期有效管理。

Description

一种基于数字孪生技术的工业涂胶机器人运行状态监测方法
技术领域
本发明涉及工业设备运行监测技术领域,具体来说是一种基于数字孪生技术的工业涂胶机器人运行状态监测方法。
背景技术
智能制造行业是立国之本、强国之基,从根本上决定着一个国家的综合实力和国际竞争力,我国经济进入新常态,制造业转型升级进入关键时期。运用互联网、物联网、机器人等新技术、新业态、新模式,大力改造提升传统制造业,将是我国从制造大国走向制造强国的重要一步。在利用智能机器人作业时,设备的稳定运行与准确的设备状态预测是保证安全生产的关键环节,对产业升级与企业经济效益的提升具有重大意义。
工业涂胶机器人在作业时,由于需要进行重物抓取、涂胶界面自我识别、喷涂、对接安装、放回、归位等复杂环节,长时间连续作业可能引起作业手臂的物理偏移。而目前的作业臂故障预警,往往是在出现较大设备故障或出现产线质量残次时的报警,这对产线正常生产及智能工厂效率总体提升产生不良影响。同时,复杂环节的作业臂操作所受的影响因素较多,难以利用简单的分析手段实现耦合数据的综合性能分析。
在已有的寿命预测技术中,多采用运行时间作为评价指标,量化指标单一,对于具有复杂、多维状态的工业机器人作业臂,难以实现准确的预判;且寿命预测后的健康指标评价对产线设备稳定生产意义重大。通过智能化软件服务系统及其辅助装置,可以实现设备在运行工况生命周期内的高效管控,实现产品生产的自适应与动态响应。
数字孪生技术是充分利用物理模型、传感器更新、运行历史等数据,集成多学科、多物理量、多尺度、多概率的仿真过程,在虚拟空间中完成映射,从而反映相对应的实体装备的全生命周期过程。数字孪生思想是高维空间中的复杂系统,在数字孪生思想的背景下,需要解决耦合数据融合,量化不确定性,设计和优化结构,同时处理大而嘈杂数据等诸多问题。
目前国内外对产品数字孪生体的系统性研究成果较少,产品数字孪生体的实现主要分布在产品设计、工艺设计、产品制造、产品服务、产品报废回收等阶段,用于模拟、监控、诊断、预测和控制产品在现实环境中的形成过程和行为。尽管已经初步探索了数字孪生的相关应用,但数字孪生的关键技术在与智能制造生产系统相结合的技术研究方面尚不成熟,主要表现在:1、在模型构建、信息物理数据融合、交互与协同等方面的理论与技术比较缺乏;2、针对独立的生产、装配设备,难以适应不同产品生产周期的变化等问题,构建完善的实物与装配模型的关联机制方面缺乏深度融合与应用;3、故障捕捉,虚实自主交互方面,由于实际干扰因素的影响,难以获得精准判断。
因此如何利用数字孪生思想,实现设备与模型的虚实交互,预测设备在生产运行阶段生命周期内的疲劳状态、运维状态已经成为急需解决的技术问题。
发明内容
本发明的目的是为了解决现有技术中难以针对工业涂胶机器人进行运行状态预测的缺陷,提供一种基于数字孪生技术的工业涂胶机器人运行状态监测方法来解决上述问题。
为了实现上述目的,本发明的技术方案如下:
一种基于数字孪生技术的工业涂胶机器人运行状态监测方法,所述的工业涂胶机器人包括往复运动关节组件,往复运动关节组件的关联部位安装有用于获取位置信息三维数值的传感器,往复运动关节组件的作业臂前端安装有姿态传感器,图像采集装置的图像采集范围位于作业臂前端,位差监测终端的监测范围位于作业臂前端,所述的工业涂胶机器人运行状态监测方法包括以下步骤:
11)获取传感器所采集的工业涂胶机器人运行状态数据:所述的运行状态数据包括设备位置三维数值、姿态图片、数值、信号,反映设备环境状态的温度、噪声数值信息,反映设备运行状态的数值信息,反映设备生产运行状态的质量信息;
12)工业涂胶机器人运行状态数据的预处理:对工业涂胶机器人的运行状态数据进行初步比对,剔除干扰数据;
13)构建数字孪生设备驱动模型:根据工业涂胶机器人的分析需要构建数字孪生设备驱动模型;
14)设备运行状态的初期评价:将预处理后的工业涂胶机器人运行状态数据输入数字孪生设备驱动模型,数字孪生设备驱动模型以虚拟方式实现工业涂胶机器人作业,得到工业涂胶机器人随运行时间变化而演化的运行状态数值;
15)数字孪生设备驱动模型的修正:将数字孪生设备驱动模型演化后的运行状态数值与实时传输的设备状态进行对比匹配,实现数字孪生设备驱动模型的修正;
16)工业涂胶机器人运行状态的监测:获取传感器采集的工业涂胶机器人运行状态实时数据,将实时数据预处理后输入修正后的数字孪生设备驱动模型,预测出工业涂胶机器人运行的远期趋势,判定出运行疲劳点。
所述工业涂胶机器人运行状态数据的预处理使用经验模式分析与辅助机械报警相结合的处理方法,其具体步骤如下:
21)获取利用经验判断建立的专家库,输入数据符合专家库规则的判定为稳定数据,将辅助报警装置提供的故障报警规则作为排他判定,接近排他判定的运行数据,作为干扰数据进行剥离;
22)将经过预处理的数据进行多源数据融合:使用带有权重分析的特征融合法对影响模拟性能的关键数据进行权重分析,并抽取不同来源数据的数据特征开展模型训练,对不同结果进行合并,使得数字孪生模型中的数据特征能够准确反映工业机器人作业臂的实际工作状态用于下一步的状态预测,其表达式如下:
Figure BDA0002405835630000031
式中:k为影响因素的权重个数,u为影响因素代号,reck(u)为各项权重值,rec(u)为综合权重分析值。
所述构建数字孪生设备驱动模型包括以下步骤:
31)利用三维模型构建数字孪生设备驱动模型,实现通过逻辑控制构建出基于数字孪生的工业机器人作业臂运行工作特征;
32)利用虚拟引擎开发软件,开展针对三维模型的深度个性化定制,使其准确反映工业涂胶机器人作业臂的运行工况特性;自定义三维模型及适用于机器人产线工业作业水平的工艺流程;建立适用于工业涂胶机器人作业特征的控制逻辑规划,利用计算机语言及虚拟引擎的控制函数,完成逻辑控制触发,实现虚拟引擎软件开发环境下的适用于工业涂胶机器人真实产线设备的数字孪生体,用于运行特征的实时展现与关键部位、关键时间节点的个性化信息定制与评估;
33)将预处理后的数据信息,输入构建的数字孪生设备驱动模型当中,利用计算机语言及虚拟引擎的控制函数,调整逻辑控制,完成满足工业涂胶机器人实际个性化工作需求的点对点位置控制,实现实时数据输入与响应运行状态的耦合同步;通过在数字孪生虚拟开发环境中添加虚拟传感器,获得虚拟操作环境的工艺状态特征,使得数字孪生驱动系统具备在任意目标状态时点、任意目标操作位置实现工业涂胶机器人实时状态监控与个性化信息提取的目的;利用自定义控制函数,将预处理后的数据信息,在数字孪生设备驱动模型中建立数据偏差分析,结合耦合权重获得随记录时间变化的点对点偏差累积函数η(t),用于远期偏差预测;
34)初步完善的数字孪生设备驱动模型,结合已采集信号的已知偏差,利用工业涂胶机器人实际运行工况的运行经验组建经验库,建立基于已知偏差的置信区间;利用自定义控制函数,根据实际工业涂胶机器人的工艺控制需要,选择相应的控制手段或编译自定义控制特征,采用记忆累加方式对偏差进行累加触发并满足△S(uii)=△S(ui0)+η(t),使得数字孪生驱动模型具备适用于工业涂胶机器人个性化运行需求的随时间变化的远期演化分析能力;
其中△S(u)为系统触发故障偏差,△S(uii)为不同采集特征的末次触发偏差,△S(ui0)为新触发前的原始偏差,u为影响因素代号,η(t)为随时间变化偏差累积函数;
35)设定构建数字孪生设备驱动模型利用持续采集的设备运行特征与演化结果进行比较,符合运行需求的结果直接进入下一轮迭代,不符合运行需求的结果,根据不同设备、部件特征及其允许的特征阈值偏差进行系统调整与自修正;
36)设定疲劳点判定:设定数字孪生设备驱动模型利用疲劳分析与算法优化实现机械臂疲劳点判定。
所述数字孪生设备驱动模型的修正包括以下步骤:
41)获取指标标准状态库,指标标准状态库中存放所采集信号的模糊症状、异常信息、故障规则、诊断专家初步的诊断逻辑;
42)预处理后的现场实时数据在导入数字孪生模型后形成状态子集S={N1,N2,N3,…,Nm},与标准库标准状态库子集P={T1,T2,T3,…,Tn}之间建立模糊推理,
构建状态子集S的归属集合SX=(X1/N1,X2/N2,X3/N3,…,Xm/Nm},根据X的拟合程度,判定状态的接近程度,其中Xi(i=1,2..,m,)表示在状态表征过程中第i类采集数据值,Ni(i=1,2..,m,)表示在状态表征过程中第i类状态数据值,Xi/Ni(i=1,2..,m,)表示在状态表征过程中第i类采集数据对第i类状态数据的呈现度,SXn(n=1,2,..m)=0表示不接近,SXn(n=1,2,..m)=1表示完全接近,S与P子集之间通过置信度Cm建立置信区间与偏差因子;
自修正模糊推理表示为:
Figure BDA0002405835630000051
所述设定疲劳点判定包括以下步骤:
51)将数字孪生模型中累积触发的工作臂运行状态,根据运行监控时点的需要,利用MCD软件结合UG或ANSYS建立工业涂胶机器人关键易受损部件的热力学分析模型计算应力场,或者持续运行状态下的设备作业演化分析云图,获得关键时点的疲劳点状态图;
52)在数字孪生驱动模型中,建立算法优化模块,采集并保存不同关键时点的数据特征,输入算法优化模块;模型采用神经网络Viterbi算法,将由数字孪生模型内演化训练得到的疲劳值进行聚类分析与疲劳状态评估,获得疲劳点数值精准分类,使疲劳评估结果按照所具有的最大可能性的疲劳状态序列进行排列分布,分布情况采用疲劳评价指标r(α,i)折合的0-1的占比值进行直观展示;
53)将Viterbi算法推测得到的最可能的状态序列与常规疲劳经验分析专家库的状态序列进行对比分析,得出疲劳评价指标r(α,i)用于疲劳指标的量化评价与准确率评估,其表达式如下:
r(α,i)=(recentstates(α,i)-states(α,i))/states(α,i)
进而获得疲劳指标量化评分,其表达式如下:
Figure BDA0002405835630000061
式中:i为疲劳要素的权重个数,α为疲劳要素的代号,recentstates(α,i)为当前预测值,states(α,i)疲劳经验值,Score(α,i)为疲劳指标量化评分;
Figure BDA0002405835630000062
式中:m1为健康值下限评分,m1为健康值上限评分;
当Score(α,i)≥m1时,系统输出等级为三级,显示为健康;当m2<Score(α,i)≤m1时,输出等级为二级,显示为疲劳;当Score(α,i)≤m2时,系统输出等级为一级,显示为极度疲劳。
有益效果
本发明的一种基于数字孪生技术的工业涂胶机器人运行状态监测方法,与现有技术相比通过利用数字孪生驱动的设备运行模型,实现了工业涂胶机器人运行状态的实时评价,从而获得更为精准的操作建议、运维建议,实现设备运行的生命周期有效管理。
本发明实现设备在运行工况生命周期内的高效管控,实现产品生产的自适应与动态响应;利用数字孪生思想,实现设备与模型的虚实交互,预测设备在生产运行阶段生命周期内的疲劳状态、运维状态,能够为设备的稳定运行与高效生产提供巨大保障。
本发明借助数字孪生技术探索耦合因素与虚实传递、数据融合、干扰评判、状态演化之间的准确性关联规律,获得状态监测的最优评价方法,有效验证数字孪生体演化分析的准确性,为真实复杂工况下,工业涂胶机器人作业臂的安全稳定运行与高效使用,提供有效增益。
附图说明
图1为本发明的方法顺序图;
图2为本发明所涉及的模块连接原理框图;
图3为本发明的方法原理图。
具体实施方式
为使对本发明的结构特征及所达成的功效有更进一步的了解与认识,用以较佳的实施例及附图配合详细的说明,说明如下:
本发明所应用的工业涂胶机器人包括往复运动关节组件,往复运动关节组件包括作业臂、抓手、摄像头、喷胶装置等组件。往复运动关节组件的关联部位安装有用于获取位置信息三维数值的传感器,往复运动关节组件的作业臂前端安装有姿态传感器,图像采集装置(摄像头、高速拍照相机)的图像采集范围位于作业臂前端,位差监测终端(误差分析仪、红外观测仪)的监测范围位于作业臂前端,在实际应用时,可以根据工业涂胶机器人的具体型号根据传统分析在相应位置安装传感或视频摄像装置。
如图1、图2和图3所示,本发明所述的一种基于数字孪生技术的工业涂胶机器人运行状态监测方法包括以下步骤:
第一步,获取传感器所采集的工业涂胶机器人运行状态数据。
在此运行状态数据可以包括设备位置三维数值、姿态图片、数值、信号,反映设备环境状态的温度、噪声数值信息,反映设备运行状态的数值信息,反映设备生产运行状态的质量信息等。
第二步,工业涂胶机器人运行状态数据的预处理。对工业涂胶机器人的运行状态数据进行初步比对,剔除干扰数据。针对工业涂胶机器人运行状态数据的预处理使用经验模式分析与辅助机械报警相结合的处理方法,其具体步骤如下:
(1)获取利用经验判断建立的专家库,输入数据符合专家库规则的判定为稳定数据,将辅助报警装置提供的故障报警规则作为排他判定,接近排他判定的运行数据,作为干扰数据进行剥离;
(2)将经过预处理的数据进行多源数据融合:使用带有权重分析的特征融合法对影响模拟性能的关键数据进行权重分析,并抽取不同来源数据的数据特征开展模型训练,对不同结果进行合并,使得数字孪生模型中的数据特征能够准确反映工业机器人作业臂的实际工作状态用于下一步的状态预测,其表达式如下:
Figure BDA0002405835630000081
式中:k为影响因素的权重个数,u为影响因素代号,reck(u)为各项权重值,rec(u)为综合权重分析值。
第三步,构建数字孪生设备驱动模型:根据工业涂胶机器人的分析需要构建数字孪生设备驱动模型。
产品数字孪生体将促进建模、仿真与优化技术无缝集成到产品全生命周期中的各个阶段,有效呈现机器人作业臂的运行演化趋势,获得在任何必要工况点的状态特征,有效降低故障率和故障的现实表现。
尽管数字孪生在产品开发与状态管理方面取的了一定成效,但针对复杂环境,耦合因素影响下的状态监测模拟及预测准确性研究依然是一大难题,主要表现在:1)针对工业涂胶机器人特殊运行条件和工况,耦合因素影响下的数据信息提取与数据融合准确性难以准确评价;2)故障状态捕捉,虚实自主交互方面,由于实际干扰因素的影响,难以获得精准判断;3)受耦合因素综合影响,单一流程评定,难以实现作业臂疲劳状态判定与运行状态演化的准确性关联。
因此,本发明在数据采集、提取、筛选、分析、算法优化、疲劳判断等诸多方面开展针对性方法优化,有效提高仿真技术应用于工业涂胶机器人实际作业需求的准确率。其具体步骤如下:
(1)利用三维模型构建数字孪生设备驱动模型,实现通过逻辑控制构建出基于数字孪生的工业机器人作业臂运行工作特征。
(2)利用虚拟引擎开发软件,开展针对三维模型的深度个性化定制,使其准确反映工业涂胶机器人作业臂的运行工况特性;自定义三维模型及适用于机器人产线工业作业水平的工艺流程;建立适用于工业涂胶机器人作业特征的控制逻辑规划,利用计算机语言及虚拟引擎的控制函数,完成逻辑控制触发,实现虚拟引擎软件开发环境下的适用于工业涂胶机器人真实产线设备的数字孪生体,用于运行特征的实时展现与关键部位、关键时间节点的个性化信息定制与评估。
目前特征构建主要用于3d游戏软件,在此将其应用于工业机器人作业臂运行,通过逻辑驱动实现后,可以模拟出工业机器人作业臂运行过程、工作特征。使得企业可利用虚拟环境有效呈现机器人作业臂的运行演化趋势,获得在任何必要工况点的状态特征,有效降低故障率和故障的现实表现,避免了故障带来的经济与能耗损失,便于企业开展早期故障排查与监测服务。
将预处理后的数据信息,输入构建的数字孪生设备驱动模型当中,利用计算机语言及虚拟引擎的控制函数,调整逻辑控制,完成满足工业涂胶机器人实际个性化工作需求的点对点位置控制,实现实时数据输入与响应运行状态的耦合同步。
通过在数字孪生虚拟开发环境中添加虚拟传感器,获得虚拟操作环境的工艺状态特征,使得数字孪生驱动系统具备在任意目标状态时点、任意目标操作位置实现工业涂胶机器人实时状态监控与个性化信息提取的目的。
再利用自定义控制函数,将预处理后的数据信息,在数字孪生设备驱动模型中建立数据偏差分析,结合耦合权重获得随记录时间变化的点对点偏差累积函数η(t),用于远期偏差预测。
(3)初步完善的数字孪生设备驱动模型,结合已采集信号的已知偏差,利用工业涂胶机器人实际运行工况的运行经验组建经验库,建立基于已知偏差的置信区间。
利用自定义控制函数,根据实际工业涂胶机器人的工艺控制需要,选择相应的控制手段或编译自定义控制特征,采用记忆累加方式对偏差进行累加触发并满足△S(uii)=△S(ui0)+η(t),使得数字孪生驱动模型具备适用于工业涂胶机器人个性化运行需求的随时间变化的远期演化分析能力;
其中△S(u)为系统触发故障偏差,△S(uii)为不同采集特征的末次触发偏差,△S(ui0)为新触发前的原始偏差,u为影响因素代号,η(t)为随时间变化偏差累积函数。
使用偏差累加触发结合系统自修正,有助于准确反映和修正设备故障及疲劳累加的发展趋势,提高系统触发准确性。
(4)设定构建数字孪生设备驱动模型利用持续采集的设备运行特征与演化结果进行比较,预测过程可以采用基于神经网络的回归分析法,采用时间序列的实时递推演化思路,将每步演化结果经过递推演化自动进行自迭代,将迭代结果结合不确定性迭代值进行耦合输入后进入下一轮,系统将符合运行需求的结果直接进入下一轮迭代,不符合运行需求的结果,根据不同设备、部件特征及其允许的特征阈值偏差进行系统调整与自修正。通过算法优化,有效避免了单一迭代造成的理论偏差物理失真等问题,有效提高系统预测准确性。
(5)设定疲劳点判定:设定数字孪生设备驱动模型利用疲劳分析与算法优化实现机械臂疲劳点判定。
其具体步骤如下:
A1)将数字孪生模型中累积触发的工作臂运行状态,根据运行监控时点的需要,利用MCD软件结合UG或ANSYS建立工业涂胶机器人关键易受损部件的热力学分析模型计算应力场,或者持续运行状态下的设备作业演化分析云图,获得关键时点的疲劳点状态图。
A2)在数字孪生驱动模型中,建立算法优化模块,采集并保存不同关键时点的数据特征,输入算法优化模块,用于疲劳分析值的有效分类与疲劳点判定的准确划分,大幅降低常规采用的简单数值分类对系统疲劳判定的错误率;模型采用神经网络Viterbi算法,将由数字孪生模型内演化训练得到的疲劳值进行聚类分析与疲劳状态评估,获得疲劳点数值精准分类,使疲劳评估结果按照所具有的最大可能性的疲劳状态序列进行排列分布,分布情况采用疲劳评价指标r(α,i)折合的0-1的占比值进行直观展示,结果表明基于数字孪生仿真的改进算法优于传统故障预测算法。
A3)将Viterbi算法推测得到的最可能的状态序列与常规疲劳经验分析专家库的状态序列进行对比分析,得出疲劳评价指标r(α,i)用于疲劳指标的量化评价与准确率评估,其表达式如下:
r(α,i)=(recentstates(α,i)-states(α,i))/states(α,i)
进而获得疲劳指标量化评分,其表达式如下:
Figure BDA0002405835630000111
式中:i为疲劳要素的权重个数,α为疲劳要素的代号,recentstates(α,i)为当前预测值,states(α,i)疲劳经验值,Score(α,i)为疲劳指标量化评分;
Figure BDA0002405835630000112
式中:m1为健康值下限评分,m1为健康值上限评分;
当Score(α,i)≥m1时,系统输出等级为三级,显示为健康;当m2<Score(α,i)≤m1时,输出等级为二级,显示为疲劳;当Score(α,i)≤m2时,系统输出等级为一级,显示为极度疲劳。
第四步,设备运行状态的初期评价。
将预处理后的工业涂胶机器人运行状态数据输入数字孪生设备驱动模型,数字孪生设备驱动模型以虚拟方式实现工业涂胶机器人作业,得到工业涂胶机器人随运行时间变化而演化的运行状态数值。
第五步,数字孪生设备驱动模型的修正。采用基于数据驱动的自修正优化算法,将常规基于误差的比较过程变更为基于算法的修正优化过程,将数字孪生设备驱动模型演化后的运行状态数值与实时传输的设备状态进行对比匹配,实现数字孪生设备驱动模型的修正。
实际的算法优化过程,将实时传输的状态数据与初期标准库进行模糊映射,获得运行状态初步判断,使得高误差数据被优先排畸;随后将基于数字孪生模型的运行状态演化数值与处理后的末一期运行数值进行比较,构成拟合子集SX,根据拟合子集的数据接近呈现度,在有效置信区间内,确认偏差因子;偏差因子用于模型自修复,并计入累加触发的递归演化,实现作业臂状态的准确预测演化。其具体步骤如下:
(1)获取指标标准状态库,指标标准状态库中存放所采集信号的模糊症状、异常信息、故障规则、诊断专家初步的诊断逻辑。
(2)预处理后的现场实时数据在导入数字孪生模型后形成状态子集S={N1,N2,N3,…,Nm},与标准库标准状态库子集P={T1,T2,T3,…,Tn}之间建立模糊推理,
构建状态子集S的归属集合SX=(X1/N1,X2/N2,X3/N3,…,Xm/Nm},根据X的拟合程度,判定状态的接近程度,其中Xi(i=1,2..,m,)表示在状态表征过程中第i类采集数据值,Ni(i=1,2..,m,)表示在状态表征过程中第i类状态数据值,Xi/Ni(i=1,2..,m,)表示在状态表征过程中第i类采集数据对第i类状态数据的呈现度,SXn(n=1,2,..m)=0表示不接近,SXn(n=1,2,..m)=1表示完全接近,S与P子集之间通过置信度Cm建立置信区间与偏差因子;
自修正模糊推理表示为:
Figure BDA0002405835630000121
指标标准状态库利用预处理后的现场实时数据以及模型演化的计算数据,实现标准库的自我修正与完善,修正后的状态库为数字孪生模型提供反馈修正,以提高数字模型准确率;数字孪生设备驱动模型,利用已采集信号的已知偏差,结合状态库的修正反馈,采用记忆累加方式对偏差进行累加触发,帮助数字孪生设备驱动模型开展模型训练与演化预测。综合修正后的数字孪生模型得出疲劳评测分级指标,用于工业涂装机器人手臂的远期运行评价与运维建议。
第六步,工业涂胶机器人运行状态的监测。
获取传感器采集的工业涂胶机器人运行状态实时数据,将实时数据预处理后输入修正后的数字孪生设备驱动模型,预测出工业涂胶机器人运行的远期趋势,判定出运行疲劳点。
本发明相关算法曾应用于多变环境的发动机状态监控与寿命预测中,利用提出的疲劳分析模型与优化算法能够准确预测疲劳寿命的分布,预测结果与实验结果具有较好的吻合性。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是本发明的原理,在不脱离本发明精神和范围的前提下本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明的范围内。本发明要求的保护范围由所附的权利要求书及其等同物界定。

Claims (3)

1.一种基于数字孪生技术的工业涂胶机器人运行状态监测方法,所述的工业涂胶机器人包括往复运动关节组件,往复运动关节组件的关联部位安装有用于获取位置信息三维数值的传感器,往复运动关节组件的作业臂前端安装有姿态传感器,图像采集装置的图像采集范围位于作业臂前端,位差监测终端的监测范围位于作业臂前端,其特征在于,
所述的工业涂胶机器人运行状态监测方法包括以下步骤:
11)获取传感器所采集的工业涂胶机器人运行状态数据:所述的运行状态数据包括设备位置三维数值、姿态图片、数值、信号,反映设备环境状态的温度、噪声数值信息,反映设备运行状态的数值信息,反映设备生产运行状态的质量信息;
12)工业涂胶机器人运行状态数据的预处理:对工业涂胶机器人的运行状态数据进行初步比对,剔除干扰数据;
13)构建数字孪生设备驱动模型:根据工业涂胶机器人的分析需要构建数字孪生设备驱动模型;所述构建数字孪生设备驱动模型包括以下步骤:
131)利用三维模型构建数字孪生设备驱动模型,实现通过逻辑控制构建出基于数字孪生的工业机器人作业臂运行工作特征;
132)利用虚拟引擎开发软件,开展针对三维模型的深度个性化定制,使其准确反映工业涂胶机器人作业臂的运行工况特性;自定义三维模型及适用于机器人产线工业作业水平的工艺流程;建立适用于工业涂胶机器人作业特征的控制逻辑规划,利用计算机语言及虚拟引擎的控制函数,完成逻辑控制触发,实现虚拟引擎软件开发环境下的适用于工业涂胶机器人真实产线设备的数字孪生体,用于运行特征的实时展现与关键部位、关键时间节点的个性化信息定制与评估;
133)将预处理后的数据信息,输入构建的数字孪生设备驱动模型当中,利用计算机语言及虚拟引擎的控制函数,调整逻辑控制,完成满足工业涂胶机器人实际个性化工作需求的点对点位置控制,实现实时数据输入与响应运行状态的耦合同步;通过在数字孪生虚拟开发环境中添加虚拟传感器,获得虚拟操作环境的工艺状态特征,使得数字孪生驱动系统具备在任意目标状态时点、任意目标操作位置实现工业涂胶机器人实时状态监控与个性化信息提取的目的;利用自定义控制函数,将预处理后的数据信息,在数字孪生设备驱动模型中建立数据偏差分析,结合耦合权重获得随记录时间变化的点对点偏差累积函数η(t),用于远期偏差预测;
134)初步完善的数字孪生设备驱动模型,结合已采集信号的已知偏差,利用工业涂胶机器人实际运行工况的运行经验组建经验库,建立基于已知偏差的置信区间;利用自定义控制函数,根据实际工业涂胶机器人的工艺控制需要,选择相应的控制手段或编译自定义控制特征,采用记忆累加方式对偏差进行累加触发并满足△S(uii)=△S(ui0)+η(t),使得数字孪生驱动模型具备适用于工业涂胶机器人个性化运行需求的随时间变化的远期演化分析能力;
其中△S(u)为系统触发故障偏差,△S(uii)为不同采集特征的末次触发偏差,△S(ui0)为新触发前的原始偏差,u为影响因素代号,η(t)为随时间变化偏差累积函数;
135)设定构建数字孪生设备驱动模型利用持续采集的设备运行特征与演化结果进行比较,符合运行需求的结果直接进入下一轮迭代,不符合运行需求的结果,根据不同设备、部件特征及其允许的特征阈值偏差进行系统调整与自修正;
136)设定疲劳点判定:设定数字孪生设备驱动模型利用疲劳分析与算法优化实现机械臂疲劳点判定;所述设定疲劳点判定包括以下步骤:
1361)将数字孪生模型中累积触发的工作臂运行状态,根据运行监控时点的需要,利用MCD软件结合UG或ANSYS建立工业涂胶机器人关键易受损部件的热力学分析模型计算应力场,或者持续运行状态下的设备作业演化分析云图,获得关键时点的疲劳点状态图;
1362)在数字孪生驱动模型中,建立算法优化模块,采集并保存不同关键时点的数据特征,输入算法优化模块;模型采用神经网络Viterbi算法,将由数字孪生模型内演化训练得到的疲劳值进行聚类分析与疲劳状态评估,获得疲劳点数值精准分类,使疲劳评估结果按照所具有的最大可能性的疲劳状态序列进行排列分布,分布情况采用疲劳评价指标r(α,i)折合的0-1的占比值进行直观展示;
1363)将Viterbi算法推测得到的最可能的状态序列与常规疲劳经验分析专家库的状态序列进行对比分析,得出疲劳评价指标r(α,i)用于疲劳指标的量化评价与准确率评估,其表达式如下:
r(α,i)=(recentstates(α,i)-states(α,i))/states(α,i)
进而获得疲劳指标量化评分,其表达式如下:
Figure FDA0002814392190000031
式中:i为疲劳要素的权重个数,α为疲劳要素的代号,recentstates(α,i)为当前预测值,states(α,i)疲劳经验值,Score(α,i)为疲劳指标量化评分;
Figure FDA0002814392190000032
式中:m1为健康值下限评分,m1为健康值上限评分;
当Score(α,i)≥m1时,系统输出等级为三级,显示为健康;当m2<Score(α,i)≤m1时,输出等级为二级,显示为疲劳;当Score(α,i)≤m2时,系统输出等级为一级,显示为极度疲劳;
14)设备运行状态的初期评价:将预处理后的工业涂胶机器人运行状态数据输入数字孪生设备驱动模型,数字孪生设备驱动模型以虚拟方式实现工业涂胶机器人作业,得到工业涂胶机器人随运行时间变化而演化的运行状态数值;
15)数字孪生设备驱动模型的修正:将数字孪生设备驱动模型演化后的运行状态数值与实时传输的设备状态进行对比匹配,实现数字孪生设备驱动模型的修正;
16)工业涂胶机器人运行状态的监测:获取传感器采集的工业涂胶机器人运行状态实时数据,将实时数据预处理后输入修正后的数字孪生设备驱动模型,预测出工业涂胶机器人运行的远期趋势,判定出运行疲劳点。
2.根据权利要求1所述的一种基于数字孪生技术的工业涂胶机器人运行状态监测方法,其特征在于,所述工业涂胶机器人运行状态数据的预处理使用经验模式分析与辅助机械报警相结合的处理方法,其具体步骤如下:
21)获取利用经验判断建立的专家库,输入数据符合专家库规则的判定为稳定数据,将辅助报警装置提供的故障报警规则作为排他判定,接近排他判定的运行数据,作为干扰数据进行剥离;
22)将经过预处理的数据进行多源数据融合:使用带有权重分析的特征融合法对影响模拟性能的关键数据进行权重分析,并抽取不同来源数据的数据特征开展模型训练,对不同结果进行合并,使得数字孪生模型中的数据特征能够准确反映工业机器人作业臂的实际工作状态用于下一步的状态预测,其表达式如下:
Figure FDA0002814392190000051
式中:k为影响因素的权重个数,u为影响因素代号,reck(u)为各项权重值,rec(u)为综合权重分析值。
3.根据权利要求1所述的一种基于数字孪生技术的工业涂胶机器人运行状态监测方法,其特征在于,所述数字孪生设备驱动模型的修正包括以下步骤:
31)获取指标标准状态库,指标标准状态库中存放所采集信号的模糊症状、异常信息、故障规则、诊断专家初步的诊断逻辑;
32)预处理后的现场实时数据在导入数字孪生模型后形成状态子集S={N1,N2,N3,…,Nm},与标准库标准状态库子集P={T1,T2,T3,…,Tn}之间建立模糊推理,
构建状态子集S的归属集合SX=(X1/N1,X2/N2,X3/N3,…,Xm/Nm},根据X的拟合程度,判定状态的接近程度,其中Xi(i=1,2..,m,)表示在状态表征过程中第i类采集数据值,Ni(i=1,2..,m,)表示在状态表征过程中第i类状态数据值,Xi/Ni(i=1,2..,m,)表示在状态表征过程中第i类采集数据对第i类状态数据的呈现度,SXn(n=1,2,..m)=0表示不接近,SXn(n=1,2,..m)=1表示完全接近,S与P子集之间通过置信度Cm建立置信区间与偏差因子;
自修正模糊推理表示为:
Figure FDA0002814392190000052
CN202010161152.5A 2020-03-10 2020-03-10 一种基于数字孪生技术的工业涂胶机器人运行状态监测方法 Active CN111230887B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010161152.5A CN111230887B (zh) 2020-03-10 2020-03-10 一种基于数字孪生技术的工业涂胶机器人运行状态监测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010161152.5A CN111230887B (zh) 2020-03-10 2020-03-10 一种基于数字孪生技术的工业涂胶机器人运行状态监测方法

Publications (2)

Publication Number Publication Date
CN111230887A CN111230887A (zh) 2020-06-05
CN111230887B true CN111230887B (zh) 2021-01-26

Family

ID=70868387

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010161152.5A Active CN111230887B (zh) 2020-03-10 2020-03-10 一种基于数字孪生技术的工业涂胶机器人运行状态监测方法

Country Status (1)

Country Link
CN (1) CN111230887B (zh)

Families Citing this family (41)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111240209B (zh) * 2020-03-16 2020-10-09 广东工业大学 构型动型控型优型联动响应的自适应组态方法及系统
CN111680749B (zh) * 2020-06-08 2023-11-07 北京百度网讯科技有限公司 用于获取点胶机输出结果的方法及装置
CN111964575B (zh) * 2020-07-06 2022-02-01 北京卫星制造厂有限公司 一种面向移动机器人铣削加工的数字孪生建模方法
CN111950133B (zh) * 2020-07-24 2022-04-26 华东交通大学 一种基于数字孪生的发动机可靠寿命预测方法
CN112085261B (zh) * 2020-08-19 2022-06-21 浙江工业大学 基于云端融合和数字孪生技术的企业生产现状诊断方法
CN112149969B (zh) * 2020-09-08 2023-06-27 国网湖南省电力有限公司 特高压直流控制保护合并单元运行状态评价方法及系统
CN112231966B (zh) * 2020-09-08 2023-04-07 合肥学院 一种基于数字孪生的协作机器人可装配性预测系统及方法
CN112085393A (zh) * 2020-09-11 2020-12-15 青岛恒拓环保科技有限公司 一种用于工业生产的大数据交互方法
CN112370035A (zh) * 2020-10-15 2021-02-19 同济大学 一种基于数字孪生平台的人机协作疲劳检测系统
CN112348251B (zh) * 2020-11-05 2024-02-09 傲林科技有限公司 一种决策辅助方法、装置、电子设备及存储介质
CN112346440B (zh) * 2020-11-20 2022-02-11 深圳优地科技有限公司 机器人健康监控方法、装置、设备及可读存储介质
CN112634203A (zh) * 2020-12-02 2021-04-09 富泰华精密电子(郑州)有限公司 图像检测方法、电子设备及计算机可读存储介质
CN112560263B (zh) * 2020-12-11 2023-02-03 太原理工大学 一种基于数字孪生的移动机器人状态监控与保养系统
CN112731894A (zh) * 2021-02-02 2021-04-30 深圳华龙讯达信息技术股份有限公司 一种基于数字孪生技术的包装机智能控制系统及方法
CN112947173B (zh) * 2021-02-03 2023-06-02 北京理工大学 一种数字孪生车间运行状态的预测方法、控制器及系统
CN113033055B (zh) * 2021-03-29 2022-07-01 武汉理工大学 一种基于数字孪生的船用发动机状态评估方法及系统
CN113139659A (zh) * 2021-04-09 2021-07-20 智科云创(北京)科技有限公司 基于数字孪生的水利监测方法及系统
CN113259486B (zh) * 2021-06-24 2021-09-17 国网天津市电力公司营销服务中心 一种基于数字孪生的计量设备自动化检定线运维系统
CN113534679B (zh) * 2021-07-06 2023-04-28 上海新氦类脑智能科技有限公司 系统监测模型生成方法、处理器芯片以及工业系统
CN113485295A (zh) * 2021-07-07 2021-10-08 西北工业大学 基于数字孪生的四足机器人故障预测方法、装置及设备
CN113400353B (zh) * 2021-07-07 2024-02-02 西北工业大学 基于数字孪生的多指灵巧手状态监测方法、装置及设备
CN113408764B (zh) * 2021-08-17 2023-01-13 能科科技股份有限公司 基于数字孪生的产品在线技术状态管控系统
CN114136670A (zh) * 2021-10-26 2022-03-04 中国石油化工股份有限公司 基于管道检测机器人测试平台的管道检测机器人评价方法
CN113902222B (zh) * 2021-11-17 2022-04-15 北京协合运维风电技术有限公司 一种新能源智慧运维系统
CN114323725A (zh) * 2022-01-26 2022-04-12 深圳市信润富联数字科技有限公司 点胶机健康度检测方法、装置、设备及存储介质
CN114663060B (zh) * 2022-03-25 2022-11-18 北京清大宏源环保科技有限公司 一种基于数字孪生技术的产品制造生产线协同智能管理系统
CN114770499B (zh) * 2022-03-30 2023-09-12 清华大学 工业机器人的能耗高效建模预测方法及装置
CN114529547B (zh) * 2022-04-24 2022-07-12 常州捷仕特机器人科技有限公司 一种机器人涂胶工作站的涂胶质量检测系统及方法
CN114918976B (zh) * 2022-06-16 2022-12-02 慧之安信息技术股份有限公司 基于数字孪生技术的关节机器人健康状态评估方法
CN114918927A (zh) * 2022-07-22 2022-08-19 南京励业智能科技有限公司 一种工业数字孪生的仿真运维系统及方法
CN115577488A (zh) * 2022-08-17 2023-01-06 国网山东省电力公司营销服务中心(计量中心) 一种关口电能计量装置性能诊断系统
CN115115289B (zh) * 2022-08-29 2022-11-11 合肥市满好科技有限公司 一种基于xr技术的交互监控管理系统
US20240078442A1 (en) * 2022-09-07 2024-03-07 International Business Machines Corporation Self-development of resources in multi-machine environment
CN115329446B (zh) * 2022-10-13 2023-01-31 江苏航运职业技术学院 一种装配式建筑预制构件智能吊装过程数字孪生建模方法
CN116117811A (zh) * 2023-02-07 2023-05-16 北京机科国创轻量化科学研究院有限公司 一种基于数字孪生的工业机器人综合性能分析方法及系统
CN116071053B (zh) * 2023-04-07 2023-11-28 沃德传动(天津)股份有限公司 基于数字孪生的往复式压缩机状态预警系统
CN116168164B (zh) * 2023-04-25 2023-07-04 深圳墨影科技有限公司 基于机器人产业链的数字孪生可视化方法
CN116502182B (zh) * 2023-06-25 2023-10-27 金锐同创(北京)科技股份有限公司 基于测绘数据的孪生平台构建方法、装置、设备和介质
CN116720415B (zh) * 2023-08-09 2023-12-05 中国人民解放军火箭军工程大学 一种基于数字孪生的目标红外特性建模方法
CN117140485A (zh) * 2023-09-21 2023-12-01 苏州诺克汽车工程装备有限公司 一种基于数字孪生的协作机器人示教装配方法与系统
CN117670864A (zh) * 2023-12-28 2024-03-08 北汽利戴工业技术服务(北京)有限公司 基于工业相机的图像识别系统

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108919760A (zh) * 2018-07-05 2018-11-30 长安大学 一种基于数字孪生的智能车间自治生产过程动态联动控制方法
CN109571476A (zh) * 2018-12-14 2019-04-05 南京理工大学 工业机器人数字孪生实时作业控制、监控与精度补偿方法
CN109719730A (zh) * 2019-01-25 2019-05-07 温州大学 一种断路器柔性化装配过程的数字孪生机器人
KR101989982B1 (ko) * 2019-03-08 2019-06-17 이봉규 디지털 트윈 기반의 실내 공간 분석을 통한 모델링 시스템과 모델링 방법
CN110181519A (zh) * 2019-06-25 2019-08-30 广东希睿数字科技有限公司 基于数字孪生机器人的地铁台站门故障检测方法及系统
CN110399642A (zh) * 2019-06-21 2019-11-01 浙江大学 一种针对生产流水线的数字孪生体及其构建方法和应用
CN110531846A (zh) * 2018-05-24 2019-12-03 明日基金知识产权控股有限公司 在实时3d虚拟世界代表真实世界的范围内的实时3d虚拟对象的双向实时3d交互操作
CN110605709A (zh) * 2019-09-25 2019-12-24 西南交通大学 数字孪生与精密滤波驱动机器人集成系统及其使用方法
WO2020040772A1 (en) * 2018-08-23 2020-02-27 Siemens Aktiengesellschaft Automatic setup of digital twins for industrial controllers
CN110850818A (zh) * 2019-10-21 2020-02-28 温州大学 一种实现断路器自动化装配车间的数字孪生建模的方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110531846A (zh) * 2018-05-24 2019-12-03 明日基金知识产权控股有限公司 在实时3d虚拟世界代表真实世界的范围内的实时3d虚拟对象的双向实时3d交互操作
CN108919760A (zh) * 2018-07-05 2018-11-30 长安大学 一种基于数字孪生的智能车间自治生产过程动态联动控制方法
WO2020040772A1 (en) * 2018-08-23 2020-02-27 Siemens Aktiengesellschaft Automatic setup of digital twins for industrial controllers
CN109571476A (zh) * 2018-12-14 2019-04-05 南京理工大学 工业机器人数字孪生实时作业控制、监控与精度补偿方法
CN109719730A (zh) * 2019-01-25 2019-05-07 温州大学 一种断路器柔性化装配过程的数字孪生机器人
KR101989982B1 (ko) * 2019-03-08 2019-06-17 이봉규 디지털 트윈 기반의 실내 공간 분석을 통한 모델링 시스템과 모델링 방법
CN110399642A (zh) * 2019-06-21 2019-11-01 浙江大学 一种针对生产流水线的数字孪生体及其构建方法和应用
CN110181519A (zh) * 2019-06-25 2019-08-30 广东希睿数字科技有限公司 基于数字孪生机器人的地铁台站门故障检测方法及系统
CN110605709A (zh) * 2019-09-25 2019-12-24 西南交通大学 数字孪生与精密滤波驱动机器人集成系统及其使用方法
CN110850818A (zh) * 2019-10-21 2020-02-28 温州大学 一种实现断路器自动化装配车间的数字孪生建模的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于数字孪生仿真建模的机器人状态监测技术研究;耿琦琦;《中国优秀硕士学位论文全文数据库 信息科技辑》;20200115;全文 *
车间生产过程数字孪生系统构建及应用;柳林燕等;《计算机集成制造系统》;20190630;全文 *

Also Published As

Publication number Publication date
CN111230887A (zh) 2020-06-05

Similar Documents

Publication Publication Date Title
CN111230887B (zh) 一种基于数字孪生技术的工业涂胶机器人运行状态监测方法
CN111539583B (zh) 一种基于数字孪生的生产过程仿真优化方法
CN111496781B (zh) 数字孪生驱动的机械臂建模、控制与监测一体化集成系统
CN113887016A (zh) 一种基于数字孪生的船舶数字化车间仿真方法与系统
CN112085261B (zh) 基于云端融合和数字孪生技术的企业生产现状诊断方法
CN110571792A (zh) 一种电网调控系统运行状态的分析评估方法及系统
CN112070279A (zh) 一种基于数字孪生技术的产品加工控制方法
CN110757510B (zh) 一种机器人剩余寿命预测方法及系统
CN107644297B (zh) 一种电机系统节能量计算和验证方法
CN111638707B (zh) 基于som聚类和mpca的间歇过程故障监测方法
CN113695109B (zh) 一种用于自动化喷涂过程的程序匹配和迭代优化方法
CN114118673A (zh) 一种基于数字孪生技术的车间智能化故障诊断预警方法
CN110750455B (zh) 基于系统日志分析的智能在线自更新故障诊断方法和系统
CN111190429A (zh) 一种基于强化学习的无人机主动容错控制方法
CN114330026A (zh) 一种数字孪生系统仿真方法及装置
CN115438726A (zh) 一种基于数字孪生技术的设备寿命与故障类型预测方法及系统
CN112231966B (zh) 一种基于数字孪生的协作机器人可装配性预测系统及方法
CN111765075B (zh) 一种锻造液压机泵源故障预测方法及系统
CN117391674A (zh) 一种基于可靠性的电气设备预防性维修优化方法及装置
CN112732541A (zh) 一种用于复杂装备故障诊断的智能判据挖掘系统
CN116663419A (zh) 一种基于优化的Elman神经网络的无传感器设备故障预测方法
CN109800866B (zh) 一种基于GA-Elman神经网络的可靠性增长预测方法
CN114239938A (zh) 一种基于状态的能源数字孪生体的构建方法
CN115952888A (zh) 基于多变量灰色模型的能源碳排放预测方法及系统
WO2021111936A1 (ja) 予測システム、情報処理装置および情報処理プログラム

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
GR01 Patent grant
GR01 Patent grant