基于共存理論的熔渣組元活度的預(yù)測(cè)畢業(yè)論文-資料下載頁(yè)
2025-06-27 19:23本頁(yè)面
【正文】 參考文獻(xiàn)參考文獻(xiàn)[1]張鑒. 冶金熔體和溶液的計(jì)算熱力學(xué)[M]. 北京:冶金工業(yè)出版社, 2007: 297~298.[2]文永才, 趙海雷, 公茂秀等. 煉鋼用錳球團(tuán)礦直接合金化動(dòng)力學(xué)[J]. 北京科技大學(xué)學(xué)報(bào)1999,10: 24~28.[3]黃希祜. 鋼鐵冶金原理[Ｍ].北京:冶金工業(yè)出版社, 2008.[4]德國(guó)鋼鐵工程師協(xié)會(huì), 王儉譯. 渣圖集[Ｍ]. 北京:冶金工業(yè)出版社, 1989.[5]鄂加強(qiáng), 王春華, 龔金科等. 銅火法冶煉熱動(dòng)力學(xué)系統(tǒng)實(shí)測(cè)數(shù)據(jù)EMD處理[J]. 中國(guó)有色金屬學(xué)報(bào), 2008, 18(5): 946~951.[6]汪金良, 盧宏, 曾青云, 張傳福. 基于遺傳算法的銅閃速熔煉過(guò)程控制優(yōu)化[J]. 中國(guó)有色金屬學(xué)報(bào), 2007,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, 2002, 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, 2004, 35(2): 203~215.[9][J].重慶大學(xué)，2003.[10]梁英教, 車(chē)蔭昌. 無(wú)機(jī)物熱力學(xué)數(shù)據(jù)手冊(cè)[M]. 沈陽(yáng):東北大學(xué)出版社, 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, 2002, 33(1): 55~59.[12]李正邦. 礦物直接合金化冶煉合金鋼－理論與實(shí)踐[Ｍ]. 北京:冶金工業(yè)出版社, 2009.[13]文永才, 趙海雷,公茂秀等. 煉鋼用錳球團(tuán)礦直接合金化動(dòng)力學(xué)[J]. 北京科技大學(xué)學(xué)報(bào)1990,1:24~28.[14] 周蘭花. MIVM參數(shù)的算法研究[J]. 材料科學(xué)與工程學(xué)報(bào)[J].2011,3: 45~49.26攀枝花學(xué)院本科畢業(yè)設(shè)計(jì)（論文） 致謝致謝經(jīng)過(guò)一個(gè)學(xué)期的努力終于有了部分成果，本次畢業(yè)設(shè)計(jì)已經(jīng)接近尾聲，一份沉甸甸的畢業(yè)論文終于呈現(xiàn)在了我的面前。作為一個(gè)本科生的畢業(yè)設(shè)計(jì)，由于經(jīng)驗(yàn)的匱乏，難免有許多考慮不周全的地方，如果沒(méi)有導(dǎo)師的督促指導(dǎo)以及一起學(xué)習(xí)的同學(xué)們的支持，想要完成本次設(shè)計(jì)是十分困難的。 在這里首先要感謝導(dǎo)師周蘭花老師。周老師平日里工作繁多，但在我做畢業(yè)設(shè)計(jì)的每個(gè)階段，她都兢兢業(yè)業(yè)，為我們的論文做仔細(xì)的指導(dǎo)，特別是對(duì)于我課題中的軟件使用，給予了細(xì)心地指導(dǎo)。同時(shí)周老師還在思想、生活上也給我們無(wú)微不至的關(guān)懷，在此謹(jǐn)向周老師致以誠(chéng)摯的謝意和崇高的敬意，這些美好的瞬間將永遠(yuǎn)記載在我的心中。 大學(xué)四年一晃而過(guò)，但是這其中的美好回憶是一生的，在此感謝大學(xué)四年來(lái)所有的老師，為我們打下良好的冶金專(zhuān)業(yè)知識(shí)基礎(chǔ)，這些知識(shí)部分運(yùn)用到了我的畢業(yè)論文當(dāng)中特別是冶金原理的運(yùn)用。其次要感謝周老師課題組其余同學(xué)，他們?cè)谖艺撐牡膶?xiě)作過(guò)程中給予了我極大的幫助?！ ∽詈蟾兄x我的母?！手▽W(xué)院四年來(lái)對(duì)我的大力栽培，四年的學(xué)習(xí)生活過(guò)程讓我學(xué)會(huì)了如何做人，如何做一個(gè)對(duì)社會(huì)有用的人。30攀枝花學(xué)院本科畢業(yè)設(shè)計(jì)（論文） 附錄I附錄Ⅰ預(yù)測(cè)CaOSiO2中CaO作用濃度用MATLAB計(jì)算程序： %excel表中導(dǎo)入x,a1_exp,a2_exp。%a1_exp為組元CaO實(shí)驗(yàn)活度值。%K1平衡常數(shù)1。%K2平衡常數(shù)2。%K3平衡常數(shù)3。%x為平衡常數(shù)KKK3對(duì)應(yīng)由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 xx=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)攀枝花學(xué)院本科畢業(yè)設(shè)計(jì)（論文） 附錄II附錄Ⅱ預(yù)測(cè)PbOSiO2中PbO作用濃度用Matlab計(jì)算程序：%excel表中導(dǎo)入x,a1_exp,a2_exp。%a1_exp為組元PbO實(shí)驗(yàn)活度值。%a2_exp為組元PbO實(shí)驗(yàn)活度值。%K1平衡常數(shù)1。%K2平衡常數(shù)2。%K3平衡常數(shù)3。%xdwc3為組元PbO活度的計(jì)算值與實(shí)驗(yàn)值之間的相對(duì)誤差的計(jì)算結(jié)果。%xdwc4為組元SiO2活度的計(jì)算值與實(shí)驗(yàn)值之間的相對(duì)誤差的計(jì)算結(jié)果。%bc1為組元PbO活度的計(jì)算值與實(shí)驗(yàn)值之間的偏差的計(jì)算結(jié)果。%bc2為組元SiO2活度的計(jì)算值與實(shí)驗(yàn)值之間的偏差的計(jì)算結(jié)果。%nn為組元PbO、SiO2的作用濃度計(jì)算結(jié)果。%x為平衡常數(shù)KKK3對(duì)應(yīng)由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=::1N=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)。 xx(2,1)=N(2,1)。 %以下以與實(shí)驗(yàn)活度作對(duì)比進(jìn)行偏差達(dá)最小為目標(biāo)繼續(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_temp2)。 if bcbc_temp amp。 N(1,1)0 amp。 N(2,1)0 bc=bc_temp。 zynd=N。 end end end end nn(k,1)=zynd(1,1)。 nn(k,2)=zynd(2,1)。 k=k+1。end xdwc3=0。 xdwc4=0。 bc1=0。 bc2=0。for ii=1:9 xdwc3=xdwc3+abs(nn(ii,1)a1_exp(ii,1))/a1_exp(ii,1)。 xdwc4=xdwc4+abs(nn(ii,2)a2_exp(ii,1))/a2_exp(ii,1)。 bc1=bc1+(nn(ii,1)a1_exp(ii,1))^2。 bc2=bc2+(nn(ii,2)a2_exp(ii,1))^2。end bc1=(bc1/9)^。 xdwc3=xdwc3/9*100。 bc2=(bc2/9)^。 xdwc4=xdwc4/9*100。
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1