CN105703680A - 基于脉冲注入的三相开关磁阻电机无位置传感器控制方法 - Google Patents

基于脉冲注入的三相开关磁阻电机无位置传感器控制方法 Download PDF

Info

Publication number
CN105703680A
CN105703680A CN201610143963.6A CN201610143963A CN105703680A CN 105703680 A CN105703680 A CN 105703680A CN 201610143963 A CN201610143963 A CN 201610143963A CN 105703680 A CN105703680 A CN 105703680A
Authority
CN
China
Prior art keywords
phase
inductance
rotor
theta
current
Prior art date
Application number
CN201610143963.6A
Other languages
English (en)
Inventor
张辑
陈丽安
魏荣宗
陈天翔
彭彦卿
郑雪钦
Original Assignee
厦门理工学院
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 厦门理工学院 filed Critical 厦门理工学院
Priority to CN201610143963.6A priority Critical patent/CN105703680A/zh
Publication of CN105703680A publication Critical patent/CN105703680A/zh

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P2203/00Indexing scheme relating to controlling arrangements characterised by the means for detecting the position of the rotor
    • H02P2203/01Motor rotor position determination based on the detected or calculated phase inductance, e.g. for a Switched Reluctance Motor

Abstract

本发明涉及开关磁阻电机无位置传感器控制技术领域,特别地涉及一种基于脉冲注入的三相开关磁阻电机无位置传感器控制方法。本发明提出了一种基于脉冲注入的三相开关磁阻电机无位置传感器控制方法,适用于三相开关磁阻电机的控制,该方法采用对开关磁阻电机非导通相注入脉冲,导通相进行斩波控制,通过电流传感器检测相电流,计算出电流斜率差值,再根据电流斜率差得到三相等效电感,合成电感矢量,最后从电感矢量中提取转子位置信息,进而得到转子转速。本发明结构简单易于实现,计算量大大降低,可靠性和精度高。

Description

