熱連軋輥系變形三維建模及有限元分析

1、<p>　　需要圖紙聯(lián)系QQ97666224</p><p><b>　　1、 緒論</b></p><p>　　1.1 四輥軋機發(fā)展情況概論</p><p>　　近年來我國軋鋼行業(yè)得到了飛速發(fā)展，鋼材年產量突破了2億噸，已連續(xù)多年成為世界鋼產量第一大國。板帶材的軋制生產能力逐步提升到了一個較高的水平，各種板帶產品也得以廣泛的應用

2、于生產和生活中的方方面面。但是我國目前軋鋼生產的技術水平與國際先進水平相比還有相當大的差距，軋制產品的主要技術指標與國際先進水平相比仍有相當大的差距，我國已經入世，國外鋼材生產技術強國的行業(yè)沖擊愈發(fā)明顯起來，要想在空前激烈的競爭中得以生存、獲得發(fā)展，我們就必須在軋機精度控制等方面多做工作。</p><p>　　四輥軋機以其較高的生產能力和良好的產品質量廣泛應用于板帶生產中，近年來隨著國民經濟的不斷發(fā)展以及工業(yè)生產

3、需求的不斷增長，用戶對板帶產品的平直度等指標要求越來越高，這就對板帶軋制中輥縫的控制精度提出了更高的要求。對四輥軋機輥系變形進行分析，是關乎板帶材質量的決定性因素。如何提高軋機輥系變形分析的水平，對各個工廠來說是要亟待解決的，傳統(tǒng)的分析方法，繁雜且精度不高。本課題采用基于ANSYS軟件的有限元分析法對四輥軋機輥系變形進行研究，是近年來一種正在被逐步廣泛應用的方法。</p><p>　　1.2 輥系變形計算的常用理

4、論與計算方法</p><p>　　1.2.1軋輥變形模型的分類</p><p>　　關于板形的軋輥變形模型的研究發(fā)展可追溯到1958年，那時薩克斯爾(Saxl)第一次對四輥軋機做了全面深入的研究。此后由于引進了數學模型，這一領域得到了更進一步的拓展。這些模型的分類如下：</p><p>　　(1) 二輥軋機的簡支梁模型；</p><p>　

5、　(2) 四輥軋機的簡支梁模型；</p><p>　　(3) 分割梁模型；</p><p>　　(4) 有限元分析模型。</p><p>　　1.2.2 二輥軋機的簡支梁模型</p><p>　　在二輥軋機簡支梁模型中，將工作輥視為線彈性應力梁。在推導梁的撓曲公式時，我們做了以下假定：</p><p>　　梁的材

6、質均勻，在拉伸與壓縮時的彈性模量相同；</p><p><b>　　梁的橫斷面相同;</b></p><p>　　梁至少關于一個軸向平面對稱；</p><p>　　所有的加載和反作用力都與梁的軸線垂直；</p><p>　　對于具有緊湊斷面的金屬梁，其寬高比等于或大于8。 </p><p>　

7、　板帶材的板形可以通過對以下的軋輥的兩類撓曲進行疊加來確定：</p><p>　　由于軋制力引起的彎曲力使軋輥產生的撓曲；</p><p>　　由于軋制力引起的剪切力使軋輥產生的撓曲。</p><p>　　由彎曲力產生的撓曲可由如下的微分方程描述：</p><p><b>　　(1-1)</b></p>&

8、lt;p>　　式中 ——軋輥彈性模量；</p><p>　　——在距離x處軋輥斷面的慣性矩(圖1-1)</p><p><b>　　P—— 軋制力；</b></p><p>　　——在距離為x處的軋輥撓度；</p><p>　　L—— 軋制力作用點的間距；</p><p>　　——

9、帶材的寬度。 </p><p>　　在軋輥與帶材的接觸區(qū)中，x的變化范圍為：</p><p><b>　　(1-2)</b></p><p>　　在這一范圍內，方程1-1的解適用于二輥軋機。此解由拉克(Larke)給出如下：</p><p>　　(1-3)

10、 </p><p><b>　　圖1-1</b></p><p><b>　　其中：</b></p><p>　　因剪切力產生的軋輥撓度由拉克計算得出，這一結果是在假定軋輥垂直斷面上的剪切應力呈均勻分布條件下獲得的。軋輥和帶材

11、的接觸區(qū)內的軋輥撓度可由如下微分方程給出：</p><p><b>　　(1-4)</b></p><p>　　式中： ——軋輥彈性剪切模量；</p><p>　　——在x處的軋輥撓度。</p><p>　　拉克給出了方程1-4的解如下：</p><p><b>　　(1-5)<

12、/b></p><p>　　由彎曲力和剪切力產生的軋輥總的撓度為：</p><p><b>　　(1-6)</b></p><p>　　1.2.3 四輥軋機的簡支梁模型</p><p>　　在四輥軋機的簡直梁模型中，工作輥和支撐輥都被認為是完全彈性應力梁。</p><p>　　斯通(Sto

13、ne)和戈雷(Gray)采用的模型中，四輥軋機的軋輥撓曲可以看成一個置于彈性基礎上的簡支梁的撓曲情況，如圖1-2所示。鐵木辛克(Timoshenko)推導了作為這一模型依據的微分方程：</p><p><b>　　(1-7)</b></p><p>　　式中 ——工作輥的彈性模量；</p><p>　　——工作輥的慣性矩；</p>

14、;<p>　　—— 工作輥在x處的撓度。</p><p>　　此撓度曲線方程的通解為：</p><p><b>　　(1-8)</b></p><p>　　錯誤!未指定書簽。錯誤!未指定書簽。</p><p>　　式中 A、B、C、D——積分常量，取決于載荷類型和邊界條件。</p>

15、<p><b>　　圖1-2</b></p><p>　　彈性量的邊界條件和載荷情況如圖1-2所示，其中k1為代表支撐輥和相鄰工作輥之間彈性條件的簡支梁常數，k2是代表相互接觸的工作輥和軋材彈性的常數。根據斯通模型：</p><p>　　k= k1+k2, 其中k1和k2的值為：</p><p>　　i=1,2

