Numerical Simulation of Wind-Induced Vibration Response Characteristics of High-Rise Buildings with Podiums

1. Introduction

As the name suggests, wind-induced vibration is the term used to describe the shaking of structures. The effect of wind on building vibration was not immediately apparent in an era when buildings were typically not very tall. It was only later, with the development of high-rise structures and wide-span bridges, etc., that wind-induced vibration events started to become more common and even resulted in wind-damaged buildings, like the Tacoma Suspension Bridge in the United States in 1940[1], that people started to investigate this effect. Numerous academics have investigated wind vibration in buildings. To study the wind loads of supertall multi-tower conjoined buildings from the perspectives of mean and pulsating wind pressure characteristics and vortex volume distribution, researchers Shitang Ke and Hao Wang [2,3] performed wind tunnel tests based on large vortex simulation for a supertall three-tower conjoined building. In order to study the spatial distribution characteristics of downwind pulsating wind loads and pulsating wind speed coherence functions, Zeng Jadong and Li Mingshui[4] et al. conducted surface pressure simultaneous measurement experiments on a rectangular rigid model. They then proposed a model expression for the downwind pulsating wind load coherence function of a rectangular building that can involve flow field and structural parameters like wind field category, vertical spacing, and others. In order to study the interrelationships reflecting the statistical coupling between modal responses under spatial-temporal varying dynamic wind excitation, M. F. Huang and C. M. Chan[5] et al. use a 60-story asymmetric full-size steel-concrete hybrid building as an example. They then propose a framework for analyzing lateral torsional wind vibration forces in high-rise buildings based on three-dimensional modal vibration patterns. In order to study wind-induced torsional loads on rectangular high-rise buildings, Vladimir Guzmán-Sols, Adrian Pozos-Estrada[6], et al. conducted wind tunnel tests on five rigid models with the same dimensions and various aspect ratios. Based on the results, they proposed a parameter on the variation of foundation torque coefficients with aspect ratio and wind direction for rectangular high-rise buildings that agrees well with the experimental results equations. A popular area of study in recent years has been the damage theory analysis of high-rise building structures. The numerical simulation study of these structures’ wind vibration response characteristics has some value for engineering.

A survey of existing domestic and foreign research shows that few studies of this kind have been carried out, despite the fact that the main podium one-piece high-rise building structure is a common structure in a class of projects. In this study, the pulsating wind velocity at various heights is solved using the AR approach, and the pulsating wind load is calculated using a formula. Finally, the software simulates the building’s response to wind vibration under pulsing wind stress[7].

2. Overview of AR Method

The pulsating wind speed time range must be resolved first in order to get the building pulsating wind load time range by self-programming. The wind tunnel test, harmonic synthesis method, and linear filtering method are the current three methods for determining the pulsing wind speed. A series of artificially generated random numbers with zero mean and a white spectrum are input into the filter for processing and are then output as a random process with a specific spectrum. This technique is known as the linear filtering method, also known as the white noise filtering method. The linear filter method can be further broken down into the Autoregressive (AR), Sliding Average, etc. Due to its benefits of low computational effort, quick computation, and high accuracy in determining the pulsing wind speed of buildings, the AR method has been extensively employed by a significant number of academics[8,9,10].

The basic idea behind the AR method’s steps is to first locate the stochastic process, then determine the regression coefficients given the stochastic process. The M-dimensional AR model can be used to represent the pulsing wind velocity field.

$$\mathit{v}\left(t\right)=-\sum _{k=1}^q{\mathit{\psi }}_k\mathit{v}\left(t-k∆t\right)+\mathit{N}\left(t\right)$$

(1)

where: q is the order of the M-dimensional AR model [11], ${\mathit{\psi }}_k$ is the regression coefficient matrix of order M × M, $k$ is the surface resistance coefficient, $∆t$ is the time step, and $\mathit{N}\left(t\right)$ is a random vector with zero mean for a given variance. The matrix of autoregressive coefficients can be written as follows:

$$\mathit{R}=-{\mathit{R}}^{\ast }{\mathit{\psi }}_k$$

(2)

