王永逵，王健，耿加會，余春榮，馮立艷

（1. 上海奇齊科技開發(fā)實業(yè)有限公司，上?！?01600；2. 河南舞陽縣惠達公路工程公司，河南　舞陽　462400；3. 建筑材料工業(yè)技術(shù)情報研究所，北京　100024；4. 四川德陽明鴻商品混凝土有限公司，四川　德陽　618099）

從混凝土歷史發(fā)展的長河中學(xué)習(xí)、理解A定則B公式之二（II）

王永逵1，王健1，耿加會2，余春榮3，馮立艷4

（1. 上海奇齊科技開發(fā)實業(yè)有限公司，上海201600；2. 河南舞陽縣惠達公路工程公司，河南舞陽462400；3. 建筑材料工業(yè)技術(shù)情報研究所，北京100024；4. 四川德陽明鴻商品混凝土有限公司，四川德陽618099）

阿勃萊姆斯的 A 定則是在費雷特普適公式基礎(chǔ)上的發(fā)展和提高，對混凝土強度理論具有里程碑的意義。保羅米的 B公式是在 A 定則基礎(chǔ)上的簡化，使數(shù)理統(tǒng)計對混凝土中的應(yīng)用進入新階段。但與A定則沒有質(zhì)的區(qū)別，都是對混凝土強度期望值估計值的計算，其結(jié)果應(yīng)該很接近。A 定則—B 公式成了現(xiàn)代混凝土的基礎(chǔ)理論之一，有很強的現(xiàn)實指導(dǎo)意義。

A 定則—B 公式；期望值；預(yù)報；控制

（接前節(jié)）

3.4B 公式的電腦擬合

上述列表計算的方法可以作為學(xué)習(xí)、掌握一元線性回歸分析的基本訓(xùn)練，可以看出計算還是很繁瑣的，僅子樣為27對﹙Xi，Yi﹚數(shù)據(jù)的小樣本，計算數(shù)據(jù)約300個次上下，而且計算中不能出任何差錯，假使有一個錯誤，就不可能得到正確的結(jié)果?？梢娫谀莻€僅有手搖計算機的年代，擬合一個用于生產(chǎn)和科研的公式多么不易。如果你能動手親自算一算，你就體會到我們今天重溫歷史，學(xué)習(xí)當(dāng)年的成果，對費雷特、阿勃萊姆斯和保羅米他們所做的開創(chuàng)性工作是多么值得敬重?，F(xiàn)在好了，幾百個幾千幾甚至上萬個數(shù)據(jù)，甚至更大的樣本，用身邊的電腦，僅用幾分鐘的時間就準(zhǔn)確無誤的得到需要的結(jié)果。

3.4.1電腦擬合 B 公式的步驟

（1）在 Excel 文檔中列表，將膠水比﹙Xi﹚和混凝土實測強度﹙Yi﹚按由小到大順序排列；

（2）選中 Xi和 Yi后，點擊“插入”中的“圖表”，點擊相應(yīng)的“散點圖”，點“完成”；

（3）選擇“散點圖”里面的一個“點”鼠標(biāo)右擊，在顯示界面中選擇【添加趨勢線】，如本例的直線方程；

（4）在【添加趨勢線】界面里面點直線方程，在顯示界面中點【選項】進入界面；

（5）在界面【選項】里面勾選【顯示公式】和【顯示 R均方差】，即可。

完成圖中除直接標(biāo)出散點和線型外，同時標(biāo)出一元一次的線性方程。

電腦擬合公式為：R28＝26.206－18.612(B/W) 與列表計算的十分接近，但兩者的計算過程相比后者要簡單多了，即使對幾千或數(shù)萬對的（Xi，Yi）試驗數(shù)據(jù)的回歸方程擬合，也就是上面的幾個步驟，就可在很短時間內(nèi)，得到你要擬合的一元回歸方程，并能同時得到你要檢驗的回歸方程與理論曲線的相關(guān)系數(shù)檢驗方差值 R2，大小與相關(guān)系數(shù)檢驗r值相同，可視為相關(guān)關(guān)系很好。

圖3　電腦計算擬合 B 公式小直線圖

