CN110045608A - 基于数字孪生的机械设备零部件结构参数动态优化方法 - Google Patents

基于数字孪生的机械设备零部件结构参数动态优化方法 Download PDF

Info

Publication number
CN110045608A
CN110045608A CN201910261337.0A CN201910261337A CN110045608A CN 110045608 A CN110045608 A CN 110045608A CN 201910261337 A CN201910261337 A CN 201910261337A CN 110045608 A CN110045608 A CN 110045608A
Authority
CN
China
Prior art keywords
mechanical equipment
model
value
parameter
component structural
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.)
Granted
Application number
CN201910261337.0A
Other languages
English (en)
Other versions
CN110045608B (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.)
Taiyuan University of Technology
Original Assignee
Taiyuan University 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 Taiyuan University of Technology filed Critical Taiyuan University of Technology
Priority to CN201910261337.0A priority Critical patent/CN110045608B/zh
Publication of CN110045608A publication Critical patent/CN110045608A/zh
Application granted granted Critical
Publication of CN110045608B publication Critical patent/CN110045608B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance

Abstract

本发明公开了一种基于数字孪生的机械设备零部件结构参数动态优化方法,通过构建高保真模型,在虚拟空间实现物理空间对应设备的数字化镜像,方便后期结构参数修改和超写实仿真;通过进行超写实仿真,在虚拟空间实时动态的同步反映物理空间对应实体的状态,在虚拟空间实现物理设备真实情况的写实运动;同时利用深度学习理论,构建神经网络结构,借助其强大的数字挖掘和映射能力,挖掘建立结构参数和疲劳寿命之间的关系,通过结合高保真模型和超写实仿真环境,实现结构参数的动态优化。通过本发明,实现虚拟空间对物理空间的结构参数的动态优化与反向指导,提高了优化效率和真实性。

Description

