胡文龍, 羅喜東, 楊語, 駱彬, 程暉

(1.西北工業(yè)大學(xué) 機電學(xué)院, 陜西 西安 710072; 2.中航工業(yè)陜西飛機工業(yè)(集團)有限公司, 陜西 漢中 723099)

碳纖維增強復(fù)合材料機械連接的可靠性高,承載能力強,是主承力結(jié)構(gòu)的重要連接手段[1-2]。孔壁作為材料斷口和連接表面,是連接件上的應(yīng)力集中部位和薄弱環(huán)節(jié)?？妆诒砻嫘蚊苍谖⒂^層面是粗糙的,具有接觸的不均勻性,雖然結(jié)構(gòu)尺寸較小,但在很大程度上影響著機械系統(tǒng)的性能及機器的服役壽命等[3-5]。此外,由于復(fù)合材料各向異性和纖維束導(dǎo)向結(jié)構(gòu),加工形成的孔壁表面幾何形貌具有不確定性,表面幾何形貌特征復(fù)雜且伴隨著大量斷崖式突降區(qū)域和不規(guī)則凸起[6-7]。由于復(fù)合材料連接孔壁微觀結(jié)構(gòu)尺寸屬于微米級且具有不確定性,在連接強度研究中難以使用現(xiàn)有的模型,需要根據(jù)實測微觀數(shù)據(jù)進行精準表征重構(gòu)。

目前的重構(gòu)方法大多是基于逆向建模技術(shù),根據(jù)實體模型,掃描數(shù)據(jù),重新生成數(shù)字模型用于后續(xù)的工程問題研究。Li等[8]利用三維數(shù)據(jù)掃描儀測量物理樣品的點云數(shù)據(jù),借助逆向工程軟件Geomagic Studio實現(xiàn)零件的逆向建模和數(shù)字模型的誤差檢測分析。Wang等[9]利用實驗測量得到的數(shù)據(jù)點云重構(gòu)具有隨機形貌的柔性關(guān)節(jié)界面,建立2個粗糙隨機表面之間的接觸模型。但對于復(fù)合材料孔壁形貌的重構(gòu),上述方法并不適用?，F(xiàn)有的逆向建模重構(gòu)技術(shù)是基于插值函數(shù)對采集點云進行擬合重構(gòu)得到連續(xù)曲面[10]。相較于傳統(tǒng)插值方法,分形插值根據(jù)整體與局部相似的原理,將插值數(shù)據(jù)點的變化特征映射到相鄰點之間的局部區(qū)域。復(fù)合材料孔壁復(fù)雜微觀表面上相鄰2個信息點之間的連線并不是線性或者光滑的,而是存在局部的變化特征[11]。傳統(tǒng)插值重構(gòu)方法難以對其進行準確的擬合表征,從而影響后續(xù)物理模型的準確性。然而這種2個信息點之間更精細一級的特征,用分形插值方法得到的結(jié)果更加符合實際。

分形理論作為研究自然界復(fù)雜形狀和結(jié)構(gòu)的工具,在復(fù)雜結(jié)構(gòu)研究中起著至關(guān)重要的作用[12-13]。隨著分形理論研究的深入,針對分形插值應(yīng)用于粗糙表面的研究,出現(xiàn)了一些實用性的方法。Yang等[14]分別使用了中點位移法(MD)和Weierstrass-Mandelbrot分形函數(shù)法(WM)對一維和二維分形表面進行模擬。該研究突出了2種分形模型的形態(tài)特征、統(tǒng)計特性和適用性,為今后分形模型的發(fā)展與應(yīng)用提供了科學(xué)依據(jù)。翟曉雅等[15]提出了一種基于分形模型的3D打印路徑規(guī)劃算法,保證邊緣細節(jié)的精度條件下最大限度地提高打印效率。趙林等[16]以實測激光切割表面輪廓數(shù)據(jù)為基礎(chǔ),應(yīng)用三維表面分形插值算法,對激光切割表面形貌進行分形插值模擬,分析模擬精度。此外,分形模型在海底地形[17]、地震研究[18]、風(fēng)速預(yù)測[19]等領(lǐng)域均有廣泛應(yīng)用。由于周期性的鉆孔路徑,復(fù)合材料連接孔壁具有周期性的形貌特性,較適合使用自仿射的分形插值方法來研究。然而在現(xiàn)有的分形理論中,給定縱向壓縮比的分形插值曲面在任意小的局部都與整體的結(jié)構(gòu)相似,只能得到具有周期變化形態(tài)的自相似結(jié)構(gòu)[20]。而復(fù)合材料粗糙表面的粗糙度變量是一種區(qū)域化的變量,既有隨機性,也有相關(guān)性。因此需要對現(xiàn)有的自仿射分形插值方法進行改進。

