CN105302070A - 一种非正交摆头转台类五轴机床后置处理方法 - Google Patents

一种非正交摆头转台类五轴机床后置处理方法 Download PDF

Info

Publication number
CN105302070A
CN105302070A CN201510839453.8A CN201510839453A CN105302070A CN 105302070 A CN105302070 A CN 105302070A CN 201510839453 A CN201510839453 A CN 201510839453A CN 105302070 A CN105302070 A CN 105302070A
Authority
CN
China
Prior art keywords
machine tool
axis
shaft
lathe
axle
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
CN201510839453.8A
Other languages
English (en)
Other versions
CN105302070B (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.)
Shandong Sma Pharmatech Co ltd
Original Assignee
Shandong 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 Shandong University of Technology filed Critical Shandong University of Technology
Priority to CN201510839453.8A priority Critical patent/CN105302070B/zh
Publication of CN105302070A publication Critical patent/CN105302070A/zh
Application granted granted Critical
Publication of CN105302070B publication Critical patent/CN105302070B/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
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
    • G05B19/404Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by control arrangements for compensation, e.g. for backlash, overshoot, tool offset, tool wear, temperature, machine construction errors, load, inertia
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/36Nc in input of data, input key till input tape
    • G05B2219/36573X, y, z and tooloffset values or direction values

Landscapes

  • Engineering & Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Numerical Control (AREA)

Abstract

本发明公开了一种非正交摆头转台类五轴机床后置处理方法。本发明以平动轴Y和回转轴B作为倾斜轴的摆头转台类五轴机床作为研究对象,分析该类型机床的结构特点。然后,分别将倾斜轴Y轴与OXZ平面之间夹角和刀具摆长作为变量,建立前置刀位数据与该类机床各坐标轴之间的运动变换方程,并推导出各坐标轴的计算公式。在此基础上,利用C++语言开发出适用于该类型五轴数控机床的后置处理软件。该方法不仅能够将刀位文件转换为机床可以识别的数控程序,而且可以在数控程序中进行五轴刀具长度补偿、工件原点偏置补偿功能和非正交轴角度可以修改。

Description