基于数字孪生的机械设备零部件结构参数动态优化方法
技术领域
本发明涉及深度学习领域,特别是涉及一种基于数字孪生的机械设备零部件结构参数动态优化方法。
背景技术
数字孪生为虚拟空间和物理空间的实时交互与双向连接提供了重要的理论依据和技术支撑,近几年在理论和应用层面均取得了快速发展,目前数字孪生技术被应用于航天设备和车间生产控制。对于通用机械这种工作环境复杂,条件多变的大型设备,其关键零部件结构参数优化手段落后,基本停留在经验设计、参数静态仿真设计,无法对实际运行情况进行高逼真度模拟仿真,这使得对真实工作面环境中机械设备零部件结构参数的优化准确性差、效率低、智能化程度滞后。在结构参数设计阶段,很少考虑到设备的实际工作环境。而机械设备的实际工作条件多变,存在诸多的不确定因素,因此,对机械设备关键零部件的结构参数设计,更应该实时动态的考虑其实际工作环境的影响。现有机械设备零部件结构参数优化手段落后,优化效率低,真实性有待提升,对于复杂工况条件下的零部件关键结构参数优化困难。
发明内容
本发明的目的是为解决上述现有技术的不足之处而提供一种基于数字孪生的机械设备零部件结构参数动态优化方法。
为解决上述技术问题,本发明采用的一个技术方案是:提供一种基于数字孪生的机械设备零部件结构参数动态优化方法,包括:
构建机械设备全部或部分结构的三维模型,研究三维模型动力学与运动状态的关联性,确定反映运动状态的参数,并根据参数分析确定机械设备全部或部分结构中所需的传感器类型及安装位置,按照分析结果在物理空间物理实体中布置传感器,通过传感器监测数据驱动虚拟空间参数更新,完成高保真建模;
依据物理实体的历史监测数据及物理空间传感器实时更新的监测数据,进行超写实仿真,获取高保真三维模型零件结构参数与对应疲劳寿命数据集作为数据分析样本;
构建深度卷积神经网络模型,用高保真三维模型零件结构参数与对应的疲劳寿命数据样本集训练模型,保存训练好的模型;
将目标疲劳寿命数据作为预测模型的输入,通过传感器实时监测数据动态更新虚拟空间参数和预测模型动态训练,得到动态优化后的零部件结构参数。
其中,在进行高保真建模的步骤中,包括步骤:
依据物理空间机械设备零部件结构的受力大小、力的类型、结构、工作时的姿态,轮廓尺寸、相对位置信息,用三维软件构建机械设备零部件结构的三维模型;
导入ANSYS软件中进行动力学分析,研究动力学与运动状态的关联性,确定反映运动状态的关键参数,其中,反映运动状态的关键参数至少包括应变、振型、节点位移、节点位置;
依据得到的关键参数布置需要的传感器的类型,以及传感器在机械设备零部件结构位置的布置;依据物理空间机械设备零部件结构的材料类型、几何尺寸,利用UG参数化建模在虚拟空间构建机械设备零部件结构的三维模型,实现高保真建模。
其中,所述传感器类型至少包括加速度、温度传感器和力传感器。
其中,传感器布置在机械设备零部件结构上应变大的位置。
其中,通过VR技术获得的机械设备零部件结构物理空间多物理参数,其中,机械设备零部件结构的物理空间多物理参数至少包括受力、速度、振动、温度、环境温度、位置参数。
其中,对机械设备零部件结构进行超写实仿真的步骤包括:
依据有限元参数化技术,通过传感器实时监测数据实时动态更新虚拟空间相关参数设置,实现动态疲劳仿真,得到对应的疲劳寿命值;其中,传感器实时监测机械设备零部件结构的监测数据至少包括力、速度、温度及位置参数;
通过调整虚拟空间高保真模型对应的结构参数,继续进行仿真,得到该结构参数下机械设备零部件对应的疲劳寿命值,构建结构参数与机械设备零部件疲劳寿命的对应关系;
通过调整不同的结构参数,得到不同的疲劳寿命值,产生机械设备零部件结构参数与对应的疲劳寿命数据样本集。
其中,构建深度卷积神经网络模型的步骤包括:
利用python语言实现DCNN模型,用测试集进行DCNN模型的训练,进行参数优化,调整模型的层数、模型的具体结构及卷积核的大小,直至训练集中的预测值和实际值的均方误差达到最小,训练预测结果达到最优确定训练次数、模型层数、模型具体结构以及卷积核大小,保存模型;
其中,均方误差表达式为
N表示参加训练的数据量,ypi表示对第i个输入的预测值,yti表示第i个输入对应的实际值;
基于训练好的模型,利用测试集测试模型预测值,将模型的预测值和实际值做对比,得到模型预测的准确性;
其中,对模型预测结果进行评判时,采用均方根误差(RMSE)、拟合优度检验(2)、平均绝对百分误差(MAPE)和Score_function三个指标进行评判;
表达式分别为:
预测分析中RMSE越接近0,代表预测结果越精确;
R2值越接近1,代表预测结果越好;
平均绝对百分误差(MAPE)不仅仅考虑预测值与真实值之间的误差,还考虑该误差占真实值的比例;
Score值越接近0,代表预测结果越精确;
上述公式中,yi表示实际值,表示预测值,表示预测的均值,yi表示第i 组的寿命对应参数。
其中,在将目标疲劳寿命数据作为训练模型的输入,对高保真三维模型的零件结构参数进行预测的步骤中,包括步骤:
预设机械设备零部件结构的目标寿命值,并将目标寿命值(目标寿命值)作为预测模型的输入,得到预测的机械设备零部件结构参数值;
以预测的机械设备零部件结构参数值作为指导,更新虚拟空间模型的对应的尺寸,通过有限元参数化进行模型的仿真,得到对应的疲劳寿命值;
比较仿真得到的疲劳寿命值与目标寿命值,若两者相等,则预测得到的参数就是优化后的参数;若不相等,进行训练集样本补充,继续训练模型,再将目标寿命值作为输入进行参数预测,依据预测的结构参数更新虚拟空间模型对应的参数,进行新一轮超写实仿真,比较目标寿命值与仿真寿命值;循环上述步骤,直至两个寿命值相等。
区别于现有技术,本发明的基于数字孪生的机械设备零部件结构参数动态优化方法,通过构建高保真模型,在虚拟空间实现物理空间对应设备的数字化镜像,方便后期结构参数修改和超写实仿真;通过进行超写实仿真,在虚拟空间实时动态的同步反映物理空间对应实体的状态,在虚拟空间实现物理设备真实情况的写实运动;同时利用深度学习理论,构建神经网络结构,借助其强大的数字挖掘和映射能力,挖掘建立结构参数和疲劳寿命之间的关系,通过结合高保真模型和超写实仿真环境,实现结构参数的动态优化。通过本发明,实现虚拟空间对物理空间的结构参数的动态优化与反向指导,提高了优化效率和真实性。
附图说明
图1是本发明提供的一种基于数字孪生的机械设备零部件结构参数动态优化方法的流程示意图。
图2是本发明提供的一种基于数字孪生的机械设备零部件结构参数动态优化方法中依据目标寿命得到模型对齿宽参数的动态优化效果图。
图3是本发明提供的一种基于数字孪生的机械设备零部件结构参数动态优化方法中依据目标寿命进行变位系数的动态优化效果图。
图4是本发明提供的一种基于数字孪生的机械设备零部件结构参数动态优化方法中依据目标寿命进行压力角值动态优化效果图。
图5是本发明提供的一种基于数字孪生的机械设备零部件结构参数动态优化方法中依据目标寿命值进行重合度值动态优化效果图。
具体实施方式
在下面的描述中阐述了很多具体细节以便于充分理解本发明。但是本发明能够以很多不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施的限制。
其次,本发明利用示意图进行详细描述,在详述本发明实施例时,为便于说明,示意图只是实例,其在此不应限制本发明保护的范围。
参阅图1,图1是本发明提供的一种基于数字孪生的机械设备零部件结构参数动态优化方法的流程示意图。该方法的步骤包括:
构建机械设备全部或部分结构的三维模型,研究三维模型动力学与运动状态的关联性,确定反映运动状态的参数,并根据参数分析确定机械设备全部或部分结构中所需的传感器类型及安装位置,按照分析结果在物理空间物理实体中布置传感器,通过传感器监测数据驱动虚拟空间参数更新,完成高保真建模;
依据物理实体的历史监测数据及物理空间传感器实时更新的监测数据,进行超写实仿真,获取高保真三维模型零件结构参数与对应疲劳寿命数据集作为数据分析样本;
构建深度卷积神经网络模型,用高保真三维模型零件结构参数与对应的疲劳寿命数据样本集训练模型,保存训练好的模型;
将目标疲劳寿命数据作为预测模型的输入,通过传感器实时监测数据动态更新虚拟空间参数和预测模型动态训练,得到动态优化后的零部件结构参数。
其中,在进行高保真建模的步骤中,包括步骤:
依据物理空间机械设备零部件结构的受力大小、力的类型、结构、工作时的姿态,轮廓尺寸、相对位置信息,用三维软件构建机械设备零部件结构的三维模型;
导入ANSYS软件中进行动力学分析,研究动力学与运动状态的关联性,确定反映运动状态的关键参数,其中,反映运动状态的关键参数至少包括应变、振型、节点位移、节点位置;
依据得到的关键参数布置需要的传感器的类型,以及传感器在机械设备零部件结构位置的布置;依据物理空间机械设备零部件结构的材料类型、几何尺寸,利用UG参数化建模在虚拟空间构建机械设备零部件结构的三维模型,实现高保真建模。
其中,所述传感器类型至少包括加速度、温度传感器和力传感器。
其中,传感器布置在机械设备零部件结构上应变大的位置。
其中,通过VR技术获得的机械设备零部件结构物理空间多物理参数,其中,机械设备零部件结构的物理空间多物理参数至少包括受力、速度、振动、温度、环境温度、位置参数。对矿井设备而言,还包括煤尘的物理参数。
其中,对机械设备零部件结构进行超写实仿真的步骤包括:
依据有限元参数化技术,通过传感器实时监测数据实时动态更新虚拟空间相关参数设置,实现动态疲劳仿真,得到对应的疲劳寿命值;其中,传感器实时监测机械设备零部件结构的监测数据至少包括力、速度、温度及位置参数;
通过调整机械设备零部件结构对应的结构参数,继续进行仿真,得到该参数下机械设备零部件结构对应的疲劳寿命值,构建结构参数与机械设备零部件疲劳寿命的对应关系;
通过调整不同的参数,得到不同的疲劳寿命值,产生机械设备零部件结构参数与对应的疲劳寿命数据样本集。
其中,构建深度卷积神经网络模型的步骤包括:
利用python语言实现DCNN模型,用测试集进行DCNN模型的训练,使模型学习到不同阶段的特征,进行参数优化,调整模型的层数、模型的具体结构及卷积核的大小,直至训练集中的预测值和实际值的均方误差达到最小,训练预测结果达到最优确定训练次数、模型层数、模型具体结构以及卷积核大小,保存模型;
其中,均方误差表达式为
N表示参加训练的数据量,ypi表示对第i个输入的预测值,yti表示第i个输入对应的实际值;
基于训练好的模型,利用测试集测试模型预测值,将模型的预测值和实际值做对比,得到模型预测的准确性;
其中,对模型预测结果进行评判时,采用均方根误差(RMSE)、拟合优度检验(R2)、平均绝对百分误差(MAPE)和Score_function三个指标进行评判;
表达式分别为:
预测分析中RMSE越接近0,代表预测结果越精确;
R2值越接近1,代表预测结果越好;
平均绝对百分误差(MAPE)不仅仅考虑预测值与真实值之间的误差,还考虑该误差占真实值的比例;
Score值越接近0,代表预测结果越精确;
上述公式中,yi表示实际值,表示预测值,表示预测的均值,yi表示第i 组的寿命对应参数。
其中,在将目标疲劳寿命数据作为训练模型的输入,对高保真三维模型的零件结构参数进行预测的步骤中,包括步骤:
预设机械设备零部件结构的目标寿命值,并将目标寿命值(目标寿命值)作为训练模型的输入,得到预测的机械设备零部件结构参数值;
以预测的机械设备零部件结构参数值作为指导,更新虚拟空间模型的对应的尺寸,通过有限元参数化进行模型的仿真,得到对应的疲劳寿命值;其中,对齿轮而言,其结构参数值至少包括齿宽、重合度、压力角、变位系数;
比较仿真得到的疲劳寿命值与目标寿命值,若两者相等,则预测得到的参数就是优化后的参数;若不相等,进行训练集样本补充,继续训练模型,再将目标寿命值作为输入进行参数预测,依据预测参数更新虚拟空间模型结构,进行新一轮超写实仿真,比较目标寿命值与仿真寿命值。循环上述步骤,直至两个寿命值相等。
具体以采煤机摇臂为例进行说明。
首先,依据物理空间采煤机摇臂的受力大小、力的类型、结构、工作时的姿态,轮廓尺寸、相对位置等信息,用UG三维软件构建采煤机摇臂的三维模型,导入ANSYS软件中进行动力学分析,研究动力学与运动状态的关联性,确定反映运动状态的关键参数,如应变、振型、节点位移、节点位置等信息。依据得到的参数布置需要的传感器的类型(加速度、温度传感器、力传感器等),以及传感器在采煤机摇臂位置的布置(传感器布置在应变大的位置);依据物理空间采煤机的材料类型、几何尺寸(长、宽、高、孔尺寸,齿轮几何参数等),利用UG参数化建模在虚拟空间构建采煤机摇臂的三维模型,实现采煤机摇臂的高保真建模。
通过VR技术获得的采煤机摇臂物理空间多物理参数(采煤机摇臂受力、速度、振动、温度、环境温度、煤尘、位置等),将UG中构建的采煤机摇臂三维模型导入有限元分析软件中,进行仿真分析。
依采煤机摇臂齿轮为例,依据有限元参数化技术,通过传感器实时监测数据 (力,速度,温度、位置等)实时动态更新虚拟空间相关参数设置,实现动态疲劳仿真,得到对应的疲劳寿命值。通过调整齿轮对应的结构参数(齿宽、重合度、压力角、变位系数等),继续进行仿真,得到该参数下齿轮对应的疲劳寿命值,构建参数与齿轮疲劳寿命的对应关系。通过调整不同的参数,得到不同的疲劳寿命值,产生齿轮结构参数与对应的疲劳寿命数据样本集。为后期结构参数预测提供数据基础。
依据上述得到的齿轮结构参数与对应的疲劳寿命数据样本集进行模型训练,保存训练好的模型。将齿轮应该具有的寿命值(目标寿命值)作为模型的输入,得到预测的齿轮参数值。以预测的齿轮结构参数值(齿宽、重合度、压力角、变位系数等)作为指导,更新虚拟空间模型的对应的尺寸,通过有限元参数化进行模型的仿真,得到对应的疲劳寿命值。比较仿真得到的疲劳寿命与目标寿命值,若两者相等,则预测得到的参数就是优化后的参数;若不相等,进行训练集样本补充,继续返回去训练模型,再将目标寿命值作为输入进行参数预测,依据预测参数更新虚拟空间模型结构,进行新一轮超写实仿真,比较目标寿命值与仿真寿命值。循环上述步骤,直至两个寿命值相等。最后,输出对应的预测的齿轮结构参数,此时的结构参数为动态优化后的结构参数,实现齿轮结构参数的动态优化。
现有机械设备关键零部件的结构参数一般是采用经验值进行选取,并通过后期的强度校核验证设计的正确性。但在实际应用中,由于在设计阶段无法考虑到实际的工作环境。因此,无论是通过经验值设计经过校核的还是依据理论公式推导计算得到的结构参数值,都会和实际运行环境中要达到的寿命值对应的结构参数存在较大的偏差。
实例验证中选取减速器齿轮的齿宽、变位系数、压力角和重合度作为变量,验证提出方法的有效性。
依据目标寿命得到模型对齿宽参数的动态优化效果图如图2所示,齿宽优化指标值如表1所示。
RMSE R<sup>2</sup> MAPE Score
指标值 0.04317 0.98192 0.02158 0.23316
表1齿宽优化指标值表
由图2和表1可知,模型对齿宽的优化指标值接近理想值,若当优化值与实际值完全相等时,RMSE,MAPE应该等于0,R2值应该等于1,Score应该接近 0。观察实际值发现,各个指标值均接近理想值。因此,模型对齿宽的动态优化是有效的。
依据目标寿命进行变位系数的动态优化,效果如图3所示,变位系数动态优化指标值如表2所示。
RMSE R<sup>2</sup> MAPE Score
指标值 0.01621 0.99171 0.01047 0.08366
表2变位系数动态优化指标值
由图3和表2分析,相比于模型对齿宽优化指标值,模型对变位系数的优化指标值更接近于理想值,也就是说模型对变位系数的优化更接近实际值,而相对百分误差为MAPE=1.047%。
依据目标寿命进行压力角值动态优化,效果如4所示,压力角优化指标值如表3所示。
RMSE R<sup>2</sup> MAPE Score
指标值 0.07517 0.97689 0.09155 0.67519
表3压力角优化指标值
由图4和表3分析可知,模型对压力角的动态优化误差比对齿宽和变位系数的优化误差大,相对百分误差为MAPE=9.155%,但是在误差允许的范围内,该优化结果是可以采取的。由于齿轮压力角对齿轮寿命的影响比齿宽和变位系数的影响小,所以导致对压力角的优化结构误差比较大,但是,该优化结果仍然可以指导实际生产。
依据目标寿命值进行重合度值动态优化,效果如图5所示,重合度优化指标值如表4所示。
RMSE R<sup>2</sup> MAPE Score
指标值 0.05933 0.98745 0.04788 0.14859
表4重合度优化指标值
由图5和表4分析,模型对重合度的动态优化误差比压力角优化误差小,重合度的动态优化相对百分误差值MAPE=4.788%,故该动态优化结果是理想的。而且,重合度对寿命的影响比压力角对寿命的影响大。
综合分析以上动态优化结果,该模型对减速器齿轮结构参数的动态优化结果是接近实际值的,各个参数的动态优化结果相对百分误差为,齿宽动态优化误差 MAPE=2.158%,变位系数动态优化误差MAPE=1.047%,压力角动态优化误差 MAPE=9.155%,重合度动态优化误差MAPE=4.788%.故在误差允许的范围内,该动态优化结果是可取的,则该方法对于机械设备关键零部件结构参数的动态优化是有效的。
区别于现有技术,本发明的基于数字孪生的机械设备零部件结构参数动态优化方法,通过构建高保真模型,在虚拟空间实现物理空间对应设备的数字化镜像,方便后期结构参数修改和超写实仿真;通过进行超写实仿真,在虚拟空间实时动态的同步反映物理空间对应实体的状态,在虚拟空间实现物理设备真实情况的写实运动;同时利用深度学习理论,构建神经网络结构,借助其强大的数字挖掘和映射能力,挖掘建立结构参数和疲劳寿命之间的关系,通过结合高保真模型和超写实仿真环境,实现结构参数的动态优化。通过本发明,实现虚拟空间对物理空间的结构参数的动态优化与反向指导,提高了优化效率和真实性。
本发明虽然已以较佳实施例公开如上,但其并不是用来限定本发明,任何本领域技术人员在不脱离本发明的精神和范围内,都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改,因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰,均属于本发明技术方案的保护范围。

