馮 君,孫巍巍,張武滿,劉志林,王曉鳴

(1.南京理工大學(xué) 瞬態(tài)物理國家重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210094;2.南京理工大學(xué) 機(jī)械工程學(xué)院,江蘇 南京 210094; 3.南京理工大學(xué) 土木工程系,江蘇 南京 210094;4.北京航空航天大學(xué) 土木工程系,北京 100191)

離散元模擬鋼筋混凝土靶板貫穿響應(yīng)研究

馮 君1,2,孫巍巍3,張武滿4,劉志林2,王曉鳴2

(1.南京理工大學(xué) 瞬態(tài)物理國家重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210094;2.南京理工大學(xué) 機(jī)械工程學(xué)院,江蘇 南京 210094; 3.南京理工大學(xué) 土木工程系,江蘇 南京 210094;4.北京航空航天大學(xué) 土木工程系,北京 100191)

提出了圓截面梁單元高斯積分點(diǎn)分布方案以及一種鋼筋-基體接觸作用模型,基于LDPM離散單元建立了鋼筋混凝土侵徹?cái)?shù)值計(jì)算模型。模擬剛性彈貫穿48 MPa和140 MPa壓縮強(qiáng)度的鋼筋混凝土靶板,通過對比彈體剩余速度和靶板破壞形態(tài),驗(yàn)證了模型對于鋼筋混凝土貫穿問題的適用性。仿真結(jié)果表明,140 MPa強(qiáng)度混凝土靶板內(nèi)鋼筋對于彈體作用更強(qiáng),對于出靶剩余速度影響更大。對比不同彈著點(diǎn)和鋼筋尺寸的貫穿仿真,彈著點(diǎn)在一根鋼筋位置和兩根鋼筋交叉處,出靶速度分別降低了約12 m/s和45 m/s;彈著點(diǎn)位置在鋼筋交叉處時(shí),通過增大鋼筋尺寸提高配筋率能夠明顯降低貫穿剩余速度。

鋼筋混凝土;彈著點(diǎn);貫穿;離散單元

隨著鉆地武器發(fā)展及防護(hù)結(jié)構(gòu)評估需要,鋼筋混凝土靶的侵徹貫穿問題成為當(dāng)前沖擊工程領(lǐng)域的研究熱點(diǎn)。彈體對鋼筋混凝土侵徹和貫穿是一個(gè)伴隨著靶體材料大變形、高應(yīng)變率和高壓的過程,涉及斷裂、擠壓、破碎等復(fù)雜響應(yīng)[1]。

關(guān)于有限厚鋼筋混凝土板的抗侵徹貫穿問題研究,前人已取得大量成果。吳飚等[2]通過石砌鋼筋混凝土結(jié)構(gòu)板侵徹實(shí)驗(yàn)研究,發(fā)現(xiàn)鋼筋在混凝土內(nèi)部的排布結(jié)構(gòu)相關(guān);馬愛娥等[1]結(jié)合TCK和HJC模型提出了混合型連續(xù)動態(tài)損傷積累模型,通過共節(jié)點(diǎn)的方式建立了鋼筋混凝土仿真模型,該模型能夠較好得到?jīng)_擊貫穿導(dǎo)致的混凝土板內(nèi)部拉壓損傷分布和破壞特征。劉志林等[3]將侵徹過程中彈丸與鋼筋的相互作用簡化為接觸點(diǎn)在阻尼介質(zhì)中受到?jīng)_擊載荷作用的響應(yīng)問題,得到了鋼筋對于彈丸的侵徹阻力,結(jié)合空腔膨脹分析得到了較為完整、合理的鋼筋混凝土侵徹阻力模型。

