吳長(zhǎng)， 楊佑佩， 茍寶龍， 張瑞芹

(1.蘭州理工大學(xué) 土木工程學(xué)院， 甘肅 蘭州 730050； 2.甘肅省土木工程防災(zāi)減災(zāi)重點(diǎn)實(shí)驗(yàn)室， 甘肅 蘭州 730050)

網(wǎng)殼結(jié)構(gòu)在沖擊荷載作用下會(huì)發(fā)生動(dòng)力失穩(wěn)甚至倒塌，一直以來(lái)是空間結(jié)構(gòu)動(dòng)力學(xué)研究領(lǐng)域中具有挑戰(zhàn)性的課題。沖擊荷載本身具有復(fù)雜性，其作為一種短暫而強(qiáng)烈的三角形脈沖荷載作用在結(jié)構(gòu)上，在極短的時(shí)間內(nèi)引起結(jié)構(gòu)的動(dòng)力響應(yīng)甚至倒塌。再者，由于網(wǎng)殼結(jié)構(gòu)本身具有復(fù)雜性和多樣性，即使是沖擊試驗(yàn)的縮尺模型也很難精確加工，且試驗(yàn)成本較高。因此，關(guān)于網(wǎng)殼結(jié)構(gòu)在沖擊荷載作用下的動(dòng)力響應(yīng)及失效模式的研究方法大多數(shù)仍然采用理論推導(dǎo)與數(shù)值模擬，而試驗(yàn)研究較少。

李海旺等[1-2]展開了K8型單層球面網(wǎng)殼的單點(diǎn)沖擊試驗(yàn)研究，對(duì)網(wǎng)殼在單點(diǎn)沖擊荷載作用下的動(dòng)力穩(wěn)定性進(jìn)行了研究，并總結(jié)了桿件的動(dòng)力響應(yīng)特點(diǎn)。沈世釗、范峰、支旭東、王多智等[3-8]對(duì)單層球面網(wǎng)殼做了大量較為系統(tǒng)的研究，分析總結(jié)了網(wǎng)殼結(jié)構(gòu)的沖擊失效模式，提出沖擊荷載下網(wǎng)殼結(jié)構(gòu)的動(dòng)力響應(yīng)分析方法，并對(duì)2個(gè)相同的K6型單層球面網(wǎng)殼進(jìn)行了單點(diǎn)沖擊試驗(yàn)研究，試驗(yàn)結(jié)果證明了網(wǎng)殼結(jié)構(gòu)在沖擊荷載作用下的2類典型失效模式的存在以及有限元分析結(jié)果的可靠性。此后，王秀麗、吳長(zhǎng)、馬肖彤等[9-11]對(duì)帶下部支撐的單層球面網(wǎng)殼進(jìn)行了沖擊試驗(yàn)，分析總結(jié)了結(jié)構(gòu)的動(dòng)力響應(yīng)規(guī)律，同時(shí)也通過(guò)試驗(yàn)驗(yàn)證了4類典型失效模式的客觀存在性。

文獻(xiàn)[12-13]對(duì)網(wǎng)殼結(jié)構(gòu)的動(dòng)力穩(wěn)定性和低速?zèng)_擊荷載下的動(dòng)力特性進(jìn)行了研究。Mohammad等[14]在分析波紋構(gòu)件沖擊性能的基礎(chǔ)上，研究了超高強(qiáng)度鋼管在沖擊荷載作用下對(duì)HC構(gòu)件性能的貢獻(xiàn)。

目前對(duì)于網(wǎng)殼結(jié)構(gòu)承載力和穩(wěn)定性的研究較為成熟[15-18]，網(wǎng)殼結(jié)構(gòu)在沖擊荷載作用下動(dòng)力響應(yīng)的有限元分析與試驗(yàn)研究上取得了較為豐碩的成果，但是無(wú)論是有限元分析還是試驗(yàn)研究，對(duì)于考慮節(jié)點(diǎn)剛度和沖擊點(diǎn)數(shù)量等因素對(duì)網(wǎng)殼結(jié)構(gòu)整體動(dòng)力響應(yīng)及失效模式的影響的研究很少，本文基于前人研究，對(duì)采用焊接空心球節(jié)點(diǎn)的K6型球面網(wǎng)殼進(jìn)行多點(diǎn)沖擊試驗(yàn)研究，以此來(lái)研究沖擊點(diǎn)數(shù)量和沖擊位置對(duì)單層球面網(wǎng)殼動(dòng)力響應(yīng)的影響。

1 試驗(yàn)概況

1.1 模型設(shè)計(jì)

