大連理工大學邵龍潭教授腦出血

Ⅰ 邵龍潭的學術著作

[1] 邵龍潭，李紅軍著，土工結構穩定分析—有限元極限平衡方法及其應用，科學出版社，2011.1.
[2] 邵龍潭著，土力學研究和探索，科學出版社，2011.3.
1. 土力學基本理論、應力應變性質及本構關系方面
[3] 邵龍潭，飽和土的土骨架應力方程，2011-2-20
[4] 邵龍潭，土力學的幾個基本問題，2011-2-20
[5] 邵龍潭，土體極限平衡條件和穩定安全系數定義，2011-2-20
[6] 邵龍潭，劉士乙，基於極限平衡條件的土體漸進破壞分析，2011-2-20
[7] 姚濤，邵龍潭. 基於圖像測量技術的非飽和黃土固結變形試驗研究. 岩土力學，2009，30(7)：2011-2015.
[8] 熊保林，邵龍潭. 考慮中主應力效應Gudehus-Bauer亞塑性模型的改進. 哈爾濱工業大學學報，2009，41(6)：127-130.
[9] 邵龍潭，熊保林，王洪波. Wu-Bauer亞塑性模型參數確定方法的改進及其驗證. 岩土力學，2008，29(3)：707-710.
[10] 熊保林，邵龍潭. 考慮主應力軸旋轉Gudehus-Bauer亞塑性模型的改進. 工程力學, 2008，25(1)：127-132.
[11] 熊保林，邵龍潭. 土亞塑性理論研究進展. 第十屆土力學及岩土工程學術會議，重慶, 2007：582-587.
[12] 董建軍，邵龍潭. 非飽和土三軸壓縮應力路徑試驗的數值模擬研究. 岩土力學, 2006，27(S1)：95-98.
[13] 王洪波，邵龍潭，熊保林. 確定亞塑性模型參數n、hs的一種改進方法. 岩土工程學報，2006，28(9)：1173-1176.
[14] 孫益振，邵龍潭. 基於局部與整體變形測量的粉土泊松比試驗研究. 岩土工程學報，2006，28(8)：1033-1038.
[15] 王洪波，邵龍潭，張學增. 基於亞塑性理論的無粘性土壓縮試驗應力應變的研究. 岩土工程學報, 2006，28(6)：780-783.
[16] Xiong B.L, Shao L.T. Strain Path Research of a Hypo-plasticConstitutive Model for Cohesion-less Soil., Geo-Shanghai Internationalconference Shanghai China, 2006 June6-8 298-303.
[17] 孫益振，邵龍潭，李根華. 基於數字圖象測量的Duncan-Chang模型參數研究. 岩土力學，2006，27(8)：1335-1338.
[18] 邵龍潭,孫益振.考慮孔隙變形的孔隙介質本構關系初探.岩土力學，2006，
27(4)：561-565.
[19] 邵龍潭，洪帥，鄭衛鋒. 循環孔隙水壓力作用下飽和砂土變形的試驗研究. 岩土工程學報，2006，28(4)：428-431.
[20] 孫益振，邵龍潭. 三軸循環加卸載條件下砂性土變形特性研究. 岩土工程學報，2005，27(11)：1353-1357.
2、土壤滲流方面
[21] 梁愛民，邵龍潭.土壤中空氣對土扒滑結構和入滲過程的影響. 水科學進展，2009，20(4)：春滾臘502-506.
[22] 馬國濤，邵龍潭，龍飛，馮恩民. 水流入滲和氣體驅替過程二維問題的耦合求解. 水科學進展，2002，13(6)：741-746.
[23] 邵龍潭. 相間相互作用原理與土壤水動力學基本方程.水科學進展, 2002，
13(5)：78-82.
[24] 邵龍潭，王助貧, 關立軍，許志強. 非飽和土中水流入滲和氣體排出過程的求解. 水科學進展，2000，11(1)：8-13.
3、試驗測試技術及儀器方面
[25] 姚濤，邵龍潭. 三軸數字圖像測量技術在黃土力學特性研究中的應用.大連理工大學學報，2009，49(1)：92-97.
[26] 桑勇，邵龍潭. 土工試驗中土樣動態變形測試系統的設計. 液壓與氣動，2008，(7)：23-25.
[27] 邵龍潭，董建軍，劉永祿，孫益振. 基於亞像素角點檢測的試樣變形圖像測量方法備純. 岩土力學，2008，29(5)：1329-1333.
[28] 董建軍，邵龍潭，劉永祿，姚濤. 基於圖像測量方法的非飽和壓實土三軸試樣變形測量. 岩土力學，2008，29(6)：1618-1622.
[29] 孫益振，邵龍潭. 三軸試樣變形數字圖像測量系統影響因素分析. 第十屆土力學及岩土工程學術會議，重慶，2007：543-548.
[30] 董建軍，邵龍潭，劉永祿. 基於圖像測量的非飽和土p-s平面屈服軌跡研究. 水利水運工程學報，2007，9 (3)：6-11.
[31] 王助貧，邵龍潭，孫益振. 基於數字圖像測量技術的粉煤灰三軸試樣剪切帶研究. 岩土工程學報，2006，28(9)：1163-1167.
[32] 邵龍潭，孫益振，王助貧，劉永祿. 數字圖像測量技術在土工三軸試驗中的應用研究. 岩土力學，2006，27(1)：29-34.
[33] 劉永祿，邵龍潭，潘石. 三軸實驗土樣變形數字圖像測量（DDIM）的一種改進演算法. 第十一屆全國實驗力學會議，大連, 2005：1107-1112.
[34] 邵龍潭，梁愛民，王助貧，孫健. 非飽和土穩態滲流試驗裝置的研製與應用. 岩土工程學報，2005，27(11)：1338-1340.
[35] 王助貧，邵龍潭. 三軸試驗土樣的端部影響問題研究. 岩土力學，2003，24(3)：363-368.
[36] 鄭衛鋒,洪帥,邵龍潭. 精密水壓力發生器的研製與應用.實驗力學，2003，
18(3)：389-392.
[37] 邵龍潭，王助貧. 三軸土樣局部變形的數字圖像測量方法. 岩土工程學，
2002，24(2)：159-163.
[38] 邵龍潭，王助貧，劉永祿. 三軸試樣變形數字圖像測量的誤差和精度分析. 大連理工大學學報，2002，42(1)：98-103.
[39] 邵龍潭，王助貧，韓國城，燕宏實. 三軸試驗土樣徑向變形的計算機圖像測量. 岩土工程學報，2001，23(3)：337-341.
4、穩定分析方面
[40] 趙傑，邵龍潭. 深基坑土釘支護的有限元數值模擬及穩定性分析. 岩土力學，2008，29(4)：983-988.
[41] 趙傑，邵龍潭. 基坑土釘支護邊坡有限元穩定性分析方法探討. 岩土力，
2008，29(6)：1654-1658.
[42] 邵龍潭. 邊坡穩定性分析極限平衡法的討論. 第十屆土力學及岩土工程學術會議，重慶, 2007：292-297.
[43] 趙傑，邵龍潭. 土體結構極限承載力的有限元分析. 岩石力學與工程學，
2007，26(s1)：3183-3189.
[44] 趙傑，邵龍潭. 平面應變條件下兩類有限元邊坡穩定分析方法比較研究. 大連理工大學學報， 2007，47(6)：873-879.
[45] 趙傑，邵龍潭. 填築和開挖邊坡的穩定性分析. 岩土力學，2007，28(5)：944-950.
[46] 張士林，邵龍潭. 非飽和土坡的穩定分析，西部探礦工程，2006，(12)：281-283.
[47] ZhaoJie, Shao L.T. Finite Element Method for Stability Analysis with Its
Applicationto Limit Load of Soil Mass. International conference
ShanghaiChina,2006.June 6-8, 336-342.
[48] 趙傑，邵龍潭. 有限元穩定分析法在確定土體結構極限承載力中的應用. 水利學報，2006，37(6)： 668-673.
[49] 唐洪祥，邵龍潭. 地震動力作用下有限元土石壩邊坡穩定性分析. 岩石力學與工程學報，2004，23(8)：1318-1324.
[50] 邵龍潭，唐洪祥，韓國城. 有限元邊坡穩定分析方法及其應用. 計算力學學報，2001，18(1)：81-87.
[51] L.T.Shao and Z. P. Wang. On the Stability of Unsaturated SoilSlopes. Asian Conference on Unsaturated Soils, Sigapore，2000：825-829.
[52] 邵龍潭，唐洪祥，孔憲京，韓國城. 隨機地震作用下土石壩邊坡的穩定性分析. 水利學報，1999，11：65-71.
[53] 邵龍潭，韓國城. 水流作用下堆石邊坡的穩定分析方法. 水利學報，1997，1：51-55.
[54] 楊蘊明，邵龍潭.土工加筋結構的有限元邊坡穩定分析.地基處理，1996，7(4)：1-7.
[55] 邵龍潭，韓國城.堆石邊坡穩定分析的一種方法.大連理工大學學報，1994，31(3)：365-369.

