軟土地鐵車(chē)站地震響應(yīng)數(shù)值計(jì)算方法的研究
楊超 楊林德 季倩倩
摘要:對(duì)飽和軟粘土采用粘彈塑性動(dòng)力本構(gòu)模型,利用拉格朗日差分法對(duì)典型軟土地鐵車(chē)站結(jié)構(gòu)建立地鐵車(chē)站地震響應(yīng)的數(shù)值計(jì)算方法。并進(jìn)一步利用該方法對(duì)軟土地鐵車(chē)站結(jié)構(gòu)振動(dòng)臺(tái)模型試驗(yàn)進(jìn)行數(shù)值擬合分析,結(jié)果表明土體和結(jié)構(gòu)模型的加速度響應(yīng)、結(jié)構(gòu)模型表面的動(dòng)土壓力以及結(jié)構(gòu)構(gòu)件的應(yīng)變規(guī)律的計(jì)算結(jié)果與試驗(yàn)結(jié)果基本吻合。
關(guān)鍵詞:軟土;地鐵車(chē)站結(jié)構(gòu);振動(dòng)臺(tái)試驗(yàn);數(shù)值計(jì)算方法 1 引言 神戶(hù)地震和歷史上發(fā)生的大震一再表明,對(duì)軟土地基中的地鐵車(chē)站等地下結(jié)構(gòu)的抗震設(shè)計(jì)開(kāi)展研究有重要的意義。對(duì)地下結(jié)構(gòu)地震響應(yīng)的計(jì)算,迄今已提出多種算法[1],然而由于對(duì)其涉及的各類(lèi)復(fù)雜因素的影響尚認(rèn)識(shí)不足,不同的計(jì)算方法或模型得出的結(jié)果存在很大的差異,且很難鑒別各自的合理性。本文建立軟土地鐵車(chē)站地震響應(yīng)的分析理論與計(jì)算方法,并通過(guò)對(duì)模型試驗(yàn)進(jìn)行擬合分析驗(yàn)證了所建立的車(chē)站結(jié)構(gòu)動(dòng)力響應(yīng)的計(jì)算方法的正確性和合理性,以便工程設(shè)計(jì)實(shí)踐參考。 2 軟土地鐵車(chē)站結(jié)構(gòu)的振動(dòng)臺(tái)試驗(yàn) 軟土地鐵車(chē)站結(jié)構(gòu)的振動(dòng)臺(tái)模型試驗(yàn)分自由場(chǎng)振動(dòng)臺(tái)模型試驗(yàn)、典型地鐵車(chē)站結(jié)構(gòu)振動(dòng)臺(tái)模型試驗(yàn)和地鐵車(chē)站接頭結(jié)構(gòu)振動(dòng)臺(tái)模型試驗(yàn)等三種。試驗(yàn)開(kāi)展過(guò)程中遇到的技術(shù)難題包括對(duì)地鐵車(chē)站縱向長(zhǎng)度的模擬,場(chǎng)地土的動(dòng)力特性與地震響應(yīng)的模擬,模型箱的構(gòu)造與邊界效應(yīng)的模擬,以及量測(cè)元件設(shè)置位置的優(yōu)選等。筆者對(duì)這些技術(shù)難題逐一進(jìn)行了研究,并都提出了行之有效的解決方法, 使試驗(yàn)取得了可靠的數(shù)據(jù)[2][3]。 試驗(yàn)過(guò)程中,首先進(jìn)行了自由場(chǎng)振動(dòng)臺(tái)模型試驗(yàn),用以模擬自由場(chǎng)地土層的地震反應(yīng),據(jù)以獲得模型箱內(nèi)不同位置處的土的加速度響應(yīng),確定“邊界效應(yīng)”的影響程度和鑒別模型箱構(gòu)造的合理性;然后通過(guò)典型地鐵車(chē)站結(jié)構(gòu)振動(dòng)臺(tái)模型試驗(yàn)了解地鐵車(chē)站結(jié)構(gòu)與土共同作用時(shí)地震動(dòng)反應(yīng)的規(guī)律與特征,為建立地鐵車(chē)站地震響應(yīng)的分析理論和計(jì)算方法提供試驗(yàn)數(shù)據(jù)。振動(dòng)臺(tái)模型試驗(yàn)記錄了在不同荷載級(jí)別的EI-Centro波、上海人工波和正弦波激振下,加速度測(cè)點(diǎn)傳感器的反應(yīng);由動(dòng)土壓力傳感器,得到了各測(cè)點(diǎn)在不同加載工況下的動(dòng)土壓力反應(yīng)時(shí)程;根據(jù)結(jié)構(gòu)模型構(gòu)件上布置的應(yīng)變片,測(cè)得了構(gòu)件應(yīng)變的變化。 3 軟土地鐵車(chē)站計(jì)算方法 將自由場(chǎng)土體簡(jiǎn)化為多自由度體系,其動(dòng)力平衡方程可表示為: [M]{ü}+[C]{u}+[K]{u}={f}(1) 式中[M]、[C]、[K]分別為體系的質(zhì)量矩陣,阻尼矩陣及剛度矩陣,{櫣}、{敶}{u}分別為相對(duì)加速度向量,相對(duì)速度向量和相對(duì)位移向量,{f}為荷載向量,對(duì)于非周期性的地震作用,初始時(shí)刻的結(jié)構(gòu)體系的速度和位移一般為零,求解式(1)可得結(jié)構(gòu)體系的瞬態(tài)反應(yīng)。 本文采用拉格朗日差分法對(duì)式(1)求解,特點(diǎn)為在時(shí)域內(nèi)將動(dòng)力平衡方程轉(zhuǎn)化為運(yùn)動(dòng)方程和應(yīng)力-應(yīng)變關(guān)系,即將計(jì)算區(qū)域離散為二維單元,單元之間由節(jié)點(diǎn)聯(lián)結(jié),并將運(yùn)動(dòng)方程:
采用如圖1的過(guò)程求解,直到不平衡力足夠小為止。
軟粘土在卸載再加載及反向加載的過(guò)程中同時(shí)伴有彈性和塑性變形[4];同時(shí)對(duì)軟粘土進(jìn)行的動(dòng)三軸試驗(yàn)表明,應(yīng)變趨向于零時(shí),其阻尼并不趨向于零,即在應(yīng)變趨于零時(shí),仍存在能量耗散[5]。筆者結(jié)合軟土地鐵車(chē)站振動(dòng)臺(tái)模型試驗(yàn)中對(duì)飽和軟粘土進(jìn)行的動(dòng)三軸試驗(yàn),將軟土在動(dòng)荷載作用下的能量耗散分為粘性和非粘性?xún)刹糠郑眠吔缑婺P屠碚摻④浾惩恋恼硰椝苄詣?dòng)力本構(gòu)模型,并通過(guò)對(duì)自由場(chǎng)振動(dòng)臺(tái)試驗(yàn)進(jìn)行數(shù)值模擬計(jì)算,驗(yàn)證了該模型的有效性[6]。本文中對(duì)軟粘土的動(dòng)力本構(gòu)模型采用粘彈塑型模型,將軟粘土在動(dòng)荷載作用下能量耗散分為粘性和非粘性?xún)刹糠郑渲姓承圆糠种慌c應(yīng)變率有關(guān)而與應(yīng)變的大小無(wú)關(guān);非粘性部分為塑性變形的加卸載過(guò)程中的能量耗散。 4 地鐵車(chē)站結(jié)構(gòu)的振動(dòng)臺(tái)試驗(yàn)的擬合分析 4.1 計(jì)算簡(jiǎn)圖 對(duì)地鐵車(chē)站結(jié)構(gòu)進(jìn)行的三維計(jì)算與分析表明,橫向激振條件下離端部較遠(yuǎn)的地鐵車(chē)站結(jié)構(gòu)可簡(jiǎn)化為平面應(yīng)變問(wèn)題進(jìn)行分析。本文擬對(duì)離端部較遠(yuǎn)的主觀測(cè)斷面按平面應(yīng)變問(wèn)題計(jì)算,方向與激振方向平行,并與車(chē)站結(jié)構(gòu)模型的縱軸垂直。計(jì)算區(qū)域以模型箱為界,底部邊界在豎直方向固定,側(cè)向邊界在水平方向固定,上表面為自由變形邊界。振動(dòng)過(guò)程中,模型箱發(fā)生的變形,可略去不計(jì),故側(cè)向和底部邊界在水平方向的加速度始終與臺(tái)面輸入波一致。計(jì)算網(wǎng)格劃分如圖2所示。 模型箱內(nèi)襯厚17.5cm的泡沫塑料板,用以模擬場(chǎng)地土易于變形的特性,劃分網(wǎng)格時(shí)泡沫塑料板和模型土均被離散為四邊形單元,車(chē)站結(jié)構(gòu)模型離散為梁?jiǎn)卧⒃谂菽芰习迮c土體、土體與車(chē)站結(jié)構(gòu)之間設(shè)置了接觸面單元,接觸面單元由法向彈簧、切向彈簧、抗拉元件和滑片組成,滑片剪切強(qiáng)度采用莫爾-庫(kù)侖準(zhǔn)則。 振動(dòng)臺(tái)模型試驗(yàn)中模型土的參數(shù)值示于表1。車(chē)站結(jié)構(gòu)材料的動(dòng)力特性參數(shù),擬按常規(guī)方法由將混凝土材料的靜彈性模量提高給出,研究表明動(dòng)彈性模量比靜彈性模量約高出3050%,微粒混凝土試樣的試驗(yàn)表明本次試驗(yàn)中微粒混凝土的靜彈性模量可取為Es=7.0GPa,則其動(dòng)彈性模量值為Es=7.0×1.4=9.8GPa。
4.2 計(jì)算結(jié)果與試驗(yàn)結(jié)果的擬合分析 自由場(chǎng)振動(dòng)臺(tái)模型試驗(yàn)表明,模型箱結(jié)構(gòu)合理,其邊界效應(yīng)的影響未波及到地鐵車(chē)站結(jié)構(gòu)模型所處的位置,鑒于典型地鐵車(chē)站結(jié)構(gòu)振動(dòng)臺(tái)模型試驗(yàn)中,用于接受激振響應(yīng)信息的傳感器有加速度傳感器、動(dòng)土壓力傳感器和應(yīng)變片等多種,以下擬對(duì)其分別作出擬合分析。 4.2.1 加速度反應(yīng)的擬合分析 (1)加速度反應(yīng)的放大系數(shù) 放大系數(shù)是指測(cè)點(diǎn)加速度反應(yīng)的峰值與振動(dòng)臺(tái)臺(tái)面輸入的峰值之比。地鐵車(chē)站結(jié)構(gòu)振動(dòng)臺(tái)模型試驗(yàn)中,土體一半厚度處測(cè)點(diǎn)和車(chē)站結(jié)構(gòu)模型上部測(cè)點(diǎn)的放大系數(shù)的計(jì)算結(jié)果、試驗(yàn)結(jié)果及相對(duì)誤差分別如表2和表3所示。由表可見(jiàn)各加載工況下土體與車(chē)站結(jié)構(gòu)模型加速度反應(yīng)放大系數(shù)的計(jì)算結(jié)果與試驗(yàn)結(jié)果均吻合較好,且上海人工波各工況的擬合程度更好。 (2)加速度反應(yīng)時(shí)程與富氏譜 對(duì)地鐵車(chē)站結(jié)構(gòu)振動(dòng)臺(tái)模型試驗(yàn),圖3、4給出了SH-4工況下土體一半厚度處測(cè)點(diǎn)的加速度反應(yīng)時(shí)程及其富氏譜的計(jì)算結(jié)果及相應(yīng)的試驗(yàn)結(jié)果,圖5、6給出了SH-4工況下車(chē)站結(jié)構(gòu)上部測(cè)點(diǎn)的加速度反應(yīng)時(shí)程及其富氏譜的計(jì)算結(jié)果與試驗(yàn)結(jié)果。由圖可見(jiàn)土體內(nèi)及結(jié)構(gòu)上測(cè)點(diǎn)的計(jì)算結(jié)果的波形、幅值與試驗(yàn)結(jié)果均基本吻合,兩者在各頻段的頻率組成也均基本吻合,表明文中的計(jì)算方法可較好地模擬地鐵車(chē)站結(jié)構(gòu)的地震加速度響應(yīng)。
4.2.2 車(chē)站結(jié)構(gòu)模型的動(dòng)土壓力的擬合分析 (1)動(dòng)土壓力的幅值 典型地鐵車(chē)站結(jié)構(gòu)模型試驗(yàn)中,側(cè)墻動(dòng)土壓力幅值的計(jì)算結(jié)果、試驗(yàn)結(jié)果及相對(duì)誤差如表4所列。由表4可見(jiàn)計(jì)算結(jié)果與實(shí)測(cè)結(jié)果基本吻合,且隨著輸入地震波荷載的增強(qiáng),兩者的相對(duì)誤差趨向增大,并在SH-10工況達(dá)到近20%,原因主要為隨著輸入荷載的增強(qiáng),土體的應(yīng)變?cè)龃螅狗蔷€(xiàn)性特征更加明顯,計(jì)算分析時(shí)產(chǎn)生的誤差逐漸增大。鑒于上海地區(qū)的地震設(shè)防烈度為7度,可認(rèn)為車(chē)站結(jié)構(gòu)模型側(cè)墻的動(dòng)土壓力的計(jì)算結(jié)果的幅值與實(shí)測(cè)結(jié)果基本吻合。
(2)動(dòng)土壓力的時(shí)程 圖7給出了在SH-4工況下,結(jié)構(gòu)模型側(cè)墻中部測(cè)點(diǎn)的動(dòng)土壓力時(shí)程的計(jì)算結(jié)果與試驗(yàn)結(jié)果。由圖可見(jiàn)地鐵車(chē)站結(jié)構(gòu)振動(dòng)臺(tái)模型試驗(yàn)中,結(jié)構(gòu)模型側(cè)墻不同部位測(cè)點(diǎn)的動(dòng)土壓力時(shí)程的計(jì)算結(jié)果的波形與試驗(yàn)結(jié)果基本吻合,也表明文中的計(jì)算方法可較好地模擬地鐵車(chē)站結(jié)構(gòu)與周?chē)馏w間的動(dòng)力相互作用。
4.2.3 車(chē)站結(jié)構(gòu)模型的動(dòng)應(yīng)變 車(chē)站結(jié)構(gòu)模型構(gòu)件的動(dòng)應(yīng)變幅值的實(shí)測(cè)結(jié)果表明,結(jié)構(gòu)構(gòu)件在各級(jí)荷載下均處于彈性受力狀態(tài)。鑒于下中柱下端的應(yīng)變最大,擬將各構(gòu)件的動(dòng)應(yīng)變與相同工況下下中柱下端的動(dòng)應(yīng)變相比較,并將比值稱(chēng)為構(gòu)件的相對(duì)應(yīng)變。計(jì)算結(jié)果和實(shí)測(cè)結(jié)果的相對(duì)應(yīng)變及其相對(duì)誤差如表5所示。由表可見(jiàn)車(chē)站結(jié)構(gòu)模型各構(gòu)件相對(duì)應(yīng)變的計(jì)算結(jié)果與實(shí)測(cè)結(jié)果基本吻合,本文采用的計(jì)算方法也可較好地模擬地鐵車(chē)站結(jié)構(gòu)的動(dòng)力變形特性。 5 結(jié)論 本文的軟土地鐵車(chē)站結(jié)構(gòu)的振動(dòng)臺(tái)模型試驗(yàn)為建立地鐵車(chē)站地震響應(yīng)的分析理論和計(jì)算方法提供了試驗(yàn)數(shù)據(jù)。采用本文的計(jì)算方法對(duì)振動(dòng)臺(tái)模型試驗(yàn)進(jìn)行擬合分析,結(jié)果表明該計(jì)算模型可較好地模擬軟土的動(dòng)力特性、地鐵車(chē)站結(jié)構(gòu)與土體的動(dòng)力相互作用,及地鐵車(chē)站結(jié)構(gòu)的動(dòng)力響應(yīng)特點(diǎn)。該數(shù)值計(jì)算方法較好地模擬了軟土地鐵車(chē)站的地震響應(yīng),可供工程設(shè)計(jì)實(shí)踐參考。