韓志騰+秦洪遠(yuǎn)+劉一鳴+黃丹

摘要： 基于非局部近場動力學(xué)（Peridynamics，PD）理論，對含預(yù)置裂紋的混凝土巴西圓盤劈裂破壞問題進(jìn)行建模分析.將結(jié)構(gòu)離散為包含混凝土材料信息的粒子，引入動態(tài)松弛、分級加載和失衡力守恒等粒子系統(tǒng)數(shù)值算法，構(gòu)建可以自然模擬脆性裂紋擴(kuò)展的PD算法體系.對含不同角度單預(yù)置中心裂紋巴西圓盤的裂紋擴(kuò)展過程進(jìn)行數(shù)值模擬，所得結(jié)果與試驗(yàn)結(jié)果吻合較好，驗(yàn)證所提出的模型和算法正確.進(jìn)一步采用該方法模擬雙預(yù)置裂紋巴西圓盤劈裂破壞過程中的裂紋擴(kuò)展、交匯、貫通過程，通過將所得模擬結(jié)果與試驗(yàn)結(jié)果進(jìn)行比較分析，探究該方法處理多裂紋擴(kuò)展問題的可行性.

關(guān)鍵詞： 巴西圓盤； 劈裂破壞； 近場動力學(xué)； 非局部理論； 微極模型； 多裂紋； 脆性材料

中圖分類號： TU37文獻(xiàn)標(biāo)志碼： A

Abstract： The concrete cracked chevron notched Brazilian disks with split are modeled and analyzed based on the nonlocal Peridynamic（PD） theory. The structure is discretized as material points containing concrete material information， and a series of numerical algorithms of particle system including dynamic relaxation， stepbystep loading and nonequilibrium criterion are introduced into the framework of PD for naturally simulating crack growth. The crack propagation process of Brazilian disks with single preexisting central crack under different initial angles is numerically simulated， and the result is in good consistence with the test result， which verifies the proposed model and algorithm. The crack propagation， intersection and coalescence process in doublenotched Brazilian disks are simulated by the method. The feasibility to study multiple crack growth problems using the proposed method is investigated through the comparison and analysis of the simulation results and the test results.

Key words： Brazilian disk； split failure； peridynamics； nonlocal theory； micropolar model； multiple cracks； brittle material

0引言

混凝土是典型的準(zhǔn)脆性拉壓異性材料和主要的工程材料，其裂紋的出現(xiàn)和擴(kuò)展往往由局部或整體承受拉應(yīng)力所引起.目前主要存在2類試驗(yàn)方法分析材料的拉伸特性：其一為直接拉伸試驗(yàn)，由于混凝土材料自身的非均勻性，直接拉伸試驗(yàn)往往存在偏心力，常導(dǎo)致試驗(yàn)結(jié)果不能正確反映材料的拉伸特性；其二為間接拉伸試驗(yàn)，該方法可避免直接拉伸試驗(yàn)中存在的限制，能夠更為準(zhǔn)確地分析混凝土類材料的拉伸特性，其中以巴西圓盤劈裂拉伸試驗(yàn)[1]最受關(guān)注.

為對巴西圓盤的靜態(tài)和動態(tài)破壞進(jìn)行分析，國內(nèi)外學(xué)者進(jìn)行大量的試驗(yàn)研究和數(shù)值模擬，如：HAERI等[2]分別研究含單、雙預(yù)置裂紋的巴西圓盤試樣受對徑壓縮載荷作用下的裂紋擴(kuò)展過程，得到不同角度預(yù)置裂紋作用下巴西圓盤試件的裂紋擴(kuò)展路徑；錢瑩等[3]采用二階流形元法對沖擊載荷作用下巴西圓盤的劈裂拉伸破壞進(jìn)行數(shù)值模擬，所得模擬結(jié)果與試驗(yàn)結(jié)果吻合較好.陳楓等[4]運(yùn)用權(quán)函數(shù)方法對中心裂紋巴西圓盤試樣的對徑加載復(fù)合斷裂進(jìn)行分析，得到較為準(zhǔn)確的應(yīng)力強(qiáng)度因子公式.

