基于均匀设计及BP神经网络的大体积混凝土热学参数反分析

张玉平　马超　李传习　高树威

摘 要：为了解决桥梁工程大体积混凝土热学参数失真的问题，提出一种基于均匀设计理论与BP神经网络的大体积混凝土热学参数反分析方法。该方法通过BP神经网络建立大体积混凝土温度场与热学参数的非线性关系;BP神经网络的训练样本由均匀设计方法确定;在BP神经网络训练阶段，采用附加动量法对网络结构进行优化;对优化前后的误差曲线及多次训练过程的分析结果表明：附加动量法可明显缩短网络训练时间，多次训练过程的平均绝对百分比误差值及均方根误差值稳定。在太洪长江大桥散索鞍支墩承台大体积混凝土施工中反演了绝热温升、反应速率常数及导热系数，基于反演值的温度计算值与现场实测值吻合较好，温度峰值最大误差仅为1.1 ℃。基于均匀设计理论与BP神经网络的大体积混凝土热学参数反分析方法可行，且反演过程稳定收敛，反演精度高，应用于指导温控施工能降低大体积混凝土开裂风险。

关键词：桥梁工程;大体积混凝土;BP神经网络;均匀设计;热学参数;反分析

中图分类号：U444;U445.57   文献标志码：A   文章编号：2096-6717（2021）02-0148-10

Abstract： In order to solve the problem of distortion of thermal parameters of mass concrete in bridge engineering， an inverse analysis method of thermal parameters of mass concrete based on uniform design theory and BP neural network was proposed. This method uses BP neural network to establish the non-linear relationship between the temperature field and thermal parameters of large-volume concrete; training samples of BP neural network are determined by uniform design method; during the training phase of BP neural network， additional momentum method is used to optimize the network structure; The error curve before and after optimization and the analysis results of multiple training processes show that the additional momentum method can significantly shorten the network training time， and the average absolute percentage error value and root mean square error value of the multiple training processes are stable. The adiabatic temperature rise， reaction rate constant， and thermal conductivity were inverted during the construction of the bulk concrete of the supporting platform of the saddle pier of the Taihong Yangtze River Bridge. The calculated temperature based on the inversion value agrees well with the actual measured value， and the maximum temperature peak error is only 1.1 ℃. Therefore， the inverse analysis method for thermal parameters of large-volume concrete based on uniform design theory and BP neural network is feasible and the inversion process is stable and convergent， and the inversion accuracy is high. It can be used to guide the temperature-controlled construction to reduce the risk of large-scale concrete cracking.

Keywords：bridge engineering;mass concrete;BP neural network;uniform design;thermal parameters;inverse analysis

大體积混凝土作为现代桥梁的重要组成部分，其开裂问题一直是工程技术人员最为关注的问题之一。为了指导大体积混凝土施工，最大限度降低施工阶段大体积混凝土的开裂风险，可根据现场实际施工状况和现有的温度场理论对其进行有限元分析。但是，在对大体积混凝土温度场和应力场进行有限元分析时，所采用的热学参数主要通过经验公式或试验得到，由于经验公式难免有误差，而试验方法花费高、耗时长而较少使用，并且这些参数在施工期受气象条件、时空、荷载、施工条件等多种因素影响，往往使得所采用的热学参数失真，大体积温控计算偏离甚至严重偏离实际情况，导致一些大体积混凝土结构出现或多或少的裂缝[1-2]，轻则影响结构的耐久性，重则产生工程质量事故，危及人民生命财产安全。因此，准确确定热学参数是大体积混凝土温控能否成功的基础和前提，除试验确定外，通过反演分析确定大体积混凝土的热学参数也是行之有效的方法之一[3]。

目前，对大体积混凝土热学参数的反演分析研究不多。王振红等[4]基于施工现场的混凝土立方体非绝热温升试验，对混凝土热学参数进行了反分析。Amirfakhrian等[5]利用径向基函数来获得内热源反分析的数值解。李守巨等[6]基于模糊理论建立了目标函数，将热传导反问题作为非线性优化问题处理。还有许多学者基于蚁群算法、粒子群算法、遗传算法等优化算法提出热学参数反演方法。但是，其中大部分研究还集中在坝体等水利工程中[7-9]，对桥梁工程中大体积混凝土热学参数反分析较少[10-11]。喻正富等[10]基于施工现场测得的温度数据，采用遗传算法对大体积混凝土的热力学参数进行了反演分析，但在反演过程中采取的是分步反演，而非同步反演。文豪等[11]对遗传算法进行优化改进，提出在MATLAB中调用ANSYS温度场数据，反演得到真实的热学参数，但运用MATLAB调用数据，需另写程序代码，比较繁琐。

