CN113916432A - 拉索的索力测量方法 - Google Patents

拉索的索力测量方法 Download PDF

Info

Publication number
CN113916432A
CN113916432A CN202111060435.1A CN202111060435A CN113916432A CN 113916432 A CN113916432 A CN 113916432A CN 202111060435 A CN202111060435 A CN 202111060435A CN 113916432 A CN113916432 A CN 113916432A
Authority
CN
China
Prior art keywords
cable
value
cable force
vibration frequency
vibration
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.)
Pending
Application number
CN202111060435.1A
Other languages
English (en)
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.)
Shanghai Construction Group Co Ltd
Original Assignee
Shanghai Construction Group Co Ltd
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 Shanghai Construction Group Co Ltd filed Critical Shanghai Construction Group Co Ltd
Priority to CN202111060435.1A priority Critical patent/CN113916432A/zh
Publication of CN113916432A publication Critical patent/CN113916432A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L5/00—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
    • G01L5/04—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes for measuring tension in flexible members, e.g. ropes, cables, wires, threads, belts or bands
    • G01L5/042—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes for measuring tension in flexible members, e.g. ropes, cables, wires, threads, belts or bands by measuring vibrational characteristics of the flexible member

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Force Measurement Appropriate To Specific Purposes (AREA)

Abstract

本发明公开一种拉索的索力测量方法，包括：在已建或者在建的拉索中，预估拉索的索力预估值，代入以下索力理论公式F＝af2+bf+d计算出索力的各阶振动频率估计值；以振动频率估计值的倍数作为振动法测量拉索的采样频率；根据采样频率对拉索进行振动加速度测量，获取不同采样频率下的加速度时程数据；将各组采样频率下测量获取的加速度时程数据分别进行移动傅里叶变换，计算得出各组振动频率实测值；将多阶振动频率实测值代入索力理论公式，计算得出拉索的多组索力计算值；将多组索力计算值进行失真剔除，获取置信度较高的索力实测值。本发明能够提高拉索的索力测量的准确性。

Description

