冲击载荷下平板玻璃裂纹扩展和破坏形态的数值模拟

王木飞; 李志强

doi:10.11858/gywlxb.20220558

冲击载荷下平板玻璃裂纹扩展和破坏形态的数值模拟

doi: 10.11858/gywlxb.20220558

王木飞^{1, 2, 3,},
李志强^{1, 2, 3,*, ,}

1.
太原理工大学机械与运载工程学院应用力学研究所, 山西太原　030024
2.
山西省材料强度与结构冲击重点实验室, 山西太原　030024
3.
力学国家级实验教学示范中心, 山西太原　030024

基金项目: 国家自然科学基金（11972244）

详细信息

作者简介:
王木飞（1996-），男，硕士研究生，主要从事冲击动力学和近场动力学研究.E-mail：1470979047@qq.com

通讯作者:
李志强（1973-），男，博士，教授，主要从事冲击动力学研究. E-mail：lizhiqiang@tyut.edu.cn

中图分类号: O346.1; O521.2
计量
- 文章访问数: 315
- HTML全文浏览量: 224
- PDF下载量: 53
出版历程
- 收稿日期: 2022-04-04
- 修回日期: 2022-05-27
- 网络出版日期: 2022-09-24
- 刊出日期: 2022-10-11

Numerical Simulation of Crack Propagation and Damage Behavior of Glass Plates under Impact Loading

WANG Mufei^{1, 2, 3
,},
LI Zhiqiang^{1, 2, 3,*
, ,}

1.
Institute of Applied Mechanics, College of Mechanical and Vehicle Engineering, Taiyuan University of Technology, Taiyuan 030024, Shanxi, China
2.
Shanxi Key Laboratory of Material Strength and Structural Impact, Taiyuan 030024, Shanxi, China
3.
National Demonstration Center for Experimental Mechanics Education, Taiyuan 030024, Shanxi, China

摘要

摘要: 裂纹的萌生和扩展是计算力学领域长期存在的难点和热点，也是玻璃、岩石、混凝土等脆性材料中常见的工程实际问题。为了探究冲击载荷下平板钠钙玻璃的损伤破坏行为和细观裂纹扩展规律，分别采用单元删除法、不连续伽辽金近场动力学法（discontinuous Galerkin peridynamic, DG-PD）和无网格粒子近场动力学法（meshless peridynamic, M-PD），研究其裂纹扩展行为。单元删除法采用JH-2材料模型，同时添加最大主应力和最大主应变失效准则；DG-PD法采用节点分离操作，并施加临界能量释放率准则；M-PD法采用自编程序粒子离散方法，选择合适的计算域，并施加临界伸长率准则。模拟结果表明：(1) 单元删除法可大致模拟出玻璃在冲击载荷下的损伤形貌，但在捕捉裂纹分叉和贯通等方面略显不足，未见玻璃碎片的飞溅，无法通过碎片飞溅速度评估其安全性能；(2) DG-PD法中环状裂纹和径向裂纹明显，裂纹具有很高的对称性，冲击点和边框处有大量玻璃碎片飞溅；(3) M-PD法中能捕捉到径向裂纹和环向裂纹，且裂纹的对称性较好，近场域和冲击速度对平板玻璃的动态响应有着重要的影响，就损伤形态而言，M-PD法和DG-PD法具有很高的一致性。
- 冲击荷载 /
- 平板玻璃 /
- 裂纹扩展 /
- 单元删除法 /
- 近场动力学 /
- 无网格法
Abstract: Crack initiation and propagation is a long-standing difficult problem in solid mechanics, especially for elastic-brittle material. To explore the damage and crack propagation behavior of glass plates under impact loading, the element deletion, discontinuous Galerkin peridynamic (DG-PD), and meshless peridynamic (M-PD) methods are used to conduct numerical simulations, respectively. The JH-2 material model, and the maximum principal stress and maximum principal strain failure criteria are adopted in the element deletion method. The node separation operation and the critical energy release rate criterion are used in the DG-PD method. In the M-PD method, a self-programmed particle discretization method is utilized along with an appropriate computational domain, and a critical elongation criterion is imposed. The simulation results show that: (1) the element deletion method can roughly simulate the damage morphology of glass under impact loading, but it is insufficient in capturing crack bifurcation and penetration. (2) In the DG-PD method, circumferential cracks and radial cracks are observed, and the cracks are of high symmetry. In addition, there are a lot of glass fragments splashing at the impact point and the frame. (3) Radial cracks and circumferential cracks can be captured in the M-PD method, and the symmetry of the cracks is good. The size of horizon and the impact velocity show great influence on the dynamic responses of the glass plates. As far as the damage form is concerned, the M-PD method and the DG-PD method yield consistent results.
- impact loading /
- glass plates /
- crack propagation /
- element deletion method /
- peridynamic /
- meshless method