为获得混凝土真实的热学参数，笔者将BP神经网络引入到温度场热学参数反分析中，BP神经网络特有的强大非线性映射能力可准确映射各热学参数与温度场之间的复杂非线性关系，结合实测数据即可快速、准确地反演混凝土热学参数。此方法有效避免了采用经验公式求解热学参数时产生的误差，还解决了因试验设备昂贵或试验过程漫长而未及时获取热学参数的问题。并且，在建立神经网络的训练样本时引入均匀设计理论，可在一定程度上减少训练样本，提高效率。研究表明：基于均匀设计及BP神经网络的反分析，得到了能反映大体积混凝土真实性能的多种热学参数，使有限元仿真分析结果能准确地反映实际工程情况，及时采取的相应温控措施避免了温度裂缝的产生。

1 均匀设计与BP神经网络的基本原理

1.1 均匀设计理论

均匀设计又称均匀设计试验法，是基于试验点在整个试验范围内均匀散布的，从均匀性角度出发提出的一种试验设计方法[12]，由方开泰教授和数学家王元在1978年共同提出。均匀设计表一般用符号Un（qs）表示，其中U表示均匀表，n表示试验次数，q表示因素水平数，s表示最多可安排的因素数。其中，n=q，s=q-1，即均匀设计表的试验次数n等于水平数q，最多可安排的因素数比水平数少1，而正交设计表中试验次数n等于水平数q的平方，相比于正交设计试验，均匀设计试验可大大减少试验次数。例如，某试验有6个因素，每个因素取31个水平，其全部组合有316=887 503 681，若用正交设计至少需要312=961次试验，而用均匀设计只需31次试验，由此可见，均匀设计较正交设计更适用于多因素多水平试验。

均匀试验设计还具有一个鲜明的特点，即能从全面试验点中挑选出部分代表性的试验点，保证每个因素的每个水平做且仅做一次试验。例如，做2因素11水平的试验，应选用U11（1110）表，表中共有10列，最多可安排10个因素，現在有2个因素，根据U11（1110）的使用表，应取1、7列安排试验，如图1所示。当有4个因素时，应取1、2、5、7列安排试验。笔者采用均匀设计理论确定热学参数样本，在保持样本均匀性不变的前提下，可大幅提高计算效率。

1.2 BP神经网络的基本原理

BP（Back Propagation）神经网络[13]是目前应用最广泛的一种前馈型神经网络模型。在函数逼近、模式识别、分类问题、数据压缩等领域均有广泛应用。BP神经网络由信息的正向传播和误差的反向传播两个过程组成。正向传播中，输入层接收外界信息并向隐含层传播，隐含层负责信息变换，最终传至输出层。当实际输出与期望输出不符时，进入误差的反向传播阶段。反向传递中，按误差梯度下降的方式逐层反传至隐含层、输入层，直到预测输出无限逼近期望输出。BP神经网络结构一般由输入层、一个以上的隐含层、输出层组成。理论证明，3层的BP神经网络（单隐含层）可以实现从输入到输出的任意非线性映射。由图2可见，若输入层有m个神经元，输出层有n个神经元，则可实现m维至n维的映射。

BP神经网络最大的优点就是可以避免函数的具体形式，不必像显式函数那样需要提前确定定义域和值域，所以，在工程领域多应用于岩石力学中的岩石行为预测[14]、边坡位移反分析[15]等问题中。混凝土热学参数反分析同边坡位移反分析一样，均是复杂的非线性问题，难以用显式的函数来描述，而这种复杂的非线性问题可以通过BP神经网络得到较好的映射。

2 BP神经网络反演分析方法

2.1 概述

