Geometrically nonlinear static problems sometimes involve buckling or collapse behavior, where the load-displacement response shows a negative stiffness, and the structure must release strain energy to remain in equilibrium. The modified Riks method allows you to find static equilibrium states during the unstable phase of the response.
You can use this method for cases where the load magnitudes are governed by a single scalar parameter. It is also useful for solving ill-conditioned problems such as limit load problems or almost unstable problems that exhibit softening. For more information, see Unstable collapse and postbuckling analysis, Section 6.2.4 of the ABAQUS Analysis User's Manual.
When you configure a static, Riks procedure, the step editor displays the Basic, Incrementation, and Other tabs. Settings you can configure with these tabbed pages include stopping criteria, the maximum number of increments, the arc increment length, and whether to account for geometric nonlinearity.
To configure settings on the Basic tabbed page:
In the Edit Step dialog box, display the Basic tabbed page.
(For information on displaying the Edit Step dialog box, see Creating a step, Section 14.9.2, or Editing a step, Section 14.9.3.)
In the Description field, enter a short description of what occurs during this analysis step. ABAQUS stores the text that you enter in the output database, and the text is displayed in the state block by the Visualization module.
Select an Nlgeom option:
Turn Nlgeom Off to perform a geometrically linear analysis during the current step.
Turn Nlgeom On to indicate that ABAQUS/Standard should account for geometric nonlinearity during the step. Once you have turned Nlgeom on, it will be active during all subsequent steps in the analysis.
Toggle on Include adiabatic heating effects if you are performing an adiabatic stress analysis. This option is relevant only for isotropic metal plasticity materials with a Mises yield surface. For more information, see Adiabatic analysis, Section 6.5.5 of the ABAQUS Analysis User's Manual.
Since the loading magnitude is part of the solution, you need a method to specify when the step is completed. Choose one or both of the following options:
Toggle on Maximum load proportionality factor to enter a maximum value for the load proportionality factor, . ABAQUS/Standard uses this value to terminate the step when the load exceeds a certain magnitude. For more information, see Proportional loading” in “Unstable collapse and postbuckling analysis, Section 6.2.4 of the ABAQUS Analysis User's Manual
Toggle on Maximum displacement to enter a maximum displacement value at a specific degree of freedom (DOF). You must also specify the Node Region that ABAQUS/Standard will monitor for finishing displacement. If this maximum displacement is exceeded, ABAQUS/Standard terminates the step.
To configure settings on the Incrementation tabbed page:
In the Edit Step dialog box, display the Incrementation tabbed page.
(For information on displaying the Edit Step dialog box, see Creating a step, Section 14.9.2, or Editing a step, Section 14.9.3.)
Select the Type option of your choice:
Select Automatic to allow ABAQUS/Standard to choose the arc length increment sizes based on computational efficiency.
Select Fixed to specify direct user control of the incrementation. ABAQUS/Standard uses an arc length increment that you specify as the constant increment size throughout the step. This method is not recommended for a Riks analysis since it prevents ABAQUS/Standard from reducing the arc length when a severe nonlinearity is encountered.
In the Maximum number of increments field, enter the upper limit to the number of increments in the step. The analysis stops if this maximum is exceeded before ABAQUS/Standard arrives at the complete solution for the step.
If you selected Automatic in Step 2, enter values for Arc length increment:
In the Initial field, enter the initial increment in arc length along the static equilibrium path in scaled load-displacement space, .
In the Minimum field, enter the minimum arc length increment, . If you enter zero, ABAQUS assumes a default value of the smaller of the suggested initial arc length or 105 times the total arc length.
In the Maximum field, enter the maximum arc length increment, . If this value is not specified, no upper limit is imposed.
In the Estimated total arc length field, enter the total arc length scale factor associated with this step, . If this entry is zero or is unspecified, ABAQUS/Standard assumes a default value of .
If you selected Fixed in Step 2, enter a value for the constant arc length increment in the Arc length increment field.
To configure settings on the Other tabbed page:
In the Edit Step dialog box, display the Other tabbed page.
(For information on displaying the Edit Step dialog box, see Creating a step, Section 14.9.2, or Editing a step, Section 14.9.3.)
Select the Matrix storage option of your choice:
Select Use solver default to allow ABAQUS/Standard to decide whether a symmetric or unsymmetric matrix storage and solution scheme is needed.
Select Unsymmetric to restrict ABAQUS/Standard to the unsymmetric storage and solution scheme.
Select Symmetric to restrict ABAQUS/Standard to the symmetric storage and solution scheme.
Toggle on Apply contact iteration solution technique to use contact iterations instead of regular severe discontinuity iterations to speed up computations. Contact iterations are effective for the solution of large, geometrically linear, small-sliding, frictionless static problems with many severe discontinuity iterations.
If you select this technique, enter the following values:
Adjustment factor for the number of solutions in any iteration. This value is a correction factor on the maximum number of right-hand-side solutions during any contact iteration.
Maximum number of contact iterations. This value specifies the maximum number of contact iterations allowed before new global matrix assemblage and factorization.
Click the arrow to the right of the Convert severe discontinuity iterations field, and select an option for dealing with severe discontinuities during nonlinear analysis:
Select Off to force a new iteration if severe discontinuities occur during an iteration.
Select On to estimate residual forces associated with severe discontinuities and check whether the equilibrium tolerances are satisfied. A solution may converge if the severe discontinuities (such as penetrations or tensile contact forces) are small. This option also changes some time incrementation parameters and uses different criteria to determine whether to do another iteration or to make a new attempt with a smaller increment size.
Select Propagate from previous step to use the value specified in the previous general analysis step. This value appears in parenthesis to the right of the field.
Click the arrow to the right of the Extrapolation of previous state at start of each increment field, and select a method for determining the first guess to the incremental solution:
Select Linear to indicate that the process is essentially monotonic, and ABAQUS/Standard should use a 1% linear extrapolation of the previous incremental solution to begin the nonlinear equation solution for the current increment.
Select None to suppress any extrapolation.
Toggle on Stop when region region name is fully plastic if “fully plastic” analysis is required with deformation theory plasticity. If you toggle on this option, enter the name of the region being monitored for fully plastic behavior.
The step ends when the solutions at all constitutive calculation points in the element set are fully plastic (defined by the equivalent strain being 10 times the offset yield strain). However, the step can end before this point if the maximum number of increments that you specified on the Incrementation tabbed page is exceeded.
If you selected Fixed time incrementation in the Incrementation tabbed page, you can toggle on Accept solution after reaching maximum number of iterations. This option directs ABAQUS/Standard to accept the solution to an increment after the maximum number of iterations allowed has been completed, even if the equilibrium tolerances are not satisfied. Very small increments and a minimum of two iterations are usually necessary if you use this option.
This approach is not recommended; you should it only in special cases when you have a thorough understanding of how to interpret results obtained in this way.
Toggle on Obtain long-term solution with time-domain material properties to obtain the fully relaxed long-term elastic solution with time-domain viscoelasticity or the long-term elastic-plastic solution for two-layer viscoplasticity. This parameter is relevant only for time-domain viscoelastic and two-layer viscoplastic materials.
Once you have finished configuring settings for the static, Riks step, click OK to close the Edit Step dialog box.