基于Newmark隐式时间积分方案的裂纹动态扩展的数值计算方法

doi:10.16183/j.cnki.jsjtu.2020.021

上海交通大学学报 ›› 2021, Vol. 55 ›› Issue (6): 689-697.doi: 10.16183/j.cnki.jsjtu.2020.021

所属专题：《上海交通大学学报》2021年“土木建筑工程”专题；《上海交通大学学报》2021年12期专题汇总专辑

基于Newmark隐式时间积分方案的裂纹动态扩展的数值计算方法

郭德平¹^,², 李铮³(), 彭森林⁴, 曾志凯³, 吴岱峰³

1.叙镇铁路有限责任公司, 云南昭通 657900
2.武汉大学土木建筑工程学院, 武汉 430072
3.重庆市城市建设投资(集团)有限公司, 重庆 400023
4.重庆大学土木工程学院, 重庆 400045

收稿日期:2020-01-14 出版日期:2021-06-28 发布日期:2021-06-30
通讯作者: 李铮 E-mail:2756232300@qq.com
作者简介:郭德平(1983-),男,高级工程师,主要从事隧道工程和岩土工程等方面的科研工作
基金资助:
国家自然科学基金重点项目(51679017);国家自然科学基金重点项目(51839009)

Numerical Calculation Method for Crack Dynamic Propagation Based on Newmark Implicit Time Integration Scheme

GUO Deping¹^,², LI Zheng³(), PENG Senlin⁴, ZENG Zhikai³, WU Daifeng³

1. Xuzhen Railway Co., Ltd., Zhaotong 657900, Yunnan, China
2. School of Civil Engineering, Wuhan University, Wuhan 430072, China
3. Chongqing City Construction Investment (Group) Co., Ltd., Chongqing 400023, China
4. Department of Civil Engineering, Chongqing University, Chongqing 400045, China

Received:2020-01-14 Online:2021-06-28 Published:2021-06-30
Contact: LI Zheng E-mail:2756232300@qq.com

摘要/Abstract

摘要：

扩展有限单元法(XFEM)是基于单位分解的思想,在常规有限元的位移模式中加入能够反映裂纹面不连续性的跳跃函数和裂纹尖端的渐近位移场函数,避免了常规有限单元法计算断裂问题时需要对裂纹尖端重新划分网格的不便以及繁重的计算量,并且裂纹的扩展独立于网格.标准有限元在处理时间积分时,在裂纹不断扩展的过程中整体刚度矩阵的自由度也会不断地增大,从而导致迭代计算无法进行.本文基于扩展有限单元法模拟动态裂纹扩展的方法,提出了新的Newmark隐式时间积分方案.此方法将所有节点都富集Heaviside函数和裂纹尖端的渐近位移场函数,即每个节点都有12个自由度,从而使得总体刚度矩阵式保持一致,避免迭代计算无法进行.同时,提出了一种稀疏矩阵技术来解决矩阵所占内存大和计算时间长的问题.

关键词: 扩展有限元, 动态裂纹, Newmark隐式时间积分, 稀疏矩阵技术

Abstract:

Extended finite element method (XFEM) is based on the idea of unit decomposition. The jump function that can reflect the discontinuity of the crack surface and the progressive displacement field function of the crack tip is added to the conventional finite element displacement mode, which avoids the inconvenience of remeshing the crack tip and the heavy calculation. Then the conventional finite element method calculates the fracture problem, and the crack propagation is independent of the mesh. When the standard finite element deals with time integration, the degree of freedom of the overall stiffness matrix will continue to increase in the process of crack propagation, which makes iterative calculation impossible. This paper proposes a novel Newmark implicit time integration scheme based on the XFEM to simulate dynamic crack growth. This method enriches all the nodes with the Heaviside function and the asymptotic displacement field function at the crack tip, that is, each node has 12 degrees of freedom, so that the overall stiffness matrix is consistent without making iterative calculation impossible. At the same time, a sparse matrix technology is proposed to solve the problems of large memory and long calculation time occupied by the matrix.

Key words: extended finite element method (XFEM), dynamic cracks, Newmark implicit time integration scheme, sparse matrix technique

中图分类号:

TU43

郭德平, 李铮, 彭森林, 曾志凯, 吴岱峰. 基于Newmark隐式时间积分方案的裂纹动态扩展的数值计算方法[J]. 上海交通大学学报, 2021, 55(6): 689-697.

GUO Deping, LI Zheng, PENG Senlin, ZENG Zhikai, WU Daifeng. Numerical Calculation Method for Crack Dynamic Propagation Based on Newmark Implicit Time Integration Scheme[J]. Journal of Shanghai Jiao Tong University, 2021, 55(6): 689-697.

图/表 12

图1

图2

图3

图4

图5

图6

图7

图8

图9

图10

图11

图12

参考文献 23

