二次曲线指数函数两点式拟合回归方程在混凝土凝结时间算法中的应用.pdfVIP

二次曲线、指数函数、两点式拟合回归方程 在混凝土凝结时间算法中的应用 张 龙 安映东 曾庆文 （中国水电三局有限公司勘测设计研究院） [摘要] 混凝土凝结时间的测定，试验技术人员按照规范步骤实施，但是对实验数据 的后续处理，目前国内规程规范或行业标准没有统一的数学方法，因此会导致不同的数 据处理结果。本文介绍了二次曲线、指数函数、两点式等拟合回归方程在混凝土凝结时 间算法中的应用，并且用Microsoft Excel中LINEST、CORREL、MINVERSE等函数进行 示例解算。进一步促进实验数据科学化、规范化处理。 [[[关键词]] 混凝土 凝结时间 函数 回归方程 数据处理[ ]] 一、引言 混凝土试验过程中，应严格按照相关规程操作。但是对实验数据的处理，以求平均 值、加权平均值或读趋势线等作简单化处理，这样往往会使严谨的试验数据，经粗放式 的数据处理，得到粗略结果。那么回归分析、方差分析、插值法等数学方法的应用，会 使问题较好解决。 混凝土凝结时间利用贯入阻力法测定，绘制初凝时段和终凝时段贯入阻力与时间关 系曲线，以3.5MPa及28MPa划两条平行横坐标的直线，直线与曲线交点对应的横坐标 值即为初凝时间与终凝时间。那么受作图与读图影响，结果自由度较大。为此，我们可 以对所测数据回归分析，建立凝结时间与贯入阻力关系的回归方程，利用解析法直接求 解。本文分别利用线性方程、指数函数、二次曲线等数学模型，对所测数据进行拟合， 得到回归方程。在解算过程中涉及矩阵求逆等问题，因此又采用MicrosoftExcel函数， 使得运算简便，可操作性强。 二、混凝土拌合物凝结时间计算 1、线性回归计算 所谓线性回归方法，就是将贯入阻力p与时间t建立线性关系。 Lnp=A+Blnt （2.1.1） 式中 t—时间，min； p—贯入阻力，MPa A、B—线性回归系数 由于凝结时间与贯入阻力呈曲线图，不能直接建立线性关系。因此对贯入阻力p及时间 t求自然对数，这样会使其线性关系明显，在建立直线方程。 根据式（2.1.1）求得当贯入阻力3.5MPa时为初凝时间t，贯入阻力28MPa时为终s 凝时间t，求解得e ln3.5 -A ts = e B （2.1.2） ln28 -A te = e B （2.1.3） 线性回归计算实例： 已知贯入阻力(MPa)及累计时间t（min），对应的自然对数，lnp及lnt，见表1； 表1 混凝土拌合物凝结时间试验结果 分别绘制了凝结时间与贯入阻力关系曲线图1，凝结时间自然对数（lnt）与贯入阻力 自然对数（lnp）关系曲线图2。前着明显呈曲线状，不能直接建立直线方程，后者则 近似于直线，可建立直线方程。 图1 凝结时间与贯入阻力关系曲线图 图2 凝结时间lnt与贯入阻力lnp关系曲线图 将lnt作为横坐标，lnp作为纵坐标，代入方程（2.1.1），用最小二乘法求回归系数A 及B。利用Microsoft Excel中函数求解。 B11=INDEX(LINEST(C2:C10,D2:D10),1),即B=6.6359 A11=INDEX(LINEST(C2:C10,D2:D10),2),即A=-43.2757 将A，B代入式（2.1.2）、式（2.1.3）中； A12=EXP((LN(3.5)-A11)/B11)，即t=821mins B12=EXP((LN(28)-A11)/B11)，即t=1123mine C11=CORREL(C2:C10,D2:D10),即相关系数R=0.99702 从求解结果看，对混凝土凝结时间及贯入阻力求自然对数，然后利用最小二乘法拟 合为最佳直线方程，其相关系数接近1，说明建立的直线方程比较