CN104992038B - 一种刚柔-机电耦合进给***的动态性能优化设计方法 - Google Patents

一种刚柔-机电耦合进给***的动态性能优化设计方法 Download PDF

Info

Publication number
CN104992038B
CN104992038B CN201510452283.8A CN201510452283A CN104992038B CN 104992038 B CN104992038 B CN 104992038B CN 201510452283 A CN201510452283 A CN 201510452283A CN 104992038 B CN104992038 B CN 104992038B
Authority
CN
China
Prior art keywords
feed system
model
design
evaluation index
mechanical
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.)
Expired - Fee Related
Application number
CN201510452283.8A
Other languages
English (en)
Other versions
CN104992038A (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.)
Tianjin University
Original Assignee
Tianjin 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 Tianjin University filed Critical Tianjin University
Priority to CN201510452283.8A priority Critical patent/CN104992038B/zh
Publication of CN104992038A publication Critical patent/CN104992038A/zh
Application granted granted Critical
Publication of CN104992038B publication Critical patent/CN104992038B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Numerical Control (AREA)
  • Feedback Control In General (AREA)

Abstract

一种刚柔‑机电耦合进给***的动态性能优化设计方法:建立刚柔‑机电耦合进给***模型;确定进给***动态性能优化的设计变量、设计范围以及评价指标;采用中心复合试验设计方法,在变量设计空间内选取试验样本点;利用Matlab/Simulink软件对试验样本点进行分析计算,获取相应的响应值;建立进给***的机械结构参数和控制参数与动态性能评价指标的响应面模型;分析设计变量对进给***动态性能的灵敏度影响,确定所述影响规律,并进行进给***的机械结构参数和控制参数的优化改进。本发明通过对进给***动态性能响应面模型分析计算,得到机械结构和控制参数对进给***动态性能的影响规律,指导进给***结构及控制参数的优化改进。

Description

