一種精軋帶鋼局部高點的控制方法

            文檔序號:3144081閱讀:217來源:國知局

            專利名稱::一種精軋帶鋼局部高點的控制方法
            技術領域
            :本發明涉及一種熱軋帶鋼成型加工
            技術領域
            ,尤其涉及一種精軋帶鋼局部高點的控制方法。
            背景技術
            :熱軋帶鋼橫斷面的局部高點是指該點處的帶鋼厚度比相鄰區域厚,它是板形質量的重要指標。產生局部高點的因素有多個方面,但關鍵是由于精軋后段機架工作輥的不均勻磨損及帶鋼凸度偏小造成。因此通過降低后段機架的不均勻磨損、采用邊部潤滑,提高帶鋼凸度是解決局部高點的有效途徑。軋輥橫移是實現軋輥磨損均勻化的重要技術。它可以均勻軋輥磨損、分散軋輥熱凸度、擴大帶鋼凸度控制范圍,減少局部高點的產生。它是進行逆寬軋制、增加同寬軋制長度,實現自由程序軋制的重要手段之一。周期橫移方式是實現自由程序軋制的一種有效方法,周期橫移分定步長橫移和變步長橫移兩種方式。周期定步長橫移是指在最大橫移范圍內軋輥采用恒定步長移動。周期變步長橫移是指軋輥在最大橫移范圍內,力求在整個軋制計劃中,軋制每卷帶鋼時均有恰當的工作橫移位置,并使軋輥輥型在完成軋制計劃的最后一卷帶鋼時仍保持平滑,其與周期定步長橫移方式相比,軋輥磨損均勻,但控制過程復雜。目前主要有以下幾種精軋帶鋼局部高點相關的控制技術(1)"ROLLINGMETHODFORPREVENTINGHIGHSPOTOFHOTSTRIP,,(JP56009006),其技術關鍵通過安裝在精軋出口的多功能儀測量出帶鋼實際斷面形狀,利用過程計算機判斷是否有必要控制帶鋼局部高點,當需要時移動工作輥。本質上它是一種反饋控制方法,通過檢測的手段來決定是否需要進行工作輥橫移。(2)"OPTIMISEDSHIFTSTRATEGYASAFUNCTIONOFSTRIPWIDTH"(W02005058517),其技術關鍵提出一種優化的工作輥橫移策略,通過以帶鋼寬度為自變量的函數來決定工作輥的橫移位置,以此來控制帶鋼的邊部質量,更充分地發揮CVC軋機技術的優勢。上述現有技術存在如下問題技術(1)考慮的是通過反饋控制來抑制帶鋼局部高點,它沒有建立一個工作輥橫移方案的優化模型;技術(2)考慮的是僅通過帶鋼寬度來優化工作輥橫移策略,而且只是針對CVC軋機的橫移策略。另外,兩者都沒有以整個軋制計劃的角度來建立一個評價工作輥橫移方案優劣的優化模型。在熱連軋過程中,軋輥磨損非常嚴重,它不但影響帶鋼的表面質量,而且使帶鋼橫截面產生局部高點。因此,要對軋制計劃的安排及一個軋制周期內總軋制長度和同寬軋制長度進行限制。軋輥橫移是實現軋輥磨損均勻化的重要技術,以往人們通常采用固定模式結合人工經驗的方法對橫移方案實現調整。為滿足客戶對熱軋帶鋼產品的規格、品種、批量等方面日益多變的要求,使熱連軋機組能靈活自如地適應鋼材市場的變化,提高市場競爭力,固定模式結合人工經驗調整的方法已顯得力不從心。
            發明內容本發明的目的是為了解決現有技術存在的問題,提供一種通過工作輥橫移優化模型,確定軋制計劃每卷帶鋼的工作輥橫移位置,從而使得軋輥磨損均勻,并能分散軋輥熱凸度,擴大帶鋼凸度控制范圍以及減少帶鋼局部高點的精軋帶鋼局部高點的控制方法。本發明的目的是通過以下技術方案實現的一種精軋帶鋼局部高點的控制方法,其特征在于包括第一步驟對精軋帶鋼的熱連軋機組軋輥磨損進行計算,橫移位置決定了磨損增量在軋輥長度方向的分布位置,為得到橫移位置,先對精軋帶鋼的軋輥磨損進行計算,軋制一卷帶鋼產生的磨損增量的數學模型為Aff(i,j,iw)={aw(i,iw)‘Dwe(i,iw)+bw(i,iw)}‘A(i,j,iw)a(i'iw)‘B(i,iw)Mi'iw)·C(i,iw)式1,式中Aff(i,j,iw)在第i機架、第iw工作輥、第j點的的磨損增量,其中i機架號,熱連軋機組為末三機架F5、F6、F7配置工作輥橫移,i取值1、2、3,分別表示F5、F6、F7,j為1-N,是工作輥長度方向每隔一固定距離e取一個點的各點序號,N為軋輥長度為L的輥身長度方向上總的計算點數,N=L/e,iw上、下工作輥索引號,iw=1為上工作輥,iw=2為下工作輥,A(i,j,iw)為軋制力作用項,B(i,iw)為變形區幾何形狀影響項,C(i,iw)為軋制帶鋼長度影響項,Dwe(i,iw)第i機架第iw工作輥的直徑,aw(i,iw),bw(i,iw),α(i,iw),β(i,iw)第i機架第iw工作輥的模型系數,上述A(i,j,iw)、B(i,iw)、C(i,iw)、DWK(i,iw)在一個軋制計劃中根據機架工作輥的實測磨損值確定;第二步驟計算一個軋制計劃產生的總磨損量一個軋制計劃結束后,軋輥輥身各點的總的磨損量等于該計劃中所有帶鋼磨損增量在輥身對應點的累加,其數學式為式2中AWk(i,j,iw)為第k卷帶鋼在i機架、第iw工作輥、第j點的的磨損增量,ρ為一個軋制計劃的帶鋼卷數,軋輥橫移后,橫移位置決定了輥身1,2,...,N計算點的磨損增量;第三步驟提供工作輥橫移的優化方案,使軋輥輥型在完成軋制計劃最后一卷帶鋼時仍保持最平滑的磨損輪廓曲線,所述工作輥橫移的優化方案為一個軋制計劃內各卷帶鋼工作輥橫移位置構成的序列,所述軋輥磨損輪廓曲線為軋輥長度方向上各計算點總磨損量構成的曲線;工作輥橫移的優化方案的優劣由評價函數表示,采用軋制計劃結束時表示磨損輪廓曲線上凸起的尖峰高低與表示曲線左右對稱性的綜合指標作為橫移方案優劣的評價函數,該評價函數用下列式子表示式3其中,fitness(i,iw)為第i機架第iw工作輥橫移方案的評價值,ff(i,j,iw)為一個計劃在第i機架、第iw工作輥、第j點產生的總磨損量,max(ff(i,j,iw),20)為當i與iw相同時取出某一工作輥磨損輪廓曲線上總磨損量排在前20個點的總磨損量,mean{max(W(i,j,iw),20)}為對上述20個點的總磨損量取平均值,std{W(i,jaw)"W(i,(N_j),iw)}為對稱性指標,用磨損輪廓曲線左半部分與右半部分對稱位置相減后的序列取標準差表示,為尖峰指標,用磨損輪廓曲線上總磨損量排在前20個點的磨損均值與中間平坦區域20個點的磨損均值之差表示,C1,C2為加權系數,C1=0.6,C2=0.4;對同一個機架i,對上下工作輥評價值取平均值,此時其中fitness(1,1)、fitness(2,1)、fitness(3,1)分別為第F5、F6、F7機架上工作輥橫移方案評價值,fitness(1,2),fitness(2,2),fitness(3,2)分別為第F5、F6、F7機架下工作輥橫移方案評價值,fitnessUfitness2,fitness3分別為F5、F6、F7機架橫移方案的評價值;在一個軋制計劃內,式3中表示的磨損輪廓曲線上的尖峰與對稱性綜合指標反映該計劃內各卷帶鋼產生局部高點大小,工作輥橫移優化方案表述為一個帶約束條件的多目標優化函數,其數學描述為Min:fitnessl=Kx1,x2,…,xp)Min:fitness2=f2(xi;x2,xp)式5Min:fitness3=f3(X1,x2,···,xp)約束條件xix+1-xix<=ASmax|Xxi+1-Xix|>=ASmin式6IxixI<=M式中,fitnessUfitness2,fitness3分別為F5、F6、F7機架橫移方案的評價值,在滿足式6的前提下式5求得橫移位置序列X1,x2,…,xp,使fitness1、fitness2、fitness3達到最小,以找到工作輥橫移的優化方案,^(x1X2,…,Xp)、f2(Xl,x2,…,Xp)、f3(Xl,x2,…,xp)分別為F5、F6、F7機架由橫移位置序列映射到橫移方案評價值的函數關系,Xix為軋制第ix卷帶鋼工作輥的橫移位置,ρ為整個軋制計劃的帶鋼卷數,M為工作輥橫移的最大行程,ΔSfflin,ΔSfflax分別為工作輥橫移的最小允許步長和最大允許步長;第四步驟對工作輥橫移優化方案進行求解,求解策略1,采用橫移位置生成包減少式5中待優化變量的數目,所述橫移位置生成包用4個變量參數dpe^NXj來控制一個計劃內各卷帶鋼橫移位置的生成,其中Cl1-第一階段工作輥橫移步長,是第一階段中從前一卷帶鋼到后一卷帶鋼工作輥橫移步長,d2-第二階段工作輥橫移步長,是第二階段中從前一卷帶鋼到后一卷帶鋼工作輥橫移步長,NX-階段切換點,第一階段、第二階段的切換點,稱為階段切換點,A-—個軋制計劃內的工作輥橫移位置構成的序列表現為一個振蕩型曲線,并在第二階段中該振蕩型曲線的邊界是衰減的,A為振蕩型曲線的衰減控制率;求解策略2采用滿足指導性規則排除大部分可行的橫移方案的可行解,以獲得工作輥橫移的優化方案,所述指導性規則第一條規則,工作輥橫移第一階段移動范圍要大;第二條規則,工作輥橫移第二階段位置構成的序列表現為一個振蕩型曲線,該振蕩型曲線的邊界要以一速率衰減;第三條規則,計劃內所有帶鋼的工作輥橫移位置左右對稱;上述第一條規則的評價值fit_rulel,用下面的表達式描述Γμfit_rule\=M--式9表示對i}取整,M為所述工作輥橫移的最大行程,上述第二條規則的評d、價值fit_rule2,用下面的表達式描述(dfit—rule2=3——式10VA0J其中-A0為衰減率控制參數的參考值,取0.1,A為所述振蕩型曲線的衰減控制率,上述第三條規則的評價值fit_rule3,用下面的表達式描述1^pfit_rule>=—Υχ,ζ11ρ其中表示橫移位置序列χ”Χ2,…,Xp圍繞X軸的對稱性,P為所述一個計劃/=1的帶鋼卷數;由此,構造出一個規則適應度評價函數fit—rule=P1‘fit_rulel+p2‘fit_rule2+p3‘fit_rule3,艮PΣχ'式12式中,fit_rule為規則適應度評價值,Pl,p2,P3為加權系數,取值=P1=P2=P3=1;設共有U組可行的橫移方案的可行解,每組可行解記為<,<,ΛΧ\,其中l^k^U,相應的規則適應度評價值為fit_rulek,橫移方案三個機架的評價值分別為fitnesslk、fitness2k及fitneSS3k,將其平均值作為綜合目標評價值fitnessk=~{fitness\k+fitness2k+fitness3")^13第五步驟利用上述第1-第4步驟,工作輥進行橫移方案優化求解流程S31分別給式12中的規則適應度評價值fit_rule中的最優規則適應度評價值fit_rule_best及給由式3、4、13計算得到的橫移方案綜合目標評價值中的最優綜合目標評價值fitness_best賦一個初值10000;根據所述4個待優化參數屯、d2、NX、A的取值范圍及取值間隔生成可行的橫移方案的可行解集合<,<,ΛΤ4,#(k=1,2,"·υ),同時將可行解計數器k置初值1;S32從可行解集合中取出當前解f<,<,iVAMhS33將當前解<,<,ΛΤ,Υ送入橫移位置生成包中,由橫移位置生成包計算出當前解{^,《,ΛΟ;}對應的橫移位置序列<(ix=1,2,…,ρ);S34將當前解f<,<,ALT,#f及橫移位置序列<(ix=1,2,…,ρ)代入規則適應度評價函數式12,計算出當前解的規則適應度評價值fit_rUlek;S35:判斷規則適應度評價值是否小于安全系數與最優規則適應度評價值的乘積,即,fit_rulek<p_threshold·fitness_best是否成立,如果成立則往下執行S36,如果不成立,則執行步驟S44;其中Mfit_rule=pr\M-M_■dA+P23-A+/V一.PS36帶鋼卷數計數器ix置初值1;S37將當前帶鋼對應的橫移位置<代入上述磨損計算模型式1,得到該卷帶鋼在軋輥長度方向上各點產生的磨損增量;S38判斷帶鋼卷數計數器iχ是否大于軋制計劃長度P,如果成立則執行S40,不成立則執行S39;S39帶鋼卷數計數器ix累加1,往上執行S37;S40用上述式2對整個軋制計劃所有帶鋼在軋輥長度方向上各點的磨損增量在對應點相加得到軋制計劃結束時軋輥總的磨損量;S41將當前解對應的軋輥總的磨損量代入式3和4給出的橫移方案評價函數計算出三機架的評價值,再按式13做平均值處理后得到綜合目標評價值fitness15;S42判斷fitnessk<fitness_best是否成立,如果成立,則執行S43,否則執行S44;S43將fitnessk更新到fitness_best,同時將fit_rulek更新到fit_rule_best,并將當前解k,<,NXk,Ak;}更新到最優解,dT,NXbes',Abes'};S44將可行解計數器k累加1;S45判斷可行解計數器k是否大于所有可行解的數目U,如果成立則執行S46,否則回到S32;S46輸出最優解對應的橫移位置序列Xl,X2,…,V本發明的有益效果本發明突破了基于固定模式結合人工經驗對于軋制計劃多變難以快速調整的障礙,以整個軋制計劃的角度建立一種工作輥橫移優化模型,對一個軋制計劃,模型能夠根據不同的橫移方案計算軋制計劃結束時軋輥的磨損輪廓曲線,基于對軋輥磨損輪廓曲線尖峰及對稱性指標綜合最優的評價對不同軋制計劃的橫移方案作動態調整,確定軋制計劃每卷帶鋼的工作輥橫移位置,以適應熱軋帶鋼生產軋制計劃多變的情況,充分發揮了利用工作輥橫移均勻軋輥磨損的優勢,從而達到分散軋輥熱凸度,擴大帶鋼凸度控制范圍以及減少帶鋼局部高點的目的。本發明是一種應用于熱連軋過程的工作輥橫移方案的優化控制方法,對提高帶鋼產品的質量起了積極作用。為進一步說明本發明的上述目的、結構特點和效果,以下將結合附圖對本發明進行詳細說明。圖1為軋輥磨損計算的不同分區示意圖;圖2為本發明中橫移位置生成包的示意圖;圖3為本發明中工作輥橫移方案優化的流程示意圖;圖4為繼續圖3的流程示意具體實施例方式下面結合實施例的附圖對本發明的具體實施方式進行詳細說明。首先對精軋帶鋼的熱連軋機組軋輥磨損計算模型(僅計算軋制一卷帶鋼產生的13磨損增量)進行說明。工作輥與帶鋼間及工作輥與支撐輥間的摩擦,都會使工作輥產生磨損。磨損變化是一個緩慢積累的過程,它的計算精度對板凸度設定、板平直度控制、帶鋼表面質量改善產生影響。磨損分中部定常磨損區和邊部集中磨損區兩部分。形成邊部磨損加劇的主要原因是邊部溫度降低、負荷在邊部區域作用的增強以及軋件金屬三維流動。參見圖1,圖1為軋輥磨損計算的不同分區的示意圖,如圖所示,根據帶鋼與軋輥接觸位置的不同,將軋輥分為7個部分中部14;左、右第1邊部11、1Γ;左、右第2邊部12,12';左、右第3邊部13、13'。針對7個部分磨損狀況的不同,分別計算各個區域的磨損量。圖1中的L為帶鋼長度,S、T為帶鋼兩邊部,0帶鋼中心線。圖1中dwO,dwl,dw2,dw3為用于分區的長度;kwl,kw2為軋制力沿帶鋼寬度分布的增益系數。本實施例中dwO=162mm,dwl=28.5mm,dw2=0,dw3=30mm,kwl=1.14,kw2=1,通過dwO,dwl,dw2,dw3四個參數可以定出左右第1邊部、第2邊部、第3邊部的距離。計算軋制每卷帶鋼后造成的磨損增量的計算模型Aff(i,j,iw)={aw(i,iw)‘Dwe(i,iw)+bw(i,iw)}‘A(i,j,iw)a(i'iw)‘B(i,iw)Mi'iw)·C(i,iw)式(1)式(1)中的符號說明Aff(i,j,iw)在第i機架、第iw工作輥、第j點的的磨損增量,其中i機架號(熱連軋機組為末三機架F5、F6、F7配置工作輥橫移,共需計算末三個機架F5、F6、F7的磨損增量),i取值1、2、3,分別表示F5、F6、F7);j為1-N,是工作輥長度方向各點序號,N為磨損模型在軋輥的輥身長度方向上總的計算點數,例沿輥身每隔一固定距離e(例5mm)取一個點,如果軋輥長L為2280mm,每隔5mm取一個點,那么共需計算N=L/e=457個點,從左至右依次j的編號為1,2,...,457);iw上、下工作輥索引號(iw=1為上工作輥,iw=2為下工作輥);A(i,j,iw)為軋制力作用項;B(i,iw)為變形區幾何形狀影響項;C(i,iw)為軋制帶鋼長度影響項;Dwe(i,iw)第i機架第iw工作輥的直徑;aw(i,iw),bw(i,iw),α(i,iw),β(i,iw)第i機架第iw工作輥的模型系數。上述A(i,j,iw)、B(i,iw)、C(i,iw)、DWK(i,iw)項與軋制計劃及軋制工況有關,即,在給定的軋制規程(軋制計劃)下,每卷帶鋼磨損的基本影響因素即軋制長度、軋制力及接觸弧長等都是可以計算的,軋輥磨損計算模型中的模型系數根據機組不同機架軋輥的實測磨損值確定。上述符號在下面公式中表示的含義是相同的,不贅述。式(1)就是磨損增量的計算模型,磨損增量主要受軋制力、軋制長度及接觸弧長等因素的影響,軋制力與軋制長度乘積越大,則磨損增量越大。因為不同分區的磨損狀況不同,故在輥身各點上的磨損值不同,其沿輥身長度分布如圖1中的10所示。另外,本發明定義AW(i,j,iw)僅在帶鋼與軋輥接觸部分的帶鋼邊部S、T兩側dw3(通常為30mm)范圍內的點上才有值,其余點都是0。下面說明軋輥(工作輥)橫移位置是如何影響磨損的工作輥橫移后,必須根據每次的橫移位置確定與帶鋼實際的接觸單元,并計算其磨損增量。也就是說,橫移位置決定了磨損增量在軋輥長度方向的分布位置。一個軋制計劃產生的總磨損量計算(疊加各卷帶鋼的磨損增量)。每軋制一卷帶鋼會在軋輥產生一定的磨損,磨損大小由式1計算,在一個軋制計劃中軋制多卷帶鋼后軋輥的累積磨損就是各卷帶鋼磨損增量在輥身長度方向相應點上的疊加,磨損的累積使得軋輥表面變得凹凸不同,這是產生局部高點的原因。一個軋制計劃結束后,軋輥輥身各點的總的磨損量等于該計劃中所有帶鋼磨損增量在輥身對應點的累加,即存在式中ΔWk(i,j,iw)為第k卷帶鋼在i機架、第iw工作輥、第j點的的磨損增量;ρ為一個軋制計劃的帶鋼卷數,N為輥身長度方向上總的計算點數;i分別取1、2、3,iw分別取1、2,共需計算6個工作輥的總磨損量;工作輥橫移后,橫移位置決定了輥身1,2,...,N點上的磨損增量,因而會影響到式(2)疊加的效果。由上面的分析可知,不同的橫移方案(一個計劃內各卷帶鋼工作輥橫移位置構成的序列)決定了不同的軋輥磨損輪廓曲線(軋輥長度方向上各點總磨損量構成的曲線),即,不同的橫移方案得到不同的磨損輪廓曲線,本發明就是要提供工作輥橫移的優化方案。要得到工作輥橫移的優化方案,首先應定義橫移方案優劣與否的工作輥橫移方案的評價函數。定義工作輥橫移方案的評價函數中首先定義使軋輥輥型在完成軋制計劃最后一卷帶鋼時仍保持最平滑的磨損輪廓曲線所對應的橫移方案為最優方案。因此,可選用軋制計劃結束時磨損輪廓曲線上凸起的尖峰高低(俗稱“貓耳”,以下稱為尖峰指標)與曲線左右對稱性(以下稱為對稱性指標)的綜合指標作為橫移方案優劣的評價函數,該評價函數用下列式子表示式(3)其中,fitness(i,iw)為第i機架第iw工作輥的評價值;ff(i,j,iw)為一個計劃在第i機架、iw工作輥的第j點產生的總磨損量;max(ff(i,j,iw),20)為取出某一工作輥(當i與iw相同時)磨損輪廓曲線上總磨損量排在前20個點的總磨損量,mean{max(W(i,j,iw),20)}為對所述20個點的總磨損量取平均值;為尖峰指標,用磨損輪廓曲線上磨損量排在前20個點的磨損均值與中間平坦區域20個點的磨損均值之差表示;std{W(i,jaw)-W(i,(N_j),iw)}為對稱性指標,用磨損輪廓曲線左半部分與右半部分對稱位置相減后的序列取標準差表示(在i與iw相等時,即對同一個工作輥,計算點j與(N-j)所在的輥身位置關于軋輥中心線對稱);C1,C2為加權系數;(取值一取決于兩個指標的重要程度,可取C1=0.6,C2=0.4);上述已提到N為磨損模型在軋輥長度方向上計算的總點數。考慮到對稱性,對同一個機架i,可對上下工作輥評價值取平均值,此時其中:fitness(l,1)、fitness(2,1)、fitness(3,1)—第F5、F6、F7機架上工作輥橫移方案評價值;fitness(1,2),fitness(2,2),fitness(3,2)-—第F5、F6、F7機架下工作輥橫移方案評價值;fitness1、fitness2、fitness3分別為F5、F6、F7機架橫移方案的評價值;熱連軋機組通常為末三機架F5、F6、F7配置工作輥橫移式軋機,為保證軋制穩定性,軋制每卷帶鋼時,通常設置F5與F7橫移位置相同,F6則與F5或F7相反。在一個軋制計劃內,如果式(3)中表示的磨損輪廓曲線上的尖峰與對稱性綜合指標越好,那么該計劃內各卷帶鋼產生局部高點的可能性就越小。因此,工作輥橫移優化方案可表述為一個帶約束條件的多目標優化問題,其數學描述為Min:fitnessl=ι(X1,χ2,…,χρ)Min:fitness2=f2(X1,x2,xp)式(5)Min:fitness3=f3(X1,x2,·..,xp)約束條件χ…-XiI<=AS,>=ASmin式(6)Xi<=M式中,fitness1、fitness2、fitness3分別為F5、F6、F7機架橫移方案的評價值;fi(χι;x2,…,Xp)、f2(Xl,x2,…,Xp)、f3(Xl,x2,…,xp)分別為F5、F6、F7機架由橫移位置序列映射到橫移方案評價值的函數關系,先根據式(1)分別計算各卷帶鋼在輥身各點產生的磨損增量,再根據式(2)將對應點的磨損增量累加得到磨損輪廓曲線,最后根16據式(3)、(4)將磨損輪廓曲線代入評價函數得到評價值;Xix為軋制第ix卷帶鋼工作輥的橫移位置;ρ為整個軋制計劃的帶鋼卷數。M為工作輥橫移的最大行程,通常為200mm;ΔSfflin,ΔSfflax分別為工作輥橫移的最小允許步長和最大允許步長(本實施例最小允許橫移步長ASmin=10mm,最大允許橫移步長ASmax=50mm)。在滿足式(6)的前提下,通過式(5)找合適的橫移位置序列X1,X2,…,Xp,使fitness1、fitness2、fitness3達到最小,以找到工作輥橫移的優化方案。下面對工作輥橫移方案優化問題的求解進行說明。求解策略1,根據周期橫移的特點減少式5中待優化變量的數目。由上面的分析可知,工作輥橫移方案是一個帶約束條件的多目標優化問題,待優化變量就是軋制各卷帶鋼工作輥的橫移位置。通常一個計劃的帶鋼卷數在40個以上,如果直接求解問題,則需設計非常復雜的優化算法。根據工作輥橫移的周期性特點,本發明設計一個如圖2所示的橫移位置生成包20,減少式5中待優化變量的數目,該橫移位置生成包20用4個變量參數(dpdyNXJ)來控制一個計劃內各卷帶鋼橫移位置的生成。下面描述橫移位置生成包的計算過程。圖2中的軋制歷程為一個計劃的帶鋼卷數(如圖2所示本實施例為62卷),定義在通常一個計劃的帶鋼卷數(例,62卷)分為兩個階段進行工作輥橫移(本實施例中第一階段為軋制第1-12卷,第二階段為軋制第13-62卷),兩個階段中從前一卷帶鋼到后一卷帶鋼工作輥橫移步長分別為屯、d2,兩個階段的切換點稱為階段切換點NX(第12卷結束),即Cl1-第一階段工作輥橫移步長,第一階段從前一卷帶鋼到后一卷帶鋼工作輥橫移步長(每步移動的大小,如圖2所示,Cl1=50mm);d2-第一階段工作輥橫移步長,第二階段從前一卷帶鋼到后一卷帶鋼工作輥橫移步長(每步移動的大小,如圖2所示為,d2=30mm);NX-階段切換點,決定在軋制多少卷帶鋼后由第一階段過渡到第二階段(第12卷結束)的切換點;A-一個軋制計劃的工作輥橫移位置序列表現為為一個振蕩型曲線,并在第二階段中該振蕩型曲線的邊界(峰點位置)是衰減的,定義A為振蕩型曲線的衰減控制率,其決定振蕩型曲線邊界衰減的快慢;所述橫移位置生成包20的作業流程如下S21,第一階段工作輥橫移軋制卷數ix為1到NX;S211,軋制第1卷帶鋼時工作輥橫移位置X1=Cl1,此后工作輥橫移步長為Cl1;S212,軋制2卷到第NX卷,卷數ix從2變化到NX,其中每步都是在前一步橫移位置的基礎上加Cl1(即Xix=Xix-Jd1),但如果加Cl1后越過了定義的第一階段工作輥橫移位置邊界M(本實施例M為200mm)(即|Ud11>200),則將橫移反向進行(即令(I1=-(I1);S22,第二階段工作輥橫移軋制卷數ix為(NX+1)到ρ;S221,軋制完第NX卷帶鋼后,工作輥橫移步長變為d2(本實施例d2=30mm);S222,軋制第NX+1卷到第ρ卷(ix從NX+1變化到ρ),每步都是在前一步橫移位置的基礎上加d2(xix=XixJd2),但如果加d2后越過了橫移邊界xLimit,則將橫移反向(即令d2=-d2)。上述第二階段的橫移邊界xLimit可以表述為其中,sp-衰減起始點,決定橫移邊界衰減開始的帶鋼卷數;ep_衰減結束點,決定橫移邊界衰減結束的帶鋼卷數;即當ix大于sp并且不大于印時,橫移邊界在不斷減小;當ix大于印時,橫移邊界保持不變。只要按上述方法計算各卷帶鋼工作輥的橫移位置,并滿足ΓΔ5<d<AS那么自然就滿足約束條件式(6)。上述橫移位置生成包將工作輥橫移過程分解為前后兩個階段,并在第二階段通過指數函數B*exp(-A·(ix-sp))來控制橫移范圍的衰減。只要給定dpdpNXj參數值,就可以生成一個軋制計劃內每卷帶鋼的工作輥橫移位置。因此,上述橫移位置生成包使工作輥的橫移優化方案只與4個變量有關,即從式5中ρ個待優化變量Xl,x2,…,xp減少為4個待優化變量屯、d2、NX、Α,縮短了橫移方案優化問題的求解時間。求解策略2采用滿足指導性規則排除大部分可行的橫移方案,以獲得工作輥橫移的優化方案,所述指導性規則是在大量計算和分析的基礎上,歸納出的第一條規則,工作輥橫移第一階段移動范圍要盡可能大;第二條規則,工作輥橫移第二階段位置構成的序列表現為一個振蕩型曲線,該振蕩型曲線的邊界要以一速率衰減;第三條規則,計劃內所有帶鋼的工作輥橫移位置左右對稱性要好;分析表明,待優化變量滿足三條規則的程度越好,獲得最優解的可能性就越大。在第一階段橫移步長Cl1、第二階段橫移步長d2、階段切換點NX、衰減控制率參數A四個待優化變量中,Cl1影響橫移第一階段的移動范圍,A影響橫移第二階段移動邊界的衰減快慢,而NX、d2、Cl1及A共同影響橫移位置序列的左右對稱性。上述第一條規則的評價值fit_rulel,可用下面的表達式描述其中#表示對#取整,M為所述工作輥橫移的最大行程;式9意義第一階段橫移移動范圍越大,則(9)式越小;上述第二條規則的評價值fit_rule2,可用下面的表達式描述(Αλfit—rule2=3——式(10)VAoJ其中=Atl為衰減率控制參數的參考值,取0.1,A為所述振蕩型曲線的衰減控制率,18式10意義A越大,(10)式越小,橫移邊界衰減越快;上述第三條規則fit_rule3,可用下面的表達式描述其中Σχ·表示橫移位置序列X1,X2,…,Xp圍繞X軸的對稱性,,P為所述一個計劃的帶鋼卷數。式11意義橫移位置序列Xl,x2,…,Xp圍繞χ軸對稱性越好,式(11)越〃結論①Cl1滿足第一條規則的程度越好,等價于(9)式越小;②A滿足第二條規則的程度越好,等價于(10)式越小;③NX、C^d1及A滿足第三條規則的程度越好,等價于(11)式越小。由此,可構造出一個規則適應度評價函數fit_rule式(12)式中,Pl,p2,p3為加權系數,取值大小決定于每條規則的重要程度,通常取Pl=P2=p3=1)。結論戲、d2、NX及A滿足三條指導性規則的程度越好,等價于式(12)取值越小。求解橫移方案優化問題前,首先要確定4個待優化參數C^dpNXj的取值范圍及取值間隔(舉例假設dl取值范圍為40到50,取值間隔為2,則dl可取到40,42,44,46,48,50),見表1。表1待優化變量的取值范圍及取值間隔根據表1,可以得到4個變量的不同組合,所有的組合構成可行的橫移方案的全搜索空間。假設共有U組可行的橫移方案的可行解,每組可行解為其中1彡k彡U,相應的規則適應度評價值為fit_rulek,三個機架F5、F6、F7的工作輥橫移方案評價值分別為fitneSSlk、fitneSS2k及fitneSS3k,將其平均值作為工作輥橫移方案綜合目標評價值fitnessk~^{fitnessV+fitness2k+fitness3k、式(丄3)下面根據圖3、圖4對本發明的工作輥橫移方案優化求解流程進行說明S31分別給式12中的規則適應度評價值fit_rule中的最優規則適應度評價值fit_rule_best及給由式3、4、13計算得到的橫移方案綜合目標評價值中的最優綜合目標評價值fitness_best賦一個很大的初值10000;根據所述4個待優化參數屯、d2、NX、A的取值范圍及取值間隔生成可行的橫移方案的可行解集合<,<,Λ「Ζ\Υ;[(k=1,2,…U);同時將可行解計數器k置初值1;S32從可行解集合中取出當前解丨,Λ^Μ;};S33將當前解<,<,7νΖ\#;}送入橫移位置生成包中,由橫移位置生成包計算出當前解{<,<,旭〃,#;|對應的橫移位置序列<(ix=1,2,…,ρ);S34將當前解<,<,ΤνΤ,#;[及橫移位置序列<(ix=1,2,…,ρ)代入規則適應度評價函數式12,計算出當前解的規則適應度評價值fit_rUlek;S35:判斷規則適應度評價值是否小于安全系數與最優規則適應度評價值的乘積,即,fit_rulek<p_threshold·fitness_best是否成立,如果成立則往下執行S36,如果不成立,則執行步驟S44。所述安全系數p_threshold是為擴大最優解的可能范圍而設置的系數,取值為1.2,使得規則適應度評價值除了滿足fit_rulek<fitness_best的可行解外,fitness一best(fit_rulek<p_threshold·fitness一best白勺"SJif角軍。S36帶鋼卷數計數器ix置初值1;S37將當前帶鋼對應的橫移位置<代入上述磨損計算模型式1,得到該卷帶鋼在軋輥長度方向上各點產生的磨損增量;S38判斷帶鋼卷數計數器ix是否大于軋制計劃長度P,如果成立則執行S40,不成立則執行S39;S39帶鋼卷數計數器ix累加1,往上執行S37;S40用上述式2對整個軋制計劃所有帶鋼在軋輥長度方向上各點的磨損增量在對應點相加得到軋制計劃結束時軋輥總的磨損量;S41將當前解對應的軋輥總的磨損量代入式3和4給出的橫移方案評價函數計算出三機架的評價值,再按式13做平均值處理后得到綜合目標評價值fitness15;S42判斷fitnessk<fitness_best是否成立,如果成立,則執行S43,否則執行S44;S43將fitnessk更新到fitness_best,同時將fit_rulek更新到fit_rule_best,并將當前解{<,<,NXk,Ak1更新到最優解fi/廣,db;s',NX"-',Ah-'};S44將可行解計數器k累加1;S45判斷可行解計數器k是否大于所有可行解的數目U,如果成立則執行S46,否則回到S32;S46輸出最優解,db;“,NX-,Abesl;}對應的橫移位置序列X1,X2,…,Xp。上述步驟S36到S41相當于完成了從橫移位置序列x/(i=1,2,…,ρ)到三機架評價值的計算過程。由于綜合目標評價值的計算需要經歷每卷帶鋼的磨損計算、軋輥對應點磨損增量疊加、磨損輪廓曲線評價值計算等多個過程,數據處理量大,計算時間長;而規則適應度評價值計算得相對簡單,不涉及到磨損計算,數據處理量小,計算時間短。通過先計算規則適應度評價值,并與當前最佳的規則適應度評價值對比,將大部分可行解排除在綜合目標評價值計算之前,這樣的優化算法設計可以極大節省找到橫移方案最優解的計算時間。本
            技術領域
            中的普通技術人員應當認識到,以上的實施例僅是用來說明本發明的目的,而并非用作對本發明的限定,只要在本發明的實質范圍內,對以上所述實施例的變化、變型都將落在本發明的權利要求的范圍內。權利要求一種精軋帶鋼局部高點的控制方法,其特征在于包括第一步驟對精軋帶鋼的熱連軋機組軋輥磨損進行計算,橫移位置決定了磨損增量在軋輥長度方向的分布位置,為得到橫移位置,先對精軋帶鋼的軋輥磨損進行計算,軋制一卷帶鋼產生的磨損增量的數學模型為ΔW(i,j,iw)={aw(i,iw)·DWR(i,iw)+bw(i,iw)}·A(i,j,iw)a(i,iw)·B(i,iw)β(i,iw)·C(i,iw)式1,式中ΔW(i,j,iw)在第i機架、第iw工作輥、第j點的的磨損增量,其中i機架號,熱連軋機組為末三機架F5、F6、F7配置工作輥橫移,i取值1、2、3,分別表示F5、F6、F7,j為1N,是工作輥長度方向每隔一固定距離e取一個點的各點序號,N為軋輥長度為L的輥身長度方向上總的計算點數,N=L/e,iw上、下工作輥索引號,iw=1為上工作輥,iw=2為下工作輥,A(i,j,iw)為軋制力作用項,B(i,iw)為變形區幾何形狀影響項,C(i,iw)為軋制帶鋼長度影響項,DWR(i,iw)第i機架第iw工作輥的直徑,aw(i,iw),bw(i,iw),α(i,iw),β(i,iw)第i機架第iw工作輥的模型系數,上述A(i,j,iw)、B(i,iw)、C(i,iw)、DWR(i,iw)在一個軋制計劃中根據機架工作輥的實測磨損值確定;第二步驟計算一個軋制計劃產生的總磨損量一個軋制計劃結束后,軋輥輥身各點的總的磨損量等于該計劃中所有帶鋼磨損增量在輥身對應點的累加,其數學式為;<mfencedopen='{'close=''><mtable><mtr><mtd><mi>W</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mn>1</mn><mo>,</mo><mi>iw</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msup><mi>&Delta;W</mi><mi>k</mi></msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mn>1</mn><mo>,</mo><mi>iw</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>W</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mn>2</mn><mo>,</mo><mi>iw</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msup><mi>&Delta;W</mi><mi>k</mi></msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mn>2</mn><mo>,</mo><mi>iw</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><mi>W</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>N</mi><mo>,</mo><mi>iw</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msup><mi>&Delta;W</mi><mi>k</mi></msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>N</mi><mo>,</mo><mi>iw</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>式2,式2中ΔWk(i,j,iw)為第k卷帶鋼在i機架、第iw工作輥、第j點的的磨損增量,p為一個軋制計劃的帶鋼卷數,軋輥橫移后,橫移位置決定了輥身1,2,...,N計算點的磨損增量;第三步驟提供工作輥橫移的優化方案,使軋輥輥型在完成軋制計劃最后一卷帶鋼時仍保持最平滑的磨損輪廓曲線,所述工作輥橫移的優化方案為一個軋制計劃內各卷帶鋼工作輥橫移位置構成的序列,所述軋輥磨損輪廓曲線為軋輥長度方向上各計算點總磨損量構成的曲線;工作輥橫移的優化方案的優劣由評價函數表示,采用軋制計劃結束時表示磨損輪廓曲線上凸起的尖峰高低與表示曲線左右對稱性的綜合指標作為橫移方案優劣的評價函數,該評價函數用下列式子表示<mrow><mi>fitness</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>iw</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>c</mi><mn>1</mn></msub><mo>&CenterDot;</mo><mi>std</mi><mo>{</mo><mi>W</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>iw</mi><mo>)</mo></mrow><mo>-</mo><mi>W</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mi>j</mi><mo>)</mo></mrow><mo>,</mo><mi>iw</mi><mo>)</mo></mrow><mo>}</mo><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><mo>&CenterDot;</mo><mo>{</mo><mi>mean</mi><mo>{</mo><mi>max</mi><mrow><mo>(</mo><mi>W</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>iw</mi><mo>)</mo></mrow><mo>,</mo><mn>20</mn><mo>)</mo></mrow><mo>}</mo><mo>-</mo><mfrac><mn>1</mn><mn>20</mn></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mrow><mo>(</mo><mi>N</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mo>-</mo><mn>10</mn></mrow><mrow><mrow><mo>(</mo><mi>N</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mo>+</mo><mn>9</mn></mrow></munderover><mi>W</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>iw</mi><mo>)</mo></mrow><mo>}</mo></mrow>式3其中,fitness(i,iw)為第i機架第iw工作輥橫移方案的評價值,W(i,j,iw)為一個計劃在第i機架、第iw工作輥、第j點產生的總磨損量,max(W(i,j,iw),20)為當i與iw相同時取出某一工作輥磨損輪廓曲線上總磨損量排在前20個點的總磨損量,mean{max(W(i,j,iw),20)}為對上述20個點的總磨損量取平均值,std{W(i,j,iw)W(i,(Nj),iw)}為對稱性指標,用磨損輪廓曲線左半部分與右半部分對稱位置相減后的序列取標準差表示,為尖峰指標,用磨損輪廓曲線上總磨損量排在前20個點的磨損均值與中間平坦區域20個點的磨損均值之差表示,c1,c2為加權系數,c1=0.6,c2=0.4;對同一個機架i,對上下工作輥評價值取平均值,此時<mfencedopen='{'close=''><mtable><mtr><mtd><mi>fitness</mi><mn>1</mn><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>&CenterDot;</mo><mo>[</mo><mi>fitness</mi><mrow><mo>(</mo><mn>1,1</mn><mo>)</mo></mrow><mo>+</mo><mi>fitness</mi><mrow><mo>(</mo><mn>1,2</mn><mo>)</mo></mrow><mo>]</mo></mtd></mtr><mtr><mtd><mi>fitness</mi><mn>2</mn><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>&CenterDot;</mo><mo>[</mo><mi>fitness</mi><mrow><mo>(</mo><mn>2,1</mn><mo>)</mo></mrow><mo>+</mo><mi>fitness</mi><mrow><mo>(</mo><mn>2,2</mn><mo>)</mo></mrow><mo>]</mo></mtd></mtr><mtr><mtd><mi>fitness</mi><mn>3</mn><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>&CenterDot;</mo><mo>[</mo><mi>fitness</mi><mrow><mo>(</mo><mn>3,1</mn><mo>)</mo></mrow><mo>+</mo><mi>fitness</mi><mrow><mo>(</mo><mn>3,2</mn><mo>)</mo></mrow><mo>]</mo></mtd></mtr></mtable></mfenced>式4其中fitness(1,1)、fitness(2,1)、fitness(3,1)分別為第F5、F6、F7機架上工作輥橫移方案評價值,fitness(1,2)、fitness(2,2)、fitness(3,2)分別為第F5、F6、F7機架下工作輥橫移方案評價值,fitness1、fitness2、fitness3分別為F5、F6、F7機架橫移方案的評價值;在一個軋制計劃內,式3中表示的磨損輪廓曲線上的尖峰與對稱性綜合指標反映該計劃內各卷帶鋼產生局部高點大小,工作輥橫移優化方案表述為一個帶約束條件的多目標優化函數,其數學描述為Minfitness1=f1(x1,x2,…,xp)Minfitness2=f2(x1,x2,…,xp)式5Minfitness3=f3(x1,x2,…,xp)約束條件|xix+1xix|<=ΔSmax|xxi+1xix|>=ΔSmin式6|xix|<=M式中,fitness1、fitness2、fitness3分別為F5、F6、F7機架橫移方案的評價值,在滿足式6的前提下式5求得橫移位置序列x1,x2,…,xp,使fitness1、fitness2、fitness3達到最小,以找到工作輥橫移的優化方案,f1(x1,x2,…,xp)、f2(x1,x2,…,xp)、f3(x1,x2,…,xp)分別為F5、F6、F7機架由橫移位置序列映射到橫移方案評價值的函數關系,xix為軋制第ix卷帶鋼工作輥的橫移位置,p為整個軋制計劃的帶鋼卷數,M為工作輥橫移的最大行程,ΔSmin、ΔSmax分別為工作輥橫移的最小允許步長和最大允許步長;第四步驟對工作輥橫移優化方案進行求解,求解策略1,采用橫移位置生成包減少式5中待優化變量的數目,所述橫移位置生成包用4個變量參數d1、d2、NX、A來控制一個計劃內各卷帶鋼橫移位置的生成,其中d1第一階段工作輥橫移步長,是第一階段中從前一卷帶鋼到后一卷帶鋼工作輥橫移步長,d2第二階段工作輥橫移步長,是第二階段中從前一卷帶鋼到后一卷帶鋼工作輥橫移步長,NX階段切換點,第一階段、第二階段的切換點,稱為階段切換點,A一個軋制計劃內的工作輥橫移位置構成的序列表現為一個振蕩型曲線,并在第二階段中該振蕩型曲線的邊界是衰減的,A為振蕩型曲線的衰減控制率;求解策略2采用滿足指導性規則排除大部分可行的橫移方案的可行解,以獲得工作輥橫移的優化方案,所述指導性規則第一條規則,工作輥橫移第一階段移動范圍要大;第二條規則,工作輥橫移第二階段位置構成的序列表現為一個振蕩型曲線,該振蕩型曲線的邊界要以一速率衰減;第三條規則,計劃內所有帶鋼的工作輥橫移位置左右對稱;上述第一條規則的評價值fit_rule1,用下面的表達式描述<mrow><mi>fit</mi><mo>_</mo><mi>rule</mi><mn>1</mn><mo>=</mo><mi>M</mi><mo>-</mo><mo>[</mo><mfrac><mi>M</mi><msub><mi>d</mi><mn>1</mn></msub></mfrac><mo>]</mo><mo>&CenterDot;</mo><msub><mi>d</mi><mn>1</mn></msub></mrow>式9其中表示對取整,M為所述工作輥橫移的最大行程,上述第二條規則的評價值fit_rule2,用下面的表達式描述<mrow><mi>fit</mi><mo>_</mo><mi>rule</mi><mn>2</mn><mo>=</mo><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mfrac><mi>A</mi><msub><mi>A</mi><mn>0</mn></msub></mfrac><mo>)</mo></mrow></mrow>式10其中A0為衰減率控制參數的參考值,取0.1,A為所述振蕩型曲線的衰減控制率,上述第三條規則的評價值fit_rule3,用下面的表達式描述<mrow><mi>fit</mi><mo>_</mo><mi>rule</mi><mn>3</mn><mo>=</mo><mfrac><mn>1</mn><mi>p</mi></mfrac><mo>|</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msub><mi>x</mi><mi>i</mi></msub><mo>|</mo></mrow>式11其中表示橫移位置序列x1,x2,…,xp圍繞x軸的對稱性,p為所述一個計劃的帶鋼卷數;由此,構造出一個規則適應度評價函數fit_rule=p1·fit_rule1+p2·fit_rule2+p3·fit_rule3,即<mrow><mi>fit</mi><mo>_</mo><mi>rule</mi><mo>=</mo><msub><mi>p</mi><mn>1</mn></msub><mo>&CenterDot;</mo><mo>{</mo><mi>M</mi><mo>-</mo><mo>[</mo><mfrac><mi>M</mi><msub><mi>d</mi><mn>1</mn></msub></mfrac><mo>]</mo><mo>&CenterDot;</mo><msub><mi>d</mi><mn>1</mn></msub><mo>}</mo><mo>+</mo><msub><mi>p</mi><mn>2</mn></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mfrac><mi>A</mi><msub><mi>A</mi><mn>0</mn></msub></mfrac><mo>)</mo></mrow><mo>+</mo><msub><mi>p</mi><mn>3</mn></msub><mo>&CenterDot;</mo><mfrac><mn>1</mn><mi>p</mi></mfrac><mo>&CenterDot;</mo><mo>|</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msub><mi>x</mi><mi>i</mi></msub><mo>|</mo></mrow>式12式中,fit_rule為規則適應度評價值,p1,p2,p3為加權系數,取值p1=p2=p3=1;設共有U組可行的橫移方案的可行解,每組可行解記為其中1≤k≤U,相應的規則適應度評價值為fit_rulek,橫移方案三個機架的評價值分別為fitness1k、fitness2k及fitness3k,將其平均值作為綜合目標評價值<mrow><msup><mi>fitness</mi><mi>k</mi></msup><mo>=</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mrow><mo>(</mo><msup><mrow><mi>fitness</mi><mn>1</mn></mrow><mi>k</mi></msup><mo>+</mo><msup><mrow><mi>fitness</mi><mn>2</mn></mrow><mi>k</mi></msup><mo>+</mo><msup><mrow><mi>fitness</mi><mn>3</mn></mrow><mi>k</mi></msup><mo>)</mo></mrow></mrow>式13第五步驟利用上述第1第4步驟,工作輥進行橫移方案優化求解流程S31分別給式12中的規則適應度評價值fit_rule中的最優規則適應度評價值fit_rule_best及給由式3、4、13計算得到的橫移方案綜合目標評價值中的最優綜合目標評價值fitness_best賦一個初值10000;根據所述4個待優化參數d1、d2、NX、A的取值范圍及取值間隔生成可行的橫移方案的可行解集合(k=1,2,...U),同時將可行解計數器k置初值1;S32從可行解集合中取出當前解S33將當前解送入橫移位置生成包中,由橫移位置生成包計算出當前解對應的橫移位置序列(ix=1,2,…,p);S34將當前解及橫移位置序列(ix=1,2,…,p)代入規則適應度評價函數式12,計算出當前解的規則適應度評價值fit_rulek;S35判斷規則適應度評價值是否小于安全系數與最優規則適應度評價值的乘積,即,fit_rulek<p_threshold·fitness_best是否成立,如果成立則往下執行S36,如果不成立,則執行步驟S44;S36帶鋼卷數計數器ix置初值1;S37將當前帶鋼對應的橫移位置代入上述磨損計算模型式1,得到該卷帶鋼在軋輥長度方向上各點產生的磨損增量;S38判斷帶鋼卷數計數器ix是否大于軋制計劃長度p,如果成立則執行S40,不成立則執行S39;S39帶鋼卷數計數器ix累加1,往上執行S37;S40用上述式2對整個軋制計劃所有帶鋼在軋輥長度方向上各點的磨損增量在對應點相加得到軋制計劃結束時軋輥總的磨損量;S41將當前解對應的軋輥總的磨損量代入式3和4給出的橫移方案評價函數計算出三機架的評價值,再按式13做平均值處理后得到綜合目標評價值fitnessk;S42判斷fitnessk<fitness_best是否成立,如果成立,則執行S43,否則執行S44;S43將fitnessk更新到fitness_best,同時將fit_rulek更新到fit_rule_best,并將當前解更新到最優解S44將可行解計數器k累加1;S45判斷可行解計數器k是否大于所有可行解的數目U,如果成立則執行S46,否則回到S32;S46輸出最優解對應的橫移位置序列x1,x2,…,xp。F200910053608XC0000022.tif,F200910053608XC0000042.tif,F200910053608XC0000043.tif,F200910053608XC0000046.tif,F200910053608XC0000052.tif,F200910053608XC0000054.tif,F200910053608XC0000055.tif,F200910053608XC0000056.tif,F200910053608XC0000057.tif,F200910053608XC0000058.tif,F200910053608XC0000059.tif,F200910053608XC00000510.tif,F200910053608XC00000511.tif,F200910053608XC0000061.tif,F200910053608XC0000062.tif,F200910053608XC0000063.tif2.如權利要求1所述的精軋帶鋼局部高點的控制方法,其特征在于所述式1中,Aff(i,j,iw)在帶鋼與軋輥接觸部分的帶鋼邊部兩側為30mm范圍內的點上有值,其余點都是為0。3.如權利要求1所述的精軋帶鋼局部高點的控制方法,其特征在于所述第一步驟中熱連軋機組的末三機架F5、F6、F7配置工作輥橫移是指,在軋制每卷帶鋼時,設置F5與F7橫移位置相同,F6則與F5或F7相反。4.如權利要求1所述的精軋帶鋼局部高點的控制方法,其特征在于所述第三步驟中工作輥橫移的最大行程M為200mm,最小允許橫移步長ΔSmin=10mm,最大允許橫移步長ΔSmax=50_。5.如權利要求1或5所述的精軋帶鋼局部高點的控制方法,其特征在于所述第四步驟中橫移位置生成包的作業流程如下S21,第一階段工作輥橫移軋制卷數ix為1到NX;S211,軋制第1卷帶鋼,工作輥橫移位置X1=Cl1,此后工作輥橫移步長為Cl1;S212,軋制第2卷到第NX卷,卷數ix從2變化到NX,其中每卷都是在前一卷橫移位置的基礎上加屯,Xix=XixJd1,當Xix=XixJd1大于工作輥橫移的最大行程M時,則將橫移反向進行,即令Cl1=-Cl1;S22,第二階段工作輥橫移軋制卷數ix為(NX+1)到ρ;S221,軋制第NX卷帶鋼后,工作輥橫移步長為d2;S222,從第NX+1卷到第ρ卷(ix從NX+1變化到ρ),每卷都是在前一卷帶鋼橫移位置的基礎上加d2,Xix=XixJd2,但如果加d2后越過了橫移邊界xLimit,則將橫移反向進行,即令d2=-d2;所述橫移位置生成包在第二階段的橫移邊界xLimit表述為根據上述作業流程計算各卷帶鋼工作輥的橫移位置,并滿足Wm<d<AS式此時,達到第3步驟中式6的約束條件;sp-衰減起始點,決定橫移邊界衰減開始的帶鋼卷數;ep-衰減結束點,決定橫移邊界衰減結束的帶鋼卷數。6.如權利要求1所述的精軋帶鋼局部高點的控制方法,其特征在于所述步驟S35中的安全系數p_threshold是為擴大最優解的可能范圍而設置的系數,使得規則適應度評價值除了滿足fit_rUlek<fitness_best的可行解外,還允許滿足fitness_best^fit_rulek<p_threshold·fitness_best的可行角軍。7.如權利要求6所述的精軋帶鋼局部高點的控制方法,其特征在于所述p_threShold取值為1.2。全文摘要本發明涉及一種精軋帶鋼局部高點的控制方法,是以整個軋制計劃的角度建立一種工作輥橫移優化模型,對一個軋制計劃,模型能夠根據不同的橫移方案計算軋制計劃結束時軋輥的磨損輪廓曲線,基于對軋輥磨損輪廓曲線尖峰及對稱性綜合指標最優的評價對不同軋制計劃的橫移方案作動態調整,確定軋制計劃每卷帶鋼的工作輥橫移位置,以適應熱軋帶鋼生產軋制計劃多變的情況,充分發揮了利用工作輥橫移均勻軋輥磨損的優勢,從而達到分散軋輥熱凸度,擴大帶鋼凸度控制范圍以及減少帶鋼局部高點的目的。本發明的方法對提高熱連軋帶鋼產品的質量起了積極作用。文檔編號B21B37/00GK101927264SQ200910053608公開日2010年12月29日申請日期2009年6月23日優先權日2009年6月23日發明者張健民,張曉峰,李維剛,陳水宣申請人:寶山鋼鐵股份有限公司
            網友詢問留言 已有0條留言
            • 還沒有人留言評論。精彩留言會獲得點贊!
            1
            婷婷六月激情在线综合激情,亚洲国产大片,久久中文字幕综合婷婷,精品久久久久久中文字幕,亚洲一区二区三区高清不卡,99国产精品热久久久久久夜夜嗨 ,欧美日韩亚洲综合在线一区二区,99国产精品电影,伊人精品线视天天综合,精品伊人久久久大香线蕉欧美
            亚洲精品1区 国产成人一级 91精品国产欧美一区二区 亚洲精品乱码久久久久久下载 国产精品久久久久久久伊一 九色国产 国产精品九九视频 伊人久久成人爱综合网 欧美日韩亚洲区久久综合 欧美日本一道免费一区三区 夜夜爽一区二区三区精品 欧美日韩高清一区二区三区 国产成人av在线 国产精品对白交换绿帽视频 国产视频亚洲 国产在线欧美精品 国产精品综合网 国产日韩精品欧美一区色 国产日韩精品欧美一区喷 欧美日韩在线观看区一二 国产区精品 欧美视频日韩视频 中文字幕天天躁日日躁狠狠躁97 视频一二三区 欧美高清在线精品一区二区不卡 国产精品揄拍一区二区久久 99久久综合狠狠综合久久aⅴ 亚洲乱码视频在线观看 日韩在线第二页 亚洲精品无码专区在线播放 成人亚洲网站www在线观看 欧美三级一区二区 99久久精品免费看国产高清 91麻豆国产在线观看 最新日韩欧美不卡一二三区 成人在线观看不卡 日韩国产在线 在线亚洲精品 亚洲午夜久久久久中文字幕 国产精品成人久久久久久久 精品国产一区二区在线观看 欧美精品国产一区二区三区 中文在线播放 亚洲第一页在线视频 国产午夜精品福利久久 九色国产 精品国产九九 国产永久视频 久久精品人人做人人综合试看 国产一区二区三区免费观看 亚洲精品国产电影 9999热视频 国产精品资源在线 麻豆久久婷婷国产综合五月 国产精品免费一级在线观看 亚洲国产一区二区三区青草影视 中文在线播放 国产成人综合在线 国产在线观看色 国产亚洲三级 国产片一区二区三区 久久99精品久久久久久牛牛影视 亚洲欧美日韩国产 四虎永久免费网站 国产一毛片 国产精品视频在 九九热在线精品 99精品福利视频 色婷婷色99国产综合精品 97成人精品视频在线播放 精品久久久久久中文字幕 亚洲欧美一区二区三区孕妇 亚洲欧美成人网 日韩高清在线二区 国产尤物在线观看 在线不卡一区二区 91网站在线看 韩国精品福利一区二区 欧美日韩国产成人精品 99热精品久久 国产精品免费视频一区 高清视频一区 精品九九久久 欧美日韩在线观看免费 91欧美激情一区二区三区成人 99福利视频 亚洲国产精品91 久热国产在线 精品久久久久久中文字幕女 国产精品久久久久久久久99热 成人自拍视频网 国产精品视频久久久久久 久久影院国产 国产玖玖在线观看 99精品在线免费 亚洲欧美一区二区三区导航 久久久久久久综合 国产欧美日韩精品高清二区综合区 国产精品视频自拍 亚洲一级片免费 久久久久久九九 国产欧美自拍视频 视频一区二区在线观看 欧美日韩一区二区三区久久 中文在线亚洲 伊人热人久久中文字幕 日韩欧美亚洲国产一区二区三区 欧美亚洲国产成人高清在线 欧美日韩国产码高清综合人成 国产性大片免费播放网站 亚洲午夜综合网 91精品久久一区二区三区 国产无套在线播放 国产精品视频网站 国产成人亚洲精品老王 91在线网站 国产视频97 欧美黑人欧美精品刺激 国产一区二区三区免费在线视频 久久久国产精品免费看 99re6久精品国产首页 久久精品91 国产成人一级 国产成人精品曰本亚洲 日本福利在线观看 伊人成综合网 久久综合一本 国产综合久久久久久 久久精品成人免费看 久久福利 91精品国产91久久久久久麻豆 亚洲精品成人在线 亚洲伊人久久精品 欧美日本二区 国产永久视频 国产一区二 一区二区福利 国产一毛片 亚洲精品1区 毛片一区二区三区 伊人久久大香线蕉综合影 国产欧美在线观看一区 亚洲国产欧洲综合997久久 国产一区二区免费视频 国产91精品对白露脸全集观看 久久亚洲国产伦理 欧美成人伊人久久综合网 亚洲性久久久影院 久久99国产精一区二区三区! 91精品国产欧美一区二区 欧美日韩亚洲区久久综合 日韩精品一二三区 久久久夜色精品国产噜噜 国产在线精品福利91香蕉 久久久久久久亚洲精品 97se色综合一区二区二区 91国语精品自产拍在线观看性色 91久久国产综合精品女同我 日韩中文字幕a 国产成人亚洲日本精品 久久国产精品-国产精品 久久国产经典视频 久久国产精品伦理 亚洲第一页在线视频 国产精品久久久久三级 日韩毛片网 久久免费高清视频 麻豆国产在线观看一区二区 91麻豆国产福利在线观看 国产成人精品男人的天堂538 一区二区三区中文字幕 免费在线视频一区 欧美日韩国产成人精品 国产综合网站 国产资源免费观看 亚洲精品亚洲人成在线播放 精品久久久久久中文字幕专区 亚洲人成人毛片无遮挡 国产一起色一起爱 国产香蕉精品视频在 九九热免费观看 日韩亚洲欧美一区 九九热精品在线观看 精品久久久久久中文字幕专区 亚洲欧美自拍偷拍 国产精品每日更新 久久久久国产一级毛片高清板 久久天天躁狠狠躁夜夜中文字幕 久久精品片 日韩在线毛片 国产成人精品本亚洲 国产成人精品一区二区三区 九九热在线观看 国产r级在线观看 国产欧美日韩精品高清二区综合区 韩国电影一区二区 国产精品毛片va一区二区三区 五月婷婷伊人网 久久一区二区三区免费 一本色道久久综合狠狠躁篇 亚洲综合色站 国产尤物在线观看 亚洲一区亚洲二区 免费在线视频一区 欧洲精品视频在线观看 日韩中文字幕a 中文字幕日本在线mv视频精品 91精品在线免费视频 精品国产免费人成在线观看 精品a级片 中文字幕日本在线mv视频精品 日韩在线精品视频 婷婷丁香色 91精品国产高清久久久久 国产成人精品日本亚洲直接 五月综合视频 欧美日韩在线亚洲国产人 精液呈暗黄色 亚洲乱码一区 久久精品中文字幕不卡一二区 亚洲天堂精品在线 激情婷婷综合 国产免费久久精品久久久 国产精品亚洲二区在线 久久免费播放视频 五月婷婷丁香综合 在线亚洲欧美日韩 久久免费精品高清麻豆 精品久久久久久中文字幕 亚洲一区网站 国产精品福利社 日韩中文字幕免费 亚洲综合丝袜 91精品在线播放 国产精品18 亚洲日日夜夜 伊人久久大香线蕉综合影 亚洲精品中文字幕乱码影院 亚洲一区二区黄色 亚洲第一页在线视频 一区二区在线观看视频 国产成人福利精品视频 亚洲高清二区 国内成人免费视频 精品亚洲性xxx久久久 国产精品合集一区二区三区 97av免费视频 国产一起色一起爱 国产区久久 国产资源免费观看 99精品视频免费 国产成人一级 国产精品九九免费视频 欧美91精品久久久久网免费 99热国产免费 久久精品色 98精品国产综合久久 久久精品播放 中文字幕视频免费 国产欧美日韩一区二区三区在线 精品久久蜜桃 国产小视频精品 一本色道久久综合狠狠躁篇 91在线免费观看 亚洲精品区 伊人成综合网 伊人热人久久中文字幕 伊人黄色片 99国产精品热久久久久久夜夜嗨 久久免费精品视频 亚洲一区二区三区高清不卡 久久久久国产一级毛片高清板 国产片一区二区三区 久久狠狠干 99久久婷婷国产综合精品电影 国产99区 国产精品成人久久久久 久久狠狠干 青青国产在线观看 亚洲高清国产拍精品影院 国产精品一区二区av 九九热在线免费视频 伊人久久国产 国产精品久久久久久久久久一区 在线观看免费视频一区 国产精品自在在线午夜区app 国产精品综合色区在线观看 国产毛片久久久久久国产毛片 97国产免费全部免费观看 国产精品每日更新 国产尤物视频在线 九九视频这里只有精品99 一本一道久久a久久精品综合 久久综合给会久久狠狠狠 国产成人精品男人的天堂538 欧美一区二区高清 毛片一区二区三区 国产欧美日韩在线观看一区二区三区 在线国产二区 欧美不卡网 91在线精品中文字幕 在线国产福利 国内精品91久久久久 91亚洲福利 日韩欧美国产中文字幕 91久久精品国产性色也91久久 亚洲性久久久影院 欧美精品1区 国产热re99久久6国产精品 九九热免费观看 国产精品欧美日韩 久久久久国产一级毛片高清板 久久国产经典视频 日韩欧美亚洲国产一区二区三区 欧美亚洲综合另类在线观看 国产精品自在在线午夜区app 97中文字幕在线观看 视频一二三区 精品国产一区在线观看 国产欧美日韩在线一区二区不卡 欧美一区二三区 伊人成人在线观看 国内精品91久久久久 97在线亚洲 国产在线不卡一区 久久久全免费全集一级全黄片 国产精品v欧美精品∨日韩 亚洲毛片网站 在线不卡一区二区 99re热在线视频 久久激情网 国产毛片一区二区三区精品 久久亚洲综合色 中文字幕视频免费 国产视频亚洲 婷婷伊人久久 国产一区二区免费播放 久久99国产精品成人欧美 99国产在线视频 国产成人免费视频精品一区二区 国产不卡一区二区三区免费视 国产码欧美日韩高清综合一区 久久精品国产主播一区二区 国产一区电影 久久精品国产夜色 国产精品国产三级国产 日韩一区二区三区在线 久久97久久97精品免视看 久久国产免费一区二区三区 伊人久久大香线蕉综合电影网 99re6久精品国产首页 久久激情网 亚洲成人高清在线 国产精品网址 国产成人精品男人的天堂538 香蕉国产综合久久猫咪 国产专区中文字幕 91麻豆精品国产高清在线 久久国产经典视频 国产精品成人va在线观看 国产精品爱啪在线线免费观看 日本精品久久久久久久久免费 亚洲综合一区二区三区 久久五月网 精品国产网红福利在线观看 久久综合亚洲伊人色 亚洲国产精品久久久久久网站 在线日韩国产 99国产精品热久久久久久夜夜嗨 国产综合精品在线 国产区福利 精品亚洲综合久久中文字幕 国产制服丝袜在线 毛片在线播放网站 在线观看免费视频一区 国产精品久久久精品三级 亚洲国产电影在线观看 最新日韩欧美不卡一二三区 狠狠综合久久综合鬼色 日本精品1在线区 国产日韩一区二区三区在线播放 欧美日韩精品在线播放 亚洲欧美日韩国产一区二区三区精品 久久综合久久网 婷婷六月激情在线综合激情 亚洲乱码一区 国产专区91 97av视频在线观看 精品久久久久久中文字幕 久久五月视频 国产成人福利精品视频 国产精品网址 中文字幕视频在线 精品一区二区三区免费视频 伊人手机在线视频 亚洲精品中文字幕乱码 国产在线视频www色 色噜噜国产精品视频一区二区 精品亚洲成a人在线观看 国产香蕉尹人综合在线 成人免费一区二区三区在线观看 国产不卡一区二区三区免费视 欧美精品久久天天躁 国产专区中文字幕 久久精品国产免费中文 久久精品国产免费一区 久久无码精品一区二区三区 国产欧美另类久久久精品免费 欧美精品久久天天躁 亚洲精品在线视频 国产视频91在线 91精品福利一区二区三区野战 日韩中文字幕免费 国产精品99一区二区三区 欧美成人高清性色生活 国产精品系列在线观看 亚洲国产福利精品一区二区 国产成人在线小视频 国产精品久久久久免费 99re热在线视频 久久久久久久综合 一区二区国产在线播放 成人国产在线视频 亚洲精品乱码久久久久 欧美日韩一区二区综合 精品久久久久免费极品大片 中文字幕视频二区 激情粉嫩精品国产尤物 国产成人精品一区二区视频 久久精品中文字幕首页 亚洲高清在线 国产精品亚洲一区二区三区 伊人久久艹 中文在线亚洲 国产精品一区二区在线播放 国产精品九九免费视频 亚洲二区在线播放 亚洲狠狠婷婷综合久久久久网站 亚洲欧美日韩网站 日韩成人精品 亚洲国产一区二区三区青草影视 91精品国产福利在线观看 国产精品久久久久久久久99热 国产一区二区精品尤物 久碰香蕉精品视频在线观看 亚洲日日夜夜 在线不卡一区二区 国产午夜亚洲精品 九九热在线视频观看这里只有精品 伊人手机在线视频 91免费国产精品 日韩欧美中字 91精品国产91久久久久 国产全黄三级播放 视频一区二区三区免费观看 国产开裆丝袜高跟在线观看 国产成人欧美 激情综合丝袜美女一区二区 国产成人亚洲综合无 欧美精品一区二区三区免费观看 欧美亚洲国产日韩 日韩亚州 国产欧美日韩精品高清二区综合区 亚洲午夜国产片在线观看 精品久久久久久中文字幕 欧美精品1区 久久伊人久久亚洲综合 亚洲欧美日韩精品 国产成人精品久久亚洲高清不卡 久久福利影视 国产精品99精品久久免费 久久久久免费精品视频 国产日产亚洲精品 亚洲国产午夜电影在线入口 精品无码一区在线观看 午夜国产精品视频 亚洲一级片免费 伊人久久大香线蕉综合影 国产精品久久影院 久碰香蕉精品视频在线观看 www.欧美精品 在线小视频国产 亚洲国产天堂久久综合图区 欧美一区二区三区不卡 日韩美女福利视频 九九精品免视频国产成人 不卡国产00高中生在线视频 亚洲第一页在线视频 欧美日韩在线播放成人 99re视频这里只有精品 国产精品91在线 精品乱码一区二区三区在线 国产区久久 91麻豆精品国产自产在线观看一区 日韩精品成人在线 九九热在线观看 国产精品久久不卡日韩美女 欧美一区二区三区综合色视频 欧美精品免费一区欧美久久优播 国产精品网址 国产专区中文字幕 国产精品欧美亚洲韩国日本久久 日韩美香港a一级毛片 久久精品123 欧美一区二区三区免费看 99r在线视频 亚洲精品国产字幕久久vr 国产综合激情在线亚洲第一页 91免费国产精品 日韩免费小视频 亚洲国产精品综合一区在线 国产亚洲第一伦理第一区 在线亚洲精品 国产精品一区二区制服丝袜 国产在线成人精品 九九精品免视频国产成人 亚洲国产网 欧美日韩亚洲一区二区三区在线观看 在线亚洲精品 欧美一区二区三区高清视频 国产成人精品男人的天堂538 欧美日韩在线观看区一二 亚洲欧美一区二区久久 久久精品中文字幕首页 日本高清www午夜视频 久久精品国产免费 久久999精品 亚洲国产精品欧美综合 88国产精品视频一区二区三区 91久久偷偷做嫩草影院免费看 国产精品夜色视频一区二区 欧美日韩导航 国产成人啪精品午夜在线播放 一区二区视频在线免费观看 99久久精品国产自免费 精液呈暗黄色 久久99国产精品 日本精品久久久久久久久免费 精品国产97在线观看 99re视频这里只有精品 国产视频91在线 999av视频 亚洲美女视频一区二区三区 久久97久久97精品免视看 亚洲国产成人久久三区 99久久亚洲国产高清观看 日韩毛片在线视频 综合激情在线 91福利一区二区在线观看 一区二区视频在线免费观看 激情粉嫩精品国产尤物 国产成人精品曰本亚洲78 国产成人精品本亚洲 国产精品成人免费视频 国产成人啪精品视频免费软件 久久精品国产亚洲妲己影院 国产精品成人久久久久久久 久久大香线蕉综合爱 欧美一区二区三区高清视频 99热国产免费 在线观看欧美国产 91精品视频在线播放 国产精品福利社 欧美精品一区二区三区免费观看 国产一区二区免费视频 国产午夜精品一区二区 精品视频在线观看97 91精品福利久久久 国产一区福利 国产综合激情在线亚洲第一页 国产精品久久久久久久久久久不卡 九色国产 在线日韩国产 黄网在线观看 亚洲一区小说区中文字幕 中文字幕丝袜 日本二区在线观看 日本国产一区在线观看 欧美日韩一区二区三区久久 欧美精品亚洲精品日韩专 国产日产亚洲精品 久久综合九色综合欧美播 亚洲国产欧美无圣光一区 欧美视频区 亚洲乱码视频在线观看 久久无码精品一区二区三区 九九热精品免费视频 久久99精品久久久久久牛牛影视 国产精品成久久久久三级 国产一区福利 午夜国产精品视频 日本二区在线观看 99久久网站 国产亚洲天堂 精品国产一区二区三区不卡 亚洲国产日韩在线一区 国产成人综合在线观看网站 久久免费高清视频 欧美在线导航 午夜精品久久久久久99热7777 欧美久久综合网 国产小视频精品 国产尤物在线观看 亚洲国产精品综合一区在线 欧美一区二区三区不卡视频 欧美黑人欧美精品刺激 日本福利在线观看 久久国产偷 国产手机精品一区二区 国产热re99久久6国产精品 国产高清啪啪 欧美亚洲国产成人高清在线 国产在线第三页 亚洲综合一区二区三区 99r在线视频 99精品久久久久久久婷婷 国产精品乱码免费一区二区 国产在线精品福利91香蕉 国产尤物视频在线 五月婷婷亚洲 中文字幕久久综合伊人 亚洲精品一级毛片 99国产精品电影 在线视频第一页 久久99国产精品成人欧美 国产白白视频在线观看2 成人精品一区二区www 亚洲成人网在线观看 麻豆91在线视频 色综合合久久天天综合绕视看 久久精品国产免费高清 国产不卡一区二区三区免费视 欧美国产中文 99精品欧美 九九在线精品 国产中文字幕在线免费观看 国产一区中文字幕在线观看 国产成人一级 国产精品一区二区制服丝袜 国产一起色一起爱 亚洲精品成人在线 亚洲欧美精品在线 国产欧美自拍视频 99精品久久久久久久婷婷 久99视频 国产热re99久久6国产精品 视频一区亚洲 国产精品视频分类 国产精品成在线观看 99re6久精品国产首页 亚洲在成人网在线看 亚洲国产日韩在线一区 久久国产三级 日韩国产欧美 欧美在线一区二区三区 国产精品美女一级在线观看 成人午夜免费福利视频 亚洲天堂精品在线 91精品国产手机 欧美日韩视频在线播放 狠狠综合久久综合鬼色 九一色视频 青青视频国产 亚洲欧美自拍一区 中文字幕天天躁日日躁狠狠躁97 日韩免费大片 996热视频 伊人成综合网 亚洲天堂欧美 日韩精品亚洲人成在线观看 久久综合给会久久狠狠狠 日韩精品亚洲人成在线观看 日韩国产欧美 亚洲成aⅴ人片在线影院八 亚洲精品1区 99久久精品免费 国产精品高清在线观看 国产精品久久久免费视频 在线亚洲欧美日韩 91在线看视频 国产精品96久久久久久久 欧美日韩国产成人精品 91在线亚洲 热久久亚洲 国产精品美女免费视频观看 日韩在线毛片 亚洲永久免费视频 九九免费在线视频 亚洲一区网站 日本高清二区视频久二区 精品国产美女福利在线 伊人久久艹 国产精品久久久久三级 欧美成人精品第一区二区三区 99久久精品国产自免费 在线观看日韩一区 国产中文字幕一区 成人免费午夜视频 欧美日韩另类在线 久久99国产精品成人欧美 色婷婷中文网 久久天天躁夜夜躁狠狠躁2020 欧美成人伊人久久综合网 国产精品福利资源在线 国产伦精品一区二区三区高清 国产精品亚洲综合色区韩国 亚洲一区欧美日韩 色综合视频 国语自产精品视频在线区 国产高清a 成人国内精品久久久久影 国产在线精品香蕉综合网一区 国产不卡在线看 国产成人精品精品欧美 国产欧美日韩综合精品一区二区三区 韩国电影一区二区 国产在线视频www色 91中文字幕在线一区 国产人成午夜免视频网站 亚洲综合一区二区三区 色综合视频一区二区观看 久久五月网 九九热精品在线观看 国产一区二区三区国产精品 99久热re在线精品996热视频 亚洲国产网 在线视频亚洲一区 日韩字幕一中文在线综合 国产高清一级毛片在线不卡 精品国产色在线 国产高清视频一区二区 精品日本久久久久久久久久 亚洲国产午夜精品乱码 成人免费国产gav视频在线 日韩欧美一区二区在线观看 欧美曰批人成在线观看 韩国电影一区二区 99re这里只有精品6 日韩精品一区二区三区视频 99re6久精品国产首页 亚洲欧美一区二区三区导航 欧美色图一区二区三区 午夜精品视频在线观看 欧美激情在线观看一区二区三区 亚洲热在线 成人国产精品一区二区网站 亚洲一级毛片在线播放 亚洲一区小说区中文字幕 亚洲午夜久久久久影院 国产自产v一区二区三区c 国产精品视频免费 久久调教视频 国产成人91激情在线播放 国产精品欧美亚洲韩国日本久久 久久亚洲日本不卡一区二区 91中文字幕网 成人国产在线视频 国产视频91在线 欧美成人精品第一区二区三区 国产精品福利在线 久久综合九色综合精品 欧美一区二区三区精品 久久国产综合尤物免费观看 久久99青青久久99久久 日韩精品免费 久久国产精品999 91亚洲视频在线观看 国产精品igao视频 色综合区 在线亚洲欧国产精品专区 国产一区二区三区在线观看视频 亚洲精品成人在线 一区二区国产在线播放 中文在线亚洲 亚洲精品第一国产综合野 国产一区二区精品久久 一区二区三区四区精品视频 99热精品久久 中文字幕视频二区 国产成人精品男人的天堂538 99精品影视 美女福利视频一区二区 久久午夜夜伦伦鲁鲁片 综合久久久久久久综合网 国产精品国产欧美综合一区 国产99视频在线观看 国产亚洲女在线精品 婷婷影院在线综合免费视频 国产亚洲3p一区二区三区 91成人爽a毛片一区二区 亚洲一区二区高清 国产欧美亚洲精品第二区首页 欧美日韩导航 亚洲高清二区 欧美激情观看一区二区久久 日韩毛片在线播放 亚洲欧美日韩高清中文在线 亚洲日本在线播放 国产精品一区二区制服丝袜 精品国产一区二区三区不卡 国产不卡在线看 国产欧美网站 四虎永久在线观看视频精品 国产黄色片在线观看 夜夜综合 一本色道久久综合狠狠躁篇 欧美亚洲综合另类在线观看 国产91在线看 伊人久久国产 欧美一区二区在线观看免费网站 国产精品久久久久三级 久久福利 日韩中文字幕a 亚洲午夜久久久久影院 91在线高清视频 国产亚洲一区二区三区啪 久久人精品 国产精品亚洲午夜一区二区三区 综合久久久久久 久久伊人一区二区三区四区 国产综合久久久久久 日韩一区精品视频在线看 国产精品日韩欧美制服 日本精品1在线区 99re视频 无码av免费一区二区三区试看 国产视频1区 日韩欧美中文字幕一区 日本高清中文字幕一区二区三区a 亚洲国产欧美无圣光一区 国产在线视频一区二区三区 欧美国产第一页 在线亚洲欧美日韩 日韩中文字幕第一页 在线不卡一区二区 伊人久久青青 国产精品一区二区在线播放 www.五月婷婷 麻豆久久婷婷国产综合五月 亚洲精品区 久久国产欧美另类久久久 99在线视频免费 伊人久久中文字幕久久cm 久久精品成人免费看 久久这里只有精品首页 88国产精品视频一区二区三区 中文字幕日本在线mv视频精品 国产在线精品成人一区二区三区 伊人精品线视天天综合 亚洲一区二区黄色 国产尤物视频在线 亚洲精品99久久久久中文字幕 国产一区二区三区免费观看 伊人久久大香线蕉综合电影网 国产成人精品区在线观看 日本精品一区二区三区视频 日韩高清在线二区 久久免费播放视频 一区二区成人国产精品 国产精品免费精品自在线观看 亚洲精品视频二区 麻豆国产精品有码在线观看 精品日本一区二区 亚洲欧洲久久 久久中文字幕综合婷婷 中文字幕视频在线 国产成人精品综合在线观看 91精品国产91久久久久福利 精液呈暗黄色 香蕉国产综合久久猫咪 国产专区精品 亚洲精品无码不卡 国产永久视频 亚洲成a人片在线播放观看国产 一区二区国产在线播放 亚洲一区二区黄色 欧美日韩在线观看视频 亚洲精品另类 久久国产综合尤物免费观看 国产一区二区三区国产精品 高清视频一区 国产精品igao视频 国产精品资源在线 久久综合精品国产一区二区三区 www.五月婷婷 精品色综合 99热国产免费 麻豆福利影院 亚洲伊人久久大香线蕉苏妲己 久久电影院久久国产 久久精品伊人 在线日韩理论午夜中文电影 亚洲国产欧洲综合997久久 伊人国产精品 久草国产精品 欧美一区精品二区三区 亚洲成人高清在线 91免费国产精品 日韩精品福利在线 国产一线在线观看 国产不卡在线看 久久99青青久久99久久 亚洲精品亚洲人成在线播放 99久久免费看国产精品 国产日本在线观看 青草国产在线视频 麻豆久久婷婷国产综合五月 国产中文字幕一区 91久久精品国产性色也91久久 国产一区a 国产欧美日韩成人 国产亚洲女在线精品 一区二区美女 中文字幕在线2021一区 在线小视频国产 久久这里只有精品首页 国产在线第三页 欧美日韩中文字幕 在线亚洲+欧美+日本专区 精品国产一区二区三区不卡 久久这里精品 欧美在线va在线播放 精液呈暗黄色 91精品国产手机 91在线免费播放 欧美视频亚洲色图 欧美国产日韩精品 日韩高清不卡在线 精品视频免费观看 欧美日韩一区二区三区四区 国产欧美亚洲精品第二区首页 亚洲韩精品欧美一区二区三区 国产精品视频免费 在线精品小视频 久久午夜夜伦伦鲁鲁片 国产无套在线播放 久热这里只精品99re8久 欧美久久久久 久久香蕉国产线看观看精品蕉 国产成人精品男人的天堂538 亚洲人成网站色7799在线观看 日韩在线第二页 一本色道久久综合狠狠躁篇 国产一区二区三区不卡在线观看 亚洲乱码在线 在线观看欧美国产 久久福利青草精品资源站免费 国产玖玖在线观看 在线亚洲精品 亚洲成aⅴ人在线观看 精品91在线 欧美一区二三区 日韩中文字幕视频在线 日本成人一区二区 日韩免费专区 国内精品在线观看视频 久久国产综合尤物免费观看 国产精品系列在线观看 一本一道久久a久久精品综合 亚洲免费播放 久久精品国产免费 久久人精品 亚洲毛片网站 亚洲成a人一区二区三区 韩国福利一区二区三区高清视频 亚洲精品天堂在线 一区二区三区中文字幕 亚洲国产色婷婷精品综合在线观看 亚洲国产成人久久笫一页 999国产视频 国产精品香港三级在线电影 欧美日韩一区二区三区四区 日韩国产欧美 国产精品99一区二区三区 午夜国产精品理论片久久影院 亚洲精品中文字幕麻豆 亚洲国产高清视频 久久免费手机视频 日韩a在线观看 五月婷婷亚洲 亚洲精品中文字幕麻豆 中文字幕丝袜 www国产精品 亚洲天堂精品在线 亚洲乱码一区 国产日韩欧美三级 久久999精品 伊人热人久久中文字幕 久热国产在线视频 国产欧美日韩在线观看一区二区三区 国产一二三区在线 日韩国产欧美 91精品国产91久久久久 亚洲一区小说区中文字幕 精品一区二区免费视频 国产精品视频免费 国产精品亚洲综合色区韩国 亚洲国产精品成人午夜在线观看 欧美国产日韩精品 中文字幕精品一区二区精品