拉索的索力测量方法
技术领域
本发明涉及建筑工程技术领域，特别涉及一种拉索的索力测量方法。
背景技术
拉索是索结构的关键受力构件，通过自身索体预应力的大小实现整体索结构内力分布、保证结构稳定性、传递结构荷载、控制结构线形形态等作用。因此，无论在索结构施工过程阶段还是运营使用阶段，准确的测量以及控制拉索的索力大小对索结构的整体安全性和稳定性都至关重要。目前，拉索的索力测量方法主要包括：压力传感器测试法、磁通量法、振动法。其中压力传感器测试法：通过在拉索锚固端设置压力传感器，将压力传感器所测压力作为拉索的索力，该方法需在建造过程中事前将压力传感器安装至拉索锚固端，无法测量已有索结构的拉索索力。其中磁通量法：在拉索索体上安装电磁传感器，通过拉索受拉后索体横截面积变化引起的磁通量变化来推算拉索的索力。该方法同样需事前安装磁通量传感器，仅适用于在建索结构的拉索索力测量。其中振动法：首先利用加速度传感器采集测量拉索在环境激励或人工激励下的自振频率，然后通过自振频率与索力的关系计算拉索的索力，该方法可适用于在建和已建的索结构，应用范围较广，但该方法的测量精度受限于理论公式的准确性以及拉索边界条件的影响。
发明内容
本发明所要解决的技术问题是，提供一种拉索的索力测量方法，以解决测量准确性较差的问题。
为了解决上述技术问题，本发明提供的技术方案是：一种拉索的索力测量方法，包括：
步骤S10，在已建或者在建的拉索中，预估拉索的索力预估值，代入以下索力理论公式(1)计算出索力的各阶振动频率估计值；
F＝af2+bf+d (1)；
其中，F为索力值，f为多阶振动频率，a、b、d为常数系数；
步骤S20，以振动频率估计值的倍数作为振动法测量拉索的采样频率；
步骤S30，根据采样频率对拉索进行振动加速度测量，获取不同采样频率下的加速度时程数据；
步骤S40，将各组采样频率下测量获取的加速度时程数据分别进行移动傅里叶变换，计算得出各组振动频率实测值；
步骤S50，将多阶振动频率实测值代入索力理论公式(1)，计算得出拉索的多组索力计算值。
步骤S60，将多组索力计算值进行失真剔除，获取置信度较高的索力实测值。其中索力是指拉力。
进一步地，本发明提供的拉索的索力测量方法，在步骤S40中，通过振动频率实测值的离散度来剔除失真的振动频率实测值，获取置信度较高的振动频率实测值。
进一步地，本发明提供的拉索的索力测量方法，在步骤S60中，所述将多组索力计算值进行失真剔除获取置信度较高的索力实测值的方法包括：
计算多组索力计算值的均方差σ1，计算每个索力计算值与均值的绝对偏差δs(i) j，将绝对偏差δs(i) j＞3σ1的索力计算值认定为置信度较低，将该置信度较低的索力计算值作为失真值剔除；计算失真值剔除完毕后的多组索力计算值的作为置信度较高的索力实测值。
进一步地，本发明提供的拉索的索力测量方法，每组所述加速度时程数据进行移动傅里叶变换方法为：
首先，获取加速度时程数据的序列数量p，自定义移动次数C和移动间隔时间t，确定傅里叶变换计算的数据序列长度l＝p-Ct；
然后，依次针对加速度时程数据的第1～l+1个、第1+t～l+1+t个、第1+2t～l+1+2t个……第1+Ct～l+1+Ct个数据进行傅里叶变换，共计C+1次；
最后，统计每次傅里叶变换计算获取的振动频率值，同时统计同一个振动频率值或者相近振动频率值的重现次数，以重现次数作为振动频率值的置信度；当重现次数小于预定值时，则该振动频率值的置信度低，确定该振动频率值为外界噪声干扰信号，判定该振动频率值作为失真值，当重现次数大于或者等于预定值时判定该振动频率值的置信度越高，该振动频率值作为可靠的振动频率实测值。
进一步地，本发明提供的拉索的索力测量方法，根据采样频率对拉索进行振动加速度测量，测量方向选定x向或者y向。
进一步地，本发明提供的拉索的索力测量方法，所述步骤S10中的索力理论公式(1)的确定方法为：
步骤S11，通过有限元分析软件结合拉索的规格尺寸、锚固边界刚度和减震器边界刚度建立拉索的梁振动模型，以梁振动模型模拟拉索；
步骤S12，在拉索的梁振动模型上施加多个索力值，计算拉索在多个索力值作用下对应的多阶振动频率理论值；
步骤S13，根据索力值和对应的多阶振动频率理论值，拟合出各阶振动频率与其相对应的索力值之间的索力理论公式：
F＝af2+bf+d (1)；
其中，F为索力值，f为多阶振动频率，a、b、d为常数系数。
进一步地，本发明提供的拉索的索力测量方法，所述索力理论公式(1)中的常数系数的数值的确定方法为：将索力值及其对应的振动频率理论值代入索力理论公式(1)，计算常数系数的数值。
进一步地，本发明提供的拉索的索力测量方法，根据振动向量将多阶振动频率理论值分类为：x向多阶振动频率理论值和y向多阶振动频率。
进一步地，本发明提供的拉索的索力测量方法，所述振动向量是指振动方向。
进一步地，本发明提供的拉索的索力测量方法，所述步骤S20中，以振动频率估计值的倍数为10倍作为振动法测量拉索的采样频率。
与现有技术相比，本发明的有益效果如下：
本发明提供的拉索的索力测量方法，根据拉索的索力预估值通过索力理论公式F＝af2+bf+d计算出准确性较高的索力的各阶振动频率估计值并以振动频率估计值的倍数作为振动法测量拉索的索力的采样频率，获取不同采样频率下的加速度时程数据并对加速度时程数据进行移动傅里叶变换，计算得出准确性较高的各组振动频率实测值；然后将多阶振动频率实测值代入索力理论公式F＝af2+bf+d计算得出准确性较高的拉索的多组索力计算值，将组索力计算值进行失真剔除，获取置信度较高的索力实测值，从而提高了振动法对拉索的索力测量的准确性。即通过拉索的多阶振动频率来计算拉索的索力，提高了振动法测索力的准确性。
具体实施方式
本发明实施例提供一种拉索的索力测量方法，可以包括以下步骤：
步骤S10，在已建或者在建的拉索中，预估拉索的索力预估值Fg，代入以下索力理论公式(1)计算出索力的各阶振动频率估计值fg(1，2，3…m)，其中m为正整数；
F＝af2+bf+d (1)；
其中，F为索力值，f为多阶振动频率，a、b、d为常数系数。本步骤S10中，将索力预估值Fg代入索力理论公式(1)中根据a、b、d为已知的常数系数求得各阶振动频率估计值fg(1 ，2，3…m)，此时公式中的f等于各阶振动频率估计值fg(1，2，3…m)。
步骤S20，以振动频率估计值fg的倍数作为振动法测量拉索的采样频率fc；其中振动频率估计值fg倍数可以是2-16倍，例如：10倍。其目的是获得准确值较高的采样频率fc。当然倍数可以根据需要进行调整。
步骤S30，根据采样频率fc对拉索进行振动加速度测量，获取不同采样频率fc下的加速度时程数据。本步骤S30中的加速度里程数据测量为振动法测量拉索的公知技术，此处不作具体阐述。其中对拉索进行振动加速度测量的测量方向选定x向或者y向。
步骤S40，将各组采样频率fc下测量获取的加速度时程数据分别进行移动傅里叶变换，计算得出各组振动频率实测值fs(1，2，3…m)1，2，3…n，其中m、n为正整数；其中各组振动频率实测值fs(1，2，3…m1)1，2，3…n即为多阶振动频率实测值fs(1，2，3…m)1，2，3…n。在本步骤S40中，为了提高准确性较高的振动频率实测值，可以通过振动频率实测值fs(1，2,3…m)1，2，3…n的离散度R1，2，3…n来剔除失真的振动频率实测值，获取置信度较高的振动频率实测值fs(1 ，2，3…m)1，2，3…n。
其中每组所述加速度时程数据进行移动傅里叶变换方法可以为：
首先，获取加速度时程数据的序列数量p，自定义移动次数C和移动间隔时间t，确定傅里叶变换计算的数据序列长度l＝p-Ct。
然后，依次针对加速度时程数据的第1～l+1个、第1+t～l+1+t个、第1+2t～l+1+2t个……第1+Ct～l+1+Ct个数据进行傅里叶变换，共计C+1次。
最后，统计每次傅里叶变换计算获取的振动频率值，同时统计同一个振动频率值或者相近振动频率值的重现次数，以重现次数作为振动频率值的置信度；当重现次数小于预定值时，则该振动频率值的置信度低，确定该振动频率值为外界噪声干扰信号，判定该振动频率值作为失真值，当重现次数大于或者等于预定值时判定该振动频率值的置信度越高，该振动频率值作为可靠的振动频率实测值。
步骤S50，将多阶振动频率实测值fs(1，2，3…m)1，2，3…n代入索力理论公式(1)，计算得出拉索的多组索力计算值Fj(1，2，3…m)1，2，3…n，其中m、n为正整数。本步骤S50中的索力计算值Fj(1，2，3…m)1，2，3…n包括失真值。
步骤S60，将多组索力计算值Fj(1，2，3…m)1，2，3…n进行失真剔除，获取置信度较高的索力实测值Fs(1，2，3…m)1，2，3…n，其中m、n为正整数。
其中步骤S60中，所述将多组索力计算值Fj(1，2，3…m)1，2，3…n进行失真剔除获取置信度较高的索力实测值Fs(1，2，3…m)1，2，3…n的方法可以包括：
计算多组索力计算值Fj(1，2，3…m)1，2，3…n的均方差σ1，计算每个索力计算值Fs(i) j与均值
Figure BDA0003256280830000061
的绝对偏差δs(i) j，将绝对偏差δs(i) j＞3σ1的索力计算值认定为置信度较低，将该置信度较低的索力计算值作为失真值剔除；然后，再次计算失真值剔除后的多组索力计算值Fj(1，2，3…k)1，2，3…k的均方差σ2，其中k为正正数且k<m，再次计算每个索力计算值Fs (i) j与均值
Figure BDA0003256280830000062
的绝对偏差δs(i) j，将绝对偏差δs(i) j＞3σ2的索力计算值认定为置信度较低，将该索力计算值作为失真值剔除……，以此循环，计算失真值剔除完毕后的多组索力计算值的作为置信度较高的索力实测值。为了便于理解，每个索力计算值用Fs(i) j来表示，其中Fs(i) j中的i等于1至m中的具体一个。
为了提供测量的准确性，本发明实施例提供的拉索的索力测量方法，所述步骤S10中的索力理论公式(1)的确定方法为：
步骤S11，通过有限元分析软件结合拉索的规格尺寸、锚固边界刚度和减震器边界刚度建立拉索的梁振动模型，以梁振动模型模拟拉索。由于有限元分析软件没办法用拉索去计算振动频率，所以用梁振动模型来等效拉索。
步骤S12，在拉索的梁振动模型上施加多个索力值，计算拉索在多个索力值作用下对应的多阶振动频率理论值。
步骤S13，根据索力值和对应的多阶振动频率理论值，拟合出各阶振动频率与其相对应的索力值之间的索力理论公式：
F＝af2+bf+d (1)；
其中，F为索力值，f为多阶振动频率，a、b、d为常数系数。
步骤S14，将索力值及其对应的振动频率理论值代入索力理论公式(1)，计算常数系数a、b、d的数值，从而获得常数系数a、b、d的数值，以使索力理论公式(1)能够在后续步骤中进行使用。
为了使索力理论公式(1)在已建或者在建的拉索中能够在x向和y向兼容性测量以计算x向或y向的索力计算值，使公式(1)具有通用性，本发明实施例提供的拉索的索力测量方法，根据振动向量将多阶振动频率理论值分类为：x向多阶振动频率理论值和y向多阶振动频率。其中振动向量是指振动方向。
本发明实施例利用有限元模拟分析结合现场实测的方法，在考虑拉索锚固及减震器的边界刚度影响的情况下，通过拉索多阶振动频率来计算拉索的索力，提高了振动法测索力的准确性。
本发明不限于上述具体实施方式，显然，上述所描述的实施例是本发明实施例的一部分实施例，而不是全部的实施例。基于所描述的本发明的实施例，本领域普通技术人员所获得的所有其他实施例，都属于本发明保护的范围。本领域的技术人员可以对本发明进行其他层次的修改和变动。如此，若本发明的这些修改和变动属于本发明权利要求书的范围之内，则本发明也意图包括这些改动和变动在内。