基于脉冲注入的H相开关磁阻电机无位置传感器控制方法
技术领域
[0001] 本发明属于开关磁阻电机无位置传感器控制技术领域,具体地设及一种结构简单 易于实现,可靠性和精度高的基于脉冲注入的=相开关磁阻电机无位置传感器控制方法。
背景技术
[0002] 开关磁阻电机(Switched Reluctance Moto;r,SRM),又被称为"磁阻电机",其结构 简单,控制灵活,功率变换器简单,工作可靠,具有广阔的应用前景。对于开关磁阻电动机调 速系统(Switched Reluc1:ance Motor Drives,S畑),准确的转子位置是电机稳定运行的必 要条件。对于某些环境较为恶劣的工业应用场合,比如矿井采煤机,传统位置传感器的检测 精度和可靠性受到严重影响。近年来开发开关磁阻电机无位置传感器控制系统已成为研究 执占之一。
[0003] 现有无位置传感器控制方法主要有电流磁链法、观测器法和智能控制法。A Position Sensorless Control System of SRM over Wide Speed Range. The 7th International Conference on Power Elect;ronics:640~643采用电流磁链法,其基本思 想是根据电机凸极性效应,不同转子位置与不同的磁链-电流曲线相对应,通过电流查表得 到参考磁链,但该方法需要大量查表数据,可移植性差。Rotor position observer for sensorless control of SRM,XIX Symposium Electromagnetic Phenomena in Nonlinear Circuits ,Maribor ,Slovenia,28-3 Ojune 2006,pp. 161 ~162 中采用 Luenberger观测器通过开关磁阻电机运行中的电流和电压信号计算转子角度,该方法计算 量比较大,而且观测器线性模型对参数敏感,鲁棒性较差。ANN based sensorless rotor position estimation for the Switched Reluctance Motor[C].Nirma University International Engineering(MJiC0肥),2011:1-6中采用神经网络智能控制法,该方法理 论上可W逼近任何非线性有理函数,但是需要大量训练样本,离实际运用仍有较大差距
发明内容
[0004] 本发明的目的在于为解决上述问题而提供一种结构简单易于实现,计算量大大降 低,可靠性和精度高的基于脉冲注入的=相开关磁阻电机无位置传感器控制方法。
[0005] 为此,本发明公开了一种基于脉冲注入的=相开关磁阻电机无位置传感器控制方 法,对开关磁阻电机的非导通相注入脉冲,对导通相进行电流斩波控制,测量各相的相电 流,计算出电流斜率差值,通过电流斜率差值得出=相等效电感,合成电感矢量,通过电感 矢量估算转子位置角度0,采用如下公式进行估算
Figure CN105703680AD00031
Figure CN105703680AD00041
[000引其中,目诚电感旋转角度,Nr为转子极数,La、Lb、Lc分别为立相的各相电感。
[0009] 进一步的,根据转子位置角度,计算出转子转速。
[0010] 更进一步的,转子转速采用W下公式进行计算
Figure CN105703680AD00042
[001^ 其中,0n和0n+l为两个相邻采样周期内通过电感矢量计算得到的转子位置角度,ft 是系统采样频率,《为转子速度。
[0013] 进一步的,采用实验法获取最佳的脉冲宽度。
[0014] 进一步的,对电流斜率差值进行数字滤波后,再计算得出=相等效电感。
[0015] 更进一步的,所述数字滤波采用最大值滤波法。
[0016] 本发明的有益效果:
[0017] 本发明通过对开关磁阻电机的非导通相注入脉冲,导通相采用斩波控制,根据相 电流斜率差值得到等效电感,从电感矢量中提取转子位置信息,结构简单易于实现,计算量 大大降低,易于实现实时控制,可靠性和精度高,鲁棒性好。
附图说明
[0018] 图1为本发明实施例的开关磁阻电机S相电感波形图;
[0019] 图2为本发明实施例的脉冲宽度与电流响应峰值表图;
[0020] 图3(a)为本发明实施例的电机运行时对空闲相注入脉冲,导通相进行斩波控制时 候的相电流波形图;
[0021] 图3(b)为本发明实施例的将相电流采样后计算得到的电流斜率差值图;
[0022] 图3(c)为本发明实施例的计算得到等效电感值图;
[0023] 图4为本发明实施例的实测转子位置与根据电感矢量估计的转子位置角波形图;
[0024] 图5为本发明实施例的转速波形图;
[0025] 图6为本发明实施例的实验平台结构图;
[0026] 图7为本发明实施例的位置辨识中断程序流程图;
[0027] 图8为本发明实施例的转速给定12化/min时实验波形;
[00%]图9为本发明实施例的转速给定25化/min时实验波形;
[0029] 图10为本发明实施例的无位置传感器启动实验波形。
具体实施方式
[0030] 现结合附图和具体实施方式对本发明进一步说明。
[0031] 基于脉冲注入的=相开关磁阻电机无位置传感器控制方法,通过对开关磁阻电机 的非导通相注入脉冲,对导通相进行电流斩波控制,然后测量各相的相电流,计算出各相电 流斜率差值,再通过电流斜率差值得出=相等效电感,合成电感矢量,最后通过电感矢量估 算出转子位置角度0,进而计算出转子转速。本发明适用于=相开关磁阻电机的控制。下面 将W =相(分别命名为a相、b相和C相)12/8结构的开关磁阻电机为例来对本发明进一步详 细说明。
[0032]对样机的非导通相注入激励脉冲,导通相采用电流斩波控制,关于脉冲的注入和 电流斩波控制可W参考现有技术,此不再细说。采用单相激励电流斩波控制模式,驱动器功 率管导通时,开关磁阻电机处于励磁状态,电压方程可写为:
Figure CN105703680AD00051
[0034]当驱动器功率管关断时,电流通过续流二极管流过,电机处于退磁状态,此时电压 方程可写为:
Figure CN105703680AD00052
[0036] 式(1)和(2)中,Udc, Rj, ij分别为开关磁阻电机的相电压,相电阻和相电流山(0)为 相电感。dijMtU和dij/dtUff分别为相电流的上升率和下降率,j表示第j相,j = a、b、c。两 模式下相电流和反电势的瞬时值并不相同,但是i神日k(e)变化很小,dL(0)/d0基本上不变, 式(1)减去式(2)整理后可得:
Figure CN105703680AD00053
[0038] 综上分析,通过对样机的非导通相注入激励脉冲,导通相采用斩波控制,测量S相 电流(可W采用电流传感器进行测量),计算出各相电流斜率差,然后通过式(3)可W计算出 S相电感La、Lb和Lc。
[0039] 为了进一步准确地得到电感,必须考虑注入激励脉冲的持续时间A t,本实施例采 用实验法得到最佳的脉冲宽度。样机对于不同脉冲宽度激励,电流响应的峰值关系如图2所 示。样机磁链特性饱和电流约20A,实验结果表明选择脉冲宽度为15化S和20化S比较合适, 而且对应的脉冲宽度越大,电流响应峰值越大,采样的电流信息信噪比含量越高,精度越 高。因而在实际系统中,选择脉冲宽度200US,脉冲周期500US。
[0040] 设样机的定转子凸极对齐位置为0 = 22.5°处,定子凸极和转子凹槽对齐位置为0 =0°处,相电感的电周期为360°,则样机的S相电感空间波形如图1所示。由图1可知,相电 感与标准正弦函数仍有些差别,但可W等效为=角函数叠加一个直流偏量,定义电感旋转 角度如式(4)所示
Figure CN105703680AD00054
[0042]其中,Nr为转子极数,Nr = 8,根据SM电感解析模型,可W近似得到毎相电感的数学 表达式,如式(5)所示
Figure CN105703680AD00061
[0044]其中,山=化祉。)/2心=山-以)/2山为转子与定子磁极对其时的相电感心为转 子与定子磁极没对齐时的相电感,由于=相开关磁阻电机相邻相绕组在空间上互差120°电 角度,所W =相电感可W看成相位互差120°的空间向量,其合成的旋转的电感矢量L,如式 (6)所示
Figure CN105703680AD00062
[0046]其中,将合成的电感矢量L正交分解分解成La和Le可得
Figure CN105703680AD00063
[004引 WLa作为基准电感,假设a相电感La达到峰值时,合成电感矢量L和a相电感La电感 方向相同,当目0 = 0°时,La达到最大值,电角度目elec = 3T,因此,00和0ele的关系可W表示为:
[0049]目 0=目 eiec-JT......(8)
[(K)加]又因为:
Figure CN105703680AD00064
[0052]实际转子位置角度(机械角度)为
Figure CN105703680AD00065
[0054]而空间电感的旋转角度00为
Figure CN105703680AD00066
[0056] 其中,目〇e(-ji/2,3T/2)。
[0057] 综上所述,通过=相电感合成电感矢量,再通过式(10)和(11)可W计算得出转子 位置角度。
[0058] 转子速度可由下式估算得到
Figure CN105703680AD00071
[0060] 式中0n和0n+功两个相邻采样周期内通过电感矢量法计算得到的转子位置,ft是系 统采样频率,CO为转子速度。
[0061] 得到转子位置角度和转子速度信息后,就可W进行SRD控制。
[0062] 为验证本发明的可行性,在MATLAB/SIMULINK软件平台搭建仿真模型,对上述方案 进行验证。本实施例W18.5kWS相12/8极的SRM电机为模型,将Matlab中自带的S相SRM电 机进行改进,仿真的具体参数如下:功率W = 18.5Kw,额定电压U = 380V,额定转速n= 1000 r/ min,给定转速 n*= 200;r/min,给定转矩 T* = 30N. m。
[0063] 图3所示为采用本发明得到的电感信息仿真结果。其中图3(a)所示为电机运行时 对非导通相注入脉冲,导通相进行斩波控制时候的相电流波形,图3(b)所示为将相电流采 样后计算得到的电流斜率差值,进而经过滤波后由公式(3)计算得到等效电感值,如图3(c) 所示。由图3可知,仿真得到的等效电感波形与理论分析吻合,故而能够通过电感矢量辨识 出转子位置信息。
[0064] 图4和图5所示为给定转速20化/min,负载30N*m时,根据电感矢量得到的转子位置 与转速仿真结果。其中,图4所示为开关磁阻电机实测转子位置与根据电感矢量估计的转子 位置角波形。从图中可W看出,估计得到转子位置与实际转子位置相吻合,精度较高,且实 现了静止启动到稳定运行状态下的跟踪。图5所示为由转子位置信息根据公式(12)计算得 到的转速波形,调速系统的稳态转速达到了给定值,且稳态误差小,采用该方法得到的转子 位置和转速信息能够满足调速系统的要求,验证了脉冲法用于实现无位置传感器系统的有 效性和可行性。
[0065] 为验证本发明的可行性,搭建了 SRD调速系统实验平台。实验平台结构图如图6所 示,包括:功率变换器、微处理器控制板、驱动电路、电压电流检测、保护电路和SM本体W及 直流负载电机组成。
[0066] 功率变换器为S相不对称半桥结构,主开关器件为英飞凌的IGBTdSRM为S相12/8 结构,额定转速lOOOr/min,额定电压380V,主控微处理器是TI公司的电机控制领域专用32 位DSP控制忍片TMS2812。
[0067] 图7所示为中断流程图,首先对非导通相注入激励脉冲,利用电流传感器对电机相 电流进行采样,通过AD转化后,用当前电流采样值减去上一采样周期储存的电流采样值,计 算得到到各相电流斜率,然后判断电流斜率差的符号,当斜率为正时存为LPslope,当斜率 为负时存为LNslope。用数字滤波得到电流斜率差的包络线,数字滤波选用最大值滤波方 法,再计算出=相等效电感,最后通过公式(7)和(8)计算出转子位置角度。
[0068] 图8所示为开关磁阻电机调速系统的给定转速为12化/min时,采用本发明算法时, 电机相电流、根据电流采样得到的电流斜率差、获取的等效相电感、根据电感信息辨识出的 转子位置波形。由图可知,采用注入法时在非导通相注入脉冲,相电流中可W明显表现出 来,电流斜率差波形和等效相电感波形变化趋势一致,得到的估计位置精度也较高。
[0069] 图9所示为开关磁阻电机调速系统的给定转速为25化/min时,电机a相电流波形, 等效相电感波形W及转子实际位置与估计位置波形。由图可知,转子位置估计值与实际值 吻合,辨识精度较高,能够满足SRD调速系统的要求,从电机相电流中可W看出电机运行平 稳,验证了脉冲注入法无位置传感器控制方法的有效性。
[0070] 图10所示为采用本发明方法,电机从静止状态启动到稳定运行时的动态波形。从 图中可知,在开关磁阻电机启动瞬间电机保持静止,对=相绕组注入脉冲,通过采样=相电 流,计算出电流斜率差值,得到=相等效电感进而估计出的转子初始位置,精度较高,电机 运行平稳。验证了脉冲注入法无位置传感器控制方法的动态性能。
[0071] 本发明通过对开关磁阻电机的非导通相注入脉冲,导通相采用斩波控制,根据相 电流斜率差值得到=相等效电感,从电感矢量中提取转子位置信息,进而可计算得到转子 速度,从而实现SRD控制,结构简单易于实现,计算量大大降低,易于实现实时控制,可靠性 和精度局,鲁棒性好。
[0072] 基于MATLAB/SIMULINK软件的仿真分析及硬件平台上的实验结果证明了本发明的 准确性和可行性。
[0073] 尽管结合优选实施方案具体展示和介绍了本发明,但所属领域的技术人员应该明 白,在不脱离所附权利要求书所限定的本发明的精神和范围内,在形式上和细节上可W对 本发明做出各种变化,均为本发明的保护范围。

Claims (6)

1. 基于脉冲注入的Ξ相开关磁阻电机无位置传感器控制方法,其特征在于:对开关磁 阻电机的非导通相注入脉冲,对导通相进行电流斩波控制,测量各相的相电流,计算出电流 斜率差值,通过电流斜率差值得出Ξ相等效电感,合成电感矢量,通过电感矢量估算转子位 置角度Θ,采用如下公式进行估算
Figure CN105703680AC00021
其中,目0为电感旋转角度,Nr为转子极数,La、Lb、Lc分别为Ξ相的各相电感。
2. 根据权利要求1所述的基于脉冲注入的Ξ相开关磁阻电机无位置传感器控制方法, 其特征在于:根据转子位置角度,计算出转子转速。
3. 根据权利要求2所述的基于脉冲注入的Ξ相开关磁阻电机无位置传感器控制方法, 其特征在于:转子转速采用W下公式进行计算
Figure CN105703680AC00022
其中,θη和θη+1为两个相邻采样周期内通过电感矢量计算得到的转子位置角度,ft是系 统采样频率,ω为转子速度。
4. 根据权利要求1所述的基于脉冲注入的Ξ相开关磁阻电机无位置传感器控制方法, 其特征在于:采用实验法获取最佳的脉冲宽度。
5. 根据权利要求1所述的基于脉冲注入的Ξ相开关磁阻电机无位置传感器控制方法, 其特征在于:对电流斜率差值进行数字滤波后,再计算得出Ξ相等效电感。
6. 根据权利要求5所述的基于脉冲注入的Ξ相开关磁阻电机无位置传感器控制方法, 其特征在于:所述数字滤波采用最大值滤波法。
CN201610143963.6A 2016-03-14 2016-03-14 基于脉冲注入的三相开关磁阻电机无位置传感器控制方法 CN105703680A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610143963.6A CN105703680A (zh) 2016-03-14 2016-03-14 基于脉冲注入的三相开关磁阻电机无位置传感器控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610143963.6A CN105703680A (zh) 2016-03-14 2016-03-14 基于脉冲注入的三相开关磁阻电机无位置传感器控制方法

Publications (1)

Publication Number Publication Date
CN105703680A true CN105703680A (zh) 2016-06-22

Family

ID=56221352

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610143963.6A CN105703680A (zh) 2016-03-14 2016-03-14 基于脉冲注入的三相开关磁阻电机无位置传感器控制方法

Country Status (1)

Country Link
CN (1) CN105703680A (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107359838A (zh) * 2017-08-28 2017-11-17 南京理工大学 一种基于有限元分析法的超高速永磁同步电机无传感器转速和位置估计方法
CN108075701A (zh) * 2017-05-25 2018-05-25 烟台仙崴机电有限公司 一种无位置传感器下开关磁阻发电机系统效率提升控制方法
US10079566B1 (en) 2017-07-14 2018-09-18 Caterpillar Inc. Pulse injection current pattern for rotor position estimation
CN109039208A (zh) * 2018-08-10 2018-12-18 南京信息工程大学 一种开关磁阻电机增量电感特性在线检测方法
CN110022098A (zh) * 2019-02-12 2019-07-16 湖南科技大学 基于相电感非饱和区定位的开关磁阻电机无位置传感器控制方法及装置
CN110120766A (zh) * 2019-05-13 2019-08-13 中国矿业大学 基于虚拟单位电感矢量的开关磁阻电机转子位置估计方法
WO2020151370A1 (zh) * 2018-10-26 2020-07-30 湖南科技大学 三相开关磁阻电机无位置传感器控制方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101699757A (zh) * 2009-10-28 2010-04-28 南京航空航天大学 一种适合低速的开关磁阻电机无位置传感器方法
JP2012228083A (ja) * 2011-04-20 2012-11-15 Yaskawa Electric Corp 交流電動機の制御装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101699757A (zh) * 2009-10-28 2010-04-28 南京航空航天大学 一种适合低速的开关磁阻电机无位置传感器方法
JP2012228083A (ja) * 2011-04-20 2012-11-15 Yaskawa Electric Corp 交流電動機の制御装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张旭隆等: "基于全周期电感空间矢量的开关磁阻电动机无位置传感器控制", 《工矿自动化》 *
蔡骏等: "基于全周期电感法的开关磁阻电机无位置传感器控制技术", 《电工技术学报》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108075701B (zh) * 2017-05-25 2020-06-26 烟台仙崴机电有限公司 一种无位置传感器下开关磁阻发电系统效率提升控制方法
CN108075701A (zh) * 2017-05-25 2018-05-25 烟台仙崴机电有限公司 一种无位置传感器下开关磁阻发电机系统效率提升控制方法
US10079566B1 (en) 2017-07-14 2018-09-18 Caterpillar Inc. Pulse injection current pattern for rotor position estimation
CN107359838A (zh) * 2017-08-28 2017-11-17 南京理工大学 一种基于有限元分析法的超高速永磁同步电机无传感器转速和位置估计方法
CN109039208A (zh) * 2018-08-10 2018-12-18 南京信息工程大学 一种开关磁阻电机增量电感特性在线检测方法
WO2020151370A1 (zh) * 2018-10-26 2020-07-30 湖南科技大学 三相开关磁阻电机无位置传感器控制方法及装置
CN110022098A (zh) * 2019-02-12 2019-07-16 湖南科技大学 基于相电感非饱和区定位的开关磁阻电机无位置传感器控制方法及装置
CN110120766A (zh) * 2019-05-13 2019-08-13 中国矿业大学 基于虚拟单位电感矢量的开关磁阻电机转子位置估计方法
CN110120766B (zh) * 2019-05-13 2020-11-27 中国矿业大学 基于虚拟单位电感矢量的开关磁阻电机转子位置估计方法

Similar Documents

Publication Publication Date Title
Luo et al. PMSM sensorless control by injecting HF pulsating carrier signal into estimated fixed-frequency rotating reference frame
US9270220B2 (en) Circuits and methods of determining position and velocity of a rotor
Tang et al. PMSM sensorless control by injecting HF pulsating carrier signal into ABC frame
Wang et al. A simple startup strategy based on current regulation for back-EMF-based sensorless control of PMSM
Xiao et al. Dynamic permanent magnet flux estimation of permanent magnet synchronous machines
Mizutani et al. Current model-based sensorless drives of salient-pole PMSM at low speed and standstill
CN102868350B (zh) 无位置传感器无刷直流电机准闭环启动方法
CN102170261B (zh) 交流电动机的驱动装置以及电动机车辆
Stirban et al. Motion-sensorless control of BLDC-PM motor with offline FEM-information-assisted position and speed observer
Dutta et al. A comparative analysis of two test methods of measuring $ d $-and $ q $-axes inductances of interior permanent-magnet machine
Consoli et al. A new zero-frequency flux-position detection approach for direct-field-oriented-control drives
US8400086B2 (en) Initial rotor position detection for permanent magnet synchronous motors
Damodharan et al. Sensorless brushless DC motor drive based on the zero-crossing detection of back electromotive force (EMF) from the line voltage difference
Rahman et al. Problems associated with the direct torque control of an interior permanent-magnet synchronous motor drive and their remedies
Rahman et al. Identification of machine parameters of a synchronous motor
Schmidt et al. Initial rotor angle detection of a nonsalient pole permanent magnet synchronous machine
CN101242154B (zh) 一种无位置传感器的内嵌式永磁无刷直流电机控制系统
US6441572B2 (en) Detection of rotor angle in a permanent magnet synchronous motor at zero speed
JP3419725B2 (ja) 位置センサレスモータ制御装置
Hu et al. Eddy current effects on rotor position estimation and magnetic pole identification of PMSM at zero and low speeds
JP2009055748A (ja) 電流検出ユニット及びモータ制御装置
CN101355337B (zh) 基于磁场正交控制的永磁同步电动机的驱动控制方法
CN101436843B (zh) 一种无位置传感器开关磁阻电机中转子位置的检测方法
CN106655942B (zh) 永磁同步电机无位置传感器控制方法
CN101286727B (zh) 电动机控制装置

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20160622

WD01 Invention patent application deemed withdrawn after publication