針對復(fù)合材料孔壁復(fù)雜表面形貌難以快速表征重構(gòu)的問題,本文以碳纖維增強復(fù)合材料連接表面為研究對象,基于自仿射分形理論,提出了一種分形插值重構(gòu)復(fù)雜表面的方法。利用連接表面少量測量單元,對復(fù)合材料孔壁表面微觀形貌進行了重構(gòu);同時利用大量測量單元建立了準確的原始實測表面模型,通過對比實測表面模型與分形插值重構(gòu)表面模型特征,討論了采樣間隔對插值重構(gòu)結(jié)果的影響,驗證重構(gòu)表面的幾何有效性;最后開展了有限元仿真驗證重構(gòu)表面的力學(xué)有效性。

1 分形曲面與分形插值

1.1 分形曲面與分形插值

分形是指具有自相似性的幾何圖形,它們在工業(yè)設(shè)計領(lǐng)域和自然界中是普遍存在的,例如車輪輪廓、閃電外形等幾何形狀都具有分形元素。分形曲面在計算機圖形學(xué)、地震學(xué)、巖土學(xué)等許多實際領(lǐng)域有廣泛的應(yīng)用。常見的分形曲面主要有:Koch分形曲面、布朗分形曲面、隨機網(wǎng)格分形曲面等[20]。

分形插值曲面是根據(jù)分形曲面的自相似性原理和迭代函數(shù),將已知數(shù)據(jù)插值成具有自相似結(jié)構(gòu)的曲面。與分形插值曲線類似,分形插值曲面任何一個局部都與整體自相似。

1.2 分形插值數(shù)學(xué)模型

根據(jù)分形曲面的理論,學(xué)者已經(jīng)研究了通過迭代函數(shù)系統(tǒng)(IFS)或遞歸迭代函數(shù)系統(tǒng)(RIFS)來構(gòu)造分形曲面[21]。迭代函數(shù)系統(tǒng)以仿射壓縮變換為框架,根據(jù)幾何對象的整體與局部具有自相似的結(jié)構(gòu),將總體形狀以一定的概率按不同的仿射壓縮變換迭代下去,直至得到滿意的分形圖形。在迭代函數(shù)系統(tǒng)的基礎(chǔ)上,分形插值曲面有多種構(gòu)造方法。根據(jù)復(fù)合材料表面幾何形貌自仿射特征,本文選用矩形區(qū)域分形插值曲面的自仿射分形插值方法。

令I(lǐng)=[a,b],J=[c,d],設(shè)區(qū)域D=I×J={(x,y)|a≤x≤b,c≤y≤d},以Δx,Δy為步長,將D剖分為網(wǎng)格

(1)

給定一組網(wǎng)格點上的數(shù)據(jù)(xi,yj,zi,j)(i=0,1,…,N;j=0,1,…,M),使用常見的二元分形插值函數(shù)f:D→R,并且滿足f(xi,yj)=zi,j(i=0,1,…,N;j=0,1,…,M)。

在區(qū)域(xn-1,ym-1,zn-1,m-1)→(xn,ym,zn,m)中,令x方向的壓縮變換為φn(x)=anx+bn,y方向的壓縮變換為ψm(y)=cmy+dm,并且需要滿足條件

(2)

由此可得

(3)

將結(jié)果帶入壓縮變換中,得

(4)

令z方向的壓縮變換為

(5)

式中,en,m,fn,m,gn,m,sn,m,kn,m均為待求解的參數(shù)且滿足

(6)

令sn,m(n∈{1,2,…,N},m∈{1,2,…,M})為自由參數(shù),稱為縱向比例因子,表示大范圍內(nèi)的數(shù)據(jù)變化特征映射到小范圍區(qū)間的縮小倍數(shù),故需要滿足|sn,m|<1。聯(lián)立(5)～(6)式得到方程組,求解得

(7)

將(7)式代入(5)式中得到z方向的壓縮變換表達式。其中參數(shù)en,m,fn,m,gn,m,kn,m根據(jù)(7)式求得,sn,m為縱向壓縮比,控制分形插值曲面形狀,根據(jù)插值形態(tài)人為給定。

2 微觀表面分形插值重構(gòu)

2.1 微觀表面分形插值

2.1.1 整體框架

本文算法總體框架如圖1所示,主體分為三部分。主程序中讀取實驗測量輸出的原始數(shù)據(jù),并對原始數(shù)據(jù)進行填補,再轉(zhuǎn)換坐標矩陣,將三維數(shù)據(jù)文件轉(zhuǎn)換為對應(yīng)x,y坐標的絕對z值矩陣;再通過對分形插值參數(shù)進行計算,得到三維壓縮比;最后對z值矩陣進行插值,得到插值矩陣,完成分形插值過程。

圖1 分形插值算法總體框架圖

2.1.2 幾何特征提取

實驗使用的測量儀器是由奧地利Alicona公司生產(chǎn)的全自動刀具測量儀,如圖2a)所示,能夠測量零部件的微觀表面,其水平光學(xué)分辨率為0.44 μm,垂直光學(xué)分辨率為10 nm,最小測量高度0.01 μm。實驗所用螺栓為航空級鈦合金螺栓,材料為Ti6Al4V,規(guī)格M6。T700級復(fù)合材料層合板連接孔名義直徑5.98 mm,鋪層順序[-45/90/45/0]4s,纖維體積分數(shù)為60%。實驗前將層合板沿著孔的中軸線進行橫向剖切,便于后期的測量觀察。

圖2 微觀幾何形貌測量儀器及測量對象

首先對螺栓桿桿壁進行測量,測量區(qū)域如圖2b)所示。由于金屬的各向同性,鈦合金螺栓只需選取桿壁的一個區(qū)域進行測量,命名為Ti-region。對復(fù)合材料層合板的孔壁進行測量,由于層合板的疊層結(jié)構(gòu),需要在孔的入口處至中間層選取4個具有代表性區(qū)域進行多次測量,分別命名為CFRP-region-1,CFRP-region-2,CFRP-region-3,CFRP-region-mid。具體測量區(qū)域如圖2c)所示。測量區(qū)域為284 μm×216 μm的類矩形區(qū)域,每個測量區(qū)域選取采樣間隔為4個數(shù)據(jù)點,共7 815個測量點。

2.1.3 丟失數(shù)據(jù)填補

主程序中首先讀取原始測量文件,將測量數(shù)據(jù)變換為矩陣形式。由于在進行原始測量數(shù)據(jù)采集時,所使用的光學(xué)顯微鏡容易受到樣本反光率、環(huán)境光和景深等因素的影響,導(dǎo)致表面形貌的個別數(shù)據(jù)丟失。由于丟失數(shù)據(jù)點的位置隨機,需要對整體數(shù)據(jù)進行遍歷,尋找丟失數(shù)據(jù)點位置并對其進行填補。當丟失點位于中心位置時,如圖3a)所示,利用周圍8個點的z值數(shù)據(jù)取平均進行缺失點的填補,即缺失點z值為

圖3 不同位置缺失點的填補方法

(8)

式中,z(i,j)表示第i行j列的數(shù)據(jù)點數(shù)值。同理,當丟失點位于邊角位置時,如圖3b)所示,對丟失的z值進行填補。這種方法符合粗糙表面變化規(guī)律,較好地解決光學(xué)測量丟失少量數(shù)據(jù)點的問題。

將矩陣中三維測量數(shù)據(jù)轉(zhuǎn)換為對應(yīng)x,y坐標的z值。由于光學(xué)顯微鏡測量得到的數(shù)據(jù)是基于全局坐標系,無法直接進行使用,需要得到相對坐標后對其進行全局坐標系與局部坐標系轉(zhuǎn)換。

2.1.4 改進的自仿射分形插值重構(gòu)

在全局統(tǒng)一縱向壓縮比的分形插值曲面中,在任意小的局部都與整體的結(jié)構(gòu)相似,只能得到具有周期變化形態(tài)的自相似結(jié)構(gòu)。復(fù)合材料微觀表面的粗糙度變量是一種區(qū)域化的變量,即既有隨機性,也有相關(guān)性。因此需要使用改進的自仿射分形插值方法,即讓每個局部子區(qū)域總生成自己的縱向壓縮比。

每個子區(qū)間上的壓縮比可以通過該方法得到:首先根據(jù)觀測值數(shù)據(jù),運用最小二乘法,得到一次擬合面方程,該擬合面反映區(qū)域觀測值的平均平面;再利用此方程,求出每個插值點位置上的預(yù)測值;用各點上的觀測值減去相應(yīng)的預(yù)測值,得到偏差值,偏差值反映觀測值偏離平均平面的程度。根據(jù)變差函數(shù)原理,數(shù)據(jù)具有相關(guān)性,某個觀測點的變異特征能夠反映附近較小領(lǐng)域內(nèi)的變化特征。因此以觀測點偏差值作為縱向壓縮比的依據(jù)是合理的。

每個子區(qū)間上的一次擬合面方程表達式為

z=b0+b1x+b2y

(9)

式中,b0,b1,b2為待求系數(shù)。根據(jù)最小二乘法,使得觀測值與預(yù)測值的離差平方和Q最小,即

(10)

分別求Q對b0,b1,b2的偏導(dǎo)數(shù),并令其為0。求解整理可得

(11)

(12)

(13)

取每個子區(qū)間上壓縮比sn,m為

(14)

按照改進的自仿射分形插值方法,對網(wǎng)格點數(shù)據(jù)進行分形插值重構(gòu)。對于給定一組網(wǎng)格點的數(shù)據(jù){xi,yj,zi,j}(i=0,1,…,N,j=0,1,…,M),網(wǎng)格剖分為(1)式所示。x方向和y方向的插值公式由(4)式求出;z方向的插值公式由(5)式求出,每個子區(qū)間的縱向壓縮比sn,m由觀測值結(jié)合(9)～(14)式求出。最后對每組網(wǎng)格數(shù)據(jù)進行插值重構(gòu),得到插值重構(gòu)數(shù)據(jù)。根據(jù)x,y,z的值,重構(gòu)出表面輪廓。

2.2 微觀表面分形插值重構(gòu)結(jié)果

分別對螺栓桿表面區(qū)域與復(fù)合材料孔壁表面4個特征區(qū)域進行分形插值,初始采樣間隔均為4個數(shù)據(jù)點,得到多組分形插值數(shù)據(jù)。其中CFRP-region-mid的表面輪廓重構(gòu)結(jié)果如圖4所示,圖4a)為實驗采集的原始表面輪廓,圖4b)為分形插值重構(gòu)表面。分形插值表面與實測表面的幾何結(jié)構(gòu)相似且重要幾何特征基本相同,如圖4中的紅圈所示。分形插值重構(gòu)表面更平滑,相較于原始實測表面更加的簡化。

圖4 表面輪廓分形插值結(jié)果圖

通過上述分析發(fā)現(xiàn),分形插值后的粗糙表面幾何形貌和原始實測粗糙表面幾何形貌基本相似,但在某些細節(jié)特征上略有不同。圖5左側(cè)為原始實測表面,中間為采樣表面,右側(cè)為插值重構(gòu)表面。分析發(fā)現(xiàn)局部細微差異主要為兩部分,首先如圖5a)所示,當采樣間隔設(shè)置過大時,部分包含在采樣間隔內(nèi)的細節(jié)特征數(shù)據(jù)點會被篩選出去,出現(xiàn)細節(jié)特征丟失。其次如圖5b)所示,當采樣間隔過大時,采樣點與采樣點間的距離過大,會導(dǎo)致插值點所依據(jù)的采樣點超出鄰域范圍,不能很好地表征插值點的原始形貌,導(dǎo)致細節(jié)特征畸變。因此,選取適當?shù)牟蓸娱g隔對于分形插值重構(gòu)的結(jié)果至關(guān)重要,將在3.1.2節(jié)進行詳細說明。

圖5 分形插值表面微觀細節(jié)特征

3 分形插值重構(gòu)表面有效性

3.1 幾何有效性驗證

3.1.1 三維表面輪廓表征與誤差分析

三維表面輪廓表征參數(shù)主要分為幾何參數(shù)和性能參數(shù)[22]。由于本研究主要是對復(fù)合材料連接孔壁進行表征,幾何參數(shù)選用常用的三維輪廓算術(shù)平均偏差Sa;性能參數(shù)選用表面接觸面積比Sdr。其中,三維輪廓算術(shù)平均偏差Sa表示被測表面上所有像素點到基準面上對應(yīng)點的直線距離的算數(shù)平均值。

(15)

式中:l(i,j)表示像素點到基準面的直線距離;M和N分別為基準面的邊長。表面接觸面積比Sdr描述采樣面積內(nèi)表面在水平面上的投影面積和表面積的相差程度。

(16)

式中:A表示總展開面積;Δx,Δy分別為x和y方向相鄰像素點的間隔距離。表征結(jié)果如圖6所示。

圖6 實測表面與重構(gòu)表面三維輪廓參數(shù)表征

對比原始實測表面與插值重構(gòu)表面的三維表面輪廓參數(shù)表征結(jié)果發(fā)現(xiàn),復(fù)合材料孔壁4個代表性區(qū)域的三維輪廓算術(shù)平均偏差Sa的誤差均在2%以內(nèi),接觸面積比Sdr的誤差均在10%以內(nèi)。通過對重構(gòu)模型與實測模型表面幾何參數(shù)分析,認為分形插值模型具有較好的幾何精確性。復(fù)合材料表面的粗糙度起伏較大,并且在分形插值過程中出現(xiàn)細節(jié)特征丟失與細節(jié)特征畸變現(xiàn)象,從而導(dǎo)致粗糙表面總展開面積A產(chǎn)生較大波動,帶來接觸面積比Sdr的誤差。通過上述分析,結(jié)果符合精度要求。

3.1.2 采樣間隔對分形插值的影響

由上文發(fā)現(xiàn)采樣間隔對分形插值結(jié)果影響較大。針對此問題進行研究,又分別進行了相同區(qū)域,采樣間隔為2～8的原始采樣表面測量實驗,使用上文相同的方法,得到采樣間隔分別為2～8的原始采樣表面形貌與分形插值結(jié)果并分析。在某些微觀細節(jié)特征上,原始實測表面與插值重構(gòu)表面有細微差別,并且隨著采樣間隔的增大,相似性降低,會出現(xiàn)原始形貌所不具有的表面特征。使用三維輪廓算術(shù)平均偏差Sa進行表征,如圖7所示。

圖7 采樣間隔對插值重構(gòu)結(jié)果影響

隨著采樣間隔的增大,原始實測表面與插值重構(gòu)表面三維輪廓算術(shù)平均偏差之間的誤差增大。當采樣間隔超過5個數(shù)據(jù)點時,誤差量超過預(yù)估范圍,分形插值結(jié)果與原始表面產(chǎn)生較大差異,精度無法保證。當采樣間隔小于等于5個數(shù)據(jù)點時,誤差值均小于5%,符合精度要求。采樣間隔越小,原始實測表面與分形插值表面的差異越小,難以顯示出等效表征的效果。因此,取采樣間隔為3～5個數(shù)據(jù)點,既可以保證分形插值精度,也可以對原始實測表面進行有效的等效表征。

綜上,通過對比重構(gòu)結(jié)果發(fā)現(xiàn),選取合適的采樣間隔時,分形插值后三維模型的粗糙表面幾何形貌和原始三維模型的粗糙表面幾何形貌基本相似,分形插值重構(gòu)結(jié)果具有幾何有效性。

3.2 力學(xué)有效性驗證

3.2.1 有限元模型

