Load Stepping

For cases with nonlinear geometry (large deformations) and nonlinear materials (plastic or hyperelastic solids), the governing equations can be highly nonlinear, causing the solution to diverge. To decrease the nonlinearity of the system and improve convergence, it is convenient to apply the external loads progressively over several steps.

Considering the static nonlinear problem:

Figure 1. EQUATION_DISPLAY
r ( u ) : = f int f e x t = 0
(4627 4628 4631)

With load stepping, Eqn. (4627) can be written as:

Figure 2. EQUATION_DISPLAY
r ( u s ) : = f int λ s f e x t = 0
(4627 4628 4631)

where the stepping factor λ s accounts for the ramping of the load over the steps s . Currently, Simcenter STAR-CCM+ uses the load loop approach. With this approach, the external load is split into parts which are incrementally applied in separate load steps. In each load step, the nonlinear sytem Eqn. (4627) is solved with the Newton-Raphson method:

Figure 3. EQUATION_DISPLAY
u k + 1 s = u s 1 + Δ u k + 1 s
(4629 4630)
Figure 4. EQUATION_DISPLAY
Δ u k + 1 s = Δ u k s + d u
(4629 4630)
Figure 5. EQUATION_DISPLAY
r ( u k + 1 s ) = f int ( u k + 1 s ) λ s f e x t = 0
(4627 4628 4631)

Linear Load Stepping

The Linear Load Step Control applies the load increment at each step based on the convergence of the previous load steps. If the current load step is:
  • The first load step of the simulation, then the user-defined Initial load increment is used.
  • The first load step of the time-step, then all load steps of the first load loop of the previous time-step are considered.
  • The first load step of a load loop, then all load steps of the previous load loop are considered.
  • A subsquent step of a load loop, then all load steps of the current load loop are considered.
The following diagram illustrates the above conditions:

At each load step the solver increments the load linearly using a least squares method:
a = i x i y i i x i x i
(4632)
x = y a
(4633)
Where y is the optimal number of newton steps, x is the load step size of the next load step, y i is the number of iterations of a previous load step i , and, x i is the load step size of a previous load step i .