16、 (1-9)</p><p>　　式中 ——兩接觸軋輥沿水平軸方向的接觸長度；</p><p>　　——兩接觸軋輥中心線徑向接近量。</p><p>　　根據弗普爾(Foppl)的研究，當支撐輥的彈性模量與工作輥的相等，即= =時，一對軋輥的中心線徑向接近量可由如下的方程確定：</p><p><b>　　

17、(1-10)</b></p><p>　　式中 、——在壓力狀態(tài)下的工作輥和支撐輥直徑；</p><p>　　b——工作輥和支撐輥的壓扁接觸寬度，等于：</p><p><b>　　(1-11)</b></p><p>　　式中 ——軋輥材質的泊松比。</p><p><

18、b>　　簡支梁模型的局限性</b></p><p>　　簡支梁撓曲模型考慮了許多影響板形的重要因素，但這些模型仍然存在著以下的不足：</p><p>　　通常軋機支撐輥的徑長比遠遠小于8，而8時從上述方程獲得精確結果所要求的最小值。</p><p>　　針對二輥軋機推導出的方程1-1和1-2的解若要用于四輥軋機，應假定在四輥軋機中，工作輥與支撐輥間

19、的壓力是沿帶材寬度傳遞的。而事實上，此壓力的傳遞是通過整個接觸區(qū)進行的。</p><p>　　簡支梁模型不能模擬在工作輥和軋件及工作輥與支撐輥間載荷橫向不均勻分布，因此，該模型沒有考慮一些重要的因素，如軋輥凸度、隨后的板凸度、材質硬度沿帶寬方向的分布、軋輥的磨損等。</p><p>　　上述局限性使得我們有理由去研究發(fā)展更為復雜的模型來計算板形。潑普勞斯基(Poplawski)和麥克迪爾摩

20、特(McDeermott)提出了一種可以模擬軋輥和板帶凸度，以及沿軋輥長度方向的載和 不均勻分布影響的方法。在該模型中，工作輥和 圖1-3 </p><p>　　支撐輥的接觸面積工作輥和帶材的接觸面杯模擬</p><p>　　成一系列的彈簧，但軋輥仍被視為簡支的彈性梁。</p>

21、<p>　　在由王國棟等提出的模型中，將簡支梁布置換成了由兩個懸臂梁代表輥身 的布置形式，如圖2-3所示，對于采用了橫移軋輥的四輥和六輥軋機，計算軋輥的彈性變形和出口板凸度時，采用了矩陣的方法。 </p><p>　　1.2.4分割法模型

22、 </p><p>　　在分割法中梁的撓曲模型最先由紹特(Shohet)和湯森德(TowpseHd)提出，在此模型中，板凸度通過確定以下三個未知量來求出：</p><p>　　(1)工作輥和板帶問的橫向載荷分布；</p><p>　　(2)工作輥和支承輥間的橫向載荷分布；</p><p>　　(3)工作輥的剛體移動。</p>

23、<p>　　因為軋機是關于軋輥中心對稱的，所以計算過程僅需考慮軋輥的一半。所用數值方法是將軋輥切分成m個單元，如圖1-4所示，并將軋輥的分布載荷代之以施加在每個單元中心的集中載荷。因為軋件寬度小于輥身長度，所以與軋件接觸的軋輥單元數n小于m。</p><p>　　支撐輥的變形。在分割法模型中，支撐輥的變形可以表示為在位置i處的軋輥表面的垂直位移：</p><p><b>

24、;　　(1-12)</b></p><p>　　式中 ——工作輥和支撐輥間第j個單元每單位寬度上的載荷；</p><p><b>　　——影響系數；</b></p><p><b>　　——單元的寬度；</b></p><p>　　——在i處的支撐輥表面的局部接觸變形；</p&

25、gt;<p>　　——支撐輥的剛體移動。</p><p>　　工作輥的變形工作輥變形可表示為工作輥表面在位全i處的垂直位移：</p><p><b>　　(1-13)</b></p><p>　　式中 ——工作輥和軋件間第i個單元上每單位寬度的載荷；</p><p>　　——工作輥的影響系數；</

26、p><p>　　——I點處的系數因子(在有軋件處,=1,在無軋件處,=0)；</p><p>　　——工作輥表面在i處的接觸變形；</p><p><b>　　——工作輥的剛度。</b></p><p>　　影響系數的計算。影響系數是由第j個單元重點的載荷在第i個單元中點產生的撓曲(圖1-4)。

27、 </p><p>　　圖1-4 </p><p>　　當時：

28、 </p><p><b>　　(1-14)</b></p><p><b>　　當時：</b></p><p><b>　　(1-15)</b></p><p>　　將工作輥和支撐輥的參數對應值代入方程1-14和1-15，可得系數和的值。</p>&

29、lt;p>　　工作輥和支撐輥接觸匹配關系。工作輥和支撐輥接觸匹配可在假定無載荷時得出，此時凸工作輥和支撐輥只是點接觸，在此點以外存在輥縫γ(i)。在有軋值力作用時，此輥縫會減小。那么，相互接觸的工作輥和支撐輥的撓度換算方程為：</p><p><b>　　(1-16)</b></p><p>　　工作輥和軋件接觸匹配關系。工作輥和軋件接觸匹配考慮了受負載軋輥任意

30、點的輥縫高度，還有軋件在相應點的出口厚度隨著軋輥的彈性壓扁和工作輥撓區(qū)的總和而變化。由此，其協(xié)調方程為：</p><p><b>　　(1-17)</b></p><p>　　式中 ——在i點處軋件出口厚度的一半；</p><p>　　——在i點處空載輥縫的高度。</p><p>　　軋簡的出口厚度可以采用軋值理論

31、的線性方程計算出來，其中是軋件入口厚度、軋件力、帶鋼張力和變形抗力的函數。</p><p>　　(6)工作輥的靜平衡。工作輥的靜平衡時將軋輥間的垂直載荷、軋件和工作輥間的載荷及彎輥設備施加給工作輥的載荷進行疊加，則方程為：</p><p>　　(1-18) </p><p>　　式中