反分析问题的求解方式一般分为解析反分析法和数值反分析法[16]。解析反分析法主要是通过找出现场监测值与待反演参数之间的显式关系，建立数学模型，进而求得待反演参数。但是，由于实际工程问题的复杂多变性，一般难以确定其数学模型，显式函数关系亦很难建立。因此，针对上述问题，有研究人员提出，建立反分析问题的目标函数，将参数求解问题转化为目标函数最优解问题[17]，采用有限元法等数值方法进行迭代计算，并逐步修正待反参数，直至函数寻到最优解。但当反演参数比较多时操作起来将非常费时，因为寻优过程需要反复迭代，且对复杂问题搜索到最优解的概率较低。BP神经网络为反分析问题提供了一条有效的途径。其一，BP神经网络由其特有的拓扑结构来表述反分析中的复杂非线性问题，无需建立数学模型;其二，BP神经网络解决问题时，只需事先提供训练样本，完成训练即可，无需反复迭代。

2.2 大体积混凝土温度场影响因素分析

在温度场分析中，主要热学参数包括：混凝土的比热c、密度ρ、导温系数α、导热系数λ、绝热温升θ0以及常数m。上述参数中混凝土的比热c、密度ρ可直接测得，计算精度可满足计算要求。绝热温升θ0受到水泥品种、用量、骨料粒径以及实验室环境与施工环境差异等多种因素的影响，如果不通过试验确定就应进行反分析;导热系数λ受混凝土密实性、材料特性以及骨料岩性的影响，难以确定，应予以反演。导热系数λ=αcρ，因此，只需反演λ、α任意一个参数即可;水泥水化热计算表达式主要有指数式、双曲线式、复合指数式3种[1]。计算所用的水化热表达式采用指数式函数θ（t）=θ0（1-e-mt），其常数m随水泥品种、比表面及浇筑温度不同而不同，更能体现施工现场实际浇筑环境，也确保所选表达式与有限元分析软件热源函数的定义保持一致。除绝热温升θ0外，还应对常数m进行反演。因此，选取绝热温升θ0、常数m、导热系数λ作为反演参数，根据相关文献，确定取值范围见表1。

3 算例分析

3.1 工程概况

太洪长江大桥是主跨为808 m的钢箱梁悬索桥，桥梁全长为1 436 m，为重庆南川至两江新区高速公路NL5标段控制性工程。大桥南岸锚碇为隧道式锚碇，隧道锚混凝土施工包括后锚室、锚塞体、前锚室以及散索鞍支墩承台，隧道锚立面布置见图4。散索鞍支墩承台分左右两幅，每幅尺寸为15 m（长） ×13 m（宽）×6 m（高），混凝土浇筑总量约2 340 m3，是典型的大体积混凝土结构。支墩承臺分3层浇筑，每层浇筑厚度均为2 m，见图5。采用C40混凝土，混凝土设计配合比的材料用量如表2所示。

3.2 温度测点布置及监测方法

考虑到结构的对称性和温度的变化规律，选取1/4结构进行温度测点布置，测点布置如图6所示[18]。采用实时无线温度综合测试系统进行监测，可根据需要设置温度采集频率，温度监测频率为2 h采集一次。太洪长江大桥南岸散索鞍支墩承台于2019年3月18日21：00开始浇筑，3月19日19：00浇筑完成，浇筑历时22 h，测得混凝土入模温度为16.9 ℃。为获得真实的桥址环境，环境温度从混凝土开盘前一星期开始监测，在施工现场的背阴处布置测点，把实测环境温度输入有限元模型，排除昼夜温差干扰。温度曲线如图7所示。

3.3 有限元模型

采用通用有限元软件 Midas FEA 建立有限元模型，共有87 035个单元，16 848个节点。如图8所示。模型考虑基岩对混凝土水化热的吸收作用，建立地基扩大模型，固定温度取为20 ℃，基岩侧面和底面给予固定约束。根据现场实际的施工过程，分为3个施工阶段，即2 m一层。混凝土抗压强度通过试验确定，根据对试拌混凝土的各项技术指标检测及抗压强度试验，得到7、28 d抗压强度分别为44.2、49.7 MPa。

3.4 BP神经网络的建立

3.4.1 隐含层的确定

