CN107191483B - 一种三自由度三极混合磁轴承的设计方法 - Google Patents

一种三自由度三极混合磁轴承的设计方法 Download PDF

Info

Publication number
CN107191483B
CN107191483B CN201710284370.6A CN201710284370A CN107191483B CN 107191483 B CN107191483 B CN 107191483B CN 201710284370 A CN201710284370 A CN 201710284370A CN 107191483 B CN107191483 B CN 107191483B
Authority
CN
China
Prior art keywords
radial
axial
magnetic
rmax
bearing capacity
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
Application number
CN201710284370.6A
Other languages
English (en)
Other versions
CN107191483A (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.)
Yangzhong Guanjie Technology Innovation Co ltd
Zhenjiang Top Management Consulting Co ltd
Original Assignee
Jiangsu 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 Jiangsu University filed Critical Jiangsu University
Priority to CN201710284370.6A priority Critical patent/CN107191483B/zh
Publication of CN107191483A publication Critical patent/CN107191483A/zh
Application granted granted Critical
Publication of CN107191483B publication Critical patent/CN107191483B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16CSHAFTS; FLEXIBLE SHAFTS; ELEMENTS OR CRANKSHAFT MECHANISMS; ROTARY BODIES OTHER THAN GEARING ELEMENTS; BEARINGS
    • F16C32/00Bearings not otherwise provided for
    • F16C32/04Bearings not otherwise provided for using magnetic or electric supporting means
    • F16C32/0406Magnetic bearings
    • F16C32/044Active magnetic bearings
    • F16C32/0474Active magnetic bearings for rotary movement
    • F16C32/0485Active magnetic bearings for rotary movement with active support of three degrees of freedom
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16CSHAFTS; FLEXIBLE SHAFTS; ELEMENTS OR CRANKSHAFT MECHANISMS; ROTARY BODIES OTHER THAN GEARING ELEMENTS; BEARINGS
    • F16C32/00Bearings not otherwise provided for
    • F16C32/04Bearings not otherwise provided for using magnetic or electric supporting means
    • F16C32/0406Magnetic bearings
    • F16C32/044Active magnetic bearings
    • F16C32/0459Details of the magnetic circuit
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16CSHAFTS; FLEXIBLE SHAFTS; ELEMENTS OR CRANKSHAFT MECHANISMS; ROTARY BODIES OTHER THAN GEARING ELEMENTS; BEARINGS
    • F16C2202/00Solid materials defined by their properties
    • F16C2202/30Electric properties; Magnetic properties

Landscapes

  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Magnetic Bearings And Hydrostatic Bearings (AREA)

Abstract

本发明公开一种三自由度三极混合磁轴承的设计方法,根据磁轴承在实际应用场合确定出轴向承载力FZmax和径向承载力Frmax,比较轴向承载力FZmax和径向承载力Frmax;若3FZmax>8Frmax,则先确定出轴向饱和磁感应强度BZS和轴向偏置磁感应强度,再确定出径向偏置磁感应强度;若3FZmax<8Frmax,则先确定出径向饱和磁感应强度和径向偏置磁感应强度,然后确定轴向偏置磁感应强度;最后计算出径向磁极面积、轴向磁极面积和轴向饱和磁感应强度;径向和轴向承载力能独立选取,相互之间没有影响的,使得磁轴承的设计能够完全符合实际工况的要求,避免了控制线圈、空间体积的浪费。

Description

