1.3.11 Flat rolling: transient and steady-state

Product: ABAQUS/Explicit  

This example illustrates the use of adaptive meshing to simulate a rolling process using both transient and steady-state approaches, as shown in Figure 1.3.11–1. A transient flat rolling simulation is performed using three different methods: a “pure” Lagrangian approach, an adaptive meshing approach using a Lagrangian domain, and a mixed Eulerian-Lagrangian adaptive meshing approach in which material upstream from the roller is drawn from an Eulerian inflow boundary but the downstream end of the blank is handled in a Lagrangian manner. In addition, a steady-state flat rolling simulation is performed using an Eulerian adaptive mesh domain as a control volume and defining inflow and outflow Eulerian boundaries. Solutions using each approach are compared.

Problem description

For each analysis case quarter symmetry is assumed; the model consists of a rigid roller and a deformable blank. The blank is meshed with C3D8R elements. The roller is modeled as an analytical rigid surface using the *SURFACE, TYPE=CYLINDER and *RIGID BODY options. The radius of the cylinder is 175 mm. Symmetry boundary conditions are prescribed on the right (z=0 plane) and bottom (y=0 plane) faces of the blank.

Coulomb friction with a friction coefficient of 0.3 is assumed between the roller and the plate. All degrees of freedom are constrained on the roller except rotation about the z-axis, where a constant angular velocity of 6.28 rad/sec is defined. For each analysis case the blank is given an initial velocity of 0.3 m/s in the x-direction to initiate contact.

The blank is steel and is modeled as a von Mises elastic-plastic material with isotropic hardening. The Young's modulus is 150 GPa, and the initial yield stress is 168.2 MPa. The Poisson's ratio is 0.3; the density is 7800 kg/m3. The *FIXED MASS SCALING option is used to scale the masses of all the blank elements by a factor of 2750 so that the analysis can be performed more economically. This scaling factor represents an approximate upper bound on the mass scaling possible for this problem, above which significant inertial effects would be generated.

The *STEADY STATE DETECTION option is used to define the criteria for stopping the rolling analyses based on the achievement of a steady-state condition. The criteria used require the satisfaction of the steady-state detection norms of equivalent plastic strain, spread, force, and torque within the default tolerances. The exit plane for each norm is defined as the plane passing through the center of the roller with the normal to the plane coincident with the rolling direction. The SAMPLING parameter is set to PLANE BY PLANE for Case 1 through Case 3 for the steady-state detection norms to be evaluated as each plane of elements passes the exit plane. Case 4 requires that the SAMPLING parameter is set to UNIFORM since the initial mesh is roughly stationary due to the initial geometry and the inflow and outflow Eulerian boundaries.

The finite element models used for each analysis case are shown in Figure 1.3.11–2. A description of each model and the adaptive meshing techniques used follows:

Case 1: Transient simulation—pure Lagrangian approach

The blank is initially rectangular and measures 224 × 20 × 50 mm. No adaptive meshing is performed. The analysis is run until steady-state conditions are achieved.

Case 2: Transient simulation—Lagrangian adaptive mesh domain

The finite element model for this case is identical to that used for Case 1, with the exception that a single adaptive mesh domain that incorporates the entire blank is defined to allow continuous adaptive meshing. Symmetry planes are defined as Lagrangian surfaces (the default), and the contact surface on the blank is defined as a sliding surface (the default). The analysis is run until steady-state conditions are achieved.

Case 3: Transient simulation—mixed Eulerian-Lagrangian approach

This analysis is performed on a relatively short initial blank measuring 65 × 20 × 50 mm. Material is continuously drawn by the action of the roller on the blank through an inflow Eulerian boundary defined on the upstream end. The blank is meshed with the same number of elements as in Cases 1 and 2 so that similar aspect ratios are obtained as the blank lengthens and steady-state conditions are achieved.

An adaptive mesh domain is defined that incorporates the entire blank. Because it contains at least one Eulerian surface, this domain is considered Eulerian for the purpose of setting parameter defaults. However, the analysis model has both Lagrangian and Eulerian aspects. The amount of material flow with respect to the mesh will be large at the inflow end and small at the downstream end of the domain. To account for the Lagrangian motion of the downstream end, the MESHING PREDICTOR option on the *ADAPTIVE MESH CONTROLS option is changed from the default of PREVIOUS to CURRENT for this problem. To mesh the inflow end accurately and to perform the analysis economically, the FREQUENCY parameter is set to 5 and the MESH SWEEPS parameter is set to 5.

As in Case 2, symmetry planes are defined as Lagrangian boundary regions (the default), and the contact surface on the blank is defined as a sliding boundary region (the default). In addition, an Eulerian boundary region is defined on the upstream end using the *SURFACE, REGION TYPE=EULERIAN option. Adaptive mesh constraints are defined on the Eulerian surface using the *ADAPTIVE MESH CONSTRAINT option to hold the inflow surface mesh completely fixed while material is allowed to enter the domain normal to the surface. The *EQUATION option is used to ensure that the velocity normal to the inflow boundary is uniform across the surface. The velocity of nodes in the direction tangential to the inflow boundary surface is constrained.