Ⅱ 大連理工大學博導名單

大連理工大帶脊學博導名單：

1、邵龍潭：1963年7月生，博士，二級教授，博士生導師，國務院特殊津貼專家。2005-2013年先後擔任大連理工大學黨委副書記、紀委書記、副校長。現為國際土力學學會會員，中國力學學會實驗力學專業委員會委員，非飽和土與特殊土專業委員會常委。《大連理工大學學報》、《實驗力學》等期刊編委，國家自然科學基金項目評審專家。

長期從事孔隙材料的實驗力學特性、土力學基本理論、土工測試技術與試驗設備研發、土工結構穩定分析等方面的研究工作。主持參加國家自然科學基金項目、國家科技攻關項目、企業委託項目共70餘項。發表學術論文150餘篇，申請授權國家發明專利37項，出版學術著作5部，教材1部。獲國家水利電力部科技進步二等獎1項，海洋工程科學技術獎特等獎1項。

Ⅲ 李俊傑的出版著作和論文

部分代表性論文如下： [1] Fei Kang, Junjie Li. Artificial bee colony algorithm optimized support vector regression for system reliability analysis of slopes. Journal of Computing in Civil Engineering, ASCE, 2015, Accepted. (SCI&EI)
[2] Fei Kang, Shaoxuan Han, Rodrigo Salgado, Junjie Li. System probabilistic stability analysis of soil slopes using Gaussian process regression with Latin hypercube sampling. Computers and Geotechnics
[3] Haojin Li, Junjie Li, Fei Kang. Application of the artificial bee colony algorithm-based projection pursuit method in statistical rock mass stability estimation. Environmental Earth Sciences
[4] Fei Kang, Junjie Li, Haojin Li, Artificial bee colony algorithm and pattern search hybridized for global optimization, Applied Soft Computing Top 25 Hottest Articles
[5] Fei Kang, Junjie Li, Zhenyue Ma. An artificial bee colony algorithm for locating the critical slip surface in slope stability analysis. Engineering Optimization
[6] Fei Kang, Junjie Li, Qing Xu. Damage detection based on improved particle swarm optimization using vibration data. Applied Soft Computing
[7] Fei Kang, Junjie Li, Zhenyue Ma. Rosenbrock artificial bee colony algorithm for accurate global optimization of numerical functions. Information Sciences,(SCI&EI) Most Cited Articles since 2010; Top 25 Hottest Articles
[8] Haojin Li, Junjie Li, Fei Kang. Risk analysis of dam based on artificial bee colony algorithm with fuzzy c-means clustering. Canadian Journal of Civil Engineering
[9] Zhou Hui, Li Jun-jie, Kang Fei. Distribution of acceleration and empirical formula for calculating maximum acceleration of rockfill dams. Journal of Central South University of Technology
[10] Fei Kang, Junjie Li, Qing Xu. Structural inverse analysis by hybrid simplex artificial bee colony algorithms. Computers & Structures, (SCI&EI) Most Cited Articles since
[11] Fei Kang, Junjie Li, Qing Xu. Virus coevolution partheno-genetic algorithms for optimal sensor placement. Advanced Engineering Informatics
[12] Wei Zeng, Junjie Li, and Fei Kang, Numerical Manifold Method with Endochronic Theory for Elastoplasticity Analysis, Mathematical Problems in Engineering
[13] Fei Kang, Junjie Li, Sheng Liu. Combined data with particle swarm optimization for structural damage detection. Mathematical Problems in Engineering, Volume
[14] Xu Wang, Fei Kang, Junjie Li, Xin Wang. Inverse parametric analysis of seismic permanent deformation for earth-rockfill dams using artificial neural networks. Mathematical Problems in Engineering, Volume [15] 仝宗良, 曾偉, 李俊傑. 基於數值流形法的土質邊坡動力穩定性分析. 岩土工程學報
[16] 汪旭, 康飛, 李俊傑. 土石壩地震永久變形參數反演方法研究. 岩土力學
[17] 李浩瑾, 李俊傑, 康飛. 基於 LSSVM 的重力壩地震穩定易損性分析. 振動與沖擊
[18] 李浩瑾, 李俊傑, 康飛, 張勇. 重力壩縱縫非連續接觸的地震反應分析. 大連理工大學學報
[19] 周暉, 李俊傑, 康飛.面板堆石壩壩頂加速度沿壩軸線分布規律.岩土力學
[20] 康飛, 李俊傑, 許青. 混合蜂群演算法及其在混凝土壩動力參數反演中的應用. 水利學報
[21] 康飛, 李俊傑, 許青. 堆石壩材料參數反演的蟻群聚類RBF網路模型. 岩石力學與工程學報
[22] 宋志宇，李俊傑，汪宏宇. 混沌人工魚群演算法在重力壩材料參數反演中的應用. 岩土力學
[23] 宋志宇，李俊傑. 基於微粒群演算法的大壩材料參數反分析研究. 岩土力學
[24] 胡軍, 李俊傑, 劉德志. 考慮剪切抗力的修正土釘單元及其應用. 岩土力學
[25] 劉德志, 李俊傑. 大壩安全監測資料的非線性檢驗, 應用基礎與工程科學學報
[26] 劉德志, 李俊傑. 土石壩安全監測軟體系統設計與實現，大連理工大學學報
[27] 楊清平, 李俊傑. 重力壩壩踵主拉應力區分布規律的探討，水利學報
[28] 李俊傑，邵龍潭，邵宇. 面板堆石壩永久變形研究，大連理工大學學報
[29] 李俊傑, 馬恆春. 蓄水期面板堆石壩動力特性研究. 岩土工程學報
[30] 李俊傑, 韓國城, 林皋. 混凝土面板堆石壩自振周期簡化公式. 振動工程學報
[31] 李俊傑, 韓國城, 孔憲京. 關門山面板堆石壩三維地震反應分析. 水利學報
[32] 李俊傑, 孔憲京, 韓國城. 面板堆石壩動力破壞計算方法研究. 大連理工大學學報
[33] 李俊傑, 韓國城, 林皋. 混凝土面板堆石壩地震加速度反應規律的幾點研究. 水利學報
[34] 孔憲京, 韓國城, 李俊傑. 關門山面板堆石壩二維地震反應分析. 大連理工大學學報
[35] 韓國城, 孔憲京, 李俊傑. 面板堆石壩動力破壞性態及抗震措施試驗研究. 水利學報
[36] 孔憲京，韓國城，李俊傑，林皋. 防滲面板對堆石壩體自振特性的影響，大連理工大學學報 [37] 楊春雨, 李俊傑. 改進的 SSOR-PCG 快速求解法在高面板堆石壩求解效率和節約內存中的實踐. 水電能源科學,
[38] 曾偉, 李俊傑. 基於 NMM-DDA 的直剪試驗數值模擬. 水電能源科學
[39] 劉景, 李俊傑. 不同開度時溢流壩弧形閘門水流三維數值模擬. 水電能源科學
[40] 曾偉, 李俊傑. 基於數值流形法的土石壩靜力計算數值模擬. 水利水電技術
[41] 宋宜祥, 李俊傑, 康飛. 虹吸井對尾礦壩地震液化的影響分析. 水電能源科學
[42] 康飛, 李俊傑, 馬震岳. 基於人工蜂群演算法的邊坡最危險滑動面搜索. 防災減災工程學報
[43] 楊秀萍, 李俊傑, 康飛. 基於 ACC-RBF 的水布埡面板堆石壩參數反演分析. 水電自動化與大壩監測
[44] 李浩瑾, 李俊傑, 康飛. 基於 ABCA-LSSVM 的復雜工程結構可靠度計算. 水電能源科學
[45] 康飛, 李俊傑, 馬震岳. 邊坡穩定分析的差分進化全局求解. 水電能源科學
[46] 杜文才, 李俊傑. 貯灰壩安全預警模型研究. 水電能源科學
[47] 李浩瑾, 李俊傑, 康飛. 基於 PSO-AHP 的大壩致災因子權重計算. 防災減災工程學報
[48] 胡崢嶸, 李俊傑. 面板堆石壩三維非線性有限元並行計算. 力學與實踐
[49] 康飛，李俊傑，許青，張運花. 改進人工蜂群演算法及其在反演分析中的應用. 水電能源科學
[50] 張運花，李俊傑，康飛. 西龍池面板堆石壩應力變形三維有限元分析. 水電能源科學，第
[51] 宋志宇，李俊傑. 基於模擬退火神經網路模型的岩質邊坡穩定性評價方法. 長江科學院院報
[52] 李俊傑，胡軍，康飛，王誼. 大頂子山溢流壩長閘墩溫度應力模擬計算分析. 水電能源科學
[53] 康飛，馬妹英，李俊傑. 支持向量回歸在貯灰壩滲流監測中的應用. 水電自動化與大壩監測
[54] 宋志宇，李俊傑. 最小二乘支持向量機在大壩變形預測中的應用. 水電能源科學
[55] 張振國，李俊傑，楊曉明. 基於變分原理的三維土坡穩定分析方法研究及應用. 水電能源科學
[56] 劉德志,李俊傑,許青. 基於Internet-Intranet的火電廠貯灰壩自動化安全監測系統. 水電能源科學
[57] 李俊傑，馬妹英，許青. RBF網路在貯灰壩浸潤線預測中的應用. 水電能源科學
[58] 譚志軍， 李俊傑. 混合遺傳演算法在貯灰壩監測系統上的應用. 水電能源科學
[59] 譚志軍，李俊傑. BP 演算法在貯灰壩監測系統中的應用. 水電自動化與大壩監測
[60] 費璟昊，李俊傑，李輝，楊建林. 利用圖像處理實現隧洞測量. 測繪通報
[61] 李俊傑, 李黎, 許勁松等. 中遠船塢抽水工程監測成果分析. 港口工程