一种三自由度三极混合磁轴承的设计方法
技术领域
本发明涉及混合磁悬浮领域,具体是三自由度混合磁悬浮轴承的设计方法,用于径向和轴向共用偏置磁通结构的三自由度磁轴承。
背景技术
三自由度三极混合磁轴承的结构如图1和图2所示,三自由度三极混合磁轴承由转子1、轴向定子2、径向定子3、径向控制线圈4、轴向控制线圈5和永磁体6构成,最中间是转子1,径向定子3同轴套在转子1外部,轴向定子2位于转子1的轴向两端并且轴向定子2内部容纳了径向定子3。径向定子3有A相磁极31、B相磁极32、C相磁极33三个磁极,A、B、C相磁极31、32、33上分别缠绕了对应的A、B、C相控制线圈41、42、43。A、B、C相磁极31、32、33与转子1之间分别形成A相径向气隙81、B相径向气隙82和C相径向气隙83三个气隙。轴向定子2的两端部与转子1的两端部之间形成左侧轴向气隙71和右侧轴向气隙72。
对于磁轴承来说,由饱和磁感应强度和磁极面积共同决定磁轴承的承载力。在传统的磁轴承设计方法中,饱和磁感应强度的设定取决于铁磁材料的饱和值,偏置磁感应强度为饱和磁感应强度的一半,因此磁轴承的磁极面积就能由承载力计算得到。但对于三自由度三极混合磁轴承来说,其轴向偏置磁通和径向偏置磁通都由永磁体提供,因此相互之间存在一定的关系,导致了径向和轴向磁极面积之间有着比例关系,因此轴向承载力和径向承载力之间有一个固定的比例关系。然而,三自由度三极混合磁轴承在实际的不同工况下,其径向承载力和轴向承载力之间的关系并不能满足这一固定比例,因此,会造成磁轴承体积的浪费、功耗的增加和成本的提高,
发明内容
本发明的目的为了解决传统磁轴承设计方法用于三自由度三极混合磁轴承时径向与轴向承载力是固定比例关系因而存在不符合实际工况的问题,提出一种三自由度三极混合磁轴承的设计方法,径向、轴向承载力可以独立选取,相互之间没有影响。
本发明一种三自由度三极混合磁轴承的设计方法采用的技术方案为:根据磁轴承在实际应用场合确定出轴向承载力FZmax和径向承载力Frmax,比较轴向承载力FZmax和径向承载力Frmax;若3FZmax>8Frmax,则先根据铁磁材料特性确定出轴向饱和磁感应强度BZS和轴向偏置磁感应强度BZ0,再根据条件0.4<Br0<0.8确定出径向偏置磁感应强度Br0;若3FZmax<8Frmax,则先根据铁磁材料特性确定出径向饱和磁感应强度BrS和径向偏置磁感应强度Br0,然后根据条件0.4<BZ0<0.8确定轴向偏置磁感应强度BZ0;最后计算出径向磁极面积Sr、轴向磁极面积SZ和轴向饱和磁感应强度BZS
与现有技术比,本发明的优点在于:由于将传统设计方法中导致轴向承载力和径向承载力之间固定比例关系的约束条件取消,径向和轴向承载力能独立选取,相互之间没有影响的,使得三自由度三极混合磁轴承的设计能够完全符合实际工况的要求,避免了控制线圈、空间体积的浪费。
附图说明
图1是三自由度三极混合磁轴承的主视图;
图2是图1的右视图;
图3是三自由度三极混合磁轴承的等效磁路图;
图4是本发明一种三自由度三极混合磁轴承的设计方法的流程图。
图中:1.转子;2.轴向定子;3.径向定子;4.径向控制线圈;5.轴向控制线圈;6.永磁体;31.A相磁极;32.B相磁极;33.C相磁极;41.A相控制线圈;42.B相控制线圈;43.C相控制线圈;71.左侧轴向气隙;72.右侧轴向气隙;81.A相气隙;82.B相气隙;83.C相气隙。
具体实施方式
如图3所示,根据等效磁路法,忽略图1-2中三自由度三极混合磁轴承的转子1、轴向定子2、径向定子3的磁阻,忽略漏磁、涡流效应,将永磁体6、径向控制线圈4、轴向控制线圈5等效为磁动势,将左侧轴向气隙71、右侧轴向气隙72、A、B、C相气隙81、82、83等效为磁导,得到图3所示的三自由度三极混合磁轴承的等效磁路图,建立磁轴承的数学模型。图3中,GA、GB、GC分别为A、B、C相径向气隙81、82、83的磁导,GZ1、GZ2分别为左侧轴向气隙71和右侧轴向气隙72的磁导,Nr为A、B、C相磁极31、32、33上缠绕的控制线圈匝数,NZ为轴向控制线圈5的总匝数,iA为A相径向控制线圈41中通入的电流,iB为B相径向控制线圈42中通入的电流,iC为C相控制线圈43中通入的电流,iZ为轴向控制线圈5中通入的电流,ΦA、ΦB、ΦC分别为A、B、C相气隙81、82、83中流过的磁通,ΦZ1、ΦZ2分别为左侧轴向气隙71和右侧轴向气隙72中流过的磁通,Fm为永磁体6产生的对外磁动势,Φm为永磁体6产生的总磁通。A、B、C相中的每相径向气隙等效的磁导GA、GB、GC分别与对应相的径向控制线圈41、42、43等效的磁动势相串接后再并联,形成径向并联磁路。左右侧轴向气隙71、72等效的磁导GZ1、GZ2分别与对应侧的轴向控制线圈5等效的磁动势相串接后再并联,形成轴向并联磁路。径向并联磁路和轴向并联磁路串接后与永磁体6产生的磁动势Fm相接形成磁轴承的等效磁路。根据等效磁路能推导出径向、轴向承载力与径向磁极面积、轴向磁极面积和轴向饱和磁感应强度之间的关系式。具体如下:
根据图3中的等效磁路,推导出各气隙中磁通的表达式:
设计承载力时以转子1位于中心位置为前提,因此,轴向两个气隙71、72的磁导相等:GZ1=GZ2=GZ0=μ0SZZ,μ0为真空磁导率,SZ为轴向定子2与转子1所构成轴向气隙部分的截面积,δZ为转子1位于中心位置时左侧轴向气隙71和右侧轴向气隙72的长度。径向三个气隙81、82、83的磁导相等:GA=GB=GC=Gr0=μ0Srr,Sr为A、B、C相径向磁极31、32、33的截面积,δr为转子1位于中心位置时A、B、C相气隙81、82、83的长度,因此式(1)可以简化为:
式(2)中,ΦZ0为永磁体6在轴向气隙中产生的磁通,Φr0为永磁体6在径向气隙中产生的磁通,kZ为轴向电流的刚度,kr为径向控制电流的刚度,具体表达式为:
再结合式(2)中左侧轴向气隙71的磁通ΦZ1可以计算出:
ΦZ1=BZ0SZ+kZNZiZmax=BZSSZ (3)
式(3)中,iZmax为轴向控制线圈5中能通入的最大电流,BZ0为永磁体6在轴向气隙中产生的饱和磁感应强度,由于ΦZ0为永磁体6在轴向气隙中产生的磁通,因此有BZ0SZ=ΦZ0,BZS为轴向气隙中能达到的饱和磁感应强度,因轴向控制线圈5中通入轴向最大控制电流iZmax,所以有iZ=iZmax,左侧轴向气隙71中的磁感应强度为轴向饱和磁感应强度BZS
由式(3)可以得到:
kZNZiZmax=BZSSZ-BZ0SZ (4)
因此右侧轴向气隙72中的磁通ΦZ2可以通过式(2)和式(4)计算得到:
ΦZ2=BZ0SZ-kZNZiZmax=2BZ0SZ-BZSSZ (5)
令FZmax为轴向承载力,则轴向承载力FZmax与轴向饱和磁感应强度BZS、轴向偏置磁感应强度BZ0和轴向磁极面积SZ之间关系的表达式可以将式(3)和式(5)代入下式得到:
令irmax为径向控制线圈5中能通入的最大控制电流,Br0为永磁体6在径向气隙中产生的磁感应强度,由于Φr0为永磁体6在径向气隙中产生的磁通,因此有Br0Sr=Φr0,令BrS为径向气隙中能达到的饱和磁感应强度。以A相磁极31为例,令A相控制线圈41中通入径向最大控制电流irmax,有iA=irmax,A相气隙81中的磁感应强度达到径向饱和磁感应强度BSr,A相气隙81中的磁通ΦA为BrSSr,则可以表示为:
ΦA=Br0Sr+krNrirmax=BrSSr (7)
则由式(7)可以得到径向最大控制电流irmax与径向饱和磁感应强度BrS和径向偏置磁感应强度Br0之间的关系为:
krNrirmax=BrSSr-Br0Sr (8)
令A相控制线圈41中通入负径向最大控制电流-irmax,即iA=-irmax,A相气隙81中的磁通ΦA可以根据式(2)和式(8)计算得到:
ΦA=Br0Sr-krNrirmax=2Br0Sr-BrSSr (9)
令A相控制线圈41中通入负径向最大控制电流,即iA=-irmax,令B相控制线圈42、C相控制线圈43中通入径向最大控制电流的一半0.5irmax,即iB=iC=0.5irmax,产生径向承载力Frmax,此时B相气隙82、C相气隙83中的磁通ΦB、ΦC可以由式(2)和式(8)计算得到:
ΦB=ΦC=Br0Sr+0.5krNrirmax=0.5Br0Sr+0.5BrSSr (10)
Frmax为径向承载力,则径向承载力Frmax与径向饱和磁感应强度BrS、径向偏置磁感应强度Br0和径向磁极面积Sr之间关系的表达式可以将式(7)、式(9)和式(10)代入下式得到:
由于A、B、C三相径向气隙81、82、83与左侧轴向气隙71、右侧轴向气隙72共用一个偏置磁通,因此,A、B、C三相气隙81、82、83中的偏置磁通之和(3Br0Sr)等于左侧轴向气隙71、右侧轴向气隙72中的偏置磁通之和(2BZ0SZ),则有:
3Br0Sr=2BZ0SZ (12)
由式(6)、(11)和(12)可以看出,三个方程中有轴向承载力FZmax、径向承载力Frmax、轴向偏置磁感应强度BZ0、轴向饱和磁感应强度BZS、径向偏置磁感应强度Br0、径向饱和磁感应强度BrS和轴向磁极面积Sz、径向磁极面积Sr,总共8个变量,而有3个等式,说明有5个变量可以自由设定。如果根据铁磁材料的特性先确定轴向和径向的磁感应强度,即令BZS=BrS=0.8,BZ0=Br0=0.4,则根据式(6)、式(11)和式(12)可以计算出轴向承载力和径向承载力之间有一个固定的关系:3FZmax=8Frmax。而实际工况中轴向承载力和径向承载力的要求不一定满足这个关系,因此这样设计的磁轴承会造成体积浪费、功耗的增加和成本的提高。
本发明根据三自由度三极混合磁轴承在实际应用场合的要求,先确定轴向承载力FZmax和径向承载力Frmax的具体值,然后再需要确定轴向偏置磁感应强度BZ0、轴向饱和磁感应强度BZS、径向偏置磁感应强度Br0和径向饱和磁感应强度BrS中的三个变量,其具体方法为:比较轴向承载力FZmax和径向承载力Frmax的比例大小:若3FZmax>8Frmax,先根据铁磁材料的本质特性确定轴向饱和磁感应强度BZS和轴向偏置磁感应强度BZ0,然后再确定径向偏置磁感应强度Br0的值,并且需要满足条件0.4<Br0<0.8。若3FZmax<8Frmax,先根据铁磁材料的本质特性确定径向饱和磁感应强度BrS的值和径向偏置磁感应强度Br0的值,然后再确定轴向偏置磁感应强度BZ0的值,并且需要满足条件0.4<BZ0<0.8。
然后根据式(11)计算出径向磁极的面积Sr
再根据式(12)计算出轴向磁极的面积SZ
最后根据公式(6)计算出轴向饱和磁感应强度BZS的值BZS
基于以上步骤,就可以设计出三自由度三极混合磁轴承的基本参数:径向饱和磁感应强度BrS、径向偏置磁感应强度Br0、轴向饱和磁感应强度BZS、轴向偏置磁感应强度BZ0以及径向磁极面积Sr和轴向磁极面积SZ。并且其径向承载力Frmax与轴向承载力FZmax可以独立选取,相互之间没有影响。
以FZmax=100N,Frmax=200N为例,由于所选取的径向和轴向承载力满足3FZmax<8Frmax,因此可以先确定径向饱和磁感应强度BrS=0.8和径向偏置磁感应强度Br0=0.4,再选取轴向偏置磁感应强度,这里取BZ0=0.5,计算出Sr=1396mm2、SZ=1675.2mm2以及BZS≈0.575。以这些基本参数设计出的三自由度三极混合磁轴承完全符合实际工况的要求。

