， ，

1 引 言

薄膜結(jié)構(gòu)的斷裂和穿透破壞過程是一類重要且同時(shí)存在幾何、材料及接觸非線性的復(fù)雜不連續(xù)行為[1]。具體實(shí)例主要有強(qiáng)臺(tái)風(fēng)或雨雪作用下的膜材撕裂、物體撞擊作用下的穿透過程以及充氣氣球的充氣破裂等。

目前，主要通過理論、試驗(yàn)及統(tǒng)計(jì)等方法分析膜材的撕裂強(qiáng)度。吳明兒等[2]針對(duì)梯形形狀的薄膜試件進(jìn)行了拉伸破壞試驗(yàn)，研究各向異性膜材在經(jīng)緯向撕裂時(shí)的破壞強(qiáng)度。陳國(guó)華[3]基于薄膜材料的拉伸束鏈模型，模擬并預(yù)估了其斷裂破壞過程及相應(yīng)的斷裂破壞強(qiáng)度。Aldrich-Smith等[4]對(duì)拉伸應(yīng)力作用下脆性膜/基材料的預(yù)先開裂這一典型失效破壞進(jìn)行了研究。楊班權(quán)等[5]研究了殘余應(yīng)力作用對(duì)脆性薄膜材料撕裂行為的影響，并推導(dǎo)了該膜材裂紋密度隨殘余應(yīng)力變化情況的理論關(guān)系式。至于暴雨、冰雹及物體等撞擊作用下膜材的穿透破壞，目前主要考慮撞擊后的膜材彈性動(dòng)力響應(yīng)和變形。宋維舉[6]利用 LS -DYNA 對(duì)剛性小球撞擊薄膜的動(dòng)力過程進(jìn)行研究，模擬了其在沖擊作用下的大撓度振動(dòng)響應(yīng)。膜材的斷裂和穿透過程具有復(fù)雜不連續(xù)性，已有文獻(xiàn)尚未有通過有限單元法對(duì)其破壞全過程進(jìn)行的數(shù)值研究，本文研究方法將為該類型問題的模擬計(jì)算提供一種新的解決方案。

向量式有限元VFIFE(Vector Form Intrinsic Finite Element)[7-9]通過點(diǎn)值描述和向量力學(xué)分析來獲得結(jié)構(gòu)體系的行為。各質(zhì)點(diǎn)運(yùn)動(dòng)方程在每個(gè)循環(huán)子步中解耦并獨(dú)立逐步求解，解耦后的顯式差分方程不存在集成的大規(guī)模剛度矩陣，也無需進(jìn)行迭代計(jì)算，在結(jié)構(gòu)大變形大轉(zhuǎn)動(dòng)、斷裂及穿透等復(fù)雜不連續(xù)行為分析中具有較好的應(yīng)用。對(duì)于結(jié)構(gòu)斷裂和穿透問題，在統(tǒng)一的主分析程序中，前者引入斷裂失效準(zhǔn)則和質(zhì)點(diǎn)分裂實(shí)現(xiàn)模塊，后者同時(shí)引入碰撞和斷裂機(jī)制模塊；在不增加計(jì)算難度和大量計(jì)算工作的基礎(chǔ)上，通過結(jié)構(gòu)構(gòu)形的步步更新可有效跟蹤獲得其變形、斷裂及穿透破壞全過程。

本文基于向量式有限元薄膜單元理論，研究了薄膜結(jié)構(gòu)的斷裂失效和穿透破壞行為。前者采用失效應(yīng)力判斷準(zhǔn)則，提出了向量式有限元質(zhì)點(diǎn)分裂的斷裂實(shí)現(xiàn)方式，后者則同時(shí)結(jié)合碰撞和斷裂實(shí)現(xiàn)機(jī)制，有效模擬了薄膜結(jié)構(gòu)的變形、斷裂及穿透破壞全過程。編制了薄膜結(jié)構(gòu)的斷裂和穿透求解程序，并通過算例分析驗(yàn)證了理論和程序的可靠性和有效性。

2 向量式膜單元理論

向量式有限元法通過質(zhì)點(diǎn)系的運(yùn)動(dòng)過程來模擬結(jié)構(gòu)的動(dòng)力響應(yīng)，結(jié)構(gòu)質(zhì)量分布于各質(zhì)點(diǎn)，質(zhì)點(diǎn)內(nèi)力則通過質(zhì)點(diǎn)間的單元來計(jì)算并集成。求解過程包括質(zhì)點(diǎn)全位移、節(jié)點(diǎn)純變形和節(jié)點(diǎn)內(nèi)力計(jì)算，詳細(xì)推導(dǎo)參考文獻(xiàn)[10]。

