Search results

The purpose of this work is to apply a recently proposed constitutive model for mechanically induced martensitic transformations to the prediction of transformation loci. Additionally, this study aims to elucidate if a stress-assisted criterion can account for transformations in the so-called strain-induced regime.

The model is derived by generalising the stress-based criterion of Patel and Cohen (1953), relying on lattice information obtained using the Phenomenological Theory of Martensite Crystallography. Transformation multipliers (cf. plastic multipliers) are introduced, from which the martensite volume fraction evolution ensues. The associated transformation functions provide a variant selection mechanism. Austenite plasticity follows a classical single crystal formulation, to account for transformations in the strain-induced regime. The resulting model is incorporated into a fully implicit RVE-based computational homogenisation finite element code.

Results show good agreement with experimental data for a meta-stable austenitic stainless steel. In particular, the transformation locus is well reproduced, even in a material with considerable slip plasticity at the martensite onset, corroborating the hypothesis that an energy-based criterion can account for transformations in both stress-assisted and strain-induced regimes.

A recently developed constitutive model for mechanically induced martensitic transformations is further assessed and validated. Its formulation is fundamentally based on a physical metallurgical mechanism and derived in a thermodynamically consistent way, inheriting a consistent mechanical dissipation. This model draws on a reduced number of phenomenological elements and is a step towards the fully predictive modelling of materials that exhibit such phenomena.

An implicit type algorithm for the integration of hypoelastic constitutive equations is proposed for large strain and large rotation conditions. Constitutive relations are derived in a deformation‐neutralized form. This provides the basis for integration in time resulting in an incremental tensor relation. Proposed algorithm can be considered as a generalization of the closest‐point‐projection method in the sense that the projection property applies to a ‘midstep’ rather than the final stress state. Hill's yield criterion under plane stress conditions suitable for metal‐forming applications is used in presented benchmark problems. Numerical results are discussed illustrating the accuracy of the algorithm for different values of the midstep parameter.

Under restriction of an isotropic elastic response of deformed lattice, develops a covariant theory of finite elastoplasticity in principal axes of a pair of deformation tensors. In material description, the tensor pair consists of the plastic deformation tensor and the total deformation Cauchy‐Green tensor. Applies the proposed theory to elastoplastic membrane shells, whose references and current configurations can be arbitrary space‐curved surfaces. Pressure‐insensitive von Mises yield criterion with isotropic hardening and a quadratic form of the strain energy function given in terms of elastic principal stretches are considered as a model problem. Through an explicit enforcement of the plane stress condition we arrive at a reduced two‐dimensional problem representation, which is set in the membrane tangent plane. Numerical implementation of the presented theory relies crucially on the operator split methodology to simplify the state update computation. Presents a set of numerical examples in order to illustrate the performance of the presented methodology and indicate possible applications in the area of sheet metal forming.

The initial stiffness method has been extensively adopted for elasto‐plastic finite element analysis. The main problem associated with the initial stiffness method, however, is its slow convergence, even when it is used in conjunction with acceleration techniques. The Newton‐Raphson method has a rapid convergence rate, but its implementation resorts to non‐symmetric linear solvers, and hence the memory requirement may be high. The purpose of this paper is to develop more advanced solution techniques which may overcome the above problems associated with the initial stiffness method and the Newton‐Raphson method.

In this work, the accelerated symmetric stiffness matrix methods, which cover the accelerated initial stiffness methods as special cases, are proposed for non‐associated plasticity. Within the computational framework for the accelerated symmetric stiffness matrix techniques, some symmetric stiffness matrix candidates are investigated and evaluated.

Numerical results indicate that for the accelerated symmetric stiffness methods, the elasto‐plastic constitutive matrix, which is constructed by mapping the yield surface of the equivalent material to the plastic potential surface, appears to be appealing. Even when combined with the Krylov iterative solver using a loose convergence criterion, they may still provide good nonlinear convergence rates.

Compared to the work by Sloan et al., the novelty of this study is that a symmetric stiffness matrix is proposed to be used in conjunction with acceleration schemes and it is shown to be more appealing; it is assembled from the elasto‐plastic constitutive matrix by mapping the yield surface of the equivalent material to the plastic potential surface. The advantage of combining the proposed accelerated symmetric stiffness techniques with the Krylov subspace iterative methods for large‐scale applications is also emphasized.

