方建豐,陳 玲
(1.天津理工大學,天津市復雜系統(tǒng)控制理論及應用重點實驗室,天津300384;2.天津理工大學機械工程學院,天津300384)
鋼帶纏繞增強塑料復合軟管的有限元分析
方建豐1,2,陳 玲1,2
(1.天津理工大學,天津市復雜系統(tǒng)控制理論及應用重點實驗室,天津300384;2.天津理工大學機械工程學院,天津300384)
提出了一種行之有效的利用有限元法計算螺旋鋼帶彈性模量的方法,在此基礎上,對復合軟管整體的單向拉伸、彎曲、扭轉載荷下的力學性能進行了數(shù)值模擬分析,并與復合軟管的實驗結果進行了對比,結果表明,數(shù)值模擬與實驗結果具有較強的一致性,從而驗證了數(shù)值模擬結果的正確性,為解決復雜結構數(shù)值模擬分析中材料屬性的確定問題提供了一種有效的方法。
塑料復合軟管;螺旋鋼帶;增強;有限元法
復合軟管一般是以加強骨架、合成樹脂黏結劑和塑料制品經(jīng)過一定的成型方法制作而成的新型管道[1]。目前,國外一些發(fā)達國家的復合軟管工業(yè)已經(jīng)發(fā)展得相當成熟,已經(jīng)廣泛應用于化工、石油、城市建設等相關領域[2]。但是由于復合軟管結構的復雜性,導致無法對復雜工作環(huán)境下的復合軟管的力學行為進行精確預測,復合軟管在使用過程中經(jīng)常出現(xiàn)由于管內(nèi)泄露,端部脫扣甚至發(fā)生軟管爆裂等質(zhì)量問題[3-4]。目前的傳統(tǒng)理論與經(jīng)驗方法尚無法解決此類問題。
本文提出了一種利用有限元法計算螺旋鋼帶彈性模量的方法,解決了用傳統(tǒng)的實驗方法難以確定此類螺旋鋼帶彈性模量的問題?;谏鲜龇椒ㄓ嬎愠龅穆菪搸У膹椥阅A?對鋼帶纏繞增強塑料復合軟管的單向拉伸、彎曲、扭轉載荷下的力學性能進行了數(shù)值模擬分析。為了對數(shù)值模擬結果進行驗證,分別對復合軟管進行單向拉伸、彎曲、扭轉實驗。
鋼帶纏繞增強塑料復合軟管是由某公司研制的一種新型的鋼塑復合軟管。這種復合軟管以高強度鋼帶螺旋纏繞形成增強體,以高密度聚乙烯(PE-HD)為基體,并在鋼帶纏繞時施加預緊力,使得外圍的鋼帶與內(nèi)部的PE-HD管緊密地連接在一起。鋼帶纏繞增強塑料復合軟管的內(nèi)部結構共分為5層,包括內(nèi)管、螺旋鋼帶、中層包覆、螺旋鋼帶、外層包覆。復合軟管的結構簡圖如圖1所示。
圖1 鋼帶纏繞增強塑料復合軟管結構簡圖Fig.1 The structure diagram of the steel strip winding reinforced plastic composite pipe
由于復合軟管的螺旋鋼帶的厚度小于1 mm,而且采用雙向螺旋纏繞結構,導致其彈性模量低于常規(guī)結構低碳鋼材料的彈性模量。因此,需要對螺旋鋼帶的彈性模量進行重新確定。由于采用傳統(tǒng)的實驗方法不僅耗時長、成本高,而且其精度難以得到保證,因此本文采用耗時較少、成本較低且精度較高的數(shù)值模擬方法。
為了確定螺旋鋼帶的彈性模量,本文特別設計了一種數(shù)值模擬方法,其實驗原理如下:
由材料力學可知,單向拉伸實驗的變形計算公式為:
式中 A——材料的橫截面積,mm2
P——施加的拉力,N
l——模型長度 ,mm
E——材料的彈性模量,MPa
Δl——材料的伸長量,mm
由式(1)可以看出,彈性模量(E)越大,鋼帶的拉伸變形越小。只要測出一定長度(l)的鋼帶在軸向力 P的作用下的伸長量Δl(彈性變形范圍內(nèi)),即可計算出鋼帶的彈性模量。
由式(1)可導出彈性模量計算公式:
根據(jù)以上原理,本文利用有限元計算軟件ANSYS對螺旋鋼帶的彈性模量進行數(shù)值模擬計算。其有限元模型長度取 l=350 mm,橫截面面積 A=1580 mm2;材料為低碳鋼,其力學性能參數(shù)為 E=2.1×105MPa,泊松比μ=0.3;采用shell181單元進行自由網(wǎng)格劃分,其單元數(shù)為13808,節(jié)點數(shù)36137,該單元是4節(jié)點單元,每個節(jié)點擁有6個自由度,適用于線性、大轉動、大應變等非線性問題的求解。螺旋鋼帶的有限元模型如圖2所示。
圖2 螺旋鋼帶有限元模型Fig.2 The finite element model of the spiral steel strips
由于各層鋼帶之間相互接觸,為了使模擬更加接近實際,提高計算精度,在每一層螺旋鋼帶之間均加入面面接觸單元,其中,目標面為 TARGE170單元,接觸面為CONTA174單元。6層螺旋鋼帶共有5個接觸曲面,劃分13203個接觸單元。圖3為加入接觸單元后的螺旋鋼帶有限元模型。
圖3 加入接觸單元后的螺旋鋼帶有限元模型Fig.3 The finite element model within contact elements
由于各層鋼帶之間尺寸、結構及旋向各有不同,導致每層鋼帶的變形均不相同,在此采用施加摩擦因數(shù)的方法模擬各層間產(chǎn)生的相對滑動。為避免產(chǎn)生剛體位移,在測定彈性模量時,將鋼帶的一端全部自由度約束,另一端施加沿軸線方向的均布載荷 P,該載荷 P以分步加載方式施加。經(jīng)過ANSYS求解,可得出各載荷步的相應變形,載荷 P及變形Δl的具體數(shù)值如表1所示。
表1 各載荷步載荷 P及相應變形ΔlTab.1 Pand displacementΔl of each load step
將表1中各載荷步下的變形Δl代入式(2)即可求得螺旋鋼帶的彈性模量。圖4是螺旋鋼帶在各載荷步下計算出的應力-應變曲線。求得該曲線的斜率,即可以確定螺旋鋼帶的彈性模量 E=0.7×105MPa。
圖4 螺旋鋼帶的應力-應變曲線Fig.4 The stress-strain curve of the spiral steel strip
復合軟管的整體結構主要采用低碳鋼及 PE-HD制成,內(nèi)管、中層包覆、外層包覆均為 PE-HD,其彈性模量 E=108 MPa,泊松比μ=0.4;復合軟管中,螺旋鋼帶的彈性模量已由本文前述計算得出,即 E=0.7×105MPa,泊松比μ=0.3。為減小模擬計算誤差,在建立鋼帶螺旋纏繞增強復合軟管的整體模型時,嚴格按照該復合軟管的設計結構及其外形尺寸精確建立,選擇 solid 186單元作為內(nèi)管、中層包覆、外層包覆模型的單元,選擇shell 181單元作為螺旋鋼帶的單元。網(wǎng)格劃分過程中,內(nèi)管、中層包覆、外層包覆采用映射網(wǎng)格劃分;螺旋鋼帶采用自由網(wǎng)格劃分。該有限元模型將用于復合軟管的單向拉伸、彎曲、扭轉的力學性能數(shù)值模擬。復合軟管整體有限元模型如圖5所示。
圖5 復合軟管整體有限元模型Fig.5 The finite element model of the whole composite pipe
2.2.1 單向拉伸數(shù)值模擬的邊界條件
為了保證單向拉伸計算結果的準確性,在施加邊界條件時,盡量選擇與實際情況相一致的邊界條件。具體做法如下:
(1)考慮軟管各層之間的接觸作用,除上述計算彈性模量時,在各層鋼帶之間加入接觸單元外,在軟管整體結構中的塑料管與鋼帶之間也引入了面面接觸單元,目標面為 TARGE170單元,接觸面為 CONTA174單元。并綜合考慮單元厚度、時間步長、法向接觸剛度以及接觸容差因子等因素對計算結果的影響,選擇出一組最為合理的參數(shù)組合;
(2)為了模擬軟管的固定情況,將模型一端所有節(jié)點的自由度全部約束;
(3)在模型的另一端施加沿軸線方向的均布載荷。按表2的具體載荷數(shù)值,采用分步加載的方式進行加載。
表2 各載荷步的載荷數(shù)值Tab.2 The value of each load steps
2.2.2 單向拉伸模擬結果分析及驗證
圖6 復合軟管各層載荷-變形曲線Fig.6 The load-displacement curve of each layer of the composite pipe
圖6為選取復合軟管內(nèi)部各層中間位置截面全載荷步的變形響應曲線。從圖6可以看出,復合軟管在單向拉伸載荷作用下,其內(nèi)部各層的變形響應趨勢基本一致,符合單向拉伸的變形規(guī)律;由PE-HD制成的內(nèi)管(HDPE_layer1)、中層包覆(HDPE_layer2)和外層包覆(HDPE_layer3)的變形比由多層鋼帶螺旋纏繞制成的螺旋鋼帶(Steel_i1~i6)及(Steel_ii1~ii2)的變形小,這說明復合軟管內(nèi)部各層所使用的材料不同會導致其抗拉剛度有所差異;此外同種材料制成的不同層之間的變形也有所不同,如:同為 PE-HD的內(nèi)管、中層包覆、外層包覆,橫截面積較小的內(nèi)管變形較大;而橫截面積較大的外層包覆變形較小。這說明同種材料各層的抗拉剛度還與其橫截面積的大小有關。
為了驗證復合軟管單向拉伸模擬結果,選用該型號的復合軟管進行單向拉伸實驗。使用拉力試驗機對復合軟管加載,記錄符合軟管端部的變形值。圖7為復合軟管整體結構的端部模擬計算的載荷—變形曲線與實驗結果的對比圖。由圖7可知,復合軟管在軸向拉力載荷作用下,變形比較均勻,計算值與實驗值相比其平均誤差僅為4.81%,計算結果與實驗結果吻合程度較高,數(shù)值模擬結果和實驗結果趨勢基本一致。
圖7 單向拉伸載荷-變形曲線的模擬結果與實驗結果的對比Fig.7 Comparison of simulation and experiment of tension load-displacement results
除單向拉伸載荷的數(shù)值模擬之外,還對模型施加了彎曲及扭轉載荷,其有限元模型的處理方法同上,鑒于篇幅所限,僅給出載荷—變形曲線的模擬結果與實驗結果對比。
從圖8可以看出,彎曲載荷數(shù)值模擬結果與實驗結果基本吻合,除第二載荷步計算誤差較大以外,其余載荷步計算誤差均比較小,數(shù)值模擬結果和實驗結果趨勢基本一致。說明本文所采用的數(shù)值模擬方法能夠較好地預測復合軟管在彎曲載荷下的變形情況。
圖8 彎曲載荷-變形曲線的模擬結果與實驗結果的對比Fig.8 Comparison of simulation and experiment of bending load-displacement results
從圖9可以看出,扭轉載荷數(shù)值模擬結果與實驗結果的趨勢基本一致,其中除第9至第18載荷步的結果與實驗結果誤差相對較大外,其他載荷步的結果與實驗結果誤差均較小。這是由于為了保證數(shù)值模擬計算結果比較容易收斂,在給定材料的本構模型時,只考慮其線彈性階段的本構模型所致。
圖9 扭轉載荷-變形曲線的模擬結果與實驗結果的對比Fig.9 Comparison of simulation and experiment of torsion load-displacement results
(1)通過結果分析以及與實驗結果的對比,證明有限元法可以較好地模擬鋼帶纏繞增強塑料復合軟管在單向拉伸,彎曲及扭轉載荷作用下的變形,并且說明本文所采用的數(shù)值模擬方法能夠較好地預測軟管的變形情況;
(2)采用有限元法模擬鋼帶纏繞增強塑料復合軟管單向拉伸、彎曲、扭轉實驗的精度較高、成本較小,避免了大量重復性的實驗。
[1] 程小平.鋼骨架增強塑料復合管性能研究[D].杭州:浙江大學機械與能源工程學院,2003.
[2] 鄭光明,李秉海,孫曉光.國外復合管的制造和施工技術(一)[J].國外油田工程,2001,17(1):50-54.
[3] 劉輝義,蔡永生,王忠華.國內(nèi)幾種主要給水鋼塑復合管材的對比及選擇[J].焊管,2005,28(3):24-26.
[4] 林秀鋒.鋼絲纏繞增強塑料復合管強度分析與優(yōu)化設計[D].杭州:浙江大學材料與化學工程學院,2006.
FEM Analysis of Steel Strip Winding Reinforced Plastic Composite Pipes
FAN GJianfeng1,2,CHEN Ling1,2
(1.Tianjin Key Laboratory for Control Theory&Applications in Complicated Systems,Tianjin University of Technology,Tianjin 300384,China;2.School of Mechanical Engineering,Tianjin University of Technology,Tianjin 300384,China)
An efficient FEM method for calculating the spiral steel strip′s Young′s modulus was proposed in this paper.Based on the calculation,mechanical properties in three modes,tension,bending and torsion were simulated.The simulation results were in good consistency with the experiment.This method could be applied to determine the material properties of more complex structures.
plastic composite pipe;spiral steel strip;reinforcement;finite element method
TQ320.72+4
B
1001-9278(2011)06-0077-05
2011-01-30
聯(lián)系人,unver@163.com