近年來,隨著多尺度模擬等計(jì)算力學(xué)技術(shù)的發(fā)展,出現(xiàn)了越來越多的數(shù)值算法以及更加準(zhǔn)確的本構(gòu)模型,因此數(shù)值模擬成為研究混凝土貫穿的重要手段。混凝土等準(zhǔn)脆性材料在沖擊載荷作用下的破壞,是材料從連續(xù)到不連續(xù)的過程,而離散元法對不連續(xù)體破壞過程的模擬往往具有一定的優(yōu)勢[4]。2011年,Cusatis等[4-5]開發(fā)出了一種模擬水泥基體內(nèi)部骨料相互作用的細(xì)觀離散元模型(LDPM)。根據(jù)混凝土骨料級配信息在幾何體內(nèi)部隨機(jī)生成骨料顆粒,并在其周圍構(gòu)建潛在裂紋微面以及細(xì)觀本構(gòu)關(guān)系。LDPM能夠模擬混凝土材料斷裂、多軸響應(yīng)[5],并在結(jié)構(gòu)受到侵徹與爆炸等動態(tài)載荷工況的模擬中獲得較好結(jié)果[6-7]。對比有限元模型中采用刪除單元來模擬侵徹、裂紋等現(xiàn)象,離散元模型避免了單元刪除帶來的質(zhì)量、能量不守恒。本文以LDPM為基礎(chǔ),提出梁單元圓截面高斯積分點(diǎn)分布方案和一種鋼筋-基體滑移接觸模型,建立了鋼筋混凝土侵徹?cái)?shù)值模型。通過對鋼筋混凝土靶板貫穿響應(yīng)的仿真驗(yàn)證了模型的準(zhǔn)確性,并對彈著點(diǎn)、鋼筋等影響因素進(jìn)行了進(jìn)一步的數(shù)值模擬研究。

1 LDPM模型簡介

1.1 LDPM離散元基本假設(shè)

首先,根據(jù)Fuller曲線將骨料顆粒等效成不同直徑的圓球,隨機(jī)投放在混凝土體內(nèi)部。相比于骨料,水泥基體材料以及骨料與基體交界面的強(qiáng)度和剛度都要弱很多,LDPM假設(shè)混凝土裂紋都發(fā)生在基體和基體-骨料接觸面,通過Delaunay技術(shù)等生成三角形潛在裂紋微面包絡(luò)骨料顆粒圓球的LDPM單元[4,8],如圖1所示。

采用剛體動力學(xué)描述細(xì)觀模型系統(tǒng)的變形,將細(xì)觀應(yīng)變ε用潛在裂紋微面中心位置的位移uc來表示,在法向N和剪切方向L、M的投影量εN,εL,εM分別為

(1)

1.2LDPM本構(gòu)關(guān)系

在彈性響應(yīng)階段,法向和切向細(xì)觀應(yīng)力與對應(yīng)的細(xì)觀應(yīng)變的線性關(guān)系:σN=ENεN;σM=ETεM;σL=ETεL,式中:EN=E0,ET=αE0,E0為等效法向彈性模量,α為切向-法向耦合系數(shù)。對于非線性響應(yīng)階段,LDPM主要考慮以下3個(gè)方面的細(xì)觀響應(yīng):斷裂、空隙坍塌和壓實(shí)以及壓剪耦合內(nèi)摩擦[4]。

法向壓縮強(qiáng)度滿足:-σbc(εD,εV)≤σN≤0。式中σbc為體應(yīng)變相關(guān)的邊界條件,偏應(yīng)變εD=εN-εV,εV為細(xì)觀體應(yīng)變。彈性響應(yīng)過后,由于孔隙坍塌導(dǎo)致法向應(yīng)力以較小的剛度Hc隨著體應(yīng)變增大而增大,σbc=σc0+〈-εV-εc0〉Hc(rDV),σc0和εc0分別為細(xì)觀壓縮強(qiáng)度及其相應(yīng)的應(yīng)變值,Hc(rDV)=Hc0/(1+κc2〈rDV-κc1〉),rDV=εD/εV,κc1和κc2為材料參數(shù)。當(dāng)-εV大于臨界值εc1時(shí),密實(shí)硬化現(xiàn)象發(fā)生,σbc=σc1(rDV)exp{(-εV-εc1)×Hc(rDV)/σc1(rDV)},考慮剪脹效應(yīng)的壓縮強(qiáng)度σc1(rDV)=σc0(rDV)+(εc1-εc0)×Hc(rDV)。

基于粘聚斷裂理論,LDPM中裂紋開裂速率相關(guān)的動態(tài)屈服準(zhǔn)則為

(2)

2 鋼筋混凝土建模

2.1 圓截面梁單元高斯積分點(diǎn)分布方案

