Riks Method

5/24/12

Abaqus Analysis User's Manual (6.10)

6.2.4 Unstable Unstable collapse and postbucklin postbuckling g analysis analys is Products: Abaqus/Standard Abaqus/CAE References

Procedures: overview, Section 6.1.1

“

”

Static stress analysis procedures: overview, Section 6.2.1

“

”

Introducing a geometric imperfection into a model, Section 11.3.1

“

”

*STATIC *IMPERFECTION Configuring a static, Riks procedure” in “Configuring general analysis procedures, Section 14.11.1 of the Abaqus/CAE User's Manual “

”

Overview

The Riks method: is generally used to predict unstable, geometrically nonlinear collapse of a structure; can include nonlinear materials and boundary conditions; often follows an eigenvalue buckling analysis to provide complete information about a structure's collapse; and and can be used to speed convergence of ill-conditioned or snap-through problems that do not exhibit instability. Unstable response

Geometrically nonlinear static problems sometimes involve buckling or collapse behavior, where the loaddisplacement response shows a negative stiffness and the structure must release strain energy to remain in equilibrium. Several approaches are possible for modeling such behavior. One is to treat the buckling response dynamically, thus actually modeling the response with inertia effects included as the structure snaps. This approach is easily accomplished by restarting the terminated static procedure ( Restarting an analysis, Section 9.1.1) 9.1.1 ) and switching to a dynamic procedure ( Implicit dynamic analysis using direct integration, Section 6.3.2) 6.3.2) when the static solution becomes unstable. In some simple cases displacement control can provide a solution, even when the conjugate load (the reaction force) is decreasing as the displacement increases. Another approach would be to use dashpots to stabilize the structure during a static analysis. Abaqus/Standard offers an automated version of this stabilization approach for the static analysis procedures (see Static stress analysis, Section 6.2.2; 6.2.2 ; Quasi-static analysis, Section 6.2.5; 6.2.5; Fully coupled thermal-stress analysis, Section 6.5.4; 6.5.4; or Coupled pore fluid diffusion and stress analysis, Section 6.8.1). 6.8.1). “

”

“

”

“

”

“

”

”

“

“

”

Alternatively, static equilibrium states during the unstable phase of the response can be found by using the prabhakar-pc:2080/v6.10/books/usb/ prabhakar-pc:2080/v6.10/books/usb/default.htm?sta default.htm?startat=pt03ch06s02a rtat=pt03ch06s02at03.html#usb-anl-ap t03.html#usb-anl-apostbuckling ostbuckling

1/7

5/24/12


“modified Riks method.” This method is used for cases where the loading is proportional; that is, where the load magnitudes are governed by a single scalar parameter. The method can provide solutions even in cases of complex, unstable response such as that shown in Fi Figu gure re 6.2. 6.2.4–1 4–1.. Figure Figu re 6.2.4–1 6 .2.4–1 Proportional loading with unstable response.

The Riks method is also useful for solving ill-conditioned problems such as limit load problems or almost unstable problem problemss that that exh exhibit bit sof softeni tening. The Riks method

In simple cases linear eigenvalue analysis ( Eigenvalue buckling prediction, Section 6.2.3) 6.2.3 ) may be sufficient for design evaluation; but if there is concern about material nonlinearity, geometric nonlinearity prior to buckling, or unstable postbuckling response, a load-deflection (Riks) analysis must be performed to investigate the problem further. “

”

The Riks method uses the load magnitude as an additional unknown; it solves simultaneously for loads and displacements. Therefore, another quantity must be used to measure the progress of the solution; Abaqus/Standard uses the “arc length,” l , along the static equilibrium path in load-displacement space. This approach provides solutions regardless of whether the response is stable or unstable. See the Modified Riks algorithm, Section 2.3.2 of the Abaqus Theory Manual, Manual, for a detailed description of the method. “

”

Proportional loading

If the Riks step is a continuation of a previous history, any loads that exist at the beginning of the step and are not redefined are treated as “dead” loads with constant magnitude. A load whose magnitude is defined in the Riks step is referred to as a “reference” load. All prescribed loads are ramped from the initial (dead load) value to the reference values values spec s peciified. fied. prabhakar-pc:2080/v6.10/books/usb/ prabhakar-pc:2080/v6.10/books/usb/default.htm?sta default.htm?startat=pt03ch06s02a rtat=pt03ch06s02at03.html#usb-anl-ap t03.html#usb-anl-apostbuckling ostbuckling

