|
|
|
| ALGORITHM AND NUMERICAL ANALYSIS METHOD FOR NONLINEAR CONTACT PROBLEMS IN UNDERGROUND ENGINEERING |
| ZHOU Mozhen1,QIAN Xiaoxiang1,2,ZHANG Bingyin1 |
| (1. State Key Laboratory of Hydroscience and Engineering,Tsinghua University,Beijing 100084,China;2. Beijing Rail and Transit Design and Research Institute Company Limited,Beijing 100089,China) |
|
|
|
|
Abstract The node-to-segment approach using penalty method was widely applied in engineering field due to its great simplicity. To deal with the problem in which a large penalty factor may cause ill-condition,the degree of freedom relevant to contact is described by relative displacement in local coordinate system,which enlarges significantly the upper-limit of penalty factor and provides an efficient solution for large-scale contact problem without using any special preconditions. A general 3D contact smoothing method based on radial point interpolation is proposed to avoid probable projecting singularity and jump of contact force in node-to-segment. This method reproduces smooth surfaces which passes exactly through the surface nodes even for coarse or hybrid meshes. The doubly connected edge list is applied to the local searching for contact information and effectively reduces the computing cost. A nonlinear contact analysis program is written and the results of numerical examples indicate that the program effectively describes the phenomenon of displacement discontinuity,solve the contact nonlinear equation,and tackle the large-scale contact problems in underground engineering.
|
|
Received: 21 October 2013
|
|
|
|
| [1] 谢和平,陈忠辉,易 成,等. 基于工程体–地质体相互作用的接触面变形破坏研究[J]. 岩石力学与工程学报,2008,27(9):1 767–1 780.(XIE Heping,CHEN Zhonghui,YI Cheng,et al. Research on deformation and failure of interface based on interaction between structural body and geo-body[J]. Chinese Journal of Rock Mechanics and Engineering,2008,27(9):1 767–1 780.(in Chinese))
[2] 张奇华,邬爱清,石根华. 关键块体理论在百色水利枢纽地下厂房岩体稳定性分析中的应用[J]. 岩石力学与工程学报,2004,23(15):2 609–2 614.(ZHANG Qihua,WU Aiqing,SHI Genhua. Application of key block theory to analysis of rock stability for underground plant in Baise hydraulic project[J]. Chinese Journal of Rock Mechanics and Engineering,2004,23(15):2 609–2 614.(in Chinese))
[3] 余华中,阮怀宁,褚卫江. 岩石节理剪切力学行为的颗粒流数值模拟[J]. 岩石力学与工程学报,2013,32(7):1 482–1 490.(YU Huazhong,RUAN Huaining,CHU Weijiang. Particle flow code modeling of shear behavior of rock joints[J]. Chinese Journal of Rock Mechanics and Engineering,2013,32(7):1 482–1 490.(in Chinese))
[4] 李宇杰,何 平,秦东平,等. 基于三维非连续几何模型的盾构施工模拟[J]. 北京交通大学学报,2012,36(4):38–43.(LI Yujie,HE Ping,QIN Dongping,et al. Simulation of shield tunnel construction process base on 3D discontinuous geometry model[J]. Journal of Beijing Jiaotong University,2012,36(4):38–43.(in Chinese))
[5] GOODMAN R E,TAYLOR R L,BREKKE T L. A model for the mechanics of jointed rock[J]. Journal of Soil Mechanics and Foundations Division,1968,94(3):637–659.
[6] DESAI C S,ZAMAN M M,LIGHTNER J G,et al. Thin-layer element for interfaces and joints[J]. International Journal for Numerical and Analytical Methods in Geomechanics,1984,8(1):19–43.
[7] 张丙印,师瑞锋,王 刚. 高面板堆石坝面板脱空问题的接触力学分析[J]. 岩土工程学报,2003,25(3):361–364.(ZHANG Bingyin,SHI Ruifeng,WANG Gang. Contact analysis of separation between concrete slab and cushion layer in high concrete face rock-fill dam[J]. Chinese Journal of Geotechnical Engineering,2003,25(3):361–364. (in Chinese))
[8] HALLQUIST J O,GOUDREAU G L,BENSON D J. Sliding interfaces with contact-impact in large-scale Lagrangian computations[J]. Computer Methods in Applied Mechanics and Engineering,1985,51(1):107–137.
[9] ZAVARISE G,WRIGGERS P. A segment-to-segment contact strategy[J]. Mathematical and Computer Modelling,1998,28(4):497–515.
[10] BELYTSCHKO T,NEAL M O. Contact-impact by the pinball algorithm with penalty and Lagrangian methods[J]. International Journal for Numerical Methods in Engineering,1991,31(3):547–572.
[11] ZAVARISE G,WRIGGERS P,SCHREFLER B A. On augmented Lagrangian algorithms for thermomechanical contact problems with friction[J]. International Journal for Numerical Methods in Engineering,1995,38(17):2 929–2 949.
[12] WOHLMUTH B I. A mortar finite element method using dual spaces for the Lagrange multiplier[J]. SIAM Journal on Numerical Analysis,2000,38(3):989–1 012.
[13] POPP A,SEITZ A,GEE M W,et al. Improved robustness and consistency of 3D contact algorithms based on a dual mortar approach[J]. Computer Methods in Applied Mechanics and Engineering,2013,264:67–80.
[14] 王勖成. 有限单元法[M]. 北京:清华大学出版社,2003:669–675.(WANG Xucheng. Finite element method[M]. Beijing:Tsinghua University Press,2003:669–675.(in Chinese))
[15] RAHMAN M U,ROWLANDS R E,COOK R D,et al. An iterative procedure for finite-element stress analysis of frictional contact problems[J]. Computers and Structures,1984,18(6):947–954.
[16] JANNA C,FERRONATO M,GAMBOLATI G. Two-level block preconditioners for contact problems[M]. [S.l.]:Springer Berlin Heidelberg,2011:211–226.
[17] QIAN X,YUAN H,ZHOU M,et al. A general 3D contact smoothing method based on radial point interpolation[J]. Journal of Computational and Applied Mathematics,2014,257:1–13 |
| [1] |
CHEN Weizhong1*, LIU Xinyu1, 2, YANG Jianping1, WANG Wei1, 2, ZANG Zhonghai3, DING Hongyuan3, ZHANG Zheyuan3, WANG Xiaogang3, SHI Zhengrong1. Development of a large-scale 3D physical model test system for underground energy storage caverns and its model experimental study[J]. , 2026, 45(6): 1615-1628. |
|
|
|
|