Fracture analysis employing equivalent domain integral method and nodal integration techniques based on reproducing kernel particle method

Springer Science and Business Media LLC - Tập 9 - Trang 1265-1278 - 2022
Satoyuki Tanaka1, Akihiro Takata1, Ming-Jyun Dai1, Hanlin Wang1, Shota Sadamoto1,2
1Graduate School of Advanced Science and Engineering, Hiroshima University, Higashi-Hiroshima, Japan
2Fujitsu Limited, Kawasaki city, Japan

Tóm tắt

A novel technique to evaluate fracture mechanics parameters is investigated employing the equivalent domain integral (EDI) method and nodal integration (NI) techniques. Galerkin-based meshfree method is adopted. Reproducing kernel (RK) is chosen for the meshfree interpolant. Stabilized conforming nodal integration (SCNI) and sub-domain stabilized conforming integration (SSCI) are adopted for numerically integrating the stiffness matrix. Voronoi diagram is employed to compute volume of each NI domain. The EDI method is addressed to evaluate the fracture mechanics parameters, i.e., energy release rate and stress intensity factors (SIFs). Because the displacement and its derivatives are computed based on SCNI/SSCI, the EDI can be discretized by summing up the physical quantities and volume of each cell/sub-cell. No special quadrature rule is required. To separate the energy release rate into the mixed-mode SIFs, interaction integral method is chosen. Efficient and accurate fracture parameter computation is achieved. Some numerical examples are demonstrated for mixed-mode fracture parameter evaluation and crack propagation analysis. Accuracy and effectiveness of the presented approach are studied.

Tài liệu tham khảo

Kishimoto K, Aoki S, Sakata M (1980) Dynamic stress intensity factors using \(\hat{J}\)-integral and finite element method. Eng Fract Mech 13(2):387–394

Fleming M, Chu YA, Moran B, Belytschko T (1997) Enriched element-free galerkin methods for crack tip fields. Int J Numer Meth Eng 40:1483–1504

Ansys Mechanical APDL Verification Manual (2021), VM267: Inclined cracked in 2-D plate under uniform tension loading. 757-758

MSC software (2014) MSC.Marc2014 User’s manual

https://www.code-aster.org/spip.php?rubrique2 [accessed 21.11.21]

Geniaut S, eXtended finite element method: general information, Code_Aster Reference manual R7.02.12, 2017

Geniaut S, Massin P, Mo\(\ddot{{\rm e}}\)s N (2005) Evaluation of stress intensity factors with G-theta method and level sets in \(Code_Aster\). 11th Int Conf of Fract, Turin, Italy