Claims (5)

1.一种三自由度三极混合磁轴承的设计方法,根据磁轴承在实际应用场合确定出轴向承载力FZmax和径向承载力Frmax,其特征是:比较轴向承载力FZmax和径向承载力Frmax;若3FZmax>8Frmax,则先根据铁磁材料特性确定出轴向饱和磁感应强度BZS和轴向偏置磁感应强度BZ0,再根据条件0.4<Br0<0.8确定出径向偏置磁感应强度Br0;若3FZmax<8Frmax,则先根据铁磁材料特性确定出径向饱和磁感应强度BrS和径向偏置磁感应强度Br0,然后根据条件0.4<BZ0<0.8确定轴向偏置磁感应强度BZ0;最后计算出径向磁极面积Sr、轴向磁极面积SZ和轴向饱和磁感应强度BZS
2.根据权利要求1所述的一种三自由度三极混合磁轴承的设计方法,其特征是:由式计算出径向磁极面积Sr
3.根据权利要求2所述的一种三自由度三极混合磁轴承的设计方法,其特征是:由式计算出轴向磁极面积SZ
4.根据权利要求3所述的一种三自由度三极混合磁轴承的设计方法,其特征是:由式计算出轴向饱和磁感应强度BZS
5.根据权利要求1所述的一种三自由度三极混合磁轴承的设计方法,其特征是:将磁轴承的每相径向气隙等效的磁导分别与对应相的径向控制线圈等效的磁动势相串接后再并联,形成径向并联磁路;将左右侧轴向气隙等效的磁导分别与对应侧的轴向控制线圈等效的磁动势相串接后再并联,形成轴向并联磁路;径向并联磁路和轴向并联磁路串接后与永磁体产生的磁动势相接形成磁轴承的等效磁路;根据等效磁路推导出径向承载力Frmax、轴向承载力FZmax与径向磁极面积Sr、轴向磁极面积SZ和轴向饱和磁感应强度BZS之间的关系式。
CN201710284370.6A 2017-04-27 2017-04-27 一种三自由度三极混合磁轴承的设计方法 Active CN107191483B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710284370.6A CN107191483B (zh) 2017-04-27 2017-04-27 一种三自由度三极混合磁轴承的设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710284370.6A CN107191483B (zh) 2017-04-27 2017-04-27 一种三自由度三极混合磁轴承的设计方法

