工程力学 ›› 2016, Vol. 33 ›› Issue (8): 8-14.doi: 10.6052/j.issn.1000-4750.2015.01.0003

• 基本方法 • 上一篇    下一篇

基于水平集算法的扩展比例边界有限元法研究

李建波1, 陈白斌1,2, 林皋1   

  1. 1. 大连理工大学工程抗震研究所, 辽宁, 大连 116024;
    2. 中铁第四勘察设计院集团有限公司, 湖北, 武汉 430063
  • 收稿日期:2015-01-04 修回日期:2015-03-10 出版日期:2016-08-25 发布日期:2016-08-25
  • 通讯作者: 陈白斌(1991-),男,湖南人,博士生,主要从事工程安全与评价研究(E-mail:cbbdlut@126.com). E-mail:cbbdlut@126.com
  • 作者简介:李建波(1977-),男,河北人,副教授,博士,主要从事核电与水工结构抗震安全评价研究(E-mail:jianboli@dlut.edu.cn);林皋(1929-),男,江西人,教授,博导,中科院院士,主要从事核电与水工结构抗震研究(E-mail:gaolin@dlut.edu.cn).
  • 基金资助:

    国家自然科学基金重点项目(51138001);国家重大科技专项项目(2011ZX06002-10);中央高校基本科研业务费专项资金项目(DUT13LK16)

STUDY ON THE EXTENED SCALED BOUNDARY FINITE ELEMENT METHOD BASED ON LEVEL SET ALGORITHM

LI Jian-bo1, CHEN Bai-bin1,2, LIN Gao1   

  1. 1. Earthquake Engineering Research Division, Dalian University of Technology, Dalian, Liaoning 116024, China;
    2. China Railway Siyuan Survey and Design Group Co. Ltd, Wuhan, Hubei 430063, China
  • Received:2015-01-04 Revised:2015-03-10 Online:2016-08-25 Published:2016-08-25
  • Contact: 10.6052/j.issn.1000-4750.2015.01.0003 E-mail:cbbdlut@126.com

摘要:

扩展比例边界有限元法在裂纹贯穿单元采用Heaviside阶跃函数描述裂纹面两侧的不连续位移,在裂尖则采用半解析的比例边界有限元描述奇异应力场。该方法具有无需预先知道裂尖渐进场的形式,无需采用特殊的数值积分技术直接生成裂尖刚度阵,对多种应力奇异类型可根据定义直接求解广义应力强度因子的特点。该文将扩展比例边界有限元法与水平集方法相结合,进一步发展了扩展比例边界有限元法,并将其应用于解决裂纹扩展的问题。在数值算例中,通过编写完整的MATLAB分析计算程序,求解了单边缺口的三点弯曲梁和四点剪切梁的裂纹扩展问题,计算结果显示扩展比例边界有限元法能有效地预测裂纹轨迹和荷载-位移曲线。通过参数敏感性分析,还可得出该方法具有较低的网格依赖性,且对裂纹扩展步长不敏感。

关键词: 扩展比例边界有限元法, 线弹性断裂力学, 裂纹扩展, 应力强度因子, 水平集法

Abstract:

In the extended scaled boundary finite element method (X-SBFEM), the Heaviside enrichment is used to represent the jump across the discontinuity surface in split element while the nonsmooth behavior around the crack tip is described by the semi-analytical SBFEM. This method does not require a priori knowledge of the asymptotic fields. It also does not require special numerical integration technique so that the stiffness of the region containing the crack tip is computed directly, while the generalized stress intensity factors of many kinds of singularities could be obtained directly from definition by the consistent formula. In this paper, the X-SBFEM is combined with the level set method so that it can be applied to solve the problem of crack propagation. During its numerical evaluation, the three-point bending specimen and the four-point shear specimen is modeled by a complete MATLAB code. The results show that the X-SBFEM is capable of predicting crack trajectories and load-displacement relations accurately. Based on the sensitivity analysis of parameters, it was found that different mesh densities and crack propagation lengths led to consistent results.

Key words: extended scale boundary finite element method, linear elastic fracture mechanics, crack propagation, stress intensity factors, level set method

中图分类号: 

  • O342