以上是以筆者在商品混凝土近二十幾年技術(shù)工作中，積累的部分膠水比與強度的樣本的子樣數(shù)據(jù)（Xi，Yi），并用列表和電腦分別擬合的 B 公式。僅27組數(shù)據(jù)對，樣本很小，有可能受積累數(shù)據(jù)的局限性，影響公式的合理性，應(yīng)盡可能接近事件的母體，從更大范圍收集數(shù)據(jù)，用更大的樣本擬合B公式似乎更具合理些。列寧還說過“要從事實的整體，從事實之間的聯(lián)系中去掌握事實……如果不是從整體上，不是從聯(lián)系中去掌握事實，如果事實是零碎的和隨意挑出來的，那么它們就只能是一種兒戲，或者連兒戲也不如?！本褪菑娬{(diào)要從現(xiàn)實生產(chǎn)的整體﹙稱為“母體”﹚或更大數(shù)據(jù)中去尋求事實的規(guī)律才更有說服力。為此，筆者以十年來《商品混凝土》雜志發(fā)表的論文為基礎(chǔ)，收集水膠比與強度數(shù)據(jù)子樣（Xi，Yi）249對。這其中有你、有我、也有他（她），可以說是業(yè)內(nèi)各位專家學(xué)者同行十幾年的共同積累，雖不能稱之為總體，甚至也說不成是大樣本，但在現(xiàn)有條件下與27對子樣的樣本相比，是相對的較大樣本了?，F(xiàn)在以此249個的膠水比為自變量（X）28天抗壓強度為因變量（Y），用電腦重新擬合的大樣本保羅米公式，試看其與小樣本和規(guī)范及實際統(tǒng)計的差別（為了區(qū)別標(biāo)注 B 公式大與 B 公式?。Y(jié)果如下：

3.4.2“大樣本”的 B 公式電腦擬合

對比小樣本中的式（22）小、（24）、（25）數(shù)據(jù)，與用較大樣本建立的B公式大相比基本一致。常數(shù)項回歸系數(shù)K1，K2的偏差都不足2%。對 Y28估算和 B/W 的估算值計算沒有影響。從大小兩 B 公式對估計值計算的比較，自然是大樣本的 B 公式大具有更大的覆蓋率。兩者的相關(guān)系數(shù)相比，B 公式大的 R2＝0.9419；B 公式小的 R2＝0.9587，數(shù)字相比有少許差別，查自由度 γ＝249－2＝247，臨界值 γn＝0.1615?，F(xiàn)｜γ｜＝0.9419＞＞0.1615，∴方程式（26）大與式（22）小都屬于特別顯著。

3.4.3以3.2C30為例分析混凝土 R28隨機變量分布的特征

從表5中，列出249個以水膠比為自變量（X）和28天抗壓強度為因變量（Y）原始數(shù)據(jù)的目的，借以觀察混凝土抗壓強度的不確定性數(shù)據(jù)的分布特點：在給定水膠比條件下，對應(yīng)的并非一個且并不相同。如表中給定 W/B＝0.46的條件下，序號123～145中為水膠比皆為0.46的混凝土抗壓強度值有24個，，與 C30用式（26）計算的 C30＝38.2MPa 相差甚少，且為正態(tài)分布，曲線略有點偏，這是因為樣本小，遠(yuǎn)沒達到接近母體要求的原因，抗壓強度值均方差為 σ＝2.94MPa。

圖4　電腦計算擬合 B 公式大直線圖

表5　“大樣本”的采集

序號W/BB/W強度序號W/BB/W強度序號W/BB/W強度序號W/BB/W強度330.6001.66724.5960.4702.12838.61590.4402.27342.52220.3402.94166.4340.6001.66728.0970.4702.12836.51600.4302.32639.52230.3303.03068.5350.6001.66725.4980.4702.12838.21610.4302.32644.02240.3303.03072.3360.5901.69523.8990.4702.12836.81620.4302.32643.52250.3303.03059.6370.5801.72428.21000.4702.12837.51630.4302.32642.62260.3203.12558.2380.5701.75427.41010.4702.12836.81640.4302.32642.12270.3203.12557.6390.5601.78631.91020.4702.12839.71650.4302.32640.22280.3203.12558.4400.5501.81831.21030.4702.12840.91660.4302.32638.22290.3203.12559.7410.5501.81832.51040.4702.12836.51670.4302.32642.72300.3203.12565.8420.5501.81828.91050.4702.12840.11680.4302.32642.52310.3103.22668.6430.5401.85233.41060.4702.12837.81690.4202.38143.72320.3103.22670.5440.5401.85230.71070.4702.12835.91700.4202.38144.02330.3103.22662.8450.5401.85232.31080.4702.12849.51710.4202.38148.52340.3103.22665.8460.5401.85235.11090.4702.12836.81720.4202.38145.12350.3103.22672.8470.5401.85232.71100.4702.12838.71730.4202.38144.82360.3103.22663.2480.5401.85234.91110.4702.12836.81740.4202.38145.62370.3103.22662.2490.5301.88728.71120.4702.12835.71750.4202.38147.22380.3003.33367.9500.5301.88729.51130.4702.12838.41760.4202.38141.22390.3003.33365.8510.5301.88732.11140.4702.12833.41770.4202.38140.62400.3003.33365.8520.5301.88734.71150.4702.12838.01780.4102.43942.92410.3003.33363.1530.5201.92329.21160.4702.12842.11790.4102.43946.02420.3003.33365.8540.5201.92329.61170.4702.12841.81800.4002.50044.42430.2903.44865.9550.5201.92327.81180.4702.12840.41810.4002.50045.32440.2903.44868.2560.5101.96133.21190.4702.12834.71820.4002.50047.62450.2903.44866.2570.5101.96130.81200.4702.12840.31830.4002.50044.92460.2903.44874.1580.5101.96129.51210.4702.12835.21840.4002.50045.22470.2903.44872.8590.5101.96133.81220.4702.12838.11850.4002.50046.72480.2703.70478.2600.5101.96133.81230.4602.17442.21860.3902.56446.82490.2703.70472.6610.5101.96135.11240.4602.17435.01870.3902.56445.2620.5101.96132.11250.4602.17435.51880.3902.56446.6630.5101.96128.31260.4602.17434.01890.3902.56450.0

圖5　C30強度正態(tài)分布

圖5平均值37.8MPa 出現(xiàn)的概率最高約為50%。

3.5B 公式的擬合可以用水膠（灰）比近似成反比嗎？

瑞典科學(xué)家保羅米在170多年前將美國伊利諾斯大學(xué)的達夫?阿勃萊姆斯﹙Duff·Abrams﹚教授建立的混凝土強度與其水灰比成雙曲線函數(shù)關(guān)系的A定則，簡化為成的線性關(guān)系，從而使A定則的應(yīng)用更為簡單，這無疑是個進步。但自然會聯(lián)想到，把阿波萊姆斯用的強度因子水灰比改為灰水比的做法有這個必要嗎？這與習(xí)慣應(yīng)用的水灰比，似乎多一層麻煩。因此，B 公式也許還可以用另一種形式。那就是混凝土強度與水灰比是否可成反比的一元線性回歸方程？方法同上，直接以水膠比為自變量（Xi）與對應(yīng)強度為因變量（Yi）列表，用上述電腦擬合的步驟得到。

混凝土強度與水膠比成反比的公式擬合結(jié)果見圖6。

圖6　混凝土強度與水膠比成反比的數(shù)學(xué)模型

從電腦繪制出的散點圖看，斜直線穿過水膠比散點的較少，兩端的散點在斜直線的以上（低強和高強混凝土的強度估計值較實際偏低），相關(guān)性檢驗值也從0.9419降到0.9149，說明混凝土強度與水膠比近似成反比的公式覆蓋率相對較低，與直線的相關(guān)性沒有與膠水比成正比的高。為驗證其用性，試用3.2C30泵送混凝土用式（17）估算的 W/ B＝0.46比較，式（30）估算 C30水膠比：W/B＝(103.05-38.22)/134.72＝0.48，用式（17）和式（30）估算 C30水膠比，顯然兩者的水膠比估算（W/B）相差較大，表明：

（1）混凝土強度與膠（灰）水比近似成正比的一元線性回歸B公式比用強度與水膠（灰）比近成反比的一元線性回歸更貼近實際。保羅米之所以認(rèn)為混凝土強度與膠（灰）水比近似成正比，而沒用混凝土強度與水膠（灰）比近似成反比是經(jīng)過慎重驗證的，不是別出心裁。

（2）強度—水膠（灰）比近似成反比的一元線性回歸方程，相關(guān)性有少許差別，但仍不失為相關(guān)性特別顯著。因此，可稱為混凝土強度與灰（膠）水比近似成正比，也可稱之為混凝土強度與水灰（膠）比近似成反比，但在配合比設(shè)計應(yīng)用中，混凝土強度與仍以水膠比成正比更貼近實際。另外，從歷史的源頭回顧這一結(jié)論，深感前人對科學(xué)的嚴(yán)緊態(tài)度很值得尊敬、學(xué)習(xí)。

3.6B 公式仍可用于現(xiàn)代混凝土

A 定則—B 公式提出的100多年來，是被混凝土實踐證明了的符合生產(chǎn)實際的，并為已往混凝土的生產(chǎn)、技術(shù)發(fā)展和提高帶來應(yīng)有的理論依據(jù)。但對普遍使用高效減水劑和大摻量礦物細(xì)粉料的現(xiàn)代混凝土，原來的水灰比已改為水膠比的今天，它還能用于指導(dǎo)現(xiàn)代混凝土研究和配合比設(shè)計的嗎？這是值得深入探討的問題。為此，我們以 B 公式（26）估算各強度等級混凝土的水膠比與現(xiàn)代普通商品混凝土實際生產(chǎn)常用的水膠比作比較，試看表6。

表6　B 公式 ﹙26) 估算水膠比與實際生產(chǎn)的水膠比的對比

