J-integral evaluation for 2D mixed-mode crack problems employing a meshfree stabilized conforming nodal integration method

Published on Aug 1, 2016in Computational Mechanics3.459
· DOI :10.1007/S00466-016-1288-9
Satoyuki Tanaka22
Estimated H-index: 22
(Hiroshima University),
Hirotaka Suzuki2
Estimated H-index: 2
(Hiroshima University)
+ 3 AuthorsTinh Quoc Bui54
Estimated H-index: 54
(TITech: Tokyo Institute of Technology)
Two-dimensional (2D) in-plane mixed-mode fracture mechanics problems are analyzed employing an efficient meshfree Galerkin method based on stabilized conforming nodal integration (SCNI). In this setting, the reproducing kernel function as meshfree interpolant is taken, while employing the SCNI for numerical integration of stiffness matrix in the Galerkin formulation. The strain components are smoothed and stabilized employing Gauss divergence theorem. The path-independent integral (J-integral) is solved based on the nodal integration by summing the smoothed physical quantities and the segments of the contour integrals. In addition, mixed-mode stress intensity factors (SIFs) are extracted from the J-integral by decomposing the displacement and stress fields into symmetric and antisymmetric parts. The advantages and features of the present formulation and discretization in evaluation of the J-integral of in-plane 2D fracture problems are demonstrated through several representative numerical examples. The mixed-mode SIFs are evaluated and compared with reference solutions. The obtained results reveal high accuracy and good performance of the proposed meshfree method in the analysis of 2D fracture problems.
📖 Papers frequently viewed together
65 Citations
36 Citations
4,430 Citations
Abstract Accurate numerical modeling of multifield piezoelectric materials is challenging because of the inherent electro-mechanical coupling effect and material anisotropic behaviors. The modeling becomes even more difficult especially for problems with non-smooth solutions like crack under dynamic loading. We present in this paper an extension of the extended isogeometric analysis (XIGA) for simulation of two-dimensional fracture mechanics problems in piezoelectric materials under dynamic and ...
124 CitationsSource
#1Zuoyi Kang (TITech: Tokyo Institute of Technology)H-Index: 3
#2Tinh Quoc Bui (TITech: Tokyo Institute of Technology)H-Index: 54
Last. Sohichi Hirose (TITech: Tokyo Institute of Technology)H-Index: 15
view all 5 authors...
A novel extended 4-node quadrilateral finite element (XCQ4) based on a consecutive-interpolation procedure (CIP) with continuous nodal stress for accurately modeling singular stress fields near crack tips of two-dimensional (2D) cracks in solids is presented. In contrast to the traditional approaches, the approximation functions constructed based on the CIP involve both nodal values and averaged nodal gradients as interpolation conditions. Our objective is to exhibit a pioneering extension of th...
64 CitationsSource
#1Satoyuki Tanaka (Hiroshima University)H-Index: 22
#2Shogo Sannomaru (Hiroshima University)H-Index: 4
Last. Hiroshi Okada (Ridai: Tokyo University of Science)H-Index: 18
view all 6 authors...
Abstract Two-dimensional (2D) dynamic stress concentration problems are analyzed using the wavelet Galerkin method (WGM). Linear B-spline scaling/wavelet functions are employed. We introduce enrichment functions for the X-FEM to represent a crack geometry. In the WGM, low-resolution scaling functions are periodically located across the entire analysis domain to approximate deformations of a body. High-resolution wavelet functions and enrichment functions including crack tip singular fields are s...
18 CitationsSource
#1Satoyuki Tanaka (Hiroshima University)H-Index: 22
#2Hirotaka Suzuki (Hiroshima University)H-Index: 2
Last. Tinh Quoc Bui (TITech: Tokyo Institute of Technology)H-Index: 54
view all 5 authors...
Abstract A novel technique to evaluate the path-independent integral using nodal integration is developed for cracked plate problems. The formulation is based on the Mindlin–Reissner plate theory, and the reproducing kernel is used as a meshfree interpolant. A visibility criterion, diffraction method, and enriched basis are included. To integrate the stiffness matrix, the stabilized conforming nodal integration and the sub-domain conforming integration are adopted. A moment intensity factor is a...
65 CitationsSource
#1Satoyuki Tanaka (Hiroshima University)H-Index: 22
#2H. Suzuki (Hiroshima University)H-Index: 1
Last. Shogo Sannomaru (Hiroshima University)H-Index: 4
view all 4 authors...
Two-dimensional (2D) crack problems are solved employing a novel technique based on a combination of wavelet Galerkin method and X-FEM with a high-order interpolant. Multiresolution analysis of the wavelet basis functions (scaling/wavelet functions) plays an important role in the numerical simulation. High-order B-spline scaling/wavelet functions are chosen as the basis functions. Severe stress concentration near a crack tip is represented by superposing the multiresolution wavelet functions. In...
22 CitationsSource
#1Nha Thanh Nguyen (VNUHCM-UT: Ho Chi Minh City University of Technology)H-Index: 5
#2Tinh Quoc Bui (University of Siegen)H-Index: 54
Last. Thien Tich Truong (VNUHCM-UT: Ho Chi Minh City University of Technology)H-Index: 7
view all 4 authors...
Abstract We present a new approach based on local partition of unity extended meshfree Galerkin method for modeling quasi-static crack growth in two-dimensional (2D) elastic solids. The approach utilizing the local partition of unity as a priori knowledge on the solutions of the boundary value problems that can be added into the approximation spaces of the numerical solutions. It thus allows for extending the standard basis functions by enriching the asymptotic near crack-tip fields to accuratel...
78 CitationsSource
#1Shota Sadamoto (Hiroshima University)H-Index: 10
#2Satoyuki Tanaka (Hiroshima University)H-Index: 22
Last. Shigenobu Okazawa (Hiroshima University)H-Index: 9
view all 3 authors...
A meshfree approach for plate buckling/post-buckling problems in the case of uniaxial thrust is presented. A geometrical nonlinear formulation is employed using reproducing kernel approximation and stabilized conforming nodal integration. The bending components are represented by Mindlin---Reissner plate theory. The formulation has a locking-free property in imposing the Kirchhoff mode reproducing condition. In addition, in-plane deformation components are approximated by reproducing kernels. Th...
21 CitationsSource
#1Jan Sladek (SAV: Slovak Academy of Sciences)H-Index: 53
#2V. Sladek (SAV: Slovak Academy of Sciences)H-Index: 13
Last. M. Wünsche (University of Siegen)H-Index: 8
view all 5 authors...
A meshless method based on the local Petrov-Galerkin approach is proposed to solve initial-boundary-value crack problems in decagonal quasicrystals. These quasicrystals belong to the class of two-dimensional (2-d) quasicrystals, where the atomic arrangement is quasiperiodic in a plane, and periodic in the perpendicular direction. The ten-fold symmetries occur in these quasicrystals. The 2-d crack problem is described by a coupling of phonon and phason displacements. Both stationary governing equ...
21 CitationsSource
#1Satoyuki Tanaka (Hiroshima University)H-Index: 22
#2Hiroshi Okada (Ridai: Tokyo University of Science)H-Index: 18
Last. Masahiko Fujikubo (Osaka University)H-Index: 14
view all 4 authors...
SUMMARY This paper presents fracture mechanics analysis using the wavelet Galerkin method and extended finite element method. The wavelet Galerkin method is a new methodology to solve partial differential equations where scaling/wavelet functions are used as basis functions. In solid/structural analyses, the analysis domain is divided into equally spaced structured cells and scaling functions are periodically placed throughout the domain. To improve accuracy, wavelet functions are superposed on ...
33 CitationsSource
#1Peng Liu (Hohai University)H-Index: 12
#2Tiantang Yu (Hohai University)H-Index: 32
Last. Chuanzeng Zhang (University of Siegen)H-Index: 38
view all 4 authors...
Abstract We study transient dynamic fracture behaviors of stationary cracked functionally graded piezoelectric materials (FGPMs) under impact loading by using the extended finite element method (X-FEM). The material properties are assumed to be varied exponentially along one direction. A dynamic X-FEM model associated with the stable implicit integration technique is developed to serve that purpose, while the contour interaction integral technique is employed to accurately evaluate the relevant ...
58 CitationsSource
Cited By46
#1Hanlin Wang (Hiroshima University)H-Index: 5
#2Satoyuki Tanaka (Hiroshima University)H-Index: 22
Last. Erkan Oterkus (University of Strathclyde)H-Index: 26
view all 4 authors...
In this article, several thermoelastic benchmark cases are studied within the framework of ordinary state based peridynamic theory (OSPD). By using OSPD, the limitations of geometrical discontinuity in fracture analysis can be overcome. Meanwhile, double nodes can also be avoided during crack definition. A domain integral method with thermal effect is applied in calculating the thermal stress intensity factors (TSIFs). Meanwhile, peridynamic differential operators are utilized to rewrite the spa...
#1Ming-Jyun Dai (Hiroshima University)H-Index: 1
#2Satoyuki Tanaka (Hiroshima University)H-Index: 22
Last. Erkan Oterkus (University of Strathclyde)H-Index: 26
view all 4 authors...
Abstract A recently developed peridynamic (PD) shell model based on the ordinary state-based peridynamic theory is adopted to evaluate stress intensity factors (SIFs) under in-plane loading. Strain and stress components are obtained by introducing the peridynamic differential operator. In order to evaluate mixed-mode SIFs, the domain form of the interaction integral is employed. The adaptive dynamic relaxation technique is utilized to obtain steady-state solutions, and the energy method is appli...
3 CitationsSource
#1Septia Hardy Sujiatanti (Hiroshima University)H-Index: 1
#2Satoyuki Tanaka (Hiroshima University)H-Index: 22
Last. Daisuke Yanagihara (Kyushu University)H-Index: 6
view all 5 authors...
Abstract The buckling and collapse behaviors of a cracked thin steel panel subjected to sequential tensile and compressive loading were investigated. An experimental apparatus was developed. Three different test specimens were employed, i.e., an intact panel and cracked panels with different crack lengths. The load-displacement curve and crack opening displacement (COD) were measured during the loading. The experimental results were examined via finite element (FE) computation. The results revea...
#1P.P. Peng (SHU: Shanghai University)H-Index: 1
#2Yumin Cheng (SHU: Shanghai University)H-Index: 38
Abstract In this paper, to improve the computational speed of the reproducing kernel particle method (RKPM), the dimension splitting reproducing kernel particle method (DSRKPM) is presented for solving three-dimensional (3D) transient heat conduction problems. The idea of the proposed method is transforming the 3D problem into a series of two-dimensional (2D) ones by using the finite difference method in the splitting direction. Since the shape function of the RKPM for 2D problem is much simpler...
2 CitationsSource
#1Zheng Liu (Qilu University of Technology)H-Index: 3
#2Gaofeng Wei (Qilu University of Technology)H-Index: 6
Last. Zhiming Wang (Qilu University of Technology)H-Index: 3
view all 3 authors...
Abstract In this paper, the radial basis function (RBF) is introduced into the reproducing kernel particle method (RKPM), and the radial basis reproducing kernel particle method (RRKPM) is proposed for solving geometrically nonlinear problem of functionally graded materials (FGM). Compared with the RKPM, the advantages of the proposed method are that it can eliminate the negative effect of different kernel functions on the computational accuracy, and has higher computational accuracy and stabili...
5 CitationsSource
#1Ming-Jyun Dai (Hiroshima University)H-Index: 1
#2Satoyuki Tanaka (Hiroshima University)H-Index: 22
Last. Tinh Quoc Bui (TITech: Tokyo Institute of Technology)H-Index: 54
view all 5 authors...
Abstract A meshfree approach for analyzing the fracture mechanics parameters in cracked curved shells is presented. The reproducing kernel (RK) meshfree method and mapping technique are employed to approximate cracked curvilinear surfaces and field variables. In order to model the crack segment, the meshfree discretization techniques are used. The stabilized conforming nodal integration (SCNI) and sub-domain stabilized conforming integration (SSCI) techniques are adopted to accurately integrate ...
6 CitationsSource
#1Wentao Ma (Ningxia University)H-Index: 2
#1Wentao Ma (Ningxia University)H-Index: 16
Last. Wenshuai Wang (Ningxia University)H-Index: 1
view all 3 authors...
Abstract A coupled extended meshfree–smoothed meshfree method (CXSMM) is presented for simulating crack growth in 2D solids. In terms of the relationship between the crack surface and the background integration cells, the analysis domain is easy to be divided into discontinuous and continuous regions, in which the meshfree Galerkin method and the smoothed meshfree method are employed to construct the discretized system equations, respectively. To precisely capture the singular stress field near ...
4 CitationsSource
#1Ramy Gadallah (Osaka University)H-Index: 5
#2Seiichiro Tsutsumi (Osaka University)H-Index: 4
Last. Naoki Osawa (Osaka University)H-Index: 8
view all 4 authors...
Abstract An efficient numerical analysis procedure based on FEM was implemented for evaluating the effect of welding residual stress (WRS) on mixed-mode stress intensity factors (MM-SIFs) for a non-planar surface-cracked tubular T-joint with a rounded weld toe. A multi-pass welding simulation was carried out using real welding conditions by means of an effective in-house welding simulation code. A fully compatible 3D finite element mesh generation system was implemented to model non-planar perpe...
6 CitationsSource
#1Kyaw Phyo Myat (Osaka University)H-Index: 1
#2Naoki Osawa (Osaka University)H-Index: 8
Last. Satoyuki Tanaka (Hiroshima University)H-Index: 22
view all 4 authors...
Abstract An effective technique to compute mixed-mode stress intensity factors (MM-SIFs) for inclined surface cracks in semi-infinite bodies under arbitrary stress fields is proposed and validated. The proposed technique in this study is accurate and handiness for calculating MM-SIFs of different crack angles by using only one non-inclined cracked model. The benefit of this technique is achieved by using the influence function method (IFM) to evaluate MM-SIFs. A fully-automated influence coeffic...
1 CitationsSource
#1Pramod Yallappa Kumbhar (Indian Institute of Technology Madras)H-Index: 2
#2Amrita Francis (Indian Institute of Technology Madras)H-Index: 5
Last. Sundararajan Natarajan (Indian Institute of Technology Madras)H-Index: 31
view all 5 authors...
In this paper, we discuss the implementation of a cell-based smoothed finite element method (CSFEM) within the commercial finite element software Abaqus. The salient feature of the CSFEM is that it...
9 CitationsSource