(1) 質(zhì)點(diǎn)全位移

質(zhì)點(diǎn)受力情況如圖1所示，滿足運(yùn)動(dòng)微分方程：

(1)

在每個(gè)時(shí)間步內(nèi)，采用中央差分公式求解質(zhì)點(diǎn)運(yùn)動(dòng)微分方程(1)，得到該時(shí)間步末時(shí)刻質(zhì)點(diǎn)的新位置，進(jìn)而得到該時(shí)間步的質(zhì)點(diǎn)全位移。

(2) 節(jié)點(diǎn)純變形

質(zhì)點(diǎn)全位移相對(duì)應(yīng)的單元節(jié)點(diǎn)全位移包括平移和轉(zhuǎn)動(dòng)(圖2)以及純變形。為避免剛體位移引起的誤差，采用逆向運(yùn)動(dòng)進(jìn)行去除后獲得節(jié)點(diǎn)純變形，從而將質(zhì)點(diǎn)空間運(yùn)動(dòng)轉(zhuǎn)換到單元平面上。剛體平移取為參考節(jié)點(diǎn)的全位移，剛體轉(zhuǎn)動(dòng)則通過單元面外和面內(nèi)的逆向剛體轉(zhuǎn)動(dòng)矩陣來?yè)Q算求解。

(3) 節(jié)點(diǎn)內(nèi)力

通過節(jié)點(diǎn)純變形的單元虛功方程來計(jì)算節(jié)點(diǎn)內(nèi)力。先引入變形坐標(biāo)系和單元形函數(shù)分布，將空間問題轉(zhuǎn)化到平面上，并以變形坐標(biāo)系下單元節(jié)點(diǎn)的內(nèi)力向量為坐標(biāo)軸方向，如圖3所示；再利用坐標(biāo)系轉(zhuǎn)換和面外、面內(nèi)正向轉(zhuǎn)動(dòng)，將其轉(zhuǎn)換到整體坐標(biāo)系下；節(jié)點(diǎn)內(nèi)力反向作用于質(zhì)點(diǎn)并集成即為質(zhì)點(diǎn)內(nèi)力。

圖1 質(zhì)點(diǎn)受力圖

Fig.1 Particle force diagram

圖2 單元的剛體運(yùn)動(dòng)

Fig.2 Rigid body motion of element

圖3 單元節(jié)點(diǎn)純變形和內(nèi)力

Fig.3 Element node pure deformation and internal force

3 斷裂準(zhǔn)則及實(shí)現(xiàn)

薄膜結(jié)構(gòu)的斷裂破壞包括薄膜單元的破裂和破碎情況，本文基于失效應(yīng)力準(zhǔn)則和質(zhì)點(diǎn)分裂方式來實(shí)現(xiàn)薄膜結(jié)構(gòu)的斷裂破壞全過程，以下給出了基本原理和推導(dǎo)過程。

采用von Mises應(yīng)力瞬時(shí)變量達(dá)到失效應(yīng)力限值作為膜材單元的斷裂失效判據(jù)，見式(2)。薄膜單元的向量式有限元基本斷裂過程即可通過該斷裂判據(jù)來完成。

(斷裂失效) (2)

斷裂行為是通過質(zhì)點(diǎn)分裂方式使得相連單元的對(duì)應(yīng)節(jié)點(diǎn)斷開連接，斷裂前后的相關(guān)變量如下。

(1) 單元節(jié)點(diǎn)的von Mises應(yīng)力等狀態(tài)參量是通過該節(jié)點(diǎn)對(duì)應(yīng)單元最近積分點(diǎn)處的狀態(tài)參量獲得，且單元僅在節(jié)點(diǎn)處出現(xiàn)斷裂。節(jié)點(diǎn)Mises應(yīng)力達(dá)到失效應(yīng)力值時(shí)斷開連接。

(2) 質(zhì)點(diǎn)分裂形成的新質(zhì)點(diǎn)屬性根據(jù)斷裂后情況進(jìn)行重新定義。質(zhì)點(diǎn)質(zhì)量m根據(jù)斷裂后新的單元連接進(jìn)行重分配；斷裂后新質(zhì)點(diǎn)的位移x和速度v則取保持?jǐn)嗔亚百|(zhì)點(diǎn)狀態(tài)不變。

以質(zhì)點(diǎn)相連單元中僅一個(gè)單元斷開為例，質(zhì)點(diǎn)分裂示意如圖4所示，其中E1～E4為單元編號(hào)，N1和Ni(i=2,3)分別為分裂前和分裂后的質(zhì)點(diǎn)編號(hào)。質(zhì)點(diǎn)分裂前后各質(zhì)點(diǎn)的狀態(tài)變量滿足

