摘要
对于气动性能,钝体断面的气动外形非常重要,采用传统风洞试验及CFD模拟计算得到钝体断面气动性能需消耗大量时间,大大影响钝体断面气动外形的气动性能评估效率。通过卷积神经网络深度学习技术实现对气动性能的快速预测,深度学习模型训练完成后,输入形状信息和与形状相关的流场信息,即可输出不同几何形状下的阻力系数,进而得到钝体断面的气动性能。为寻找性能最优的深度学习模型,通过综合判定误差和参数量大小对卷积神经网络结构的深度和宽度进行优化。对深度学习模型输出阻力系数与CFD计算结果进行对比发现,误差符合预期要求,并且相较于传统方法,基于深度学习网络的预测所需时间达到数量级的提升,未来可作为钝体断面气动外形优化的关键方法。
平均风产生的静荷载简称静力风荷载,通常由无量纲参数三分力系数来描述,而过大的阻力和升力矩会让钝体断面产生横向位移和扭转,严重时甚至会出现风致静力失稳,进而影响结构的安全性和稳定性,由此可见,三分力系数对钝体断面抗风设计至关重
目前测定三分力系数的常用方法主要有两种,1)将节段模型放入风洞进行试验,得到三分力系
近年来,随着神经网络技术的不断发展,人工智能预测技术逐渐应用在抗风领
深度学习是人工智能领域近年来取得的最重大的研究成果之一,由于出色的数据拟合能力、GPU等硬件的出现,其进入高速发展阶段。卷积神经网络(CNN)可以对图像特征进行自动识别,且卷积核平移不变性,可以极大地降低训练难度,非常适合处理不同形状钝体断面图
深度学习中的神经网络模型输入需能够表达钝体断面的气动外形信息,且该信息应该具有一致的表达结构。对于钝体断面气动外形的表达问题,使用坐标点的传统方式会使输入数据的结构和大小与外形相关联,引发数据异构问题。针对该问题,提出一种图像形式的气动外形描述方法,可以有效地与卷积神经网络结合,实现输入数据一致性表达的同时,提高了信息传递效率。基于此方法,以3个场对气动外形进行描述,分别是距离场、顺风向坐标场和横风向坐标场。
外形影响气动力的最主要因素之一是流场的边界层。采用距离场的形式描述外形,将外形以周围空间各点与外形的距离进行表达,可以强化输入信息与边界层的逻辑联系。具体表达如
(1) |
式中:指数函数内部表示流场空间位置与气动外形最近边界之间的无量纲距离;B为气动外形特征长度;U为平均风速,m/s。因为更加关注靠近气动外形边界周围的流场,忽略了离气动外形较远的流场,需要输入更多的边界周围信息,为此,通过负指数形式来增加边界层附近的权重。采样位置在边界上或者边界内时,为0,反之为1。

图1 距离场
Fig. 1 Distance field
单纯采用距离场作为输入存在信息缺失的问题,为了方便描述,考虑如

图2 气动外形示意图
Fig. 2 Schematic diagram of aerodynamic shape
为兼顾输入信息获取的便携性和流场信息的必要性,设置一个流场流动方向的图像数据和一个描述横风向方向的图像数据作为输入数据,即顺风向坐标场和横风向坐标场。
其中,顺风向坐标场见
(2) |
式中:X为空间点的X方向坐标,原点为形状在顺风向投影的中心;采样位置在边界上或者边界内时,为0,反之为1。

图3 顺风向坐标场
Fig. 3 Coordinate field of streamwise wind
此外,横风向坐标场见
(3) |
式中:Y为空间点的Y方向坐标;采样位置在边界上或者边界内时,为0,反之为1。

图4 横风向坐标场
Fig. 4 Cross-wind direction coordinate field
输入为3个图像,按照式(1)~

