CN102663195A - 一种降低旋翼桨叶振动和动应力水平的设计方法 - Google Patents

一种降低旋翼桨叶振动和动应力水平的设计方法 Download PDF

Info

Publication number
CN102663195A
CN102663195A CN2012101130179A CN201210113017A CN102663195A CN 102663195 A CN102663195 A CN 102663195A CN 2012101130179 A CN2012101130179 A CN 2012101130179A CN 201210113017 A CN201210113017 A CN 201210113017A CN 102663195 A CN102663195 A CN 102663195A
Authority
CN
China
Prior art keywords
blade
vibration
rotor blade
design
level
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.)
Granted
Application number
CN2012101130179A
Other languages
English (en)
Other versions
CN102663195B (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.)
Beihang University
Original Assignee
Beihang 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 Beihang University filed Critical Beihang University
Priority to CN201210113017.9A priority Critical patent/CN102663195B/zh
Publication of CN102663195A publication Critical patent/CN102663195A/zh
Application granted granted Critical
Publication of CN102663195B publication Critical patent/CN102663195B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Turbine Rotor Nozzle Sealing (AREA)
  • Structures Of Non-Positive Displacement Pumps (AREA)

Abstract

本发明公开了一种降低旋翼桨叶振动和动应力水平的设计方法,属于直升机技术领域。本发明首先根据直升机旋翼桨叶提取设计变量,根据直升机旋翼强度要求确定设计变量的上限和下限;然后以直升机旋翼桨叶振动水平指数作为目标函数,旋翼桨叶的振型节点位置要求为约束函数,建立降低直升机旋翼桨叶振动和动应力水平的桨叶优化设计模型;最后对得到的桨叶优化设计模型进行求解,输出桨叶设计参数,得到桨叶振动水平指数。通过本发明提出的振动水平指数设计的桨叶能够明显降低桨根切力,减小直升机响应振动水平。本发明提出的设计方法加入了固有频率和振型节点的要求,因此通过此优化设计方法得到的桨叶还能够改善桨叶的疲劳状况。

Description