Claims (8)

1.一种基于数字孪生的机械设备零部件结构参数动态优化方法,其特征在于,包括:
构建机械设备全部或部分结构的三维模型,研究三维模型动力学与运动状态的关联性,确定反映运动状态的参数,并根据参数分析确定机械设备全部或部分结构中所需的传感器类型及安装位置,按照分析结果在物理空间物理实体中布置传感器,通过传感器监测数据驱动虚拟空间参数更新,完成高保真建模;
依据物理实体的历史监测数据及物理空间传感器实时更新的监测数据,进行超写实仿真,获取高保真三维模型零件结构参数与对应疲劳寿命数据集作为数据分析样本;
构建深度卷积神经网络模型,用高保真三维模型零件结构参数与对应的疲劳寿命数据样本集训练模型,保存训练好的模型;
将目标疲劳寿命数据作为预测模型的输入,通过传感器实时监测数据动态更新虚拟空间参数和预测模型动态训练,得到动态优化后的零部件结构参数。
2.根据权利要求1所述的基于数字孪生的机械设备零部件结构参数动态优化方法,其特征在于,在进行高保真建模的步骤中,包括步骤:
依据物理空间机械设备零部件结构的受力大小、力的类型、结构、工作时的姿态,轮廓尺寸、相对位置信息,用三维软件构建机械设备零部件结构的三维模型;
导入ANSYS软件中进行动力学分析,研究动力学与运动状态的关联性,确定反映运动状态的关键参数,其中,反映运动状态的关键参数至少包括应变、振型、节点位移、节点位置;
依据得到的关键参数布置需要的传感器的类型,以及传感器在机械设备零部件结构位置的布置;依据物理空间机械设备零部件结构的材料类型、几何尺寸,利用UG参数化建模在虚拟空间构建机械设备零部件结构的三维模型,实现高保真建模。
3.根据权利要求2所述的基于数字孪生的机械设备零部件结构参数动态优化方法,其特征在于,所述传感器类型至少包括加速度、温度传感器和力传感器。
4.根据权利要求2所述的基于数字孪生的机械设备零部件结构参数动态优化方法,其特征在于,传感器布置在机械设备零部件结构上应变大的位置;其中,应变大小通过分析结果的应变图进行判断。
5.根据权利要求1所述的基于数字孪生的机械设备零部件结构参数动态优化方法,其特征在于,通过VR技术获得的机械设备零部件结构物理空间多物理参数,其中,机械设备零部件结构的物理空间多物理参数至少包括受力、速度、振动、温度、环境温度、位置参数。
6.根据权利要求1所述的基于数字孪生的机械设备零部件结构参数动态优化方法,其特征在于,对机械设备零部件结构进行超写实仿真的步骤包括:
依据有限元参数化技术,通过传感器实时监测数据实时动态更新虚拟空间相关参数设置,实现动态疲劳仿真,得到对应的疲劳寿命值;其中,传感器实时监测机械设备零部件结构的监测数据至少包括力、速度、温度及位置参数;
通过调整机械设备零部件结构对应的结构参数,继续进行仿真,得到该参数下机械设备零部件结构对应的疲劳寿命值,构建参数与机械设备零部件的疲劳寿命的对应关系;
通过调整不同的参数,得到不同的疲劳寿命值,产生机械设备零部件结构参数与对应的疲劳寿命数据样本集。
7.根据权利要求1所述的基于数字孪生的机械设备零部件结构参数动态优化方法,其特征在于,构建深度卷积神经网络模型的步骤包括:
利用python语言实现DCNN模型,用测试集进行DCNN模型的训练,进行参数优化,调整模型的层数、模型的具体结构及卷积核的大小,直至训练集中的预测值和实际值的均方误差达到最小,训练预测结果达到最优确定训练次数、模型层数、模型具体结构以及卷积核大小,保存模型;
其中,均方误差表达式为
N表示参加训练的数据量,ypi表示对第i个输入的预测值,yti表示第i个输入对应的实际值;
基于训练好的模型,利用测试集测试模型预测性能,将模型的预测值和实际值做对比,得到模型预测的准确性;
其中,对模型预测结果进行评判时,采用均方根误差(RMSE)、拟合优度检验(R2)、平均绝对百分误差(MAPE)和Score_function三个指标进行评判;
表达式分别为:
预测分析中RMSE越接近0,代表预测结果越精确;
R2值越接近1,代表预测结果越好;
平均绝对百分误差(MAPE)不仅仅考虑预测值与真实值之间的误差,还考虑该误差占真实值的比例;
Score值越接近0,代表预测结果越精确;
上述公式中,yi表示实际值,表示预测值,表示预测的均值,yi表示第i组的寿命对应参数。
8.根据权利要求1所述的基于数字孪生的机械设备零部件结构参数动态优化方法,其特征在于,在将目标疲劳寿命数据作为结构参数预测模型的输入,对高保真三维模型的零件结构参数进行预测的步骤中,包括步骤:
预设机械设备零部件结构的目标寿命值,并将目标寿命值(目标寿命值)作为预测模型的输入,得到预测的机械设备零部件结构参数值;
以预测的机械设备零部件结构参数值作为指导,更新虚拟空间模型的对应的尺寸,通过有限元参数化进行模型的仿真,得到对应的疲劳寿命值;
比较仿真得到的疲劳寿命值与目标寿命值,若两者相等,则预测得到的参数就是优化后的参数;若不相等,进行训练集样本补充,继续训练模型,再将目标寿命值作为输入进行参数预测,依据预测参数更新虚拟空间模型结构,进行新一轮超写实仿真,比较目标寿命值与仿真寿命值;循环上述步骤,直至两个寿命值相等。
CN201910261337.0A 2019-04-02 2019-04-02 基于数字孪生的机械设备零部件结构参数动态优化方法 Active CN110045608B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910261337.0A CN110045608B (zh) 2019-04-02 2019-04-02 基于数字孪生的机械设备零部件结构参数动态优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910261337.0A CN110045608B (zh) 2019-04-02 2019-04-02 基于数字孪生的机械设备零部件结构参数动态优化方法