用 B 公式估算的普通商品混凝土各強度等級的 W/B，在C20～C55（W/B＝0.57～0.31）范圍內(nèi)與王元[8]對遼寧省百余家商品混凝土企業(yè)調(diào)查統(tǒng)計后，用1000多組水膠比與強度數(shù)據(jù)統(tǒng)計的結(jié)果 十分接近。也與筆者在近二十年中，從事商品混凝土技術(shù)實踐的經(jīng)驗數(shù)據(jù)基本一致。說明：

（1）B 公式雖誕生于170余年前的干硬性、塑性的四組分混凝土，對于今天的流動性、大流動性的多組分混凝土的配合比設(shè)計與生產(chǎn)技術(shù)管理，仍有現(xiàn)實的理論指導(dǎo)意義，可直接用于現(xiàn)代普通商品混凝土配合比設(shè)計，不代入水泥強度及粉煤灰礦粉影響系數(shù)，亦可求出所要求配制強度的水膠比估計值；

（2）混凝土強度在 C20～C50范圍內(nèi)估算值與試驗值基本－致，只是在 C10～C15和＞C50時，B 公式計算的W/B和 R28略有偏差，但均未超過5%。根據(jù)這一特點，在試驗時對其進行適當(dāng)?shù)恼{(diào)整即可。另外，B公式對現(xiàn)代多組分混凝土應(yīng)用范圍不是縮小，而是進一步擴大，不再是早年認(rèn)為的（C20～C40）B/W＝1.5～2.5，而是1.43～3.45。說明 A 定則—B 公式隨時代的科技發(fā)展不是失去實用意義，而是進一步擴大適用范圍，對現(xiàn)代混凝土發(fā)展仍有較強的理論指導(dǎo)意義。究其原理是：保羅米是把灰水比—混凝土強度的關(guān)系，近似簡化成直線關(guān)系的數(shù)學(xué)摸型，建立在大量試驗建立的線型數(shù)模和對不確定性數(shù)據(jù)運用最小二乘法數(shù)學(xué)原理基礎(chǔ)上的，在理論上是嚴(yán)謹(jǐn)?shù)摹?/p>

