熱連軋輥系變形三維建模及有限元分析畢業(yè)設(shè)計(jì)

1、<p><b>　　1、 緒論</b></p><p>　　1.1 四輥軋機(jī)發(fā)展情況概論</p><p>　　近年來我國軋鋼行業(yè)得到了飛速發(fā)展，鋼材年產(chǎn)量突破了2億噸，已連續(xù)多年成為世界鋼產(chǎn)量第一大國。板帶材的軋制生產(chǎn)能力逐步提升到了一個(gè)較高的水平，各種板帶產(chǎn)品也得以廣泛的應(yīng)用于生產(chǎn)和生活中的方方面面。但是我國目前軋鋼生產(chǎn)的技術(shù)水平與國際先進(jìn)水平相比還有

2、相當(dāng)大的差距，軋制產(chǎn)品的主要技術(shù)指標(biāo)與國際先進(jìn)水平相比仍有相當(dāng)大的差距，我國已經(jīng)入世，國外鋼材生產(chǎn)技術(shù)強(qiáng)國的行業(yè)沖擊愈發(fā)明顯起來，要想在空前激烈的競(jìng)爭(zhēng)中得以生存、獲得發(fā)展，我們就必須在軋機(jī)精度控制等方面多做工作。</p><p>　　四輥軋機(jī)以其較高的生產(chǎn)能力和良好的產(chǎn)品質(zhì)量廣泛應(yīng)用于板帶生產(chǎn)中，近年來隨著國民經(jīng)濟(jì)的不斷發(fā)展以及工業(yè)生產(chǎn)需求的不斷增長，用戶對(duì)板帶產(chǎn)品的平直度等指標(biāo)要求越來越高，這就對(duì)板帶軋制中輥縫

3、的控制精度提出了更高的要求。對(duì)四輥軋機(jī)輥系變形進(jìn)行分析，是關(guān)乎板帶材質(zhì)量的決定性因素。如何提高軋機(jī)輥系變形分析的水平，對(duì)各個(gè)工廠來說是要亟待解決的，傳統(tǒng)的分析方法，繁雜且精度不高。本課題采用基于ANSYS軟件的有限元分析法對(duì)四輥軋機(jī)輥系變形進(jìn)行研究，是近年來一種正在被逐步廣泛應(yīng)用的方法。</p><p>　　1.2 輥系變形計(jì)算的常用理論與計(jì)算方法</p><p>　　1.2.1軋輥?zhàn)冃文?/p>

4、型的分類</p><p>　　關(guān)于板形的軋輥?zhàn)冃文Ｐ偷难芯堪l(fā)展可追溯到1958年，那時(shí)薩克斯?fàn)?Saxl)第一次對(duì)四輥軋機(jī)做了全面深入的研究。此后由于引進(jìn)了數(shù)學(xué)模型，這一領(lǐng)域得到了更進(jìn)一步的拓展。這些模型的分類如下：</p><p>　　(1) 二輥軋機(jī)的簡(jiǎn)支梁模型；</p><p>　　(2) 四輥軋機(jī)的簡(jiǎn)支梁模型；</p><p>　

5、　(3) 分割梁模型；</p><p>　　(4) 有限元分析模型。</p><p>　　1.2.2 二輥軋機(jī)的簡(jiǎn)支梁模型</p><p>　　在二輥軋機(jī)簡(jiǎn)支梁模型中，將工作輥視為線彈性應(yīng)力梁。在推導(dǎo)梁的撓曲公式時(shí)，我們做了以下假定：</p><p>　　梁的材質(zhì)均勻，在拉伸與壓縮時(shí)的彈性模量相同；</p><p>

6、;<b>　　梁的橫斷面相同;</b></p><p>　　梁至少關(guān)于一個(gè)軸向平面對(duì)稱；</p><p>　　所有的加載和反作用力都與梁的軸線垂直；</p><p>　　對(duì)于具有緊湊斷面的金屬梁，其寬高比等于或大于8。 </p><p>　　板帶材的板形可以通過對(duì)以下的軋輥的兩類撓曲進(jìn)行疊加來確定：</p>

7、;<p>　　由于軋制力引起的彎曲力使軋輥產(chǎn)生的撓曲；</p><p>　　由于軋制力引起的剪切力使軋輥產(chǎn)生的撓曲。</p><p>　　由彎曲力產(chǎn)生的撓曲可由如下的微分方程描述：</p><p><b>　　(1-1)</b></p><p>　　式中 ——軋輥彈性模量；</p>&l

8、t;p>　　——在距離x處軋輥斷面的慣性矩(圖1-1)</p><p><b>　　P—— 軋制力；</b></p><p>　　——在距離為x處的軋輥撓度；</p><p>　　L—— 軋制力作用點(diǎn)的間距；</p><p>　　—— 帶材的寬度。 </p><p>　　在軋輥

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

10、 </p><p><b>　　圖1-1</b></p><p><b>　　其中：</b></p><p>　　因剪切力產(chǎn)生的軋輥撓度由拉克計(jì)算得出，這一結(jié)果是在假定軋輥垂直斷面上的剪切應(yīng)力呈均勻分布條件下獲得的。軋輥和帶材的接觸區(qū)內(nèi)的軋輥撓度可由如下微分方程給出：</p><p&

11、gt;<b>　　(1-4)</b></p><p>　　式中： ——軋輥彈性剪切模量；</p><p>　　——在x處的軋輥撓度。</p><p>　　拉克給出了方程1-4的解如下：</p><p><b>　　(1-5)</b></p><p>　　由彎曲力和剪切力產(chǎn)生

12、的軋輥總的撓度為：</p><p><b>　　(1-6)</b></p><p>　　1.2.3 四輥軋機(jī)的簡(jiǎn)支梁模型</p><p>　　在四輥軋機(jī)的簡(jiǎn)直梁模型中，工作輥和支撐輥都被認(rèn)為是完全彈性應(yīng)力梁。</p><p>　　斯通(Stone)和戈雷(Gray)采用的模型中，四輥軋機(jī)的軋輥撓曲可以看成一個(gè)置于彈性基

13、礎(chǔ)上的簡(jiǎn)支梁的撓曲情況，如圖1-2所示。鐵木辛克(Timoshenko)推導(dǎo)了作為這一模型依據(jù)的微分方程：</p><p><b>　　(1-7)</b></p><p>　　式中 ——工作輥的彈性模量；</p><p>　　——工作輥的慣性矩；</p><p>　　—— 工作輥在x處的撓度。</p>

14、<p>　　此撓度曲線方程的通解為：</p><p><b>　　(1-8)</b></p><p>　　式中 A、B、C、D——積分常量，取決于載荷類型和邊界條件。</p><p><b>　　圖1-2</b></p><p>　　彈性量的邊界條件和載荷情況如圖1-2所示，其中

15、k1為代表支撐輥和相鄰工作輥之間彈性條件的簡(jiǎn)支梁常數(shù)，k2是代表相互接觸的工作輥和軋材彈性的常數(shù)。根據(jù)斯通模型：</p><p>　　k= k1+k2, 其中k1和k2的值為：</p><p>　　i=1,2 (1-9)</p><p>　　式中 ——兩接觸軋輥沿水平軸方向的接觸長度；</p><

16、p>　　——兩接觸軋輥中心線徑向接近量。</p><p>　　根據(jù)弗普爾(Foppl)的研究，當(dāng)支撐輥的彈性模量與工作輥的相等，即= =時(shí)，一對(duì)軋輥的中心線徑向接近量可由如下的方程確定：</p><p><b>　　(1-10)</b></p><p>　　式中 、——在壓力狀態(tài)下的工作輥和支撐輥直徑；</p><