[1]	SHEN J H, WHEELER C, ILIC D, et al. Application of open source FEM and DEM simulations for dynamic belt deflection modelling[J]. Powder Technology, 2019, 357:171-185. doi: 10.1016/j.powtec.2019.08.068 URL
[2]	ELGUEDJ T, JAN Y, COMBESCURE A, et al. X-FEM Analysis of dynamic crack growth under transient loading in thick shells[J]. International Journal of Impact Engineering, 2018, 122:228-250. doi: 10.1016/j.ijimpeng.2018.08.013 URL
[3]	SHARMA V, FUJISAWA K, MURAKAMI A. Space-time FEM with block-iterative algorithm for nonlinear dynamic fracture analysis of concrete gravity dam[J]. Soil Dynamics and Earthquake Engineering, 2020, 131:105995. doi: 10.1016/j.soildyn.2019.105995 URL
[4]	SUN J S, LEE K H, LEE H P. Comparison of implicit and explicit finite element methods for dynamic problems[J]. Journal of Materials Processing Technology, 2000, 105(1/2):110-118. doi: 10.1016/S0924-0136(00)00580-X URL
[5]	NILSSON K, LIDSTRÖM P. Simulation of ductile fracture of slabs subjected to dynamic loading using cohesive elements[J]. International Journal of Da-mage Mechanics, 2012, 21(6):871-892.
[6]	BELYTSCHKO T, BLACK T. Elastic crack growth in finite elements with minimal remeshing[J]. International Journal for Numerical Methods in Engineering, 1999, 45(5):601-620. doi: 10.1002/(ISSN)1097-0207 URL
[7]	ZHOU X P, ZHANG J Z, BERTO F. Fracture ana-lysis in brittle sandstone by digital imaging and AE techniques: Role of flaw length ratio[J]. Journal of Materials in Civil Engineering, 2020, 32(5):04020085. doi: 10.1061/(ASCE)MT.1943-5533.0003151 URL
[8]	ZHOU X P, CHENG H. Multidimensional space method for geometrically nonlinear problems under total Lagrangian formulation based on the extended finite-element method[J]. Journal of Engineering Mechanics, 2017, 143(7):04017036. doi: 10.1061/(ASCE)EM.1943-7889.0001241 URL
[9]	CHEN J W, ZHOU X P. The enhanced extended finite element method for the propagation of complex branched cracks[J]. Engineering Analysis With Boundary Elements, 2019, 104:46-62. doi: 10.1016/j.enganabound.2019.03.028 URL
[10]	STOLARSKA M, CHOPP D L, MOËS N, et al. Modelling crack growth by level sets in the extended finite element method[J]. International Journal for Numerical Methods in Engineering, 2001, 51(8):943-960. doi: 10.1002/nme.201 URL
[11]	ZHOU X P, CHEN J W. Extended finite element simulation of step-path brittle failure in rock slopes with non-persistent en-echelon joints[J]. Engineering Geology, 2019, 250:65-88. doi: 10.1016/j.enggeo.2019.01.012 URL
[12]	CHEN J W, ZHOU X P, BERTO F. The improvement of crack propagation modelling in triangular 2D structures using the extended finite element method[J]. Fatigue & Fracture of Engineering Materials & Structures, 2019, 42(2):397-414.
[13]	黄宏伟, 刘德军, 薛亚东, 等. 基于扩展有限元的隧道衬砌裂缝开裂数值分析[J]. 岩土工程学报, 2013, 35(2):266-275.
	HUANG Hongwei, LIU Dejun, XUE Yadong, et al. Numerical analysis of cracking of tunnel linings based on extended finite element[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(2):266-275.
[14]	阮滨, 陈国兴, 王志华. 基于扩展有限元法的均质土坝裂纹模拟[J]. 岩土工程学报, 2013, 35(Sup.2):49-54.
	RUAN Bin, CHEN Guoxing, WANG Zhihua. Numerical simulation of cracks of homogeneous earth dams using an extended finite element method[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(Sup.2):49-54.
[15]	MENOUILLARD T, BELYTSCHKO T. Dynamic fracture with meshfree enriched XFEM[J]. Acta Mechanica, 2010, 213(1/2):53-69. doi: 10.1007/s00707-009-0275-z URL
[16]	WEN L F, TIAN R. Improved XFEM: Accurate and robust dynamic crack growth simulation[J]. Compu-ter Methods in Applied Mechanics and Engineering, 2016, 308:256-285.
[17]	ZHOU X P, ZHANG J Z, QIAN Q H, et al. Expe-rimental investigation of progressive cracking processes in granite under uniaxial loading using digital imaging and AE techniques[J]. Journal of Structural Geology, 2019, 126:129-145. doi: 10.1016/j.jsg.2019.06.003 URL
[18]	WANG Y T, ZHOU X P, XU X. Numerical simulation of propagation and coalescence of flaws in rock materials under compressive loads using the extended non-ordinary state-based peridynamics[J]. Engineering Fracture Mechanics, 2016, 163:248-273. doi: 10.1016/j.engfracmech.2016.06.013 URL
[19]	WU Z J, FAN L F, LIU Q S, et al. Micro-mechanical modeling of the macro-mechanical response and fracture behavior of rock using the numerical manifold method[J]. Engineering Geology, 2017, 225:49-60. doi: 10.1016/j.enggeo.2016.08.018 URL
[20]	CHEN S, HANSEN J M, TORTORELLI D A. Unconditionally energy stable implicit time integration: Application to multibody system analysis and design[J]. International Journal for Numerical Methods in Engineering, 2000, 48(6):791-822. doi: 10.1002/(ISSN)1097-0207 URL
[21]	WANG J R, WU J C, WANG D D. A quasi-consis-tent integration method for efficient meshfree analysis of Helmholtz problems with plane wave basis functions[J]. Engineering Analysis With Boundary Elements, 2020, 110:42-55. doi: 10.1016/j.enganabound.2019.10.002 URL
[22]	WANG D D, WU J C. An inherently consistent reproducing kernel gradient smoothing framework toward efficient Galerkin meshfree formulation with explicit quadrature[J]. Computer Methods in Applied Mechanics and Engineering, 2019, 349:628-672. doi: 10.1016/j.cma.2019.02.029 URL
[23]	WANG D D, WANG J R, WU J C, et al. A three-dimensional two-level gradient smoothing meshfree method for rainfall induced landslide simulations[J]. Frontiers of Structural and Civil Engineering, 2019, 13(2):337-352. doi: 10.1007/s11709-018-0467-5 URL

基于Newmark隐式时间积分方案的裂纹动态扩展的数值计算方法

Numerical Calculation Method for Crack Dynamic Propagation Based on Newmark Implicit Time Integration Scheme

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

图/表 12

参考文献 23

相关文章 15

编辑推荐

Metrics

本文评价

[1]	杨其润, 李明广, 陈锦剑, 吴航. 同步实施的相邻基坑相互作用机理[J]. 上海交通大学学报, 2022, 56(6): 722-729.
[2]	邱杰凯, 丁肇伟, 宋春雨, 陈龙珠. 成层广义Gibson地基中桩的水平动力响应[J]. 上海交通大学学报, 2022, 56(4): 431-442.
[3]	胡亚元, 王啊强. 饱和孔隙-裂隙黏土双层地基的一维固结分析[J]. 上海交通大学学报, 2022, 56(4): 443-453.
[4]	潘上, 刘谨豪, 张琪, 叶冠林. 三轴仪K₀系数测量与应力路径试验功能的研发[J]. 上海交通大学学报, 2021, 55(4): 372-379.
[5]	应宏伟, 许鼎业, 王迪, 章丽莎. 波动承压水下基坑底部弱透水层的非Darcy渗流分析[J]. 上海交通大学学报, 2020, 54(12): 1300-1306.
[6]	韩红桂, 杨士恒, 张璐, 乔俊飞. 城市污水处理过程出水氨氮优化控制[J]. 上海交通大学学报, 2020, 54(9): 916-923.
[7]	吴威, 李明广, 史玉金, 王建华. 上海地层长期变形特性及其对基础设施的影响[J]. 上海交通大学学报, 2018, 52(11): 1429-1436.
[8]	章红兵，范凡，胡昊. 基坑群施工对邻近隧道影响与隧道保护[J]. 上海交通大学学报（自然版）, 2016, 50(05): 803-809.
[9]	武朝军, 叶冠林, 王建华. 上海莲花路浅部土层超固结特性试验研究[J]. 上海交通大学学报, 2016, 50(03): 331-335.
[10]	程演, 张璐璐, 张磊, 王建华. 基于随机场的非饱和土固结分析[J]. 上海交通大学学报, 2014, 48(11): 1528-1535.
[11]	朱龙, 徐震, 白占伟, 陈锦剑, 王建华. 大顶力作用下圆形沉井后背土体受力变形特性分析[J]. 上海交通大学学报, 2014, 48(11): 1510-1516.
[12]	刘翔, 白海梅, 陈晓晨, 陈锦剑, 王建华. 软土中大直径顶管管道受力特性测试[J]. 上海交通大学学报, 2014, 48(11): 1503-1509.
[13]	张校通, 郑东生, 夏小和, 王建华. 非均质饱和-非饱和土的变形特性[J]. 上海交通大学学报, 2012, 46(10): 1544-1547.
[14]	刘映晶, 陈锦剑, 王建华, 徐安军, 李耀良. 逆作土方通道式施工工艺对基坑变形的影响[J]. 上海交通大学学报, 2012, 46(01): 89-93.
[15]	向国威, 叶冠林, 王建华. 冲击振动在砂土中传播的模型试验研究[J]. 上海交通大学学报, 2012, 46(01): 136-141.