CN102520279A - 空间电子设备加速寿命试验中温度加速基准应力确定方法 - Google Patents

空间电子设备加速寿命试验中温度加速基准应力确定方法 Download PDF

Info

Publication number
CN102520279A
CN102520279A CN2011104038146A CN201110403814A CN102520279A CN 102520279 A CN102520279 A CN 102520279A CN 2011104038146 A CN2011104038146 A CN 2011104038146A CN 201110403814 A CN201110403814 A CN 201110403814A CN 102520279 A CN102520279 A CN 102520279A
Authority
CN
China
Prior art keywords
ambient temperature
operating ambient
temperature
electronic equipment
reliability
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
CN2011104038146A
Other languages
English (en)
Other versions
CN102520279B (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.)
510 Research Institute of 5th Academy of CASC
Original Assignee
510 Research Institute of 5th Academy of CASC
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 510 Research Institute of 5th Academy of CASC filed Critical 510 Research Institute of 5th Academy of CASC
Priority to CN201110403814.6A priority Critical patent/CN102520279B/zh
Publication of CN102520279A publication Critical patent/CN102520279A/zh
Application granted granted Critical
Publication of CN102520279B publication Critical patent/CN102520279B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Tests Of Electronic Circuits (AREA)

Abstract

本发明为一种空间电子设备加速寿命试验中温度加速基准应力确定方法,属于电子产品技术领域。首先针对工作环境温度的周期性变化,将环境温度进行等步阶分割,再分别计算每个台阶温度下的对应时间段的可靠度;然后计算一个周期和整个寿命末期的可靠度,最后给出改进的可靠度指数计算公式;当选择工作环境温度作为加速应力即通过提高工作环境温度进行空间电子设备的加速寿命试验时,应采用改进的可靠度指数计算公式,在可靠性不变的前提下将周期性变化工作环境温度等效变换为恒定型工作环境温度,从而确定基准应力。本发明解决了通过提高工作环境温度进行空间电子设备的加速寿命试验时,如何准确确定温度加速基准应力这一技术问题。

Description