(3)

式中x1和v1為分裂前的質(zhì)點(diǎn)位移和速度，xi和vi(i=2,3)為分裂后的質(zhì)點(diǎn)位移和速度。

(3) 單元的節(jié)點(diǎn)斷開方式包括部分?jǐn)嚅_和所有節(jié)點(diǎn)斷開，前者一般先出現(xiàn)，對(duì)應(yīng)為單元的破裂行為，后者則對(duì)應(yīng)單元的破碎行為。

圖4 質(zhì)點(diǎn)分裂示意圖

Fig.4 Particle splitting schematic diagram

4 穿透行為實(shí)現(xiàn)

膜結(jié)構(gòu)的穿透破壞行為同時(shí)結(jié)合了碰撞和斷裂情況，在每個(gè)循環(huán)分析子步中，首先出現(xiàn)碰撞檢測(cè)并進(jìn)行碰撞響應(yīng)，接著通過斷裂判據(jù)來判斷斷裂的出現(xiàn)情況，并通過質(zhì)點(diǎn)分裂來實(shí)現(xiàn)斷裂模擬，跟蹤獲得膜材變形及穿透破壞的過程。

碰撞檢測(cè)通過質(zhì)點(diǎn)和三角形單元的點(diǎn)-三角形碰撞檢測(cè)方法來處理，碰撞響應(yīng)則基于結(jié)合中央差分與罰函數(shù)的碰撞接觸反力來實(shí)現(xiàn)。時(shí)間子步中，首先通過中央差分運(yùn)動(dòng)方程求解步末的質(zhì)點(diǎn)新位移；接著對(duì)質(zhì)點(diǎn)進(jìn)行碰撞檢測(cè)處理，對(duì)其中發(fā)生實(shí)際碰撞的質(zhì)點(diǎn)新位置通過響應(yīng)處理來更新替換。薄膜結(jié)構(gòu)碰撞接觸的推導(dǎo)過程參考文獻(xiàn)[11]。

采用von Mises應(yīng)力狀態(tài)變量達(dá)到失效應(yīng)力限值作為膜材斷裂失效判據(jù)，而薄膜結(jié)構(gòu)的斷裂行為則通過質(zhì)點(diǎn)分裂方式，將各相連單元對(duì)應(yīng)節(jié)點(diǎn)進(jìn)行斷開操作。單元部分節(jié)點(diǎn)斷開和單元所有節(jié)點(diǎn)斷開對(duì)應(yīng)為破裂和破碎過程。

5 程序?qū)崿F(xiàn)

在進(jìn)行結(jié)構(gòu)分析計(jì)算時(shí)，傳統(tǒng)有限元需事先區(qū)分類型(如斷裂及碰撞等)并采用專用的公式和程序來實(shí)現(xiàn)，在結(jié)構(gòu)的復(fù)雜不連續(xù)領(lǐng)域應(yīng)用不便。向量式有限元對(duì)各類復(fù)雜不連續(xù)力學(xué)行為均可采用相同概念和流程實(shí)現(xiàn)，即用一個(gè)完整而通用的理論進(jìn)行廣義的結(jié)構(gòu)模擬和行為預(yù)測(cè)。對(duì)于結(jié)構(gòu)斷裂和穿透問題，只需在基于統(tǒng)一牛頓運(yùn)動(dòng)方程的主分析程序中，前者引入斷裂失效準(zhǔn)則和質(zhì)點(diǎn)分裂實(shí)現(xiàn)模塊，后者同時(shí)引入碰撞和斷裂機(jī)制模塊；在不增加計(jì)算難度和大量計(jì)算工作的基礎(chǔ)上，通過結(jié)構(gòu)構(gòu)形的步步更新可跟蹤獲得其變形、斷裂及穿透破壞全過程。

本文采用Matlab編制了薄膜結(jié)構(gòu)的斷裂和穿透求解計(jì)算程序，并進(jìn)行了算例分析驗(yàn)證，表明了理論和程序的可靠性和有效性，實(shí)現(xiàn)流程如圖5所示。

6 算例分析

