CN107505918B - 一种切割机的速度规划方法 - Google Patents
一种切割机的速度规划方法 Download PDFInfo
- Publication number
- CN107505918B CN107505918B CN201710694486.7A CN201710694486A CN107505918B CN 107505918 B CN107505918 B CN 107505918B CN 201710694486 A CN201710694486 A CN 201710694486A CN 107505918 B CN107505918 B CN 107505918B
- Authority
- CN
- China
- Prior art keywords
- acceleration
- speed
- velocity
- deceleration
- corner
- 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
Links
- 238000005520 cutting process Methods 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 title claims abstract description 38
- 230000001133 acceleration Effects 0.000 claims abstract description 76
- 238000009795 derivation Methods 0.000 claims abstract description 5
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 230000003044 adaptive effect Effects 0.000 claims description 3
- 210000003127 knee Anatomy 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 230000003247 decreasing effect Effects 0.000 claims description 2
- 239000004576 sand Substances 0.000 claims description 2
- 230000001788 irregular Effects 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/18—Numerical 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/416—Numerical 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 of velocity, acceleration or deceleration
- G05B19/4163—Adaptive control of feed or cutting velocity
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/36—Nc in input of data, input key till input tape
- G05B2219/36521—Select by combination of detected force, acceleration, speed, work rate
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
一种切割机的速度规划方法,所述方法包括以下步骤:1)速度规划方案如下:采用类S型五段加减速控制方法,确定相应的阶梯步长表;2)可变加速度的实现如下:确定拐点速度、拐角大小以及加速度三者之间的关系,在经过轨迹的不同拐角时,根据当前实时速度调整加速度;3)自适应拐点速度的计算,根据角速度求解公式和速度反向推导法确定最终的拐角速度。本发明提供一种在切割过程中有效缩短运动加减速时间和提高拐角速度的切割机的速度规划方法。
Description
技术领域
本发明涉及一种运动速度规划的方法,尤其是一种切割机的速度规划方法,属于运动控制领域。
背景技术
中国有着巨大的消费市场,人们对各种产品的需求不断扩大,在对产品的包装盒、标签以及其他材料加工时,切割机是最为理想的选择,因为切割机在稳定性、精度以及工作效率上能较好满足生产的需求。大部分的切割机在刀具等相关工具更换上很方面,能够根据实际情况满足多种材料的裁切。
目前,国内市场上存在大量中低端切割机,在切割不规则物件时,由于运行轨迹存在不同大小的夹角,切割机在加减速时间以及拐角转弯的速度受到较为严重的影响。而现有市场上的切割机在处理这个拐角问题时,仅仅考虑了角度大小的问题,而未结合当前运行速度,并且在加减速时间上没有做到很好的优化。因此在这些方面还存在很大的提升空间。
发明内容
为了解决现有嵌入式切割机在对不规则物件切割时运行速度较慢的问题,本发明提供一种类S型五段曲线速度控制和可变加速度的控制方法,可有效提升切割机在对不规则物件切割时的速度和稳定性。
本发明解决其技术问题所采用的技术方案如下:
一种切割机的速度规划方法,所述方法包括以下步骤:
1)采用类S型五段加减速控制方法,确定相应的阶梯步长表
类S型的五段速度分别为加加速段、匀加速段、匀速段、匀减速段、减减速段,结合加速度表和电机转动的物理特性得到相应的阶梯步长表;
2)可变加速度的实现如下:确定拐点速度、拐角大小以及加速度三者之间的关系,在经过轨迹的不同拐角时,根据当前实时速度调整加速度;
3)自适应拐点速度的计算
上位机下发要切割的目标轨迹文件由离散后的n个坐标点组成,相应坐标点为(x1,y1),(x2,y2),…,(xn-1,yn-1),(xn,yn),相邻坐标点相连形成的密集线段即为切割机需要运行切割的小线段,相邻线段之间的点即为拐点,n-1条小线段分别表示为:p1p2,p3p4,…,pn-2,pn-1pn,计算出它们与X轴方向的夹角分别为:θ1,θ2,…,θn-1;结合步骤(2)中所述,根据当前的实时速度V、加速度a以及拐角之间的关系,以计算得到在切割轨迹各个坐标拐点对应X和Y轴的最大速度(v1x,v1y),(v2x,v2y),…,(vn-1x,vn-1y);将以上求得的信息带入角度度求解公式即可求得各坐标点的拐角速度Vj;
将求得的X和Y轴的最大速度作为初始目标速度,结合步骤(1)中的阶梯步长表,进行反向推导,对轨迹重新规划,求得最终的拐角速度为V′j。
进一步,所述速度规划方法还包括以下步骤:
4)对长度小于最小可规划的线段处理
在步骤3)的速度规划中,若是遇到切割线段小于设定的最小可加减速规划的线段时,Lmin表示设定的最小线段长度,在这种情况下的速度规划应做如下处理:首先是修改当前加速度,同时将下一个离散点的目标速度作为当前线段加减速过程中的最大运行速度;若在此线段的运动过程中达到了这一最大速度,后续就进行匀速运动,不再进行加减速,直到运动到下一离散点;
5)建立加速度表和速度表
根据步骤2)、3)所述,由速度和加速度的关系以及电机的升降频特性,并且通过对系统的测试,能够确定相应的加速度表和速度表。
本发明的有益效果表现在:(1)采用类S型的五段曲线加减速控制算法,相比标准七段S型控制方法,略去减加速段和加减速段,能明显缩短加减速时间。(2)引入可变加速度控制方法,在轨迹拐点角速度计算时,根据实时速度、加速度以及拐角大小三者的关系,并且结合速度反向推导法,能快速确定最终的拐角速度。(3)根据加速度表、速度表、阶梯步长表以及各拐点角度大小,切割机能够较好的进行自适应轨迹切割运动,保证系统的稳定运行,避免了过切的现象,从而提高了生产效率。
附图说明
图1为类S型五段加减速度控制算法的速度曲线图。
图2为相邻小线段拐点速度变化及在坐标轴上的速度分量图。
具体实施方式
下面结合附图对本发明的实施方式作进一步描述。
参照图1和图2,一种切割机的速度规划方法,包括以下步骤:
1)采用类S型五段加减速控制方法,确定相应的阶梯步长表
本系统加减速的控制采用的是类S型的五段速度控制方法,分为加加速段、匀加速段、匀速段、匀减速段、减减速段。与七段标准S型速度曲线相比,略去减加速段和加减速段。这种控制方法能够平稳起步,并且以最短时间加速到最大速度,从而减少了加速(减速)时间,提高了运行效率。
为了避免切割机在对不规则物件切割时,因惯性引起的切割误差,所以需要一定的脉冲数来适应电机的速度变化,速度越大所需脉冲数越多。结合加速度表和电机转动的物理特性可得到相应的阶梯步长表。
2)可变加速度的实现如下
本发明提出一种可变加速度计算方法,主要是针对不规则小线段轨迹的裁切,在对运行轨迹的速度规划时,确定当前实时速度V和加速度a之间的函数关系,并将该函数关系式带入拐角速度计算公式即可求得相应的拐角速度。相对于固定加速度控制而言,能很好的提高运行速度。
已知刀头最高运行速度为Vmax时,可测得在一个脉冲周期T内最大加速度为aT。同理在刀头运行速度为Vmin时,可测得最大加速度为a′T,可知a′T>aT。再结合切割机平台的机械结构以及电机的型号通过现场调试,可以建立运行速度v和加速度a之间的函数关系。这样就能得到在各个时刻对应X轴和Y轴方向上电机能允许运行的最大加速度大小atx、aty。
a(v)=k*v+b (1)
如上公式所示,公式1即为拐角速度与对应合加速度的函数关系式。其中a(v)表示合加速度与当前运行速度v的关系,K是该函数关系式的斜率,b是该函数关系式的常数项。(根据拐角处速度不可突变的原理,即v=ve=vs=vj,表示上一条线段的末速度ve和下一条线段的起始速度vs以及拐角速度Vj三者在数值上相等并且等于运行合速度v)。
3)自适应拐点速度的计算
上位机下发要切割的目标轨迹文件由离散后的n个坐标点组成,相应坐标点为(x1,y1),(x2,y2),…,(xn-1,yn-1),(xn,yn),相邻坐标点相连形成的密集线段即为切割机需要运行切割的小线段。相邻线段之间的点即为拐点。n-1条小线段(此处以切割不封闭的图形为例)分别表示为:p1p2,p3p4,…,pn-2,pn-1pn,可计算出它们与X轴方向的夹角分别为:θ1,θ2,…,θn-1。结合步骤2中所述,根据当前的实时速度V、加速度a以及拐角之间的关系,可以计算得到在切割轨迹各个坐标拐点对应X和Y轴的最大速度(v1x,v1y),(v2x,v2y),…,(vn-1x,vn-1y)。将以上求得的信息带入角度度求解公式即可求得各坐标点的拐角速度Vj。
将公式(1)带入公式(6)、(7)可计算得到相应坐标点的角速度,对求得的角速度值取决对值为确保切割稳定,可求得实际的角速度为根据以上计算得到各坐标点的拐角速度Vj_1,Vj_2,...,Vj_(n-2)。
进一步,为了切割机系统的稳定和防止过切的现象,将前面求得的X和Y轴的最大速度作为初始目标速度,进行反向推导,对轨迹重新规划,求得最终的拐角速度为V′j_1,V′j_2,…,V′j_(n-2)。
4)对长度小于设定的最小可规划线段的处理
在步骤3)的速度规划中,若是遇到切割线段小于设定的最小可加减速规划的线段时,Lmin表示设定的最小线段长度,在这种情况下的速度规划应做如下处理:首先是修改当前加速度,同时将下一个离散点的目标速度作为当前线段加减速过程中的最大运行速度。若在此线段的运动过程中达到了这一最大速度,后续就进行匀速运动,不再进行加减速,直到运动到下一离散点。
5)建立加速度表和速度表
根据步骤2)、3)所述,可以由速度和加速度的关系确定,并且通过对系统的测试,能够确定相应的加速度表和速度表。
进一步,所述的运动速度规划方法还包含以下步骤:
6)在步骤4)的速度规划中,若是遇到切割线段小于设定的最小可加减速规划线段时,Lmin表示设定的最小线段长度,在这种情况下本发明对速度规划的处理如下:通过改变原运动过程中的加速度来控制,即将当前使用的加速度a改成a’,公式如下所示:
其中加速度α是根据原类S型五段加减速度规划算法得到的,t为常系数。同时将下一个离散点的目标速度作为当前加减速过程中的最大速度,如果在此线段的运动过程中达到了这一最大速度,后续就进行匀速运动,不再进行加减速,直到运动到下一离散点。
Claims (2)
1.一种切割机的速度规划方法,其特征在于:所述方法包括以下步骤:
1)采用类S型五段加减速控制方法,确定相应的阶梯步长表;
类S型的五段速度分别为加加速段、匀加速段、匀速段、匀减速段、减减速段,结合加速度表和电机转动的物理特性得到相应的阶梯步长表;
2)可变加速度的实现如下:确定拐点速度、拐角大小以及加速度三者之间的关系,在经过轨迹的不同拐角时,根据当前实时速度调整加速度;
已知刀头最高运行速度为Vmax时,可测得在一个脉冲周期T内最大加速度为aT,同理在刀头运行速度为Vmin时,可测得最大加速度为a′T,可知a′T>aT,再结合切割机平台的机械结构以及电机的型号通过现场调试,建立运行速度v和加速度a之间的函数关系,这样就能得到在各个时刻对应X轴和Y轴方向上电机能允许运行的最大加速度大小atx、aty;
a(v)=k*v+b (1)
如上公式所示,公式(1)即为拐角速度与对应合加速度的函数关系式,其中a(v)表示合加速度与当前运行速度v的关系,K是该函数关系式的斜率,b是该函数关系式的常数项,根据拐角处速度不可突变的原理,即v=ve=vs=vj,表示上一条线段的末速度ve和下一条线段的起始速度vs以及拐角速度Vj三者在数值上相等并且等于运行合速度v;
3)自适应拐点速度的计算
上位机下发要切割的目标轨迹文件由离散后的n个坐标点组成,相应坐标点为(x1,y1),(x2,y2),...,(xn-1,yn-1),(xn,yn),相邻坐标点相连形成的密集线段即为切割机需要运行切割的小线段,相邻线段之间的点即为拐点,n-1条小线段分别表示为:p1p2,p3p4,...,pn-2,pn-1pn,计算出它们与X轴方向的夹角分别为:θ1,θ2,...,θn-1;结合步骤(2)中所述,根据当前的实时速度V、加速度a以及拐角之间的关系,以计算得到在切割轨迹各个坐标拐点对应X和Y轴的最大速度(v1x,v1y),(v2x,v2y),...,(vn-1x,vn-1y);将以上求得的信息带入角度求解公式即可求得各坐标点的拐角速度Vj;
将公式(1)带入公式(6)、(7)计算得到相应坐标点的角速度,对求得的角速度值取决对值为确保切割稳定,求得实际的角速度为根据以上计算得到各坐标点的拐角速度Vj_1,Vj_2,...,Vj_(n-2);
将求得的X和Y轴的最大速度作为初始目标速度,结合步骤(1)中的阶梯步长表,进行反向推导,对轨迹重新规划,求得最终的拐角速度为V′j。
2.如权利要求1所述的一种切割机的速度规划方法,其特征在于:所述速度规划方法还包括以下步骤:
4)对长度小于最小可规划的线段处理
在步骤3)的速度规划中,若是遇到切割线段小于设定的最小可加减速规划的线段时,Lmin表示设定的最小线段长度,在这种情况下的速度规划应做如下处理:首先是修改当前加速度,通过改变原运动过程中的加速度来控制,即将当前使用的加速度a改成a′,公式如下所示:
其中,加速度α是根据原类S型五段加减速度规划算法得到的,t为常系数;
同时将下一个离散点的目标速度作为当前线段加减速过程中的最大运行速度;若在此线段的运动过程中达到了这一最大速度,后续就进行匀速运动,不再进行加减速,直到运动到下一离散点;
5)建立加速度表和速度表
根据步骤2)、3)所述,由速度和加速度的关系以及电机的升降频特性,并且通过对系统的测试,能够确定相应的加速度表和速度表。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710694486.7A CN107505918B (zh) | 2017-08-15 | 2017-08-15 | 一种切割机的速度规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710694486.7A CN107505918B (zh) | 2017-08-15 | 2017-08-15 | 一种切割机的速度规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107505918A CN107505918A (zh) | 2017-12-22 |
CN107505918B true CN107505918B (zh) | 2020-01-10 |
Family
ID=60690920
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710694486.7A Active CN107505918B (zh) | 2017-08-15 | 2017-08-15 | 一种切割机的速度规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107505918B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108153246B (zh) * | 2017-12-26 | 2020-07-10 | 哈工大机器人(合肥)国际创新研究院 | 一种基于指定速度的参数自适应s形速度规划插补方法 |
CN108388206B (zh) * | 2018-03-07 | 2020-11-24 | 深圳市汇川技术股份有限公司 | 进给速度实时动态规划方法及系统 |
CN108817695A (zh) * | 2018-07-17 | 2018-11-16 | 大族激光科技产业集团股份有限公司 | 激光切割方法、装置及数控系统 |
CN109048091B (zh) * | 2018-07-17 | 2020-12-25 | 大族激光科技产业集团股份有限公司 | 激光切割速度规划方法、装置、存储介质及计算机设备 |
CN109300158B (zh) * | 2018-08-01 | 2021-05-18 | 浙江工业大学 | 一种基于Mark点定位功能切割PVC板的方法 |
CN109634219A (zh) * | 2018-12-24 | 2019-04-16 | 杭州澳星科技有限公司 | 一种有效保护电机的平面双轴协同切割方法 |
CN111977571A (zh) * | 2019-05-21 | 2020-11-24 | 北京京东尚科信息技术有限公司 | 一种升降机构速度控制方法和装置 |
CN110286653A (zh) * | 2019-06-14 | 2019-09-27 | 杭州爱科科技股份有限公司 | 用于任意曲线运动s加减速控制的速度计算方法 |
CN111015785B (zh) * | 2019-12-27 | 2021-05-18 | 湖南鼎一致远科技发展有限公司 | 一种切刀减速方法及装置 |
CN111679633B (zh) * | 2020-06-19 | 2023-06-09 | 重庆大学 | 一种基于自抗扰的材料追剪控制方法 |
CN112327756B (zh) * | 2020-11-19 | 2021-08-17 | 杭州爱科科技股份有限公司 | 一种柔性材料轨迹数据处理方法、装置、设备及存储介质 |
CN113441848B (zh) * | 2021-06-29 | 2022-05-06 | 苏州科韵激光科技有限公司 | 一种偏光片的切割方法及切割装置 |
CN113325807B (zh) * | 2021-08-02 | 2021-10-08 | 杭州爱科科技股份有限公司 | 一种切割运动速度的控制方法、装置、设备及介质 |
CN115581617A (zh) * | 2022-10-21 | 2023-01-10 | 海南先端医疗科技有限公司 | 一种智能模拟人工捻转提插的电针灸治疗仪系统 |
CN118219055B (zh) * | 2024-05-23 | 2024-08-20 | 宁波经纬数控股份有限公司 | 一种异形图案切割速度控制方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103326646A (zh) * | 2013-05-17 | 2013-09-25 | 浙江工业大学 | 一种基于步进电机的运动控制器速度控制方法 |
CN103324141B (zh) * | 2013-06-14 | 2015-04-29 | 浙江工业大学 | 一种高精度变插补周期的多轴联动运动控制方法 |
CN106444645A (zh) * | 2016-08-17 | 2017-02-22 | 义乌朝晖智能科技有限公司 | 一种基于嵌入式裁床控制器的多轴联动运动控制方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3370040B2 (ja) * | 2000-03-02 | 2003-01-27 | 日本リライアンス株式会社 | 速度制御装置 |
ES2683177T3 (es) * | 2014-05-07 | 2018-09-25 | Fives Oto S.P.A. | Máquina para cortar un objeto en movimiento |
-
2017
- 2017-08-15 CN CN201710694486.7A patent/CN107505918B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103326646A (zh) * | 2013-05-17 | 2013-09-25 | 浙江工业大学 | 一种基于步进电机的运动控制器速度控制方法 |
CN103324141B (zh) * | 2013-06-14 | 2015-04-29 | 浙江工业大学 | 一种高精度变插补周期的多轴联动运动控制方法 |
CN106444645A (zh) * | 2016-08-17 | 2017-02-22 | 义乌朝晖智能科技有限公司 | 一种基于嵌入式裁床控制器的多轴联动运动控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107505918A (zh) | 2017-12-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107505918B (zh) | 一种切割机的速度规划方法 | |
CN107850887B (zh) | S型曲线规划方法、装置及数控机床 | |
CN104678894B (zh) | 数控加工路径的规划方法、数控加工系统和方法 | |
CN102681487B (zh) | 数控系统中操作设备的轨迹平滑方法和装置及数控机床 | |
CN101976060B (zh) | 基于机床动力学和曲线特性的nurbs插补方法 | |
CN106094737B (zh) | 一种指定加工误差条件下的数控加工速度优化控制方法 | |
CN107077126A (zh) | 刀具路径的生成方法及机床 | |
CN110618659B (zh) | 一种面向五轴线性插补的平动轴旋转轴协调运动规划控制方法 | |
CN109799778A (zh) | 一种电子凸轮曲线生成方法及相关装置 | |
CN106444645A (zh) | 一种基于嵌入式裁床控制器的多轴联动运动控制方法 | |
CN112965443B (zh) | 一种裁床拐角轨迹跟踪高精度插值控制方法 | |
CN112486101B (zh) | Nurbs曲线自适应前瞻插补方法 | |
CN105717873A (zh) | 一种基于模板缝纫机控制器的自动送料速度控制方法 | |
CN105629882A (zh) | 一种用于样条插补的三角函数速度规划方法 | |
CN111722591B (zh) | 一种商标模切机高精度联动插补的方法 | |
CN116300698A (zh) | 一种基于动态规划的螺纹切削加工方法 | |
CN115616983A (zh) | 一种五轴联动同步刀具路径插补方法和系统 | |
CN105717874B (zh) | 一种五轴数控加工奇异区域刀位点优化方法 | |
CN114415598B (zh) | 加工路径的过渡方法、装置、存储介质及计算机设备 | |
CN108098453A (zh) | 带刀长有旋转轴运动的数控机床的速度控制方法 | |
CN112068488B (zh) | 点胶机防堆胶的运动控制方法 | |
CN104750925B (zh) | 一种关于伺服冲床主轴非匀速运动曲线的分析方法 | |
TWI541623B (zh) | 五軸加工數值控制系統及其數值控制方法 | |
CN113733099A (zh) | 一种机器人平滑轨迹规划方法、计算机及存储介质 | |
CN108628255B (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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20210303 Address after: No.206, Lane 777, Qingfeng Road, Cicheng Town, Jiangbei District, Ningbo City, Zhejiang Province Patentee after: Ningbo Yongshi precision tools Co.,Ltd. Address before: 310014 Zhejiang University of Technology, 18, Chao Wang Road, Xiacheng District, Hangzhou, Zhejiang Patentee before: ZHEJIANG University OF TECHNOLOGY |