2/7

5/24/12


The loading during during a Riks Riks step step is always propo pro portional. rtional. The current curre nt load magnitu magnitude de,,

, is defi de fined ned by

where is the “dead load,” oad ,” is the reference load vector, and is the “load “load proportional prop ortionaliity factor.” The load proporti proportion onal aliity factor actor is is fou foun nd as part part of the the solu solution tion.. Abaqu Abaqus/St s/Stan andard dard pri prin nts out out th the curre curren nt val valu ue of the the load load proporti proportion onal aliity factor actor at each incremen crement. t. Incrementation

Abaqus/Standard Abaqus/Standard uses Newton Ne wton's 's method (as described in Static stress analysis, Section 6.2.2) 6.2.2) to solve the nonlinear equilibrium equations. The Riks procedure uses only a 1% extrapolation of the strain increment. “

You provi pro vide de an initi initial al incre increment ment in arc length length along the static static equi eq uili librium brium path, path, initi initial al load oa d propo pro portionali rtionality ty factor, factor, , is computed computed as

”

, when when you defi de fine ne the step. ste p. The

where is a user-spe user- specif cified ied total tota l arc length length scale scale factor factor (typically (typically set set equal eq ual to 1). This This value value of is used during during the first first iteration of a Riks Riks step. ste p. For Fo r subsequent subseq uent iterations terations and incre increments ments the value value of is computed co mputed automatically, automatically, so you have no no control co ntrol over the loa load d magni magnitude. The value value of is part of the solution. solution. Mini Minimu mum m and maximu maximum m arc length incre increments, ments, and , can ca n be used to control co ntrol the automatic incrementati ncrementation. on. *STATIC,, RIKS *STATIC

Input Inp ut File File Usage :

Abaqus/CAE Usage: Step module: Create Step: General: Static, Riks

Direct use userr control of the increment increment size size is also provided; provided ; in this this case case the the incremental incremental arc length, length, , is is kept kept constant. 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. *STATIC,, RIKS, DIRECT *STATIC

Input Inp ut File File Usage :

Abaqus/CAE Usage: Step module: Create Step: General: Static, Riks : Incrementation: Type: Fixed

Ending a Riks analysis step

Since the loading magnitude is part of the solution, you need a method to specify when the step is completed. You can ca n specif spe cify y a maximu maximum m value value of the load propo pro portionali rtionality ty factor, factor, , or a maximu maximum m displacement displace ment value value at a specified degree of freedom. The step will terminate when either value is crossed. If neither of these finishing conditions is specified, the analysis will continue for the number of increments specified in the step definition (see Procedures: overview, Section 6.1.1). 6.1.1). “

”

Bifurcation

prabhakar-pc:2080/v6.10/books/usb/ prabhakar-pc:2080/v6.10/books/usb/default.htm?sta default.htm?startat=pt03ch06s02a rtat=pt03ch06s02at03.html#usb-anl-ap t03.html#usb-anl-apostbuckling ostbuckling

3/7

5/24/12


The Riks method works well in snap-through problems—those in which the equilibrium path in load-displacement space is smooth and does not branch. Generally you do not need take any special precautions in problems that do not exhibit exhibit branching b ranching (bifurcat (bifurcatiion). Snap-through buckling analysis of circular arches, Section 1.2.1 of the Abaqus Example Problems Manual, Manual, is an example of a smooth snap-through problem. “

”

The Riks method can also be used to solve postbuckling problems, both with stable and unstable postbuckling behav behaviior. Howev However, er, th the exact exact postb postbu uckli ckling problem problem cann cannot be anal analy yzed direct directlly due due to th the discon disconti tin nuous ous respon response se at the point of buckling. To analyze a postbuckling problem, it must be turned into a problem with continuous response instead of bifurcation. This effect can be accomplished by introducing an initial imperfection into a “perfect” geometry so that there is some response in the buckling mode before the critical load is reached. Introducing geometric imperfections

Imperfections are usually introduced by perturbations in the geometry. Unless the precise shape of an imperfection is known, an imperfection consisting of multiple superimposed buckling modes must be introduced ( Eigenvalue buckl bu ckliing predi predicti ction on,, Section 6.2.3). 6.2.3 ). Abaqus allows you to define imperfections; see Introducing a geometric imperfection into a model, Section 11.3.1. 11.3.1 . “