Classical continuum models, i.e. continuum models that do not incorporate an internal length scale, suffer from excessive mesh dependence when strain‐softening models are used in numerical analyses and cannot reproduce the size effect commonly observed in quasi‐brittle failure. In this contribution three different approaches will be scrutinized which may be used to remedy these two intimately related deficiencies of the classical theory, namely (i) the addition of higher‐order deformation gradients, (ii) the use of micropolar continuum models, and (iii) the addition of rate dependence. By means of a number of numerical simulations it will be investigated under which conditions these enriched continuum theories permit localization of deformation without losing ellipticity for static problems and hyperbolicity for dynamic problems. For the latter class of problems the crucial role of dispersion in wave propagation in strain‐softening media will also be highlighted.

Contact problems are among the most difficult ones in mechanics. Due to its practical importance, the problem has been receiving extensive research work over the years. The finite element method has been widely used to solve contact problems with various grades of complexity. Great progress has been made on both theoretical studies and engineering applications. This paper reviews some of the main developments in contact theories and finite element solution techniques for static contact problems. Classical and variational formulations of the problem are first given and then finite element solution techniques are reviewed. Available constraint methods, friction laws and contact searching algorithms are also briefly described. At the end of the paper, a bibliography is included, listing about seven hundred papers which are related to static contact problems and have been published in various journals and conference proceedings from 1976.

This paper proposes a unified original general framework, designed to theoretically develop and to extremely easily implement elastoplastic constitutive laws defined in the so called two-invariants space, both in small and finite strain regime.

A general return mapping algorithm is proposed, and particularly a standard procedure is developed to compute the two algorithmic tangent operators, required to solve the Newton–Raphson scheme at the local and global level and thus cast the elastoplastic algorithm within a FEM code.

This work demonstrates that the proposed procedure is fully general and can be applied whatever is the elastic law, the yield surface, the plastic potential function and the hardening law. Several numerical examples are reported, not only to demonstrate the accuracy and robustness of the algorithm, but also explain how to use this general algorithm also in other applications.

The proposed algorithm and its numerical implementation into a FEM code is new and original. The usefulness and the value of the algorithm is twofold: (1) it can be implemented in a small and finite strain simulation FEM code, in order to handle different types of constitutive laws in the same modular way, thus fully leveraging on modern object-oriented coding approach; (2) it can be used as a framework to develop (and then to implement) new constitutive models, since the researcher can simply define the relevant functions (and its main derivatives) and automatically get the numerical algorithm.

Presents a robust and unconditionally stable return‐mapping algorithm based on the discrete counterpart of the principle of maximum plastic dissipation. Develops the explicit expression for the consistent elasto‐plastic tangent modulus. All expressions are derived via tensor formulation showing the advantage over the classical matrix notation. The integration algorithm is implemented in the formulation of the four‐node isoparametric assumed‐strain finite‐rotation shell element employing the Mindlin‐Reissner‐type shell model. By applying the layered model, plastic zones can be displayed through the shell thickness. Material non‐linearity described by the von Mises yield criterion and isotropic hardening is combined with a geometrically non‐linear response assuming finite rotations. Numerical examples illustrate the efficiency of the present formulation in conjunction with the standard Newton iteration approach, in which no line search procedures are required. Demonstrates the excellent performance of the algorithm for large time respective load steps.

The purpose of this paper is to compare the performance of implicit and explicit integration schemes for simulating the metal rolling process using commercial software packages ANSYS™ and LS-DYNA™.

For the industrial application of finite element method, the time discretization is one of the most important factors that determine the stability and efficiency of the analysis. An iterative approach, which is unconditionally stable in linear analyses, is the obvious choice for a quasi-static problem such as metal rolling. However, this approach may be challenging in achieving convergence with non-linear material behavior and complicated contact conditions. Therefore, a non-iterative method is usually adopted, in order to achieve computational accuracy through very small time steps. Models using both methods were constructed and compared for computational efficiency.

The results indicate that the explicit method yields higher levels of efficiency compared to the implicit method as model complexity increases. Furthermore, the implicit method displayed instabilities and numerical difficulties in certain load conditions further disfavoring the solver’s performance.

Comparison of the implicit and explicit procedures for time stepping was applied in 3D finite element analysis of the plate rolling process in order to evaluate and quantify the computational efficiency.

In this paper, an interpolation method called the quasi‐higher order element method (QHOEM) is proposed for solving elastoplastic problems using symmetric Galerkin boundary element method (SGBEM). At the initial stage, it uses higher order elements to interpolate the coordinates and the field variables. Then, for the numerical integration involved, it further uses interpolation to decompose the higher order elements into lower order elements so that the existing analytical integration formulas can be applied. By adopting this procedure, the proposed method yields good adaptability and reduces the computational cost. Its accuracy and efficiency are demonstrated by analyzing three practical examples.