李志強(qiáng)， 倪 莉， 趙則祥， 張 洪， 梁 穎

(中原工學(xué)院 機(jī)電學(xué)院， 河南 鄭州 450007)

梳理機(jī)工作輥結(jié)合件是梳理機(jī)的關(guān)鍵件之一，2個工作輥的工作間隙理想狀況是保持一致性，從而保證梳理通過的纖維層的厚度一致，進(jìn)而達(dá)到高質(zhì)量的梳理效果，但實際狀況是，工作輥在重力載荷和針布載荷作用下產(chǎn)生撓度，導(dǎo)致工作輥間的隔距不均勻，對梳理機(jī)整體性能產(chǎn)生很大影響。目前國外同行業(yè)工作輥最大撓度值和間距調(diào)整方面做得較好。影響工作輥性能的指標(biāo)主要有材料和結(jié)構(gòu)2方面，因此，探索工作輥最佳結(jié)構(gòu)成為提高其性能的關(guān)鍵之一。作為結(jié)構(gòu)優(yōu)化概念設(shè)計層次的拓?fù)鋬?yōu)化，可從根本上定義結(jié)構(gòu)的形式，是真正意義上的最優(yōu)設(shè)計[1]。1988年，Bendsoe等基于“微結(jié)構(gòu)”思想提出了均勻化方法，由此拓?fù)鋬?yōu)化進(jìn)入繁榮時期。隨后發(fā)展的變密度法，使得拓?fù)鋬?yōu)化的工程應(yīng)用成為可能[2]。

國內(nèi)外學(xué)者在減少結(jié)構(gòu)撓度和結(jié)構(gòu)優(yōu)化方面做過很多研究：Yi等[3]提出如何在拓?fù)鋬?yōu)化的結(jié)果上建立參數(shù)化CAD模型；郭鵬等[4]推導(dǎo)了魚腹型中梁的優(yōu)化，降低了中梁撓度；王琴等[5]提出基于加載狀態(tài)下，動態(tài)特性預(yù)測的電主軸結(jié)構(gòu)參數(shù)優(yōu)化的思想；陳蘭等[6]對大型船用卷板機(jī)上的輥系統(tǒng)進(jìn)行力學(xué)分析和優(yōu)化，降低了上輥撓度；夏建芳等[7]以主梁幾種主要型鋼結(jié)構(gòu)參數(shù)為設(shè)計變量，運用ANSYS APDL語言編制了結(jié)構(gòu)優(yōu)化程序，主梁結(jié)構(gòu)優(yōu)化結(jié)果表明，其應(yīng)力分布更加合理，最大應(yīng)力下降4.8%，整體質(zhì)量下降16%，最大撓度上升50%，仍符合設(shè)計標(biāo)準(zhǔn)；李兆凱等[8]分析了二級層級褶皺結(jié)構(gòu)的板模型，在剪切或壓縮載荷下對應(yīng)的6種失效模式及對應(yīng)的名義應(yīng)力。通過比較名義應(yīng)力值，控制發(fā)生的失效模式類型,并結(jié)合多目標(biāo)優(yōu)化算法，同時進(jìn)行最小質(zhì)量和最小撓度的性能優(yōu)化；Yoshihiro等[9]研究了桁架在自身重力下的全局優(yōu)化；李志強(qiáng)等[10]研究了以位移約束、應(yīng)力約束、應(yīng)變能約束質(zhì)量最小化為目標(biāo)的拓?fù)鋬?yōu)化關(guān)系。通過對目標(biāo)函數(shù)與約束函數(shù)的解析敏度推導(dǎo),證明了他們之間的等價性。得到的準(zhǔn)則方程表明，在最優(yōu)結(jié)構(gòu)中,單元質(zhì)量與該單元應(yīng)變能之比等于結(jié)構(gòu)總質(zhì)量與結(jié)構(gòu)總應(yīng)變能之比；Zhu等[11]研究了動力響應(yīng)下的結(jié)構(gòu)優(yōu)化，提高了動力性能。

