孫慶利
(濟南大學工程訓練中心,山東濟南250022)
代數法是根據物料平衡建立多元一次聯立方程,可以分別求出各種原料在生料中的配合比。文獻[1]介紹了簡化代數法的配料計算方法,該方法省略了諸如物料成分的干燥基準向灼燒基準的換算等,是一種較為簡化的方法。這里采用該方法為例,借助Excel輔助計算,向大家展示Excel在代數配料法中的應用。
將已知的原料、煤灰化學成分分析(干燥基準)數據(見表1)及配料要求的其它數據石灰飽和系數KH、硅率SM、鋁率IM和熟料中煤灰摻入量gA(%)帶入下列方程就可以直接求出干燥基準的原料1~3的百分配合比X干、Y干和Z干。

表1 已知的原料、煤灰干燥基準化學成分 %

式中的行列式為:


上述7式中:ΔX、ΔY、ΔZ分別是X、Y、Z的系數的行列式值;ΔXYZ是方程常數項的行列式值。它們都沒有單位。
用前一組方程,即用率值KH、SM可求得X干、Y干、Z干。
行列式中的各個元素為:

用后一組方程,即用率值KH、IM也可求得X干、Y干、Z干,其它均與以上元素相同,只是以下元素不同:

已知條件[2][1]:原料和煤灰的干燥基準化學成分見表1,熟料熱耗q為3762kJ/kg。熟料率值要求:KH=0.89、SM=2.0、IM=1.2。煤灰的沉降率S為100%。試計算各原料百分配合比:石灰石X干、黏土Y干、鐵粉Z干。

表2 原料和煤灰的干燥基準化學成分 %

表3 煤的工業分析 %
(1)首先啟動Excel,將工作簿命名為“Excel在代數配料法中的輔助計算”,保存在指定的磁盤中。
(2)原始數據輸入。
在Excel工作表Sheet1中輸入說明文本、原料和煤灰的化學成分分析數據,建立“原料和煤灰的干燥基準化學成分表”。在該表下方輸入“煤的工業分析表”。在“煤的工業分析表”下方輸入配料計算采用的其它參數,如圖1所示(注:煤工業分析中各個量的單位見表3)。
(3)計算煤灰摻入量gA。
100kg熟料中的煤灰摻入量gA按下式計算:

用鼠標選中工作表Sheet1的單元格D17,然后在公式編輯欄中輸入煤灰摻入量的計算公式“=F14*C10*F13/(100*D10)”。再用鼠標點擊“確定”按鈕或者點擊公式按鈕左邊的“√”,Excel就自動在“D17”中顯示煤灰摻入量的計算結果。
(4)計算行列式中的各個元素。
這里是采用率值KH、SM來求X干、Y干、Z干的方法。求a1:

用鼠標選中工作表Sheet1的單元格B21,然后在公式編輯欄中輸入計算公式“=2.8*B13*B3+1.65* C3+0.35*D3-E3”。再用鼠標點擊“確定”按鈕或者點擊公式按鈕左邊的“√”,Excel就自動在“B21”中顯示a1的計算結果。

圖1 原始數據的輸入
求b1:

用鼠標選中工作表Sheet1的單元格C21,然后在公式編輯欄中輸入計算公式“=2.8*B13*B4+ 1.65*C4+0.35*D4-E4”。再用鼠標點擊“確定”按鈕或者點擊公式按鈕左邊的“√”,Excel就自動在“C21”中顯示b1的計算結果。
求c1:

用鼠標選中工作表Sheet1的單元格D21,然后在公式編輯欄中輸入計算公式“=2.8*B13*B5+ 1.65*C5+0.35*D5-E5”。再用鼠標點擊“確定”按鈕或者點擊公式按鈕左邊的“√”,Excel就自動在“D21”中顯示c1的計算結果。求d1:

用鼠標選中工作表Sheet1的單元格E21,然后在公式編輯欄中輸入計算公式“=D17*[E6-(2.8*B13*B6+1.65*C6+0.35*D6)]”。再用鼠標點擊“確定”按鈕或者點擊公式按鈕左邊的“√”,Excel就自動在“E21”中顯示d1的計算結果。
求a2:

用鼠標選中工作表Sheet1的單元格B22,然后在公式編輯欄中輸入計算公式“=B14*(C3+D3) -B3”。再用鼠標點擊“確定”按鈕或者點擊公式按鈕左邊的“√”,Excel就自動在“B22”中顯示a2的計算結果。
求b2:

用鼠標選中工作表Sheet1的單元格C22,然后在公式編輯欄中輸入計算公式“=B14*(C4+D4)-B4”。再用鼠標點擊“確定”按鈕或者點擊公式按鈕左邊的“√”,Excel就自動在“C22”中顯示b2的計算結果。
求c2:

用鼠標選中工作表Sheet1的單元格D22,然后在公式編輯欄中輸入計算公式“=B14*(C5+D5)-B5”。再用鼠標點擊“確定”按鈕或者點擊公式按鈕左邊的“√”,Excel就自動在“D22”中顯示c2的計算結果。
求d2:

