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

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

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

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

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

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

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

7、　板帶材的板形可以通過(guò)對(duì)以下的軋輥的兩類(lèi)撓曲進(jìn)行疊加來(lái)確定：</p><p>　　由于軋制力引起的彎曲力使軋輥產(chǎn)生的撓曲；</p><p>　　由于軋制力引起的剪切力使軋輥產(chǎn)生的撓曲。</p><p>　　由彎曲力產(chǎn)生的撓曲可由如下的微分方程描述：</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—— 軋制力作用點(diǎn)的間距；</p><p>　　——

9、帶材的寬度。 </p><p>　　在軋輥與帶材的接觸區(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)力呈均勻分布條件下獲得的。軋輥和帶材

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

12、/b></p><p>　　由彎曲力和剪切力產(chǎn)生的軋輥總的撓度為：</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>　　斯通(Sto

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

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

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

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

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

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

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

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

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

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

23、<p>　　因?yàn)檐垯C(jī)是關(guān)于軋輥中心對(duì)稱(chēng)的，所以計(jì)算過(guò)程僅需考慮軋輥的一半。所用數(shù)值方法是將軋輥切分成m個(gè)單元，如圖1-4所示，并將軋輥的分布載荷代之以施加在每個(gè)單元中心的集中載荷。因?yàn)檐埣挾刃∮谳伾黹L(zhǎng)度，所以與軋件接觸的軋輥單元數(shù)n小于m。</p><p>　　支撐輥的變形。在分割法模型中，支撐輥的變形可以表示為在位置i處的軋輥表面的垂直位移：</p><p><b>

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

25、gt;<p>　　——支撐輥的剛體移動(dòng)。</p><p>　　工作輥的變形工作輥?zhàn)冃慰杀硎緸楣ぷ鬏伇砻嬖谖蝗玦處的垂直位移：</p><p><b>　　(1-13)</b></p><p>　　式中 ——工作輥和軋件間第i個(gè)單元上每單位寬度的載荷；</p><p>　　——工作輥的影響系數(shù)；</

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

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

28、 </p><p><b>　　(1-14)</b></p><p><b>　　當(dāng)時(shí)：</b></p><p><b>　　(1-15)</b></p><p>　　將工作輥和支撐輥的參數(shù)對(duì)應(yīng)值代入方程1-14和1-15，可得系數(shù)和的值。</p>&

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

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

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

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

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

34、(Hollander)和萊茵(Reinen)曾經(jīng)對(duì)另一種模型做了闡述。</p><p>　　分割梁模型的局限性 </p><p>　　分割梁撓曲模型的研究發(fā)展對(duì)于在軋制過(guò)程中提高帶材板形的模擬能力邁出了很大一步。然而，這類(lèi)模型也有自身的不足，因?yàn)榉指盍耗Ｐ褪墙⒃诩僭O(shè)在軋制力作用下工作輥和支承輥完全接觸之上的。實(shí)際上，當(dāng)采用特殊的輥型系統(tǒng)如CVC軋輥、

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

36、 </p><p>　　1.2.5 有限元分析理論 </p><p>　　在有限元分析中采用的是矩陣結(jié)構(gòu)分析方法。這種方法采用了一種直接的物理方法來(lái)建立和求解梁及框架結(jié)構(gòu)問(wèn)題。</p><p&

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

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

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

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

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

42、結(jié)構(gòu)的幾何形狀及模型可能移動(dòng)的方向的個(gè)數(shù)，即自由度。</p><p>　　簡(jiǎn)單的彈簧單元可承受張力并有兩個(gè)自由度，因?yàn)榱旱膬啥丝梢匝貜椈傻闹鬏S方向自由移動(dòng)。然而，具有更多自由度的單元體可以作其他運(yùn)動(dòng)如彎曲運(yùn)動(dòng)。單元類(lèi)型的選擇也決定了在一個(gè)模型中的自由度。</p><p>　　(2) 應(yīng)變-位移關(guān)系的建立。應(yīng)變-位移關(guān)系可通常表示為：</p><p><b>

43、;　　(1-21)</b></p><p>　　式中 ——單元應(yīng)變矩陣；</p><p>　　——單元應(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)力矩陣；<

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

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

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

47、]構(gòu)成的總剛度矩陣[K]和由各個(gè)節(jié)點(diǎn)的力矢量{P}構(gòu)成的總的力矢量{R}的聯(lián)立。</p><p>　　聯(lián)立方法的根據(jù)是由于每個(gè)節(jié)點(diǎn)的相互聯(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)求解未

48、知參數(shù)?？紤]物體的幾何和力的邊界條件，求解代數(shù)方程1-16，可得所有未知的位移。在線性平衡問(wèn)題中，可直接應(yīng)用矩陣代數(shù)技術(shù)來(lái)求解：</p><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

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

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

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

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

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

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