Claims (10)

1.一种拉索的索力测量方法，其特征在于，包括：
步骤S10，在已建或者在建的拉索中，预估拉索的索力预估值，代入以下索力理论公式(1)计算出索力的各阶振动频率估计值；
F＝af2+bf+d (1)；
其中，F为索力值，f为多阶振动频率，a、b、d为常数系数；
步骤S20，以振动频率估计值的倍数作为振动法测量拉索的采样频率；
步骤S30，根据采样频率对拉索进行振动加速度测量，获取不同采样频率下的加速度时程数据；
步骤S40，将各组采样频率下测量获取的加速度时程数据分别进行移动傅里叶变换，计算得出各组振动频率实测值；
步骤S50，将多阶振动频率实测值代入索力理论公式(1)，计算得出拉索的多组索力计算值。
步骤S60，将多组索力计算值进行失真剔除，获取置信度较高的索力实测值。
2.根据权利要求1所述的拉索的索力测量方法，其特征在于，在步骤S40中，通过振动频率实测值的离散度来剔除失真的振动频率实测值，获取置信度较高的振动频率实测值。
3.根据权利要求1所述的拉索的索力测量方法，其特征在于，在步骤S60中，所述将多组索力计算值进行失真剔除获取置信度较高的索力实测值的方法包括：
计算多组索力计算值的均方差σ1，计算每个索力计算值与均值的绝对偏差δs(i) j，将绝对偏差δs(i) j＞3σ1的索力计算值认定为置信度较低，将该置信度较低的索力计算值作为失真值剔除；计算失真值剔除完毕后的多组索力计算值的作为置信度较高的索力实测值。
4.根据权利要求1所述的拉索的索力测量方法，其特征在于，每组所述加速度时程数据进行移动傅里叶变换方法为：
首先，获取加速度时程数据的序列数量p，自定义移动次数C和移动间隔时间t，确定傅里叶变换计算的数据序列长度l＝p-Ct；
然后，依次针对加速度时程数据的第1～l+1个、第1+t～l+1+t个、第1+2t～l+1+2t个……第1+Ct～l+1+Ct个数据进行傅里叶变换，共计C+1次；
最后，统计每次傅里叶变换计算获取的振动频率值，同时统计同一个振动频率值或者相近振动频率值的重现次数，以重现次数作为振动频率值的置信度；当重现次数小于预定值时，则该振动频率值的置信度低，确定该振动频率值为外界噪声干扰信号，判定该振动频率值作为失真值，当重现次数大于或者等于预定值时判定该振动频率值的置信度越高，该振动频率值作为可靠的振动频率实测值。
5.根据权利要求1所述的拉索的索力测量方法，其特征在于，根据采样频率对拉索进行振动加速度测量，测量方向选定x向或者y向。
6.根据权利要求1所述的拉索的索力测量方法，其特征在于，所述步骤S10中的索力理论公式(1)的确定方法为：
步骤S11，通过有限元分析软件结合拉索的规格尺寸、锚固边界刚度和减震器边界刚度建立拉索的梁振动模型，以梁振动模型模拟拉索；
步骤S12，在拉索的梁振动模型上施加多个索力值，计算拉索在多个索力值作用下对应的多阶振动频率理论值；
步骤S13，根据索力值和对应的多阶振动频率理论值，拟合出各阶振动频率与其相对应的索力值之间的索力理论公式：
F＝af2+bf+d (1)；
其中，F为索力值，f为多阶振动频率，a、b、d为常数系数。
7.根据权利要求6所述的拉索的索力测量方法，其特征在于，所述索力理论公式(1)中的常数系数的数值的确定方法为：将索力值及其对应的振动频率理论值代入索力理论公式(1)，计算常数系数的数值。
8.根据权利要求6所述的拉索的索力测量方法，其特征在于，根据振动向量将多阶振动频率理论值分类为：x向多阶振动频率理论值和y向多阶振动频率。
9.根据权利要求8所述的拉索的索力测量方法，其特征在于，所述振动向量是指振动方向。
10.根据权利要求1所述的拉索的索力测量方法，其特征在于，所述步骤S20中，以振动频率估计值的倍数为10倍作为振动法测量拉索的采样频率。
CN202111060435.1A 2021-09-10 2021-09-10 拉索的索力测量方法 Pending CN113916432A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111060435.1A CN113916432A (zh) 2021-09-10 2021-09-10 拉索的索力测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111060435.1A CN113916432A (zh) 2021-09-10 2021-09-10 拉索的索力测量方法