32、 ——總的軋制彎曲力(見圖1-4)。</p><p>　　(7)方程的解。求解方程1-16、1-17、1-18在于如何找出m+n+1個未知量，即表示工作輥和支撐輥之間載荷分布的m個力的值(),表示工作輥和軋件間載荷分布的n個力的值()和工作輥的剛度系數。有了上述三個基本方程，可以得到求解、和的m+n+1個方程。這些方程可以用代數矩陣的方法求解。</p><p>　　分割梁撓曲模型還考慮了帶

33、材沿寬度方向的張力，這里另外增加了n個未知量。這一問題可以通過疊加的方式來解決。 圖1-5 </p><p>　　在由郭任明(Guo)研制的模型中，將軋輥劃分為一系列彈簧單元(圖1-5)。這些彈簧變形的相互關系通過線性相關方法計算出來。這種方法考慮了由彎矩和剪切力產生的梁的撓曲變形。這一模型也將帶材視為一系列的彈簧單元，并認為這些彈簧單元的剛度是帶材軋制特性的函數?；籼m德

34、(Hollander)和萊茵(Reinen)曾經對另一種模型做了闡述。</p><p>　　分割梁模型的局限性 </p><p>　　分割梁撓曲模型的研究發(fā)展對于在軋制過程中提高帶材板形的模擬能力邁出了很大一步。然而，這類模型也有自身的不足，因為分割梁模型是建立在假設在軋制力作用下工作輥和支承輥完全接觸之上的。實際上，當采用特殊的輥型系統(tǒng)如CVC軋輥、

35、UPC軋輥和錐形軋輥時，情況并非如此。此時還需考慮可能存在的接觸面不吻合(如圖1-6)。</p><p>　　此模型的另一不足是該模型對影響系數的計算是根據簡支架撓曲方程而來的。但是，如前所述當軋輥的徑長比小時，這些方程的實用性是值得懷疑的，而且分割梁模型用二維問題代替三維 圖1-6</p><p>　　問題，因此在某些情況下，不可能獲得良好的計算精度。

36、 </p><p>　　1.2.5 有限元分析理論 </p><p>　　在有限元分析中采用的是矩陣結構分析方法。這種方法采用了一種直接的物理方法來建立和求解梁及框架結構問題。</p><p&

37、gt;　　求解工程問題的解要用到以下三個條件：</p><p><b>　　力平衡方程；</b></p><p><b>　　變形協(xié)調方程；</b></p><p>　　材料行為的本構關系。</p><p>　　利用這三個條件可以建立未知應力(力法)或未知位移(位移法)組成的方程。在有限元技術中，常

38、用位移法。</p><p>　　有限元分析模型是對如下一些連續(xù)物體問題進行近似求解的一種方法：</p><p>　　連續(xù)物體被劃分成有限個單元，各個單元的行為由有限個參數給定；</p><p>　　整個系統(tǒng)用它的單元集合體求解同樣精確地遵循適用于標準離散問題的那些原則。</p><p>　　由有限個單元構成的集合體通常是指網格，這些單元通過節(jié)

39、點相互連接起來。在平面問題中，單元可以是三角形或是四邊形的。在三維問題中單元可以是三棱柱或長方體或六棱柱。</p><p>　　有限元分析采用了多年發(fā)展起來的應用于離散問題的標準方法。這套方法包含了對物體或結構的每一個單元的力-位移關系的計算。按既定的步驟，對每個有限元的每個節(jié)點都能建立起局部的力平衡方程來。求解這些方程可得到未知位移的解。有限元分析方法一般分為以下六個步驟：</p><p&g

40、t;　　(1)選擇位移模型。位移模型將各個單元的位移表示成各個節(jié)點位移的函數，通常表示為：</p><p><b>　　(1-19)</b></p><p>　　式中 ——單元的位移矩陣；</p><p>　　——單元的形函數矩陣；</p><p>　　——單元的節(jié)點位移矩陣(未知的變量)。</p>

41、<p>　　建立位移模型時，沿給定方向的位移分布，通常用一個簡單的函數表示，例如下面的多項式：</p><p><b>　　(1-20)</b></p><p>　　式中 ——沿x方向的位移。</p><p>　　多項式的系數被稱為廣義位移系數，他們確定了位移模型的形狀。通過邊界條件可求得上述系數。有限元方程的個數依賴于要被模型化的

42、結構的幾何形狀及模型可能移動的方向的個數，即自由度。</p><p>　　簡單的彈簧單元可承受張力并有兩個自由度，因為梁的兩端可以沿彈簧的主軸方向自由移動。然而，具有更多自由度的單元體可以作其他運動如彎曲運動。單元類型的選擇也決定了在一個模型中的自由度。</p><p>　　(2) 應變-位移關系的建立。應變-位移關系可通常表示為：</p><p><b>

43、;　　(1-21)</b></p><p>　　式中 ——單元應變矩陣；</p><p>　　——單元應變-位移關系矩陣。</p><p>　　應變-應變關系的建立。應變-應變關系可表示為：</p><p>　　(1-22) </p><p>　　式中 ——單元應力矩陣；<

44、/p><p>　　——單元材料的應力-應變關系矩陣。</p><p>　　在彈性變形時，材料的應力—應變關系矩陣[C]常根據虎克(Hooke)定律得到。這一矩陣建立起了線應變、剪應變、正應力和剪應力之間的關系(如圖1-7所示)：</p><p><b>　　(1-23)</b></p><p>　　矩陣[C]中的系數常用楊氏

45、模量E和泊松比表示。</p><p>　　(4) 建立有限元剛度矩陣。有限元剛度將節(jié)點位移與節(jié)點處的力聯(lián)系起來，即：</p><p><b>　　(1-24)</b></p><p>　　式中 ——節(jié)點處的力矢量(已知量)；</p><p>　　——單元剛度矩陣(剛度影響系數)。</p>&l