目前，國內(nèi)在工作輥撓度優(yōu)化效果上始終沒有達(dá)到國外先進(jìn)水平，為設(shè)計出高性能梳理機(jī)工作輥，借鑒前人研究優(yōu)化思路，本文提出采取二步法對工作輥結(jié)構(gòu)進(jìn)行優(yōu)化：首先，基于ANSYS APDL平臺應(yīng)用自編優(yōu)化程序?qū)ぷ鬏佭M(jìn)行拓?fù)鋬?yōu)化，得到工作輥合理拓?fù)錁?gòu)型；然后，基于該拓?fù)錁?gòu)型，對拓?fù)錁?gòu)型建立數(shù)學(xué)模型，再采用尺寸優(yōu)化方法對工作輥構(gòu)型做進(jìn)一步尺寸優(yōu)化，最終設(shè)計出滿足強(qiáng)度與撓度要求的高性能梳理機(jī)工作輥。

1 梳理機(jī)工作輥數(shù)學(xué)模型及優(yōu)化方法

1.1 工作輥拓?fù)鋬?yōu)化數(shù)學(xué)模型

工作輥拓?fù)鋬?yōu)化數(shù)學(xué)模型表達(dá)式為

(1)

式中：ρ為設(shè)計變量；ρi為連續(xù)體第i單元的偽密度；u(ρ)為最大位移；σj(ρ)和[σj]分別為第j個節(jié)點的應(yīng)力及其上限。

由于模型(1)的約束條件隨著有限元模型網(wǎng)格的增多，應(yīng)力對偽密度的敏度計算會急劇增加，因此，在優(yōu)化過程中會耗費大量機(jī)時，文獻(xiàn)[10]論證了位移、應(yīng)力、應(yīng)變能約束的等價性，因此，模型(1)可轉(zhuǎn)化為應(yīng)變能最小化問題，更有利于問題的解決，對應(yīng)的數(shù)學(xué)模型為

(2)

1.2 敏度分析

在優(yōu)化過程中采用理性數(shù)學(xué)準(zhǔn)則法需要進(jìn)行敏度分析，針對模型(2)的敏度推導(dǎo)如下

(3)

式中，針布力F2為定值，所以對偽密度的敏度為0。

1.3 工作輥尺寸優(yōu)化數(shù)學(xué)模型

在拓?fù)鋬?yōu)化的基礎(chǔ)上，建立結(jié)構(gòu)特征的數(shù)字參數(shù)化數(shù)學(xué)優(yōu)化模型：

(4)

式中：設(shè)計變量xi為工作輥各部分結(jié)構(gòu)尺寸；u(X)為最大位移；σj(X)和[σj]分別為第j個節(jié)點的應(yīng)力及其上限。

在尺寸優(yōu)化過程中，設(shè)計變量大大減少，從而敏度分析計算量也相應(yīng)減少，因此，通過目標(biāo)函數(shù)對設(shè)計變量的敏度分析推導(dǎo)出優(yōu)化準(zhǔn)則，然后基于ANSYS仿真平臺進(jìn)行優(yōu)化設(shè)計，計算耗費時間較少。