試驗(yàn)?zāi)Ｐ蜑榭缍? 500 mm、矢高375 mm、頻數(shù)5環(huán)的K6型單層球面網(wǎng)殼，網(wǎng)殼通過(guò)小立柱焊接在底部環(huán)梁上，環(huán)梁為方鋼管，下部焊接2塊帶有螺栓孔的鋼板，整體試驗(yàn)?zāi)Ｐ屯ㄟ^(guò)6個(gè)M24普通螺栓連接在試驗(yàn)平臺(tái)的鋼板上，鋼板尺寸為1.9 m×1.9 m×0.02 m。試驗(yàn)?zāi)Ｐ蜆?gòu)件尺寸如表1所示，加工完成后的試驗(yàn)?zāi)Ｐ腿鐖D1所示。

表1 桿件及節(jié)點(diǎn)參數(shù)表Table 1 Parameters of rods and nodes mm

圖1 試驗(yàn)?zāi)Ｐ虵ig.1 Text model

1.2 多點(diǎn)沖擊試驗(yàn)內(nèi)容與方案

網(wǎng)殼多點(diǎn)沖擊試驗(yàn)在沖擊試驗(yàn)臺(tái)上進(jìn)行，升降平臺(tái)最高可升至8 m，平臺(tái)兩側(cè)分別架設(shè)2個(gè)12號(hào)工字鋼梁，每個(gè)鋼梁配有電動(dòng)葫蘆和電磁吸附裝置，可通過(guò)電磁吸附裝置控制沖擊物的自由下落，試驗(yàn)平臺(tái)如圖2所示。各試驗(yàn)工況方案如表2所示，本文主要研究沖擊點(diǎn)數(shù)量及沖擊位置對(duì)結(jié)構(gòu)的影響，因而每個(gè)沖擊點(diǎn)的沖擊能量應(yīng)保持一致，鋼球選用此參數(shù)的原因是所有的試驗(yàn)設(shè)計(jì)工況均為彈性沖擊試驗(yàn)，不能對(duì)結(jié)構(gòu)造成破壞，具體沖擊位置如圖3～圖6所示。

表2 動(dòng)力響應(yīng)試驗(yàn)方案Table 2 Test scheme of dynamic response

圖2 沖擊加載裝置Fig.2 Impact loading equipment

圖3 第1組工況沖擊點(diǎn)布置Fig.3 Impact point arrangement of the first group

圖4 第2組工況沖擊點(diǎn)布置Fig.4 Impact point arrangement of the second group

圖5 第3組工況沖擊點(diǎn)布置Fig.5 Impact point arrangement of the third group

圖6 第4組工況沖擊點(diǎn)布置Fig.6 Impact point arrangement of the fourth group

1.3 測(cè)點(diǎn)布置與數(shù)據(jù)采集

本次試驗(yàn)采集的數(shù)據(jù)有關(guān)鍵桿件的動(dòng)應(yīng)變及關(guān)鍵節(jié)點(diǎn)的動(dòng)位移和加速度，桿件上共布置10個(gè)應(yīng)變片，編號(hào)為S1～S10；位移傳感器共布置2個(gè)，編號(hào)分別為D1和D2；加速度傳感器共布置3個(gè)，編號(hào)依次為A1、A2、A3，具體布置方案如圖7所示，各測(cè)試元件固定安裝完成后，如圖8所示。

圖7 測(cè)點(diǎn)布置方案Fig.7 Arrangement of measure point

圖8 安裝完成后的模型Fig.8 Model after installation

1.4 精細(xì)化有限元分析模型

在ANSYS/LS-DYNA中建立試驗(yàn)?zāi)Ｐ偷木?xì)化有限元模型，焊接球與桿件均采用SHELL163單元，施加沖擊荷載的鋼球采用SOLID164單元。金屬材料承受短時(shí)超強(qiáng)載荷時(shí)，應(yīng)變率對(duì)其材料性能有較大的影響，因此，鋼材的材料模型采用分段線性塑性模型，該材料模型可考慮應(yīng)變率的影響，沖擊物材料模型采用剛性材料模型。精細(xì)化有限元分析模型如圖9所示。

圖9 數(shù)值模型Fig.9 Numerical mode

2 沖擊試驗(yàn)結(jié)果分析

2.1 應(yīng)變分析