46、t;p>　　剛度矩陣可由最小余能原理求得，通常表示為：</p><p><b>　　(1-25)</b></p><p>　　式中 ——[B]的轉置矩陣；</p><p><b>　　V——單元的體積。</b></p><p>　　(5)聯(lián)立代數方程。這一步包括對由每個單元的剛度矩陣[k

47、]構成的總剛度矩陣[K]和由各個節(jié)點的力矢量{P}構成的總的力矢量{R}的聯(lián)立。</p><p>　　聯(lián)立方法的根據是由于每個節(jié)點的相互聯(lián)系要求所有于該節(jié)點相鄰的節(jié)點在該節(jié)點位移應相同。總的剛度矩陣、總的力矢量和總的位移矢量{r}的平衡關系可有一套聯(lián)立方程表示：</p><p><b>　　(1-26)</b></p><p>　　(6)求解未

48、知參數?？紤]物體的幾何和力的邊界條件，求解代數方程1-16，可得所有未知的位移。在線性平衡問題中，可直接應用矩陣代數技術來求解：</p><p><b>　　(1-27)</b></p><p>　　1.2.6 二維有限元分析模型</p><p>　　圖1-7時一個典型的四輥軋機軋輥系統(tǒng)的二維有限元網格圖，位于X-Y平面內。在最簡單的情況中，Z

49、方向的厚度是一個常量，這很大程度上減少了模型中有限元的數量，使我們無需額外增加大量的計算時間就可設計出一個十分準確的有限元網格。但是，只用二維作分析而忽略第三維的影響，我們很難得到高精度的結果來。</p><p>　　這一問題在一定程度上被陳先霖(XianLin)和鄒家祥(Jiaxiang)解 圖1-7</p><p>　　決了,他們采用了一種在

50、Z方向上具有不同厚度的二維模型。在這一模型中，每個單 </p><p>　　元根據其到軋輥軸線的縱坐標軸距離和軋輥的半徑對應不同的厚度，第i個單元對應的厚度為：</p><p><b>　　(1-28)</b></p><p>　　式中 R——軋輥半徑。<

51、;/p><p>　　參數和由下面的式子定義：</p><p><b>　　(1-29)</b></p><p>　　式中 ——第i個單元的縱坐標；</p><p>　　——第i個單元的高度。</p><p>　　盡管二維模型只是物理模型的一種簡化表示，當其在分析許多因素如軋輥壓扁和帶材張力對板形的影響

52、時仍不失為一個極其方便的工具。</p><p>　　1.2.7 三維有限元分析模型</p><p>　　三維有限元網格為一個所研究系統(tǒng)的物理模型提供了最確切的表示。但是，在確定網格中有限元的數量和類型時要費點勁。</p><p>　　圖1-8描述了由一對軋輥和軋制帶材構成的體系的三維有限元網格，曾應用于由聯(lián)合工程公司和國際軋鋼咨詢公司聯(lián)合研制的ROLL-FLEX離線

53、模型中。</p><p>　　在研究網格時，單元的數量或向格的尺寸可由所要求的板凸度計算精度而定。但是，選擇單元的數量還要考慮計算時間和花費的成本。在評價計算精度時，要用到如下的參數：</p><p>　　(1)中心凸度的相對誤差；</p><p>　　(2)軋輥壓扁影響的相對誤差。</p><p>　　軋輥中心凸度的相對誤差可表示為：<

54、;/p><p>　　(1-30) </p><p>　　圖1-8

55、 </p><p>　　式中 、——在有限元網格中，當網格數為n和∞時分別對應的中心凸度；</p><p>　　——帶材的出口厚度。</p><p>　　因為中心凸度的相對誤差隨著網格數量的增加而減小，所以我們可以根據指數定律來確定的值。</p>

56、;<p>　　如圖1-9所示，隨著軋材剛度的增加，為了達到相同的精度，有限元網格的數量也增加了。當中心凸度的相對誤差為0.25%時，表示F1、F4、F6機架的網格數分別減少為2500、5000和7000個。軋輥壓扁影響的相對誤差由如下方程計算得出：</p><p><b>　　(1-31)</b></p><p>　　式中 ——有限元分析中有n個網

57、格時工作輥和支撐輥的中心線徑向接近量；</p><p>　　——由軋輥壓扁方程計算得到的工作輥和支撐輥的中心線徑向接近量。</p><p>　　根據弗浦爾的研究，當工作輥和支撐輥具有相同的彈性模量E時，一對軋輥如工作輥和支撐輥的中心接近量 可由方程1-10和1-11確定。若要將F1、F4、F6機架的軋輥的壓扁影響的相對誤差減小到0.25%,那么所要求的有限元網格數量必須分別等于或大于350

58、0、4600和5200個，如圖1-10所示。</p><p><b>　　圖1-9</b></p><p><b>　　圖1-10 </b></p><p>　　1.3 ANSYS軟件</p><p>　　1.3.1 ANSYS的工作原理——有限單元法</p><p>　　

59、有限單元法的基本思想是將連續(xù)的結構離散成有限個單元，并在每一個單元中設定有限個節(jié)點，將連續(xù)體看作是只在節(jié)點處相連接的一組單元的集合體；同時選定場函數的節(jié)點值作為基本未知量，并在每一單元中假設一近似差值函數以表示單元中場函數的分布規(guī)律；進而利用力學中的某些變分原理去建立用于求解節(jié)點未知量的有限元法方程，從而將一個連續(xù)域中的無限自由度問題化為離散域中的有限自由度問題。一經求解就可以利用解得的節(jié)點值和設定的插值函數確定單元上以至整個集合體上的

60、場函數。有限元求解程序的內部過程可從圖1-11中看出。</p><p><b>　　↓</b></p><p><b>　　↓</b></p><p><b>　　↓</b></p><p><b>　　↓</b></p><p>

61、<b>　　↓</b></p><p><b>　　↓</b></p><p><b>　　圖1-11</b></p><p>　　近40年來，隨著計算機的飛速發(fā)展和廣泛應用，各種行之有效的數值計算方法得到了巨大的發(fā)展。而有限元方法則是計算機誕生以后，在計算數學、計算力學、和計算工程科學領域里誕生的最

