王占山， 呂利葉， 張佳麗， 楊育林， 張芳芳

(1. 燕山大學(xué) 機(jī)械工程學(xué)院， 河北 秦皇島 066004； 2. 燕山大學(xué) 自潤滑關(guān)節(jié)軸承共性技術(shù)航空科技重點(diǎn)實(shí)驗(yàn)室，河北 秦皇島 066004； 3. 全興工業(yè)研發(fā)(廊坊)有限公司， 河北 廊坊 065000)

聚四氟乙烯/芳綸破斜紋織物襯墊拉伸性能數(shù)值仿真

王占山1，2， 呂利葉1，2， 張佳麗3， 楊育林1，2， 張芳芳1

(1. 燕山大學(xué) 機(jī)械工程學(xué)院， 河北 秦皇島 066004； 2. 燕山大學(xué) 自潤滑關(guān)節(jié)軸承共性技術(shù)航空科技重點(diǎn)實(shí)驗(yàn)室，河北 秦皇島 066004； 3. 全興工業(yè)研發(fā)(廊坊)有限公司， 河北 廊坊 065000)

為了解自潤滑關(guān)節(jié)軸承織物襯墊的拉伸性能，對其拉伸性能進(jìn)行了數(shù)值仿真與試驗(yàn)驗(yàn)證。通過進(jìn)行切片試驗(yàn)，利用激光共聚焦顯微鏡觀測纖維截面形態(tài)，對截面形貌輪廓尺寸進(jìn)行測量計(jì)算，得到了纖維束的幾何參數(shù)，利用ANSYS有限元軟件建立了織物襯墊增強(qiáng)相單胞模型；基于復(fù)合材料細(xì)觀力學(xué)方法和周期性邊界條件，研究了織物襯墊拉伸性能，利用ANSYS得出了聚四氟乙烯(PTFE)/芳綸(Kevlar)混織破斜紋襯墊的拉伸彈性常數(shù)；并將數(shù)值仿真分析結(jié)果與拉伸試驗(yàn)結(jié)果對比，結(jié)果相差約20%，在誤差允許范圍內(nèi)。

織物襯墊； 細(xì)觀模型； 拉伸試驗(yàn)； 有限元仿真

自潤滑關(guān)節(jié)軸承的內(nèi)外圈之間黏結(jié)一層自潤滑材料，該自潤滑材料的結(jié)構(gòu)主要為織物襯墊形式。自潤滑關(guān)節(jié)軸承是關(guān)節(jié)軸承的一種，其關(guān)鍵功能件為自潤滑襯墊，自潤滑襯墊的壽命與可靠性對自潤滑關(guān)節(jié)軸承的重要性不言而喻[1-3]。

在研究計(jì)算中，隨著織物結(jié)構(gòu)的改變，襯墊的基礎(chǔ)數(shù)據(jù)相應(yīng)改變。獲取這些數(shù)據(jù)的傳統(tǒng)方式為物理試驗(yàn)，但是該方式耗時(shí)、耗力、成本高[4]。利用有限元法可以獲得材料力學(xué)性能的更多信息，能夠分析材料非線性過程的歷史，便于分析材料摩擦磨損過程。目前，針對自潤滑材料拉伸性能的研究已取得一些成果。Jaehyon Choi等[5]研究了機(jī)織復(fù)合材料單胞內(nèi)裂紋演化機(jī)制，對單胞模型施加了面內(nèi)的單軸拉伸及剪切載荷。沈雪瑾等[6]根據(jù)具體纖維結(jié)構(gòu)建立了襯墊的彈性性能分析模型，但是僅僅適用于斜紋織物。

織物襯墊的細(xì)觀結(jié)構(gòu)決定了其彈性性能和強(qiáng)度，即決定了宏觀性能。為了實(shí)現(xiàn)織物襯墊的性能優(yōu)化，預(yù)測織物的細(xì)觀結(jié)構(gòu)對宏觀性能的影響程度，需要采用細(xì)觀力學(xué)的分析方法[7]。本文建立了接近真實(shí)狀態(tài)的襯墊細(xì)觀模型，預(yù)測了混合材料性能參數(shù)，并獲得了緯向拉伸情況下的應(yīng)力場和應(yīng)變場，為后續(xù)磨損過程有限元模型的建立奠定了理論基礎(chǔ)。

