13.2 Loading rates

The actual time taken for a physical process is called its natural time. Generally, it is safe to assume that performing an analysis in the natural time for a quasi-static process will produce accurate static results. After all, if the real-life event actually occurs in a natural time scale in which velocities are zero at the conclusion, a dynamic analysis should be able to capture the fact that the analysis has, in fact, achieved a steady state. You can increase the loading rate so that the same physical event occurs in less time as long as the solution remains nearly the same as the true static solution and dynamic effects remain insignificant.


13.2.1 Smooth amplitude curves

For accuracy and efficiency quasi-static analyses require the application of loading that is as smooth as possible. Sudden, jerky movements cause stress waves, which can induce noisy or inaccurate solutions. Applying the load in the smoothest possible manner requires that the acceleration changes only a small amount from one increment to the next. If the acceleration is smooth, it follows that the changes in velocity and displacement are also smooth.

ABAQUS has a simple, built-in smooth step amplitude curve that automatically creates a smooth loading amplitude. When you define a smooth step amplitude curve, ABAQUS automatically connects each of your data pairs with curves whose first and second derivatives are smooth and whose values are zero at each of your data points. Since both of these derivatives are smooth, you can apply a displacement loading with a smooth step amplitude curve using only the initial and final data points, and the intervening motion will be smooth. Using this type of loading amplitude allows you to perform a quasi-static analysis without generating waves due to discontinuity in the rate of applied loading. An example of a smooth step amplitude curve is shown in Figure 13–2.

Figure 13–2 Amplitude definition using a smooth step amplitude curve.


13.2.2 Structural problems

In a static analysis the lowest mode of the structure usually dominates the response. Knowing the frequency and, correspondingly, the period of the lowest mode, you can estimate the time required to obtain the proper static response. To illustrate the problem of determining the proper loading rate, consider the deformation of a side intrusion beam in a car door by a rigid cylinder, as shown in Figure 13–3. The actual test is quasi-static.

Figure 13–3 Rigid cylinder impacting beam.

The response of the beam varies greatly with the loading rate. At an extremely high impact velocity of 400 m/s, the deformation in the beam is highly localized, as shown in Figure 13–4. To obtain a better quasi-static solution, consider the lowest mode.

Figure 13–4 Impact velocity of 400 m/s.

The frequency of the lowest mode is approximately 250 Hz, which corresponds to a period of 4 milliseconds. The natural frequencies can easily be calculated using the eigenfrequency extraction procedure in ABAQUS/Standard. To deform the beam by the desired 0.2 m in 4 milliseconds, the velocity of the cylinder is 50 m/s. While 50 m/s still seems like a high impact velocity, the inertial forces become secondary to the overall stiffness of the structure, and the deformed shape—shown in Figure 13–5—indicates a much better quasi-static response. While the overall structural response appears to be what we expect as a quasi-static solution, it is usually desirable to increase the loading time to 10 times the period of the lowest mode to be certain that the solution is truly quasi-static. To improve the results even further, the velocity of the rigid cylinder could be ramped up gradually—for example, using a smooth step amplitude curve—thereby easing the initial impact.

Figure 13–5 Impact velocity of 50 m/s.


13.2.3 Metal forming problems

Artificially increasing the speed of forming events is necessary to obtain an economical solution, but how large a speedup can we impose and still obtain an acceptable static solution? If the deformation of the sheet metal blank corresponds to the deformed shape of the lowest mode, the time period of the lowest structural mode can be used as a guideline for forming speed. However, in forming processes the rigid dies and punches can constrain the blank in such a way that its deformation may not relate to the structural modes. In such cases a general recommendation is to limit punch speeds to less than 1% of the sheet metal wave speed. For typical processes the punch speed is on the order of 1 m/s, while the wave speed of steel is approximately 5000 m/s. This recommendation, therefore, suggests a factor of 50 as an upper bound on the speedup of the punch velocity.

The suggested approach to determining an acceptable punch velocity involves running a series of analyses at various punch speeds in the range of 3 to 50 m/s. Perform the analyses in the order of fastest to slowest since the solution time is inversely proportional to the punch velocity. Examine the results of the analyses, and get a feel for how the deformed shapes, stresses, and strains vary with punch speed. Some indications of excessive punch speeds are unrealistic, localized stretching and thinning as well as the suppression of wrinkling. If you begin with a punch speed of, for example, 50 m/s, and decrease it from there, at some point the solutions will become similar from one punch speed to the next—an indication that the solutions are converging on a quasi-static solution. As inertial effects become less significant, differences in simulation results also become less significant.

As the loading rate is increased artificially, it becomes more and more important to apply the loads in a gradual and smooth manner. For example, the simplest way to load the punch is to impose a constant velocity throughout the forming step. Such a loading causes a sudden impact load onto the sheet metal blank at the start of the analysis, which propagates stress waves through the blank and may produce undesired results. The effect of any impact load on the results becomes more pronounced as the loading rate is increased. Ramping up the punch velocity from zero using a smooth step amplitude curve minimizes these adverse effects.

Springback

Springback is often an important part of a forming analysis because the springback analysis determines the shape of the final, unloaded part. While ABAQUS/Explicit is well-suited for forming simulations, springback poses some special difficulties. The main problem with performing springback simulations within ABAQUS/Explicit is the amount of time required to obtain a steady-state solution. Typically, the loads must be removed very carefully, and damping must be introduced to make the solution time reasonable. Fortunately, the close relationship between ABAQUS/Explicit and ABAQUS/Standard allows a much more efficient approach.

Since springback involves no contact and usually includes only mild nonlinearities, ABAQUS/Standard can solve springback problems much faster than ABAQUS/Explicit can. Therefore, the preferred approach to springback analyses is to import the completed forming model from ABAQUS/Explicit into ABAQUS/Standard.