1、基于成桥荷载试验的襄阳汉江五桥有限元模型修正导师：*答辩人：*学 院：理 学 院专 业：工程力学 汉江五桥是位于襄阳市的一座典型的大跨度连续刚构梁拱组合桥构造特点：拱肋：拱肋：1.1.变变高度八边形箱型高度八边形箱型截面，截面，内部设加劲肋、横隔板。内部设加劲肋、横隔板。2.2.净宽：净宽：2.47m,2.47m,高度：高度：2.473.97m2.473.97m。3.3.共共8484根吊杆，间隔根吊杆，间隔6m6m。跨径：跨径：77+138+138+77=430m77+138+138+77=430m 主梁主梁:1.:1.大大悬臂斜腹板单箱三室悬臂斜腹板单箱三室截面截面，桥面横向坡度桥面横向坡度
2、2 2%。2.2.桥墩处箱梁高度达桥墩处箱梁高度达7.8m7.8m，跨中跨中3.0m3.0m。3.3.桥墩与箱梁之间采用固结形式。桥墩与箱梁之间采用固结形式。汉江五桥是内环线闭合环的主要过江通道，作为市区内环线控制性工程，对其全桥有限元模型进行“整体-局部”修正，对于长期健康监测具有重要价值。背景介绍模型修正局部分析刚域探讨根据工程实际所建立的MIDAS全桥模型，对于整体有限元模型修正中几个关键因素进行分析研究在局部分析的基础上利用成桥荷载试验结果和人工神经网络方法，重点对拱桥拱脚部位的刚域效应及结构的杨氏模量合理取值进行探讨由“整体到局部”的分析思路，利用ANSYS软件掌握拱脚处应力大小及分
3、布规律，从而对拱脚处构造做出综合评价处理方式的不同处理方式的不同：梁单元模型：拱肋隔板以节点荷载的形式施加于拱肋节点处；板单元模型：按设计对应截面处增设加强隔板、标准隔板。梁与板单元模型边墩对比梁与板单元模型边墩对比梁与板单元模型中墩对比梁与板单元模型中墩对比挠度对比挠度对比箱梁剪力对比箱梁剪力对比箱梁弯矩对比箱梁弯矩对比非常接近，一部分节点存在0.010.02mm差别，跨中0.41mm模型二最大竖向位移83.9mm，相差1%拱脚附近：10%15%拱肋拱肋箱梁箱梁箱梁应力除墩顶外，其余位置均比较接近拱肋应力：模型一 模型二 由于模型一中没有考虑隔板参加受力，使得顶底板分担的应力比模型二稍有增加
4、，拱脚部位相差稍大。拱脚处于拱肋与箱梁交接处，是结构传力的重要部位，因此底板存在应力集中现象，梁单元模型则无法体现这一特点。由于无加强横隔板，梁单元模型在吊杆连接处出现了明显的应力集中现象，在板单元模型中由于横隔板的存在较好地消除了这一现象。在加劲肋与横隔板连接处出现较大拉应力，这些部位在梁单元模型中都是被简化忽略掉的。有无桩土作用模型对比“m”法：zxzm z x()()pzszxspzzza bm z xPAKa bm zxxx工况一：恒载(1.0)+活载(1.0)+风荷载(1.2)工况二：恒载(1.0)+活载(0.85)+风荷载(0.75)+整体升温(1.2)工况三：恒载(1.0)+活载
5、(0.85)+风荷载(0.75)+整体降温(1.2)试验工况：工况（不考虑）轴力Fx 剪力Fy剪力Fz 扭矩Mx 弯矩My弯矩Mz 应力(kN)(kN)(kN)(kN*m)(kN*m)(kN*m)(MPa)115671.72-354.0855234.9021083.56-1501738.545739.09-8.96216541.29-68.9236803.2416435.04-1458615.956883.38-8.68312569.27-68.9266989.6716435.04-1550818.956883.38-9.28工况（考虑）轴力Fx 剪力Fy剪力Fz 扭矩Mx 弯矩My弯矩Mz
6、应力(kN)(kN)(kN)(kN*m)(kN*m)(kN*m)(MPa)115679.54-206.5951717.0119257.90-1497744.466998.48-8.98216625.92-88.0436789.5616106.39-1454937.407023.18-8.70312465.21-88.0461503.9716106.39-1546669.327023.18-9.30箱梁拱脚内力对比表：1.考虑桩土作用模型中的土弹簧分担了很大一部分横向荷载。2.基底桩为摩擦桩，因此剪力相差较大，达到6.8%。3.在考虑桩土作用时未考虑桩身转动弹性刚度。工况（不考虑）轴力Fx 剪力
7、Fy剪力Fz 扭矩Mx 弯矩My弯矩Mz 应力(kN)(kN)(kN)(kN*m)(kN*m)(kN*m)(MPa)113108.4961.202040.289831.82281650.46-1500.7515.1211573.85101.401188.158461.53258090.162769.2413.8311296.96101.402344.418461.53298843.732769.2515.9工况（考虑）轴力Fx 剪力Fy剪力Fz 扭矩Mx 弯矩My弯矩Mz 应力(kN)(kN)(kN)(kN*m)(kN*m)(kN*m)(MPa)112970.70149.962144.6299
8、03.30282111.60-1900.2415.0211664.15139.701229.518489.45258313.24432.9213.7310998.63139.702510.828489.45299425.64432.9215.8箱梁跨中内力对比表：1.与箱梁拱脚类似，横向剪力与水平弯矩增幅对应拱脚部位的降幅。在全桥结构要保持平衡的前提条件下，拱脚部位因为桩土效应而减小的内力值转移到了跨中。2.结构内部由于膨胀和收缩产生变形，内力累积效应在跨中部位体现最为明显。工况（不考虑）轴力Fx 剪力Fy剪力Fz 扭矩Mx 弯矩My弯矩Mz 应力(kN)(kN)(kN)(kN*m)(kN*m
9、)(kN*m)(MPa)1-17369.45463.091347.07600.299892.4013826.26-1392-22614.51289.651519.32376.9811111.818649.09-1493-11686.76289.651036.05376.986914.018649.09-111工况（考虑）轴力Fx 剪力Fy剪力Fz 扭矩Mx 弯矩My弯矩Mz 应力(kN)(kN)(kN)(kN*m)(kN*m)(kN*m)(MPa)1-17199.76698.641307.94902.679442.5420852.43-1492-22641.83289.601510.62376
10、.8111007.938645.71-1493-11488.29289.601009.27376.816717.718645.71-110拱肋拱脚内力对比表：1.工况1增幅近50%，基底的水平方向的弹簧刚度增加，对上部的拱肋结构的抵抗荷载能力也有了更高的要求。2.摩擦桩的设置使得桩底弹性系数增加，竖向弯矩减小5%，变化幅度较小。3.侧重考察桩土作用对全桥整体内力影响方面，观察发现数据变化幅度很小。未知系数法未知系数法：手工迭代法手工迭代法：AXD A 12,inAAAA111212122212mmnnnmAAAAAAAAA=未知系数法是MIDAS软件的调索功能，理论基础来自于影响矩阵法：=得到
11、响应影响矩阵，通过影响矩阵计算，得到每个阶段的施调向量 ，满足结构平衡状态条件下使吊杆力逐渐趋向目标值。A1=()nnniniinidesresSSSS 式中，Sini为索力输入值，Sres为索力输出值，Sdes为成桥索力设计值，为下次迭代输入值。1niniS初始张拉力：750KN设置上下限：751KN，749KN不同点：1.初始值 2.简易程度拱肋与箱梁弯矩对比图拱肋与箱梁弯矩对比图1.相对相对于手工迭代法，未知荷载系数法的拱肋计算弯矩平均减小于手工迭代法，未知荷载系数法的拱肋计算弯矩平均减小4%；2.与与之相比箱梁弯矩在左右拱脚处变化较大，相差在之相比箱梁弯矩在左右拱脚处变化较大，相差在1
12、7%左右左右。3.1/4跨跨 1/2跨与跨与1/23/4跨之间吻合，差值跨之间吻合，差值5%以下。以下。两种方法在不同区段有不同的特征两种方法在不同区段有不同的特征拱肋与箱梁剪力对比图拱肋与箱梁剪力对比图1.两种计算方法所得拱肋剪力与箱梁剪力均比较吻合两种计算方法所得拱肋剪力与箱梁剪力均比较吻合；2.未知荷载系数法的结果要小于手工迭代法，差值在未知荷载系数法的结果要小于手工迭代法，差值在3%左右左右。拱肋与箱梁挠度对比图拱肋与箱梁挠度对比图1.挠度差值均出现两头小、中间大的分布态势挠度差值均出现两头小、中间大的分布态势；2.拱肋对箱梁的位移约束能力不如手工迭代法拱肋对箱梁的位移约束能力不如手工
13、迭代法。1.横隔板的存在对抑制结构的变形很有必要横隔板的存在对抑制结构的变形很有必要2.板单元应力内力计算结果更理想3.两种模型总体差别不大1.未知系数法未知系数法在弯矩、轴力、正应力三个方面的参数表现在弯矩、轴力、正应力三个方面的参数表现要优于手工迭代法要优于手工迭代法2.未知系数法的成桥索力分布很不均匀3.手工迭代法对吊杆力的精确控制更加高效，但费时间1.桩土效应对结构的影响主要表现在水平方向桩土效应对结构的影响主要表现在水平方向2.桩土作用对梁拱组合桥型的静力荷载作用影响较小3.有限元模型整体静力分析可不考虑桩土作用梁、板单元对比梁、板单元对比吊杆力优化吊杆力优化桩土效应桩土效应ANSY
14、S拱脚局部模型 全桥模型无法对结构局部的受力性能做出准确判断，因此有必要采用实体单元有限元法对其进行局部分析，掌握拱脚处的应力大小及分布规律。疏密网格过渡区疏密网格过渡区金字塔单元金字塔单元壳壳单元嵌入单元嵌入箱梁：Solid95拱肋：Shell181结构分析采用1/4对称边界条件：纵向对称面：X向平动自由度 Y、Z向转动自由度横向对称面：Z向平动自由度 X、Y向转动自由度 利用平截面假定，主节点设置于界面的几何中心，主节点处添加相对刚度无限大的虚拟梁单元，用来施加节点荷载。结构部位最大竖向位移(cm)结果分析MIDAS/CivilANSYS箱梁-0.49-0.44-0.66混凝土结构位移绝对
15、数值较小，结果相近拱肋-2.53-1.58-1.73钢结构位移相差较大，原因是拱脚处建模简单，刚度偏小 箱梁整体第一主应力处于箱梁整体第一主应力处于-3.791.92MPa之间，整体之间，整体第三第三应力处于应力处于-15.8-5.39MPa之之间，两者都存在明显的应力集中现象，应力集中点位于预应力束与混凝土实体单元耦合间，两者都存在明显的应力集中现象，应力集中点位于预应力束与混凝土实体单元耦合处。处。拱拱脚脚箱箱梁梁应应力力分分布布图图箱梁整体应力分布图拱拱脚脚箱箱梁梁应应力力分分布布图图14号块整体应力分布图 1.1.一号块第一主应力处于一号块第一主应力处于-2.312.12MPa之间，之
16、间，第三第三应力处于应力处于-11.92.32MPa之间，之间，边边室室底板底板、箱梁箱梁与横隔板与横隔板交界处吊杆交界处吊杆锚固处第一主应力值较高，局部超出混凝土设计锚固处第一主应力值较高，局部超出混凝土设计抗抗拉强度拉强度。2.2.二号块第一主应力处于二号块第一主应力处于-2.242.38MPa之间，之间，第三第三应力处于应力处于-10.060.201MPa之间之间，箱梁与横隔板交界处第一主应力值较高，局部超出混凝土设计箱梁与横隔板交界处第一主应力值较高，局部超出混凝土设计抗拉强度抗拉强度，边室隔板、边室隔板、中室隔板与腹板交界处存在应力集中中室隔板与腹板交界处存在应力集中现象现象，接近设
17、计强度。，接近设计强度。3.3.三号号块第一主应力处于块第一主应力处于-2.423.06MPa之间，之间，第三第三应力处于应力处于-13.40.234MPa之间之间，三号块中室底板、箱梁与横隔板交界处三号块中室底板、箱梁与横隔板交界处、吊杆、吊杆锚固处第一主应力值较高，局部超出混凝锚固处第一主应力值较高，局部超出混凝土设计土设计抗拉强度抗拉强度。4.4.四四号号块第一主应力处于块第一主应力处于-2.353.02MPa之间，之间，第三第三应力处于应力处于-15.8-5.39MPa之间之间，四四号块中室底板第一主应力值较高，局部超出混凝土设计抗拉强度，需引起重视。此外，在号块中室底板第一主应力值较
18、高，局部超出混凝土设计抗拉强度，需引起重视。此外，在箱梁与横隔板交界处拉应力值也较大。箱梁与横隔板交界处拉应力值也较大。锚固区锚固区第一主应力处于第一主应力处于-0.9422.25MPa之间，整体第三主应力处于之间，整体第三主应力处于-16.70.11MPa之间之间。钢板嵌入区域外侧第一主应力偏高，局部接近混凝土设计抗拉强度，钢拱肋底部第三钢板嵌入区域外侧第一主应力偏高，局部接近混凝土设计抗拉强度，钢拱肋底部第三主应力值偏高，需引起重视。主应力值偏高，需引起重视。拱拱脚脚箱箱梁梁应应力力分分布布图图锚固区整体应力分布图 拱肋整体应力分布均匀，拱肋整体应力分布均匀，整体整体Von MisesVo
19、n Mises应力处于应力处于0.055144MPa0.055144MPa之间，局部最大当量之间，局部最大当量应力达到应力达到258MPa258MPa。在在GL1bGL1b上面板、竖向横隔板以及拱肋底部，与混凝土相衔接处存在应上面板、竖向横隔板以及拱肋底部，与混凝土相衔接处存在应力集中，应力梯度较大力集中，应力梯度较大。拱拱脚脚拱拱肋肋应应力力分分布布图图拱肋整体应力分布图 箱梁整体箱梁整体Von Mises应力处于应力处于0.055115MPa之间，局部最大当量应力达到之间，局部最大当量应力达到258MPa，拱脚拱脚底部底部结构结构受力复杂，结构处于高应力状态，在拱肋底面板处可受力复杂，结构
20、处于高应力状态，在拱肋底面板处可以明显看到纵肋的以明显看到纵肋的变形变形，同时根部的横隔板处于高应力状态，同时根部的横隔板处于高应力状态。拱拱脚脚拱拱肋肋应应力力分分布布图图拱肋根部应力分布图1.1.拱拱脚在三向预应力束体系作用下，局部混凝土压应力大大高脚在三向预应力束体系作用下，局部混凝土压应力大大高于悬臂梁部分于悬臂梁部分。由于上部钢拱肋的倒角，结构构造以及外力分布由于上部钢拱肋的倒角，结构构造以及外力分布复杂，而导致局部混凝土结构存在明显应力集中处，使结构处于复杂，而导致局部混凝土结构存在明显应力集中处，使结构处于高高应力状态应力状态。2.2.与与拱脚连接处的拱肋处于高应力状态，在拱肋底
21、面板处纵肋拱脚连接处的拱肋处于高应力状态，在拱肋底面板处纵肋变形较大。除部分应力集中点外，拱脚处拱肋应力均处于合理变形较大。除部分应力集中点外，拱脚处拱肋应力均处于合理范范围围。3.3.传统以杆系理论为基础的数值分析方法难以准确反映拱脚处的传统以杆系理论为基础的数值分析方法难以准确反映拱脚处的结构应力状态，在用于荷载试验时容易造成误判，此时应采用比结构应力状态，在用于荷载试验时容易造成误判，此时应采用比较细致的有限元局部分析较细致的有限元局部分析。4.拱脚处箱梁与拱肋之间为固结形式，由于采用混合结构，两者拱脚处箱梁与拱肋之间为固结形式，由于采用混合结构，两者刚度均比较大，因此应考虑其产生刚域效
22、应，对此将在下节讨论。刚度均比较大，因此应考虑其产生刚域效应，对此将在下节讨论。231456789101112134*19.25=774*19.25=774*34.5=1384*34.5=138L1L2L3L4测点分布：按照桥面线形测量的惯例，在左航道桥上每隔L/4设置一个测点，测点设在左右防撞墙内侧桥面上。在不同试验工况下，对大桥的关键部位如各跨的L/4处、跨中、3L/4处及墩顶截面处的试验前后挠度值进行测量，得到试验荷载作用下结构的挠度值。本次试验以汉江五桥主桥左航道桥为试验对象。利用检测仪器测试结构的控制部位和控制截面在荷载作用下的挠度等特性的变化，从而对结构在试验荷载作用下的刚度、承载
23、能力进行评价。2#3#4#36,59,677138L1L22#3#4#36,59,69,635,877138L1L22#3#4#77138L1L264,62#3#4#36,5L1L2L3698,917,2771381381.试验工况1：检验边跨最大正弯矩截面A承受正弯矩的荷载。2.试验工况2：检验边跨墩顶截面B承受负弯矩的承受能力。3.试验工况3：检验中跨跨中截面C承受正弯矩的能力。4.试验工况4：检验拱脚截面D承受（对称）负弯矩的能力。试验选用350KN双后轴载重车进行最不利加载，试验车辆前轴重70KN，中后轴重均为140KN，前轴距中轴4m，中轴距后轴1.4m，荷载效率系数（加上冲击荷载系
24、数影响）为1.05。测点测点编号编号工况三工况三工况四工况四实测值理论值 校验系数 实测值 理论值 校验系数15.496.060.912.593.870.67 23.735.390.691.582.620.60 31.282.210.580.210.430.49 4-1.33-1.860.724.585.370.85 5-1.32-1.420.938.809.200.96 6-0.12-0.250.481.352.640.51 70.030.050.60-5.60-7.900.71 80.030.050.60-0.25-0.620.40 90.020.030.67-0.05-0.130.381
25、11.250.04-1.50-0.39 -12-1.00-1.400.713.003.800.79 13-0.65-0.820.796.006.700.90 静载试验主要测点挠度值静载试验主要测点挠度值表表：测点测点编号编号工况三工况三工况四工况四实测值理论值 校验系数 实测值 理论值 校验系数1-1.09-1.230.895.116.100.842-1.07-1.560.695.965.451.093-1.04-1.001.041.162.250.5241.104.620.24-2.59-2.151.20511.212.200.92-2.70-2.970.9162.904.670.621.3
26、02.250.587-1.00-1.500.673.485.560.638-1.00-1.180.855.487.730.719-0.25-0.241.041.132.060.55112.002.720.74-1.00-1.250.80126.008.650.69-1.00-1.700.59131.002.750.36-1.00-1.330.75静载试验主要测点挠度值静载试验主要测点挠度值表表：结论：结论：多数测点在静载试验工况下表现正常，除个别测点以外，挠度值校验系数均未超过1.0。但也反映出另一方面的问题，理论值相对实测值偏大，导致挠度校验系数普遍偏小，说明有限元建模过程中对材料参数的选取
27、，以及结构刚度、质量等参数与工程实际结构不相吻合，在此条件下得到的结构校验系数不能准确反映大桥的承载能力及安全储备。修正因素：刚域弹性模量修正措施：在箱梁与拱肋端部设置一定区段的刚域，同时对混凝土、拱肋、预应力束等弹性模量上下调整10%进行合理性探讨。拱脚部位因为拱脚部位因为构造构造原因原因会形成一段刚度趋近于无限大会形成一段刚度趋近于无限大的区域，即为刚域，在此部位的区域，即为刚域，在此部位区段的材料刚度系数如果不相区段的材料刚度系数如果不相应作以调整，将使模型模拟实应作以调整，将使模型模拟实际结构出现际结构出现失真失真。弹性模量精确弹性模量精确与否直接影响到分与否直接影响到分析结果析结果精
28、度。精度。在工程实际中，在工程实际中，由于由于混混凝土凝土徐徐变、变、收缩收缩，钢拱肋和预应力束，钢拱肋和预应力束老化老化等等原因原因，结构结构的的弹性模量相比设弹性模量相比设计规范值会有一定计规范值会有一定的的偏差和出入。偏差和出入。杆系计算内力真实内力(试验及细分有限元计算内力)(1/6-1/7)bhb 刚刚域是指构件与构件相交部分弯曲刚度按照无限大考虑的区域，域是指构件与构件相交部分弯曲刚度按照无限大考虑的区域，刚域一般产生于框架中梁柱重叠部分，其轴向刚度与抗弯刚度都很刚域一般产生于框架中梁柱重叠部分，其轴向刚度与抗弯刚度都很大，基本上可以认为不发生轴向、弯曲变形，从而在整体上看作一大，
29、基本上可以认为不发生轴向、弯曲变形，从而在整体上看作一个刚性个刚性块体块体。在在有限元有限元软件中模拟节点刚域的方式主要有两种软件中模拟节点刚域的方式主要有两种:1.主从主从节点自由度节点自由度法法 2.刚性刚性材料法材料法。优缺点：第一种第一种：需要将一个主节点与多个从属节点耦合在一起，耦合节点需要将一个主节点与多个从属节点耦合在一起，耦合节点之间的单元刚度趋于无穷大，节点区域实现了刚域效应，缺点是需之间的单元刚度趋于无穷大，节点区域实现了刚域效应，缺点是需要重新划分节点，工作量稍大。要重新划分节点，工作量稍大。第二种第二种：采用刚度很大的单元模拟节点附近的刚性区域。采用刚度很大的单元模拟节
30、点附近的刚性区域。a1a2c1c2lc1lc2lb1lb2bchb112211220.250.250.250.25bbbbcccclhlhlcblcb 日本奥田勇教授提出下列刚域尺寸计算公式：=-0.25=-0.25zfzfLBhLhB梁柱高规高规对于梁柱交汇点处的刚性节点区域的具体刚域尺寸进行了对于梁柱交汇点处的刚性节点区域的具体刚域尺寸进行了规定规定：其中，hf是梁高，Bz是柱宽。由于梁柱斜交，交汇处箱梁高度为7.9m,拱肋高度为4.0m，得到：L梁刚域=0.025m L拱刚域=2.95m，梁端刚域计算尺寸过小。采用奥田勇教授的理论公式估算，得到箱梁刚域长度：L梁刚域=2.205m L拱刚
31、域=2.95m刚刚域对测点域对测点挠度的挠度的影响影响弹模弹模对对测点测点挠度的挠度的影响影响绝对差值：0.021.98mm相对差值：7%差值：混凝土：0.520.91mm拱肋：0.120.36mm预应力束：0.110.14mmBP网络模型BP神经网络神经网络算法步骤算法步骤如下：如下：（1）对阈值进行初始赋值。（2）提供输入样本集与期望输出值。（3）计算输出值：()jji iyfW x1()1 exp()f xx式中，f(*)为Sigmoid函数，即：（4）网络调整。从输出层开始反馈，将误差反向传播对权值进行调整：(1)()ijijjjW tW ty 式中，j按下式计算：(1)()(1)()
32、iiijiikjkkyyTyjyyWj当 为输出节点当 为隐节点(1)()()(1)ijijjjijijWtWtyWtWt使用冲量值时，权值按下式计算：通过这种方式，不仅确定了最佳试验组合，并且建立起各相关因素之间的对应模型，利用均匀设计方法的的特点，使得试验样本具有分布均匀，代表性强的特点。用BP神经网络法挑选出表现最好的试验组合，结束模型修正试验。建立起描述指标与对应因素的关系模型，在一个较大的试验范围内筛选出最优试验条件组合。缩小范围，以对建立的模型在新一轮试验中进行验证，从而达到对模型进一步修正的目的。测点测点编号编号刚域区段刚域区段0.10.20.30.40.50.60.70.80.
33、91.01-1.21-1.19-1.18-1.17-1.15-1.13-1.09-1.03-0.95-0.912-1.52-1.49-1.46-1.43-1.41-1.39-1.37-1.35-1.34-1.333-1.01-1.02-1.03-1.04-1.06-1.08-1.11-1.14-1.18-1.2344.494.384.294.224.154.104.064.034.013.99512.1111.9811.8311.6811.5111.3311.1310.9110.6610.3864.604.544.494.454.414.374.344.314.294.277-1.46-1.4
34、3-1.40-1.37-1.35-1.32-1.29-1.26-1.24-1.228-1.18-1.19-1.20-1.21-1.23-1.25-1.27-1.30-1.33-1.379-0.24-0.24-0.23-0.22-0.22-0.21-0.21-0.20-0.20-0.20112.732.752.772.792.812.832.862.902.943.00128.478.318.178.007.867.757.677.617.567.54132.652.572.482.422.362.322.282.252.222.19各区段刚域挠度各区段刚域挠度多数测点考虑刚域多数测点考虑刚域效
35、效应应后的挠度值更接近实测值后的挠度值更接近实测值。某些测点的某些测点的挠度值在中挠度值在中不断增不断增大中逐渐大中逐渐越越偏离实测值偏离实测值。一些测点理论值在刚域一些测点理论值在刚域扩展扩展后期后期减小至实测挠度值减小至实测挠度值以下以下。与实测值对比：与实测值对比：1.2.3.效验系数对比效验系数对比1.通过节点刚域的模拟，修正模型中的校验系数精度整体提高2%10%，平均增幅为6%。2.9号测点的校验系数修正值超出了规定校验系数1.05，出现不利因素。3.在4、6、11号等测点处，实测值与初始理论值相差较大，经过理论值修正后，仍然没有较大改观。分析：一方面在于挠度实际测量值的误差，另一方
36、面存在于荷载试验中车辆加载具体位置的偏差。测点测点编号编号刚域范围刚域范围0.10.20.30.40.50.60.70.80.91.016.216.296.366.426.476.516.546.566.576.5825.315.185.095.024.954.894.834.784.744.7232.192.142.112.082.052.032.011.991.971.964-2.15-2.14-2.13-2.13-2.12-2.12-2.11-2.11-2.10-2.095-2.94-2.89-2.83-2.76-2.68-2.60-2.52-2.44-2.36-2.2862.202.1
37、62.132.112.082.052.032.011.991.9775.385.295.225.175.125.075.034.994.974.9587.617.487.357.237.116.986.876.776.686.6192.021.991.961.931.911.891.871.861.851.8411-1.26-1.27-1.28-1.28-1.29-1.30-1.30-1.31-1.32-1.3312-1.63-1.58-1.53-1.49-1.46-1.43-1.41-1.39-1.37-1.3513-1.24-1.17-1.10-1.05-1.01-0.98-0.95-0.
38、93-0.92-0.91各区段刚域挠度各区段刚域挠度多数测点考虑刚域多数测点考虑刚域效效应应后的挠度值更接近实测值后的挠度值更接近实测值。某些测点的某些测点的挠度值在挠度值在不断增不断增大中逐渐大中逐渐越越偏离实测值偏离实测值。一些测点理论值在刚域一些测点理论值在刚域扩展扩展后期后期减小至实测挠度值减小至实测挠度值以下以下。与实测值对比：与实测值对比：1.2.3.效验系数对比效验系数对比1.7号和8号测点出现较大测量误差，其误差在合理范围之内，通过节点刚域的模拟，理论修正值与实测值在多数测点实现了更高的吻合度。2.校验系数相比初始值精确度提高6%至27%，平均值为10.1%，与跨中截面工况情况
39、类似。分析：多数测点的挠度校验系数提高至0.70.8范围内，与工程实际结构高度吻合的有限元模型还有一定差距。均匀设计分布均匀设计分布样本样本：将弹性模量与刚域尺寸同时作为可调结构参数。将三种材料的弹性模量参数分别上下调整10%，将节点刚域尺寸范围缩小至0.40.5区段，将弹性模量作为模型修正的主考虑因素，将刚域长度作为副考虑因素。计计算算值值与与实实测测值值对对比比工况三工况三工况四工况四分析：可以认为考虑刚域效应的有限元模型已能较好的反映桥梁的内部真实构造。确定最终有限元模型合理参数组合：确定最终有限元模型合理参数组合：阶数阶数计算频率计算频率实测频率实测频率绝对误差绝对误差相对误差相对误差
40、振型（完全吻合）振型（完全吻合）11.291.3050.0151.15%主梁一阶反对称竖弯21.551.52-0.03-1.97%主梁一阶对称竖弯32.011.98-0.03-1.52%主梁一阶侧弯42.222.24-0.02-0.89%主梁二阶反对称竖弯52.562.40-0.16-6.67%主梁二阶对称竖弯一阶振型图二阶振型图三阶振型图四阶振型图五阶振型图阶数阶数计算频率计算频率实测频率实测频率绝对误差绝对误差相对误差相对误差振型（完全吻合）振型（完全吻合）11.271.305-0.035-2.68%主梁一阶反对称竖弯21.431.52-0.09-5.92%主梁一阶对称竖弯31.951.9
41、8-0.03-1.52%主梁一阶侧弯42.022.24-0.22-9.82%主梁二阶反对称竖弯52.152.4-0.25-10.42%主梁二阶对称竖弯初始模型自振频率计算结果：初始模型自振频率计算结果：修正修正模型自振频率计算结果：模型自振频率计算结果：初始模型平均相对误差为6.072%，修正模型为 2.44%，平均减小3.632%，有了大幅减小。襄阳汉江五桥建模拱肋采用梁单元形式或板单元形式，若单纯做静力分析，两种建模方式均可。索力优化手工迭代法虽然较为麻烦，但是其吊杆力的精确控制更加充分有效，因此建议采用。桩土作用对梁拱组合桥型的静力荷载作用影响较小，因此仅做静力分析时可不予考虑。结论一ANSYS模型计算结果对拱脚处混凝土整体应力水平做了精确的评估，由此结果可以判断，拱脚局部应力水平处于安全状态。并可知以杆系理论为基础的数值分析方法在反映拱脚处的细部应力状态时容易造成误判，因此应采用实体有限云模型进行分析。结论二 类似襄阳汉江五桥的大跨度连续刚构组合体系拱桥中，梁、拱交汇处存在着明显的刚域效应。经过人工神经网络确定的刚域尺寸修正的有限元模型可以对实际大桥结构进行有效评估。由于施工原因、材料老化、外部环境等因素的影响，在对有限元模型中混凝土、拱肋与预应力束的杨氏模量修正后，可作为大桥长期健康监测系统的基准计算模型。结论三 谢谢！
