×

A frictional spring and cohesive contact model for accurate simulation of contact forces in numerical manifold method. (English) Zbl 07841923

Summary: “Open-close iteration” is a crucial algorithm for handling complex contacts in numerical manifold method (NMM) and discontinuous deformation analysis (DDA). This algorithm has proved to be robust and efficient for decades. However, as some researchers have pointed out, the original open-close iteration may involve errors in sliding tests, especially in critical sliding tests with cohesive contacts. In this study, two major problems in the original algorithm are found to be nonconvergent contact force and early removed cohesive strength. The modifications are the following: (a) a frictional spring. By avoiding the trail value of normal contact force, we added a new frictional spring to the iteration scheme. This spring can apply accurate friction and can help ensure the convergence of contact forces. (b) A cohesive contact model. The original scheme can encounter an “early failure” in cohesive contacts. After investigating how contacts provide shear resistance, we found the cause and then provided a simple correction of the cohesive issue. The new algorithms in this article are essential for accurately simulating contacts by NMM/DDA.
{© 2020 John Wiley & Sons, Ltd.}

MSC:

74Sxx Numerical and other methods in solid mechanics
74Mxx Special kinds of problems in solid mechanics
65Nxx Numerical methods for partial differential equations, boundary value problems
Full Text: DOI

References:

