李楊，蔡業(yè)彬，張興芳

（1.廣東石油化工學(xué)院機(jī)電工程學(xué)院，廣東茂名525000；2.太原理工大學(xué)化學(xué)化工學(xué)院，山西太原030024）

焊接滾輪架包角對(duì)筒體影響的有限元分析

李楊1，2，蔡業(yè)彬1，張興芳2

（1.廣東石油化工學(xué)院機(jī)電工程學(xué)院，廣東茂名525000；2.太原理工大學(xué)化學(xué)化工學(xué)院，山西太原030024）

首先介紹了ANSYS求解接觸問(wèn)題的基本類(lèi)型和接觸單元的類(lèi)型，確定了焊接滾輪架與筒體的接觸是柔體—柔體的面—面接觸問(wèn)題。利用ANSYS軟件分別建立了包角取50°、60°、70°、80°、90°、100°、110°的滾輪和筒體的有限元幾何模型，分別進(jìn)行接觸分析，并對(duì)筒體的應(yīng)力和位移進(jìn)行分析對(duì)比。結(jié)果表明，若只考慮筒體的應(yīng)力和變形，焊接滾輪架包角應(yīng)選在60°～70°范圍內(nèi)，對(duì)包角的選取具有一定的指導(dǎo)價(jià)值。

ANSYS；焊接滾輪架；包角；筒體

焊接滾輪架是一種常用的焊接輔助裝置，在焊接生產(chǎn)過(guò)程中被廣泛應(yīng)用。主要借助焊件與主動(dòng)滾輪間的摩擦力帶動(dòng)圓筒類(lèi)焊件旋轉(zhuǎn)，與焊接操作機(jī)、焊接變位機(jī)等組合可實(shí)現(xiàn)對(duì)焊件的內(nèi)外環(huán)縫、角焊縫和內(nèi)外縱縫的自動(dòng)焊接。焊接滾輪架與筒體的結(jié)構(gòu)簡(jiǎn)圖如圖1所示，其中α表示兩滾輪對(duì)筒體的包角。近年來(lái)，對(duì)焊接滾輪架的研究主要集中在滾輪架的軸向防竄動(dòng)和自調(diào)式滾輪架的研制［1-4］，對(duì)包角α的研究也僅限于傳動(dòng)過(guò)程，即當(dāng)α＜45°或α＞120°時(shí)，都會(huì)使傳動(dòng)過(guò)程不穩(wěn)定［5］，并未涉及靜態(tài)分析。因此，分析不同包角α對(duì)筒體的影響，對(duì)焊接滾輪架包角的選取具有重要意義。

筒體和滾輪屬于接觸問(wèn)題，接觸問(wèn)題是一種狀態(tài)非線性行為，傳統(tǒng)方法是運(yùn)用赫茲接觸理論進(jìn)行求解計(jì)算，但這種方法計(jì)算精度不高而且局限性較大［6］。隨著有限元分析軟件的發(fā)展，ANSYS可以解決所有赫茲接觸問(wèn)題而且計(jì)算精度較高，應(yīng)用廣泛。因此，本文利用ANSYS分析焊接滾輪架包角對(duì)筒體的影響。

圖1　焊接滾輪架與筒體的結(jié)構(gòu)簡(jiǎn)圖

1　ANSYS接觸分析

1.1接觸問(wèn)題類(lèi)型

接觸問(wèn)題有兩種基本類(lèi)型：剛體—柔體的接觸和柔體—柔體的接觸。一般情況下，一種軟材料和一種硬材料接觸可以被假定為剛體—柔體的接觸，例如許多金屬的成形問(wèn)題；柔體—柔體的接觸是一種更普遍的類(lèi)型，兩個(gè)接觸體都是變形體［7］。筒體和滾輪的接觸就是柔體—柔體的接觸。

1.2接觸單元類(lèi)型

ANSYS支持三種接觸方式：點(diǎn)—點(diǎn)接觸、點(diǎn)—面接觸和面—面接觸。每種接觸方式可選用不同的接觸單元以適應(yīng)于不同類(lèi)型的問(wèn)題。