基于新興的非局部近場動力學(xué)（Peridynamics， PD）理論[56]分析材料和結(jié)構(gòu)的破壞問題已成為當(dāng)前計(jì)算力學(xué)和相關(guān)領(lǐng)域的熱點(diǎn).十余年來，國內(nèi)外研究人員將PD方法應(yīng)用于混凝土類脆性材料的破壞分析中，如：劉一鳴等[7]基于態(tài)型PD理論對混凝土板的裂紋擴(kuò)展過程進(jìn)行數(shù)值模擬，分析預(yù)置裂紋角度對于裂紋擴(kuò)展路徑的影響；郁楊天等[8]采用改進(jìn)的PMB模型計(jì)算沖擊載荷作用下單邊缺口混凝土梁的破壞全過程；GERSTLE等[9]，HA等[10]和KILIC等[11]分別運(yùn)用PD方法對脆性材料混凝土構(gòu)件的沖擊破壞、失穩(wěn)等問題進(jìn)行分析；筆者所在課題組早期也曾模擬混凝土板的拉壓和沖擊破壞過程等[1214].

在已有研究基礎(chǔ)上，本文基于PD理論，通過引入動態(tài)松弛、分級加載、失衡力守恒等粒子系統(tǒng)數(shù)值算法，構(gòu)建能夠自然模擬結(jié)構(gòu)由彈性變形到脆性斷裂的PD算法體系，并對含不同角度單預(yù)置裂紋巴西圓盤受對徑壓縮載荷作用下的破壞過程進(jìn)行數(shù)值模擬，所得裂紋擴(kuò)展路徑與試驗(yàn)結(jié)果吻合較好，驗(yàn)證提出的模型和算法正確.在此基礎(chǔ)上，進(jìn)一步采用本文方法模擬雙預(yù)置裂紋巴西圓盤受對徑壓縮載荷作用的裂紋擴(kuò)展、交匯、貫通過程，并將模擬結(jié)果與試驗(yàn)結(jié)果進(jìn)行比較分析，驗(yàn)證所提出的模型和算法在分析多裂紋擴(kuò)展問題時(shí)的可引性.

1PD基本框架

1.1PD基本思想和模型

設(shè)物體占據(jù)空間域R，見圖1.對某一時(shí)刻t，域內(nèi)任一物質(zhì)點(diǎn)x與其鄰近一定范圍δ（稱為近場范圍尺寸）內(nèi)其他物質(zhì)點(diǎn)存在相互作用力f（點(diǎn)對力），則根據(jù)牛頓第二定律可得到關(guān)于物質(zhì)點(diǎn)空間積分形式的運(yùn)動方程為

2.2含中心單裂紋巴西圓盤劈裂破壞

單預(yù)置裂紋巴西圓盤幾何參數(shù)和加載條件見圖4.分析對徑壓縮載荷作用下的單初始裂紋巴西圓盤，其材料參數(shù)和模型尺寸與文獻(xiàn)[2]一致，具體為：彈性模量E=15 GPa，泊松比ν=0.21，密度ρ=2 400 kg/m3，抗拉強(qiáng)度ft=3.81 MPa，抗壓強(qiáng)度fc=20 MPa，圓盤直徑為0.1 m，預(yù)置中心裂紋長度2b=0.03 m.物質(zhì)點(diǎn)間距Δx=0.000 5 m，近場范圍δ=3Δx=0.001 5 m，迭代時(shí)間步長Δt=0.1 μs.按照文獻(xiàn)[2]的試驗(yàn)方案，預(yù)置裂紋角度β分別取為0°，15°，30°，45°，60°和90°.

不同角度預(yù)置裂紋的圓盤裂紋擴(kuò)展路徑不同.以β=60°為例，對圓盤的裂紋擴(kuò)展過程進(jìn)行分析，劈裂破壞過程見圖5.當(dāng)t=1.78 ms時(shí)，預(yù)置裂紋兩個(gè)尖端出現(xiàn)損傷累積；當(dāng)t=1.94 ms時(shí)，裂紋開始擴(kuò)展；當(dāng)t=2.06 ms時(shí)，裂紋貫通.