1 織物襯墊細(xì)觀力學(xué)模型

試驗(yàn)所用混織破斜紋織物襯墊為外購，其經(jīng)向纖維束是Kevlar49，緯向纖維束由PTFE與Kevlar49混織而成，固化基體是酚醛縮醛樹脂。

織物襯墊本身的編織結(jié)構(gòu)具有周期性，最小重復(fù)體積單元(單胞，RVE)反映襯墊的整體性能，單胞結(jié)構(gòu)示意圖如圖1所示。

1.1 織物組織截面形貌

1.2 織物襯墊幾何模型

根據(jù)織物襯墊截面的實(shí)際形貌特征，依據(jù)Peirce模型[8]，使用文獻(xiàn)[9]提出的任意壓扁系數(shù)下的橢圓形截面建模方法，并結(jié)合1.1切片試驗(yàn)測量數(shù)據(jù)，得到織物襯墊纖維束相關(guān)參數(shù)，圖4示出織物緯向截面幾何結(jié)構(gòu)示意圖。圖5示出破斜紋襯墊增強(qiáng)相單胞模型。相關(guān)幾何參數(shù)詳見表1。

Tab.1 Section parameters of warp and weft fibers of broken twill liner μm

注：aj、bj、aw、bw分別為經(jīng)向長短半軸長和緯向長短半軸長，lja、ljb分別為經(jīng)向紗線交叉區(qū)域長度和非交叉區(qū)域長度，lwa、lwb分別為緯向紗線交叉區(qū)域長度和非交叉區(qū)域長度，hj、hw分別為經(jīng)緯紗的屈曲波高。

由圖5可以看出，在經(jīng)緯纖維束交叉處，纖維束之間并沒有相互接觸。這是由于建立模型時(shí)，接觸會(huì)產(chǎn)生干涉問題，因而進(jìn)行了簡化，在誤差允許范圍內(nèi)對參數(shù)進(jìn)行適當(dāng)處理，獲得近似模型。

2 織物襯墊拉伸力學(xué)性能

2.1 有限元計(jì)算結(jié)果及分析

纖維束的性質(zhì)取決于單絲的性質(zhì)，不同廠家生產(chǎn)的單絲的性能不盡相同。根據(jù)廠家提供的參數(shù)，運(yùn)用Chamis模型[10]，使用式(1)計(jì)算纖維束的性能參數(shù)。纖維束及基體性能參數(shù)如表2所示。

(1)

式中：Ef1、Gf1分別為纖維方向的彈性模量和剪切模量，GPa；Ef2、Gf2分別為垂直于纖維方向的彈性模量和剪切模量，GPa；μf12為纖維方向泊松比；vf為紗線填充系數(shù)，萃取法測定襯墊中樹脂百分比得到，vf=80%；Em、Gm分別為基體的彈性模量和剪切模量，GPa；μm為基體的泊松比。

表2 纖維束及基體性能參數(shù)

最小重復(fù)單元單胞構(gòu)成宏觀織物襯墊，為使單胞的性能代表襯墊的整體性能，單胞邊界處需滿足位移連續(xù)條件和應(yīng)力連續(xù)條件[11]。根據(jù)文獻(xiàn)[10]提供的周期性邊界條件公式，導(dǎo)出襯墊拉伸預(yù)施加的周期性條件，本文僅針對緯向拉伸進(jìn)行探討，對織物襯墊施加緯向方向的拉伸，即x方向的拉伸，預(yù)施加的周期性邊界條件如式(2)所示。

(2)

式中k為位移系數(shù)。

