[發(fā)明專利]一種玻璃光纖拉絲設備優(yōu)化生產(chǎn)方法

有效

• 申請?zhí)枺?/td>: 201310560642.2
• 申請日: 2013-11-12
• 公開（公告）號: CN103613272A
• 公開（公告）日: 2014-03-05
• 發(fā)明（設計）人: 王春林
• 申請（專利權(quán)）人: 杭州電子科技大學
• 主分類號: C03B37/07
• 分類號: C03B37/07
• 代理公司: 杭州求是專利事務所有限公司 33200
• 代理人: 杜軍

1.一種玻璃光纖拉絲機優(yōu)化生產(chǎn)方法，其特征在于該方法的步驟包括：

步驟(1)針對給定的玻璃光纖拉絲設備及光纖直徑的要求，設定分為N級拉伸最終達到光纖產(chǎn)品的直徑要求,每一級的拉出光纖直徑為d_i∈(d_imin,d_imax),i=1,2…N,N≥1，其中d_imin為第i級拉出光纖直徑允許最小值,d_imax為第i級拉出光纖直徑允許最大值；

根據(jù)生產(chǎn)工藝要求設定產(chǎn)品光纖拉絲速度允許的范圍：V∈(V_min,V_max)，其中V為光纖拉絲速度，V_min為允許光纖拉絲最小速度,V_max為允許光纖拉絲最大速度；設定加熱爐加熱溫度的允許范圍T∈(T_min,T_max)及預制棒的饋送速度允許范圍V_y∈(V_ymin,V_ymax)，其中T為加熱爐加熱溫度，T_min為允許加熱爐最低加熱溫度,T_max為允許加熱爐最高加熱溫度，V_y為光纖預制棒饋送速度，V_ymin為允許光纖預制棒最小饋送速度,V_ymax為允許光纖預制棒最大饋送速度；以上所述：光纖拉伸級數(shù)(N),每級對應的光纖拉出直徑(di)，光纖拉絲速度(V)，加熱爐加熱溫度(T)及光纖預制棒饋送速度(Vy)為玻璃光纖生產(chǎn)操作參數(shù)；在參數(shù)允許速度范圍內(nèi)，將各參數(shù)的允許范圍均勻分為M個區(qū)間,M≥1，各參數(shù)允許范圍內(nèi)每個區(qū)間取一個數(shù)據(jù)作為代表，做不同拉絲速度、不同加熱爐的加熱溫度和不同預制棒的饋送速度的組合試驗，并記錄相應參數(shù)組合下光纖產(chǎn)品品質(zhì)，所述光纖產(chǎn)品品質(zhì)指標包括：芯包比(α)、同心率(β)、折射率分布(γ)、直徑比(δ)；將以上試驗數(shù)據(jù)采集入試驗數(shù)據(jù)庫中備用；

步驟(2)采集玻璃光纖拉絲設備生產(chǎn)操作參數(shù)及相關的表征光纖產(chǎn)品品質(zhì)的指標，建立生產(chǎn)實時數(shù)據(jù)庫；具體的玻璃光纖拉絲設備生產(chǎn)操作參數(shù)和產(chǎn)品品質(zhì)指標，通過玻璃光纖拉絲設備運行數(shù)控系統(tǒng)獲取，或直接通過儀器測量采集，將以上生產(chǎn)實時數(shù)據(jù)采集入實時數(shù)據(jù)庫中備用；

步驟(3)建模數(shù)據(jù)選擇，針對給定的玻璃光纖拉絲設備，選取試驗數(shù)據(jù)庫中全部數(shù)據(jù)，另外，再在實時數(shù)據(jù)庫中選取試驗數(shù)據(jù)庫中數(shù)據(jù)量的2倍到5倍的數(shù)據(jù)作為建模數(shù)據(jù)，數(shù)據(jù)要包括步驟(1)和步驟(2)中所述的生產(chǎn)操作參數(shù)及光纖產(chǎn)品品質(zhì)指標的情況；

步驟(4)采用多目標支持向量機方法針對建模數(shù)據(jù)建模，建立玻璃光纖的產(chǎn)品品質(zhì)指標與生產(chǎn)操作參數(shù)間的模型；具體建模方法如下：

用于建模樣本的輸入?yún)?shù)及輸出參數(shù)表示為其中x_i表示第i組作為輸入數(shù)據(jù)向量的玻璃光纖拉絲設備生產(chǎn)操作參數(shù)，y_i表示第i組作為輸出參數(shù)向量的玻璃光纖產(chǎn)品品質(zhì)的指標，N為樣本數(shù)量，以建模數(shù)據(jù)為基礎，建立各玻璃光纖拉絲設備的產(chǎn)品品質(zhì)指標與各生產(chǎn)操作參數(shù)間的關系模型；