梳理機(jī)工作輥最大撓度最小化依據(jù)上述模型和設(shè)計變量對敏度的分析，就可以展開二步優(yōu)化法的實施：第1步是對工作輥的拓?fù)鋬?yōu)化；第2步是對工作輥的尺寸優(yōu)化。這2個階段的優(yōu)化目標(biāo)函數(shù)相同，設(shè)計變量不同，因此，優(yōu)化迭代準(zhǔn)則也不同，拓?fù)鋬?yōu)化階段的設(shè)計變量為單元的偽密度，尺寸優(yōu)化有限元模型(FEM)參數(shù)化建模中的設(shè)計變量為拓?fù)鋬?yōu)化結(jié)果中提取的幾何特征，在拓?fù)鋬?yōu)化結(jié)果向尺寸優(yōu)化有限元模型轉(zhuǎn)變過程中很關(guān)鍵，直接影響尺寸優(yōu)化結(jié)果。拓?fù)鋬?yōu)化結(jié)果的幾何特征提取應(yīng)該與優(yōu)化對象的生產(chǎn)工藝相聯(lián)系，避免設(shè)計出來而生產(chǎn)不出來的困境，也應(yīng)該與制造經(jīng)濟(jì)成本聯(lián)系起來，盡量做到結(jié)構(gòu)簡單、生產(chǎn)便捷、制造成本低、產(chǎn)品性能優(yōu)的效果。拓?fù)鋬?yōu)化結(jié)果只是提供了材料大體的分布情況，但不能代替第2步的尺寸優(yōu)化。工作輥優(yōu)化方案流程圖如圖1所示。

圖1 工作輥優(yōu)化基本流程Fig.1 Basic procedure of work roll optimization

2 算 例

梳理機(jī)工作輥軸承間距為3 891 mm，外徑為 220 mm，滾筒軸長為3 780 mm，軸頭材料為Q235-A圓鋼，滾筒材料為45鋼，彈性模量2.06×1011Pa，泊松比為0.3，密度為7.85 g/cm3，針布質(zhì)量為 30 kg。工作輥示意圖如圖2所示。

圖2 工作輥示意圖Fig.2 Sketch of work roll

2.1 工作輥構(gòu)型的拓?fù)鋬?yōu)化

優(yōu)化前，對國內(nèi)算例一指定的工作輥進(jìn)行試驗與仿真，試驗與仿真結(jié)果較吻合，工作輥最大撓度為0.3～0.4 mm，國外相同尺寸和材料的工作輥在相同工況下最大撓度稍大于0.1 mm。國內(nèi)外的區(qū)別在于工作輥結(jié)構(gòu)的差異，由此可見合理改變工作輥質(zhì)量分布可解決最大撓度最小化問題。拓?fù)鋬?yōu)化可打破人們以往根據(jù)經(jīng)驗預(yù)判提出最佳結(jié)構(gòu)，然后對提出結(jié)構(gòu)進(jìn)行尺寸優(yōu)化的不足，實現(xiàn)真正意義上的結(jié)構(gòu)優(yōu)化。根據(jù)拓?fù)鋬?yōu)化數(shù)學(xué)模型(2)對工作輥進(jìn)行拓?fù)鋬?yōu)化，結(jié)果如圖3所示。

圖3 工作輥簡化模型剖視圖Fig.3 Section view of work-roll simplified model

工作輥分為2個區(qū)域，拓?fù)鋬?yōu)化區(qū)域與非拓?fù)鋬?yōu)化區(qū)域：非拓?fù)鋬?yōu)化部分是為了保持工作輥功能特性的外形，不參與優(yōu)化；優(yōu)化區(qū)域是通過拓?fù)鋬?yōu)化對質(zhì)量進(jìn)行最佳分布，從而使工作輥最大撓度最小化。按照數(shù)學(xué)模型(2)，應(yīng)用ANSYS APDL平臺和自編程序?qū)ぷ鬏佭M(jìn)行拓?fù)鋬?yōu)化，找到材料最佳分布，并對拓?fù)錁?gòu)型結(jié)果提取數(shù)字特征，為第2步尺寸優(yōu)化建立參數(shù)化數(shù)學(xué)模型做準(zhǔn)備。

