【導讀】熔體中組作用的原理。通過收集有關活度計算用文獻，調(diào)研計算用原始數(shù)據(jù)，根。據(jù)共存理論建立CaO-SiO2、MnO-SiO2、PbO-SiO2、Na2O-SiO2、MgO-SiO2、一個組元作用濃度，并對預測結(jié)果與實驗測定值進行了對比。結(jié)果發(fā)現(xiàn)，模型計。論建立健的模型較為正確。

　　

【正文】 212421 42 , nnKnnnnK ?? () TG 62???? ? J/mol 質(zhì)量平衡為： 14321 ???? nnnn () ? ?? ??? 432 nnnxA () ? ?? ??? 431 nnnxB () 通過上述方程整理得： ? ? ? ? 4321 ?????? nbanbabnan () 通過方程式（ ）和式（ ）聯(lián)立可建立模型。 MnOTiO2渣系計算結(jié)果及結(jié)論 在 1823K下，反應 （ ）、（ ） 的平衡常數(shù) K K2由對應的 ?G? =RTlnK計算得出，計算結(jié)果列入表 。采用模型（ ）、（ ）求解模型不同 A、B值下 MnO作用濃度 nMnO，采用 Matlab編程求解，其程序附錄 Ⅱ 相似 ，計算結(jié)果列入表 。對模型預測的 nMnO與對應 MnO實驗測定 的活度進行比較，結(jié)果見圖 ，計算的相對誤差、偏差結(jié)果列入表 。 計算發(fā)現(xiàn)， 在 1823K下 ， 通過建立的模型預測的 MnOSiO2二元體系中 MnO作用濃度 與實測活度值之間相對誤差為 %，偏差為 ，結(jié)果表明了基于熔渣共存理論建立 MnOSiO2二元體系中 MnO作用濃度預測較為合理 。 攀枝花學院本科畢業(yè)設計（論文） 2 基于共存理論的熔渣 中 組元活度的預測 23 表 1823K下 MnOSiO2熔體組元 MnO作用濃度 預測用數(shù)據(jù) 平衡常數(shù)代號 平衡常數(shù)值 K1 K2 溫度 /K 1823 表 1823K下 MnOSiO2熔體組元 MnO作用濃度 預測結(jié)果 MnO 的摩爾濃度 MnO 的實驗測定活度 MnO 的預測定的作用濃度 MnO 作用濃度計算值與活度實驗測定值之間的相對誤差 /% MnO 作用濃度計算值與活度實驗測定值之間的偏差 01 x M n OaMnO/NMnO 圖 在 1823K下 MnOTiO2熔體中組元的預測作用濃度與實驗測定的活度對比 (MnO計算的作用濃度 N（ — ）、實測活度 a（ ● ） ) 攀枝花學院本科畢業(yè)設計（論文） 結(jié)論 24 結(jié)論 采用共存理論，預測了 CaOSiO2渣系、 MnOSiO2渣系、 PbOSiO2渣系、Na2OSiO2渣系 、 MgOSiO2渣系 、 MnOTiO2渣系 六個二元熔渣體系中組元的活度，同時應用 Matlab編程求解了模型中關鍵參數(shù)。計算得到： （ 1）由爐渣結(jié)構(gòu)共存理論和 CaOSiO2渣系的相圖從而制定了在 1500℃ 溫度的作用濃度計算模型，由模型預測的 CaO作用濃度與實驗測定的 CaO活度之間相對誤差為 %，偏差為 。 （ 2） 由爐渣結(jié)構(gòu)共存理論和 MnOSiO2渣系的相圖從而制定了在 1873K溫度的作用濃度計算模型，由模型預測的 MnO作用濃度與實驗測定的 MnO活度之間相對誤差為 %，偏差為 。 （ 3） 由爐渣結(jié)構(gòu)共存理論和 PbOSiO2渣系的相圖從而制定了在 1273K溫度的作用濃度計算模型，由模型預測的 PbO作用濃度與實驗測定的 PbO活度之間相對誤差為 %，偏差為 。 （ 4） 由爐渣結(jié)構(gòu)共存理論和 Na2OSiO2渣系的相圖從而制定了在 1373K溫度的作用濃度計算模型，由模型預測的 Na2O作用濃度與實驗測定的 Na2O活度之間相對誤差為 %，偏差為 。 （ 5） 由爐渣結(jié)構(gòu)共存理論和 MgOSiO2渣系的相圖從而制定了在 1873K溫度的 作用濃度計算模型，由模型預測的 MgO作用濃度與實驗測定的 MgO活度之間相對誤差為 %，偏差為 。 （ 6） 由爐渣結(jié)構(gòu)共存理論和 MnOTiO2渣系的相圖從而制定了在 1823K溫度的作用濃度計算模型，由模型預測的 MnO作用濃度與實驗測定的 MnO活度之間相對誤差為 %，偏差為 。 除 Na2OSiO2渣 系外，基于共存模型 計算所得的 5組熔渣中一組元的 作用濃度與 對應 實測活度對比相差并不大， 表明基于共存理論建立的爐渣 模型 較為 正確 。攀枝花學院本科畢業(yè)設計（論文） 參考文獻 25 參考文獻 [1]張鑒 . 冶金熔體和溶液的計 算熱力學 [M]. 北京 :冶金工業(yè)出版社 , 20xx: 297~298. [2]文永才 , 趙海雷 , 公茂秀等 . 煉鋼用錳球團礦直接合金化動力學 [J]. 北京科技大學學報1999,10: 24~28. [3]黃希祜 . 鋼鐵冶金原理 [Ｍ ].北京 :冶金工業(yè)出版社 , 20xx. [4]德國鋼鐵工程師協(xié)會 , 王儉譯 . 渣圖集 [Ｍ ]. 北京 :冶金工業(yè)出版社 , 1989. [5]鄂加強 , 王春華 , 龔金科等 . 銅火法冶煉熱動力學系統(tǒng)實測數(shù)據(jù) EMD處理 [J]. 中國有色金屬學報 , 20xx, 18(5): 946~951. [6]汪金良 , 盧宏 , 曾青云 , 張傳福 . 基于遺傳算法的銅閃速熔煉過程控制優(yōu)化 [J]. 中國有色金屬學報 , 20xx,17(1): 156~160. [7] NICHITA D V, GOMEZ S, LUNA E. Multiphase equilibria calculation by direct minimization of Gibbs free energy with a global optimization method[J]. Computers and Chemical Engineering, 20xx, 26(12): 1703~1724. [8]ILYUSHECHKIN A, HAYES P C, JAK E. Liquids temperatures in calcium ferrite slag in equilibrium with molten copper[J]. Metallurgical and Materials Transactions B, 20xx, 35(2): 203~215. [9]賈娟魚 .含鈦爐渣中 TiO2作用濃度及爐渣粘度研究 [J].重慶大學， 20xx. [10]梁英教 , 車蔭昌 . 無機物熱力學數(shù)據(jù)手冊 [M]. 沈陽 :東北大學出版社 , 1993: 458. [11] SANG H L, SEOK M M, DONG J M, JOD H P. Thermodynamic behavior of nickel in CaOSiO2FetO slag[J]. Metallurgical and Materials Transactions B, 20xx, 33(1): 55~59. [12]李正邦 . 礦物直接合金化冶煉合金鋼－理論與實踐 [Ｍ ]. 北京 :冶金工業(yè)出版社 , 20xx. [13]文永才 , 趙海雷 ,公茂秀等 . 煉鋼用錳球團礦直接合金化動力學 [J]. 北京科技大學學報1990,1:24~28. [14] 周蘭花 . MIVM參數(shù)的算法研究 [J]. 材料科學與工程學報 [J].20xx,3: 45~49.攀枝花學院本科畢業(yè)設計（論文） 致謝 26 致謝 經(jīng)過一個學期的努力終于有了部分成果，本次畢業(yè)設計已經(jīng)接近尾聲，一份沉甸甸的畢業(yè)論文終于呈現(xiàn)在了我的面前。作為一個本科生的畢業(yè)設計，由于經(jīng)驗的匱乏，難免有許多考慮不周全的地方，如果沒有導師的督促指導以及一起學習的同學們的支持，想要完成本次設 計是十分困難的。 在這里首先要感謝導師周蘭花老師。周老師平日里工作繁多，但在我做畢業(yè)設計的每個階段，她都兢兢業(yè)業(yè)，為我們的論文做仔細的指導，特別是對于我課題中的軟件使用，給予了細心地指導。同時周老師還在思想、生活上也給我們無微不至的關懷，在此謹向周老師致以誠摯的謝意和崇高的敬意，這些美好的瞬間將永遠記載在我的心中。 大學四年一晃而過，但是這其中的美好回憶是一生的，在此感謝大學四年來所有的老師，為我們打下良好的冶金專業(yè)知識基礎，這些知識部分運用到了我的畢業(yè)論文當中特別是冶金原理的運用。其次要感謝周老 師課題組其余同學，他們在我論文的寫作過程中給予了我極大的幫助。 最后感謝我的母校 —攀枝花學院四年來對我的大力栽培，四年的學習生活過程讓我學會了如何做人，如何做一個對社會有用的人。 攀枝花學院本科畢業(yè)設計（論文） 附錄 I 27 附錄 Ⅰ 預測 CaOSiO2中 CaO 作用濃度用 MATLAB計算程序： %excel 表中導入 x,a1_exp,a2_exp。 %a1_exp 為組元 CaO 實驗活度值 。 %K1 平衡常數(shù) 1。 %K2 平衡常數(shù) 2。 %K3 平衡常數(shù) 3。 %x 為平衡常數(shù) K K K3對應由 excel 輸入值 。 function [y]=cal_x(x,xx,a1_exp) K1=x(1,1)。 K2=x(2,1)。 K3=x(3,1)。 %nn=zeros(9,2)。 %k=1。 zynd=zeros(2,1)。 %for i=:: a=xx(1,1)。 b=1a。 syms x x=sym2poly(K3*(3*ab+x)*x^3+K2*(2*ab+x)*(x+*a)*x^2+K1*(ab+x)*(x+*a)^2*x(bx)*(x+*a)^3)。 y=roots(x)攀枝花學院本科畢業(yè)設計（論文） 附錄 II 28 附錄 Ⅱ 預測 PbOSiO2中 PbO 作 用濃度用 Matlab 計算程序 ： %excel表中導入 x,a1_exp,a2_exp。 %a1_exp 為組元 PbO 實驗活度值 。 %a2_exp 為組元 PbO 實驗活度值 。 %K1 平衡常數(shù) 1。 %K2 平衡常數(shù) 2。 %K3 平衡常數(shù) 3。 %xdwc3 為組元 PbO活度的計算值與實驗值之間的相對誤差的計算結(jié)果 。 %xdwc4 為組元 SiO2活度的計算值與實驗值之間的相對誤差的計算結(jié)果 。 %bc1 為組元 PbO 活度的計算值與實驗值之間的偏差的計算結(jié)果 。 %bc2 為組元 SiO2活度的計算值與實驗值之間的偏差的計算結(jié)果 。 %nn 為組 元 PbO、 SiO2的作用濃度計算結(jié)果 。 %x 為平衡常數(shù) K K K3對應由 excel輸入值 。 function [xdwc3,xdwc4,bc1,bc2,nn]=cal_N1_N2(x,a1_exp,a2_exp) K1=x(1,1)。 K2=x(2,1)。 K3=x(3,1)。 nn=zeros(9,2)。 k=1。 zynd=zeros(2,1)。 for i=:: a=1i。 b=i。 bc=10000。 for t=::1 for tt=::1 N=fsolve( @(N)[*a*N(1)b*N(2)+(ab)*K1*N(1)*N(2)+(2*ab)*K2*N(1)^2*N(2)+(4*ab)*K3*N(1)^4*N(2)。 1(*a+1)*N(1)(1b)*N(2)(ab+1)*K1*N(1)*N(2)(2*ab+1)*K2*N(1)^2*N(2)(4*ab+1)*K3*N(1)^4*N(2)],[t,tt]39。,optimset(39。Display39。,39。off39。))。 %N 作 用濃度優(yōu)化求解 if isreal(N) xx(1,1)=N(1,1)。 攀枝花學院本科畢業(yè)設計（論文） 附錄 II 29 xx(2,1)=N(2,1)。 %以下以與實驗活度作對比進行偏差達最小為目標繼續(xù)優(yōu)化求 N bc_temp1=abs(xx(1,1)a1_exp(k,1))。 bc_temp2=abs(xx(2,1)a2_exp(k,1))。 bc_temp=(bc_temp1+bc