Ⅳ 我今年大三，想考大連理工的工程力學研究生，但是不知道初試和復試的科目，還望哪位前輩指點一下，不勝感

大工研究生院網站裡面是不是沒有啊？我找了半天沒有看到，就看到一個「大連理工大學有碩士學位授權學科、專業一覽表（學術型）」。http://gs1.dlut.e.cn/newVersion/Front/sszs/xsxylb.html
一般研究生院都會有相關信息的，可能是沒找到地方吧。。。
中國考研網上面倒是有些信息，你參考下吧：
招生目錄返回 大連理工大學 首頁
大連理工大學 工程力學 2010年 招生目錄（一般情況是不會有大變動的。。）
專業代碼：080104
研究方向 01 工程結構多學科優化與反問題
02 結構與多學科耦合系統模擬軟體與應用
03 先進材料與結構性能表徵與現代設計理論
04 工程結構動力分析與控制
05 工程結構可靠性分析與健康診斷
初試科目 ① 101思想政治理論
② 201英語一
③ 301數學一
④ 816材料力學
復試科目 考試形式：面試加筆試
專業綜合筆試內容：包括工程力學基礎知識、動力學及數學
參考書目 816材料力學
《材料力學》，編者：劉鴻文，高等教育出版社，第三版上下冊
復試參考書目：
《高等數學》，四川大學數學系編，高等教育出版社
《常微分方程講義》，編者：葉彥謙，人民教育出版社，第二版，其他教材也適用
《理論力學》，哈爾濱工業大學編，高等教育出版社，第七版
《材料力學》，編者：孫訓方、方孝淑、關來泰，高等教育出版社
導師名單 該碩士點指導教師：程耿東**、歐進萍**、陳浩然*、郭杏林*、郭旭*、李剛*、李曉傑*、林家浩*、劉書田*、劉迎曦*、邵龍潭*、王希誠*、吳承偉*、岳前進*、胡平*、張洪武*、關振群*、亢戰*、王躍方*、季順迎*、周霞*、張文首、武湛君、白瑞祥、畢祥軍、陳金濤、董守華、李雲鵬、李守巨、王平、奚進一、閆鴻浩、楊迪雄、趙國忠、張小鵬、王博、孫國、張昭、閻軍、阮詩倫、馬國軍、任明法、張永存、王小紅、高強、梅玉林、王健
註：標「**」者為中科院院士或工程院院士；標「*」者為博士生導師
同等學力 同等學力考生復試期間必須加試兩門報考專業本科主幹課程。
嘻嘻~~我今年也大三，也准備考研啦~~加油吧~~~~

