CN112632820A - 一种用于壁板结构减振降噪的动力吸振器参数设计方法 - Google Patents
一种用于壁板结构减振降噪的动力吸振器参数设计方法 Download PDFInfo
- Publication number
- CN112632820A CN112632820A CN202011508230.0A CN202011508230A CN112632820A CN 112632820 A CN112632820 A CN 112632820A CN 202011508230 A CN202011508230 A CN 202011508230A CN 112632820 A CN112632820 A CN 112632820A
- Authority
- CN
- China
- Prior art keywords
- modal
- vibration absorber
- dynamic
- vibration
- space
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Buildings Adapted To Withstand Abnormal External Influences (AREA)
Abstract
本发明公开了一种用于壁板结构减振降噪的动力吸振器(DVA)参数设计方法,获取壁板结构上DVA最优位置及材料参数。首先,对壁板结构模态空间开展基变换,构建新的n维正交完备空间,使得外力在新建n维空间中前n‑1维空间的投影为零,从而将壁板结构DVA参数设计问题转换为单模态振动系统的DVA参数设计问题;进一步,通过等效参数代换将单模态振动系统的DVA参数设计问题转换为单自由度振动系统的DVA参数设计问题;最后,利用单自由度振动系统DVA参数优化设计中的定点理论完成壁板结构DVA最优材料参数设计,通过单模态极值确定DVA布置位置。本发明解决了复杂壁板结构由于响应参与模态较多而难以确定DVA最优位置及材料参数的难题,克服了现有DVA参数设计中难以适用于多模态响应结构的局限性。
Description
技术领域
本发明涉及一种用于壁板结构减振降噪的动力吸振器参数设计方法,属于结构减振技术领域。
背景技术
动力吸振器由于其理论相对成熟、成本低、安装方便等优点,被广泛应用于各类工程减振降噪中,但适合于壁板结构减振的理论研究并不多。目前较多采用的方案为采用模态参数识别的方法,将主系统等效为单自由度系统进行优化设计,但该方法受限于识别算法的精度且针对耦合模态识别精度不高,很难应用到壁板结构上。基于理论推导的动力吸振器参数设计方法目前仅适用板、梁等简单结构,难以推广到壁板结构。因此,为了更好的将动力吸振器应用到壁板结构,亟待发展一种适用于壁板结构的动力吸振器参数设计方法。
发明内容
发明目的:本发明的目的是提供一种用于壁板结构减振降噪的动力吸振器参数设计方法。
技术方案:为实现以上目的,本发明公开了一种用于壁板结构减振降噪的动力吸振器参数设计方法,该方法包含如下步骤:
(1)、建立壁板结构有限元模型,进行模态分析,获取其模态参数;
(2)、对壁板结构施加外力,获取结构响应,将结构响应和外力投影到模态空间得到广义坐标和广义力;根据研究频段内模态的参与度,即广义坐标的大小建立模态空间内n维动力学模型;
(3)、对壁板结构在模态空间内的n维动力学模型进行基变换,使得在变换后的n维空间中,外力在前n-1维空间的投影为零,将模态空间内的n维动力学模型的动力吸振器参数设计问题降维成新基向量张成空间内单模态振动系统的动力吸振器参数设计问题;
(4)、通过等效参数替换将降维后的单模态振动系统动力吸振器参数设计问题转化为等效单自由度振动系统的动力吸振器参数设计问题,通过模态振型函数极值处所对应空间坐标来确定动力吸振器布置位置;
(5)、根据定点理论进行等效单自由度振动系统减振的动力吸振器参数设计,完成壁板结构减振降噪的动力吸振器参数设计。所述壁板结构为加筋板、机翼蒙皮等壁板结构。
进一步的,所述步骤(1)中,对结构进行有限元建模并利用试验数据修正初始模型,开展模态分析获取结构准确的模态参数。
进一步的,所述步骤(2)中广义力的计算方法及模态空间的维度n选取包括以下步骤:
fg(x,y)=∫sf(x,y)φj(x,y)ds (1)
式中,fg(x,y)为广义力;f(x,y)为分布力激励;φj(x,y)为第j阶结构模态振型;∫s·ds表示在壁板结构面积域s内进行积分,则第j阶模态下结构响应的广义坐标qj可以表示为:
式中,ω为外激励频率,Mj、Cj和Kj分别为第j阶模态质量、模态阻尼和模态刚度;j为虚数单位,广义坐标qj代表了该阶模态对响应的参与度,根据广义坐标的值由大到小挑选模态,使得选取的前n阶模态对该响应的参与度总和大于90%,即确定n的数值,将其它阶模态舍弃,建立模态空间内的n维动力学模型。
进一步的,所述步骤(3)中基变换包括以下步骤:
对于模态空间内的n维动力学模型,利用其模态振型φi(x,y)的线性组合构造新的n维完备正交基ψj(x,y),则有:
式中,aji表示线性组合中对应φi(x,y)的权重系数,使得分布力激励f(x,y)在ψj(x,y),j=1,2,…,n-1上的投影为零,即满足:
此时,分布力激励f(x,y)只在ψn(x,y)上的投影为非零,定义ψn(x,y)为与激励相关的特征模态,同时,针对新构造的正交基,应满足如下正交性条件:
∫sρhψi(x,y)ψj(x,y)ds=δij,i,j=1,2,…n (5)
式中,ρ为密度,h为壁板结构的厚度,δij为狄拉克函数,由式(3)和式(5)可得:
此时,模态空间内n维动力学模型通过基变换降维成新基向量张成空间内的单模态振动系统,该系统的特征模态ψn(x,y)的表达式为:
式中,Ki为降维前模态空间内的n维动力学模型的第i阶模态刚度。
进一步的,所述步骤(4)中等效参数替换包括以下步骤:
基于单模态振动系统,在壁板结构空间位置(x1,y1)处附加质量、刚度和阻尼分别为m1、k1和c1的动力吸振器后,其振动响应复振幅X的表达式为:
式(12)中,等式右边乘积的第一项为关于动力吸振器放置位置的表达式,要使被减振结构振动量级最小,应使ψ2 n(x1,y1)最大,即动力吸振器应放置在振型幅值最大处。
进一步的,所述步骤(5)中定点理论包括以下步骤:
对一单自由度振动系统,其质量、刚度和阻尼也分别等于Mref、Kref、Cref,附加质量、刚度和阻尼分别为m1、k1和c1的动力吸振器后,其振动响应复振幅Xs的表达式为:
式中,F为外激励的幅值,取复振幅Xs的模可得:
式中
式中,为动力吸振器固有频率,为单自由度振动系统固有频率;为动力吸振器阻尼比;为动力吸振器与单自由度振动系统的质量比;单自由度振动系统的阻尼比,式(14)所示单自由度振动系统响应复振幅幅值表达式存在两个不随吸振器阻尼比改变的定点,定点的位置可由下式得出:
令这两个定点的|Xs|取值相等,并使曲线|Xs|在这两定点的导数为零定为最优设计条件,根据最优设计条件可得单自由度振动系统动力吸振器最优材料参数,式(13)与式(12)仅相差一个倍数关系,此时单自由度振动系统动力吸振器最优材料参数即为壁板结构上附加动力吸振器的最优材料参数:
式中
有益效果:与现有技术相比,本发明的技术方案具有以下有益技术效果:
本发明通过基变换将模态空间内n维结构降维成新基坐标下一维结构,再通过等效参数替换将一维结构降维成零维结构,使得针对模态空间内n维结构的动力吸振器优化设计问题降维成零维单自由度主系统动力吸振器优化设计问题,大大降低模型的复杂度。本发明可以针对壁板结构进行动力吸振器优化设计,解决了现有动力吸振器理论在壁板结构中应用的局限性。
附图说明
图1为本发明的流程图;
图2为本发明算例中舱体结构有限元模型;
图3为本发明算例模型第三阶振型/37Hz;
图4为本发明算例模型第四阶振型/37Hz;
图5为本发明算例模型第五阶振型/46Hz;
图6为本发明算例模型第六阶振型/46Hz;
图7为本发明算例模型广义坐标;
图8为本发明算例模型特征模态振型1/40Hz;
图9为本发明算例模型特征模态振型2/46Hz;
图10为本发明算例模型特征模态振型1幅值处振动响应;
图11为本发明算例模型特征模态振型2幅值处振动响应;
图12为本发明算例模型辐射声功率曲线。
具体实施方式
下面对本发明技术方案进行详细说明,但是本发明的保护范围不局限于所述实施例。
实施例:一种用于壁板结构减振降噪的动力吸振器参数设计方法,选取复杂舱体结构为分析对象,具体操作如下:
(1)、对结构进行模态分析,获取模态参数。采用壳单元离散舱体结构模型,舱体结构有限元模型如图2所示;3-6阶模态振型如图3-6所示。
(2)、将结构响应和外力投影到模态空间得到广义坐标和广义力;根据研究频段内广义坐标的大小建立模态空间内n维结构模型。广义力的计算方法及空间维度的选取包括以下步骤:
fg(x,y)=∫sf(x,y)φj(x,y)ds (1)
式中,fg(x,y)为广义力;f(x,y)为分布力激励;φj(x,y)为第j阶模态振型;∫s·ds表示对域内进行积分,则第j阶模态下的广义坐标可以表示为:
式中,ω为外激励频率;Mj、Cj和Kj分别为第j阶模态质量、模态阻尼和模态刚度;j为虚数单位。将步骤1中模态参数代入上式可得广义坐标随频率变换曲线,如图7所示。当研究频段为35.9Hz附近频段时,添加动力吸振器会使声响应在50.5Hz出现较大响应峰值,因此将研究频段选为35.9Hz和50.5Hz附近频段。图7显示35.9Hz附近频段3、4、5阶模态参与度较高,且其对响应的贡献占比达90%以上,因此建立模态空间内三维结构模型;50.5Hz频段内第5阶模态参与度较高,因此建立模态空间内一维结构模型。
(3)、对n维模态空间基坐标进行基变换,使得外力在变换后n维空间中的n-1维投影为零;从而将模态空间内n维结构降维成新基向量张成空间内的一维结构。基变换包括以下步骤:
对于模态空间内n维结构,构造新的n维正交完备基表示该n维结构,则有:
式中,φi(x,y)为结构的第i阶振型函数,ψj(x,y)为新构造的完备正交基函数,aji表示对应系数。使得分布力f(x,y)在新构造的前n-1阶完备正交基上投影为零,即满足:
此时分布力f(x,y)只在ψn(x,y)上的投影非零,此时的ψn(x,y)定义为特征模态。同时,针对新构造的正交基,应满足如下正交性条件:
∫sρhψi(x,y)ψj(x,y)ds=δij (5)
式中,δij为狄拉克函数。将ψj(x,y)代入上式进一步可得
式中:k≠l;Mi为结构的第i阶模态质量。联立可得特征模态ψn(x,y)系数ani的表达式:
此时模态空间内n维结构通过基变换降维成新基向量张成空间内的一维结构,该一维空间内所对应的特征模态振型的表达式为:
一维空间内对应的特征模态质量和特征模态刚度可表示为:
式中,Ki为结构第i阶模态刚度。
将步骤2建立的模态空间内n维结构模型代入上式可得特征模态振型1、2,如图8、图9所示。
(4)、通过等效参数替换将一维结构优化问题降维成零维单自由度主系统优化设计问题。等效参数替换包括以下步骤:
根据单模态工况下的包含动力吸振器的响应表达式:
式中,ψn(x,y)为主系统的第n阶特征模态振型;分别为主系统第n阶特征模态所对应的模态质量、模态刚度和模态阻尼;m1、k1和c1分别为动力吸振器的质量、刚度和阻尼;(x1,y1)表示动力吸振器放置位置;j为虚数单位。将式 代入上式进行等效参数替换可得:
上式中等式右边乘积的第一项为关于动力吸振器放置位置的表达式,要使主系统振动量级最小,应使ψ2 n(x1,y1)最大,即动力吸振器应放置在振型幅值最大处;第二项与单自由度主系统动力吸振器设计方法表达式形式完全一致。
将步骤3建立的特征模态代入上式可得动力吸振器最优位置及等效单自由度主系统动力吸振器优化设计模型。
(5)、根据单自由度主系统动力吸振器优化理论中的定点理论完成对动力吸振器参数的优化。定点理论指单自由度主系统附加动力吸振器结构的传递函数存在两个定点,通过设计使得这两定点取值相等且为域内极大值。将步骤4建立的单自由度主系统动力吸振器优化设计模型代入定点理论可得动力吸振器最优材料参数如表1所示。
表1动力吸振器最优材料参数
图10-12结果可知,添加动力舱体模型在35.9Hz力学响应及声学响应均大幅降低。
Claims (6)
1.一种用于壁板结构减振降噪的动力吸振器参数设计方法,其特征在于,该方法包含如下步骤:
(1)、建立壁板结构有限元模型,进行模态分析,获取其模态参数;
(2)、对壁板结构施加外力,获取结构响应,将结构响应和外力投影到模态空间得到广义坐标和广义力;根据研究频段内模态的参与度,即广义坐标的大小建立模态空间内n维动力学模型;
(3)、对壁板结构在模态空间内的n维动力学模型进行基变换,使得在变换后的n维空间中,外力在前n-1维空间的投影为零,将模态空间内的n维动力学模型的动力吸振器参数设计问题降维成新基向量张成空间内单模态振动系统的动力吸振器参数设计问题;
(4)、通过等效参数替换将降维后的单模态振动系统动力吸振器参数设计问题转化为等效单自由度振动系统的动力吸振器参数设计问题,通过模态振型函数极值处所对应空间坐标来确定动力吸振器布置位置;
(5)、根据定点理论进行等效单自由度振动系统减振的动力吸振器参数设计,完成壁板结构减振降噪的动力吸振器参数设计。
2.根据权利要求1所述的一种用于壁板结构减振降噪的动力吸振器参数设计方法,其特征在于:所述步骤(1)中,对结构进行有限元建模并利用试验数据修正初始模型,开展模态分析获取结构准确的模态参数。
3.根据权利要求1所述的一种用于壁板结构减振降噪的动力吸振器参数设计方法,其特征在于,步骤(2)中广义力的计算方法及模态空间的维度n选取包括以下步骤:
fg(x,y)=∫sf(x,y)φj(x,y)ds (1)
式中,fg(x,y)为广义力;f(x,y)为分布力激励;φj(x,y)为第j阶结构模态振型;∫s·ds表示在壁板结构面积域s内进行积分,则第j阶模态下结构响应的广义坐标qj可以表示为:
式中,ω为外激励频率,Mj、Cj和Kj分别为第j阶模态质量、模态阻尼和模态刚度;j为虚数单位,广义坐标qj代表了该阶模态对响应的参与度,根据广义坐标的值由大到小挑选模态,使得选取的前n阶模态对该响应的参与度总和大于90%,即确定n的数值,将其它阶模态舍弃,建立模态空间内的n维动力学模型。
4.根据权利要求3所述的一种用于壁板结构减振降噪的动力吸振器参数设计方法,其特征在于,所述步骤(3)中基变换包括以下步骤:
对于模态空间内的n维动力学模型,利用其模态振型φi(x,y)的线性组合构造新的n维完备正交基ψj(x,y),则有:
式中,aji表示线性组合中对应φi(x,y)的权重系数,使得分布力激励f(x,y)在ψj(x,y),j=1,2,…,n-1上的投影为零,即满足:
此时,分布力激励f(x,y)只在ψn(x,y)上的投影为非零,定义ψn(x,y)为与激励相关的特征模态,同时,针对新构造的正交基,应满足如下正交性条件:
∫sρhψi(x,y)ψj(x,y)ds=δij,i,j=1,2,…n (5)
式中,ρ为密度,h为壁板结构的厚度,δij为狄拉克函数,由式(3)和式(5)可得:
此时,模态空间内n维动力学模型通过基变换降维成新基向量张成空间内的单模态振动系统,该系统的特征模态ψn(x,y)的表达式为:
式中,Ki为降维前模态空间内的n维动力学模型的第i阶模态刚度。
6.根据权利要求5所述的一种用于壁板结构减振降噪的动力吸振器参数设计方法,其特征在于:所述步骤(5)中定点理论包括以下步骤:
对一单自由度振动系统,令其质量、刚度和阻尼分别等于Mref、Kref、Cref,附加质量、刚度和阻尼分别为m1、k1和c1的动力吸振器后,其振动响应复振幅Xs的表达式为:
式中,F为外激励的幅值,取复振幅Xs的模可得:
式中
式中,为动力吸振器固有频率,为单自由度振动系统固有频率;为动力吸振器阻尼比;为动力吸振器与单自由度振动系统的质量比;单自由度振动系统的阻尼比,式(14)所示单自由度振动系统响应复振幅幅值表达式存在两个不随动力吸振器阻尼比改变的定点,定点的位置可由下式得出:
令这两个定点的|Xs|取值相等,并使曲线|Xs|在这两定点的导数为零定为最优设计条件,根据最优设计条件可得单自由度振动系统动力吸振器最优材料参数,式(13)与式(12)仅相差一个倍数关系,此时单自由度振动系统动力吸振器最优材料参数即为壁板结构上附加动力吸振器的最优材料参数:
式中
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011508230.0A CN112632820B (zh) | 2020-12-18 | 2020-12-18 | 一种用于壁板结构减振降噪的动力吸振器参数设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011508230.0A CN112632820B (zh) | 2020-12-18 | 2020-12-18 | 一种用于壁板结构减振降噪的动力吸振器参数设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112632820A true CN112632820A (zh) | 2021-04-09 |
CN112632820B CN112632820B (zh) | 2022-11-18 |
Family
ID=75318123
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011508230.0A Active CN112632820B (zh) | 2020-12-18 | 2020-12-18 | 一种用于壁板结构减振降噪的动力吸振器参数设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112632820B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116090257A (zh) * | 2023-03-07 | 2023-05-09 | 中大智能科技股份有限公司 | 一种压路机-土基材料系统的理论模型 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107480322A (zh) * | 2017-06-23 | 2017-12-15 | 中国工程物理研究院总体工程研究所 | 自由体多点相关脉动压力随机振动分析计算方法 |
CN110083994A (zh) * | 2019-05-31 | 2019-08-02 | 吉林大学 | 一种行星混联式混合动力系统扭振优化方法 |
CN111209668A (zh) * | 2020-01-06 | 2020-05-29 | 北京工业大学 | 针对复杂薄壁结构件多模态时变切削振动的动力吸振器参数的优化方法 |
-
2020
- 2020-12-18 CN CN202011508230.0A patent/CN112632820B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107480322A (zh) * | 2017-06-23 | 2017-12-15 | 中国工程物理研究院总体工程研究所 | 自由体多点相关脉动压力随机振动分析计算方法 |
CN110083994A (zh) * | 2019-05-31 | 2019-08-02 | 吉林大学 | 一种行星混联式混合动力系统扭振优化方法 |
CN111209668A (zh) * | 2020-01-06 | 2020-05-29 | 北京工业大学 | 针对复杂薄壁结构件多模态时变切削振动的动力吸振器参数的优化方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116090257A (zh) * | 2023-03-07 | 2023-05-09 | 中大智能科技股份有限公司 | 一种压路机-土基材料系统的理论模型 |
CN116090257B (zh) * | 2023-03-07 | 2023-06-16 | 中大智能科技股份有限公司 | 一种压路机-土基材料系统的理论模型 |
Also Published As
Publication number | Publication date |
---|---|
CN112632820B (zh) | 2022-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Van Genechten et al. | A direct hybrid finite element–wave based modelling technique for efficient coupled vibro-acoustic analysis | |
CN107885945B (zh) | 面向3d打印的多尺度孔洞结构轻量化建模方法 | |
Muravyov et al. | Determination of nonlinear stiffness with application to random vibration of geometrically nonlinear structures | |
Van Genechten et al. | An efficient Wave Based Method for solving Helmholtz problems in three-dimensional bounded domains | |
CN109918712B (zh) | 一种基于p型有限元法和围线积分法求解复合型应力强度因子的计算方法 | |
CN108846192B (zh) | 一种结构任意阻尼处理的船舶三维声弹性分析方法 | |
Wang et al. | A coupled smoothed finite element method (S-FEM) for structural-acoustic analysis of shells | |
Xie et al. | A variational formulation for vibro-acoustic analysis of a panel backed by an irregularly-bounded cavity | |
Panagiotopoulos et al. | Krylov subspaces recycling based model order reduction for acoustic BEM systems and an error estimator | |
Liew et al. | Mesh‐free radial point interpolation method for the buckling analysis of Mindlin plates subjected to in‐plane point loads | |
Woo et al. | Free vibration of skew Mindlin plates by p-version of FEM | |
CN112632820B (zh) | 一种用于壁板结构减振降噪的动力吸振器参数设计方法 | |
CN112528411A (zh) | 一种基于模态减缩的几何非线性结构噪声振动响应计算方法 | |
Vergote et al. | An efficient wave based approach for the time-harmonic vibration analysis of 3D plate assemblies | |
Zheng et al. | Topology optimization of PCLD on plates for minimizing sound radiation at low frequency resonance | |
Zhang et al. | Acoustic topology optimization of sound power using mapped acoustic radiation modes | |
CN115062500A (zh) | 分布式随机激励下结构振动响应分析方法 | |
Mi et al. | An interpolation method for coupling non-conforming patches in isogeometric analysis of vibro-acoustic systems | |
Park et al. | A mode selection criterion based on flexibility approach in component mode synthesis | |
Roque et al. | Transient analysis of composite and sandwich plates by radial basis functions | |
Aumann et al. | Robust error assessment for reduced order vibro-acoustic problems | |
Ng et al. | On the accuracy and stability of a variety of differential quadrature formulations for the vibration analysis of beams | |
Shang et al. | Stochastic dynamic analysis of acoustic-structural coupled systems under non-stationary random excitations | |
Schmid et al. | Three-dimensional stability of a Burgers vortex | |
CN108959776B (zh) | 一种声压的耦合建模方法、装置、设备及介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |