1. Introduction
In order to understand the essence of the work of steel bar structures, it is worth going back to the beginnings related to their mechanics. In the case of a slender bar, the most typical failure mechanism is loss of stability. In general, the stability of a structure can be defined as a phenomenon in which a slight disturbance causes a very large change. Depending on the effect of the disturbance, the condition of the structure can be defined as stable or unstable. Usually, such a definition is referred to as stability in Lyapunov’s meaning. Another approach to the analysis of the equilibrium states of the structure is based on the potential energy analysis. It is assumed that the structure is in equilibrium when potential energy reaches a minimum. The necessary condition for the occurrence of minimum potential energy is the disappearance of its first variation. The nature of the equilibrium path is determined by higher order variations of potential energy. The type of bifurcation point can be defined based on the analysis of the second variation of potential energy. This approach was proposed by Stepan Timoshenko [
1], who was considered to be the father of modern engineering mechanics. In the classical theory of stability, the type of equilibrium path can be illustrated by a ball moving on a differently shaped surface as stable, unstable and neutral. In problems with many degrees of freedom, there is also another case, called the metastable point of equilibrium. This state is characterised by temporary equilibrium and local minimum of potential energy. Note that the nature of a bifurcation point can be defined by the higher-order variation of potential energy and can be symmetric or asymmetric.
In engineering practice, i.e., in the current codes [
2], the most common concept is the one of critical stress, which was proposed originally by the Swiss mathematician Leonhard Euler, who in 1757 derived the formula of critical load. This formula describes the compressive load at which, due to initial geometrical imperfection, a slender column suddenly bends. This phenomenon is called buckling and can take different forms with respect to the bar cross-sectional characteristic. In case of slender, thin-walled bar members, a spatial loss of stability is observed. In the profile, warping occurs, combined with elastic couplings involving bending and torsion. Vasilii Zakharovich Vlasov developed a torsion theory in which restrained warping is included [
3]. This theory is also called “warping torsion” or “non-uniform torsion”. Note that the loss of stability is limited by the occurrence of initial geometrical imperfections. In typical structural elements, geometrical imperfections can take global and local forms. Global imperfections refer mainly to the bending or twisting of rod axes. On the other hand, local imperfections may occur in the form of deviations of the dimensions of individual cross-sectional walls, or the shape of the contour of the cross-section associated with the local deformation of walls, or the symmetric and asymmetric “opening” and “closing” of the cross-section. These imperfections are particularly dangerous when their shape coincides with particular forms of the local or distortion buckling mode [
4]. Taking into account the actual geometrical imperfections, both global and local, is possible using the finite element method [
5,
6]. Issues connected with the accuracy of the computational models which take into account the initial geometrical imperfections and residual stresses were dealt with in [
7]. The influence of local geometrical imperfections of thin-walled sections on the stability of columns and short beams was analysed in [
8]. There are many methods of introducing geometrical imperfections into the FEM numerical model. One of the simplest methods is the equivalent load method in which the force induces in the numerical model deformation corresponding to initial geometrical imperfections [
9,
10]. This method is especially suitable for generating global geometrical imperfections, but its disadvantage is the creation of additional stresses, which in fact are not accompanied by the formation of geometrical imperfections. An alternative method of imperfection modelling is a concept based on the equivalent bar stiffness proposed by [
11]. In [
12] the author shows that geometrically exact beam finite elements can be successfully employed to determine, accurately, the behaviour of steel I-section beams undergoing large displacements and finite rotations. For this purpose, a two-node, geometrically exact, beam element was presented and validated, which can handle arbitrary initial configurations (e.g., curved configurations), plasticity, geometric imperfections and residual stresses. A very realistic model of construction imperfections, both global and local, can be obtained by introducing a distorted geometry of the FE mesh by measuring the actual imperfections “in situ”. This approach requires the use of highly advanced measuring devices (such as 3D scanners) to create 3D imperfection maps and convert them to FEM programs. Unfortunately, in the case of civil engineering structures, this method is still rarely used because of a relatively high cost. It has one more drawback, namely it allows the creation of an initial imperfections pattern only for a given structural element. Therefore, it only covers one particular case. Due to the high randomness of imperfections, the results obtained on the basis of deterministically defined imperfections are not always reliable and do not allow potentially dangerous configurations to be identified. An alternative method of taking into account imperfections may be a method of shaping deformed geometry as a linear superposition of buckling modes obtained from a linear solution of the eigenvalue problem [
13]. In this approach, the amplitude of buckling modes can be scaled by proportional coefficient factors determined arbitrarily based on standard tolerances or on a limited number of measurements of real imperfections.
Nonlinear FE formulations based on von Kármán kinematics with initial geometrical imperfections are taken into account in [
14]. The authors proposed a predictor–corrector reduced-order modelling method for the nonlinear buckling analysis of thin-walled structures with initial geometrical imperfections based on Koiter theory. The main innovation of this method is the implementation of initial geometrical imperfection, instead of equivalent loads. This makes the Koiter–Newton reduced-order modelling method applicable for pre- and post-buckling, even for relatively large imperfection amplitudes.
The idea of increasing buckling capacity with modal geometrical imperfections was discussed by [
15]. They note that not all imperfections are detrimental to buckling problems. A recent design trend called “modal nudging” takes advantage of geometrical imperfections to change the primary structure to equilibrium paths with greater bearing capacity. In this article, the authors proposed an automated procedure based on a Koiter-inspired reduced-order model that allows the equilibrium paths of structures with geometrical imperfections to be estimated rapidly according to the first linearized buckling modes. Fast random search is used to select a geometrical imperfection that pushes the structure onto a favourable equilibrium path. In this case, a geometrical imperfection may play a beneficial role. “Modal nudging” aims to improve the buckling response of slender elastic structures by introducing a small change to the base geometry, capable of pushing the structure onto equilibrium paths with higher bearing capacity, both for unstable and stable equilibrium pathways. However, the authors note that it is not enough to guarantee a response that is insensitive to imperfections. In fact, designers should be aware that small random imperfections could push the structure back onto unstable paths. The “modal nudging” was also used by [
16], as a design strategy that intends to improve the load-carrying capacity of a structure, by nudging the structure to follow a specific equilibrium path. They proved that “modal nudging” by small deterministically chosen geometric changes lead the steel frames onto their higher load-carrying paths. This is why it is so important in engineering practice to study different imperfection configurations.
The nonlinear finite element model used to simulate equilibrium paths is presented in [
17]. Special interest is focused on the identification of critical points which contribute to the understanding and quantification of stability processes in a nonlinear system. The paper discusses the methods of stochastic and sensibility analyses which are frequently applied to the assessment of the safety and reliability of supporting structural systems. It is pointed out that the failures of the slender structural system caused by stability loss very often causes its decrease in safety and reliability. Therefore, reliable methods of stability analysis are important instruments to ensure the general effectiveness of modern structural systems.
The snap-through phenomenon is very important in the stability analyses of steel members with initial imperfection. It is also observed in real engineering structures as well as during experiments.
This issue is studied by [
18] on the example of geometrically nonlinear truss structures. By tracing their equilibrium paths, using an arc length approach within a finite element analysis, snap-through behaviour is demonstrated. The snap-through phenomenon of the Mises truss taking into account spatial asymmetric and symmetric buckling with an out-of-plane lateral linear spring is also investigated in [
19]. In the paper, the tangent stiffness matrix, for the primary equilibrium trajectory, is obtained in a closed-form and eigenvalues of the stiffness matrix are also derived analytically. They are used for the detection, classification and sequencing of critical points, in the primary equilibrium path of the Mises truss. The kinematic approach was used to formulate a spatial truss finite element, in which the expressions of the internal forces and the tangent stiffness are dependent on the direction and the current length of the bar element. The main innovation of this paper is the introduction of an additional DOF in the Mises truss, which allows a more complex scenario of equilibrium paths and instabilities to be studied. It is important to highlight, in particular, the simplicity and efficiency of the test function based on monitoring the eigenvalue signals. Finally, the remarkably accurate numerical results obtained, in perfect agreement with those of the derived analytical results, make the developed truss element very reliable and robust.
Recently, machine learning algorithms have been used successfully in various fields of structural engineering. This modern technological tool utilization plays a key role in developing robust frameworks that can enhance the decision-making process within the construction industry. For example, machine and deep learning-based technological tools were used in [
20] for providing a robust, fast, accurate, and flexible forecasting framework for the prediction influence of concrete mix properties on the shear strength of slender structured concrete beams without stirrups. In addition, in [
21], machine learning algorithms were used in order to predict the shear strength accurately, which enable the automation to be applied in steel buildings.
2. Problem Formulation
The research problem posed in this work concerns the influence of the loss of stability on the load-bearing capacity of the bar structure, widely used in civil engineering, with respect to a complex geometrical imperfection pattern. Due to the high complexity of the loss-of-stability phenomenon, the routine approach to designing slender structures, as proposed by professional computer programs, can be dangerous or can lead to non-optimal solutions. Therefore, the model structure proposed in this paper can be very useful in order to understand the essence of the influence of complex imperfection configuration on the post-critical behaviour of steel bar structures. It is worth noting that the ultimate load-bearing capacity of the structure caused by the loss of stability depends on many factors, above all, on the unavoidable geometrical deviations of the structure, understood as a deviation from the ideal calculation model.
In order to study the pre- and post-buckling behaviour of bar structures with geometrical imperfections, nonlinear algebraic equilibrium equations were developed. The potentially most dangerous configurations of imperfections were taken into account using a method of shaping deformed geometry as a linear superposition of the buckling modes obtained from a linear solution of the eigenvalue problem. Next, the analysis of the equilibrium states of the structure based on the concept of the minimum of potential energy was performed. Moreover, the so-called snap-through phenomenon based on transition from the local to the global minimum of potential energy, commonly associated with the release of internal energy due to a sufficiently large impulse [
22], was investigated. It was assumed that this phenomenon is caused by the specific configuration of initial geometric imperfections and is associated with the so-called metastable bifurcation point. Despite the fact that snap-through is a fully reversible phenomenon, multiple snap-throughs in the case of variable loads may cause material fatigue and, consequently, its brittle failure.
The advantage of the proposed model is the possibility of tracing various forms of deformation corresponding to the buckling modes and it enables the easy and clear description of the equilibrium equations. At the same time, the model enables the control of imperfections, which are developed in accordance with modal forms obtained from the solution of the linear eigenvalue problem. Therefore, it seems that the proposed model has full equivalents in engineering practice, and performs the ability to track various configurations of imperfections and multi-bifurcation post-critical equilibrium paths.
4. Numerical Examples
From among many cases studied by the author, the examples shown in
Figure 1 are considered. This model consists of four and three elastic nodes. In this example, the length of the bar was assumed
L = 1, and compression load
P = 1. Numerical examples are solved for various amplitudes of initial imperfections
(Equation (15)) introduced as a linear combination of eigenvectors
obtained from the linear eigenvalue problem (Equation (7)) scaled by buckling mode proportionality factor
. Keeping in mind that the eigenvector is normalised, thus the displacements are not real, the buckling mode proportionality factor should be assumed so as to reflect real imperfections. This can be achieved by the measurement of imperfection in situ, based on tolerances given in codes or assumed as arbitrary amplitudes.
In this example, equilibrium paths in the model considered are analysed, taking into account different configurations of geometrical imperfections. The increasing values of imperfections associated with three successive forms of the eigenvector are considered. The amplitudes of the mode proportionality factor
, used to create different configurations of the initial vector of geometrical imperfections are presented in
Table 1.
Moreover, based on the analysis of potential energy, the nature of equilibrium states is described. It should be noted that the structure is in equilibrium when potential energy reaches a minimum. The necessary condition for the occurrence of minimum potential energy is the disappearance of its first variation. This condition can be written as:
In the case of the example considered, potential energy can be described according to the formula:
where
λ is a load multiplier.
The nature of the equilibrium path is determined by higher order variations of potential energy. The type of bifurcation point can be determined based on the analysis of the second variation of potential energy as stable, neutral or unstable, according to:
Numerical examples are described by symbols composed of small letters of the Latin alphabet and numbers. The letter means the amplitude of the adopted buckling mode proportionality factor (m—small, s—medium, d—large), and the number—the form of a buckling mode, according to which the imperfection vector is determined, shown in
Figure 3. Thus, symbol m1 denotes an example in which imperfections are introduced according to the shape of the first buckling mode with a small amplitude
α1 = 0.001.
Figure 5 presents the equilibrium paths obtained for the structure with small amplitudes of initial geometrical imperfections developed in the shape of the first, second and third buckling modes, respectively (
m1,
m2 and
m3). It is worth noting that in all three cases, a stable equilibrium path is observed, where there are points of bifurcation and the accompanying unstable equilibrium path.
When imperfection related to the first buckling, the bifurcation point corresponds almost exactly to the first eigenvalue obtained for the ideal structure. However, when imperfections are developed in accordance with the second buckling mode, the bifurcation point corresponds to the second eigenvalue, etc. It can be seen that the shape of the initial geometrical imperfection strongly determines the nature of the equilibrium path.
Figure 6 shows the equilibrium paths obtained for examples in which geometrical imperfections are determined according to the small, medium and large amplitude of the first buckling mode (
m1,
s1 and
d1).
In all three cases, the bifurcation point associated only with the first critical load is observed. However, it can be observed that an increase in the amplitude of the initial geometrical imperfection results in an increase in the value of the critical load. Because the structure under consideration is primarily characterised by a stable equilibrium path, an increase in the imperfection amplitude somehow defies the snap-through of the structure to the configuration described by the unstable equilibrium path, and in this situation plays a positive role. In the case considered, the displacement both before and after the snap-through has the shape of the first buckling mode (
Figure 7). Moreover, based on the analysis of the second-order variation of potential energy, it is found that before snap-through the equilibrium path is stable, and afterwards it turns into an unstable one. This is illustrated in the figures by a ball moving on a differently shaped surface.
To determine the contribution of particular buckling mode in the displacement vector for each load increment, factor
is calculated according to the following relationship:
Figure 8 shows the contribution of an individual buckling mode in the displacement vector. Coefficients
,
and
illustrate the influence of the first, second and third buckling mode on the displacement vector, respectively. It is worth noting that in the example considered, in the whole range of the equilibrium path, the first buckling mode has a significant share in the displacement vector. Thus, the initial imperfection, determined by the first buckling mode, strongly influences the shape of the total displacement vector.
In the subsequent examples, initial geometric imperfections with three different amplitudes, developed according to the second or third buckling mode, are taken into account. For all imperfection configurations, only one bifurcation point is observed. However, the imperfection in the shape of the second buckling mode is associated with a second eigenvalue (
Figure 9). In turn, for imperfections developed in accordance with the third buckling mode, the bifurcation point is associated with the third eigenvalue (
Figure 10). As in the first case, an increase in the amplitude of the initial geometrical imperfection results in an increase in the value of the critical load. The displacements before and after a snap-through for all the initial amplitudes of the geometrical imperfections developed according to the second and third buckling forms have the shape of the second (
Figure 11) and the third buckling form, respectively (
Figure 12).
This phenomenon is confirmed by the analysis of coefficient
ρ, which obtains non-zero values before and after snap-through only for the second buckling mode (
Figure 13) for both small and large initial amplitudes of the geometrical imperfections developed in accordance with the second buckling mode.
Similarly, when imperfections in the form of the third buckling mode are introduced into the structure model, non-zero values of coefficient
ρ are obtained only for the third buckling mode (
Figure 14).
Next, analyses were carried out taking into account examples in which initial geometrical imperfections developed in the form of two different buckling modes.
Figure 15 shows equilibrium paths obtained for geometrical imperfections developed according to the first low-amplitude buckling mode and second or third medium amplitude buckling mode, respectively (
m1
s2 or
m1
s3). In turn,
Figure 16 shows the equilibrium paths for geometrical imperfections corresponding with a large amplitude of the first buckling mode and medium amplitude of the second or third buckling mode (
d1
s2 or
d1
s3). In all cases, two bifurcation points are observed: the first is associated with the first and the second with the second or third eigenvalue. It is observed that the introduction of geometrical imperfections developed according to two different buckling mode results in two or three snap-throughs associated with transition from the local to the global minimum of potential energy. The nature of a bifurcation point is defined by the second-order variation of potential energy and is observed to be stable before the first snap-through, and next, in case of a higher snap-through, unstable or metastable. A metastable bifurcation point is characterised by local minimum of potential energy and temporary equilibrium state. This phenomenon is illustrated by a ball moving on a differently shaped surface. The analysis of displacements before and after a snap-through shows that initial imperfection, which is a superposition of two different buckling modes, strongly influences the shape of the displacement vector.
However, it is still visible that at the first bifurcation point, the most influential is the first buckling mode, at the second bifurcation point—the second buckling mode, and at third bifurcation point—the third buckling mode (
Figure 17 and
Figure 18). Thus, as in the previous examples, the type of initial geometrical imperfection strongly affects the whole critical path, which is also confirmed by the analysis of the coefficient of the buckling mode contribution in displacement vector
ρ (
Figure 19 and
Figure 20). It can also be seen that in this case an increase in the amplitude of imperfection results in an increase in the value of the critical force.
Another group are examples in which initial geometrical imperfections are developed according to three different buckling modes. In these cases, the structural response is more complicated.
Figure 21 shows equilibrium paths obtained for selected configurations of initial geometrical imperfections corresponding to a small amplitude of the first buckling mode, medium or large amplitude of the second or third buckling mode
m1
d2
s3 lub
m1
s2
d3. In turn, in
Figure 22, the pre- and post-snap through displacements observed in the above-mentioned example were compared. Three bifurcation points are observed which, depending on the type of an initial geometrical imperfection pattern, differ both in the amplitude of the displacements and the value of the eigenvalue.
It can be observed that the introduction of geometrical imperfections developed according to the three different buckling modes results in several snap-throughs which are stable, unstable or metastable. In such cases, the nature of the equilibrium path in the post-buckling range is very complex and associated with multi-transition from the local to global minimum of potential energy. Several temporary states are observed in the form of metastable bifurcation points.
However, the snap-through is associated with the first, second, and third buckling mode, respectively. Increasing amplitudes of geometrical imperfections in the shape of the second or third form of buckling do not significantly affect the eigenvalue associated with the first buckling mode. In turn, a significant increase in the second eigenvalue is observed for a large amplitude of geometrical imperfection in the shape of the second buckling mode and, similarly, an increase in the third eigenvalue is observed for large amplitudes of imperfections developed according to the third buckling mode. The analysis of coefficient
ρ leads to the conclusion that we are dealing with the interaction of three buckling modes, because in the whole range of the considered equilibrium path, all these modes are involved in the displacement vector (
Figure 23). The interaction of the analysed phenomenon is also confirmed by the analysis of the deformation shape of the bar considered, before and after a snap-through. It is clearly visible that the initial geometrical imperfections, which in this case is a superposition of various amplitudes of three successive buckling modes, determine the shape of the bar displacement before and after a snap-through. There are no pure forms of deformation here that would perfectly suit subsequent buckling modes, but disrupted and asymmetrical ones. It is interesting that the maximum value of displacement before a snap-through is greater than afterwards, which indicates that the potential energy of the structure configuration changes from higher to lower.