Publications (2)

Publication Number Publication Date
CN110045608A true CN110045608A (zh) 2019-07-23
CN110045608B CN110045608B (zh) 2022-04-05

Family

ID=67275793

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910261337.0A Active CN110045608B (zh) 2019-04-02 2019-04-02 基于数字孪生的机械设备零部件结构参数动态优化方法

Country Status (1)

Country Link
CN (1) CN110045608B (zh)

Cited By (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110442936A (zh) * 2019-07-24 2019-11-12 中国石油大学(北京) 基于数字孪生模型的设备故障诊断方法、装置及系统
CN110532626A (zh) * 2019-07-31 2019-12-03 西安交通大学 基于数字孪生的航空发动机主轴承剩余寿命预测方法
CN110605709A (zh) * 2019-09-25 2019-12-24 西南交通大学 数字孪生与精密滤波驱动机器人集成系统及其使用方法
CN110795874A (zh) * 2019-09-30 2020-02-14 武汉大学 一种用于柔性电路板制造工艺过程的数字孪生模型
CN110851966A (zh) * 2019-10-30 2020-02-28 同济大学 一种基于深度神经网络的数字孪生模型修正方法
CN110941238A (zh) * 2019-12-18 2020-03-31 广东工业大学 一种基于机器学习的平面运动直线度误差补偿系统及方法
CN111008502A (zh) * 2019-11-25 2020-04-14 北京航空航天大学 一种数字孪生驱动的复杂装备故障预测方法
CN111177942A (zh) * 2020-01-06 2020-05-19 中国矿业大学(北京) 矿井无人化综掘工作面数字孪生智能监控系统
CN111445081A (zh) * 2020-04-01 2020-07-24 浙江大学 产品作业动态调度的数字孪生虚实自适应迭代优化方法
CN111523260A (zh) * 2020-03-18 2020-08-11 上海发电设备成套设计研究院有限责任公司 汽轮机转子数字孪生体构建方法及监控系统
CN111523274A (zh) * 2020-03-18 2020-08-11 上海发电设备成套设计研究院有限责任公司 汽轮机高温承压部件数字孪生体构建方法及监控系统
CN111652406A (zh) * 2020-04-10 2020-09-11 安徽继远软件有限公司 云-边协同下电网设备运行状态预测方法及系统
CN111859566A (zh) * 2020-07-17 2020-10-30 重庆大学 基于数字孪生的表面粗糙度稳定方法
CN112115649A (zh) * 2020-09-29 2020-12-22 郑州轻工业大学 基于数字孪生的立磨机多场耦合系统工艺参数优化方法
CN112163325A (zh) * 2020-09-10 2021-01-01 浙江工业大学 基于数字孪生的rv减速器使用寿命预测方法
CN112355713A (zh) * 2020-09-24 2021-02-12 北京航空航天大学 基于图像预测切削力的方法、装置及设备
CN112417742A (zh) * 2021-01-22 2021-02-26 浙江中自庆安新能源技术有限公司 一种基于数字孪生模型的齿轮箱寿命动态评估方法及系统
WO2021036635A1 (zh) * 2019-08-31 2021-03-04 深圳市广宁股份有限公司 基于振动检测的数字孪生智能健康预测方法及装置
CN112784343A (zh) * 2021-01-29 2021-05-11 中国农业银行股份有限公司 一种基于数字孪生模型的机房设计方法、装置及设备
CN113420465A (zh) * 2021-07-30 2021-09-21 太原理工大学 一种基于数字孪生模型的液压支架全寿命周期管理方法
CN113515848A (zh) * 2021-05-13 2021-10-19 三一重型装备有限公司 矿山设备的结构设计优化方法、系统和可读存储介质
CN113706037A (zh) * 2021-08-31 2021-11-26 天津大学 基于虚拟传感器的耙吸挖泥船的施工安全辅助监控方法
CN113858146A (zh) * 2021-12-02 2021-12-31 徐州安普瑞特能源科技有限公司 一种工作台装置
CN113883869A (zh) * 2021-12-08 2022-01-04 江苏东方恒基通用航空有限公司 一种烘干装置
CN115037590A (zh) * 2022-03-25 2022-09-09 烽火通信科技股份有限公司 一种网络虚拟化体系结构以及虚拟化方法
CN115046575A (zh) * 2021-03-08 2022-09-13 Vega格里沙贝两合公司 通过虚拟孪生对传感器的自动参数化
CN115345034A (zh) * 2022-10-18 2022-11-15 中煤科工开采研究院有限公司 一种液压支架群的数字孪生体的管理方法及系统
CN115730401A (zh) * 2022-11-18 2023-03-03 中复连众风电科技有限公司 一种风电叶片模具参数化设计系统及方法
CN115859508A (zh) * 2022-11-23 2023-03-28 北京百度网讯科技有限公司 流场分析方法、元件模型生成方法、训练方法及装置
CN116205101A (zh) * 2023-02-10 2023-06-02 上海益实智能科技有限公司 一种基于数字孪生的卷烟工厂设备参数化设计方法及装置
CN116611190A (zh) * 2023-07-20 2023-08-18 宁波东力传动设备有限公司 一种轻量化的多级减速器设计方法
CN116822380A (zh) * 2023-08-29 2023-09-29 昆明理工大学 一种基于数字孪生在铜冶炼过程中尾气回收利用的协同优化方法

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2004104917A1 (ja) * 2003-05-20 2004-12-02 Fujitsu Limited センサおよびセンサプログラム
CN101308523A (zh) * 2008-06-24 2008-11-19 中国科学院力学研究所 一种虚拟激光加工过程几何仿真方法
EP2226411A1 (en) * 2009-03-06 2010-09-08 Recherche 2000 Inc. Method for ensuring and monitoring electrolyzer safety and performances
US20130317629A1 (en) * 2012-05-25 2013-11-28 Statistics & Control, Inc. Method of large scale process optimization and optimal planning based on real time dynamic simulation
CN105738109A (zh) * 2016-02-22 2016-07-06 重庆大学 基于稀疏表示与集成学习的轴承故障分类诊断方法
CN106325072A (zh) * 2016-10-12 2017-01-11 浙江理工大学 一种直线伺服系统机械残余振动控制方法
CN106980822A (zh) * 2017-03-14 2017-07-25 北京航空航天大学 一种基于选择性集成学习的旋转机械故障诊断方法
US20170286572A1 (en) * 2016-03-31 2017-10-05 General Electric Company Digital twin of twinned physical system
US20180054376A1 (en) * 2016-08-22 2018-02-22 General Electric Company Internet of things associate
CN108256778A (zh) * 2018-01-31 2018-07-06 北京仿真中心 一种基于机器学习和平行仿真的高级计划与排程方法
CA3051241A1 (en) * 2017-01-24 2018-08-02 Siemens Aktiengesellschaft System and method for cognitive engineering technology for automation and control of systems
CN108806364A (zh) * 2018-06-08 2018-11-13 潍坊护理职业学院 一种基于双网互联的智慧仿真型护理教学系统及构建方法
CN109270899A (zh) * 2018-09-03 2019-01-25 江苏科技大学 一种基于数字孪生的船用柴油机关重件制造过程管控方法
CN109343496A (zh) * 2018-11-14 2019-02-15 中国电子工程设计院有限公司 应用于工业生产的数字孪生系统及其形成方法
WO2019055538A1 (en) * 2017-09-12 2019-03-21 Arconic Inc. SYSTEMS AND METHODS OF ADDITIVE MANUFACTURING

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2004104917A1 (ja) * 2003-05-20 2004-12-02 Fujitsu Limited センサおよびセンサプログラム
CN101308523A (zh) * 2008-06-24 2008-11-19 中国科学院力学研究所 一种虚拟激光加工过程几何仿真方法
EP2226411A1 (en) * 2009-03-06 2010-09-08 Recherche 2000 Inc. Method for ensuring and monitoring electrolyzer safety and performances
US20130317629A1 (en) * 2012-05-25 2013-11-28 Statistics & Control, Inc. Method of large scale process optimization and optimal planning based on real time dynamic simulation
CN105738109A (zh) * 2016-02-22 2016-07-06 重庆大学 基于稀疏表示与集成学习的轴承故障分类诊断方法
US20170286572A1 (en) * 2016-03-31 2017-10-05 General Electric Company Digital twin of twinned physical system
US20180054376A1 (en) * 2016-08-22 2018-02-22 General Electric Company Internet of things associate
WO2018038902A1 (en) * 2016-08-22 2018-03-01 General Electric Company Internet of things associate
CN106325072A (zh) * 2016-10-12 2017-01-11 浙江理工大学 一种直线伺服系统机械残余振动控制方法
CA3051241A1 (en) * 2017-01-24 2018-08-02 Siemens Aktiengesellschaft System and method for cognitive engineering technology for automation and control of systems
CN106980822A (zh) * 2017-03-14 2017-07-25 北京航空航天大学 一种基于选择性集成学习的旋转机械故障诊断方法
WO2019055538A1 (en) * 2017-09-12 2019-03-21 Arconic Inc. SYSTEMS AND METHODS OF ADDITIVE MANUFACTURING
CN108256778A (zh) * 2018-01-31 2018-07-06 北京仿真中心 一种基于机器学习和平行仿真的高级计划与排程方法
CN108806364A (zh) * 2018-06-08 2018-11-13 潍坊护理职业学院 一种基于双网互联的智慧仿真型护理教学系统及构建方法
CN109270899A (zh) * 2018-09-03 2019-01-25 江苏科技大学 一种基于数字孪生的船用柴油机关重件制造过程管控方法
CN109343496A (zh) * 2018-11-14 2019-02-15 中国电子工程设计院有限公司 应用于工业生产的数字孪生系统及其形成方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
CHEN XIUYU;GUAN TIANYI: "Research on the predicting model of convenience store model based on digital twins", 《2018 INTERNATIONAL CONFERENCE ON SMART GRID AND ELECTRICAL AUTOMATION》 *
Q. QI AND F. TAO: "Digital twin and big data towards smart manufacturing and industry 4.0: 360 degree comparison", 《IEEE ACCESS》 *
SCHLUSE M ETAL.: "Experimentable Digital Twins-Streamlining Simulation-based Systems Engineering for Industry", 《IEEE TRANSACTION ON INDUSTRY INFORMATION》 *
于勇等: "数字孪生在工艺设计中的应用探讨", 《航空制造技术》 *
张新生: "基于数字孪生的车间管控系统的设计与实现", 《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》 *
李柏松等: "数字孪生体及其在智慧管网应用的可行性", 《油气储运》 *
谢嘉成: "VR环境下综采工作面"三机"监测与动态规划方法研究", 《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅰ辑》 *

Cited By (49)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110442936B (zh) * 2019-07-24 2021-02-23 中国石油大学(北京) 基于数字孪生模型的设备故障诊断方法、装置及系统
CN110442936A (zh) * 2019-07-24 2019-11-12 中国石油大学(北京) 基于数字孪生模型的设备故障诊断方法、装置及系统
CN110532626B (zh) * 2019-07-31 2021-01-19 西安交通大学 基于数字孪生的航空发动机主轴承剩余寿命预测方法
CN110532626A (zh) * 2019-07-31 2019-12-03 西安交通大学 基于数字孪生的航空发动机主轴承剩余寿命预测方法
WO2021036635A1 (zh) * 2019-08-31 2021-03-04 深圳市广宁股份有限公司 基于振动检测的数字孪生智能健康预测方法及装置
CN110605709B (zh) * 2019-09-25 2022-11-08 西南交通大学 数字孪生与精密滤波驱动机器人集成系统及其使用方法
CN110605709A (zh) * 2019-09-25 2019-12-24 西南交通大学 数字孪生与精密滤波驱动机器人集成系统及其使用方法
CN110795874A (zh) * 2019-09-30 2020-02-14 武汉大学 一种用于柔性电路板制造工艺过程的数字孪生模型
CN110851966A (zh) * 2019-10-30 2020-02-28 同济大学 一种基于深度神经网络的数字孪生模型修正方法
CN111008502A (zh) * 2019-11-25 2020-04-14 北京航空航天大学 一种数字孪生驱动的复杂装备故障预测方法
CN111008502B (zh) * 2019-11-25 2021-07-13 北京航空航天大学 一种数字孪生驱动的复杂装备故障预测方法
CN110941238A (zh) * 2019-12-18 2020-03-31 广东工业大学 一种基于机器学习的平面运动直线度误差补偿系统及方法
CN110941238B (zh) * 2019-12-18 2020-09-29 广东工业大学 一种基于机器学习的平面运动直线度误差补偿系统及方法
CN111177942B (zh) * 2020-01-06 2023-04-18 中国矿业大学(北京) 矿井无人化综掘工作面数字孪生智能监控系统
CN111177942A (zh) * 2020-01-06 2020-05-19 中国矿业大学(北京) 矿井无人化综掘工作面数字孪生智能监控系统
CN111523260A (zh) * 2020-03-18 2020-08-11 上海发电设备成套设计研究院有限责任公司 汽轮机转子数字孪生体构建方法及监控系统
CN111523260B (zh) * 2020-03-18 2022-12-02 上海发电设备成套设计研究院有限责任公司 汽轮机转子数字孪生体构建方法及监控系统
CN111523274A (zh) * 2020-03-18 2020-08-11 上海发电设备成套设计研究院有限责任公司 汽轮机高温承压部件数字孪生体构建方法及监控系统
CN111445081A (zh) * 2020-04-01 2020-07-24 浙江大学 产品作业动态调度的数字孪生虚实自适应迭代优化方法
CN111652406A (zh) * 2020-04-10 2020-09-11 安徽继远软件有限公司 云-边协同下电网设备运行状态预测方法及系统
CN111859566A (zh) * 2020-07-17 2020-10-30 重庆大学 基于数字孪生的表面粗糙度稳定方法
CN111859566B (zh) * 2020-07-17 2023-11-17 重庆大学 基于数字孪生的表面粗糙度稳定方法
CN112163325A (zh) * 2020-09-10 2021-01-01 浙江工业大学 基于数字孪生的rv减速器使用寿命预测方法
CN112163325B (zh) * 2020-09-10 2024-03-29 浙江工业大学 基于数字孪生的rv减速器使用寿命预测方法
CN112355713A (zh) * 2020-09-24 2021-02-12 北京航空航天大学 基于图像预测切削力的方法、装置及设备
CN112115649A (zh) * 2020-09-29 2020-12-22 郑州轻工业大学 基于数字孪生的立磨机多场耦合系统工艺参数优化方法
CN112115649B (zh) * 2020-09-29 2021-09-24 郑州轻工业大学 基于数字孪生的立磨机多场耦合系统工艺参数优化方法
CN112417742B (zh) * 2021-01-22 2021-04-23 浙江中自庆安新能源技术有限公司 一种基于数字孪生模型的齿轮箱寿命动态评估方法及系统
CN112417742A (zh) * 2021-01-22 2021-02-26 浙江中自庆安新能源技术有限公司 一种基于数字孪生模型的齿轮箱寿命动态评估方法及系统
CN112784343A (zh) * 2021-01-29 2021-05-11 中国农业银行股份有限公司 一种基于数字孪生模型的机房设计方法、装置及设备
CN115046575A (zh) * 2021-03-08 2022-09-13 Vega格里沙贝两合公司 通过虚拟孪生对传感器的自动参数化
CN113515848A (zh) * 2021-05-13 2021-10-19 三一重型装备有限公司 矿山设备的结构设计优化方法、系统和可读存储介质
CN113420465A (zh) * 2021-07-30 2021-09-21 太原理工大学 一种基于数字孪生模型的液压支架全寿命周期管理方法
CN113706037A (zh) * 2021-08-31 2021-11-26 天津大学 基于虚拟传感器的耙吸挖泥船的施工安全辅助监控方法
CN113858146A (zh) * 2021-12-02 2021-12-31 徐州安普瑞特能源科技有限公司 一种工作台装置
CN113883869B (zh) * 2021-12-08 2022-08-05 江苏东方恒基通用航空有限公司 一种烘干装置
CN113883869A (zh) * 2021-12-08 2022-01-04 江苏东方恒基通用航空有限公司 一种烘干装置
CN115037590A (zh) * 2022-03-25 2022-09-09 烽火通信科技股份有限公司 一种网络虚拟化体系结构以及虚拟化方法
CN115345034A (zh) * 2022-10-18 2022-11-15 中煤科工开采研究院有限公司 一种液压支架群的数字孪生体的管理方法及系统
CN115345034B (zh) * 2022-10-18 2023-02-03 中煤科工开采研究院有限公司 一种液压支架群的数字孪生体的管理方法及系统
CN115730401A (zh) * 2022-11-18 2023-03-03 中复连众风电科技有限公司 一种风电叶片模具参数化设计系统及方法
CN115730401B (zh) * 2022-11-18 2024-05-07 中复连众风电科技有限公司 一种风电叶片模具参数化设计系统及方法
CN115859508B (zh) * 2022-11-23 2024-01-02 北京百度网讯科技有限公司 流场分析方法、元件模型生成方法、训练方法及装置
CN115859508A (zh) * 2022-11-23 2023-03-28 北京百度网讯科技有限公司 流场分析方法、元件模型生成方法、训练方法及装置
CN116205101A (zh) * 2023-02-10 2023-06-02 上海益实智能科技有限公司 一种基于数字孪生的卷烟工厂设备参数化设计方法及装置
CN116611190A (zh) * 2023-07-20 2023-08-18 宁波东力传动设备有限公司 一种轻量化的多级减速器设计方法
CN116611190B (zh) * 2023-07-20 2023-10-03 宁波东力传动设备有限公司 一种轻量化的多级减速器设计方法
CN116822380A (zh) * 2023-08-29 2023-09-29 昆明理工大学 一种基于数字孪生在铜冶炼过程中尾气回收利用的协同优化方法
CN116822380B (zh) * 2023-08-29 2023-11-14 昆明理工大学 一种基于数字孪生在铜冶炼过程中尾气回收利用的协同优化方法

Also Published As

Publication number Publication date
CN110045608B (zh) 2022-04-05

Similar Documents

Publication Publication Date Title
CN110045608A (zh) 基于数字孪生的机械设备零部件结构参数动态优化方法
CN104615840B (zh) 一种数字仿真模型的修正方法及系统
CN104866525A (zh) 一种虚拟核电站仿真教学系统与构建方法
CN107238388A (zh) 多无人机任务分配与航迹规划联合优化方法及装置
CN107590853A (zh) 一种城市建筑群震害高真实度展示方法
CN103331756A (zh) 一种机械臂运动控制方法
CN105069826A (zh) 弹性物体变形运动的建模方法
CN109918726A (zh) 一种力学结构异常状态快速识别方法、存储介质
CN110188417A (zh) 一种基于多级超单元的复杂结构有限元模型修正方法
CN104570736A (zh) 一种星-臂耦合系统的动力学参数在轨辨识方法和装置
KR20230065343A (ko) 메쉬 표현 및 그래프 신경망을 사용한 물리적 환경 시뮬레이션
CN103914872A (zh) 一种基于简化模态分析法的树动画模拟方法
CN113654756A (zh) 一种海上浮式风机的主动式实时混合模型试验方法
CN115659756A (zh) 一种输电杆塔防风性能分析方法
CN111680823A (zh) 一种风向信息预测方法及系统
CN104834790B (zh) 一种空间尺寸链公差分析方法
Zhang et al. Genetic algorithms for optimal design of vehicle suspensions
CN116629775B (zh) 一种基于物联网技术的装配式建筑建材智能生产控制系统
Du et al. Industrial robot digital twin system motion simulation and collision detection
CN114818128B (zh) 船体局部曲面优化神经网络建模方法及优化方法
CN115293069B (zh) 一种用于飞行器外流场仿真控制参数智能优化的系统
CN103529803A (zh) 一种基于云端虚拟化的烧结偏析布料先进控制实验系统
CN115455754A (zh) 一种基于数字孪生的矿山液压支架设计方法
CN104820728A (zh) 一种水田激光平地机调平系统动力学建模方法
CN107967387A (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
GR01 Patent grant
GR01 Patent grant