先来看计算方法,在双对数坐标纸上,以含水量w为横坐标,锥入深度h为纵坐标,点绘a、b、c三点含水量的h-w图,连接此三点,应呈一条直线。如三点不在同一条直线上,要通过a点与b、c两点连成两条直线,根据液限(a点含水量)在hp-wl图上查得Ph,以此hP再在h-w图上的ab及ac直线上求出相应的两个含水量,当两个含水量的差值小于2%时,以该两点的平均值与a点连成一条直线.当两个含水量的差值大于2%时,应重做试验。ywP检测VBA
再来看一张图:ywP检测VBA
ywP检测VBA
ywP检测VBA
ywP检测VBA
ywP检测VBA
ywP检测VBA
上图中紧挨的三条线是三点不在同一直线上但a点与b、c两点连成两条直线含水量的差值小于等于2%时重新画了一条直线,而另一条单独的直线是三点在一条直线上时画出的线。ywP检测VBA
首先我们要根据规范先求出ab及ac直线上求出相应的两个含水量,如下图,假设这是我们试验得出的含水率与锥入深度平均值,ywP检测VBA
根据平面几何的理论,平面上两点确定一直线,由a、b或a、c两点必有一条直线且只有一条直线通过两点。我们可以先分别求得这两条直线的线性回归方程,直线的回归方程一般都是y=ax+b这样的情况,但这里需要注意一个问题,我们的h-w图是一个双对数坐标图,常规的线性回归方法是不适用的,例如我直接用 LOPE函数(斜率)和INTERCEPT函数(截距)求出回归方程中的a和b,再通过在hp-wl图计算出塑限,这时候的塑限锥入深度就是方程的y值,塑限含水率就是方程的x值,代入y值,求出x。这样直接计算得出的结果一定是错误的,我们需要先将前面的锥入深度和含水率以及塑限对应的锥入深度的数值全部转换成对数,ywP检测VBA
如上图,用LOG10函数(返回给定值以10为底的对数)把B2、C2、D2、B3、C3、D4以及利用hp-wl图计算出塑限对应的锥入深度5.42转换成对数, 5.42转换成对数为0.73,利用常规回归计算代入y=0.73得出x=1.22,但显然1.22不可能是塑限含水率,1.22只是代表log10为底x的对数等于1.22,代表x=10的1.22次方,利用Excel计算功能求出10^1.22=16.63这就是塑限含水率。ywP检测VBA
利用上面的方法分别求出由a、b和a、c两条直线对应的塑限含水率,若差值大于2%,则试验应重做,若相等,则三点在一条直线上,若不等(可修约至整数)则计算两条直线对应的塑限含水率的平均值,重新做回归计算。ywP检测VBA
前面的方法比较繁琐,这里再介绍一种比较简洁的方法,如图:ywP检测VBA
我们画出双对数坐标图趋势线时勾选显示公式,可以看到双对数图中直线的回归方程,他的格式是y=axb,这样我们就可以通过两个已知的坐标点得出两个二元一次的指数方程,求解这两个方程,得出a、b的值,再代入塑限对应的锥入深度y值求出x,就是塑限含水率。ywP检测VBA
假设有两个指数方程如下:y1=ax1b ; y2=ax2bywP检测VBA
两式相除,解得:b等于 log以y1/y2为底x1/x2的对数,Excel函数中对数函数LOG计算起来很简单,代入b,求出a。ywP检测VBA
还有一种直接用函数求得a值和b值得方法,如图:ywP检测VBA
求a点和b点之间的系数,公式是这样的,ywP检测VBA
a值=EXP(INDEX(LINEST(LN(B2:C2),LN(B3:C3)),2))ywP检测VBA
b值=INDEX(LINEST(LN(B2:C2),LN(B3:C3)),1)ywP检测VBA
用函数公式算出来的值和解方程求出的值是一样的。ywP检测VBA
求出a、b的值,再求塑限含水率跟前面一样,把塑限锥入深度代入方程求解。ywP检测VBA |