一种非正交摆头转台类五轴机床后置处理方法
技术领域
本发明涉及一种数控机床后置处理方法,尤其涉及一种非正交摆头转台类五轴机床后置处理方法,属于五轴数控加工机床领域。
背景技术
五轴数控机床通常由三平动轴和两回转轴组成,通过平动轴和回转轴的不同配置而构成各种不同类型的五轴机床。一般可将其分为三种基本类型:1)双摆头类型,2)双转台类型,3)摆头转台类型,而其他类型的五轴机床都是在其基础上演化而来。因此,许多学者针对各种类型五轴数控机床后置处理方法开展了大量研究工作,如双摆头类五轴机床、双转台类五轴机床、摆头转台类五轴机床、非正交回转轴双转台五轴机床等。现有技术中,主要是针对非正交回转轴的双转台五轴机床开发的专用后置处理程序,其中非正交轴仅指回转轴,而没有考虑平动轴也为非正交轴的情况。因此,目前文献对于非正交摆头转台类五轴机床,尤其平动轴和回转轴同时为非交情况研究较少。
发明内容
本发明的目的是提供一种具有在数控程序中进行五轴刀具长度补偿、工件原点偏置补偿和非正交轴角度修改功能的非正交摆头转台类五轴机床后置处理方法。
本发明的目的是通过以下技术方案实现的:
一种非正交摆头转台类五轴机床后置处理方法,该方法包括以下步骤:
A、分析以平动轴Y和回转轴B作为倾斜轴的摆头转台类五轴机床的结构特征,分别将倾斜轴Y轴与OXZ平面之间夹角和刀具摆长作为变量,建立前置刀位数据与该类机床各坐标轴之间的运动变换方程,并推导出各坐标轴的计算公式;
B、在机床各平动轴和回转轴的计算公式基础上,结合西门子数控系统或其他数控系统提供的宏变量(如R参数)功能,基于VC++6.0平台利用C++语言开发出一种适用于非正交摆头转台类五轴控机床的后置处理软件;其中在该软件界面中,刀具摆长可以作为宏变量或数值输入,倾斜轴Y轴与OXZ平面之间夹角作为数值输入,工件原点在机床坐标系中的偏置量作为数值输入。
由上述发明技术方案可以看出,该方法不仅能够将刀位文件转换为机床可以识别的数控程序,而且具有在数控程序中进行五轴刀具长度补偿、工件原点偏置补偿和非正交轴角度修改功能。
附图说明
图1为非正交摆头转台类五轴机床;
图2为非正交摆头转台类五轴机床运动链;
图3为非正交摆头转台类五轴机床中的坐标系;
图4为后置处理软件界面。
具体实施方式
本发明的非正交摆头转台类五轴机床后置处理方法,以平动轴Y和回转轴B作为倾斜轴的摆头转台类五轴机床为例进行阐述,如图1所示,其较佳的具体实施方式是,包括:
A、分析以平动轴Y和回转轴B作为倾斜轴的摆头转台类五轴机床的结构特征,分别将倾斜轴Y轴与OXZ平面之间夹角和刀具摆长作为变量,建立前置刀位数据与该类机床各坐标轴之间的运动变换方程,并推导出各坐标轴的计算公式。
首先,以平动轴Y和回转轴B作为倾斜轴的摆头转台类五轴数控机床作为研究对象,如图1所示;该机床由三平动轴XYZ和两回转轴AB组成,其中X轴和YZ轴夹角都等于90°,而Y轴与Z轴夹角小于90°,且回转轴B的中心平行于Y轴,与Z轴夹角亦小于90°,回转轴A的中心平行于X轴,机床主轴与回转轴B夹角大于90°。XZ轴方向的平动和绕A轴的转动由机床转台实现,而Y轴方向的平动和绕B轴的转动由机床主轴实现;根据该机床的具体结构,可以得到该机床各坐标轴之间的运动关系,即机床运动链,它由回转工作台、平动工作台、床身、主轴和刀具等单元按顺序串联而成,如图2所示。
其次,为将由CAM得到的刀位数据转换为机床数控加工程序,需要根据该机床的结构特点来确定机床的运动变化关系,然后再将刀位数据变换分解到机床各坐标轴上,从而获得机床各坐标轴的数值。为描述该机床运动,在该五轴机床上建立图3所示各相关坐标系,其中O m X m Y m Z m为机床坐标系,其原点O m位于A轴转台右端面中心,X m轴方向平行于A轴回转中心,Z m轴方向为竖直向上,Y m轴方向则垂直于O m X m Z m平面;O m1 X m1 Y m1 Z m1为与回转轴B固连的坐标系,其原点为主轴回转中心与B轴回转中心的交点,各坐标轴方向与机床坐标系一致;O m2 X m2 Y m2 Z m2为与回转轴B固连的坐标系,是机床Y轴、B轴实际运动的坐标系,其原点为主轴回转中心与B轴回转中心的交点,各坐标轴方向由机床坐标系绕X m轴旋转α角得到;O w X w Y w Z w为与工件固连的工件坐标系,其原点在机床坐标系中的位置矢量=(d x ,d y ,d z ),且各坐标轴方向与机床坐标系一致,而前置刀位数据是在该坐标系下表示;O t X t Y t Z t为与刀具固连的刀具坐标系,其原点设在刀位点上,其各坐标轴方向与机床坐标系一致。在机床初始状态时,假设刀具轴线平行于Z轴,坐标系O m2 X m2 Y m2 Z m2与机床坐标系O m X m Y m Z m的原点重合。假设坐标系O m1 X m1 Y m1 Z m1的原点到刀具坐标系的原点O t的距离为L(即刀具摆长),则刀位点O t在坐标系O m1 X m1 Y m1 Z m1中的位置为(0,0,-L);在刀具坐标系中,刀位点的位置矢量和刀轴矢量分别为(0,0,0)和(0,0,1)。假设机床平动轴相对于初始状态的移动量为(X,Y,Z),回转轴相对于初始状态的转动量分别为BA(其正方向如图3所示),工件坐标系中刀位点和刀轴矢量分别为(x,y,z)和(i,j,k)。
根据该机床运动链和各坐标系,可分别建立回转轴和平动轴的运动变换方程:
(1)
(2)
式中: T R 分别为平移和回转运动的齐次变换矩阵:
是由编程坐标系O w X w Y w Z w到机床坐标系O m X m Y m Z m的变换矩阵:
由式(1)可以得到:
(3)
若该机床B轴运动范围(-0.5π,0.5π),A轴运动范围[-π,π],则。由式(3)可解得
(4)
式中:
由式(2)可以得到各平动轴的计算公式:
(5)
因此,由式(4)和式(5)就可以确定机床各回转轴和平动轴的数值,即对前置刀位数据进行后置处理而得到机床可以识别执行的数控加工程序,其中机床平动轴计算公式中包含刀具摆长变量L、工件原点偏置变量d x d y d z 和倾斜轴Y m2O m X m Y m平面之间夹角α。
B、在机床各平动轴和回转轴的计算公式基础上,结合西门子数控系统或其他数控系统提供的宏变量(如R参数)功能,基于VC++6.0平台利用C++语言开发出一种适用于非正交摆头转台类五轴控机床的后置处理软件;其中在该软件界面中,刀具摆长可以作为宏变量或数值输入,倾斜轴Y轴与OXZ平面之间夹角作为数值输入,工件原点在机床坐标系中的偏置量作为数值输入。
在步骤A所得到的机床各平动轴和回转轴的计算公式基础上,结合西门子数控系统或其他数控系统所提供的宏变量(如R参数)功能,基于VC++6.0平台利用C++语言开发出一种适用于非正交摆头转台类五轴控机床的后置处理软件。该软件界面主要包括文件输入输出设置、刀具摆长的设置、倾斜轴YOXZ平面夹角、工件原点在机床坐标系中的位置、回转轴的运动范围、后置处理命令按钮等,如图4所示。此外,刀具摆长既可以输入数值也可以输入宏变量。对于不带“RTCP/RPCP”功能的五轴数控机床,若工件坐标原点不在A转台回转中心时,则可以根据工件原点在机床坐标系中的实际位置进行相应的补偿,避免了重新返回CAM软件中进行数控编程;若刀具长度发生变化,通过采用“宏变量后置处理”命令则可将刀具摆长作为宏变量写入数控加工程序,通过直接修改数控程序中宏变量的数值即可实现五轴刀具长度补偿。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。