Publications (1)

Publication Number Publication Date
CN113916432A true CN113916432A (zh) 2022-01-11

Family

ID=79234473

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111060435.1A Pending CN113916432A (zh) 2021-09-10 2021-09-10 拉索的索力测量方法

Country Status (1)

Country Link
CN (1) CN113916432A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116269266A (zh) * 2023-05-22 2023-06-23 广州培生智能科技有限公司 基于ai的老年人健康监测方法和系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102735386A (zh) * 2012-07-14 2012-10-17 福州大学 考虑弯曲刚度的斜拉索索力数值计算方法
CN104933254A (zh) * 2015-06-25 2015-09-23 哈尔滨大金工程试验检测有限公司 基于频比偏移系数的改进线性模型振动法索力测量方法
CN110285909A (zh) * 2019-07-05 2019-09-27 河海大学 基于同步压缩变换的索承桥梁瞬时索力计算方法
CN111928890A (zh) * 2020-07-14 2020-11-13 宁波大学 一种实时测量拉索自振频率与索力的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102735386A (zh) * 2012-07-14 2012-10-17 福州大学 考虑弯曲刚度的斜拉索索力数值计算方法
CN104933254A (zh) * 2015-06-25 2015-09-23 哈尔滨大金工程试验检测有限公司 基于频比偏移系数的改进线性模型振动法索力测量方法
CN110285909A (zh) * 2019-07-05 2019-09-27 河海大学 基于同步压缩变换的索承桥梁瞬时索力计算方法
CN111928890A (zh) * 2020-07-14 2020-11-13 宁波大学 一种实时测量拉索自振频率与索力的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
况中华: "基于振动法的拉索索力测试研究", vol. 43, no. 7, pages 1401 - 1404 *
胡宏祥等: "《环境土壤学》", pages: 119 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116269266A (zh) * 2023-05-22 2023-06-23 广州培生智能科技有限公司 基于ai的老年人健康监测方法和系统
CN116269266B (zh) * 2023-05-22 2023-08-04 广州培生智能科技有限公司 基于ai的老年人健康监测方法和系统