空间电子设备加速寿命试验中温度加速基准应力确定方法
技术领域
本发明用于航天器的空间电子设备在变化工作环境温度下加速寿命试验中温度加速基准应力的确定,属于电子产品技术领域。
背景技术
可靠性预计是为了估计产品在给定工作条件下的可靠性而进行的工作,其主要目的之一是为可靠性增长试验、鉴定试验及费用核算等方面的研究提供依据。从此目的可以看出,提供产品准确的可靠性预计结果显得至关重要。对电子产品而言,可靠性预计通常是计算产品的可靠度。
作为电子产品,空间电子设备采用应力分析法进行可靠性预计,其工作期间的应力值应当恒定才能按照可靠度指数公式计算,实际情况更多是:空间电子设备工作期间的电应力值是恒定,环境应力值在变化,随着航天器姿态、轨道位置的周期性变化,空间电子设备经受的环境应力亦呈现周期性的变化,尤以对空间电子设备在轨可靠性产生重要影响的工作环境温度的周期性变化为显著特点。在空间电子设备采用应力分析法进行可靠性预计时,目前国内外航天领域的通行方法是:电应力由于恒定而取其实际值,工作环境温度由于变化而取其最大值,亦即将工作环境温度取其最大值并视为恒温。很显然,这样得到的可靠性预计结果不准确,很保守,尤其当工作环境温度的变化范围很大时,采用此法的缺陷更是显而易见。因此,需要寻求更为科学的可靠性预计方法并基此获得更为准确的可靠性预计结果。空间电子设备的可靠性要求都很高,有的品种要求其可靠度达到R≥0.999,按照可靠度的指数函数分布规律,要达到如此高的可靠度,工作失效率将很低,因而其预期寿命将极长。对于这样的高可靠空间电子设备,如果采用在正常的工作条件下,通过寿命试验的方法来评估其可靠性,往往需要耗费很长的时间,付出昂贵的经费和很高的人员投入,是很不现实的,长期寿命试验方法已经不能适应迅速发展的空间电子设备研制的需要。因此,通过加速寿命试验,可以在相对短的时间内评估其可靠性。
由于工作环境温度是影响空间电子设备寿命以及可靠性的一个重要因素,因而通常选择工作环境温度作为加速应力亦即通过提高工作环境温度进行空间电子设备的加速寿命试验,此时必须先确定空间电子设备正常工作条件下的工作环境温度应力水平,且此工作环境温度应力应为恒定应力或等效为恒定应力,亦称基准应力。
目前的方法是将变化工作环境温度的最大值作为基准应力。很显然,这样得到的加速系数以及随之获得的可靠性评估结果都小于实际情况,还会增加不必要的试验成本,尤其当工作环境温度的变化范围很大时,采用此法的缺陷更是不容忽视。因此,如何准确确定变化工作环境温度下的温度加速基准应力,这不但关系到加速系数的准确获得以及最终的可靠性评估结果,而且关系到包括子样数量、试验时间、人员安排等的试验成本问题,是进行温度加速寿命试验前首先需要解决的关键问题。
发明内容
本发明的目的是解决通过提高工作环境温度进行空间电子设备的加速寿命试验时,如何准确确定温度加速基准应力这一技术问题。
本发明的目的是通过以下技术方案实现的。
该空间电子设备加速寿命试验中温度加速基准应力确定方法,首先针对工作环境温度的周期性变化,将环境温度进行等步阶分割,再分别计算每个台阶温度下的对应时间段的可靠度;然后计算一个周期和整个寿命末期的可靠度,最后给出改进的可靠度指数计算公式;当选择工作环境温度作为加速应力即通过提高工作环境温度进行空间电子设备的加速寿命试验时,应采用改进的可靠度指数计算公式,在可靠性不变的前提下将周期性变化工作环境温度等效变换为恒定型工作环境温度,从而确定基准应力。
所述的工作环境温度分为恒定型、三角波型、组合型;
其中恒定型工作环境温度下的空间电子设备的可靠度R计算公式为:
Figure BDA0000117213220000021
λ0是对应恒定工作环境温度T0的工作失效率,t是工作时间;
其中三角波型工作环境温度下的空间电子设备的可靠度R计算公式为:
Figure BDA0000117213220000022
λi是对应工作环境温度Ti=TL+n·ΔT的工作失效率,n为步阶数,n=1,2,…,(TU-TL)/ΔT;ΔT是步阶温度,取ΔT=5℃;当n=1时,T1=TL+5℃;当n=2时,T2=TL+10℃;…;当n=(TU-TL)/ΔT时,Tn=TU;t是工作时间;TU、TL分别是上限工作环境温度、下限工作环境温度,tc是工作环境温度变化一个周期的时间;
其中组合型工作环境温度是恒定型和三角波型的组合,该空间电子设备的可靠度R计算公式为:λ0是对应恒定工作环境温度TL的工作失效率,λi是对应工作环境温度Ti=TL+n·ΔT的工作失效率,n为步阶数,n=1,2,…,(TU-TL)/ΔT;tc是工作环境温度变化一个周期的时间,ΔT是步阶温度,取ΔT=5℃;当n=1时,T1=TL+5℃;当n=2时,T2=TL+10℃;…;当n=(TU-TL)/ΔT时,Tn=TU;t是工作时间。
本发明的有益效果:
1)采用等步阶分割方法由于步阶温度小,步阶数多,因而得到的空间电子设备可靠度结果比现有算法更准确,此法可以推广应用到更一般的变化环境温度下的电子设备的可靠性预计工作中;
2)基于改进的可靠度指数计算公式可以确定变化工作环境温度下的加速寿命试验的温度加速基准应力;
3)由于准确确定了温度加速基准应力,因而可以准确确定加速系数以及最终的可靠性评估结果;
4)温度加速基准应力的准确确定为加速寿命试验中包括子样数量、试验时间、人员安排等的试验成本确定问题提供了科学依据。
附图说明
图1是空间电子设备工作环境温度的种类图。
具体实施方式
下面结合实施例对本发明做进一步说明。
首先针对工作环境温度的周期性变化,将环境温度进行等步阶分割,其中分图(a)中的T0是恒定工作环境温度;分图(b)中的TU、TL分别是上限工作环境温度、下限工作环境温度,tc是工作环境温度变化一个周期的时间;分图(c)中的TU、TL分别是上限工作环境温度、下限工作环境温度,此时TL亦代表恒定工作环境温度,t0、tc分别是对应工作环境温度变化一个周期的时间(t0+tc)中的恒定工作环境温度的时间、变化工作环境温度的时间;按照技术要求,上述T0、TU、TL均是5℃的整倍数;需要说明的是:组合型是恒定型和三角波型的组合,这里一个周期中的恒定段可以是TL,也可以是TU,也可以是TL、TU或含有TL、TU以及介于二者之间的某一个或多个恒定段的组合;一个周期中的三角波型则代表温度上升段和下降段的组合,可以是一组或多组上升段和下降段的组合,且上升段和下降段的温度变化速率可以不一样;
恒定型工作环境温度下的空间电子设备的可靠度R计算公式为:
R = e - λ 0 t - - - ( 1 )
λ0是对应恒定工作环境温度T0的工作失效率,t是工作时间;在恒定型工作环境温度下进行加速寿命试验,此恒定工作环境温度即为温度加速基准应力,即是图1中分图(a)中的T0
三角波型工作环境温度下的空间电子设备的可靠度R计算公式为:
R = e - ( 1 n Σ i = 1 n λ i ) t - - - ( 2 )
λi是对应工作环境温度Ti=TL+n·ΔT的工作失效率,n为步阶数,n=1,2,…,(TU-TL)/ΔT;ΔT是步阶温度,取ΔT=5℃;当n=1时,T1=TL+5℃;当n=2时,T2=TL+10℃;…;当n=(TU-TL)/ΔT时,Tn=TU;t是工作时间;在三角波型工作环境温度下进行加速寿命试验,令公式(2)等号右边等于公式(1)等号右边,化简后则
λ 0 = 1 n Σ i = 1 n λ i - - - ( 4 )
公式(4)中λ0是三角波型等效为恒定型工作环境温度下的等效工作失效率,再按照就近就高的原则在λi中查出对应此等效工作失效率的实际工作失效率,则此实际工作失效率对应的工作环境温度Ti即为温度加速基准应力;
组合型工作环境温度下的空间电子设备的可靠度R计算公式为
Figure BDA0000117213220000043
λ0是对应恒定工作环境温度TL的工作失效率,λi是对应工作环境温度Ti=TL+n·ΔT的工作失效率,n为步阶数,n=1,2,…,(TU-TL)/ΔT;ΔT是步阶温度,取ΔT=5℃;当n=1时,T1=TL+5℃;当n=2时,T2=TL+10℃;…;当n=(TU-TL)/ΔT时,Tn=TU;t是工作时间;在组合型工作环境温度下进行加速寿命试验,令公式(3)等号右边等于公式(1)等号右边,为区别公式(3)中的λ0这里用λ0′代替公式(1)中的λ0,化简后则
λ 0 ′ = [ λ 0 t 0 + ( 1 n Σ i = 1 n λ i ) t c ] 1 t 0 + t c - - - ( 5 )
公式(5)中λ0′是组合型等效为恒定型工作环境温度下的等效工作失效率,再按照就近就高的原则在λi中查出对应此等效工作失效率的实际工作失效率,则此实际工作失效率对应的工作环境温度Ti即为温度加速基准应力。
再分别计算每个台阶温度下的对应时间段的可靠度;然后计算一个周期和整个寿命末期的可靠度,最后给出改进的可靠度指数计算公式;
工作环境温度曲线如图1中的分图(c)所示,其中t0=4380h,tc=4380h,t=131400h,TU=70℃,TL=0℃,其温度加速基准应力确定的具体步骤为:
1)将变化工作环境温度按照ΔT=5℃进行等步阶分割,步阶数n=(TU-TL)/ΔT=(70-0)℃/5℃=14;
2)采用应力分析法进行可靠度计算,获得DC-DC变换器的步阶数(n)-工作环境温度(Ti)-工作失效率(λi)数值之间的对应关系,见表1;
表1DC-DC变换器的步阶数-工作环境温度-工作失效率数值
Figure BDA0000117213220000051
注1:n=0表示不分割,此时为恒定工作环境温度。
注2:GJB/Z 299C-2006、MIL-HDBK-217F分别是GJB/Z 299C-2006,电子设备可靠性预计手册[S].北京:总装备部军标出版发行部,2006和MIL-HDBK-217F,ReliabilityPrediction of Electronic Equipment[S].Department of Defense,Washington D.C.,USA,1991。
注3:Fit(菲特)表示失效率的概念,其含义是1Fit=1×10-9/h,即103个器件工作106h,其中有1个失效。
3)代入表1中的数据,则 1 n Σ i = 1 n λ i = 1 14 Σ i = 1 14 λ i = 98.74 Fit ;
4)将以上已知数据代入公式(5)中,得到DC-DC变换器组合型等效为恒定型工作环境温度下的等效工作失效率λ0′:
λ 0 ′ = [ λ 0 t 0 + ( 1 n Σ i = 1 n λ i ) t c ] 1 t 0 + t c
= [ 43.36 × 10 - 9 × 4380 + ( 1 14 Σ i = 1 14 λ i ) × 4380 ] × 1 4380 + 4380
= [ 43.36 × 10 - 9 × 4380 + 89.74 × 10 - 9 × 4380 ] × 1 4380 + 4380
= 66.55 Fit ;
5)按照就近就高的原则在表1中查出λ6=67.68Fit与λ0′=66.55Fit最为接近且高于它,λ6即是对应此等效工作失效率λ0′的实际工作失效率,则λ6对应的工作环境温度T6=30℃即为温度加速基准应力。
当选择工作环境温度作为加速应力即通过提高工作环境温度进行空间电子设备的加速寿命试验时,应采用改进的可靠度指数计算公式,在可靠性不变的前提下将周期性变化工作环境温度等效变换为恒定型工作环境温度,从而确定基准应力。