點(diǎn)—點(diǎn)接觸單元主要用于點(diǎn)—點(diǎn)的接觸行為，使用這些單元，需預(yù)先知道接觸位置，而且只適用于有較小相對(duì)滑動(dòng)的接觸面。對(duì)于大多數(shù)點(diǎn)—點(diǎn)接觸問(wèn)題，CONTA178單元是最好的選擇。

點(diǎn)—面接觸單元主要用于點(diǎn)—面的接觸行為，使用這些單元，不需預(yù)先知道確切的接觸位置，允許有大的相對(duì)滑動(dòng)和大變形，且接觸件間不需要網(wǎng)格保持一致。點(diǎn)—面接觸單元有CONTA175、CONTA26、CONTA48和CONTA49，通常推薦使用CONTA175，因?yàn)樗褂玫膯卧獢?shù)較少，且適用于大的摩擦滑動(dòng)問(wèn)題。

面—面接觸單元主要用于面—面的接觸行為，這些單元通過(guò)接觸面和目標(biāo)面形成接觸對(duì)來(lái)模擬實(shí)體間的接觸，且接觸對(duì)共享一個(gè)實(shí)常號(hào)。目標(biāo)面單元有TARGE169（2D）和TARGE170（3D），接觸面單元有CONTA171（2D）、CONTA172（2D）、CONTA173（3D）、CONTA174（3D），這些單元可用于處理裝配零件間的過(guò)盈配合等問(wèn)題［7-9］。筒體和滾輪的接觸選用面—面接觸單元模擬。

2　問(wèn)題描述

筒體材料為Q345R，焊接滾輪架材料為ZG35 CrMnSi，具體參數(shù)如表1所示。主要對(duì)不同包角α下的筒體進(jìn)行有限元分析，討論包角α對(duì)筒體的變形和應(yīng)力的影響。JB/T 9187-1999《焊接滾輪架》中規(guī)定：焊接滾輪架每對(duì)滾輪的中心距必須能根據(jù)筒體類(lèi)工件的直徑作相應(yīng)的調(diào)整，保證兩滾輪對(duì)筒體的包角大于45°，小于110°.因此，包角α分別取50°、60°、70°、80°、90°、100°、110°進(jìn)行分析。

表1　筒體和滾輪的參數(shù)表

3　分析步驟

3.1定義單元和材料

由于滾輪和筒體的結(jié)構(gòu)具有對(duì)稱(chēng)性，而且僅受重力影響，故取模型的二分之一進(jìn)行有限元分析。采用SOLID45單元，該單元主要用于三維實(shí)體結(jié)構(gòu)模型，具有8個(gè)節(jié)點(diǎn)，且每個(gè)節(jié)點(diǎn)有x、y、z三個(gè)方向的3個(gè)自由度。同時(shí)，具有塑性、徐變、膨脹、應(yīng)力強(qiáng)化、大變形和大應(yīng)變的能力。

選用各項(xiàng)同性材料，分別定義筒體和滾輪的材料常數(shù)，筒體和滾輪所用材料的材料特性如表2所示。

表2　筒體和滾輪的材料特性表

3.2創(chuàng)建幾何模型和網(wǎng)格劃分

根據(jù)表1中筒體和滾輪的尺寸建立幾何模型（如圖2所示），分別估計(jì)筒體和滾輪的接觸區(qū)域，并用工作平面將筒體切割為三部分，為網(wǎng)格劃分和定義接觸對(duì)做準(zhǔn)備。分別指定幾何模型中各線段的單元等分?jǐn)?shù)，然后用SWEEP方式進(jìn)行網(wǎng)格劃分，為保證計(jì)算精度同時(shí)減少系統(tǒng)的資源消耗，對(duì)筒體與滾輪接觸區(qū)域進(jìn)行網(wǎng)格細(xì)化（如圖3所示）。

圖2　筒體和滾輪架的幾何模型（以α=60°為例）

圖3　筒體和滾輪架的有限元網(wǎng)格模型（以α=60°為例）

