Search results

The purpose of this study is to simultaneously determine the constitutive parameters and boundary conditions by solving inverse mechanical problems of power hardening elastoplastic materials in three-dimensional geometries.

The power hardening elastoplastic problem is solved by the complex variable finite element method in software ABAQUS, based on a three-dimensional complex stress element using user-defined element subroutine. The complex-variable-differentiation method is introduced and used to accurately calculate the sensitivity coefficients in the multiple parameters identification method, and the Levenberg–Marquardt algorithm is applied to carry out the inversion.

Numerical results indicate that the complex variable finite element method has good performance for solving elastoplastic problems of three-dimensional geometries. The inversion method is effective and accurate for simultaneously identifying multi-parameters of power hardening elastoplastic problems in three-dimensional geometries, which could be employed for solving inverse elastoplastic problems in engineering applications.

The constitutive parameters and boundary conditions are simultaneously identified for power hardening elastoplastic problems in three-dimensional geometries, which is much challenging in practical applications. The numerical results show that the inversion method has high accuracy, good stability, and fast convergence speed.

The constitutive equations for the deformation of elastoplastic, viscoplastic or compressible materials are presented for the small strain approximation and for the large strain theory of Hill. A velocity approach is proposed for time discretization, which leads to a second order approximation for small strain, and an incrementally objective second order approximation for large deformation processes. Two other quasi second order formulations are discussed. The finite element space discretization is outlined and the solution procedure is described.

The finite element quasi‐static analysis of elastoplastic systems is studied by making use of a generalized variable approach for the spatial discretization and a generalized mid‐point rule for the time integration. Both the classical form of the constitutive law and the convex analysis formulation are presented. The relation between the mid‐point time integration and the extremal path theory is discussed. Extremal properties for the finite‐step solution are formulated.

The present contribution deals with the sensitivity analysis and optimization of structures for path‐dependent structural response. Geometrically as well as materially non‐linear behavior with hardening and softening is taken into account. Prandtl‐Reuss‐plasticity is adopted so that not only the state variables but also their sensitivities are path‐dependent. Because of this the variational direct approach is preferred for the sensitivity analysis. For accuracy reasons the sensitivity analysis has to be consistent with the analysis method evaluating the structural response. The proposed sensitivity analysis as well as its application in structural optimization is demonstrated by several examples.

The purpose of this paper is to present an efficient return mapping algorithm for elastoplastic constitutive problems of ductile metals with an exact closed form solution of the local constitutive problem in the small strain regime. A Newton Raphson iterative method is adopted for the solution of the boundary value problem.

An efficient return mapping algorithm is illustrated which is based on an elastic predictor and a plastic corrector scheme resulting in an implicit and accurate numerical integration method. Nonlinear kinematic hardening rules and linear isotropic hardening rules are used to describe the components of the hardening variables. In the adopted algorithmic approach the solution of the local constitutive equations reduces to only one straightforward nonlinear scalar equation.

The presented algorithmic scheme naturally leads to a particularly simple form of the nonlinear scalar equation which ultimately scales down to an algebraic (polynomial) equation with a single variable. The straightforwardness of the present approach allows to find the analytical solution of the algebraic equation in a closed form. Further, the consistent tangent operator is derived as associated with the proposed algorithmic scheme and it is shown that the proposed computational procedure ensures a quadratic rate of asymptotic convergence when used with a Newton Raphson iterative method for the global solution procedure.

In the present approach the solution of the algebraic nonlinear equation is found in a closed form and accordingly no iterative method is required to solve the problem of the local constitutive equations. The computational procedure ensures a quadratic rate of asymptotic convergence for the global solution procedure typical of computationally efficient solution schemes. In the paper it is shown that the proposed algorithmic scheme provides an efficient and robust computational solution procedure for elastoplasticity boundary value problems. Numerical examples and computational results are reported which illustrate the effectiveness and robustness of the adopted integration algorithm for the finite element analysis of elastoplastic structures also under elaborate loading conditions.

The importance of contact and friction problems in different application areas is discussed. Methods and algorithms for numerical solutions using the finite element method are presented. Both elastic and elastic plastic materials are included as well as combination of contact and crack problems. The methods are applied to practical applications such as bolted joints, lugs and roller bearings.

The boundary element method (BEM) and the finite element method (FEM) may be computationally expensive if complex problems are to be solved; thus there is the need of implementing them on fast computer architectures, especially parallel computers. Because these methods are complementary to each other, the coupling of FEM and BEM is widely used. In this paper, the coupling of displacement‐based FEM and collocation BEM and its implementation on a distributed memory system (Parsytec MultiCluster2) is described. The parallelization is performed by data partitioning which leads to a very high efficiency. As model problems, we assume linear elasticity for the boundary element method and elastoplasticity for the finite element method. The efficiency of our implementation is shown by various test examples. By numerical examples we show that a multiplicative Schwarz method for coupling BEM with FEM is very well suited for parallel implementation.

In the framework of the finite element method, the problem of elasto‐plastic consolidation gives rise to a system of non‐linear, coupled residual equations which satisfy the conditions of balance of momentum and balance of mass. In determining the roots of these equations it is necessary that the coupled equations be linearized. To this end, the concept of ‘consistent linearization’ proposed by Simo and Taylor for a single‐phase system is applied to the two‐phase soil‐water system. The roots of the coupled residual equations are solved iteratively by employing Newton's method. It is shown that in non‐linear consolidation analyses, the use of a tangent coefficient matrix derived consistently from the integrated constitutive equation defining the characteristics of the solid skeletal phase results in an iterative solution scheme which preserves the asymptotic rate of quadratic convergence of Newton's method. Numerical examples involving combined radial and vertical flows through an elasto‐plastic soil medium are presented to demonstrate the computational superiority of the above technique over the method based on standard ‘elasto‐plastic continuum formulations’ adopted in most finite element codes.

Large strain model can be formulated in terms of the Lagrangian or the Eulerian frame. In this paper, the Eulerian type large strain models are studied. Numerical examples on the Lagrangian and Eulerian types large strain models are investigated and compared. It is found that the differences in the choice of large strain model under large strain and rotation problems are noticeable but not significant if small load step is used for analysis. Furthermore, we have also found that unsymmetrical formulation instead of symmetrical formulation should be adopted for Eulerian type large strain models.

This work addresses the computational aspects of a model for elastoplastic damage at finite strains. The model is a modification of a previously established model for large strain elastoplasticity described by Perić et al. which is here extended to include isotropic damage and kinematic hardening. Within the computational scheme, the constitutive equations are numerically integrated by an algorithm based on operator split methodology (elastic predictor—plastic corrector). The Newton—Raphson method is used to solve the discretized evolution equations in the plastic corrector stage. A numerical assessment of accuracy and stability of the integration algorithm is carried out based on iso‐error maps. To improve the stability of the local N—R scheme, the standard elastic predictor is replaced by improvedinitial estimates ensuring convergence for large increments. Several possibilities are explored and their effect on the stability of the N—R scheme is investigated. The finite element method is used in the approximation of the incremental equilibrium problem and the resulting equations are solved by the standard Newton—Raphson procedure. Two numerical examples are presented. The results are compared with those obtained by the original elastoplastic model.