17、p>　　b——工作輥和支撐輥的壓扁接觸寬度，等于：</p><p><b>　　(1-11)</b></p><p>　　式中 ——軋輥材質(zhì)的泊松比。</p><p><b>　　簡(jiǎn)支梁模型的局限性</b></p><p>　　簡(jiǎn)支梁撓曲模型考慮了許多影響板形的重要因素，但這些模型仍然存在

18、著以下的不足：</p><p>　　通常軋機(jī)支撐輥的徑長比遠(yuǎn)遠(yuǎn)小于8，而8時(shí)從上述方程獲得精確結(jié)果所要求的最小值。</p><p>　　針對(duì)二輥軋機(jī)推導(dǎo)出的方程1-1和1-2的解若要用于四輥軋機(jī)，應(yīng)假定在四輥軋機(jī)中，工作輥與支撐輥間的壓力是沿帶材寬度傳遞的。而事實(shí)上，此壓力的傳遞是通過整個(gè)接觸區(qū)進(jìn)行的。</p><p>　　簡(jiǎn)支梁模型不能模擬在工作輥和軋件及工作輥與

19、支撐輥間載荷橫向不均勻分布，因此，該模型沒有考慮一些重要的因素，如軋輥凸度、隨后的板凸度、材質(zhì)硬度沿帶寬方向的分布、軋輥的磨損等。</p><p>　　上述局限性使得我們有理由去研究發(fā)展更為復(fù)雜的模型來計(jì)算板形。潑普勞斯基(Poplawski)和麥克迪爾摩特(McDeermott)提出了一種可以模擬軋輥和板帶凸度，以及沿軋輥長度方向的載和 不均勻分布影響的方法。在該模型中，工作輥和

20、 圖1-3 </p><p>　　支撐輥的接觸面積工作輥和帶材的接觸面杯模擬</p><p>　　成一系列的彈簧，但軋輥仍被視為簡(jiǎn)支的彈性梁。</p><p>　　在由王國棟等提出的模型中，將簡(jiǎn)支梁布置換成了由兩個(gè)懸臂梁代表輥身

21、 的布置形式，如圖2-3所示，對(duì)于采用了橫移軋輥的四輥和六輥軋機(jī)，計(jì)算軋輥的彈性變形和出口板凸度時(shí)，采用了矩陣的方法。 </p><p>　　1.2.4分割法模型 </p><p>　　在分割法中梁的撓曲模型最先由紹特(Shohet)和湯森德(TowpseHd)提出，在此模型中，板凸度通過確定以

22、下三個(gè)未知量來求出：</p><p>　　(1)工作輥和板帶問的橫向載荷分布；</p><p>　　(2)工作輥和支承輥間的橫向載荷分布；</p><p>　　(3)工作輥的剛體移動(dòng)。</p><p>　　因?yàn)檐垯C(jī)是關(guān)于軋輥中心對(duì)稱的，所以計(jì)算過程僅需考慮軋輥的一半。所用數(shù)值方法是將軋輥切分成m個(gè)單元，如圖1-4所示，并將軋輥的分布載荷代之以

23、施加在每個(gè)單元中心的集中載荷。因?yàn)檐埣挾刃∮谳伾黹L度，所以與軋件接觸的軋輥單元數(shù)n小于m。</p><p>　　支撐輥的變形。在分割法模型中，支撐輥的變形可以表示為在位置i處的軋輥表面的垂直位移：</p><p><b>　　(1-12)</b></p><p>　　式中 ——工作輥和支撐輥間第j個(gè)單元每單位寬度上的載荷；</p&g

24、t;<p><b>　　——影響系數(shù)；</b></p><p><b>　　——單元的寬度；</b></p><p>　　——在i處的支撐輥表面的局部接觸變形；</p><p>　　——支撐輥的剛體移動(dòng)。</p><p>　　工作輥的變形工作輥?zhàn)冃慰杀硎緸楣ぷ鬏伇砻嬖谖蝗玦處的垂直位移

25、：</p><p><b>　　(1-13)</b></p><p>　　式中 ——工作輥和軋件間第i個(gè)單元上每單位寬度的載荷；</p><p>　　——工作輥的影響系數(shù)；</p><p>　　——I點(diǎn)處的系數(shù)因子(在有軋件處,=1,在無軋件處,=0)；</p><p>　　——工作輥表面在i

26、處的接觸變形；</p><p><b>　　——工作輥的剛度。</b></p><p>　　影響系數(shù)的計(jì)算。影響系數(shù)是由第j個(gè)單元重點(diǎn)的載荷在第i個(gè)單元中點(diǎn)產(chǎn)生的撓曲(圖1-4)。 </p><p>　　圖1-4

27、 </p><p>　　當(dāng)時(shí)： </p><p><b>　　(1-14)</b></p><p><

28、b>　　當(dāng)時(shí)：</b></p><p><b>　　(1-15)</b></p><p>　　將工作輥和支撐輥的參數(shù)對(duì)應(yīng)值代入方程1-14和1-15，可得系數(shù)和的值。</p><p>　　工作輥和支撐輥接觸匹配關(guān)系。工作輥和支撐輥接觸匹配可在假定無載荷時(shí)得出，此時(shí)凸工作輥和支撐輥只是點(diǎn)接觸，在此點(diǎn)以外存在輥縫γ(i)。在有軋值

29、力作用時(shí)，此輥縫會(huì)減小。那么，相互接觸的工作輥和支撐輥的撓度換算方程為：</p><p><b>　　(1-16)</b></p><p>　　工作輥和軋件接觸匹配關(guān)系。工作輥和軋件接觸匹配考慮了受負(fù)載軋輥任意點(diǎn)的輥縫高度，還有軋件在相應(yīng)點(diǎn)的出口厚度隨著軋輥的彈性壓扁和工作輥撓區(qū)的總和而變化。由此，其協(xié)調(diào)方程為：</p><p><b&g

30、t;　　(1-17)</b></p><p>　　式中 ——在i點(diǎn)處軋件出口厚度的一半；</p><p>　　——在i點(diǎn)處空載輥縫的高度。</p><p>　　軋簡(jiǎn)的出口厚度可以采用軋值理論的線性方程計(jì)算出來，其中是軋件入口厚度、軋件力、帶鋼張力和變形抗力的函數(shù)。</p><p>　　(6)工作輥的靜平衡。工作輥的靜平衡時(shí)將

31、軋輥間的垂直載荷、軋件和工作輥間的載荷及彎輥設(shè)備施加給工作輥的載荷進(jìn)行疊加，則方程為：</p><p>　　(1-18) </p><p>　　式中 ——總的軋制彎曲力(見圖1-4)。</p><p>　　(7)方程的解。求解方程1-16、1-17、1-18在于如何找出m+n+1個(gè)未

32、知量，即表示工作輥和支撐輥之間載荷分布的m個(gè)力的值(),表示工作輥和軋件間載荷分布的n個(gè)力的值()和工作輥的剛度系數(shù)。有了上述三個(gè)基本方程，可以得到求解、和的m+n+1個(gè)方程。這些方程可以用代數(shù)矩陣的方法求解。</p><p>　　分割梁撓曲模型還考慮了帶材沿寬度方向的張力，這里另外增加了n個(gè)未知量。這一問題可以通過疊加的方式來解決。 圖1-5 </p>&

