CN112883667A - 一种基于大涡模拟在竖向脉动风荷载作用下的膜结构抗风优化方法 - Google Patents
一种基于大涡模拟在竖向脉动风荷载作用下的膜结构抗风优化方法 Download PDFInfo
- Publication number
- CN112883667A CN112883667A CN202110335340.XA CN202110335340A CN112883667A CN 112883667 A CN112883667 A CN 112883667A CN 202110335340 A CN202110335340 A CN 202110335340A CN 112883667 A CN112883667 A CN 112883667A
- Authority
- CN
- China
- Prior art keywords
- wind
- membrane structure
- large vortex
- pulsating
- model
- 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
Links
- 239000012528 membrane Substances 0.000 title claims abstract description 121
- 238000004088 simulation Methods 0.000 title claims abstract description 57
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000005457 optimization Methods 0.000 title claims abstract description 11
- 230000000694 effects Effects 0.000 title description 2
- 230000009471 action Effects 0.000 claims abstract description 25
- 238000013461 design Methods 0.000 claims abstract description 24
- 238000011160 research Methods 0.000 claims abstract description 16
- 238000013433 optimization analysis Methods 0.000 claims abstract description 4
- 230000004044 response Effects 0.000 claims description 45
- 238000005192 partition Methods 0.000 claims description 18
- 230000008878 coupling Effects 0.000 claims description 10
- 238000010168 coupling process Methods 0.000 claims description 9
- 238000005859 coupling reaction Methods 0.000 claims description 9
- 239000007787 solid Substances 0.000 abstract description 7
- 230000001808 coupling effect Effects 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 30
- 238000004364 calculation method Methods 0.000 description 20
- 239000010408 film Substances 0.000 description 7
- 238000006073 displacement reaction Methods 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 230000010349 pulsation Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000012821 model calculation Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 241000486761 Chilobrachys Species 0.000 description 1
- 241001077262 Conga Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004215 lattice model Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000012466 permeate Substances 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
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/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提供了一种基于大涡模拟在竖向脉动风荷载作用下的膜结构抗风优化方法，在UG中建立模型并导入ICEM中进行网格划分；设置结构模型参数；设置边界条件；将大涡模拟模型数据与RNGk‑ε模型对比，验证大涡模拟的正确性及精确性，利用大涡模拟方法研究不同风速、不同膜结构模型在竖向脉动风荷载作用下脉动风压分布特性规律，针对鞍形膜结构和拱形膜结构进行优化分析研究，进而进行膜结构抗风设计。本发明考虑流固耦合作用基于大涡模拟方法，对于膜结构在竖向脉动风荷载作用下的脉动风压特性进行研究，对于鞍形和拱形膜结构进行优化研究，对于完善膜结构的抗风设计理论有着重要的意义，并为膜结构建筑在实际工程中的抗风设计及优化设计提供了依据。
Description
技术领域
本发明涉及建筑行业结构抗风分析的技术领域，尤其涉及一种基于大涡模拟在竖向脉动风荷载作用下的膜结构抗风优化方法。
背景技术
目前大跨度膜结构的风压特性研究方法主要有现场实测，风洞试验及数值模拟三种方法。现场实测是最为直接有效的方法，依据对结构的实时监测从而得到真实的数据结果，但此方法因周期长、费时费钱、人力物力需求量大，且不能预前进行，现场实测有一定的局限性。风洞实验试验费用高，周期长，很难保证试验模型是抗风性能最好的模型，对于监测点的布置以及风压的同步测定等问题上有不足之处。数值模拟法与上述两种方法相比，具有周期短、费用低、控制方便的优点，并且可以进行反复模拟。该方法不但克服了诸多局限性，并且可以为建筑结构风场提供准确、高效的模拟结果。
目前膜结构风压荷载特性的数值模拟研究依然是数值计算中非常具有挑战性的工作，诸多问题尚处于初探阶段。国内外的许多学者做了很多行之有效的研究工作，为大跨度膜结构的风压特性研究奠定了基础。
1993年，Murakami等人通过大涡模拟方法模拟了德克萨斯科技大学的TTU标准低矮房屋建筑的风压和风流场，并同现场实测和风洞试验所得的数据进行了数据对比。潘毅等人采用简化气弹模型对考虑流固耦合效应的薄膜结构振动频率进行研究，通过对开敞式薄膜和封闭式薄膜在来流风下的振动频率特性进行研究，分析两种膜结构在不同条件下振动频率的变化规律和影响因素。王彬、杨庆山、朱伟亮等人通过运用弱耦合算法，实现了对一系列三维张拉膜结构的流固耦合数值模拟。孙晓颖基于弱耦合分区求解策略，搭建了薄膜结构三维流固耦合效应的数值模拟平台沈世钊、武岳等人采用分区法对张拉结构的风致流固耦合作用进行了研究，对湍流的模拟采用了LES模拟方法编制了二维结构的流固耦合程序，模拟实现了二维情况下对结构及其周围流场运动的实时模拟。康楠对采用弱耦合算法求解流固耦合问题时的求解策略进行了分析，基于大涡模拟方法对湍流控制方程进行了推导，研究了湍流控制方程的离散格式以及对离散方程进行数值求解的SIMPLE算法。李启基于弱耦合算法，将亚格子模型模块和动网格模块与求解任意拉格朗日欧拉动网格-大涡模拟(ALE-LES)方程的程序结合；采用九点网格重构法优化了动网格更新策略，解决了网格畸变的问题；开发了适用于模拟动边界钝体绕流风场的CFD计算程序，实现了二维封闭式平屋盖膜结构和四周封闭的三维大跨平屋盖膜结构的数值模拟。
目前对大跨度膜结构的抗风设计研究中研究的风荷载特性普遍针对于水平方向或是平均风荷载，对竖向脉动风荷载作用下膜结构的风压特性及膜结构优化却鲜有涉及，而竖向风荷载对于膜结构来讲，与水平风荷载有着同样重要的。
发明内容
针对上述技术问题，本发明的目的在于提供一种基于大涡模拟在竖向脉动风荷载作用下的膜结构抗风优化方法，研究膜结构在竖向风荷载作用下的脉动风压特性，为柔性膜结构的抗风设计提供参考依据。
为实现上述目的，本发明提供了一种基于大涡模拟在竖向脉动风荷载作用下的膜结构抗风优化方法，包括如下步骤：
步骤1：根据模型尺寸在UG软件中建立基准模型并导入ICEM中进行网格划分；
步骤2：设置边界条件参数；
步骤3：将大涡模拟模型数据与RNG k-ε模型数据对比，证明大涡模拟的正确性及精确性；
步骤4：利用大涡模拟方法研究不同风速下不同膜结构模型在竖向脉动风荷载作用下脉动风压特性规律；
步骤5：在步骤4的基础上针对鞍形膜结构和拱形膜结构进行优化分析研究，进而进行膜结构抗风设计。
可选的，所述步骤1中建立模型选取波浪形、鞍形及连续拱形三种不同形式的大跨度膜结构计算模型；网格划分采用混合网格划分，膜结构外部区域采用结构化网格划分。在划分膜结构网格时利用ANSYS ICEM局部加密、边界层加密等技术对流动变化较剧烈和关键位置进行网格加密。计算域离散均匀，重要部件的y+值控制在1左右，以达到湍流模型的最佳要求。
可选的，所述步骤2的边界条件参数的设置中，采用基于压力的求解方式，求解使用压力-速度场耦合，并且差分方式均选用高阶精度进行计算。
可选的，所述步骤3中，在于选取的RNG k-ε是在标准的k-ε模型上改进得到的，改进后的湍流模型为Prandtl数提供解析公式，相较于之前的常数有效的增加了模拟仿真的精确性，并且模型考虑了湍流涡结构影响。RANS方程计算得到的一般是时均流场，所以流场中是没有流动细节的，大多数研究者均采用此湍流模型，方法可行行真实可靠。
进一步的，所述步骤4中，本发明采用大涡模拟方法对竖向脉动风荷载进行时程模拟，研究三种不同形状的膜结构在竖向脉动风作用下的抖振响应及脉动风压分布规律，对结构进行脉动风速分别为10m/s、15m/s及20m/s下的膜结构脉动风压的数值模拟并得出分区系数值。
所述步骤5中，从多角度、系统化将流线设计理念巧妙渗透到膜结构外观设计各环节，包括流线设计原则、流线设计方法手段等，对鞍形膜结构和拱形膜结构进行优化设计。
由上，本发明基于计算流体力学中的大涡模拟方法，选取实际工程中较为典型的波浪形、鞍形及连续拱形大跨度膜结构形式，研究了竖向脉动风荷载作用下，各不同形式膜结构的脉动风压特性。研究了不同风速工况下，各种膜结构在竖向脉动风作用下结构表面的脉动风压特性，分析得到竖向脉动风作用下不同形状膜结构的脉动风压分布规律。研究了基于流线型设计理念，对鞍形膜结构和拱形膜结构进行优化。研究成果对膜结构在竖向风荷载下的抗风设计及优化设计提供参考。
附图说明
图1为本发明具体实施方式中基于大涡模拟方法在竖向脉动风荷载作用的膜结构抗风设计及优化方法的流程图；
图2为本发明数值模拟中涉及到的模型；
其中(a)为鞍形膜结构模型；(b)为波浪形膜结构模型；(c)为连续拱形膜结构模型；
图3为三种模型的网格划分图；
其中(a)为鞍形膜结构网格划分图；(b)为连续拱形膜结构网格划分图；(c)为波浪形膜结构网格划分图；
图4为本发明数值模拟中各模型相应的计算模型尺寸图；
图5为不同湍流模型计算结果对比图；
其中(a)为LES计算风振响应图；(b)为RNG k-ε计算风振响应图；
图6为不同湍流模型位移计算结果对比图；
其中(a)为LES计算位移图；(b)为RNG k-ε计算位移图
图7为不同形状膜结构屋盖脉动风速时程对比图；
其中(a)鞍形膜结构屋盖竖向脉动风速时程图；(b)为波浪形膜结构屋盖竖向脉动风速时程图；(c)为连续拱形膜结构屋盖竖向脉动风速时程图；
图8为鞍形膜结构屋盖测点分布示意图；
图9为鞍形膜结构屋盖在不同风速下的风振响应对比图；
其中(a)为风速10m/s的脉动风时的风振响应图；(b)为风速15m/s的脉动风时的风振响应图；(c)为风速20m/s的脉动风时的风振响应图；
图10为鞍形膜结构在不同风速下的脉动风压分区对比图；
其中(a)为鞍形膜结构在10m/s的脉动风速下的脉动风压分区图；(b)为鞍形膜结构在15m/s的脉动风速下的脉动风压分区图；(c)为鞍形膜结构在20m/s的脉动风速下的脉动风压分区图；
图11为波浪形膜结构屋盖测点分布示意图；
图12为波浪形膜结构屋盖在不同风速下的风振响应对比图；
其中(a)为风速10m/s的脉动风时的风振响应图；(b)为风速15m/s的脉动风时的风振响应图；(c)为风速20m/s的脉动风时的风振响应图；
图13为波浪形膜结构在不同风速下的脉动风压分区对比图；
其中(a)为波浪形膜结构在10m/s的脉动风速下的脉动风压分区图；(b)为波浪形膜结构在15m/s的脉动风速下的脉动风压分区图；(c)为波浪形膜结构在20m/s的脉动风速下的脉动风压分区图；
图14为连续拱形膜结构屋盖测点分布示意图；
图15为连续拱形膜结构屋盖在不同风速下的风振响应对比图；
其中(a)为风速10m/s的脉动风时的风振响应图；(b)为风速15m/s的脉动风时的风振响应图；(c)为风速20m/s的脉动风时的风振响应图；
图16为连续拱形膜结构在不同风速下的脉动风压分区对比图；
其中(a)为连续拱形膜结构在10m/s的脉动风速下的脉动风压分区图；(b)为连续拱形膜结构在15m/s的脉动风速下的脉动风压分区图；(c)为连续拱形膜结构在20m/s的脉动风速下的脉动风压分区图。
具体实施方式
下面参见附图对本发明所述基于大涡模拟在竖向脉动风荷载作用下的膜结构抗风优化方法进行详细说明。
如图1至16所示，基于大涡模拟在竖向脉动风荷载作用下的膜结构抗风优化方法，本发明的大涡模拟方法的竖向脉动风荷载作用下的膜结构抗风设计及优化方法，包括以下步骤：
步骤1：根据模型尺寸在UG软件中建立基准模型并导入ICEM中进行网格划分。
本实施方式中，建立模型选取波浪形、鞍形及连续拱形三种不同形式的大跨度膜结构作为计算模型；流体网格的划分方式为混合网格划分方式。如图3所示，其中马鞍形结构网格数351万，655060个节点；连续拱形结构网格数335万，3441565个节点；波浪形结构网格数338万，3470026个节点。
步骤2：设置边界条件参数；
本实施方式中，模型尺寸缩放比为1：100，结构模型如图2，不同结构模型的具体参数见图4所示，其中假设计算模型宽度为L1，长度为L2，最大高度为Hmax，选取计算域精度的准确性直接影响其数值模拟结果评估程度，为减小误差，本发明选取的计算域尺寸为20、L1×10、L2×10Hmax，分别对应x、y和z轴。因满足阻塞率小于3％，本发明将计算模型放置距入风口垂直距离W1大于3Hmax，计算模型的侧面和顶面距计算域边界的距离要大于4Hmax；计算模型的背风面距出风口的距离W2通常取9-10Hmax。基本方向角为沿z轴负方向，将模型凹面作为迎风面。本发明选取对网格边界层要求很高的LES湍流模型进行计算，采用基于压力的求解方式，求解使用压力-速度场耦合，并且差分方式均选用高阶精度进行。整个模拟过程使用残差收敛，以数值为10的均方残差做为迭代收敛标准，当求解器收敛到设置的精度后自动停止计算，并将自动进入后处理界面。
步骤3：将大涡模拟模型数据与RNG k-ε模型对比，证明大涡模拟的正确性及精确性。
本实施方式中，将大涡模拟模型与RNG k-ε模型模拟对比，图5为不同湍流模型风振响应结果对比图，根据其中a、b对比可以明显看出随着计算时间的推进，风振相应呈现逐渐增大趋势，但采用RNG k-ε计算时其最大值在300左右进行波动，而LES计算的结果则是在350左右进行波动，并且波动的幅值明显大于RNG模型计算结果。由此可见，采用精度更高的LES模型可以更加精确的捕捉到膜表面的细节流动，同时LES模型对流场的变化也更为敏感。
图6为不同湍流模型位移计算对比，根据其中a、b膜结构位移对比可以发现，对于不同的湍流模型，其位移变化趋势没有明显区别。RNG湍流模型也可以在一定程度上计算出较为靠谱的结果，并且计算速度快，如果只观察流动分布趋势，可以采用此湍流模型计算。在鞍形膜结构的迎风面上，由于气流发生偏转流速升高，此处压力系数下降至-1893，膜结构表面压力系数最大处位于中部达到757，对应于图6中的最小变形区域。由于RNG湍流模型计算误差较大，在提取详细数据时明显存在无法分辨的情况，而精度更高的LES模型则可以很精细的分辨出不同风压系数分区，有利于研究膜结构表面更为细致的流动现象，验证了大涡模拟的正确性及精确性。
步骤4：利用大涡模拟方法研究不同风速下不同膜结构模型在竖向脉动风荷载作用下脉动风压特性规律。
本实施方式中，脉动风速共有三种，分别为10m/s、15m/s、20m/s。大跨度膜结构形式共有三种为鞍形膜结构，连续拱形膜结构和波浪形结构。
图9为鞍形膜结构屋盖在不同风速下的风振响应对比图，其中a为10m/s风速的风振响应图，b为15m/s风速的风振响应图，c为20m/s风速的风振响应图；根据其中a、b、c三个对比图可知，风振响应频率并不均匀，并且尾部频率开始变得剧烈，从模拟结果中可以得出，第一耦合作用与竖向脉动风风振响应有直接关系，第二这些现象符合脉动紊流风的随机性和无序性，流层不再是成层流动，而是呈现不规则紊动，并且质点间也有相互混掺问题出现。以上两条结论可以大致总结出，竖向脉动风引起的风振响应对膜结构表面影响较为严重，所以必须考虑该情况对结构进行抗风设计时的影响。
图10为鞍形膜结构屋盖在不同风速下的脉动平均风压分区对比图，其中a为10m/s风速的脉动风压分区图，b为15m/s风速的脉动风压分区图，c为20m/s风速的脉动风压分区图；根据其中a、b、c三个对比图可知，采用欧洲规范中的风压分区排布方式，对三种风速下的鞍形膜结构屋面进行脉动风进行时成模拟，并提取数据，生成风压系数平均分布图，由脉动风谱生成的脉动风速图上可以看出，速度时程毛刺较多，曲线振动非常强烈，证明受竖向风作用的影响极其明显。整体竖向风脉动风考虑风振响应时，脉动风压值比只考虑平均风风压时，风压数值增大明显，可见竖向脉动风在对大跨度膜结构做抗风研究时是必须考虑的影响因素。
图12为波浪形膜结构屋盖在不同风速下的风振响应对比图，其中a为10m/s风速的风振响应图，b为15m/s风速的风振响应图，c为20m/s风速的风振响应图；根据其中a、b、c三个对比图可知，当风速提升到20m/s时，响振频率变得极为剧烈，这些现象的出现跟波浪形屋盖形状有关，因为屋盖表面高度不断的变化，结构自由频率降低，波浪形屋盖对风振响应变得更加敏感，说明竖向脉动风引起的风振响应对波浪形膜结构屋盖表面影响较为严重，所以必须考虑风振响应。
图13为波浪形膜结构屋盖在不同风速下的脉动平均风压分区对比图，其中a为10m/s风速的脉动平均风压分区图，b为15m/s风速的脉动平均风压分区图，c为20m/s风速的脉动平均风压分区图；根据其中a、b、c三个对比图可知，平均风压大致呈对称分布，最大风压值主要集中在膜结构屋盖中央区域，但由于是脉动风，不是完全对称，与大涡模拟模型有关，模拟数值结果更贴合实际工程。
图15为连续拱形膜结构屋盖在不同风速下的风振响应对比图，其中a为10m/s风速的风振响应图，b为15m/s风速的风振响应图，c为20m/s风速的风振响应图；根据其中a、b、c三个对比图可知，其响应频率大多较为平稳，末尾部位出现衰减，这符合随机性振动现象。但当风速达到20m/s时，初始的峰值频率和尾部的峰值频率接近，并且尾部频率开始变得剧烈，曲线振动强烈，这表明随着模态的光滑运动，风振响应逐渐接近并且稳定于起始点，这些结论充分体现了竖向脉动风引起的风振响应对膜结构表面影响是巨大的，应当予以考虑。
图16为连续拱形膜结构在不同风速下的脉动风压分区对比图，其中a为10m/s风速的脉动平均风压分区图，b为15m/s风速的脉动平均风压分区图，c为20m/s风速的脉动平均风压分区图；根据其中a、b、c三个对比图可知，平均风压最大值部分均集中在一侧，这可能是大跨度膜结构形状导致，风压起伏较大，曲线振动非常强烈，证明受竖向风作用的影响极其明显。因此竖向脉动风在对大跨度膜结构做抗风研究时，应当重视。
步骤5：在步骤5的基础上针对鞍形膜结构和拱形膜结构进行优化分析研究，进而进行膜结构抗风设计。
本实施方式中，针对于鞍形膜结构和拱形膜结构基于流线型设计理念进行优化设计。本发明从多角度、系统化剖析空气动力学因素对膜结构外部形状的影响，深化把握环境风速、空气阻力大小二者之间的关系，包括膜结构所受空气阻力类型，即迎风阻力、形状阻力，迎风阻力大小和膜结构横截面积大小相关，形状阻力大小和膜结构外形相关，结合流线设计要求以及重点，将流线设计理念巧妙渗透到膜结构外观设计各环节，包括流线设计原则、流线设计方法手段等，最大化发挥流线设计作用，有效解决空气阻力问题的同时优化膜结构使用性能。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (4)
1.一种基于大涡模拟在竖向脉动风荷载作用下的膜结构抗风优化方法，其特征在于，包括如下步骤：
步骤1：根据模型尺寸在UG软件中建立基准模型并导入ICEM中进行网格划分；
步骤2：设置边界条件参数；
步骤3：将大涡模拟模型数据与RNG k-ε模型数据对比，证明大涡模拟的正确性及精确性；
步骤4：利用大涡模拟方法研究不同风速下不同膜结构模型在竖向脉动风荷载作用下脉动风压特性规律；
步骤5：在步骤4的基础上针对鞍形膜结构和拱形膜结构进行优化分析研究，进而进行膜结构抗风设计。
2.根据权利要求1所述的基于大涡模拟在竖向脉动风荷载作用下的膜结构抗风优化方法，其特征在于，所述步骤1中建立的模型选取波浪形、鞍形及连续拱形三种不同形式的大跨度膜结构；网格划分采用混合网格划分，膜结构外部区域采用结构化网格划分。
3.根据权利要求1所述的基于大涡模拟在竖向脉动风荷载作用下的膜结构抗风优化方法，其特征在于，所述步骤2的边界条件参数的设置中，采用基于压力的求解方式，求解使用压力-速度场耦合，并且差分方式均选用高阶精度进行计算。
4.根据权利要求1所述的基于大涡模拟在竖向脉动风荷载作用下的膜结构抗风优化方法，其特征在于，所述步骤4中，采用大涡模拟方法对竖向脉动风荷载进行时程模拟，研究三种不同形状的膜结构在竖向脉动风作用下的抖振响应及脉动风压分布规律，对结构进行脉动风速分别为10m/s、15m/s及20m/s下的膜结构脉动风压的数值模拟并得出分区系数值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110335340.XA CN112883667A (zh) | 2021-03-29 | 2021-03-29 | 一种基于大涡模拟在竖向脉动风荷载作用下的膜结构抗风优化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110335340.XA CN112883667A (zh) | 2021-03-29 | 2021-03-29 | 一种基于大涡模拟在竖向脉动风荷载作用下的膜结构抗风优化方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112883667A true CN112883667A (zh) | 2021-06-01 |
Family
ID=76039900
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110335340.XA Pending CN112883667A (zh) | 2021-03-29 | 2021-03-29 | 一种基于大涡模拟在竖向脉动风荷载作用下的膜结构抗风优化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112883667A (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113607373A (zh) * | 2021-07-05 | 2021-11-05 | 北京航空航天大学 | 基于离散多点测量数据的大跨度屋面风压分布快速重构算法 |
| CN114417750A (zh) * | 2022-01-20 | 2022-04-29 | 重庆大学 | 基于主动-被动混合试验技术的三维气动导纳识别方法及系统 |
| CN114492251A (zh) * | 2022-04-18 | 2022-05-13 | 国家超级计算天津中心 | 超算环境的低速流场发散处理方法、装置、设备及介质 |
| CN114608788A (zh) * | 2022-03-02 | 2022-06-10 | 广州大学 | 一种钢骨架支承式膜结构风致动力响应测量方法及系统 |
| CN116150893A (zh) * | 2023-04-19 | 2023-05-23 | 北京航空航天大学 | 一种航空发动机篦齿封严机构设计方法及篦齿封严机构 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111079330A (zh) * | 2019-12-09 | 2020-04-28 | 桂林理工大学 | 一种考虑流固耦合效应的大跨度膜结构抗风设计方法 |
-
2021
- 2021-03-29 CN CN202110335340.XA patent/CN112883667A/zh active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111079330A (zh) * | 2019-12-09 | 2020-04-28 | 桂林理工大学 | 一种考虑流固耦合效应的大跨度膜结构抗风设计方法 |
Non-Patent Citations (2)
| Title |
|---|
| 张强: "基于动边界技术的单向张拉膜结构流固耦合效应研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》, pages 2 * |
| 赵欣: "考虑流固耦合作用的膜结构在竖向风荷载作用下的风压特性研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》, pages 3 - 4 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113607373A (zh) * | 2021-07-05 | 2021-11-05 | 北京航空航天大学 | 基于离散多点测量数据的大跨度屋面风压分布快速重构算法 |
| CN114417750A (zh) * | 2022-01-20 | 2022-04-29 | 重庆大学 | 基于主动-被动混合试验技术的三维气动导纳识别方法及系统 |
| CN114417750B (zh) * | 2022-01-20 | 2024-05-24 | 重庆大学 | 基于主动-被动混合试验技术的三维气动导纳识别方法及系统 |
| CN114608788A (zh) * | 2022-03-02 | 2022-06-10 | 广州大学 | 一种钢骨架支承式膜结构风致动力响应测量方法及系统 |
| CN114492251A (zh) * | 2022-04-18 | 2022-05-13 | 国家超级计算天津中心 | 超算环境的低速流场发散处理方法、装置、设备及介质 |
| CN114492251B (zh) * | 2022-04-18 | 2022-07-15 | 国家超级计算天津中心 | 超算环境的低速流场发散处理方法、装置、设备及介质 |
| CN116150893A (zh) * | 2023-04-19 | 2023-05-23 | 北京航空航天大学 | 一种航空发动机篦齿封严机构设计方法及篦齿封严机构 |
| CN116150893B (zh) * | 2023-04-19 | 2023-06-16 | 北京航空航天大学 | 一种航空发动机篦齿封严机构设计方法及篦齿封严机构 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112883667A (zh) | 一种基于大涡模拟在竖向脉动风荷载作用下的膜结构抗风优化方法 | |
| Toja-Silva et al. | A review of computational fluid dynamics (CFD) simulations of the wind flow around buildings for urban wind energy exploitation | |
| CN104298869B (zh) | 一种弹性水翼的流固耦合特性数值预测方法 | |
| Hong et al. | Validation of an open source CFD code to simulate natural ventilation for agricultural buildings | |
| King et al. | Optimization of wind plant layouts using an adjoint approach | |
| CN103969010B (zh) | 桥梁风浪流耦合场、弹性模型及动态响应试验测试系统 | |
| Ravensbergen et al. | A variational multiscale framework for atmospheric turbulent flows over complex environmental terrains | |
| CN104573269B (zh) | 一种基于强耦合整体技术的索膜结构抗风设计方法 | |
| CN109726433B (zh) | 基于曲面边界条件的三维无粘低速绕流的数值模拟方法 | |
| Jafari et al. | Aerodynamic shape optimization of rectangular and elliptical double-skin façades to mitigate wind-induced effects on tall buildings | |
| CN109885864A (zh) | 一种三维钢桥塔涡激振动计算方法 | |
| CN104200102A (zh) | 一种基于cfd的光伏阵列风荷载预测方法 | |
| Ran et al. | Three-dimensional numerical simulation of flow in trapezoidal cutthroat flumes based on FLOW-3D | |
| CN105260806A (zh) | 一种管路系统的流固耦合动力学特性预测方法 | |
| Li et al. | Investigation of the flow structure of single-and dual-celled tornadoes and their wind effects on a dome structure | |
| CN104200003A (zh) | 基于流固耦合的自锁阻尼叶片气流弯应力设计方法 | |
| CN113360979B (zh) | 钢箱梁抑涡除湿一体流动智能控制系统及使用方法 | |
| CN110543677A (zh) | 一种涡特征驱动的旋转湍流pans模型 | |
| CN102163263B (zh) | 风机叶片振动位移及其威布尔分布拟合方法 | |
| CN109141542B (zh) | 基于临界水深槽工作原理的矩形渠道无槽量水方法 | |
| CN111079330A (zh) | 一种考虑流固耦合效应的大跨度膜结构抗风设计方法 | |
| CN115906682A (zh) | 一种基于水动力模型的鱼礁设计优化方法 | |
| Amiri et al. | A review of physical and numerical modeling techniques for horizontal-axis wind turbine wakes | |
| CN108052692A (zh) | 垃圾在气力输运管线中运移规律的计算方法 | |
| Sacedo | Computational Study of the Unsteady Pressure Around a 2D Circular Cylinder Undergoing Forced Motion |
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 |