CN109039166B - 一种永磁同步直线伺服系统速度环pi-ip控制参数自校正方法 - Google Patents
一种永磁同步直线伺服系统速度环pi-ip控制参数自校正方法 Download PDFInfo
- Publication number
- CN109039166B CN109039166B CN201810828582.0A CN201810828582A CN109039166B CN 109039166 B CN109039166 B CN 109039166B CN 201810828582 A CN201810828582 A CN 201810828582A CN 109039166 B CN109039166 B CN 109039166B
- Authority
- CN
- China
- Prior art keywords
- speed
- controller
- permanent magnet
- magnet synchronous
- servo system
- 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
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P6/00—Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
- H02P6/006—Controlling linear motors
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
- H02P21/18—Estimation of position or speed
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P25/00—Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
- H02P25/02—Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details characterised by the kind of motor
- H02P25/06—Linear motors
- H02P25/064—Linear motors of the synchronous type
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Ac Motors In General (AREA)
Abstract
本发明公开了一种永磁同步直线伺服系统速度环PI‑IP控制参数自校正方法,采用PI‑IP控制器,并实时自动校正PI‑IP控制器的参数,实现对永磁同步直线伺服系统的高性能速度控制,该方法包括以下步骤:S1、提取永磁同步直线伺服系统的推力电流指令与直线速度反馈,实时辨识速度环被控模型参数;S2、基于速度环被控模型,预测k+j时刻永磁同步直线伺服系统的速度输出;建立李雅普诺夫评价指标,判断速度跟踪性能;S3、简化李雅普诺夫评价指标增量函数,在稳定条件下得到PI‑IP控制器的参数在线优化结果,实现速度环PI‑IP控制器的控制参数自校正。本发明的方法利用PI‑IP控制器取代传统的PI或IP控制器,并实时校正控制器参数,具有控制结构简单、抗扰动能力强和速度响应快等优点。
Description
技术领域
本发明涉及高频响永磁同步直线伺服系统的技术领域,尤其涉及一种永磁同步直线伺服系统速度环PI-IP控制参数自校正方法。
背景技术
永磁同步直线伺服系统结构简单,不需要中间传动环节就能直接实现直线运动,具有相对小的负载惯量及高动态响应的优点,广泛应用于数控机床、半导体芯片制造及精密仪器等领域。永磁同步直线伺服系统的控制性能不仅取决于其硬件制造水平,而且决定于直线伺服驱动中所采用的控制策略和设置的控制参数。只有当直线伺服驱动采用的控制策略、设置的控制参数和永磁同步直线伺服系统自身固有特性之间形成良好匹配时,永磁同步直线伺服系统才能处于最优工作状态。
在永磁同步直线电机运行过程中,速度指令可能需要做出调整,针对这样的调整,直线伺服驱动需要具备良好的瞬态响应跟踪。当速度指令恒定时,针对不同的运行工况,直线伺服驱动需要具备较强的抗扰动能力。这样直线伺服驱动使用PI或IP控制器很难同时满足瞬态响应和抗扰动能力的要求。大量研究表明,PI-IP控制器综合了PI和IP控制器的优缺点,能够在不影响永磁同步直线伺服系统闭环稳定性的基础上,有效地提高永磁同步直线伺服系统的动态性能。但是PI-IP控制器存在更多的控制参数需要实时调整,为了满足永磁同步直线伺服系统高速高精的发展趋势,需要探求一种高效的永磁同步直线伺服系统速度环PI-IP控制器的参数自校正方法和途径。
一般而言,可将控制器参数自校正方法分为以下两类:一类是基于规则的自校正方法,如模糊PID,神经网络等,文献(A.Rohan,F.Asghar,S.H.Kim,Design of Fuzzy LogicTuned PID Controller for Electric Vehicle based on IPMSM Using Flux-weakening[J],Journal of Electrical Engineering and Technology,2018,13(1):451-459提出一种随机模糊PID的自校正方法。这类算法的计算量较大,不能满足伺服系统实时性需求,同时不恰当的控制参数初始值会使得在线校正过程陷入局部最优,无法保证最优的实时控制效果。另一类是基于模型的自校正方法,文献(Yasuki Kansha,Li Jia,Min-SenChiu.Self-tuning PID controllers based on the Lyapunov approach[J].ChemicalEngineering Science,2008,63(10):2732-2740)提出一种基于模型的PID控制参数自校正方法。基于模型的自校正方法算法简单,稳定性好,但是依赖被控模型结构和参数的辨识精度。鉴于上述两类自校正方法的优缺点,本发明拟采用改进的递推经验频率参数估计法对速度环被控模型参数进行在线高精辨识,同时提出一种李雅普诺夫控制方法来实现速度环PI-IP控制器的参数自动校正。
发明内容
本发明要解决的技术问题在于针对现有技术中的缺陷,提供一种永磁同步直线伺服系统速度环PI-IP控制参数自校正方法,该控制方法能适应永磁同步直线电机高频响特性,快速跟踪系统指令,也能适应具有负载质量、负载力等非线性特点的高速高精应用场合。
本发明解决其技术问题所采用的技术方案是:
本发明提供一种永磁同步直线伺服系统速度环PI-IP控制参数自校正方法,通过在永磁同步直线伺服系统中采用PI-IP控制器,并实时自动校正PI-IP控制器的参数,实现对永磁同步直线伺服系统的高性能速度控制,该方法包括以下步骤:
S1、提取永磁同步直线伺服系统的推力电流指令与直线速度反馈,建立速度环被控模型,并实时辨识速度环被控模型参数;
S2、基于速度环被控模型,预测k+j时刻永磁同步直线伺服系统的速度输出;建立李雅普诺夫评价指标,判断速度跟踪性能;
S3、简化李雅普诺夫评价指标增量函数,在稳定条件下得到PI-IP控制器的参数在线优化结果,实现速度环PI-IP控制器的控制参数自校正。
进一步地,本发明的该方法中实时辨识速度环被控模型参数的方法为:
永磁同步直线伺服系统中,速度环被控模型的参数通过改进的递推经验频率参数估计法来在线辨识,其中,速度环被控对象模型的离散表达式为:
其在线辨识过程可以通过如下方程组来进行:
H(k)=H(k-1)-m-1(k)H(k-1)
×[ψ(k)βT(k)Φ(k-1)
+ΦT(k-1)β(k)ψT(k)]H(k-1)
+l-1(k)m-1(k)H(k-1)
×[ψ(k)H(k-1)ψT(k)ΦT(k-1)β(k)
×βT(k)Φ(k-1)-σ(k)ψ(k)ψT(k)]H(k-1)
Φ(k)=[ψ(k-1),ψ(k)]T
Y(k)=[ωf(k-1),ωf(k)]T
β(k)=[1,1]T×[x(k)x(k-1)+x(k-1)x(k-2)+x(k)x(k-2)]
σ(k)=x*(k)-βT(k)Φ(k-1)H(k-1)ΦT(k-1)β(k)
l(k)=1+ψT(k)H(k-1)ΦT(k-1)β(k)
m(k)=l(k)+l-1(k)σ(k)ψT(k)H(k-1)ψ(k)
x*(k)=x2(k)+x2(k-1)+x2(k-2)
其中,Φ(k)为输入输出历史与当前数据,Y(k)为输出历史与当前数据,β(k)为加权向量,σ(k)、l(k)、m(k)和x*(k)均为参与递推运算的中间变量。
进一步地,本发明的该方法中采用李雅普诺夫来实现PI-IP控制器的参数最优校正,其具体步骤如下:
(1)预测输出:结合永磁同步直线伺服系统速度环被控模型参数,预测k+1时刻系统的速度输出,得到预测速度误差,其推导过程如下:
e(k+1)=ωr(k+1)-ωf(k+1)
(2)建立李雅普诺夫评价指标:通过更新速度环PI-IP控制参数,使得速度预测输出和速度指令保持一致,其评价指标函数表示为:
其中,λ是一个正实数;
(3)PI-IP控制参数在线自校正:永磁同步直线伺服系统的加速度定义为α,优化李雅普诺夫评价指标增量函数,得到PI-IP控制参数在线学习过程;
其中,η(k)=[η1(k),η2(k),η3(k)]T为参数矢量,其余中间变量如下:
eu(k)=[e(k),ωr(k)-ωr(k-1),ωf(k-1)-ωf(k)]
进一步地,本发明的该方法中PI-IP控制参数在线校正结果通过如下公式计算得到:
进一步地,本发明的该方法中的PI-IP控制器表达为以下增量模式:
其中,kv,ki和kα为PI-IP控制器的控制参数且kα∈[0,1];当kα=1,PI-IP控制器退化为PI控制器,当kα=0,PI-IP控制器退化为IP控制器。
本发明产生的有益效果是:本发明的永磁同步直线伺服系统速度环PI-IP控制参数自校正方法,1、在被控对象模型结构已知的情况下,采用改进的递推经验频率参数估计法,直接根据当前和过去的输入输出数据估计被控模型动态参数,算法实时性强,辨识精度高。2、本发明进一步扩大李雅普诺夫方法的适用范围,针对PI-IP控制器独特的控制结构,在系统全局稳定情况下,李雅普诺夫法能有效地调整控制参数,同时也保持了PI和IP控制器的优良特性。3、本发明可满足高频响永磁同步直线伺服系统的指令快速调整,也能适应负载质量和负载力等非线性特性应用场合。工程人员不需要根据手动设定和调节控制参数,系统自动完成速度环PI-IP控制参数自校正。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1是本发明实施例的永磁同步直线伺服系统矢量控制结构示意图。
图2是本发明实施例的PI-IP控制器结构图。
图3是本发明实施例的控制参数自校正原理结构示意图。
图4是本发明实施例的控制参数自校正流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
图1为本发明永磁同步直线伺服系统矢量控制结构示意图。在实际工程应用中,通常采用id=0来实现电流的近似解耦。在图2中,通过对永磁同步直线伺服系统的速度环被控对象模型进行离散化处理,可以得到速度环被控对象的二阶离散模型:
其中,kv,ki和kα为PI-IP控制器的控制参数且kα∈[0,1]。当kα=1,PI-IP控制器退化为PI控制器,当kα=0,PI-IP控制器退化为IP控制器。e(k)=ωr(k)-ωf(k),ωr(k)为直线速度指令,e(k)为实际速度误差。
自校正方法采用李雅普诺夫控制方法,其基本原理如图3所示。在得到速度环被控模型参数后,对k+1时刻的速度输出进行预测,从而对系统预测输出误差进行在线评价,并根据评价结果推导出PI-IP控制参数的迭代公式,从而实现PI-IP控制参数的在线校正,满足永磁同步直线伺服系统瞬态响应和抗扰动能力的要求。
基于李雅普诺夫的永磁同步直线伺服速度环PI-IP控制参数的自校正流程图如图4所示,主要有以下几步:
第一步,首先需要实时提取永磁同步直线伺服速度环中的实际直线速度ωf和推力电流指令作为改进后的递推经验频率参数估计方法的输入输出数据。通过实时在线辨识得到所需的被控模型参数和改进后的递推经验频率参数估计方法如下:
H(k)=H(k-1)-m-1(k)H(k-1)
×[ψ(k)βT(k)Φ(k-1)
+ΦT(k-1)β(k)ψT(k)]H(k-1)
+l-1(k)m-1(k)H(k-1)
×[ψ(k)H(k-1)ψT(k)ΦT(k-1)β(k)
×βT(k)Φ(k-1)-σ(k)ψ(k)ψT(k)]H(k-1)
Φ(k)=[ψ(k-1),ψ(k)]T
Y(k)=[ωf(k-1),ωf(k)]T
β(k)=[1,1]T×[x(k)x(k-1)+x(k-1)x(k-2)+x(k)x(k-2)]
σ(k)=x*(k)-βT(k)Φ(k-1)H(k-1)ΦT(k-1)β(k)
l(k)=1+ψT(k)H(k-1)ΦT(k-1)β(k)
m(k)=l(k)+l-1(k)σ(k)ψT(k)H(k-1)ψ(k)
x*(k)=x2(k)+x2(k-1)+x2(k-2)
其中,其中,Φ(k)为输入输出历史与当前数据,Y(k)为输出历史与当前数据,β(k)为加权向量,σ(k)、l(k)、m(k)和x*(k)均为参与递推运算的中间变量,ωf为直线速度反馈,为推力电流指令。
第二步,结合速度环被控模型参数,预测k+1时刻系统的速度输出,判断速度是否跟踪良好,设其李雅普诺夫评价指标如下:
ωr(k+1)为设定速度;
定义er(k+1)=er(k)+Δer(k+1),李雅普诺夫评价指标增量函数ΔV(k)可通过如下方式获得:
第三步,当永磁同步直线电机以加速度α运行的时候,上式中的Δer(k+1)可以表示为:
其中,上述式子中关键函数可表示如下:
在李雅普诺夫评价指标增量函数基础上,可得到全局收敛的PI-IP控制参数在线学习过程。
当调整参数η(k)=[η1(k),η2(k),η3(k)]T满足上面的式子时,可以保证ΔV(k)在全局范围内总为负数,即表示速度的跟踪误差会趋近于零,系统是稳定的。得到PI-IP控制器的参数在线校正结果如下:
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (3)
1.一种永磁同步直线伺服系统速度环PI-IP控制参数自校正方法,其特征在于,通过在永磁同步直线伺服系统中采用PI-IP控制器,并实时自动校正PI-IP控制器的参数,实现对永磁同步直线伺服系统的高性能速度控制,该方法包括以下步骤:
S1、提取永磁同步直线伺服系统的推力电流指令与直线速度反馈,建立速度环被控模型,并实时辨识速度环被控模型参数;
S2、基于速度环被控模型,预测k+j时刻永磁同步直线伺服系统的速度输出;建立李雅普诺夫评价指标,判断速度跟踪性能;
S3、简化李雅普诺夫评价指标增量函数,在稳定条件下得到PI-IP控制器的参数在线优化结果,实现速度环PI-IP控制器的控制参数自校正;
该方法中实时辨识速度环被控模型参数的方法为:
永磁同步直线伺服系统中,速度环被控模型的参数通过改进的递推经验频率参数估计法来在线辨识,其中,速度环被控对象模型的离散表达式为:
其在线辨识过程可以通过如下方程组来进行:
H(k)=H(k-1)-m-1(k)H(k-1)×[ψ(k)βT(k)Φ(k-1)+ΦT(k-1)β(k)ψT(k)]H(k-1)+l-1(k)m-1(k)H(k-1)×[ψ(k)H(k-1)ψT(k)ΦT(k-1)β(k)×βT(k)Φ(k-1)-σ(k)ψ(k)ψT(k)]H(k-1)
Φ(k)=[ψ(k-1),ψ(k)]T
Y(k)=[ωf(k-1),ωf(k)]T
β(k)=[1,1]T×[x(k)x(k-1)+x(k-1)x(k-2)+x(k)x(k-2)]
σ(k)=x*(k)-βT(k)Φ(k-1)H(k-1)ΦT(k-1)β(k)
l(k)=1+ψT(k)H(k-1)ΦT(k-1)β(k)
m(k)=l(k)+l-1(k)σ(k)ψT(k)H(k-1)ψ(k)
x*(k)=x2(k)+x2(k-1)+x2(k-2)
该方法中采用李雅普诺夫来实现PI-IP控制器的参数最优校正,其具体步骤如下:
(1)预测输出:结合永磁同步直线伺服系统速度环被控模型参数,预测k+1时刻系统的速度输出,得到预测速度误差,其推导过程如下:
e(k+1)=ωr(k+1)-ωf(k+1)
(2)建立李雅普诺夫评价指标:通过更新速度环PI-IP控制参数,使得速度预测输出和速度指令保持一致,其评价指标函数表示为:
其中,λ是一个正实数;
(3)PI-IP控制参数在线自校正:永磁同步直线伺服系统的加速度定义为α,优化李雅普诺夫评价指标增量函数,得到PI-IP控制参数在线学习过程;
其中,η(k)=[η1(k),η2(k),η3(k)]T为参数矢量,其余中间变量如下:
eu(k)=[e(k),ωr(k)-ωr(k-1),ωf(k-1)-ωf(k)]
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810828582.0A CN109039166B (zh) | 2018-07-25 | 2018-07-25 | 一种永磁同步直线伺服系统速度环pi-ip控制参数自校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810828582.0A CN109039166B (zh) | 2018-07-25 | 2018-07-25 | 一种永磁同步直线伺服系统速度环pi-ip控制参数自校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109039166A CN109039166A (zh) | 2018-12-18 |
CN109039166B true CN109039166B (zh) | 2022-03-22 |
Family
ID=64646221
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810828582.0A Active CN109039166B (zh) | 2018-07-25 | 2018-07-25 | 一种永磁同步直线伺服系统速度环pi-ip控制参数自校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109039166B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109951115B (zh) * | 2019-04-08 | 2021-08-31 | 东莞市诺必然智能科技有限公司 | 一种直流无刷电机转速自适应控制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1623273A (zh) * | 2002-03-20 | 2005-06-01 | 株式会社安川电机 | 控制常数调整装置 |
CN103124158A (zh) * | 2012-12-06 | 2013-05-29 | 华中科技大学 | 基于分数阶的永磁同步电机速度环控制参数的自整定方法 |
CN104270055A (zh) * | 2014-09-24 | 2015-01-07 | 武汉科技大学 | 一种交流伺服系统位置环ip控制器的参数自校正方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101988088B1 (ko) * | 2016-10-31 | 2019-06-12 | 현대자동차주식회사 | 모터 구동 제어 방법, 시스템 및 이를 적용한 연료전지 시스템의 압축기 구동 제어 방법 |
-
2018
- 2018-07-25 CN CN201810828582.0A patent/CN109039166B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1623273A (zh) * | 2002-03-20 | 2005-06-01 | 株式会社安川电机 | 控制常数调整装置 |
CN103124158A (zh) * | 2012-12-06 | 2013-05-29 | 华中科技大学 | 基于分数阶的永磁同步电机速度环控制参数的自整定方法 |
CN104270055A (zh) * | 2014-09-24 | 2015-01-07 | 武汉科技大学 | 一种交流伺服系统位置环ip控制器的参数自校正方法 |
Non-Patent Citations (1)
Title |
---|
基于RBF神经网络的永磁同步电机速度PI-IP控制;周佳等;《组合机床与自动化加工技术》;20170120(第01期);第116-118、122页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109039166A (zh) | 2018-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sun et al. | Neural network-based adaptive controller design of robotic manipulators with an observer | |
CN108303885A (zh) | 一种基于干扰观测器的电机位置伺服系统自适应控制方法 | |
CN109491248B (zh) | 基于rbf-arx模型和拉盖尔函数的磁悬浮球位置预测控制方法 | |
CN111965976B (zh) | 基于神经网络观测器的机器人关节滑模控制方法及系统 | |
Yu et al. | Stable neural PID anti-swing control for an overhead crane | |
Bartolini et al. | Higher-order sliding modes for the output-feedback control of nonlinear uncertain systems | |
CN109885077A (zh) | 一种四旋翼飞行器姿态控制方法及控制器 | |
CN110209122B (zh) | 一种多轴运动平台的控制方法、装置、介质及设备 | |
CN115890668A (zh) | 一种机器人关节模组分散式优化学习控制方法及系统 | |
CN112769364A (zh) | 一种直流电机伺服系统的快速自适应抗扰控制方法 | |
CN111459031A (zh) | 一种面向学习的扰动观测器设计方法 | |
CN111930010A (zh) | 一种基于lstm网络的通用mfa控制器设计方法 | |
CN109039166B (zh) | 一种永磁同步直线伺服系统速度环pi-ip控制参数自校正方法 | |
Lin et al. | Data‐driven set‐point tuning of model‐free adaptive control | |
Perng et al. | Design of robust PI control systems based on sensitivity analysis and genetic algorithms | |
Wang et al. | Discrete-time sliding mode control with adaptive reaching law via implicit euler method | |
Gao et al. | Event-triggered based model-free adaptive sliding mode constrained control for nonlinear discrete-time systems | |
CN116442223A (zh) | 一种机械手系统轨迹跟踪的非线性动态控制器设计方法 | |
CN105739311B (zh) | 基于预设回声状态网络的机电伺服系统受限控制方法 | |
Menées et al. | Fuzzy model predictive control for nonlinear processes | |
CN110209055B (zh) | 基于参考模型和扰动观测的二阶系统控制器及控制方法 | |
Abdelouaheb et al. | Synergetic adaptive fuzzy control for a class of nonlinear discrete-time systems | |
Zhang et al. | An online control approach for forging machine using reinforcement learning and taboo search | |
Huang et al. | Performance optimisation of discrete time linear active disturbance rejection control approach | |
CN111679578A (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: 20220817 Address after: Building 12, No. 6, Zhongnan South Road, Shangsha, Chang'an Town, Dongguan City, Guangdong Province, 523846 Patentee after: Guangdong Samson Technology Co.,Ltd. Address before: 430081 No. 947 Heping Avenue, Qingshan District, Hubei, Wuhan Patentee before: WUHAN University OF SCIENCE AND TECHNOLOGY |