邵冬冬章 青夏曉舟

(河海大學(xué)力學(xué)與材料學(xué)院,江蘇南京 210098)

混凝土等顆粒型材料的偶極效應(yīng)分析

邵冬冬,章 青,夏曉舟

(河海大學(xué)力學(xué)與材料學(xué)院,江蘇南京 210098)

考慮物質(zhì)點的偶極效應(yīng),基于Cosserat連續(xù)體介質(zhì)力學(xué)模型,把球形微顆粒推廣到方形微顆粒,推導(dǎo)出方形微顆粒對應(yīng)的偶應(yīng)力和曲率應(yīng)變之間的線彈性關(guān)系,并提出了一種能考慮微顆粒傾覆破壞的廣義屈服準則。編寫了基于Cosserat連續(xù)體介質(zhì)力學(xué)模型的有限元分析程序,進行了混凝土試樣單軸壓縮試驗和土體邊坡穩(wěn)定性分析的數(shù)值模擬。研究結(jié)果表明,物質(zhì)點的偶極效應(yīng)確實存在,且是材料發(fā)生應(yīng)變局部化的直接原因。

偶極效應(yīng);方形微顆粒;應(yīng)變局部化;數(shù)值模擬

傳統(tǒng)連續(xù)介質(zhì)力學(xué)廣泛應(yīng)用于當(dāng)今各個領(lǐng)域,但在微細觀領(lǐng)域,由于尺度效應(yīng)廣泛存在,其適用性受到制約。為了正確而有效地解釋尺度效應(yīng),偶應(yīng)力逐漸成為力學(xué)研究的熱點。早在19世紀法國Cosserat兄弟[1]就提出了偶應(yīng)力理論,即每個顆粒形變時不僅產(chǎn)生位移,還伴隨著轉(zhuǎn)動。在基于Cosserat連續(xù)體理論的有限元方法中,除引入線位移自由度外,還另加入獨立于線位移的相對轉(zhuǎn)動自由度?；贑osserat連續(xù)體理論進行的數(shù)值模擬計算,能有效模擬變形局部化效果并消除了網(wǎng)格依賴性[2]。該理論在分析特定材料[3-5]和特定物理現(xiàn)象[6]時,表現(xiàn)出很高的求解精度和良好的模擬效果。

混凝土、巖體等高度非均質(zhì)材料形成的顆粒夾雜相形狀大小各異,只用球形或圓形來模擬混凝土顆粒等過于單一。本文在Cosserat連續(xù)體介質(zhì)力學(xué)模型的框架下,把球形微顆粒推廣到方形微顆粒,建立考慮方形微顆粒偶極效應(yīng)的偶應(yīng)力與曲率應(yīng)變之間的線彈性關(guān)系,在此基礎(chǔ)上,提出一種能考慮微顆粒傾覆破壞的廣義屈服準則。

1 考慮物質(zhì)點偶極效應(yīng)的基本方程

物質(zhì)點被賦予內(nèi)稟尺度,其廣義應(yīng)力狀態(tài)如圖1所示,圖中σij為應(yīng)力張量分量,μij為偶應(yīng)力張量分量,因此平衡方程還包括旋轉(zhuǎn)平衡,即

式中:bj——體積力分量;Ψi——體力偶分量;eijk——置換符。

與經(jīng)典連續(xù)體相比,剪應(yīng)力互等定理已不再成立[7],由于與曲率能量共軛的偶應(yīng)力的出現(xiàn),角動量的守恒將產(chǎn)生不對稱的應(yīng)力張量。

幾何協(xié)調(diào)方程:

式中:εij——應(yīng)變張量分量;ωk——旋轉(zhuǎn)向量分量;χij——曲率張量分量;ωj,i——對旋轉(zhuǎn)向量求偏導(dǎo)。當(dāng)不考慮旋轉(zhuǎn)向量的影響時,協(xié)調(diào)方程退化為連續(xù)介質(zhì)力學(xué)的幾何方程。應(yīng)變張量ε不僅與位移有關(guān),同時與旋轉(zhuǎn)向量也存在關(guān)系。

圖1 細觀物質(zhì)點的廣義應(yīng)力狀態(tài)Fig.1 Generalized stress state of meso material point

本構(gòu)方程:

式中:D*——Cosserat物質(zhì)點的彈性矩陣。寫成分量形式為

式中:Dijkl——應(yīng)力與應(yīng)變之間的彈性矩陣分量形式;Hijkl——偶應(yīng)力與曲率(扭率)應(yīng)變之間的彈性矩陣分量形式。

2 物質(zhì)點的彈性矩陣