33、lt;p>　　在由郭任明(Guo)研制的模型中，將軋輥劃分為一系列彈簧單元(圖1-5)。這些彈簧變形的相互關(guān)系通過線性相關(guān)方法計(jì)算出來。這種方法考慮了由彎矩和剪切力產(chǎn)生的梁的撓曲變形。這一模型也將帶材視為一系列的彈簧單元，并認(rèn)為這些彈簧單元的剛度是帶材軋制特性的函數(shù)?；籼m德(Hollander)和萊茵(Reinen)曾經(jīng)對(duì)另一種模型做了闡述。</p><p>　　分割梁模型的局限性

34、 </p><p>　　分割梁撓曲模型的研究發(fā)展對(duì)于在軋制過程中提高帶材板形的模擬能力邁出了很大一步。然而，這類模型也有自身的不足，因?yàn)榉指盍耗Ｐ褪墙⒃诩僭O(shè)在軋制力作用下工作輥和支承輥完全接觸之上的。實(shí)際上，當(dāng)采用特殊的輥型系統(tǒng)如CVC軋輥、UPC軋輥和錐形軋輥時(shí)，情況并非如此。此時(shí)還需考慮可能存在的接觸面不吻合(如圖1-6)。</p><p>　　此模型的另一不足是該模型對(duì)

35、影響系數(shù)的計(jì)算是根據(jù)簡(jiǎn)支架撓曲方程而來的。但是，如前所述當(dāng)軋輥的徑長比小時(shí)，這些方程的實(shí)用性是值得懷疑的，而且分割梁模型用二維問題代替三維 圖1-6</p><p>　　問題，因此在某些情況下，不可能獲得良好的計(jì)算精度。 </p><p>　　1.2.5 有限元分析理論

36、 </p><p>　　在有限元分析中采用的是矩陣結(jié)構(gòu)分析方法。這種方法采用了一種直接的物理方法來建立和求解梁及框架結(jié)構(gòu)問題。</p><p>　　求解工程問題的解要用到以下三個(gè)條件：</p><p><b>　　力平衡方程；</b></p>

37、;<p><b>　　變形協(xié)調(diào)方程；</b></p><p>　　材料行為的本構(gòu)關(guān)系。</p><p>　　利用這三個(gè)條件可以建立未知應(yīng)力(力法)或未知位移(位移法)組成的方程。在有限元技術(shù)中，常用位移法。</p><p>　　有限元分析模型是對(duì)如下一些連續(xù)物體問題進(jìn)行近似求解的一種方法：</p><p>

38、　　連續(xù)物體被劃分成有限個(gè)單元，各個(gè)單元的行為由有限個(gè)參數(shù)給定；</p><p>　　整個(gè)系統(tǒng)用它的單元集合體求解同樣精確地遵循適用于標(biāo)準(zhǔn)離散問題的那些原則。</p><p>　　由有限個(gè)單元構(gòu)成的集合體通常是指網(wǎng)格，這些單元通過節(jié)點(diǎn)相互連接起來。在平面問題中，單元可以是三角形或是四邊形的。在三維問題中單元可以是三棱柱或長方體或六棱柱。</p><p>　　有限元分

39、析采用了多年發(fā)展起來的應(yīng)用于離散問題的標(biāo)準(zhǔn)方法。這套方法包含了對(duì)物體或結(jié)構(gòu)的每一個(gè)單元的力-位移關(guān)系的計(jì)算。按既定的步驟，對(duì)每個(gè)有限元的每個(gè)節(jié)點(diǎn)都能建立起局部的力平衡方程來。求解這些方程可得到未知位移的解。有限元分析方法一般分為以下六個(gè)步驟：</p><p>　　(1)選擇位移模型。位移模型將各個(gè)單元的位移表示成各個(gè)節(jié)點(diǎn)位移的函數(shù)，通常表示為：</p><p><b>　　(1-

40、19)</b></p><p>　　式中 ——單元的位移矩陣；</p><p>　　——單元的形函數(shù)矩陣；</p><p>　　——單元的節(jié)點(diǎn)位移矩陣(未知的變量)。</p><p>　　建立位移模型時(shí)，沿給定方向的位移分布，通常用一個(gè)簡(jiǎn)單的函數(shù)表示，例如下面的多項(xiàng)式：</p><p><b&g

41、t;　　(1-20)</b></p><p>　　式中 ——沿x方向的位移。</p><p>　　多項(xiàng)式的系數(shù)被稱為廣義位移系數(shù)，他們確定了位移模型的形狀。通過邊界條件可求得上述系數(shù)。有限元方程的個(gè)數(shù)依賴于要被模型化的結(jié)構(gòu)的幾何形狀及模型可能移動(dòng)的方向的個(gè)數(shù)，即自由度。</p><p>　　簡(jiǎn)單的彈簧單元可承受張力并有兩個(gè)自由度，因?yàn)榱旱膬啥丝梢匝貜椈?/p>

42、的主軸方向自由移動(dòng)。然而，具有更多自由度的單元體可以作其他運(yùn)動(dòng)如彎曲運(yùn)動(dòng)。單元類型的選擇也決定了在一個(gè)模型中的自由度。</p><p>　　(2) 應(yīng)變-位移關(guān)系的建立。應(yīng)變-位移關(guān)系可通常表示為：</p><p><b>　　(1-21)</b></p><p>　　式中 ——單元應(yīng)變矩陣；</p><p>　　

43、——單元應(yīng)變-位移關(guān)系矩陣。</p><p>　　應(yīng)變-應(yīng)變關(guān)系的建立。應(yīng)變-應(yīng)變關(guān)系可表示為：</p><p>　　(1-22) </p><p>　　式中 ——單元應(yīng)力矩陣；</p><p>　　——單元材料的應(yīng)力-應(yīng)變關(guān)系矩陣。</p><p>　　在彈性變形時(shí)，材料的應(yīng)力—應(yīng)變關(guān)系矩陣[C

44、]常根據(jù)虎克(Hooke)定律得到。這一矩陣建立起了線應(yīng)變、剪應(yīng)變、正應(yīng)力和剪應(yīng)力之間的關(guān)系(如圖1-7所示)：</p><p><b>　　(1-23)</b></p><p>　　矩陣[C]中的系數(shù)常用楊氏模量E和泊松比表示。</p><p>　　(4) 建立有限元?jiǎng)偠染仃嚒Ｓ邢拊獎(jiǎng)偠葘⒐?jié)點(diǎn)位移與節(jié)點(diǎn)處的力聯(lián)系起來，即：</p>

45、<p><b>　　(1-24)</b></p><p>　　式中 ——節(jié)點(diǎn)處的力矢量(已知量)；</p><p>　　——單元?jiǎng)偠染仃?剛度影響系數(shù))。</p><p>　　剛度矩陣可由最小余能原理求得，通常表示為：</p><p><b>　　(1-25)</b>&l

46、t;/p><p>　　式中 ——[B]的轉(zhuǎn)置矩陣；</p><p><b>　　V——單元的體積。</b></p><p>　　(5)聯(lián)立代數(shù)方程。這一步包括對(duì)由每個(gè)單元的剛度矩陣[k]構(gòu)成的總剛度矩陣[K]和由各個(gè)節(jié)點(diǎn)的力矢量{P}構(gòu)成的總的力矢量{R}的聯(lián)立。</p><p>　　聯(lián)立方法的根據(jù)是由于每個(gè)節(jié)點(diǎn)的相互

