随着国家对基建项目的大量投入,箱梁桥由于良好的受力性能而广泛应用。本文主要研究了箱梁的约束扭转问题,首先将变截面连续箱梁离散化,即变成等截面的梁段单元,采用能量变分法建立微分方程,并确定边界的约束条件。然后用初参数法求解微分方程,求出箱梁任意截面的内力与位移。再根据刚度系数的一般定义,推导了单元刚度矩阵和等效节点荷载。接着参考平面杆系结构的一般有限元分析思路,运用Fortran语言自编计算程序,求出箱梁任意截面的扭矩与双力矩。最后结合团结河大桥主桥部分,运用自编程序计算分析,计算出弯曲应力、翘曲正应力以及约束扭转剪应力,并进一步计算出应力放大系数。与有限元软件相比,自编程序可求出梁全截面上的内力,而有限元软件只能求出梁的局部应力。结合数值算例,本文的主要结论如下:（1）考虑悬臂板和顶板的厚度不同,推导了主扇性坐标、扇性静矩等扭转几何特性的计算公式。通过引入剪应力系数研究总剪应力受二次剪应力的影响。用本文解析法和ANSYS实体单元计算了简支箱梁算例的翘曲正应力和剪应力,并研究了悬臂板宽度和箱室高度变化对广义扇性静矩、二次剪应力及剪应力系数的影响规律。研究结果表明:解析解与ANSYS解吻合良好,验证了公式的正确性;当悬臂板宽度比约为0.3时,腹板剪应力最大点处的剪应力系数可达到最大值2.4;随着箱室高宽比的增大,特征点处的二次剪应力逐渐减小;随着悬臂板宽度比和箱室高宽比的增大,顶、底板中点处的广义扇性静矩呈现出完全不同的变化规律。（2）通过研究三种不同荷载工况下的扭转效应,可得双力矩、扭矩以及扭转角的分布规律为:双力矩在支座处和集中扭矩处达到最大值;扭矩在集中扭矩处会发生突变;扭转角在集中扭矩处达到最大值,在支座处为零。（3）同一荷载工况下,双车道的应力放大系数普遍小于单车道的应力放大系数。不同荷载工况下,正应力放大系数最大值的位置不固定。工况一,在中跨1/4截面处有最大值;工况二,在中跨的3/4截面处有最大值;工况三,在左支座截面处有最大值。（4）剪应力放大系数在支座截面上顶板与腹板交点处达到最大值。工况一,单车道（1.6022）,双车道（1.4127）;工况二,单车道（2.1741）,双车道（1.8046）;工况三,单车道（1.6919）,双车道（1.4741）。