混料設計的均勻設計條件下的組織試驗研究

混料設計的均勻設計條件下的組織試驗研究

這種設計在化工、材料工業(yè)、食品、低溫等領域非常重要。為了提高各種用途，應添加抗氧劑、分散劑、粘度指數(shù)改進劑、分凝劑等不同類型的添加劑。為了提高石油和天然氣的用量，應添加抗氧劑、分散劑、粘度指數(shù)、抗凝劑等不同類型的添加劑。添加量和基本油的比例關系決定了油罐的優(yōu)缺點。在實驗設計中,各種添加劑和基礎油稱為配方的成分(ingredient),決定諸成分比例的設計方法稱為配方設計,混料設計(designofexperimentswithmixtures)是其中的一種。設在一個油品配方中有s個成分X1,X2,X3,…,Xs等,則混料設計是區(qū)間Ts={(X1,X2,X3,…,Xs)=Xy≥0,j=1,2,…,s,X1+X2+X3+…+Xs=1}上的一個點,油品配方基于均勻設計思想條件下的混料設計是在Ts上選擇具有代表性的幾個點,目的是減少工作量,降低實驗費用,使有限的數(shù)據(jù)可以很好地擬合實際情況,并通過統(tǒng)計建模,找到一個最好的油品配方。1有約束條件的混料設計算法由于油品成分之間有各種成分之和等于1的約束條件,因此其混料設計大大難于以前的s個無約束條件下的混料設計,下面介紹無約束條件的混料設計和有約束條件的混料設計的算法。1.1ts上的n個點王元和方開泰提出了混料均勻設計,其思想是將n個實驗點(n種不同的實驗配方)均勻地分配在Ts上,采用遞變換的方法,若在單位立方體Cs-1=s-1上有一個均勻設計{Ck=(Ck,1,Ck,2,…,Ck,s-1),k=1,2,…,n},然后變化成為Ts上的n個點,理論上可以闡明這n個點在Ts上是散布均勻的,具體的算法如下:(1)設U(Ns-1)為均勻設計表,其中一個元素記為Ui,j;(2)令Ck,i=(Uk,i-0.5)/n,i=1,2,…,s-1;k=1,2,…,n。則{Ck=(Ck,1,Ck,2,…,Ck,s-1),k=1,2,…,n}為Ts上的一個均勻散布點集;(3)計算:Xk?i=(1?C1/(s?i)k,j)∏j=1i?1C1/(s?i)k,j,k=1,2,?,s?1Xk?s=∏j=1s?1C1/(s?i)k,j,k=1,2???nXk?i=(1-Ck,j1/(s-i))∏j=1i-1Ck,j1/(s-i),k=1,2,?,s-1Xk?s=∏j=1s-1Ck,j1/(s-i),k=1,2???n1.2方開泰和楊振海采用條件分布法ts在油品配方中,通常對成分都有一定的限制,有限制的混料設計可以表示為以下的區(qū)域:T(a,b)={(X1,X2,X3,…,Xs):0≤ai≤Xi≤bj≤1,j=1,2,…,s,X1+X2+X3+…+Xs=1},這里邊界a=(a1,a2,…,as),b=(b1,b2,…,bs),給定邊界的合理條件為,約束條件0≤ai≤Xi≤bj≤1。在設計之前要對邊界條件進行修改:a1=max(a1,b1+1-b),b1=min(b1,a1+1-a),用上式修改過的(a1,b1)來定義Ts(a,b)比較準確,如何在Ts(a,b)上產(chǎn)生均勻設計有不同的方法,王元和方開泰采用變換的方法,但該方法的缺點是不能精確地給出要求的實驗數(shù)n,于是方開泰和楊振海改用條件分布法在Ts(a,b)上來生成均勻設計,其算法如下:令b*=(b-a)/(1-a),y=(y1,y2,…,ys)=(x-a)/(1-a),ds=max(0,1+bs-b)根據(jù)ys的實驗分布:Fs=???????(1?ds)s?1?(1?y)s?1(1?ds)s?1?(1?bs)s?110ds≤y≤bsy≥bsy≤dsFs={(1-ds)s-1-(1-y)s-1(1-ds)s-1-(1-bs)s-1ds≤y≤bs1y≥bs0y≤ds運用逆變化可知:ys=1-[u(1-bs)s-1+(1-u)(1-ds)s-1]1/(s-1),ds≤y≤bs式中:u為均勻分布u(0,1)的一個隨機樣本。(1)令Δs=1?Δk=1?∑i=k+1syi?k=s?1,?,2dk=max{ak/Δk,1?∑i=1k?1bi/Δk},k=s,s?1,?,2?k=max{bk/Δk,1?∑i=1k?1ai/Δk},k=s,s?1,?,2(1)令Δs=1?Δk=1-∑i=k+1syi?k=s-1,?,2dk=max{ak/Δk,1-∑i=1k-1bi/Δk},k=s,s-1,?,2?k=max{bk/Δk,1-∑i=1k-1ai/Δk},k=s,s-1,?,2(2)在Ts(a,b)上產(chǎn)生均勻分布的算法:①產(chǎn)生s-1個相互獨立的隨機數(shù)us…u2;②xk=G(uk,ds,?k,Δk,k-1),k=s,s-1,…,2;③x1=(x1,…,xs)為Ts(a,b)上的一個樣本。(3)在以上的基礎上不難給出有混料均勻設計的算法:①在Cs-1上產(chǎn)生一個均勻設計un(ns-1),記為U=(Ui,j);②計算xi,j=G(uk,ds,?k,Δk,k-1),k=s,s?1,?,2xi,1=1?∑k=2sxkk=s,s-1,?,2xi,1=1-∑k=2sxk則{x1=(x1,x2,…,xs),i=1,2,…,n}為一個在限制Ts(a,b)上的混料設計。2配方設計的限制本文中的例子是研究可生物降解潤滑劑的抗氧性能?？缮锝到獾臐櫥瑒μ砑觿┯兄車栏竦囊?例如,德國的“藍色天使”組織要求可生物降解潤滑劑的不可生物降解添加劑的加量在7%以下。因此,油品的配方設計是一個有限制的混料設計,上文第二種方法理論要求高、計算量大,目前在沒有專業(yè)軟件的情況下難以推廣。因此,本文根據(jù)油品配方的具體情況,采用了第一種算法,把各種不同類型添加劑看為一個整體,限量固定值,應用無限制混料條件來確定配方,其它的部分用基礎油來填充,在本文中應用清凈劑復劑(T102,T106,T113)加量4%,分散劑復劑(T151,T152,L117,T154)加量5%,來確定復劑中的單劑比例。2.1i的計算選用U6(62)正交表(表1)。計算:Ck,i=(Uk,i-0.5)/n,i=1,2;k=1,2,…,6,得出表2所示的結果。應用Xk?i=(1?C1/(s?i)k,j)∏j=1i?1C1/(s?i)k,j(k=1,2???6Xk?i=(1-Ck,j1/(s-i))∏j=1i-1Ck,j1/(s-i)(k=1,2???6;i=1,2)和Xk?s=∏j=1s?1C1/(s?i)k,j(k=1,2,?,6Xk?s=∏j=1s-1Ck,j1/(s-i)(k=1,2,?,6;s=3)得出表3所示的結果。2.2ii,j1,2,3e選用U8(83)正交表(表4)。計算:Ck,i=(Uk,i-0.5)/n,i=1,2,3;k=1,2,…,8,得出表5所示的結果。應用Xk?i=(1?C1/(s?i)k,j)∏j=1i?1C1/(s?i)k,j(k=1,2,?,6Xk?i=(1-Ck,j1/(s-i))∏j=1i-1Ck,j1/(s-i)(k=1,2,?,6;i=1,2,3)和Xk?s=∏j=1s?1C1/(s?i)k,j(k=1?2???8Xk?s=∏j=1s-1Ck,j1/(s-i)(k=1?2???8;s=4),得出表6所示的結果。2.3實驗加量結果按表3來組織實驗,清凈劑復劑(T102,T106,T113)總加量為4%,實驗結果如表7所示。按表6來組織實驗,分散劑復劑(T151,T152,L117,T154)總加量為5%,實驗結果如表8所示。2.4觀法其為最佳方案表7中,清凈劑復劑(T102,T106,T113)在方案6表現(xiàn)最好,而且根據(jù)直觀法其也是最佳方案。本實驗證明了該實驗設計覆蓋了最佳方案。表8中,分散劑復劑(T151,T152,L117,T154)在方案4中表現(xiàn)最好(甚至比理論最佳方案好),本實驗證明了該實驗設計覆蓋了最佳方案,以上的兩個實驗的復劑表現(xiàn)都比單劑強。3混料設計的實驗進展(1)應用混