（3）筆者擬合的 B 公式計算出的 C20～C50混凝土 R28和 W/B 估計值，之所以與遼寧省建科院王元根據(jù)1000多組調(diào)查得出的統(tǒng)計值殊途同歸，完全一致絕不是偶然的巧合，它是反映當(dāng)前我國南北方商品混凝土的整體水平，具有一定的客觀現(xiàn)實性，不以人的主觀意志為轉(zhuǎn)移。這也說明 B 公式有其應(yīng)用上的普適性。

（4）對最常用的 C25、C30的引氣型混凝土用式（29）估算 W/B 與表中 ACI 推薦的 W/B 十分接近（特別是最常用的 C25、C30），差別主要在于樣本取于不同的國家和年代以及 B 公式與 A 定則本來就是近似關(guān)系有關(guān)。

4　小結(jié)

（1）B 公式是對 A 定則的簡化，兩者沒有質(zhì)的區(qū)別。時代的發(fā)展、混凝土的科技進步、高效減水劑和礦物摻合料等的應(yīng)用，都可在其樣本子樣（Xi，Yi）中得以反映，改變的是兩回歸系數(shù)，而 A 定則—B 公式的數(shù)學(xué)表達形式?jīng)]有變化，仍可用于現(xiàn)代混凝土的配合比設(shè)計，有很強的現(xiàn)實指導(dǎo)意義；

（2）A 定則—B 公式是人類對混凝土認(rèn)識從感性到理性進步的重要標(biāo)志，他們的貢獻應(yīng)給予充份肯定。但隨實踐、認(rèn)識、再實踐、再認(rèn)識的不斷循環(huán)往復(fù)與發(fā)展，人們對混凝土強度的認(rèn)識總是希望有所發(fā)現(xiàn)發(fā)展，認(rèn)識永遠(yuǎn)不會停留在一個層面上，A 定則—B 公式出現(xiàn)新的演變是必然的。

[8] 王元．減水劑的應(yīng)用對使用水灰比定則的影響[J]．商品混凝土，2013(3):1-6．

［通訊地址］上海松江區(qū)榮樂中路228弄北九峰小區(qū)152號302室（201600）

王永逵（1934—），前蘭州鐵道學(xué)院土木系建材室主任。