[1] Yang Z J, Chen J F. Finite element modelling of multiple discrete cohesive crack propagation in reinforced concrete beams[J]. Engineering Fracture Mechanics, 2005, 72(14):2280-97.
[2] Belytschko T, Gracie R. On XFEM applications to dislocations and interfaces[J]. International Journal of Plasticity, 2007, 23(10):1721-1738.
[3] 张晓东, 丁勇, 任旭春. 混凝土裂纹扩展过程模拟的扩展有限元法研究[J]. 工程力学, 2013, 30(7):14-21. Zhang Xiaodong, Ding Yong, Ren Xuchun. Simulation of the concrete crack propagation process with the extend finite element method[J]. Engineering Mechanics, 2013, 30(7):14-21. (in Chinese)
[4] Sundararajan N, Song C M. Representation of singular fields without asymptotic enrichment in the extended finite element method[J]. International Journal for Numerical Methods in Engineering, 2013, 96(13):813-841.
[5] Belytschko T, Gracie R, Ventura G. A review of extended/generalized finite element methods for material model[J]. Modelling and Simulation in Materials Science and Engineering, 2009, 17(4):1-24.
[6] 刘钧玉, 林皋, 胡志强. 裂纹面荷载作用下多裂纹应力强度因子计算[J]. 工程力学, 2011, 28(4):7-12. Liu Junyu, Lin Gao, Hu Zhiqiang. The calculation of stress intensity factors of multiple cracks under surface tractions[J]. Engineering Mechanics, 2011, 28(4):7-12. (in Chinese)
[7] Song C M, Francis T, Wei G. A definition and evaluation procedure of generalized stress intensity factors at cracks and multi-material wedges[J]. Engineering Fracture Mechanics, 2010, 10(1):2316-2336.
[8] 陈白斌, 李建波, 林皋. 基于X-SBFEM的裂纹体非网格重剖分耦合模型研究[J]. 工程力学, 2015, 32(3):15-21. Chen Baibin, Li Jianbo, Lin Gao. Study on the coupling model of crack without remeshing based on X-SBFEM[J]. Engineering Mechanics, 2015, 32(3):15-21. (in Chinese)
[9] 茹忠亮, 朱传锐, 赵洪波. 裂纹扩展问题的改进XFEM算法[J]. 工程力学, 2012, 29(7):12-16. Ru Zhongliang, Zhu Chuanrui, Zhao Hongbo. An improved algorithm of XFEM for the crack propagation problems[J]. Engineering Mechanics, 2012, 29(7):12-16. (in Chinese)
[10] Deeks A J, Wolf J P. A virtual work derivation of the scaled boundary finite-element method for elastostatics[J]. Computational Mechanics, 2002, 28(6):489-504.
[11] Yang Z J. Fully automatic modelling of mixed-mode crack propagation using scaled boundary finite element method[J]. Engineering Fracture Mechanics, 2006, 73(12):1711-1731.
[12] Stolarska M, Chopp D L, Moes N, Belytschko T. 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.
[13] Petersson P E. Crack growth and development of fracture zone in plain concrete and similar materials, Report TVBM-1006[R]. Lund:Lund Institute of Technology, 1981.
[14] Yang Z J, Deeks A J. Modelling cohesive crack growth using a two-step finite element-scaled boundary finite element coupled method[J]. International Journal of Fracture. 2007, 143(4):333-354.
[15] Arrea M, Ingraffea A R. Mixed-mode crack propagation in mortar and concrete[R]. America, State of New York:Department of Structural Engineering, Cornell University, 1982:Report No. 81-13.
[16] Huang R, Sukumar N, Prevost J H. Modeling quasi-static crack growth with the extended finite element method Part II:Numerical applications[J]. International Journal Solids Structure, 2003, 40(26):7539-7552.