[1] ShiGH. Contact theory. Sci China Technol Sci. 2015;58(9):1450‐1496.
[2] BoonCW, HoulsbyGT, UtiliS. A new algorithm for contact detection between convex polygonal and polyhedral particles in the discrete element method. Comput Geotech. 2012;44:73‐82.
[3] ZhengF, JiaoY, SitarN. Generalized contact model for polyhedra in three‐dimensional discontinuous deformation analysis. Int J Numer Anal Methods Geomech. 2018;42(13):1471‐1492.
[4] BatheK, ChaudharyA. A solution method for planar and axisymmetric contact problems. Int J Numer Methods Eng. 1985;21(1):65‐88. · Zbl 0551.73099
[5] ArmeroF, PetoczE. Formulation and analysis of conserving algorithms for frictionless dynamic contact impact problems. Comput Methods Appl Mech Eng. 1998;158(3):269‐300. · Zbl 0954.74055
[6] ArmeroF, PetoczE. A new dissipative time‐stepping algorithm for frictional contact problems: formulation and analysis. Comput Methods Appl Mech Eng. 1999;179(1):151‐178. · Zbl 0980.74044
[7] JonesRE, PanayiotisP. A novel three‐dimensional contact finite element based on smooth pressure interpolations. Int J Numer Methods Eng. 2000;51(7):791‐811. · Zbl 1024.74043
[8] LinXC, LiX, WangXG, WangYJ. A compact 3D block cutting and contact searching algorithm. Sci China Technol Sci. 2019;62(8):1438‐1454.
[9] WriggersP, SimoJC. A note on tangent stiffness for fully nonlinear contact problems. Commun Appl Num Methods. 1985;1(5):199‐203. · Zbl 0582.73110
[10] CundallPA. A computer model for simulating progressive, large‐scale movement in blocky rock system. Proc Int Symp Rock Mech. Nancy, France; 1971:128‐132.
[11] CundallPA, StrackOD. A discrete numerical model for granular assemblies. Geotechnique. 1979;29(1):47‐65.
[12] NocedalJ, WrightSJ. Numerical Optimization. Beijing. China: China Science Publishing & Media Ltd; 2016:492‐532.
[13] ShiGH. Discontinuous Deformation Analysis—A New Numerical Model for the Statics and Dynamics of Block Systems. [Ph. D thesis]. Berkeley, CA: Department of Civil Engineering University of California; 1988.
[14] ShiGH. Discontinuous deformation analysis: a new numerical model for the statics and dynamics of deformable block structures. Eng Comput. 1992;9(2):157‐168.
[15] ShiGH. Manifold Method of Material Analysis. Research Triangle Park, NC: Army Research Office; 1992.
[16] MaGW, AnXM, HeL. The numerical manifold method: a review. Int J Comput Methods. 2010;7(01):1‐32. · Zbl 1267.74129
[17] LinCT, AmadeiB, StureS, et al. Using an augmented Lagrangian method and block fracturing in the DDA Method. Albuquerque, NM: Sandia National Labs; 1994.
[18] YuY, YinJ. Some modifications to the process of discontinuous deformation analysis. J Rock Mech Geotech Eng. 2015;7(1):95‐100.
[19] LiX, ZhengH. Condensed form of complementarity formulation for discontinuous deformation analysis. Sci China Technol Sci. 2015;58(9):1509‐1519.
[20] ZhengH, ZhangP, DuX. Dual form of discontinuous deformation analysis. Comput Methods Appl Mech. 2016;305:196‐216. · Zbl 1425.74507
[21] PusoMA. A 3D mortar method for solid mechanics. Int J Numer Methods Eng. 2004;59(3):315‐336. · Zbl 1047.74065
[22] YangB, LaursenTA, MengX. Two dimensional mortar contact methods for large deformation frictional sliding. Int J Numer Methods Eng. 2005;62(9):1183‐1225. · Zbl 1161.74497
[23] LaursenTA, PusoMA, SandersJ. Mortar contact formulations for deformable-deformable contact: past contributions and new extensions for enriched and embedded interface formulations. Comput Methods Appl Mech Eng. 2012;205‐208:3‐15. · Zbl 1239.74070
[24] WriggersP. Computational Contact Mechanics. 2nd ed.Berlin, New York: Springer; 2006. · Zbl 1104.74002
[25] LaursenTA. Computational Contact and Impact Mechanics: Fundamentals of Modeling Interfacial Phenomena in Nonlinear Finite Element Analysis. 2nd ed.Berlin, New York: Springer; 2011.
[26] MacLaughlinMM, DoolinDM. Review of validation of the discontinuous deformation analysis (DDA) method. Int J Numer Anal Methods Geomech. 2006;30(4):271‐305. · Zbl 1140.74574
[27] YeungMCR. Application of Shi’s discontinuous deformation analysis to the study of rock behaviour. [Ph.D. Thesis]. Berkeley, CA: University of California; 1991.
[28] DoolinDM, SitarN. Displacement accuracy of discontinuous deformation analysis method applied to sliding block. J Eng Mech. 2002;128(11):1158‐1168.
[29] WangLZ, JiangHY, YangZX, XuYC, ZhuXB. Development of discontinuous deformation analysis with displacement‐dependent interface shear strength. Comput Geotech. 2013;47:91‐101.
[30] ZhangY, XuQ, ChenG, ZhaoJX, ZhengL. Extension of discontinuous deformation analysis and application in cohesive‐frictional slope analysis. Int J Rock Mech Min Sci. 2014;70:533‐545.
[31] CaiY, HeT, WangR. Numerical simulation of dynamic process of the Tangshan earthquake by a new method—LDDA. Pure Appl Geophys. 2000;157(11):2083‐2104.
[32] WriggersP, ZavariseG. Application of augmented Lagrangian techniques for non‐linear constitutive laws in contact interfaces. Commun Numer Methods Eng. 1993;9(10):815‐824. · Zbl 0788.73077
[33] ZhengH, LiJ. A practical solution for KKT systems. Numer Algorithms. 2007;46(2):105‐119. · Zbl 1128.65026
[34] DoolinDM, SitarN. Time integration in discontinuous deformation analysis. J Eng Mech. 2004;130(3):249‐258.
[35] JiangW, ZhengH. An efficient remedy for the false volume expansion of DDA when simulating large rotation. Comput Geotech. 2015;70:18‐23.
[36] ZhangN, LiX, JiangQH, LinXC. Rotation errors in numerical manifold method and a correction based on large deformation theory. Front Struct Civ Eng. 2019;13(5):1036‐1053.
[37] ZienkiewiczOC, TaylorRL, FoxDD. The finite element method for solid & structural mechanics. Singapore: Elsevier; 2015.
[38] DoolinDM. Unified displacement boundary constraint formulation for discontinuous deformation analysis (DDA). Int J Numer Anal Methods Geomech. 2005;29(12):1199‐1207. · Zbl 1140.74569
[39] BartonNR. A model study of rock‐joint deformation. Int J Rock Mech Mining Sci Geomech Abstracts. 1972;9(5):579‐582.
[40] BandisS, LumsdenAC, BartonNR. Experimental studies of scale effects on the shear behaviour of rock joints. Int J Rock Mech Mining Sci Geomech Abstracts. 1981;18(1):1‐21.
[41] PapaliangasT, HencherSR, LumsdenAC, ManolopoulouS. The effect of frictional fill thickness on the shear strength of rock discontinuities. Int J Rock Mech Mining Sci Geomech Abstracts. 1993;30(2):81‐91.
[42] ZhengH, ChengSG, LiuQS. A decomposition procedure for nearly‐symmetric matrices with applications to some nonlinear problems. Mech Res Commun. 2010;37(1):78‐84. · Zbl 1272.74512
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.