3.3定義接觸單元和接觸對(duì)

接觸面采用CONTA174單元，目標(biāo)面采用TARGE170單元。根據(jù)接觸面和目標(biāo)面的選擇準(zhǔn)則［10］，由于主要分析筒體的變形和應(yīng)力，故將筒體接觸部分視作接觸面，將滾輪接觸部分視作目標(biāo)面，這兩個(gè)面合起來(lái)稱(chēng)為接觸對(duì)（如圖4所示）。

圖4　筒體和滾輪架的接觸對(duì)（以α=60°為例）

3.4設(shè)置邊界條件

滾輪下端各節(jié)點(diǎn)施加全部約束，筒體對(duì)稱(chēng)面上各節(jié)點(diǎn)施加對(duì)稱(chēng)約束，同時(shí)施加Y方向重力加速度g=9.8 m/s2（如圖5所示）。

圖5　筒體和滾輪架的邊界條件（以α=60°為例）

3.5求解及結(jié)果分析

設(shè)置求解選項(xiàng)，打開(kāi)自動(dòng)時(shí)間步長(zhǎng)和牛頓拉普森選項(xiàng)，開(kāi)啟線性搜索方式。選中所有節(jié)點(diǎn)和單元，求解當(dāng)前載荷。求解完成后，可通過(guò)ANSYS的后處理功能顯示筒體的等效應(yīng)力云圖（如圖6所示，由于包角不同時(shí)等效應(yīng)力云圖相似，故只以α=60°為例）和總位移分布圖（如圖7所示），可以直觀地看出筒體的受力與變形，并依此可繪制出不同包角下筒體的最大應(yīng)力值和最大位移值的折線圖（如圖8所示）。

圖6　筒體的等效應(yīng)力云圖（以α=60°為例）

（續(xù)下圖）

圖7　筒體的總位移分布圖

圖8　不同中心包角下筒體的最大應(yīng)力值和最大位移值

由圖6和圖8可知，筒體的最大應(yīng)力出現(xiàn)在筒體和滾輪的接觸區(qū)域，但接觸處的相鄰區(qū)域應(yīng)力值相對(duì)較小，越靠近接觸區(qū)域應(yīng)力值就越大，可見(jiàn)筒體和滾輪的接觸區(qū)域存在較大的擠壓和磨損。而且，隨著包角α的增大，筒體的最大應(yīng)力值先減小后增大，且當(dāng)α=60°時(shí)，最大應(yīng)力值達(dá)到最小。

由圖7和圖8可知，筒體的最大位移也出現(xiàn)在筒體和滾輪的接觸區(qū)域，接觸處的相鄰區(qū)域位移值也相對(duì)較高。而且，隨著包角α的增大，筒體的最大位移值先減小后增大再減小，且當(dāng)α=70°時(shí)，最大位移值達(dá)到最??；當(dāng)α＞90°時(shí)，最大位移值急劇增大，且α =100°時(shí)達(dá)到最大。此外，隨著包角α的增大，筒體的最大位移值的范圍也在增大，即滾輪對(duì)筒體變形的影響在增大，遠(yuǎn)離接觸處的區(qū)域的位移值依然相對(duì)較高。

綜上比較分析，若只考慮滾輪對(duì)筒體變形和應(yīng)力的影響，滾輪架的包角α選在60°～70°范圍內(nèi)為宜，因?yàn)樵谶@個(gè)范圍內(nèi)，筒體的最大變形和最大應(yīng)力值相對(duì)較小。

4　結(jié)束語(yǔ)

綜上分析得出：

（1）利用ANSYS軟件，建立了不同包角的筒體和滾輪的有限元模型，進(jìn)行了接觸分析。

（2）筒體的最大應(yīng)力和最大位移都出現(xiàn)在筒體和滾輪的接觸區(qū)域，接觸處的相鄰區(qū)域應(yīng)力值相對(duì)較小，但位移值依然較大。