一种降低旋翼桨叶振动和动应力水平的设计方法
技术领域
本发明属于直升机技术领域,具体涉及一种降低直升机旋翼桨叶振动和动应力水平的设计方法,可用于直升机低振动和动应力水平旋翼桨叶的设计。
背景技术
旋翼是直升机上最为关键部件,并为直升机提供升力、推力和操作。而由于旋翼的周期性运动,使得直升机出现了严重的振动问题。乘员舒适性与任务工作效能、机体与构件的疲劳寿命、机上系统与设备的功能与可靠性和维护费用等均会受到振动的巨大影响。此外,振动还会限制直升机前飞速度与机动性能等。因此减振是进行直升机设计时需要考虑的重要问题之一。为了保障直升机的整体性能,各国航空组织对其减振的要求也越来越高,减振技术也成为了直升机设计的重点和难点问题。直升机减振技术的研究主要集中在三个方面,首先是控制和减小机身结构的振动;其次是减小旋翼传递到机身结构上的激振力;最后是控制旋翼本身的激振力,由于旋翼是直升机产生振动问题的根源,所以降低旋翼本身的激振力是减振最理想的目标。
如何降低旋翼的激振力是一个极其复杂的问题,从80年代开始,各国学者对此开展了大量的研究,并得到了桨叶动力学调频、桨叶模态优化和极小化浆毂载荷等方法。在这些方法中大部分都围绕在旋翼桨叶设计中如何降低浆毂切力和力矩进行。通过研究表明合理选择频率和阵型能够使根切力减小20%-40%,但是到目前为止,仍然难以得到一个同时考虑挥舞、摆振和扭转及其耦合运动的结构动力学模型和一个能够准确模拟复杂空气动力环境的气动模型,因此还无法在实际桨叶设计中应用。考虑到固有频率是桨叶的特有属性,并不涉及到气动力和桨叶振型,因此大量研究从如何避免共振的角度,通过合理控制桨叶的固有频率以达到减振的目的,并且已经成为了目前重要的设计准则,但是在一般情况下固有频率控制不能有效的控制浆毂力大小,也就无法有效的达到减振和降低动应力水平的目的。
发明内容
本发明针对现有直升机减振设计方法的不足,考虑到直升机旋翼桨叶减振和降低动应力水平问题的复杂性,直接从工程应用的角度,引入了振动水平指数和振型节点位置控制等方法,提出了在桨叶外形参数和浆毂型式确定情况下降低直升机振动和动应力水平的设计方法。应用此方法设计桨叶展向质量和刚度分布,不但使桨根交变切力最小,达到降低直升机振动水平的目的,还能够控制动应力大小改善桨叶疲劳情况。
本发明的技术解决方案:一种降低旋翼桨叶振动和动应力水平的设计方法,通过如下步骤实现:
(1)根据直升机旋翼桨叶提取设计变量,包括桨叶沿展向的质量和刚度分布,根据直升机旋翼强度要求确定设计变量的上限和下限;
(2)以本发明提出的直升机旋翼桨叶振动水平指数作为目标函数,旋翼桨叶的振型节点位置要求为约束函数,建立降低直升机旋翼桨叶振动和动应力水平的桨叶优化设计模型;
(3)对得到的桨叶优化设计模型进行求解,输出桨叶设计参数,得到桨叶振动水平指数。
本发明的特别之处在于基于工程应用的目的提出了振动水平指数,并且在采用此指数对旋翼桨叶进行减振设计的同时,还考虑了振型节点位置的要求,从而建立了一套新的低振动和动应力水平旋翼桨叶的优化设计方法。这样就可以达到既能降低桨根切力,减小直升机旋翼的振动水平,又能改善疲劳状况的目的。
本发明的优点在于:
(1)通过本发明提出的振动水平指数设计的桨叶能够明显降低桨根切力,减小直升机响应振动水平。
(2)本发明提出的设计方法加入了振型节点位置的要求,因此通过此优化设计方法得到的桨叶还能够改善桨叶的疲劳状况。
附图说明
图1是本发明的降低直升机旋翼振动和动应力水平设计方法流程图;
图2是原始桨叶和本发明设计桨叶的振型节点位置对比示意图。
具体实施方式
下面结合附图和实施例对本发明提供的设计方法进行详细说明。
如图1所示为本发明的降低旋翼桨叶振动和动应力水平设计方法流程图,具体设计步骤如下:
第一步,提取直升机旋翼桨叶的设计变量,包括桨叶沿展向的质量和刚度分布,根据直升机旋翼桨叶的强度要求确定设计变量的上限和下限。
第二步,以直升机旋翼桨叶振动水平指数作为目标函数,旋翼桨叶的振型节点位置要求为约束函数,建立降低直升机旋翼桨叶振动和动应力水平的桨叶优化设计模型。
本发明在低振动旋翼动力学优化设计方法中,首次提出并采用了振动水平指数做为目标函数,并以此指数为目标函数,振型节点位置作为约束函数建立了一套新的降低旋翼振动和动应力水平的设计方法。下面对本发明提出的旋翼桨叶的优化设计方法中的目标函数和约束函数进行详细说明。
(1)本发明采用了反映直升机旋翼振动水平的振动水平指数作为目标函数。旋翼桨叶中,各阶气动力在桨根位置均会引起桨根切力,为了达到降低旋翼桨叶振动水平的目的,必须减小kΩ谐波气动力所引起的桨根切力分量值。要得到根切力大小,首先必须知道气动力的情况,此处的核心思想是通过模态参数近似表示广义气动力,进而推导出一个新的衡量旋翼桨叶振动的参数。首先得到各阶气动力引起的桨根切力,在桨叶外形参数和浆毂形式确定的情况下,kΩ谐波气动力中的第i阶广义气动力Fi与模态参数可以采用近似的正比例关系表示,即下式:
Figure BDA0000153975070000031
其中,r为桨叶的展向距离,
Figure BDA0000153975070000032
是桨叶的第i阶固有振型。
得到气动力之后,通过叠加的方法可以得到kΩ谐波气动力分量引起的桨根切力分量的近似比例,采用振动水平指数Vindex来定义,其数学模型如下:
Figure BDA0000153975070000033
其中模态阻尼比:
Figure BDA0000153975070000034
模态质量Mi的表达式为:
Figure BDA0000153975070000035
式中ωi为第i阶固有圆频率,为固有频率的2π倍,
Figure BDA0000153975070000036
是桨叶的第i阶固有振型,m为桨叶质量,R为旋翼桨盘半径,r为桨叶的展向距离,Ω为旋翼转速,n为所用模态阶数,k为谐波气动力阶数,ρ为空气密度,c为桨叶弦长,a为桨叶二元升力线的斜率。
通过此振动水平指数Vindex,使得旋翼振动水平成为了模态参数的函数,并取决于桨叶的模态参数。通过引入振动水平指数,在进行低振动水平桨叶设计过程中除了调整固有圆频率和模态阻尼比,还可以通过改变固有振型减小广义气动力。而展向质量和刚度分布的设计参数等会直接影响旋翼桨叶的固有振型,因此通过调整设计参数能够改变振型并减小气动力,从而达到降低振动水平的目的。
(2)桨叶优化设计模型中采用了振型节点位置要求作为约束函数,固有频率及固有振型是桨叶结构的固有属性,当激振力频率与桨叶结构固有频率接近时,结构振动达到最大值。而振型节点位置是结构动力学设计中的重要内容,通过调整振型节点位置能够有效降低振动水平,改善疲劳状况。在桨叶优化设计模型中采用了振型节点位置作为约束函数,以改善旋翼桨叶的疲劳状况。如图2为振型节点位置示意图,O为原点位置,P为原始桨叶第i阶振型节点位置,P′为设计要求的第i阶振型节点位置,原始桨叶振型节点位置P与原点O之间的距离为
Figure BDA0000153975070000041
而设计要求的第i阶振型节点位置P′与原点O之间的距离为
Figure BDA0000153975070000042
则满足设计要求时,两者之间的距离di应该满足下式:
d i = | X i u - X i o | < &epsiv;
式中ε为设计要求的设计误差,当两者之间的距离di小于设计误差时,可以认为满足设计要求。
根据以上的设计参数、目标函数和约束函数得到的优化模型如下:
目标函数:
Figure BDA0000153975070000044
振型节点位置约束函数: d i = | X i u - X i o | < &epsiv; , ( i = 1,2,3 . . . . . . s )
质量约束函数: W ( D &OverBar; ) &le; W 0
自旋惯量约束函数: I ( D &OverBar; ) &GreaterEqual; cI 0
设计变量约束函数: D i L &le; D i &le; D i u , ( i = 1,2,3 . . . . . . n )
其中为目标函数,即振动水平指数;di<ε为振型节点位置要求约束,s为考虑的模态振型的阶数;
Figure BDA00001539750700000410
为设计变量,即展向质量和刚度分布;
Figure BDA00001539750700000411
为质量约束函数,即桨叶总质量不得大于原始桨叶的总质量W0
Figure BDA00001539750700000412
为自旋惯量约束函数,即自旋惯量不得小于原始桨叶自旋惯量值I0的c倍,c可以根据设计需要而定;最后为设计变量约束函数Di,通过强度要求可以得到设计变量上限
Figure BDA00001539750700000413
和下限
Figure BDA00001539750700000414
第三步,对上述得到的桨叶优化设计模型进行求解,即可得到满足设计要求的桨叶结构,输出桨叶设计参数,得到桨叶振动水平指数。即得到了能够同时降低旋翼桨叶振动和动应力水平的桨叶结构。其中,求解方法可以采用基于灵敏度分析的优化设计方法,也可以采用现有的任何其它方法来求解。
实施例
本例中针对一种铰接式复合材料旋翼桨叶采用本发明提出的降低旋翼桨叶振动和动应力水平的优化设计方法进行设计。所采用的旋翼中共有4片桨叶,直径D=2m,水平铰外伸量Lsp=0.0495m,无几何扭转,转速Ω=800r/min,翼型为NACA0015。这里将本发明提出来的设计方法与其他设计得到的桨叶振动水平进行对比分析,能够更好的反映出本发明的设计方法的有效性。除了本发明提出的方法之外,还考虑了控制固有频率、控制挥舞一阶振型节点位置和最小垂直切力振动水平指数三种设计方法。桨叶沿展向的质量和刚度分布设计变量的设计区间为初始值的75%到250%,自旋惯量约束中c取为90%。由此得到的原始桨叶及满足各种设计要求时桨叶前四阶挥舞频率及相应的垂直切力的振动水平指数。
表1各设计情况下的桨叶振动水平指数比较
Figure BDA00001539750700000415
Figure BDA0000153975070000051
由表1可知,采用振动水平指数方法和本发明提出设计方法所得到桨叶前四阶模态的振动水平有较大程度的下降。比较振动水平指数方法和本发明的设计方法可以发现,尽管本发明提出的设计方法得到第四阶模态的振动水平指数较采用振动水平指数方法得到的结果要大,但是前三阶均有明显的降低,因此本发明方法明显优于其它三种方法。通过验证,采用本发明提出方法设计的桨叶能够使桨根切力幅值降低约20%,而单纯的频率设计方法则没有改进。此外由于本发明的方法还控制了振型节点位置,由此设计得到的桨叶还能改善疲劳状况。