一种刚柔-机电耦合进给***的动态性能优化设计方法
技术领域
本发明涉及一种进给***的动态性能优化设计方法。特别是涉及一种刚柔-机电耦合进给***的动态性能优化设计方法。
背景技术
目前的机床进给***机电耦合模型中机械结构的建模方法大多为基于集中质量的刚体参数化建模方法,虽然建模效率较高,但无法真实反映机床各关键结构件的柔性变形。同时,进给***的运动精度与机电耦合特性密切相关,其机械结构、伺服控制性能及其两者间的相互影响,共同制约着进给***的动态性能。
因此,寻求一种刚柔-机电耦合进给***建模方法,通过对所建立响应面模型仿真分析,揭示机械结构参数、伺服控制参数对进给***的动态性能的影响规律,并指导设计人员进行优化设计,对改善进给***性能加快设计进程具有重要意义。
发明内容
本发明所要解决的技术问题是,提供一种能够更为全面的反映机床进给***的动态性能的刚柔-机电耦合进给***的动态性能优化设计方法。
本发明所采用的技术方案是:一种刚柔-机电耦合进给***的动态性能优化设计方法,包括如下步骤:
1)建立刚柔-机电耦合进给***模型;
2)确定进给***动态性能优化的设计变量、设计范围以及评价指标;
3)采用中心复合试验设计方法,在变量设计空间内选取试验样本点;
4)利用Matlab/Simulink软件对试验样本点进行分析计算,获取相应的响应值;
5)建立进给***的机械结构参数和控制参数与动态性能评价指标的响应面模型;
6)分析设计变量对进给***动态性能的灵敏度影响,确定所述影响规律,并进行进给***的机械结构参数和控制参数的优化改进。
步骤1)所述的建立刚柔-机电耦合进给***模型,包括:
(1)建立机床进给***的CAD模型,并导入ANSYS软件中建立进给***有限元模型;
(2)在Matlab/Simulink软件中建立进给***机械结构的刚柔耦合模型;
(3)将进给***刚柔耦合模型与控制进给***的伺服控制***进行耦合得到刚柔-机电耦合进给***模型。
第(2)步所建立的进给***机械结构的刚柔耦合模型中,床身和工作台模型是基于模态缩减法建立的柔性动力学模型,而滚珠丝杠模型是基于集中参数法建立的刚体动力学模型,所述的刚柔耦合模型是通过床身、工作台和滚珠丝杠之间的作用关系,将所建立的柔性动力学模型和刚体动力学模型耦合得到进给***刚柔耦合模型。
步骤2)中:
所述的设计变量包括进给***的机械结构参数和控制参数,其中,所述的机械结构参数有丝杠导程、螺母副刚度和工作台质量,所述的控制参数有位置环增益和速度环增益;
所述的变量设计范围,是将所述的设计变量的初始值作为中间值,将初始值增大20%作为设计变量的最大值,将初始值减小20%作为设计变量的最小值而构成的范围;
所述的评价指标是指进给***的响应效率评价指标、定位精度评价指标和延迟特性评价指标。
所述的进给***的响应效率评价指标为:
其中,tp为进给***走完行程的理论时刻,tc为实际位移到达理想位移误差ε范围内的时刻,te为测试结束时刻;
所述的定位精度评价指标,是选取进给***走完行程的理论时刻作为特征点,将在理论时刻的实际位移与理想位移的比作为定位精度的评价指标;
所述的延迟特性评价指标,是将进给***实际达到峰值的时刻与进给***到达峰值的理论时刻的比作为延迟特性的评价指标。
步骤5)中所述的响应面模型,是以进给***的机械结构参数和控制参数作为设计变量,进给***的响应效率评价指标、定位精度评价指标和延迟特性评价指标作为输出响应值,建立的n个设计变量的二阶响应面模型表示为:
式中:y为输出变量;xi为设计变量;n为设计变量的个数;β为待定系数,由最小二乘回归法拟合得到。
步骤6)中所述的分析设计变量对进给***动态性能的灵敏度影响,确定所述影响规律,并进行进给***的机械结构参数和控制参数的优化改进,是通过对进给***的响应效率评价指标、定位精度评价指标和延迟特性评价指标的响应面模型的分析计算,所述的分析计算是采用局部灵敏度分析方法,具体是对机械结构参数和控制参数中的一个设计变量在变量设计范围之内进行参数调整,其它的设计变量取初始值,计算模型在所述的进行参数调整的设计变量每次参数调整前后响应值的相对变化量,就是进给***结构参数和控制参数对评价指标的灵敏度。按所述的方法依次分析其它设计变量,根据得到的灵敏度分析结果,确定设计变量对进给***动态性能影响规律。针对影响规律,对原始设计方案的结构参数和控制参数进行相应的调整改进,提高进给***的动态性能。
本发明的一种刚柔-机电耦合进给***的动态性能优化设计方法,是基于模态缩减法的刚柔耦合建模方法,并与伺服控制***耦合建立进给***机电耦合模型仿真分析,能够更为全面的反映机床进给***的动态性能。本发明通过对进给***动态性能响应面模型分析计算,得到机械结构和控制参数对进给***动态性能的影响规律,指导进给***结构及控制参数的优化改进。
附图说明
图1是本发明刚柔-机电耦合进给***的动态性能优化设计方法的流程图。
具体实施方式
下面结合实施例和附图对本发明的一种刚柔-机电耦合进给***的动态性能优化设计方法做出详细说明。
如图1所示,本发明的一种刚柔-机电耦合进给***的动态性能优化设计方法,其特征在于,包括如下步骤:
1)建立刚柔-机电耦合进给***模型,包括:
(1)建立机床进给***的CAD模型,并导入ANSYS软件中建立进给***有限元模型;
在Pro/E软件中建立机床进给***的实体模型,并进行简化处理,删除对分析结果近乎无影响的小结构信息。
将机床进给***简化模型导入到有限元分析软件ANSYS中,定义材料属性,建立实体单元和Combin14弹簧阻尼单元,在各结构件的结合面处建立相应的硬点对,添加边界条件和载荷、划分有限元网格,其中,导轨、丝杠作为机床关键运动部件,其连接刚度通过导轨丝杠产品手册得到,阻尼参数参考相关文献获得,根据查得数据设置弹簧阻尼单元的实常数,建立进给***有限元模型。
(2)在Matlab/Simulink软件中建立进给***机械结构的刚柔耦合模型;
应用ANSYS命令流来提取所述进给***有限元模型各结构件的full文件、各结合面所有硬点位置坐标参数以及刚度阻尼值等参数,并应用模态缩减软件MORforANSYS,对各结构件分别进行模态缩减处理,获取各结构件的质量矩阵、刚度矩阵、输入和输出矩阵,产生的缩减模型几乎包含了原始有限元模型的所有参数信息。在本发明实施例中将各结构件缩减后的矩阵设为100阶,以保证缩减的模型的精度。
将各结构件的缩减矩阵导入到Matlab软件中,建立多自由度***的运动微分方程,可表示为:
式中,MRi、DRi、KRi为缩减后进给***的结构件i的质量矩阵、阻尼矩阵和刚度矩阵;B、C分别为进给***的输入和输出矩阵。
在Matlab/Simulink环境下,采用状态空间建模方法,把各结构件的运动微分方程转化成状态空间方程,如下式所示:
应用集中参数法建立滚珠丝杠结构的刚体动力学模型,同样用上述方法转化成状态空间方程。
所建立的进给***机械结构的刚柔耦合模型中,床身和工作台模型是基于模态缩减法建立的柔性动力学模型,而滚珠丝杠模型是基于集中参数法建立的刚体动力学模型,所述的刚柔耦合模型是通过床身、工作台和滚珠丝杠之间的作用关系,将所建立的柔性动力学模型和刚体动力学模型耦合得到进给***刚柔耦合模型。
(3)将进给***刚柔耦合模型与控制进给***的伺服控制***进行耦合得到刚柔-机电耦合进给***模型。
本发明实施例中采用典型的三环PID伺服控制***,在Matlab/Simulink中建立了包含位置环、速度环、电流环与交流同步伺服电机的伺服控制***模型,根据伺服控制***与进给***机械结构间的输入输出关系,将伺服控制***模型与进给***机械结构的刚柔耦合模型进行耦合,建立刚柔-机电耦合进给***仿真模型。
2)确定进给***动态性能优化的设计变量、设计范围以及评价指标,其中:
所述的设计变量包括进给***的机械结构参数和控制参数,其中,所述的机械结构参数有丝杠导程、螺母副刚度和工作台质量,所述的控制参数有位置环增益和速度环增益;
所述的变量设计范围,是将所述的设计变量的初始值作为中间值,将初始值增大20%作为设计变量的最大值,将初始值减小20%作为设计变量的最小值而构成的范围;
所述的评价指标是指进给***的响应效率评价指标、定位精度评价指标和延迟特性评价指标。其中:
(1)所述的进给***的响应效率评价指标为:
其中,tp为进给***走完行程的理论时刻,tc为实际位移到达理想位移误差ε范围内的时刻,te为测试结束时刻;
(2)所述的定位精度评价指标,是选取进给***走完行程的理论时刻作为特征点,将在理论时刻的实际位移与理想位移的比作为定位精度的评价指标;
(3)所述的延迟特性评价指标,是将进给***实际达到峰值的时刻与进给***到达峰值的理论时刻的比作为延迟特性的评价指标。
3)采用中心复合试验设计方法,在变量设计范围内选取试验样本点;
所述试验样本点的选择对响应面的建立至关重要,不仅会影响所述响应面模型的精度,甚至导致响应面模型难以建立。本发明的试验设计选择中心复合试验设计方法,是根据实验设计理论来确定合理的试验样本点以保证试验的准确性和高效性。
所述中心复合试验设计方法只在中心点和扩展点处做数值分析试验,能够从所设定的试验区域中选取最具代表性的试验样本点。本发明实施例中进给***动态性能优化问题,共有5个设计变量,对于响应面模型,所需的最少试验样本点的个数为10。本发明通过中心复合试验选取70个试验样本点,以保证所建立的响应面模型的精确度。
4)利用Matlab/Simulink软件对试验样本点进行分析计算,获取相应的响应值;
在Matlab/Simulink环境下,输入测试信号为位移斜坡信号,即进给***Z轴在静止状态下,通过位移斜坡信号激励,从静止状态匀速运动到指定位置,以获得进给***对位移斜坡信号的位移响应。根据动态性能评价指标对比输入信号与输出响应,可以得到反映进给***性能的响应值。根据试验样本点的设计变量及设计范围,改变所述模型的分析参数,分别求出所有样本点对应的响应值。
5)建立进给***的机械结构参数和控制参数与动态性能评价指标的响应面模型;
响应面模型是用显性的函数表达式来表示优化问题中隐性的所述设计输入变量与所述输出响应值之间的关系。由于二阶响应面模型能够拟合多种函数形式,且参数容易估计,本发明选择具有较高准确性的二阶响应面模型。
本发明所述的响应面模型,是以进给***的机械结构参数和控制参数作为设计变量,进给***的响应效率评价指标、定位精度评价指标和延迟特性评价指标作为输出响应值,建立的n个设计变量的二阶响应面模型可表示为:
式中:y为输出变量;xi为设计变量;n为设计变量的个数;β为待定系数,由最小二乘回归法拟合得到。
本发明实施例中,用编码变量值代替实际变量值进行响应面拟合,通过响应面模型的建立可以得到各设计变量对***性能指标的影响度,以及响应面模型的预测能力,为下一步所述响应面模型的优化奠定了基础。
6)分析设计变量对进给***动态性能的灵敏度影响,确定所述影响规律,并进行进给***的机械结构参数和控制参数的优化改进。
所述的分析设计变量对进给***动态性能的灵敏度影响,确定所述影响规律,并进行进给***的机械结构参数和控制参数的优化改进,是通过对进给***的响应效率评价指标、定位精度评价指标和延迟特性评价指标的响应面模型的分析计算,所述的分析计算是采用局部灵敏度分析方法,具体是对机械结构参数和控制参数中的一个设计变量在变量设计范围之内进行参数调整,其它的设计变量取初始值,计算模型在所述的进行参数调整的设计变量每次参数调整前后响应值的相对变化量,就是进给***结构参数和控制参数对评价指标的灵敏度。按所述的方法依次分析其它设计变量,根据得到的灵敏度分析结果,确定设计变量对进给***动态性能影响规律。针对影响规律,对原始设计方案的结构参数和控制参数进行相应的调整改进,提高进给***的动态性能。
以上实施例仅用于说明本发明的技术而非限制,尽管参照实例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的进行修改或者同等替换,而不脱离本发明技术的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (1)