62、有效的計算方法。隨著有限元理論基礎的日益完善，出現(xiàn)了很多通用和專用的有限元計算軟件。在國際畫的市場經濟中，企業(yè)間的競爭日益加劇。為取得競爭優(yōu)勢，企業(yè)迫切需要以高質量低成本的產品迅速搶占市場，因此企業(yè)迫切需要高技術、高速度、低成本的設計方法。ANSYS程序即是應此要求而發(fā)展起來的計算機仿真設計工具。</p><p>　　1.3.2 ANSYS軟件簡介</p><p>　　ANSYS軟件由成立

63、于1970年的美國ANSYS公司開發(fā)完成，是融結構、熱、流體、電磁、聲學于一體的大型通用有限元分析軟件，可廣泛應用于核工業(yè)、鐵道、石油化工、航空航天、機械制造、能源、汽車交通、國防軍工、電子、土木工程、造船、生物醫(yī)學、輕工、地礦、水利、日用家電等一般工業(yè)及科學研究。</p><p>　　ANSYS程序是一個功能強大的、靈活的設計分析及優(yōu)化軟件包。該軟件可浮動運行于從PC機、NT工作站、UNIX工作站甚至巨型機的各

64、類計算機及操作系統(tǒng)中，數據文件在其所有的產品系列和工作品臺上均兼容。其多物理場耦合的功能，允許在同一模型上進行各式各樣的耦合計算，如：熱—結構耦合、磁—結構耦合、以及電—磁—流體—熱耦合。在PC機上生成的模型同樣可以運行于巨型機上，這樣就保證了所有的ANSYS用戶的多領域多變工程的求解。</p><p>　　1.3.2.1 ANSYS軟件的組成</p><p>　　ANSYS軟件 主要包括

65、三個部分：前處理模塊、分析計算模塊和后處理模塊。</p><p><b>　　前處理模塊</b></p><p>　　它為用戶提供了一個強大的實體建模及網絡劃分工具，用戶可以方便地構造有限元模型，軟件提供了100種以上的單元類型，用來模擬工程中的各種結構和材料。</p><p><b>　　實體建模</b></p&g

66、t;<p><b>　　參數化建模</b></p><p><b>　　體素庫及布爾運算</b></p><p>　　拖拉、旋轉、拷貝、蒙皮、倒角等。</p><p>　　多種網格自動劃分工具，自動進行單元形態(tài)、求解精度檢查及修正。</p><p>　　自動/映射網格劃分、智能網格劃分

67、、自適應網格劃分。</p><p>　　復雜幾何體Sweep映射網格生成。</p><p>　　六面體向四面體自動過渡網格：金字塔形。</p><p><b>　　邊界層網格劃分</b></p><p>　　在幾何模型或FE模型上加載：點載荷、分布載荷、體載荷、函數載荷。</p><p>　　可擴

68、展的標準梁截面形狀庫</p><p><b>　　分析計算模塊</b></p><p>　　包括結構分析（可進行線性分析、非線性分析和高度非線性分析）、流體動力學分析、電磁場分析、聲場分析、壓電分析及多物理場的耦合分析，可模擬多種物理介質的相互作用，具有靈敏度分析及優(yōu)化分析能力。</p><p><b>　　后處理模塊</b&g

69、t;</p><p>　　可將計算結果以彩色等值線顯示、梯度顯示、矢量顯示、力子流跡顯示、立體切片顯示、透明及半透明顯示（可看到內部結構）等圖形方式顯示出來，也可以將計算結果以圖表、曲線形式顯示或輸出。具體如下所示：</p><p>　　計算報告自動生成及定制工具，自動生成符合要求格式的計算報告。</p><p>　　結果顯示菜單：圖形顯示、抓圖、結果列表。<

70、/p><p>　　圖形：云圖、等值線、矢量顯示、粒子流跡顯示、切片、透明及半透明顯示、紋理。</p><p>　　鋼筋混凝土單元可顯示單元內的鋼筋、開裂情況以及壓碎部位。</p><p>　　梁、管、板、復合材料單元及結果按實際形狀顯示，顯示橫截面結果；顯示梁單元彎矩圖。</p><p>　　顯示優(yōu)化靈敏度及優(yōu)化變量曲線。</p>

71、<p>　　各種結果動畫顯示，可獨立保存及重放。</p><p><b>　　3D圖形注釋功能。</b></p><p>　　直接生成BMP、JPG、VRML、WMF、PNG、PS、TIFF、HPGL等格式的圖形。</p><p>　　計算結果排序、檢索、列表及再組合。</p><p>　　提供對計算結果的加、

72、減、積分、微分等計算。</p><p>　　顯示沿任意路徑的結果曲線，并可進行路徑的數學計算。</p><p>　　1.3.2.2 ANSYS主要的技術特點</p><p>　　● 唯一能實現(xiàn)多場及多耦合分析的軟件</p><p>　　唯一實現(xiàn)前后處理、求解及多場分析統(tǒng)一數據庫的一體化大型FEA分析軟件</p><p>

73、;　　唯一具有物理場優(yōu)化功能的FEA軟件</p><p>　　唯一具有中文界面的大型通用有限元分析軟件</p><p>　　具有強大的非線性分析功能</p><p>　　具有適用于不同的問題和硬件配置的多種求解器</p><p>　　支持異種異構功能網絡浮動，在異種、異構平臺上支持界面統(tǒng)一，數據文件通用</p><p>

74、;　　強大的并行計算功能，支持分布式并行和共享內存式并行</p><p>　　多種用戶網格劃分技術</p><p><b>　　完善的用戶開發(fā)環(huán)境</b></p><p>　　同時，ANSYS軟件擁有豐富和完善的單元庫、材料模型庫和求解器，保證了他能高效的求解各類結構的靜力、動力、振動、線性和非線性問題，壓縮和不可壓縮的流體問題。其友好的圖形界

