CN101571888B - 一种继电器用多有限元软件联合仿真分析方法 - Google Patents
一种继电器用多有限元软件联合仿真分析方法 Download PDFInfo
- Publication number
- CN101571888B CN101571888B CN2009100723118A CN200910072311A CN101571888B CN 101571888 B CN101571888 B CN 101571888B CN 2009100723118 A CN2009100723118 A CN 2009100723118A CN 200910072311 A CN200910072311 A CN 200910072311A CN 101571888 B CN101571888 B CN 101571888B
- Authority
- CN
- China
- Prior art keywords
- finite element
- relay
- simulation
- parameter
- 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.)
- Expired - Fee Related
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 34
- 238000004458 analytical method Methods 0.000 title claims abstract description 20
- 238000013461 design Methods 0.000 claims abstract description 11
- 238000013519 translation Methods 0.000 claims abstract description 6
- 238000004364 calculation method Methods 0.000 claims abstract description 4
- 101150041750 Surf4 gene Proteins 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 3
- 238000009434 installation Methods 0.000 abstract 1
- 230000004907 flux Effects 0.000 description 14
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000005672 electromagnetic field Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000001052 transient effect Effects 0.000 description 3
- 238000012938 design process Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
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 00
转换成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和BuildVolume可以建立模型的面和体单元,从而完成继电器的实体建模。接着,采用assignMeshPoint命令控制分网的精细程度,以及Material和assignRegionToVolumes命令定义模型材料,Flux即可按照建模脚本对继电器进行建模、分网和材料属性设置。
步骤三:Flux仿真模块根据用户输入的衔铁转角参数和通过系统注册表得到的Flux安装状况,生成仿真脚本文件F3D_INI.SPI;
步骤四:Flux仿真模块再根据线圈电流信息控制Flux按照仿真脚本文件循环仿真出各种电流情况下的吸力特性;
步骤五:Flux仿真模块生成仿真结果文件CHARACTERISTIC-ANALYSIS.PRT,通过对结果文件进行读取得到需要的继电器吸力特性。
Claims (1)
1.一种继电器用多有限元软件联合仿真分析方法,包括继电器描述语言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为面拉伸方向的矢量控制点;
模型描述部分还包括将几何坐标位置信息转换为工程图纸参数信息的转换语句,描述格式如下:
SetMap C Value K S1 S2 S3......SetMapEnd
其中,C为工程图纸中的参数名称,Value为这个参数的原始值,K为相关系数,S1、S2、S3为参数的编号,即满足式(1)的关系;
Sn=Sn0+K(C0-C1) (1)
其中,Sn为描述语言参数的当前值,Sn0为描述语言参数的原始值,C0为模型参数的原始值,C1为模型参数的当前值;
有限元仿真软件的联合调用,具体步骤如下:
步骤一:采用RML语言对继电器进行建模,将其工程图纸中的尺寸信息转化为RML语言描述的脚本文件;
步骤二:针对有限元仿真软件开发出相应的翻译模块,将RML语言描述的继电器模型转化为对应有限元仿真软件的建模脚本;
步骤三:结合有限元软件安装情况、仿真分析类型和仿真参数生成可以控制仿真软件工作的脚本文件;
步骤四:根据继电器设计中不同的分析需要控制有限元软件进行仿真计算;
步骤五:有限元仿真软件生成仿真结果文件,对结果文件进行读取得到需要的继电器特性信息。
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 CN101571888A (zh) | 2009-11-04 |
CN101571888B true 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 (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101763447B (zh) * | 2009-12-28 | 2012-01-11 | 中国农业大学 | 三维造型软件与有限元分析软件之间数据自动转换的方法 |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102012955B (zh) * | 2010-11-29 | 2014-11-05 | 中兴通讯股份有限公司 | 获取背板连接器针头钻阵列的散射矩阵的方法和装置 |
CN103309726B (zh) * | 2013-05-23 | 2016-07-06 | 南京南瑞继保工程技术有限公司 | 二次装置中液晶部件的仿真实现方法 |
CN104990706A (zh) * | 2015-06-26 | 2015-10-21 | 江苏润源控股集团有限公司 | 基于adams脚本的复杂机构实际输出动程分析方法 |
CN105893701B (zh) * | 2016-04-25 | 2019-01-18 | 哈尔滨工业大学 | 一种基于温度补偿矩阵的电磁机构动态特性的计算方法 |
CN106202794A (zh) * | 2016-07-21 | 2016-12-07 | 哈尔滨工业大学 | 基于多场有限元仿真的继电器类单机输出特性获取方法 |
CN108536882A (zh) * | 2017-03-06 | 2018-09-14 | 无锡飞谱电子信息技术有限公司 | 一种高频电磁场有限元仿真系统 |
CN107748826B (zh) * | 2017-11-08 | 2018-09-25 | 哈尔滨工业大学 | 一种继电器耐力学性能贮存退化分析方法 |
CN109858123B (zh) * | 2019-01-22 | 2022-09-27 | 安徽大学 | 基于Python和Maxwell的永磁球形电机三维电磁转矩自动化分析方法 |
CN111782289B (zh) * | 2020-05-20 | 2023-09-29 | 北京仿真中心 | 基于Adams软件的协同方法、系统、计算机设备及存储介质 |
CN111695258B (zh) * | 2020-06-11 | 2023-05-26 | 哈尔滨工业大学 | 一种电磁继电器动态特性仿真动能注入仿真方法 |
Citations (4)
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 | 哈尔滨工业大学 | 电磁继电器静态特性测试分析装置 |
CN1976154A (zh) * | 2006-12-22 | 2007-06-06 | 北京四方继保自动化股份有限公司 | 采用脚本描述实现的继电器仿真模型 |
CN101377670A (zh) * | 2008-10-07 | 2009-03-04 | 深圳创维-Rgb电子有限公司 | 一种插件程序的转换方法、系统 |
-
2009
- 2009-06-17 CN CN2009100723118A patent/CN101571888B/zh not_active Expired - Fee Related
Patent Citations (4)
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 | 哈尔滨工业大学 | 电磁继电器静态特性测试分析装置 |
CN1976154A (zh) * | 2006-12-22 | 2007-06-06 | 北京四方继保自动化股份有限公司 | 采用脚本描述实现的继电器仿真模型 |
CN101377670A (zh) * | 2008-10-07 | 2009-03-04 | 深圳创维-Rgb电子有限公司 | 一种插件程序的转换方法、系统 |
Non-Patent Citations (2)
Title |
---|
肖唐杰,杨文英,翟国富.航天电磁继电器综合仿真分析系统的设计.机电元件.2009,29(2),17-23. * |
靳现林,彭慧滔,王毅.采用脚本描述的继电器仿真模型及其软件实现.2008中国电力系统保护与控制学术研讨会论文集.2008,624-629. * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101763447B (zh) * | 2009-12-28 | 2012-01-11 | 中国农业大学 | 三维造型软件与有限元分析软件之间数据自动转换的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101571888A (zh) | 2009-11-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101571888B (zh) | 一种继电器用多有限元软件联合仿真分析方法 | |
Bézivin et al. | Bridging the ms/dsl tools and the eclipse modeling framework | |
CN100435049C (zh) | 一种面向生产现场的半沉浸式装配工艺规划方法 | |
CN105260177A (zh) | 基于SiPESC平台的Python扩展模块开发方法 | |
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 | |
Gao et al. | Research on the Application of Path Planning in Augmented Reality-Assisted Assembly | |
Rantzau et al. | The integration of immersive Virtual Reality applications into Catia V5 | |
Lee et al. | A study on iNterface and Interaction using Chatgpt System in Virtual Reality Space | |
Roßkopf et al. | Aspects of flight control software—a software engineering point of view | |
Guiling et al. | Virtual human control technology for immersed virtual maintenance system in DELMIA environment | |
Mannuß et al. | From scene graph centered to entity centered virtual environments | |
KR20080011500A (ko) | 스스로 추론하고 반응하며 행동하는 3차원 가상객체의설계방법과 구동 시스템 | |
Xiang et al. | Multi-disciplinary modeling and collaborative simulation of multi-robot systems based on HLA | |
Jin-Liang et al. | Parametric Design of Visual Simulation Based on Vega Prime | |
Figueroa et al. | A framework for 3D interaction techniques | |
Su et al. | Research on Satellite Simulation for Mobile Terminals | |
Kirov | An Object-Oriented HLA Simulation Study | |
Qin et al. | Research on the Application of 3D Modeling Technology in Mechanical Structure Teaching | |
Janyška et al. | Hydrodynamical Picture of QM | |
Hernández-Ortiz et al. | Constraining the position of the CEP through the Speed of Sound in the LSMq | |
Gao et al. | Study on the Development and Application of osgART Based on MFC | |
Zenor et al. | Software support tools for high-speed real-time simulations | |
CN118194544A (zh) | 一种采用双系统进行磁悬物建模的方法及仿真装置 |
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 |