基于ANSYS的2根拉索平行鋼絲接觸問題有限元分析

陳 璨,李 紅

(1.四川外國語學(xué)院重慶南方翻譯學(xué)院管理學(xué)院,重慶 401120;2.內(nèi)蒙古科技大學(xué)建筑與土木工程學(xué)院,內(nèi)蒙古 包頭 014010)

基于ANSYS的2根拉索平行鋼絲接觸問題有限元分析

陳 璨1,李 紅2

(1.四川外國語學(xué)院重慶南方翻譯學(xué)院管理學(xué)院,重慶 401120;2.內(nèi)蒙古科技大學(xué)建筑與土木工程學(xué)院,內(nèi)蒙古 包頭 014010)

由于2根鋼絲接觸問題可以作為整個拉索斷面平行鋼絲間接觸問題的簡化,故在用ANSYS有限元方法分析整個拉索斷面平行鋼絲排列模型的接觸問題之前,運(yùn)用經(jīng)典理論解和ANSYS有限元方法對2根拉索平行鋼絲間接觸問題進(jìn)行了較為詳盡的分析,確定出了兩根鋼絲間的接觸寬度、在接觸面上產(chǎn)生的最大接觸應(yīng)力及分布等主要參數(shù).同時,將2種分析方法得出的結(jié)果進(jìn)行了對比,表明用ANSYS有限元方法求解接觸問題是可行的,其求解的準(zhǔn)確性和精度能夠滿足要求.

拉索平行鋼絲;接觸問題;ANSYS;有限元分析

0 引 言

近年來,平行鋼絲斜拉索因其獨(dú)特的優(yōu)點(diǎn)而廣泛應(yīng)用在斜拉橋的承重上.平行鋼絲斜拉索內(nèi)鋼絲斷面排列情況是復(fù)雜多樣的(見圖1),這對分析拉索因鋼絲扭絞產(chǎn)生的鋼絲間正壓力的傳遞與分布帶來了相當(dāng)大的難度[1],目前尚不能進(jìn)行精確的理論計算或參數(shù)化編程計算.由于2根鋼絲接觸問題是拉索平行鋼絲間接觸問題的最簡化形式,只有將2根鋼絲接觸問題分析準(zhǔn)確之后,才能在此基礎(chǔ)上對多根鋼絲間接觸問題進(jìn)行分析,或者對一根鋼絲外包n圈鋼絲簡化模型接觸問題進(jìn)行分析,進(jìn)而對整個拉索斷面平行鋼絲排列模型的接觸問題進(jìn)行分析.兩根鋼絲接觸問題分析結(jié)果作為最基本的結(jié)果,也可用來檢驗(yàn)對復(fù)雜問題的解決是否有效.故本研究運(yùn)用經(jīng)典理論解和ANSYS有限元方法對2根拉索平行鋼絲間接觸問題進(jìn)行較為詳盡地分析.

1 2根鋼絲接觸問題的描述及求解

1.1 問題描述

2根拉索平行鋼絲接觸問題描述如下:2根單位長度(這里取1 mm)、半徑為3.5 mm的斜拉索鍍鋅平行鋼絲(見圖2),當(dāng)在平行鋼絲的軸線上施加均布荷載P=F/l=4 000 N/mm時,2根鋼絲均會在壓力的作用下產(chǎn)生相應(yīng)的應(yīng)力應(yīng)變,特別是在2根鋼絲的接觸表面將會發(fā)生局部接觸變形.那么,其接觸的面積以及在接觸面上產(chǎn)生的最大接觸應(yīng)力的大小可以根據(jù)接觸問題經(jīng)典理論解和有限元軟件來求解.其中,平行鋼絲的彈性模量 E1=E2=E=1.95×105MPa,v1=v2=0.3.

圖2 單位長度的2根鋼絲的接觸

1.2 經(jīng)典理論解