75、面和程序結構，交互式的前后處理和圖形軟件，大大的減輕了用戶在實際工程問題中創(chuàng)建模型、有限元求解以及結果分析和評價的工作量。他的統(tǒng)一集中式的數據庫保證了個模塊之間的有效可靠的集成，并實現(xiàn)了與多個CAD/CAE軟件的友好鏈接。</p><p>　　1.4 本文主要研究內容及創(chuàng)新</p><p>　　在深入理解有限元分析實質以及ANSYS應用的基礎上，即可對四輥軋機工作輥的變形進行計算。本課題要

76、求在給定輥系的幾何尺寸的條件下(工作輥/支撐輥直徑及長度，輥頸的位置及尺寸，彎輥力的大小等)，對相應工況的輥系變形進行計算分析。由于工作輥和與支撐輥之間、工作輥與軋件間都存在接觸區(qū)，所以必須在建模的時候專門界定接觸區(qū)并指定合適的接觸單元。作為畢業(yè)設計主要的工作是應用這一方法求解實際輥系變形問題，求出不同工況下的輥縫曲線。</p><p>　　本課題采用的有限元分析方法是近年來被廣泛使用的一種求解四輥軋機輥系變形的

77、方法，其主要特點在前一節(jié)中已有交待。較之傳統(tǒng)的材料力學及分割梁的方法優(yōu)勢顯而易見，適用于大中型鋼鐵企業(yè)的熱軋生產線輥系變形的分析及控制。本課題所分析所針對的直接對象就是梅山1420熱連軋機組，具有很強的實用性。</p><p><b>　　模型設計與計算</b></p><p>　　2.1輥系模型的建立</p><p>　　2.1.1 設計原理

78、</p><p>　　圖2-1為本課題要分析的軋輥輥系的尺寸圖。 </p><p>　　圖2-1梅山1420熱連軋機輥系尺寸圖</p><p>　　在利用ANSYS軟件對此四輥軋機模型進行建模分析之前我們發(fā)現(xiàn)，由于結構和受力上的對稱的特點，故在確定模型的幾何尺寸的時候，只需對整個輥系的一半進行建模分析，就可以完成對整個模型的研究。這樣既不會影響分析結果，又

79、可以因模型的簡化而縮短一半的運行時間。對我們的設計是很有利的。</p><p>　　2.1.2 模型設計</p><p>　　根據上圖尺寸及上述設計思想，定義軋機輥系一半模型的各個關鍵點，并經由各點生成所需的實體模型。此步驟的建模思路是：生成關鍵點、由點生成面、由面旋轉成體。主要步驟的ANSYS建模程序如下：</p><p><b>　　K,1,,,, &

80、lt;/b></p><p>　　K,2,1510,,, </p><p>　　K,3,1510,-494.5955,, </p><p>　　K,4,1510,-630,, </p><p>　　K,5,800,-630,, </p><p>　　K,6,800,-494.5955,, </p>

81、<p>　　K,7,650,-440,, </p><p>　　K,8,0,-385,,</p><p>　　K,9,1510,-680,, </p><p>　　K,10,1510,-725,, </p><p>　　K,11,1510,-950,, </p><p>　　K,12,343,-950,

82、, </p><p>　　K,13,343,-725,,</p><p>　　K,14,615,-725,, </p><p>　　K,15,615,-680,, </p><p>　　K,16,800,-680,,</p><p>　　K,1001,1510,-630,, !!!與點4的坐標值相同，但代表的是不同的體

83、上的點</p><p>　　K,1002,800,-630,, !!!與點5的坐標值相同，但代表的是不同的體上的點</p><p>　　a,1,2,3,6,7,8,</p><p>　　a,3,4,5,6,</p><p>　　a,1001,9,16,1002</p><p>　　a,9,10,14,15</p

84、><p>　　a,10,11,12,13</p><p><b>　　alls</b></p><p>　　vrotat,1,,,,,,1,2,360</p><p>　　vrotat,2,,,,,,1,2,360</p><p>　　vrotat,3,,,,,,11,12,360</p>

85、<p>　　vrotat,4,,,,,,11,12,360</p><p>　　vrotat,5,,,,,,11,12,360</p><p>　　在進行此部分的建模示要特別注意k1001、k1002兩點的定義，單從坐標值看來，似乎這兩個點的定義是與點4和點5重復的，實則不然，此兩點連線所生成的直線表示的是兩輥之間初始接觸的那條線。在ANSYS中，若將同一條線賦予兩個實體，那

86、么兩個實體在此條線上即具有了相同的材料特性。但是本課題中兩輥材料不同，而且兩個體共用一條線的情況不能保證接觸對的正確創(chuàng)建。所以此步驟中在相同的空間位置定義出兩條分別隸屬于兩個不同的體是在進行體模型創(chuàng)建時的關鍵環(huán)節(jié)。</p><p>　　模型的生成方式直接決定了單元類型的選取和求解精度，在設計該模型時我做了大量的實踐，先后使用了自由(free)、掃掠(sweep)兩種不同的體單元劃分方式劃分。采用何種體單元劃分方式

87、是一定要事先明確的，因為不同的劃分方式所要求的前期模型創(chuàng)建方式也是不同的。上述程序的特點是創(chuàng)建出的實體模型便于通過掃掠的方式劃分單元網格。相比較而言掃掠獲得的單元形狀較為規(guī)范美觀，在相當精度的情況下運算起來較為省時，而且便于后處理操作。故在以后的分析中，均采用掃掠出的模型進行分析，采用自由網格劃分的模型另附于說明書后。</p><p>　　2.2接觸區(qū)域的劃分</p><p>　　2.2.

88、1 設計原理</p><p>　　2.2.1.1軋件與工作輥的接觸區(qū)</p><p>　　按照適用于熱軋帶鋼的Sims公式，可知在軋件與工作輥的接觸區(qū)中存在著前滑和后滑兩個變形區(qū)。兩個區(qū)域由中性面分開，各自滿足不同的公式的作用情況。所以要對這兩個區(qū)域進行劃分。</p><p>　　本課題所列出的初始條件如下：</p><p><b>

89、;　　k=85MPa;</b></p><p><b>　　h0=65mm;</b></p><p><b>　　h1=35mm;</b></p><p><b>　　b=1000mm</b></p><p>　　于是我們由SIMS公式可知：</p>