Claims (3)

1.一种非正交摆头转台类五轴机床后置处理方法,其特征在于,包括以下步骤:
A、分析以平动轴Y和回转轴B作为倾斜轴的摆头转台类五轴机床的结构特征,分别将倾斜轴Y轴与OXZ平面之间夹角和刀具摆长作为变量,建立前置刀位数据与该类机床各坐标轴之间的运动变换方程,并推导出各坐标轴的计算公式;
B、在机床各平动轴和回转轴的计算公式基础上,结合西门子数控系统或其他数控系统提供的宏变量(如R参数)功能,基于VC++6.0平台利用C++语言开发出一种适用于非正交摆头转台类五轴控机床的后置处理软件;其中在该软件界面中,刀具摆长可以作为宏变量或数值输入,倾斜轴Y轴与OXZ平面之间夹角作为数值输入,工件原点在机床坐标系中的偏置量作为数值输入。
2.根据权利要求1所述的一种非正交摆头转台类五轴机床后置处理方法,其特征在于,步骤A中,首先分析该非正交摆头转台类五轴机床的结构特点,进而得到该类五轴机床各坐标轴之间的运动关系,即机床运动链,它由回转工作台、平动工作台、床身、主轴和刀具等单元按顺序串联而成;
其次,在该非正交摆头转台五轴机床上建立各相关坐标系,如工件坐标系、机床坐标系、倾斜轴的坐标系等,根据该机床运动链和各坐标系,建立回转轴和平动轴的运动变换方程
(1)
(2)
式中, T R 分别为平移和回转运动的齐次变换矩阵
是由编程坐标系O w X w Y w Z w到机床坐标系O m X m Y m Z m的变换矩阵
最后,通过求解回转轴和平动轴的运动变换方程式(1)和(2)来确定机床各回转轴和平动轴的数值,由式(1)可以得到
(3)
若该机床B轴运动范围(0.5π,π],A轴运动范围[-π,π],则;由式(3)可以得到
(4)
式中,
结合式(4),通过求解式(2)可以得到各平动轴的计算公式
(5)
因此,由式(4)和式(5)就可以确定机床各回转轴和平动轴的数值,即对前置刀位数据进行后置处理而得到机床可以识别执行的数控加工程序,其中机床平动轴计算公式中包含刀具摆长变量L、工件原点偏置变量d x d y d z 和倾斜轴Y m2O m X m Y m平面之间夹角α。
3.根据权利要求1所述的一种非正交摆头转台类五轴机床后置处理方法,其特征在于,步骤B中,刀具摆长既可以输入数值也可以输入宏变量;对于不带“RTCP/RPCP”功能的五轴数控机床,若工件坐标原点不在A转台回转中心时,则可以根据工件原点在机床坐标系中的实际位置进行相应的补偿,避免了重新返回CAM软件中进行数控编程;若刀具长度发生变化,通过采用“宏变量后置处理”命令则可将刀具摆长作为宏变量写入数控加工程序,通过直接修改数控程序中宏变量的数值即可实现五轴刀具长度补偿。
CN201510839453.8A 2015-11-27 2015-11-27 一种非正交摆头转台类五轴机床后置处理方法 Active CN105302070B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510839453.8A CN105302070B (zh) 2015-11-27 2015-11-27 一种非正交摆头转台类五轴机床后置处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510839453.8A CN105302070B (zh) 2015-11-27 2015-11-27 一种非正交摆头转台类五轴机床后置处理方法