鋼筋混凝土材料模擬中鋼筋通常采用梁單元進(jìn)行模擬,由于鋼筋通常為圓形截面,為了更準(zhǔn)確地描述復(fù)雜受力鋼筋梁單元,需要在圓形截面上定義若干高斯積分點(diǎn)。傳統(tǒng)梁單元高斯積分點(diǎn)分布不夠合理,往往集中在面心位置,當(dāng)梁單元受到剪力時(shí)往往精度不夠。每個(gè)積分點(diǎn)代表其附近區(qū)域的應(yīng)力狀態(tài),因此需要將圓形截面離散成若干接近于圓形或正方形的區(qū)域,選擇形心位置為高斯積分點(diǎn)位置。本文提出的方案:將圓形截面區(qū)域的高斯積分點(diǎn)分布做如下設(shè)定,將半徑為R的圓形徑向等分成n等份,其中心處圓等分成4份,如圖2所示,則第i個(gè)圓環(huán)等分成4(2i-1)份,則每個(gè)劃分區(qū)域有相等的面積且保證每個(gè)區(qū)域最接近于正方形。每個(gè)劃分區(qū)域的形心位置為高斯積分點(diǎn)位置,此處的應(yīng)力、應(yīng)變代表整個(gè)劃分區(qū)域的均值。第i個(gè)圓環(huán)上每個(gè)劃分域的積分點(diǎn)極坐標(biāo)表示為

(3)

式中:A=πR2/(4n2)為每個(gè)劃分區(qū)域的面積,αi=2π/[4(2i-1)]為每個(gè)劃分區(qū)域的弧度。

2.2 鋼筋-混凝土接觸模型

數(shù)值仿真鋼筋混凝土往往通過混凝土與梁單元共節(jié)點(diǎn)的方式進(jìn)行建模[1],該方法雖然將鋼筋和混凝土材料進(jìn)行了一定程度的耦合,但是無法描述鋼筋與基體的脫粘、滑移等接觸關(guān)系。為更加準(zhǔn)確描述鋼筋對混凝土基體的作用,本文根據(jù)鋼筋拉拔機(jī)理,提出一種考慮鋼筋與基體間相對滑移、分離現(xiàn)象的接觸模型。

鋼筋與混凝土之間作用力在彈性階段時(shí)軸向應(yīng)力σa和徑向應(yīng)力σr滿足關(guān)系:σa=Kasa,σr=Krsr。式中:sa和sr分別為軸向和徑向位移,軸向剛度系數(shù)Ka=λb10/s10,Kr為徑向剛度系數(shù),λ為粘結(jié)強(qiáng)度形狀系數(shù),b10為初始黏結(jié)強(qiáng)度,s10為黏結(jié)力彈性極限對應(yīng)的滑移量。在鋼筋與基體脫離接觸之前,滑移量相關(guān)的接觸力模型:徑向應(yīng)力對于軸向黏結(jié)強(qiáng)度影響關(guān)系為b1=〈b10-μσr〉,滑移強(qiáng)度b3=〈b30-μσr〉,b30為初始滑移強(qiáng)度,μ為摩擦系數(shù)。黏結(jié)應(yīng)力σb定義為當(dāng)前滑移量和當(dāng)前軸向作用強(qiáng)度的函數(shù):

(4)

式中:s1,s2,s3分別為徑向力作用下實(shí)際的接觸關(guān)系特征滑移量,s1=λb1/Ka,s2=2s20-s1-s10,s3=s30+s20-s1-s10,s20為無側(cè)向作用時(shí)脫粘開始階段對應(yīng)的滑移量,s30為無側(cè)向作用時(shí)完全脫粘對應(yīng)的滑移量。

當(dāng)前軸向作用強(qiáng)度σb=0時(shí),表示鋼筋和基體材料脫黏,此后兩者接觸力軸向分量純粹為庫倫摩擦作用,σa=μσr。

參考Leibovich等[10]通過鋼筋拉拔實(shí)驗(yàn)獲得的拉拔相應(yīng)情況特征,本研究取b10=12 MPa,b30=7 MPa,λ=2。結(jié)合Leibovich等[10]所進(jìn)行的鋼筋拉拔實(shí)驗(yàn)中拉拔載荷與鋼筋滑移量的關(guān)系,本模型其他參數(shù)選擇如下:s10=0.5 mm,s20=1.2 mm,s30=2.8 mm,Kr=1.0 GPa。

3 LDPM模擬Hanchak貫穿實(shí)驗(yàn)

3.1 彈體侵徹模型

[11]中小直徑動能彈貫穿48 MPa和140 MPa強(qiáng)度鋼筋混凝土靶板的實(shí)驗(yàn),如圖3所示,鋼筋混凝土靶板尺寸為610 mm×610 mm×178 mm,配有三層正交鋼筋, 混凝土上、下層保護(hù)厚度均為12.7 mm,鋼筋層間距76.2 mm,鋼筋直徑為5.69 mm。彈體總長為145.7 mm,彈體直徑為25.4 mm(1英寸),尖卵形頭部彈丸的頭部曲率半徑比為3。