含不同角度預(yù)置裂紋巴西圓盤劈裂破壞過程中裂紋擴(kuò)展路徑的PD模擬結(jié)果與試驗(yàn)結(jié)果對比見圖6.由此可見，本文結(jié)果與試驗(yàn)結(jié)果高度吻合.

2.3雙裂紋巴西圓盤受對徑壓縮破壞

在第2.2節(jié)基礎(chǔ)上進(jìn)一步采用本文所提出的模型和算法分析含雙預(yù)置裂紋巴西圓盤的劈裂破壞過程.材料參數(shù)與第2.2節(jié)一致，模型見圖7.圓盤直徑0.1 m，預(yù)置裂紋長度2a=0.02 m，兩條預(yù)置裂紋中點(diǎn)間距2b=0.03 m.物質(zhì)點(diǎn)離散間距Δx=0.000 5 m，近場范圍δ=3Δx=0.001 5 m，迭代時(shí)間步長Δt=0.1 μs.分別取下方裂紋角度為0°，30°，60°和90°進(jìn)行數(shù)值模擬.含雙裂紋巴西圓盤劈裂破壞的裂紋擴(kuò)展路徑數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果[2]對比見圖8.當(dāng)下方預(yù)置裂紋角度為0°和90°時(shí)，本文模擬結(jié)果與試驗(yàn)結(jié)果高度吻合；當(dāng)下方預(yù)置裂紋角度為30°和60°時(shí)，本文模擬結(jié)果與試驗(yàn)結(jié)果在裂紋向下擴(kuò)展的路徑上稍有差異，但整體基本吻合.由此可以表明本文所提出的模型和算法能夠準(zhǔn)確模擬2條裂紋的起裂、擴(kuò)展、交匯和貫通過程.

3結(jié)論

基于新興的近場動力學(xué)理論，引入動態(tài)松弛、分級加載和失衡力準(zhǔn)則等粒子系統(tǒng)數(shù)值算法，構(gòu)建具有較高定量計(jì)算精度且能夠模擬準(zhǔn)脆性破壞問題的PD算法體系，模擬含不同角度單預(yù)置裂紋巴西圓盤的劈裂破壞過程，所得裂紋擴(kuò)展路徑與試驗(yàn)結(jié)果高度吻合，驗(yàn)證本文的模型和算法.進(jìn)一步采用本文方法模擬含不同角度雙裂紋巴西圓盤的劈裂破壞過程，所得2條裂紋的擴(kuò)展、交匯、貫通過程與試驗(yàn)結(jié)果吻合較好，表明本文方法適用于多裂紋脆性擴(kuò)展問題，為脆性多裂紋擴(kuò)展問題的高效數(shù)值模擬提供新的選擇.參考文獻(xiàn)：

[1]OUINAS D， BOUIADJRA B B， SERIER B， et al. Numerical analysis of Brazilian bioceramic discs under diametrical compression loading[J]. Computational Materials Science， 2009， 45（2）： 443448. DOI： 10.1016/j.commatsci.2008.11.004.

[2]HAERI H， SHAHRIAR K， MARJI M F， et al. Experimental and numerical study of crack propagation and coalescence in precracked rocklike disks[J]. International Journal of Rock Mechanics & Mining Sciences， 2014， 67（4）： 2028. DOI： 10.1016/j.ijrmms.2014.01. 008.

[3]錢瑩， 楊軍， 王兵臣. 巖石巴西圓盤動態(tài)劈裂的流形元法模擬[J]. 爆炸與沖擊， 2009， 29（1）： 2328. DOI： 10.3321/j.issn：10011455.2009.01.005

QIAN Y， YANG J， WANG B C. Numerical simulation on dynamic split of a rock Brazilian disc using the manifold method[J]. Explosion and Shock Waves， 2009， 29（1）： 2328. DOI： 10.3321/j.issn：10011455.2009.01.005

