CN101571888A - 一种继电器用多有限元软件联合仿真分析方法 - Google Patents

一种继电器用多有限元软件联合仿真分析方法 Download PDF

Info

Publication number
CN101571888A
CN101571888A CNA2009100723118A CN200910072311A CN101571888A CN 101571888 A CN101571888 A CN 101571888A CN A2009100723118 A CNA2009100723118 A CN A2009100723118A CN 200910072311 A CN200910072311 A CN 200910072311A CN 101571888 A CN101571888 A CN 101571888A
Authority
CN
China
Prior art keywords
finite element
relay
parameter
simulation
software
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
CNA2009100723118A
Other languages
English (en)
Other versions
CN101571888B (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN2009100723118A priority Critical patent/CN101571888B/zh
Publication of CN101571888A publication Critical patent/CN101571888A/zh
Application granted granted Critical
Publication of CN101571888B publication Critical patent/CN101571888B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Devices For Executing Special Programs (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明的目的在于提供一种用于继电器设计的综合仿真分析系统,实现多有限元软件联合仿真的继电器用多有限元软件联合仿真分析方法。包括继电器描述语言和有限元仿真软件的联合调用,继电器描述语言分为两大部分:参数部分和模型描述部分;有限元仿真软件的联合调用,采用RML语言对继电器进行建模;针对有限元仿真软件开发出相应的翻译模块;结合有限元软件安装情况、仿真分析类型和仿真参数生成可以控制仿真软件工作的脚本文件;根据继电器设计中不同的分析需要控制有限元软件进行仿真计算;本发明提出了一种用于航天电磁继电器设计的建模语言,通过对Ansys、Flux和Adams等有限元软件进行二次开发,从而实现了多有限元软件的联合仿真。

Description

一种继电器用多有限元软件联合仿真分析方法
(一)技术领域
本发明涉及继电器技术,具体说就是一种继电器用多有限元软件联合仿真分析方法。
(二)背景技术
航天电磁继电器是用于导弹、运载火箭、人造卫星、宇宙飞船、航天飞机及其配套地面测控设备中完成信号传递、执行控制、系统配电等功能的继电器,是国防电子系统中主要电子元器件之一。
电磁继电器由电磁系统、触点系统和释放弹簧等组成,通过控制线圈电流所产生的电磁吸力驱动衔铁而实现触点开闭或转换功能。因此,针对电磁继电器的设计过程,涉及电磁场求解、机械动力学、热学等多个领域。而且由于磁系统的非线性,使得这一过程的定量计算和分析变得十分复杂。
目前,针对机、电、磁、热多场域的研究,多采用数值分析方法,其中以有限元方法应用最为广泛。而且,随着计算机技术的不断发展,不同领域的有限元分析软件也都日益成熟,如Ansys、Adams、Flux、Patran/Nastran等。虽然每种软件各有所长,但又都有一定的局限性,如Ansys能对静态的电磁问题进行求解,但对于瞬态电磁问题的求解,则需要用户自己编程实现;而Flux可以求解瞬态电磁场问题,但其求解动力学的模型相对简单,远不如Adams在处理运动学问题上得心应手。而在进行航天电磁继电器的设计过程中,其动态特性的分析就既要对电磁系统的瞬态电磁场问题进行研究,同时还要考虑接触系统的运动状态。对于继电器设计人员而言,每采用一种分析软件,就要重新进行一次建模,十分不方便,更不利于产品的设计和开发。
(三)发明内容
本发明的目的在于提供一种用于继电器设计的综合仿真分析系统,实现多有限元软件联合仿真的继电器用多有限元软件联合仿真分析方法。
本发明的目的是这样实现的:本发明一种继电器用多有限元软件联合仿真分析方法,包括继电器描述语言RML和有限元仿真软件的联合调用,继电器描述语言RML分为两大部分:参数部分和模型描述部分;参数部分主要包含变量的定义,包括变量名和变量值;模型描述部分包含采用规定语法对继电器模型进行的描述;在语言中参数部分和模型描述部分分别采用Variable和Geometry两个关键字进行标识;在参数部分,参数表达格式如下:
参数名参数值
在模型描述部分,按照点、线、面、体层次定义了一系列的语句进行描述;
首先定义的是坐标系,其格式如下:
coord ID x1 y1 z1 x2 y2 z2 ID2
其中,ID为建立坐标系编号,x1、y1、z1分别是新坐标系原点在直角坐标系下三个轴的数值,x2、y2、z2分别是新坐标系相对于原坐标系在三个轴向的旋转值;ID2为原坐标系的编号;
点的定义为:
point ID x y z
其中,ID描述的点的编号,x、y、z为点在坐标中的三轴向数值;
线的定义为:
lines ID P1 P2
其中,ID为线的编号,P1、P2为线段两端点的编号;
面的定义为:
surf4 ID L1 L2 L3 L4
其中,ID为面的编号,L1、L2、L3、L4为组成面的四条线的编号;
体的定义为:
solidz ID S1 P1 P2
其中,ID为体的编号,S1为基准面,P1、P2为面拉伸方向的矢量控制点;
通过这些描述定义,可以描述继电器的所有几何信息,但考虑到实际工程中采用的工程图以长、宽、角度等信息作为参数的表达方法,与RML中采用的几何坐标为参数的表达方法不同;因此,定义了将这些点的几何坐标位置信息,转换为工程图纸中的参数信息,描述格式如下:
SetMap C Value K S1 S2 S3……SetMapEnd
其中,C为工程图纸中的参数名称,Value为这个参数的原始值,K为相关系数,S1、S2、S3为参数的编号,即满足式(1)的关系;
Sn=Sn0+K(C0-C)                (1)
其中,Sn为描述语言参数的当前值,Sn0为描述语言参数的原始值,
C0为模型参数的原始值,C为模型参数的当前值;
有限元仿真软件的联合调用,具体步骤如下:
步骤一:采用RML语言对继电器进行建模,将其工程图纸中的尺寸信息转化为RML语言描述的脚本文件;
步骤二:针对有限元仿真软件开发出相应的翻译模块,将RML语言描述的继电器模型转化为对应有限元仿真软件的建模脚本;
步骤三:结合有限元软件安装情况、仿真分析类型和仿真参数生成可以控制仿真软件工作的脚本文件;
步骤四:根据继电器设计中不同的分析需要控制有限元软件进行仿真计算;
步骤五:有限元仿真软件生成仿真结果文件,对结果文件进行读取。
本发明一种继电器用多有限元软件联合仿真分析方法,提出了一种用于航天电磁继电器设计的建模语言——Relay Modeling Language(RML),并通过对Ansys、Flux和Adams等有限元软件进行二次开发,使该语言所建的继电器模型可以被上述任一有限元软件识别。在此基础上,采用C++Builder开发工具开发了用于继电器设计的综合仿真分析系统,从而实现了多有限元软件的联合仿真。
(四)附图说明
图1为本发明的多有限元软件调用方框图;
图2为本发明的有限元仿真调用流程图;
图3为Flux有限元仿真软件调用流程图。
(五)具体实施方式
下面结合附图举例对本发明作进一步说明。
实施例1:结合图1-图3,本发明一种继电器用多有限元软件联合仿真分析方法,包括以下内容:
1.RML建模语言
要使Ansys、Flux和Adams几个有限元仿真软件可以配合仿真,首要的是它们必须有共同的仿真模型来源,但是这几个有限元仿真软件的模型格式和脚本语言都不相同,而且它们之间不易转换。因此需要一个能够描述继电器模型的通用而且可以比较容易转换成需要的有限元仿真脚本语言。本发明制定了一种继电器描述语言(RML)实现这一功能。因为继电器描述语言需要转化为有限元仿真脚本,所以继电器描述语言参考有限元仿真脚本进行开发。同时参考计算机图形学的基本理论,按点、线、面、体的层次对继电器进行描述。
为了使继电器模型便于调整参数,为继电器设计提供极大的方便,继电器描述语言分为两大部分:参数部分和模型描述部分。参数部分主要包含变量的定义,包括变量名和变量值。模型描述部分包含采用规定语法对继电器模型进行的描述。在语言中参数部分和模型描述部分分别采用Variable和Geometry两个关键字进行标识。在参数部分,参数表达格式如下:
参数名参数值
在模型描述部分,按照点、线、面、体层次定义了一系列的语句进行描述。
首先定义的是坐标系,其格式如下:
coord ID x1 y1 z1 x2 y2 z2 ID2
其中,ID为建立坐标系编号,x1、y1、z1分别是新坐标系原点在直角坐标系下三个轴的数值,x2、y2、z2分别是新坐标系相对于原坐标系在三个轴向的旋转值。ID2为原坐标系的编号。
点的定义为:
point ID x y z
其中,ID描述的点的编号,x、y、z为点在坐标中的三轴向数值。
线的定义为:
lines ID P1 P2
其中,ID为线的编号,P1、P2为线段两端点的编号。
面的定义为:
surf4 ID L1 L2 L3 L4
其中,ID为面的编号,L1、L2、L3、L4为组成面的四条线的编号。
体的定义为:
solidz ID S1 P1 P2
其中,ID为体的编号,S1为基准面,P1、P2为面拉伸方向的矢量控制点。
通个这些描述定义,可以描述继电器的所有几何信息,但考虑到实际工程中采用的工程图以长、宽、角度等信息作为参数的表达方法,与RML中采用的几何坐标为参数的表达方法不同。因此,定义了将这些点的几何坐标位置信息,转换为工程图纸中的参数信息,描述格式如下:
SetMap C Value K S1 S2 S3……SetMapEnd
其中,C为工程图纸中的参数名称,Value为这个参数的原始值,K为相关系数,S1、S2、S3为参数的编号,即满足式(1)的关系。
Sn=Sn0+K(C0-C)                    (1)
其中,Sn为描述语言参数的当前值,Sn0为描述语言参数的原始值,C0为模型参数的原始值,C为模型参数的当前值。
2.有限元仿真软件的联合调用
具体步骤如下:
步骤一:采用RML语言对继电器进行建模,将其工程图纸中的尺寸信息转化为RML语言描述的脚本文件;
步骤二:针对有限元仿真软件开发出相应的翻译模块,将RML语言描述的继电器模型转化为对应有限元仿真软件的建模脚本;
步骤三:结合有限元软件安装情况、仿真分析类型和仿真参数生成可以控制仿真软件工作的脚本文件;
步骤四:根据继电器设计中不同的分析需要控制有限元软件进行仿真计算;
步骤五:有限元仿真软件生成仿真结果文件,对结果文件进行读取得到需要的继电器特性信息。
实施例2:结合图1、图2和图3,本发明是一种用于继电器设计的多有限元软件联合仿真方法,具体有限元软件调用仿真实例如下:
步骤一:采用RML语言对继电器进行建模,将其工程图纸中的尺寸信息转化为RML语言描述的脚本文件;
步骤二:Flux调用首先需要开发Flux的翻译模块,将RML模型文件转换为Flux可以使用的PYTHON语言的建模脚本。
首先是参数的转换,例如在RML语言中的参数
S5002-2.500E-03
其中,S5002是参数名,-2.500E-03是参数值
转换成PYTHON语言后是
ParameterGeom(name=′S5002′,expression=′-2.500E-03′)
接下来是坐标系的转换。
RML语言中坐标系为:
coord 5001 S5001 S5002 S5003 S5220 S5220 S5220 0 0
转换成PYTHON语言是
CoordSysCartesian(name=′5001′,parentCoordSys=Local(coordSys=CoordSys[′XYZ1′]),origin=[′S5001′,′S5002′,′S5003′],rotationAngles=RotationAngles(angleX=′0′,angleY=′0′,angleZ=′S5220′))
点为:
point 5001 S5004 S5005 S5006 5001
转化为:
PointCoordinates(color=Color[′White′],visibility=Visibility[′VISIBLE′],coordSys=CoordSys[′5001′],uvw=[′S5004′,′S5005′,′S5006′],nature=Nature[′STANDARD′])
线为:
lines 5001 5003 5007
转化为:
LineSegment(color=Color[′White′],visibility=Visibility[′VISIBLE′],defPoint=[Point[3],Point[7]],nature=Nature[′STANDARD′])
到此RML语言描述的继电器模型就可以全部转换为PYTHON语言的建模脚本。然后使用PYTHON语言中BuildFace和Build Volume可以建立模型的面和体单元,从而完成继电器的实体建模。接着,采用assignMeshPoint命令控制分网的精细程度,以及Material和assignRegionToVolumes命令定义模型材料,Flux即可按照建模脚本对继电器进行建模、分网和材料属性设置。
步骤三:Flux仿真模块根据用户输入的衔铁转角参数和通过系统注册表得到的Flux安装状况,生成仿真脚本文件F3D_INI.SPI;
步骤四:Flux仿真模块再根据线圈电流信息控制Flux按照仿真脚本文件循环仿真出各种电流情况下的吸力特性;
步骤五:Flux仿真模块生成仿真结果文件CHARACTERISTIC-ANALYSIS.PRT,通过对结果文件进行读取得到需要的继电器吸力特性。

Claims (1)

1.一种继电器用多有限元软件联合仿真分析方法,包括继电器描述语言RML和有限元仿真软件的联合调用,其特征在于:继电器描述语言RML分为两大部分:参数部分和模型描述部分;参数部分主要包含变量的定义,包括变量名和变量值;模型描述部分包含采用规定语法对继电器模型进行的描述;在语言中参数部分和模型描述部分分别采用Variable和Geometry两个关键字进行标识;在参数部分,参数表达格式如下:
参数名参数值
在模型描述部分,按照点、线、面、体层次定义了一系列的语句进行描述;
首先定义的是坐标系,其格式如下:
coord ID x1 y1 z1 x2 y2 z2 ID2
其中,ID为建立坐标系编号,x1、y1、z1分别是新坐标系原点在直角坐标系下三个轴的数值,x2、y2、z2分别是新坐标系相对于原坐标系在三个轴向的旋转值;ID2为原坐标系的编号;
点的定义为:
point ID x y z
其中,ID描述的点的编号,x、y、z为点在坐标中的三轴向数值;线的定义为:
lines ID P1 P2
其中,ID为线的编号,P1、P2为线段两端点的编号;
面的定义为:
surf4 ID L1 L2 L3 L4
其中,ID为面的编号,L1、L2、L3、L4为组成面的四条线的编号;
体的定义为:
solidz ID S1 P1 P2
其中,ID为体的编号,S1为基准面,P1、P2为面拉伸方向的矢量控制点;
通过这些描述定义,可以描述继电器的所有几何信息,但考虑到实际工程中采用的工程图以长、宽、角度等信息作为参数的表达方法,与RML中采用的几何坐标为参数的表达方法不同;因此,定义了将这些点的几何坐标位置信息,转换为工程图纸中的参数信息,描述格式如下:
SetMap C Value K S1 S2 S3……SetMapEnd
其中,C为工程图纸中的参数名称,Value为这个参数的原始值,K为相关系数,S1、S2、S3为参数的编号,即满足式(1)的关系;
Sn=Sn0+K(C0-C)                            (1)
其中,Sn为描述语言参数的当前值,Sn0为描述语言参数的原始值,
C0为模型参数的原始值,C为模型参数的当前值;
有限元仿真软件的联合调用,具体步骤如下:
步骤一:采用RML语言对继电器进行建模,将其工程图纸中的尺寸信息转化为RML语言描述的脚本文件;
步骤二:针对有限元仿真软件开发出相应的翻译模块,将RML语言描述的继电器模型转化为对应有限元仿真软件的建模脚本;
步骤三:结合有限元软件安装情况、仿真分析类型和仿真参数生成可以控制仿真软件工作的脚本文件;
步骤四:根据继电器设计中不同的分析需要控制有限元软件进行仿真计算;
步骤五:有限元仿真软件生成仿真结果文件,对结果文件进行读取得到需要的继电器特性信息。
CN2009100723118A 2009-06-17 2009-06-17 一种继电器用多有限元软件联合仿真分析方法 Expired - Fee Related CN101571888B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100723118A CN101571888B (zh) 2009-06-17 2009-06-17 一种继电器用多有限元软件联合仿真分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100723118A CN101571888B (zh) 2009-06-17 2009-06-17 一种继电器用多有限元软件联合仿真分析方法

Publications (2)

Publication Number Publication Date
CN101571888A true CN101571888A (zh) 2009-11-04
CN101571888B CN101571888B (zh) 2010-09-08

Family

ID=41231244

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100723118A Expired - Fee Related CN101571888B (zh) 2009-06-17 2009-06-17 一种继电器用多有限元软件联合仿真分析方法

Country Status (1)

Country Link
CN (1) CN101571888B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102012955A (zh) * 2010-11-29 2011-04-13 中兴通讯股份有限公司 获取背板连接器针头钻阵列的散射矩阵的方法和装置
CN103309726A (zh) * 2013-05-23 2013-09-18 南京南瑞继保工程技术有限公司 二次装置中液晶部件的仿真实现方法
CN104990706A (zh) * 2015-06-26 2015-10-21 江苏润源控股集团有限公司 基于adams脚本的复杂机构实际输出动程分析方法
CN105893701A (zh) * 2016-04-25 2016-08-24 哈尔滨工业大学 一种基于温度补偿矩阵的电磁机构动态特性的计算方法
CN106202794A (zh) * 2016-07-21 2016-12-07 哈尔滨工业大学 基于多场有限元仿真的继电器类单机输出特性获取方法
CN107748826A (zh) * 2017-11-08 2018-03-02 哈尔滨工业大学 一种继电器耐力学性能贮存退化分析方法
CN108536882A (zh) * 2017-03-06 2018-09-14 无锡飞谱电子信息技术有限公司 一种高频电磁场有限元仿真系统
CN109858123A (zh) * 2019-01-22 2019-06-07 安徽大学 基于Python和Maxwell的永磁球形电机三维电磁转矩自动化分析方法
CN111695258A (zh) * 2020-06-11 2020-09-22 哈尔滨工业大学 一种电磁继电器动态特性仿真动能注入仿真方法
CN111782289A (zh) * 2020-05-20 2020-10-16 北京仿真中心 基于Adams软件的协同方法、系统、计算机设备及存储介质

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101763447B (zh) * 2009-12-28 2012-01-11 中国农业大学 三维造型软件与有限元分析软件之间数据自动转换的方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6236956B1 (en) * 1996-02-16 2001-05-22 Avant! Corporation Component-based analog and mixed-signal simulation model development including newton step manager
CN1609625A (zh) * 2004-11-12 2005-04-27 哈尔滨工业大学 电磁继电器静态特性测试分析装置
CN100470981C (zh) * 2006-12-22 2009-03-18 北京四方继保自动化股份有限公司 基于采用脚本描述的继电器仿真模型的实现方法
CN101377670B (zh) * 2008-10-07 2010-09-01 深圳创维-Rgb电子有限公司 一种插件程序的转换方法、系统

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102012955A (zh) * 2010-11-29 2011-04-13 中兴通讯股份有限公司 获取背板连接器针头钻阵列的散射矩阵的方法和装置
CN102012955B (zh) * 2010-11-29 2014-11-05 中兴通讯股份有限公司 获取背板连接器针头钻阵列的散射矩阵的方法和装置
CN103309726A (zh) * 2013-05-23 2013-09-18 南京南瑞继保工程技术有限公司 二次装置中液晶部件的仿真实现方法
CN103309726B (zh) * 2013-05-23 2016-07-06 南京南瑞继保工程技术有限公司 二次装置中液晶部件的仿真实现方法
CN104990706A (zh) * 2015-06-26 2015-10-21 江苏润源控股集团有限公司 基于adams脚本的复杂机构实际输出动程分析方法
CN105893701A (zh) * 2016-04-25 2016-08-24 哈尔滨工业大学 一种基于温度补偿矩阵的电磁机构动态特性的计算方法
CN105893701B (zh) * 2016-04-25 2019-01-18 哈尔滨工业大学 一种基于温度补偿矩阵的电磁机构动态特性的计算方法
CN106202794A (zh) * 2016-07-21 2016-12-07 哈尔滨工业大学 基于多场有限元仿真的继电器类单机输出特性获取方法
CN108536882A (zh) * 2017-03-06 2018-09-14 无锡飞谱电子信息技术有限公司 一种高频电磁场有限元仿真系统
CN107748826A (zh) * 2017-11-08 2018-03-02 哈尔滨工业大学 一种继电器耐力学性能贮存退化分析方法
CN109858123A (zh) * 2019-01-22 2019-06-07 安徽大学 基于Python和Maxwell的永磁球形电机三维电磁转矩自动化分析方法
CN109858123B (zh) * 2019-01-22 2022-09-27 安徽大学 基于Python和Maxwell的永磁球形电机三维电磁转矩自动化分析方法
CN111782289A (zh) * 2020-05-20 2020-10-16 北京仿真中心 基于Adams软件的协同方法、系统、计算机设备及存储介质
CN111782289B (zh) * 2020-05-20 2023-09-29 北京仿真中心 基于Adams软件的协同方法、系统、计算机设备及存储介质
CN111695258A (zh) * 2020-06-11 2020-09-22 哈尔滨工业大学 一种电磁继电器动态特性仿真动能注入仿真方法
CN111695258B (zh) * 2020-06-11 2023-05-26 哈尔滨工业大学 一种电磁继电器动态特性仿真动能注入仿真方法

Also Published As

Publication number Publication date
CN101571888B (zh) 2010-09-08

Similar Documents

Publication Publication Date Title
CN101571888B (zh) 一种继电器用多有限元软件联合仿真分析方法
Ivanenko et al. The gauge treatment of gravity
Hees et al. Combined cosmological and solar system constraints on chameleon mechanism
Bézivin et al. Bridging the ms/dsl tools and the eclipse modeling framework
CN100435049C (zh) 一种面向生产现场的半沉浸式装配工艺规划方法
Ibragimov et al. A simple method for group analysis and its application to a model of detonation
Xiang et al. Agent-based composable simulation for virtual prototyping of fluid power system
Canfora et al. General relativity from Einstein-Gauss-Bonnet gravity
Krösche et al. The xPOI concept
Suk et al. A study on the effect of the use of augmented reality on students' quantity take-off performance
Zhu et al. Design and realization of virtual classroom
Gao et al. Research on the Application of Path Planning in Augmented Reality-Assisted Assembly
Lee et al. A study on iNterface and Interaction using Chatgpt System in Virtual Reality Space
Rantzau et al. The integration of immersive Virtual Reality applications into Catia V5
Guiling et al. Virtual human control technology for immersed virtual maintenance system in DELMIA environment
KR20080011500A (ko) 스스로 추론하고 반응하며 행동하는 3차원 가상객체의설계방법과 구동 시스템
Xiang et al. Multi-disciplinary modeling and collaborative simulation of multi-robot systems based on HLA
Mannuß et al. From scene graph centered to entity centered virtual environments
Roßkopf et al. Aspects of flight control software—a software engineering point of view
Jin-Liang et al. Parametric Design of Visual Simulation Based on Vega Prime
Figueroa et al. A framework for 3D interaction techniques
Kirov An Object-Oriented HLA Simulation Study
Su et al. Research on Satellite Simulation for Mobile Terminals
CN118194544A (zh) 一种采用双系统进行磁悬物建模的方法及仿真装置
Hernández-Ortiz et al. Constraining the position of the CEP through the Speed of Sound in the LSMq

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

Granted publication date: 20100908

CF01 Termination of patent right due to non-payment of annual fee