由于本文模擬的侵徹貫穿實(shí)驗(yàn)后彈體頭部幾乎沒有形狀變化和質(zhì)量損失,為了提高計(jì)算效率,本文選擇剛性彈體對下文貫穿進(jìn)行仿真分析。鋼筋采用彈塑性模型描述,密度為7.8g/cm3,泊松比為0.3,楊氏模量為200GPa,拉伸屈服強(qiáng)度為500MPa,剪切屈服強(qiáng)度為300MPa。

彈體與混凝土和鋼筋之間接觸設(shè)置為罰函數(shù)接觸,其接觸剛度Ksif為時(shí)間步長Δt相關(guān)的函數(shù)。

(5)

式中:mmin為發(fā)生接觸的顆?；蛘邌卧匈|(zhì)量最小值,Δt為顯示計(jì)算步中的時(shí)間步長,計(jì)算步長越小罰函數(shù)剛度越大,對一定相互作用力則產(chǎn)生更小的相對位移量。結(jié)合式(6),給出罰函數(shù)接觸力Fn=Ksifp,p為接觸體相互侵入量。

表1 LDPM參數(shù)取值

3.2 48MPa強(qiáng)度混凝土板貫穿仿真

對48MPa強(qiáng)度鋼筋混凝土靶板貫穿進(jìn)行模擬,結(jié)果如表2所示。彈體撞擊速度vs在434～1 058m/s范圍內(nèi),N-3號實(shí)驗(yàn)彈著點(diǎn)在鋼筋交叉位置,其他工況彈著點(diǎn)都在鋼筋網(wǎng)眼處。對比數(shù)值模型仿真得到的彈體剩余速度vrs與實(shí)驗(yàn)測得的出靶速度vr,在434m/s撞擊速度時(shí)模擬結(jié)果較大,隨著彈體撞擊速度繼續(xù)增大,貫穿出靶剩余速度仿真值與實(shí)驗(yàn)值相對誤差δ在4.8%以內(nèi)。

表2 48 MPa強(qiáng)度混凝土板貫穿剩余速度對比

本文給出的模型參數(shù),只對模擬出靶速度大于470m/s的工況有較高精度。實(shí)驗(yàn)N-3和N-4對比了相近撞擊速度下不同彈著點(diǎn)位置對于出靶速度的影響,貫穿實(shí)驗(yàn)結(jié)果表明彈體出靶速度相差10m/s;而本文數(shù)值仿真結(jié)果也顯示,在746m/s撞擊速度的情況下,彈著點(diǎn)位置引起的出靶速度相差14m/s。

N-2號實(shí)驗(yàn)彈體貫穿過程仿真結(jié)果如圖4所示,圖4(a)～圖4(d)分別為撞擊后0.1ms,0.2ms,0.4ms和0.6ms時(shí)刻,鋼筋混凝土靶體內(nèi)部剖面裂紋破壞情況。在0.1ms時(shí),彈體進(jìn)入靶板內(nèi)部并在周圍產(chǎn)生開坑和隧道;在0.2ms時(shí),彈體仍然在靶板內(nèi)部侵徹,但由于靶背自由面無約束,剖面上觀察到靶背附近產(chǎn)生了部分裂紋;0.4ms時(shí),彈體頭部開始露出靶背,整個(gè)侵徹隧道區(qū)形成;在0.6ms時(shí),靶前、靶后呈現(xiàn)出開坑后崩落區(qū)域,但由于鋼筋網(wǎng)的約束作用,基體碎片崩落不明顯。

N-4號工況下,靶前的實(shí)驗(yàn)與仿真結(jié)果對比如圖5所示,靶背的對比如圖6所示。由于沖擊產(chǎn)生的壓縮波和邊界反射的拉伸波共同作用,在靶板迎彈面留下方形開坑,仿真結(jié)果顯示其開坑面積基本一致。而靶背由于主要受拉伸反射波以及彈體沖塞作用,仿真結(jié)果顯示彈著點(diǎn)附近骨料顆粒被彈體沖塞帶走脫離靶板。對比圖5(a)和圖6(a)中貫穿實(shí)驗(yàn)后靶板表面的破壞情況,圖5(b)和圖6(b)的仿真結(jié)果呈現(xiàn)更多的輻射狀裂紋。