Ⅳ 邵龍潭的介紹

邵龍潭，1963年生，吉林省睜核梨樹人。1984年畢業於大連理工大學土木工程系。1987年獲碩士學位，1996年獲博士學位。1998-1999年赴德國留學。1987年始在大連姿瞎理工大學任教，1998年任教授，2000年受聘為博士生導師。跡早空現任大連理工大學校長助理。主持或參加完成國家自然科學基金項目、科技攻關項目以及企業委託項目多項。發表學術論文30餘篇。曾獲水利水電規劃設計總院科技進步獎一項，國家水利電力部科技進步二等獎一項。申請國家發明專利3項。

Ⅵ 地震動力穩定性資料

摘 要 本文考慮地震動的隨機性，在土石壩隨機地震反應分析和有限元邊坡穩定分析方法的基礎上，建立了隨機地震作用下土石壩邊坡的穩定性分析方法，並通過對土壩動力模型試驗的數值驗證及—理想土石壩邊坡的動力穩定性分析，證明這種方法是合理的、有效的。

關鍵詞 隨機地震反應，有限元，邊坡穩定分析，土石壩。

本文於1998年10月13日收到，系國家自然科學基金資助項目。

地震作用下邊坡的穩定性問題一直沒有得到很好地解決，以往慣用的極限平衡法及擬靜力法分析邊坡的地震動力穩定性存在著不少缺陷。本文在考慮輸入地震動荷載的平穩隨機特性進行壩體隨機動力反應分析[1、2、3]的基礎上，應用有限元邊坡穩定分析方法[4、5]，分析了邊坡的地震動力穩定性。文中對如何考慮隨機動應力作了處理，並對其合理性進行了論證。通過對模型壩和高土石壩兩個算例的計算分析，可以看出本文的方法還是很成功的。

