Finite element modeling of crack growth in thin-wall structures by method of combining sub-partition and substructure

Published on May 15, 2018in Engineering Fracture Mechanics3.426
· DOI :10.1016/J.ENGFRACMECH.2018.03.023
Chen Xing1
Estimated H-index: 1
(NUAA: Nanjing University of Aeronautics and Astronautics),
Chuwei Zhou7
Estimated H-index: 7
(NUAA: Nanjing University of Aeronautics and Astronautics)
Abstract A combination of sub-partition and substructure methods in finite element is developed to simulate the crack propagation in plate and shell structures. This method allows modeling arbitrary shape crack independent of element mesh and only the elements overlapping crack need to be operated. An element cut apart by a crack is sub-partitioned into ordinary sub-elements while an element enveloping crack tip is sub-partitioned into several singular sub-elements. The whole sub-elements constitute a substructure and the additional nodal freedom degrees introduced by element sub-partition are condensed to nodes of original mesh. In this way global re-meshing is avoided when the crack grows or new crack nucleates. Good accuracy of the crack tip fields prediction and good adaptability of moving crack simulation by the proposed method are proven through designed examples.
📖 Papers frequently viewed together
5 Citations
5 Citations
#1Yongxiang Wang (Columbia University)H-Index: 13
#2Haim Waisman (Columbia University)H-Index: 26
24 CitationsSource
#1Yongxiang Wang (Columbia University)H-Index: 13
#2Haim Waisman (Columbia University)H-Index: 26
Last. Isaac Harari (TAU: Tel Aviv University)H-Index: 33
view all 3 authors...
17 CitationsSource
#1Silvia Corbani (UFRJ: Federal University of Rio de Janeiro)H-Index: 2
#2J.T.P. Castro (PUC-Rio: Pontifical Catholic University of Rio de Janeiro)H-Index: 5
Last. Anthony R. Ingraffea (Cornell University)H-Index: 61
view all 6 authors...
Abstract In this paper, crack shape evolution under out-of-plane pure bending loads in plates with an initial through-thickness edge crack is investigated experimentally and computationally. Such cracks were numerically simulated by a step-by-step 3D finite element technique as an attempt to reproduce the evolving geometry of the crack front for the present problem, considering contact on the closed part of the crack faces. The purpose of this numerical analysis was to identify the particular is...
6 CitationsSource
#1Yongxiang Wang (Columbia University)H-Index: 13
#2Chiara Cerigato (Columbia University)H-Index: 1
Last. Elena Benvenuti (University of Ferrara)H-Index: 17
view all 4 authors...
Abstract This paper presents a high-order extended finite element method (XFEM) with a novel set of material-dependent enrichment functions for the linear elastic analysis of bimaterial interface cracks. With the proposed material-dependent high-order enrichment functions, highly accurate near-tip displacement and stress fields can be obtained for arbitrary material combinations. The aim of this contribution is the direct evaluation of the complex stress intensity factors (SIFs) of interface cra...
21 CitationsSource
#1Qinglei Zeng (THU: Tsinghua University)H-Index: 8
#2Zhanli Liu (THU: Tsinghua University)H-Index: 20
Last. Zhuo Zhuang (THU: Tsinghua University)H-Index: 22
view all 5 authors...
Summary In this paper, the extended finite element method (X-FEM) formulation for the modeling of arbitrary crack propagation in coupled shell/solid structures is developed based on the large deformation continuum-based (CB) shell theory. The main features of the new method are as follows: (1) different kinematic equations are derived for different fibers in CB shell elements, including the fibers enriched by shifted jump function or crack tip functions and the fibers cut into two segments by th...
25 CitationsSource
#1Yongxiang Wang (Columbia University)H-Index: 13
#2Haim Waisman (Columbia University)H-Index: 26
Abstract Failure of quasi-brittle materials is governed by crack formation and propagation which can be characterized by two phases: (i) diffuse material degradation process with initial crack formation and (ii) severe localization of damage leading to the propagation of large cracks and fracture. While continuum damage mechanics provides an excellent framework to describe the first failure phase, it is unable to represent discontinuous displacement fields. In sharp contrast, cohesive zone model...
73 CitationsSource
#1D.H. Li (THU: Tsinghua University)H-Index: 13
#2Y. Liu (THU: Tsinghua University)H-Index: 4
Last. Xiong Zhang (THU: Tsinghua University)H-Index: 26
view all 3 authors...
Summary For the delamination and matrix crack prediction of composite laminated structures, the methods based on the damage mechanics and fracture mechanics are most commonly used. However, there are very few methods that can accurately simulate the delaminations together with matrix cracks, although the in-plane matrix cracks always exist alongside the delaminations under impact loading. In this work, an extended layerwise method is developed to model the composite laminated beam with multiple ...
37 CitationsSource
#1Yang Jiao (Columbia University)H-Index: 5
#2Jacob Fish (Columbia University)H-Index: 52
Two hierarchical approaches, the \(s\)-method and the extended finite element method (XFEM), are compared to the classical ply-by-ply discretization approach in terms of their effectiveness in modeling delamination in laminated composites. In the two hierarchical approaches, a smooth approximation field based on the mesh made of through-the-thickness solid laminated elements is first introduced to resolve a delamination-free response of the composite structure. The initiation and propagation of ...
11 CitationsSource
#1Yongxiang Wang (Columbia University)H-Index: 13
#2Haim Waisman (Columbia University)H-Index: 26
The modeling of progressive delamination by means of a discrete damage zone model within the extended finite element method is investigated. This framework allows for both bulk and interface damages to be conveniently traced, regardless of the underlying mesh alignment. For discrete interfaces, a new mixed-mode force---separation relation, which accounts for the coupled interaction between opening and sliding modes, is proposed. The model is based on the concept of Continuum Damage Mechanics and...
72 CitationsSource
24 CitationsSource
Cited By7
Abstract In this study, a systematic three-dimensional finite element analyses (3D FEA) are performed to evaluate the Mode II stress intensity factor (SIF) for circumferential through-wall cracked pipe under torsion moment. A wide range of crack lengths and pipe radius over thickness ratios are considered. Empirical solution of the mid-thickness Mode II SIF is developed as functions of crack lengths and pipe sizes. The deviating angle and the effective SIF regarding a through-wall crack (TWC) un...
#1Hongjun Yu (HIT: Harbin Institute of Technology)H-Index: 15
#2Meinhard Kuna (Freiberg University of Mining and Technology)H-Index: 31
Abstract A critical overview of past research on the interaction integral (I-integral) method is presented. The I-integral is the two-state mutual energy release rate decided by a designable auxiliary field and the actual field. Due to the designability of the auxiliary field, the I-integral has been used extensively in extracting the individual stress intensity factors (SIFs) and the T-stress of a crack in single materials or at a bi-material interface. The practical implementation of the I-int...
#1M.T. Mohammadi Anaei (Shiraz University)
#2A. Khosravifard (Shiraz University)H-Index: 12
Last. Tinh Quoc Bui (TITech: Tokyo Institute of Technology)H-Index: 54
view all 3 authors...
Abstract This paper presents a robust meshfree technique for computation of moment and stress intensity factors, and analysis of fatigue crack growth in cracked plates under the application of in- and out of-plane loads; a methodology for efficient and accurate computation of parameters of fracture mechanics problems is thus developed. Although having small computational labor and few complications of formulation, the proposed method has a high accuracy and convergence rate. By adding a few node...
#1Thuan Ho-Nguyen-Tan (Seoul National University of Science and Technology)H-Index: 3
#2Hyun Gyu Kim (Seoul National University of Science and Technology)H-Index: 28
In this study, we present a novel finite element approach to simulate crack propagation in shell structures. A local spider-web mesh is placed at the tip of a crack propagating through a global background mesh. Interface shell elements with assumed natural strains are used to connect a non-matching interface between the background mesh and the local spider-web mesh. Interface shell elements are also employed for trimmed shell elements created by cutting shell elements with the crack line. Numeri...
An increasing number of the offshore structures in the North Sea are being operated in a life extension phase, which means that their original design life has been exceeded. As the structures age, they are deteriorating, mainly due to fatigue and corrosion. In order to ensure that the structural integrity is maintained during the life extension phase, the remaining fatigue life of the structures needs to be assessed. While simple and conservative approximations are used in fatigue design, it wil...
#1Qing Shuai Cao (ZUCC: Zhejiang University City College)H-Index: 1
#2Yang Zhao (ZJU: Zhejiang University)H-Index: 10
Last. Ben Yue Li (ZJU: Zhejiang University)H-Index: 1
view all 5 authors...
Abstract Large steel silos are typical kinds of thin-walled structures that are usually constructed with a vertical, circular form and are used as vessels for storing particulate granular solids in the industrial and agricultural sectors. In the present analyses, the buckling behaviours of circular steel silos with circumferential and meridional fabrication cracks are evaluated in accordance with Eurocodes EN1990, 1991, and 1993 [1–6]. A steel silo with a diameter of 10 m and a height of 24 m is...
1 CitationsSource
4 CitationsSource
#1Yongxiang Wang (WHU: Wuhan University)H-Index: 13
#2Haim Waisman (Columbia University)H-Index: 26
Abstract Numerical modeling of cohesive crack growth in quasi-brittle materials is challenging, primarily due to the combination of (i) nonlinearity associated with the fracture process zone (FPZ), (ii) arbitrary directions to which a crack may propagate, and (iii) snap-back or snap-through instabilities encountered in the response of the structure. To address these challenges, we propose a novel arc-length method that can follow the equilibrium path of cohesive crack propagation. The proposed a...
10 CitationsSource