CN114778048A - 一种飞机减振构件的动刚度和一致性测试方法 - Google Patents

一种飞机减振构件的动刚度和一致性测试方法 Download PDF

Info

Publication number
CN114778048A
CN114778048A CN202210686028.XA CN202210686028A CN114778048A CN 114778048 A CN114778048 A CN 114778048A CN 202210686028 A CN202210686028 A CN 202210686028A CN 114778048 A CN114778048 A CN 114778048A
Authority
CN
China
Prior art keywords
vibration
end region
matrix
triaxial acceleration
test
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
CN202210686028.XA
Other languages
English (en)
Other versions
CN114778048B (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.)
AVIC Aircraft Strength Research Institute
Original Assignee
AVIC Aircraft Strength Research Institute
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 AVIC Aircraft Strength Research Institute filed Critical AVIC Aircraft Strength Research Institute
Priority to CN202210686028.XA priority Critical patent/CN114778048B/zh
Publication of CN114778048A publication Critical patent/CN114778048A/zh
Application granted granted Critical
Publication of CN114778048B publication Critical patent/CN114778048B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M7/00Vibration-testing of structures; Shock-testing of structures
    • G01M7/02Vibration-testing by means of a shake table
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64FGROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
    • B64F5/00Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
    • B64F5/60Testing or inspecting aircraft components or systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M5/00Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
    • G01M5/0066Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by exciting or detecting vibration or acceleration
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M5/00Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
    • G01M5/0075Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by means of external apparatus, e.g. test benches or portable test systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M7/00Vibration-testing of structures; Shock-testing of structures
    • G01M7/02Vibration-testing by means of a shake table
    • G01M7/025Measuring arrangements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M7/00Vibration-testing of structures; Shock-testing of structures
    • G01M7/02Vibration-testing by means of a shake table
    • G01M7/027Specimen mounting arrangements, e.g. table head adapters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M7/00Vibration-testing of structures; Shock-testing of structures
    • G01M7/02Vibration-testing by means of a shake table
    • G01M7/06Multidirectional test stands
    • 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

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Manufacturing & Machinery (AREA)
  • Transportation (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种飞机减振构件的动刚度和一致性测试方法,包括步骤:一、金属丝网减振器组的安装;二、测试区域的划分;三、三轴加速度传感器组的安装;四、测量各测试区域的频响函数;五、金属丝网减振器组的一致性判断;六、各减振器的动刚度计算;七、动刚度设计矩阵的评估。本发明通过在底座板和测试板之间安装多个减振器,再利用简单的频响函数测试,即可同时测试多个减振器,使用的装置结构简单,安装拆卸方便,测试速度快,测试数据处理方式统一,可大幅提高减振器动刚度和一致性测试及动刚度设计矩阵的评估效率,有利于满足一致性要求的减振器产品应用于飞机构件,提高系统的动力学性能一致性。

Description

一种飞机减振构件的动刚度和一致性测试方法
技术领域
本发明属于结构动力学测试技术领域,具体涉及一种飞机减振构件的动刚度和一致性测试方法。
背景技术
飞机减振构件常用金属减振器,如金属丝网减振器、弹簧复合减振器等的工作原理都是通过自身变形吸收振动能量,再用自身阻尼消耗能量。以金属丝网减振器为例,其是通过特殊工艺制得具有疏松多孔特点的功能型阻尼材料,成型产品内的金属丝之间呈现相互啮合类似橡胶的大分子孔隙结构。振动过程中,金属丝网减振器弹性变形的同时金属丝之间相互摩擦起到阻尼减振作用。
金属丝网减振器同款产品出厂时,在安装到飞机上前,需要测试产品性能的一致性和影响减振性能核心的动刚度,另外,也需要检查动刚度初始设计值是否符合要求。传统采用的测试方式往往是将产品安装在动刚度测试平台上逐个试验测试,然后对试验数据处理来获取同型号产品的动刚度和一致性并评估动刚度设计值,需要进行大量的重复测试作业,测试效率低,且不同批次的测试对测试结果也存在影响,因此,现需要一种飞机减振构件的动刚度和一致性测试装置及方法,来测试金属丝网减振器组的动刚度和一致性,并评估动刚度设计矩阵。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足,提供一种飞机减振构件的动刚度和一致性测试方法,通过在底座板和测试板之间安装多个减振器,再利用简单的频响函数测试,即可同时测试多个减振器,使用的装置结构简单,安装拆卸方便,测试速度快,测试数据处理方式统一,可大幅提高减振器动刚度和一致性测试及动刚度设计矩阵的评估效率,有利于满足一致性要求的减振器产品应用于飞机构件,提高系统的动力学性能一致性。
为解决上述技术问题,本发明采用的技术方案是:一种飞机减振构件的动刚度和一致性测试方法,所述飞机减振构件包括金属丝网减振器组,该方法采用的装置包括用于安装金属丝网减振器组的底座板和连接在所述金属丝网减振器组顶部的测试板;金属丝网减振器组包括多个减振器,所述减振器包括减振器主体、减振器安装座和减振器安装轴;
其特征在于:该方法包括以下步骤:
步骤一、金属丝网减振器组的安装:将N个所述减振器沿底座板的上表面周侧均匀安装,减振器安装轴的顶部穿过测试板并与测试板连接;N为不小于3的正整数;
步骤二、测试区域的划分:
将减振器安装座的上表面记做主动端区;将测试板上表面位于减振器安装轴周侧的环形区域记做被动端区,减振器安装轴的轴线穿过被动端区的环心,被动端区的环宽为10mm~15mm;
步骤三、三轴加速度传感器组的安装:
在各所述被动端区内均粘贴一个第一三轴加速度传感器,在各所述主动端区内均粘贴一个第二三轴加速度传感器,在测试板上表面的几何中心粘贴中心三轴加速度传感器,将粘贴中心三轴加速度传感器的区域记做中心测试区,第二三轴加速度传感器、第一三轴加速度传感器和中心三轴加速度传感器的轴向方向均平行、切向方向均平行、径向方向均平行;
步骤四、测量各测试区域的频响函数:
以第二三轴加速度传感器的安装中心点作为主动端区的响应点,以第一三轴加速度传感器的安装中心点作为被动端区的响应点,以中心三轴加速度传感器的安装中心点作为中心测试区的响应点;测量主动端区的响应点、被动端区的响应点和中心测试区的响应点之间的频响函数,组成三轴方向上的耦合频响函数矩阵;
步骤五、金属丝网减振器组的一致性判断:根据步骤四中获取的耦合频响函数矩阵,首先分别计算各个减振器在三轴方向上的振动传递率,随后分别计算N个减振器在三轴方向上的振动传递率均值;比较各轴方向上各个减振器的振动传递率和N个减振器的振动传递率均值之间的误差,当三轴方向上误差均小于等于误差允许值时,则表明当前测试的金属丝网减振器组满足一致性要求,否则不满足一致性要求;
步骤六、各减振器的动刚度计算:根据步骤四中获取的三轴方向上的耦合频响函数矩阵计算三轴方向上的解耦频响函数矩阵,根据三轴方向上的解耦频响函数矩阵和各个减振器在三轴方向上的振动传递率矩阵计算各个减振器在三轴方向上的动刚度矩阵;
步骤七、动刚度设计矩阵的评估:根据各个减振器在三轴方向上的动刚度矩阵评估各个减振器在三轴方向上的动刚度设计矩阵。
上述的方法,其特征在于,步骤四的具体步骤为:
步骤401、使用力锤在各主动端区分别施加振动激励,N个第二三轴加速度传感器接收三轴加速度信号,得到主动端区的响应点在三轴方向上的耦合频响函数矩阵
Figure 145059DEST_PATH_IMAGE001
;其中,a为各主动端区的编号且a=1,2,...,N;A=1,2,3,当A为1时表示沿切向方向,当A为2时表示沿径向方向,当A为3时表示沿轴向方向;
步骤402、使用力锤在各主动端区分别施加振动激励,N个第一三轴加速度传感器接收三轴加速度信号,得到主动端区的激励点到被动端区的响应点在三轴方向上的耦合频响函数矩阵
Figure 745804DEST_PATH_IMAGE002
;其中,p为各被动端区的编号且p=1,2,...,N;
步骤403、使用力锤在各主动端区分别施加振动激励,中心三轴加速度传感器接收三轴加速度信号,得到主动端区的激励点到中心测试区的响应点在三轴方向上的耦合频响函数矩阵
Figure 271464DEST_PATH_IMAGE003
;其中,t为中心测试区的编号且t=1;
步骤404、使用力锤在各被动端区分别施加振动激励,N个第一三轴加速度传感器接收三轴加速度信号,得到被动端区的响应点在三轴方向上的耦合频响函数矩阵
Figure 615857DEST_PATH_IMAGE004
步骤405、使用力锤在各被动端区分别施加振动激励,中心三轴加速度传感器接收三轴加速度信号,得到被动端区的激励点到中心测试区的响应点在三轴方向上的耦合频响函数矩阵
Figure 113835DEST_PATH_IMAGE005
步骤406、使用力锤在各被动端区分别施加振动激励,N个第二三轴加速度传感器接收三轴加速度信号,得到被动端区的激励点到主动端区的响应点在三轴方向上的耦合频响函数矩阵
Figure 619902DEST_PATH_IMAGE006
上述的方法,其特征在于,步骤五中进行金属丝网减振器组的一致性判断的具体步骤为:
步骤501、根据公式
Figure 836120DEST_PATH_IMAGE007
,计算N个所述减振器在三轴方向上各自的振动传递率矩阵
Figure 780942DEST_PATH_IMAGE008
,
Figure 133426DEST_PATH_IMAGE009
,
步骤502、计算矩阵
Figure 810395DEST_PATH_IMAGE010
中所有元素的均值
Figure 248330DEST_PATH_IMAGE011
步骤503、当N个所述减振器在三轴方向上各自的振动传递均满足不等式
Figure 934526DEST_PATH_IMAGE012
时,则表明当前测试的金属丝网减振器组在三轴方向上均满足一致性要求,否则不满足一致性要求;其中,i=1,2,...,N;η为误差允许值,0≤η≤0.1。
上述的方法,其特征在于,步骤六中进行各减振器的动刚度计算的具体步骤为:
步骤601、根据公式
Figure 705298DEST_PATH_IMAGE013
,计算在三轴方向上的解耦频响函数矩阵
Figure 818748DEST_PATH_IMAGE014
步骤602、根据公式
Figure 743979DEST_PATH_IMAGE015
,计算编号为p的被动端区对应的所述减振器在三轴方向上的动刚度矩阵
Figure 968287DEST_PATH_IMAGE016
,其中,
Figure 295363DEST_PATH_IMAGE017
Figure 579713DEST_PATH_IMAGE018
的广义逆矩阵。
上述的方法,其特征在于:步骤七中根据计算出的第p个被动端区对应的减振器在三轴方向上的动刚度矩阵
Figure 788978DEST_PATH_IMAGE019
评估该减振器在三轴方向上的动刚度设计矩阵
Figure 816977DEST_PATH_IMAGE020
,具体步骤为:
S1、动刚度设计矩阵
Figure 998559DEST_PATH_IMAGE021
为对角矩阵,将
Figure 453812DEST_PATH_IMAGE021
中为0的元素的位置集合记为
Figure 353634DEST_PATH_IMAGE022
,将
Figure 919745DEST_PATH_IMAGE023
中位置为
Figure 18151DEST_PATH_IMAGE024
的元素设置为0后得到
Figure 644304DEST_PATH_IMAGE025
,
Figure 765844DEST_PATH_IMAGE025
也为对角矩阵;
S2、提取对角矩阵
Figure 135646DEST_PATH_IMAGE021
中的对角元素的值记为
Figure 291820DEST_PATH_IMAGE026
,提取对角矩阵
Figure 151192DEST_PATH_IMAGE025
中的对角元素的值记为
Figure 760028DEST_PATH_IMAGE027
,当不等式
Figure 667941DEST_PATH_IMAGE028
成立时,该减振器在三轴方向上的动刚度设计矩阵
Figure 678622DEST_PATH_IMAGE021
合理,否则不合理。
本发明与现有技术相比的区别特征为:本发明通过在底座板和测试板之间安装多个减振器,再利用简单的频响函数测试,即可同时测试多个减振器,使用的装置结构简单,安装拆卸方便,测试速度快,测试数据处理方式统一,可大幅提高减振器动刚度和一致性测试及动刚度设计矩阵的评估效率,有利于满足一致性要求的减振器产品应用于飞机构件,提高系统的动力学性能一致性。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1为本发明的方法流程框图。
图2为本发明使用的装置的结构示意图。
图3为本发明使用的装置中被动端区的位置示意图。
附图标记说明:
1—底座板; 2—测试板; 3—减振器主体;
4—减振器安装座; 5—螺母; 6—凸台;
7—被动端区; 8—减振器安装轴;
10—第一三轴加速度传感器; 11—第二三轴加速度传感器;
12—中心三轴加速度传感器。
具体实施方式
如图1、图2和图3所示,本发明的一种飞机减振构件的动刚度和一致性测试方法,所述飞机减振构件包括金属丝网减振器组,该方法采用的装置包括用于安装金属丝网减振器组的底座板1和连接在所述金属丝网减振器组顶部的测试板2;金属丝网减振器组包括多个减振器,所述减振器包括减振器主体3、减振器安装座4和减振器安装轴8;
该方法包括以下步骤:
步骤一、金属丝网减振器组的安装:将N个所述减振器沿底座板1的上表面周侧均匀安装,减振器安装轴8的顶部穿过测试板2并与测试板2连接;N为不小于3的正整数;
本实施例中,底座板1为正N边形板或圆形板,所述测试板2为金属板。
本实施例中,所述减振器安装座4连接在底座板1上,减振器安装轴8的顶部穿过测试板2并通过螺母5与测试板2连接;
本实施例中,在实际安装时,底座板1需要固定至刚性地面或刚度较大的夹具平台上。
步骤二、测试区域的划分:
将减振器安装座4的上表面记做主动端区;将测试板2上表面位于减振器安装轴8周侧的环形区域记做被动端区7,减振器安装轴8的轴线穿过被动端区7的环心,被动端区7的环宽为10mm~15mm;
需要说明的是,所述主动端区和被动端区7组成所述测试区域;
本实施例中,所述被动端区7的内径等于螺母5的外径;
需要说明的是,所述主动端区的数量和被动端区7的数量与所述减振器的数量相同且一一对应;
需要说明的是,所述主动端区为直接受到减振作用的区域,所述被动端区7为非直接受到减振作用的区域;
步骤三、三轴加速度传感器组的安装:
在各所述被动端区7内均粘贴一个第一三轴加速度传感器10,在各所述主动端区内均粘贴一个第二三轴加速度传感器11,在测试板2上表面的几何中心粘贴中心三轴加速度传感器12,将粘贴中心三轴加速度传感器12的区域记做中心测试区,第二三轴加速度传感器11、第一三轴加速度传感器10和中心三轴加速度传感器12的轴向方向均平行、切向方向均平行、径向方向均平行;
优选的,N=4,所述底座板1为正方形板,中心三轴加速度传感器12的轴向方向垂直于底座板1,中心三轴加速度传感器12的切向方向与底座板1上表面的一条边平行。
本实施例中,多个第一三轴加速度传感器10组成第一三轴加速度传感器单元,多个第二三轴加速度传感器11组成第二三轴加速度传感器单元;第一三轴加速度传感器单元、第二三轴加速度传感器单元和中心三轴加速度传感器12组成三轴加速度传感器组;
步骤四、测量各测试区域的频响函数:
以第二三轴加速度传感器11的安装中心点作为主动端区的响应点,以第一三轴加速度传感器10的安装中心点作为被动端区7的响应点,以中心三轴加速度传感器12的安装中心点作为中心测试区的响应点;测量主动端区的响应点、被动端区7的响应点和中心测试区的响应点之间的频响函数,组成三轴方向上的耦合频响函数矩阵;
步骤五、金属丝网减振器组的一致性判断:根据步骤四中获取的耦合频响函数矩阵,首先分别计算各个减振器在三轴方向上的振动传递率,随后分别计算N个减振器在三轴方向上的振动传递率均值;比较各轴方向上各个减振器的振动传递率和N个减振器的振动传递率均值之间的误差,当三轴方向上误差均小于等于误差允许值时,则表明当前测试的金属丝网减振器组满足一致性要求,否则不满足一致性要求;
步骤六、各减振器的动刚度计算:根据步骤四中获取的三轴方向上的耦合频响函数矩阵计算三轴方向上的解耦频响函数矩阵,根据三轴方向上的解耦频响函数矩阵和各个减振器在三轴方向上的振动传递率矩阵计算各个减振器在三轴方向上的动刚度矩阵;
步骤七、动刚度设计矩阵的评估:根据各个减振器在三轴方向上的动刚度矩阵评估各个减振器在三轴方向上的动刚度设计矩阵。
本实施例中,第二三轴加速度传感器11、第一三轴加速度传感器10和中心三轴加速度传感器的型号相同,第二三轴加速度传感器11的最大量程不大于50g、灵敏度10mv/g,第二三轴加速度传感器11的尺寸不超过10mm*10mm,重量不超过10g。
本实施例中,所述减振器安装轴8上具有供测试板2安装的错台,或所述减振器安装轴8为圆台状安装轴,测试板2上开设有供所述圆台状安装轴穿过的通孔,所述通孔的孔径大于所述圆台状安装轴的大径端直径且小于所述圆台状安装轴的小径端直径。目的是使测试板2能够稳定的安装在减振器安装轴8上。
本实施例中,所述测试板2上表面的几何中心设置有用于安装所述中心三轴加速度传感器12的凸台6,所述凸台6的厚度小于测试板2的厚度。
需要说明的是,所述凸台6用于限制中心加速度传感器的粘贴位置和方向,保证中心加速度传感器安装准确。
本实施例中,优选的,所述减振器的数量为4个,所述测试板2的横截面为正方形,所述测试板2与底座板1的横截面尺寸相同,所述测试板2的横截面的边长为60mm~90mm,所述测试板2的厚度为2mm~4mm。
本发明通过在底座板1和测试板2之间安装多个减振器,再利用简单的频响函数测试,即可同时测试多个减振器,整个流程步骤简单,安装拆卸方便,测试速度快,测试数据处理方式统一,可大幅提高减振器动刚度和一致性测试效率,有利于满足一致性要求的减振器产品应用于飞机构件,提高系统的动力学性能一致性。
本实施例中,步骤四的具体步骤为:
步骤401、使用力锤在各主动端区分别施加振动激励,N个第二三轴加速度传感器11接收三轴加速度信号,得到主动端区的响应点在三轴方向上的耦合频响函数矩阵
Figure 646578DEST_PATH_IMAGE001
;其中,a为各主动端区的编号且a=1,2,...,N;A=1,2,3,当A为1时表示沿切向方向,当A为2时表示沿径向方向,当A为3时表示沿轴向方向;
需要说明的是,
Figure 742710DEST_PATH_IMAGE029
为主动端区的响应点在切向方向上的耦合频响函数矩阵,
Figure 29815DEST_PATH_IMAGE030
为主动端区的响应点在径向方向上的耦合频响函数矩阵,
Figure 160582DEST_PATH_IMAGE031
为主动端区的响应点在轴向方向上的耦合频响函数矩阵;
步骤402、使用力锤在各主动端区分别施加振动激励,N个第一三轴加速度传感器10接收三轴加速度信号,得到主动端区的激励点到被动端区7的响应点在三轴方向上的耦合频响函数矩阵
Figure 299439DEST_PATH_IMAGE002
;其中,p为各被动端区7的编号且p=1,2,...,N;
步骤403、使用力锤在各主动端区分别施加振动激励,中心三轴加速度传感器12接收三轴加速度信号,得到主动端区的激励点到中心测试区的响应点在三轴方向上的耦合频响函数矩阵
Figure 617288DEST_PATH_IMAGE003
;其中,t为中心测试区的编号且t=1;
步骤404、使用力锤在各被动端区7分别施加振动激励,N个第一三轴加速度传感器10接收三轴加速度信号,得到被动端区7的响应点在三轴方向上的耦合频响函数矩阵
Figure 132583DEST_PATH_IMAGE004
步骤405、使用力锤在各被动端区7分别施加振动激励,中心三轴加速度传感器12接收三轴加速度信号,得到被动端区7的激励点到中心测试区的响应点在三轴方向上的耦合频响函数矩阵
Figure 180174DEST_PATH_IMAGE005
步骤406、使用力锤在各被动端区7分别施加振动激励,N个第二三轴加速度传感器11接收三轴加速度信号,得到被动端区7的激励点到主动端区的响应点在三轴方向上的耦合频响函数矩阵
Figure 489932DEST_PATH_IMAGE006
需要说明的是,
Figure 295077DEST_PATH_IMAGE032
为主动端区的响应点在三轴方向上的耦合频响函数矩阵,
Figure 348484DEST_PATH_IMAGE033
为主动端区的激励点到被动端区7的响应点在三轴方向上的耦合频响函数矩阵,
Figure 188264DEST_PATH_IMAGE034
为主动端区的激励点到中心测试区的响应点在三轴方向上的耦合频响函数矩阵,
Figure 668924DEST_PATH_IMAGE035
为被动端区7的响应点在三轴方向上的耦合频响函数矩阵,
Figure 289261DEST_PATH_IMAGE036
为被动端区7的激励点到中心测试区的响应点在三轴方向上的耦合频响函数矩阵,
Figure 880779DEST_PATH_IMAGE037
为被动端区7的激励点到主动端区的响应点在三轴方向上的耦合频响函数矩阵,其中X表示矩阵行号,Y表示矩阵列号;本实施例中,X取a,p,t;Y取a,p。
本实施例中,在执行步骤四之前,还需要对第二三轴加速度传感器11、第一三轴加速度传感器10和中心三轴加速度传感器12的灵敏度进行标定,其灵敏度均为10mv/g;
本实施例中,将步骤四中的耦合频响函数矩阵展开为:
Figure 840645DEST_PATH_IMAGE038
Figure 492206DEST_PATH_IMAGE039
Figure 271943DEST_PATH_IMAGE040
Figure 729470DEST_PATH_IMAGE041
Figure 543842DEST_PATH_IMAGE042
Figure 366304DEST_PATH_IMAGE043
本实施例中,考虑被动端区7的响应点的三个自由度,考虑主动端区的响应点的三个自由度;
需要说明的是,三轴加速度传感器组与数据采集装置进行连接,进行数据的采集和处理。
本实施例中,步骤五中进行金属丝网减振器组的一致性判断的具体步骤为:
步骤501、根据公式
Figure 633338DEST_PATH_IMAGE007
,计算N个所述减振器在三轴方向上各自的振动传递率矩阵
Figure 832238DEST_PATH_IMAGE008
,
Figure 501116DEST_PATH_IMAGE009
,步骤502、计算矩阵
Figure 58262DEST_PATH_IMAGE044
中所有元素的均值
Figure 812591DEST_PATH_IMAGE011
步骤503、当N个所述减振器在三轴方向上各自的振动传递均满足不等式
Figure 815182DEST_PATH_IMAGE012
时,则表明当前测试的金属丝网减振器组在三轴方向上均满足一致性要求,否则不满足一致性要求;其中,i=1,2,...,N;η为误差允许值,0≤η≤0.1。
其中,
Figure 338568DEST_PATH_IMAGE045
为编号为1的被动端区7对应的所述减振器在A方向上的振动传递率,
Figure 502833DEST_PATH_IMAGE046
为编号为2的被动端区7对应的所述减振器在A方向上的振动传递率,...,
Figure 744458DEST_PATH_IMAGE047
为编号为N的被动端区7对应的所述减振器在A方向上的振动传递率;
本实施例中,分析频率范围为5Hz~2000Hz,在测试过程中会存在多个共振频率,当前测试的金属丝网减振器组满足一致性要求时,则各共振频率处的振动传递率矩阵
Figure 347478DEST_PATH_IMAGE048
均能使不等式
Figure 725370DEST_PATH_IMAGE049
成立时,当前测试的金属丝网减振器组在三轴方向上均满足一致性要求,否则不满足一致性要求。
本实施例中,步骤六中进行各减振器的动刚度计算的具体步骤为:
步骤601、根据公式
Figure 326115DEST_PATH_IMAGE013
,计算在三轴方向上的解耦频响函数矩阵
Figure 789458DEST_PATH_IMAGE014
步骤602、根据公式
Figure 133851DEST_PATH_IMAGE015
,计算编号为p的被动端区7对应的所述减振器在三轴方向上的动刚度矩阵
Figure 428566DEST_PATH_IMAGE016
,其中,
Figure 200213DEST_PATH_IMAGE017
Figure 150852DEST_PATH_IMAGE018
的广义逆矩阵。
本实施例中,步骤七中根据计算出的第p个被动端区7对应的减振器在三轴方向上的动刚度矩阵
Figure 33357DEST_PATH_IMAGE019
评估该减振器在三轴方向上的动刚度设计矩阵
Figure 385841DEST_PATH_IMAGE020
,具体步骤为:
S1、动刚度设计矩阵
Figure 328389DEST_PATH_IMAGE021
为对角矩阵,将
Figure 828641DEST_PATH_IMAGE021
中为0的元素的位置集合记为
Figure 514837DEST_PATH_IMAGE022
,将
Figure 721827DEST_PATH_IMAGE023
中位置为
Figure 835277DEST_PATH_IMAGE024
的元素设置为0后得到
Figure 760508DEST_PATH_IMAGE025
,
Figure 984816DEST_PATH_IMAGE025
也为对角矩阵;
S2、提取对角矩阵
Figure 875673DEST_PATH_IMAGE021
中的对角元素的值记为
Figure 894445DEST_PATH_IMAGE026
,提取对角矩阵
Figure 306972DEST_PATH_IMAGE025
中的对角元素的值记为
Figure 334971DEST_PATH_IMAGE027
,当不等式
Figure 516553DEST_PATH_IMAGE028
成立时,该减振器在三轴方向上的动刚度设计矩阵
Figure 34122DEST_PATH_IMAGE021
合理,否则不合理。
本实施例中,当该减振器的动刚度设计矩阵
Figure 668366DEST_PATH_IMAGE021
不合理时对其进行改进,随后根据改进后的动刚度设计矩阵
Figure 234477DEST_PATH_IMAGE021
进行金属丝网减振器的加工生产,重复评估当前动刚度设计矩阵是否合理,实现快速迭代,便于优化动刚度设计矩阵。
本实施例中,当减振器的动刚度设计矩阵
Figure 536145DEST_PATH_IMAGE021
合理时,则各共振频率处的
Figure 162298DEST_PATH_IMAGE021
Figure 283838DEST_PATH_IMAGE025
均能使
Figure 715956DEST_PATH_IMAGE050
成立。
以上所述,仅是本发明的较佳实施例,并非对本发明作任何限制,凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化,均仍属于本发明技术方案的保护范围内。

Claims (5)

1.一种飞机减振构件的动刚度和一致性测试方法,所述飞机减振构件包括金属丝网减振器组,该方法采用的装置包括用于安装金属丝网减振器组的底座板(1)和连接在所述金属丝网减振器组顶部的测试板(2);金属丝网减振器组包括多个减振器,所述减振器包括减振器主体(3)、减振器安装座(4)和减振器安装轴(8);
其特征在于:该方法包括以下步骤:
步骤一、金属丝网减振器组的安装:将N个所述减振器沿底座板(1)的上表面周侧均匀安装,减振器安装轴(8)的顶部穿过测试板(2)并与测试板(2)连接;N为不小于3的正整数;
步骤二、测试区域的划分:
将减振器安装座(4)的上表面记做主动端区;将测试板(2)上表面位于减振器安装轴(8)周侧的环形区域记做被动端区(7),减振器安装轴(8)的轴线穿过被动端区(7)的环心,被动端区(7)的环宽为10mm~15mm;
步骤三、三轴加速度传感器组的安装:
在各所述被动端区(7)内均粘贴一个第一三轴加速度传感器(10),在各所述主动端区内均粘贴一个第二三轴加速度传感器(11),在测试板(2)上表面的几何中心粘贴中心三轴加速度传感器(12),将粘贴中心三轴加速度传感器(12)的区域记做中心测试区,第二三轴加速度传感器(11)、第一三轴加速度传感器(10)和中心三轴加速度传感器(12)的轴向方向均平行、切向方向均平行、径向方向均平行;
步骤四、测量各测试区域的频响函数:
以第二三轴加速度传感器(11)的安装中心点作为主动端区的响应点,以第一三轴加速度传感器(10)的安装中心点作为被动端区(7)的响应点,以中心三轴加速度传感器(12)的安装中心点作为中心测试区的响应点;测量主动端区的响应点、被动端区(7)的响应点和中心测试区的响应点之间的频响函数,组成三轴方向上的耦合频响函数矩阵;
步骤五、金属丝网减振器组的一致性判断:根据步骤四中获取的耦合频响函数矩阵,首先分别计算各个减振器在三轴方向上的振动传递率,随后分别计算N个减振器在三轴方向上的振动传递率均值;比较各轴方向上各个减振器的振动传递率和N个减振器的振动传递率均值之间的误差,当三轴方向上误差均小于等于误差允许值时,则表明当前测试的金属丝网减振器组满足一致性要求,否则不满足一致性要求;
步骤六、各减振器的动刚度计算:根据步骤四中获取的三轴方向上的耦合频响函数矩阵计算三轴方向上的解耦频响函数矩阵,根据三轴方向上的解耦频响函数矩阵和各个减振器在三轴方向上的振动传递率矩阵计算各个减振器在三轴方向上的动刚度矩阵;
步骤七、动刚度设计矩阵的评估:根据各个减振器在三轴方向上的动刚度矩阵评估各个减振器在三轴方向上的动刚度设计矩阵。
2.根据权利要求1所述的方法,其特征在于,步骤四的具体步骤为:
步骤401、使用力锤在各主动端区分别施加振动激励,N个第二三轴加速度传感器(11)接收三轴加速度信号,得到主动端区的响应点在三轴方向上的耦合频响函数矩阵
Figure 220336DEST_PATH_IMAGE001
;其中,a为各主动端区的编号且a=1,2,...,N;A=1,2,3,当A为1时表示沿切向方向,当A为2时表示沿径向方向,当A为3时表示沿轴向方向;
步骤402、使用力锤在各主动端区分别施加振动激励,N个第一三轴加速度传感器(10)接收三轴加速度信号,得到主动端区的激励点到被动端区(7)的响应点在三轴方向上的耦合频响函数矩阵
Figure 231017DEST_PATH_IMAGE002
;其中,p为各被动端区(7)的编号且p=1,2,...,N;
步骤403、使用力锤在各主动端区分别施加振动激励,中心三轴加速度传感器(12)接收三轴加速度信号,得到主动端区的激励点到中心测试区的响应点在三轴方向上的耦合频响函数矩阵
Figure 198973DEST_PATH_IMAGE003
;其中,t为中心测试区的编号且t=1;
步骤404、使用力锤在各被动端区(7)分别施加振动激励,N个第一三轴加速度传感器(10)接收三轴加速度信号,得到被动端区(7)的响应点在三轴方向上的耦合频响函数矩阵
Figure 295105DEST_PATH_IMAGE004
步骤405、使用力锤在各被动端区(7)分别施加振动激励,中心三轴加速度传感器(12)接收三轴加速度信号,得到被动端区(7)的激励点到中心测试区的响应点在三轴方向上的耦合频响函数矩阵
Figure 803447DEST_PATH_IMAGE005
步骤406、使用力锤在各被动端区(7)分别施加振动激励,N个第二三轴加速度传感器(11)接收三轴加速度信号,得到被动端区(7)的激励点到主动端区的响应点在三轴方向上的耦合频响函数矩阵
Figure 934214DEST_PATH_IMAGE006
3.根据权利要求2所述的方法,其特征在于,步骤五中进行金属丝网减振器组的一致性判断的具体步骤为:
步骤501、根据公式
Figure 73071DEST_PATH_IMAGE007
,计算N个所述减振器在三轴方向上各自的振动传递率矩阵
Figure 656499DEST_PATH_IMAGE008
,
Figure 906215DEST_PATH_IMAGE009
,步骤502、计算矩阵
Figure 891489DEST_PATH_IMAGE010
中所有元素的均值
Figure 765029DEST_PATH_IMAGE011
步骤503、当N个所述减振器在三轴方向上各自的振动传递均满足不等式
Figure 835753DEST_PATH_IMAGE012
时,则表明当前测试的金属丝网减振器组在三轴方向上均满足一致性要求,否则不满足一致性要求;其中,i=1,2,...,N;η为误差允许值,0≤η≤0.1。
4.根据权利要求3所述的方法,其特征在于,步骤六中进行各减振器的动刚度计算的具体步骤为:
步骤601、根据公式
Figure 889160DEST_PATH_IMAGE013
,计算在三轴方向上的解耦频响函数矩阵
Figure 728940DEST_PATH_IMAGE014
步骤602、根据公式
Figure 209600DEST_PATH_IMAGE015
,计算编号为p的被动端区(7)对应的所述减振器在三轴方向上的动刚度矩阵
Figure 502041DEST_PATH_IMAGE016
,其中,
Figure 421455DEST_PATH_IMAGE017
Figure 115742DEST_PATH_IMAGE018
的广义逆矩阵。
5.根据权利要求4所述的方法,其特征在于:步骤七中根据计算出的第p个被动端区(7)对应的减振器在三轴方向上的动刚度矩阵
Figure 767303DEST_PATH_IMAGE019
评估该减振器在三轴方向上的动刚度设计矩阵
Figure 812620DEST_PATH_IMAGE020
,具体步骤为:
S1、动刚度设计矩阵
Figure 207829DEST_PATH_IMAGE021
为对角矩阵,将
Figure 818939DEST_PATH_IMAGE021
中为0的元素的位置集合记为
Figure 906980DEST_PATH_IMAGE022
,将
Figure 174014DEST_PATH_IMAGE023
中位置为
Figure 107335DEST_PATH_IMAGE022
的元素设置为0后得到
Figure 776213DEST_PATH_IMAGE024
,
Figure 35156DEST_PATH_IMAGE024
也为对角矩阵;
S2、提取对角矩阵
Figure 586223DEST_PATH_IMAGE021
中的对角元素的值记为
Figure 588814DEST_PATH_IMAGE025
,提取对角矩阵
Figure 112200DEST_PATH_IMAGE024
中的对角元素的值记为
Figure 276465DEST_PATH_IMAGE026
,当不等式
Figure 518090DEST_PATH_IMAGE027
成立时,该减振器在三轴方向上的动刚度设计矩阵
Figure 622575DEST_PATH_IMAGE021
合理,否则不合理。
CN202210686028.XA 2022-06-17 2022-06-17 一种飞机减振构件的动刚度和一致性测试方法 Active CN114778048B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210686028.XA CN114778048B (zh) 2022-06-17 2022-06-17 一种飞机减振构件的动刚度和一致性测试方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210686028.XA CN114778048B (zh) 2022-06-17 2022-06-17 一种飞机减振构件的动刚度和一致性测试方法

Publications (2)

Publication Number Publication Date
CN114778048A true CN114778048A (zh) 2022-07-22
CN114778048B CN114778048B (zh) 2022-09-02

Family

ID=82420640

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210686028.XA Active CN114778048B (zh) 2022-06-17 2022-06-17 一种飞机减振构件的动刚度和一致性测试方法

Country Status (1)

Country Link
CN (1) CN114778048B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115750647A (zh) * 2022-11-03 2023-03-07 中国海洋大学 一种用于轨道交通的金属丝网减振器、设计方法及系统

Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3693402A (en) * 1971-07-19 1972-09-26 Caterpillar Tractor Co Vibration damper test machine
JP2001091398A (ja) * 1999-09-24 2001-04-06 Takenaka Komuten Co Ltd 免震部材の特性試験方法及び特性試験装置
CN102607788A (zh) * 2012-03-27 2012-07-25 西安航空技术高等专科学校 超大尺度光学平台隔振性能一致性检测装置与检测方法
CN103308333A (zh) * 2013-05-17 2013-09-18 同济大学 一种隔振器动刚度测试方法
CN103808499A (zh) * 2014-01-21 2014-05-21 江苏科技大学 一种隔振器动刚度测试方法及其装置
CN104897354A (zh) * 2015-06-23 2015-09-09 中国航空工业集团公司西安飞机设计研究所 一种飞机发动机隔振器刚度及阻尼的测试装置及测试方法
CN105115690A (zh) * 2015-08-18 2015-12-02 上海交通大学 一种隔振器多向阻抗矩阵及刚度测试试验装置和试验方法
CN105243172A (zh) * 2015-07-22 2016-01-13 湖南大学 一种考虑动态刚度和阻尼的发动机悬置系统优化方法
RU2607361C1 (ru) * 2015-08-10 2017-01-10 Олег Савельевич Кочетов Способ испытаний многомассовых систем виброизоляции
CN106840559A (zh) * 2015-12-03 2017-06-13 株洲时代新材料科技股份有限公司 一种可调式单层及双层隔振系统振动试验台及试验方法
CN108170643A (zh) * 2017-12-26 2018-06-15 上海交通大学 基于原位测量频响函数的免拆分子结构频响函数识别方法
CN108827457A (zh) * 2018-06-12 2018-11-16 太原理工大学 基于传递函数的振动幅值迭代控制方法
CN109557186A (zh) * 2018-12-24 2019-04-02 暨南大学 一种机械装配动态质量检测方法
CN110210179A (zh) * 2019-06-28 2019-09-06 合肥工业大学 一种弹性元件动刚度及其两端载荷的获取方法
CN111766030A (zh) * 2020-07-03 2020-10-13 中国飞机强度研究所 一种用于柔性结构的模态测试装置及模态测试方法
US20210172829A1 (en) * 2019-12-04 2021-06-10 Pukyong National University Industry-University Cooperation Foundation Device for analyzing sensitivity of object using frequency response and analyzing method using the same
CN215785038U (zh) * 2021-09-10 2022-02-11 中国航空工业集团公司西安飞行自动控制研究所 一种减振器刚度一致性的筛选工装
CN114264460A (zh) * 2021-12-15 2022-04-01 渤海造船厂集团有限公司 一种橡胶减振器动刚度测试装置及安装方法

Patent Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3693402A (en) * 1971-07-19 1972-09-26 Caterpillar Tractor Co Vibration damper test machine
JP2001091398A (ja) * 1999-09-24 2001-04-06 Takenaka Komuten Co Ltd 免震部材の特性試験方法及び特性試験装置
CN102607788A (zh) * 2012-03-27 2012-07-25 西安航空技术高等专科学校 超大尺度光学平台隔振性能一致性检测装置与检测方法
CN103308333A (zh) * 2013-05-17 2013-09-18 同济大学 一种隔振器动刚度测试方法
CN103808499A (zh) * 2014-01-21 2014-05-21 江苏科技大学 一种隔振器动刚度测试方法及其装置
CN104897354A (zh) * 2015-06-23 2015-09-09 中国航空工业集团公司西安飞机设计研究所 一种飞机发动机隔振器刚度及阻尼的测试装置及测试方法
CN105243172A (zh) * 2015-07-22 2016-01-13 湖南大学 一种考虑动态刚度和阻尼的发动机悬置系统优化方法
RU2607361C1 (ru) * 2015-08-10 2017-01-10 Олег Савельевич Кочетов Способ испытаний многомассовых систем виброизоляции
CN105115690A (zh) * 2015-08-18 2015-12-02 上海交通大学 一种隔振器多向阻抗矩阵及刚度测试试验装置和试验方法
CN106840559A (zh) * 2015-12-03 2017-06-13 株洲时代新材料科技股份有限公司 一种可调式单层及双层隔振系统振动试验台及试验方法
CN108170643A (zh) * 2017-12-26 2018-06-15 上海交通大学 基于原位测量频响函数的免拆分子结构频响函数识别方法
CN108827457A (zh) * 2018-06-12 2018-11-16 太原理工大学 基于传递函数的振动幅值迭代控制方法
CN109557186A (zh) * 2018-12-24 2019-04-02 暨南大学 一种机械装配动态质量检测方法
CN110210179A (zh) * 2019-06-28 2019-09-06 合肥工业大学 一种弹性元件动刚度及其两端载荷的获取方法
US20210172829A1 (en) * 2019-12-04 2021-06-10 Pukyong National University Industry-University Cooperation Foundation Device for analyzing sensitivity of object using frequency response and analyzing method using the same
CN111766030A (zh) * 2020-07-03 2020-10-13 中国飞机强度研究所 一种用于柔性结构的模态测试装置及模态测试方法
CN215785038U (zh) * 2021-09-10 2022-02-11 中国航空工业集团公司西安飞行自动控制研究所 一种减振器刚度一致性的筛选工装
CN114264460A (zh) * 2021-12-15 2022-04-01 渤海造船厂集团有限公司 一种橡胶减振器动刚度测试装置及安装方法

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
LIN CHANG-LIANG 等: "Adjustable stiffness structure application in helicopter vibration reduction", 《JOURNAL OF VIBRATION ENGINEERING》 *
MALY, JR 等: "Complex stiffness measurement of vibration-damped structural elements", 《PROCEEDINGS OF THE SOCIETY OF PHOTO-OPTICAL INSTRUMENTATION ENGINEERS (SPIE)》 *
任佳 等: "装甲车座椅抗地雷爆炸冲击模拟试验方法研究", 《应用力学学报》 *
任明可 等: "新型橡胶-电磁复合主被动隔振器研究", 《振动与冲击》 *
段宇星 等: "无人机捷联惯组隔振系统动力学分析与优化设计", 《中国惯性技术学报》 *
沈良杰 等: "复合隔振器隔振装置的性能研究", 《装备环境工程》 *
董万元 等: "大功率涡桨发动机隔振系统设计与试验", 《振动、测试与诊断》 *
董万元 等: "金属丝网非线性隔振器的动静态特性研究", 《实验力学》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115750647A (zh) * 2022-11-03 2023-03-07 中国海洋大学 一种用于轨道交通的金属丝网减振器、设计方法及系统

Also Published As

Publication number Publication date
CN114778048B (zh) 2022-09-02

Similar Documents

Publication Publication Date Title
CN104198197A (zh) 一种车身接附点动刚度测试方法和装置
CN114778048B (zh) 一种飞机减振构件的动刚度和一致性测试方法
CN109632159B (zh) 一种六维力和力矩传感器
CN106468623A (zh) 一种整车状态下动力总成悬置系统刚体模态参数测试方法
CN201306407Y (zh) 动力吸振器
CN110895195A (zh) 一种前下摆臂台架耐久试验载荷获取方法
CN111241706A (zh) 基于动态k&c试验系统的汽车悬架动载性能评价方法
WO2023024968A1 (zh) 一种变行程自适应调整准零刚度装置及参数校核方法
CN113340526B (zh) 一种六维力传感器静、动态标定装置及其标定方法
CN112611579B (zh) 一种两通道台架耐久测试方法
CN114689265B (zh) 用于弹性机翼气动弹性风洞试验中结构弯矩动态标定方法
CN112595528A (zh) 车辆行驶工况下动力总成悬置主动侧激励力的识别方法
CN115186377A (zh) 一种商用车悬置系统的稳健性优化方法
CN110210179A (zh) 一种弹性元件动刚度及其两端载荷的获取方法
CN114778047B (zh) 一种飞机减振构件的动刚度和一致性测试装置
CN104156550A (zh) 车辆钢板弹簧悬架系统阻尼比的分析计算方法
CN111114225B (zh) 一种基于阵列式减振器的车用悬架控制臂及参数计算方法
CN111780993B (zh) 一种车轮力测量系统及方法
US20060075810A1 (en) Devices and methods for simulating tire non-uniformity forces for vehicle vibration sensitivity measurements and tuning
CN115081116A (zh) 车辆底盘零部件最大应力计算方法、系统及存储介质
CN110160736B (zh) 一种耦合弹性模态非定常气动力测量装置及方法
CN113155452A (zh) 一种多节传动轴的模态测试装置及方法
CN219265650U (zh) 一种飞机apu减振系统振动试验装置
CN116929702B (zh) 一种飞行器气动特性测试天平装置与测试方法
CN210293701U (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