HTML全文

图 1 有限元模型及弹头放大图

Figure 1. Finite element model and enlarged view of the bullet

下载: 全尺寸图片幻灯片

图 2 迎撞面和背撞面的损伤演化

Figure 2. Damage evolution of the top face and bottom face

下载: 全尺寸图片幻灯片

图 3 不同冲击速度下平板玻璃的最终损伤模态

Figure 3. Final damage of glass plates at different impact velocities

下载: 全尺寸图片幻灯片

图 4 不同冲击速度下损伤参数的对比

Figure 4. Comparison of damage parameters at different impact velocities

下载: 全尺寸图片幻灯片

图 5 二维有限单元断裂前后的节点分离示意图

Figure 5. Node separation before and after fracture of 2D finite elements

下载: 全尺寸图片幻灯片

图 6 5 m/s冲击速度下迎撞面的损伤演化

Figure 6. Damage evolution of the impact face at the impact velocity of 5 m/s

下载: 全尺寸图片幻灯片

图 9 50 m/s冲击速度下迎撞面的损伤演化

Figure 9. Damage evolution of the impact face at the impact velocity of 50 m/s

下载: 全尺寸图片幻灯片

图 7 15 m/s冲击速度下迎撞面的损伤演化

Figure 7. Damage evolution of the top face at the impact velocity of 15 m/s

下载: 全尺寸图片幻灯片

图 8 30 m/s冲击速度下迎撞面的损伤演化

Figure 8. Damage evolution of the top face at the impact velocity of 30 m/s

下载: 全尺寸图片幻灯片

图 10 采用DG-PD法得到的不同冲击速度下中心点单元的速度变化

Figure 10. Velocity change of center element obtained by DG-PD method at different impact velocities

下载: 全尺寸图片幻灯片

图 11 M-PD非局部相互作用示意图

Figure 11. Schematic diagram of M-PD non-local interaction

下载: 全尺寸图片幻灯片

图 12 M-PD粒子模型及弹头放大图

Figure 12. M-PD particle model and enlarged view of the bullet

下载: 全尺寸图片幻灯片

图 13 30 m/s冲击速度下迎爆面的损伤演化

Figure 13. Damage evolution images of the top face at the impact velocity of 30 m/s

下载: 全尺寸图片幻灯片

图 14 不同近场域大小下的平板玻璃损伤

Figure 14. Damage of glass plate under different horizon sizes

下载: 全尺寸图片幻灯片

图 15 不同冲击速度下平板玻璃的最终损伤破坏

Figure 15. Final damage of glass plate at different impact velocities

下载: 全尺寸图片幻灯片

表 1 玻璃材料JH-2材料模型参数^[6]

Table 1. JH-2 material model parameters of glass material^[6]

$\;\rho$ /(kg·m⁻³) G/GPa A B C M N $\dot{\varepsilon }_{\rm s}$ T/GPa

2530 24 0.93 0.2 0.003 1 0.77 1 0.05042

σ_f,max σ_e/GPa p_HEL/GPa β D₁ D₂ k₁/GPa k₂/GPa k₃/GPa
0.5 5.95 2.92 1 0.043 0.85 45.4 −138 290