The expression is where $-{\mathit{R}}^{\ast }$ is the autocorrelation matrix of order I × I (I=q×M), and The elements in $-{\mathit{R}}^{\ast }$ can be obtained using the Wiener Khinchin formula:

$${\mathit{R}}_{ij}={\int }_0^{\infty }{\mathit{S}}_{ij}\cos \left(2\pi nj∆t\right)dn,i,j=\mathrm{1,2}\dots ,M$$

(3)

where ${\mathit{S}}_{ij}$ is the pulsating wind speed mutual spectrum density function, $n$ stands for the pulsating wind speed frequency, ${\mathit{S}}_{ij}$ can have a known pulsating wind speed spectrum $\mathit{S}$, and the spatial correlation coefficient $\mathrm{c}\mathrm{o}\mathrm{h}(r,n)$ is used to find. The following section will introduce a number of readily available pulsing wind speed spectra. The expression for the spatial correlation coefficient $\mathrm{c}\mathrm{o}\mathrm{h}(r,n)$ is:

$$\mathrm{c}\mathrm{o}\mathrm{h}\left(r,n\right)=\exp \left\{{\displaystyle \frac{-2n\sqrt{C_x^2{\left(x_i-x_j\right)}^2+C_y^2{\left(y_i-y_j\right)}^2+C_z^2{\left(z_i-z_j\right)}^2}}{v_{z_i}+v_{z_j}}}\right\}$$

(4)

where: $v_{z_i}$ and $v_{z_j}$ signify the average wind speed at $z_i$ and $z_j$ on the space, and $C_x$, $C_y$, $C_z$ denote the left-right, front-back, and top-bottom attenuation coefficients of a space, respectively.

By multiplying both sides of equation (1) simultaneously right by v(t)=[v_1 (t),v_2 (t),…,v_m (t)], and then substituting the derived autoregressive coefficient matrix ψ_k, we can obtain the expression:

$${\mathit{R}}_N=\sum _{k=1}^q{\mathit{\psi }}_k\mathit{R}\left(k∆t\right)+\mathit{R}\left(0\right)$$

(5)

${\mathit{R}}_N$ is subject to Cholesky decomposition, and the decomposition expression is:

$${\mathit{R}}_N=\mathit{L}{\mathit{L}}^T$$

(6)

Thus the random vector N(t) with a zero mean given the variance may be calculated from L, which stands for the lower triangular matrix. The expression is:

$$\mathit{N}\left(t\right)=\mathit{L}\mathit{N}\left(t\right)$$

(7)

By inserting the resulting $\mathit{N}\left(t\right)$ into equation (1), it is possible to find the pulsing wind velocity process vector with time interval t and spatial correlation.

The expression employed in the study, which uses the Davenport[12] pulsating wind speed spectrum, is:

$$S_v\left(n\right)={\displaystyle \frac{4k{x}^2{v}_{10}^2}{n{\left(1+x^2\right)}^{\frac{4}{3}}}}$$

(8)

$$x={\displaystyle \frac{1200n}{v_{10}}}$$

(9)

where $v_{10}$ defines the reference wind speed at 10 m height, $k$ denotes the ground roughness coefficient, $S_v\left(n\right)$ denotes the pulsing wind speed power spectrum, $n$ denotes the pulsating wind frequency.

3. Pulsating Wind Load Solution

The program is created using the AR method, imported into the Matlab software, and run to output the pulsating wind speed time history data and the pulsating wind speed time history graph at each point. Finally, the calculated simulated spectrum curve is compared with the target spectrum curve to confirm the accuracy of the simulated wind speed time history data.

3.1. Computational Model

The study conditions are shown in Table 1. The main building is 40 meters in length, width, and height. The model size is shown in Figure 1. Measurement point locations are 10, 28, 56, 84, 112, 140, 168, and 196 meters on the center line of the main building’s windward side. The measurement point map is shown in Figure 2. (podium area working condition measurement point distribution and 24m podium working condition is the same).

3.2. Time History Simulation of Fluctuating Wind Speed