神经网络隐含层神经元个数的选择对网络的性能影响很大，所以，对层内神经元数目需要进行恰当的选择。隐含层神经元的选择十分复杂，目前并没有统一的理论指导。隐含层神经元数目与待解决问题的特征、输入与输出单元数直接相关。通常情况下，隐含层神经元太少，网络可能训练不出来或网络性能较差;若隐含层神经元太多，可能使训练时间延长，训练容易陷入局部极小值点，误差也不一定达到最佳。一般情况下，仅能靠经验和多次试验来确定隐含层单元数。根据经验，可以参照公式l=n+m+a进行设计。式中：l为隐含层神经元数目，n为输入层单元数，m为输出层节点数;a为1～10之间的调节常数。在经验公式的基础上，采用逐步增加隐含层单元数的变结构法，即开始放入比较少的隐含层单元数，学习一定次数后，性能不佳就增加隐含层单元数，一直达到比较理想的隐含层单元数为止。根据训练数据，确定隐含层单元数为15，网络允许最大误差设为0.001，学习速率为0.05，进行网络训练。

3.4.2 样本设计及数据归一化处理

采用均匀设计方法确定参数样本，参数的水平数取为31，选取均匀设计表为U31（3110），样本值见表3。将表3的参数样本数据输入有限元计算模型，得到支墩承台特征点的温度计算值，见表4。

考虑到各参数之间的量纲影响以及小数值信息被大数值信息淹没现象的发生，在处理输入与输出数据时，要用到归一化方法。神经网络模型采用tan-sigmod型传递函数，该函数的值域为[-1， 1]，因此，在计算过程中用式（1）预处理数据，通过实际样本试算，使其全部归一化到[-1， 1]区间内，归一化公式为

归一化用MATLAB语言实现。

[input_train，inputps]=mapminmax（P1）;

[output_train，outputps]=mapminmax（T1）;

其中：input_train和output_train为归一化返回的值，结构体inputps和outputps是进行归一化时所用的参数。

3.4.3 网络训练及效果评估

为了提高网络训练速度，使学习时间更短，同时保证训练过程的稳定性，采用附加动量法对神经网络进行优化，附加动量法是在每个加权调节量上加上一项正比例于前次加权变化量的值，使网络在修正其权值时不仅考虑误差在梯度上的作用，而且考虑在误差曲面上变化趋势的影响，在没有附加动量的作用下，网络可能陷入浅的局部极小值，利用附加动量可以带动梯度下降过程，冲过狭窄的局部极小值，从而提高训练速度。具体表达为

在网络训练前还需对一部分参数期望值先进行设置。网络最大训练次数设置为1 000次，训练期望精度设置为1×10-4。网络误差梯度目标值设为默认值0，网络训练过程中，误差梯度实际值约为0.004 76，与目标值0接近，认为误差梯度符合要求。神经网络训练过程中，为了防止网络的过度训练而使泛化能力降低，设置有终止训练功能的有效性检查步数，即确认样本的误差曲线不再下降的连续迭代次数，在网络训练之前，对有效性检查步数的值进行不同设定后分别进行试验，最终确定为6。训练过程误差曲线见图9，可以看出，优化后的训练过程随着训练次数的递增稳定收敛，比优化前更快、更早收敛。BP神经网络训练过程结果具有随机性，为评估网络训练效果是否达到预期，用12组归一化

后的数据多次反演来检验测试样本训练效果，以平均绝对百分比误差（MAPE）和均方根误差（RMSE）作为误差评价指标。

式中：observes为预测值;predicted为设计值;亦称真实值;N是样本数。平均绝对百分比误差又可称为相对误差的绝对值平均值，由于离差被绝对值化，不会出现正负相抵消的情况，因此，更能反映训练效果的可信程度。与相对误差类似，它是一个百分比值，即如果MAPE值为5，则表示预测值较真实值平均偏离5%。均方根误差是用来衡量观测值同真实值之间的偏差，RMSE值越小，表示精度越高。每次训练的平均绝对百分比误差、均方根误差见表5，第

1次训练的热学参数设计值与预测值的趋势见图10，可以看出：平均绝对百分比误差均小于10%，均方根误差值均小于1.2，预测值与设计值拟合度较高，说明神经网络模型对大体积混凝土预测精度较高，可以用来反演大体积混凝土热学参数。