90、<p><b>　　咬入角</b></p><p>　　=0.307 (2-1)</p><p><b>　　中性角</b></p><p>　　=0.094 (2-2)</p>

91、;<p>　　出口至中性面之間的單位壓力</p><p><b>　　(2-3)</b></p><p>　　特殊地，當θ=0時，</p><p>　　= 66.75Mpa</p><p>　　所求的便是出口處的單位壓力</p><p>　　入口至中性面之間的單位壓力</p&g

92、t;<p><b>　　(2-4)</b></p><p>　　特殊地，當θ=γ時，</p><p><b>　　當θ=α時， </b></p><p><b>　　軋制總壓力</b></p><p>　　P=pm·A= (1/2模型上的軋制力，真實的軋

93、制力為該值的2倍)</p><p>　　依據中性角和咬入角的大小我們確定工作輥變形區(qū)域線比例為0.0598(中性面分界線)和0.1954。</p><p>　　2.2.1.2支承輥與工作輥之間的接觸區(qū)</p><p>　　四輥軋機支承輥和工作輥之間承載時有很大的接觸應力，Hertz認為：兩圓柱體在接觸區(qū)內產生局部的彈性壓扁，存在呈半橢圓形分布的壓應力。半徑方向產生的

94、法向正應力在接觸面的中部最大。最大壓應力及接觸區(qū)寬度由下式計算：</p><p>　　(2-5)式中 ——加在接觸表面單位長度上的負荷；</p><p>　　、及、——相互接觸的兩個軋輥的直徑及半徑；</p><p>　　——與軋輥材料有關的系數， ，。</p><p>　　其中，、及、為兩軋輥材料的泊松比及彈性模數。&l

95、t;/p><p><b>　　(2-6)</b></p><p>　　本課題中q是未知量，是要通過ANSYS求出的，所以兩輥接觸區(qū)厚度只能先估計出大概范圍，然后從對結果的分析中返回該步驟分析以確定q及相關的其他參數。</p><p><b>　　2.2.2模型設計</b></p><p>　　由于接觸區(qū)

96、寬度b在本課題中也屬未知，故因由分析得出。體現(xiàn)在ANSYS分析中即為接觸區(qū)域面積。接觸區(qū)域面積過小，不能模擬出兩輥間的完全接觸；而接觸面積過大，又會直接影響到求解速度。選定b=1000mm(軋件寬度)；J=1200KN(彎輥力)時的工況闡述編程原理，數次實踐后選定支承輥徑向線比例為0.15，工作輥徑向線比例為0.3。確定接觸區(qū)線段過程的ANSYS建模程序如下：</p><p>　　kl,58,0.15!!!創(chuàng)建接

97、觸面的輪廓線</p><p>　　kl,78,0.85!!!比例劃分線指令，用于創(chuàng)建界定接觸區(qū)的線上的點</p><p><b>　　kl,85,0.3</b></p><p>　　kl,105,0.7</p><p>　　kgen,2,71,74,1,-710,,,,0!!!將生成的一邊的邊界點拷貝到另一邊</p

98、><p>　　lstr,71,75!!!連接邊界點，生成邊界線</p><p>　　lstr,72,76</p><p>　　lstr,73,77</p><p>　　lstr,74,78</p><p><b>　　alls</b></p><p>　　kl,101,0.19

99、54,,!!!創(chuàng)建軋制區(qū)的輪廓線</p><p>　　kl,101,0.0598,,</p><p>　　lgen,2,101,,,-500,,,,0</p><p>　　kgen,2,79,80,1,-500,,,,0</p><p>　　lstr,79,83</p><p>　　lstr,80,84</p&g

100、t;<p>　　lsbl,158,159!!!用線劃分線指令，用于生成軋制區(qū)邊緣</p><p><b>　　alls</b></p><p>　　kl,101,0.1954,,</p><p>　　lstr,79,83</p><p>　　lsbl,161,160</p><p>

101、　　kl,101,0.0598,,</p><p>　　lstr,80,84</p><p><b>　　ldele,162</b></p><p>　　該過程對應的模型圖如圖2所示。</p><p>　　圖2-2 模型的線輪廓圖</p><p>　　上一步驟中的難點在于軋制區(qū)的劃分，因為軋制區(qū)是

102、曲面，所以線的創(chuàng)建有一定的難度。完成對線的創(chuàng)建以后，采用程序提供的面劃分指令對接觸面進行創(chuàng)建，此過程的ANSYS建模程序如下：</p><p>　　lsel,s,,,158,159,1</p><p>　　lsel,a,,,163</p><p>　　asbl,61,all</p><p><b>　　alls</b>&

103、lt;/p><p>　　asbl,101,160</p><p><b>　　alls</b></p><p>　　asbl,45,155!!!用現(xiàn)有直線劃分出加載面</p><p>　　asbl,66,157</p><p>　　asbl,30,154</p><p>　　a

104、sbl,51,156</p><p><b>　　alls</b></p><p>　　vsel,s,,,1,8!!!將支承輥包含的8個體塊粘在一起</p><p><b>　　vglue,all</b></p><p><b>　　alls</b></p>&l

105、t;p>　　vsel,s,,,9,20!!!將工作輥包含的12個體塊粘在一起</p><p><b>　　vglue,all</b></p><p><b>　　alls</b></p><p>　　在上述程序中將工作輥和支承輥的體塊分別粘貼(AGLUE)起來是必須的，一般來說我們常采用的合并體塊的方法是采用布爾加(

106、ADD)運算，將所有的體塊合并成一個整體，但本模型正是要通過對不同的體賦予不同的網格精度來構建模型，所以在此處采用布爾加命令不符合設計思路要求。而采用粘貼命令后的體塊個數沒有減少，但由于粘貼在一起的各個體之間是通過公共面聯(lián)系在一起的，內部的節(jié)點同樣被搭接在了一起，一樣可以很好的建立起平衡方程。</p><p>　　2.3分析模型的建立</p><p>　　2.3.1 單元的選取</p