”

“

”

In this way the Riks method can be used to perform postbuckling analyses of structures that show linear behavior prior prior to to (bif (bifurcati rcation) on) buckl buckliing. An exam exampl plee of thi this meth method od of introdu troduci cin ng geometr eometriic im imperfecti perfection onss is is presen presented ted in in Buckling of a cylindrical shell under uniform axial pressure, Section 1.2.3 of the Abaqus Benchmarks Manual. Manual. “

”

By performing a load-displacement analysis, other important nonlinear effects, such as material inelasticity or contact, can be included. In contrast, all inelastic effects are ignored in a linear eigenvalue buckling analysis and all contact conditions are fixed in the base state. Imperfections based on linear buckling modes can also be useful for the analysis of structures that behave inelastically prior to reaching peak load. Introd In troducing ucing load loading ing imp imperfe erfections ctions

Perturbations in loads or boundary conditions can also be used to introduce initial imperfections. In this case fictitious “trigger” loads can be used to initiate the instability. The trigger loads should perturb the structure in the expected buckling modes. Typically, these loads are applied as dead loads prior to the Riks step so that they have fixed magnitudes. The magnitudes of trigger loads must be sufficiently small so that they do not affect the overall postbu postbuckl ckliing solu solution tion.. It is you yourr respon responsi sibi billity to choose choose appropr appropriiate mag magn nitudes tudes and and locati location onss for for such such fictiti ctitiou ouss loads; Abaqus/Standard Abaqus/Standard does not check that they they are reasonable. Obtaining a solution at a particular load or displacement value

The Riks algorithm cannot obtain a solution at a given load or displacement value since these are treated as unknowns—termination occurs at the first solution that satisfies the step termination criterion. To obtain solutions at exact values of load or displacement, the solution must be restarted at the desired point in the step ( Restarting an analysis, Section 9.1.1) 9.1.1) and a new, non-Riks step must be defined. Since the subsequent step is a continuation of the Riks analysis, the load magnitude in that step must be given appropriately so that the step begins with the loading continuing to increase or decrease according to its behavior at the point of restart. For example, if the load was increasing at the restart point and was positive, a larger load magnitude than the current magnitude should be given in the restart step to continue this behavior. If the load was decreasing but positive, a smaller magnitude than the current magnitude should be specified. “

”


4/7

5/24/12


Restrictions

A Riks analysis is subject to the following restrictions: A Riks step cannot be followed by another step in the same analysis. Subsequent steps must be analyzed by using the restart capability. If a Riks analysis includes irreversible deformation such as plasticity and a restart using another Riks step is attempted while the magnitude of the load on the structure is decreasing, Abaqus/Standard will find the elastic unloading solution. Therefore, restart should occur at a point in the analysis where the load magnitude is increasing if plasticity is present. For postbuckling problems involving loss of contact, the Riks method will usually not work; inertia or viscous damping forces (such as those provided by dashpots) must be introduced in a dynamic or static analysis to stabilize the solution. Initi In itial al cond c onditi itions ons

Initial values of stresses, temperatures, field variables, solution-dependent state variables, etc. can be specified; Initial conditions in Abaqus/Standard and Abaqus/Explicit, Section 30.2.1, 30.2.1 , describes all of the available initial conditions. “

”

Boundary conditions

Boundary conditions can be applied to any of the displacement or rotation degrees of freedom (1–6) or to warping degree of freedom 7 in open-section beam elements ( Boun Boundary dary conditi conditions ons in Abaqus/Standard Abaqus/Sta ndard and Abaqus/Explicit, Section 30.3.1). 30.3.1 ). Amplitude definitions ( Amplitude curves, Section 30.1.2) 30.1.2) cannot be used to vary the magnitudes of prescribed boundary conditions during a Riks analysis. “

”

“

”

Loads