經(jīng)典理論認(rèn)為:對于2個長為 l,半徑分別為R1、R2的圓柱體的接觸(見圖3),當(dāng)沒有任何壓力作用于它們上時,2圓柱體只在其切平面上的一條線上有所接觸.當(dāng)沿圓柱體軸線上的均布荷載 P=F/l作用于圓柱體上時,2圓柱體間將形成面積為2 b×L的長方形接觸面.其赫茲接觸半寬b、最大接觸應(yīng)力 σHmax(接觸面上的最大單位壓力q0稱為最大接觸應(yīng)力 σHmax)、接觸相對位移 δ可按下列公式[3]計算:

圖3 2圓柱體接觸

接觸半寬b為,

接觸面上的最大接觸應(yīng)力 σHmax為,

接觸相對位移δ為,

其中,E1、E2為2圓柱體材料各自的彈性模量;v1、v2分別為2圓柱體材料各自的泊松比.

當(dāng)2圓體的材料相同,且泊松比均為0.3,即E1=E2=E,v1=v2=0.3時,接觸半寬b為,

接觸面上的最大接觸應(yīng)力 σHmax為,

接觸面上產(chǎn)生的接觸應(yīng)力的分布曲線近似為平滑的拋物線(見圖3).

由于平行鋼絲的幾何特性,可近似求解2根鋼絲的接觸問題[3].根據(jù)式(4)、(5)得,在荷載 P=F/l=4 000 N/mm的作用下,2根鋼絲間的接觸半寬b=0.288367536 mm;最大 接 觸 應(yīng)力 σHmax=8 824.793644MPa.

1.3 ANSYS有限元解

1.3.1 有限元模型的建立.

由于實(shí)體單元建模劃分網(wǎng)格時網(wǎng)格較多,故可將該問題簡化為平面問題,并建立平面模型[4].所建模型的尺寸、材料參數(shù)和力學(xué)參數(shù)為:模型厚度取單位厚度(1 mm),圓的半徑取3.5 mm,所加的荷載為4 000 N,材料的彈模取1.95×105MPa,泊松比取0.3.

在對模型進(jìn)行網(wǎng)格劃分時,為了確保計算精度,將模型劃分成四邊形網(wǎng)格.劃分網(wǎng)格后模型的單元數(shù)為640,如圖4所示(顯示1mm厚度),圖5為接觸區(qū)域的網(wǎng)格劃分情況.

網(wǎng)格劃分完畢后,將2圓面接觸區(qū)域建立接觸,將此接觸處理為面—面接觸.同時,指定目標(biāo)面和接觸面,定義目標(biāo)面、接觸面并生成接觸單元(見圖6).其中,用ESURE命令建立接觸對,用TARGE169單元在目標(biāo)面上生成目標(biāo)單元,用CONTA171單元在接觸面上生成接觸單元.建立接觸后,模型單元數(shù)增加到了672個.接著,在上面一個圓面頂部的節(jié)點(diǎn)施加向下的豎向荷載F=4 000 N,同時,在下面一個圓面底部處施加X、Y方向的約束,如圖7所示.

1.3.2 結(jié)果與分析.

合理地選擇單元的接觸剛度對于求解的收斂十分重要,即在ANSYS有限元分析程序中,需合理地選擇罰剛度比例系數(shù)(FKN).FKN取值通常將在0.01～10之間,可嘗試去尋找一個適當(dāng)?shù)腇KN值,起始估計值FKN=0.1.FKN 分別取0.1、1、10、100.圓柱體材料的彈模為1.95×105MPa,對應(yīng)的接觸剛度分別為 1.95×104、1.95×105、1.95×106、1.95×107N/mm2.以最大接觸應(yīng)力為監(jiān)控變量,分析結(jié)果見表1.

表1 不同接觸剛度,2根鋼絲模型ANSYS方法分析結(jié)果與理論值對比

表1中,FKN為罰剛度比例系數(shù);K接觸為接觸剛度;σAmax為用ANSYS有限元方法得出的最大接觸應(yīng)力;σHmax為用赫茲經(jīng)典理論得出的最大接觸應(yīng)力;bA為用ANSYS有限元方法得出的接觸半寬;bH為用赫茲經(jīng)典理論得出的接觸半寬.用ANSYS有限元方法求解時,接觸半寬bA可以通過公式(6)[5]得出,

