›› 2014, Vol. 35 ›› Issue (6): 1787-1794.

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

基于线性互补的非连续变形分析

李小凯,郑 宏   

  1. 中国科学院武汉岩土力学研究所 岩土力学与工程国家重点实验室,武汉 430071
  • 收稿日期:2013-03-03 出版日期:2014-06-10 发布日期:2014-06-20
  • 作者简介:李小凯,男,1975年生,博士研究生,主要从事岩土工程及数值分析方面的研究工作。
  • 基金资助:

    国家重点基础研究发展计划(973)项目(No. 2011CB013505);国家杰出青年基金项目(No.50925933)。

Discontinuous deformation analysis based on linear complementarity theory

LI Xiao-kai, ZHENG Hong   

  1. State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanic, Chinese Academy of Sciences, Wuhan 430071, China
  • Received:2013-03-03 Online:2014-06-10 Published:2014-06-20

摘要: 非连续变形分析(DDA)方法是一种新的用来分析块体系统运动和变形的非连续介质数值计算方法。研究的核心工作是致力于对现有DDA接触问题处理方法的改进。DDA主要采用罚函数法和Lagrange乘子法处理接触问题,合理设定罚参数很困难,此外,因开闭迭代而引起的刚度矩阵的不连续变化也会导致收敛方面的困难。为避免引入罚参数及传统意义上的开闭迭代,用混合线性互补模型(LCDDA)对DDA方法进行了重新描述。在此基础上,综合基于非光滑分析的Newton法的局部平方收敛和最速下降法的全局线性收敛的优势,提出求解LCDDA模型的有效算法。根据上述思想及理论研究成果编制了完整的计算程序,算例计算结果证明了方法的精度及可行性。

关键词: 非连续变形分析, 接触问题, 线性互补理论, 半光滑函数

Abstract: Discontinuous deformation analysis (DDA) is a newly developed discontinuum numerical method for simulating large deformation and displacement of block systems. The core research of this thesis concentrates on improving the method to enforce contact constraints in DDA. The penalty function method and Lagrange multiplier method or its variants are generally utilized to enforce contact constraints in DDA. But it is difficult to set a reasonable penalty parameter. Furthermore, discontinuous change in stiffness matrix due to the open-close iteration frequently causes poor convergence. To avoid troubles in setting artificial penalty parameters and open-close iteration, we reformulate discontinuous deformation analysis method as mixed linear complementarity model (LCDDA). On the basis of these, both advantages of local quadratic convergence rate of the Newton method based on nonsmooth analysis and global linear convergence rate of steepest descent method are integrated to set up the algorithm. According to the above achievement, the complete code for LCDDA is accomplished. Some examples are analyzed to verify the precision and feasibility.

Key words: discontinuous deformation analysis (DDA), contact problem, linear complementarity theory, semismooth function

中图分类号: 

  • O 242.21
[1] 徐栋栋, 邬爱清, 李 聪, 汪 斌, 蒋昱州, 曾 平, 杨永涛, . 破裂全过程模拟的改进非连续变形分析方法[J]. 岩土力学, 2019, 40(3): 1169-1178.
[2] 赵强, 焦玉勇, 张秀丽, 谢壁婷, 王龙, 黄刚海, . 基于显式时间积分的球颗粒DDA计算方法[J]. 岩土力学, 2019, 40(11): 4515-4522.
[3] 林 姗,李春光,孙冠华,王水林,杨永涛,. 二维接触问题的互补算法及工程应用[J]. , 2018, 39(10): 3863-3874.
[4] 刘泉声,蒋亚龙,何 军. 非连续变形分析的精度改进方法及研究趋势[J]. , 2017, 38(6): 1746-171.
[5] 黄小福,张迎宾,赵兴权,余鹏程,邢 昊,张 珏,陈岩岩,. 地震条件下危岩崩塌运动特性的初步探讨[J]. , 2017, 38(2): 583-592.
[6] 付晓东,盛 谦,张勇慧,冷先伦, . 非连续变形分析(DDA)线性方程组的高效求解算法[J]. , 2016, 37(4): 1171-1178.
[7] 付晓东,盛 谦,张勇慧. 非连续变形分析方法中的阻尼问题研究[J]. , 2015, 36(7): 2057-2062.
[8] 邬爱清,冯细霞,卢 波. 非连续变形分析中时间步及弹簧刚度取值研究[J]. , 2015, 36(3): 891-897.
[9] 虞 松 ,朱维申 ,张云鹏,. 基于DDA方法一种流-固耦合模型的建立及裂隙体渗流场分析和应用[J]. , 2015, 36(2): 555-560.
[10] 江 贝 ,李术才 ,王 琦 ,朱维申 ,王德超 , 王富奇 ,王洪涛 ,阮国强 ,邵 行,. 基于非连续变形分析方法的深部沿空掘巷围岩变形破坏及控制机制对比研究[J]. , 2014, 35(8): 2353-2360.
[11] 付晓东,盛 谦,张勇慧. 基于OpenMP的非连续变形分析并行计算方法[J]. , 2014, 35(8): 2401-2407.
[12] 李术才 ,陈云娟 ,朱维申 ,赵成龙 ,李利平 ,尹福强,. 岩石破坏非连续变形分析弹脆性本构模型的二次开发[J]. , 2014, 35(8): 2144-2149.
[13] 王 文,朱维申,马海萍,张 磊. 不同倾角节理组和锚固效应对岩体特性的影响[J]. , 2013, 34(3): 887-893.
[14] 付晓东,盛 谦,张勇慧. 水电站地下洞室群分步开挖的非连续变形分析[J]. , 2013, 34(2): 568-574.
[15] 刘永茜,杨 军. 一种改进步长自调的非连续变形分析法[J]. , 2011, 32(8): 2544-2548.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 孙树林,李 方,谌 军. 掺石灰黏土电阻率试验研究[J]. , 2010, 31(1): 51 -55 .
[2] 李英勇,张顶立,张宏博,宋修广. 边坡加固中预应力锚索失效机制与失效效应研究[J]. , 2010, 31(1): 144 -150 .
[3] 李 晶,缪林昌,钟建驰,冯兆祥. EPS颗粒混合轻质土反复荷载下变形和阻尼特性[J]. , 2010, 31(6): 1769 -1775 .
[4] 梁健伟,房营光,谷任国. 极细颗粒黏土渗流的微电场效应分析[J]. , 2010, 31(10): 3043 -3050 .
[5] 王丽艳,姜朋明,刘汉龙. 砂性地基中防波堤地震残余变形机制分析与液化度预测法[J]. , 2010, 31(11): 3556 -3562 .
[6] 李秀珍,王成华,邓宏艳. DDA法和Fisher判别法在潜在滑坡判识中的应用比较[J]. , 2011, 32(1): 186 -192 .
[7] 陈力华 ,林 志 ,李星平. 公路隧道中系统锚杆的功效研究[J]. , 2011, 32(6): 1843 -1848 .
[8] 陈立文,孙德安. 不同应力路径下水土耦合超固结黏土分叉分析[J]. , 2011, 32(10): 2922 -2928 .
[9] 赵明华,雷 勇,张 锐. 岩溶区桩基冲切破坏模式及安全厚度研究[J]. , 2012, 33(2): 524 -530 .
[10] 王松鹤,齐吉琳. 高温冻土松弛特性试验研究[J]. , 2012, 33(6): 1660 -1666 .