為了研究單層球面網(wǎng)殼結(jié)構(gòu)在不同沖擊物數(shù)量及不同沖擊位置下的動(dòng)力響應(yīng)，試驗(yàn)記錄了每個(gè)測(cè)點(diǎn)的應(yīng)變隨時(shí)間變化的曲線，現(xiàn)對(duì)第1組彈性沖擊試驗(yàn)的應(yīng)變?cè)囼?yàn)數(shù)據(jù)進(jìn)行分析，各測(cè)點(diǎn)應(yīng)變片的拉應(yīng)變及壓應(yīng)變的峰值變化如圖10所示。由工況1-1至工況1-3的應(yīng)變峰值變化分析可知，隨著沖擊點(diǎn)數(shù)量的增多，各測(cè)點(diǎn)所在位置的應(yīng)變值增大，由工況1-3至工況1-4的應(yīng)變峰值曲線分析可知，各測(cè)點(diǎn)所在位置的應(yīng)變值減小。

圖10 工況拉壓應(yīng)變峰值Fig.10 The peak strain of condition

隨著工況1-1到工況1-3沖擊點(diǎn)數(shù)量增多，作用到網(wǎng)殼結(jié)構(gòu)的總能量也逐漸增大，導(dǎo)致測(cè)點(diǎn)應(yīng)變逐漸增大。由于工況1-1、1-2、1-4只有測(cè)點(diǎn)S4附近有一個(gè)沖擊點(diǎn)，可以看出測(cè)點(diǎn)S4的拉應(yīng)變明顯大于其他測(cè)點(diǎn)。工況1-3在測(cè)點(diǎn)S4、S7附近均有沖擊點(diǎn)，可以看出測(cè)點(diǎn)S4及測(cè)點(diǎn)S7處的拉應(yīng)變明顯大于其他測(cè)點(diǎn)。其中測(cè)點(diǎn)S4處的拉應(yīng)變相較于前一個(gè)工況分別增長(zhǎng)了9.58%、41.68%和-21.21%。工況1-3與工況1-4同為3點(diǎn)沖擊，但工況1-3的沖擊位置處于網(wǎng)殼結(jié)構(gòu)的同一側(cè)，沖擊點(diǎn)較為密集，而工況1-4的沖擊位置分布均勻，導(dǎo)致工況1-3沖擊區(qū)的測(cè)點(diǎn)應(yīng)變大于工況1-4的沖擊區(qū)測(cè)點(diǎn)應(yīng)變。

2.2 位移分析

結(jié)構(gòu)的位移可以直接反映結(jié)構(gòu)在沖擊過(guò)程中的變形規(guī)律，試驗(yàn)中記錄了網(wǎng)殼結(jié)構(gòu)部分關(guān)鍵節(jié)點(diǎn)的在沖擊荷載下的動(dòng)態(tài)位移。各工況下2個(gè)位移測(cè)點(diǎn)的最大位移值變化如圖11所示，圖11(a)為第1組試驗(yàn)工況下測(cè)點(diǎn)最大位移變化曲線，從工況1-1至工況1-3可以看出，隨著沖擊點(diǎn)數(shù)量的增多，D1、D2測(cè)點(diǎn)的位移逐漸增大，而工況1-4中測(cè)點(diǎn)D2位于沖擊點(diǎn)處，導(dǎo)致D2位移劇增。圖11(b)為第2組試驗(yàn)工況下測(cè)點(diǎn)最大位移變化曲線，工況2-1至工況2-2的變化趨勢(shì)與第1組相同，但工況2-3由于測(cè)點(diǎn)距沖擊點(diǎn)距離較遠(yuǎn)，位移與工況2-2中的位移相比并未出現(xiàn)明顯變化，工況2-4由于D2測(cè)點(diǎn)距一個(gè)沖擊點(diǎn)較近，導(dǎo)致測(cè)點(diǎn)位移劇增。與第1組相比，第2組由于沖擊點(diǎn)外移，與測(cè)點(diǎn)相距較遠(yuǎn)，導(dǎo)致第2組測(cè)點(diǎn)位移整體小于第1組測(cè)點(diǎn)位移。

圖11(c)為第3組試驗(yàn)工況下測(cè)點(diǎn)最大位移變化曲線，測(cè)點(diǎn)D2距沖擊區(qū)距離比測(cè)點(diǎn)D1遠(yuǎn)，導(dǎo)致測(cè)點(diǎn)D2的位移波動(dòng)較??；隨著沖擊同一徑向桿的點(diǎn)數(shù)的增多，測(cè)點(diǎn)D1的位移逐漸增大。圖11(d)為第4組試驗(yàn)工況下測(cè)點(diǎn)最大位移變化曲線，由于測(cè)點(diǎn)D1位于沖擊點(diǎn)下，導(dǎo)致其位移明顯大于測(cè)點(diǎn)D2的位移；其中工況4-3的3個(gè)沖擊點(diǎn)位于測(cè)點(diǎn)D1的一側(cè)，工況4-4的3個(gè)沖擊點(diǎn)均勻分布在測(cè)點(diǎn)D1的兩側(cè)，較為集中，因而工況4-4測(cè)點(diǎn)D1的位移比工況4-3大。從圖11可以看出，有限元的模擬結(jié)果與試驗(yàn)結(jié)果稍有差異，這是由試驗(yàn)誤差引起的，但兩者的整體變化趨勢(shì)相同。

