藥曉江,盧華濤,尚捷,王清華,李輝,丁元皓,鄧月,賈建波
中海油田服務(wù)股份有限公司油田技術(shù)研究院,河北 廊坊 065201
隨鉆測量儀器的鉆鋌具有承載電子器件、傳感器或執(zhí)行機(jī)構(gòu)等功能,機(jī)械結(jié)構(gòu)比較復(fù)雜。在井下作業(yè)的過程中,鉆鋌承受拉伸、壓縮、彎曲、扭轉(zhuǎn)等作用力,是隨鉆測量儀器的關(guān)鍵受力部件,其力學(xué)性能的優(yōu)劣關(guān)系到隨鉆測量儀器測量功能的可靠性和井下作業(yè)的安全性。因此,開展隨鉆測量儀器鉆鋌的力學(xué)分析方法研究十分必要[1-4]。
力學(xué)試驗是驗證隨鉆測量儀器鉆鋌力學(xué)性能的重要手段之一[5-8]。對于復(fù)雜的大型鉆鋌結(jié)構(gòu),如果直接研制全尺寸模型進(jìn)行力學(xué)試驗及分析,不僅成本大、周期長,而且對試驗設(shè)備與試驗條件的要求較高。為降低技術(shù)風(fēng)險與成本,縮短儀器研發(fā)周期,通?;诳s比方法,研制縮比模型進(jìn)行力學(xué)試驗。縮比模型試驗?zāi)壳巴ǔV泛用于大型飛艇、運載火箭、水面艦船各種力學(xué)特性分析上??s比模型試驗具有成本低,尺寸小,試驗數(shù)據(jù)測量快速、簡單等特點,可有效降低力學(xué)試驗的技術(shù)風(fēng)險。通過縮比模型試驗,可獲取全尺寸模型力學(xué)性能數(shù)據(jù),以預(yù)測、驗證全尺寸模型的力學(xué)性能,同時還可用來修正全尺寸模型的數(shù)值模型,進(jìn)而驗證全尺寸模型的可行性。為了驗證縮比模型試驗的等價性與可靠性,通常會借助有限元方法對試驗工況進(jìn)行有限元分析[9-15]。下面,筆者針對某型骨架型鉆鋌機(jī)械結(jié)構(gòu),建立了縮比模型,對縮比模型進(jìn)行了拉伸、壓縮、彎曲、扭轉(zhuǎn)4種工況的有限元模擬,并進(jìn)行了4種工況的單載荷力學(xué)試驗驗證縮比模型試驗的可行性,旨在為類似結(jié)構(gòu)的工程設(shè)計及試驗提供科學(xué)依據(jù)和參考。
骨架型鉆鋌結(jié)構(gòu),即將電路骨架與鉆鋌做成一體結(jié)構(gòu)。骨架型鉆鋌結(jié)構(gòu)一般分為上側(cè)段、中間段、下側(cè)段3個部分,通常在中間段開蓋板槽,將需要經(jīng)常裝拆電子器件安裝在蓋板槽中;在上側(cè)段與下側(cè)段上開電路板槽,將電子器件、傳感器均安裝在電路板槽中。圖1為675型骨架型鉆鋌結(jié)構(gòu),全尺寸模型外徑為0.18m,長3.5m。
圖1 675型骨架型鉆鋌結(jié)構(gòu)Fig. 1 Structure of 675 skeleton drill collar
全尺寸模型考察的應(yīng)力應(yīng)變情況為受拉力、壓力、彎矩、扭矩作用下的應(yīng)力及應(yīng)變情況。假設(shè)縮比模型與全尺寸模型直徑比為1∶n,基于相似理論,即在應(yīng)力和應(yīng)變相同的情況下,通過力學(xué)公式,得出全尺寸模型與縮比模型受拉力、壓力、彎矩、扭矩的比例。
1)拉壓工況。拉應(yīng)力(壓應(yīng)力)計算公式為:
(1)
拉應(yīng)力(壓應(yīng)力)與截面積A成反比,即與外徑D的平方成反比,與拉力(壓力)F成正比。為保證應(yīng)力和應(yīng)變相同,則全尺寸模型與縮比模型的拉力(壓力)的比例為n2∶1。
2)彎曲工況。彎曲應(yīng)力計算公式為:
(2)
式中:σ2為彎曲應(yīng)力, MPa;M為彎矩,kN·m;Wz為抗彎截面系數(shù),m3。
彎曲應(yīng)力σ2與抗彎截面系數(shù)Wz成反比,與外徑D的三次方成反比。為保證彎曲應(yīng)力和彎曲角度相同,則全尺寸模型與縮比模型的彎矩比例為n3∶1。
3)扭轉(zhuǎn)工況。扭轉(zhuǎn)應(yīng)力計算公式為:
(3)
式中:τ為扭轉(zhuǎn)應(yīng)力, MPa;T為扭矩,kN·m;Wt為抗扭截面系數(shù), m3。
扭轉(zhuǎn)應(yīng)力與抗扭截面系數(shù)Wt成反比,與外徑D的三次方成反比。為保證應(yīng)力和應(yīng)變相同,則全尺寸模型與縮比模型的扭矩比例為n3∶1。
由于上側(cè)段與下側(cè)段形狀結(jié)構(gòu)相似,均僅取中間段和下側(cè)段作為試驗研究的對象,這2段包含了試件所有的結(jié)構(gòu)要素。由于下側(cè)段包含結(jié)構(gòu)要素比中間段更全面、更典型,同時承受上側(cè)段與中間段所有載荷,故以下側(cè)段為例進(jìn)行試驗研究。
由于全尺寸模型長度與外徑過大,不能滿足試驗機(jī)的工作要求,全尺寸模型不能直接用來做力學(xué)試驗。因此,采用了縮比方法,將全尺寸模型等比縮小二分之一,在縮比模型上施加等比例的力與扭矩,進(jìn)行力學(xué)試驗。縮比模型外形如圖2(a)所示,縮比模型實物如圖2(b)所示。
圖2 縮比模型外形與實物Fig. 2 Appearance and physical object of scaled model
縮比模型材料與全尺寸模型一致,采用15-15HS MAX材料。表1列出了縮比模型的材料屬性。
表1 縮比模型材料屬性
由于全尺寸模型與縮比模型直徑比為2∶1,對于拉壓試驗,由式(1)可知,全尺寸模型與縮比模型的拉力和壓力的比例為4∶1;對于彎曲試驗而言,由式(2)可知,全尺寸模型與縮比模型的彎矩比例為8∶1;對于扭轉(zhuǎn)試驗而言,由式(3)可知, 全尺寸模型與縮比模型的扭矩比例為8∶1。
全尺寸模型所采用的拉壓試驗載荷、彎曲試驗載荷與扭轉(zhuǎn)試驗載荷,按上述比例進(jìn)行縮放,得到縮比試驗的加載載荷,如表2所示。
表2 全尺寸模型和縮比模型的載荷對比
試驗過程中,由于試驗條件的影響,只能使用300kN的試驗機(jī)做拉壓試驗,因而將拉伸試驗的拉力和壓縮試驗的壓力調(diào)整為299kN,為了得到更大的應(yīng)力應(yīng)變,彎曲試驗按照純彎曲試驗加載彎矩。表3列出了調(diào)整后的試驗載荷數(shù)值。
表3 縮比試驗實際加載載荷
有限元模型由骨架型鉆鋌縮比模型和試驗工裝組成。采用ANSYS Workbench中的SCDM模塊對有限元模型進(jìn)行簡化處理,刪除一些對剛度、強(qiáng)度沒有影響的結(jié)構(gòu),將處理后的有限元模型導(dǎo)入ANSYS Workbench,進(jìn)行全四面體網(wǎng)格劃分。
根據(jù)縮比模型應(yīng)變片測量位置,采用2mm×1mm的矩形在幾何上進(jìn)行映射,并提取各矩形上相關(guān)方向的應(yīng)變平均結(jié)果,與試驗結(jié)果進(jìn)行比對。針對每個工況,分為10個載荷步進(jìn)行計算,每個應(yīng)變片位置提取10個計算結(jié)果。
3.2.1 拉伸與壓縮試驗工況
圖3(a)為拉力試驗網(wǎng)格模型,試件網(wǎng)格尺寸為4mm,節(jié)點個數(shù)為811305。拉力施加在右端的連接工裝上,在左端螺紋面的連接工裝端部施加固定約束,如圖3(b)所示。壓力與拉力模型相同,只是加載力的方向相反。
圖3 縮比模型拉力試驗網(wǎng)格模型及荷載與邊界Fig. 3 Scaled model tensile test mesh model and external load and boundary
3.2.2 彎曲試驗工況
圖4(a)為縮比模型彎曲試驗網(wǎng)格模型,試件網(wǎng)格尺寸為4mm,節(jié)點個數(shù)為810526。彎曲工況采用延長梁上遠(yuǎn)程加載拉力方式,兩端分別施加拉力34kN,在右端工裝連接孔上放開繞Z軸轉(zhuǎn)動自由度,在左端工裝連接孔上放開繞Z軸轉(zhuǎn)動與沿X軸平動自由度,載荷與邊界情況如圖4(b)所示。
圖4 縮比模型彎曲試驗網(wǎng)格模型及荷載與邊界Fig. 4 Scaled model bending test mesh model and external load and boundary
3.2.3 扭轉(zhuǎn)試驗工況
圖5(a)為縮比模型扭轉(zhuǎn)試驗網(wǎng)格模型,試件網(wǎng)格尺寸為4mm,節(jié)點個數(shù)為673825。扭矩施加在右端螺紋錐形面的連接工裝上,并使工裝只有沿軸線轉(zhuǎn)動的自由度,在左端螺紋面的連接工裝端部施加固定約束,載荷與邊界情況如圖5(b)所示。
圖5 縮比模型扭轉(zhuǎn)試驗網(wǎng)格模型及荷載與邊界Fig. 5 Scaled model torsion test mesh model and external load and boundary
拉伸試驗、壓縮試驗、彎曲試驗與扭轉(zhuǎn)試驗均在300kN電子萬能試驗機(jī)上開展,均需要設(shè)計專用加載工裝。通過試驗機(jī)加載橫梁向上或向下運行從而施加拉力、壓力、純彎矩以及純扭矩載荷,通過試驗機(jī)和應(yīng)變分析儀、計算機(jī)采集和處理數(shù)據(jù)。
在應(yīng)力較大的地方及典型的幾何形狀處布置應(yīng)變片;較大的面上、較寬溝槽中,以及溝槽的較深位置多布置應(yīng)變片,在一些拐角處,也增加一些應(yīng)變片,以考察扭轉(zhuǎn)和彎曲工況下拐角處的應(yīng)力。
圖6為縮比模型的應(yīng)變片分布情況。圖6(a)為應(yīng)變片總體分布圖。圖6(b)為短槽內(nèi)S7應(yīng)變片實物照片,圖6(c)為槽口端部短槽內(nèi)S8應(yīng)變片實物照片,圖6(d)槽口邊緣S9、S10應(yīng)變片實物照片,圖6(e)為槽內(nèi)與槽口拐角S6、S11應(yīng)變片實物照片。
圖6 縮比模型的應(yīng)變片分布Fig. 6 Strain gauge distribution of the scaled model
縮比模型在拉壓試驗設(shè)備上的連接如圖7(a)所示,圖7(b)為縮比模型試驗現(xiàn)場彎曲試驗的加載情況。彎曲試驗不僅要采集加載數(shù)據(jù)和應(yīng)變數(shù)據(jù),而且要測量彎曲位移,彎曲位移通過數(shù)值千分表測得,數(shù)據(jù)通過應(yīng)變分析儀的通道采集千分表數(shù)值。圖7(c)為縮比模型扭轉(zhuǎn)試驗現(xiàn)場。扭轉(zhuǎn)試驗除了采集加載數(shù)據(jù)和應(yīng)變數(shù)據(jù)外,還需要測定扭轉(zhuǎn)角位移,角位移通過數(shù)值千分表測得,數(shù)據(jù)通過應(yīng)變分析儀的通道采集。
4.1.1 拉伸與壓縮工況
試件的拉壓承載能力比較大,但由于試驗條件所限,只加載到299kN,接近試驗機(jī)最大載荷,此時試驗機(jī)夾頭位移4mm左右,由于試驗系統(tǒng)間隙的影響,加載初始段曲線平緩,當(dāng)間隙消除后,幾乎成線性加載。由于試驗系統(tǒng)的滯后作用,卸載曲線低于加載曲線??s比模型的拉伸試驗加載曲線如圖8(a)所示。壓縮試驗的加載力同于拉伸試驗,為299kN,位移不足5mm。加載曲線為線性??s比模型的壓縮試驗加載曲線如圖8(b)所示。
圖7 拉伸、壓縮、彎曲、扭轉(zhuǎn)試驗圖Fig. 7 Diagram of tensile, compression, bending and torsion tests
圖8 縮比模型拉伸與壓縮試驗加載曲線Fig. 8 Tensile and compression test loading curve of scaled model
4.1.2 彎曲工況
在不同的彎曲方向,抗彎截面系數(shù)不同,取某一槽口向上做純彎曲試驗。加載曲線如圖9(a)所示。根據(jù)千分表的示值計算得到了縮比模型的撓度,可見在彎矩10kN·m的情況下,撓度基本保持在1400μm左右,如圖9(b)所示。
圖9 縮比模型彎曲試驗加載曲線及撓度曲線Fig. 9 Bending test loading curve and deflection curve of scaled model
4.1.3 扭轉(zhuǎn)工況
縮比模型在進(jìn)行扭轉(zhuǎn)加載試驗時,通過對加載力臂施加向上的拉力實現(xiàn),加載到110kN,即18.7kN·m時,加載曲線不再上升,縮比模型進(jìn)入塑性失穩(wěn)階段,加載曲線如圖10(a)所示。
圖10(b)為縮比模型的扭轉(zhuǎn)角曲線。扭轉(zhuǎn)角由位移傳感器測得的位移轉(zhuǎn)換得到,2個測點分別在縮比模型的兩端,相距400mm。因為測點在轉(zhuǎn)動,當(dāng)角度大的時候測量值稍大一點。固定端的轉(zhuǎn)動主要由夾持端的彈性應(yīng)變引起,角度比較小,扭轉(zhuǎn)的驅(qū)動端轉(zhuǎn)角較大,二者的差就是縮比模型在400mm內(nèi)的扭轉(zhuǎn)角。扭轉(zhuǎn)角在縮比模型彈性扭轉(zhuǎn)時為線性,當(dāng)縮比模型進(jìn)入塑性應(yīng)變時,縮比模型扭轉(zhuǎn)剛度降低,轉(zhuǎn)角增大。扭矩10kN· m時,扭轉(zhuǎn)角約2.4°。在試驗最后,縮比模型進(jìn)入塑性失穩(wěn)階段。
圖10 縮比模型扭轉(zhuǎn)加載曲線與加載轉(zhuǎn)角曲線Fig. 10 Torsion loading curve and loading corner curve of scaled model
4.2.1 拉伸與壓縮工況
圖11 縮比模型拉伸試驗結(jié)果與仿真結(jié)果對比Fig.11 Comparison of tensile test results and simulation results of scaled model
因為拉力為單向加載,縮比模型的應(yīng)力狀態(tài)大部分都是單向應(yīng)力狀態(tài),因而只考慮軸向應(yīng)變。圖11為縮比模型拉伸試驗S6~S11結(jié)果和仿真結(jié)果的對比。應(yīng)變都與拉力成正比,大部分值在(500~600)×10-6左右,最大的是S10,達(dá)到841×10-6,對比前述S10的位置可知,S10處于槽口縱向的邊緣,是典型的應(yīng)力集中處,且應(yīng)力集中系數(shù)較大。最小的為S6,槽口中間的臺階上,該處截面較大,且中心對稱,無應(yīng)力集中,因而應(yīng)力較小。應(yīng)力集中點基本上都在槽的邊角以及截面過渡處,槽縱向的中間在拉力作用下不會應(yīng)力集中。
以仿真結(jié)果減去試驗結(jié)果定義絕對誤差,以絕對誤差除以試驗結(jié)果定義相對誤差。試驗和仿真差別的來源,除了試驗誤差、仿真結(jié)果理論誤差外,還包含試驗和仿真所確定的點的位置不同而引起的偶然誤差。拉伸的仿真結(jié)果大于試驗結(jié)果,相對誤差在10%以下,壓縮的試驗結(jié)果幾乎完全反對稱于拉伸的試驗結(jié)果。
4.2.2 彎曲工況
圖12 縮比模型彎曲試驗結(jié)果與仿真結(jié)果對比Fig. 12 Comparison of bending test results and simulation results of scaled model
彎曲試驗為純彎曲加載,軸向拉伸為其主應(yīng)變方向,因而在彎曲試驗中只考慮軸向拉伸應(yīng)變。圖12為S6~S11縮比模型彎曲試驗結(jié)果和仿真結(jié)果的對比。彎曲試驗中軸向應(yīng)變的大小不僅取決于測點與距離中性層的距離,也取決于測點處的幾何情況,即該處的應(yīng)力集中程度。彎曲引起的軸向應(yīng)變隨著加載線性增加,最大正應(yīng)變?yōu)镾10處的3081×10-6,是典型的應(yīng)力集中點。仿真結(jié)果接近試驗結(jié)果,且變化趨勢相同,相對誤差在10%以下。
4.2.3 扭轉(zhuǎn)工況
圖13 縮比模型S6~S8扭轉(zhuǎn)試驗結(jié)果與仿真結(jié)果對比Fig.13 Comparison of torsion test results and simulation results of scaled model S6~S8
扭轉(zhuǎn)試驗不僅要測得對應(yīng)于扭矩的應(yīng)變,而且要作為破壞試驗考察試件的破壞形式。圖13為縮比模型扭轉(zhuǎn)試驗S6~S8的應(yīng)變和仿真結(jié)果的對比。在彈性階段,應(yīng)變與扭矩成線性關(guān)系,10kN·m時S6處應(yīng)變曲線出現(xiàn)轉(zhuǎn)折,應(yīng)變急劇增大,意味著該點開始進(jìn)入塑性狀態(tài),接著S7、S8都進(jìn)入塑性狀態(tài)。塑性變形發(fā)生于外表面,抗扭截面模量最小的部位,從最外層,離扭轉(zhuǎn)中心最遠(yuǎn)的地方開始,槽底部最后塑性變形,最大的是S6,S7和S8 處于同樣的抗扭截面上,應(yīng)變接近。塑性變形發(fā)生于抗扭截面模量最小的區(qū)域,然后向周圍擴(kuò)展,最后導(dǎo)致整體塑性變形。仿真結(jié)果按照彈性計算,不考慮塑性應(yīng)變,仿真結(jié)果成線性,其值接近試驗結(jié)果,相對誤差在10%以下。
1)建立了鉆鋌結(jié)構(gòu)的1∶2縮比模型,基于相似理論,全尺寸模型與縮比模型拉/壓應(yīng)力、應(yīng)變的比例為4∶1,彎曲應(yīng)力、應(yīng)變的比例為8∶1,扭轉(zhuǎn)應(yīng)力、應(yīng)變的比例為8∶1。
2)進(jìn)行了縮比模型拉伸、壓縮、彎曲、扭轉(zhuǎn)4種單載荷工況的力學(xué)性能試驗,結(jié)果證明,采用縮比模型試驗的方法進(jìn)行分析和驗證是可行的。
3)進(jìn)行了拉伸、壓縮、彎曲、扭轉(zhuǎn)4種單載荷工況的有限元分析,試驗結(jié)果與仿真結(jié)果變化趨勢相同,相對誤差在10%以下,表明采用縮比模型進(jìn)行試驗驗證較合理,可為全尺寸模型的力學(xué)分析提供科學(xué)依據(jù)與參考。