{'fraction_non_alphanumeric': 0.025979721885400275, 'fraction_numerical': 0.01467970382601068, 'mean_word_length': 4.460205662769404, 'pattern_counts': {'":': 0, '<': 7, '<?xml version=': 0, '>': 7, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Parametric Analysis of the Shear Lag Effect in Tube Structural Systems of Tall Buildings

Horizontal loads such as earthquake and wind are considered dominant loads for the design of tall buildings. One of the most efficient structural systems in this regard is the tube structural system. Even though such systems have a high resistance when it comes to horizontal loads, the shear lag effect that is characterized by an incomplete and uneven activation of vertical elements may cause a series of problems such as the deformation of internal panels and secondary structural elements, which cumulatively grow with the height of the building. In this paper, the shear lag effect in a typical tube structure will be observed and analyzed on a series of different numerical models. A parametric analysis will be conducted with a great number of variations in the structural elements and building layout, for the purpose of giving recommendations for an optimal design of a tube structural system.


Introduction
Today, most of the mid to high tall buildings are built with tube structural systems. This type of system is known for its circumferential frame structure composed of closely placed columns and high beams, forming a tubelike box cross section, whose main purpose is to resist horizontal loads such as winds and earthquakes. On the inside of the building, additional columns may be found for bearing only vertical loads. Additionally, a combination of a tube structure with a central core is quite often, thus further enhancing the stiffness of the structure. The greatest advantages of the tube structural system are the ability to build very tall buildings, the enhanced stiffness of the building, the possibility of construction in various materials such as reinforced concrete and steel, possibility to use prefabricated elements for the purpose of repeating the floor dimensions and plans, and finally, a much-desired openness for the interior design of the building. Aside from very small openings that are a direct result of the column spacing on the exterior of the building, the most significant shortcoming of the tube system is the shear lag effect [1]. This effect is also very common in other structural systems of tall buildings [2][3][4] which will be discussed in detail in chapter 2 of this paper.

Shear Lag Effect in Tube Structural Systems
As mentioned previously, this type of structural system is derived from the rigid frame system, which implies that the exterior columns, which form the tube, are closely spaced and connected with high beams on every floor. So, it can be concluded that the tube system consists of multiple very rigid and connected frame panels which, while bending, act as webs and flanges of a tube cross section. During a horizontal load, due to their flexibility, beams in the web and flange frames deform in bending. This results in the increase of axial forces in the columns at the end of the frame panels and the decrease of axial forces in the columns towards the middle of the panel. Consequently, deformations occur in the 2 of 16 floor slabs and the interior secondary elements which leads to the Euler-Bernoulli beam theory bending of cross sections being no longer valid. Hence, the distribution of stress becomes nonlinear and very complex [5]. The nonlinearity is found in the flange frames (frame panels perpendicular to the direction of the horizontal load) as well as in the web frames (frame panels parallel to the direction of the horizontal load). This phenomenon can be seen in Figure 1 and is known as the positive shear lag [6].
Appl. Sci. 2021, 11, x FOR PEER REVIEW 2 of 17 of axial forces in the columns towards the middle of the panel. Consequently, deformations occur in the floor slabs and the interior secondary elements which leads to the Euler-Bernoulli beam theory bending of cross sections being no longer valid. Hence, the distribution of stress becomes nonlinear and very complex [5]. The nonlinearity is found in the flange frames (frame panels perpendicular to the direction of the horizontal load) as well as in the web frames (frame panels parallel to the direction of the horizontal load). This phenomenon can be seen in Figure 1 and is known as the positive shear lag [6]. On the other hand, the negative shear lag effect is characterized with increased axial forces in the columns in the middle of the panels ( Figure 2). To explain this phenomenon, a simple box section 5000 mm high cantilever model is used, loaded at mid-height by point load (Figure 3a) and uniform continuous load (Figure 3b). This model is introduced as an illustrative example, and as such it comprises a rectangular thin-walled steel section consisting of 1000 mm wide and 100 mm thick panels. These panels are modelled each with 10 shell elements along the panel width allowing an in-plane shear deformation which is responsible for shear lag effect. The stresses of the flange opposite to the load are shown in the graphic, ranging from zero (dark blue) to maximum compression (red). Isoareas of the same color show matching stress levels on the panel surface, thus allowing sections with positive or negative shear lag to be recognized. Since the value of the bending moment and the shear force at the free end of the cantilever equal zero, the flange stays undeformed in this position. On the fixed end of the cantilever, positive shear lag occurs due to the longitudinal shear in the flange. Due to the compatibility of deformations for all adjacent sections between the beginning and the end of the cantilever, stresses of the opposite distribution occur along the height as we move away from the fixed edge (negative shear lag). On the quarter of the cantilever length, closer to the fixed end, the axial stresses are almost equal along the whole flange width for the point load, while for continuous load some negative shear lag occurs at that section. So, in conclusion, the negative shear lag effect is a result of the compatibility of displacements (the cross section at the end of the cantilever must be straight and in-plane) [7]. On the other hand, the negative shear lag effect is characterized with increased axial forces in the columns in the middle of the panels ( Figure 2). To explain this phenomenon, a simple box section 5000 mm high cantilever model is used, loaded at mid-height by point load (Figure 3a) and uniform continuous load ( Figure 3b). This model is introduced as an illustrative example, and as such it comprises a rectangular thin-walled steel section consisting of 1000 mm wide and 100 mm thick panels. These panels are modelled each with 10 shell elements along the panel width allowing an in-plane shear deformation which is responsible for shear lag effect. The stresses of the flange opposite to the load are shown in the graphic, ranging from zero (dark blue) to maximum compression (red). Iso-areas of the same color show matching stress levels on the panel surface, thus allowing sections with positive or negative shear lag to be recognized. Since the value of the bending moment and the shear force at the free end of the cantilever equal zero, the flange stays undeformed in this position. On the fixed end of the cantilever, positive shear lag occurs due to the longitudinal shear in the flange. Due to the compatibility of deformations for all adjacent sections between the beginning and the end of the cantilever, stresses of the opposite distribution occur along the height as we move away from the fixed edge (negative shear lag). On the quarter of the cantilever length, closer to the fixed end, the axial stresses are almost equal along the whole flange width for the point load, while for continuous load some negative shear lag occurs at that section. So, in conclusion, the negative shear lag effect is a result of the compatibility of displacements (the cross section at the end of the cantilever must be straight and in-plane) [7].    The first researches of the shear lag and the effect it has on different structural systems date back to the 1980s [8]. In the very beginning this phenomenon was examined in steel and concrete girders, mostly bridges. Researches involving steel box girders and bridges were the most prominent ones [8][9][10][11]. In [9,10] theoretical foundations for the shear lag effect are given with a special emphasis on steel box girders. At about the same time [8,10,11] basic mathematical models and formulas for the calculation of the shear lag effect in box girders were given. While in [8,11] the focus is set on the positive shear lag effect and the shear lag effect in general, in [10] the focus is set on the negative shear lag effect and the calculation methods that would best approximate it. However, this effect is quite a problem in tube structural systems of tall buildings as well. Although this effect was mostly examined in detail in box girders and bridges, the first mathematical expressions for tall buildings were set in the 1990s [12]. In [12] a mathematical model for the preliminary design of the shear lag effect in tall buildings is given. In the proposed method, independent distributions of axial displacements are used for the web and flange panels and thus the shear lag in each panel is individually allowed for. Lastly, it is deduced that the shear lag effect is dependent on the type of the horizontal load and that the positive shear lag effect is greater at the base of the building than the negative shear lag is at its top. A great advancement in theoretical knowledge is seen mostly in the complexity and preciseness of analytical methods like the one shown in [13]. In this research the shear lag effect in tall tube structures is shown regarding the axial stresses and displacements of columns. At the same time, the first computer-based models for this phenomenon are developed. According to [14] a simple model is developed for the purpose of the calculation and design of a tube structure under horizontal loads. Additionally, an expression to assess the horizontal displacement of a tall building is given with the shear lag effect considered while developing the said expression. In [5] a much more complex computerbased program is developed using a 40 storey tall tube structure. It is also confirmed that The first researches of the shear lag and the effect it has on different structural systems date back to the 1980s [8]. In the very beginning this phenomenon was examined in steel and concrete girders, mostly bridges. Researches involving steel box girders and bridges were the most prominent ones [8][9][10][11]. In [9,10] theoretical foundations for the shear lag effect are given with a special emphasis on steel box girders. At about the same time [8,10,11] basic mathematical models and formulas for the calculation of the shear lag effect in box girders were given. While in [8,11] the focus is set on the positive shear lag effect and the shear lag effect in general, in [10] the focus is set on the negative shear lag effect and the calculation methods that would best approximate it. However, this effect is quite a problem in tube structural systems of tall buildings as well. Although this effect was mostly examined in detail in box girders and bridges, the first mathematical expressions for tall buildings were set in the 1990s [12]. In [12] a mathematical model for the preliminary design of the shear lag effect in tall buildings is given. In the proposed method, independent distributions of axial displacements are used for the web and flange panels and thus the shear lag in each panel is individually allowed for. Lastly, it is deduced that the shear lag effect is dependent on the type of the horizontal load and that the positive shear lag effect is greater at the base of the building than the negative shear lag is at its top. A great advancement in theoretical knowledge is seen mostly in the complexity and preciseness of analytical methods like the one shown in [13]. In this research the shear lag effect in tall tube structures is shown regarding the axial stresses and displacements of columns. At the same time, the first computer-based models for this phenomenon are developed. According to [14] a simple model is developed for the purpose of the calculation and design of a tube structure under horizontal loads. Additionally, an expression to assess the horizontal displacement of a tall building is given with the shear lag effect considered while developing the said expression. In [5] a much more complex computer-based program is Appl. Sci. 2021, 11, 278 5 of 16 developed using a 40 storey tall tube structure. It is also confirmed that the positive shear lag effect is related to the lower storeys and that the negative shear lag is related to the higher storeys of the building. Finally, it is concluded that the positive shear lag is the source of the negative shear lag. At the beginning of the 2000s the analysis of this phenomenon is expanded to more complex structural systems like tube in tube structures [15][16][17] where it is deduced that the negative shear lag effect is more prominent at the lower part of the interior tube than at the lower parts of the exterior tube. Furthermore in [17] exact analytical expressions are given for the calculation of the shear lag effect in this type of structural system. On the other hand, in diagrid structural systems the shear lag effect is mainly influenced by two parameters, the angle under which the exterior diagonal columns cross each other and the way the shear lag effect is taken into account in the preliminary design of structural elements [2,18]. The next step in the research of the shear lag effect was the pursuit of ideas and ways to minimize its impact. In [6] the optimal angle under which the exterior bracings cross to minimize the shear lag effect is given. In [19] a similar analysis is conducted with the tube structural system and diagonal bracing where the shear lag effect in an ordinary tube structure is compared to its effect in a tube structure with diagonal bracing. It is deduced that the implementation of diagonal bracing has a great impact on the shear lag effect and that the angle at which the diagonals are placed is the most important factor to consider. Furthermore, according to [20] the use of two diagonals in an "X" like fashion is much more efficient than the use of one diagonal in any configuration regarding the decrease of the shear lag effect. The same principle of bracing can be used for tube in tube structural systems [21]. In recent years, the research conducted is mostly based on parametric analysis of different structural systems [22][23][24][25]. In [22] a parametric analysis is conducted where the main parameter being varied are the layout dimensions of the building. The main idea is to see how the ratio of the layout dimensions of a rectangular shaped building affects the buildings behavior when subjected to strong winds. Even though the shear lag effect is not addressed in this paper per se, an insight is given into the optimal ratio of layout dimensions for a tall tube structure. In [23] besides the ratio of layout dimensions of a tube structure, the effect of column spacing is also taken into account. In their research it is concluded that the layout dimensions of a tube structure have no or minimal impact on the shear lag effect. Furthermore, it is stated that the column spacing has a greater impact in the negative shear lag area of the building than in the positive shear lag area of the building. In the end, the impact that diagonal bracing placement on the exterior of the structure has on the shear lag effect is examined. A similar analysis is conducted in [24] where the use of mega bracing in tall buildings is studied with the conclusion that the X pattern diagonal bracing shows the most promising results regarding the reduction of the shear lag effect.
In this paper the impact of the negative and positive shear lag effect on a tall tube structure will be analyzed in detail. Using a computer-based model a parametric analysis will be conducted with a series of different parameters being varied with the main purpose of detecting the critical height of the tall building where the transfer from positive to negative shear lag occurs. Besides that, the impact of every parameter variation on the shear lag effect will be shown. The parameters that will be examined include the ratio of column spacing and column dimensions, the ratio of storey height and the height of connecting beams, the ratio of the layout dimensions of the entire building, and the type (profile) of horizontal load being used. In the end, the main goal is to find the optimal layout dimensions of a tall building with a tube structural system and the optimal dimensions of the structural elements that are used in it.

Overview of the Existing Tall Buildings with Tube Structural System
To set the boundaries of parameters for the computer-based model, an overview of the most significant buildings with the tube structural system is given. The parameters that were included are: material used, storey height, building height, layout dimensions, column spacing and column dimensions (Table 1). According to most used parameter values range, a computer-based template model is developed.  The columns and beams are modelled with one dimensional finite elements (6 degrees of freedom), and assigned a corresponding cross section properties (moment of inertia for both axis, sectional area, shear deformation area, torsional moment of inertia, warping modulus and concrete material linear properties) according to column or beam member dimensions. They are connected with a fixed connection at the intersection of their axes, allowing the transfer of moment, shear force and axial force. Floor slabs are modeled as rigid diaphragms ( Figure 5), by coupling of the nodes in the same storey to joint in-plane movements and rotations (flexible slab-rigid disk behavior), as allowed and recommended by Eurocode design code A linear statical analysis is performed.

Computer-Based Model of the Tube Structure
The building is loaded with the dead load (no additional vertical loads were considered besides the dead load) of all structural elements and wind load uniformly distributed along the height. The value of the uniformly distributed load is determined as the maximum expected value based on EN 1991-4 [26] while considering the top of the structure as a reference height. In addition, an actual wind profile will be used in the parametric analysis that is also calculated by EN 1991-1-4 [26]. The horizontal loads are concentrated at each intersection of beams and columns ( Figure 6).       The building is loaded with the dead load (no additional vertical loads were considered besides the dead load) of all structural elements and wind load uniformly distributed along the height. The value of the uniformly distributed load is determined as the maximum expected value based on EN 1991-4 [26] while considering the top of the structure as a reference height. In addition, an actual wind profile will be used in the parametric analysis that is also calculated by EN 1991-1-4 [26]. The horizontal loads are concentrated at each intersection of beams and columns ( Figure 6). The shear lag effect is observed throughout the building's height. In Figure 7 the shear lag effect is shown in seven characteristic storeys ranging from the 2nd to the 36th storey. On the horizontal axis from left to right eleven columns are placed and on the vertical axis the coefficient α is given. This coefficient shows the ratio of the compressive force in a given column to the compressive force in the middle column. Conclusively, the value α = 1 applies to the column in the very middle, that is to column number 6 ( Figure  7). By the values shown in Figure 7 it is clear that in the lower storeys (2nd and 6th) the positive shear lag is prominent and by moving towards the building's top its impact diminishes with the negative shear lag rapidly increasing. At the 12th storey it is noticeable that the shear lag effect is practically non-existent so at this position it is considered that the positive shear lag morphs into the negative shear lag. The shear lag effect is observed throughout the building's height. In Figure 7 the shear lag effect is shown in seven characteristic storeys ranging from the 2nd to the 36th storey. On the horizontal axis from left to right eleven columns are placed and on the vertical axis the coefficient α is given. This coefficient shows the ratio of the compressive force in a given column to the compressive force in the middle column. Conclusively, the value α = 1 applies to the column in the very middle, that is to column number 6 ( Figure 7). By the values shown in Figure 7 it is clear that in the lower storeys (2nd and 6th) the positive shear lag is prominent and by moving towards the building's top its impact diminishes with the negative shear lag rapidly increasing. At the 12th storey it is noticeable that the shear lag effect is practically non-existent so at this position it is considered that the positive shear lag morphs into the negative shear lag. Appl. Sci. 2021, 11, x FOR PEER REVIEW 9 of 17

Parametric Analysis
The parametric analysis will be conducted with four main parameters ( Figure 4): (1) s/c-the ratio of column spacing perpendicular to the horizontal load direction and column dimension, ranging from 4 to 6, in every next model variation one column is removed to maintain a constant building layout) (2) H1/h-the ratio of storey height and beam height (ranging from 2.  (Figure 8, in red) that is located on the opposite side of the horizontal load. In Figure 8 the case of the positive shear lag is shown. The number of columns in this case is N = 5.
Coefficient β is defined as follows: where: n-actual number of columns in half of the flange panel (in the case of uneven number of columns, a half of the column is used) nβ-number of columns involved in the transfer of loads calculated in Equation (2)

Parametric Analysis
The parametric analysis will be conducted with four main parameters ( Figure 4): (1) s/c-the ratio of column spacing perpendicular to the horizontal load direction and column dimension, ranging from 4 to 6, in every next model variation one column is removed to maintain a constant building layout) Additionally, it is then clear that the location with the value of coefficient β = 1 is the location in which no shear lag effect exists (the positive shear lag ends and the negative one begins). Accordingly, β < 1 applies to positive shear lag and β > 1 applies to negative shear lag. Finally, as the value of β moves away from the value equal to 1, either less or more, the shear lag positive or negative, increases. For the positive shear lag, it means that as the value of β decreases (moves away from the value equal to 1), the positive shear lag increases. For the negative shear lag, it means that as the value of β increases (moves away from the value equal to 1) the negative shear lag increases.

The Impact of Column Spacing on the Shear Lag Effect
In this chapter, the impact of exterior column spacing in a tube structure on the shear lag effect will be analyzed. As already stated, the column dimension perpendicular to the horizontal load direction will stay the same throughout the analysis (c = 0.5 m). The column spacing will be varied five times with values of s = 2.0 m, s = 2.18 m, s = 2.4 m, s = 2.67 m and s = 3.0 m while the layout dimensions stay constant A/B = 24 m/24 m (number of columns is varied). The entire building height and the storey height are constant. The same seven storeys as in chapter 4 will be analyzed. It is very important to note that as the column spacing is changed, the referent node area for point load changes. That means that for different column spacing the horizontal force in the nodes where beams and columns meet is different. The results of the analysis are shown in Figure 9. Coefficient β is defined as follows: where: n-actual number of columns in half of the flange panel (in the case of uneven number of columns, a half of the column is used) n β -number of columns involved in the transfer of loads calculated in Equation (2) for half of the flange panel where: ΣF i -sum of axial compressive forces in columns (F 1 to F 5 ) F 1 -compressive force in the edge (first or last) column of the flange panel Therefore, if the axial force on the edge of the flange panel (F 1 ) is the one with maximal value, then the flange panel is under positive shear lag. On the other hand, if the same force is given with its minimal value, then the flange panel is under negative shear lag. Additionally, it is then clear that the location with the value of coefficient β = 1 is the location in which no shear lag effect exists (the positive shear lag ends and the negative one begins). Accordingly, β < 1 applies to positive shear lag and β > 1 applies to negative shear lag. Finally, as the value of β moves away from the value equal to 1, either less or more, the shear lag positive or negative, increases. For the positive shear lag, it means that as the value of β decreases (moves away from the value equal to 1), the positive shear lag increases. For the negative shear lag, it means that as the value of β increases (moves away from the value equal to 1) the negative shear lag increases.

The Impact of Column Spacing on the Shear Lag Effect
In this chapter, the impact of exterior column spacing in a tube structure on the shear lag effect will be analyzed. As already stated, the column dimension perpendicular to the horizontal load direction will stay the same throughout the analysis (c = 0. According to the results ( Figure 9) the following can be concluded: (1) All storeys bellow the 12th storey are in the positive shear lag effect area (β < 1). It is also clear that the column spacing for these storeys has no meaningful impact on the positive shear lag effect (the value of β for the 2nd, 6th and 12th storey is not changing). (2) The transition from the positive to the negative shear lag effect area for all column spacing variants is at the 12th storey (β = 1), which is approximately at 33% of the building height measured from the ground. (3) All storeys above the 12th storey are in the negative shear lag effect area (β > 1). It is also clear that the column spacing for these storeys has some impact on the negative shear lag effect. As the building rises, so does the coefficient β. Furthermore, the increase of coefficient β in one storey due to the variation of column spacing is more and more prominent as the building rises. While on the 18th storey this increase is mild, on the 36th storey the increase of coefficient β with the increase of column spacing is quite notable.
It can be deduced that the column spacing has practically no impact on the positive shear lag effect, while in the case of negative shear lag effect, the increase in column spacing results in the increase of this effect with this phenomenon being more prominent as the building rises.

The Impact of Beam Height on the Shear Lag Effect
In this chapter, the impact of the beam height in a tube structure on the shear lag effect will be analyzed. While the beam height will be varied, the storey height will stay the same with the value of H1 = 3.   According to the results ( Figure 9) the following can be concluded: (1) All storeys bellow the 12th storey are in the positive shear lag effect area (β < 1). It is also clear that the column spacing for these storeys has no meaningful impact on the positive shear lag effect (the value of β for the 2nd, 6th and 12th storey is not changing). (2) The transition from the positive to the negative shear lag effect area for all column spacing variants is at the 12th storey (β = 1), which is approximately at 33% of the building height measured from the ground. (3) All storeys above the 12th storey are in the negative shear lag effect area (β > 1).
It is also clear that the column spacing for these storeys has some impact on the negative shear lag effect. As the building rises, so does the coefficient β. Furthermore, the increase of coefficient β in one storey due to the variation of column spacing is more and more prominent as the building rises. While on the 18th storey this increase is mild, on the 36th storey the increase of coefficient β with the increase of column spacing is quite notable.
It can be deduced that the column spacing has practically no impact on the positive shear lag effect, while in the case of negative shear lag effect, the increase in column spacing results in the increase of this effect with this phenomenon being more prominent as the building rises.

The Impact of Beam Height on the Shear Lag Effect
In this chapter, the impact of the beam height in a tube structure on the shear lag effect will be analyzed. While the beam height will be varied, the storey height will stay the same with the value of H 1 = 3. According to the results ( Figure 10) the following can be concluded: (1) All storeys bellow the 12th storey are in the positive shear lag effect area (β < 1). For these storeys, by decreasing the beam height, the positive shear lag effect increases. Examining the curves in Figure 10 it can be concluded that the positive shear lag effect is more prominent for lower storeys of the building and that coefficient β is smallest for the ground storey. (2) The transition from the positive to the negative shear lag effect area for all beam height variants is at the 12th storey (β = 1), which is approximately at 33% of the building's height measured from the ground. (3) All storeys above the 12th storey can be found in the negative shear lag effect area (β > 1). As the building rises, so does the coefficient β. Furthermore, the increase of coefficient β in one storey due to the variation of beam height is more and more prominent as the building rises. While on the 18th storey this increase is mild, on the 36th storey the increase of coefficient β with the decrease of beam height is quite notable.
It can be deduced that the decrease in beam height results in the increase of positive and negative shear lag effect in tall tube buildings.

The Impact of Layout Dimensions on the Shear Lag Effect
In this chapter, the impact of building layout dimensions on the shear lag effect will be analyzed. The parameters that are going to stay constant include the beams cross section dimensions b/h = 0.5 m/1.0 m, column cross section dimension c = 0.5 m, column spacing s = 2.4 m, storey height H1 = 3.5 m and the entire buildings height H = 126 m. The layout ratio of the building will be changed from a rectangular shape to a square shape starting from A/B = 24/43.2 = 0.56 and ending with the final ratio of A/B = 24/24 = 1.0. The value of A is the side of the building parallel to the horizontal load direction (web panel) and the value of B is the side of the building perpendicular to the horizontal load direction (flange panel). The results of the analysis are shown in Figure 11.   According to the results ( Figure 10) the following can be concluded: (1) All storeys bellow the 12th storey are in the positive shear lag effect area (β < 1). For these storeys, by decreasing the beam height, the positive shear lag effect increases. Examining the curves in Figure 10 it can be concluded that the positive shear lag effect is more prominent for lower storeys of the building and that coefficient β is smallest for the ground storey. (2) The transition from the positive to the negative shear lag effect area for all beam height variants is at the 12th storey (β = 1), which is approximately at 33% of the building's height measured from the ground. (3) All storeys above the 12th storey can be found in the negative shear lag effect area (β > 1). As the building rises, so does the coefficient β. Furthermore, the increase of coefficient β in one storey due to the variation of beam height is more and more prominent as the building rises. While on the 18th storey this increase is mild, on the 36th storey the increase of coefficient β with the decrease of beam height is quite notable.
It can be deduced that the decrease in beam height results in the increase of positive and negative shear lag effect in tall tube buildings.

The Impact of Layout Dimensions on the Shear Lag Effect
In this chapter, the impact of building layout dimensions on the shear lag effect will be analyzed. The parameters that are going to stay constant include the beams cross section dimensions b/h = 0. It can be deduced that the positive and negative shear lag effect are increased with the increase of the building side on which the horizontal loads acts.

The Impact of Wind Load Type on the Shear Lag Effect
In this chapter, the impact of the type of the horizontal load profile on the shear lag effect will be analyzed. The constant parameters include the beams cross section dimensions b/h = 0.5 m/1.0 m, column cross section dimensions c = 0.5 m, column spacing s = 2.4 m, the entire buildings height H = 126 m, storey height H1 = 3.5 m and building layout dimensions A/B = 24 m/43.2 m. The building layout dimensions were chosen differently from the previous analyses, so that the difference in effect of different horizontal load profiles is most apparent. Two horizontal load profiles shall be considered: the continuous wind load with a constant value throughout the building's height and the actual wind load profile calculated using EN 1991-1-4 [26]. The constant load profile includes the use of the maximum wind load force that is calculated for the top of the building after which that load is used throughout the building height (Figure 12b)  It can be deduced that the positive and negative shear lag effect are increased with the increase of the building side on which the horizontal loads acts.

The Impact of Wind Load Type on the Shear Lag Effect
In this chapter, the impact of the type of the horizontal load profile on the shear lag effect will be analyzed. The constant parameters include the beams cross section dimensions b/h = 0.5 m/1.0 m, column cross section dimensions c = 0.5 m, column spacing s = 2.4 m, the entire buildings height H = 126 m, storey height H 1 = 3.5 m and building layout dimensions A/B = 24 m/43.2 m. The building layout dimensions were chosen differently from the previous analyses, so that the difference in effect of different horizontal load profiles is most apparent. Two horizontal load profiles shall be considered: the continuous wind load with a constant value throughout the building's height and the actual wind load profile calculated using EN 1991-1-4 [26]. The constant load profile includes the use of the maximum wind load force that is calculated for the top of the building after which that load is used throughout the building height (Figure 12b). The actual wind load profile ( Figure 12a) is calculated using EN 1991-1-4 [26]. As it was the case in prior analysis, the same seven storeys will be considered. The results of the analysis are shown in Figure 13.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 14 of 17 ( Figure 12a) is calculated using EN 1991-1-4 [26]. As it was the case in prior analysis, the same seven storeys will be considered. The results of the analysis are shown in Figure 13.
(a) (b)     ( Figure 12a) is calculated using EN 1991-1-4 [26]. As it was the case in prior analysis, the same seven storeys will be considered. The results of the analysis are shown in Figure 13.
(a) (b)    According to the results the following can be concluded: (1) The 2nd, 6th and 12th storey are found in the positive shear lag effect area (β < 1). In this area the shear lag effect is more significant (smaller β) for the constant value (C) of the wind load which is most clearly seen on the 2nd storey of the building. As the building rises the gap between two wind load profiles regarding the positive shear lag effect diminishes with the positive shear effect being more prominent for the actual wind load profile on the 12th storey. (2) All storeys above the 12th storey are found in the negative shear lag effect area (β > 1).
In this area the shear lag effect is more significant (larger β) for the constant value (C) of the wind load than for the actual wind load profile. As the building rises, the gap between two wind load profiles regarding the negative shear lag effect increases until the 30th storey after which the gap slightly diminishes to the buildings top.
It can be deduced that the positive and negative shear lag effect are more significant for the continuous wind load with a constant value in most storeys of the structure.

Conclusions
The prototype structure of a 126 m high (36 storeys) tube structural system tall building was analyzed, comprising of a reinforced concrete column and beam panel members. For a varying set of parameters defining layout dimensions from 24 to 43 m, column spacing from 2 to 3 m, beam height from 0.8 to 1.2 m, and two different horizontal load profiles, following conclusions can be made regarding changes in the shear lag effect: (1) The positive shear lag effect is characteristic for the lower storeys of the structure and as the structure rises, the positive shear lag decreases and gradually transforms into negative shear lag which then increases to the buildings top. This transformation usually happens at around 33% of the building height measured from the bottom to the top except in the case of the variation of the layout dimensions where the transformation happens at 50% of the building height. (2) The increase of column spacing, with the layout dimensions of the building staying the same, has no major impact on the positive shear lag effect that is characteristic for the lower storeys of the building. The position of the change from the positive to the negative shear lag is also not conditioned by the column spacing. In upper storeys, where negative shear lag appears, with the increase in column spacing the negative shear lag increases. Please turn to the CRediT taxonomy for the term explanation. Authorship must be limited to those who have contributed substantially to the work reported. All authors have read and agreed to the published version of the manuscript.
Funding: This research received no external funding.
Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy restrictions.