1 分析方法簡述

1.1 隨機地震反應分析 本文將地震過程看作零均值平穩高斯過程。由隨機振動理論可知，對於高阻尼體系在平穩運動激勵下的初始非平穩響應段很短，可近似忽略，而按平穩響應處理。土工建築物可以當作高阻尼體系考慮，因此可以按平穩輸入平穩響應來進行分析。

在隨機荷載作用下，決定土層反應的一般二維等價線性方程為

(1)

其中{JX}、{JY}為水平與豎向荷載指示向量，、為水平與豎向地震加速度輸入過程，阻尼陣[ceq]按單元變阻尼法形成。在頻域上示解需對一個個的頻率離散空洞點分別進行，求解上式時可輸入加速度功率譜、在一系列ωj處離散。對第j個離散點，假定系統受幅值為、的虛擬簡諧運動激勵，這時問題的求解式(1)變為下述確定性線性方程

(2)

然後用振型分解法降價可迅速求解得到位移反應幅值，此即平穩隨機響應的確定性演算法，該法計算簡便且精度較高[9]。將穩態反應的位移幅值作為結構的靜變位，計算出各單元的正應變與剪應變幅值，由平面應變狀態下的應力-應變關系

(3)

即可求得每一單元的動應如段力幅值，幅值的平方即得功率值。對每一頻率渣虧譽離散點進行上述計算即得位移、應力反應的功率譜。在得到動應力反應的功率譜後，通過積分可得反應量的方差；應用直接插值等價線性化法[10]，可得等價的動應力平均幅值。同時，從此可求出最大動應力反應的中值(也即平均最大值)[1]。可以證明，按上面計算的穩態反應功率譜已計入了各階振型互相耦合的影響，結果是比較准確的[9]。

1.2 隨機地震作用下邊坡的穩定性分析 應用有限元邊坡穩定分析方法[4，5]，取土體的抗剪強度為莫爾-庫侖強度准則，那麼曲面上任一點土體的抗剪強度為

τf=σntgφ+c
(4)

式中σn為法向應力，φ和c分別是土體的內摩擦角和粘滯力。邊坡穩定分析的目的是要在計算區域內找到這樣一個曲面(平面問題為一條曲線)，沿這個面的抗滑穩定安全系數為最小。用有限元方法計算出壩體區域的應力場，並將平面問題的曲線離散後，問題的求解可以表示為

(5)

e為離散後曲線上的一個單元。上式可進一步寫成

(6)

式中|J|為雅可比行列式。上式可以用高斯數值積分計算。在靜力條件下，一點沿曲線方向的法向應力和切向應力用下式計算

(7)

式中

(8)

其中y n'是沿曲線方向的法線斜率。

在隨機地震作用下，式(7)中的各應力分量應為靜應力與隨機動應力分量之和。由於動應力是由隨機動力反應分析得到的，只能得到動應力的平均幅值與平均最大值的大小，而方向是不確定的，所以不能簡單地迭加上隨機動應力後進行最危險滑裂面的搜索。

如前所述，本文所考慮的邊坡穩定分析方法是一種在有限元應力分析基礎之上的、假定初始滑裂面、採用虎克 捷夫(Hooke-Jeeves)方法逐點、逐步搜索求解的數學規劃方法。為了迭加上隨機動應力又不至於增加太多的計算量，在每一計算點考慮3個動應力(σdx、σdy、σdxy)的隨機組合，則法向應力和切向應力可表示為