55、 </p><p>　　式中 、——在有限元網(wǎng)格中，當(dāng)網(wǎng)格數(shù)為n和∞時(shí)分別對(duì)應(yīng)的中心凸度；</p><p>　　——帶材的出口厚度。</p><p>　　因?yàn)橹行耐苟鹊南鄬?duì)誤差隨著網(wǎng)格數(shù)量的增加而減小，所以我們可以根據(jù)指數(shù)定律來(lái)確定的值。</p>

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

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

58、0、4600和5200個(gè)，如圖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é)構(gòu)離散成有限個(gè)單元，并在每一個(gè)單元中設(shè)定有限個(gè)節(jié)點(diǎn)，將連續(xù)體看作是只在節(jié)點(diǎn)處相連接的一組單元的集合體；同時(shí)選定場(chǎng)函數(shù)的節(jié)點(diǎn)值作為基本未知量，并在每一單元中假設(shè)一近似差值函數(shù)以表示單元中場(chǎng)函數(shù)的分布規(guī)律；進(jìn)而利用力學(xué)中的某些變分原理去建立用于求解節(jié)點(diǎn)未知量的有限元法方程，從而將一個(gè)連續(xù)域中的無(wú)限自由度問(wèn)題化為離散域中的有限自由度問(wèn)題。一經(jīng)求解就可以利用解得的節(jié)點(diǎn)值和設(shè)定的插值函數(shù)確定單元上以至整個(gè)集合體上的

60、場(chǎng)函數(shù)。有限元求解程序的內(nèi)部過(guò)程可從圖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年來(lái)，隨著計(jì)算機(jī)的飛速發(fā)展和廣泛應(yīng)用，各種行之有效的數(shù)值計(jì)算方法得到了巨大的發(fā)展。而有限元方法則是計(jì)算機(jī)誕生以后，在計(jì)算數(shù)學(xué)、計(jì)算力學(xué)、和計(jì)算工程科學(xué)領(lǐng)域里誕生的最

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

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

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

65、三個(gè)部分：前處理模塊、分析計(jì)算模塊和后處理模塊。</p><p><b>　　前處理模塊</b></p><p>　　它為用戶提供了一個(gè)強(qiáng)大的實(shí)體建模及網(wǎng)絡(luò)劃分工具，用戶可以方便地構(gòu)造有限元模型，軟件提供了100種以上的單元類(lèi)型，用來(lái)模擬工程中的各種結(jié)構(gòu)和材料。</p><p><b>　　實(shí)體建模</b></p&g

66、t;<p><b>　　參數(shù)化建模</b></p><p><b>　　體素庫(kù)及布爾運(yùn)算</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)格劃分

67、、自適應(yīng)網(wǎng)格劃分。</p><p>　　復(fù)雜幾何體Sweep映射網(wǎng)格生成。</p><p>　　六面體向四面體自動(dòng)過(guò)渡網(wǎng)格：金字塔形。</p><p><b>　　邊界層網(wǎng)格劃分</b></p><p>　　在幾何模型或FE模型上加載：點(diǎn)載荷、分布載荷、體載荷、函數(shù)載荷。</p><p>　　可擴(kuò)

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

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

70、/p><p>　　圖形：云圖、等值線、矢量顯示、粒子流跡顯示、切片、透明及半透明顯示、紋理。</p><p>　　鋼筋混凝土單元可顯示單元內(nèi)的鋼筋、開(kāi)裂情況以及壓碎部位。</p><p>　　梁、管、板、復(fù)合材料單元及結(jié)果按實(shí)際形狀顯示，顯示橫截面結(jié)果；顯示梁?jiǎn)卧獜澗貓D。</p><p>　　顯示優(yōu)化靈敏度及優(yōu)化變量曲線。</p>

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

72、減、積分、微分等計(jì)算。</p><p>　　顯示沿任意路徑的結(jié)果曲線，并可進(jìn)行路徑的數(shù)學(xué)計(jì)算。</p><p>　　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ù)庫(kù)的一體化大型FEA分析軟件</p><p>

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

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

75、面和程序結(jié)構(gòu)，交互式的前后處理和圖形軟件，大大的減輕了用戶在實(shí)際工程問(wèn)題中創(chuàng)建模型、有限元求解以及結(jié)果分析和評(píng)價(jià)的工作量。他的統(tǒng)一集中式的數(shù)據(jù)庫(kù)保證了個(gè)模塊之間的有效可靠的集成，并實(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ì)算。本課題要

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

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

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

79、可以因模型的簡(jiǎn)化而縮短一半的運(yùn)行時(shí)間。對(duì)我們的設(shè)計(jì)是很有利的。</p><p>　　2.1.2 模型設(shè)計(jì)</p><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,,,, &

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,, !!!與點(diǎn)4的坐標(biāo)值相同，但代表的是不同的體