圖4示出工作輥的拓?fù)鋬?yōu)化過程中結(jié)構(gòu)材料隨著迭代步數(shù)增加，材料分布趨于最優(yōu)的進(jìn)化過程。圖中白色區(qū)域代表有材料分布，灰色區(qū)域代表工作輥沒有材料分布，白色區(qū)域外部的黑色區(qū)域為優(yōu)化過程中ANSYS窗口背景。圖4(a)～(f)中工作輥的壁厚從端部向中間逐漸變厚，從而增加了工作輥的剛度，工作輥縱截面類似魚骨型結(jié)構(gòu)，與肋板效果一樣也加強(qiáng)了工作輥剛度，但隨著迭代步數(shù)的增加，工作輥內(nèi)部魚骨狀結(jié)構(gòu)逐漸消失；在圖4(g)～(i)中隨著迭代的收斂，工作輥內(nèi)部魚骨狀結(jié)構(gòu)已經(jīng)完全消失，這是因為魚骨狀結(jié)構(gòu)一方面增強(qiáng)了工作輥剛度，有利于減少工作輥最大撓度，但另一方面其存在增加了工作輥質(zhì)量，這對減少工作輥最大撓度是不利的，這種兩難狀況依賴設(shè)計經(jīng)驗是無法完成的，從拓?fù)浣Y(jié)構(gòu)的漸變過程來看，利用拓?fù)鋬?yōu)化能清晰穩(wěn)定地得到材料的最佳分布。在優(yōu)化結(jié)束時，工作輥端部縱截面類似倒放著的M狀結(jié)構(gòu)得到保留，因為工作輥工作時上部承受壓應(yīng)力，下部承受拉應(yīng)力，這樣的結(jié)構(gòu)有利于減少縱向應(yīng)變，從而減少了工作輥的最大撓度。在拓?fù)鋬?yōu)化迭代過程中，材料分布漸變圖反映了迭代過程中材料分布變化的穩(wěn)定、清晰、快速的收斂過程，收斂的拓?fù)鋬?yōu)化結(jié)果圖(見圖4(i))為工作輥第2步尺寸優(yōu)化提供了結(jié)構(gòu)幾何特征，根據(jù)這些幾何特征，建立參數(shù)化工作輥幾何模型后就進(jìn)入了第2步優(yōu)化階段，即工作輥結(jié)構(gòu)尺寸優(yōu)化過程。

圖4 迭代過程中工作輥縱剖面結(jié)構(gòu)的演化過程Fig.4 Evolution of work roll rertical section structure in iterative process. (a) Iteration 5; (b) Iteration 6; (c) Iteration 7; (d) Iteration 8; (e) Iteration 9; (f) Iteration 10; (g) Iteration 11; (h) Iteration 12; (i) Iteration 13

2.2 工作輥的尺寸優(yōu)化

根據(jù)拓?fù)鋬?yōu)化結(jié)果對工作輥進(jìn)行參數(shù)化建模，如圖5所示。把優(yōu)化區(qū)域中的材料分布用參數(shù)化數(shù)學(xué)模型描述，中間部分可用分段線性函數(shù)描述，端部利用曲線擬合，具體選取的擬合函數(shù)與零部件生產(chǎn)水平與工藝相關(guān)，應(yīng)與相關(guān)企業(yè)協(xié)商，否則會出現(xiàn)設(shè)計出來制造不出來的狀況。算例中的設(shè)計變量有工作輥端面厚度、軸頸直徑、空腔兩端用圓臺擬合，中部用階梯錐孔擬合，工作輥壁厚等，從而得到參數(shù)化模型。

圖5 工作輥參數(shù)化幾何模型Fig.5 Work-roll geometric model with parameters

隨著3D打印的發(fā)展，制造業(yè)水平得到前所未有的提升，很多復(fù)雜結(jié)構(gòu)零件的生產(chǎn)工藝得以解決。在工作輥拓?fù)鋬?yōu)化結(jié)果的基礎(chǔ)上，提取數(shù)字特征建立參數(shù)化數(shù)學(xué)模型，進(jìn)入到工作輥尺寸優(yōu)化階段，按照數(shù)學(xué)模型(4)，基于ANSYS平臺，應(yīng)用自編程序?qū)ぷ鬏佭M(jìn)行結(jié)構(gòu)優(yōu)化，也可將模型(4)的問題轉(zhuǎn)化為應(yīng)變能約束下的質(zhì)量最小化問題。優(yōu)化過程中設(shè)計變量的選取很關(guān)鍵，直接影響優(yōu)化效果，這里軸頭直徑也作為設(shè)計變量。