1.一种刚柔-机电耦合进给***的动态性能优化设计方法,其特征在于,包括如下步骤:
1)建立刚柔-机电耦合进给***模型;所述的建立刚柔-机电耦合进给***模型,包括:
(1)建立机床进给***的CAD模型,并导入ANSYS软件中建立进给***有限元模型;
(2)在Matlab/Simulink软件中建立进给***机械结构的刚柔耦合模型;
所建立的进给***机械结构的刚柔耦合模型中,床身和工作台模型是基于模态缩减法建立的柔性动力学模型,而滚珠丝杠模型是基于集中参数法建立的刚体动力学模型,所述的刚柔耦合模型是通过床身、工作台和滚珠丝杠之间的作用关系,将所建立的柔性动力学模型和刚体动力学模型耦合得到进给***刚柔耦合模型;
(3)将进给***刚柔耦合模型与控制进给***的伺服控制***进行耦合得到刚柔-机电耦合进给***模型;
2)确定进给***动态性能优化的设计变量、设计范围以及评价指标;其中:
所述的设计变量包括进给***的机械结构参数和控制参数,其中,所述的机械结构参数有丝杠导程、螺母副刚度和工作台质量,所述的控制参数有位置环增益和速度环增益;
所述的变量设计范围,是将所述的设计变量的初始值作为中间值,将初始值增大20%作为设计变量的最大值,将初始值减小20%作为设计变量的最小值而构成的范围;
所述的评价指标是指进给***的响应效率评价指标、定位精度评价指标和延迟特性评价指标;
所述的进给***的响应效率评价指标为:
其中,tp为进给***走完行程的理论时刻,tc为实际位移到达理想位移误差ε范围内的时刻,te为测试结束时刻;
所述的定位精度评价指标,是选取进给***走完行程的理论时刻作为特征点,将在理论时刻的实际位移与理想位移的比作为定位精度的评价指标;
所述的延迟特性评价指标,是将进给***实际达到峰值的时刻与进给***到达峰值的理论时刻的比作为延迟特性的评价指标;
3)采用中心复合试验设计方法,在变量设计空间内选取试验样本点;
4)利用Matlab/Simulink软件对试验样本点进行分析计算,获取相应的响应值;
5)建立进给***的机械结构参数和控制参数与动态性能评价指标的响应面模型;
所述的响应面模型,是以进给***的机械结构参数和控制参数作为设计变量,进给***的响应效率评价指标、定位精度评价指标和延迟特性评价指标作为输出响应值,建立的n个设计变量的二阶响应面模型表示为:
式中:y为输出变量;xi为设计变量;n为设计变量的个数;β为待定系数,由最小二乘回归法拟合得到;
6)分析设计变量对进给***动态性能的灵敏度影响,确定所述影响规律,并进行进给***的机械结构参数和控制参数的优化改进;
所述的分析设计变量对进给***动态性能的灵敏度影响,确定所述影响规律,并进行进给***的机械结构参数和控制参数的优化改进,是通过对进给***的响应效率评价指标、定位精度评价指标和延迟特性评价指标的响应面模型的分析计算,所述的分析计算是采用局部灵敏度分析方法,具体是对机械结构参数和控制参数中的一个设计变量在变量设计范围之内进行参数调整,其它的设计变量取初始值,计算模型在所述的进行参数调整的设计变量每次参数调整前后响应值的相对变化量,就是进给***结构参数和控制参数对评价指标的灵敏度;按上述的方法依次分析其它设计变量,根据得到的灵敏度分析结果,确定设计变量对进给***动态性能影响规律;针对影响规律,对原始设计方案的结构参数和控制参数进行相应的调整改进,提高进给***的动态性能。
CN201510452283.8A 2015-07-27 2015-07-27 一种刚柔-机电耦合进给***的动态性能优化设计方法 Expired - Fee Related CN104992038B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510452283.8A CN104992038B (zh) 2015-07-27 2015-07-27 一种刚柔-机电耦合进给***的动态性能优化设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510452283.8A CN104992038B (zh) 2015-07-27 2015-07-27 一种刚柔-机电耦合进给***的动态性能优化设计方法