47、聯(lián)系要求所有于該節(jié)點(diǎn)相鄰的節(jié)點(diǎn)在該節(jié)點(diǎn)位移應(yīng)相同?？偟膭偠染仃嚒⒖偟牧κ噶亢涂偟奈灰剖噶縶r}的平衡關(guān)系可有一套聯(lián)立方程表示：</p><p><b>　　(1-26)</b></p><p>　　(6)求解未知參數(shù)。考慮物體的幾何和力的邊界條件，求解代數(shù)方程1-16，可得所有未知的位移。在線性平衡問題中，可直接應(yīng)用矩陣代數(shù)技術(shù)來求解：</p><

48、p><b>　　(1-27)</b></p><p>　　1.2.6 二維有限元分析模型</p><p>　　圖1-7時(shí)一個(gè)典型的四輥軋機(jī)軋輥系統(tǒng)的二維有限元網(wǎng)格圖，位于X-Y平面內(nèi)。在最簡(jiǎn)單的情況中，Z方向的厚度是一個(gè)常量，這很大程度上減少了模型中有限元的數(shù)量，使我們無需額外增加大量的計(jì)算時(shí)間就可設(shè)計(jì)出一個(gè)十分準(zhǔn)確的有限元網(wǎng)格。但是，只用二維作分析而忽略第三維

49、的影響，我們很難得到高精度的結(jié)果來。</p><p>　　這一問題在一定程度上被陳先霖(XianLin)和鄒家祥(Jiaxiang)解 圖1-7</p><p>　　決了,他們采用了一種在Z方向上具有不同厚度的二維模型。在這一模型中，每個(gè)單 </p>

50、<p>　　元根據(jù)其到軋輥軸線的縱坐標(biāo)軸距離和軋輥的半徑對(duì)應(yīng)不同的厚度，第i個(gè)單元對(duì)應(yīng)的厚度為：</p><p><b>　　(1-28)</b></p><p>　　式中 R——軋輥半徑。</p><p>　　參數(shù)和由下面的式子定義：</p><p><b>　　(1-29)</b>

51、</p><p>　　式中 ——第i個(gè)單元的縱坐標(biāo)；</p><p>　　——第i個(gè)單元的高度。</p><p>　　盡管二維模型只是物理模型的一種簡(jiǎn)化表示，當(dāng)其在分析許多因素如軋輥壓扁和帶材張力對(duì)板形的影響時(shí)仍不失為一個(gè)極其方便的工具。</p><p>　　1.2.7 三維有限元分析模型</p><p>　　三維有限

52、元網(wǎng)格為一個(gè)所研究系統(tǒng)的物理模型提供了最確切的表示。但是，在確定網(wǎng)格中有限元的數(shù)量和類型時(shí)要費(fèi)點(diǎn)勁。</p><p>　　圖1-8描述了由一對(duì)軋輥和軋制帶材構(gòu)成的體系的三維有限元網(wǎng)格，曾應(yīng)用于由聯(lián)合工程公司和國際軋鋼咨詢公司聯(lián)合研制的ROLL-FLEX離線模型中。</p><p>　　在研究網(wǎng)格時(shí)，單元的數(shù)量或向格的尺寸可由所要求的板凸度計(jì)算精度而定。但是，選擇單元的數(shù)量還要考慮計(jì)算時(shí)間和

53、花費(fèi)的成本。在評(píng)價(jià)計(jì)算精度時(shí)，要用到如下的參數(shù)：</p><p>　　(1)中心凸度的相對(duì)誤差；</p><p>　　(2)軋輥壓扁影響的相對(duì)誤差。</p><p>　　軋輥中心凸度的相對(duì)誤差可表示為：</p><p>　　(1-30)

54、 </p><p>　　圖1-8 </p><p>　　式中

55、 、——在有限元網(wǎng)格中，當(dāng)網(wǎng)格數(shù)為n和∞時(shí)分別對(duì)應(yīng)的中心凸度；</p><p>　　——帶材的出口厚度。</p><p>　　因?yàn)橹行耐苟鹊南鄬?duì)誤差隨著網(wǎng)格數(shù)量的增加而減小，所以我們可以根據(jù)指數(shù)定律來確定的值。</p><p>　　如圖1-9所示，隨著軋材剛度的增加，為了達(dá)到相同的精度，有限元網(wǎng)格的數(shù)量也增加了。當(dāng)中心凸度的相對(duì)誤差為0.25%時(shí)，表示F1、F4

56、、F6機(jī)架的網(wǎng)格數(shù)分別減少為2500、5000和7000個(gè)。軋輥壓扁影響的相對(duì)誤差由如下方程計(jì)算得出：</p><p><b>　　(1-31)</b></p><p>　　式中 ——有限元分析中有n個(gè)網(wǎng)格時(shí)工作輥和支撐輥的中心線徑向接近量；</p><p>　　——由軋輥壓扁方程計(jì)算得到的工作輥和支撐輥的中心線徑向接近量。</p&

57、gt;<p>　　根據(jù)弗浦爾的研究，當(dāng)工作輥和支撐輥具有相同的彈性模量E時(shí)，一對(duì)軋輥如工作輥和支撐輥的中心接近量 可由方程1-10和1-11確定。若要將F1、F4、F6機(jī)架的軋輥的壓扁影響的相對(duì)誤差減小到0.25%,那么所要求的有限元網(wǎng)格數(shù)量必須分別等于或大于3500、4600和5200個(gè)，如圖1-10所示。</p><p><b>　　圖1-9</b></p>

58、<p><b>　　圖1-10 </b></p><p>　　1.3 ANSYS軟件</p><p>　　1.3.1 ANSYS的工作原理——有限單元法</p><p>　　有限單元法的基本思想是將連續(xù)的結(jié)構(gòu)離散成有限個(gè)單元，并在每一個(gè)單元中設(shè)定有限個(gè)節(jié)點(diǎn)，將連續(xù)體看作是只在節(jié)點(diǎn)處相連接的一組單元的集合體；同時(shí)選定場(chǎng)函數(shù)的節(jié)點(diǎn)值作為

59、基本未知量，并在每一單元中假設(shè)一近似差值函數(shù)以表示單元中場(chǎng)函數(shù)的分布規(guī)律；進(jìn)而利用力學(xué)中的某些變分原理去建立用于求解節(jié)點(diǎn)未知量的有限元法方程，從而將一個(gè)連續(xù)域中的無限自由度問題化為離散域中的有限自由度問題。一經(jīng)求解就可以利用解得的節(jié)點(diǎn)值和設(shè)定的插值函數(shù)確定單元上以至整個(gè)集合體上的場(chǎng)函數(shù)。有限元求解程序的內(nèi)部過程可從圖1-11中看出。</p><p><b>　　↓</b></p>

60、<p><b>　　↓</b></p><p><b>　　↓</b></p><p><b>　　↓</b></p><p><b>　　↓</b></p><p><b>　　↓</b></p><

61、;p><b>　　圖1-11</b></p><p>　　近40年來，隨著計(jì)算機(jī)的飛速發(fā)展和廣泛應(yīng)用，各種行之有效的數(shù)值計(jì)算方法得到了巨大的發(fā)展。而有限元方法則是計(jì)算機(jī)誕生以后，在計(jì)算數(shù)學(xué)、計(jì)算力學(xué)、和計(jì)算工程科學(xué)領(lǐng)域里誕生的最有效的計(jì)算方法。隨著有限元理論基礎(chǔ)的日益完善，出現(xiàn)了很多通用和專用的有限元計(jì)算軟件。在國際畫的市場(chǎng)經(jīng)濟(jì)中，企業(yè)間的競(jìng)爭(zhēng)日益加劇。為取得競(jìng)爭(zhēng)優(yōu)勢(shì)，企業(yè)迫切需要以高