對于球形物質(zhì)點,給定其內(nèi)稟特征尺度,即直徑ld,則偶應(yīng)力μ與曲率(扭率)應(yīng)變χ之間的彈性矩陣H為[8]

式中:G——細觀物質(zhì)點的剪切模量。

對于方形物質(zhì)點,應(yīng)在其局部坐標(biāo)系下建立本構(gòu)方程。在局部坐標(biāo)系下,對于三維問題,物質(zhì)點的廣義應(yīng)力分量為

廣義應(yīng)變分量為

則局部坐標(biāo)系下的本構(gòu)方程可寫為

式中:Gc——Cosserat剪切模量;ν——泊松比。

對于方形物質(zhì)點,給定其內(nèi)稟尺度,即長la、寬lb和高lc,則偶應(yīng)力ˉ與曲率(扭率)應(yīng)變之間的彈性矩陣H為

式中:E——細觀物質(zhì)點的彈性模量;β1、β2、β3——方形截面桿的扭轉(zhuǎn)形狀參數(shù),根據(jù)該截面長寬(長高或?qū)捀?之比來定,具體數(shù)值可由彈性力學(xué)[9]中方形扭轉(zhuǎn)桿的解答給出。

3 考慮物質(zhì)點偶極效應(yīng)的廣義屈服準則

可直接采用截斷型摩爾-庫倫屈服準則[10-11]描述細觀物質(zhì)點拉裂與滑移2種極限狀態(tài),但由于偶應(yīng)力的存在導(dǎo)致剪應(yīng)力不互等,其莫爾圓如圖2所示(σS與σN分別代表切向應(yīng)力和法向應(yīng)力)。

圖2 細觀物質(zhì)點的莫爾圓Fig.2 Mohr's circle of meso material point

為了方便,把應(yīng)力張量分成對稱部分和反對稱部分[12],其中對稱部分就是傳統(tǒng)意義下的應(yīng)力張量,又可分解成偏量和球量,因此反映拉裂和滑移的截斷型莫爾-庫倫屈服準則直接利用對稱部分套用即可。為了數(shù)學(xué)上的方便處理,在對稱應(yīng)力張量空間下,可采用 DP (Drucker-Prager)屈服準則[13]來反映;對于傾覆這種極限狀態(tài),在偶應(yīng)力空間下描述要比反對稱應(yīng)力張量空間下方便,因為達到傾覆屈服后,偶應(yīng)力與扭率、曲率應(yīng)變之間的關(guān)系可直接套用應(yīng)力-應(yīng)變關(guān)系來處理。

對于傾覆或滾動的極限狀態(tài)可引入靜摩阻系數(shù)δmax,如果對于平面應(yīng)變問題,球形顆粒的極限平衡狀態(tài)應(yīng)為

式中:τvy、τvx——平面內(nèi)的剪應(yīng)力項;l——內(nèi)稟尺度;c——材料的黏聚力;δmax——砂漿強相顆粒相對于弱項顆粒的靜摩阻系數(shù)。

空間問題的長方體顆粒相對復(fù)雜些,考慮到棱角帶來的肋性效應(yīng),對摩阻系數(shù)加以修正,將應(yīng)力空間下DP屈服準則推廣到偶應(yīng)力空間,即:

式中:λ——實常數(shù)。s表示對相應(yīng)變量求偏微分。

對于對稱應(yīng)力空間,彈塑性部分的應(yīng)力增量計算同連續(xù)介質(zhì)力學(xué)處理方式一樣。對于偶應(yīng)力空間,采用理想塑性模型,同樣采用廣義中點法[14]算出該迭代步的真實偶應(yīng)力,然后代入平衡方程(2),算出真實的反對稱應(yīng)力部分。

4 算例分析

通過編寫的基于Cosserat連續(xù)體介質(zhì)力學(xué)模型的有限元分析程序,進行混凝土試樣單軸壓縮試驗的數(shù)值模擬,并應(yīng)用于邊坡穩(wěn)定性分析計算。

4.1 混凝土試樣單軸壓縮

某一平面應(yīng)變條件下的混凝土板,大小為40 cm×40 cm,上下端放置剛性板,施加位移荷載,即接觸面上各節(jié)點水平位移為0,垂直位移相同,取1/4(20 cm×20 cm)進行模擬。材料參數(shù)為:E=30 GPa,μ=0.17,c= 3 MPa,δmax=1 mm。圖3、圖4分別為頂部受垂直指定位移0.04 cm時試樣的變形圖及微顆粒旋轉(zhuǎn)位移圖。

圖3 試樣變形Fig.3 Deformation patterns of specimens

圖4 旋轉(zhuǎn)位移等效云圖Fig.4 Equivalent image of rotational displacement

