Simulation of springback: Through-thickness integration

• Published in: International journal of plasticity Vol. 23; no. 3; pp. 345 - 360
• Main Authors: Wagoner, R.H., Li, M.
• Format: Journal Article Conference Proceeding
• Language: English
• Published: Oxford Elsevier Ltd 01.03.2007; Elsevier Science
• Subjects: Accuracy; Applied sciences; Errors; Exact sciences and technology; Forming; Fundamental areas of phenomenology (including applications); Inelasticity (thermoplasticity, viscoplasticity...); Integration; Integration points; IP (Internet Protocol); Mathematical analysis; Mathematical models; Metals. Metallurgy; Numerical integration; Other forming methods; Permissible error; Physics; Production techniques; Shell elements; Solid mechanics; Spring-back; Springback; Structural and continuum mechanics; Shell elements; Integration points; Spring-back; Integration; Springback; Numerical integration; Forming; Modeling; Exact solution; Finite element method; Inelasticity; Elastoplasticity; Bending moment; Mechanical tension
• ISSN: 0749-6419; 1879-2154
• DOI: 10.1016/j.ijplas.2006.04.005

The number of through-thickness integration points ( N IP) required for accurate springback analysis following sheet forming simulation using shell elements is a subject of confusion and controversy. Li and Wagoner recommended, in 1999, based on a finite element analysis (FEA) of draw-bending springback, the use of 25 integration points (IP), with up to 51 IP required to ensure accuracies of 1%. Several researchers have since reported that N IP between 5 and 11 are adequate, or even that 7 or 9 IP are optimal, with reduced accuracy for more IP. These apparent contradictions are addressed with an analytical model of elasto-plastic bending under tension, followed by elastic springback. The fractional error in the evaluated bending moment, which is equal to the fractional error in springback, was determined by comparing three numerical integration schemes, with various N IP, to the closed-form result. The results illustrate the oscillatory nature of numerical integration error with small parametric changes, such that fortuitous agreement can be obtained in isolated simulations where the number of integration points is inadequate. The concept of an assured error limit is introduced as well as a maximum error limit (for a range of generally unknown sheet tensions). The assured error limit varies with the integration scheme, N IP, bending ratio ( R/ t), and sheet tension. Guidelines for the number of integration points required for given error tolerances are reported to allow practitioners to choose numerical parameters appropriately.