Claims (2)

1.空间电子设备加速寿命试验中温度加速基准应力确定方法,其特征在于:首先针对工作环境温度的周期性变化,将环境温度进行等步阶分割,再分别计算每个台阶温度下的对应时间段的可靠度;然后计算一个周期和整个寿命末期的可靠度,最后给出改进的可靠度指数计算公式;当选择工作环境温度作为加速应力即通过提高工作环境温度进行空间电子设备的加速寿命试验时,应采用改进的可靠度指数计算公式,在可靠性不变的前提下将周期性变化工作环境温度等效变换为恒定型工作环境温度,从而确定基准应力。
2.如权利要求1所述的空间电子设备加速寿命试验中温度加速基准应力确定方法,其特征在于:所述的工作环境温度分为恒定型、三角波型、组合型;
其中恒定型工作环境温度下的空间电子设备的可靠度R计算公式为:
Figure FDA0000117213210000011
λ0是对应恒定工作环境温度T0的工作失效率,t是工作时间;
其中三角波型工作环境温度下的空间电子设备的可靠度R计算公式为:
Figure FDA0000117213210000012
λi是对应工作环境温度Ti=TL+n·ΔT的工作失效率,n为步阶数,n=1,2,…,(TU-TL)/ΔT;ΔT是步阶温度,取ΔT=5℃;当n=1时,T1=TL+5℃;当n=2时,T2=TL+10℃;…;当n=(TU-TL)/ΔT时,Tn=TU;t是工作时间;TU、TL分别是上限工作环境温度、下限工作环境温度,tc是工作环境温度变化一个周期的时间;
其中组合型工作环境温度是恒定型和三角波型的组合,该空间电子设备的可靠度R计算公式为:
Figure FDA0000117213210000013
λ0是对应恒定工作环境温度TL的工作失效率,λi是对应工作环境温度Ti=TL+n·ΔT的工作失效率,n为步阶数,n=1,2,…,(TU-TL)/ΔT;tc是工作环境温度变化一个周期的时间,ΔT是步阶温度,取ΔT=5℃;当n=1时,T1=TL+5℃;当n=2时,T2=TL+10℃;…;当n=(TU-TL)/ΔT时,Tn=TU;t是工作时间。
CN201110403814.6A 2011-12-07 2011-12-07 空间电子设备加速寿命试验中温度加速基准应力确定方法 Active CN102520279B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110403814.6A CN102520279B (zh) 2011-12-07 2011-12-07 空间电子设备加速寿命试验中温度加速基准应力确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110403814.6A CN102520279B (zh) 2011-12-07 2011-12-07 空间电子设备加速寿命试验中温度加速基准应力确定方法