增強(qiáng)相單胞由纖維束和基體組成，其中纖維束是PTFE纖維和Kevlar纖維，為橫觀各向同性材料?；w是酚醛縮醛樹脂，為各向同性材料。纖維束形狀較為規(guī)則，宜選用六面體單元用體積掃略的方式劃分網(wǎng)格。纖維束單元為Solid45單元，根據(jù)Solid45單元特性，設(shè)定網(wǎng)格單元坐標(biāo)系沿纖維屈曲方向，實(shí)現(xiàn)纖維束的橫觀各向同性?；w的形狀非常不規(guī)則，無法選用六面體單元，宜選用帶中間節(jié)點(diǎn)的四面體單元(Solid92單元)采用自由網(wǎng)格劃分方式?？椢镆r墊增強(qiáng)相單胞網(wǎng)格圖如圖6所示。

根據(jù)式(2)施加邊界條件，建立位移約束方程，進(jìn)行有限元運(yùn)算，x向拉伸邊界條件下經(jīng)緯纖維束應(yīng)力云圖如圖7所示。

增強(qiáng)相單胞中，纖維束主要承受載荷，基體起固定的作用。從纖維束的應(yīng)力云圖中可以看出，經(jīng)緯纖維束交織的位置為最大應(yīng)力處，即最容易發(fā)生拉伸破壞。

利用APDL中ETABLE命令分別求出有限元模型每個(gè)網(wǎng)格單元的平均應(yīng)力(平均應(yīng)變)與單元體積，并賦予相應(yīng)單元表中；將單元平均應(yīng)力表(平均應(yīng)變表)與單元體積表相乘，賦予相應(yīng)數(shù)據(jù)表；并對數(shù)據(jù)表進(jìn)行求和計(jì)算。等效平均應(yīng)力與平均應(yīng)變分別根據(jù)式(3)、(4)計(jì)算得到。

(3)

(4)

依照式(3)、(4)計(jì)算，得出x方向平均應(yīng)力為2.260×107Pa，x方向平均應(yīng)變?yōu)?.979 2%，y方向平均應(yīng)變?yōu)?.597 3%，計(jì)算得出x方向拉伸彈性模量值為2.310GPa，泊松比為0.610。

2.2 拉伸試驗(yàn)

拉伸試驗(yàn)主要參考GB/T 3923.1—1997《紡織品織物拉伸性能：斷裂強(qiáng)力和斷裂伸長率的測定(條樣法)》和GB/T 1447—2005《纖維增強(qiáng)塑料拉伸性能試驗(yàn)方法》進(jìn)行測試。

織物拉伸試驗(yàn)設(shè)備如表3所示。同時(shí)拉伸試驗(yàn)機(jī)應(yīng)具有一定的穩(wěn)定性，載荷波動(dòng)較小，小于+1%；速度相對恒定，當(dāng)速度較小(≤10mm/min) 時(shí)誤差小于等于20%，當(dāng)速度較大(>10mm/min)時(shí)誤差小于等于10%。

表3 織物拉伸試驗(yàn)設(shè)備

準(zhǔn)備試樣：對滿足要求的試樣編號、劃線，測量試樣上工作段內(nèi)的厚度和寬度，任取3處，求平均值。緯向試樣尺寸如表4所示；將夾具和試樣中心線應(yīng)重合；施加破壞載荷5%的載荷消除間隙，保證整個(gè)系統(tǒng)處于正常工作狀態(tài)；緯向拉伸，加載速度為20 mm/min，測量試樣的拉伸載荷和斷裂伸長量。

表4 緯向試樣尺寸

2.3 試驗(yàn)結(jié)果及分析

5組緯向拉伸試驗(yàn)加載拉力與拉伸長度曲線如圖8所示。

從圖8可看出，拉伸初始階段，載荷與位移為線性關(guān)系，此時(shí)襯墊完全處于彈性變形還未開始塑性變形。定義此直線段載荷增量值為△F，△F對應(yīng)的標(biāo)距L0內(nèi)變形增量值記為△L，斷裂時(shí)斷裂伸長量記為△Lb，表5示出緯向試樣拉伸彈性模量測量數(shù)據(jù)。表6示出緯向試樣泊松比測量數(shù)據(jù)。

表5 緯向試樣拉伸彈性模量測量數(shù)據(jù)

表6 緯向試樣泊松比測量數(shù)據(jù)