3.3 140MPa強(qiáng)度混凝土板貫穿仿真

為了研究靶體強(qiáng)度對于貫穿響應(yīng)的影響,Hanchak等[11]對同樣內(nèi)部結(jié)構(gòu)的140MPa強(qiáng)度鋼筋混凝土靶板進(jìn)行443～998m/s撞擊速度的侵徹實(shí)驗(yàn),實(shí)驗(yàn)中通過控制炮口和靶板位置,保證彈體著靶位置都在鋼筋網(wǎng)眼中心處。本文模擬結(jié)果如表3所示,與上一節(jié)仿真規(guī)律類似,彈體撞擊速度為443m/s時(shí),仿真所得出靶速度比實(shí)驗(yàn)結(jié)果略小,而彈體撞擊速度在552～998m/s之間時(shí),仿真獲得的彈體剩余速度較準(zhǔn)確,相對于實(shí)驗(yàn)值的誤差在5.1%以下。本文采用的仿真參數(shù)只對出靶速度超過260m/s的工況的預(yù)測較為精確。這可能是因?yàn)槭Ｓ嗨俣容^小的工況下,靶背沖塞過程響應(yīng)對混凝土拉伸特性較敏感,由于缺少實(shí)驗(yàn)數(shù)據(jù),本文無法進(jìn)行相關(guān)參數(shù)標(biāo)定。

表3 140 MPa強(qiáng)度混凝土板貫穿剩余速度對比

彈體以不同撞擊速度貫穿140MPa強(qiáng)度鋼筋混凝土靶板的侵徹過程中,彈體的減速度ap時(shí)程曲線如圖7所示。彈體撞擊速度越大,彈體過載曲線峰值越大,侵徹貫穿時(shí)程越短。998m/s撞擊速度工況下,彈體速度衰減變化不大,因此與速度相關(guān)的侵徹阻力變化不大,過載曲線平臺段占主要侵徹時(shí)程。而撞擊速度較小的工況下,由于速度衰減較明顯且存在靶背沖塞現(xiàn)象,過載曲線在侵徹時(shí)程后半段明顯減小。特別對于443m/s撞擊速度工況,仿真中彈體出靶后繼續(xù)和靶背沖塞飛出的骨料相互作用,彈體過載隨著時(shí)間緩慢減小到0。

H-4號實(shí)驗(yàn)仿真結(jié)果對比如圖8和圖9所示,本文模型仿真得到的140MPa強(qiáng)度混凝土靶板迎彈面開坑和靶背沖塞破壞形態(tài)與實(shí)驗(yàn)基本吻合。由于140MPa強(qiáng)度混凝土材料更脆,模擬得到迎彈面開坑比N-4號實(shí)驗(yàn)仿真結(jié)果更大。靶體迎彈面和背面破壞形態(tài)相比,仿真結(jié)果顯示前者主要是由于彈體擠壓靶塊剪切破壞,而后者由于拉伸波和彈體沖塞作用呈現(xiàn)粉碎破壞。

4 140 MPa強(qiáng)度混凝土板貫穿響應(yīng)仿真分析

4.1 彈著點(diǎn)位置的影響

由于侵徹實(shí)驗(yàn)中彈著點(diǎn)與鋼筋的相對位置很難控制,彈著點(diǎn)對鋼筋混凝土侵徹貫穿響應(yīng)影響一般很難用實(shí)驗(yàn)驗(yàn)證。孫其然等[12]采用HJC模型仿真分析了48MPa強(qiáng)度鋼筋混凝土板貫穿問題中不同彈著點(diǎn)位置對于出靶速度的影響,其結(jié)論與實(shí)驗(yàn)結(jié)果[11]類似,該侵徹工況下,是否撞擊鋼筋對于彈體終點(diǎn)性能影響不大。本文結(jié)合前文140MPa強(qiáng)度混凝土靶板貫穿實(shí)驗(yàn)工況,對圖10中3種典型彈著點(diǎn)位置的情況進(jìn)行數(shù)值模擬分析。彌補(bǔ)前人鋼筋混凝土貫穿研究中高強(qiáng)度混凝土靶板彈著點(diǎn)相對位置影響規(guī)律這一空白。