Wang Xiuqiong[13] conducted a pertinent study regarding the Davenport wind speed spectrum and provided more reasonable values for the ground roughness coefficient k according to different roughness types, including 0.00129 for Class A landforms, 0.00215 for Class B landforms, 0.00464 for Class C landforms, and 0.01291 for Class D landforms(in the Code for Structural Loads of Buildings (GB50009-2012)[14] the ground roughness is divided into four categories A, B, C and D). Since B is the studied landform type, 0.00215 is used as the k value. Only height correlation is taken into account for spatial correlation, hence $C_z=10$ is the only number used (empirical value). Equation (4) uses the wind speed exponential law to calculate the variables $v_{z_i}$ and $v_{z_j}$. The wind speed exponential law’s expression is:

$$V_z=V_0{\left({\displaystyle \frac{Z}{Z_0}}\right)}^{\alpha }$$

(10)

where: $Z_0$ is the standard height; $V_0$ is the standard wind speed at $Z_0$; $V_z$ is the average wind speed at $Z$ above the ground; $\alpha$ is the ground roughness index.

The reference wind pressure in this study is 0.25 at a height of 10 meters, based on the equation between wind pressure and wind speed:

$$p={\displaystyle \frac{v^2}{1600}}$$

(11)

20 m/s can be calculated as the reference wind speed at 10 meters, and the other necessary values are established as given in Table 2.

According to the calculation of MATLAB program, the fluctuating wind speed at each height is shown in Figure 3.

According to Figure 3, each height measurement point’s peak pulsing wind speed is maintained at roughly 10 m/s, and they all oscillate up and down at 0. The simulated spectrum curve (black straight line on the right Figure) produced by Matlab software for each height measurement point is more consistent with the target spectrum curve (red dashed line on the right Figure), as can be seen from the power spectrum comparison graph, demonstrating that there is some degree of confidence in the simulated pulsating wind speed time range data in this section.

The pulsating wind pressure time range data can be calculated using the obtained pulsating wind speed time range data, and the pulsating wind load time range data can be calculated using a different formula. The wind pressure includes both average and pulsing wind pressure since the natural wind in the atmospheric environment is made up of both average and pulsating wind. The equation for the link between wind speed and wind pressure is:

$$P={\displaystyle \frac{1}{2}}\rho V^2$$

(12)

where $\rho$ is the air density, and $V$ is the total wind speed. Since the total wind speed is composed of average wind speed and fluctuating wind speed, equation (12) can be expressed as:

$$P={\displaystyle \frac{1}{2}}\rho {\left(\overline{v}+v\right)}^2$$

(13)

where $\overline{v}$ is the average wind speed, $v$ is the fluctuating wind speed. The wind speed exponential rule is used to get the average wind speed at each location. Table 3 displays the average wind speed at each height measurement point, and Figure 4 displays the pulsing wind pressure time range at each height measurement point that was determined by Equation (13).

Based on the calculated fluctuating wind pressure time history data, the fluctuating wind load time history can be directly calculated through the formula, which is:

$$F=PA{\mu }_s$$

(14)

In this study, the windward area A of the simulated height measurement point is calculated as the product of the measurement point’s distance from two neighboring measurement sites divided by half, then multiplied by the breadth of the windward surface. ${\mu }_s$ indicates the building wind carrier type coefficient, which is specified in the Code for Structural Loads of Buildings(GB50009-2012)[14].

To be conservative, the wind load coefficient on the windward side and the wind load coefficient on the leeward side are typically superimposed to take 1.4 when designing and calculating the wind resistance performance of the building. This is because in the wind resistance design of the building structure, the wind pressure on the windward side of the building and the wind suction force on the leeward side of the building should be in the same direction. The wind load coefficient on the windward side is taken to be 1.25 in this study, and Figure 5 displays the wind load time range for each height measurement point derived from Equation (14).

4. Wind Induced Vibration Analysis

4.1. Comparison of Wind Vibration Results Between High Rise Buildings with and Without Podiums

Prior to performing the time analysis, the modal analysis of the finite element model’s first 16 orders of vibration patterns is performed using the eigenvalue analysis method in Midas Gen, which may guarantee the model’s plausibility. Under the condition of calculating the vibration mode to the 16^th order, the Table4 shows the first six orders of the vibration mode diagram for the high building without podium.