將試驗(yàn)結(jié)果按照式(5)、(6)進(jìn)行處理，得到彈性模量、泊松比的試驗(yàn)結(jié)果，如表7所示。設(shè)織物緯紗方向?yàn)橹鬏S1方向，經(jīng)紗方向?yàn)橹鬏S2方向，垂直紗線方向?yàn)橹鬏S3方向，根據(jù)表6中的結(jié)果得到平均彈性模量、泊松比。將有限元結(jié)果與實(shí)驗(yàn)結(jié)果對比，如表8所示。

(5)

(6)

式中：Et為拉伸彈性模量，MPa；b、d分別為試樣寬度與厚度，mm；μ為泊松比；ε1為與載荷增量△F對應(yīng)的軸向應(yīng)變，ε1=△L1/L1；ε2為與載荷增量△F對應(yīng)的橫向應(yīng)變，ε2=△L2/L2；L1、L2分別為軸向和橫向的測量標(biāo)距，mm；△L1，△L2分別為與增量△F對應(yīng)的標(biāo)距L1、L2的變形增量，mm。

表7 緯向拉伸試驗(yàn)結(jié)果

表8 有限元模擬結(jié)果與試驗(yàn)結(jié)果對比

由表8可看出，有限元結(jié)果與試驗(yàn)結(jié)果不完全一致，存在可接受的誤差，總結(jié)原因有以下幾點(diǎn)：1)襯墊幾何模型建立時(shí)存在誤差，為防止經(jīng)緯纖維束之間產(chǎn)生干涉，模型中的纖維束不能接觸，與實(shí)際情況存在差距；2)經(jīng)緯纖維束的力學(xué)性能參數(shù)如彈性模量等是通過經(jīng)驗(yàn)公式和將基體平均化推導(dǎo)出來的，與實(shí)際情況存在差異；3)試驗(yàn)過程與真實(shí)情況也存在一定誤差，如試樣測量、試驗(yàn)操作等。綜上所述，在竭力減少誤差存在的可能性時(shí)，有限元結(jié)果與試驗(yàn)結(jié)果存在一定誤差是可以接受的。

3 結(jié) 論

針對織物襯墊采用度驗(yàn)手段觀測襯墊形貌，并依據(jù)Peirce模型，建立織物襯墊的橢圓形截面模型，并建立有限元模型，利用有限元技術(shù)進(jìn)行靜力拉伸仿真；利用試驗(yàn)方法檢測襯墊拉伸性能參數(shù)，實(shí)驗(yàn)與仿真實(shí)驗(yàn)結(jié)果對比，得到破斜紋Kevlar/PTFE混織襯墊緯紗拉伸彈性模量仿真結(jié)果與試驗(yàn)結(jié)果相差約20%；泊松比仿真結(jié)果與試驗(yàn)結(jié)果相差5.20%。本文用混織破斜紋織物襯墊進(jìn)行仿真，驗(yàn)證了織物襯墊細(xì)觀力學(xué)模型的正確性，該模擬方法可推廣到其他類型的織物襯墊，可為以后進(jìn)一步研究襯墊性能提供基礎(chǔ)，為研究織物襯墊摩擦磨損性能奠定了基礎(chǔ)。

FZXB

[1] SAUTTER S, HADEN H, KOTTWITZ B. Spherical plain bearings for on and off road vehicles[J]. Bearing Technology: Analysis, Development and Testing, 1985: 11-26.

[2] 楊育林, 祖大磊, 黃世軍. 自潤滑關(guān)節(jié)軸承現(xiàn)狀及發(fā)展[J]. 軸承, 2009(1):58-61. YANG Yulin, ZU Dalei, HUANG Shijun. Status and development of self-lubricating spherical plain bearings[J]. Bearing, 2009(1): 58-61.

[3] 沈雪瑾, 曹磊, 陳有光, 等. 織物襯墊自潤滑關(guān)節(jié)軸

承的研究現(xiàn)狀與展望[J]. 軸承, 2009(3):57-61. SHEN Xuejin, CAO Lei, CHEN Youguang, et al. Research Status and Prospect of Spherical Plain Bearings with Self-Lubricating Fabric Liner[J]. Bearing, 2009(3):57-61.

