岩土力学 ›› 2020, Vol. 41 ›› Issue (11): 3810-3822.doi: 10.16285/j.rsm.2020.0189

• 数值分析 • 上一篇    下一篇

非连续变形分析中块体大转动问题研究

巩师林1, 2,凌道盛1, 2, 3,胡成宝1, 2,钮家军1, 2   

  1. 1. 浙江大学 岩土工程研究所,浙江 杭州 310058;2. 浙江大学 软弱土与环境土工教育部重点实验室,浙江 杭州 310058; 3. 浙大宁波理工学院 土木建筑工程学院,浙江 宁波 315100
  • 收稿日期:2020-01-09 修回日期:2020-04-13 出版日期:2020-11-11 发布日期:2020-12-25
  • 通讯作者: 凌道盛,男,1968年生,博士,教授,博士生导师,主要从事土动力学、计算土力学方面的研究工作。E-mail: dsling@zju.edu.cn E-mail:slgong@zju.edu.cn
  • 作者简介:巩师林,男,1993年生,博士研究生,主要从事岩石力学与非连续变形分析方面的研究工作
  • 基金资助:
    国家重点研发计划项目(No. 2016YFC0800200);国家自然科学基金资助项目(No. 51578502)

Study on large rotation problem of a block in discontinuous deformation analysis

GONG Shi-lin1, 2, LING Dao-sheng1, 2, 3, HU Cheng-bao1, 2, NIU Jia-jun1, 2   

  1. 1. Institute of Geotechnical Engineering, Zhejiang University, Hangzhou, Zhejiang 310058, China; 2. MOE Key Laboratory of Soft Soils and Geoenvironmental Engineering, Zhejiang University, Hangzhou, Zhejiang 310058, China; 3. School of Civil Engineering and Architecture, NingboTech University, Ningbo, Zhejiang 315100, China
  • Received:2020-01-09 Revised:2020-04-13 Online:2020-11-11 Published:2020-12-25
  • Supported by:
    The work was supported by the National Key Research and Development Program of China (2016YFC0800200) and the National Natural Science Foundation of China (51578502).

摘要: 原有非连续变形分析(DDA)采用一阶近似后的位移增量表达式更新块体构形,推导相关子矩阵,且对不同时步计算出的应变增量直接叠加,当模拟的块体发生大转动时往往会产生较大误差。为考虑块体转动与变形的耦合作用,引入先变形、后转动的块体位移增量表达式。重新推导了惯性力子矩阵,将块体转动时的离心力与科氏力加到荷载矩阵中。计算时对应变分量及其相关变量进行坐标变换与修正,并采用新引入的位移增量表达式计算块体顶点位移,进行后接触修正与更新块体构形。数值算例表明,改进后的程序能够消除转动带来的误差,自动考虑了块体转动时离心力和科氏力引起的变形,应变计算精度更高。改进方法克服了块体体积自由膨胀、应变场畸变等问题,给出了合理的块体应变。

关键词: 非连续变形分析(DDA), 转动误差, 坐标变换, 离心力, 科氏力

Abstract: The original discontinuous deformation analysis (DDA) takes the first-order approximation of the incremental displacement to update block configuration, deduces related sub-matrices, and directly superimposes the incremental strains calculated at different time steps, which generates significant errors when simulating large rotation of a block. In order to consider the coupling effect of block rotation and deformation, a modified incremental displacement formulation is introduced, which calculates the displacements caused by the strains before the rigid-body rotation term. The sub-matrix of the inertial force is then modified, and the centrifugal and Coriolis forces are added to the load matrix. In the calculation, the coordinate transformation is carried out for the strain components and its related variables, and the newly introduced displacement formulation is used to calculate the displacements of block vertices. After that, the post-adjustment and post-contact adjustment are used to update the block configuration. The numerical examples indicate that the modified DDA can eliminate the errors caused by large rotation, automatically consider the deformation caused by centrifugal and Coriolis forces of the rotating block, and the precision of strain calculation is higher than that of original DDA. The modified codes overcome the free volume expansion and strain distortion problems, and give a reasonable block strain.

Key words: discontinuous deformation analysis (DDA), rotation error, coordinate transformation, centrifugal force, Coriolis force

中图分类号: O 241
[1] 刘泉声,蒋亚龙,何 军. 非连续变形分析的精度改进方法及研究趋势[J]. , 2017, 38(6): 1746-171.
[2] 黄小福,张迎宾,赵兴权,余鹏程,邢 昊,张 珏,陈岩岩,. 地震条件下危岩崩塌运动特性的初步探讨[J]. , 2017, 38(2): 583-592.
[3] 付晓东,盛 谦,张勇慧. 非连续变形分析方法中的阻尼问题研究[J]. , 2015, 36(7): 2057-2062.
[4] 邬爱清,冯细霞,卢 波. 非连续变形分析中时间步及弹簧刚度取值研究[J]. , 2015, 36(3): 891-897.
[5] 魏永权 ,罗 强 ,张 良 ,刘 钢 ,张立祥 ,翁骏祺,. 离心力场中微型土压力传感器非线性响应分析[J]. , 2015, 36(1): 286-292.
[6] 桂乐乐,周宏伟,夏昌敬,杜盛浩,刘 艳. 基于尺度律的裂隙网络生成及其对地下厂房稳定性的影响分析[J]. , 2010, 31(12): 3934-3938.
[7] 王永甫,王 成. 曲线桥下斜坡受邻近基坑开挖的影响研究[J]. , 2008, 29(S1): 421-425.
[8] 刘 君,黄盛铨,孔宪京. 强度折减DDA方法及其在稳定分析中的应用[J]. , 2006, 27(S2): 286-290.
[9] 沈振中,林伟斌,杨海林. 考虑降雨入渗的非连续性岩体边坡稳定分析[J]. , 2006, 27(S2): 269-274.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!