回归分析在界限含水率试验数据处理中的应用.docxVIP

第 27 卷第 3 期2010 年 6 月华 东 交 通 大 学 学 报Journal of East China Jiaotong UniversityVol. 27No. 3Jun. ,2010文章编号 :100520523 (2010) 0320027206回归分析在界限含水率试验数据处理中的应用赵秀绍1 ,莫林利2(1. 华东交通大学 铁路环境振动与噪声教育部工程研究中心 ,江西 南昌 330013 ; 2. 华东交通大学 软件学院 ,江西 南昌 330013)摘要 :常用的界限含水率数据处理方法中 ,作图法不够精确 ,公式法仅能同时处理三组试验数据 。针对以上特点 ,提出了线 性回归计算界限含水率的优化算法 。该算法可以充分利用平行试验的多组数据 ,使误差达到最小 。重点阐述了用最小二 乘法计算界限含水率的基本原理 ,给出了程序设计流程图 ,并用 Visual Basic 语言实现了这种算法 ,最后利用典型工程实例 分析了该算法的优缺点 。关 键 词 :界限含水率 ;回归分析 ;优化算法 ;Visual Basic中图分类号 : TP311 + TU411文献标识码 :A界限含水率试验所得结果液限 wL和塑限 wp 是确定细粒土的承载能力 、对细粒土进行填料分组的重 要指标 ,所以准确地确定土的液塑限指标对工程建设具有重要的意义 。界限含水率试验测定完成后 ,目前常用的数据处理方法有作图法和公式法 。作图法受人为因素影响 , 估读数据时一般会参照常规坐标 ,造成读数错误 ,其次每次试验都靠手工完成 ,不能实现数据处理自动化 。 近年来许多学者提出了一些新的算法 ,例如 : 殷春娟提出了用 Excel 表格处理界限含水率的方法[ 1 ] ,彭 意[ 2 ] ,冯朝辉[ 3 ] ,张传邦[ 4 ]提出了公式法计算界限含水率的方法 ,高盟[ 5 ]提出了用样条函数求解的方法 ,这 些方法大大的提高了计算的效率与精度 。以上方法能通过塑限差法 (Δwp 2 %) 控制试验的数据偏差与 精度 ,但仅能处理三组数据 ,部分试验数据被迫舍去而得不到充分利用 ,往往导致试验要制样 5～9 次进行 测试 ,而有效的仅有 3 次 。在前人工作方法的基础上 ,提出了基于线性回归的界限含水率计算方法 ,充分 利用了各数据点对试验结果的贡献 ,用 Visual Basic 语言实现了这种算法 ,并在实践中检验了该算法的合 理与优越性 。1作图法计算界限含水率根据文献[ 6 ] 可得 : 在液塑限联合测定试验中 , 锥尖入土深度 h (mm) 与含水率 w ( %) 之间存在如 式 (1) 之间的双对数线性关系 :lg h = alg w + b(1) 式中 : a 和 b 是由试验确定的常数 , 根据以上原理 , 界限含水率试验的作图法应采用双对数坐标系 。例 如《TB10102 - 2004 铁路工程土工试验规程》[7 ] 或土 工试验方法标准[8 ] 法 , 以 w 为横坐标 , h 为纵坐标 , 在双对数坐标纸上绘制 , 三组数据 ( A 、B 、C) 应连成 一条直线 。三点明显不在一条直线上时 , 连接 AB和 AC 并延长 , 与 h = 2 mm 直线有两个交点 D1 、D2 ,图 1 作图法原理图收稿日期 :2010201212基金项目 :江西省青年科学基金项目 ( GJJ09503) ,华东交通大学博士科研启动基金项目 ,华东交通大学科研 基金 (09RJ05)作者简介 :赵秀绍 (1978 - ) ,男 ,博士 ,副教授 ,研究方向 :计算机在岩土工程中的应用 。28华 东 交 通 大 学 学 报2010 年两个交点横坐标 ( wAC - wAB) 之差称为塑限差 , 当Δwp 2 %时 , 取两个交点的中点 D , 则 AD 就是表达式(1) 的线性关系所绘制的直线 。当Δwp ≥2 %时说明试验差距过大 , 应重做试验 (见图 1) 。 作图法理论上可以计算和控制试验的精度 , 但实际操作起来相当困难 , 因为Δwp ≥2 %在双对数坐标上很难估读准确 。作图法的缺点 : ①受作图精度的影响 ; ②受人为因素的影响 ; ③作图法要求三个点之间的距离应尽 量拉开 , 例如规程[7 ]规定 , C 点的范围在 3～5 mm , B 点范围在 9～11 mm , A 点范围在 16～18 mm。但实际 测定时 , 可能多次测试均不在以上范围内而需重新配水测定 , 造成试验数据浪费与人力资源浪费 。2公式法计算界限含水率依据作图法的原理 (见图 1) , 分别求出 AB 与 AC 直线的斜率 KAB和 KAC , 则 K = ( KAB + KAC) /