Similar Documents

Publication Publication Date Title
CN103078640B (zh) 一种用于adc的rls自适应滤波校准方法
CN104375976B (zh) 基于张量正则分解的欠定盲源分离中的混合矩阵识别方法
CN110376894B (zh) 一种适用于实时混合试验的两级时滞补偿方法
CN106768574B (zh) 基于磁通量法修正的拉索锚固后线性模型索力测量方法
CN113916432A (zh) 拉索的索力测量方法
CN113820003B (zh) 一种适用于桥梁振动监测的加速度实时重构动态位移方法
CN110987769A (zh) 液体颗粒计数器的标定方法
CN110879302A (zh) 一种石英谐振差动式加速度计温度补偿方法
Lagerblad et al. Dynamic response identification based on state estimation and operational modal analysis
CN106326530A (zh) 一种基于右矩阵分式模型的时变结构模态参数辨识方法
CN110186482B (zh) 一种提高惯性制导航天器的落点精度的方法
CN110348157A (zh) 动态比较器的噪声仿真方法及系统
CN112989966B (zh) 一种改进的模拟电路信号降噪方法
CN115642913A (zh) 模数转换器adc校准方法、装置、设备和存储介质
CN112747773B (zh) 基于Allan方差和随机多项式提高陀螺仪精度的方法
CN111240198B (zh) 压电陶瓷执行器迟滞分析方法
CN113155384A (zh) 用于减小结构阻尼比识别的不确定性的传感器布置方法
CN111308124B (zh) 一种激波管测速传感器时间差确定方法
CN110470925B (zh) 一种基于可拓关联函数的电驱动力总成可靠度测试方法
CN112667952A (zh) 一种结构动态位移非积分重构方法
CN1645067A (zh) 电子天平
CN104200002A (zh) 一种粘性阻尼振动信号中的模态参数提取方法
CN112541157A (zh) 一种信号频率精确估计方法
CN116839780B (zh) 无线电能传输耦合器的微弱电磁力测量装置及方法
CN110186484B (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