Sharif Digital Repository

In this paper, an automated adaptive remeshing procedure is presented for simulation of arbitrary shape crack growth in a 2D finite element mesh. The Zienkiewicz-Zhu error estimator is employed in conjunction with a modified SPR technique based on the recovery of gradients using analytical crack-tip fields in order to obtain more accurate estimation of errors. The optimization of crack-tip singular finite element size is achieved through the adaptive mesh strategy. Finally, several numerical examples are illustrated to demonstrate the effectiveness, robustness and accuracy of computational algorithm in calculation of fracture parameters and prediction of crack path pattern. © 2008 Elsevier... 

In this paper, the three-dimensional automatic adaptive mesh refinement is presented in modeling the crack propagation based on the modified superconvergent patch recovery technique. The technique is developed for the mixed mode fracture analysis of different fracture specimens. The stress intensity factors are calculated at the crack tip region and the crack propagation is determined by applying a proper crack growth criterion. An automatic adaptive mesh refinement is employed on the basis of modified superconvergent patch recovery (MSPR) technique to simulate the crack growth by applying the asymptotic crack tip solution and using the collapsed quarter-point singular tetrahedral elements... 

In this paper, the data transfer operators are developed in 3D large plasticity deformations using superconvergent patch recovery (SPR) method. The history-dependent nature of plasticity problems necessitates the transfer of all relevant variables from the old mesh to new one, which is performed in three main stages. In the first step, the history-dependent internal variables are transferred from the Gauss points of old mesh to nodal points. The variables are then transferred from nodal points of old mesh to nodal points of new mesh. Finally, the values are computed at the Gauss points of new mesh using their values at nodal points. As the solution procedure, in general, cannot be...