圖11 工況位移峰值Fig.11 Displacement peak of the first group

2.3 加速度分析

加速度與結(jié)構(gòu)的慣性力密切相關(guān)，是反映結(jié)構(gòu)振動(dòng)狀態(tài)重要因素，為了研究加速度對(duì)網(wǎng)殼結(jié)構(gòu)動(dòng)力響應(yīng)的影響，試驗(yàn)記錄了網(wǎng)殼結(jié)構(gòu)部分關(guān)鍵節(jié)點(diǎn)的加速度時(shí)程曲線。現(xiàn)對(duì)沖擊點(diǎn)位于同一主肋桿的工況3-1至工況3-4的加速度時(shí)程曲線進(jìn)行分析，其中A3加速度傳感器安裝于沖擊區(qū)，A1、A2位于非沖擊區(qū)。

圖12～15為各工況加速度傳感器的時(shí)程曲線圖。隨著沖擊點(diǎn)逐漸增多，測(cè)點(diǎn)加速度值逐漸增大。工況3-2和工況3-3同為2點(diǎn)沖擊，但工況3-2沖擊點(diǎn)較為集中，并且在試驗(yàn)過(guò)程中，由于沖擊球沒(méi)有同時(shí)下落到結(jié)構(gòu)上，導(dǎo)致A3加速度時(shí)程曲線出現(xiàn)2個(gè)明顯峰值。

圖12 工況3-1加速度時(shí)程曲線Fig.12 Acceleration time-history curve under condition 3-1

圖13 工況3-2加速度時(shí)程曲線Fig.13 Acceleration time-history curve under condition 3-2

圖14 工況3-3加速度時(shí)程曲線Fig.14 Acceleration time-history curve under condition 3-3

圖15 工況3-4加速度時(shí)程曲線Fig.15 Acceleration time-history curve under condition 3-4

其中工況3-3加速度時(shí)程曲線中A1、A2、A3的第1次加速度峰值為4 404.9、4 533.7、5 583.3 m/s2，其中由于A3位于沖擊區(qū)，A3加速度傳感器第1次峰值加速度出現(xiàn)在0.013 s，而A1、A2位于非沖擊區(qū)，導(dǎo)致A1、A2加速度傳感器的第1次峰值加速度相對(duì)于A1滯后0.003 s左右；工況3～4加速度時(shí)程曲線中A1、A2、A3的第1次加速度峰值為5 423.9、5 100.7、11 078.7 m/s2，A3加速度傳感器第1次峰值加速度出現(xiàn)在0.013 s，A1、A2加速度傳感器的第1次峰值加速度相對(duì)與A1滯后0.005 s左右。從工況3-3到工況3-4沖擊點(diǎn)由2個(gè)增加為3個(gè)，A1、A2、A3加速度峰值的增幅分別為23.13%、12.51%、98.40%，可以看出沖擊點(diǎn)的增加會(huì)導(dǎo)致沖擊區(qū)域的加速度發(fā)生較大增幅，非沖擊區(qū)加速度增幅較小；而加速度傳感器A1、A2之間的誤差主要是由試驗(yàn)中傳感器的安裝位置發(fā)生偏移以及受空氣阻力的影響導(dǎo)致沖擊位置偏離沖擊節(jié)點(diǎn)中心產(chǎn)生的。

2.4 試驗(yàn)現(xiàn)象與數(shù)值模擬

試驗(yàn)共4組，每組4個(gè)工況，均使用直徑為100 mm的鋼球從距離沖擊點(diǎn)4 m高處釋放，試驗(yàn)通過(guò)增加沖擊點(diǎn)個(gè)數(shù)和變換沖擊位置來(lái)研究K6型網(wǎng)殼結(jié)構(gòu)的動(dòng)力響應(yīng)。試驗(yàn)結(jié)果顯示，鋼球沖擊完成后，網(wǎng)殼結(jié)構(gòu)在沖擊點(diǎn)處有局部輕微凹陷，并且桿件與焊接球的連接部位有輕微變形，從而導(dǎo)致漆皮脫落現(xiàn)象，結(jié)構(gòu)整體處于彈性階段，無(wú)明顯變形，能夠繼續(xù)承載。以工況1-1為例，圖16中看出沖擊點(diǎn)的最終變形，有限元模擬結(jié)果與實(shí)驗(yàn)結(jié)果較為吻合。

