CN103809454B - 纵向飞行模型簇人机闭环复合频率鲁棒控制器设计方法 - Google Patents
纵向飞行模型簇人机闭环复合频率鲁棒控制器设计方法 Download PDFInfo
- Publication number
- CN103809454B CN103809454B CN201410070275.2A CN201410070275A CN103809454B CN 103809454 B CN103809454 B CN 103809454B CN 201410070275 A CN201410070275 A CN 201410070275A CN 103809454 B CN103809454 B CN 103809454B
- Authority
- CN
- China
- Prior art keywords
- frequency
- flight
- envelope
- loop
- aircraft
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000013461 design Methods 0.000 title claims abstract description 30
- 238000012360 testing method Methods 0.000 claims abstract description 6
- 238000012546 transfer Methods 0.000 claims description 13
- 238000007476 Maximum Likelihood Methods 0.000 claims description 6
- 239000002131 composite material Substances 0.000 claims description 5
- 230000003068 static effect Effects 0.000 claims description 4
- 241000764238 Isis Species 0.000 claims description 3
- 230000010355 oscillation Effects 0.000 description 5
- 238000004088 simulation Methods 0.000 description 4
- 238000013178 mathematical model Methods 0.000 description 3
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241001633106 Lithocarpus Species 0.000 description 1
- 101100348848 Mus musculus Notch4 gene Proteins 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000013643 reference control Substances 0.000 description 1
- 238000013112 stability test Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Landscapes
- Feedback Control In General (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明提供了一种纵向飞行模型簇人机闭环复合频率鲁棒控制器设计方法,该方法在给定不同高度、马赫数条件下通过扫频飞行试验直接确定获得全包线内的幅频和相频特性构成的模型簇;根据飞行包线内的幅频特性直接确定开环截止频率区间;根据飞行包线内的相频特性直接确定与截止频率区间所对应的相位裕度区间;通过加入多级串联滞后—超前补偿环节控制器并在飞行器全包线内的相位裕度指标和系统辨识中的模型辨识方法确定补偿环节个数和参数值;在飞行器全飞行包线内的幅值裕度指标
Description
技术领域
本发明涉及一种飞行器控制器设计方法,特别涉及纵向飞行模型簇人机闭环复合频率鲁棒控制器设计方法,属于测控技术和飞行力学等范畴。
背景技术
飞行器起降过程的控制对飞行安全有重要作用;由于飞行器起降过程中飞行速度变化大,即使按照纵向模型也会面临强非线性问题;另一方面,飞行器的操纵舵存在饱和、死区等现象;从飞行安全考虑,超低空飞行(如飞机起飞/着陆)时,控制器必须保证系统具有一定的稳定裕度、无超调和平稳性,这样,就使得超低空飞行控制器设计非常复杂,不能直接套用现有控制理论进行飞行器控制的设计。
在现代实际飞行控制器的设计中,一少部分采用状态空间法进行设计,而大多数仍然采用以PID为代表的经典频域法和逆Nyquist阵列法为代表的现代频率法进行控制器设计。现代控制理论以状态空间法为特征、以解析计算为主要手段、以实现性能指标为最优的现代控制理论,而后又展了最优控制方法、模型参考控制方法、自适应控制方法、动态逆控制方法,反馈线性化方法、直接非线性优化控制、变增益控制法、神经网络控制方法,模糊控制方法,鲁棒控制法以及多种方法组合控制等一系列控制器设计方法,发表的学术论文数以万计,例如2011年GhasemiA设计了自适应模糊滑模控制的再入飞行器(GhasemiA,MoradiM,MenhajMB.AdaptiveFuzzySlidingModeControlDesignforaLow-LiftReentryVehicle[J].JournalofAerospaceEngineering,2011,25(2):210-216),2013年BabaeiAR为非最小相位和非线性飞行器设计了模糊滑模控制自动驾驶仪(BabaeiAR,MortazaviM,MoradiMH.FuzzyslidingmodeautopilotdesignfornonminimumphaseandnonlinearUAV[J].JournalofIntelligentandFuzzySystems,2013,24(3):499-509),很多研究仅仅停留在理想化的仿真研究阶段;而且这种设计存在三个问题:(1)由于无法进行飞行器超低空操纵稳定性试验,难以得到精确的被控对象的数学模型;(2)对于军标规定的稳定裕度等评价飞行控制系统的重要性能指标,状态空间法远不像经典频率法那样能以明显的形式表达出来;(3)控制器结构过于复杂、没有考虑实际控制器和飞行状态的约束,设计的控制器物理上不可实现。
英国的学者Rosenbrock系统地、开创性地研究了如何将频域法推广到多变量系统的设计中去,利用矩阵对角优势概念,把多变量问题转化为能用人们熟知的古典方法的单变量系统的设计问题,以后相继出现了Mayne序列回差法,MacFarlane特征轨迹法、Owens并矢展开法等方法,共同特点是把多输入一多输出、回路间严重关联的多变量系统的设计,化为一系列单变量系统的设计问题,进而可选用某一种古典方法(Nyquist和Bode的频率响应法,Evans的根轨迹法等)完成系统的设计,上述这些方法保留和继承了古典图形法的优点,不要求特别精确的数学模型,容易满足工程上的限制。特别是当采用有图形显示终端的人一机对话式的计算机辅助设计程序实现时,可以充分发挥设计者的经验和智慧,设计出既满足品质要求,又是物理上可实现的、结构简单的控制器;国内外对多变量频率法进行了改进研究(高大远,罗成,沈辉,胡德文,挠性卫星姿态解藕控制器多变量频率域设计方法,宇航学报,2007,Vol.28(2),pp442-447;熊柯,夏智勋,郭振云,倾斜转弯高超声速巡航飞行器多变量频域法解耦设计,弹箭与制导学报,2011,Vol.31(3),pp25-28)但是,这种设计方法当考虑系统不确定问题时保守性过大,在飞行器操纵舵限制情况下不能得到合理的设计结果。
在高性能飞机的研制中,评价一架飞机飞行品质的好坏,不仅取决于飞机本身和驾驶员操纵的动力学特性,还取决于驾驶员与飞机之间高度一致地配合,以及驾驶员与先进飞机飞行控制系统之间功能分配的合理性。直到1980年后,评价飞机飞行品质的美军标仍然存在着严重缺陷,即没有考虑驾驶员在操纵回路中的作用,因而由此所得的评价与试飞员试飞后所得的结果仍有一定差距。近年来,正在发展一种有驾驶员参与系统的闭环准则(尼尔2史密斯准则),但如何实现至今仍无有效算法。尼尔—史密斯准则是在1970年提出的,它是一个闭环俯仰跟踪准则。它考虑问题的方法是:在驾驶员驾驶飞机时与飞机构成一个闭环系统,驾驶员轻松地操纵就能达到特定的飞行指标,则飞行品质是好的。为了获得与驾驶员一致的对飞机的评价意见,则在理论分析中必须把驾驶员包括在内。通常,驾驶员行为的数学模型是非线性的,可能是离散的,但当研究具有稳定性的操纵对象时,有用的近似模型仍为线性。由大量飞行实践和仿真研究表明,驾驶员的行为由他的心理特性、生理特性、周围环境、操纵系统、操纵对象所要完成的任务来决定,尽管驾驶员有着各自的特点,但在完成单一的飞行任务中,大多数驾驶员的动作可由完全确定的数学模型来描述,它是驾驶员行为的大量试验的平均状态,与实际情况很接近,因此现在大都采用如下形式的人机闭环特性驾驶员模型:
和系统开环传递函数或频率特性来估计人机闭环特性。
其中:为驾驶员环节的静态增益、为驾驶员的固有延时特性、为驾驶员超前补偿时间常数、为驾驶员滞后补偿时间常数;此模型加入后,在飞行控制器中要考虑反应较慢的驾驶员带来的诱发振荡问题,在全部飞行包线内,对反应较慢的驾驶员容许控制器设计还没有系统的方法,只是对单一飞行状态具有部分研究。
综上所述,目前的控制方法还不能在飞行器模型变化、按照全飞行包线内的稳定裕度指标设计出无驾驶员诱发振荡、超调量小、平稳的低空飞行控制器。
发明内容
为了克服现有方法不能在飞行器在全飞行包线内模型变化大的情况下设计出符合全飞行包线内的稳定裕度指标的无驾驶员诱发振荡、超调量小、平稳低空飞行控制器的技术缺陷,本发明提供了一种纵向飞行模型簇人机闭环复合频率鲁棒控制器设计方法,该方法在给定不同高度、马赫数条件下通过扫频飞行试验直接确定获得全包线内的幅频和相频特性构成的模型簇;根据飞行包线内的幅频特性直接确定开环截止频率区间;根据飞行包线内的相频特性直接确定与截止频率区间所对应的相位裕度区间;通过加入多级串联滞后—超前补偿环节控制器并在飞行器全包线内的相位裕度指标和系统辨识中的模型辨识方法确定补偿环节个数和参数值;在飞行器全飞行包线内的幅值裕度指标分贝数给定情况下进行控制器效果验证;从相位裕度和幅值裕度的概念出发设计出符合全飞行包线的无驾驶员诱发振荡、超调量小、平稳的低空飞行鲁棒控制器。
本发明解决其技术问题所采用的技术方案:一种纵向飞行模型簇人机闭环复合频率鲁棒控制器设计方法,其特征包括以下步骤:
1、给定不同高度、马赫数下通过扫频飞行试验直接由允许飞行的全包线内的幅频和相频特性构成飞行器全包线内的升降舵与飞行高度之模型簇,对应的飞行器升降舵与飞行高度之间开环传递函数簇描述为:
其中
、
为多项式,为传递函数中常用的拉普拉斯变化后的变量,分别为飞行高度和马赫数,是俯仰回路的延迟时间,为随变化的增益,为多项式中随变化的系数簇,为多项式中随变化的系数簇,为模型中的不确定项;
考虑人机闭环特性时驾驶员模型:
和系统开环传递函数或频率特性来估计人机闭环特性;
其中:为驾驶员环节的静态增益、为驾驶员的固有延时特性、为驾驶员超前补偿时间常数、为驾驶员滞后补偿时间常数;
这样,考虑驾驶员特性的人机系统的开环模型就变为:
其中:;
2、判断在已知模型不确定部分时,根据飞行包线内的幅频特性直接确定开环截止频率区间确定方法为:
从即中,近似为,得到开环截止频率解的最大值和最小值,开环截止频率区间为;
式中,为正实数,为频率特性中的变量,为虚部表示,为角频率;
3、判断在已知模型不确定部分时,根据飞行包线内的相频特性,计算包线内最大相位裕度:
和包线内最小相位裕度:
,
直接确定与截止频率区间所对应的相位裕度区间为:
;
其中,为正实数角度数;
4、候选多级串联滞后—超前补偿环节的传递函数为:
式中,为待确定的常数增益,N为整数,表示待确定的滞后—超前补偿环节的级数,、、、为待确定的时间常数,为待确定的参数;
加入多级串联滞后—超前补偿环节后,
从即
中,得到开环截止频率解的最大值和最小值,开环截止频率区间为,
在飞行器全包线内的相位裕度指标给定情况下,加入多级串联滞后—超前补偿环节后系统的相位裕度应该满足:
即满足:
(1)
在指标(1)式和极大似然准则共同约束下,根据系统模型结构辨识中的极大似然方法确定滞后—超前补偿环节的级数N、常数增益、时间常数、、、,待确定的参数;
5、在飞行器全包线内的幅值裕度指标分贝数给定情况下,
根据,采用近似表达:
中,得到频率解的最大值和最小值,区间为,
判断:
即满足:
若满足,则飞行控制器设计完成,若不满足,再增加补偿环节级数或减小常数增益。
本发明的有益效果是:从相位裕度和幅值裕度的概念出发,通过加入多级串联滞后—超前补偿环节控制器,在全飞行包线内按照符合给定相位裕度和幅值裕度的要求和模型辨识方法确定多级串联滞后—超前补偿环节鲁棒控制器的参数,设计出符合全飞行包线的无驾驶员诱发振荡、超调量小、平稳的低空飞行鲁棒控制器。
下面结合实施例对本发明作详细说明。
具体实施方式
1、给定不同高度、马赫数下使用线性扫频信号(为起始频率,为截止频率,,为扫频时间)或对数扫频信号
(为起始频率,为截止频率,,T为扫频时间)对飞机激励,可直接得到允许飞行的全包线内的幅频和相频特性,构成飞行器全包线内的升降舵与飞行高度之模型簇,对应的飞行器升降舵与飞行高度之间开环传递函数簇描述为:
其中
、
为多项式,为传递函数中常用的拉普拉斯变化后的变量,分别为飞行高度和马赫数,是俯仰回路的延迟时间,为随变化的增益,为多项式中随变化的系数簇,为多项式中随变化的系数簇,为模型中的不确定项;
考虑人机闭环特性时驾驶员模型:
和系统开环传递函数或频率特性来估计人机闭环特性;
其中:为驾驶员环节的静态增益、为驾驶员的固有延时特性、为驾驶员超前补偿时间常数、为驾驶员滞后补偿时间常数;
这样,考虑驾驶员特性的人机系统的开环模型就变为:
其中:;
2、判断在已知模型不确定部分时,根据飞行包线内的幅频特性直接确定开环截止频率区间确定方法为:
从即中,近似为,得到开环截止频率解的最大值和最小值,开环截止频率区间为;
式中,为正实数,为频率特性中的变量,为虚部表示,为角频率;
3、判断在已知模型不确定部分时,根据飞行包线内的相频特性,计算包线内最大相位裕度:
和包线内最小相位裕度:
,
直接确定与截止频率区间所对应的相位裕度区间为:
;
其中,为正实数角度数;
4、候选多级串联滞后—超前补偿环节的传递函数为:
式中,为待确定的常数增益,N为整数,表示待确定的滞后—超前补偿环节的级数,、、、为待确定的时间常数,为待确定的参数;
加入多级串联滞后—超前补偿环节后,
从即
中,得到开环截止频率解的最大值和最小值,开环截止频率区间为,
在飞行器全包线内的相位裕度指标给定情况下,加入多级串联滞后—超前补偿环节后系统的相位裕度应该满足:
即满足:
(1)
在指标(1)式和极大似然准则共同约束下,根据系统模型结构辨识中的极大似然方法确定滞后—超前补偿环节的级数N、常数增益、时间常数、、、,待确定的参数;
5、在飞行器全包线内的幅值裕度指标分贝数给定情况下,
根据,采用近似表达:
中,得到频率解的最大值和最小值,区间为,
判断:
即满足:
若满足,则飞行控制器设计完成,若不满足,再增加补偿环节级数或减小常数增益。
Claims (1)
1.一种纵向飞行模型簇人机闭环复合频率鲁棒控制器设计方法,其特征包括以下步骤:
(1)给定不同高度、马赫数下通过扫频飞行试验直接由允许飞行的全包线内的幅频和相频特性构成飞行器全包线内的升降舵与飞行高度之模型簇,对应的飞行器升降舵与飞行高度之间开环传递函数簇描述为:
其中
、
为多项式,为传递函数中常用的拉普拉斯变化后的变量,分别为飞行高度和马赫数,是俯仰回路的延迟时间,为随变化的增益,为多项式中随变化的系数簇,为多项式中随变化的系数簇,为模型中的不确定项;
考虑人机闭环特性时驾驶员模型:
和系统开环传递函数或频率特性来估计人机闭环特性;
其中:为驾驶员环节的静态增益、为驾驶员的固有延时特性、为驾驶员超前补偿时间常数、为驾驶员滞后补偿时间常数;
这样,考虑驾驶员特性的人机系统的开环模型就变为:
其中:;
(2)判断在已知模型不确定部分时,根据飞行包线内的幅频特性直接确定开环截止频率区间确定方法为:
从即中,近似为,得到开环截止频率解的最大值和最小值,开环截止频率区间为;
式中,为正实数,为频率特性中的变量,为虚部表示,为角频率;
(3)判断在已知模型不确定部分时,根据飞行包线内的相频特性,计算包线内最大相位裕度:
和包线内最小相位裕度:
,
直接确定与截止频率区间所对应的相位裕度区间为:
;
其中,为正实数角度数;
(4)候选多级串联滞后—超前补偿环节的传递函数为:
式中,为待确定的常数增益,N为整数,表示待确定的滞后—超前补偿环节的级数,、、、为待确定的时间常数,为待确定的参数;
加入多级串联滞后—超前补偿环节后,
从即
中,得到开环截止频率解的最大值和最小值,开环截止频率区间为,
在飞行器全包线内的相位裕度指标给定情况下,加入多级串联滞后—超前补偿环节后系统的相位裕度应该满足:
即满足:
(1)
在指标(1)式和极大似然准则共同约束下,根据系统模型结构辨识中的极大似然方法确定滞后—超前补偿环节的级数N、常数增益、时间常数、、、,待确定的参数;
(5)在飞行器全包线内的幅值裕度指标分贝数给定情况下,
根据,采用近似表达:
中,得到频率解的最大值和最小值,区间为,
判断:
即满足:
若满足,则飞行控制器设计完成,若不满足,再增加补偿环节级数或减小常数增益。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410070275.2A CN103809454B (zh) | 2014-02-28 | 2014-02-28 | 纵向飞行模型簇人机闭环复合频率鲁棒控制器设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410070275.2A CN103809454B (zh) | 2014-02-28 | 2014-02-28 | 纵向飞行模型簇人机闭环复合频率鲁棒控制器设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103809454A CN103809454A (zh) | 2014-05-21 |
CN103809454B true CN103809454B (zh) | 2016-06-15 |
Family
ID=50706404
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410070275.2A Expired - Fee Related CN103809454B (zh) | 2014-02-28 | 2014-02-28 | 纵向飞行模型簇人机闭环复合频率鲁棒控制器设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103809454B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2151905A2 (en) * | 2008-08-05 | 2010-02-10 | The Boeing Company | Optimizing usage of powered systems |
CN102707624A (zh) * | 2012-05-31 | 2012-10-03 | 西北工业大学 | 基于飞行器常规模型的纵向控制器区域设计方法 |
CN102929130A (zh) * | 2012-10-10 | 2013-02-13 | 西北工业大学 | 一种鲁棒飞行控制器设计方法 |
CN103116280A (zh) * | 2013-01-16 | 2013-05-22 | 北京航空航天大学 | 一种存在变分布网络随机延迟的微小型无人飞行器纵向控制方法 |
-
2014
- 2014-02-28 CN CN201410070275.2A patent/CN103809454B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2151905A2 (en) * | 2008-08-05 | 2010-02-10 | The Boeing Company | Optimizing usage of powered systems |
CN102707624A (zh) * | 2012-05-31 | 2012-10-03 | 西北工业大学 | 基于飞行器常规模型的纵向控制器区域设计方法 |
CN102929130A (zh) * | 2012-10-10 | 2013-02-13 | 西北工业大学 | 一种鲁棒飞行控制器设计方法 |
CN103116280A (zh) * | 2013-01-16 | 2013-05-22 | 北京航空航天大学 | 一种存在变分布网络随机延迟的微小型无人飞行器纵向控制方法 |
Non-Patent Citations (3)
Title |
---|
ADRC FRACTIONAL ORDER PID CONTROLLER DESIGN OF HYPERSONIC FLIGHT VEHICLE;Qin Changmao,et;《Transactions of Nanjing University of Aeronautics & Astronautics》;20110915;第28卷(第3期);第240-245页 * |
基于飞行品质评估的无人直升机鲁棒控制器设计;刘鹏等;《航空学报》;20120520;第33卷(第9期);第1587-1597页 * |
飞行器模型簇描述及鲁棒控制器设计;史忠科;《控制与决策》;20040830;第19卷(第8期);第911-914、926页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103809454A (zh) | 2014-05-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105607473B (zh) | 小型无人直升机的姿态误差快速收敛自适应控制方法 | |
Durmaz et al. | Sliding mode control for non-linear systems with adaptive sliding surfaces | |
Qu et al. | Adaptive output-feedback control with closed-loop reference models for very flexible aircraft | |
Bierling | Comparative analysis of adaptive control techniques for improved robust performance | |
Cordeiro et al. | Robustness of incremental backstepping flight controllers: The boeing 747 case study | |
CN103809446B (zh) | 飞行器多回路模型簇颤振抑制复合频率鲁棒控制器设计方法 | |
CN103809442B (zh) | 飞行器多回路模型簇复合频率鲁棒控制器设计方法 | |
Moin et al. | State space model of an aircraft using Simulink | |
CN103777523B (zh) | 飞行器多回路模型簇复合pid鲁棒控制器设计方法 | |
CN103809456B (zh) | 纵向飞行模型簇人机闭环pid鲁棒控制器设计方法 | |
Zhao et al. | Prescribed performance tracking control for hypersonic flight vehicles with model uncertainties | |
CN103853048B (zh) | 飞行器多回路模型簇人机闭环复合频率鲁棒控制器设计方法 | |
CN103809433B (zh) | 飞行器多回路模型簇复合根轨迹多级pid鲁棒控制器设计方法 | |
CN103809448B (zh) | 飞行器多回路模型簇颤振抑制复合根轨迹补偿鲁棒控制器设计方法 | |
CN103809444B (zh) | 飞行器多回路模型簇人机闭环pid鲁棒控制器设计方法 | |
CN103809454B (zh) | 纵向飞行模型簇人机闭环复合频率鲁棒控制器设计方法 | |
CN103853049B (zh) | 纵向飞行模型簇复合频率鲁棒控制器设计方法 | |
CN103823367B (zh) | 纵向飞行模型簇颤振抑制复合频率鲁棒控制器设计方法 | |
CN103809449B (zh) | 飞行器多回路模型簇颤振抑制复合pid鲁棒控制器设计方法 | |
CN103809445B (zh) | 飞行器多回路模型簇复合pid控制器设计方法 | |
Canal et al. | Ziegler–Nichols customization for quadrotor attitude control under empty and full loading conditions | |
CN103792848B (zh) | 纵向飞行模型簇人机闭环复合根轨迹多级pid鲁棒控制器设计方法 | |
CN103823364B (zh) | 飞行器多回路模型簇复合根轨迹补偿鲁棒控制器设计方法 | |
CN103823374B (zh) | 飞行器多回路模型簇复合根轨迹补偿控制器设计方法 | |
CN103809443B (zh) | 飞行器多回路模型簇人机闭环复合根轨迹多级pid鲁棒控制器设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160615 |