Publications (2)

Publication Number Publication Date
CN102520279A true CN102520279A (zh) 2012-06-27
CN102520279B CN102520279B (zh) 2014-06-25

Family

ID=46291265

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110403814.6A Active CN102520279B (zh) 2011-12-07 2011-12-07 空间电子设备加速寿命试验中温度加速基准应力确定方法

Country Status (1)

Country Link
CN (1) CN102520279B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103678866A (zh) * 2012-09-06 2014-03-26 国际商业机器公司 用于计算系统产品可靠度估计的方法和系统
CN103852672A (zh) * 2014-03-31 2014-06-11 中国航天科工运载技术研究院北京分院 一种以步进应力对电子设备实施环境应力筛选的方法
CN105093028A (zh) * 2015-08-21 2015-11-25 北京航天长征飞行器研究所 一种电子类产品加速贮存的试验方法
CN105572498A (zh) * 2015-12-11 2016-05-11 中国航空工业集团公司西安飞机设计研究所 一种电子产品可靠性加速试验方法
CN106990011A (zh) * 2017-05-15 2017-07-28 中国人民解放军海军航空工程学院 一种加速寿命试验的温度加速应力谱的确定方法
CN109655683A (zh) * 2018-12-13 2019-04-19 广州广电计量检测股份有限公司 一种汽车电子产品热疲劳寿命加速测试方法
CN110889083A (zh) * 2018-09-10 2020-03-17 湖南银杏可靠性技术研究所有限公司 基于窗谱估计的加速贮存与自然贮存退化数据一致性检验法
CN113092911A (zh) * 2021-04-07 2021-07-09 西安苏试广博环境可靠性实验室有限公司 空间电子设备加速寿命试验中温度加速基准应力确定方法
CN113984833A (zh) * 2021-10-29 2022-01-28 江苏徐工工程机械研究院有限公司 一种环境温度等效及加速试验方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101795060A (zh) * 2010-03-22 2010-08-04 北京工业大学 混合集成电路dc/dc变换器可靠性快速评价方法
CN102023636A (zh) * 2010-12-03 2011-04-20 浙江理工大学 一种机床数控系统的加速寿命试验方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101795060A (zh) * 2010-03-22 2010-08-04 北京工业大学 混合集成电路dc/dc变换器可靠性快速评价方法
CN102023636A (zh) * 2010-12-03 2011-04-20 浙江理工大学 一种机床数控系统的加速寿命试验方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《Proceedings of ISUMA-NAFIPS'95 IEEE》 19951231 Huamin, Liu 等 Reliability Assessment of Systems Operating in Variable Conditons 第6页 1 , *
HUAMIN, LIU 等: "Reliability Assessment of Systems Operating in Variable Conditons", 《PROCEEDINGS OF ISUMA-NAFIPS’95 IEEE》, 31 December 1995 (1995-12-31), pages 6 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103678866A (zh) * 2012-09-06 2014-03-26 国际商业机器公司 用于计算系统产品可靠度估计的方法和系统
CN103678866B (zh) * 2012-09-06 2017-06-09 国际商业机器公司 用于计算系统产品可靠度估计的方法和系统
CN103852672A (zh) * 2014-03-31 2014-06-11 中国航天科工运载技术研究院北京分院 一种以步进应力对电子设备实施环境应力筛选的方法
CN105093028B (zh) * 2015-08-21 2018-08-07 北京航天长征飞行器研究所 一种电子类产品加速贮存的试验方法
CN105093028A (zh) * 2015-08-21 2015-11-25 北京航天长征飞行器研究所 一种电子类产品加速贮存的试验方法
CN105572498A (zh) * 2015-12-11 2016-05-11 中国航空工业集团公司西安飞机设计研究所 一种电子产品可靠性加速试验方法
CN105572498B (zh) * 2015-12-11 2018-08-24 中国航空工业集团公司西安飞机设计研究所 一种电子产品可靠性加速试验方法
CN106990011A (zh) * 2017-05-15 2017-07-28 中国人民解放军海军航空工程学院 一种加速寿命试验的温度加速应力谱的确定方法
CN110889083A (zh) * 2018-09-10 2020-03-17 湖南银杏可靠性技术研究所有限公司 基于窗谱估计的加速贮存与自然贮存退化数据一致性检验法
CN109655683A (zh) * 2018-12-13 2019-04-19 广州广电计量检测股份有限公司 一种汽车电子产品热疲劳寿命加速测试方法
CN113092911A (zh) * 2021-04-07 2021-07-09 西安苏试广博环境可靠性实验室有限公司 空间电子设备加速寿命试验中温度加速基准应力确定方法
CN113092911B (zh) * 2021-04-07 2023-08-15 西安苏试广博环境可靠性实验室有限公司 空间电子设备加速寿命试验中温度加速基准应力确定方法
CN113984833A (zh) * 2021-10-29 2022-01-28 江苏徐工工程机械研究院有限公司 一种环境温度等效及加速试验方法
CN113984833B (zh) * 2021-10-29 2024-03-01 江苏徐工工程机械研究院有限公司 一种环境温度等效及加速试验方法