用鼠標選中工作表Sheet1的單元格E22,然后在公式編輯欄中輸入計算公式“=D17*[B6-B14* (C6+D6)]”。再用鼠標點擊“確定”按鈕或者點擊公式按鈕左邊的“√”,Excel就自動在“E22”中顯示d2的計算結果。
求a3:

用鼠標選中工作表Sheet1的單元格B23,然后在公式編輯欄中輸入計算公式“=(100-F3)/100”。再用鼠標點擊“確定”按鈕或者點擊公式按鈕左邊的“√”,Excel就自動在“B23”中顯示a3的計算結果。
求b3:

用鼠標選中工作表Sheet1的單元格C23,然后在公式編輯欄中輸入計算公式“=(100-F4)/100”。再用鼠標點擊“確定”按鈕或者點擊公式按鈕左邊的“√”,Excel就自動在“C23”中顯示b3的計算結果。
求c3:

用鼠標選中工作表Sheet1的單元格D23,然后在公式編輯欄中輸入計算公式“=(100-F5) /100”。再用鼠標點擊“確定”按鈕或者點擊公式按鈕左邊的“√”,Excel就自動在“D23”中顯示c3的計算結果。
求d3:

用鼠標選中工作表Sheet1的單元格E23,然后在公式編輯欄中輸入計算公式“=100-D17”。再用鼠標點擊“確定”按鈕或者點擊公式按鈕左邊的“√”,Excel就自動在“E23”中顯示d3的計算結果。
行列式中各元素所對應的計算值如圖2所示。

圖2 行列式各元素所對應的值
Excel數學和三角函數中提供的MDETERM函數(返回數組的矩陣行列式值)可以用來計算矩陣的行列式值,它的精度可達16位有效數字,因此運算結果可能由于位數的取舍而導致某些微小誤差。
行列式常被用來求解多元聯立方程式組。這正合適我們本文的需要。
以(4)~(7)和(1)~(3)公式為依據,將行列式計算相關的數據匯總到Excel工作表中(見圖3)。

圖3 分子、分母計算和最終結果的計算
在E25單元格中輸入公式“MDETERM(B25:D27)”,再用鼠標點擊“確定”按鈕或者點擊公式按鈕左邊的“√”即可。并將該公式復制(方法是:對準E25單元格,按鼠標右鍵,出現菜單,再按“復制”,然后對準J25等單元格,按鼠標右鍵,出現菜單,再按“粘貼”即出現結果。)到“J25”、“E32”單元格中。此時已經分別求出X干、Y干、Z干的分子。分母的計算是在B35單元格中輸入公式“E25+J25+E32”。再用鼠標點擊“確定”按鈕或者點擊公式按鈕左邊的“√”,Excel就自動在“B35”中顯示分母的計算結果。
在“G32”單元格中輸入公式“100*E25/B35”,在“G33”單元格中輸入公式“100*J25/B35”,在“G34”單元格中輸入公式“100*E32/B35”,再分別用鼠標點擊“確定”按鈕或者點擊公式按鈕左邊的“√”,Excel就自動在“G32”、“G33”、“G34”中最終計算結果:石灰石X干=82.63%,黏土Y干=13.85%,鐵粉Z干=3.52%。如圖3所示。該結果與文獻[1]、[2]的結果基本相同。
有關熟料率值和礦物組成的計算請讀者自行進行,這里不再贅述。
用Excel計算水泥配合料,可以將計算方法、計算結果及原始數據保存在一個或多個工作簿內,這有利于對科研和生產資料進行有效的管理。
對于一個系列的水泥配料計算,不需要改變配料的計算方法,只需要修改水泥的某個率值、原料化學成分中的一個或多個原始數據,在這種情況下,可以應用已編輯好的電子表格計算水泥的新配料。實際上所編制好Excel工作簿相當于一個編制好了的程序,只需要將需要修改的數據輸入進去,而其它的數據不用改變,就可以立即得到新的配料計算結果,既迅速又正確。不必每次都從頭開始。注意:在完成數據輸入后,鼠標要在空白單元格處點擊一下。
為了方便管理,計算時,將該計算表格復制,再在復制的表格中修改數據,這樣可以使我們的工作井然有序,不至于發生混亂。
另外,代數方程的求解還可利用逆矩陣和矩陣乘積的方法。具體到Excel中就是采用數學和三角函數中提供的MINVERSE(返回數組矩陣的逆矩陣)、MMULT(返回矩陣的乘積)函數來求解方程。
[1]孫慶利.簡化代數法和試湊法的配料計算.天津:水泥技術.1992年第4期.
[2]南京化工學院等編.水泥工藝原理.第一版,北京:中國建筑工業出版社,1980年.
[3]蘇華等編著.Excel 2003函數在辦公中的應用.第一版,北京:人民郵電出版社.2007年.