Publications (2)

Publication Number Publication Date
CN107191483A CN107191483A (zh) 2017-09-22
CN107191483B true CN107191483B (zh) 2019-06-28

Family

ID=59872563

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710284370.6A Active CN107191483B (zh) 2017-04-27 2017-04-27 一种三自由度三极混合磁轴承的设计方法

Country Status (1)

Country Link
CN (1) CN107191483B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111089116B (zh) * 2020-01-17 2021-08-17 淮阴工学院 悬浮力对称六极混合磁轴承的设计方法
CN111173838B (zh) * 2020-01-17 2023-05-26 淮阴工学院 径向无耦合三自由度直流混合磁轴承
CN112332709B (zh) * 2020-09-24 2022-03-22 江苏大学 一种圆柱转子三自由度磁悬浮开关磁阻电机

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100432461C (zh) * 2005-05-18 2008-11-12 江苏大学 三自由度交直流径向-轴向混合磁轴承及其控制方法
CN105259785A (zh) * 2015-10-28 2016-01-20 江苏大学 三自由度混合磁轴承变饱和柔性变结构控制方法

Also Published As

Publication number Publication date
CN107191483A (zh) 2017-09-22

Similar Documents

Publication Publication Date Title
CN107191484B (zh) 一种径向六极的三自由度混合磁轴承的设计方法
US8354767B2 (en) Permanent magnet rotor with flux concentrating pole pieces
Bellara et al. Two-dimensional exact analytical solution of armature reaction field in slotted surface mounted PM radial flux synchronous machines
CN107191483B (zh) 一种三自由度三极混合磁轴承的设计方法
NO20013869L (no) Elektrisk multipolmotor/generator med aksial magnetisk fluks
WO2013078861A1 (zh) 磁场调制式横向磁通多相永磁电机
Zheng et al. Analysis and optimization of a novel tubular staggered-tooth transverse-flux PM linear machine
US10447127B2 (en) Electric motor comprising solenoid cores having coil slot
CN107289003B (zh) 同极式永磁偏置径向磁轴承
CN105990922B (zh) 转子及具有其的切向式永磁同步电机
CN106981937B (zh) 一种转子错位结构电机
CN106655553B (zh) 一种复合结构电机
CN204103628U (zh) 一种混合励磁双气隙爪极电机
US20190165628A1 (en) Interior Permanent Magnet Synchronous Machine
Yamazaki et al. Modeling of salient-pole synchronous machines assisted by permanent magnets
Li et al. A novel E-core hybrid-excited flux switching machine based on biased flux
CN207625413U (zh) 一种集成式永磁复合电机
CN109995162A (zh) 一种一体多用永磁同步电机转子
Al-Ani et al. Switched flux permanent magnet machine with segmanted magnets
Bastawade et al. No-load magnetic field analysis of double stator double Rotor radial flux permanent magnet generator for low power wind turbines
CN105449959A (zh) 同轴双转子超能电动机
Bernatt et al. Experimental Validation of Hybrid Excited Permanent Magnet Synchronous Generator
CN104009557B (zh) 一种产生轴向力的短磁路开关磁阻电机
Kim et al. The influence of stator pole shape and its arrangements on cogging torque for double-sided AFPM generator
Dirba et al. Optimization of the magnetic circuit of an axial inductor machine based on the calculation and analysis of magnetic field

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: 20211215

Address after: 212200 GANGLONG road science and technology innovation center, Yangzhong Development Zone, Zhenjiang City, Jiangsu Province

Patentee after: Yangzhong Guanjie Technology Innovation Co.,Ltd.

Address before: 212200 No. 239 Chunliu North Road, Sanmao street, Yangzhong City, Zhenjiang City, Jiangsu Province

Patentee before: Zhenjiang Top Management Consulting Co.,Ltd.

Effective date of registration: 20211215

Address after: 212200 No. 239 Chunliu North Road, Sanmao street, Yangzhong City, Zhenjiang City, Jiangsu Province

Patentee after: Zhenjiang Top Management Consulting Co.,Ltd.

Address before: Zhenjiang City, Jiangsu Province, 212013 Jingkou District Road No. 301

Patentee before: JIANGSU University