從圖3可以看出模擬過程中隨著壓力的逐漸增加,混凝土在受壓方向(面內(nèi))發(fā)生壓縮變形的同時產(chǎn)生橫向膨脹,混凝土顆粒發(fā)生微旋轉(zhuǎn)(圖4),物質(zhì)點的偶極效應(yīng)確實存在。

從圖5和圖6可以看出在云圖發(fā)展區(qū)域塑性應(yīng)變相對集中在試樣對角線?？紤]壓應(yīng)力與剪應(yīng)力的關(guān)系,一般工程材料滿足τs=(0.56-0.6)σs,σs為壓應(yīng)力,τs為剪應(yīng)力。算例中得到的剪切塑性應(yīng)力的值超過上述范圍,因此剪切破壞優(yōu)先發(fā)生。

圖5 等效塑性應(yīng)變Fig.5 Equivalent plastic strain

圖6 等效塑性應(yīng)力Fig.6 Equivalent plastic stress

由于物質(zhì)點存在偶極效應(yīng),即剪切導(dǎo)致混凝土顆粒旋轉(zhuǎn),產(chǎn)生局部化變形。等效塑性應(yīng)變在加載過程中不斷發(fā)展,在由屈服到破壞過程中達到峰值,最終導(dǎo)致試樣發(fā)生破壞產(chǎn)生局部變形。從細觀層面考慮:一方面,隨著剪切應(yīng)力的逐漸變大,打破了微顆粒的平衡狀態(tài),使之發(fā)生旋轉(zhuǎn),激活了顆粒的旋轉(zhuǎn)自由度;另一方面,顆粒的微旋轉(zhuǎn)又加速了局部化變形。由此可推斷出物質(zhì)點的偶極效應(yīng)存在,且是材料發(fā)生應(yīng)變局部化的直接原因。

4.2 邊坡穩(wěn)定性分析

土體材料具有明顯的峰值和殘余強度,存在應(yīng)變軟化行為,導(dǎo)致土體的邊坡穩(wěn)定性分析十分復(fù)雜。在峰值和殘余強度之間,抗剪強度會隨著應(yīng)變的增長而迅速降低,產(chǎn)生非均勻的應(yīng)變,導(dǎo)致抗剪強度非均勻變化,最終發(fā)生漸進破壞。為合理分析這種由應(yīng)變軟化引起的漸進破壞現(xiàn)象,采用基于Cosserat連續(xù)體介質(zhì)力學(xué)模型的有限元數(shù)值方法。

取某土質(zhì)邊坡,底寬40m,頂寬30m,高10m,材料參數(shù)為:E=50MPa,μ=0.3,Gc=G/3,c=50kPa,內(nèi)稟尺度l分別為1mm、2mm、5mm3種顆粒級別且隨機分布。假定材料進入塑性后服從線性應(yīng)變軟化規(guī)律[15]:σ=σy=σy0+hpεp,軟化模量hp取-150kPa。在邊坡頂部添加剛性單元,使施加的荷載壓力與實際相符(圖7)。

基于Cosserat連續(xù)體介質(zhì)力學(xué)模型的有限元分析可以追蹤漸進破壞過程,直到出現(xiàn)明顯的破壞現(xiàn)象(圖8),相比于傳統(tǒng)有限元計算,效率更高,結(jié)果更為準確。對于應(yīng)變軟化土體,在接近結(jié)構(gòu)的極限破壞狀態(tài)時,采用方形微顆粒模型能夠模擬出局部化破壞區(qū)域的位置及漸進破壞過程,進一步驗證了物質(zhì)點偶極效應(yīng)的存在與作用。

圖7 邊坡模型Fig.7 Slope model

圖8 等效塑性應(yīng)變演化過程Fig.8 Evolution process of equivalent plastic strain

5 結(jié) 論

數(shù)據(jù)模擬證明,考慮物質(zhì)點的偶極效應(yīng),在Cosserat連續(xù)體介質(zhì)力學(xué)模型框架下推導(dǎo)出的基于方形微顆粒模型下偶應(yīng)力和曲率應(yīng)變之間的線彈性關(guān)系,以及提出的考慮微顆粒傾覆破壞的廣義屈服準則是正確可行的。通過編寫基于Cosserat連續(xù)體介質(zhì)力學(xué)模型的有限元分析程序,能夠得到更加準確合理的結(jié)果。

[1]COSSERAT E,COSSERAT F.Theorie des corps deformables[M].Paris:A Hermann&Fils,1909.

[2]LI Xikui,TANG Hongxiang.A consistent return mapping algorithm for pressure-dependent elastoplastic Cosserat continua and modeling of strain localization[J].Computers and Structures,2005,83(1):1-10.