表4給出了不同撞擊速度vs侵徹貫穿鋼筋混凝土靶板3種不同彈著點(diǎn)位置工況下(如圖10所示)彈體的出靶速度。相比于彈著點(diǎn)1的貫穿工況,彈體撞擊一根鋼筋(彈著點(diǎn)2)仿真得到的出靶速度分別減小了16m/s,10m/s,5m/s,14m/s和17m/s;若彈體撞擊鋼筋交叉點(diǎn),出靶速度相對彈著點(diǎn)在鋼筋網(wǎng)眼的工況則分別減小46m/s,55m/s,44m/s,46m/s和32m/s。考慮到離散元模型仿真中帶來的一定隨機(jī)性,本文認(rèn)為彈著點(diǎn)對于彈體剩余速度影響的變化規(guī)律與撞擊速度無關(guān)。由于鋼筋的作用,彈著點(diǎn)在一根鋼筋和兩根鋼筋交叉位置與鋼筋網(wǎng)眼位置相比,出靶速度分別降低了約12m/s和45m/s。對比48MPa強(qiáng)度鋼筋混凝土貫穿實(shí)驗(yàn)[11]以及孫其然等[12]的仿真研究,本文數(shù)值模擬結(jié)果表明,彈著點(diǎn)位置對于貫穿剩余速度的影響較明顯。參考劉志林等[3]提出的鋼筋混凝土侵徹工程分析模型,彈體與鋼筋相互作用可以簡化為固支梁在阻尼介質(zhì)中受到的沖擊載荷作用,則鋼筋對于彈體侵徹阻力作用中一部分來自混凝土對鋼筋的阻力。隨著混凝土強(qiáng)度的增大,鋼筋對于彈體阻力也一定程度增大,這與本文的模擬結(jié)果相符合。

表4 不同彈著點(diǎn)位置剩余速度對比

針對743m/s撞擊速度工況的仿真,圖11給出不同彈著點(diǎn)位置對應(yīng)的彈體減速度時(shí)程曲線。由于鋼筋和混凝土骨料之間的相互作用,彈體撞擊第一層鋼筋時(shí),其減速度變化沒有特別突出的體現(xiàn),而在0.2ms和0.3ms附近,第二層、第三層鋼筋作用于彈體的阻力對其減速度產(chǎn)生明顯的擾動,這是由于隨著彈體速度減小,彈體受到鋼筋的阻力占彈體受到的總阻力比重越來越大。

彈體以998m/s速度撞擊彈著點(diǎn)3位置,貫穿行程仿真結(jié)果如圖12所示。鋼筋受到彈體撞擊、擠壓以及混凝土粘結(jié)的共同作用,發(fā)生塑性變形,成為彈體截面圓形,第一層鋼筋向靶前方向移動,第二層、第三層鋼筋向靶后方向產(chǎn)生形變。與文獻(xiàn)[3]中鋼筋混凝土侵徹實(shí)驗(yàn)回收進(jìn)行對比分析,本文模擬鋼筋變形與圖13實(shí)驗(yàn)觀察現(xiàn)象類似。

4.2 鋼筋對于混凝土靶板貫穿響應(yīng)影響

前文通過LDPM對鋼筋混凝土靶板貫穿問題中彈著點(diǎn)相對位置進(jìn)行了討論。本節(jié)對于彈著點(diǎn)在鋼筋網(wǎng)眼位置和鋼筋交叉位置,分別進(jìn)一步仿真分析鋼筋的影響規(guī)律。圖14為鋼筋混凝土靶和素混凝土靶貫穿后破壞形態(tài)對比圖。

1)彈著點(diǎn)1位置鋼筋的影響。 為研究鋼筋在彈著點(diǎn)1位置工況中的作用,對比仿真彈著點(diǎn)在鋼筋網(wǎng)眼位置以及素混凝土侵徹的響應(yīng)規(guī)律。