Publications (2)

Publication Number Publication Date
CN104992038A CN104992038A (zh) 2015-10-21
CN104992038B true CN104992038B (zh) 2018-10-16

Family

ID=54303852

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510452283.8A Expired - Fee Related CN104992038B (zh) 2015-07-27 2015-07-27 一种刚柔-机电耦合进给***的动态性能优化设计方法

Country Status (1)

Country Link
CN (1) CN104992038B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105893709B (zh) * 2016-04-29 2018-11-06 天津大学 一种数控机床进给***滚珠丝杠预紧力的试验优化方法
CN106777820A (zh) * 2017-01-22 2017-05-31 西安交通大学 一种机床基础大件动态性能优化设计方法
CN107102620B (zh) * 2017-04-26 2019-06-21 南京理工大学 一种数控机床双驱进给***模态建模方法
CN108920882B (zh) * 2018-08-15 2020-01-10 东南大学 高速机床滚珠丝杠进给***动态设计方法
CN109815525B (zh) * 2018-12-06 2023-04-07 天津大学 一种进给***动态特性多目标优化设计方法
CN109613824A (zh) * 2018-12-13 2019-04-12 广东工业大学 一种滚珠丝杠驱动的刚柔耦合运动平台及控制方法
CN111752213B (zh) * 2019-12-26 2021-11-19 西安交通大学 一种应用于高速高精直驱进给***的机电耦合分析方法
JP2023505617A (ja) * 2020-02-28 2023-02-09 スリーエム イノベイティブ プロパティズ カンパニー 高度モデル予測制御のための深層因果学習
CN112000011B (zh) * 2020-08-12 2022-05-24 深圳市烨嘉为技术有限公司 小龙门数控加工中心的机电耦合分析与优化方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103366043A (zh) * 2013-05-31 2013-10-23 天津大学 一种基于响应面模型的机床床身三点支撑优化设计方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103366043A (zh) * 2013-05-31 2013-10-23 天津大学 一种基于响应面模型的机床床身三点支撑优化设计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于刚柔-机电耦合的机床直线电机进给***动态性能仿真研究;姚延风 等;《振动与冲击》;20110131;第30卷(第1期);第191-196页 *
数控机床机电***刚柔耦合建模技术;熊雪峰 等;《装备制造技术》;20150630(第6期);第1-4页 *