（3）隨著包角α的增大，筒體的最大應(yīng)力值先減小后增大，且當(dāng)α=60°時(shí)，最大應(yīng)力值達(dá)到最小；而最大位移值先減小后增大再減小，且當(dāng)α=70°時(shí)，最大位移值達(dá)到最小。

（4）隨著包角α的增大，遠(yuǎn)離接觸處的區(qū)域的位移值依然相對(duì)較高。

（5）若只考慮滾輪對(duì)筒體變形和應(yīng)力的影響，滾輪架的包角α應(yīng)選在60°～70°范圍內(nèi)，對(duì)焊接滾輪架中心包角的選取具有指導(dǎo)意義。

［1］羅生梅，陳利.焊接滾輪架上工件軸向竄動(dòng)防竄機(jī)理研究［J］.機(jī)床與液壓，2007，35（8）：74-75.

［2］顏平錄.如何制造焊件軸向零竄動(dòng)的焊接滾輪架［J］.機(jī)械工人（熱加工），2003，（3）：6-7.

［3］陳利.智能防軸向竄動(dòng)滾輪架系統(tǒng)的研究與開(kāi)發(fā)［D］.蘭州:蘭州理工大學(xué)，2007.

［4］李廣周，劉玉虎，王剛.基于虛擬樣機(jī)技術(shù)的自調(diào)心焊接滾輪架的研制與應(yīng)用［J］.天津科技，2015，42（10）：52-54.

［5］于有生，張旭明，鄒澤明.自調(diào)式焊接滾輪架的一個(gè)幾何約束條件［J］.焊接技術(shù)，2001，30（4）：31-32.

［6］李峰，趙杰.應(yīng)用ANSYS求解赫茲接觸問(wèn)題［J］.機(jī)電產(chǎn)品開(kāi)發(fā)與創(chuàng)新，2006，19（6）：76-78.

［7］劉相新，孟憲頤.ANSYS基礎(chǔ)與應(yīng)用教程［M］.北京：科學(xué)出版社，2006.

［8］錢(qián)俊梅，江曉紅，仲小冬，等.淺談基于ANSYS軟件的接觸分析問(wèn)題［J］.煤礦機(jī)械，2007，27（7）：62-64.

［9］李一耕，陳錦江，龍超，等.基于ANSYS的球軸承赫茲接觸問(wèn)題有限元分析［J］.機(jī)械研究與應(yīng)用，2011，（5）：35-37.

［10］姜衛(wèi)忠.壓力容器筒體與補(bǔ)強(qiáng)圈間的接觸特性［J］.壓力容器，2003，（02）：20-23.

Finite Element Analysis for the Influence of the Wrap Angle of Welding Roller Bed on the Barrel

LI Yang1，2，CAI Ye-bin1，ZHANG Xing-fang2
（1.College of Mechanical and Electrical Engineering，Guangdong University of Petrochemical Technology，Maoming Guangdong 525000，China；2.College of Chemistry and Chemical Engineering，TaiyuanUniversity of Technology，Taiyuan 030024，China）

The fundamental types and types of contact elements of ANSYS software in solution of contact problem were introduced.Then the contact of welding roller bed and barrel was determined to the contact problem of flexible body－to－flexible body surface－to－surface.And this paper used analysis software ANSYS to establish the finite element models of the roller and barrel of different wrap angle and executed the contact analysis respectively．Through compared and analyzed the stress and displacement of barrel，results showed that the wrap angle of welding roller bed is better in a range between 60°and 70°if only consider the stress and deformation of the barrel.This is valuable in the select of wrap angle.

ANSYS；welding roller bed；wrap angle；barrel

TG431；TH123.4

A

1672-545X（2016）08-0007-04

2016-05-16

廣東省自然科學(xué)基金（編號(hào)：9152500002000003）；廣東省教育部產(chǎn)學(xué)研結(jié)合項(xiàng)目（編號(hào)：2010B090400237）；廣東省教育廳科技創(chuàng)新項(xiàng)目（編號(hào)：2012KJCX0076）

李楊（1991-），女，山西臨汾人，碩士研究生，研究方向?yàn)榛ぱb備制造與檢測(cè)。