圖14中,彈著點(diǎn)在鋼筋網(wǎng)眼處,鋼筋對于整個(gè)結(jié)構(gòu)仍然起到粘結(jié)作用。對比素混凝土裂紋集中整體結(jié)構(gòu)破壞,鋼筋混凝土靶板貫穿后裂紋集中不明顯。由于鋼筋的約束作用,靶前、靶后開坑和崩落的情況也極大地改善。圖15是5種撞擊速度工況中彈體速度vp隨著彈體位移sp的變化曲線,實(shí)心線代表鋼筋混凝土工況,其他線是素混凝土貫穿響應(yīng)。不同撞擊速度下,彈體貫穿鋼筋混凝土和素混凝土靶體的過程中,彈體速度變化幾乎重合。因此,小直徑彈丸侵徹鋼筋混凝土靶的彈著點(diǎn)在鋼筋網(wǎng)眼時(shí),貫穿過程中彈體速度變化時(shí)程曲線與貫穿相應(yīng)的素混凝土靶情況基本一致。劉志林等[3]提出鋼筋混凝土侵徹分析模型,認(rèn)為鋼筋混凝土的侵徹阻力來自混凝土阻力與鋼筋阻力之和,其中鋼筋對于混凝土的阻力沒有貢獻(xiàn)。因此該模型應(yīng)用在小直徑彈丸侵徹彈著點(diǎn)位置在鋼筋網(wǎng)眼處的情況,得到的侵徹阻力與侵徹素混凝土彈體的響應(yīng)一致。這與本文模擬結(jié)果基本相符。

2)彈著點(diǎn)3位置鋼筋尺寸的影響。 彈著點(diǎn)在鋼筋交叉處,鋼筋會明顯降低彈體的出靶剩余速度。本節(jié)針對這一工況進(jìn)一步仿真研究鋼筋尺寸對于鋼筋混凝土貫穿彈體出靶速度的影響。

表5給出了鋼筋混凝土靶板中配筋直徑尺寸dj為4mm,5.7mm和7mm的3種情況下,貫穿仿真得到的彈體出靶速度對比。隨著鋼筋直徑的增大,不同撞擊速度工況下出靶速度都明顯降低。將配筋直徑從4mm提高到5.7mm,出靶速度平均減小了17m/s,繼續(xù)提高到7mm時(shí),出靶速度大概降低了25m/s。當(dāng)彈體撞擊速度較低時(shí),鋼筋直徑增大對于出靶速度影響較大。

表5 不同鋼筋尺寸下鋼筋混凝土靶貫穿剩余速度對比

5 結(jié)束語

基于LDPM離散元,本文建立了鋼筋混凝土滑移接觸模型,通過對兩組鋼筋混凝土貫穿進(jìn)行仿真,驗(yàn)證了所研究的數(shù)值模型對于模擬侵徹問題的準(zhǔn)確性。對不同工況下影響因素進(jìn)行仿真分析,得到如下結(jié)論:

①對于本文研究的48MPa和140MPa強(qiáng)度鋼筋混凝土貫穿問題,在出靶速度較大的情況下,LDPM能夠較準(zhǔn)確地模擬彈體出靶速度以及靶體破壞形態(tài)。

②相比48MPa強(qiáng)度靶,140MPa強(qiáng)度鋼筋混凝土靶板貫穿工況下,彈著點(diǎn)位置對于出靶速度影響更大。彈著點(diǎn)位置在一根鋼筋位置和兩根鋼筋交叉位置,彈體剩余速度分別降低了約12m/s和45m/s。

③相比素混凝土靶,鋼筋混凝土靶能夠有效地避免彈體貫穿導(dǎo)致的靶體結(jié)構(gòu)破壞,并能減小開坑和崩落區(qū)域面積;而當(dāng)彈著點(diǎn)位置在鋼筋網(wǎng)眼處時(shí),貫穿178mm厚素混凝土靶板和鋼混凝土靶板,出靶速度幾乎沒有區(qū)別。

④彈體著靶位置在鋼筋交叉處時(shí),通過增大鋼筋尺寸提高配筋率能夠明顯降低貫穿剩余速度。

參考文獻(xiàn)

[1] 馬愛娥,黃風(fēng)雷,李金柱,等.鋼筋混凝土抗貫穿數(shù)值模[J].北京理工大學(xué)學(xué)報(bào),2007,27(2):103-107. MA Ai-e,HUANG Feng-lei,LI Jin-zhu,et al.Simulation on the perforation of reforced concrete[J].Transactions of Beijing Institute of Technology,2007,27(2):103-107.(in Chinese)

[2] 吳飚,楊建超,劉瑞朝.有限厚塊石砌體鋼筋混凝土結(jié)構(gòu)板抗貫穿性能的實(shí)驗(yàn)研究[J].爆炸與沖擊,2013,33(1):73-78. WU Biao,YANG Jian-chao,LIU Rui-chao.Experimental study on perforation resistance of composite targets composed by granite block masonry and reinforced concrete plates[J].Explosion and Shock Waves,2013,33(1):73-78.(in Chinese)