3.5 反演分析结果与检验

将特征点实测温度值输入网络，输出即为3个热学参数的反演值，见表6。反分析完成后，还需对热学参数反演值进行检验，具体方法是，将反演值输入有限元模型，计算出第2个施工阶段特征点处溫度计算值，并与该施工阶段测得的实际温度值进行比较，检验二者的拟合程度，计算温度值与实测温度值拟合曲线如图11～图14所示。由图可知，温度计算值与实测值之间的误差较小，变化规律基本一致，表明基于BP神经网络反演得到的热学参数符合混凝土的实际施工环境，该组反演值真实可靠。

4 结论

以太洪长江大桥散索鞍支墩承台为工程背景，为得到混凝土真实的热学参数，基于均匀设计和 BP神经网络对大体积混凝土绝热温升、反应速率、导热系数等热学参数进行反演，揭示了温度场与热学参数之间的内在联系和规律，得到如下结论：

1）施工期大体积混凝土的温度峰值与混凝土热学参数之间的复杂非线性关系可由BP神经网络表述完成，运用BP神经网络可避免主观调整热学参数所造成的误差，有较强的实用性。

2）在BP神经网络训练阶段，采用附加动量法优化网络结构能明显减少网络训练时间，提高训练效率。

3）利用均匀设计和BP神经网络相结合的反分析方法，可以大大减少网络学习的样本数量。同时，可使有限元正分析与反分析过程分离，大大减少了反分析时间，提高了反分析效率和准确性。

4）通过BP神经网络反分析得到的大体积混凝土热学参数分别为：绝热温升θ0=46.834 2 ℃，常数m=0.979 504，导热系数λ=9.742 07 kJ/（m·h·℃），结合有限元正分析对后续施工的混凝土温度场进行预测，得到的特征点温度计算值与温度实测值较为接近，在变化规律上基本吻合，温度峰值最大误差仅为1.1 ℃，表明基于均匀设计及BP神经网络方法可较准确地反演大体积混凝土热学参数。

参考文献：

[1] 朱伯芳. 大体积混凝土温度应力与温度控制[M]. 北京： 中国水利水电出版社， 2012： 738.

ZHU B F. Thermal stresses and temperature control of mass concrete [M]. Beijing： China Water Power Press， 2012： 738.

[2] HATTEL J H， THORBORG J. A numerical model for predicting the thermomechanical conditions during hydration of early-age concrete [J]. Applied Mathematical Modelling， 2003， 27（1）： 1-26.

[3] 朱岳明， 刘勇军， 谢先坤. 确定混凝土温度特性多参数的试验与反演分析[J]. 岩土工程学报， 2002， 24（2）： 175-177.

ZHU Y M， LIU Y J， XIE X K. Determination of thermal parameters of concrete by reverse analysis of test results [J]. Chinese Journal of Geotechnical Engineering， 2002， 24（2）： 175-177. （in Chinese）

[4] 王振红， 朱岳明， 武圈怀， 等. 混凝土热学参数试验与反分析研究[J]. 岩土力学， 2009， 30（6）： 1821-1825， 1830.

WANG Z H， ZHU Y M， WU Q H， et al. Thermal parameters of concrete by test and back analysis [J]. Rock and Soil Mechanics， 2009， 30（6）： 1821-1825， 1830. （in Chinese）

[5] AMIRFAKHRIAN M， ARGHAND M， KANSA E J. A new approximate method for an inverse time-dependent heat source problem using fundamental solutions and RBFs [J]. Engineering Analysis with Boundary Elements， 2016， 64： 278-289.

[6] 李守巨， 刘迎曦. 基于模糊理论的混凝土热力学参数识别方法[J]. 岩土力学， 2004， 25（4）： 570-573.

LI S J， LIU Y X. Identification of concrete thermal parameters based on fuzzy theory [J]. Rock and Soil Mechanics， 2004， 25（4）： 570-573. （in Chinese）

[7] 孙伟， 何蕴龙. 基于改进遗传算法的垫层区混凝土热学参数反演和早期裂缝成因分析[J]. 武汉大学学报（工学版）， 2014， 47（2）： 193-200.