[4]陳楓， 孫宗頎， 徐紀(jì)成. 單軸壓縮下中心裂紋巴西試樣的權(quán)函數(shù)分析[J]. 巖石力學(xué)與工程學(xué)報(bào)， 2000， 19（5）： 599603. DOI： 10.3321/j.issn：10006915.2000.05.011.

CHEN F， SUN Z Q， XU J C. Weight function solution for centercracked Brazilian disk specimen subjected to uniaxial compression[J]. Chinese Journal of Rock Mechanics and Engineering， 2000， 19（5）： 599603. DOI： 10.3321/j.issn：10006915.2000.05.011.

[5]SILLING S A. Reformulation of elasticity theory for discontinuities and longrange forces[J]. Journal of Mechanics and Physics of Solids， 2000， 48（1）： 175209. DOI： 10.1016/S00225096（99）000290.

[6]黃丹， 章青， 喬丕忠， 等. 近場動力學(xué)方法及其應(yīng)用[J]. 力學(xué)進(jìn)展， 2010， 40（4）： 448459.

HUANG D， ZHANG Q， QIAO P Z， et al. A review on peridynamics method and its application [J]. Advance in Mechanics， 2010， 40（4）： 448459.

[7]劉一鳴， 黃丹， 秦洪遠(yuǎn). 混凝土板裂紋擴(kuò)展的態(tài)型近場動力學(xué)模擬[J]. 計(jì)算機(jī)輔助工程， 2016， 25（5）： 5359. DOI： 10.13340/j.cae. 2016.05.011.

LIU Y M， HUANG D， QIN H Y. Statebased peridynamics simulation for crack propagation in concrete slab[J]. Computer Aided Engineering， 2016， 25（5）： 5359. DOI： 10.13340/j.cae. 2016.05.011.

[8]郁楊天， 章青， 顧鑫. 含單邊缺口混凝土梁沖擊破壞的近場動力學(xué)模擬[J]. 工程力學(xué)， 2016， 33（12）： 8085. DOI： 10.6052/j.issn.10004750.2015.05.0396.

YU Y T， ZHANG Q， GU X. Impact failure simulation of a single edged notched concrete beam based on peridynamics[J]. Engineering Mechanics， 2016， 33（12）： 8085. DOI： 10.6052/j.issn.10004750.2015.05.0396.

[9]GERSTLE W， SAU N， SILLING S. Peridynamic modeling of concrete structures[J]. Nuclear Engineering and Design， 2007， 237（1213）： 12501258. DOI： 10.1016/j.nucengdes.2006.10.002.

[10]HA Y D， BOBARU F. Characteristics of dynamic brittle fracture captured with peridynamics[J]. Engineering Fracture Mechanics， 2011， 78（6）： 11561168. DOI： 10.1016/j.engfracmech.2010.11.020.

[11]KILIC B， MADENCI E. Structural stability and failure analysis using peridynamic theory[J]. International Journal of Nonlinear Mechanics， 2009， 44（8）： 845854. DOI： 10.1016/j.ijnonlinmec.2009. 05.007.

[12]沈峰， 章青， 黃丹， 等. 基于近場動力學(xué)理論的混凝土軸拉破壞過程模擬[J]. 計(jì)算力學(xué)學(xué)報(bào)， 2013， 30（S1）： 7983.

SHEN F， ZHANG Q， HUANG D， et al. Damage and failure process of concrete structure under uniaxialtension based on peridynamics modeling[J]. Chinese Journal of Computational Mechanics， 2013， 30（S1）： 7983.

[13]HUANG D， ZHANG Q， QIAO P Z. Damage and progressive failure of concrete structures using nonlocal peridynamic modeling[J]. Science China： Technological Science， 2011， 54（3）： 591596. DOI： 10.1007/s1143101143063.

[14]HUANG D， LU G， QIAO P. An improved peridynamic approach for quasistatic elastic deformation and brittle fracture analysis[J]. International Journal of Mechanical Sciences， 2015， 9495： 111122. DOI：10.1016/j.ijmecsci.2015.02.018.

[15]GERSTLE W H， SAU N， SAKHAVAND N. On peridynamic computational simulation of concrete structures[J]. ACI Special Publication， 2009， 265： 245264.（編輯武曉英）