大连理工大学邵龙潭教授脑出血

Ⅰ 邵龙潭的学术著作

[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至今， 大连理工大学工程力学系，教授，博士生导师。