圖7給出了幾個(gè)典型時(shí)刻含裂縫矩形膜片的拉伸變形和膜材撕裂。可以看出，隨著拉伸位移Ux的遞增，在t=0.6 ms時(shí)首先在膜片兩側(cè)出現(xiàn)拉伸變形和應(yīng)力；接著拉伸變形和應(yīng)力逐漸往中間區(qū)域擴(kuò)展，在t=1.2 ms時(shí)膜片中間區(qū)域達(dá)到較大變形和應(yīng)力，并在裂縫端部開始出現(xiàn)應(yīng)力集中效應(yīng)；在t=1.6 ms時(shí),在裂縫端部由于拉伸變形和應(yīng)力繼續(xù)增大而使其Mises應(yīng)力達(dá)到斷裂失效應(yīng)力并開始出現(xiàn)破裂現(xiàn)象；隨后,膜片沿著裂縫方向的破裂范圍逐漸擴(kuò)大，即發(fā)生撕裂現(xiàn)象(t=2.0 ms,t=2.4 ms)；在t=2.8 ms時(shí),矩形膜片出現(xiàn)完全撕裂。

圖5 斷裂和穿透分析流程圖

Fig.5 Fracture and penetration analysis flow chart

圖6 矩形膜片模型

Fig.6 Rectangular membrane model

圖7 矩形膜片拉伸撕裂典型應(yīng)力變形圖

Fig.7 Typical stress and deformation diagrams of tensile tear for rectangular membrane

本算例矩形膜片的斷裂破壞過程與文獻(xiàn)[2]的試驗(yàn)結(jié)果基本一致，驗(yàn)證了本文方法在薄膜結(jié)構(gòu)斷裂行為模擬中的有效性。由于僅采用了簡(jiǎn)單的Mises應(yīng)力斷裂準(zhǔn)則且并未考慮各向異性，與試驗(yàn)結(jié)果的定量分析比較還需進(jìn)一步研究。

圖9給出了幾個(gè)典型時(shí)刻方形氣囊膨脹直至斷裂破壞的變形圖。可以看出，隨著內(nèi)壓p的增大，t=0.02 s 時(shí)，首先在氣囊周邊出現(xiàn)膨脹變形；繼而往中間區(qū)域變形擴(kuò)展，在t=0.05 s時(shí),整個(gè)氣囊均已充滿氣并膨脹起；在t=0.06 s時(shí)，氣囊中心區(qū)域由于von Mises應(yīng)力達(dá)到失效應(yīng)力閾值而開始出現(xiàn)破裂；在t=0.07 s 時(shí)，氣囊破裂往周邊大范圍擴(kuò)展，并開始出現(xiàn)破碎；隨后，膜材的破裂破碎變化區(qū)域明顯，并大范圍擴(kuò)展(t=0.08 s)。

本例采用了簡(jiǎn)單的Mises失效應(yīng)力判斷準(zhǔn)則，并未涉及膜材各向異性及氣囊膨脹破裂后氣體泄漏引起的內(nèi)壓降低問題，因而氣囊膨脹破裂變形與實(shí)際情況有所差異。盡管如此，氣囊出現(xiàn)斷裂的位置及發(fā)展模式均基本符合實(shí)際情況，證明本文方法可有效模擬薄膜結(jié)構(gòu)的斷裂仿真過程。

圖8 氣囊模型

Fig.8 Airbag model

圖9 氣囊膨脹破裂典型變形圖

Fig.9 Typical deformation diagrams of airbag inflating rupture

圖11給出了幾個(gè)典型時(shí)刻馬鞍面薄膜的穿透變形圖?？梢钥闯?，首先薄膜在碰撞位置出現(xiàn)凹陷；隨著凹陷增大，單元節(jié)點(diǎn)應(yīng)力也增大，直至達(dá)到失效應(yīng)力并開始節(jié)點(diǎn)分裂；然后，薄膜在碰撞位置區(qū)域沿著最小應(yīng)力方向(經(jīng)向)呈現(xiàn)撕裂效果；最終，撕裂趨于明顯，剛球逐漸穿透薄膜結(jié)構(gòu)。本算例中馬鞍面薄膜的穿透破壞過程與實(shí)際情況基本一致，驗(yàn)證了本文方法在薄膜結(jié)構(gòu)穿透行為模擬中的有效性。

圖10 剛球-馬鞍面薄膜模型

Fig.10 Rigid ball-saddle membrane model

圖11 剛球-馬鞍面薄膜穿透典型變形圖

Fig.11 Typical deformation diagrams of rigid ball-saddle membrane penetration

7 結(jié) 論

本文將向量式有限元法引入薄膜結(jié)構(gòu)的斷裂和穿透過程分析。提出了膜材的斷裂失效應(yīng)力判據(jù)和質(zhì)點(diǎn)分裂實(shí)現(xiàn)方式，模擬獲得了薄膜變形和斷裂失效全過程；結(jié)合碰撞檢測(cè)和響應(yīng)機(jī)制、斷裂失效和向量式實(shí)現(xiàn)機(jī)制，模擬獲得了薄膜變形和穿透破壞全過程。

