CN110761946A - 一种风机塔架振动控制方法及装置 - Google Patents
一种风机塔架振动控制方法及装置 Download PDFInfo
- Publication number
- CN110761946A CN110761946A CN201911007511.5A CN201911007511A CN110761946A CN 110761946 A CN110761946 A CN 110761946A CN 201911007511 A CN201911007511 A CN 201911007511A CN 110761946 A CN110761946 A CN 110761946A
- Authority
- CN
- China
- Prior art keywords
- vibration
- wind turbine
- preset frequency
- tower
- signal
- 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
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000001914 filtration Methods 0.000 claims description 29
- 230000009467 reduction Effects 0.000 claims description 21
- 230000001133 acceleration Effects 0.000 claims description 15
- 230000006870 function Effects 0.000 claims description 12
- 238000012546 transfer Methods 0.000 claims description 12
- 238000013016 damping Methods 0.000 claims description 7
- 230000005284 excitation Effects 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 abstract description 11
- 239000000284 extract Substances 0.000 abstract description 3
- 230000005764 inhibitory process Effects 0.000 abstract description 3
- 238000010248 power generation Methods 0.000 abstract description 3
- 238000012544 monitoring process Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 7
- 230000007246 mechanism Effects 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000000819 phase cycle Methods 0.000 description 3
- 238000000926 separation method Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000000919 ceramic Substances 0.000 description 1
- 230000002860 competitive effect Effects 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D7/00—Controlling wind motors
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/04—Automatic control; Regulation
- F03D7/042—Automatic control; Regulation by means of an electrical or electronic controller
- F03D7/043—Automatic control; Regulation by means of an electrical or electronic controller characterised by the type of control logic
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D7/00—Controlling wind motors
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/022—Adjusting aerodynamic properties of the blades
- F03D7/0224—Adjusting blade pitch
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D7/00—Controlling wind motors
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/0296—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor to prevent, counteract or reduce noise emissions
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2270/00—Control
- F05B2270/30—Control parameters, e.g. input parameters
- F05B2270/334—Vibration measurements
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/72—Wind turbines with rotation axis in wind direction
Abstract
本发明具体涉及一种风机塔架振动控制方法及装置，属于风力发电控制技术领域。方法包括以下步骤：1)获取风机塔架振动信号；2)根据风机塔架振动信号提取至少一个预设频率下的振动能量；3)判断每个预设频率下的振动能量是否超过设定阈值，若超过，则进行该预设频率下的加阻控制。本发明对预设频率下的振动能量进行了提取，在某预设频率的振动能量超过设定阈值时进行该预设频率下的加阻控制，针对性强，负面效果小，对风机塔架振动抑制效果较好。
Description
技术领域
本申请涉及一种风机塔架振动控制方法及装置，属于风力发电控制技术领域。
背景技术
风力发电在中国经历高速发展之后，目前优质风场已经被开发殆尽，风电厂家竞争场地已转移到3类低风速大湍流风场。为提高风机发电量，增加塔架高度成为首选技术手段。风机塔高的增加使风电机组塔架振动控制技术的提升需求越发重要，目前各大风机厂商都在开发柔塔技术，原因是风机塔架高度的增加，导致塔架振动模态进一步降低，塔架一阶振动模态无法避免的与风机叶片穿越频率的1P相交了。这就进一步增加了塔架的振动危险性，风机塔架振动控制技术的提升与升级成为必须进行的课题之一。
近年来对风机塔架振动控制的研究，国内外很多整机商及学者都开展了大量工作。例如授权公告号为CN103321854B的中国专利公开了一种风力发电机组塔架振动的控制方法，该方法利用塔架振动信号为输入，进行时域加阻的方法减小塔架的振动，这对刚性塔架、半刚性塔架都是有效的，是时域信号上的塔架减振控制算法；但是此等方法严格依赖信号相序，在相序混乱情况下会短时失效甚至起到负阻尼效果，在柔塔风机中会造成严重后果。
期刊号为《机电工程》2014 31(3)：325-329，作者为陶学军，名称为“基于控制方法的风机塔架减振研究”的文章提出了在相序混乱情况下利用停机等手段避免风机进一步振动的方法，但风机大风时紧急停机本身就能引入较大的振动冲击，对柔性性塔架控制效果较差。
发明内容
本申请的目的在于提供一种风机塔架振动控制方法及装置，以解决现有控制方法对风机塔架振动抑制效果较差的问题。
为了实现上述目的，本发明提出一种风机塔架振动控制方法，包括以下步骤：
1)获取风机塔架振动信号；
2)根据风机塔架振动信号提取至少一个预设频率下的振动能量；
3)判断每个预设频率下的振动能量是否超过设定阈值，若超过，则进行该预设频率下的加阻控制。
有益效果是：本发明对预设频率下的振动能量进行了提取，在振动能量超过设定阈值时进行该预设频率下的加阻控制，针对性强，效果明确，负面效果小，控制精度较高，对风机塔架振动抑制效果较好。
进一步的，所述至少一个预设频率为至少一个模态频率和/或至少一个激励频率。
为了实现对风机塔架振动的实时控制，进一步的，步骤2)提取一个预设频率下的振动能量的过程包括以下步骤：
A：对风机塔架振动信号进行低通降噪滤波；
B：对步骤A中经低通降噪滤波后的信号进行陷波滤波，陷波滤波的传递函数
其中ω1为所述预设频率的值；ξ1为陷波滤波器阻尼系数；
C：将步骤A中进行低通降噪滤波后的信号减去步骤B中进行陷波滤波后的信号，差值记为振动能量。
进一步的，低通降噪滤波的传递函数为
其中k为滤波修正增益，τ为滤波实现常数。
进一步的，所述风机塔架振动信号为风机振动加速度信号。
另外，本发明还提出一种风机塔架振动控制装置，包括数据采集器、存储器、处理器以及存储在所述存储器中并可在处理器上运行的计算机程序，数据采集器的输入端用于采集风机塔架振动信号，数据采集器的输出端与处理器连接，用于将采集到的风机塔架振动信号传输给处理器；所述处理器在执行所述计算机程序时实现以下步骤：
1)获取风机塔架振动信号；
2)根据风机塔架振动信号提取至少一个预设频率下的振动能量；
3)判断每个预设频率下的振动能量是否超过设定阈值，若超过，则进行该预设频率下的加阻控制。
有益效果是：本发明对预设频率下的振动能量进行了提取，在振动能量超过设定阈值时进行该预设频率下的加阻控制，针对性强，效果明确，负面效果小，控制精度较高，对风机塔架振动抑制效果较好。
进一步的，所述至少一个预设频率为至少一个模态频率和/或至少一个激励频率。
为了实现对风机塔架振动的实时控制，进一步的，步骤2)提取一个预设频率下的振动能量的过程包括以下步骤：
A：对风机塔架振动信号进行低通降噪滤波；
B：对步骤A中经低通降噪滤波后的信号进行陷波滤波，陷波滤波的传递函数
其中ω1为所述预设频率的值；ξ1为陷波滤波器阻尼系数；
C：将步骤A中进行低通降噪滤波后的信号减去步骤B中进行陷波滤波后的信号，差值记为振动能量。
进一步的，低通降噪滤波的传递函数
其中k为滤波修正增益，τ为滤波实现常数。
进一步的，所述风机塔架振动信号为风机振动加速度信号。
附图说明
图1是本发明风机塔架振动控制方法实施例中的风机塔架振动控制方法流程图。
具体实施方式
以下结合实施例对本发明作进一步的详细描述。
风机塔架振动控制方法实施例：
本实施例涉及一种对风机塔架振动信号进行频域细分提取从而监测塔架振动性质，检测塔架振动状况及依据监测检测结果进行针对性的减振控制应对的方法。本实例选取某叶片直径为116m的机型进行举例。
第一步：确定风机塔架振动相关的特征频率，以及每个特征频率对应的监测阈值。风机塔架特征频率的确定可以借助bladed软件中的风机模型，直接由软件计算获得；然后在风场根据风机振动数据进行辨识修正，进行微调并确认。
该步骤一般在风机运行之前进行，得到的特征频率为风机的固有频率。通过软件仿真的方式获取特征频率是风电机组设计开发过程及机组控制环节中的常规手段。为论述简便，本实例例只举例一种特征频率，即风机塔架的一阶模态频率，大小为0.305Hz。实际应用中会根据重要程度对包括塔架其他模态频率和主要激励频率在的多种塔架重要特征频率进行辨识，比如1p、3p激励频率。
第二步：风机运行时，提取一阶模态特征频率下的风机塔架振动能量。如图1所示，塔架加速度传感器信号以V(t)表示，经正常低通降噪过滤，输出V2(t)。低通降噪传递函数以G1(s)表示，其传递函数可表示为：
其中：k为滤波修正增益，τ为滤波实现常数。
然后V2(t)经过一个以0.305Hz为极点的陷波滤波传递函数变换，输出不含塔架一阶模态能量的频率信号V3(t)。陷波滤波器的传递函数可表示为：
其中：ω1为0.305Hz；ξ1为陷波滤波器阻尼系数，本实施例中取0.4。
设V4(t)为时域信号中频率为塔架模态频率0.305Hz的分量；则其能量值为：
V4(t)＝V2(t)-V3(t) (3)
至此，实时提取一阶模态特征频率下的风机塔架振动能量完成，过程中没有进行耗时的数据傅里叶变换等频谱计算算法，而是利用陷波滤波这一实时控制方法进行分离，使该算法应用于实时控制成为可能。
第三步：特定频率振动的监测检测：利用提取出的相应特征频率下能量值V4(t)与设计好的此频率下监测阈值比较，当检测到相应特征能量超过监测阈值时启动对应频率的控制措施，否则继续监控。
第四步：有条件塔架控制措施。依据监测检测结果进行有条件塔架控制，塔架总体振动保护停机值相对原值可成倍放大，以达不到极限设计载荷为限；当检测到塔架非特定频率较大冲击振动时风机可继续运行，不做干预；当检测到相应特定频率的振动超过相应阈值时，进行特定频率加阻控制干预。加阻输入为相应特征频率这一单一频率。这样输入信号变换是相位信息明了，在对加速度进行积分变换时，相位情况易于处理，避免了加阻失效甚至负阻尼现象。
塔架加速度传感器安装于机舱内主机架上的专用传感器安装板上，既可以是安装方向与风机主轴垂直的、检测风机面内振动加速度的加速度传感器，也可以是安装方向与风机主轴平行的、检测风机面外振动加速度的加速度传感器。加阻传递函数计算输出结果加到相应执行机构中，如果塔架加速度传感器采集信号为风机面内振动加速度，则相应执行机构指的是变流器执行机构，通过变流器转矩进行受力干预及加阻；如果塔架加速度传感器采集信号为风机面外振动加速度，则相应执行机构指的是变桨控制执行机构，通过变桨进行受力干预及加阻。
本实施例中的步骤一至步骤四说明的是一个特征频率即塔架一阶模态频率的分离及监测、控制过程。塔架其他重要模态频率及激励频率同样需要类似处理过程，其关系式并行关系，在图1中以第i个特征频率分离监测表示。
当监测的特征频率的个数大于等于2时，分别判断各特征频率的特征能量是否超过相应的监测阈值；当某特征频率的特征向量超过相应的监测阈值时，相应进行该特征频率的加阻控制干预。
风机塔架振动控制装置实施例：
风机塔架振动控制装置包括数据采集器、存储器、处理器以及存储在所述存储器中并可在处理器上运行的计算机程序，数据采集器的输入端用于采集风机塔架振动信号，数据采集器的输出端与处理器连接，用于将采集到的风机塔架振动信号传输给处理器；所述处理器在执行所述计算机程序时实现风机塔架振动控制方法实施例中的风机塔架振动控制方法。
风机塔架振动控制方法的具体实施过程在上述风机塔架振动控制方法实施例中已经介绍，这里不做过多赘述。
Claims (10)
1.一种风机塔架振动控制方法，其特征在于，包括以下步骤：
1)获取风机塔架振动信号；
2)根据风机塔架振动信号提取至少一个预设频率下的振动能量；
3)判断每个预设频率下的振动能量是否超过设定阈值，若超过，则进行该预设频率下的加阻控制。
2.根据权利要求1所述的风机塔架振动控制方法，其特征在于，所述至少一个预设频率为至少一个模态频率和/或至少一个激励频率。
3.根据权利要求1所述的风机塔架振动控制方法，其特征在于，步骤2)提取一个预设频率下的振动能量的过程包括以下步骤：
A：对风机塔架振动信号进行低通降噪滤波；
B：对步骤A中经低通降噪滤波后的信号进行陷波滤波，陷波滤波的传递函数
其中ω1为所述预设频率的值；ξ1为陷波滤波器阻尼系数；
C：将步骤A中进行低通降噪滤波后的信号减去步骤B中进行陷波滤波后的信号，差值记为振动能量。
4.根据权利要求3所述的风机塔架振动控制方法，其特征在于，低通降噪滤波的传递函数为
其中k为滤波修正增益，τ为滤波实现常数。
5.根据权利要求1所述的风机塔架振动控制方法，其特征在于，所述风机塔架振动信号为风机振动加速度信号。
6.一种风机塔架振动控制装置，包括数据采集器、存储器、处理器以及存储在所述存储器中并可在处理器上运行的计算机程序，数据采集器的输入端用于采集风机塔架振动信号，数据采集器的输出端与处理器连接，用于将采集到的风机塔架振动信号传输给处理器；其特征在于，所述处理器在执行所述计算机程序时实现以下步骤：
1)获取风机塔架振动信号；
2)根据风机塔架振动信号提取至少一个预设频率下的振动能量；
3)判断每个预设频率下的振动能量是否超过设定阈值，若超过，则进行该预设频率下的加阻控制。
7.根据权利要求6所述的风机塔架振动控制装置，其特征在于，所述至少一个预设频率为至少一个模态频率和/或至少一个激励频率。
8.根据权利要求6所述的风机塔架振动控制装置，其特征在于，步骤2)提取一个预设频率下的振动能量的过程包括以下步骤：
A：对风机塔架振动信号进行低通降噪滤波；
B：对步骤A中经低通降噪滤波后的信号进行陷波滤波，陷波滤波的传递函数
其中ω1为所述预设频率的值；ξ1为陷波滤波器阻尼系数；
C：将步骤A中进行低通降噪滤波后的信号减去步骤B中进行陷波滤波后的信号，差值记为振动能量。
9.根据权利要求8所述的风机塔架振动控制装置，其特征在于，低通降噪滤波的传递函数为
其中k为滤波修正增益，τ为滤波实现常数。
10.根据权利要求6所述的风机塔架振动控制装置，其特征在于，所述风机塔架振动信号为风机振动加速度信号。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911007511.5A CN110761946B (zh) | 2019-10-22 | 2019-10-22 | 一种风机塔架振动控制方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911007511.5A CN110761946B (zh) | 2019-10-22 | 2019-10-22 | 一种风机塔架振动控制方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110761946A true CN110761946A (zh) | 2020-02-07 |
| CN110761946B CN110761946B (zh) | 2021-08-03 |
Family
ID=69332588
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911007511.5A Active CN110761946B (zh) | 2019-10-22 | 2019-10-22 | 一种风机塔架振动控制方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110761946B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113824426A (zh) * | 2020-06-18 | 2021-12-21 | 新疆金风科技股份有限公司 | 风力发电机组的滤波器控制方法及其装置 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105403402A (zh) * | 2015-11-12 | 2016-03-16 | 科诺伟业风能设备（北京）有限公司 | 风力发电机组传动链系统扭振状态的监测方法 |
| CN107850050A (zh) * | 2015-12-04 | 2018-03-27 | 远景能源（江苏）有限公司 | 风力涡轮机和减少摆振振动的运行风力涡轮机的方法 |
| WO2018077530A1 (en) * | 2016-10-28 | 2018-05-03 | Siemens Aktiengesellschaft | Damping wind turbine tower oscillations |
| US20190203692A1 (en) * | 2016-08-17 | 2019-07-04 | Vestas Wind Systems A/S | Dynamic controlled wind turbine shutdown |
-
2019
- 2019-10-22 CN CN201911007511.5A patent/CN110761946B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105403402A (zh) * | 2015-11-12 | 2016-03-16 | 科诺伟业风能设备（北京）有限公司 | 风力发电机组传动链系统扭振状态的监测方法 |
| CN107850050A (zh) * | 2015-12-04 | 2018-03-27 | 远景能源（江苏）有限公司 | 风力涡轮机和减少摆振振动的运行风力涡轮机的方法 |
| US20190203692A1 (en) * | 2016-08-17 | 2019-07-04 | Vestas Wind Systems A/S | Dynamic controlled wind turbine shutdown |
| WO2018077530A1 (en) * | 2016-10-28 | 2018-05-03 | Siemens Aktiengesellschaft | Damping wind turbine tower oscillations |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113824426A (zh) * | 2020-06-18 | 2021-12-21 | 新疆金风科技股份有限公司 | 风力发电机组的滤波器控制方法及其装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110761946B (zh) | 2021-08-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102980651B (zh) | 风电机组状态监测方法、装置及系统 | |
| AU2004259426B2 (en) | Method and apparatus for wind turbine rotor load control based on shaft radial displacement | |
| US8033788B2 (en) | Method for determining fatigue load of a wind turbine and for fatigue load control, and wind turbines therefor | |
| CN105673322B (zh) | 实现风力机mppt控制的变参数非线性反馈控制方法 | |
| CN109296506B (zh) | 一种风电机组的振动检测方法、控制方法及装置 | |
| CN104564536A (zh) | 用于变速风力涡轮机的减震系统 | |
| CN107356384A (zh) | 一种监测风力发电设备中叶片的状态的方法、计算设备和系统 | |
| CN108644069B (zh) | 叶片不平衡检测方法及装置 | |
| CN107762730B (zh) | 一种带有尾缘襟翼的大型变桨风力机控制系统及控制方法 | |
| CN110761946B (zh) | 一种风机塔架振动控制方法及装置 | |
| CN109958576B (zh) | 控制风力发电机组的转速的方法和装置 | |
| CN110259637A (zh) | 风力发电机组的叶片气动不平衡矫正方法、装置及设备 | |
| CN114258459A (zh) | 利用叶片载荷信号减少沿边振动 | |
| CN112368477A (zh) | 用于控制风力涡轮机以减小机舱振动的方法和系统 | |
| CN112018784A (zh) | 一种基于同步相量测量数据的次同步谐振溯源方法 | |
| US20230175484A1 (en) | Load reduction control method for wind turbine generator and device | |
| CN108386316A (zh) | 风力发电机组降噪系统和方法、控制装置及其存储介质 | |
| CN113565679A (zh) | 一种基于prony算法的风电机组运行控制方法、装置及存储介质 | |
| CN102878016A (zh) | 基于智能结构的大型风力机叶片摆振抑制系统及控制方法 | |
| EP4115079B1 (en) | Control method and device of a wind park | |
| CN112600199B (zh) | 一种基于转子动能的风电机组等效转动惯量在线评估方法 | |
| CN113794198A (zh) | 抑制宽频振荡的方法、装置、终端及存储介质 | |
| Lin | Novel modified Elman neural network control for PMSG system based on wind turbine emulator | |
| CN113137335A8 (zh) | 用于使用改进的功率曲线来操作风力涡轮的系统和方法 | |
| CN109799456B (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 |