Daiki Shiozawa¹, Shiro Kubo², Takahide Sakagami², Masaaki Takagi²

CMES-Computer Modeling in Engineering & Sciences, Vol.11, No.1, pp. 27-36, 2006, DOI:10.3970/cmes.2006.011.027

Abstract The passive electric potential CT (computed tomography) method using piezoelectric film was applied to the identification of plural through cracks. The use of piezoelectric material made it possible to obtain electric potential field without applying electric current. For identification of cracks an inverse analysis scheme based on the least residual method was applied, in which square sum of residuals is evaluated between the measured electric potential distributions and those computed by using the finite element method. Akaike information criterion (AIC) was used to estimate the number of cracks. Numerical simulations were carried out on the identification of plural cracks and… More >

Jae Bong Lee¹, Jai Hak Park¹, Sung Ho Lee², Hong-Deok Kim², Han-Sub Chung²

CMES-Computer Modeling in Engineering & Sciences, Vol.11, No.1, pp. 9-16, 2006, DOI:10.3970/cmes.2006.011.009

Abstract The growth of stress corrosion cracks in steam generator tubes is predicted using the Monte Carlo simulation and statistical approaches. The statistical parameters that represent the characteristics of crack growth and crack initiation are derived from in-service inspection (ISI) non-destructive evaluation (NDE) data. Based on the statistical approaches, crack growth models are proposed and applied to predict crack distribution at the end of cycle (EOC). Because in-service inspection (ISI) crack data is different from physical crack data, a simple method for predicting the physical number of cracks from periodic in-service inspection data is proposed in this study. Actual number of… More >

U. Andreaus^1,3, R.C. Batra², M. Porfiri^{2, 3}

CMES-Computer Modeling in Engineering & Sciences, Vol.9, No.2, pp. 111-132, 2005, DOI:10.3970/cmes.2005.009.111

Abstract Structural health monitoring techniques based on vibration data have received increasing attention in recent years. Since the measured modal characteristics and the transient motion of a beam exhibit low sensitivity to damage, numerical techniques for accurately computing vibration characteristics are needed. Here we use a Meshless Local Petrov-Galerkin (MLPG) method to analyze vibrations of a beam with multiple cracks. The trial and the test functions are constructed using the Generalized Moving Least Squares (GMLS) approximation. The smoothness of the GMLS basis functions requires the use of special techniques to account for the slope discontinuities at the crack locations. Therefore, a… More >

H. Okada ¹, T. Kamibeppu ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.10, No.3, pp. 229-238, 2005, DOI:10.3970/cmes.2005.010.229

Abstract In this paper, a three-dimensional VCCM (Virtual Crack Closure-Integral Method) for evaluating the energy release rate and the stress intensity factor is presented. Many engineers and researchers believe that hexahedral finite elements should be used to perform three-dimensional fracture analyses. Previous VCCM formulations assume the use of hexahedral finite elements. In present study, the authors have been developing a VCCM that works with tetrahedral finite elements. In the field of large-scale computation, the use of tetrahedral finite elements has becoming very popular as high performance mesh generation programs became available. Therefore, building a large and complex analysis model with hexahedral… More >

Toshihisa Nishioka ¹,

CMES-Computer Modeling in Engineering & Sciences, Vol.10, No.3, pp. 209-216, 2005, DOI:10.3970/cmes.2005.010.209

Abstract Recent Advances in Numerical Simulation Technologies for Various Dynamic Fracture Phenomena are summarized. First, the basic concepts of fracture simulations are explained together with pertinent simulation results. Next, Examples of dynamic fracture simulations are presented. More >

V.V. Mykhas’kiv¹, O.I. Stepanyuk²

CMES-Computer Modeling in Engineering & Sciences, Vol.10, No.1, pp. 45-64, 2005, DOI:10.3970/cmes.2005.010.045

Abstract The 3D elastostatic problem for an infinite remotely loaded matrix containing a finite number of arbitrarily located rigid disk-inclusions and plane cracks is solved by the boundary integral equation (BIE) method. Its boundary integral formulation is achieved by the superposition principle with the subsequent integral representations of superposition terms through surface integrals, which should satisfy the displacement linearity conditions in the inclusion domains and load-free conditions in the crack domains. The subtraction technique in the conjunction with mapping technique under taking into account the structure of the solution at the edges of inhomogeneities is applied for the regularization of BIE… More >

E. Ferretti¹

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.3, pp. 227-244, 2004, DOI:10.3970/cmes.2004.006.227

Abstract Defining the crack path in brittle and non-brittle crack is not easy, due to several unknowns. If the direction of crack propagation can be computed by means of one of the existing criteria, it is not known whether this direction will remain constant during crack propagation. A crack initiation leads to an enhanced stress field at crack tip. During propagation, the enhanced tip stress field propagates into the solid, locally interacting with the pre-existing stress field. This interaction can lead to modifications of the propagation direction, depending on the domain and crack geometry. Moreover, trajectory deviation affects the length of… More >

Ch. Zhang², A. Savaidis³

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.5, pp. 603-618, 2003, DOI:10.3970/cmes.2003.004.603

Abstract A novel non-hypersingular time-domain traction BEM is presented for three-dimensional (3-D) transient elastodynamic crack analysis. The initial-boundary value problem is formulated as a set of non-hypersingular time-domain traction boundary integral equations (BIEs). To solve the time-domain traction BIEs, a time-stepping scheme based on the convolution quadrature formula of Lubich (1988a,b; 1994) for temporal discretization and a collocation method for spatial discretization is adopted. Numerical examples are given for an unbounded solid with a penny-shaped crack under a tensile and shear impact loading. A comparison of the present time-domain BEM with the conventional one shows that the novel time-domain method is… More >

E. Ferretti¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.1, pp. 51-72, 2003, DOI:10.3970/cmes.2003.004.051

Abstract A numerical code for modeling crack propagation using the cell method is proposed. The Mohr-Coulomb criterion is used to compute the direction of crack propagation, and the new crack geometry is realized by an intra-element propagation technique. Automatic remeshing is then activated. Applications in Mode I and Mixed Mode are presented to illustrate the robustness of the implementation. More >

T. Nishioka¹, J. Furutuka¹, S. Tchouikov¹, T. Fujimoto¹

CMES-Computer Modeling in Engineering & Sciences, Vol.3, No.1, pp. 129-146, 2002, DOI:10.3970/cmes.2002.003.129

Abstract The governing condition of dynamic crack bifurcation phenomena had not been fully elucidated until our recent experimental studies. We found from the experimental results that the energy flux per unit time into a propagating crack tip or into a fracture process zone governs the crack bifurcation. Regarding the numerical simulation of dynamic crack bifurcation, to the authors' knowledge, no accurate simulations have been carried out, due to several unresolved difficulties. In order to overcome the difficulties, for the analysis of dynamic crack bifurcation, we developed a moving finite element method based on Delaunay automatic triangulation. Using the moving finite element… More >