(9)

其中m=±1,n=±1,l=±1,它們的取值實際上代表了動應力的方向。當m、n、l分別取值時，在每一高斯點形成8種不同的應力組合。由每組算得的σn、τ代入目標函數，取其中對目標函數值貢獻最小的應力組合作為此點的計算應力。

從數學上講，用上述方法最終都能搜索得到最危險滑裂面並求得最小安全系數。用於最後計算對目標函數值貢獻的是其中最不利組合的一種，那麼，此時的動應力方向能否代表實際邊坡在地震作用下破壞時真實的動應力方向是一個要考慮的問題。因為這里所考慮的地震作用是隨機的，對某一點來說出現這樣那樣的應力方向是可能的，但以往大量的確定性地震反應分析表明，邊坡破壞時沿破壞面的應力分布具有一定規律性，這里考慮地震動的隨機性來分析邊坡的穩定性也應符合這種規律。在後面對模型壩的計算分析中將進一步對這一問題加以說明。

1.3 計算方法 為了分析隨機地震作用下土石壩邊坡的穩定性，需進行如下計算：(1)用隨機地震反應分析方法計算壩體的隨機動力反應，求出動應力場的平均幅值及平均最大值。(2)分別考慮隨機動應力的平均幅值及平均最大值按前述方法進行搜索求解，分別得到這兩種情況下的最小安全系數及最危險滑裂面。(3)為了用數學規劃法(Hook-Jeeves搜索法)搜索得到最危險滑裂面，先給定多條初始滑裂面進行搜索，找到各自的最小安全系數及其對應的滑裂面，將各安全系數進行比較，取其中最小的安全系數及其對應的滑裂面為問題的解。

2 應用分析方法的合理性

應用本文所述的隨機地震作用下邊坡的穩定性分析方法，可以求得平均意義上的最危險滑裂面及最小安全系數。而計算得到的滑裂面正確與否以及最後確定的動應力是否與實際相符還需要得到驗證。

由於這里所考慮的是隨機動力反應，那麼輸入的動力過程也應該是隨機的激勵荷載。較為合理的做法應該是在大量的隨機荷載激勵下，進行土石壩邊坡的破壞試驗，然後在統計出的輸入荷載數據的基礎上，進行壩體的隨機動力反應分析及壩坡的穩定性分析，再與統計的試驗結果相比較。但由於缺少這方面的試驗資料，進行試驗又有一定的困難，在目前情況下還難以做到，這里引用了董軍在日本東京大學完成的在正弦波激勵下模型砂壩破壞試驗的研究成果[6]，砂壩含水量為1.6%，粘聚力c=300Pa，內摩擦角φ=36°，壩的幾何尺寸及其破壞曲線示於圖2.在輸入假定功率譜的情況下，對它進行隨機動力反應分析及穩定性分析，以便能得到對隨機動力作用下邊坡穩定性分析方法合理性的一個印證。

為了盡可能與原試驗有一定程度的近似，這里按以下兩方面的要求選擇輸入的加速度功率譜：(1)由於原模型壩試驗輸入的是正弦波激勵，只有一個頻率分量，故選擇輸入的功率譜為一窄帶過程，且其主振頻率為正弦波的頻率，即5Hz;(2)所輸入加速度功率譜的總功率與輸入正弦波的總功率一致。對輸入的正弦波過程先進行FFT變換，再進行功率譜積分，即可得到輸入正弦波加速度的總功率值，以此作為輸入加速度功率譜的總功率值。
圖1 模型壩試驗的輸入功率譜曲線

符合以上條件的功率譜是不難找到的。這里選取了圖1所示的加速度功率譜作為輸入，用上述方法進行隨機動力反應分析，分別取動應力的平均幅值與平均最大值進行穩定性分析，將試驗結果與計算結果示於圖2.由圖可見，用本文的方法計算出的最危險滑裂面趨勢與試驗結果是比較一致的。同時，可以看出，不論是用平均幅值還是用平均最大值計算出來的最危險滑裂面位置都比靜力狀態下的要淺，試驗破壞面介於二者之間且更接近用平均幅值計算出來的最危險滑裂面。
TFL為試驗破壞曲線；SFL為靜力狀態下最危險滑裂面，K=2.25；DFL為動應力取平均幅值計算出的最危險滑裂面，K=1.14；DML為動應力取平均最大值計算出的最危險滑裂面，K=0.41.
圖2 模型壩的靜動力最危險滑裂面

同時，這里將模型壩破壞時的實際主應力分布示於圖3，將由搜索所確定的最危險滑裂面上各點的動應力示於表1，主應力沿滑裂面的分布示於圖4.可以看出，動應力方向並沒有因為考慮隨機組合而出現雜亂無章的情況，而是很有規律；接近邊坡左表面的主應力分布與實際破壞時情況也相一致。

3 隨機地震作用下高土石壩邊坡的穩定性分析

取一堆石壩的壩高為100m，壩頂寬10m，壩體上下游邊坡坡比為1∶1.4，堆石壩壩體為均質堆石材料，容重γ=2.0t/m3，靜力計算時壩體材料的應力-應變關系模型採用修改的鄧肯非線性雙曲線E—B模型。泊松比μ=0.25，最大動剪切模量Gmax=69.9(K2)max(σ0)0.5，其中σ0為平均有效應力，(K2)max=150,動摩擦角φ=42°。