62、質(zhì)量低成本的產(chǎn)品迅速搶占市場(chǎng)，因此企業(yè)迫切需要高技術(shù)、高速度、低成本的設(shè)計(jì)方法。ANSYS程序即是應(yīng)此要求而發(fā)展起來的計(jì)算機(jī)仿真設(shè)計(jì)工具。</p><p>　　1.3.2 ANSYS軟件簡(jiǎn)介</p><p>　　ANSYS軟件由成立于1970年的美國ANSYS公司開發(fā)完成，是融結(jié)構(gòu)、熱、流體、電磁、聲學(xué)于一體的大型通用有限元分析軟件，可廣泛應(yīng)用于核工業(yè)、鐵道、石油化工、航空航天、機(jī)械制造、

63、能源、汽車交通、國防軍工、電子、土木工程、造船、生物醫(yī)學(xué)、輕工、地礦、水利、日用家電等一般工業(yè)及科學(xué)研究。</p><p>　　ANSYS程序是一個(gè)功能強(qiáng)大的、靈活的設(shè)計(jì)分析及優(yōu)化軟件包。該軟件可浮動(dòng)運(yùn)行于從PC機(jī)、NT工作站、UNIX工作站甚至巨型機(jī)的各類計(jì)算機(jī)及操作系統(tǒng)中，數(shù)據(jù)文件在其所有的產(chǎn)品系列和工作品臺(tái)上均兼容。其多物理場(chǎng)耦合的功能，允許在同一模型上進(jìn)行各式各樣的耦合計(jì)算，如：熱—結(jié)構(gòu)耦合、磁—結(jié)構(gòu)耦合

64、、以及電—磁—流體—熱耦合。在PC機(jī)上生成的模型同樣可以運(yùn)行于巨型機(jī)上，這樣就保證了所有的ANSYS用戶的多領(lǐng)域多變工程的求解。</p><p>　　1.3.2.1 ANSYS軟件的組成</p><p>　　ANSYS軟件 主要包括三個(gè)部分：前處理模塊、分析計(jì)算模塊和后處理模塊。</p><p><b>　　前處理模塊</b></p>

65、;<p>　　它為用戶提供了一個(gè)強(qiáng)大的實(shí)體建模及網(wǎng)絡(luò)劃分工具，用戶可以方便地構(gòu)造有限元模型，軟件提供了100種以上的單元類型，用來模擬工程中的各種結(jié)構(gòu)和材料。</p><p><b>　　實(shí)體建模</b></p><p><b>　　參數(shù)化建模</b></p><p><b>　　體素庫及布爾運(yùn)算&

66、lt;/b></p><p>　　拖拉、旋轉(zhuǎn)、拷貝、蒙皮、倒角等。</p><p>　　多種網(wǎng)格自動(dòng)劃分工具，自動(dòng)進(jìn)行單元形態(tài)、求解精度檢查及修正。</p><p>　　自動(dòng)/映射網(wǎng)格劃分、智能網(wǎng)格劃分、自適應(yīng)網(wǎng)格劃分。</p><p>　　復(fù)雜幾何體Sweep映射網(wǎng)格生成。</p><p>　　六面體向四面體自

67、動(dòng)過渡網(wǎng)格：金字塔形。</p><p><b>　　邊界層網(wǎng)格劃分</b></p><p>　　在幾何模型或FE模型上加載：點(diǎn)載荷、分布載荷、體載荷、函數(shù)載荷。</p><p>　　可擴(kuò)展的標(biāo)準(zhǔn)梁截面形狀庫</p><p><b>　　分析計(jì)算模塊</b></p><p>　

68、　包括結(jié)構(gòu)分析（可進(jìn)行線性分析、非線性分析和高度非線性分析）、流體動(dòng)力學(xué)分析、電磁場(chǎng)分析、聲場(chǎng)分析、壓電分析及多物理場(chǎng)的耦合分析，可模擬多種物理介質(zhì)的相互作用，具有靈敏度分析及優(yōu)化分析能力。</p><p><b>　　后處理模塊</b></p><p>　　可將計(jì)算結(jié)果以彩色等值線顯示、梯度顯示、矢量顯示、力子流跡顯示、立體切片顯示、透明及半透明顯示（可看到內(nèi)部結(jié)構(gòu)

69、）等圖形方式顯示出來，也可以將計(jì)算結(jié)果以圖表、曲線形式顯示或輸出。具體如下所示：</p><p>　　計(jì)算報(bào)告自動(dòng)生成及定制工具，自動(dòng)生成符合要求格式的計(jì)算報(bào)告。</p><p>　　結(jié)果顯示菜單：圖形顯示、抓圖、結(jié)果列表。</p><p>　　圖形：云圖、等值線、矢量顯示、粒子流跡顯示、切片、透明及半透明顯示、紋理。</p><p>　　鋼

70、筋混凝土單元可顯示單元內(nèi)的鋼筋、開裂情況以及壓碎部位。</p><p>　　梁、管、板、復(fù)合材料單元及結(jié)果按實(shí)際形狀顯示，顯示橫截面結(jié)果；顯示梁?jiǎn)卧獜澗貓D。</p><p>　　顯示優(yōu)化靈敏度及優(yōu)化變量曲線。</p><p>　　各種結(jié)果動(dòng)畫顯示，可獨(dú)立保存及重放。</p><p><b>　　3D圖形注釋功能。</b>

71、</p><p>　　直接生成BMP、JPG、VRML、WMF、PNG、PS、TIFF、HPGL等格式的圖形。</p><p>　　計(jì)算結(jié)果排序、檢索、列表及再組合。</p><p>　　提供對(duì)計(jì)算結(jié)果的加、減、積分、微分等計(jì)算。</p><p>　　顯示沿任意路徑的結(jié)果曲線，并可進(jìn)行路徑的數(shù)學(xué)計(jì)算。</p><p>

72、　　1.3.2.2 ANSYS主要的技術(shù)特點(diǎn)</p><p>　　● 唯一能實(shí)現(xiàn)多場(chǎng)及多耦合分析的軟件</p><p>　　唯一實(shí)現(xiàn)前后處理、求解及多場(chǎng)分析統(tǒng)一數(shù)據(jù)庫的一體化大型FEA分析軟件</p><p>　　唯一具有物理場(chǎng)優(yōu)化功能的FEA軟件</p><p>　　唯一具有中文界面的大型通用有限元分析軟件</p><p

73、>　　具有強(qiáng)大的非線性分析功能</p><p>　　具有適用于不同的問題和硬件配置的多種求解器</p><p>　　支持異種異構(gòu)功能網(wǎng)絡(luò)浮動(dòng)，在異種、異構(gòu)平臺(tái)上支持界面統(tǒng)一，數(shù)據(jù)文件通用</p><p>　　強(qiáng)大的并行計(jì)算功能，支持分布式并行和共享內(nèi)存式并行</p><p>　　多種用戶網(wǎng)格劃分技術(shù)</p><p&g

74、t;<b>　　完善的用戶開發(fā)環(huán)境</b></p><p>　　同時(shí)，ANSYS軟件擁有豐富和完善的單元庫、材料模型庫和求解器，保證了他能高效的求解各類結(jié)構(gòu)的靜力、動(dòng)力、振動(dòng)、線性和非線性問題，壓縮和不可壓縮的流體問題。其友好的圖形界面和程序結(jié)構(gòu)，交互式的前后處理和圖形軟件，大大的減輕了用戶在實(shí)際工程問題中創(chuàng)建模型、有限元求解以及結(jié)果分析和評(píng)價(jià)的工作量。他的統(tǒng)一集中式的數(shù)據(jù)庫保證了個(gè)模塊之間