式中,X′為處于接觸狀態(tài)的點(diǎn)中X坐標(biāo)最大的點(diǎn)的X坐標(biāo)值,U′X為該節(jié)點(diǎn)的X方向位移.

由表1可以看出,當(dāng)FKN的值超過10時,結(jié)果不再發(fā)生顯著的改變,即10為尋到的合適的FKN值.同時,將最大接觸應(yīng)力、接觸半寬的ANSYS解與經(jīng)典理論解對比可知:由ANSYS方法得出的結(jié)果和經(jīng)典理論得出的結(jié)果非常接近,這表明用ANSYS方法求解接觸問題是可行的,其求解的準(zhǔn)確性和精度能夠滿足要求.

圖8為2圓面接觸區(qū)域附近的節(jié)點(diǎn)編號.由于模型的對稱性,取Y軸左側(cè),接觸區(qū)域下表面附近的節(jié)點(diǎn)作為研究對象,其節(jié)點(diǎn)編號參見表2.在提取接觸區(qū)附近節(jié)點(diǎn)的求解結(jié)果時,需判斷接觸區(qū)域附近節(jié)點(diǎn)的接觸狀態(tài).可以通過提取接觸區(qū)域附近節(jié)點(diǎn)的節(jié)點(diǎn)反力FY來判斷節(jié)點(diǎn)的接觸狀態(tài),如果節(jié)點(diǎn)反力FY＞0時,則此節(jié)點(diǎn)處于接觸狀態(tài),當(dāng)FY=0時,此節(jié)點(diǎn)為非接觸狀態(tài)[6].節(jié)點(diǎn)的接觸反力之和∑FY應(yīng)為4 000 N.

圖8 2圓面接觸區(qū)域附近的節(jié)點(diǎn)編號

表2 接觸剛度為1.95×106N·mm-2時,兩根鋼絲模型接觸區(qū)附近節(jié)點(diǎn)ANSYS分析結(jié)果

圖9為ANSYS分析得出的位移云圖.

由圖9可以看出,加力位置附近位移最大,2圓面底部加橫向約束、縱向約束的位置附近位移為零,位移值至頂向下逐漸變小,接觸區(qū)域附近相對位移值范圍為0.0808～0.1010 mm.

圖10是加載前后的變形圖,加力位置附近變形最大,變形至頂向下逐漸變小.圖11為等效應(yīng)力云圖,最大等效應(yīng)力值為8 262 MPa,最大等效應(yīng)力位置位于接觸區(qū)域附近.圖12為接觸區(qū)域的接觸應(yīng)力云圖,其接觸應(yīng)力分布曲線近似為平滑的拋物線,這與經(jīng)典理論解得出的接觸應(yīng)力分布曲線相一致.

當(dāng)加大小不同荷載時,2根鋼絲模型的最大接觸應(yīng)力σAmax、接觸半寬bA、接觸點(diǎn)節(jié)點(diǎn)反力之和以及接觸區(qū)域附近節(jié)點(diǎn)的接觸狀況進(jìn)行比較.結(jié)果如表3所示.

表3 加不同大小荷載時,2根鋼絲模型的接觸應(yīng)力、接觸半寬、接觸區(qū)域附近節(jié)點(diǎn)的節(jié)點(diǎn)反力、接觸狀態(tài)

表3表明,隨著所加荷載的增大,由ANSYS有限元方法分析得出的最大接觸應(yīng)力和接觸半寬、接觸區(qū)域附近節(jié)點(diǎn)的節(jié)點(diǎn)反力之和、節(jié)點(diǎn)進(jìn)入接觸狀態(tài)的數(shù)目都有所增大.其中所加荷載增加到4倍,最大接觸應(yīng)力、接觸半徑只近似增加到2倍,即隨著荷載的增加最大接觸應(yīng)力、接觸半寬不是線性增加的.

這與接觸問題經(jīng)典理論解得出的結(jié)論:最大接觸應(yīng)力、接觸半寬與所加載荷的1/2次方成正比(見公式(1)、(2))是相吻合的,這也表明用ANSYS有限元方法求解接觸問題是可靠的.