According to Table 4, the fact that each order’s vibration type essentially complies with the requirements listed for the first six high-rise building vibration types in the Technical Specification for Concrete Structures of High-rise Buildings[15] shows that the finite element model is properly set.

The numerical simulation results of wind-induced vibration response for two working conditions (high-rise buildings with and without podiums) are visualized and compared through post-processing, as shown in Figure 6 to Figure 9. Figure 6 shows the variation of top displacement over time and the comparison of maximum displacement for each floor. Figure 7 shows the comparison of the highest acceleration of each floor and the variation of the top floor acceleration over time. Figure 8 shows the comp`arison of maximum shear forces on each floor and the variation of bottom shear force over time. Figure 9 depicts the time history of the bottom overturning moment and the comparison of the maximum overturning moment on each floor.

When looking at Figure 6(a), it is clear that there aren’t many differences between the two working conditions, as evidenced by the displacement time history curves of the top level of the building with and without a podium. Even still, the peak of the curve shows that the curve without a podium has a slightly bigger overall numerical magnitude than the one with a podium, showing that the displacement pulsation of the top floor of the building without a podium is larger.

In both operating situations, the maximum inter-story displacement can be seen in Figure 6(b) to significantly increase as the floor height increases. This is in keeping with the general rule of wind load on high-rise buildings since as a building’s height increases, so does its flexibility, making it more susceptible to wind stress. When the displacement curves for the two working conditions are compared, it is clear that the maximum displacement between each floor in the working condition without a podium is greater than that in the working condition with a podium, and the difference between the two gradually grows as the floor increases.

In general, adding a podium at the base of a high-rise building helps to reduce the wind-induced vibration displacement response of the entire structure because the podium structure at the base increases the building’s bottom stiffness to some extent, increasing the building’s overall wind resistance.

Figure 7(a) demonstrates that the acceleration time history curves of the top level of the building with and without the podium are still very comparable in both situations, much like the displacement time history scenario. The maximum acceleration difference between the two working conditions does not increase with the height of the building, but instead begins at zero and increases until it reaches its maximum value around the tenth floor, after which it decreases until the two curves essentially coincide, and finally, the acceleration of the working condition without a podium is lower than the acceleration of the working condition with a podium, as can be seen from the detailed diagram in the lower right corner of the Figure 7 (b). The stiffness difference is what drives the curve’s development tendency. The inter-story acceleration of a building drops as its stiffness rises, and the inter-story acceleration of a building with a podium condition is lower from the bottom level than it would be without a podium condition because the podium makes the bottom of the building stiffer. The influence of the bottom podium on the building’s interstory stiffness then steadily reduces as the number of floors rises. Overall, the podium predominantly suppresses inter-story acceleration in the middle and lower floors, and the strength of the suppression weakens as the number of storeys rises.

The time history curves for bottom shear in the two operating circumstances are essentially the same, as can be seen in Figure 8 (a). The difference between the peaks of the time history curves for the cases with and without podiums demonstrates that the effect of the podium on the building’s bottom shear is not immediately apparent. The maximum acceleration between floors in both working circumstances drops noticeably with an increase in floor height, as shown in Figure 8(b). The difference in maximum acceleration between floors is primarily visible in the lower floors, and the inter-story shear force without podium condition is greater than that with podium condition. However, the difference in shear force between the two does not increase with the height of the building; instead, it starts at zero until the maximum value is reached around the tenth floor, then it gradually decreases until the two curves essentially coincide. Overall, the middle and low levels still represent the majority of the podium’s effect on the building’s inter-story shear, and the strength of the impact will continue to diminish as the number of storeys rises.

Because there is a transformation relationship between the shear force and the overturning moment, it can be seen from Figure 9(a) that the shape of the bottom overturning moment time history curve under the two working conditions is essentially the same as the bottom shear time history curve in Figure 8(a). This indicates that the podium’s influence on the building’s bottom overturning moment is as minimal as the bottom shear force.

Figure 9(b) demonstrates that the amplitude and variety of the two curves’ differences are still comparable to Figure 8(b). This means that rather than increasing as a building’s height rises, the difference between the two overturning moments instead increases from 0 until it reaches its greatest value around the 10^th story, then reduces until the two curves almost converge. The scale in the Figure is too tiny and the overturning moment values are quite great, which reduces the distance between the two, causing the curves to be very similar.

Overall, the middle and lower floors still reflect the majority of the podium’s shear stress on the overturning moment between the building’s floors, and the severity of the effect reduces as the number of stories rises.

To better analyze the data, important nodes in the building model were selected for quantitative comparison to highlight the podium effect. The data comparison is shown in Table 5. From Table 5 it can be seen that the average values of various parameters for buildings with podiums and the parameters for the middle floor are lower than those for buildings without podiums under the same conditions. But in the case of a podium on the top floor, except for the top floor displacement being smaller than that without podium, all other parameters are higher than those without a podium. It can be inferred that one of the reasons for this phenomenon is the relative weakening of the stiffness effect of the upper part of the building.

4.2. The Influence of Podium Height on Wind-Induced Vibration Results

Modal analysis was conducted on the finite element models of high-rise buildings with podium heights of 10 meters, 17 meters, and 24 meters, respectively. The analysis is shown in Table 6. This indicates that the vibration modes of each order basically meet the standard of the first six natural vibrations of high-rise buildings, proving that the finite element model is properly set.

Figure 10 shows the time history of the top floor displacement and the maximum displacement of each floor under the three working conditions. Figure 11 shows the time history of the top floor acceleration and the maximum acceleration of each floor. Figure 12 shows the time history of the bottom shear force and the maximum inter-story shear force of each floor. Figure 14 shows the time history of the bottom overturning moment and the maximum overturning moment of each floor.

It can be seen from Figure 10(a) that with the increase of the podium height, the overall value of the displacement time history curve of the top floor is gradually decreasing, and from the peak value of the curve under the three working conditions, the overall amplitude of the curve is also decreasing. This shows that the higher the podium, the smaller the displacement fluctuation amplitude of the top floor of the building. By observing Figure 10(b), it can be seen that with the increase of the building floor, the maximum inter floor displacement under the three working conditions shows an obvious increasing trend. Secondly, comparing the maximum displacement curves between the floors under the three working conditions, it can be clearly seen that the maximum displacement between the floors of the building gradually decreases with the increase of the podium height, and the curve difference between the three working conditions gradually increases with the increase of the floors.

The podium structure at the bottom of the high-rise building can boost the bottom stiffness of the building to a certain level, improving the building’s total wind resistance, according to the research of the working conditions with and without the podium. Similar to this, the research of the podium height condition revealed that as the podium height increased, the stiffness of the building’s main body increased. As a result, the building’s inter-story displacement decreased and the displacement pulsation amplitude increased.

As the podium’s height increases, it is evident from Figure 11(a) that both the amplitude and period of the pulsation of the acceleration time history curve are gradually diminishing. As can be observed from Figure 11(b), the acceleration rises noticeably as the floor height rises, and the maximum acceleration curve between floors follows the expected trend. The acceleration difference between each curve, which starts at zero and increases until it reaches its maximum value around the tenth floor, then gradually decreases until the two curves essentially coincide, can be seen in Figure 11(b). During this time, the maximum acceleration magnitude between building floors gradually decreases with the increase in podium height. The acceleration numbers are in the opposite order of magnitude as before starting at about the 45^th story of the building, and the difference in acceleration is growing. The detailed Figure in Figure 12(b)’s lower right corner makes it easy to determine the magnitude difference between the acceleration values. Following analysis of the podium height scenario, it was discovered that elevating the podium can successfully reduce inter-story acceleration in the building’s middle and lower floors, similar to the cases with and without a podium. Reaching the middle and upper stories, this suppression effect is totally reversed as the floors rise higher.

Figure 11. Inter-storey displacement under podium height condition.

Figure 11. Inter-storey displacement under podium height condition.

The bottom shear time history curves for the three working conditions are extremely similar, as can be shown in Figure 12(a). The three have different pulsation amplitudes, but there isn’t a very noticeable difference between them, so it can be concluded that altering the podium’s height doesn’t significantly affect the building’s bottom shear.

As the floors rise higher, it is clear from Figure 12(b) that the maximum inter-story shear forces for the three operating conditions only differ marginally at the lower floors. The shear force difference between the curves in the Figure is not as noticeable as it would be in the case with and without the podium because the value of the shear force at the base of the building does not vary significantly as a result of the little increase in podium height. However, as can be seen from the detail diagram in the lower right corner of Figure 12(b), the maximum inter-story shear force close to the intersection level of the main structure of the building and the podium gradually reduces with an increase in podium height.

Figure 12. Inter-storey shear under podium height condition.

Figure 12. Inter-storey shear under podium height condition.

Figure 13(a) reveals that the bottom overturning moment time history curve’s trend and period for the three working conditions are essentially identical to those of Figure 11(a). This is due to the fact that shear force and overturning moment have a transformation connection that is comparable to the situation in operating conditions with and without podium. Because the overturning moment values in Figure 13(b) are quite high and the shear forces between the layers for the three working conditions have been somewhat near. The three curves’ differences are thus barely discernible and even close to overlapping. The inter-story overturning moment of the middle and lower floors of the building does, however, gradually decrease with an increase in podium height, as shown by the detail diagram in Figure 13(b)’s right corner.

Figure 13. Overturning bending moment between building floors under podium height condition.

Figure 13. Overturning bending moment between building floors under podium height condition.

Similarly, in order to better analyze the data, important location nodes in the building model were selected for quantitative comparison to highlight the effect of podium height. The data comparison is shown in Table 7.

From Table 7, it can be seen that as the height of the podium increases, the average values of various parameters and the parameters of the middle floor decrease. But at the top level, except for displacement, all other parameters increase with the height of the podium building. It can be inferred that one of the reasons for this phenomenon is still the relative weakening of the stiffness effect on the upper part of the building.

4.3. The Influence of Podium Area on Wind-Induced Vibration Results

The finite element models of the three podium area working conditions perform modal analysis, The specific settings for the three area conditions are shown in Table 1. Perform wind vibration response calculation and post-processing analysis on finite element models under three types of podium floor area conditions.

Following the time history analysis, Figure 14, Figure 15, and Figure 16 show the time history of the top floor acceleration and the maximum acceleration of each floor under the three working conditions, the time history of the bottom shear force and the maximum inter-story shear force of each floor under the three working conditions, and Figure 19 shows the time history of the bottom shear force.

The displacement time history curve for the top floor of the building under the three operating conditions has a very similar trend and period, as can be shown in Figure 14(a). The overall value of the displacement time history curve of the top level, however, gradually declines as the podium’s surface area grows. Additionally, under the three operating conditions, the overall amplitude of the curve is likewise falling from the apex of the curve. This suggests that the displacement pulsation amplitude of the building’s top floor decreases as the podium’s area increases. The maximum displacement between floors in all three sets of working situations grows dramatically with the addition of building floors, as can be seen by first looking at Figure 14(b). Second, when comparing the three curves, it is clear that as the height of the podium rises, the maximum displacement between floors of the building progressively reduces, and as floors rise, the disparity in the curves between the three groups of working circumstances gradually grows.

Combining the two sets of displacement diagrams, it is clear that expanding the podium area in the main podium one-story high-rise building is a good way to reduce the building’s overall response to wind-induced vibration displacement. This is because the inter-story displacement of the building will decrease and the displacement pulsation amplitude will weaken as the size of the podium increases, increasing the stiffness of the main body of the building. This situation is comparable to the former working condition, and the outcomes are as anticipated.

Figure 14. Inter-storey displacement under podium area condition.

Figure 14. Inter-storey displacement under podium area condition.

As demonstrated in Figure 15(a), with an increase in podium floor area, the acceleration time history curves’ pulsation amplitude increases. As can be observed from Figure 15(b), the maximum acceleration curves between floors grow significantly as the number of floors increases, and these changes essentially fulfill the expected requirements. The acceleration difference between each curve doesn’t increase as the building height rises; instead, it increases from zero until it reaches its maximum value somewhere around the tenth floor, then it starts to decline until the two curves essentially coincide, during which point the acceleration magnitude decreases as the podium area grows. Beginning at the 45^th level, the magnitude sequence of the acceleration values is reversed from what it was, followed by a growing discrepancy in the acceleration values, which lasts all the way to the building’s top floor. The detail graphic shows the magnitude and sequence of the acceleration values close to the top floor.In general, expanding the podium’s area can successfully reduce the building’s middle and lower floors’ inter-story acceleration. Reaching the middle and upper floors, however, this suppression effect will be completely contrary to the prior one, while the difference is still minimal. Similar results to this one have been seen in previous working conditions.

Figure 15. Inter-storey displacement under podium area condition.

Figure 15. Inter-storey displacement under podium area condition.

The bottom shear time history curves for the three working conditions are quite similar to one another, as can be seen in Figure 16(a). It can be seen that the podium height of a slight modification for the building bottom shear does not have a very visible impact since the overall difference in the pulsation amplitude of the three is not very great.

The maximum acceleration between floors under the three operating circumstances only minimally varies at the lower floors, as can be shown in Figure 16(b), as the floor area increases. Due to the relatively slight change in the area of the podium, the change in the shear force at the base of the building is not obviously evident. However, as can also be seen from the detail diagram in the lower left corner of Figure 16(b), the value of the largest inter-story shear force in the area around the location where the building’s main structure and the podium meet diminishes as the area of the podium increases.

Figure 16. Inter-storey shear under podium area condition.

Figure 16. Inter-storey shear under podium area condition.

It is evident from Figure 17(a) that the bottom overturning moment time history curves under the three working conditions have a trend and period that are essentially the same as the bottom shear time history curves in Figure 16(a) due to a transformation relationship between the shear force and the overturning moment. The value of the overturning moment at the base of the building does not vary obviously because the change in the area of the podium is rather minor.

In Figure 17 (b), the distinctions between the three curves are not immediately apparent and have even come very close to overlapping since the shear forces between the layers in the three operating conditions have been extremely close and the values of the overturning moments have been quite large. The inter-story overturning moment of the building’s middle and lower floors, however, nevertheless gradually reduces with an increase in podium area, as shown by the detailed subfigure in right corner of Figure 17(b), and the law is compatible with this change in the condition of the podium height.

Figure 17. Overturning bending moment between building floors under podium area condition.

Figure 17. Overturning bending moment between building floors under podium area condition.

The following is a summary of the research results on different podium area working conditions. Nodes from important parts of each model are selected for quantitative comparison, highlighting the podium effect. The data comparison is shown in Table 8.

From Table 8, it can be seen that as the area of the podium increases, the average values of various parameters and the parameters of the middle floor decrease. But at the top level, except for displacement, all other parameters increase with the increase of the podium area, which is similar to the first two working conditions, once again verifying the reliability of the simulation results. And the three parameters of acceleration, shear force, and bending moment decrease with the increase of podium area at the lower level of the building, but increase with the increase of podium area at the middle and upper levels, showing the opposite variation. It can be inferred that one of the reasons for this phenomenon is still the weakening of the stiffness effect of the upper part of the building.

5. The Influence of Podium Layout on Wind-Induced Vibration Results

5.1. Calculation Models for Several Podium Layouts

In all layout methods, the size of the podium is 60m x 60m x 10m, and the cross-sectional size of the main body is the same as that of the building without a podium. The overall height of the building is 90m. C1, C2, and C3 are symmetrical layouts, while C4, C5, and C6 are asymmetrical layouts. The model is shown in Figure 18.

Figure 18. Finite element model of working conditions for podium layout.

Figure 18. Finite element model of working conditions for podium layout.

The top views of the six models presented in this section are shown in Figure 19.

Figure 19. Top view of architectural models with different podium layouts.

Figure 19. Top view of architectural models with different podium layouts.

The loading method of pulsating wind load (nodal dynamic load) on the attached podium building is shown in Figure 20.

Figure 20. Schematic diagram of fluctuating wind load loading on the building of conditions for podium layout.

Figure 20. Schematic diagram of fluctuating wind load loading on the building of conditions for podium layout.

5.2. Analysis of Numerical Simulation Results of Wind-Induced Vibration Response

The maximum displacement of each layer under six working conditions is shown in Figure 21, and the values of the maximum displacement, average displacement, and the ratio of maximum displacement to average displacement of key layers are shown in Table 9 Displacement of key floors in the conditions of podium layout Comparative analysis of figure and table data, it can be seen that the maximum displacement of symmetrical layout always increases with the height of the floor, while the displacement of asymmetrical layout undergoes a sudden change at the junction of the podium and the main body, and the maximum displacement of asymmetrical layout is generally greater than that of symmetrical layout.

Figure 21. Comparison diagram of displacement(S) of buildings in the conditions of podium layout.

Figure 21. Comparison diagram of displacement(S) of buildings in the conditions of podium layout.

The maximum shear force and maximum overturning moment of each floor under six working conditions are shown in Figure 22 and Figure 23, respectively. The values of maximum shear force and overturning moment of key floors are shown in Table 8 From the Figure, it can be seen that the maximum overturning moment and maximum shear force trend of the six working conditions are similar, but the difference is that the maximum shear force graph is convex on the top and concave on the bottom. Like the maximum displacement, the maximum shear force curves and maximum overturning moment curves of C1 and C3, C4 and C6 completely overlap. Except for the layout at the center, the maximum shear force curves of the other five layouts have small differences and are not smooth at the junction between the podium and the main body. The maximum shear force and maximum overturning moment of the layout at the center of the high-rise area are smaller than those of the other working conditions, while the values of the layout at the center of the low rise area are relatively large.

Figure 22. Comparison diagram of shear forces of buildings in the conditions of podium layout.

Figure 22. Comparison diagram of shear forces of buildings in the conditions of podium layout.

Figure 23. Comparison diagram of overturning moment of buildings in the conditions of podium layout.

Figure 23. Comparison diagram of overturning moment of buildings in the conditions of podium layout.

From the above analysis, it can be concluded that the overall displacement of symmetrical layout is smaller than that of asymmetrical layout. It is speculated that this is due to the relative weakening of stiffness caused by eccentric layout, resulting in an increase in displacement. By observing the maximum shear force and maximum overturning moment of each floor, it can be concluded that the design of the central podium layout can appropriately reduce the maximum shear force and overturning moment of the higher floors of the building, but the effect is opposite in lower floors.

6. Conclusions

In order to study the wind-induced response characteristics of the main podium one-piece high-rise building structure, this paper presents a process algorithm for solving the pulsating wind speed time history and determining the pulsating wind load based on the AR method. It then applies this algorithm to the study of the wind-induced response characteristics of the structure, and the following conclusions are drawn.

(1) It is discovered that the building’s response to wind-induced displacement is mostly impacted by the podium’s existence or absence. The middle and bottom floors of the building exhibit the most visible effects of the podium’s presence or absence on the acceleration response, bottom shear, and bottom overturning moment, and the influence diminishes with increasing building height.

(2) The maximum displacement between floors of the building steadily reduces as the podium height rises, which is evident from the building’s response to wind-induced displacement. The influence of wind-induced vibration acceleration is most noticeable in the middle and lower floors of the building structure, and it gradually diminishes as the number of floors rises. At levels close to the top floor of the structure, the inter-story acceleration rises with the height of the podium. With an increase in podium height, the building’s bottom shear force and bottom overturning moment somewhat decrease.

(3) Expanding the podium’s area can also successfully reduce the building’s displacement due to wind, and this effect is particularly pronounced on the upper stories. Even while expanding the podium’s surface area can reduce the rate of acceleration on the middle and lower floors, the suppression effect is less noticeable on the upper floors. In line with other working circumstances, as the podium area is increased, the shear force and overturning moment of the building’s middle and lower floors diminish.

(4) The layout of the podium has an impact on the wind-induced vibration response of the main building. When the podium and main building are symmetrically arranged in the downwind direction, the maximum displacement of each floor is small, while the maximum displacement curve of buildings with asymmetric layout at the junction of the podium and main building is not smooth. The design of the central layout of the podium and main building can effectively reduce the maximum shear force and maximum overturning moment of the higher floors of the building, but the effect is opposite at lower floors.