“…Since previous studies [31,32,10,11,26,27] have clearly demonstrated that it is essential to include the shear strain contributions into the total potential energy formulation to model the interactive buckling behaviour, Timoshenko beam theory is assumed currently. The global mode is decomposed into two components: a pure lateral displacement W and a pure rotation of the plane sections θ, see Figure 2, known as the 'sway' and 'tilt' modes [33,31] respectively.…”