表1 模型壩沿最危險滑裂面各點的動應力

--------------------------------------------------------------------------------

序號 沿滑裂線各點坐標 動應力(平均值)分量 序號 沿滑裂線各點坐標 動應力(平均值)分量

--------------------------------------------------------------------------------

--------------------------------------------------------------------------------

--------------------------------------------------------------------------------

--------------------------------------------------------------------------------

X/m Y/m σd/104Pa σdy/104Pa τxy/104Pa X/m Y/m σd/104Pa σdy/104Pa τxy/104Pa

--------------------------------------------------------------------------------

1 0.08 0.04 -0.01 -0.01 0.01 16 0.47 0.15 -0.11 -0.11 0.05
2 0.11 0.05 -0.01 -0.01 0.01 17 0.49 0.16 -0.11 -0.11 0.05
3 0.13 0.05 -0.02 -0.02 0.01 18 0.52 0.17 -0.11 -0.11 0.05
4 0.16 0.06 -0.03 -0.03 0.02 19 0.54 0.19 -0.11 -0.11 0.05
5 0.19 0.06 -0.03 -0.04 0.02 20 0.56 0.20 -0.10 -0.11 0.04
6 0.21 0.07 -0.04 -0.04 0.02 21 0.58 0.22 -0.09 -0.09 0.04
7 0.24 0.07 -0.05 -0.05 0.02 22 0.60 0.23 -0.08 -0.08 0.04
8 0.27 0.08 -0.06 -0.06 0.03 23 0.62 0.25 -0.07 -0.07 0.04
9 0.29 0.08 -0.07 -0.07 0.03 24 0.64 0.26 -0.05 -0.05 0.03
10 0.32 0.09 -0.07 -0.08 0.03 25 0.66 0.28 -0.04 -0.04 0.03
11 0.34 0.10 -0.08 -0.08 0.04 26 0.69 0.30 -0.03 0.03 0.02
12 0.37 0.10 -0.08 -0.08 0.03 27 0.71 0.33 -0.03 0.03 0.01
13 0.39 0.11 -0.08 -0.09 0.04 28 0.73 0.35 -0.02 0.02 0.01
14 0.42 0.12 -0.09 -0.10 0.05 29 0.75 0.36 -0.03 0.03 0.01
15 0.44 0.13 -0.10 -0.10 0.05

--------------------------------------------------------------------------------

(基底輸入加速度0.5g)
圖3 模型壩破壞時的主應力分布

圖4 模型壩取平均動應力幅值搜索得到的
最危險滑裂面上主應力分布

波選用了塔夫脫波的水平向分量、唐山波的水平向分量和豎直分量，而塔夫脫波又按最大加速度的不同考慮了幾種情況。隨機地震反應分析需要輸入加速度功率譜，對於已經記錄到的某加速度波形，可以將其看成是一平穩隨機地震過程的一個樣本的實現，按照Vanmarcke等介紹的尋求等價平穩運動的方法[7]，換出與歷時曲線相應的等價平穩運動的功率譜，同時也可求得這一平穩運動的持續時間。其功率譜曲線見圖5～7.

圖5 輸入塔夫脫波的加速度(0.2g)功率譜曲線

圖6 輸入水平向唐山波加速度(0.2g)功率譜曲線

由輸入加速度功率譜曲線可見，在相同的最大加速度情況下，唐山波的最大譜值明顯比塔夫脫的最大譜值小，而頻帶(約35Hz)明顯比塔夫脫波的頻帶(約10Hz)寬；塔夫脫波在最大加速度不同情況下的頻譜特性相同，只是最大加速度大的其譜值也大，圖中只表示出了最大加速度為0.2g的功率譜曲線。用前述方法進行堆石壩的隨機地震反應分析，分別取動應力的平均幅值和平均最大值與靜應力迭加進行堆石壩的隨機動力穩定性分析，求出相應的最危險滑裂面及相應的抗滑穩定安全系數，見圖8～12.同時，為了便於比較，靜力狀態下的最危險滑裂面及最小安全系數在圖8及圖12給出。

圖7 輸入豎直向唐山波加速度(0.133g)功率譜曲線

AFL 動應力取平均幅值計算出的最危險滑裂面，K=1.13
MFL 動應力取平均最大值計算出的最危險滑裂面，K=0.64
STA 靜力狀態下計算出的最危險滑裂面，K=1.96
圖8 輸入塔夫脫波加速度為0.2g時的滑裂面

AFL 動應力取平均幅值計算出的最危險滑裂面，K=0.86
MFL 動應力取平均最大值計算出的最危險滑裂面，K=0.37
圖9 輸入塔夫脫波加速度為0.4g時的滑裂面

AFL 動應力取平均幅值計算出的最危險滑裂面，K=0.49
MFL 動應力取平均最大值計算出的最危險滑裂面，K=0.17
圖10 輸入塔夫脫波加速度為0.6g時的滑裂面

SAFL 唐山波(水平向)動應力取平均幅值計算出的最危險滑裂面，K=1.36 SMFL 唐山波(水平向)動應力取平均最大值計算出的最危險滑裂面，K=0.63 TAFA 塔夫脫波動應力取平均幅值算出的最危險滑裂面，K=1.13TMFL 塔夫脫波動應力取平均最大值計算出的最危險滑裂面，K=0.64
圖11 輸入不同波(加速度均為0.2g，唐山波只考慮水平向)滑裂面的比較