編制了薄膜結(jié)構(gòu)斷裂和穿透求解計(jì)算程序，并進(jìn)行了算例驗(yàn)證。結(jié)果表明，該程序可較好完成薄膜結(jié)構(gòu)的大變形大轉(zhuǎn)動(dòng)、斷裂及穿透等復(fù)雜不連續(xù)行為分析，驗(yàn)證了理論和程序的可靠性和有效性。

本文僅對(duì)薄膜結(jié)構(gòu)的斷裂及穿透過程進(jìn)行了基本的模擬分析，要精確實(shí)現(xiàn)膜材斷裂和穿透這一高度非線性問題需要后續(xù)進(jìn)一步的研究，如各向異性材料的引入、斷裂機(jī)制的完善和碰撞摩擦的考慮等。

:

[1] 武 岳,胥傳喜.膜結(jié)構(gòu)設(shè)計(jì)(3)——膜結(jié)構(gòu)的荷載態(tài)分析與結(jié)構(gòu)設(shè)計(jì)[J].工業(yè)建筑,2004,34(8):73-77.(WU Yue,XU Chuan-xi.Structural analysis and design of membrane structures [J].IndustrialConstruction,2004,34(8):73-77.(in Chinese))

[2] 吳明兒,戴 璐,慕 仝.ePTFE膜材力學(xué)性能試驗(yàn)研究[J].建筑材料學(xué)報(bào),2010,13(5):632-635.(WU Ming-er,DAI Lu,MU Tong.Experimental studies on mechanical properties of ePTFE fabric [J].JournalofBuildingMaterials,2010,13(5):632-635.(in Chinese))

[3] 陳國(guó)華.機(jī)織物拉伸斷裂過程模擬及強(qiáng)度預(yù)測(cè)[D].東華大學(xué),2006.(CHEN Guo -hua.The Tensile Breaking Progress Simulation and Strength Prediction of Woven Fabric [D].Donghua University,2006.(in Chinese))

[4] Aldrich-Smith G,Jennett N,Housden J.Adhesion of thin coatings:the VAMAS (TWA 22-2) inter laboratory exercise [J].SurfaceandCoatingsTechno-logy,2005,197(2-3):336-344.

[5] 楊班權(quán),李 朋.殘余應(yīng)力作用下脆性膜/基材料的斷裂行為[J].裝甲兵工程學(xué)院學(xué)報(bào),2011,25(1):88-91.(YANG Ban-quan,LI Peng.Fracture behavior of a brittle thin film / substrate system under residual stress [J].JournalofAcademyofArmoredForceEngineering,2011,25(1):88-91.(in Chinese))

[6] 宋維舉.薄膜結(jié)構(gòu)在沖擊荷載作用下的動(dòng)力響應(yīng)研究[D].重慶大學(xué),2011.(SONG Wei-ju.Study on Dynamic Response of Orthotropic Membranes under Impact Loading [D].Chongqing University,2011.(in Chinese))

[7] Ting E C,Shih C,Wang Y K.Fundamentals of a vector form intrinsic finite element (Part I) basic procedure and a plane frame element [J].JournalofMechanics,2004,20(2):113-122.

[8] Ting E C,Shih C,Wang Y K.Fundamentals of a vector form intrinsic finite element (Part II) plane solid elements [J].JournalofMechanics,2004,20(2):123-132.

[9] Ting E C,Shih C,Wang Y K.Fundamentals of a vector form intrinsic finite element (Part III) convected material frame and examples [J].JournalofMechanics,2004,20(2):113-123.

[10] 趙 陽(yáng),王 震,彭 濤.向量式有限元膜單元及其在膜結(jié)構(gòu)褶皺分析中的應(yīng)用[J].建筑結(jié)構(gòu)學(xué)報(bào),2015,36(1):127-135.(ZHAO Yang,WANG Zhen,PENG Tao.Membrane element based on vector form intrinsic finite element and its application in wrinkling ana-lysis of membrane structures [J].JournalofBuil-dingStructures,2015,36(1):127-135.(in Chinese))

[11] 王 震,趙 陽(yáng).膜材碰撞接觸分析的向量式有限元法[J].計(jì)算力學(xué)學(xué)報(bào),2014,31(3):378-383.(WANG Zhen,ZHAO Yang.Vector Form Intrinsic Finite Element method for collision-contact analysis of membrane material [J].ChineseJournalofComputa-tionalMechanics,2014,31(3):378-383.(in Chinese))