Creep crack-growth: A new path-independent integral (T sub c), and computational studies

The development of valid creep fracture criteria is considered. Two path-independent integral parameters which show some degree of promise are the C* and (Delta T)sub c integrals. The mathematical aspects of these parameters are reviewed by deriving generalized vector forms of the parameters using conservation laws which are valid for arbitrary, three dimensional, cracked bodies with crack surface tractions (or applied displacements), body forces, inertial effects, and large deformations. Two principal conclusions are that (Delta T)sub c has an energy rate interpretation whereas C* does not. The development and application of fracture criteria often involves the solution of boundary/initial value problems associated with deformation and stresses. The finite element method is used for this purpose. An efficient, small displacement, infinitesimal strain, displacement based finite element model is specialized to two dimensional plane stress and plane strain and to power law creep constitutive relations. A mesh shifting/remeshing procedure is used for simulating crack growth. The model is implemented with the quartz-point node technique and also with specially developed, conforming, crack-tip singularity elements which provide for the r to the n-(1+n) power strain singularity associated with the HRR crack-tip field. Comparisons are made with a variety of analytical solutions and alternate numerical solutions for a number of problems

Creep crack-growth: A new path-independent T sub o and computational studies

Two path independent integral parameters which show some degree of promise as fracture criteria are the C* and delta T sub c integrals. The mathematical aspects of these parameters are reviewed. This is accomplished by deriving generalized vector forms of the parameters using conservation laws which are valid for arbitrary, three dimensional, cracked bodies with crack surface tractions (or applied displacements), body forces, inertial effects and large deformations. Two principal conclusions are that delta T sub c is a valid crack tip parameter during nonsteady as well as steady state creep and that delta T sub c has an energy rate interpretation whereas C* does not. An efficient, small displacement, infinitestimal strain, displacement based finite element model is developed for general elastic/plastic material behavior. For the numerical studies, this model is specialized to two dimensional plane stress and plane strain and to power law creep constitutive relations

Stress and Fracture Analyses Under Elastic-plastic and Creep Conditions: Some Basic Developments and Computational Approaches

A new hybrid-stress finite element algorith, suitable for analyses of large quasi-static deformations of inelastic solids, is presented. Principal variables in the formulation are the nominal stress-rate and spin. A such, a consistent reformulation of the constitutive equation is necessary, and is discussed. The finite element equations give rise to an initial value problem. Time integration has been accomplished by Euler and Runge-Kutta schemes and the superior accuracy of the higher order schemes is noted. In the course of integration of stress in time, it has been demonstrated that classical schemes such as Euler's and Runge-Kutta may lead to strong frame-dependence. As a remedy, modified integration schemes are proposed and the potential of the new schemes for suppressing frame dependence of numerically integrated stress is demonstrated. The topic of the development of valid creep fracture criteria is also addressed

Stress-intensity factors for small surface and corner cracks in plates

Three-dimensional finite-element and finite-alternating methods were used to obtain the stress-intensity factors for small surface and corner cracked plates subjected to remote tension and bending loads. The crack-depth-to-crack-length ratios (a/c) ranged from 0.2 to 1 and the crack-depth-to-plate-thickness ratios (a/t) ranged from 0.05 to 0.2. The performance of the finite-element alternating method was studied on these crack configurations. A study of the computational effort involved in the finite-element alternating method showed that several crack configurations could be analyzed with a single rectangular mesh idealization, whereas the conventional finite-element method requires a different mesh for each configuration. The stress-intensity factors obtained with the finite-element-alternating method agreed well (within 5 percent) with those calculated from the finite-element method with singularity elements

Elastic-Plastic Finite Element Analysis of Fatigue Crack Growth in Mode 1 and Mode 2 Conditions

Presented is an alternate cost-efficient and accurate elastic-plastic finite element procedure to analyze fatigue crack closure and its effects under general spectrum loading. Both Modes 1 and 2 type cycling loadings are considered. Also presented are the results of an investigation, using the newly developed procedure, of various factors that cause crack growth acceleration or retardation and delay effects under high-to-low, low-to-high, single overload, and constant amplitude type cyclic loading in a Mode 1 situation. Further, the results of an investigation of a centercracked panel under external pure shear (Mode 2) cyclic loading, of constant amplitude, are reported

Editorial

This article has no abstract

Quasi ALE finite element method for nonlinear water waves

This paper presents a newly developed quasi arbitrary Lagrangian-Eulerian finite element method (QALE-FEM) for simulating water waves based on fully nonlinear potential theory. The main difference of this method from the conventional finite element method developed by one of authors of this paper and others (see, e.g., [11] and [22]) is that the complex mesh is generated only once at the beginning and is moved at all other time steps in order to conform to the motion of the free surface and structures. This feature allows one to use an unstructured mesh with any degree of complexity without the need of regenerating it every time step, which is generally inevitable and very costly. Due to this feature, the QALE-FEM has high potential in enhancing computational efficiency when applied to problems associated with the complex interaction between large steep waves and structures since the use of an unstructured mesh in such a case is likely to be necessary. To achieve overall high efficiency, the numerical techniques involved in the QALE-FEM are developed, including the method to move interior nodes, technique to re-distribute the nodes on the free surface, scheme to calculate velocities and so on. The model is validated by water waves generated by a wavemaker in a tank and the interaction between water waves and periodic bars on the bed of tank. Satisfactory agreement is achieved with analytical solutions, experimental data and numerical results from other methods

A Meshless Local Petrov-Galerkin Method for Solving the Bending Problem of a Thin Plate

Meshless methods have been extensively popularized in literature in recent years, due to their flexibility in solving boundary value problems. The meshless local Petrov-Galerkin(MLPG) method for solving the bending problem of the thin plate is presented and discussed in the present paper. The method uses the moving least-squares approximation to interpolate the solution variables, and employs a local symmetric weak form. The present method is a truly meshless one as it does not need a mesh, either for the purpose of interpolation of the solution or for the integration of the energy. All integrals can be easily evaluated over regularly shaped domains (in general, spheres in three-dimensional problems) and their boundaries. The essential boundary conditions are enforced by the penalty method. Several numerical examples are presented to illustrate the implementation and performance of the present method. The numerical examples presented in the paper show that high accuracy can be achieved for arbitrary nodal distributions for clamped and simply-supported edge conditions. No post processing procedure is required to compute the strain and stress, since the original solution from the present method, using the moving least squares approximation, is of C2 type

Computational Nano-mechanics and Multi-scale Simulation

This article provides a review of the computational nanomechanics, from the ab initio methods to classical molecular dynamics simulations, and multi- temporal and spatial scale simulations. The recent improvements and developments are briefly discussed. Their applications in nanomechanics and nanotubes are also summarized