最小二乘法在预应力张拉中应用.docVIP

最小二乘法在预应力张拉中应用【摘要】预应力结构张拉前，配套的千斤顶和油压表等张拉设备必须在试验室标定，试验室所得标准测力仪显示值和压力表读数呈线性关系。在它们的相关系数满足要求时，精确求出它们的线性方程，是确保工程质量的前提。本文以宁波市庆丰桥配套的千斤顶和压力表的标定数据为例，以实际值与趋势值的偏差的平方和最小为依据，根据函数存在的极值条件，求出线性方程，此为最小二乘法。因为在试验室标定时，压力表读数一般是按照等差数列递增的，所以为了计算的简便，可用“零点法”设置等差序号，使这一列数据的总和为零，此时要注意当试验次数分别为偶数和奇数时不同地设置等差序号的方法。在最后的结果中，一定要注意参数关系的转换。最小二乘法在预应力张拉仪器标定和数据处理中加以应用，提高了计算精度，体现了最小二乘法的优越性。 【关键词】最小二乘法；预应力张拉；零点法；等差序号 1 引言 在桥梁工程中，箱梁砼结构的后张法预应力张拉是关键性工程，在张拉前，要先标定油表和千斤顶，理论上压力表读数和标准测力仪显示值呈线性关系，试验室所得的这两组数据应满足相关系数的要求。在现场施工中，有的技术员在计算线性回归方程时以相邻标准测力仪显示值差值的平均值为依据进行计算，虽然也能求出线性方程，但计算出的结果不精确。本人根据试验所得油压表和对应的千斤顶的数据，用最小二乘法来计算验证[1]，即简便，又精确。 2 最小二乘法 最小二乘法，亦称最小平方法。它是根据试验统计数据，通过建立数学模型，以求得线性回归方程的方法。具体有以下步骤： 2.1 设所求回归方程 设所求线性回归方程为=。 式中：――标准测力仪显示值 ――压力表读数 ，――常数 2.2 求出和的值 从数学意义上看要使=这一直线方程具有较高的代表性，就必须使实际值与趋势值之间的偏差（）的平方和最小[2]，即： 这种确定，使偏差平方和最小的方法就叫做最小二乘法。 上式中的和全都是常数(代表压力表读数，代表标准测力仪显示值)，所以是和的二元函数，为了求的最小值。就要求出两个偏导数： 根据极值原理，即极值存在的必要条件[3]，令： 将上式整理，即可得出关于，的方程组： 方程组简写为： （1） 由该方程就可以求出，的数值[4]。 式中：―――试验次数 ―――压力表读数 ―――总和符号 2.3 得出所求的方程并检验计算结果 把，值代入=，求得具体的一元一次线性回归方程，并对其做检验。 2.4 用“零点法”计算简便 注意，用“零点法”求值，为简便计算，引用等差序号,使得，根据方程(1)，把当成0，得出下式计算式： （2） [5] 式中：―――等差序号 得出一元一次方程 （3） 而是的一次函数，容易得到 （4） 将(4)代入(3)，求得 得到， ，（均为常数） 3 举例说明 宁波市庆丰桥某千斤顶和油压表标定试验数据如表1： 表1标定试验数据 压力表读数(Mpa) 标准测力仪显示值（KN） 5 254.3 10 480.1 15 696.3 20 916.6 25 1118.3 30 1326.0 35 1529.4 40 1736.8 45 1936.9 50 2133.6 55 2336.9 解：设所求线性方程为 = 3.1 第一种解法： 列表计算如表2： 表2 计算表 试验次数n 压力表读数(Mpa)x 标准测力仪显值（KN）y xy X― 1 5 254.3 1271.5 25 2 10 480.1 4801 100 3 15 696.3 10444.5 225 4 20 916.6 18332 400 5 25 1118.3 27957.5 625 6 30 1326.0 39780 900 7 35 1529.4 53529 1225 8 40 1736.8 69472 1600 9 45 1936.9 87160.5 2025 10 50 2133.6 106680 2500 11 55 2336.9 128529.5 3025 总计 330 14465.2 547957.5 12650 代入方程组(1)， 得， 所求方程为 3.2 第二种解法： 列表计算如表3： 表3引入等差序号的计算表 试验 次数n 压力表读数（Mpa） 标准测力仪显示值（KN） 等差序号 1