采用支持向量機方法建模，支持向量機核函數(shù)選為徑向基函數(shù)

K(xi,xj)=φ(xi)·φ(xj)=exp|(||xi-xj||22σ2)|]]>

其中，參數(shù)“σ”為徑向基函數(shù)的寬度，φ(x)為映射函數(shù)，設所求的目標函數(shù)為：f(x_i)＝w·φ(x_i)+b，f(x_i)為模型輸出的向量值，w為權(quán)重系數(shù)向量，b為截距向量；引入松弛因子ξ^*_i≥0和ξ_i≥0和允許擬合誤差ε，模型可以通過在約束：

yi-w·φ(xi)-b≤ϵ+ξiw·φ(xi)+b-yi≤ϵ+ξi*ξi≥0ξi*≥0i=1,...,N,]]>條件下，最小化：

minR(w,ξ,ξ*)=12w·w+CΣi=1kξ+ξ*]]>

獲得，其中常數(shù)C>0為懲罰系數(shù)，k為預測超出允許擬合誤差點的數(shù)量0≤k＜N；該最小化問題為一個凸二次規(guī)劃問題，引入拉格朗日函數(shù)：

L(w,b,ξ,ξ*,α,α*,γ,γ*)=12w·w+CΣi=1N(ξ+ξ*)-Σi=1Nαi[yi-(ξi+ϵ+f(xi))]-Σi=1Nαi*[ξi*+ϵ+f(xi)-yi]-Σi=1N(γiξi+γi*ξi*)]]>

其中：為拉格朗日乘數(shù)；

在鞍點處，函數(shù)L是關于w,b,ξ_i,ξ_i^*的極小點，也是極大點，最小化問題轉(zhuǎn)化為求其對偶問題的最大化問題；

拉格朗日函數(shù)L在鞍點處是關于w,b,ξ_i,ξ_i^*極小點，得：

∂∂wL=0→w=Σi=1N(αi-αi*)φ(xi)∂∂bL=0→Σi=1N(αi-αi*)=0∂∂ξiL=0→C-αi-γi=0∂∂ξi*L=0→C-αi*-γi*=0]]>

可得拉格朗日函數(shù)的對偶函數(shù)：

此時，

w=Σi=1N(αi-αi*)φ(xi)]]>

f(x)=Σi=1N(αi-αi*)K(x,xi)+b]]>

按照庫恩-塔克(KKT)條件定理，在鞍點有下式成立：

αi[ϵ+ξi-yi+f(xi)]=0αi*[ϵ+ξi+yi-f(xi)]=0i=1,...,N]]>

由上式可見，α_i和都不會同時為非零，可得：

ξiγi=0ξi*γi*=0i=1,...,N]]>

從上式可求出向量b，獲得模型；

步驟(5)利用粒子群優(yōu)化算法，結(jié)合以上步驟(4)所建各玻璃光纖拉絲設備的產(chǎn)品品質(zhì)指標與各生產(chǎn)操作參數(shù)間的關系模型，進行各生產(chǎn)任務的優(yōu)化，具體步驟如下：

a.定義粒子群算法初始群體x向量的各維分量，分別為各玻璃光纖拉絲設備的生產(chǎn)操作參數(shù)；

b.設定粒子群算法的搜索目標和迭代次數(shù)，搜索目標為各玻璃光纖產(chǎn)品品質(zhì)指標的設定值；

c.初始化群體向量x，其中各生產(chǎn)操作參數(shù)的優(yōu)化范圍均為相應允許范圍；結(jié)合步驟(4)所建模型，根據(jù)初始群體中各生產(chǎn)操作參數(shù)的初始化值，獲得各光纖產(chǎn)品品質(zhì)指標的優(yōu)化目標初始值，然后根據(jù)上一步設定的搜索目標，用粒子群算法進行迭代計算，搜索最優(yōu)生產(chǎn)生產(chǎn)操作參數(shù)配置情況；

d.當粒子群算法完成迭代次數(shù)或找到設定要求的最優(yōu)時，停止計算獲得相應最優(yōu)的群體向量，即獲得最優(yōu)的各玻璃光纖生產(chǎn)設備的生產(chǎn)操作參數(shù)及最優(yōu)的光纖產(chǎn)品品質(zhì)指標；依據(jù)該優(yōu)化結(jié)果可以指導實際生產(chǎn)，獲得最優(yōu)的產(chǎn)品品質(zhì)。