我国南水北调工程中渡槽应用十分广泛，由于该工程中渡槽多为自流输水，水头损失较小，而输水流量巨大，普通梁式渡槽很难满足要求，因此对渡槽结构提出了更高要求。为了减小槽身底板厚度，设计时将槽身分成两槽或三槽并联的形式，使其纵梁达到三或四条，改变了以往利用大侧墙承重的结构形式，减少了上部结构的重量，也减小了渡槽槽身各构件的尺寸，使渡槽槽身承载能力大大增加，纵向跨越能力也随之增大，槽身结构形式和受力趋于优化。在水荷载作用时，由于结构的整体作用，各纵梁不同程度地都要产生挠曲变形而形成一个挠曲面，显示结构变形和受力的空间性，从而造成水荷载在纵梁间的分配非常复杂，结构静力计算的难度相应增加，渡槽的这种受力特性实际上已属于空间结构力学分析范畴。
一般认为，中纵梁承担一个单槽的水荷载，而边纵梁只承担半个槽的水荷载，理论上如果边纵梁及中纵梁的刚度与其承受的荷载相对应，他们的挠度变形将相同，横梁支座及荷载作用的不同点处将不会产生相对变形，这在实际工程中很难做到的。一般情况下，挠度较大的中纵梁的（水荷载为边纵梁的2倍）部分荷载必然通过底肋以剪力形式传递给边纵梁，从而中纵梁分配的水荷载在减小，而边纵梁分配的水荷载在增大，使渡槽结构达到平衡。
目前，多槽式渡槽结构的静力分析方法还不成熟，竺慧珠等采用传统结构力学方法，在渡槽横向受力计算时，没有考虑边纵梁与中纵梁的不均匀变形，而把纵梁看成横梁的固定支承是不合适的；赵顺波等研究中侧肋和底肋组成横向受力框架，纵梁和横梁结点为固定支承，各纵梁不均匀变形按结构力学中支座沉降进行槽身横向受力分析，其支座的不均匀沉降需要通过三维实体有限元计算；王云仓考虑了中纵梁传递给边纵梁的次生荷载和横梁产生的次生弯矩；赵顺波等采用了改进的空间刚架结构计算法和三维混凝土有限元法分析了南水北调工程中双洎河渡槽。本文以南水北调某大型三槽式矩形渡槽为例，运用结构力学法和有限元法对渡槽槽身空间结构受力进行分析，获取渡槽槽身纵横向计算的简化方法和计算公式。
1渡槽槽身结构计算模型简化
多槽式矩形渡槽槽身受力分析过程中，将槽身视为带翼缘板的横梁和纵梁交叉组成的空间结构，分析时按平面问题简化为横梁和纵梁。横向简化为由侧肋、中墙和底肋组成的受力框架，承受水荷载作用，纵梁和横梁的结点为弹性支承；纵向简化为一简支梁，承受横梁弹性支承传递的荷载。渡槽槽身结构计算简化模型见图1。
图1渡槽槽身结构计算简化模型
Fig.1Simplified model of the calculation of aqueduct structure
2横梁弹性支承刚度计算
根据结构力学中位移计算原理，取单位均布荷载作用在纵梁上，求横梁位置处截面的挠度。根据弹簧变形原理，该截面的弹性支座刚度即为其挠度的倒数。由于各纵梁高跨比很大，挠度计算时要同时考虑弯矩和剪力的共同作用。
结构力学中平面受弯构件在荷载作用下的位移计算公式为：
Δp=∑∫MMpdx EI+∑∫kFsFspdx GA（1）
式中：E为材料的弹性模量；I和A分别为杆件截面惯性矩和面积；G为材料的切变模量；k为切应力沿截面分布不均匀而引用的修正系数，其值与截面形状有关，矩形截面k=6 5，工字型截面k≈A A′、A′为腹板截面面积。
纵梁上作用一单位均布荷载，其弯矩和剪力分布为：
Mp（x）=1 2xL-1 2x2Fsp（x）=1 2L-x（2）
式中：L为简支梁支座间距（m）；x为计算截面距离支座的距离（m）。
纵梁上在计算截面上作用一单位集中荷载，其弯矩和剪力分布为：
M（x）=x（1-a L），x≤a
a（1-x L），x≥a
Fs（x）=1-a L，x≤a
1-a L，x≥a
（3）
式中：a为纵梁上横梁位置距支座的距离（m）。
将（2）、（3）式代入（1）式，化简整理后得到纵梁在单位均布荷载作用下横梁位置处的竖向位移，即为纵梁在该截面的柔度系数：
Ca=1 24EI（a4+L3a-2La3）+k 2GA（La-a2）（4）
纵梁弹性支座的刚度即为：
ka=1 Ca（5）
3计算截面选取
渡槽槽身横向计算时，简化为弹性支承在各主梁上的多跨弹性支承连续框架，在跨中截面处，弹性支承刚度最小；当截面由跨中向支座处靠近时，弹性支承刚度增大，在支座处的位移为零，此时端横梁近似为固定支承连续梁。因此，在计算时沿跨径方向选取槽跨跨中、1/4、1/8和支座处为特征截面进行计算分析。
4渡槽槽身横向结构弹性支承反力计算
双槽式渡槽槽身底板受均布水荷载为q，侧墙所受水荷载作用在横梁梁端弯矩为M，中纵梁对横梁的弹性支承刚度为k2，采用力法进行分析，见图2。
图2槽身横向受力计算示意图
Fig.2Schematic diagram of aqueduct lateral force calculation
力法典型方程为：
δ11X1+Δ1p=Δ1=-（X1 k2+Δ1Δ）（6）
式中：δ11=∑∫M21ds EI=l3 6EI；Δ1p=∑∫M1Mpds EI=1 EI（1 2Ml2-5 24ql4）；
Δ1Δ=-∑FR1c=-ql k1+X1 2k1
将以上系数代入典型方程（6）中，化简计算得到：
X1=5 12ql4 EI-Ml2 EI+2ql k1 l3 3EI+1 k1+2 k2（7）
令1 k1=C1，1 k2=C2，C1、C2为弹性支座的柔度系数。
则，双槽式渡槽弹性支承反力：
R2=X1=5 12ql4 EI-El2 EI+2qlC1 l3 3EI+C1+2C2（8）
R1=R3=ql-1 2R2（9）
式中：q为作用在横向结构上的均布水荷载；l为各纵梁的间距；C1为边纵梁与横梁交叉点弹性支座的柔度系数；C2为中纵梁与横梁交叉点弹性支座的柔度系数。
同理可得三槽式渡槽的弹性支承反力为：
R2=R3=11ql4 12EI-Ml2 EI+3 2qlC1 5l3 6EI+10 9C1+C2（10）
R1=R4=3 2ql-R2（11）
5槽身纵向计算
作为纵向、横向和高度3个方向整体受力的结构
，对于其所受荷载不是在各纵梁之间进行简单平均分配，而是依靠相互协调平衡荷载的作用。所以通过横向结构分析中得到各纵梁对横梁的弹性支承反力，即可求出水荷载沿横向在各纵梁之间的分配系数。将其对应水荷载作用在各纵梁上，对纵梁进行受力分析。从支座到跨中，弹性支承的刚度逐渐变小，使弹性支承的反力在槽跨不同计算截面处产生变化，因此各纵梁上分配的水荷载沿跨径方向的分布是不均匀的。
6算例分析
6.1工程背景
南水北调中线某大型三槽式渡槽，单槽断面60 m×53 m，设计水深4792 m。边墙厚06 m，顶部设20 m宽的人行道板，中墙厚07 m，顶部设27 m宽的人行道板。槽身设侧肋和底肋，间距为25 m。侧肋宽05 m，高07 m。底肋宽05 m，高11 m。槽身纵向为多跨简支梁，跨径为30 m。槽身为预应力混凝土结构，强度为C50。
6.2渡槽有限元法计算结果
采用有限元软件MIDAS进行计算，对槽身进行单元划分，共有9 444个结点，29 275个单元。渡槽有限元模型见图3。弹性模量E=345×104MPa，泊松比λ=02。水荷载作用下槽身有限元计算结果见表1、表2。
图3有限元模型
Fig.3Aqueduct FEM model
6.3渡槽结构力学法计算结果
槽身结构中，边纵梁截面惯性矩Ib=52.06 m4，中纵梁的截面惯性矩Iz=65.91 m4，横梁的截面惯性矩为I0=1.42 m4。根据上述原理计算分析，所得结果见表3至表5。
6.4结果分析
（1）同一截面处中纵梁位移和底缘应力均大于边纵梁，说明中纵梁所受水荷载大于边纵梁所受水荷载。
（2）从支座到跨中，弹性支承的刚度逐渐变小，使弹性支承的反力在槽跨不同计算截面处产生变化，水荷载在各纵梁上的分配情况为：在跨中截面，中纵梁分配的水荷载是边纵梁的1.5倍，随着向梁端截面的靠近，中纵梁分配的水荷载在增加，而边纵梁分配的水荷载在减少，在梁端截面处，中纵梁分配的水荷载大约是边纵梁的2倍。
（3）一般认为，边纵梁承担半个槽的水荷载，中纵梁承担一个单槽的水荷载，如果这样对槽身进行纵向受力计算分析，对边纵梁的计算是偏于不安全的，对中纵梁的计算是偏于保守的。
（4）边纵梁从支点到1/4截面的荷载分配系数逐渐增大，近似呈直线分布，从1/4截面到跨中截面的荷载分配系数基本不变；中纵梁从支点到1/4截面的荷载分配系数逐渐减小，近似呈直线分布，从1/4截面到跨中截面的荷载分配系数基本不变。所以在横向计算时沿跨径选取槽跨跨中和支座处为特征截面进行计算分析就可以了。
（5）边横梁为固定支撑的连续梁，支点截面处为负弯矩，下缘受压；跨中截面处为正弯矩，下缘受拉。
（6）中横梁为弹性支承的连续梁，在靠近边纵梁处为负弯矩，下缘受压，随着远离边纵梁，弯矩由负变正，下缘由受压向受拉转变。其与固定支承的连续梁不同，不能把刚性支承连续梁上弯矩分布的概念用以说明中横梁的内力情况。
（7）结构力学法在槽身结构计算过程中，没有考虑梁的扭转问题，所以与有限元法计算稍有差异。
7结语
进行渡槽槽身空间结构计算时，将渡槽槽身简化为平面计算模型，从横向结构和纵向结构单独分析，用弹性支座使横向和纵向结构产生联系。横梁弹性支承在各纵梁上，通过计算弹性支座反力，获得水荷载横向分配系数，将其对应水荷载作用在各纵梁上，以此进行渡槽槽身结构受力分析，计算简便，结果准确，便于实际工程的设计配筋，计算精度也能满足实际工程的要求。
参考文献（References）：
竺慧珠，陈德亮，管枫年.渡槽.北京：中国水利水电出版社，2004.（ZHU Huizhu，CHEN Deliang，GUAN Fengnian.Aqueducts .Beijing：China WaterPower Press，2004.（in Chinese））
赵顺波，李晓克，赵平.混凝土叠合结构设计原理与应用.北京：中国水利水电出版社，2001.（ZHAO Shunbo，LI Xiaoke，ZHAO Ping.Design principles and Applications of Concrete Composite Structure .Beijing：China WaterPower Press，2004.（in Chinese））
王云仓.三向预应力多侧墙渡槽结构研究.天津大学，2009.（WANG Yuncang.Structural Research on Three Dimension Prestressed MultiSide Wall Aqueduct .Tianjin University，2009.（in Chinese））
赵顺波，李晓克，赵平.大型钢筋混凝土多纵梁渡槽结构设计方法的研究.水利学报，1999（4）：3539.（ZHAO Shunbo，LI Xiaoke，ZHAO Ping.On the Design Method for Large Reinforced Concrete Aqueduct with Multilongitudinal Beam .Journal of Hydraulic Engineering，1999（4）：3539.（in Chinese））
李廉锟.结构力学.北京：高等教育出版社，2004.（LI Liankun.Structural Mechanics.Beijing：Higher Education Press.（2004）.（in Chinese））
季日臣，陈尧隆.大型多纵梁矩形渡槽结构静力计算方法研究.水力发电学报，2005，24（5）：8084.（JI Richen，CHEN Yaolong.The Study on Static Structure Calculation Methods of Large Rectangle Aqueduct with Longitudinal Beams.Journal of hydroelectric Engineering，2005，24（4）.8084.（in Chinese））
季日臣，陈尧隆，王春国.大型多纵梁矩形渡槽受力计算的实用空间法.西安理工大学学报，2005，21 （3）：257260.（JI Richen，Chen Yaolong，WANG Chunguo.The Practicality Space Analysis Method of Structure Calculation for Large Rectangle Aqueduct with MultiLongitudinal Beam .Journal of Xi’an University of Technology，2005，21（3）：257260.（in Chinese））
季日臣，张永亮，王军玺.大型多纵梁矩形渡槽槽身横向结构计算.兰州交通大学学报，2004，23（4）：58.（JI Richen，ZHANG Yongliang，WANG Junxi.Orientation Structure Calculation Method of Large Rectangle Aqueduct with Multilongitudinal Beam .Journal of Lanzhou Jiaotong University，2004，23（4）：58.（in Chinese）） 本文由第一论文网摘自《南水北调与水利科技》2014年第2期，版权归原作者和期刊所有。