107、><p>　　有限元的思想是將模型分解成有限個可以體現(xiàn)其特性的微小單元，因此單元的選取尤為重要。在對軋輥模型進行體單元選取的時候，我先后比較試驗了SOLID45、SOLID92、SOLID185、SOLID186等單元，多次試用后發(fā)現(xiàn)SOLID185單元比較滿足本課題要求。它是3維8節(jié)點結構單元(3-D 8-Node Structural Solid)，既可以很好地體現(xiàn)單元在外載作用下的受力情況，又方便程序運算(因其

108、節(jié)點數較之20節(jié)點的SOLID186等單元較少，節(jié)約了大量的運算時間)。</p><p>　　在軋制區(qū)域我選用了SHELL63單元對加載面進行劃分，SHELL63單元是具有8節(jié)點的彈性殼結構單元，可以很好地和SOLID185單元形成匹配。注意在選擇殼單元進行網格劃分的時候，要定義適當的單元厚度，并且使加載區(qū)單元劃分地盡可能細致一些，因為在進行加載時程序要從模型中獲得該區(qū)域單元中心的坐標值作為參數參加計算分析。單元

109、選取越細致，載荷施加的精度就越高。此外，SHELL63單元還肩負著另外一個任務，就是在采用掃掠創(chuàng)建體單元時，需要采用它對源面進行劃分。</p><p>　　在進行掃掠的時候，首先確定兩軋輥的中間截面為源面，在這些源面上用不同的單元尺寸劃分成密度不同的四邊形面單元，然后確定是否存在目標面以使掃掠可以進行，在對體掃掠的時候可以對預期大變形的體塊進行網格細化。在掃掠工作完成之后，要記住將多余的源面上的殼單元刪除。<

110、;/p><p>　　在工作輥和支承輥的接觸區(qū)我們選用了一個被ANSYS稱為接觸對的組操作來劃分兩輥之間的接觸區(qū)。ANSYS的接觸對創(chuàng)建會將互相接觸的實體中剛度較大的那個定義為目標面，而將剛度較小的定義為接觸面。在本課題中因支承輥的剛度較大，故將支承輥的接觸面定義為TARGE170單元(3-D Target Segment)，將工作輥的接觸面定義為CONTA174單元(3-D 8-Node Surface-to-Sur

111、face Contact)，然后對這些區(qū)域進行網格劃分。網格劃分的精度由劃分體單元時的對面上得單元尺寸控制來保證，由于接觸單元的網格精度控制的質量直接決定了求解的精確性，故可對接觸面進行線或面的局部網格細化控制。</p><p>　　2.3.2 材料特性的定義</p><p>　　在劃分單元的同時，還要對這些單元的材料特性作出一定的定義。</p><p>　　本課題

112、中工作輥和支承輥的材料選取有所不同，支承輥采用鑄鋼，工作輥采用合金鑄鐵作為材料。由《機械設計課程設計手冊》可以查出兩種材料的不同的材料性質：</p><p>　　鑄鋼: 彈性模量：175Gpa,</p><p>　　切變模量：71-84Gpa,選取80Gpa</p><p>　　泊松比： 0.25-0.29，選取0.27</p><p

113、>　　合金鑄鐵： 彈性模量：150Gpa,</p><p>　　切變模量：61Gpa,</p><p>　　泊松比： 0.25-0.29，選取0.27</p><p>　　將兩種不同的材料分別定義于兩軋輥，然后再定義接觸對。所有單元創(chuàng)建完成以后的模型如圖2-3所示。</p><p>　　圖2-3單元劃分后的模型圖</p>

114、<p>　　2.3.3 模型設計</p><p>　　用于生成軋制區(qū)殼單元的程序語句如下：</p><p><b>　　type,4</b></p><p><b>　　real,3</b></p><p>　　lsel,s,,,92,101,9!!!劃分加載區(qū)</p>&

115、lt;p>　　lsel,a,,,159,163,4</p><p>　　lesize,all,3!!!控制軸向單元尺寸</p><p><b>　　alls</b></p><p>　　lsel,s,,,158,162,2</p><p>　　lesize,all,25!!! 控制徑向單元尺寸</p>

116、<p><b>　　alls</b></p><p>　　amesh,61,103,42</p><p><b>　　alls</b></p><p>　　2.用于生成源面殼單元的程序語句如下：</p><p>　　type,4!!!殼單元的單元和實常數類型</p><

117、;p><b>　　real,3</b></p><p>　　asel,s,,,6,24,6!!!支承輥的源面殼單元劃分</p><p>　　aesize,all,100</p><p><b>　　alls</b></p><p>　　mshkey,1!!!單元形狀</p>&l

118、t;p>　　asel,s,,,6,24,6</p><p><b>　　amesh,all</b></p><p>　　lsel,s,,,105,161,56</p><p>　　lesize,all,,,20!!!設置四邊形網格尺寸</p><p><b>　　alls</b></p&

119、gt;<p>　　mat,2!!!支承輥的材料特性</p><p><b>　　mshkey,2</b></p><p>　　asel,s,,,112,116,2</p><p>　　asel,a,,,109</p><p>　　aesize,all,60</p><p><b

120、>　　amesh,all</b></p><p><b>　　alls</b></p><p>　　vatt,2,1,1!!!掃掠生成支承輥</p><p>　　vsel,s,,,1,4,1</p><p>　　vsel,a,,,21,24,1</p><p>　　esize,8

121、0,0</p><p>　　vsweep,all</p><p><b>　　alls</b></p><p>　　asel,s,,,34,42,4!!!工作輥的源面殼單元劃分</p><p>　　asel,a,,,29 </p><p>　　asel,a,,,44,11

122、9,75</p><p>　　asel,a,,,122,125,3</p><p>　　aesize,all,100</p><p>　　mshkey,1!!!工作輥的材料特性</p><p><b>　　mat,1</b></p><p><b>　　amesh,all</b&g

123、t;</p><p><b>　　alls</b></p><p>　　lsel,s,,,58,78,20</p><p>　　lesize,all,,,20</p><p><b>　　alls</b></p><p><b>　　mshkey,2</b&g