Case 4: Steady-state simulation—Eulerian adaptive mesh domain

This analysis employs a control volume approach in which material is drawn from an inflow Eulerian boundary and is pushed out through an outflow boundary by the action of the roller. The blank geometry for this analysis case is defined such that it approximates the shape corresponding to the steady-state solution: this geometry can be thought of as an “initial guess” to the solution. The blank initially measures 224 mm in length and 50 mm in width and has a variable thickness such that it conforms to the shape of the roller. The surface of the blank transverse to the rolling direction is not adjusted to account for the eventual spreading that will occur in the steady-state solution. Actually, any reasonable initial geometry will reach a steady state, but geometries that are closer to the steady-state geometry often allow a solution to be obtained in a shorter period of time.

As in the previous two cases an adaptive mesh domain is defined on the blank, symmetry planes are defined as Lagrangian surfaces (the default), and the contact surface is defined as a sliding surface (the default). Inflow and outflow Eulerian surfaces are defined on the ends of the blank using the same techniques as in Case 3, except that for the outflow boundary adaptive mesh constraints are applied only normal to the boundary surface and no material constraints are applied tangential to the boundary surface.

To improve the computational efficiency of the analysis, the frequency of adaptive meshing is increased to every fifth increment because the Eulerian domain undergoes very little overall deformation and the material flow speed is much less than the material wave speed. This frequency will cause the mesh at Eulerian boundaries to drift slightly. However, the amount of drift is extremely small and does not accumulate. There is no need to increase the mesh sweeps because this domain is relatively stationary and the default MESHING PREDICTOR setting for Eulerian domains is PREVIOUS. Very little mesh smoothing is required.

Results and discussion

The final deformed configurations of the blank for each of the three transient cases are shown in Figure 1.3.11–3. The transient cases have reached a steady-state solution and have been terminated based on the criteria defined using the *STEADY STATE DETECTION option. Steady-state conditions are determined to have been reached when the reaction forces and moments on the roller have stabilized and the cross-sectional shape and distribution of equivalent plastic strain under the roller become constant over time. When using the *STEADY STATE DETECTION option, these conditions imply that the force, moment, spread, and equivalent plastic strain norms have stabilized such that the changes in the norms over three consecutive sampling intervals have fallen below the user-prescribed tolerances. See Steady-state detection, Section 11.8.1 of the ABAQUS Analysis User's Manual, for a detailed discussion on the definition of the norms. Contours of equivalent plastic strain for each of the three transient cases are in good agreement and are shown in the final configuration of each blank in Figure 1.3.11–4. Figure 1.3.11–5 shows the initial and final mesh configurations at steady state. With the exception of Case 3 all analyses were terminated using the default steady-state norm tolerances. Case 3 required that the force and torque norm tolerances be increased from .005 to .01 due to the force and torque at the roller being rather noisy.

To compare the results from the transient and steady-state approaches, the steady-state detection norms are summarized for each case in Table 1.3.11–1. The table shows a comparison of the values of the steady-state detection norms after the analyses have been terminated. The only significant difference is in the value of the spread norm for Case 4, which is higher than the others. The spread norm is defined as the largest of the second principle moments of inertia of the workpiece's cross-section. Since the spread norm is a cubic function of the lateral deformation of the workpiece, rather small differences in displacements between the test cases can lead to significant differences in the spread norms.

Time history plots of the steady-state detection norms are also shown. Figure 1.3.11–9 and Figure 1.3.11–10 show time history plots of the steady-state force and torque norms, respectively, for all cases. The force and torque norms are essentially running averages of the force and moment on the roller and show good agreement for all four test cases. Figure 1.3.11–7 and Figure 1.3.11–8 show time history plots of the steady-state equivalent plastic strain and spread norms, respectively, for all cases. The equivalent plastic strains norms are in good agreement for all cases.

Input files

Table

Table 1.3.11–1 Comparison of steady-state detection norms.

FormulationSpread normEffective plastic strain normForce norm Torque norm
Case 11.349 E–7.8037–1.43 E63.59 E4
Case 21.369 E–7.8034–1.43 E63.55 E4
Case 31.365 E–7.8018–1.43 E63.61 E4
Case 41.485 E–7.8086–1.40 E63.65 E4


Figures

Figure 1.3.11–1 Diagram illustrating the four analysis approaches used in this problem.

Figure 1.3.11–2 Initial configurations for each case.

Figure 1.3.11–3 Deformed mesh for Cases 1–3.

Figure 1.3.11–4 Contours of equivalent plastic strain for Cases 1–3.

Figure 1.3.11–5 Deformed mesh for Case 4 (shown with initial mesh for comparison).

Figure 1.3.11–6 Contours of equivalent plastic strain for Case 4.

Figure 1.3.11–7 Comparison of equivalent plastic strain norm versus time for all cases.

Figure 1.3.11–8 Comparison of spread norm versus time for all cases.

Figure 1.3.11–9 Comparison of force norm versus time for all cases.

Figure 1.3.11–10 Comparison of torque norm versus time for all cases.