E. Ferretti¹, E. Casadio, A. Di Leo¹

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.3, pp. 163-190, 2008, DOI:10.3970/cmes.2008.030.163

Abstract In this study, the Cell Method (CM) is applied in order to investigate the failure mechanisms of masonry walls under shear force. The direction of propagation is computed step-wise by the code, and the domain is updated by means of a propagation technique of intra-element nodal relaxation with re-meshing. The crack extension condition is studied in the Mohr/Coulomb plane, using the criterion of Leon. The main advantage of using the CM for numerical analyses of masonry is that the mortar, the bricks and the interfaces between mortar and bricks can be modeled without any need… More >

S.S. Xu, Y. Dong, Y. Zhang

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.1, pp. 17-26, 2008, DOI:10.3970/cmes.2008.030.017

Abstract A meshless method for arbitrary crack growths is presented. The new method is based on a local partition of unity by introducing additional degrees of freedom that determine the opening of the crack. The crack is modeled with overlapping crack segments located at the nodes. The crack segments are rotated at directional changes of the principal tensile stress such that smearing of the crack is avoided. Such smearing occurs in fixed crack method probably because of inaccurate stress state around the crack tip when the crack propagates. The key feature of our method is that More >

Abdulnaser M. Alshoaibi^1,2, M. S. A. Hadi², A. K. Ariffin²

Structural Durability & Health Monitoring, Vol.3, No.1, pp. 15-28, 2007, DOI:10.3970/sdhm.2007.003.015

Abstract An adaptive finite element method is employed to analyze two-dimensional linear elastic fracture problems. The mesh is generated by the advancing front method and the norm stress error is taken as a posteriori error estimator for the h-type adaptive refinement. The stress intensity factors are estimated by a displacement extrapolation technique. The near crack tip displacements used are obtained from specific nodes of natural six-noded quarter-point elements which are generated around the crack tip defined by the user. The crack growth and its direction are determined by the calculated stress intensity factors as the maximum circumference More >

T. Fujimoto¹ and T. Nishioka¹

CMES-Computer Modeling in Engineering & Sciences, Vol.11, No.2, pp. 91-102, 2006, DOI:10.3970/cmes.2006.011.091

Abstract In the dynamic fracture of metallic material, some cracks propagate with the incidence of plastic deformation, and distinct plastic strain remains near the post-propagation area. In order to elucidate these dynamic nonlinear fracture processes, the moving finite element method is developed for nonlinear crack propagation. The T* integral is used as the parameter to estimate crack tip condition. First, the effect of material viscosity and crack propagation velocity have been discussed based on the numerical results for fracture under pure mode I high speed loading. Under mixed mode loading, numerical simulations for fracture path prediction More >

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

CMC-Computers, Materials & Continua, Vol.3, No.3, pp. 147-154, 2006, DOI:10.3970/cmc.2006.003.147

Abstract In this study, phenomena of multiple branching of dynamically propagating crack are investigated numerically. The complicated paths of cracks propagating in a material are simulated by moving finite element method based on Delaunay automatic triangulation (MFEM BODAT), which was extended for such problems. For evaluation of fracture parameters for propagating and branching cracks switching method of the path independent dynamic J integral was used. Using these techniques the generation phase simulation of multiple dynamic crack branching was performed. Various dynamic fracture parameters, which are almost impossible to obtain by experimental technique alone, were accurately evaluated. More >

G. Dhondt¹

Structural Durability & Health Monitoring, Vol.1, No.1, pp. 21-34, 2005, DOI:10.3970/sdhm.2005.001.021

Abstract An algorithm is described which allows for the automatic calculation of crack propagation due to cyclic loading under mixed-mode conditions. The core of the procedure deals with the insertion of an arbitrarily formed crack into a virgin 20-node brick element mesh, thereby generating new quadratic bricks. One especially difficult aspect is the extension of the triangulation of the crack surface up to the boundary of the crack front elements. In the present article the technique is applied to linear elastic calculations using the stress intensity factor concept and a Paris-type law. However, other crack propagation More >

S.C. Forth¹, A. Staroselsky²

CMES-Computer Modeling in Engineering & Sciences, Vol.9, No.3, pp. 287-298, 2005, DOI:10.3970/cmes.2005.009.287

Abstract A new hybrid surface-integral-finite-element numerical scheme has been developed to model a three-dimensional crack propagating through a thin, multi-layered coating. The finite element method was used to model the physical state of the coating, and the surface integral method was used to model the fatigue crack growth. The two formulations are coupled through the need to satisfy boundary conditions on the crack and external surface. The coupling is sufficiently weak that the surface integral mesh of the crack surface and the finite element mesh of the uncracked volume can be set up independently. Thus, when More >

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

CMC-Computers, Materials & Continua, Vol.1, No.2, pp. 191-204, 2004, DOI:10.3970/cmc.2004.001.191

Abstract Phenomena of dynamic crack branching are investigated numerically from a macroscopic point of view. Repetitive branching phenomena, interaction of cracks after bifurcation and their stability, bifurcation into two and three branches were the objectives of this research. For the analysis of dynamic crack branching, recently we developed moving finite element method based on Delaunay automatic triangulation [Nishioka, Furutuka, Tchouikov and Fujimoto (2002)]. In this study this method was extended to be applicable for complicated crack branching phenomena, such as bifurcation of the propagating crack into more than two branches, multiple crack bifurcation and so on. 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 >

K. Nishimura¹, N. Miyazaki²

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.2, pp. 143-154, 2001, DOI:10.3970/cmes.2001.002.143

Abstract In this paper, we present a classical molecular dynamics algorithm and its implementation on Cray C90 and Fujitsu VPP700. The characters of this algorithm consist in a grid based on the block division of the atomic system and a neighbor list based on the use of a short range potential. The computer program is used for large scale simulations on a Cray C90 and a 32-node VPP700, and measurements of computational performance are reported. Then, we examine the interaction between a crack propagating and a tilt grain boundary under uniaxial tension using this computer program.… More >