Also Published As

Publication number Publication date
CN104992038A (zh) 2015-10-21

Similar Documents

Publication Publication Date Title
CN104992038B (zh) 一种刚柔-机电耦合进给***的动态性能优化设计方法
Zhang et al. Application framework of digital twin-driven product smart manufacturing system: A case study of aeroengine blade manufacturing
Li et al. Dual sliding mode contouring control with high accuracy contour error estimation for five-axis CNC machine tools
CN109240204B (zh) 一种基于两步法的数控机床热误差建模方法
US10120368B2 (en) Manufacturing adjustment system that adjusts manufacturing statuses of multiple machines
CN103995937A (zh) 基于响应面和遗传算法的精密机床质量匹配设计方法
CN103197596B (zh) 一种数控加工参数自适应模糊控制规则优化方法
Fujishima et al. Thermal displacement reduction and compensation of a turning center
CN102109822A (zh) 伺服电机的变结构双自由度智能集成控制方法
JP7151547B2 (ja) 予測制御開発装置、予測制御開発方法、及び予測制御開発プログラム
Brecher et al. Evaluation of toolpath quality: User-assisted CAM for complex milling processes
CN103927411B (zh) 一种龙门式激光切割机床的设计方法
CN113569453B (zh) 一种面向加工的数字孪生产品模型的迁移方法
CN105938329A (zh) 一种数控机床动态特性仿真及其与数控***集成的方法
WO2020037747A1 (zh) 一种面向动态精度的数控***插补生成运动指令评价方法
Gong et al. Tool orientation optimization method based on ruled surface using genetic algorithm
Yong et al. Dynamic characteristic optimization of ball screw feed drive in machine tool based on modal extraction of state space model
Shi et al. Online monitoring dynamic characteristics in thin-walled structure milling: A physics-constrained bayesian updating approach
Zaeh et al. Prediction of the dynamic behaviour of machine tools during the design process using mechatronic simulation models based on finite element analysis
CN109446721A (zh) 基于标识符软件线程执行顺序排列的机床工艺交互算法
CN103413049A (zh) 基于机电耦合特性的并联机床结构优化参数值的获取方法
Naumann et al. Optimization of characteristic diagram based thermal error compensation via load case dependent model updates
Lou et al. The selection of key temperature measuring points for the compensation of thermal errors of CNC machining tools
CN110989357B (zh) 一种复杂机电***的辨识控制方法和***
Yau et al. NC simulation with dynamic errors due to high-speed motion

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP02 Change in the address of a patent holder
CP02 Change in the address of a patent holder

Address after: 300350 Haijing garden, Haihe Education Park, Jinnan, Tianjin, 135, Tianjin University.

Patentee after: Tianjin University

Address before: 300072 Tianjin City, Nankai District Wei Jin Road No. 92

Patentee before: Tianjin University

CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20181016