Claims (2)

1.一种降低旋翼桨叶振动和动应力水平的设计方法,其特征在于:
(1)根据直升机旋翼桨叶提取设计变量,包括桨叶沿展向的质量和刚度分布,根据直升机旋翼强度要求确定设计变量的上限和下限;
(2)以本发明提出的直升机旋翼桨叶振动水平指数作为目标函数,旋翼桨叶的振型节点位置要求为约束函数,建立降低直升机旋翼桨叶振动和动应力水平的桨叶优化设计模型;
所述的桨叶优化设计模型如下:
目标函数:
Figure FDA0000153975060000011
振型节点位置约束函数: d i = | X i u - X i o | < &epsiv; , ( i = 1,2,3 . . . . . . s ) ;
质量约束函数: W ( D &OverBar; ) &le; W 0 ;
自旋惯量约束函数: I ( D &OverBar; ) &GreaterEqual; cI 0 ;
设计变量约束函数: D i L &le; D i &le; D i u , ( i = 1,2,3 . . . . . . n ) ;
其中
Figure FDA0000153975060000016
为目标函数,即振动水平指数;di<ε为振型节点位置要求约束,s为考虑的模态振型的阶数;
Figure FDA0000153975060000017
为设计变量,即展向质量和刚度分布;
Figure FDA0000153975060000018
为质量约束函数,即桨叶总质量不得大于原始桨叶的总质量W0为自旋惯量约束函数,即自旋惯量不得小于原始桨叶自旋惯量值I0的c倍,c根据设计需要而定;最后为设计变量约束函数Di,通过强度要求得到设计变量上限
Figure FDA00001539750600000110
和下限
Figure FDA00001539750600000111
(3)对得到的桨叶优化设计模型进行求解,输出桨叶设计参数,得到桨叶振动水平指数。
2.根据权利要求1所述的一种降低旋翼桨叶振动和动应力水平的设计方法,其特征在于:
所述的振动水平指数Vindex如下:
Figure FDA00001539750600000112
其中模态阻尼比:
Figure FDA00001539750600000113
模态质量Mi的表达式为:
Figure FDA00001539750600000114
式中ωi为第i阶固有圆频率,是桨叶的第i阶固有振型,m为桨叶质量,R为旋翼桨盘半径,r为桨叶的展向距离,Ω为旋翼转速,n为所用模态阶数,k为谐波气动力阶数,ρ为空气密度,c为桨叶弦长,a为桨叶二元升力线的斜率。
CN201210113017.9A 2012-04-17 2012-04-17 一种降低旋翼桨叶振动和动应力水平的设计方法 Expired - Fee Related CN102663195B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210113017.9A CN102663195B (zh) 2012-04-17 2012-04-17 一种降低旋翼桨叶振动和动应力水平的设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210113017.9A CN102663195B (zh) 2012-04-17 2012-04-17 一种降低旋翼桨叶振动和动应力水平的设计方法