[3] 劉志林,孫巍巍,王曉鳴,等.卵形彈體垂直侵徹鋼筋混凝土靶的工程解析模型[J].彈道學(xué)報(bào),2015,27(3):84-90. LIU Zhi-lin,SUN Wei-wei,WANG Xiao-ming,et al.Engineering analytical model of ogive-nose steel projectiles vertically penetrating reinforced concrete target[J].Journal of Ballistics,2015,27(3):84-90.(in Chinese)

[4] CUSATIS G,PELESSONE D,MENCARELLI A.Lattice discrete particle model(LDPM)for failure behavior of concrete.I:Theory[J].Cement and Concrete Composites,2011,33(9):881-890.

[5] CUSATIS G,MENCARELLI A,PELESSONE D,et al.Lattice discrete particle model(LDPM)for failure behavior of concrete.II:Calibration and validation[J].Cement and Concrete Composites,2011,33(9):891-905.

[6] SMITH J,CUSATIS G,PELESSONE D,et al.Discrete modeling of ultra-high-performance concrete with application to projectile penetration[J].International Journal of Impact Engineering,2014,65:13-32.

[7] FENG J,SUN W,LIU Z,et al.An armour-piercing projectile penetration in a double-layered target of ultra-high-performance fiber reinforced concrete and armour steel:Experimental and numerical analyses[J].Materials & Design,2016,102:131-141.

[8] WAN L,WENDNER R,LIANG B,et al.Analysis of the behavior of ultra high performance concrete at early age[J].Cement and Concrete Composites,2016,74:120-135.

[9] REZAKHANI R,CUSATIS G.Asymptotic expansion homogenization of discrete fine-scale models with rotational degrees of freedom for the simulation of quasi-brittle materials[J].Journal of the Mechanics and Physics of Solids,2016,88:320-334.

[10] LEIBOVICH O,DANCYGIER A N,YANKELEVSKY D Z.An innovative experimental procedure to study local rebar-concrete bond by direct observations and measurements[J].Experimental Mechanics,2016,56(5):673-682.

[11] HANCHAK S J,FORRESTAL M J,YOUNG E R,et al.Perforation of concrete slabs with 48 MPa(7 ksi)and 140 MPa(20 ksi)unconfined compressive strengths[J].International Journal of Impact Engineering,1992,12(1):1-7.

[12] 孫其然,李睿宇,趙亞運(yùn),等.HJC模型模擬鋼筋混凝土侵徹實(shí)驗(yàn)的參數(shù)研究[J].工程力學(xué),2016(8):1-9. SUN Qi-ran,LI Rui-yu,ZHAO Ya-yun,et al.Investigation on parameters of HJC model applied to simulations on perforation experiments of reinforced concrete[J].Engineering Mechanics,2016(8):1-9.(in Chinese)

Study on Perforation Responses of Reinforced Concrete Slab by Discrete Modeling

FENG Jun1,2,SUN Wei-wei3,ZHANG Wu-man4,LIU Zhi-lin2,WANG Xiao-ming2

(1.National Key Laboratory of Transient Physics,Nanjing University of Science & Technology,Nanjing 210094,China; 2.School of Mechanical Engineering,Nanjing University of Science & Technology,Nanjing 210094,China; 3.Department of Civil Engineering,Nanjing University of Science & Technology,Nanjing 210094,China; 4.Department of Civil Engineering,Beihang University,Beijing 100191,China)

A new distribution scheme of Gauss integration point and a new rebar-matrix interaction model were presented.A numerical model for rebar concrete(RC)plate perforation was established based on lattice discrete particle model(LDPM).The compressive strength of 48 MPa and 140 MPa of RC plates perforation was simulated respectively.The proposed model was validated for the perforation simulation in terms of residual velocity and damage mode.Under the conditions of 140 MPa compressive strength,the rebar of reinforced concrete has greater resistance on the impact projectile compared to the conditions of 48 MPa compressive strength.Numerical results indicate that the residual velocity decreases by 12 m/s and 45 m/s while projectile impacting one rebar or two rebar joints.It can significantly decrease the projectile residual-velocity to increase the rebar diameter while the impact position is at the junction of two rebars.

reinforced concrete;impact point;perforation;discrete element

2016-12-13

國家自然科學(xué)基金項(xiàng)目(11602111,51278250);中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金資助(30917011343)

馮君(1988- ),男,講師,博士,研究方向?yàn)槎喑叨饶M。E-mail:jun.feng@njust.edu.cn。

O385

A

1004-499X(2017)02-0026-08