CN111665765B - 一种基于dsp的旋转变压器解码系统 - Google Patents
一种基于dsp的旋转变压器解码系统 Download PDFInfo
- Publication number
- CN111665765B CN111665765B CN202010605738.6A CN202010605738A CN111665765B CN 111665765 B CN111665765 B CN 111665765B CN 202010605738 A CN202010605738 A CN 202010605738A CN 111665765 B CN111665765 B CN 111665765B
- Authority
- CN
- China
- Prior art keywords
- entering
- quadrant
- dsp chip
- circuit
- rotary transformer
- 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
- 230000003750 conditioning effect Effects 0.000 claims abstract description 12
- 230000005284 excitation Effects 0.000 claims abstract description 8
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000000034 method Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 3
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 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/04—Programme control other than numerical control, i.e. in sequence controllers or logic controllers
- G05B19/042—Programme control other than numerical control, i.e. in sequence controllers or logic controllers using digital processors
-
- 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/04—Programme control other than numerical control, i.e. in sequence controllers or logic controllers
- G05B19/042—Programme control other than numerical control, i.e. in sequence controllers or logic controllers using digital processors
- G05B19/0423—Input/output
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01F—MAGNETS; INDUCTANCES; TRANSFORMERS; SELECTION OF MATERIALS FOR THEIR MAGNETIC PROPERTIES
- H01F38/00—Adaptations of transformers or inductances for specific applications or functions
- H01F38/18—Rotary transformers
-
- 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/20—Pc systems
- G05B2219/21—Pc I-O input output
- G05B2219/21137—Analog to digital conversion, ADC, DAC
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Power Engineering (AREA)
- Transmission And Conversion Of Sensor Element Output (AREA)
- Measurement Of Length, Angles, Or The Like Using Electric Or Magnetic Means (AREA)
Abstract
本发明公开了一种基于DSP的旋转变压器解码系统,属于伺服控制领域,该系统由DSP芯片、旋转变压器、激励信号源、信号调理电路和象限判断电路组成,系统采用同步采样和DMA技术实现信号高速采样,采用单向量S变换和ATO算法实现高精度旋转变压器解码算法,本发明为旋转变压器解码提供了一种低成本、高精度和实时的解决方案。
Description
技术领域
本发明属于伺服控制领域,涉及一种基于DSP的旋转变压器解码系统,尤其涉及一种采用TMS320F28335PGFA芯片的DMA技术、单向量S变换和角度跟踪算法的旋转变压器解码系统。
背景技术
旋转变压器作为一种位置传感器,其精度高、结构坚固、可靠耐用,广泛应用于轨道交通、机器人、飞机等工业领域,但旋转变压器输出信号是模拟信号,从中解算出转轴的位置信息,通常采用专用的解码芯片,这会增加系统成本,降低系统灵活性。
利用软件实现旋转变压器解码替代专用解码芯片越来越受到关注,这种方案通过将解码算法嵌入到伺服控制系统,在保证解码精度的同时,几乎不增加硬件成本,但需要对解码方法做精心设计;本发明以DSP芯片为核心,结合简单的外部电路,采用强抗干扰的单向量S变换和角度跟踪算法,实现旋转变压器的高精度解码;同时利用DSP芯片的DMA技术,减少DSP芯片CPU开销,保证解码系统的实时性。
发明内容
本发明的目的在于提供一种基于DSP的旋转变压器解码系统,用于伺服控制系统的旋转变压器角度解算,降低系统成本,提高解码精度和实时性。
本发明所述的一种基于DSP的旋转变压器解码系统,包括DSP芯片、旋转变压器、激励信号源、信号调理电路和象限判断电路,所述DSP芯片为TMS320F28335PGFA,DSP芯片的引脚GPIO20和GPIO21分别连接象限判断电路的输出信号qs和qc,DSP芯片的ADC通道ADCINA0和ADCINB0分别连接信号调理电路的输出信号和/>
本发明所述的象限判断电路的3路输入信号分别为激励信号源输出信号ue,旋转变压器输出正弦信号us和余弦信号uc,象限判断电路的2路输出信号分别为qs和qc,所述象限判断电路由3个过零比较器和2个异或电路组成,us,ue和uc分别输入过零比较器1、过零比较器2和过零比较器3,过零比较器1和过零比较器2的输出端分别连接异或电路1的2个输入端,异或电路1的输出信号qs连接DSP芯片引脚GPIO20,过零比较器2和过零比较器3的输出端分别连接异或电路2的2个输入端,异或电路2的输出信号qc连接DSP芯片引脚GPIO21。
S1:启动ADC同步采样通道ADCINA0和ADCINB0,进入步骤S2;
S2:判断采样是否完成,若是则进入步骤S3,否则等待;
S3:启动DMA进行数据传输,所述DMA是指DSP芯片的直接存储器访问技术,进入步骤S4;
S4:判断暂存标志fm是否等于0,若是则进入步骤S5,否则进入步骤S7;
S5:将采集的数据传输至数据缓冲区0,进入步骤S6;
S6:令fm=1,进入步骤S9;
S7:将采集的数据传输至数据缓冲区1,进入步骤S8;
S8:令fm=0,进入步骤S9;
S9:判断DMA传输是否完成,若是则进入步骤S10,否则等待;
S10:进入DMA中断,执行旋转变压器解码算法。
本发明所述的DSP芯片完成旋转变压器解码算法,求解出旋转变压器转动角度,所述旋转变压器解码算法包括如下步骤:
T2:采用如下公式计算角度差序列Δθ(n):
进入步骤T3;
T3:判断Δθ(n)是否小于et,所述et为预设的误差容限,其设置范围为0.01~0.0001,若是则进入步骤T5,否则进入步骤T4;
T4:令θ(n)=θ(n)+Δθ(n),返回步骤T2;
T5:判断θ(n)是否处于象限1,若是则进入步骤T11,否则进入步骤T6;
T6:判断θ(n)是否处于象限2,若是则进入步骤T7,否则进入步骤T8;
T7:令θ(n)=180-θ(n),进入步骤T1l;
T8:判断θ(n)是否处于象限3,若是则进入步骤T9,否则进入步骤T10;
T9:令θ(n)=180+θ(n),进入步骤T11;
T10:令θ(n)=360-θ(n),进入步骤T11;
T11:输出角度序列θ(n)。
本发明所述的DSP芯片检测象限判断电路的2路输出信号qs和qc,并利用这2路信号判断旋转变压器转轴位置所处象限,DSP芯片采用如下步骤实现象限判断:
P1:根据检测信号qs和qc的状态组合成二进制变量qx,进入步骤P2;
P2:判断qx是否等于00b,其中b表示二进制数值,若是则进入步骤P3,否则进入步骤P4;
P3:令象限值Qu=1,进入步骤P9;
P4:判断qx是否等于01b,若是则进入步骤P5,否则进入步骤P6;
P5:令象限值Qu=2,进入步骤P9;
P6:判断qx是否等于11b,若是则进入步骤P7,否则进入步骤P8;
P7:令象限值Qu=3,进入步骤P9;
P8:令象限值Qu=4,进入步骤P9;
P9:输出象限值Qu。
本发明的有益效果是,通过提供一种基于DSP的旋转变压器解码系统,只需在DSP芯片外部增加简单的信号调理和象限判断电路,结合DSP内部的解码算法,就可实现旋转变压器的实时高精度解码。
附图说明
图1为本发明的系统整体框图。
图2为本发明的象限判断电路原理图。
图3为本发明的信号采集流程图。
图4为本发明的旋转变压器解码算法流程图。
图5为本发明的象限判断方法流程图。
具体实施方式
以下结合附图对本发明的优选实施实例进行阐述,需要说明的是,优选实施实例是为了进一步说明本发明,而不是为了限制本发明的保护范围。
图1为本发明的系统整体框图,所述的一种基于DSP的旋转变压器解码系统,包括DSP芯片、旋转变压器、激励信号源、信号调理电路和象限判断电路,所述DSP芯片为TMS320F28335PGFA,DSP芯片的引脚GPIO20和GPIO21分别连接象限判断电路的输出信号qs和qc,DSP芯片的ADC通道ADCINA0和ADCINB0分别连接信号调理电路的输出信号和/>
图2为本发明的象限判断电路原理图,其3路输入信号分别为激励信号源输出信号ue,旋转变压器输出正弦信号us和余弦信号uc,象限判断电路的2路输出信号分别为qs和qc,所述象限判断电路由3个过零比较器和2个异或电路组成,us,ue和uc分别输入过零比较器1、过零比较器2和过零比较器3,过零比较器1和过零比较器2的输出端分别连接异或电路1的2个输入端,异或电路1的输出信号qs连接DSP芯片引脚GPIO20,过零比较器2和过零比较器3的输出端分别连接异或电路2的2个输入端,异或电路2的输出信号qc连接DSP芯片引脚GPIO21。
本发明所述的信号调理电路由电压偏置电路和比例运算电路组成,作用是将旋转变压器的双极性电压信号转换为0~3V之间的电压信号,以满足DSP芯片的ADC通道对输入信号的要求。
S1:启动ADC同步采样通道ADCINA0和ADCINB0,进入步骤S2;
S2:判断采样是否完成,若是则进入步骤S3,否则等待;
S3:启动DMA进行数据传输,所述DMA是指DSP芯片的直接存储器访问技术,进入步骤S4;
S4:判断暂存标志fm是否等于0,若是则进入步骤S5,否则进入步骤S7;
S5:将采集的数据传输至数据缓冲区0,进入步骤S6;
S6:令fm=1,进入步骤S9;
S7:将采集的数据传输至数据缓冲区1,进入步骤S8;
S8:令fm=0,进入步骤S9;
S9:判断DMA传输是否完成,若是则进入步骤S10,否则等待;
S10:进入DMA中断,执行旋转变压器解码算法。
图4为本发明的旋转变压器解码算法,所述旋转变压器解码算法包括如下步骤:
所述的单向量S变换采用快速算法,首先对输入信号做FFT(快速傅里叶变换)运算,然后对FFT结果针对激励频率点移位,最后对移位结果做反FFT运算。
T2:采用如下公式计算角度差序列Δθ(n):
进入步骤T3;
T3:判断Δθ(n)是否小于et,所述et为预设的误差容限,其设置范围为0.01~0.0001,若是则进入步骤T5,否则进入步骤T4;
T4:令θ(n)=θ(n)+Δθ(n),返回步骤T2;
T5:判断θ(n)是否处于象限1,若是则进入步骤T11,否则进入步骤T6;
T6:判断θ(n)是否处于象限2,若是则进入步骤T7,否则进入步骤T8;
T7:令θ(n)=180-θ(n),进入步骤T11;
T8:判断θ(n)是否处于象限3,若是则进入步骤T9,否则进入步骤T10;
T9:令θ(n)=180+θ(n),进入步骤T11;
T10:令θ(n)=360-θ(n),进入步骤T11;
T11:输出角度序列θ(n)。
图5为本发明的象限判断方法流程图,所述DSP芯片检测象限判断电路的2路输出信号qs和qc,并利用这2路信号判断旋转变压器转轴位置所处象限,DSP芯片采用如下步骤实现象限判断:
P1:根据检测信号qs和qc的状态组合成二进制变量qx,进入步骤P2;
P2:判断qx是否等于00b,其中b表示二进制数值,若是则进入步骤P3,否则进入步骤P4;
P3:令象限值Qu=1,进入步骤P9;
P4:判断qx是否等于01b,若是则进入步骤P5,否则进入步骤P6;
P5:令象限值Qu=2,进入步骤P9;
P6:判断qx是否等于11b,若是则进入步骤P7,否则进入步骤P8;
P7:令象限值Qu=3,进入步骤P9;
P8:令象限值Qu=4,进入步骤P9;
P9:输出象限值Qu。
以上所述的本发明的实施方式,并非成为本发明保护范围的限定,倘若对本发明实施方式进行各种变形或修改,但尚在本发明的精神和原则之内,均应包含在本发明的权利要求保护范围之内。
Claims (4)
1.一种基于DSP的旋转变压器解码系统,包括DSP芯片、旋转变压器、激励信号源、信号调理电路和象限判断电路,所述DSP芯片为TMS320F28335PGFA,DSP芯片的引脚GPIO20和GPIO21分别连接象限判断电路的输出信号qs和qc,DSP芯片的ADC通道ADCINA0和ADCINB0分别连接信号调理电路的输出信号和/>所述象限判断电路的3路输入信号分别为激励信号源输出信号ue,旋转变压器输出正弦信号us和余弦信号uc,象限判断电路的2路输出信号分别为qs和qc,其特征在于,所述象限判断电路由3个过零比较器和2个异或电路组成,us、ue和uc分别输入过零比较器1、过零比较器2和过零比较器3,过零比较器1和过零比较器2的输出端分别连接异或电路1的2个输入端,异或电路1的输出信号qs连接DSP芯片引脚GPIO20,过零比较器2和过零比较器3的输出端分别连接异或电路2的2个输入端,异或电路2的输出信号qc连接DSP芯片引脚GPIO21。
2.根据权利要求1所述的一种基于DSP的旋转变压器解码系统,其DSP芯片的ADC通道ADCINA0和ADCINB0分别连接信号调理电路的输出信号和/>其特征在于,DSP芯片采集信号/>和/>包括如下步骤:
S1:启动ADC同步采样通道ADCINA0和ADCINB0,进入步骤S2;
S2:判断采样是否完成,若是则进入步骤S3,否则等待;
S3:启动DMA进行数据传输,所述DMA是指DSP芯片的直接存储器访问技术,进入步骤S4;
S4:判断暂存标志fm是否等于0,若是则进入步骤S5,否则进入步骤S7;
S5:将采集的数据传输至数据缓冲区0,进入步骤S6;
S6:令fm=1,进入步骤S9;
S7:将采集的数据传输至数据缓冲区1,进入步骤S8;
S8:令fm=0,进入步骤S9;
S9:判断DMA传输是否完成,若是则进入步骤S10,否则等待;
S10:进入DMA中断,执行旋转变压器解码算法。
3.根据权利要求1所述的一种基于DSP的旋转变压器解码系统,其DSP芯片完成旋转变压器解码算法,求解出旋转变压器转动角度,其特征在于,所述旋转变压器解码算法包括如下步骤:
T2:采用如下公式计算角度差序列Δθ(n):
进入步骤T3;
T3:判断Δθ(n)是否小于et,所述et为预设的误差容限,其设置范围为0.01~0.0001,若是则进入步骤T5,否则进入步骤T4;
T4:令θ(n)=θ(n)+Δθ(n),返回步骤T2;
T5:判断θ(n)是否处于象限1,若是则进入步骤T11,否则进入步骤T6;
T6:判断θ(n)是否处于象限2,若是则进入步骤T7,否则进入步骤T8;
T7:令θ(n)=180-θ(n),进入步骤T11;
T8:判断θ(n)是否处于象限3,若是则进入步骤T9,否则进入步骤T10;
T9:令θ(n)=180+θ(n),进入步骤T11;
T10:令θ(n)=360-θ(n),进入步骤T11;
T11:输出角度序列θ(n)。
4.根据权利要求1所述的一种基于DSP的旋转变压器解码系统,其DSP芯片检测象限判断电路的2路输出信号qs和qc,并利用这2路信号判断旋转变压器转轴位置所处象限,其特征在于,DSP芯片采用如下步骤实现象限判断:
P1:根据检测信号qs和qc的状态组合成二进制变量qx,进入步骤P2;
P2:判断qx是否等于00b,其中b表示二进制数值,若是则进入步骤P3,否则进入步骤P4;
P3:令象限值Qu=1,进入步骤P9;
P4:判断qx是否等于01b,若是则进入步骤P5,否则进入步骤P6;
P5:令象限值Qu=2,进入步骤P9;
P6:判断qx是否等于11b,若是则进入步骤P7,否则进入步骤P8;
P7:令象限值Qu=3,进入步骤P9;
P8:令象限值Qu=4,进入步骤P9;
P9:输出象限值Qu。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010605738.6A CN111665765B (zh) | 2020-06-29 | 2020-06-29 | 一种基于dsp的旋转变压器解码系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010605738.6A CN111665765B (zh) | 2020-06-29 | 2020-06-29 | 一种基于dsp的旋转变压器解码系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111665765A CN111665765A (zh) | 2020-09-15 |
CN111665765B true CN111665765B (zh) | 2023-06-13 |
Family
ID=72390315
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010605738.6A Active CN111665765B (zh) | 2020-06-29 | 2020-06-29 | 一种基于dsp的旋转变压器解码系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111665765B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113358965A (zh) * | 2021-06-24 | 2021-09-07 | 玄天智能科技(常州)有限公司 | 一种旋转变压器的象限鉴相系统及鉴相方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN202083398U (zh) * | 2011-04-22 | 2011-12-21 | 长丰集团有限责任公司 | 一种旋转变压器信号调理电路 |
CN102403938A (zh) * | 2011-12-14 | 2012-04-04 | 兰州交通大学 | 基于单fpga的旋转变压器解码处理装置和方法 |
CN105703773A (zh) * | 2014-11-28 | 2016-06-22 | 中国科学院沈阳自动化研究所 | 一种旋转变压器解码方法 |
CN106374931A (zh) * | 2016-09-27 | 2017-02-01 | 湖南工业大学 | 一种采用单频率点s变换的旋转变压器信号解码方法 |
WO2017166793A1 (zh) * | 2016-04-01 | 2017-10-05 | 北京新能源汽车股份有限公司 | 旋转变压器解码芯片的故障信息处理方法和系统 |
CN209656071U (zh) * | 2019-04-19 | 2019-11-19 | 武汉蔚铭自控科技有限公司 | 一种数字输出的分装式旋转变压器 |
-
2020
- 2020-06-29 CN CN202010605738.6A patent/CN111665765B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN202083398U (zh) * | 2011-04-22 | 2011-12-21 | 长丰集团有限责任公司 | 一种旋转变压器信号调理电路 |
CN102403938A (zh) * | 2011-12-14 | 2012-04-04 | 兰州交通大学 | 基于单fpga的旋转变压器解码处理装置和方法 |
CN105703773A (zh) * | 2014-11-28 | 2016-06-22 | 中国科学院沈阳自动化研究所 | 一种旋转变压器解码方法 |
WO2017166793A1 (zh) * | 2016-04-01 | 2017-10-05 | 北京新能源汽车股份有限公司 | 旋转变压器解码芯片的故障信息处理方法和系统 |
CN106374931A (zh) * | 2016-09-27 | 2017-02-01 | 湖南工业大学 | 一种采用单频率点s变换的旋转变压器信号解码方法 |
CN209656071U (zh) * | 2019-04-19 | 2019-11-19 | 武汉蔚铭自控科技有限公司 | 一种数字输出的分装式旋转变压器 |
Non-Patent Citations (1)
Title |
---|
张伟鹏.一种高精度旋转变压器解码算法.科学技术与工程.2019,全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN111665765A (zh) | 2020-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101836085B (zh) | 使用定向跟踪无接触地感测旋转和角度位置的方法和设备 | |
US20100250184A1 (en) | Rotation angle detection apparatus | |
CN110044388B (zh) | 一种旋转体的非接触式绝对角度位置传感器的使用方法 | |
CN111665765B (zh) | 一种基于dsp的旋转变压器解码系统 | |
CN112945288B (zh) | 一种旋转变压器全角度测量装置及方法 | |
CN112985325B (zh) | 正余弦编码器的位置解码方法、装置及计算机可读介质 | |
CN103036572A (zh) | 一种用于双通道旋转变压器的解码器及其解码方法 | |
CN113148234B (zh) | 一种低成本高精度微小卫星反作用飞轮及其控制实现方法 | |
CN110243401B (zh) | 光编码器的光电池和光磁编码器的解码装置 | |
Sun et al. | Software-based resolver-to-digital converter using the PLL tracking algorithm | |
CN209877990U (zh) | 一种旋转体的非接触式绝对角度位置传感器 | |
CN111181559B (zh) | 一种旋变软解码方法、装置、设备和存储介质 | |
RU167428U1 (ru) | Преобразователь синусно-косинусного сигнала в код положения | |
CN109739122B (zh) | 一种半波激励快速锁相方法 | |
CN107241099B (zh) | 角度传感器信号处理电路及处理方法 | |
CN218822420U (zh) | 一种编码器、电机和自动化系统 | |
CN116961511A (zh) | 一种改进型msmu硬件电流环及控制方法 | |
CN110784223B (zh) | 单芯片数字型象限分割技术的轴角-数字转换电路与方法 | |
CN112556732A (zh) | 一种磁电角度传感器、编码器和电机 | |
CN109217743B (zh) | 一种绝对式磁角度编码器输出电压控制系统及方法 | |
Kang et al. | Research on decoding algorithm of split-type magnetic encoder | |
CN215528993U (zh) | 一种轴角转换器模块 | |
CN114257151B (zh) | 一种永磁同步电机电流采样在线校准方法 | |
Sun et al. | Implementation of SVPWM Algorithm in Sensorless Motor Based on Different Microcontrollers | |
CN115189607B (zh) | 一种基于dsp的旋转变压器解码方法 |
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: 20240716 Address after: Floor 1, No. 232 Xilin Road, Chenghua District, Chengdu City, Sichuan Province, 610000 Patentee after: Sichuan Qifu Cloud Enterprise Management Group Co.,Ltd. Country or region after: China Address before: 412007 No. 88 Taishan Road, Tianyuan District, Zhuzhou City, Hunan Province Patentee before: HUNAN University OF TECHNOLOGY Country or region before: China |