专利名称：Sr电机行至相电感不随角度变化区域的无位置检测方法
技术领域：
本发明属于SR电机转子位置检测技术领域，尤其是涉及一种SR电机行至相电感 不随角度变化区域的无位置检测方法。
背景技术：
开关磁阻电动机简称(Switched Reluctance Motor简称SRD)是20世纪80年代
迅速发展起来的一种新型电机，它融SR电机与现代电力电子技术控制技术、微计算机技术
为一体，具有结构简单、坚固，易于调速，控制灵活，可靠性高、容错性强等特点，兼有异步电
动机变频调速系统和直流电动机调速系统的优点。位置闭环控制是开关磁阻电机的基本特
征之一，但位置传感器的存在使结构简单的优点逊色不少，同时也增加了成本、降低了系统
可靠性，探索实用的无位置检测器方案，成为众多科研人员十分关注的课题。 国内外对间接转子位置检测技术的研究始于八十年代初，各国学者在近二十年的
研究工作中，提出了多种间接位置检测方案。各间接转子位置检测方法的基本思想或原理
是相同的，无论是对工作相还是非工作相的处理，都是通过分析绕组电压方程，得到绕组磁
链、电感或电动势等参数的解析形式，转子位置信息可以从这些参数中提取出来，但大都普
遍在电感不随角度变化的区域内不能很好地满足检测精度的问题。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种SR电机 行至相电感不随角度变化区域的无位置检测方法，其步骤简单、实现方便、实用性强且检测 精度高、误差小，能有效解决现有SR电机运行进入相电感不随角度变化区域时所存在的检 测精度低、误差大等缺陷和不足。 为解决上述技术问题，本发明采用的技术方案是一种SR电机行至相电感不随角 度变化区域的无位置检测方法，其特征在于该方法包括以下步骤 步骤一、确定被检测SR电机的综合性能参数，并根据所确定的综合性能参数选择 适宜被检测SR电机的常规无位置检测方法所述综合性能参数包括被检测SR电机的结构 参数和基本电参数以及经检测所确定的被检测SR电机的电机特性曲线，其中所述结构参 数包括被检测SR电机的定子及转子的极对数、转子极距角和定子及转子的极弧；所述基本 电参数包括被检测SR电机的额定电压、额定电流和额定转速，所述电机特性曲线包括被检 测SR电机的电感特性曲线即电磁特性曲线和静态转矩特性曲线； 步骤二、通过处理器建立被检测SR电机的实际数学模型首先将步骤一中所确定 的被检测SR电机的综合性能参数输入至处理器，并通过所述处理器建立被检测SR电机的 实际数学模型且相应进行存储；同时，在所述处理器中建立能实现步骤一中所选择无位置 检测方法的智能无位置检测模块； 步骤三、利用处理器对被检测SR电机进行分区域角度位置检测，其分区域角度位 置检测过程如下
301、划分区域利用处理器根据步骤二中所建立被检测SR电机的实际数学模型， 确定被检测SR电机运行过程中相电感不随转子角度变化区域并相应进行存储；所述相电 感不随转子角度变化区域为被检测SR电机的定子极与转子极重叠区域和定子极与转子槽 重叠区域； 302、利用处理器对被检测SR电机运行过程中的转子角度位置进行连续检测，其 检测过程包括以下几步 3021、所述处理器调用比较判断模块比较判断运行过程中任一时刻被检测SR电 机的导通相即检测相是否进入步骤301中所划定的相电感不随转子角度变化区域，并将比 较判断结果上传至所述处理器； 对于所述检测相而言，所述相电感不随转子角度变化区域为所述检测相的定子极 与转子极重叠区域即相电感最大区域(b， c)和定子极与转子槽重叠区域即相电感最小区 域(0， a); 3022、通过处理器对所述比较判断结果进行分析处理当判断得出此时被检测SR 电机的检测相未进入所述相电感不随转子角度变化区域时，所述处理器调用步骤二中所建 立的智能无位置检测模块对此时所述检测相的角度位置进行分析并相应估算得出此时被 检测SR电机的转子角度位置后输出；之后，返回步骤3021 ; 反之，当判断得出此时被检测SR电机的检测相进入所划定的相电感不随转子角 度变化的区域时，所述处理器调用差值比较模块对此时被检测SR电机检测相的导通角e c
与该SR电机的步进角e^p作差值比较当比较得出此时检测相的导通角e。不大于一个
步进角e step时，说明检测相和与其相邻的相邻相间不存在导通重叠区域，则说明此时被检 测sr电机的检测相进入匀速检测区域(b， e^p+a)，在所述匀速检测区域(b， estep+a)内，
被检测SR电机的检测相处于匀速运动状态且匀速运行速度为所述检测相运行至b点处时 的运行速度Vb，此时所述处理器调用所述智能无位置检测模块对检测相运行至b点处时的 角度位置进行检测并相应推算得出运行速度Vb，之后所述处理器根据所述检测相运行至b
点处的角度位置和运行速度vb相应推算得出匀速检测区域(b， e step+a)内被检测sr电机
的转子角度位置后输出；之后，返回步骤3021 ; 当比较得出此时检测相的导通角e。大于一个步进角e^p时，说明检测相和与其 相邻的相邻相间存在导通重叠区域，则在所述导通重叠区域内，所述处理器首先调用所述 比较判断模块比较判断此时所述相邻相即重叠导通相是否进入所述相电感不随转子角度
变化区域，并将比较判断结果上传至所述处理器当判断得出此时被检测SR电机的重叠导
通相未进入所述相电感不随转子角度变化区域时，所述处理器调用步骤二中所建立的智能 无位置检测模块对此时所述重叠导通相的角度位置进行分析并相应估算得出此时被检测
SR电机的转子角度位置后输出；反之，当判断得出此时被检测SR电机的重叠导通相进入所 述相电感不随转子角度变化区域时，说明此时被检测SR电机的重叠导通相和被检测的导 通相均在匀速检测区域内，可视被检测SR电机为在单相导通状态下进入匀速检测区域，匀 速运行速度为所述检测相运行至b点处时的运行速度Vb，此时所述处理器调用所述智能无 位置检测模块对检测相运行至b点处时的角度位置进行检测并相应推算得出运行速度Vb， 之后所述处理器根据所述检测相运行至b点处的角度位置和运行速度Vb相应推算得出匀 速检测区域(b， estro+a)内被检测SR电机的转子角度位置。
上述步骤一所述无位置检测方法为根据基于电机特性曲线所建立的基于自适应 模糊推理系统ANFIS的无位置检测方法。
上述步骤二中所述的处理器为PC机。
本发明与现有技术相比具有以下优点 1、主要针对现有无位置检测方法在相电感不随转子位置变化的区域即定子和转 子重叠范围以及定子和转子槽重叠范围内检测时所存在的检测精度低、误差大等缺陷，不 仅原理清晰、容易理解，操作过程简单方便，实现方便。 2、相电感恒定区域内所采用的检测方法计算简单，速度快且实时性强。 3、将传无位置检测方法常用的智能算法与电机实际数学模型分析方法有机结合
在非线性区域采用非线性逼近能力强大的智能算法；在相电感不随转子变化区域内，根据
电机数学模型采用匀速检测法不同与以往的直接对位置检测的方法而是利用在该区域的
匀速特点以速度推算位置，不同区域内应用不同算法实现位置估算，能够让两种方法发挥
各自优势在不同区域的边界无缝连接。 4、设计合理，SR电机进入相电感不随转子变化区域无位置检测方法实现步骤简 单易行，首先，根据电机实际结构、数学模型建立ANFIS无位置检测方法，确立电感恒值区 域分别为最高电感区域和最低电感区域；之后，根据控制策略设定的开通角e。n、关断角
e。『导通角e。、步进角e^p，判断相邻两相是否有导通重叠区域在导通重叠区域内， 根据重叠导通相电感的状态来决定是选择从两相之中选择处于非线性区域的导通相应用
ANFIS无位置检测方法进行位置估算作为结果还是应用匀速方法检测角度位置；没有重叠 现象时，在其它非线性区域应用ANFIS方法，而在电感不随转子变化区域(b， estep+a)内应 用匀速方法检测角度位置，通过ANFIS方法在非线性边界区域b的位置推算出此点的速度， 作为可匀速检测区域的速度，根据匀速特点推算角度位置。具体而言当电机导通相进入相 电感不随转子变化区域之前，应用原有无位置检测方法进行位置估算；当电机导通相进入 相电感不随转子变化区域时，根据导通角大小判断有无重叠导通现象，结合具体控制角度 策略和该区域边界重新设定可匀速检测区域范围，从而实现以速度推算角度的检测。因为 控制SR电机要求得到较大有效转矩，并且尽量减少制动转矩，通常开通角9 。n在相电感最
小区域内，关断角在相电感最大到来之前，由此可知匀速检测区域范围在(b， a+e^p)，可
通过在最大电感区域边界b点的速度推算在此匀速检测区域的角度位置。 5、所述采用ANFIS理论的无位置检测方法具有很好非线性逼近能力在非恒值电 感区检测精度高，误差小，但在一个导通周期的相电感不随角度变化的区检测精度底，误差 大。 综上所述，本发明步骤简单、实现方便、实用性强且检测精度高、误差小，能有效解 决现有SR电机运行进入相电感不随角度变化区域时所存在的检测精度低、误差大等缺陷 和不足。 下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
图1为本发明的方法流程图。
图2为本发明一个转矩角范围内的相电感变化示意图。
具体实施例方式
如图1所示的一种SR电机行至相电感不随角度变化区域的无位置检测方法，包括 以下步骤 步骤一、确定被检测SR电机的综合性能参数，并根据所确定的综合性能参数选择 适宜被检测SR电机的常规无位置检测方法所述综合性能参数包括被检测SR电机的结构 参数和基本电参数以及经检测所确定的被检测SR电机的电机特性曲线，其中所述结构参 数包括被检测SR电机的定子及转子的极对数、转子极距角和定子及转子的极弧；所述基本 电参数包括被检测SR电机的额定电压、额定电流和额定转速，所述电机特性曲线包括被检 测SR电机的电感特性曲线即电磁特性曲线和静态转矩特性曲线。 本实施例中，被检测SR电机为一台8/6相SR电机。所述无位置检测方法为根据
基于电机特性曲线所建立的基于自适应模糊推理系统ANFIS的无位置检测方法。 步骤二、通过处理器建立被检测SR电机的实际数学模型首先将步骤一中所确定
的被检测SR电机的综合性能参数输入至处理器，并通过所述处理器建立被检测SR电机的
实际数学模型且相应进行存储；同时，在所述处理器中建立能实现步骤一中所选择无位置
检测方法的智能无位置检测模块。 本实施例中，所述的处理器为PC机。 步骤三、利用处理器对被检测SR电机进行分区域角度位置检测，其分区域角度位 置检测过程如下 301、划分区域利用处理器根据步骤二中所建立被检测SR电机的实际数学模型， 确定被检测SR电机运行过程中相电感不随转子角度变化区域并相应进行存储；所述相电 感不随转子角度变化区域为被检测SR电机的定子极与转子极的重叠区域和定子极与转子 槽的重叠区域。 302、利用处理器对被检测SR电机运行过程中的转子角度位置进行连续检测，其 检测过程包括以下几步 3021、所述处理器调用比较判断模块比较判断运行过程中任一时刻被检测SR电 机的导通相即检测相是否进入步骤301中所划定的相电感不随转子角度变化区域，并将比 较判断结果上传至所述处理器。 结合图2，对于所述检测相而言，所述相电感不随转子角度变化区域为所述检测相 的定子极与转子极重叠区域即相电感最大区域(b， c)和定子极与转子槽的重叠区域即相 电感最小区域(O，a)。 所述相电感最大区域(b，c)和相电感最小区域(O，a)内不产生有转子位置改变引 起绕组中磁链变化而感应的电动势即运动电动势，而不产生运动电动势将导致无电磁转矩 产生。 本实施例中，所述8/6相SR电机步进角e^p为3(T ，极距角、为90° ，测量 分析其电机电感特性曲线和静态转矩曲线，并通过基于自适应模糊推理系统ANFIS的无位 置检测方法实现无位置检测，该无位置检测方法在相电感随转子角度变化区域即非线性区 域具有良好精度的逼近能力，误差范围为±0.2° 。所述8/6相SR电机的一个电感周期内 相电感不随角度变化的区域，包括最大相电感区域(即定子凸极(或定子极)与转子凸极(或转子极)的重合区域)(b， c)为(33° ，57° )和最小相电感区域(即定子凸极(或定 子极)和转子槽重合区域)(0， a)为(0° ，7° )，在这些区域内无运动电动势产生，故不产 生电磁转矩，均可视为匀速运动区域。 3022、通过处理器对所述比较判断结果进行分析处理当判断得出此时被检测SR 电机的检测相未进入所述相电感不随转子角度变化区域时，所述处理器调用步骤二中所建 立的智能无位置检测模块对此时所述检测相的角度位置进行分析并相应估算得出此时被 检测SR电机的转子角度位置后输出；之后，返回步骤3021。也就是说，在检测相未进入所 述相电感不随转子角度变化区域时，采用原有的无位置检测方法进行角度位置检测，实际 操作过程中，可采用导通相检测法、非导通相检测法、基于智能控制的无位置检测法等常规 检测方法进行检测。 反之，当判断得出此时被检测SR电机的检测相进入所划定的相电感不随转子角 度变化的区域时，所述处理器调用差值比较模块对此时被检测SR电机检测相的导通角e c
与该SR电机的步进角e^p作差值比较当比较得出此时检测相的导通角e。不大于一个
步进角e step时，说明检测相和与其相邻的相邻相同不存在导通重叠区域，则说明此时被检 测sr电机的检测相进入匀速检测区域(b， e^p+a)，在所述匀速检测区域(b， estep+a)内，
被检测SR电机的检测相处于匀速运动状态且匀速运行速度为所述检测相运行至b点处时 的运行速度Vb，此时所述处理器调用所述智能无位置检测模块对检测相运行至b点处时的 角度位置进行检测并相应推算得出运行速度Vb，之后所述处理器根据所述检测相运行至b
点处的角度位置和运行速度vb相应推算得出匀速检测区域(b， e step+a)内被检测sr电机
的转子角度位置后输出；之后，返回步骤3021 ;当比较得出此时检测相的导通角e。大于一个步进角e^p时，说明检测相和与其
相邻的相邻相间存在导通重叠区域，则在所述导通重叠区域内，所述处理器首先调用所述 比较判断模块比较判断此时所述相邻相即重叠导通相是否进入所述相电感不随转子角度
变化区域，并将比较判断结果上传至所述处理器当判断得出此时被检测SR电机的重叠导
通相未进入所述相电感不随转子角度变化区域时，所述处理器调用步骤二中所建立的智能 无位置检测模块对此时所述重叠导通相的角度位置进行分析并相应估算得出此时被检测
SR电机的转子角度位置后输出；反之，当判断得出此时被检测SR电机的重叠导通相进入所 述相电感不随转子角度变化区域时，说明此时被检测SR电机的重叠导通相和被检测的导 通相均在匀速检测区域内，可视SR电机为在单相导通状态下进入匀速检测区域，匀速运行 速度为所述检测相(非重叠导通相)运行至b点处时的运行速度Vb，此时所述处理器调用 所述智能无位置检测模块对检测相运行至b点处时的角度位置进行检测并相应推算得出 运行速度Vb，之后所述处理器根据所述检测相运行至b点处的角度位置和运行速度Vb相应
推算得出匀速检测区域(b， estep+a)内被检测SR电机的转子角度位置。 本实施例中，所述8/6相SR电机的检测相运行在相电感随角度变化区域(7° ， 33° )时，采用基于自适应模糊推理系统ANFIS的无位置检测方法估算其转子的角度位置。
所述8/6相SR电机的检测相运行在相电感不随转子角度变化区域(0° ，7° )和 (33° ，57° )时，根据具体角度控制策略，判断是否出现重叠导通现象。具体判断时，所述 处理器调用差值比较模块对此时被检测SR电机检测相的导通角9 。与该SR电机的步进角
estep = 30°作差值比较当比较得出此时检测相的导通角ee不大于30。时，说明检测相和与其相邻的相邻相间即所述重叠导通相不存在导通重叠区域；反之，当导通角e 。大于 30°时，说明检测相和重叠导通相间存在导通重叠区域。 具体而言，当检测相和重叠导通相间不存在导通重叠区域时，则说明此时被检测 8/6相SR电机的检测相进入匀速检测区域(33° ，7° +30° )即(33° ，37° )，在所述匀速 检测区域(33° ，37° )内，所述8/6相SR电机的检测相处于匀速运动状态且匀速运行速度 为所述检测相运行至b点处时的运行速度Vb，此时所述处理器调用所述智能无位置检测模 块对检测相运行至b点处时的角度位置进行检测并相应推算得出运行速度Vb，之后所述处 理器根据所述检测相运行至b点处的角度位置和运行速度Vb相应推算得出匀速检测区域 (33° ，37° )内被检测8/6相SR电机的转子角度位置后输出。 当检测相和重叠导通相间存在导通重叠区域时，则在所述导通重叠区域内，所述 处理器首先调用所述比较判断模块比较判断此时8/6相SR电机的重叠导通相是否进入不 随相电感变化区域(O。 ，7° )和(33° ，57° )，并将比较判断结果上传至所述处理器当 判断得出此时被检测8/6相SR电机的重叠导通相未进入所述相电感不随转子角度变化区 域时，所述处理器调用步骤二中所建立的智能无位置检测模块对此时所述重叠导通相的角 度位置进行分析并相应估算得出此时被检测SR电机的转子角度位置后输出。
反之，当判断得出此时被检测8/6相SR电机的重叠导通相进入所述相电感不随转 子角度变化区域时，说明此时被检测SR电机的重叠导通相和导通相均在匀速检测区域内， 可视8/6相SR电机为在单相导通状态下进入匀速检测区域，且匀速运行速度为所述检测相 (非重叠导通相)运行至b点处时的运行速度Vb，此时所述处理器调用所述智能无位置检 测模块对检测相运行至b点处时的角度位置进行检测并相应推算得出运行速度Vb，之后所 述处理器根据所述检测相运行至b点处的角度位置和运行速度Vb相应推算得出匀速检测 区域(b， 9 step+a)内被检测SR电机的转子角度位置。 综上所述，在SR电机运行时，当SR电机行至相电感随角度变化区域时，采用基于
自适应模糊推理系统ANFIS的无位置检测方法进行角度位置检测；当SR电机行至相电感不
随角度变化区域时，采用匀速检测方法并划定匀速检测区域，将相电感随角度位置变化区
域到相电感不随角度位置变化区域的相交点作为基准点，并将此点的速度作为相电感不随
角度变化区域的速度，根据基准点的速度和基准点的位置估算匀速检测区域的角度位置，
匀速检测法在相电感不随角度变化区域内检测精度高误差小能够达到检测要求。 以上所述，仅是本发明的较佳实施例，并非对本发明作任何限制，凡是根据本发明
技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化，均仍属于本发明技
术方案的保护范围内。
权利要求
一种SR电机行至相电感不随角度变化区域的无位置检测方法，其特征在于该方法包括以下步骤步骤一、确定被检测SR电机的综合性能参数，并根据所确定的综合性能参数选择适宜被检测SR电机的常规无位置检测方法所述综合性能参数包括被检测SR电机的结构参数和基本电参数以及经检测所确定的被检测SR电机的电机特性曲线，其中所述结构参数包括被检测SR电机的定子及转子的极对数、转子极距角和定子及转子的极弧；所述基本电参数包括被检测SR电机的额定电压、额定电流和额定转速，所述电机特性曲线包括被检测SR电机的电感特性曲线即电磁特性曲线和静态转矩特性曲线；步骤二、通过处理器建立被检测SR电机的实际数学模型首先将步骤一中所确定的被检测SR电机的综合性能参数输入至处理器，并通过所述处理器建立被检测SR电机的实际数学模型且相应进行存储；同时，在所述处理器中建立能实现步骤一中所选择无位置检测方法的智能无位置检测模块；步骤三、利用处理器对被检测SR电机进行分区域角度位置检测，其分区域角度位置检测过程如下301、划分区域利用处理器根据步骤二中所建立被检测SR电机的实际数学模型，确定被检测SR电机运行过程中相电感不随转子角度变化区域并相应进行存储；所述相电感不随转子角度变化区域为被检测SR电机的定子极与转子极重叠区域和定子极与转子槽重叠区域；302、利用处理器对被检测SR电机运行过程中的转子角度位置进行连续检测，其检测过程包括以下几步3021、所述处理器调用比较判断模块比较判断运行过程中任一时刻被检测SR电机的导通相即检测相是否进入步骤301中所划定的相电感不随转子角度变化区域，并将比较判断结果上传至所述处理器；对于所述检测相而言，所述相电感不随转子角度变化区域为所述检测相的定子极与转子极重叠区域即相电感最大区域(b，c)和定子极与转子槽重叠区域即相电感最小区域(0，a)；3022、通过处理器对所述比较判断结果进行分析处理当判断得出此时被检测SR电机的检测相未进入所述相电感不随转子角度变化区域时，所述处理器调用步骤二中所建立的智能无位置检测模块对此时所述检测相的角度位置进行分析并相应估算得出此时被检测SR电机的转子角度位置后输出；之后，返回步骤3021；反之，当判断得出此时被检测SR电机的检测相进入所划定的相电感不随转子角度变化的区域时，所述处理器调用差值比较模块对此时被检测SR电机检测相的导通角θc与该SR电机的步进角θstep作差值比较当比较得出此时检测相的导通角θc不大于一个步进角θstep时，说明检测相和与其相邻的相邻相间不存在导通重叠区域，则说明此时被检测SR电机的检测相进入匀速检测区域(b，θstep+a)，在所述匀速检测区域(b，θstep+a)内，被检测SR电机的检测相处于匀速运动状态且匀速运行速度为所述检测相运行至b点处时的运行速度Vb，此时所述处理器调用所述智能无位置检测模块对检测相运行至b点处时的角度位置进行检测并相应推算得出运行速度Vb，之后所述处理器根据所述检测相运行至b点处的角度位置和运行速度Vb相应推算得出匀速检测区域(b，θstep+a)内被检测SR电机的转子角度位置后输出；之后，返回步骤3021；当比较得出此时检测相的导通角θc大于一个步进角θstep时，说明检测相和与其相邻的相邻相间存在导通重叠区域，则在所述导通重叠区域内，所述处理器首先调用所述比较判断模块比较判断此时所述相邻相即重叠导通相是否进入所述相电感不随转子角度变化区域，并将比较判断结果上传至所述处理器当判断得出此时被检测SR电机的重叠导通相未进入所述相电感不随转子角度变化区域时，所述处理器调用步骤二中所建立的智能无位置检测模块对此时所述重叠导通相的角度位置进行分析并相应估算得出此时被检测SR电机的转子角度位置后输出；反之，当判断得出此时被检测SR电机的重叠导通相进入所述相电感不随转子角度变化区域时，说明此时被检测SR电机的重叠导通相和被检测的导通相均在匀速检测区域内，可视被检测SR电机为在单相导通状态下进入匀速检测区域，匀速运行速度为所述检测相运行至b点处时的运行速度Vb，此时所述处理器调用所述智能无位置检测模块对检测相运行至b点处时的角度位置进行检测并相应推算得出运行速度Vb，之后所述处理器根据所述检测相运行至b点处的角度位置和运行速度Vb相应推算得出匀速检测区域(b，θstep+a)内被检测SR电机的转子角度位置。
2. 按照权利要求1所述的SR电机行至相电感不随角度变化区域的无位置检测方法，其 特征在于步骤一所述无位置检测方法为根据基于电机特性曲线所建立的基于自适应模糊 推理系统ANFIS的无位置检测方法。
3. 按照权利要求1或2所述的SR电机行至相电感不随角度变化区域的无位置检测方 法，其特征在于步骤二中所述的处理器为PC机。
全文摘要
本发明公开了一种SR电机行至相电感不随角度变化区域的无位置检测方法，包括步骤一、确定被检测SR电机的综合性能参数并选择适宜被检测SR电机的常规无位置检测方法；二、建立被检测SR电机的实际数学模型；三、利用处理器对被检测SR电机进行分区域角度位置检测首先划定相电感不随转子角度变化区域，当电机导通相进入该区域之前应用原有无位置检测方法进行位置估算；当电机导通相进入相电感不随转子变化区域时，重新设定匀速检测区域范围实现以速度推算角度的检测。本发明步骤简单、实现方便、实用性强且检测精度高、误差小，能有效解决现有SR电机运行进入相电感不随角度变化区域时所存在的检测精度低、误差大等缺陷和不足。
文档编号H02P6/16GK101741301SQ200910254670
公开日2010年6月16日 申请日期2009年12月31日 优先权日2009年12月31日
发明者张晶园, 柴钰 申请人:西安科技大学