Publications (2)

Publication Number Publication Date
CN105302070A true CN105302070A (zh) 2016-02-03
CN105302070B CN105302070B (zh) 2018-01-02

Family

ID=55199458

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510839453.8A Active CN105302070B (zh) 2015-11-27 2015-11-27 一种非正交摆头转台类五轴机床后置处理方法

Country Status (1)

Country Link
CN (1) CN105302070B (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108345266A (zh) * 2017-01-24 2018-07-31 常州数控技术研究所 一种五轴数控机床数控程序生成方法
CN108490872A (zh) * 2018-01-31 2018-09-04 深圳市拓智者科技有限公司 一种五轴rtcp测定方法
CN108549323A (zh) * 2018-04-18 2018-09-18 西南交通大学 一种五轴数控工具磨床通用后置处理方法
CN108563186A (zh) * 2018-04-02 2018-09-21 西南交通大学 一种五轴球头铣削几何误差补偿方法
CN109277658A (zh) * 2018-12-10 2019-01-29 航天科工哈尔滨风华有限公司 一种多轴电火花机床加工闭式叶轮轨迹后置处理方法
CN109696884A (zh) * 2019-01-25 2019-04-30 山东大学 一种非正交双转台型五轴机床后置处理方法
CN111367237A (zh) * 2020-04-10 2020-07-03 山东理工大学 一种非正交五轴立卧转换数控机床后处理方法
CN111413924A (zh) * 2020-04-10 2020-07-14 山东理工大学 装万向角度头非正交五轴数控机床后处理方法
CN112363453A (zh) * 2020-10-19 2021-02-12 上海航天设备制造总厂有限公司 一种利用90°侧铣头实现大型非正交角度异型深腔类零件的加工方法
WO2022110521A1 (zh) * 2020-11-30 2022-06-02 北京工业大学 一种应用于特殊七轴五联动形式机床的后置处理方法
CN116000699A (zh) * 2023-03-28 2023-04-25 中科航迈数控软件(深圳)有限公司 一种适用于双摆头的快速倾斜坐标系定义方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001162586A (ja) * 1999-12-03 2001-06-19 Nippon Reliance Kk ロータリーカッタの制御装置
CN104460515A (zh) * 2013-09-25 2015-03-25 山东理工大学 一种基于后置处理五轴刀具长度补偿方法
CN104678888A (zh) * 2013-12-02 2015-06-03 昌河飞机工业(集团)有限责任公司 基于机床运动约束环形刀多轴加工复杂曲面刀轨生成方法
CN204658052U (zh) * 2015-06-03 2015-09-23 山东永华机械有限公司 五轴联动机床摇篮转台及a/c摆头的防a轴跌落装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001162586A (ja) * 1999-12-03 2001-06-19 Nippon Reliance Kk ロータリーカッタの制御装置
CN104460515A (zh) * 2013-09-25 2015-03-25 山东理工大学 一种基于后置处理五轴刀具长度补偿方法
CN104678888A (zh) * 2013-12-02 2015-06-03 昌河飞机工业(集团)有限责任公司 基于机床运动约束环形刀多轴加工复杂曲面刀轨生成方法
CN204658052U (zh) * 2015-06-03 2015-09-23 山东永华机械有限公司 五轴联动机床摇篮转台及a/c摆头的防a轴跌落装置

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
CHAOLIN 等: "A five-axis CNC machining method of orthogonal variable transmission ratio face gear", 《JOURNAL OF ADVANCED MECHANICAL DESIGN, SYSTEMS, AND MANUFACTURING》 *
周续等: "非正交双转台五轴机床后置处理通用方法", 《机械工程学报》 *
左晓明等: "非正交双转台加工中心五轴定位铣削的后置处理与仿真验证", 《机械科学与技术》 *
葛振红: "非正交五轴联动数控机床后置处理算法", 《机械设计与研究》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108345266A (zh) * 2017-01-24 2018-07-31 常州数控技术研究所 一种五轴数控机床数控程序生成方法
CN108490872A (zh) * 2018-01-31 2018-09-04 深圳市拓智者科技有限公司 一种五轴rtcp测定方法
CN108563186B (zh) * 2018-04-02 2020-09-15 西南交通大学 一种五轴球头铣削几何误差补偿方法
CN108563186A (zh) * 2018-04-02 2018-09-21 西南交通大学 一种五轴球头铣削几何误差补偿方法
CN108549323A (zh) * 2018-04-18 2018-09-18 西南交通大学 一种五轴数控工具磨床通用后置处理方法
CN109277658A (zh) * 2018-12-10 2019-01-29 航天科工哈尔滨风华有限公司 一种多轴电火花机床加工闭式叶轮轨迹后置处理方法
CN109696884A (zh) * 2019-01-25 2019-04-30 山东大学 一种非正交双转台型五轴机床后置处理方法
CN111367237A (zh) * 2020-04-10 2020-07-03 山东理工大学 一种非正交五轴立卧转换数控机床后处理方法
CN111413924A (zh) * 2020-04-10 2020-07-14 山东理工大学 装万向角度头非正交五轴数控机床后处理方法
CN112363453A (zh) * 2020-10-19 2021-02-12 上海航天设备制造总厂有限公司 一种利用90°侧铣头实现大型非正交角度异型深腔类零件的加工方法
CN112363453B (zh) * 2020-10-19 2022-10-18 上海航天设备制造总厂有限公司 一种利用90°侧铣头实现大型非正交角度异型深腔类零件的加工方法
WO2022110521A1 (zh) * 2020-11-30 2022-06-02 北京工业大学 一种应用于特殊七轴五联动形式机床的后置处理方法
CN116000699A (zh) * 2023-03-28 2023-04-25 中科航迈数控软件(深圳)有限公司 一种适用于双摆头的快速倾斜坐标系定义方法及装置

Also Published As

Publication number Publication date
CN105302070B (zh) 2018-01-02

Similar Documents

Publication Publication Date Title
CN105302070A (zh) 一种非正交摆头转台类五轴机床后置处理方法
CN101859126B (zh) 五坐标数控加工刀具路径转换方法
CN104460515B (zh) 一种基于后置处理五轴刀具长度补偿方法
CN103048947B (zh) 五坐标联动机床数控程序后置处理方法
WO2022174657A1 (zh) 一种五轴数控机床平动轴几何误差补偿数据快速优化方法
CN108829038B (zh) 一种刀尖运动轨迹控制算法
CN101710236B (zh) 大型球体空间曲面上多个空间孔焊接坡口的数控加工方法
CN104102171A (zh) 一种球头刀多轴加工刀轴矢量优化方法
US20190271965A1 (en) Numerical controller and numerical control method
CN108549319B (zh) 一种双转台五轴数控机床通用后处理方法
CN111367237A (zh) 一种非正交五轴立卧转换数控机床后处理方法
CN106959664B (zh) 基于五轴双转台在线非线性误差补偿方法
CN111413924A (zh) 装万向角度头非正交五轴数控机床后处理方法
She et al. Development of a five-axis postprocessor system with a nutating head
CN108608425A (zh) 六轴工业机器人铣削加工离线编程方法及系统
Li et al. Development of post-processing system for three types of five-axis machine tools based on solid model
Liu et al. Analysis on influence of perpendicularity error of five axis NC machine tool error modeling accuracy and complexity
Mohamed et al. A review on intelligence STEP-NC data model and function blocks CNC machining protocol
Martinov et al. The approach of creating a particular postprocessor and using CNC measuring cycles
Solvang et al. STEP-NC based industrial robot CAM system
CN107942942A (zh) 一种机床设备应用于相交倾斜面的倾斜坐标系建立方法
Li et al. Design of a generic Postprocessor Based on UG
Lang et al. An efficient motion trajectory planning method in CNC system
Tang A five-axis CNC machine postprocessor based on inverse kinematics transformation
JPH05108134A (ja) 主軸の割出対応の座標変換方法

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220509

Address after: 255000 No. 2989, Zunxian Road, high tech Zone, Zibo City, Shandong Province

Patentee after: SHANDONG SMA PHARMATECH CO.,LTD.

Address before: 255086 room 313, block a, Gaochuang Park, high tech Zone, Zibo City, Shandong Province

Patentee before: Shandong University of Technology