CN109543302A - 一种建筑屋顶檐口的风效应数值模拟方法 - Google Patents

一种建筑屋顶檐口的风效应数值模拟方法 Download PDF

Info

Publication number
CN109543302A
CN109543302A CN201811401246.4A CN201811401246A CN109543302A CN 109543302 A CN109543302 A CN 109543302A CN 201811401246 A CN201811401246 A CN 201811401246A CN 109543302 A CN109543302 A CN 109543302A
Authority
CN
China
Prior art keywords
wind
building
building roof
numerical simulation
cornice
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
CN201811401246.4A
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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201811401246.4A priority Critical patent/CN109543302A/zh
Publication of CN109543302A publication Critical patent/CN109543302A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Civil Engineering (AREA)
  • Architecture (AREA)
  • Software Systems (AREA)
  • Computer Graphics (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明提供了一种建筑屋顶檐口的风效应数值模拟方法，包括：步骤一，建立计算模型，在数值模拟时给定有限的三维计算区域，确定计算流域；步骤二，进行网格划分并设定边界条件；步骤三，进行求解，求解计算采用3D单精度求解器为压力基隐式；步骤四，对建筑物屋顶檐口垂线位置对应的建筑物部分门窗的开敞或突然开孔导致的内压增大对具有檐口的屋顶结构进行风效应数值模拟，将模拟位置转移到建筑物屋顶檐口垂线位置对应的建筑物部分门窗处，建立内压传递方程，采用Schinozuka‑Deodatis法模拟得到脉动风速的功率谱，对开孔结构的内压传递方程进行推导。风洞试验和模拟的结果比较吻合，数值模拟可以模拟真实的风环境，为建筑屋顶檐口的抗风设计提供参考依据。

Description

一种建筑屋顶檐口的风效应数值模拟方法
技术领域
本发明涉及风效应数值模拟领域，特别是一种针对建筑屋顶檐口的风效应数值模拟方法。
背景技术
随着社会经济的发展和科技的进步，各种造型独特、结构新颖的空间结构大量涌现，具有檐口的屋顶结构因为具有优美的造型和良好的性能，被广泛应用于各种大型公共建筑中，由于这些建筑物大多具有质量轻、柔性大、阻力小的特点，对风荷载敏感，风荷载往往成为其结构设计的主要控制荷载，此外，由于建筑物部分门窗的开敞或突然开孔导致的内压增大对具有檐口的屋顶结构安全产生较大的影响，大量的风灾调查表明，屋顶结构的风致破坏在很大程度上都是由于内外压联合作用引起的，而现有技术对内压变化的产生机理以及评估尚处于研究阶段，在进行结构的抗风设计时，国内外现行载荷规范仅是提出名义封闭或者开敞时的内压系数建议值。
结构风工程学是综合空气动力学、气象学和建筑结构力学、振动工程学、试验力学、气动弹性力学等多学科的多领域交叉学科，主要研究风对结构响应、减小结构风致响应和风损风毁事故措施。主要研究方法包括现场实测、理论分析、风洞试验等，以获得结构表面风压分布情况，由此直接得到结构响应或者计算结构风致响应和等效风荷载。现场实测是风工程研究中最直接有效的方法，风洞试验是当今风工程研究中最广泛采用的一种方法。随着计算机技术的迅速发展，新的有效方法数值模拟逐渐发展起来，可同风洞试验一样来预测建筑物表面风压和建筑物周围风场，从而同实测以及风洞试验协同形成研究流体流动问题的完整体系。数值模拟方法可以看作利用计算机模拟生成一个假想的风洞，即数值风洞，数值求解问题的基本思想在于：将原来在时间与空间内有着连续坐标的物理量的场，替代成为一系列有限个离散点(称为节点)上的值的集合，然后通过这些离散点上变量之间的相互关系建立起代数方程(即离散方程)，并对这些代数方程求解来得到变量的近似值。然而，目前的数值模拟中需要改进湍流模型使其更加成熟的解决结构风工程中吨体绕流的问题，并且风载脉动效应的数值模拟精度需要更加提高从而更加真实的模拟湍流来流条件，此外，数值模拟的计算量通常很大，需要进一步优化数值模拟的方法降低对计算机CPU的要求。
发明内容
为克服现有技术的缺陷，本发明的目的在于提供一种建筑屋顶檐口的风效应数值模拟方法，避免风洞试验的部分局限性，更高效的得到更准确的建筑屋顶檐口风场模拟效果，并且对于檐口表面以及建筑物内压都能够进行模拟，对周围风速和湍流特性能够进行预测。
本发明的目的在于提供一种建筑屋顶檐口的风效应数值模拟方法，包括：
步骤一，建立计算模型，在数值模拟时给定有限的三维计算区域，确定计算流域；
步骤二，进行网格划分并设定边界条件；
步骤三，进行求解，求解计算采用3D单精度求解器为压力基隐式；
步骤四，对建筑物屋顶檐口垂线位置对应的建筑物部分门窗的开敞或突然开孔导致的内压增大对具有檐口的屋顶结构进行风效应数值模拟，将模拟位置转移到建筑物屋顶檐口垂线位置对应的建筑物部分门窗处，建立内压传递方程，采用Schinozuka-Deodatis法模拟得到脉动风速的功率谱，对开孔结构的内压传递方程进行推导。
优选的，所述步骤一中确定所述计算流域对应的数值风洞对阻塞度的要求，在建筑屋顶檐口情况下表示为：(屋顶檐口最大的迎风面积/流域横截面面积)*100％<3％，保证计算域足够大的同时，保证边界与模型之间具有足够的距离，同时考虑到随着计算域的增大网格数量增加、计算量加大。
优选的，所述步骤一中确定所述计算流域设计还包括使得流场得到充分的发展的计算流域范围。
优选的，所述步骤一重所述计算区域尺寸设计为700m*400m*200m，将建筑屋顶檐口移到计算流域沿风流向的前1/3内，出口边界和建筑物后有着远大于建筑物高度的10倍的距离。
优选的，所述步骤二的网格划分包括：建筑屋细部形状中对计算结果影响较大的位置是拐角区域及建筑屋壁面区域加密网格，根据湍流模型的不同，当计算结果随着网格的分辨率变化时，取最小网格宽度为建筑物宽度及进深的1/50。
优选的，所述步骤二的网格划分包括：对整个计算流域分块进行网格划分，对于檐口周围的比较近的风敏感区湍流复杂，且体型复杂，采用适应性较强的非结构化四面体网格进行局部加密，远离檐口的区域，采用四面体单元离散，远离建筑屋顶檐口逐渐稀疏。
优选的，所述步骤二的边界条件的设定，模拟计算中假定气流为定常的不可压缩流动，进流面采用速度入口，其速度随高度的变化为UZ＝U0(Z/Z0)a，其中，Z0和U0分别为参考高度和该高度处的平均风速，本实施例中参考高度为20m，平均风速为8m/s，Z和UZ分别为计算高度和该高度处的平均风速，a为地面粗糙度指数，本实施例处于C类地貌，因此a取值为0.22，湍流度采用与现场实测条件对应的湍流度，即对C类地形的湍流强度模拟，我国现行荷载规范：
Iu＝I10*(z/10)-a,其中，I10为10m高名义湍流度，对应C类地面粗糙度取值为0.23，同时参考日本规范中III类地貌取值：Iu＝0.23，z<＝z0，Iu＝0.1(z/zG)-a-0.05，
z0<＝z<＝zG，式中z0取值为5m，zG为梯度风高度，取值为450m。入口处的湍流动能k和耗散率d分别为：k＝1.5(VzIu)2,d＝a0.75*k1.5/l，式中，Vz表示入口处的平均风速，Iu表示湍流强度，l表示湍流积分尺度。
优选的，所述步骤三基于峰对建筑屋顶檐口的作用为定常，采用SIMPLEC算法加快收敛速度，采用压力速度耦合求解方式下的SIMPLEC方法，高精度的二阶迎风格式并配合较小的松弛因子。
优选的，所述步骤三的压力松弛因子取为0.2，其余欠松弛因子取值为0.5。
优选的，所述步骤三采用稳态计算方法，设置残差为1*10-4，计算过程中设置监控图像，同时监测计算过程中对应的物理量，当其在迭代过程中保持基本不变，确定计算流域已经进入稳定状态。
本发明的有益效果：经过与实际风洞试验结果的比对，风洞试验和模拟的结果比较吻合，数值模拟可以模拟真实的风环境，构造足尺的建筑屋顶檐口计算模型，按实际的风环境进行模拟，有着广阔的应用前景，与风洞试验相结合，为建筑屋顶檐口的抗风设计提供参考依据。
根据下文结合附图对本发明具体实施例的详细描述，本领域技术人员将会更加明了本发明的上述以及其他目的、优点和特征。
附图说明
后文将参照附图以示例性而非限制性的方式详细描述本发明的一些具体实施例。附图中相同的附图标记标示了相同或类似的部件或部分。本领域技术人员应该理解，这些附图未必是按比例绘制的。本发明的目标及特征考虑到如下结合附图的描述将更加明显，附图中：
附图1为根据本发明实施例的建筑屋顶檐口的风效应数值模拟方法流程图；
附图2为为根据本发明实施例的建筑屋顶檐口的实际风洞试验结果与风效应数值模拟结果比对图。其中左侧为建筑屋顶檐口的实际风洞试验结果，右侧为风效应数值模拟结果。
具体实施方式
参见图1，一种建筑屋顶檐口的风效应数值模拟方法，包括：
步骤一，建立计算模型，位于大气边界层中的建筑物，处于一个完全开口的流动风场中，风对建筑物的影响具有一定的范围，为方便计算，在数值模拟时给定有限的三维计算区域，在确定计算流域时，保证计算域足够大的同时，保证边界与模型之间具有足够的距离，同时考虑到随着计算域的增大网格数量增加、计算量加大，因此，数值风洞对阻塞度的要求，在建筑屋顶檐口情况下表示为：(屋顶檐口最大的迎风面积/流域横截面面积)*100％<3％。
除了保证阻塞度要求，还需要考虑流场是否得到充分的发展，对计算区域的尺寸要求既能够消除其对模型附近流场的影响，又不会对计算量产生巨大影响。因此本实施例中计算区域尺寸设计为700m*400m*200m，将建筑屋顶檐口移到计算流域沿风流向的前1/3内，出口边界和建筑物后有着远大于建筑物高度的10倍的距离，保证了充分流动发展状态。
步骤二，进行网格划分并设定边界条件，本实施例中，对建筑物表面及周围流场配置合理的计算网格，相对于理论上讨论计算效率和可靠度更加重要，建筑屋细部形状中对计算结果影响较大的位置是拐角区域及建筑屋壁面，从而在这些区域加密网格，根据湍流模型的不同，当计算结果随着网格的分辨率变化时，可取最小网格宽度为建筑物宽度及进深的1/50，对于流场变化较少的区域，可以减少网格提高计算效率，本实施例中为了保证计算的精度，同时考虑到计算效率，对整个计算流域分块进行网格划分，对于檐口周围的比较近的风敏感区湍流复杂，且体型复杂，采用适应性较强的非结构化四面体网格进行局部加密，远离檐口的区域，采用四面体单元离散，远离建筑屋顶檐口逐渐稀疏，9种工况下网格划分时的最小尺寸以及网格总数如下表1所示。
表1计算模型网格划分
对于边界条件的设定，模拟计算中假定气流为定常的不可压缩流动，进流面采用速度入口，其速度随高度的变化为UZ＝U0(Z/Z0)a，其中，Z0和U0分别为参考高度和该高度处的平均风速，本实施例中参考高度为20m，平均风速为8m/s，Z和UZ分别为计算高度和该高度处的平均风速，a为地面粗糙度指数，本实施例处于C类地貌，因此a取值为0.22。湍流度采用与现场实测条件对应的湍流度，即对C类地形的湍流强度模拟，我国现行荷载规范：Iu＝I10*(z/10)-a,其中，I10为10m高名义湍流度，对应C类地面粗糙度取值为0.23，同时参考日本规范中III类地貌取值：Iu＝0.23，z<＝z0，Iu＝0.1(z/zG)-a-0.05，z0<＝z<＝zG，式中z0取值为5m，zG为梯度风高度，取值为450m。入口处的湍流动能k和耗散率d分别为：
k＝1.5(VzIu)2,d＝a0.75*k1.5/l，式中，Vz表示入口处的平均风速，Iu表示湍流强度，l表示湍流积分尺度。入口边界处速度与湍流强度采用与Fluent接口的UDF(UseDefined Functions)编程，在Fluent中输入恰当的边界条件，数值风洞的边界有入口、出口、地面、壁面以及建筑物壁面，入口边界条件包括风速、湍流动能k和耗散率d，通过Fluent模块提供的用户自定义模块UDF自编边界条件函数myudf.c与Fluent接口实现，边界条件以及参数参见表2。
表2模型参数信息
步骤三，进行求解，求解计算采用3D单精度求解器为压力基隐式。基于峰对建筑屋顶檐口的作用为定常，采用SIMPLEC算法加快收敛速度，采用压力速度耦合求解方式下的SIMPLEC方法，并配合较小的欠松弛因子。合适的离散格式加速收敛，并提高求解速度和解的精度，采用高精度离散方法可以得到较高精度的解，但是容易使数值计算的结果失稳和发散，因而在计算中同时要选择较小的松弛因子，在需要对建筑屋顶檐口进行定量分析的情况下，选择高精度的二阶迎风格式并配合较小的松弛因子，本实施例经过试算后压力松弛因子取为0.2，其余值为0.5。
采用软件Fluent进行计算，采用稳态计算方法，设置残差为1*10-4，计算过程中设置监控图像，同时监测计算过程中对应的物理量，当其在迭代过程中保持基本不变，确定计算流域已经进入稳定状态。本实施例中当所有控制方程的迭代残差量都小于1*10-4，且建筑物屋顶檐口面某点的风压值趋于平缓时，计算结果收敛，迭代次数为1000次；
步骤四，对建筑物屋顶檐口垂线位置对应的建筑物部分门窗的开敞或突然开孔导致的内压增大对具有檐口的屋顶结构进行风效应数值模拟，将模拟位置转移到建筑物屋顶檐口垂线位置对应的建筑物部分门窗处，建立内压传递方程，采用Schinozuka-Deodatis法模拟得到脉动风速的功率谱，对开孔结构的内压传递方程进行推导。
对于固定内部容积建筑，随着开孔面积的减小，孔口等效阻尼比增大，当建筑物开孔面积一定时，结构孔口等效阻尼比随着建筑内部容积的增大而增大，当建筑的开孔面积与建筑内部容积一定时，孔口等效阻尼比随着参考风速U10以及空口处风压高度变化系数的增大而增大，当开孔结构的开孔率较大时，会发生Helmholtz共振，内压的脉动量也会随之增大，随着开孔率的减小，建筑屋顶檐口外部风压脉动的高频部分经过洞口后得到衰减，系统的阻尼增大，内压脉动量也随之减小，最终Helmholtz共振消失。对内部容积不同的建筑，当开孔率相同时，对应的孔口阻尼特性极为相似，仅Helmholtz频率有所差异。紊流中，开孔结构内压瞬态响应的极值随着时间的推移而迅速下降，但下降的趋势并不是严格随着背立面幕墙开孔率的增大而下降，在背立面与迎风面幕墙开孔率的比值约为0.25时，即表1的工况3，内压瞬态响应达到最大，当比值大于或者小于这个值时，内压振荡逐渐消失，当建筑背立面开孔率一定时，随着迎风面幕墙开孔率的增大，内压系数呈现增大的趋势，迎风面幕墙开孔率的减小能减小内压响应的峰值，内压系数时程整体向负风压移动且内压脉动量减小。
参见图2，经过与实际风洞试验结果的比对，风洞试验和模拟的结果比较吻合，数值模拟可以模拟真实的风环境，构造足尺的建筑屋顶檐口计算模型，按实际的风环境进行模拟，有着广阔的应用前景，与风洞试验相结合，为建筑屋顶檐口的抗风设计提供参考依据。
虽然本发明已经参考特定的说明性实施例进行了描述，但是不会受到这些实施例的限定而仅仅受到附加权利要求的限定。本领域技术人员应当理解可以在不偏离本发明的保护范围和精神的情况下对本发明的实施例能够进行改动和修改。

Claims (10)

1.一种建筑屋顶檐口的风效应数值模拟方法，其特征在于包括：
步骤一，建立计算模型，在数值模拟时给定有限的三维计算区域，确定计算流域；
步骤二，进行网格划分并设定边界条件；
步骤三，进行求解，求解计算采用3D单精度求解器为压力基隐式；
步骤四，对建筑物屋顶檐口垂线位置对应的建筑物部分门窗的开敞或突然开孔导致的内压增大对具有檐口的屋顶结构进行风效应数值模拟，将模拟位置转移到建筑物屋顶檐口垂线位置对应的建筑物部分门窗处，建立内压传递方程，采用Schinozuka-Deodatis法模拟得到脉动风速的功率谱，对开孔结构的内压传递方程进行推导。
2.根据权利要求1所述的一种建筑屋顶檐口的风效应数值模拟方法，其特征在于：所述步骤一中确定所述计算流域对应的数值风洞对阻塞度的要求，在建筑屋顶檐口情况下表示为：(屋顶檐口最大的迎风面积/流域横截面面积)*100％<3％，保证计算域足够大的同时，保证边界与模型之间具有足够的距离，同时考虑到随着计算域的增大网格数量增加、计算量加大。
3.根据权利要求2所述的一种建筑屋顶檐口的风效应数值模拟方法，其特征在于：所述步骤一中确定所述计算流域设计还包括使得流场得到充分的发展的计算流域范围。
4.根据权利要求2-3任一所述的一种建筑屋顶檐口的风效应数值模拟方法，其特征在于：所述步骤一重所述计算区域尺寸设计为700m*400m*200m，将建筑屋顶檐口移到计算流域沿风流向的前1/3内，出口边界和建筑物后有着远大于建筑物高度的10倍的距离。
5.根据权利要求1所述的一种建筑屋顶檐口的风效应数值模拟方法，其特征在于：所述步骤二的网格划分包括：建筑屋细部形状中对计算结果影响较大的位置是拐角区域及建筑屋壁面区域加密网格，根据湍流模型的不同，当计算结果随着网格的分辨率变化时，取最小网格宽度为建筑物宽度及进深的1/50。
6.根据权利要求5所述的一种建筑屋顶檐口的风效应数值模拟方法，其特征在于：所述步骤二的网格划分包括：对整个计算流域分块进行网格划分，对于檐口周围的比较近的风敏感区湍流复杂，且体型复杂，采用适应性较强的非结构化四面体网格进行局部加密，远离檐口的区域，采用四面体单元离散，远离建筑屋顶檐口逐渐稀疏。
7.根据权利要求1所述的一种建筑屋顶檐口的风效应数值模拟方法，其特征在于：所述步骤二的边界条件的设定，模拟计算中假定气流为定常的不可压缩流动，进流面采用速度入口，其速度随高度的变化为UZ＝U0(Z/Z0)a，其中，Z0和U0分别为参考高度和该高度处的平均风速，本实施例中参考高度为20m，平均风速为8m/s，Z和UZ分别为计算高度和该高度处的平均风速，a为地面粗糙度指数，本实施例处于C类地貌，因此a取值为0.22，湍流度采用与现场实测条件对应的湍流度，即对C类地形的湍流强度模拟，我国现行荷载规范：Iu＝I10*(z/10)-a,其中，I10为10m高名义湍流度，对应C类地面粗糙度取值为0.23，同时参考日本规范中III类地貌取值：Iu＝0.23，z<＝z0，Iu＝0.1(z/zG)-a-0.05，z0<＝z<＝zG，式中z0取值为5m，zG为梯度风高度，取值为450m。入口处的湍流动能k和耗散率d分别为：k＝1.5(VzIu)2,d＝a0.75*k1.5/l，式中，Vz表示入口处的平均风速，Iu表示湍流强度，l表示湍流积分尺度。
8.根据权利要求1所述的一种建筑屋顶檐口的风效应数值模拟方法，其特征在于：所述步骤三基于峰对建筑屋顶檐口的作用为定常，采用SIMPLEC算法加快收敛速度，采用压力速度耦合求解方式下的SIMPLEC方法，高精度的二阶迎风格式并配合较小的松弛因子。
9.根据权利要求8所述的一种建筑屋顶檐口的风效应数值模拟方法，其特征在于：所述步骤三的压力松弛因子取为0.2，其余欠松弛因子取值为0.5。
10.根据权利要求9所述的一种建筑屋顶檐口的风效应数值模拟方法，其特征在于：所述步骤三采用稳态计算方法，设置残差为1*10-4，计算过程中设置监控图像，同时监测计算过程中对应的物理量，当其在迭代过程中保持基本不变，确定计算流域已经进入稳定状态。
CN201811401246.4A 2018-11-22 2018-11-22 一种建筑屋顶檐口的风效应数值模拟方法 Pending CN109543302A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811401246.4A CN109543302A (zh) 2018-11-22 2018-11-22 一种建筑屋顶檐口的风效应数值模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811401246.4A CN109543302A (zh) 2018-11-22 2018-11-22 一种建筑屋顶檐口的风效应数值模拟方法

Publications (1)

Publication Number Publication Date
CN109543302A true CN109543302A (zh) 2019-03-29

Family

ID=65850324

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811401246.4A Pending CN109543302A (zh) 2018-11-22 2018-11-22 一种建筑屋顶檐口的风效应数值模拟方法

Country Status (1)

Country Link
CN (1) CN109543302A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110196990A (zh) * 2019-04-19 2019-09-03 中国人民解放军63921部队 大跨超高空间网格结构风工程分析方法
CN110489854A (zh) * 2019-08-14 2019-11-22 中船第九设计研究院工程有限公司 一种单墙开洞建筑的内压估算简化方法
CN111027201A (zh) * 2019-12-04 2020-04-17 北京工业大学 一种考虑潜在开孔的低矮建筑围护结构风致损失估计方法
CN112528571A (zh) * 2020-11-11 2021-03-19 武汉理工大学 适用于近地面湍流特征的数值风洞模拟的建模方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003203194A (ja) * 2001-08-17 2003-07-18 Ohbayashi Corp 風環境予測プログラム、このプラグラムを格納した媒体および風環境予測方法
CN104200102A (zh) * 2014-09-03 2014-12-10 中国建筑设计咨询有限公司 一种基于cfd的光伏阵列风荷载预测方法
WO2015122756A1 (es) * 2014-02-13 2015-08-20 Dirección General De Educación Superior Tecnológica De La Sep Extractor mixto de aire, que utiliza la energía solar y la del viento, para ventilar pasivamente naves industriales, viviendas, y/o edificaciones (chimenea solar-eólica)
CN106227951A (zh) * 2016-07-26 2016-12-14 南京航空航天大学 适用于超规范特大型冷却塔的二维风振系数计算方法
CN106250590A (zh) * 2016-07-21 2016-12-21 哈尔滨工业大学 一种基于cfd数值模拟的高空露天平台行人风环境评估方法
CN107300853A (zh) * 2017-06-26 2017-10-27 中船第九设计研究院工程有限公司 一种分区计算空间相关性的脉动风速时程输入方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003203194A (ja) * 2001-08-17 2003-07-18 Ohbayashi Corp 風環境予測プログラム、このプラグラムを格納した媒体および風環境予測方法
WO2015122756A1 (es) * 2014-02-13 2015-08-20 Dirección General De Educación Superior Tecnológica De La Sep Extractor mixto de aire, que utiliza la energía solar y la del viento, para ventilar pasivamente naves industriales, viviendas, y/o edificaciones (chimenea solar-eólica)
CN104200102A (zh) * 2014-09-03 2014-12-10 中国建筑设计咨询有限公司 一种基于cfd的光伏阵列风荷载预测方法
CN106250590A (zh) * 2016-07-21 2016-12-21 哈尔滨工业大学 一种基于cfd数值模拟的高空露天平台行人风环境评估方法
CN106227951A (zh) * 2016-07-26 2016-12-14 南京航空航天大学 适用于超规范特大型冷却塔的二维风振系数计算方法
CN107300853A (zh) * 2017-06-26 2017-10-27 中船第九设计研究院工程有限公司 一种分区计算空间相关性的脉动风速时程输入方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
于勇: "《FLUENT入门与进阶教程》", 《北京理工大学出版社》, pages: 22 *
周绪红等: "低层双坡屋面建筑三维定常风场的数值模拟", 《工程力学》, no. 03, 25 March 2010 (2010-03-25), pages 19 - 29 *
武岳: "《风工程与结构抗风设计》", 31 May 2014, pages: 34 - 38 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110196990A (zh) * 2019-04-19 2019-09-03 中国人民解放军63921部队 大跨超高空间网格结构风工程分析方法
CN110489854A (zh) * 2019-08-14 2019-11-22 中船第九设计研究院工程有限公司 一种单墙开洞建筑的内压估算简化方法
CN110489854B (zh) * 2019-08-14 2023-08-01 中船第九设计研究院工程有限公司 一种单墙开洞建筑的内压估算简化方法
CN111027201A (zh) * 2019-12-04 2020-04-17 北京工业大学 一种考虑潜在开孔的低矮建筑围护结构风致损失估计方法
CN111027201B (zh) * 2019-12-04 2023-09-19 北京工业大学 一种考虑潜在开孔的低矮建筑围护结构风致损失估计方法
CN112528571A (zh) * 2020-11-11 2021-03-19 武汉理工大学 适用于近地面湍流特征的数值风洞模拟的建模方法

Similar Documents

Publication Publication Date Title
CN109543302A (zh) 一种建筑屋顶檐口的风效应数值模拟方法
Zheng et al. CFD simulations of wind flow and mean surface pressure for buildings with balconies: Comparison of RANS and LES
Vasaturo et al. Large eddy simulation of the neutral atmospheric boundary layer: performance evaluation of three inflow methods for terrains with different roughness
Parente et al. A comprehensive modelling approach for the neutral atmospheric boundary layer: consistent inflow conditions, wall function and turbulence model
Parsons et al. Numerical modelling of flow structures over idealized transverse aeolian dunes of varying geometry
CN106682262B (zh) 一种获取飞行器流场的数值模拟方法
Balendra et al. Evaluation of flow characteristics in the NUS-HDB Wind Tunnel
CN114357571A (zh) 建成建筑环境下大气边界层风场特性的反演方法和系统
Pattanapol et al. Modeling of surface roughness for flow over a complex vegetated surface
Emmons Critique of numerical modeling of fluid-mechanics phenomena
Averkova et al. Numerical simulation of air currents at the inlet to slot leaks of ventilation shelters
Liu et al. Large-eddy simulations of the flow over an isolated three-dimensional hill
True et al. Wind driven cross-flow analysed as a catchment problem and as a pressure driven flow
Zheng et al. Evaluation and comparison of various fast fluid dynamics modeling methods for predicting airflow around buildings
Baghapour et al. Numerical simulation of wind-induced turbulence over precipitation gauges
Wright et al. Unsteady CFD simulations for natural ventilation
Kewalramani et al. Energy consistent Gaussian integral model for jet with off-source heating
Nardecchia et al. CFD analysis of urban canopy flows employing the V2F model: Impact of different aspect ratios and relative heights
Yoshikawa et al. LES for wind load estimation by unstructured grid system
Kevlahan et al. Adaptive wavelet simulation of global ocean dynamics using a new Brinkman volume penalization
Roy et al. Atmospheric boundary layer airflow through CFD simulation on pyramidal roof of square plan shape buildings
CN111046542B (zh) 一种评估二十面体的le网格计算特性方法及离散方法
Hanson et al. Numerical modelling of wind flow over buildings in two dimensions
Vasconcellos et al. Flow past a circular cylinder: a comparison between commercial finite volume and finite element codes
Djilali An investigation of two-dimensional flow separation with reattachment

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