Low cycle fatigue life prediction in shot‐peened components of different geometries—part I: residual stress relaxation

In this study, the residual stress relaxation behaviour occurring during low‐cycle fatigue in shot‐peened specimens with either a flat or a notched geometry has been studied. A representative low‐pressure steam turbine material, FV448, was used. The residual stress and strain hardening profiles caused by shot peening were measured experimentally and were then incorporated into a finite element model. By allowing for both effects of shot peening, the residual stress relaxation behaviour was successfully simulated using this model and correlated well with the experimental data. Although more modelling work may be required to simulate the interaction between shot peening effects and external loads in a range of notched geometries, the model predictions are consistent with the specimens tested in the current study. The novelty of this study lies in the development of such a modelling approach which can be used to effectively simulate the complex interaction between shot peening effects and external loads in notched regions. Compared with the un‐notched geometry, the notched geometry was found to be more effective in retaining the improvement in fatigue life resulting from shot peening, by restricting the compressive residual stress relaxation during fatigue loading.     

Probabilistic crack trajectory analysis by a dimension reduction method

This paper proposes a novel analysis method of stochastic crack trajectory based on a dimension reduction approach. The developed method allows efficiently estimating the statistical moments, probability density function and cumulative distribution function of the crack trajectory for cracked elastic structures considering the randomness of the loads, material properties and crack geometries. First, the traditional dimension reduction method is extended to calculate the first four moments of the crack trajectory, in which the responses are eigenvectors rather than scalars. Then the probability density function and cumulative distribution function of the crack trajectory can be obtained using the maximum entropy principle constrained by the calculated moments. Finally, the simulation of the crack propagation paths is realized by using the scaled boundary finite element method. The proposed method is well validated by four numerical examples performed on varied cracked structures. It is demonstrated that this method outperforms the Monte Carlo simulation in terms of computational efficiency, and in the meanwhile, it has an acceptable computational accuracy.     

Velocity‐based formulations for standard and quasi‐incompressible hypoelastic‐plastic solids

We present three velocity‐based updated Lagrangian formulations for standard and quasi‐incompressible hypoelastic‐plastic solids. Three low‐order finite elements are derived and tested for non‐linear solid mechanics problems. The so‐called V‐element is based on a standard velocity approach, while a mixed velocity–pressure formulation is used for the VP and the VPS elements. The two‐field problem is solved via a two‐step Gauss–Seidel partitioned iterative scheme. First, the momentum equations are solved in terms of velocity increments, as for the V‐element. Then, the constitutive relation for the pressure is solved using the updated velocities obtained at the previous step. For the VPS‐element, the formulation is stabilized using the finite calculus method in order to solve problems involving quasi‐incompressible materials. All the solid elements are validated by solving two‐dimensional and three‐dimensional benchmark problems in statics as in dynamics. Copyright © 2016 John Wiley & Sons, Ltd.     

The tetrahedral finite cell method: Higher‐order immersogeometric analysis on adaptive non‐boundary‐fitted meshes

The finite cell method (FCM) is an immersed domain finite element method that combines higher‐order non‐boundary‐fitted meshes, weak enforcement of Dirichlet boundary conditions, and adaptive quadrature based on recursive subdivision. Because of its ability to improve the geometric resolution of intersected elements, it can be characterized as an immersogeometric method. In this paper, we extend the FCM, so far only used with Cartesian hexahedral elements, to higher‐order non‐boundary‐fitted tetrahedral meshes, based on a reformulation of the octree‐based subdivision algorithm for tetrahedral elements. We show that the resulting TetFCM scheme is fully accurate in an immersogeometric sense, that is, the solution fields achieve optimal and exponential rates of convergence for h‐refinement and p‐refinement, if the immersed geometry is resolved with sufficient accuracy. TetFCM can leverage the natural ability of tetrahedral elements for local mesh refinement in three dimensions. Its suitability for problems with sharp gradients and highly localized features is illustrated by the immersogeometric phase‐field fracture analysis of a human femur bone. Copyright © 2016 John Wiley & Sons, Ltd.     

A new multi‐scale dispersive gradient elasticity model with micro‐inertia: Formulation and ‐finite element implementation