The following loads can be prescribed in a Riks analysis: Concentrated nodal forces can be applied to the displacement degrees of freedom (1–6); see Concentrated loads, Section 30.4.2. 30.4.2. “

”

Distributed pressure forces or body forces can be applied; see Di Distribu stributed ted load oads, s, Section 30.4.3. 30.4.3. The distributed load types available with particular elements are described in Part VI, Elements Elements.” .” “

”

“

Since Abaqus/Standard scales loading magnitudes proportionally based on the user-specified magnitudes, amplitude references are ignored when the Riks method is chosen. If follower loads are prescribed, their contribution to the stiffness matrix may be unsymmetric; the unsymmetric matrix storage and solution scheme can be used to improve computational efficiency in such cases (see Procedures: overview, Section 6.1.1). 6.1.1). “

”

Predefined fields

Nodal temperat temperatu ures can can be specif specified (see (see Predefined fields, Section 30.6.1). 30.6.1 ). Any difference between the applied “

”


5/7

5/24/12


and initial temperatures will cause thermal strain if a thermal expansion coefficient is given for the material ( Thermal expansion, Section 23.1.2). 23.1.2 ). The loads generated by the thermal strain contribute to the “reference” load specified for the Riks analysis and are ramped up with the load proportionality factor. Hence, the Riks procedure can analyze postbuckling and collapse due to thermal straining. “

”

The values of other user-defined field variables can be specified. These values affect only field-variable-dependent material properties, if any. Since the concept of time is replaced by arc length in a Riks analysis, the use of properti properties es th that ch chang ange due due to chan chang ges in in temperat temperatu ures an and/or field eld vari variabl ables es is is not not recom recomm mended. ended. Material options

Most material models that describe mechanical behavior are available for use in a Riks analysis. The following material properties are not active during a Riks analysis: acoustic properties, thermal properties (except for thermal expansion), expansion), mass diffu diffusion sion properti prope rties, es, electrical properti prope rties, es, and pore flui fluid d flow properties. prop erties. Materials with with history history dependence can be used; however, it should be realized that the results will depend on the loading history, which is not known in advance. The concept of time is replaced by arc length in a Riks analysis. Therefore, any effects involving time or strain rate (such as viscous damping or rate-dependent plasticity) are no longer treated correctly and should not be used. See Part V, Materials Materials,” ,” for details on the material models available in Abaqus/Standard. “

Elements

Any of the stress/displacement elements in Abaqus/Standard (including those with temperature or pressure degrees of freedom) can be used in a Riks analysis (see Choosing the appropriate element for an analysis type, Section 24.1.3). 24.1.3 ). Dashpots should not be used since velocities will be calculated as displacement increments divided by arc length, which is meaningless. “

”

Output

Output options are provided to allow the magnitudes of individual load components (pressure, point loads, etc.) to be prin printed or or to to be wri writt tten en to th the resu results file. Th The curre curren nt val valu ue of the the load load proport proportiional onaliity factor, actor, LPF, wil will be gi given automatically with any results or output database file output request. These output options are recommended when the Riks method is used so that load magnitudes can be seen directly. All of the output variable identifiers are outlined in Abaq Abaqus/Standard us/Standard output variable variable identi dentifi fiers, ers, Section 4.2.1. 4.2.1. “

”

Input file template

*HEADING … *INITIAL CONDITIONS

Data lines lines to define initia initiall condition conditionss *BOUNDARY

Data lines lines to specify specify zero-valued boundary conditions conditions ** *STEP, *STEP , NLGEOM *STATIC *CLOAD and/or *DLOAD and/or *TEMPERATURE prabhakar-pc:2080/v6.10/books/usb/ prabhakar-pc:2080/v6.10/books/usb/default.htm?sta default.htm?startat=pt03ch06s02a rtat=pt03ch06s02at03.html#usb-anl-ap t03.html#usb-anl-apostbuckling ostbuckling

6/7

5/24/12


Data lines lines to specify specify preload (dead load), load), *END STEP ** *STEP, *STEP , NLGEOM *STATIC, *STATIC , RIKS

Data line line to define define incrementation incrementation and stopping criteria criteria *CLOAD and/or *DLOAD and/or *TEMPERATURE

Data lines lines to specify specify reference loadi loading ng , *END STEP


7/7

Riks Method

Recommend Documents