利用有限元軟件進行力學(xué)擠壓仿真模擬,根據(jù)擠壓力學(xué)性能所關(guān)注的物理量,對比原始表面模型與重構(gòu)表面模型的應(yīng)力應(yīng)變以及相互作用面積,驗證分形插值重構(gòu)結(jié)果的力學(xué)有效性,從而驗證分形插值重構(gòu)方法的正確性。

將逆向建模得到的三維實體模型導(dǎo)入有限元軟件中,分別對復(fù)合材料板與鈦合金板進行材料屬性設(shè)置、接觸設(shè)置和邊界條件設(shè)置,給予擠壓載荷,模擬復(fù)合材料裝配連接擠壓過程。最終得到的有限元模型如圖8所示。

圖8 力學(xué)有效性驗證模型

3.2.2 結(jié)果分析

將每組原始表面模型與重構(gòu)表面模型進行有限元仿真,圖9為CFRP-region-mid區(qū)域的仿真結(jié)果。

圖9 力學(xué)有效性驗證仿真結(jié)果

根據(jù)有限元模擬結(jié)果圖可以發(fā)現(xiàn)原始表面輪廓擠壓模型和插值重構(gòu)表面輪廓擠壓模型的接觸區(qū)域大致相同,部分細節(jié)特征有所不同。對每組模型最大應(yīng)力值與最大變形量進行分析,如表1～2所示。結(jié)果發(fā)現(xiàn),4組模型的最大應(yīng)力值和最大變形量差距均在10%以內(nèi)。

表1 各區(qū)域原始模型與重構(gòu)模型最大應(yīng)力值及誤差

表2 各區(qū)域原始模型與重構(gòu)模型最大變形量及誤差

利用二值化程序分別對4組有限元變形云圖的變形區(qū)進行處理計算,得出4組變形區(qū)域與整個面積的百分比,如表3所示。通過分析發(fā)現(xiàn),CFRP-region-1的原始表面輪廓與插值表面輪廓的變形區(qū)域占比相差較大,達到9.5%,4組的差距均在10%以內(nèi),整體趨勢與圖6b)中接觸面積比誤差結(jié)果一致。分析重構(gòu)模型與實測模型力學(xué)性能,認為分形插值模型具有較好的力學(xué)預(yù)測精確性。分析原因發(fā)現(xiàn),由于CFRP-region-1區(qū)域的接觸表面粗糙度較大,加劇細節(jié)特征丟失與細節(jié)特征畸變現(xiàn)象。此外,在進行曲面逆向重構(gòu)擬合時,急劇起伏的曲面擬合誤差較大,擬合的效果不太理想,進而導(dǎo)致原始輪廓模型與插值輪廓模型的誤差。

表3 各區(qū)域原始模型與重構(gòu)模型變形區(qū)域占比及誤差

通過對比重構(gòu)前后力學(xué)擠壓模擬結(jié)果發(fā)現(xiàn),分形插值重構(gòu)預(yù)測模型的力學(xué)性能誤差均在10%以內(nèi),分形插值重構(gòu)結(jié)果具有力學(xué)有效性。

4 結(jié) 論

本文主要提出了一種利用分形插值重構(gòu)研究復(fù)合材料孔壁微觀表面等效表征的方法,并使用重構(gòu)誤差分析和有限元仿真驗證其有效性。分析得出以下結(jié)論:

1) 分形插值過程中,由于細節(jié)特征丟失與細節(jié)特征畸變,原始表面與插值表面有細微差別,并隨采樣間隔增大,相似性降低。當采樣間隔為3～5時,既可以保證分形插值誤差,又可以對原始實測表面進行有效的等效表征。

2) 通過重構(gòu)的曲面模型,對比原始表面與插值表面幾何形貌相似性,三維輪廓算術(shù)平均偏差的誤差均在2%以內(nèi),接觸面積比的誤差均在10%以內(nèi),證明了分形插值重構(gòu)結(jié)果具有幾何有效性。

3) 通過有限元仿真,在擠壓載荷下,實測原始表面與插值重構(gòu)表面的最大應(yīng)力、最大變形量以及接觸變形區(qū)域面積預(yù)測誤差均在10%以內(nèi),證明了分形插值重構(gòu)結(jié)果具有力學(xué)有效性。