图5 设计的CNN结构图
Fig. 5 Design of CNN structure diagram
注: D为初始宽度。
输入层为包含距离信息和流场信息的图像。输入图像分辨率越高,模型预测越准确,但模型复杂度和训练时间也会增加,因此,选择合适分辨率非常重要。一般输入层大小需要能够被2整除多次,以保证卷积层深度。选择64×64像素为输入分辨率,有3个场作为输入,即64×64×3的三维矩阵,每个场取值均确保在0到1之间。输出层输出数据用作回归分析,对钝体断面阻力系数进行预测,为保证输出数据范围在0到1之间,还加入了sigmoid函数,见
(4) |
隐藏层中:卷积层是卷积神经网络的核心模块,主要确定卷积核尺寸、数量、步长、零填充的数量,卷积核尺寸有越来越小的趋势,因为小尺寸卷积核可以通过堆叠替代大尺寸卷积层并且感受野大小不
CNN框架搭建完成后,输入和输出数据的获取也非常关键,输入数据为3个类图片的数据场,输出数据为阻力系数,获取不同钝体断面气动外形下的输入和输出数据,为深度学习模型提供数据支持非常重要。
按照基本形状,分为正六边形、圆形、矩形、正菱形4种,为降低模型训练难度,所有形状做到上下和左右对称。根据横向长度与竖向长度比值的不同,每种基本形状下又有不同形状,为方便起见,采用无量纲形式来确定形状,固定所有形状的横向长度为1。考虑到正六边形非常接近钝体断面,因此,正六边形输入数据占到总数据的57%,增强了模型对钝体断面的学习能力。
根据正六边形上下角点距离X=0直线垂直距离的不同,分成4组,圆形、矩形、正菱形各1组,共7组,每组15个不同形状,共105个,如
基本形状 | ![]() | ![]() | ![]() | ![]() |
---|---|---|---|---|
上角点到X=0的距离R | [0.1,0.2,0.3,0.4] | |||
B | 1 | 1 | 1 | 1 |
AR=B/H | [0.5, 0.56, 0.625, 0.71, 0.83, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10] | |||
数量 | 60 | 15 | 15 | 15 |
注: AR为横向长度与竖向长度的比值。
输入数据获取需先生成64×64的网格位置坐标,再根据式(1)~
输出数据为阻力系数,采用CFD数值模拟计算得到稳定后的压力流场。计算域为矩形,长度和宽度分别为15B和10B,B为钝体截面宽度。为了防止网格变化速度太快,设置3个计算区域,定义刚性区域,以保持网格形状不变,在刚性区外使用变形区,实现网格平滑处理,最后,使用固定区域对计算域中的其他部分进行网格划分。
采用Fluent软件,湍流模型选择RANs模型,具体选择k-ω SST模型,在近壁面采用k-ω模型,边界层外缘和自由剪切层采用k-模型,在k-ω模型和k-模型间采用一个混合函数过渡。入口边界设置为速度入口,出口设置为压力出口,形状表面采用无滑移条件,这意味着物理表面的流动状态等于形状的运动。待流场稳定后,读取数据,通过对稳定压力流场进行数据处理,得到阻力系数。
将全部105个不同形状按照上述操作进行处理,得到对应形状的输入和输出数据。综上所述,输入数据的图像张量为(84, 3, 64, 64),输出数据张量为(21,3),快速得到输入和输出数据,为之后的模型训练做好准备。
深度学习模型的深度和宽度直接影响模型性能,模型加深能一定程度提升模型性能,但网络深度也不是越深越好,深层网络容易带来梯度不稳定、网络退化等问题,且无法通过优化消除,导致模型网络深度不断加深,性能并未出现明显提升,甚至出现下降的情况,模型宽度同理,随着模型宽度的加深,模型性能并不一定显著提升。基于之前设计的卷积神经网络,寻找最优的模型深度与宽度,主要通过模型预测阻力系数的平均绝对误差MAE,见
(5) |
式中:为阻力系数真实值;为阻力系数预测值。
与模型宽度相比,模型深度对于模型性能的影响更大,因此,先确定最优的模型深
CNN层数 | MAE | 参数量 |
---|---|---|
CNN_2 | 0.314 1 | 8 270 419 |
CNN_3 | 0.238 3 | 4 469 587 |
CNN_4 | 0.291 1 | 3 601 747 |
CNN_5 | 0.246 0 | 7 297 363 |
CNN_6 | 0.271 1 | 25 661 779 |
由
在确定深度学习模型层数为3层的情况下,优化深度学习模型宽度D,主要从模型预测误差大小和模型参数量来进行评估。基于CNN采样,第1层宽度为64,选择对比的CNN网络第1层宽度分别为128和256,第2层和第3层模型宽度分别为前一层模型宽度的2倍,采用完全相同数据处理后的输入和输出数据,模型参数配置完全相同,结果见
模型宽度 | MAE | 参数量 |
---|---|---|
CNN_64 | 0.276 7 | 4 469 587 |
CNN_128 | 0.196 4 | 9 673 683 |
CNN_256 | 0.248 3 | 22 293 715 |
由
综上所述,通过MAE大小和模型参数量大小对比,选择CNN层数为3层的网络深度,网络宽度为128-256-512的网络模型,第1层全连接层中神经元数量S为32 768个。
Pytorch在科学研究方面表现非常优秀,主要体现在Pytorch风格非常Python化,降低了入门的难度,且搭建深度学习模型时可以逐层搭建,方便实时修改。鉴于搭建模型时修改细节比较多,因此,使用Pytorch来搭建深度学习模
在所使用的105个数据中选取80%作为训练数据集,剩下的作为验证数据集,训练数据集共有84个不同形状的输入和输出,验证数据集有21个不同形状的输入和输出,每次训练前都对样本数据进行随机打乱和批量读取。模型训练方法选择SGD,使用0.2作为训练模型的Dropout率,以防止过拟合。
通过卷积神经网络预测阻力系数的结果如