[1] 胡光辉, 徐涛, 陈崇枫, 杨学凯. 基于离散元法的脆性岩石细观蠕变失稳研究[J]. 工程力学, 2018, 35(9): 26-36.
[2] 王综轶, ZHANG Yun-feng, 王元清, 杜新喜, 袁焕鑫. 交通信号支撑结构疲劳裂纹扩展有限元分析[J]. 工程力学, 2018, 35(9): 180-187.
[3] 王珊. 基于Kirchhoff板中裂纹尖端辛本征解的有限元应力恢复方法[J]. 工程力学, 2018, 35(5): 10-16.
[4] 白凡, 林原胜, 代路, 柳勇, 肖颀. 船舶动力装置高温高压管道破前漏分析[J]. 工程力学, 2018, 35(3): 242-248.
[5] 牛彦泽, 徐业鹏, 黄丹. 双轴动载作用下脆性裂纹扩展问题的近场动力学建模与分析[J]. 工程力学, 2018, 35(10): 249-256.
[6] 廖小伟, 王元清, 石永久, 陈宏. 低温环境下桥梁钢Q345qD疲劳裂纹扩展行为研究[J]. 工程力学, 2018, 35(10): 85-91.
[7] 关鹏涛, 闾川阳, 唐夏焘, 朱建伟, 丁振宇, 高增梁. Q345R钢多个过载作用下疲劳裂纹扩展行为研究[J]. 工程力学, 2017, 34(7): 224-231,240.
[8] 郁大照, 陈跃良, 王允良. 含多处损伤宽板螺接搭接件疲劳寿命研究[J]. 工程力学, 2017, 34(6): 217-225.
[9] 秦洪远, 黄丹, 刘一鸣, 章青. 基于改进型近场动力学方法的多裂纹扩展分析[J]. 工程力学, 2017, 34(12): 31-38.
[10] 李戎, 梁斌, NODA Nao-Aki. 基于临界应力强度因子的胶接对接接头抗拉强度简便预测方法[J]. 工程力学, 2017, 34(11): 218-224.
[11] 黄凯, 尹函彬, 孙庆平. 红外热像技术测量I型裂纹应力强度因子的研究[J]. 工程力学, 2017, 34(11): 231-239.
[12] 李伟, 郑宏, 郭宏伟. 基于MLS-NMM的摩擦接触问题研究[J]. 工程力学, 2017, 34(11): 18-25.
[13] 刘学聪, 章青, 夏晓舟. 一种新型裂尖加强函数的显式动态扩展有限元法[J]. 工程力学, 2017, 34(10): 10-18.
[14] 傅宇光, 童乐为, 刘博. 基于Beach Marking方法的钢结构疲劳裂纹检测研究[J]. 工程力学, 2016, 33(8): 93-100,131.
[15] 宗亮, 施刚, 王元清, 廖小伟, 于兴哲. WNQ570桥梁钢及其对接焊缝疲劳裂纹扩展性能试验研究[J]. 工程力学, 2016, 33(8): 45-51.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 原 园;徐颖强;吕国志;朱贤飞. 齿轮啮合过程中安定状态残余应力的数值方法研究[J]. 工程力学, 2008, 25(10): 0 -211, .
[2] 邢德进;李忠献. 应用SMA智能阻尼器的结构模糊控制[J]. 工程力学, 2008, 25(10): 0 -228, .
[3] 周小平;杨海清;张永兴. 有限宽偏心裂纹板在裂纹面受两对集中拉力作用时裂纹线的弹塑性解析解[J]. 工程力学, 2008, 25(1): 0 -027 .
[4] 龚耀清;包世华. 超高层建筑空间巨型框架自由振动计算的新方法[J]. 工程力学, 2008, 25(10): 0 -140 .
[5] 刘金兴;邓守春;张 晶;梁乃刚. 颗粒复合材料断裂的梁网格模型[J]. 工程力学, 2008, 25(10): 0 -037 .
[6] 郎风超;邢永明;朱 静. 应用纳米压痕技术研究表面纳米化后316L 不锈钢力学性能[J]. 工程力学, 2008, 25(10): 0 -071 .
[7] 郭小刚;刘人怀;曾 娜;金 星. 子结构位移迭代法修正软管空间形态[J]. 工程力学, 2008, 25(10): 0 -032 .
[8] 邢静忠;柳春图. 线弹性土壤中埋设悬跨管道的屈曲分析[J]. 工程力学, 2008, 25(10): 0 -075 .
[9] 刘祥庆;刘晶波. 基于纤维模型的拱形断面地铁车站结构弹塑性地震反应时程分析[J]. 工程力学, 2008, 25(10): 0 -157 .
[10] 郝庆多;王言磊;侯吉林;欧进萍;. GFRP带肋筋粘结性能试验研究[J]. 工程力学, 2008, 25(10): 0 -165, .
X

近日,本刊多次接到来电,称有不法网站冒充《工程力学》杂志官网,并向投稿人收取高额费用。在此,我们郑重申明:

1.《工程力学》官方网站是本刊唯一的投稿渠道(原网站已停用),《工程力学》所有刊载论文必须经本刊官方网站的在线投稿审稿系统完成评审。我们不接受邮件投稿,也不通过任何中介或编辑收费组稿。

2.《工程力学》在稿件符合投稿条件并接收后会发出接收通知,请作者在接到版面费或审稿费通知时,仔细检查收款人是否为“《工程力学》杂志社”,千万不要汇款给任何的个人账号。请广大读者、作者相互转告,广为宣传!如有疑问,请来电咨询:010-62788648。

感谢大家多年来对《工程力学》的支持与厚爱,欢迎继续关注我们!

《工程力学》杂志社

2018年11月15日