圖6示出利用ANSYS軟件創(chuàng)建的工作輥幾何模型的縱剖面。對幾何模型劃分網(wǎng)格、加載、輸入利用敏度推導(dǎo)的優(yōu)化準(zhǔn)則進(jìn)行尺寸優(yōu)化，就得到工作輥最優(yōu)結(jié)構(gòu)。圖7示出工作輥尺寸優(yōu)化收斂時的一次優(yōu)化結(jié)果。位移云圖顯示，工作輥最大撓度為0.105 mm，約等于0.11 mm，工作輥最大撓度的最小值已經(jīng)與國外相同材料和工況的同類先進(jìn)工作輥的性能一致，優(yōu)化后的工作輥極大改善了工作輥撓度。

圖7 優(yōu)化后位移云圖Fig.7 Displacement of work-roll after optimization

圖6 工作輥幾何模型縱剖面Fig.6 Vertical section of work-roll geometric model

優(yōu)化前工作輥最大撓度為0.33 mm，最大應(yīng)力為7.12 MPa；優(yōu)化后其最大撓度約為 0.11 mm，最大應(yīng)力為9.15 MPa。說明優(yōu)化后的工作輥最大應(yīng)力有所增大，但遠(yuǎn)遠(yuǎn)低于許用應(yīng)力，最大撓度減少了67%。最大撓度發(fā)生在工作輥中間，中間壁厚的變化從直觀上看對撓度的影響是一個矛盾，壁厚加厚一方面提高了剛度，但另一方面也加大了自身質(zhì)量，這個問題無法從直觀經(jīng)驗上解決，正是某企業(yè)工作輥撓度減不下來的問題所在，工作輥端部構(gòu)形也是遠(yuǎn)非經(jīng)驗和想象可以設(shè)計出來的，通過拓?fù)鋬?yōu)化和尺寸優(yōu)化后工作輥結(jié)構(gòu)異于原結(jié)構(gòu)，材料分布更加合理，工作輥性能得到很大提高。圖8示出未優(yōu)化時原工作輥結(jié)構(gòu)圖。圖6近似認(rèn)為優(yōu)化后的工作輥結(jié)構(gòu)，與優(yōu)化前的工作輥區(qū)別是尺寸參數(shù)的變化，結(jié)構(gòu)形式不變。

圖8 原工作輥結(jié)構(gòu)圖Fig.8 Structure chart of original work-roll

3 結(jié)束語

本文研究了針布載荷和重力載荷下梳理機(jī)工作輥結(jié)構(gòu)的拓?fù)鋬?yōu)化和尺寸優(yōu)化，工作輥結(jié)構(gòu)優(yōu)化過程分為2個步驟進(jìn)行：首先對工作輥劃定優(yōu)化區(qū)域與非優(yōu)化區(qū)域，按照給定的數(shù)學(xué)模型對優(yōu)化區(qū)域進(jìn)行拓?fù)鋬?yōu)化，找到優(yōu)化區(qū)域最佳材料分布；然后在前一步優(yōu)化結(jié)果的基礎(chǔ)上，抽取材料分布特征，建立參數(shù)化數(shù)學(xué)模型，選定優(yōu)化方法?；贏NSYS軟件對工作輥進(jìn)行尺寸優(yōu)化，獲得工作輥最佳結(jié)構(gòu)，通過算例仿真對比，驗證了優(yōu)化方法有效性。

優(yōu)化過程中，采用理性數(shù)學(xué)準(zhǔn)則使拓?fù)浣Y(jié)果更加精確可靠，為工程人員提供了新的設(shè)計方法，具有一定的工程應(yīng)用價值。由于與企業(yè)合作過程中未提到工作輥的固有頻率限制，在相應(yīng)的固有頻率約束下，工作輥撓度最小化需要進(jìn)一步研究。

FZXB