75、的有效可靠的集成，并實(shí)現(xiàn)了與多個(gè)CAD/CAE軟件的友好鏈接。</p><p>　　1.4 本文主要研究?jī)?nèi)容及創(chuàng)新</p><p>　　在深入理解有限元分析實(shí)質(zhì)以及ANSYS應(yīng)用的基礎(chǔ)上，即可對(duì)四輥軋機(jī)工作輥的變形進(jìn)行計(jì)算。本課題要求在給定輥系的幾何尺寸的條件下(工作輥/支撐輥直徑及長度，輥頸的位置及尺寸，彎輥力的大小等)，對(duì)相應(yīng)工況的輥系變形進(jìn)行計(jì)算分析。由于工作輥和與支撐輥之間、工作輥

76、與軋件間都存在接觸區(qū)，所以必須在建模的時(shí)候?qū)ｉT界定接觸區(qū)并指定合適的接觸單元。作為畢業(yè)設(shè)計(jì)主要的工作是應(yīng)用這一方法求解實(shí)際輥系變形問題，求出不同工況下的輥縫曲線。</p><p>　　本課題采用的有限元分析方法是近年來被廣泛使用的一種求解四輥軋機(jī)輥系變形的方法，其主要特點(diǎn)在前一節(jié)中已有交待。較之傳統(tǒng)的材料力學(xué)及分割梁的方法優(yōu)勢(shì)顯而易見，適用于大中型鋼鐵企業(yè)的熱軋生產(chǎn)線輥系變形的分析及控制。本課題所分析所針對(duì)的直接

77、對(duì)象就是梅山1420熱連軋機(jī)組，具有很強(qiáng)的實(shí)用性。</p><p><b>　　模型設(shè)計(jì)與計(jì)算</b></p><p>　　2.1輥系模型的建立</p><p>　　2.1.1 設(shè)計(jì)原理</p><p>　　圖2-1為本課題要分析的軋輥輥系的尺寸圖。 </p><p>　　圖2-1梅山1

78、420熱連軋機(jī)輥系尺寸圖</p><p>　　在利用ANSYS軟件對(duì)此四輥軋機(jī)模型進(jìn)行建模分析之前我們發(fā)現(xiàn)，由于結(jié)構(gòu)和受力上的對(duì)稱的特點(diǎn)，故在確定模型的幾何尺寸的時(shí)候，只需對(duì)整個(gè)輥系的一半進(jìn)行建模分析，就可以完成對(duì)整個(gè)模型的研究。這樣既不會(huì)影響分析結(jié)果，又可以因模型的簡(jiǎn)化而縮短一半的運(yùn)行時(shí)間。對(duì)我們的設(shè)計(jì)是很有利的。</p><p>　　2.1.2 模型設(shè)計(jì)</p><

79、p>　　根據(jù)上圖尺寸及上述設(shè)計(jì)思想，定義軋機(jī)輥系一半模型的各個(gè)關(guān)鍵點(diǎn)，并經(jīng)由各點(diǎn)生成所需的實(shí)體模型。此步驟的建模思路是：生成關(guān)鍵點(diǎn)、由點(diǎn)生成面、由面旋轉(zhuǎn)成體。主要步驟的ANSYS建模程序如下：</p><p><b>　　K,1,,,, </b></p><p>　　K,2,1510,,, </p><p>　　K,3,1510,-494

80、.5955,, </p><p>　　K,4,1510,-630,, </p><p>　　K,5,800,-630,, </p><p>　　K,6,800,-494.5955,, </p><p>　　K,7,650,-440,, </p><p>　　K,8,0,-385,,</p><p

81、>　　K,9,1510,-680,, </p><p>　　K,10,1510,-725,, </p><p>　　K,11,1510,-950,, </p><p>　　K,12,343,-950,, </p><p>　　K,13,343,-725,,</p><p>　　K,14,615,-725,, &

82、lt;/p><p>　　K,15,615,-680,, </p><p>　　K,16,800,-680,,</p><p>　　K,1001,1510,-630,, !!!與點(diǎn)4的坐標(biāo)值相同，但代表的是不同的體上的點(diǎn)</p><p>　　K,1002,800,-630,, !!!與點(diǎn)5的坐標(biāo)值相同，但代表的是不同的體上的點(diǎn)</p>

83、<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><p>　　a,10,11,12,13</p><p><b>　　alls</b><

84、/p><p>　　vrotat,1,,,,,,1,2,360</p><p>　　vrotat,2,,,,,,1,2,360</p><p>　　vrotat,3,,,,,,11,12,360</p><p>　　vrotat,4,,,,,,11,12,360</p><p>　　vrotat,5,,,,,,11,12,3

85、60</p><p>　　在進(jìn)行此部分的建模示要特別注意k1001、k1002兩點(diǎn)的定義，單從坐標(biāo)值看來，似乎這兩個(gè)點(diǎn)的定義是與點(diǎn)4和點(diǎn)5重復(fù)的，實(shí)則不然，此兩點(diǎn)連線所生成的直線表示的是兩輥之間初始接觸的那條線。在ANSYS中，若將同一條線賦予兩個(gè)實(shí)體，那么兩個(gè)實(shí)體在此條線上即具有了相同的材料特性。但是本課題中兩輥材料不同，而且兩個(gè)體共用一條線的情況不能保證接觸對(duì)的正確創(chuàng)建。所以此步驟中在相同的空間位置定義出兩條

86、分別隸屬于兩個(gè)不同的體是在進(jìn)行體模型創(chuàng)建時(shí)的關(guān)鍵環(huán)節(jié)。</p><p>　　模型的生成方式直接決定了單元類型的選取和求解精度，在設(shè)計(jì)該模型時(shí)我做了大量的實(shí)踐，先后使用了自由(free)、掃掠(sweep)兩種不同的體單元?jiǎng)澐址绞絼澐帧２捎煤畏N體單元?jiǎng)澐址绞绞且欢ㄒ孪让鞔_的，因?yàn)椴煌膭澐址绞剿蟮那捌谀Ｐ蛣?chuàng)建方式也是不同的。上述程序的特點(diǎn)是創(chuàng)建出的實(shí)體模型便于通過掃掠的方式劃分單元網(wǎng)格。相比較而言掃掠獲得的單

87、元形狀較為規(guī)范美觀，在相當(dāng)精度的情況下運(yùn)算起來較為省時(shí)，而且便于后處理操作。故在以后的分析中，均采用掃掠出的模型進(jìn)行分析，采用自由網(wǎng)格劃分的模型另附于說明書后。</p><p><b>　　圖3-21</b></p><p>　　2) b=800mm(軋件寬度)，J=1200KN(彎輥力)；</p><p><b>　　圖3-22&l

88、t;/b></p><p>　　3) b=800mm(軋件寬度)，J=800KN(彎輥力)；</p><p><b>　　圖3-23</b></p><p>　　4) b=800mm(軋件寬度)，J=0KN(彎輥力)；</p><p><b>　　圖3-24</b></p>&l

89、t;p>　　5) b=1000mm(軋件寬度)，J=1200KN(彎輥力)，ε=1/3(相對(duì)壓下率)；</p><p><b>　　圖3-25</b></p><p>　　3.4 接觸面上的接觸壓力分布</p><p>　　接觸變形的查看方法和查看撓曲線的方法相同，先定義出路徑，然后將單位接觸應(yīng)力映射到路徑上，由于單位應(yīng)力與加在接觸表面單