下载: 导出CSV

表 2 边框和弹头刚体材料模型参数

Table 2. Rigid body material model parameters of frame and bullet

$\,\rho$ /(kg·m⁻³) $\,\nu$ E/GPa

7850 0.28 200

下载: 导出CSV

表 3 玻璃材料DG-PD材料模型参数

Table 3. DG-PD material model parameters of glass material

$\,\rho$ /(kg·m⁻³) E/GPa G_t/(J·m⁻²) HSFAC

2530 72 15.47 0.8

下载: 导出CSV

表 4 玻璃材料M-PD材料模型参数

Table 4. M-PD model parameters of glass

$\,\rho$ /(kg·m⁻³) E/GPa G/GPa S₀ $\Delta x$ /mm $\delta$ /mm

2530 72 24 $7.17\times 1{0}^{-5}$ 3 9.03

下载: 导出CSV

参考文献(18)

[1]	KRAUTHAMMER T, ALTENBERG A. Negative phase blast effects on glass panels [J]. International Journal of Impact Engineering, 2000, 24: 1–17. doi: 10.1016/S0734-743X(99)00043-3
[2]	YUAN Y, TAN P J, LI Y B. Dynamic structural response of laminated glass panels to blast loading [J]. Composite Structure, 2017, 182: 579–589. doi: 10.1016/j.compstruct.2017.09.028
[3]	HOLMQUIST T J, JOHNSON G R, LOPATIN C M, et al. High strain rate properties and constitutive modeling of glass: SAND 95-0379C [R]. Albuquerque, NM: Sandia National Laboratories, 1995.
[4]	DARYADEL S S, MANTENA P R, KIM K, et al. Dynamic response of glass under low-velocity impact and high strain-rate SHPB compression loading [J]. Journal of Non-Crystalline Solids, 2016, 432B: 432–439.
[5]	张晓颖, 李胜杰, 李志强. 爆炸载荷作用下夹层玻璃动态响应的数值模拟 [J]. 兵工学报, 2018, 39(7): 1379–1388. doi: 10.3969/j.issn.1000-1093.2018.07.016 ZHANG X Y, LI S J, LI Z Q. Numerical simulation of dynamic response of laminated glass subjected to blast load [J]. Acta Armamentarii, 2018, 39(7): 1379–1388. doi: 10.3969/j.issn.1000-1093.2018.07.016
[6]	唐剑. 冲击荷载下建筑玻璃裂纹扩展及破碎的试验与数值模拟研究 [D]. 天津: 天津大学, 2018.
[7]	DAUX C, MOES N, DOLBOW J, et al. Arbitrary branched and intersecting cracks with the extended finite element method [J]. International Journal for Numerical Methods in Engineering, 2000, 48(12): 1741−1760.
[8]	SILLING S A. Reformulation of elasticity theory for discontinuities and long-range forces [J]. Journal of Mechanics Physics of Solids, 2000, 48(1): 175–209. doi: 10.1016/S0022-5096(99)00029-0
[9]	徐泳, 孙其诚, 张凌, 等. 颗粒离散元法研究进展 [J]. 力学进展, 2003, 33(2): 251–260. doi: 10.3321/j.issn:1000-0992.2003.02.010 XU Y, SUN Q C, ZHANG L, et al. Advances in discrete element methods for particulate materials [J]. Advances in Mechanics, 2003, 33(2): 251–260. doi: 10.3321/j.issn:1000-0992.2003.02.010
[10]	张锁春. 光滑质点流体动力学(SPH)方法(综述) [J]. 计算物理, 1996, 13(4): 2–14. ZHANG S C. Smoothed particle hydrodynamics (SPH) method (a review) [J]. Chinese Journal of Computation Physics, 1996, 13(4): 2–14.
[11]	SUN X, LIU W. Modeling and characterization of dynamic failure of borosilicate glass under compression/shear loading [J]. International Journal of Impact Engineering, 2009, 36(2): 226–234. doi: 10.1016/j.ijimpeng.2008.01.014
[12]	NIE X, CHEN W. Dynamic failure of borosilicate glass under compression/shear loading experiments [J]. Journal of the American Ceramic Society, 2007, 90(8): 2556–2562. doi: 10.1111/j.1551-2916.2007.01819.x
[13]	ZANG M, CHEN H, LEI Z. Simulation on high velocity impact process of windshield by sph/fem coupling method [C]//2010 WASE International Conference on Information Engineering. Institute of Electrical and Electronics Engineers, 2010: 381–384.
[14]	DENG R B, JIN X L. Numerical simulation for blast analysis of insulating glass in a curtain wall [J]. International Journal for Computational Methods in Engineering Science and Mechanics, 2010, 11(3): 162–171. doi: 10.1080/15502281003702302
[15]	BOBARU F, HA Y D, HU W K, et al. Damage progression from impact in layered glass modeled with peridynamics [J]. Central European Journal of Engineering, 2012, 2(4): 551–561.
[16]	BLESS S, CHEN T. Impact damage in layered glass [J]. International Journal of Fracture, 2010, 162(1/2): 151–158. doi: 10.1007/s10704-009-9379-7
[17]	REN B, WU C T, SELESON P, et al. A peridynamic failure analysis of fiber-reinforced composite laminates using finite element discontinuous Galerkin approximations [J]. International Journal of Fracture, 2018, 214: 49–68. doi: https://doi.org/10.1007/s10704-018-0317-4
[18]	师燕超, 于海涛, 李忠献. 爆炸荷载作用下建筑玻璃的破碎分析 [J]. 天津大学学报(自然科学与工程技术版), 2015, 48(10): 888–893. doi: 10.11784/tdxbz201411025 SHI Y C, YU H T, LI Z X. Analysis of architecture glass breaking under explosion load [J]. Journal of Tianjin University (Science and Technology), 2015, 48(10): 888–893. doi: 10.11784/tdxbz201411025