Publications (2)

Publication Number Publication Date
CN102663195A true CN102663195A (zh) 2012-09-12
CN102663195B CN102663195B (zh) 2014-01-22

Family

ID=46772686

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210113017.9A Expired - Fee Related CN102663195B (zh) 2012-04-17 2012-04-17 一种降低旋翼桨叶振动和动应力水平的设计方法

Country Status (1)

Country Link
CN (1) CN102663195B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104063549A (zh) * 2014-07-02 2014-09-24 湖南大学 一种基于振动性能的动力与传动系统匹配方法
CN106326501A (zh) * 2015-06-15 2017-01-11 上海东浩兰生国际服务贸易(集团)有限公司 建筑物结构动力分析用自振频率和振型的计算方法
CN111523178A (zh) * 2020-04-21 2020-08-11 北京航空航天大学 一种降低复合材料旋翼桨毂振动载荷的方法
CN111950076A (zh) * 2020-07-10 2020-11-17 北京航空航天大学 一种降低复合材料旋翼动应力水平的设计方法
CN112052528A (zh) * 2020-09-25 2020-12-08 中国直升机设计研究所 一种直升机新型旋翼桨叶气动外形设计方法
CN112867671A (zh) * 2020-04-28 2021-05-28 深圳市大疆创新科技有限公司 振动模态优化方法、振动模态优化装置和无人机

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008524075A (ja) * 2004-12-20 2008-07-10 シコルスキー エアクラフト コーポレイション 航空機機体の吸音システム
US20090121081A1 (en) * 2007-10-12 2009-05-14 Abe Karem Composite Bulkhead and Skin Construction
CN101723091A (zh) * 2009-12-16 2010-06-09 李游 旋翼直升机的旋翼变距控制装置
WO2010134923A1 (en) * 2009-05-22 2010-11-25 Bell Helicopter Textron Inc. Rotor blade spacing for vibration attenuation

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008524075A (ja) * 2004-12-20 2008-07-10 シコルスキー エアクラフト コーポレイション 航空機機体の吸音システム
US20090121081A1 (en) * 2007-10-12 2009-05-14 Abe Karem Composite Bulkhead and Skin Construction
WO2010134923A1 (en) * 2009-05-22 2010-11-25 Bell Helicopter Textron Inc. Rotor blade spacing for vibration attenuation
CN101723091A (zh) * 2009-12-16 2010-06-09 李游 旋翼直升机的旋翼变距控制装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
GUO JUN-XIAN ET.AL.: "Composite Rotor Blade Design Optimization for Vibration Reduction with Aeroelastic Constraints", 《CHINESE JOURNAL OF AERONAUTICS》, vol. 17, no. 3, 31 August 2004 (2004-08-31), pages 152 - 158 *
向锦武等: "直升机减振的旋翼桨叶优化设计研究综述", 《北京航空航天大学学报》, vol. 27, no. 1, 8 February 2001 (2001-02-08), pages 32 - 35 *
徐敏等: "降低直升机振动水平的桨叶襟翼主动控制技术综述", 《直升机技术》, no. 2, 31 December 2007 (2007-12-31), pages 56 - 61 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104063549A (zh) * 2014-07-02 2014-09-24 湖南大学 一种基于振动性能的动力与传动系统匹配方法
CN104063549B (zh) * 2014-07-02 2017-05-10 湖南大学 一种基于振动性能的动力与传动系统匹配方法
CN106326501A (zh) * 2015-06-15 2017-01-11 上海东浩兰生国际服务贸易(集团)有限公司 建筑物结构动力分析用自振频率和振型的计算方法
CN106326501B (zh) * 2015-06-15 2019-06-25 上海东浩兰生国际服务贸易(集团)有限公司 建筑物结构动力分析用自振频率和振型的计算方法
CN111523178A (zh) * 2020-04-21 2020-08-11 北京航空航天大学 一种降低复合材料旋翼桨毂振动载荷的方法
CN112867671A (zh) * 2020-04-28 2021-05-28 深圳市大疆创新科技有限公司 振动模态优化方法、振动模态优化装置和无人机
CN111950076A (zh) * 2020-07-10 2020-11-17 北京航空航天大学 一种降低复合材料旋翼动应力水平的设计方法
CN112052528A (zh) * 2020-09-25 2020-12-08 中国直升机设计研究所 一种直升机新型旋翼桨叶气动外形设计方法
CN112052528B (zh) * 2020-09-25 2022-04-29 中国直升机设计研究所 一种直升机新型旋翼桨叶气动外形设计方法

