CN107084830A - 一种大气边界层风洞的非定常风的模拟方法 - Google Patents
一种大气边界层风洞的非定常风的模拟方法 Download PDFInfo
- Publication number
- CN107084830A CN107084830A CN201710537102.0A CN201710537102A CN107084830A CN 107084830 A CN107084830 A CN 107084830A CN 201710537102 A CN201710537102 A CN 201710537102A CN 107084830 A CN107084830 A CN 107084830A
- Authority
- CN
- China
- Prior art keywords
- wind
- frequency
- frequency converter
- tunnel
- wind speed
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/02—Wind tunnels
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
本发明涉及一种大气边界层风洞的非定常风的模拟方法,其通过变频器的频率值变化调节多风机转速来控制风洞内风速,在不同频率下测量风机风速,标定风机的响应关系曲线,通过组态王软件控制变频器频率变化,实现恒定风速、正弦风速调节,通过热线风速仪的风速测量结果,对变频器的频率值变化进行修正,可保证风洞内风速的高频变化,克服风机惯性影响,在高流速的前提下实现了风洞内可控的非定常风,能够高度近似模拟野外阵风,为对大气边界层空气流动状态以及流动特性的研究提供参考意义。
Description
技术领域
本发明属于流体动力学研究技术领域,具体涉及一种大气边界层风洞的非定常风的模拟方法。
背景技术
土地沙漠化、沙尘暴等是全球性的严重环境问题,荒漠化防治是国际社会关注的焦点.按照《联合国防治沙漠化公约》确立的定义,我国荒漠化的土地面积高达280多万平方公里,占我国近1/3国土面积。在华北、东北、西北的“三北”地区,有80%的漠化、土壤风蚀及草场退化所形成的荒漠化面积达2460km2.尤其是我国西部的沙漠面积达70多万平方公里,约占国土面积的1/14,全球第二大流动性沙漠新疆塔克拉玛干沙漠一直是我国风沙危害最为严重的地区。在西部干旱半干旱地区特殊生态环境下,加之以水土开发为主的人类活动,使得以土地荒漠化为典型象征的生态环境恶化趋势越演越烈,以至影响到区域生态环境安全和区域社会经济的持续发展,违背了党和政府倡导经济发展的可持续性、注重经济发展和环境保护同步的科学发展观,这一基本国情决定了荒漠化问题始终是我国现代化建设进程中必须处理好的一个重大问题,也决定了防沙治沙工作始终是关系经济社会可持续发展全局的一项重要工作。风洞实验是风沙迁移研究的基本手段,其研究结果是防风治沙工作的理论基础,具有重要的实践意义。
大气边界层风洞(Atmospherical Boundary Layer Wind-tunnel,ABL Wind-tunnel)是用来模拟大气边界层内空气流动状态的一类风洞,其最早被用于研究污染物扩散问题,进而在环境、建筑、汽车、化工、机械等领域得到了广泛应用。ABL风洞在风沙迁移研究中的应用最早可追溯到上世纪40年代,Bagnold和Chepil等利用边界层风洞开展了风沙运动和土壤侵蚀的实验研究。
目前风洞的研究常限于对风洞内定常风速流场的模拟,但是野外阵风为无规律随机变化的非定常湍流,有大尺度结构的非定常特征和小尺度结构的高频脉动特性,因此,定常风的模拟并不能全面反应出空气流动状态以及流动特性。
发明内容
为了克服现有技术所存在的不足,本发明提供了一种能够克服风机惯性影响、能够在高流速的前提下实现风洞内风速的非定常变化,模拟过程简单且对大气边界层空气流动状态以及流动特性的研究提供参考意义的大气边界层风洞内非定常风的模拟方法。
本发明实现上述目的所采用的技术方案是:
该大气边界层风洞的非定常风的模拟方法由以下步骤实现:
(1)在风洞的进风端和出风端设置有多台与变频器连接的风机,进风端的送风风机与出风端的吸风风机串联,在风洞的试验段中部自上而下纵向梯度分布多个热线探头,用于检测风洞内的风速廓线分布;
(2)通过改变变频器的电流频率,在不同频率下测量多台风机同时运行时的风速,标定吸风风机和送风风机的变频器频率变化与风速跟随性响应曲线,拟合得到频率变化与风速跟随性响应关系为:u=F(fm),u表示风速,fm表示变频器的频率值;
(3)根据步骤(2)的标定结果,通过组态王软件控制调节变频器频率实现风机的恒定风速、正弦风速的调节,并对变频器的运转实现监控;
(4)对步骤(3)变频器的频率值变化进行修正,调整多台风机叠加时吸风风机与送风风机的正弦风速周期、幅值的差异以及相位角变化叠加,使风洞内生成的阵风能量谱分布符合柯尔莫哥洛夫的-5/3次方律,即可模拟得到非定常风。
进一步说明,所述步骤(3)中变频器的频率控制方程为:
fm=Asin(2π/Tt+θ)+k,0hz≤fm≤50hz;
式中fm为变频器的频率值;A为振幅;T为周期;t为时间;θ为相位角;k为截距。
进一步说明,所述步骤(4)具体是:用热线风速仪测量风洞中的主流风速,将所得风速序列进行傅里叶变换,得到风速能量谱分布,确定由送风风机和吸风风机的生成风主含能频率w3和w4,根据野外实测阵风数据,调整送风风机、吸风风机的变频器频率变化控制方程的周期、相位角及振幅使风洞内生成的阵风能量谱分布符合柯尔莫哥洛夫的-5/3次方律且与野外实测阵风的低频能谱一致。
进一步说明,所述送风风机的变频器频率变化控制方程为fm3=A-3sin(2π/T3t3+θ3)+k3,吸风风机的变频器频率变化控制方程为fm4=A-4sin(2π/T4t4+θ4)+k4;
式中:fm3为对应送风风机的变频器的频率值,fm4为对应吸风风机的变频器的频率值;A3、A4为对应送风/吸风风机的变频器的频率值周期变化的振幅;T3、T4为对应送风/吸风风机的变频器的频率值周期变化的周期;t3、t4为时间;θ3、θ4为对应送风/吸风风机的变频器的频率值周期变化的相位角;k3、k4为对应送风/吸风风机的变频器的频率值周期变化的平均值,即曲线的截距。
进一步说明,所述的风机为串联的2~5台。
本发明的大气边界层风洞的非定常风的模拟方法,通过变频器的频率值变化调节双/多风机转速来控制风洞内风速,在不同频率下测量风机风速,标定风机的响应关系曲线,通过组态王软件控制变频器频率变化,实现恒定风速、正弦风速调节,通过热线风速仪的风速测量结果,对变频器的频率值变化进行修正,可保证风洞内风速的高频变化,克服风机惯性影响,在高流速的前提下实现了风洞内可控的非定常风,能够高度近似模拟野外阵风,为对大气边界层空气流动状态以及流动特性的研究提供参考意义。
附图说明
图1为风洞的结构示意图。
图2为热线风速仪示意图。
图3为双风机可控恒定风频率变化与风速图像。
图4为双风机可控正弦风频率变化与风速图像。
图5为风洞内所模拟的非定常风的频率变化与风速曲线。
图6为风洞内所模拟的非定常风的频率变化与风速曲线。
具体实施方式
参见附图和实施例对本发明的技术方案进行进一步说明。
参见图1的风洞结构示意图,本实施例的风洞的直径为940mm,实验段长度为9000mm,在风洞的进风端安装送风风机,出风端安装出风风机,并在风洞的实验段中部自上而下纵向梯度分布多个用于检测风洞内风速的热线探头,如图2所示。将风机与变频器以及控制柜连接之后,开始模拟大气边界层风洞的非定常风,具体步骤如下:
(1)在风洞的两端各自安装风机,在风洞的试验段中部自上而下纵向梯度分布多个热线探头之后,确定风机以及与风机连接的变频器的工作参数,以调整其在工作参数范围内工作,如最大转速3000r/min、最大风速为30m/s以及变频器的最大工作频率50hz等。
(2)通过改变变频器的电流频率,测量吸风风机和送风风机的单机风速和双台风机同时开启时的双机风速,分别获得吸风风机和送风风机的变频器频率变化与单机风速跟随性响应曲线以及变频器频率变化与双机风速跟随性响应曲线,标定风速与变频器频率之间的定常变化关系为:u=F(fm)
u表示风速,fm表示变频器的频率值;
在单风机标定时,吸气风机和吹气风机对应的变频器频率分别以5hz的增量逐步增大,待风洞内流场稳定后依次测得该频率下对应的风洞主流风速,同样步骤变频器频率下降再进行一次测量,以测得数据拟合得到风速与变频器频率的标定关系,以确定双风机工作时的干扰损耗。之后两台风机同时开启按照单风机标定方法,拟合得到双台风机与风速与变频器频率的标定关系。
由于具体实验时粗糙元铺设、模型尺寸及通流情况会有变化,因此在每次实验前需要进行风速与变频器频率的标定。
(3)根据步骤(2)的标定结果,通过组态王程序编译对变频器的运转实现监控,同时通过组态画面周期调节变频器频率实现风机的恒定风速、正弦风速的调节;组态王程序发送命令通过PLC控制器去设定送风风机和吸风风机同时开启时对应的变频器频率的周期变化,对应的变频器频率值的控制方程为:
fm=Asin(2π/Tt+θ)+k,0hz≤fm≤50hz;
式中fm为变频器的频率值;A为振幅;T为周期;t为时间;θ为相位角;k为截距。
根据步骤(2)标定结果,若在组态王中设定双风机变频器控制方程的振幅A1为0,截距k恒定为40hz,可形成恒定风速如图3所示;若在组态王中设定双风机变频器控制方程的振幅A2皆为10hz,截距k2皆为35hz,周期T2皆为60s,相位角θ2皆为3π/2,则使双风机频率值周期性变化实现正弦风速的设定如图4所示。
根据步骤(2)标定结果,改变组态王中双风机变频器控制方程的振幅、周期、相位角、截距可模拟出需要的正弦风速。
(4)对步骤(3)变频器的频率值变化进行修正,调整多台风机叠加时吸风风机与送风风机的正弦风速周期、幅值的差异以及相位角变化叠加,使风洞内生成的阵风能量谱分布符合柯尔莫哥洛夫的-5/3次方律,即可模拟得到非定常风。
设定送风风机变频器频率值的控制方程为fm3=A3sin(2π/T3t3+θ3)+k3,吸风风机变频器频率值的控制方程为fm4=A4sin(2π/T4t4+θ4)+k4
式中fm3为对应送风风机的变频器的频率值,fm4为对应吸风风机的变频器的频率值;A3、A4为对应送风/吸风风机的变频器的频率值周期变化的振幅;T3、T4为对应送风/吸风风机的变频器的频率值周期变化的周期;t3、t4为时间;θ3、θ4为对应送风/吸风风机的变频器的频率值周期变化的相位角;k3、k4为对应送风/吸风风机的变频器的频率值周期变化的平均值,即曲线的截距。
用热线风速仪测量风洞中的主流风速,将所得风速序列进行傅里叶变换,得到风速能量谱分布,确定由送风风机和吸风风机的生成风主含能频率w3和w4,根据野外实测阵风数据,调整送风风机、吸风风机的变频器频率变化控制方程的周期、相位角及振幅使风洞内生成的阵风能量谱分布符合柯尔莫哥洛夫的-5/3次方律且与野外实测阵风的低频能谱一致。
对上述实施例所得到的非定常风进行测定,送风风机变频器频率值变化的振幅A5为20hz,截距k5为20hz,周期T5为120s,相位角θ5为3π/2;吸风风机变频器频率值变化的振幅A6为10hz,截距k6为10hz,周期T6为60s,相位角θ6为3π/2;双风机振幅、截距相差10hz,周期相差60s,非定常风的频率变化与风速曲线如图5所示。送风风机变频器频率值变化的振幅A7为10hz,截距k7为30hz,周期T6为60s,相位角θ7为0;吸风风机变频器频率值变化的振幅A8为10hz,截距k8为30hz,周期T8为60s,相位角θ8为π;双风机相位角相差π,非定常风的频率变化与风速曲线如图6所示。
Claims (5)
1.一种大气边界层风洞的非定常风的模拟方法,其特征在于由以下步骤组成:
(1)在风洞的进风端和出风端设置有多台与变频器连接的风机,进风端的送风风机与出风端的吸风风机串联,在风洞的试验段中部自上而下纵向梯度分布多个热线探头,用于检测风洞内的风速廓线分布;
(2)通过改变变频器的电流频率,在不同频率下测量多台风机同时运行时的风速,标定吸风风机和送风风机的变频器频率变化与风速跟随性响应曲线,拟合得到频率变化与风速跟随性响应关系为:u=F(fm),u表示风速,fm表示变频器的频率值;
(3)根据步骤(2)的标定结果,通过组态王软件控制调节变频器的频率值实现风机的恒定风速、正弦风速的调节,并对变频器的运转实现监控;
(4)对步骤(3)变频器的频率值变化进行修正,调整多台风机叠加时吸风风机与送风风机的正弦风速周期、幅值的差异以及相位角变化叠加,使风洞内生成的阵风能量谱分布符合柯尔莫哥洛夫的-5/3次方律,即可模拟得到非定常风。
2.根据权利要求1所述的大气边界层风洞的非定常风的模拟方法,其特征在于:所述步骤(3)中变频器的频率值控制方程为:
fm=Asin((2π/T)t+θ)+k,0hz≤fm≤50hz;
式中fm为变频器的频率值;A为振幅;T为周期;t为时间;θ为相位角;k为截距。
3.根据权利要求1所述的大气边界层风洞的非定常风的模拟方法,其特征在于:所述步骤(4)具体是:用热线风速仪测量风洞中的主流风速,将所得风速序列进行傅里叶变换,得到风速能量谱分布,确定由送风风机和吸风风机的生成风主含能频率w3和w4,根据野外实测阵风数据,调整送风风机、吸风风机的变频器频率变化控制方程的周期、相位角及振幅使风洞内生成的阵风能量谱分布符合柯尔莫哥洛夫的-5/3次方律且与野外实测阵风的低频能谱一致。
4.根据权利要求3所述的大气边界层风洞的非定常风的模拟方法,其特征在于:所述送风风机的变频器频率变化控制方程为fm3=A3sin((2π/T3)t3+θ3)+k3,吸风风机的变频器频率变化控制方程为fm4=A4sin((2π/T4)t4+θ4)+k4;
式中:fm3为对应送风风机的变频器的频率值,fm4为对应吸风风机的变频器的频率值;A3、A4为对应送风/吸风风机的变频器的频率值周期变化的振幅;T3、T4为对应送风/吸风风机的变频器的频率值周期变化的周期;t3、t4为时间;θ3、θ4为对应送风/吸风风机的变频器的频率值周期变化的相位角;k3、k4为对应送风/吸风风机的变频器的频率值周期变化的平均值,即曲线的截距。
5.根据权利要求1所述的大气边界层风洞的非定常风的模拟方法,其特征在于:所述的风机为串联的2~5台。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710537102.0A CN107084830B (zh) | 2017-07-04 | 2017-07-04 | 一种大气边界层风洞的非定常风的模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710537102.0A CN107084830B (zh) | 2017-07-04 | 2017-07-04 | 一种大气边界层风洞的非定常风的模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107084830A true CN107084830A (zh) | 2017-08-22 |
CN107084830B CN107084830B (zh) | 2018-12-18 |
Family
ID=59606566
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710537102.0A Active CN107084830B (zh) | 2017-07-04 | 2017-07-04 | 一种大气边界层风洞的非定常风的模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107084830B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107764578A (zh) * | 2017-11-22 | 2018-03-06 | 吉林省农业机械研究院 | 旋翼植保无人机喷雾性能试验台负压吸风装置及其控制方法 |
CN107966263A (zh) * | 2017-11-15 | 2018-04-27 | 南京航空航天大学 | 一种风洞实验阵风发生装置 |
CN111521365A (zh) * | 2020-05-18 | 2020-08-11 | 南京航空航天大学 | 基于旋转缝栅的可调频率和相位非定常流动控制实验装置 |
CN111929027A (zh) * | 2020-08-03 | 2020-11-13 | 中交天津港湾工程研究院有限公司 | 一种实验室风谱模拟方法 |
CN112067246A (zh) * | 2020-09-21 | 2020-12-11 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种非对称圆弧形边界层测压排架 |
CN113188752A (zh) * | 2021-06-11 | 2021-07-30 | 日照坤仑智能科技有限公司 | 一种非定常立式风洞 |
CN113221311A (zh) * | 2021-02-02 | 2021-08-06 | 中国人民解放军国防科技大学 | 一种大气边界层风速的不确定性量化方法 |
CN113326624A (zh) * | 2021-06-07 | 2021-08-31 | 国网河南省电力公司电力科学研究院 | 一种沙漠地区大气边界层高度的预测方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1411165A (en) * | 1973-10-01 | 1975-10-22 | Sage Action Inc | Gas flow visualisation system and control apparatus |
JPS61246649A (ja) * | 1985-04-24 | 1986-11-01 | Mitsubishi Heavy Ind Ltd | 流速制御装置 |
CN101520365A (zh) * | 2009-03-06 | 2009-09-02 | 山东科技大学 | 矿井巷道粉尘环境控制模拟试验系统 |
CN203249751U (zh) * | 2013-04-24 | 2013-10-23 | 东华大学 | 智能风洞风速调节及数据采集系统 |
CN106441789A (zh) * | 2016-11-08 | 2017-02-22 | 浙江大学 | 一种不同速度及温度分布的风洞流场产生装置 |
-
2017
- 2017-07-04 CN CN201710537102.0A patent/CN107084830B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1411165A (en) * | 1973-10-01 | 1975-10-22 | Sage Action Inc | Gas flow visualisation system and control apparatus |
JPS61246649A (ja) * | 1985-04-24 | 1986-11-01 | Mitsubishi Heavy Ind Ltd | 流速制御装置 |
CN101520365A (zh) * | 2009-03-06 | 2009-09-02 | 山东科技大学 | 矿井巷道粉尘环境控制模拟试验系统 |
CN203249751U (zh) * | 2013-04-24 | 2013-10-23 | 东华大学 | 智能风洞风速调节及数据采集系统 |
CN106441789A (zh) * | 2016-11-08 | 2017-02-22 | 浙江大学 | 一种不同速度及温度分布的风洞流场产生装置 |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107966263A (zh) * | 2017-11-15 | 2018-04-27 | 南京航空航天大学 | 一种风洞实验阵风发生装置 |
CN107966263B (zh) * | 2017-11-15 | 2023-10-24 | 南京航空航天大学 | 一种风洞实验阵风发生装置 |
CN107764578A (zh) * | 2017-11-22 | 2018-03-06 | 吉林省农业机械研究院 | 旋翼植保无人机喷雾性能试验台负压吸风装置及其控制方法 |
CN107764578B (zh) * | 2017-11-22 | 2024-03-29 | 吉林省农业机械研究院 | 旋翼植保无人机喷雾性能试验台负压吸风装置及其控制方法 |
CN111521365A (zh) * | 2020-05-18 | 2020-08-11 | 南京航空航天大学 | 基于旋转缝栅的可调频率和相位非定常流动控制实验装置 |
CN111929027A (zh) * | 2020-08-03 | 2020-11-13 | 中交天津港湾工程研究院有限公司 | 一种实验室风谱模拟方法 |
CN112067246A (zh) * | 2020-09-21 | 2020-12-11 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种非对称圆弧形边界层测压排架 |
CN112067246B (zh) * | 2020-09-21 | 2022-11-04 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种非对称圆弧形边界层测压排架 |
CN113221311A (zh) * | 2021-02-02 | 2021-08-06 | 中国人民解放军国防科技大学 | 一种大气边界层风速的不确定性量化方法 |
CN113326624A (zh) * | 2021-06-07 | 2021-08-31 | 国网河南省电力公司电力科学研究院 | 一种沙漠地区大气边界层高度的预测方法及系统 |
CN113188752A (zh) * | 2021-06-11 | 2021-07-30 | 日照坤仑智能科技有限公司 | 一种非定常立式风洞 |
CN113188752B (zh) * | 2021-06-11 | 2022-02-22 | 日照坤仑智能科技有限公司 | 一种非定常立式风洞 |
Also Published As
Publication number | Publication date |
---|---|
CN107084830B (zh) | 2018-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107084830B (zh) | 一种大气边界层风洞的非定常风的模拟方法 | |
Murata et al. | Experimental and numerical investigation of the effect of turbulent inflow on a Horizontal Axis Wind Turbine (Part I: Power performance) | |
Shigetomi et al. | Interactive flow field around two Savonius turbines | |
Chen et al. | Effects of pitch angle and blade camber on flow characteristics and performance of small-size Darrieus VAWT | |
CN103886513B (zh) | 一种风电场微观选址用地表粗糙度变化模型的建模方法 | |
Liu et al. | Towards uncovering the structure of power fluctuations of wind farms | |
Højstrup | Spectral coherence in wind turbine wakes | |
CN103969010A (zh) | 桥梁风浪流耦合场、弹性模型及动态响应试验测试系统 | |
CN104655391B (zh) | 一种应用于地形模型风场特性风洞试验的三维渐变式边界过渡装置 | |
Ma et al. | Hydrodynamics and wake flow analysis of a Π-type vertical axis twin-rotor tidal current turbine in surge motion | |
Danao et al. | Effects on the performance of vertical axis wind turbines with unsteady wind inflow: A numerical study | |
Chaudhry et al. | The influence of structural morphology on the efficiency of Building Integrated Wind Turbines (BIWT) | |
CN103792391B (zh) | 新型矿用风表校正装置 | |
Danao et al. | The Performance of a Vertical Axis Wind Turbine in Fluctuating Wind-A Numerical Study | |
Huang et al. | Improving the self-starting and operating characteristics of vertical axis wind turbine by changing center distance in part of blades | |
Cheng et al. | Simulation of wind gust structure in the atmospheric boundary layer with Lattice Boltzmann Method | |
Wosnik et al. | Experimental investigation and similarity solution of the axisymmetric turbulent wake with rotation | |
Jafari et al. | Performance analysis of a small wind turbine mounted inside a tall building | |
Grosvenor et al. | Performance and condition monitoring of tidal stream turbines | |
Qian et al. | Variations of horizontal and vertical velocities over two-dimensional transverse dunes: A wind tunnel simulation of the effect of windward slope | |
Ozbay et al. | An experimental investigation on the wake interference of multiple wind turbines in atmospheric boundary layer winds | |
Chaudhry et al. | Computational analysis to factor wind into the design of an architectural environment | |
Zhao et al. | Simulation of downburst in a multiple fan wind tunnel and research on its load on high-rise structure by wind tunnel experiment | |
CN108491617B (zh) | 一种基于尾流边界膨胀的bp模型的简化模型 | |
CN111982453A (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 |