Abstract:In this paper, we propose an implicit staggered algorithm for crystal plasticity finite element method (CPFEM) which makes use of dynamic relaxation at the constitutive integration level. An uncoupled version of the constitutive system consists of a multi-surface flow law complemented by an evolution law for the hardening variables. Since a saturation law is adopted for hardening, a sequence of nonlinear iteration followed by a linear system is feasible. To tie the constitutive unknowns, the dynamic relaxation method is adopted. A Green-Naghdi plasticity model is adopted based on the Hencky strain calculated using a $[2/2]$ Padé approximation. For the incompressible case, the approximation error is calculated exactly. A enhanced-assumed strain (EAS) element technology is adopted, which was found to be especially suited to localization problems such as the ones resulting from crystal plasticity plane slipping. Analysis of the results shows significant reduction of drift and well defined localization without spurious modes or hourglassing.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to effectively reduce the drift problem caused by the decoupling algorithm and improve the computational efficiency in the crystal plasticity finite element method (CPFEM) analysis of aluminum. Specifically, the author proposes an implicit staggered algorithm, which adopts the dynamic relaxation method at the material constitutive integration level. This method aims to avoid drift through a strongly - coupled solution, thereby providing more accurate and stable results when simulating the anisotropic behavior of metal polycrystals. ### Background and Challenges - **Fully - coupled constitutive system**: The fully - coupled constitutive system of a single crystal involves a large number of constitutive unknowns (such as plastic strain tensors, hardening variables, etc.), which form a nonlinear and often non - smooth system. This makes the computational cost very high, especially in polycrystal simulations. - **Limitations of the decoupling algorithm**: Although the decoupling algorithm can reduce the computational complexity, it is prone to the drift phenomenon, that is, the error of the solution changes with the time step and the order of operations. ### Solutions - **Implicit staggered algorithm**: The author proposes an implicit staggered algorithm, which adopts the dynamic relaxation method at the material constitutive integration level. This method can effectively reduce the drift phenomenon by introducing sub - step lengths and appropriate predictors. - **Strongly - coupled solution**: Through the strongly - coupled solution, the accurate solution of the material constitutive equation is ensured, thereby improving the stability and accuracy of the algorithm. ### Key Technologies - **Dynamic relaxation method**: Within each time step, the hardening variables are updated through multiple iterations and linear combinations to gradually approach the real solution. - **Enhanced assumed strain technology (EAS)**: The EAS technology is adopted, which is particularly suitable for capturing strain localization problems, such as planar slip in crystal plasticity. ### Experimental Verification - **Single - crystal cylindrical tensile test**: The effectiveness of the algorithm is verified through tensile tests of single - crystal cylinders in different directions. The results show that even with a large time step, the dynamic relaxation method can significantly reduce the drift phenomenon. - **Polycrystalline rectangular block compression test**: A compression test is carried out on a polycrystalline rectangular block to show the texture evolution and strain distribution. The results are similar to those reported in the existing literature, but the treatment of edge degrees of freedom is improved. ### Conclusions - **Advantages of the algorithm**: The proposed implicit staggered algorithm performs well in reducing the drift phenomenon and improving the computational efficiency, and can still maintain high precision and stability especially with a large time step. - **Future work**: Further optimize the algorithm and expand its application in more complex materials and structures. Through these technologies and experimental verifications, this paper provides an efficient and stable solution for the crystal plasticity finite element analysis of aluminum.

An implicit staggered algorithm for CPFEM-based analysis of aluminum

A temporal stable smoothed particle finite element method for large deformation problems in geomechanics

Crystal Plasticity Modeling of Texture Evolution and Heterogeneity in Equal Channel Angular Pressing of Aluminum Single Crystal

Open-source implementations and comparison of explicit and implicit crystal-plasticity finite element methods

An Effective Semi-Implicit Integration Scheme for Rate Dependent Crystal Plasticity Using Explicit Finite Element Codes

Investigation on Plastic Flow Behaviors of FCC Polycrystalline Aluminum under Pre-Cyclic Tension-Compression Loading: Experiments and Crystal Plasticity Modeling.

A Robust Integration Algorithm for Implementing Rate Dependent Crystal Plasticity into Explicit Finite Element Method

An Insight into the Deformation and Orientation Development of Severely Plastic Deformed Aluminum

Multi-scale Modelling of Evolving Plastic Anisotropy During Al-alloy Sheet Forming

Crystal Plasticity Investigation of Texture Evolution of Al Single Crystal During Extrusion Through a Bent Channel

A trust region algorithm for finite-strain plasticity with strongly coupled hardening

A Crystalline Plasticity Finite Element Method for Simulation of the Plastic Deformation of Az31 Magnesium Alloys

Explicit incremental-update algorithm for modeling crystal elasto-viscoplastic response in finite element simulation

A Computationally Efficient Multiscale, Multi-Phase Modeling Approach Based on CPFEM to Assess the Effect of Second Phase Particles on Mechanical Properties

An Efficient Parallel-Operational Explicit Algorithm for Taylor-Type Model of Rate Dependent Crystal Plasticity

Key Problems for Rate-Dependent Crystal Plasticity Applied in Finite Element Simulation

Microstructure-based multiscale modeling of large strain plastic deformation by coupling a full-field crystal plasticity-spectral solver with an implicit finite element solver

Simulation of Texture Evolution During Plastic Deformation of FCC, BCC and HCP Structured Crystals with Crystal Plasticity Based Finite Element Method

A multiplicative finite strain crystal plasticity formulation based on additive elastic corrector rates: Theory and numerical implementation

$ρ$-CP: Open Source Dislocation Density Based Crystal Plasticity Framework for Simulating Temperature- and Strain Rate-Dependent Deformation

An asynchronous parallel high-throughput model calibration framework for crystal plasticity finite element constitutive models