83、上的點(diǎn)</p><p>　　K,1002,800,-630,, !!!與點(diǎn)5的坐標(biāo)值相同，但代表的是不同的體上的點(diǎn)</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>　　在進(jìn)行此部分的建模示要特別注意k1001、k1002兩點(diǎn)的定義，單從坐標(biāo)值看來(lái)，似乎這兩個(gè)點(diǎn)的定義是與點(diǎn)4和點(diǎn)5重復(fù)的，實(shí)則不然，此兩點(diǎn)連線所生成的直線表示的是兩輥之間初始接觸的那條線。在ANSYS中，若將同一條線賦予兩個(gè)實(shí)體，那

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

87、是一定要事先明確的，因?yàn)椴煌膭澐址绞剿蟮那捌谀Ｐ蛣?chuàng)建方式也是不同的。上述程序的特點(diǎn)是創(chuàng)建出的實(shí)體模型便于通過(guò)掃掠的方式劃分單元網(wǎng)格。相比較而言掃掠獲得的單元形狀較為規(guī)范美觀，在相當(dāng)精度的情況下運(yùn)算起來(lái)較為省時(shí)，而且便于后處理操作。故在以后的分析中，均采用掃掠出的模型進(jìn)行分析，采用自由網(wǎng)格劃分的模型另附于說(shuō)明書(shū)后。</p><p>　　2.2接觸區(qū)域的劃分</p><p>　　2.2.

88、1 設(shè)計(jì)原理</p><p>　　2.2.1.1軋件與工作輥的接觸區(qū)</p><p>　　按照適用于熱軋帶鋼的Sims公式，可知在軋件與工作輥的接觸區(qū)中存在著前滑和后滑兩個(gè)變形區(qū)。兩個(gè)區(qū)域由中性面分開(kāi)，各自滿足不同的公式的作用情況。所以要對(duì)這兩個(gè)區(qū)域進(jìn)行劃分。</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>　　特殊地，當(dāng)θ=0時(shí)，</p><p>　　= 66.75Mpa</p><p>　　所求的便是出口處的單位壓力</p><p>　　入口至中性面之間的單位壓力</p&g

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

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

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

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

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

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

98、><p>　　lstr,71,75!!!連接邊界點(diǎn)，生成邊界線</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>　　該過(guò)程對(duì)應(yīng)的模型圖如圖2所示。</p><p>　　圖2-2 模型的線輪廓圖</p><p>　　上一步驟中的難點(diǎn)在于軋制區(qū)的劃分，因?yàn)檐堉茀^(qū)是

102、曲面，所以線的創(chuàng)建有一定的難度。完成對(duì)線的創(chuàng)建以后，采用程序提供的面劃分指令對(duì)接觸面進(jìn)行創(chuàng)建，此過(guò)程的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個(gè)體塊粘在一起</p><p><b>　　vglue,all</b></p><p><b>　　alls</b></p>&l

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

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

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

108、節(jié)點(diǎn)數(shù)較之20節(jié)點(diǎn)的SOLID186等單元較少，節(jié)約了大量的運(yùn)算時(shí)間)。</p><p>　　在軋制區(qū)域我選用了SHELL63單元對(duì)加載面進(jìn)行劃分，SHELL63單元是具有8節(jié)點(diǎn)的彈性殼結(jié)構(gòu)單元，可以很好地和SOLID185單元形成匹配。注意在選擇殼單元進(jìn)行網(wǎng)格劃分的時(shí)候，要定義適當(dāng)?shù)膯卧穸?，并且使加載區(qū)單元?jiǎng)澐值乇M可能細(xì)致一些，因?yàn)樵谶M(jìn)行加載時(shí)程序要從模型中獲得該區(qū)域單元中心的坐標(biāo)值作為參數(shù)參加計(jì)算分析。單元

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

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

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

112、中工作輥和支承輥的材料選取有所不同，支承輥采用鑄鋼，工作輥采用合金鑄鐵作為材料。由《機(jī)械設(shè)計(jì)課程設(shè)計(jì)手冊(cè)》可以查出兩種材料的不同的材料性質(zhì)：</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>　　將兩種不同的材料分別定義于兩軋輥，然后再定義接觸對(duì)。所有單元?jiǎng)?chuàng)建完成以后的模型如圖2-3所示。</p><p>　　圖2-3單元?jiǎng)澐趾蟮哪Ｐ蛨D</p>

114、<p>　　2.3.3 模型設(shè)計(jì)</p><p>　　用于生成軋制區(qū)殼單元的程序語(yǔ)句如下：</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.用于生成源面殼單元的程序語(yǔ)句如下：</p><p>　　type,4!!!殼單元的單元和實(shí)常數(shù)類(lèi)型</p><

117、;p><b>　　real,3</b></p><p>　　asel,s,,,6,24,6!!!支承輥的源面殼單元?jiǎng)澐?lt;/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!!!設(shè)置四邊形網(wǎng)格尺寸</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!!!工作輥的源面殼單元?jiǎng)澐?lt;/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