HAFL 唐山波(水平向)動應力取平均幅值計算出的最危險滑裂面，K=1.36 HMFL 唐山波(水平向)動應力取平均最大值計算出的最危險滑裂面，K=0.63 VAFA 唐山波(兩向)動應力取平均幅值算出的最危險滑裂面，K=1.30 VMFL 唐山波(兩向)動應力取平均最大值計算出的最危險滑裂面，K=0.62 STA 靜力狀態下計算出的最危險滑裂面，K=1.96
圖12 輸入唐山波的水平向(0.2g)與輸入唐山波兩向
(水平向，0.2g；豎直向，0.133g)滑裂面的比較

由計算結果可以看到：

(1)無論取動應力平均幅值，還是取動應力的平均最大值計算出的最小安全系數比靜力狀態下的要小，最危險滑裂面較為接近於壩頂及壩坡的表面，這與試驗觀察得到的結果是一致的[8]；只是取動應力的平均幅值計算出的最危險滑裂面與靜力狀態下的差不多或稍淺，而取動應力的平均最大值計算出的最危險滑裂面位置卻較深一些，這在輸入加速度較大時更為明顯。(2)取動應力平均幅值計算出的最危險滑裂面接近一直線，而取動應力的平均最大值計算出的最危險滑裂面較為接近圓弧。(3)對於輸入的同一加速度功率譜，一般來說，取動應力的平均最大值計算出的最危險滑裂面比以動應力平均幅值計算出最危險滑裂面要深一些，在輸入地震動強度較大時尤其如此，見圖9及圖10.

對於輸入同樣頻譜特性的塔夫脫波加速度功率譜，其最大加速度越大，計算出的最小安全系數越小，最危險滑裂面相對也越深，見圖8～圖10.

對於輸入最大加速度相同的不同地震波，當動應力取平均幅值計算時，功率譜值越大(塔夫脫波)，最小安全系數越小，而最危險滑裂面位置相差不大；當動應力取平均最大值計算時，功率譜值越大(塔夫脫波)，最危險滑裂面位置也越深，見圖11.

從圖12可以看出，對唐山波而言，只輸入水平向地震波與同時輸入水平向和豎直向地震波，計算出來的最危險滑裂面及最小安全系數差別甚微。這說明，在一般的計算應用中，只考慮水平向地震波是可行的。

4 結語與討論

由前述分析可見，本文所進行的隨機地震作用下邊坡的穩定分析是合理和有效的。它具有以下特點：(1)隨機振動反應得出的是動應力的平均值，它包含了大量歷時曲線統計的平均，比單一的歷時曲線響應分析得出的結果更具有代表性和普遍性。(2)求出的動應力的平均幅值和平均最大值較為直觀，用來分析邊坡的穩定性得出的結果也比較直觀明了，而不象時程分析那樣繁瑣。當然，這種直觀是在統計平均的意義上的。(3)隨機振動反應輸入的是功率譜，這實際上是從能量的角度來分析問題。對於類似的動力作用過程，如果對應的功率譜能量是已知的，則可用這種方法作類似的分析。

參考文獻

1 吳再光，韓國城，林皋。隨機土動力學概論。大連：大連理工大學出版社，1992..

2 吳再光。地基土石壩隨機地震反應及動力穩定性的概率分析〔學位論文〕。大連：大連工學院，1987.

3 劉文廷。土石壩隨機地震反應分析〔學位論文〕。大連：大連理工大學，1993.

4 邵龍潭，韓國城。堆石壩邊坡穩定分析的一種方法。大連理工大學學報，1994，34(3).

5 邵龍潭，韓國城。水流作用下堆石邊坡的穩定分析方法。水利學報，1997，(1).

6 JUN DONG. STUDY ON DYNAMIC SLOPE STABILITY OF FILL-TYPE DAM MODELS.A Dissertation submitted for the Degree of Doctor of Engineering at the Graate School of Civil Engineering University of Tokyo.

7 Vanmarcke E H, Lai S S. Strong Motion Duration and RMS Amplitude of Earthquake. BSSA, 1980,70(4).

8 韓國城，孔憲京，李寇。面板堆石壩動力破壞性態及抗震措施試驗研究。水利學報，1994，(12).

9 林家浩。隨機地震響應的確定性演算法。地震工程與工程振動，1985，5(1).

10 吳再光。土層隨機地震反應的一種改進演算法。振動工程學報，1990，3(1).

http://www.cws.net.cn/Journal/slxb/199911/12.html
裡面有圖的 可以仔細看看

Ⅶ 邵龍潭的主要學歷

1980-1984，大連理工大學土木工程系（原水利工程御銀系）學習，學士學位； 1987-1989，大連理工大學土木工程系學習，碩士學位；
1991-1996，大連理工大學土木工程系學習，博士學位；
1998-1999，德國KARLSRUHE大學土力學與岩石力學研究所，訪問學者；
1984-1990，大連理工大學土木工程系，助宏拆備教；
1990-1993，大連理工大學土木工程系，講師；
1993-1998，大連理工大學土木工程系，副教授；
1998-2001，蔽毀大連理工大學土木工程系，教授，2000年增列為博士生導師。
2001至今， 大連理工大學工程力學系，教授，博士生導師。