SUN W， HE Y L. Inversion analysis of thermal parameters of cushion concrete and reasons of early stage cracks based on improved genetic algorithms [J]. Engineering Journal of Wuhan University， 2014， 47（2）： 193-200. （in Chinese）

[8] 倪智强， 周兰庭. 基于改进蚁群算法的混凝土坝热学参数反演[J]. 水电能源科学， 2018， 36（4）： 82-85.

NI Z Q， ZHOU L T. Inverse analysis of concrete dam thermal parameters based on improved ant colonyalgorithm [J]. Water Resources and Power， 2018， 36（4）： 82-85. （in Chinese）

[9] 王峰， 周宜红， 赵春菊， 等. 基于改进粒子群算法的混凝土坝热学参数反演研究[J]. 振动与冲击， 2019， 38（12）： 168-174， 181.

WANG F， ZHOU Y H， ZHAO C J， et al. Inverse analysis of concrete dam thermal parameters based on an improved particle swarm optimization method [J]. Journal of Vibration and Shock， 2019， 38（12）： 168-174， 181. （in Chinese）

[10] 喻正富， 张筱雨， 刘来君， 等. 基于遗传算法的大体积混凝土热力学参数反演分析[J]. 建筑科学与工程学报， 2015， 32（5）： 81-88.

YU Z F， ZHANG X Y， LIU L J， et al. Inverse analysis of mass concrete thermal parameter based on genetic algorithm [J]. Journal of Architecture and Civil Engineering， 2015， 32（5）： 81-88. （in Chinese）

[11] 文豪， 王高峰， 侯章偉. 基于MATLAB和ANSYS的大体积混凝土水化热效应分析[J]. 内蒙古公路与运输， 2018（3）： 5-10.

WEN H， WANG G F， HOU Z W. Hydration heat effect analysis of massive concrete based on MATLAB and ANSYS [J]. Highways & Transportation in Inner Mongolia， 2018（3）： 5-10. （in Chinese）

[12] 方开泰. 实用回归分析[M]. 北京： 科学出版社， 1988.

FANG K T. Practical regression analysis [M]. Beijing： Science Press， 1988.（in Chinese）

[13] 周品. MATLAB神经网络设计与应用[M]. 北京： 清华大学出版社， 2013.

ZHOU P. Design and application of MATLAB neural network [M]. Beijing： Tsinghua University Press， 2013.（in Chinese）

[14] 张清， 宋家蓉. 利用神经元网络预测岩石或岩石工程的力学性态[J]. 岩石力学与工程学报， 1992， 11（1）： 35-43.

ZHANG Q， SONG J R. Predicting mechanical behaviors of rock or rock engineering by neural network [J]. Chinese Journal of Rock Mechanics and Engineering， 1992， 11（1）： 35-43. （in Chinese）

[15] 李端有， 李迪， 马水山. 三峡永久船闸开挖边坡岩体力学参数反分析[J]. 长江科学院院报， 1998， 15（2）： 10-13.

LI D Y， LI D， MA S S. Backward analysis of the rock mass mechanical parameters in excavated slope of TGP permanent lock [J]. Journal of Yangtze River Scientific Research Institute， 1998， 15（2）： 10-13. （in Chinese）

[16] 宫经伟. 水工准大体积混凝土分布式光纤温度监测与智能反馈研究[D]. 武汉： 武汉大学， 2013.

GONG J W. Study on temperature monitoring by distributed optical fible and the intelligent feedback for hydraulic quasi-mass concrete structure [D]. Wuhan： Wuhan University， 2013.（in Chinese）

[17] CHEN G R， YANG B. Back analysis of thermal parameters of mass concrete and its application based on micro-genetic algorithm [J]. Advanced Materials Research， 2014， 936： 1463-1467.

[18] 张玉平， 李传习， 张磊， 等. 张花高速澧水大桥索塔基础大体积混凝土温控[J]. 科技导报， 2012， 30（34）： 55-60.

ZHANG Y P， LI C X， ZHANG L， et al. Mass concrete temperature control of Lishui bridge tower foundation in Zhanghua highway [J]. Science & Technology Review， 2012， 30（34）： 55-60. （in Chinese）

（编辑 王秀玲）