[4] 陳繼剛, 薛亞紅, 邱洪亮. 芳綸/聚四氟乙烯纖維織物襯墊力學(xué)建模與性能計(jì)算[J]. 紡織學(xué)報(bào), 2015, 36(1): 82-87. CHEN Jigang, XUE Yahong, QIU Hongliang. Mechanics modeling and performance calculation of aramid/polytetrafluoroethene fabric liner[J]. Journal of Textile Research, 2015, 36(1): 82-87.

[5] JAEHYON C, KUMAR K. Woven fabric composites: part Ⅰ: predictions of homogenized elastic properties and micromechanical damage analysis[J]. International Journal for Numerical Methods in Engineering, 2001, 50: 2285-2298.

[6] SHEN Xuejin, LIU Yunfei, CAO Lei, et al. Numerical simulation of sliding wear for self-lubricating spherical plain bearings[J]. Journal of Materials Research & Technology, 2012, 13(1):8-12.

[7] 燕瑛. 編織復(fù)合材料彈性性能的細(xì)觀力學(xué)模型[J].力學(xué)學(xué)報(bào), 1997(4):429-438. YAN Ying. A micromechanical model for elastic behavior analysis of woven fabric composites[J]. Acta Mechanica Sinica, 1997(4): 429-438.

[8] PEIRCE FT. The geometry of cloth structure[J]. Journal of the Texitle Institute Transactions, 1937, 28(3): 45-96.

[9] 谷大鵬. 織物自潤滑復(fù)合材料組元厚向分布研究及磨損試驗(yàn)表征[D]. 秦皇島:燕山大學(xué), 2013:33-37. GU Dapeng. Research on constituent distribution along with thickness of fabric self-lubricating composite and characterization of wear experiment[D]. Qinhuangdao: Yanshan University, 2013: 33-37.

[10] CHAMIS CC. Mechanics of composites materials- past, present and future[J]. Journal of Composites Technology & Research, 1989, 11(1): 3-14.

[11] 王新峰, 周光明, 周儲(chǔ)偉, 等. 基于周期性邊界條件的機(jī)織復(fù)合材料多尺度分析[J]. 南京航空航天大學(xué)學(xué)報(bào), 2005, 37(6):730-735. WANG Xinfeng, ZHOU Guangming, ZHOU Chuwei, et al. Multi-scale analyses of woven composite based on periodical boundary condition[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2005, 37(6): 730-735.

Numerical simulation of tensile properties of polytetrafluoroethylene/Kevlar broken twill fabric liner

WANG Zhanshan1,2, Lü Liye1,2, ZHANG Jiali3, YANG Yulin1,2, ZHANG Fangfang1

(1.CollegeofMechanicalEngineering,YanshanUniversity,Qinhuangdao,Hebei066004,China; 2.AviationKeyLaboratoryofScienceandTechnologyonGenericTechnologyofSelf-LubricatingSphericalPlainBearing,YanshanUniversity,Qinhuangdao,Hebei066004,China; 3.QuanxingIndustrialDevelopment(Langfang)Co.,Ltd.,Langfang,Hebei065000,China)

In order to understand the tension property of the fabric liner for self-lubricating spherical plain bearings, numerical simulations and experiment verification of the tensile properties were investigated. First of all, the section testing was performed, and laser scanning confocal micro-scopy (LSCM)was used to observe the section morphology of fibers. By measuring and calculating for the profile, the geometry parameters of the fibers were obtained. The reinforcementunit-cell model of the fabric liner was built using ANSYS finite element software. Then, based on the mesomechanics of composite materials and the periodic boundary conditions, the tensile properties of the fabric liner were studied and the tensile elastic constants of the hybrid polytetrafluoroethylene (PTFE)/Kevlar broken twill fabric liner were obtained using ANSYS. At last, comparison of the experimental results and the simulated ones shows that the error is about 20%, which is within the allowed range.

broken twill fabric liner; mesoscopic model; tensile testing; finite element simulation

10.13475/j.fzxb.20150605506

2015-06-28

2015-12-26

王占山(1986—)，男，博士生。主要從事自潤滑關(guān)節(jié)軸承織物的性能研究。楊育林，通信作者，E-mail：yang_yulin123@163.com。

TB 332

A