图6 阻力系数预测结果
Fig 6 Drag coefficient prediction results
形状 | 平均相对误差/% |
---|---|
R=0.1正六边形 | 38.75 |
R=0.2正六边形 | 6.98 |
R=0.3正六边形 | 9.39 |
R=0.4正六边形 | 21.49 |
圆形 | 12.39 |
矩形 | 31.77 |
正菱形 | 38.81 |
由
深度学习模型针对2个截面预测的平均相对误差在10%以下,针对3个截面预测的平均相对误差在30%以上,主要是模型输入和输出的逻辑关系过于薄弱,导致模型性能容易出现反复。
提出使用卷积神经网络来预测阻力系数是为了减少得到钝体断面气动性能的时间,因此,有必要评估深度学习方法的时间性能。
通过编写好的代码,输入设计约1 s即可完成,输出设计需要进行网格划分、服务器计算等一系列操作,得到稳定流场后,计算得到阻力系数约需10 s。本文卷积神经网络进行一次迭代训练的平均时间为1.08 s,进行1 000次训练迭代所需的总时间约为1 080 s,使用已经训练好的模型进行一次回归预测,即产生阻力系数所需的平均时间为1 s,一共需要1 092 s。
相比于传统基于CFD计算和风洞试验需要几周时间才能得到三分力系数,基于深度学习的预测方法在时间性能上得到了数量级的提升。
提出将不同钝体断面的气动外形信息转化为图像数据而不依赖于传统设计参数,使用修改和优化后的卷积神经网络对气动外形图像进行回归预测,建立钝体断面气动外形和阻力系数之间的联系,且预测精度相对较高。针对相对误差较大的钝体断面,后续将尝试采用新的输出,强化输入和输出逻辑关系的同时,针对深度学习模型来进行改善,进一步提升阻力系数的预测精度。
与传统的风洞试验和CFD仿真计算方法相比,在允许误差范围内,提出的卷积神经网络模型计算时间达到了数量级的提升,为钝体断面气动外形优化提供了计算基础。
参考文献
陈政清. 桥梁风工程[M]. 北京: 人民交通出版社, 2005. [百度学术]
CHEN Z Q. Bridge wind engineering [M]. Beijing: China Communications Press, 2005. (in Chinese) [百度学术]
欧阳克俭, 陈政清, 李红利, 等. 空间主缆悬索桥非线性静风稳定分析[J]. 土木建筑与环境工程, 2010, 32(3): 83-87. [百度学术]
OUYANG K J, CHEN Z Q, LI H L, et al. Nonlinear aerostatic stability of suspension bridge with spatial main cables [J]. Journal of Civil, Architectural & Environmental Engineering, 2010, 32(3): 83-87. (in Chinese) [百度学术]
廖海黎, 李明水, 马存明, 等. 桥梁风工程2019年度研究进展[J]. 土木与环境工程学报(中英文), 2020, 42(5): 56-66. [百度学术]
LIAO H L, LI M S, MA C M, et al. State-of-the-art review of bridge wind engineering in 2019 [J]. Journal of Civil and Environmental Engineering, 2020, 42(5): 56-66. (in Chinese) [百度学术]
高亮. 桥梁断面三分力系数试验研究[D]. 西安: 长安大学, 2010. [百度学术]
GAO L. Investigation into three-component force coefficients of bridge deck section [D]. Xi,an: Chang,an University, 2010. (in Chinese) [百度学术]
张丹. 流线型桥梁断面三分力系数及表面压力系数雷诺数效应[D]. 西安: 长安大学, 2013. [百度学术]
ZHANG D. Reynolds number effect of streamline-like bridge deck section static coeffcient and surface pressure [D]. Xi,an: Chang,an University, 2013. (in Chinese) [百度学术]
刘昊苏, 雷俊卿. 大跨度双层桁架主梁三分力系数识别[J]. 浙江大学学报(工学版), 2019, 53(6): 1092-1100. [百度学术]
LIU H S, LEI J Q. Identification of three-component coefficients of double deck truss girder for long-span bridge [J]. Journal of Zhejiang University (Engineering Science), 2019, 53(6): 1092-1100. (in Chinese) [百度学术]
宋锦忠, 徐建英. 桥梁节段模型风洞试验技术研讨[C]//2004全国结构风工程实验技术研讨会论文集. 长沙, 2004: 68-70. [百度学术]
SONG J Z. XU J Y. Discussion on wind tunnel test technology of bridge segment model [C]//Proceedings of the 2004 National Symposium on Structural Wind Engineering Experimental Technology. Changsha, 2004: 68-70. (in Chinese) [百度学术]
晏致涛, 李正良, 祁仁俊. 大跨度中承式拱桥节段模型风洞试验[J]. 重庆建筑大学学报, 2006, 28(6): 55-57, 98. [百度学术]
YAN Z T, LI Z L, QI R J. Wind tunnel tests of long span half-through arch bridge section model [J]. Journal of Chongqing Jianzhu University, 2006, 28(6): 55-57, 98. (in Chinese) [百度学术]
王锋, 赖笑. 基于CFD桥梁典型断面三分力系数差异研究[J]. 建筑技术开发, 2017, 44(1): 131-132. [百度学术]
WANG F, LAI X. Study on difference of typical section of three component coefficients based on CFD bridge [J]. Building Technology Development, 2017, 44(1): 131-132. (in Chinese) [百度学术]
刘婷婷, 张文首. CFD方法的大跨度桥梁抖振荷载[J]. 土木建筑与环境工程, 2011, 33(6): 59-63. [百度学术]
LIU T T, ZHANG W S. Load analysis of butteting bridges by CFD [J]. Journal of Civil, Architectural & Environmental Engineering, 2011, 33(6): 59-63. (in Chinese) [百度学术]
勾红叶, 杨彪, 华辉, 等. 桥梁信息化及智能桥梁2019年度研究进展[J]. 土木与环境工程学报(中英文), 2020, 42(5): 14-27. [百度学术]
GOU H Y, YANG B, HUA H, et al. State-of-the-art review of bridge informatization and intelligent bridge in 2019 [J]. Journal of Civil and Environmental Engineering, 2020, 42(5): 14-27. (in Chinese) [百度学术]
李乔, 杨兴旺. 桥梁主梁断面空气力学特性分析的人工神经网络方法[J]. 中南公路工程, 2000, 25(3): 56-59. [百度学术]
LI Q, YANG X W. ANN method for analysis of aerodynamic characteristics of bridge girder section [J]. Central South Highway Engineering, 2000, 25(3): 56-59. (in Chinese) [百度学术]
杨兴旺. 人工神经网络在桥梁主梁断面空气力学特性分析中的应用[D]. 成都: 西南交通大学, 2000. [百度学术]
YANG X W. Application of artificial neural network in analysis of aerodynamic characteristics of bridge main girder section [D]. Chengdu: Southwest Jiaotong University, 2000. (in Chinese) [百度学术]
黄继鸿, 苏红莲, 赵新华. 基于BP神经网络的翼型空气动力系数预测[J]. 航空工程进展, 2010, 1(1): 36-39. [百度学术]
HUANG J H, SU H L, ZHAO X H. Aerodynamic coefficient prediction of airfoil using BP neural network [J]. Advances in Aeronautical Science and Engineering, 2010, 1(1): 36-39. (in Chinese) [百度学术]
陈讷郁, 葛耀君. 基于人工神经网络的典型桥梁断面气动参数识别[J]. 土木工程学报, 2019, 52(8): 91-97, 128. [百度学术]
CHEN N Y, GE Y J. Aerodynamic parameter identification of typical bridge sections based on artificial neural network [J]. China Civil Engineering Journal, 2019, 52(8): 91-97, 128. (in Chinese) [百度学术]
LIAO H L, MEI H Y, HU G, et al. Machine learning strategy for predicting flutter performance of streamlined box girders [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2021, 209: 104493. [百度学术]
HU G, KWOK K C S. Predicting wind pressures around circular cylinders using machine learning techniques [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2020, 198: 104099. [百度学术]
HU G, LIU L B, TAO D C, et al. Deep learning-based investigation of wind pressures on tall building under interference effects [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2020, 201: 104138. [百度学术]
CHEN D L, GAO X, XU C F, et al. FlowGAN: A conditional generative adversarial network for flow prediction in various conditions [C]//2020 IEEE 32nd International Conference on Tools with Artificial Intelligence. November 9-11, 2020, Baltimore, MD, USA. IEEE, 2020: 315-322. [百度学术]
MIYANAWALA T P, JAIMAN R K. A novel deep learning method for the predictions of current forces on bluff bodies [C]//Proceedings of ASME 2018 37th International Conference on Ocean, Offshore and Arctic Engineering, June 17-22, 2018, Madrid, Spain. 2018. [百度学术]
陈海, 钱炜祺, 何磊. 基于深度学习的翼型气动系数预测[J]. 空气动力学学报, 2018, 36(2): 294-299. [百度学术]
CHEN H, QIAN W Q, HE L. Aerodynamic coefficient prediction of airfoils based on deep learning [J]. Acta Aerodynamica Sinica, 2018, 36(2): 294-299. (in Chinese) [百度学术]
廖鹏, 姚磊江, 白国栋, 等. 基于深度学习的混合翼型前缘压力分布预测[J]. 航空动力学报, 2019, 34(8): 1751-1758. [百度学术]
LIAO P, YAO L J, BAI G D, et al. Prediction of hybrid airfoil leading edge pressure distribution based on deep learning [J]. Journal of Aerospace Power, 2019, 34(8): 1751-1758. (in Chinese) [百度学术]
KHAN A, SOHAIL A, ZAHOORA U, et al. A survey of the recent architectures of deep convolutional neural networks [J]. Artificial Intelligence Review, 2020, 53(8): 5455-5516. [百度学术]
SZEGEDY C, VANHOUCKE V, IOFFE S, et al. Rethinking the inception architecture for computer vision [C]//2016 IEEE Conference on Computer Vision and Pattern Recognition. June 27-30, 2016, Las Vegas, NV, USA. IEEE, 2016: 2818-2826. [百度学术]
尹宝才, 王文通, 王立春. 深度学习研究综述[J]. 北京工业大学学报, 2015, 41(1): 48-59. [百度学术]
YIN B C, WANG W T, WANG L C. Review of deep learning [J]. Journal of Beijing University of Technology, 2015, 41(1): 48-59. (in Chinese) [百度学术]
ELDAN R, SHAMIR O. The power of depth for feedforward neural networks [C]//Conference on Learning Theory, 2016: ArXiv:1512.03965. [百度学术]
LU Z, PU H, WANG F, et al. The expressive power of neural networks: a view from the width[C]//Neural Information Processing Systems, 2017: 6232-6240. [百度学术]