90、位長度上的載荷之間在代數(shù)上有著0.049倍的關(guān)系，故只需將公式LabR=(FACT1*Lab1)+( FACT2*Lab2)中的FACT1令為0.049，將Lab2令為none defined即可。 </p><p>　　b=1000mm(軋件寬度)，J=1200KN(彎輥力)；</p><p>　　圖3-26接觸面接觸壓應(yīng)力分布示意圖</p><p>　　圖3-

91、27接觸面接觸壓力q分布</p><p>　　2) b=800mm(軋件寬度)，J=1200KN(彎輥力)；</p><p>　　圖3-28接觸面接觸壓應(yīng)力分布示意圖</p><p>　　圖3-29接觸面接觸壓力q分布</p><p>　　3) b=800mm(軋件寬度)，J=800KN(彎輥力)；</p><p>　

92、　圖3-30接觸面接觸壓應(yīng)力分布示意圖</p><p>　　圖3-31接觸面接觸壓力q分布</p><p>　　4)b=800mm(軋件寬度)，J=0KN(彎輥力)；</p><p>　　圖3-32接觸面接觸壓應(yīng)力分布示意圖</p><p>　　圖3-33接觸面接觸壓力q分布</p><p>　　5) b=1000mm

93、(軋件寬度)，J=1200KN(彎輥力)，ε=1/3(相對(duì)壓下率)；</p><p>　　圖3-34接觸面接觸壓應(yīng)力分布示意圖</p><p>　　圖3-35接觸面接觸壓力q分布</p><p>　　3.5 繪制軋件表面變形曲線</p><p>　　1) b=1000mm(軋件寬度)，J=1200KN(彎輥力)；</p><

94、;p><b>　　圖3-36</b></p><p>　　2) b=800mm(軋件寬度)，J=1200KN(彎輥力)；</p><p><b>　　圖3-37</b></p><p>　　3) b=800mm(軋件寬度)，J=800KN(彎輥力)；</p><p><b>　　圖3

95、-38</b></p><p>　　4) b=800mm(軋件寬度)，J=0KN(彎輥力)；</p><p><b>　　圖3-39</b></p><p>　　5) b=1000mm(軋件寬度)，J=1200KN(彎輥力)，ε=1/3(相對(duì)壓下率)；</p><p><b>　　圖3-40<

96、/b></p><p>　　3.6 圖形數(shù)據(jù)分析</p><p>　　3.6.1 對(duì)輥系變形情況的分析</p><p>　　在對(duì)工作輥輥系變形的分析中我們發(fā)現(xiàn)，工作輥在施加彎輥力的地方變形最大，這主要是因?yàn)閺澼伭κ鞘┘釉趩我还?jié)點(diǎn)上的，相對(duì)節(jié)點(diǎn)位移較大。伴隨著彎輥力的逐漸降低，節(jié)點(diǎn)位移變化量也將減小。由對(duì)第二種條件下的位移變化的圖形分析我們又可以看出，在彎輥力不

97、變的情況下，節(jié)點(diǎn)最大位移隨著軋制板寬的減小而增大。此外，由56頁的圖3-41和圖3-42 (b=800mm,J=1200KN)所示，除彎輥力施加處外，在支承輥輥頸位置安裝動(dòng)壓油膜軸承的位置，由于近似地將約束條件簡(jiǎn)化成在兩個(gè)節(jié)點(diǎn)上施加約束，故在圖形中約束點(diǎn)的應(yīng)力值很大。我們可以通過改善約束條件來消除這個(gè)不利影響，改用一個(gè)小區(qū)域的面約束就能很好地解決這個(gè)問題。</p><p>　　3.6.2. 支承輥和工作輥的撓度分

98、析</p><p>　　3.6.2.1與簡(jiǎn)支梁算法求出的撓度值的比較</p><p>　　一般常采用的求解撓曲線的方法是將承載軋輥看作簡(jiǎn)支梁，用材料力學(xué)中計(jì)算直短梁撓區(qū)的方法來處理。由于坐標(biāo)零位選取的不同，故本課題中不對(duì)撓度值進(jìn)行分析，而只分析撓度差值的大小。現(xiàn)以，為例計(jì)算分析如下：</p><p>　　輥身中點(diǎn)和輥身邊上的撓度差值</p><p

99、><b>　?。?-1）</b></p><p><b>　　對(duì)于支承輥來說</b></p><p>　　彎矩引起的 （3-2）</p><p>　　切力引起的 （3-3）<

100、/p><p><b>　　故 </b></p><p><b>　　對(duì)于工作輥來說</b></p><p>　　彎矩引起的 （3-4）</p><p>　　切力引起的 （3-5）&l

101、t;/p><p><b>　　故 </b></p><p>　　其中 E和G——彈性模數(shù)和剪切模數(shù),</p><p><b>　　對(duì)于支承輥來說，</b></p><p><b>　　對(duì)于工作輥來說，；</b></p><p><b>

102、　　P——軋制力,；</b></p><p><b>　　b——軋件寬度,；</b></p><p>　　a——軋輥軸承中心線之間的距離,；</p><p>　　K——截面系數(shù)，對(duì)圓截面K=10/9。</p><p>　　由采用ANSYS計(jì)算分析后的相同工況下的撓度曲線圖中可以看出,支承輥從中間截面到輥頸面處

103、的撓度差為；工作輥從中間截面到輥頸面處的撓度差為。兩種不同的求解方法得出的撓度差值不大。</p><p>　　圖3-41輥系的mises應(yīng)力圖</p><p>　　圖3-42輥系的mises應(yīng)力圖</p><p>　　3.6.2.2不同工況下的撓度差值比較</p><p>　　不同工況下兩輥的撓度差值如表3-1所示：</p>&

104、lt;p><b>　　表3-1</b></p><p>　　由此表可以看出不同的工況對(duì)于撓度的影響是很大的。彎輥力越大，工作輥的撓度值變化就越小，這說明在軋輥材料強(qiáng)度條件得以滿足的前提下施加盡可能大的彎輥力對(duì)于保持板形的良好是大有好處的；彎輥力的重要性亦正如第四種工況所示，工作輥的撓度差值遠(yuǎn)遠(yuǎn)大于了支承輥的差值，在這種情況下，工作輥的變型的不到很好的補(bǔ)償，板形的控制精度是很差的。當(dāng)單位

105、軋制壓力下降時(shí)，支承輥撓度差減小，工作輥撓度差增加，這也是比較容易理解的：對(duì)于工作輥而言，由于中間截面處的撓度在初始條件下小于輥頸處撓度，于是軋制壓力的減小使得軋輥中段的撓度變得更小，致使差值增加；對(duì)于支承輥而言，由于中間截面處的撓度在初始條件下是要大于輥頸處撓度的，所以撓度差值自然隨著軋輥中段的撓度變小而減小。</p><p>　　3.6.3 彈性壓扁曲線的計(jì)算分析</p><p>　　

106、對(duì)不同工況下兩輥之間的彈性壓扁量的分析如下表所示：</p><p><b>　　表3-2</b></p><p>　　分析表格中的數(shù)據(jù)，不難看出，當(dāng)彎輥力的大小足夠的時(shí)候，板寬邊緣的彈性壓扁量要大于中間截面的彈性壓扁量，而當(dāng)彎輥力較小時(shí)則恰恰相反。在實(shí)際軋制中，由于板帶材邊部壓下量大，壓應(yīng)力高，所以靠近邊部的彈性壓扁量大一些對(duì)于減小帶材邊部凸度是很有好處的。</