施引文献

资源附件(0)

访问统计

点击查看大图

图(15) / 表(4)

计量

文章访问数: 315
HTML全文浏览量: 224
PDF下载量: 53

$\;\rho$ /(kg·m⁻³)	G/GPa	A	B	C	M	N	$\dot{\varepsilon }_{\rm s}$	T/GPa
2530	24	0.93	0.2	0.003	1	0.77	1	0.05042

σ_f,max	σ_e/GPa	p_HEL/GPa	β	D₁	D₂	k₁/GPa	k₂/GPa	k₃/GPa
0.5	5.95	2.92	1	0.043	0.85	45.4	−138	290

$\,\rho$ /(kg·m⁻³)	E/GPa	G_t/(J·m⁻²)	HSFAC
2530	72	15.47	0.8

$\,\rho$ /(kg·m⁻³)	$\,\nu$	E/GPa
7850	0.28	200

$\,\rho$ /(kg·m⁻³)	E/GPa	G/GPa	S₀	$\Delta x$ /mm	$\delta$ /mm
2530	72	24	$7.17\times 1{0}^{-5}$	3	9.03

冲击载荷下平板玻璃裂纹扩展和破坏形态的数值模拟

doi: 10.11858/gywlxb.20220558

作者简介: 王木飞（1996-），男，硕士研究生，主要从事冲击动力学和近场动力学研究.E-mail：1470979047@qq.com

通讯作者: 李志强（1973-），男，博士，教授，主要从事冲击动力学研究. E-mail：lizhiqiang@tyut.edu.cn

计量

出版历程

Numerical Simulation of Crack Propagation and Damage Behavior of Glass Plates under Impact Loading

计量

出版历程

目录

作者简介:
王木飞（1996-），男，硕士研究生，主要从事冲击动力学和近场动力学研究.E-mail：1470979047@qq.com

通讯作者:
李志强（1973-），男，博士，教授，主要从事冲击动力学研究. E-mail：lizhiqiang@tyut.edu.cn