Also Published As

Publication number Publication date
CN102520279B (zh) 2014-06-25

Similar Documents

Publication Publication Date Title
CN102520279B (zh) 空间电子设备加速寿命试验中温度加速基准应力确定方法
Postnikov et al. Tidal Love numbers of neutron and self-bound quark stars
CN102622473A (zh) 基于贝叶斯理论的步进应力加速退化试验优化设计方法
CN104156500A (zh) 一种预测材料疲劳寿命的方法
CN103913328A (zh) 广义位移混合监测受损索载荷递进式识别方法
CN103530449A (zh) 一种弹上寿命件的多变量加速贮存试验优化设计方法
CN101782946A (zh) 有支座沉降时空间坐标监测识别松弛的支承索的递进式方法
CN105022792A (zh) 基于数据挖掘的被动雷达信号分选关联权重计算方法
CN105547325A (zh) 一种基于k均值聚类的光纤惯组温度模型系数确定方法
CN104134013A (zh) 一种风力机叶片模态分析方法
CN103257000A (zh) 一种桥梁结构日照作用分析的极值温度预测方法
CN103852317A (zh) 角位移混合监测问题索载荷递进式识别方法
CN101788404B (zh) 有支座沉降时角度监测的识别松弛的支承索的递进式方法
CN103913327A (zh) 线位移空间坐标监测问题索载荷递进式识别方法
CN103868732A (zh) 线位移应变监测问题索载荷递进式识别方法
CN104331775A (zh) 一种计量资产全寿命周期质量分析方法
CN103868707A (zh) 应变监测问题索载荷递进式识别方法
CN101629407A (zh) 一种路面结构强度预测方法
CN103852321A (zh) 角度监测问题索载荷递进式识别方法
CN103868714A (zh) 线位移应变监测受损索载荷递进式识别方法
CN103868722A (zh) 角位移空间坐标监测问题索载荷递进式识别方法
CN103852308A (zh) 混合监测问题索载荷线位移识别方法
CN103604655A (zh) 线位移混合监测受损索集中载荷递进式识别方法
CN102495715B (zh) 双精度浮点定义下的深空Doppler测速计算方法
CN105975697A (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