2 結(jié) 論

本研究先后用經(jīng)典理論解和ANSYS有限元方法對2根拉索平行鋼絲接觸問題進(jìn)行了分析,并得出以下3點(diǎn)結(jié)論.

1)用ANSYS有限元方法分析2根拉索平行鋼絲接觸問題時,在合理選擇接觸剛度的基礎(chǔ)上,得出2根鋼絲間的接觸寬度以及在接觸面上產(chǎn)生的最大接觸應(yīng)力值和經(jīng)典理論得出的結(jié)果是非常接近的.用ANSYS有限元方法分析得出的2根鋼絲模型的位移云圖、加載后變形圖、等效應(yīng)力云圖及接觸區(qū)域的接觸應(yīng)力分布曲線較為直觀、清晰.分析發(fā)現(xiàn),隨著所加荷載的增大,由ANSYS有限元方法分析得出的最大接觸應(yīng)力和接觸半寬、接觸區(qū)域附近節(jié)點(diǎn)的節(jié)點(diǎn)反力之和、節(jié)點(diǎn)進(jìn)入接觸狀態(tài)的數(shù)目都有所增大.

2)將2種分析方法的結(jié)果進(jìn)行對比,驗(yàn)證了用ANSYS有限元方法求解接觸問題是可行的,其求解接觸問題的準(zhǔn)確性和精度能夠滿足要求.

3)2根鋼絲接觸問題作為拉索平行鋼絲間接觸問題的最簡化形式,對其分析結(jié)果的較好把握為接下來解決更加復(fù)雜的拉索平行鋼絲間接觸問題奠定了實(shí)踐基礎(chǔ).同時,其分析結(jié)果也可用來檢驗(yàn)對復(fù)雜問題的解決是否有效.

[1] 謝宜.考慮分層滑移的斜拉索彎曲應(yīng)力的進(jìn)一步分析[D].重慶:重慶交通大學(xué),2009.

[2] 陸萬明,羅學(xué)富.彈性力學(xué)理論基礎(chǔ)[M].北京:清華大學(xué)出版社,1997.

[3] 李紅.斜拉索內(nèi)鋼絲間聯(lián)合工作模式對局部彎曲應(yīng)力的影響[D].重慶:重慶交通大學(xué),2010.

[4] 王新敏.ANSYS工程結(jié)構(gòu)數(shù)值分析[M].北京:人民交通出版社,2007.

[5] 李峰,趙杰.應(yīng)用ANSYS求解赫茲接觸問題[J].機(jī)電產(chǎn)品開發(fā)與創(chuàng)新,2006,19(6):76-78.

[6] 王水林,鄧建輝,葛修潤.改進(jìn)的拉氏乘子法在接觸摩擦問題中的應(yīng)用[J].巖土工程學(xué)報,1998,20(5):64-67.

Finite Element Analysis on Contact Problem between T wo Parallel Stay-cables Steel Wires Based on ANSYS

CHEN Canl,LI Hong2
(1.School of Management,Chongqing Nanfang Translators College of SISU,Chongqing 401120,China;2.School of Civil Engineering&Architecture,Inner Mongolia University of Science and Technology,Baotou 014010,China)

Contact problem of two steel wires is the simplest one of contact problems,so contact problem of two parallel steel wires of stay-cables is analyzed in detail using both the classic theoretical solution and the ANSYS finite element software before contact problem analysis of the parallel wires of the whole cable section by the ANSYS finite element software.Some major parameters,such as contact width,maximum contact stress and distribution of the maximum contact stress between two steel wires are determined.By comparing the results of both methods,we can find that ANSYS is feasible to solve the contact problem and the accuracy and precision of this solution can meet our requirements.

parallel stay-cables steel wire;contact problem;ANSYS;finite element analysis

U443.38

A

1004-5422(2013)01-0071-04

2012-12-19.

陳 璨(1987—),女,碩士,從事結(jié)構(gòu)工程力學(xué)問題研究.