許珊珊, 張營營, 徐俊豪, 楚時海
(1.中國礦業(yè)大學 江蘇省土木工程環(huán)境災變與結構可靠性重點實驗室, 江蘇 徐州 221116;2.中國礦業(yè)大學 江蘇建筑節(jié)能與建造技術協(xié)同創(chuàng)新中心, 江蘇 徐州 221116;3.江蘇誠意住宅工業(yè)科技發(fā)展有限公司, 江蘇 徐州 221131)
聚氯乙烯(PVC)涂層織物膜材是以聚酯纖維為基布、PVC為涂層的復合材料.由于其價格合理、力學性能良好,被廣泛應用于大跨空間結構中[1-2].聚酯纖維基布由相互垂直的經(jīng)向和緯向纖維編織而成,其細觀編織結構決定了膜材的各向異性特性[3].PVC涂層織物膜材的材料非線性、幾何非線性及摩擦接觸非線性,導致其在拉伸荷載作用下具有復雜的力學響應、非線性的應力-應變關系、各向異性和荷載歷史依賴性等[4-5].
目前,關于PVC涂層織物膜材本構模型的研究主要分為2種:一種是基于組分材料、纖維微觀結構和力學特性細觀結構的分析模型[4-6];另一種是將膜材視為連續(xù)均質(zhì)材料,唯象描述材料宏觀行為的宏觀模型[7-8].由于雙軸拉伸試驗與膜面實際工作狀態(tài)更為接近,有學者根據(jù)雙軸拉伸試驗,引入2個與彈性模量變化相對應的新參數(shù),補充現(xiàn)有平面應力正交異性模型,提出了1種新型非線性各向異性模型[2].隨著膜材拉伸過程的進行,變形表現(xiàn)出明顯的黏彈性,基于材料的黏彈性蠕變特性,Chung等[9]和Sun等[10]根據(jù)各向異性塑性方程和流動法則,定義了蠕變勢函數(shù),在此基礎上,Yu等[11]和Zhang等[12]由偏軸拉伸試驗,得出材料的非線性各向異性蠕變本構方程[11-12].基于聚二氟乙烯涂層織物膜材的偏軸拉伸試驗,易洪雷等[13]研究了其彈性本構關系和強度準則.
本文采用偏軸拉伸試驗研究了PVC涂層織物膜材的各向異性本構關系.同時,結合試驗結果,基于各向異性彈塑性理論,提出了1種三參數(shù)形式的彈塑性屈服函數(shù),并通過分段函數(shù)擬合了PVC涂層織物膜材的應力-塑性應變曲線,以反映其固有的非線性性質(zhì),最終得到能有效預測PVC涂層織物膜材任意角度應力-應變關系的非線性彈塑性本構方程.
材料的變形為彈塑性變形,且變形表現(xiàn)出明顯的各向異性.假設材料變形為正交各向異性,根據(jù)正交各向異性材料的塑性關系理論和流動法,得到以下表達式:
(1)
(2)
由于二次多項式形式的屈服函數(shù)通常能夠與試驗數(shù)據(jù)之間取得良好的吻合[15-16],故假定屈服函數(shù)f為σij的二次函數(shù),則有:
(3)
(4)
(5)
式中:Sp為塑性柔度.
(6)
由此,得到膜材的屈服函數(shù)表達式為:
(7)
通常,在平面應力狀態(tài)下,材料的塑性功增量dWp可由應力和塑性應變增量確定[16]:
(8)
用式(6)中的應力代替應變增量,則塑性功增量表達式為:
(9)
(10)
進一步定義等效應力為:
(11)
由式(10)、(11)得到等效塑性應變增量與等效應力的關系式為:
(12)
由式(11)、(12)可將等效塑性應變增量表示為:
(13)
對于偏軸加載試件,其應力狀態(tài)可根據(jù)材料方向上的應力轉軸公式確定[17]:
(14)
式中:θ為膜材試樣主軸偏離經(jīng)向纖維的角度;σθ為偏軸方向應力.
將式(14)代入到式(10)、(12)中,得到等效應力和等效塑性應變增量的表達式為:
(15)
其中,
H(θ)=
(16)
假設加載過程為比例加載,則等效塑性應變的表達式為:
(17)
偏軸拉伸試驗選取法拉利公司產(chǎn)PVC涂層織物膜材Ferrari 1002T2作為研究對象.由于PVC涂層織物膜材表面平滑,長條形試件在拉伸過程中經(jīng)常在夾持端處斷裂或滑落,導致試驗結果失效,因此本試驗采用啞鈴形試件,以經(jīng)向為基準方向,選取0°、5°、15°、45°、75°、85°和90°這7個偏軸角度試件進行單軸拉伸試驗,拉伸速率為100mm/min,每個方向取5個試件進行檢測.偏軸拉伸試件的應力-應變(σ-ε)曲線如圖1所示.
圖1 偏軸拉伸試件的應力-應變曲線Fig.1 Stress-strain curves of off-axial specimens
由圖1可見:膜材的σ-ε曲線呈明顯的各向異性,且曲線呈明顯階段性的先凸后凹;隨拉伸變形的增大,膜材除了彈性變形外還有明顯的塑性變形.
膜材在初始彈性階段拉伸變形很小,纖維基布和涂層共同受力,膜材剛度較為均勻.為了獲得膜材偏軸方向的彈性模量(E),選取膜材初始直線段進行線性擬合.具體方法為:從初始的5個應力-應變數(shù)據(jù)點開始,逐個增加數(shù)據(jù)點來計算應力和應變之間的相關系數(shù),取相關系數(shù)最大的1組數(shù)據(jù)點進行線性擬合,擬合曲線的斜率即為偏軸拉伸試件的彈性模量.各偏軸方向的彈性模量計算結果列于表1.
表1 偏軸拉伸試件的彈性模量Table 1 Elastic modulus of off-axial specimens
在確定偏軸方向初始彈性模量后,根據(jù)偏軸拉伸試件的應力-應變曲線得到試件偏軸方向上應力與塑性應變的關系,其表達式為:
(18)
式中:εp為試件在偏軸方向上的塑性應變;ε為試件在偏軸方向上的應變;εe為試件在偏軸方向上的彈性應變.
根據(jù)式(18)得到試件在偏軸方向上的應力-塑性應變(σ-εp)曲線,如圖2所示.
圖2 偏軸拉伸試件的應力-塑性應變曲線Fig.2 Stress-plastic strain curves of off-axial specimens
由圖2可見,偏軸拉伸試件的σ-εp曲線呈明顯的非線性,且由于膜材的細觀結構,試件在拉伸過程中σ-εp曲線呈現(xiàn)先凸后凹.針對此現(xiàn)象,在對偏軸拉伸試件σ-εp曲線進行擬合時采用分段函數(shù)形式,即假設前半段的本構關系為指數(shù)函數(shù)形式,后半段的本構關系為冪函數(shù)形式,其表達式為:
(19)
式中:a1、b1、a2、b2為待定系數(shù).
利用式(18)對偏軸拉伸試件的σ-εp曲線進行擬合,擬合效果如圖3所示.
圖3 偏軸拉伸試件的應力-塑性應變擬合曲線Fig.3 Stress-plastic strain fitting curves of off-axial specimens
由圖3可見,分段形式的本構關系表達式能夠很好地擬合PVC涂層織物膜材的應力-塑性應變關系,各擬合曲線的相關系數(shù)R2均在0.99以上,說明前半段為指數(shù)函數(shù)、后半段為冪函數(shù)的表達式能夠用來表達PVC涂層織物膜材的本構關系.
(20) 表2 H(θ)中待定系數(shù)的值Table 2 Values of the tobe determined coefficients in H(θ)
圖曲線聚合及回歸效果Fig.4 Superposition and nonlinear regression
(21)
根據(jù)表2列出的H(θ)中各待定系數(shù)的值,由式(16)擬合得到H(θ)與θ的關系式:
(22)
H(θ)與θ的關系曲線如圖5所示.
圖5 H(θ)與θ的關系曲線Fig.5 Relationship between H(θ) and θ
將偏軸角度(5°和85°)帶入式(21)中,可以得到7種工況下試件的理論應力-塑性應變曲線,如圖6所示.其中帶點曲線表示試驗曲線,淺色實線表示理論曲線.
圖6 偏軸拉伸試件的應力-塑性應變試驗與理論對比曲線Fig.6 Comparison of stress-plastic strain test and theoretical curves of off-axial specimens
采用上述σ-εp的關系式來反推σ-ε的關系式.在應力值不變的情況下,此時應變值ε為塑性應變值與彈性應變值之和,即:
(23)
由式(21)、(23)可將應變表示為:
(24)
圖7 模型預測結果與試驗曲線對比Fig.7 Comparison between model prediction results and test curves
將模型預測的σ-ε關系與試驗結果進行比較,如圖7所示.由圖7可見,上述模型能夠用來預測PVC涂層織物膜材的σ-ε曲線,預測精度存在較小偏差,最大誤差約為15.9%.產(chǎn)生這些偏差的原因是:分析時采用的是名義應力和名義應變,而且膜材在拉伸過程不僅包含彈性變形、不可恢復的塑性變形,還有部分隨著時間可恢復的黏性變形.
(1)PVC涂層織物膜材應力-應變曲線呈明顯的非線性和各向異性,且塑性變形出現(xiàn)在任意加載水平下.
(2)根據(jù)偏軸拉伸曲線的特征,定義了前半段為指數(shù)函數(shù)形式、后半段為冪函數(shù)形式的分階段形式的本構表達式.該表達式能夠預測PVC涂層織物膜材的塑性本構關系.
(3)通過定義屈服函數(shù)、等效應力和等效應變的表達式,推導得到能夠反應材料各向異性的本構關系表達式,并取得了較好的預測結果.