Motivated by nano‐scale experimental evidence on the dispersion characteristics of materials with a lattice structure, a new multi‐scale gradient elasticity model is developed. In the framework of gradient elasticity, the simultaneous presence of acceleration and strain gradients has been denoted as dynamic consistency. This model represents an extension of an earlier dynamically consistent model with an additional micro‐inertia contribution to improve the dispersion behaviour. The model can therefore be seen as an enhanced dynamic extension of the Aifantis' 1992 strain‐gradient theory for statics obtained by including two acceleration gradients in addition to the strain gradient. Compared with the previous dynamically consistent model, the additional micro‐inertia term is found to improve the prediction of wave dispersion significantly and, more importantly, requires no extra computational cost. The fourth‐order equations are rewritten in two sets of symmetric second‐order equations so that ‐continuity is sufficient in the finite element implementation. Two sets of unknowns are identified as the microstructural and macrostructural displacements, thus highlighting the multi‐scale nature of the present formulation. The associated energy functionals and variationally consistent boundary conditions are presented, after which the finite element equations are derived. Considerable improvements over previous gradient models are observed as confirmed by two numerical examples. Copyright © 2016 John Wiley & Sons, Ltd.     

A particle‐based moving interface method (PMIM) for modeling the large deformation of boundaries in soft matter systems

The mechanics of the interaction between a fluid and a soft interface undergoing large deformations appear in many places, such as in biological systems or industrial processes. We present an Eulerian approach that describes the mechanics of an interface and its interactions with a surrounding fluid via the so‐called Navier boundary condition. The interface is modeled as a curvilinear surface with arbitrary mechanical properties across which discontinuities in pressure and tangential fluid velocity can be accounted for using a modified version of the extended finite element method. The coupling between the interface and the fluid is enforced through the use of Lagrange multipliers. The tracking and evolution of the interface are then handled in a Lagrangian step with the grid‐based particle method. We show that this method is ideal to describe large membrane deformations and Navier boundary conditions on the interface with velocity/pressure discontinuities. The validity of the model is assessed by evaluating the numerical convergence for a axisymmetrical flow past a spherical capsule with various surface properties. We show the effect of slip length on the shear flow past a two‐dimensional capsule and simulate the compression of an elastic membrane lying on a viscous fluid substrate. Copyright © 2015 John Wiley & Sons, Ltd.     

A two‐scale generalized finite element approach for modeling localized thermoplasticity

Predicting localized, nonlinear, thermoplastic behavior and residual stresses and deformations in structures subjected to intense heating is a prevalent challenge in a range of modern engineering applications. The authors present a generalized finite element method targeted at this class of problems, involving the solution of intrinsically parallelizable local boundary value problems to capture localized, time‐dependent thermo‐elasto‐plastic behavior, which is embedded in the coarse, structural‐scale approximation via enrichment functions. The method accommodates approximation spaces that evolve in between time or load steps while maintaining a fixed global mesh, which avoids the need to map solutions and state variables on changing meshes typical of traditional adaptive approaches. Representative three‐dimensional examples exhibiting localized, transient, nonlinear thermal and thermomechanical effects are presented to demonstrate the advantages of the method with respect to available approaches, especially in terms of its flexibility and potential for realistic future applications in this area. Parallelism of the approach is also discussed. Copyright © 2016 John Wiley & Sons, Ltd.     

Vademecum‐based GFEM (V‐GFEM): optimal enrichment for transient problems

This paper proposes a generalized finite element method based on the use of parametric solutions as enrichment functions. These parametric solutions are precomputed off‐line and stored in memory in the form of a computational vademecum so that they can be used on‐line with negligible cost. This renders a more efficient computational method than traditional finite element methods at performing simulations of processes. One key issue of the proposed method is the efficient computation of the parametric enrichments. These are computed and efficiently stored in memory by employing proper generalized decompositions. Although the presented method can be broadly applied, it is particularly well suited in manufacturing processes involving localized physics that depend on many parameters, such as welding. After introducing the vademecum‐generalized finite element method formulation, we present some numerical examples related to the simulation of thermal models encountered in welding processes. Copyright © 2016 John Wiley & Sons, Ltd.     

An n‐material thresholding method for improving integerness of solutions in topology optimization

It is common in solving topology optimization problems to replace an integer‐valued characteristic function design field with the material volume fraction field, a real‐valued approximation of the design field that permits ‘fictitious’ mixtures of materials during intermediate iterations in the optimization process. This is reasonable so long as one can interpolate properties for such materials and so long as the final design is integer valued. For this purpose, we present a method for smoothly thresholding the volume fractions of an arbitrary number of material phases which specify the design. This method is trivial for two‐material design problems, for example, the canonical topology design problem of specifying the presence or absence of a single material within a domain, but it becomes more complex when three or more materials are used, as often occurs in material design problems. We take advantage of the similarity in properties between the volume fractions and the barycentric coordinates on a simplex to derive a thresholding, method which is applicable to an arbitrary number of materials. As we show in a sensitivity analysis, this method has smooth derivatives, allowing it to be used in gradient‐based optimization algorithms. We present results, which show synergistic effects when used with Solid Isotropic Material with Penalty and Rational Approximation of Material Properties material interpolation functions, popular methods of ensuring integerness of solutions. Copyright © 2016 John Wiley & Sons, Ltd.