Also Published As

Publication number Publication date
CN102663195B (zh) 2014-01-22

Similar Documents

Publication Publication Date Title
CN102663195B (zh) 一种降低旋翼桨叶振动和动应力水平的设计方法
Gaertner et al. Definition of the IEA 15-megawatt offshore reference wind turbine
Gaertner et al. IEA wind TCP task 37: definition of the IEA 15-megawatt offshore reference wind turbine
Liu et al. Effects of aerodynamic damping on the tower load of offshore horizontal axis wind turbines
Bak et al. Light Rotor: The 10-MW reference wind turbine
CN102722606B (zh) 一种降低直升机旋翼桨毂振动载荷的方法
Yang et al. Investigation on aerodynamics and active flow control of a vertical axis wind turbine with flapped airfoil
Worasinchai et al. The physics of H-Darrieus turbine starting behavior
Resor et al. Aeroelastic instability of very large wind turbine blades
CN109933837A (zh) 一种风电机组强度快速评估方法和系统
Pourrajabian et al. Hollow blades for small wind turbines operating at high altitudes
Soykasap et al. Performance enhancement of a composite tilt-rotor using aeroelastic tailoring
Dong et al. Comparative analysis of different criteria for the prediction of vortex ring state of floating offshore wind turbines
Meister et al. Simulation of a 5MW wind turbine in an atmospheric boundary layer
CN111950076B (zh) 一种降低复合材料旋翼动应力水平的设计方法
Zhu et al. Static and dynamic characteristics study of wind turbine blade
Young et al. Optimization of the skin thickness distribution in the composite wind turbine blade
Choi et al. Fluid–Structure Interaction Analysis of a Fabric Skin for Fabric-Covered Wind Turbine Blades
Malim et al. UNSTEADY AEROELASTIC RESPONSE SIMULATION OF A 3D PRINTED HIGH ALTITUDE PROPELLER USING FLUID-STRUCTURE INTERACTION
Quintana et al. Nonlinear characterization of the piecewise structural effects on whirl flutter of a rotor-nacelle system
Tiwari Dynamics and Analysis of Rectangular Rotor Blades
Cheng et al. Fluid-structure coupled analyses of composite wind turbine blades
Higgins et al. A Propeller Blade Design for Experimental Stall Flutter Investigations
Roscher et al. Combined structural optimization and aeroelastic analysis of a Vertical Axis Wind Turbine
BONETTI Whirl flutter stability enhancement by rotor twist optimization

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent of invention or patent application
CB03 Change of inventor or designer information

Inventor after: Xiang Jinwu

Inventor after: Ren Yiru

Inventor after: Luo Zhangping

Inventor after: Zhang Lianhong

Inventor after: Guo Junxian

Inventor after: Zhang Yajun

Inventor before: Xiang Jinwu

Inventor before: Ren Yiru

Inventor before: Luo Zhangping

Inventor before: Guo Junxian

Inventor before: Zhang Yajun

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: XIANG JINWU REN YIRU LUO ZHANGPING GUO JUNXIAN ZHANG YAJUN TO: XIANG JINWUREN YIRU LUO ZHANGPING ZHANG LIANHONG GUO JUNXIAN ZHANG YAJUN

C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140122

Termination date: 20150417

EXPY Termination of patent right or utility model