Part 1. Theoretical Background #
AutoMcr is an application used in the Steel Design module to calculate the elastic critical moment (Mcr). Mcr is required in the calculation of lateral torsional buckling resistance. AutoMcr creates an individual finite element submodel of each steel design element, for which it determines the Mcr value by solving an eigenvalue problem. The submodel is built-up of special beam finite elements only with those degrees of freedom that are relevant for lateral torsional buckling:
- v lateral displacement, in the direction of local y axis;
- θx torsion: rotation about beam axis / local x axis;
- θz rotation about weak axis / local z axis;
- w warping.
When creating the submodel, the program automatically identifies lateral supports, which can be edited by the user. The rigidity components of the support, indexed according to the local coordinate system of the submodel: Ry, Rxx, Rzz, Rw.
The AutoMcr is based on the same theory as the LTBeam program, of which further information can be read in the following article: Yvan Galea: Moment critique de deversement elastique de poutres flechies presentation du logiciel ltbeam .
This Guide has two main goals. In Part 1, examples demonstrate the possibilities and limits of AutoMcr, while helping users to properly use the program. Part 2 is a summary of verification models, in which results of AutoMcr are compared to literature and to other programs. For basics of the AutoMcr method and to learn how to use it, check AxisVM13 User’s Manual: 6.6.2. Steel beam design based on Eurocode.
The AutoMcr is capable of analysing straight elements with a cross section symmetric at least about the weak axis. Moreover, it can handle:
- elements with variable cross-section;
- cantilevers: no need to define if it is a cantilever or not, as in AxisVM12;
- eccentric load: distance from the weak axis, one value for all load cases analysed at a time;
- eccentric support conditions: defined individually for each support.
The AutoMcr method handles only continuous elements, therefore it splits up design members in the following two cases:
- tapered beam: when part of the beam has variable cross-section, the rest is constant;
- elements with intermediate pin.
- Lateral supports
With default settings, the Auto Mcr method automatically determines the lateral supports of the designed member; which will be detailed in the following. The program finds not only the supports defined earlier in the main model, but also the elements that are connected to the designed member. These connected elements may be:
- truss, beam of rib elements;
- surface elements;
- rigid elements, node-to-node interface elements.
Based on the properties of these elements, lateral support stiffness values are estimated by the program. This is detailed in Table -.
In the Design Parameters window (Fig. 1) the lateral supports may be edited after pressing the […] button which is below the Auto Mcr setting and next to the Lateral Supports caption. The Lateral supports window will appear (Fig. 2), in which the assumed lateral supports are visible. These supports are dependent on the settings of the AutoMcr method:
Automatic default setting; see Table 2-3.
Estimated from kz, kw Based on the user-defined kz and kw parameters, similarly to AxisVM 12, lateral support location and stiffness values are estimated. For details see Table 1.
Fork supports at both ends In the end of the designed member, lateral supports are assumed with rigid Ry and Rxx components. If the user-defined cantilever option is checked, then supports appear only on one end with rigid Ry, Rxx and Rzz components.
User defined Only the user-defined supports are considered defined in the Lateral supports window.
Figure 1: Design Parameters window
Figure 2: Lateral Supports window
Table 1: Lateral supports determined based on kz and kw
Table 2: Lateral supports determined by the program automatically – supports and connected line elements
* if the designed member is not braced in x-y plane; otherwise Ry = 0 kN/m
Table 3: Lateral supports determined by the program automatically – further connected elements
α – smallest angle between the axis of designed member + the axis of connected member / surface plane (0÷90°)
β – smallest angle between the major axis of designed member + the axis of connected member / surface plane (0÷90°)
For example, when designing an I beam these angles for the bracing elements:
Part 2. Examples #
In the girders below, lateral torsional buckling is prevented by using fork supports in the ends and by laterally connected beams in two intermediate points of the girder.
Figure 3: Girders with stiffening beams and connection detail (source: )
The goal of this example is to demonstrate:
- how to determine the support stiffness provided by the connected beams;
- comparing Mcr obtained by AutoMcr with those of shell models and the LTBeam program.
The structure in the following book served as a basis for this example, which gives guidance in determining the support stiffness provided by adjacent beams: Teil 2 – Stabilität und Theorie II. Ordnung .
- Cross-section [mm]:
- girder: in order to be able to compare results with shell finite element models, welded I section similar to IPE 300: web: 300*7, flanges: 150*11 mm;
- connected beam UPE 140;
- girder: l=6m;
- connected beam: a=3m;
- Loading: distributed force along the whole girder or point load in the middle of the girder; applied in the geometric centre or on top of the flange;
- Support condition: supports in the ends of the girder according to Figure 3 (either of the two girders may move laterally)
Name of AxisVM models:
- Beam finite element model with AutoMcr: I. Girder – beam finite element model.axs
- Shell finite element model as an eigenvalue problem: I. Girder – shell finite element model.axs
Lateral support stiffness
In the ends of the beams, there are fork supports. In AxisVM13, when creating the AutoMcr submodel, the program automatically adopts the supports defined earlier in Elements >> Nodal supports. These supports of the AutoMcr submodel can be seen in the table at Design Parameters >> Lateral supports. For the girder, these adopted supports can be seen in Figure 4, of which the lateral Ry and rotational Rxx stiffness components are stiff.
Figure 4: Defining lateral supports in AxisVM13
In the table above, additionally to the adopted supports (Support form model), the connected beams also provide support (Connecting element) against lateral torsional buckling. The program automatically gives approximate values for the Ry and Rxx components of such a support:
- Ry = 1010 kN/m if the analysed member is braced in local x-y plane; otherwise: Ry = 0 kN/m;
- Rxx=2*EI/a based on the length (a) and the inertia (I) of the connecting member.
It is the User’s responsibility to define this stiffness value accurately, if needed. To calculate the stiffness provided by the connected beams,  gives the following recommendation: the rotational support stiffness (Rxx) may easily be calculated based on the stiffness of the connected beam (EI/a). The stiffness values may be determined by the following two formulas, based on the deformation of the structure:
Girders exhibit lateral displacements and rotate in the same direction. The connected beams do not provide any lateral support.
Rxx = 6*EI/a =
= 6*21000kN/cm2 * 599.6cm4 / 3m
= 2520 kNm/rad
Ry = Rzz = Rw =0
Figure 5: Possible deformation of the girder structure: non-symmetric case (source: )
Girders do not exhibit lateral displacements and rotate in the opposite direction. The connected beams provide some lateral support. Auto_Mcr_girder_deformation_2.png
Rxx = 2*EI/a =
= 6*21000kN/cm2 * 599.6cm4 / 3m
= 840 kNm/rad
Ry > 0
Rzz = Rw =0
Figure 6: Possible deformation of the girder structure: symmetric case (source: )
In reality, semi-rigid connections and the distortions of the girder may lower the above support stiffness values, therefore to stay on the safe side, the program uses the second case. In the following comparison, both cases will be presented, in the second case by neglecting Ry.
Comparison of results
The obtained Mcr results are compared to results of shell models created in AxisVM13, and of the LTBeam program, which works on the same basis as AutoMcr. The models created in LTBeam (v1.0.10) have the same settings. The differences in the obtained results are due to the used numerical algorithm and to the differences in the discretisation.
The shell models in AxisVM13 were created with the help of the Edit >> Convert beams to shell model function. After defining the load, by solving an eigenvalue problem (Buckling tab), a load factor is obtained. Mcr can be calculated by multiplying the load factor with the maximal moment along the beam. Compared to beam models, shell models are capable of a more detailed and precise modelling, thus the obtained Mcr is more accurate. Another advantage of shell models is that there is no need to create a sub-model, and thus there is no error caused by defining lateral supports. The disadvantage is that the modelling is more complex and more time consuming. The calculation time for AutoMcr is about a 100 times lower that for an appropriate shell model. To avoid local deformations in the shell model, the web of the girder is stiffened by rigid elements at the intersection of the beams (a more accurate modelling of the stiffening plate is neglected). The obtained lowest eigenform is the symmetric case, while the second is the non-symmetric case (Figure 7).
Figure 7: Eigenforms of shell finite element models; left: symmetric case; right: non-symmetric case [mm]
In Table 4, the Δ columns show the difference of the AutoMcr results (MAutoMcr) compared to either of the other methods (Mcr), based on this formula: Δ = (MAutoMcr – Mcr) / Mcr.
Table 4: Comparison of results
Comparing the results to the LTBeam program, the AutoMcr method is accurate. Furthermore, it can be concluded, that the results obtained by the shell finite element model and the beam finite element model with AutoMcr correspond well, thus the applied support stiffness values are accurate enough.
Part 3. Verification #
In this part, the verification of the AutoMcr method is summarized. The calculated Mcr values are compared to those of other methods and programs, among which is the LTBeam program that is based on the same theoretical background as AutoMcr. In the first section, the LTBeam and shell models are taken from the verification documentation of LTBeam: Yvan Galea: LTBeam – Report on Validation Tests . Afterwards, comparison is made with the ENV  analytic formula. Lastly, the differences of the AutoMcr method in AxisVM12 and 13 are summarized.
The error (Δ) of the AutoMcr results (MAutoMcr) compared to either of the other methods (Mcr) was calculated based on this formula: Δ = (MAutoMcr – Mcr) / Mcr.
- Validating with LTBeam program and shell models
Ansys shell models
Based on Chapter 2 of .
This section presents simple examples of all the types of models that can be calculated with AutoMcr. Results are compared to those of the LTBeam program and of shell models in Ansys  and are presented in Table 5Table 6. Mcr values are only -4÷3% different, which is a very good result.
Name of Axis model: LTBeam Validation – Chapter 2 – #.axs (where # is the number of the example)
Table 5: Comparison of results I.
Table 6: Compariso of results II.
Based on Chapter 5 of .
The analysed beam has variable web height (hw1÷hw2), fork supports in the end points, and end moments (M1 and M2). The results are generally +2% and maximum -9% different from the results of LTBeam and Finelg , the reason of which lies in the different discretisation of the sections. These differences are negligible compared to the general uncertainty of modelling variable cross-sections.
Name of Axis model: LTBeam Validation – Chapter 5 – Variable cross-section.axs
Table 7: Comparison of results of beam with variable cross section
- Basic Cases with the analytic expression in ENV
In order to determine Mcr values, AxisVM program has long been using the so-called “3 factor formula”, which can be found in the pre-standard of the Eurocode  (in the following referred to as ENV). Additionally to the 3 C factors, the formula uses the kz and kw effective length factors. Recommended values for all these factors may be found in several literatures for basic cases only, and in some cases giving different results. To calculate the C1 factor, Lopez et al. proposed a simple analytic formula that AxisVM program implemented. This formula was calibrated by numerical results in several support conditions and load cases.
In Table 6, results are summarized and compared for the AutoMcr method and for the ENV formula based on factors of several sources. All the examples are beams supported on the ends, loaded and supported in their shear centre, with a double- or single-symmetric I cross-section and various effective length factors.
In line with , in the ENV formula, kz and kw are assumed to be equal. Additionally to pinned and fixed beams,  provides factors for a third “semi-fixed” support condition: when k values are taken as 0.7. This provides less information about the support condition, than what needs to be defined in AutoMcr. Therefore, in the following, this case was modelled with three different settings. Logically, the k=0.7 corresponds to a beam, that is fully-fixed on one end and pinned on the other; for this setting, the smallest possible Mcr value is included in the table. In the other two settings either kz or kw is 0.5, the other is 1, which are generally used in practice. Table 5 summarizes these support conditions (the support components not included in the table are assumed to be zero for the AutoMcr method).
Table 8: Lateral support conditions as defined in for the different methods
End moments only
Figure 8. End moments only
Cross-section: Symmetric: welded (same plate size as IPE 300)
Name of Axis model: Basic cases – End moments – Symmetric cross-section.axs
Table 9: Comparison of analytic and numerical results, end moments only
It can be seen in Table 9, that the various methods give significantly different results. In all cases, the results of AutoMcr and LTBeam are very close.
- For pinned beams, results are always very similar for all methods.
- For fixed beams, the results from the ENV method combined with C1 factor based on the Lopez formula  is closest to the AutoMcr results, mainly if Ψ>0.
- The differences between the methods for the “semi-fixed” cases lie in the different definition of the support condition.
Name of Axis model: Basic cases – Transverse loading – Symmetric cross-section.axs
Table 10: Comparison of analytic and numerical results, transverse loading
- Differences between AxisVM Version 12 and 13
In AxisVM12, when defining the sub-model, the support conditions are assumed based on the user defined kz and kw values. The obtained Mcr values are very similar to the results in AxisVM13 in the basic cases (k=0.5 or k=1) but are less accurate if kz≠kw.
A further important difference is, that in version 13, for a safe design, in case of a simple beam with fixed end-supports, AutoMcr automatically assumes that Ry=Rxx=Rzz=1010, while the user shall determine Rw. In version 12, if kz=kw=0.5, Rw is also assumed to be rigid.
Table 11: Lateral support conditions
The AutoMcr method of AxisVM13 is numerically more precise in version 12. The Mcr results are maximum ±10% different. When first opening a model in version 13, that was created and saved in version 12, the support conditions are the same as they were in version 12, but the Mcr values are calculated by the more precise algorithm. In the Steel Design Parameters window, such a model will appear to have the Mcr method: „AutoMcr_v12”. Conversion of such models are recommended, and the redefinition of lateral support conditions, to facilitate a more accurate design.
 Yvan Galea: Moment critique de deversement elastique de poutres flechies presentation du logiciel ltbeam, CTICM, www.cticm.com, 2003 (in French)
 Stahlbau: Teil 2 – Stabilität und Theorie II. Ordnung, 10.4 Stabilisierung durch behinderung der verdrehungen, szerző: Rolf Kindmann, Ernst&Sohn, 4th edition, pp. 336-338., 2008 (in German)
 Yvan Galea: LTBeam – Report on Validation Tests, CTICM, July 2002, www.cticm.com
 ENV 1993-1-1: Appendix F
 Ádány, Dulácska, Dunai, Fernezely, Horváth: Acélszerkezetek, 1. Általános eljárások, Tervezés az Eurocode alapján, 2006 (in Hungarian)
 López, Yong, Serna: Lateral-torsional buckling of steel beams: A general expression for the moment gradient factor. Proceedings of the International Colloquium of Stability and Ductility of Steel Structures, D. Camotim et al. Eds., Lisbon, Portugal, September 6-8, 2006.
 Access Steel: NCCI: Elastic critical moment for lateral torsional buckling, SN003a-EN-EU, 2008
 Braham M. – “Le déversement élastique des poutres en I à section monosysmétrique soumises à un gradient de moment de flexion” – Revue Construction Métallique n°1-2001 – CTICM (in French)