107、p><p>　　3.6.4 接觸面接觸應(yīng)力的分析</p><p>　　四輥軋機(jī)的支承輥和工作輥之間承載時(shí)有很大的接觸應(yīng)力，我們要對(duì)其進(jìn)行必要的校核計(jì)算。</p><p>　　對(duì)比各種工況下接觸應(yīng)力的分布情況可以知道，接觸應(yīng)力的最大值出現(xiàn)在接觸面靠近輥頸的位置(正常情況下)或者靠近軋輥中間截面的位置(彎輥力很小時(shí))。我們按照接觸壓應(yīng)力最大的一種情況進(jìn)行校核。</p&

108、gt;<p>　　在b=1000mm，J=1200KN的條件下，接觸應(yīng)力</p><p><b>　　（3-6）</b></p><p><b>　　最大反復(fù)切應(yīng)力</b></p><p>　　由《軋鋼機(jī)械》(第三版)表3-7查知支承輥材料的輥面硬度較低，HS=58~68，查表3-8知最大許用接觸應(yīng)力，，所以

109、加載條件是滿足接觸應(yīng)力強(qiáng)度要求的。</p><p><b>　　由公式推知</b></p><p>　　其中 ——加在接觸表面單位長度上的負(fù)荷；</p><p>　　、——相互接觸的兩個(gè)軋輥的半徑；</p><p>　　從圖形中可以看出，在除第四種工況的分析中，接觸壓力的變化并不是單調(diào)遞增的，在中間截面處，壓力稍大

110、一些的原因是施加邊界條件時(shí)約束方式的不盡完美。而在靠近輥頸的位置壓力值突然下降的原因則是由于在建模的時(shí)候接觸網(wǎng)格仍不算太細(xì)密使然。</p><p>　　3.6.5板帶凸度分析</p><p>　　由軋件表面的變形曲線可以看出，板帶材截面大體可以分為三個(gè)區(qū)域：在兩邊，厚度變化劇烈，稱為邊部減薄，這與板帶材邊部工作輥的彈性壓扁較小有關(guān)，而在板帶材中間90%以上的區(qū)域厚度變化較為平緩，大致具有二

111、次曲線的特征，該區(qū)域的形狀一般用中心厚度與邊部厚度之差表示，稱為板凸度或帶鋼凸度。</p><p>　　在b=800mm時(shí)的三種工況中，當(dāng)J=1200KN時(shí)，凸度；當(dāng)J=800KN時(shí)，凸度；當(dāng)J=0KN時(shí)，凸度?？梢姀澼伭Φ氖┘訉?duì)于降低板帶材的凸度是大有好處的。在強(qiáng)度允許的條件下，較大的彎輥力是有益于保證板帶的平直度的。</p><p><b>　　3.6.6綜合分析</b

112、></p><p>　　綜合各軋制過程中的影響因素不難發(fā)現(xiàn)，在軋制的過程中，如何保持合適的彎輥力與軋制力的搭配是很重要的。由于彎輥力的施加對(duì)于保持軋件的邊步板形良好具有很重要的作用，所以在軋制時(shí)盡可能使彎輥力大一些，軋制力也是如此，雖然較大的軋制力使板帶產(chǎn)生了較大的變形，但由于整個(gè)軋制過程還有后續(xù)軋機(jī)的精軋過程，所以再此不必過多考慮彎輥力和軋制力增加所帶來的不利影響。只要滿足強(qiáng)度要求即可。</p>

113、;<p><b>　　結(jié)論</b></p><p>　　本論文研究的目的是借助ANSYS軟件實(shí)現(xiàn)對(duì)軋制過程中影響因素的控制，并對(duì)工作輥的變形狀態(tài)進(jìn)行準(zhǔn)確的預(yù)報(bào)。論文從實(shí)際應(yīng)用的角度出發(fā)，研究中模型的建立和方法的應(yīng)用兼顧良好的仿真模擬和計(jì)算的快速、易實(shí)現(xiàn)性。由于精力和時(shí)間的限制，本論文難免有存在一些不足之處，比如在更高精度的建模以及約束的施加方面還存在著一些可以細(xì)化處理的環(huán)節(jié)，下

114、面對(duì)論文的主要研究分析結(jié)果作做一簡(jiǎn)單的總結(jié)：</p><p>　　論文介紹了四輥軋機(jī)輥系變形的常用理論和較為通用的相關(guān)模型，并著重介紹了本課題所采用的三位有限元建模思想。</p><p>　　論文介紹了本人應(yīng)用ANSYS通用有限元分析軟件建立模型以及分析輥系變形的過程。著重介紹了接觸問題的處理以及非線性力的加載方法。</p><p>　　論文從對(duì)工作輥和支承輥的撓曲

115、線、彈性壓扁、接觸壓力分布、軋件凸度圖形的計(jì)算分析中闡述了板寬和彎輥力的大小對(duì)于軋制過程的影響情況。</p><p><b>　　致謝</b></p><p>　　用一個(gè)學(xué)期的時(shí)間專注于一件事情，這是我在以前的學(xué)習(xí)中不曾預(yù)見過的，而且很可能是在以后的學(xué)習(xí)工作中也不會(huì)多次出現(xiàn)的機(jī)會(huì)，我很感謝在即將離開大學(xué)之前有這樣的機(jī)會(huì)可以讓我系統(tǒng)地綜合地學(xué)習(xí)運(yùn)用知識(shí)來解決工程實(shí)際問題

116、，學(xué)會(huì)做人做事做學(xué)問。</p><p>　　感謝我的指導(dǎo)老師徐致讓教授在整個(gè)畢業(yè)設(shè)計(jì)期間對(duì)我的指導(dǎo)和幫助，我的畢業(yè)設(shè)計(jì)中匯聚了徐老師濃濃的關(guān)懷，徐老師身上的學(xué)者風(fēng)范深深影響著我，使我明白了一個(gè)專業(yè)的機(jī)械人才應(yīng)該具備怎樣的素質(zhì)，謹(jǐn)祝徐老師工作順利，生活安康；感謝包家漢老師在ANSYS的應(yīng)用方面對(duì)我的指導(dǎo)，我為包老師熱心而細(xì)致的輔導(dǎo)所深深感動(dòng)，祝愿包老師在事業(yè)中取得更卓越的成就；感謝4年來機(jī)械學(xué)院所有的專業(yè)課老師對(duì)我

117、的培養(yǎng)，感謝4年來院領(lǐng)導(dǎo)的關(guān)懷和提攜。</p><p><b>　　參考文獻(xiàn)</b></p><p>　　1.徐致讓,薛家國.四輥軋機(jī)輥系變形分析的有限元方法.力學(xué)與實(shí)踐,2000, 22: 22~24</p><p>　　2.宗家福等.不銹鋼/碳鋼復(fù)合板帶可逆冷軋過程有限元模擬.軋鋼,2004, 1: 8~11</p><

118、p>　　3.徐建忠等.四輥軋機(jī)軋輥彈性變形解析模塊的開發(fā).軋鋼,2003,2: 8~11 </p><p>　　4.Guo.R-M.‘Development of a mathematical model for strip thickness profile’ Iron and Steel Engineering. Sept,1991.</p><p>　　5.Guo.R-M.‘S

119、tress analysis and life expectancy of rolling mill housings’ Iron and Steel Engineering. Aug. 1992:49~53 </p><p>　　6.曹新九.熱連軋精軋機(jī)組板形控制系統(tǒng)的研究.[學(xué)位論文].北京：北京科技大學(xué)，2002</p><p>　　7.V.B.金茲伯格. 高精度板帶軋制理論與實(shí)踐.