[3]NUBEL K,HUANG Wenxiong.A study of localized deformation pattern in granular media[J].Computer Methods in Applied Mechanics and Enginnering,2004,193:2719-2724.

[4]BEVERIDGE A J,WHEEL M A,NASH D H.Themicropolar elastic behavior of model macroscopically heterogeneous materials [J].International Journal of Solids and Structures,2013,50:246-255.

[5]RAMEZANI H,MOUNANGA P,JEONG J,et al.Role of cement paste composition on the self induced stress in early-age mortars: application of the Cosserat size number[J].Cement&Concrete Composites,2013,39:43-59.

[6]KHOEI A R,YADEGARI S,BIABANAKI S O R.3D finite element modeling of shear band localizatonin via the micro-polar Cosserat continuum theory[J].Computational Materials Science,2010,49:720-733.

[7]ZHANG Hongwu,WANG Hui,CHEN Biaosong,et al.Parametric variational principle based elastic-plastic analysis of cosserat continuum[J].Acta Mechanica Solida Sinica,2007,20(1):65-74.

[8]張洪武.參變量變分原理與材料和結(jié)構(gòu)力學(xué)分析[M].北京:科學(xué)出版社,2010:132-136.

[9]徐芝綸.彈性力學(xué)[M].4版.北京:高等教育出版社,2006:246-249.

[10]李育超,凌道盛,陳云敏.Cosserat連續(xù)介質(zhì)的Mohr-Coulomb屈服準則及其應(yīng)用[J].浙江大學(xué)學(xué)報:工學(xué)版,2005,39(2): 253-258.(LI Yuchao,LING Daosheng,CHEN Yunmin.The Cosserat continuum Mohr-Coulomb yield criterion and its application [J].Journal of Zhejiang University:Engineering Science,2005,39(2):253-258.(in Chinese))

[11]馬剛,常曉林,周偉,等.基于Cosserat理論的重力壩深層抗滑穩(wěn)定分析[J].巖土力學(xué),2012,33(5):1505-1512.(MA Gang,CHANG Xiaolin,ZHOU Wei,et al.Analysis of stability against deep sliding of gravity dam based on Cosserat theory[J].Rock and Soil Mechanics,2012,33(5):1505-1512.(in Chinese))

[12]王勝軍,陳萬吉.關(guān)于彈性偶應(yīng)力C0和C1理論及有限元[D].大連:大連理工大學(xué),2008.

[13]李錫夔,唐洪祥.壓力相關(guān)彈塑性Cosserat連續(xù)體模型與應(yīng)變局部化有限元模擬[J].巖石力學(xué)與工程學(xué)報,2005,24(9): 1497-1505.(LI Xikui,TANG Hongxiang.The pressure dependent elastoplastic Cosserat continuum model and finite element simulation of strain localization[J].Chinese Journal of Rock Mechanics and Engineering,2005,24(9):1497-1505.(in Chinese))

[14]王勖成,邵敏.有限單元法基本原理和數(shù)值方法[M].2版.北京:清華大學(xué)出版社,1997:448-450.

[15]李錫夔,唐洪祥.基于Cosserat連續(xù)體模型的應(yīng)變局部化有限元模擬[D].大連:大連理工大學(xué),2005.

Dipole effect analysis of granular materials(such as concrete)

SHAO Dongdong,ZHANG Qing,XIA Xiaozhou
(College of Mechanics and Materials,Hohai University,Nanjing 210098,China)

Within the framework of the Cosserat continuum medium mechanics model,the spherical micro particle is extended and transformed into a square micro particle with consideration of the dipole effect of material points.Based on the Cosserat continuum medium mechanics model,a linear elastic relationship between the couple stress and the curvature strain for the square micro particle is derived,and a generalized yield criterion that can consider the overturning failure of micro particles is put forward.With the finite element analysis program based on the Cosserat continuum medium mechanics model,a numerical simulation was performed in order to analyze the concrete specimens under uniaxial compression and the soil slope stability.The results show that the dipole effect of material points indeed exists,and it is a direct cause of strain localization of materials.

dipole effect;square micro particles;strain localization;numerical simulation

TV314

:A

:1000-1980(2014)06-0513-05

10.3876/j.issn.1000-1980.2014.06.009

2013-11 12

水文水資源與水利工程科學(xué)國家重點實驗室開放基金(2011490911);國家自然科學(xué)基金(10972072,51179064,11132003)

邵冬冬(1989—),男,江蘇南通人,碩士研究生,主要從事細觀本構(gòu)建摸和數(shù)值模擬研究。E-mail:sbsdd@126.com

章青,教授。E-mail:lxzhangqing@hhu.edu.cn