圖16 工況1-1沖擊點(diǎn)變形對(duì)比Fig.16 Deformation comparison of impact point under condition 1-1

沖擊試驗(yàn)值和數(shù)值模擬結(jié)果值對(duì)比分析表明，測(cè)點(diǎn)響應(yīng)的基本規(guī)律一致，誤差在允許范圍之內(nèi)，驗(yàn)證了數(shù)值模擬的有效性，同時(shí)驗(yàn)證了網(wǎng)殼結(jié)構(gòu)在多點(diǎn)彈性沖擊荷載作用下的變形模式。以第1組工況為例，具體數(shù)據(jù)如表3所示。從整體來(lái)說(shuō),與模擬值相比，試驗(yàn)值偏小，其中加速度值誤差范圍為1.99%～12.45%，位移值誤差范圍為2.67%～16.85%。造成這些誤差的原因是：

表3 第1組工況試驗(yàn)值與模擬值Table 3 Test value and simulation value of the first working condition

1)有限元模型為理想模型，但實(shí)際沖擊試驗(yàn)時(shí)，網(wǎng)殼模型在加工廠制作時(shí)就存在一定的初始誤差。

2)有限元模擬沖擊時(shí)，沖擊球能精準(zhǔn)沖擊到網(wǎng)殼結(jié)構(gòu)的沖擊點(diǎn)上，但實(shí)際沖擊試驗(yàn)中，由電磁鐵吸附的鋼球在高空中存在微小擺動(dòng)，導(dǎo)致沖擊球下落到?jīng)_擊點(diǎn)時(shí)沖擊位置發(fā)生了偏差。

3 結(jié)論

1)隨著沖擊點(diǎn)數(shù)的增加，網(wǎng)殼沖擊區(qū)域的應(yīng)變、位移及加速度逐漸增加，網(wǎng)殼非沖擊區(qū)域的應(yīng)變、位移及加速度變化無(wú)明顯規(guī)律且影響較小。在相同的沖擊點(diǎn)數(shù)量及沖擊速度下，沖擊點(diǎn)位置對(duì)網(wǎng)殼結(jié)構(gòu)的動(dòng)力響應(yīng)影響較大，沖擊點(diǎn)越集中，網(wǎng)殼結(jié)構(gòu)沖擊區(qū)的應(yīng)變、位移及加速度越大；網(wǎng)殼結(jié)構(gòu)非沖擊區(qū)域的應(yīng)變、位移及加速度變化較小。

2)綜合比較第1組試驗(yàn)工況和第2組試驗(yàn)工況，當(dāng)沖擊速度一定時(shí)，沖擊點(diǎn)距網(wǎng)殼結(jié)構(gòu)頂點(diǎn)越遠(yuǎn)，網(wǎng)殼結(jié)構(gòu)沖擊區(qū)的位移越小，頂點(diǎn)位移也越小。綜合比較第3組試驗(yàn)工況與第4組試驗(yàn)工況，當(dāng)沖擊速度一定時(shí)，沖擊網(wǎng)殼結(jié)構(gòu)頂點(diǎn)時(shí)結(jié)構(gòu)沖擊區(qū)的位移較大。

3)沖擊試驗(yàn)和數(shù)值模擬的結(jié)果基本吻合，變化規(guī)律一致，驗(yàn)證了數(shù)值模擬方法的有效性，同時(shí)驗(yàn)證了網(wǎng)殼結(jié)構(gòu)在多點(diǎn)彈性沖擊下的變形模式：沖擊點(diǎn)處有局部輕微凹陷，沖擊區(qū)桿件與焊接球連接處有輕微變形。

4)在網(wǎng)殼結(jié)構(gòu)的抗沖擊設(shè)計(jì)中，為避免結(jié)構(gòu)遭受多點(diǎn)密集沖擊和結(jié)構(gòu)頂部遭受沖擊時(shí)結(jié)構(gòu)破壞，設(shè)計(jì)中應(yīng)對(duì)結(jié)構(gòu)頂部構(gòu)件進(jìn)行加強(qiáng)，并對(duì)沖擊荷載下易破壞的關(guān)鍵桿件與節(jié)點(diǎn)的連接處做加強(qiáng)處理。