An innovative analytical approach for predicting the fundamental time period of moment-resisting frames

Eman SharafDepartment of Civil Engineering, Misr Higher Institute for Engineering and Technology, Mansoura, EgyptStructural Engineering Department, Faculty of Engineering, Zagazig University, Sharkia Governorate, Zagarig 44519, Egypt.e.sharaf22@eng.zu.edu.eg, emanzaki2050@gmail.comAtef Eraky, Abdallah Salama, Shimaa EmadStructural Engineering Department, Faculty of Engineering, Zagazig University, Sharkia Governorate, Zagazig 44519, Egypt.aeamien@eng.zu.edu.eg, ASAbdullatij@eng.zu.edu.eg, seattia@zu.edu.eg

Introduction

The fundamental period of vibration of a building is one of the important parameters appearing in the different building codes to compute the design base shear and lateral loads. The empirical formulae found in codes are mostly based on the building type (e.g., frame or shear wall), overall size, and construction material (e.g., steel or reinforced concrete). The determination of this period is critical since it has a direct impact on the calculation of seismic loads for improving both the accuracy and reliability of the building design process.
Almost all the formulas currently found in the "Uniform Building Code" (UBC) [1], the "Structural Engineers Association of California" (SEAOC) [2] recommendations, and the "Egyptian Code" (EGC) [3] are primarily derived from building vibration data recorded during the 1971 San Fernando earthquake. Later, in 1996 and 1997, Goel and Chopra [4, 5] evaluated existing empirical formulas and subsequently proposed new empirical formulas to estimate the fundamental vibration period of reinforced concrete moment-resisting frame (RC MRF) buildings. These methods were developed for use in equivalent lateral force analysis, as required by building codes, utilizing motion data from various structures recorded during multiple earthquakes. The data analyzed by [4] were obtained from earthquakes, including the 1971 San Fernando, 1984 Morgan Hill, 1986 Mt. Lewis and Palm Springs, 1987 Whittier, 1989 Loma Prieta, 1990 Upland, 1991 Sierra Madre, and 1994 Northridge events.
Salama [6] also updated formulas for the estimation of the fundamental vibration period by using regression analysis on data from eight recorded earthquakes in California. Their formulas account for both the building height and the number of stories in the moment-resisting reinforced concrete frames. Several other researchers have also proposed an adjustment to the Ct coefficient based on the number of stories, according to the formulas in American and Egyptian building codes. Similarly, Kalpan et al. [7] emphasized that the empirical formulas in building codes should be based on the particular region in which the typical design and architectural character of local construction are exhibited. In this regard, they derived a simplified formula to estimate the fundamental elastic period of mid-rise reinforced concrete buildings and compared it with the formula available in the "Turkish Earthquake Code" (TEC) [8]. Their study indicated that the empirical formula in the code yields non-conservative base shear values. They also investigated the contribution of infill walls to the lateral stiffness of buildings and offered preliminary recommendations based on their findings.
Several researchers, including [ 9 , 10 , 11 , 12 , 13 , 14 , 15 , 16 ] [ 9 , 10 , 11 , 12 , 13 , 14 , 15 , 16 ] [9,10,11,12,13,14,15,16][9,10,11,12,13,14,15,16][9,10,11,12,13,14,15,16], have derived simplified equations based on building height by using ambient vibration measurements to estimate the fundamental period of buildings. Mohamed et al. [17] conducted nonlinear dynamic analyses to estimate the fundamental period of mid-rise moment-resisting RC frames. In their study, the results were compared with the period equations given by Salama and "American Society of Civil Engineers: Minimum Design Loads for Buildings and Other Structures" (ASCE 7-16) [6, 18]. They developed several formulas from the regression analysis of their test results to estimate the period of vibration. They also supported the application of rational approaches,
like modal analysis, to calculate the fundamental period of RC frames and, in this way, remove some of the restrictive codebased limits.
Young and H. Adeli [19, 20] criticized the methodology for determining the estimated fundamental period of momentresisting frames obtained from ASCE 7-16 [18] as being overly conservative. They emphasized that this approach doesn't fully account for geometric irregularity in structures, which may greatly affect the structural responses during seismic events. Due to this deficiency, an extensive study was carried out on the fundamental period of moment-resisting frames with various types of geometric irregularities. They developed several formulas that took into consideration both vertical and horizontal irregularities; hence, these formulas became more accurate and flexible methods of estimating the fundamental period.
The degree of irregularity in a stepped building frame is estimated in [21], using the Rayleigh approach , and suggested modifications to the empirical formula used in building codes for estimating the fundamental period of such stepped frames. However, the calculation of the regularity index, as proposed by them, is complex and time-consuming. Ricci et al. [22] studied the elastic vibration period of uncracked and infilled reinforced concrete moment-resisting frame buildings. They employed a modal analysis approach to assess the structural response and subsequently offered simplified formulas to estimate the fundamental period. These formulas were based on regression analysis of numerical data, providing a practical alternative for efficient computation.
An equation was proposed for estimating the fundamental vibration period of regular frame structures with a maximum of ten stories, aiming to improve the empirical expressions widely used in building codes [23]. They utilized structural features such as stiffness and mass distribution, which enable more accurate estimations to be obtained by getting a correction factor based on effective mass and height, and the effective lateral stiffness of the structure. They also saw its applicability to fairly irregular frame structures with minimal loss of accuracy.
Alrudaini [24] critically examined the estimation of the fundamental vibration period ( T ) of reinforced concrete momentresisting frame (RC-MRF) buildings, noting the limitation in existing seismic design codes. Based on regression analysis of recorded motions from RC moment-resisting frame (RC-MRF) buildings during eight California earthquakes, an improved formula was proposed in [24] that implicitly accounts for the floor height (h) as a significant parameter by incorporating both the total building height ( H ) ( H ) (H)(\mathrm{H})(H) and the number of stories ( N ) ( N ) (N)(\mathrm{N})(N). The study suggests the modification of the coefficient Ct in the current code equations as a function of N to enhance the precision of T estimates.
Current studies have focused more on improving the precision of fundamental period estimation of steel and reinforced concrete structures, particularly in terms of complex structural properties such as height, irregularity, and bracing systems. Kumar et al. [25] performed a comparative analysis of machine learning algorithms, including artificial neural networks (ANN), genetic programming (GP), and regression trees (RT), demonstrating that ANN models exhibited superior predictive accuracy, especially in the presence of geometric irregularities.
Likewise, Rahman et al. [26] employed interpretable machine learning methods with SHapley Additive exPlanations (SHAP) to determine the significant structural parameters influencing the fundamental period in steel-braced reinforced concrete buildings and the way its configurations have significant influences on its dynamic behavior. In order to accurately capture the uncertainties associated with aging and material degradation, Shan et al. [27] developed probabilistic machine learning techniques to forecast the natural period of existing high-rise reinforced concrete buildings. Moreover, Karampinis et al. [28] created hybrid analytical-machine learning models that used SHAP values to generate interpretable equations for framed structures. This allowed them to preserve engineering transparency while achieving prediction accuracy on par with purely data-driven methods. These developments highlight a growing movement to improve period estimation, especially for complex and irregular structures, by combining data-driven techniques with conventional analytical models.
Most seismic design codes recommend empirical equations for the estimation of the fundamental vibration period of buildings. However, it has been demonstrated through several studies that these code-based equations typically underestimate significantly the actual measured periods, as stipulated by Goel and Chopra [4]. They once again highlighted the inadequacy of the code provisions available by citing that the simple equations that are commonly used can vary as much as 100 % 100 % 100%100 \%100% from observed values. These broad variations can lead to excessively conservative estimates of base shear force, thus structurally over-strengthened and uneconomical designs. Although several equations have been presented both in design codes and in the literature, many are empirical and usually derived from regression analyses of small datasets. These formulations often neglect the physical parameters most critical to dynamic behavior. Although many recent studies have also looked into data-driven and machine learning strategies can provide high predictive accuracy, their applicability and reliability in engineering practice are limited by the need for large, high-quality datasets and the possibility of overfitting or a lack of physical interpretability. Furthermore, they can be computationally intensive and require specialized knowledge to implement, which makes them less useful for everyday design. Hence, there is an evident necessity for a more logical and
physically based model, including the equivalent mass and stiffness, two of the most significant parameters that control the fundamental period, precisely to increase both accuracy and usability of period estimation techniques.
This paper proposes a novel analysis procedure for the accurate estimation of the fundamental period of vibration of reinforced concrete moment-resisting frames. The equivalent seismic mass and the building's lateral stiffness are both included in the proposed procedure for higher accuracy than the current procedure. This approach differs from the previous approaches because it provides an easier process to compute the equivalent stiffness and equivalent mass of a multi-story building and also computes its fundamental period with accuracy, which gives more accurate results, and the resulting equations are in the form of closed-form expressions. Sensitivity analysis was carried out to investigate the influence of significant design parameters such as building height, span length, beam-to-column stiffness, and material elasticity. This procedure not only serves to confirm the strength and validity of the proposed equation but also establishes the key parameters that need to be taken into consideration in the design. It confirms the proposed model as being stable and representative of a wide range of buildings, thus making it more generalizable and usable in practice.

Formulas from the codes and literature

The empirical formulas for the fundamental vibration period of moment-resisting frames (MRFs), adopted by most design codes, including U.S. building codes, UBC [1], the "Applied Technology Council" (ATC,) [29], SEAOC, [2], the "National Earthquake Hazards Reduction Program" (NEHRP) [30], and the more recent EGC [3], are generally in the form:
(1) T a = C t H 0.9 (1) T a = C t H 0.9 {:(1)T_(a)=C_(t)H^(0.9):}\begin{equation*} T_{a}=C_{t} H^{0.9} \tag{1} \end{equation*}(1)Ta=CtH0.9
In this empirical formula, H represents the total height of the building in meters, which is measured as the vertical distance from the ground level or foundation. Ct represents a numerical coefficient related to the lateral force-resisting system. Values adopted for Ct by the above codes are 0.047 for RC MRF and 0.051 for Steel MRF buildings. For some design codes, such as the NEHRP [30] regulations and earlier versions of other seismic codes, including the EGC [3], an alternative formula is available for RC MRF buildings. This formula expresses the fundamental period as:
(2) T a = 0.1 N (2) T a = 0.1 N {:(2)T_(a)=0.1N:}\begin{equation*} T_{a}=0.1 \mathrm{~N} \tag{2} \end{equation*}(2)Ta=0.1 N
The simple formula is restricted to structures with a height limit of a story of 10 feet and a height of not more than 12 stories, where N denotes the number of stories. Salama [6] provides the estimation of the fundamental period of vibration (T) of reinforced concrete moment-resisting frame buildings, one of the more important parameters in seismic design. Conventional seismic codes, such as ASCE [18] and EGC [3], estimate T using empirical equations that account for either building height ( H ) or the number of stories ( N ) individually, not both together. An improved empirical formula that includes both H and N , based on regression analysis of the data for moment-resisting frame reinforced concrete buildings experiencing eight California earthquakes was suggested in [6]. The method is intended to give an improved estimate of T by considering variations in floor heights and the number of stories. The study recommends that the coefficient Ct of the current code formulas be modified for the effect of N in order to make the correlation between estimated and observed vibration periods even better.
Salama [6] adopted a theoretical form of the period equation that included both the number of stories ( N ) and the total height ( H ) ( H ) (H)(\mathrm{H})(H) of the building
(3) T = a N b H c (3) T = a N b H c {:(3)T=a*N^(b)*H^(c):}\begin{equation*} T=a \cdot N^{b} \cdot H^{c} \tag{3} \end{equation*}(3)T=aNbHc
in this approach, the constants a , b a , b a,b\mathrm{a}, \mathrm{b}a,b, and c as variables in an unconstrained regression analysis, allowing the optimal fit of statistical data to measured data. This method minimized the standard error of the estimate, resulting in:
(3-a) T = 0.027 N 0.17 H 0.74 (3-a) T = 0.027 N 0.17 H 0.74 {:(3-a)T=0.027*N^(0.17)*H^(0.74):}\begin{equation*} T=0.027 \cdot N^{0.17} \cdot H^{0.74} \tag{3-a} \end{equation*}(3-a)T=0.027N0.17H0.74
where T is the fundamental period of vibration (in seconds), N is the number of stories, and H is the height of the building (in meters). Aninthaneni et al.[23] derived an improved analytical method to forecast the fundamental period (T) of regular
reinforced concrete (RC) moment-resisting frame buildings in the Bulletin of the New Zealand Society for Earthquake Engineering. The recommended equation estimates T more correctly than empirical equations by considering seismic weight, storey height, number of storeys, bays, and effective stiffness. The method is validated against eigenvalue analysis and experimental findings on low- to mid-rise RC frames, even with minor deviation.
The improved empirical formulas proposed are:
(4) T = 0.47 φ 3 W s b 3 n s [ ( 1 β d ) 3 + λ ( 1 β c ) 3 ] ( n b + 1 ) E c I c e f (4) T = 0.47 φ 3 W s b 3 n s 1 β d 3 + λ 1 β c 3 n b + 1 E c I c e f {:(4)T=0.47*varphi_(3)*sqrt((W_(s)*b)/(3n_(s)[(1-beta_(d))^(3)+lambda(1-beta_(c))^(3)](n_(b)+1)E_(c)I_(cef))):}\begin{equation*} T=0.47 \cdot \varphi_{3} \cdot \sqrt{\frac{W_{s} \cdot b}{3 n_{s}\left[\left(1-\beta_{d}\right)^{3}+\lambda\left(1-\beta_{c}\right)^{3}\right]\left(n_{b}+1\right) E_{c} I_{c e f}}} \tag{4} \end{equation*}(4)T=0.47φ3Wsb3ns[(1βd)3+λ(1βc)3](nb+1)EcIcef
In this equation, ϕ 3 ϕ 3 phi_(3)\phi_{3}ϕ3 is a correcting factor that takes into account the effective distribution of mass and height in the structure. The variable Ws is the total seismic weight of the frame, and h h hhh is the storey height, or the height between two consecutive floors. The storey number is denoted by n s n s n_(s)\mathrm{n}_{\mathrm{s}}ns, meanwhile, n b n b n_(b)\mathrm{n}_{\mathrm{b}}nb is the total number of bays in the frame. The parameters β d β d beta_(d)\beta_{\mathrm{d}}βd and β c β c beta_(c)\beta_{c}βc are the depiction of the distribution of stiffness over the building height for the beam and column systems, respectively. The ratio of the beam-to-column stiffness is given by λ λ lambda\lambdaλ, and the material properties are described by E c E c E_(c)\mathrm{E}_{\mathrm{c}}Ec, the modulus of elasticity of the concrete, and I cef I cef I_(cef)\mathrm{I}_{\mathrm{cef}}Icef, the effective moment of inertia of the columns. The model enables the prediction of the fundamental period of the building, considering both material and geometric properties, with better accuracy compared to standard empirical approaches for common moment-resisting frame-type structures.
Goel and Chopra [4] experimentally validated the equations provided in U.S. seismic codes for approximating the fundamental period of vibration (T) of buildings. Based on the assessment of observed periods of 106 buildings that experienced eight California earthquakes from 1971 through 1994, they discovered that the current code equations underestimate the fundamental time period. By regression analysis of this expanded database, [5] developed better empirical equations that more accurately estimate the natural periods of reinforced concrete (RC) and steel moment-resisting frame (MRF) buildings. These improved equations had a significant impact in revising seismic design provisions in, for example, the NEHRP [31] and ASCE 7-05 [32] standards, and are expressed as:
(5) T = 0.0466 H 0.9 (5) T = 0.0466 H 0.9 {:(5)T=0.0466H^(0.9):}\begin{equation*} T=0.0466 H^{0.9} \tag{5} \end{equation*}(5)T=0.0466H0.9
Here, T is the fundamental period of vibration of the building in seconds, and H is the height of the building in meters. Most current design codes specify that the design base shear should be calculated using the following equation:
(6) V = C W (6) V = C W {:(6)V=CW:}\begin{equation*} V=C W \tag{6} \end{equation*}(6)V=CW
where V is the design base shear, W is the total seismic dead load, and C is the seismic coefficient. The seismic coefficient C depends on several parameters: the soil profile, the seismic zone factor, the importance factor, the building's fundamental period T , and a numerical coefficient representing overstrength and global ductility capacity inherent in the lateral loadresisting system. The fundamental period Ta calculated by the empirical formulas from equations (1) and (2) is typically designed to be shorter than the true period in order to provide a conservative estimate of the base shear, as structures are required to be designed for larger seismic forces. Thus, code formulas are intentionally calibrated to underestimate the fundamental period by about 10-20%. Furthermore, while building codes allow the determination of the fundamental period using rational methods, such as Rayleigh's method [33], they also specify that the period computed based on these methods shall not exceed the period obtained from the empirical formula by more than a certain factor. In this way, consistency in seismic design is achieved, while there is provision for making more accurate calculations when required.
Figure 1: Comparaison of FEM and code periods for RC MRF buildings.

EFFECTIVENESS OF FORMULAS FROM CODES AND LITERATURE

To evaluate the validity of frequently used empirical formulas in identifying the fundamental period of structures, comparative study was performed between code expressions and FEM. The comparison given in Figure 1 shows a detailed assessment of accuracy of different empirical expressions to comprehensive finite element method (FEM) analysis for calculation of the fundamental period of structures for a broad range of heights. Specifically, the comparison incorporates empirical period formulae from UBC [1], ASCE 7-05 [32], and also Salama's equation [6] as written in Equations (1), (5) and (3) respectively. Figure 1 presents several important observations concerning the accuracy of these formulae relative to FEM results:
  • FEM calculations consistently forecast longer fundamental periods than all empirical formulas, particularly for structures exceeding 30 meters in height. This reflects the enhanced flexibility represented in numerical models, which is typically underestimated in empirical formulations.
  • ASCE 7-05 [32] calculates the period using the formula T a = 0.0466 h 0.9 T a = 0.0466 h 0.9 T_(a)=0.0466h^(0.9)\mathrm{T}_{\mathrm{a}}=0.0466 \mathrm{~h}^{0.9}Ta=0.0466 h0.9, with an upper limit adjusted by C u = 1.4 C u = 1.4 C_(u)=1.4\mathrm{C}_{\mathrm{u}}=1.4Cu=1.4. While the upper bound curve somewhat aligns with the FEM results for mid-heights; nevertheless, it significantly underestimates periods for high-rise buildings, indicating inadequate applicability for tall structures.
  • The ASCE 7-05 [32] code formula produces significantly shorter period estimates than those obtained from FEM analysis, particularly for buildings exceeding 45 m in height. At 90 m , the percentage error is approximately 45 % 45 % -45%-45 \%45%, reflecting a rather large underestimation of the fundamental period and that the code may not capture accurately the dynamic behaviour of high-rise buildings.
  • The UBC [1] predictions are relatively aligned with FEM results for low- to mid-rise structures (about 15 40 m 15 40 m 15-40m15-40 \mathrm{~m}1540 m ), but diverge significantly as building height increases. It is obvious that the UBC formula is relatively adequate up to a height of 40 m , but is not reliable in the case of high-rise buildings.
  • Salama [6] presents the most conservative of the three estimates, significantly underestimating the fundamental period for all heights of buildings. This shows that the model is probably missing the structural flexibility and cracking influences, and hence it overestimates the lateral stiffness.
    The discrepancy between empirical predictions and FEM highlights the restrictions of generalized code equations, especially for irregular geometries or taller structures. This warrants calibration of empirical models with newer numerical data and consideration of FEM-based strategies in performance-based seismic design.

Estimation of the Fundamental Period of Vibration

Seismic design typically begins with estimating a trustworthy seismic force, usually based on the periods of the building's vibration. These periods are primarily influenced by the mass and stiffness of the building. In this case, seismic design often becomes a process of trial and error, since the relationship between seismic capacity, mass, and stiffness must be well recognized to ensure effective seismic performance. To avoid unnecessary trial and error, most standards include starting points for element size, construction orientation, number of levels, and material properties. These provide natural period calculations, and many formulas from codes and literature are used to calculate the building's seismic load. The fundamental period of the building, during which a significant portion of the mass participates in the associated mode of vibration, is by definition the natural period of vibration. While more complex formulas are available in building codes, simpler equations based on height or the number of stories are often recommended to streamline the calculation process. Several studies have examined the effect of both structural and non-structural factors on the fundamental period of vibration in reinforced concrete buildings. Consequently, the empirical equations have revealed that the natural period is mostly determined by the building's height or number of stories. More specifically, the taller or higher a building is, the longer the period of vibration will be [34, 35]. These two properties, height and number of stories, are found to be the two most influential ones that determine a building's fundamental period. Consequently, most design codes rely on one of these two properties, preferably height or number of floors, in developing their empirical equations for the vibration period. Notable examples of such codes include the "building standard law of Japan" (BSLJ) [36], the National Building Code of Canada (NBCC, 1995) [37] the Uniform UBC, [1], and the ATC [29], all of which incorporate height or number of stories as key parameters in their period formulas.

A simplified technique for determining the fundamental period

The most practical approach to developing an analytical model of a real-frame structure is to model it as a skeleton with masses lumped at each floor level. In this section, a simplified technique is outlined for calculating the fundamental period of vibration of such structures. As previously mentioned, and as widely recognized, the period of vibration of a building depends on its mass and stiffness. For an idealized Single-Degree-of-Freedom (SDOF) system, the fundamental period is classically expressed by Equation (7), which defines the dynamic relationship between mass and stiffness.
(7) T = 2 π m K (7) T = 2 π m K {:(7)T=2pisqrt((m)/(K)):}\begin{equation*} T=2 \pi \sqrt{\frac{m}{K}} \tag{7} \end{equation*}(7)T=2πmK
Equation (7) forms the base of this study. An RC building can be modeled by a single-degree-of-freedom (SDOF) system using an equivalent lumped mass m eq m eq  m_("eq ")\mathrm{m}_{\text {eq }}meq  and equivalent stiffness k eq k eq  k_("eq ")\mathrm{k}_{\text {eq }}keq  as shown in Fig. 2. This idealization is concerned with the fundamental period of vibration and the interaction between the mass and stiffness in the system, making it much easier to analyze the dynamic behavior of the building.
Procedure for estimating the equivalent lateral stiffness of the frame ( K eq K eq  K_("eq ")K_{\text {eq }}Keq  )
The following explanation may help clarify the development of the current formula. Initially, the building is simplified to achieve a more accurate estimate of lateral stiffness. The proposed method begins by calculating the equivalent stiffness for a two-storey building, as shown in Fig. 3. Using the Rayleigh approach, the process for estimating the fundamental period of the building is formulated in this section.
Figure 2: Equivalent cantilever column.
Figure 3: Configuration of MRF with simplification procedures.
To develop the shear stiffness matrix of the two-storey MRF, considering the stiffness of the beams and anti-symmetry, the flexibility matrix is developed first by applying two horizontal forces separately and calculating the horizontal deformations. The stiffness matrix of the frame due to translational and rotational degrees of freedom takes the form:
(8) K = [ 24 E I c 1 h 1 3 + 24 E I c 2 h 2 3 12 E I c 1 h 1 2 12 E I c 2 h 2 2 24 E I c 2 h 2 3 12 E I c 2 h 2 3 12 E I c 1 h 1 2 12 E I c 2 h 2 2 8 E I b 1 L + 8 E I c 2 h 2 + 8 E I c 1 h 1 12 E I c 2 h 2 2 4 E I c 2 h 2 24 E I c 2 h 2 3 12 E I c 2 h 2 2 24 E I c 2 h 2 3 12 E I c 2 h 2 2 12 E I c 2 h 2 3 4 E I c 2 h 2 12 E I c 2 h 2 2 8 E I c 2 L + 8 E I c 2 h 2 ] (8) K = 24 E I c 1 h 1 3 + 24 E I c 2 h 2 3 12 E I c 1 h 1 2 12 E I c 2 h 2 2 24 E I c 2 h 2 3 12 E I c 2 h 2 3 12 E I c 1 h 1 2 12 E I c 2 h 2 2 8 E I b 1 L + 8 E I c 2 h 2 + 8 E I c 1 h 1 12 E I c 2 h 2 2 4 E I c 2 h 2 24 E I c 2 h 2 3 12 E I c 2 h 2 2 24 E I c 2 h 2 3 12 E I c 2 h 2 2 12 E I c 2 h 2 3 4 E I c 2 h 2 12 E I c 2 h 2 2 8 E I c 2 L + 8 E I c 2 h 2 {:(8)K=[[(24 EIc_(1))/(h_(1)^(3))+(24 EIc_(2))/(h_(2)^(3)),(12 EIc_(1))/(h_(1)^(2))-(12 EIc_(2))/(h_(2)^(2)),-(24 EIc_(2))/(h_(2)^(3)),-(12 EIc_(2))/(h_(2)^(3))],[(12 EIc_(1))/(h_(1)^(2))-(12 EIc_(2))/(h_(2)^(2)),(8EIb_(1))/(L)+(8EIc_(2))/(h_(2))+(8EIc_(1))/(h_(1)),(12 EIc_(2))/(h_(2)^(2)),(4EIc_(2))/(h_(2))],[-(24 EIc_(2))/(h_(2)^(3)),(12 EIc_(2))/(h_(2)^(2)),(24 EIc_(2))/(h_(2)^(3)),(12 EIc_(2))/(h_(2)^(2))],[-(12 EIc_(2))/(h_(2)^(3)),(4EIc_(2))/(h_(2)),(12 EIc_(2))/(h_(2)^(2)),(8EIc_(2))/(L)+(8EIc_(2))/(h_(2))]]:}\mathrm{K}=\left[\begin{array}{cccc} \frac{24 E I c_{1}}{h_{1}^{3}}+\frac{24 E I c_{2}}{h_{2}^{3}} & \frac{12 E I c_{1}}{h_{1}^{2}}-\frac{12 E I c_{2}}{h_{2}^{2}} & -\frac{24 E I c_{2}}{h_{2}^{3}} & -\frac{12 E I c_{2}}{h_{2}^{3}} \tag{8}\\ \frac{12 E I c_{1}}{h_{1}^{2}}-\frac{12 E I c_{2}}{h_{2}^{2}} & \frac{8 E I b_{1}}{L}+\frac{8 E I c_{2}}{h_{2}}+\frac{8 E I c_{1}}{h_{1}} & \frac{12 E I c_{2}}{h_{2}^{2}} & \frac{4 E I c_{2}}{h_{2}} \\ -\frac{24 E I c_{2}}{h_{2}^{3}} & \frac{12 E I c_{2}}{h_{2}^{2}} & \frac{24 E I c_{2}}{h_{2}^{3}} & \frac{12 E I c_{2}}{h_{2}^{2}} \\ -\frac{12 E I c_{2}}{h_{2}^{3}} & \frac{4 E I c_{2}}{h_{2}} & \frac{12 E I c_{2}}{h_{2}^{2}} & \frac{8 E I c_{2}}{L}+\frac{8 E I c_{2}}{h_{2}} \end{array}\right](8)K=[24EIc1h13+24EIc2h2312EIc1h1212EIc2h2224EIc2h2312EIc2h2312EIc1h1212EIc2h228EIb1L+8EIc2h2+8EIc1h112EIc2h224EIc2h224EIc2h2312EIc2h2224EIc2h2312EIc2h2212EIc2h234EIc2h212EIc2h228EIc2L+8EIc2h2]
To obtain the flexibility matrix, the horizontal deformations due to horizontal unit forces are obtained as follows:
(9) f = T T K T (9) f = T T K T {:(9)f=T^(T)KT:}\begin{equation*} f=T^{T} K T \tag{9} \end{equation*}(9)f=TTKT
where T is the transformation matrix that takes the form:
(10) T trans = [ 1 0 0 0 0 1 0 0 ] (10) T trans  = 1 0 0 0 0 1 0 0 {:(10)T_("trans ")=[[1,0],[0,0],[0,1],[0,0]]:}T_{\text {trans }}=\left[\begin{array}{ll} 1 & 0 \tag{10}\\ 0 & 0 \\ 0 & 1 \\ 0 & 0 \end{array}\right](10)Ttrans =[10000100]
Then the shear stiffness matrix of the shown frame is obtained as follows:
(11) K = f 1 (11) K = f 1 {:(11)K=f^(-1):}\begin{equation*} K=f^{-1} \tag{11} \end{equation*}(11)K=f1
The structure mass is condensed to one mass at the first-floor level, then the structure mass matrix takes the form:
(12) M = [ m 1 0 0 0 ] (12) M = m 1 0 0 0 {:(12)M=[[m_(1),0],[0,0]]:}M=\left[\begin{array}{cc} m_{1} & 0 \tag{12}\\ 0 & 0 \end{array}\right](12)M=[m1000]
where m 1 m 1 m_(1)m_{1}m1 is the equivalent mass of the structure.
Then the first natural frequency of this frame is obtained by:
(13) | K ω e q 1 2 M | = 0 (13) K ω e q 1 2 M = 0 {:(13)|K-omega_(eq1)^(2)M|=0:}\begin{equation*} \left|K-\omega_{e q 1}^{2} M\right|=0 \tag{13} \end{equation*}(13)|Kωeq12M|=0
Then the equivalent first-floor stiffness is obtained as follows:
(14) K e q = ω e q 1 2 m 1 (14) K e q = ω e q 1 2 m 1 {:(14)K_(eq)=omega_(eq1)^(2)m_(1):}\begin{equation*} K_{e q}=\omega_{e q 1}^{2} m_{1} \tag{14} \end{equation*}(14)Keq=ωeq12m1
which takes the following form:
(15) K e q = 6 E I c 1 b 1 3 ( I c 1 I c 2 L 2 + 4 I b 1 I c 2 L h 1 + 4 I b 2 I c 1 L h 2 + 4 I b 2 I c 2 L h 1 + 16 I b 1 I b 2 h 1 h 2 I c 1 I c 2 L 2 + I b 1 I c 2 L h 1 + 4 I b 2 I c 1 L h 2 + I b 2 I c 2 L h 1 + 4 I b 1 I b 2 h 1 h 2 ) (15) K e q = 6 E I c 1 b 1 3 I c 1 I c 2 L 2 + 4 I b 1 I c 2 L h 1 + 4 I b 2 I c 1 L h 2 + 4 I b 2 I c 2 L h 1 + 16 I b 1 I b 2 h 1 h 2 I c 1 I c 2 L 2 + I b 1 I c 2 L h 1 + 4 I b 2 I c 1 L h 2 + I b 2 I c 2 L h 1 + 4 I b 1 I b 2 h 1 h 2 {:(15)K_(eq)=(6EIc_(1))/(b_(1)^(3))((Ic_(1)Ic_(2)L^(2)+4Ib_(1)Ic_(2)Lh_(1)+4Ib_(2)Ic_(1)Lh_(2)+4Ib_(2)Ic_(2)Lh_(1)+16 Ib_(1)Ib_(2)h_(1)h_(2))/(Ic_(1)Ic_(2)L^(2)+Ib_(1)Ic_(2)Lh_(1)+4Ib_(2)Ic_(1)Lh_(2)+Ib_(2)Ic_(2)Lh_(1)+4Ib_(1)Ib_(2)h_(1)h_(2))):}\begin{equation*} K_{e q}=\frac{6 E I c_{1}}{b_{1}^{3}}\left(\frac{I c_{1} I c_{2} L^{2}+4 I b_{1} I c_{2} L h_{1}+4 I b_{2} I c_{1} L h_{2}+4 I b_{2} I c_{2} L h_{1}+16 I b_{1} I b_{2} h_{1} h_{2}}{I c_{1} I c_{2} L^{2}+I b_{1} I c_{2} L h_{1}+4 I b_{2} I c_{1} L h_{2}+I b_{2} I c_{2} L h_{1}+4 I b_{1} I b_{2} h_{1} h_{2}}\right) \tag{15} \end{equation*}(15)Keq=6EIc1b13(Ic1Ic2L2+4Ib1Ic2Lh1+4Ib2Ic1Lh2+4Ib2Ic2Lh1+16Ib1Ib2h1h2Ic1Ic2L2+Ib1Ic2Lh1+4Ib2Ic1Lh2+Ib2Ic2Lh1+4Ib1Ib2h1h2)
where K eq K eq K_(eq)\mathrm{K}_{\mathrm{eq}}Keq is the lateral equivalent stiffness, E is the concrete modulus of elasticity, Ic 1 Ic 1 Ic_(1)\mathrm{Ic}_{1}Ic1 and Ic 2 Ic 2 Ic_(2)\mathrm{Ic}_{2}Ic2 are the moments of inertia of the first and second floors columns, respectively, L is the frame span, Ib 1 Ib 1 Ib_(1)\mathrm{Ib}_{1}Ib1 and Ib 2 Ib 2 Ib_(2)\mathrm{Ib}_{2}Ib2 are the moments of inertia of the first and second floors girders, respectively, and h 1 h 1 h_(1)\mathrm{h}_{1}h1 and h 2 h 2 h_(2)\mathrm{h}_{2}h2 are the heights of the first and second floors, respectively. By substituting in equation (15) with Ic 1 = Ic 2 Ic 1 = Ic 2 Ic_(1)=Ic_(2)\mathrm{Ic}_{1}=\mathrm{Ic}_{2}Ic1=Ic2 and Ib 1 = Ib 2 Ib 1 = Ib 2 Ib_(1)=Ib_(2)\mathrm{Ib}_{1}=\mathrm{Ib}_{2}Ib1=Ib2 to obtain a more simplified form as follows:
(16) K e q = ( 6 E I c 1 ( I c 1 2 L 2 + 16 I b 1 2 h 1 h 2 + 8 I b 1 I c 1 L h 1 + 4 I b 1 I c 1 L h 2 ) ) ( h 1 3 ( I c 1 2 L 2 + 4 I b 1 2 h 1 h 2 + 2 I b 1 I c 1 L h 1 + 4 I b 1 I c 1 L h 2 ) ) (16) K e q = 6 E I c 1 I c 1 2 L 2 + 16 I b 1 2 h 1 h 2 + 8 I b 1 I c 1 L h 1 + 4 I b 1 I c 1 L h 2 h 1 3 I c 1 2 L 2 + 4 I b 1 2 h 1 h 2 + 2 I b 1 I c 1 L h 1 + 4 I b 1 I c 1 L h 2 {:(16)K_(eq)=((6EIc_(1)(Ic_(1)^(2)L^(2)+16 Ib_(1)^(2)h_(1)h_(2)+8Ib_(1)Ic_(1)Lh_(1)+4Ib_(1)Ic_(1)Lh_(2))))/((h_(1)^(3)(Ic_(1)^(2)L^(2)+4Ib_(1)^(2)h_(1)h_(2)+2Ib_(1)Ic_(1)Lh_(1)+4Ib_(1)Ic_(1)Lh_(2)))):}\begin{equation*} K_{e q}=\frac{\left(6 E I c_{1}\left(I c_{1}^{2} L^{2}+16 I b_{1}^{2} h_{1} h_{2}+8 I b_{1} I c_{1} L h_{1}+4 I b_{1} I c_{1} L h_{2}\right)\right)}{\left(h_{1}^{3}\left(I c_{1}^{2} L^{2}+4 I b_{1}^{2} h_{1} h_{2}+2 I b_{1} I c_{1} L h_{1}+4 I b_{1} I c_{1} L h_{2}\right)\right)} \tag{16} \end{equation*}(16)Keq=(6EIc1(Ic12L2+16Ib12h1h2+8Ib1Ic1Lh1+4Ib1Ic1Lh2))(h13(Ic12L2+4Ib12h1h2+2Ib1Ic1Lh1+4Ib1Ic1Lh2))
By substituting into equation 16 with I b 1 = α Ic 1 I b 1 = α Ic 1 I_(b1)=alpha^(**)Ic_(1)\mathrm{I}_{\mathrm{b} 1}=\alpha^{*} \mathrm{Ic}_{1}Ib1=αIc1.
(17) K e q = 6 E I b 1 α ( 16 h 1 h 2 + L 2 α 2 + 8 L α h 1 + 4 L α h 2 ) ( h 1 3 ( 4 h 1 h 2 + L 2 α 2 + 2 L α h 1 + 4 L α h 2 ) ) (17) K e q = 6 E I b 1 α 16 h 1 h 2 + L 2 α 2 + 8 L α h 1 + 4 L α h 2 h 1 3 4 h 1 h 2 + L 2 α 2 + 2 L α h 1 + 4 L α h 2 {:(17)K_(eq)=(6EIb_(1)alpha(16h_(1)h_(2)+L^(2)alpha^(2)+8L alphah_(1)+4L alphah_(2)))/((h_(1)^(3)(4h_(1)h_(2)+L^(2)alpha^(2)+2L alphah_(1)+4L alphah_(2)))):}\begin{equation*} K_{e q}=\frac{6 E I b_{1} \alpha\left(16 h_{1} h_{2}+L^{2} \alpha^{2}+8 L \alpha h_{1}+4 L \alpha h_{2}\right)}{\left(h_{1}^{3}\left(4 h_{1} h_{2}+L^{2} \alpha^{2}+2 L \alpha h_{1}+4 L \alpha h_{2}\right)\right)} \tag{17} \end{equation*}(17)Keq=6EIb1α(16h1h2+L2α2+8Lαh1+4Lαh2)(h13(4h1h2+L2α2+2Lαh1+4Lαh2))
where α α alpha\alphaα is defined as the beam-to-column moment of inertia ratio.
This method is applied on a broader scale to include multi-story frames, as shown in Fig. 4. The equivalent lateral stiffness is derived using the same equation, based on the sections and heights of the first and second floors. This indicates that the
upper floors have a negligible impact on determining the lateral stiffness, allowing for a simplified analysis focused primarily on the lower floors. This assumption helps streamline calculations while maintaining reasonable accuracy for structural behavior under lateral loads. It becomes clear that the previous equation for K eq K eq K_(eq)\mathrm{K}_{\mathrm{eq}}Keq, expressed in terms of the properties and behaviors of the first two stories, can be generalized for multi-story structures with high accuracy.
Figure 4: The simplification of a multi-storey frame into an SDOF.

Equivalent seismic mass of the multi-storey buildings

Equivalent seismic mass is a simplified representation of how each floor's weight contributes to the building's overall seismic response. It considers the distribution of mass and how it affects the structure's dynamic behavior during an earthquake. The sequence of the method is as illustrated in Fig. 5, where the following building consists of three floors. Two iterations are performed to obtain an equivalent mass concentrated at the level of the first floor. The first iteration involves shifting the mass of the third floor to the level of the second floor to obtain m eq m eq m_(eq)\mathrm{m}_{\mathrm{eq}}meq. Then, the mass obtained from the previous iteration, m eq 1 m eq  1 m_("eq "1)\mathrm{m}_{\text {eq } 1}meq 1, is shifted to the level of the first floor to calculate m eq 2 m eq  2 m_("eq "2)\mathrm{m}_{\text {eq } 2}meq 2, resulting in a single equivalent mass at the first-floor level.
Figure 5: Procedures for obtaining equivalent seismic mass.
To perform this task, return again to the two-storey frame shown in Fig. 3, but the mass of the second floor is set by m 2 m 2 m_(2)\mathrm{m}_{2}m2, and the mass of the first floor is omitted; then the structure mass matrix takes the form:
(18) M = [ 0 0 0 m 2 ] (18) M = 0 0 0 m 2 {:(18)M=[[0,0],[0,m_(2)]]:}M=\left[\begin{array}{cc} 0 & 0 \tag{18}\\ 0 & m_{2} \end{array}\right](18)M=[000m2]
Again, by substituting into the eigen equation:
(19) | K ω e q 2 2 M | = 0 (19) K ω e q 2 2 M = 0 {:(19)|K-omega_(eq2)^(2)M|=0:}\begin{equation*} \left|K-\omega_{e q 2}^{2} M\right|=0 \tag{19} \end{equation*}(19)|Kωeq22M|=0
Then the equivalent mass is expressed as:
(20) m e q = K e q ω e q 2 2 (20) m e q = K e q ω e q 2 2 {:(20)m_(eq)=(K_(eq))/(omega_(eq2)^(2)):}\begin{equation*} m_{e q}=\frac{K_{e q}}{\omega_{e q 2}^{2}} \tag{20} \end{equation*}(20)meq=Keqωeq22
To get the mass contribution factor, apply the following equation:
(21) f 0 = m e q m 2 = K e q m 2 ω e q 2 2 (21) f 0 = m e q m 2 = K e q m 2 ω e q 2 2 {:(21)f_(0)=(m_(eq))/(m_(2))=(K_(eq))/(m_(2)omega_(eq2)^(2)):}\begin{equation*} f_{0}=\frac{m_{e q}}{m_{2}}=\frac{K_{e q}}{m_{2} \omega_{e q 2}^{2}} \tag{21} \end{equation*}(21)f0=meqm2=Keqm2ωeq22
By substituting into Eqn. (21) for α = I b / I c α = I b / I c alpha=I_(b)//I_(c)\alpha=\mathrm{I}_{\mathrm{b}} / \mathrm{I}_{\mathrm{c}}α=Ib/Ic and β 2 = L / h 2 β 2 = L / h 2 beta_(2)=L//h_(2)\beta_{2}=\mathrm{L} / \mathrm{h}_{2}β2=L/h2, obtaining:
(22) f 0 = ( 4 α 2 h 1 4 + 4 α 2 h 1 h 2 3 + 3 α β 2 h 1 4 + 4 α β 2 h 1 3 h 2 + 6 α β 2 h 1 2 h 2 2 + 8 α β 2 h 1 h 2 3 + β 2 2 h 1 3 h 2 + 3 β 2 2 h 1 2 h 2 2 + 3 β 2 2 h 1 h 2 2 + β 2 2 h 2 4 ) ( h 1 3 ( 4 α 2 h 1 + β 2 2 h 2 + 2 α β 2 h 1 + 4 α β 2 h 2 ) ) (22) f 0 = 4 α 2 h 1 4 + 4 α 2 h 1 h 2 3 + 3 α β 2 h 1 4 + 4 α β 2 h 1 3 h 2 + 6 α β 2 h 1 2 h 2 2 + 8 α β 2 h 1 h 2 3 + β 2 2 h 1 3 h 2 + 3 β 2 2 h 1 2 h 2 2 + 3 β 2 2 h 1 h 2 2 + β 2 2 h 2 4 h 1 3 4 α 2 h 1 + β 2 2 h 2 + 2 α β 2 h 1 + 4 α β 2 h 2 {:(22)f_(0)=((4alpha^(2)h_(1)^(4)+4alpha^(2)h_(1)h_(2)^(3)+3alphabeta_(2)h_(1)^(4)+4alphabeta_(2)h_(1)^(3)h_(2)+6alphabeta_(2)h_(1)^(2)h_(2)^(2)+8alphabeta_(2)h_(1)h_(2)^(3)+beta_(2)^(2)h_(1)^(3)h_(2)+3beta_(2)^(2)h_(1)^(2)h_(2)^(2)+3beta_(2)^(2)h_(1)h_(2)^(2)+beta_(2)^(2)h_(2)^(4)))/((h_(1)^(3)(4alpha^(2)h_(1)+beta_(2)^(2)h_(2)+2alphabeta_(2)h_(1)+4alphabeta_(2)h_(2)))):}\begin{equation*} f_{0}=\frac{\left(4 \alpha^{2} h_{1}^{4}+4 \alpha^{2} h_{1} h_{2}^{3}+3 \alpha \beta_{2} h_{1}^{4}+4 \alpha \beta_{2} h_{1}^{3} h_{2}+6 \alpha \beta_{2} h_{1}^{2} h_{2}^{2}+8 \alpha \beta_{2} h_{1} h_{2}^{3}+\beta_{2}^{2} h_{1}^{3} h_{2}+3 \beta_{2}^{2} h_{1}^{2} h_{2}^{2}+3 \beta_{2}^{2} h_{1} h_{2}^{2}+\beta_{2}^{2} h_{2}^{4}\right)}{\left(h_{1}^{3}\left(4 \alpha^{2} h_{1}+\beta_{2}^{2} h_{2}+2 \alpha \beta_{2} h_{1}+4 \alpha \beta_{2} h_{2}\right)\right)} \tag{22} \end{equation*}(22)f0=(4α2h14+4α2h1h23+3αβ2h14+4αβ2h13h2+6αβ2h12h22+8αβ2h1h23+β22h13h2+3β22h12h22+3β22h1h22+β22h24)(h13(4α2h1+β22h2+2αβ2h1+4αβ2h2))
To simplify calculations, the formula can be modified using h 1 = n h 2 h 1 = n h 2 h_(1)=n*h_(2)\mathrm{h}_{1}=\mathrm{n} \cdot \mathrm{h}_{2}h1=nh2 in Eqn. (22), where h 2 h 2 h_(2)\mathrm{h}_{2}h2 is the floor height, h 1 h 1 h_(1)\mathrm{h}_{1}h1 is the height at which the mass has accumulated from the previous floor, and n n nnn is the floor number to which the mass was moved to simplify the calculation.
(23) f 0 = ( 4 α 2 n 4 + 4 α 2 n + 3 α β 2 n 4 + 4 α β 2 n 3 + 6 α β 2 n 2 + 8 α β 2 n + β 2 2 n 3 + 3 β 2 2 n 2 + 3 β 2 2 n + β 2 2 ) ( n 3 ( 4 α β 2 + 4 α 2 n + β 2 + 2 α β 2 n ) ) (23) f 0 = 4 α 2 n 4 + 4 α 2 n + 3 α β 2 n 4 + 4 α β 2 n 3 + 6 α β 2 n 2 + 8 α β 2 n + β 2 2 n 3 + 3 β 2 2 n 2 + 3 β 2 2 n + β 2 2 n 3 4 α β 2 + 4 α 2 n + β 2 + 2 α β 2 n {:(23)f_(0)=((4alpha^(2)n^(4)+4alpha^(2)n+3alphabeta_(2)n^(4)+4alphabeta_(2)n^(3)+6alphabeta_(2)n^(2)+8alphabeta_(2)n+beta_(2)^(2)n^(3)+3beta_(2)^(2)n^(2)+3beta_(2)^(2)n+beta_(2)^(2)))/((n^(3)(4alphabeta_(2)+4alpha^(2)n+beta^(2)+2alphabeta_(2)n))):}\begin{equation*} f_{0}=\frac{\left(4 \alpha^{2} n^{4}+4 \alpha^{2} n+3 \alpha \beta_{2} n^{4}+4 \alpha \beta_{2} n^{3}+6 \alpha \beta_{2} n^{2}+8 \alpha \beta_{2} n+\beta_{2}^{2} n^{3}+3 \beta_{2}^{2} n^{2}+3 \beta_{2}^{2} n+\beta_{2}^{2}\right)}{\left(\mathrm{n}^{3}\left(4 \alpha \beta_{2}+4 \alpha^{2} n+\beta^{2}+2 \alpha \beta_{2} \mathrm{n}\right)\right)} \tag{23} \end{equation*}(23)f0=(4α2n4+4α2n+3αβ2n4+4αβ2n3+6αβ2n2+8αβ2n+β22n3+3β22n2+3β22n+β22)(n3(4αβ2+4α2n+β2+2αβ2n))
The following equation is obtained by replacing γ = β 2 α γ = β 2 α gamma=(beta_(2))/(alpha)\gamma=\frac{\beta_{2}}{\alpha}γ=β2α, in Eqn. (23):
(24) f 0 = ( 4 n 4 + 4 n + 3 γ n 4 + 4 γ n 3 + 6 γ n 2 + 8 γ n + γ 2 n 3 + 3 γ 2 n 2 + 3 γ 2 n + γ 2 ) ( n 3 ( 4 γ + 4 n + γ 2 + 2 γ n ) (24) f 0 = 4 n 4 + 4 n + 3 γ n 4 + 4 γ n 3 + 6 γ n 2 + 8 γ n + γ 2 n 3 + 3 γ 2 n 2 + 3 γ 2 n + γ 2 n 3 4 γ + 4 n + γ 2 + 2 γ n {:(24)f_(0)=((4n^(4)+4n+3**gamman^(4)+4gamman^(3)+6gamman^(2)+8gamma n+gamma^(2)n^(3)+3gamma^(2)n^(2)+3gamma^(2)n+gamma^(2)))/((n^(3)(4gamma+4n+gamma^(2)+2gamma n)):}\begin{equation*} f_{0}=\frac{\left(4 n^{4}+4 n+3 * \gamma n^{4}+4 \gamma n^{3}+6 \gamma n^{2}+8 \gamma n+\gamma^{2} n^{3}+3 \gamma^{2} n^{2}+3 \gamma^{2} n+\gamma^{2}\right)}{\left(n^{3}\left(4 \gamma+4 n+\gamma^{2}+2 \gamma n\right)\right.} \tag{24} \end{equation*}(24)f0=(4n4+4n+3γn4+4γn3+6γn2+8γn+γ2n3+3γ2n2+3γ2n+γ2)(n3(4γ+4n+γ2+2γn)
A simplified form of Equation (22) can be obtained by setting the stiffness ratio α = 1 α = 1 alpha=1\alpha=1α=1, which corresponds to the case where the beam and column sections are identical. This assumption leads to the following simplified expression:
(25) f 0 = ( β 2 2 h 1 3 h 2 + 3 β 2 2 h 1 2 h 2 2 + 3 β 2 2 h 1 h 2 3 + β 2 2 h 2 4 + 3 β 2 h 1 4 + 4 β 2 h 1 3 h 2 + 6 β 2 h 1 2 h 2 2 + 8 β 2 h 1 h 2 3 + 4 h 1 4 + 4 h 1 h 2 3 ) ( h 1 3 ( 4 h 1 + 2 β 2 h 1 + 4 β 2 h 2 + β 2 2 h 2 ) ) (25) f 0 = β 2 2 h 1 3 h 2 + 3 β 2 2 h 1 2 h 2 2 + 3 β 2 2 h 1 h 2 3 + β 2 2 h 2 4 + 3 β 2 h 1 4 + 4 β 2 h 1 3 h 2 + 6 β 2 h 1 2 h 2 2 + 8 β 2 h 1 h 2 3 + 4 h 1 4 + 4 h 1 h 2 3 h 1 3 4 h 1 + 2 β 2 h 1 + 4 β 2 h 2 + β 2 2 h 2 {:(25)f_(0)=((beta_(2)^(2)h_(1)^(3)h_(2)+3beta_(2)^(2)h_(1)^(2)h_(2)^(2)+3beta_(2)^(2)h_(1)h_(2)^(3)+beta_(2)^(2)h_(2)^(4)+3beta_(2)h_(1)^(4)+4beta_(2)h_(1)^(3)h_(2)+6beta_(2)h_(1)^(2)h_(2)^(2)+8beta_(2)h_(1)h_(2)^(3)+4h_(1)^(4)+4h_(1)h_(2)^(3)))/((h_(1)^(3)(4h_(1)+2beta_(2)h_(1)+4beta_(2)h_(2)+beta_(2)^(2)h_(2)))):}\begin{equation*} f_{0}=\frac{\left(\beta_{2}^{2} h_{1}^{3} h_{2}+3 \beta_{2}^{2} h_{1}^{2} h_{2}^{2}+3 \beta_{2}^{2} h_{1} h_{2}^{3}+\beta_{2}^{2} h_{2}^{4}+3 \beta_{2} h_{1}^{4}+4 \beta_{2} h_{1}^{3} h_{2}+6 \beta_{2} h_{1}^{2} h_{2}^{2}+8 \beta_{2} h_{1} h_{2}^{3}+4 h_{1}^{4}+4 h_{1} h_{2}^{3}\right)}{\left(h_{1}^{3}\left(4 h_{1}+2 \beta_{2} h_{1}+4 \beta_{2} h_{2}+\beta_{2}^{2} h_{2}\right)\right)} \tag{25} \end{equation*}(25)f0=(β22h13h2+3β22h12h22+3β22h1h23+β22h24+3β2h14+4β2h13h2+6β2h12h22+8β2h1h23+4h14+4h1h23)(h13(4h1+2β2h1+4β2h2+β22h2))
By substituting in Eq. (24) with h 1 = n h 2 h 1 = n h 2 h_(1)=n^(**)h_(2)\mathrm{h}_{1}=\mathrm{n}^{*} \mathrm{~h}_{2}h1=n h2
(26) f 0 = ( β 2 2 n 3 + 3 β 2 2 n 2 + 3 β 2 2 n + β 2 2 + 3 β 2 n 4 + 4 β 2 n 3 + 6 β 2 n 2 + 8 β 2 n + 4 n 4 + 4 n ) ( n 3 ( 4 β 2 + 4 n + 2 β 2 n + β 2 2 ) ) (26) f 0 = β 2 2 n 3 + 3 β 2 2 n 2 + 3 β 2 2 n + β 2 2 + 3 β 2 n 4 + 4 β 2 n 3 + 6 β 2 n 2 + 8 β 2 n + 4 n 4 + 4 n n 3 4 β 2 + 4 n + 2 β 2 n + β 2 2 {:(26)f_(0)=((beta_(2)^(2)n^(3)+3beta_(2)^(2)n^(2)+3beta_(2)^(2)n+beta_(2)^(2)+3beta_(2)n^(4)+4beta_(2)n^(3)+6beta_(2)n^(2)+8beta_(2)n+4n^(4)+4n))/((n^(3)(4beta_(2)+4n+2beta_(2)n+beta_(2)^(2)))):}\begin{equation*} f_{0}=\frac{\left(\beta_{2}^{2} n^{3}+3 \beta_{2}^{2} n^{2}+3 \beta_{2}^{2} n+\beta_{2}^{2}+3 \beta_{2} n^{4}+4 \beta_{2} n^{3}+6 \beta_{2} n^{2}+8 \beta_{2} n+4 n^{4}+4 n\right)}{\left(n^{3}\left(4 \beta_{2}+4 n+2 \beta_{2} n+\beta_{2}^{2}\right)\right)} \tag{26} \end{equation*}(26)f0=(β22n3+3β22n2+3β22n+β22+3β2n4+4β2n3+6β2n2+8β2n+4n4+4n)(n3(4β2+4n+2β2n+β22))
The equivalent mass in an SDOF system can be calculated using the following formula:
(27) m e q [ n , ( n 1 ) ] = m ( m 1 ) + f 0 [ n , ( n 1 ) ] (27) m e q [ n , ( n 1 ) ] = m ( m 1 ) + f 0 [ n , ( n 1 ) ] {:(27)m_(eq[n,(n-1)])=m_((m-1))+f_(0[n,(n-1)]):}\begin{equation*} m_{e q[\mathrm{n},(\mathrm{n}-1)]}=m_{(m-1)}+f_{0[\mathrm{n},(\mathrm{n}-1)]} \tag{27} \end{equation*}(27)meq[n,(n1)]=m(m1)+f0[n,(n1)]
where n represents the number of the floor whose mass is to be shifted, and f 0 f 0 f_(0)\mathrm{f}_{0}f0 is the mass-shifting coefficient to the lower floor, which was derived in its original form from Eqn. (22). Meanwhile, m eq m eq m_(eq)\mathrm{m}_{\mathrm{eq}}meq is defined as the total equivalent mass of the frame.
A multi-degree-of-freedom (MDOF) frame structure may have many modes of vibration. Usually, we use the first mode of vibration in orthogonal directions to estimate seismic demand. The fundamental period of the equivalent SDOF system is given by:
(28) T = 2 π m e q K e q (28) T = 2 π m e q K e q {:(28)T=2pisqrt((m_(eq))/(K_(eq))):}\begin{equation*} T=2 \pi \sqrt{\frac{m_{e q}}{K_{e q}}} \tag{28} \end{equation*}(28)T=2πmeqKeq
A step-by-step flowchart of the proposed method for the computation of the fundamental time period of the RC MRFs is provided in Fig. 6.
Figure 6: Proposed procedure flowchart.

Numerical modeling framework

Assumptions and scope

The assumption of this study mainly examines bare frame buildings, excluding the stiffness contribution of infill walls. The building models examined range from 3 to 30 stories in height. For each number of stories studied, different storey heights were chosen, which vary from 2.5 m to 4 m . The suggested formula has mainly been derived through estimating the fundamental time period of RC moment-resisting frames under elastic conditions.

Modeling approach

The frame structure shown was modeled using frame elements in the ETABS environment. The beams and columns were idealized as linear elastic elements. Each member (beams and columns) was discretized into a single frame element between
joints, given the simple geometry, and to maintain computational efficiency. The finite element model (FEM) was used to validate the proposed analytical framework for seismic response prediction.

Element type and meshing strategy

Beam-column members were modeled by [Timoshenko beam elements] to account for both bending and shear deformations. A uniform structured mesh was employed for the model, with local refining in significant regions, particularly in beam-column joints and positions of supports, for better local stress accuracy and realistic distribution of stiffness. Mesh density was selected following initial convergence tests to ensure negligible variation in global response parameters such as the fundamental period and the maximum story drift with additional refinement.

Mesh sensitivity and convergence

A detailed mesh sensitivity analysis was not performed. The initial mesh density was selected by trial, starting with element lengths of approximately 1.0 m , and subsequently refined to about 0.6 m , at which point the variations in fundamental frequency and maximum storey drift became negligible. Local refinement was applied at beam-column connections and support areas, where element lengths were further reduced to the range of 0.20 m to 0.35 m to improve accuracy in regions with high stress concentrations. The final mesh meets the minimum convergence criteria, with global response parameters remaining stable under further refinement.

Mass and boundary condition modeling with diaphragm assumptions

In order to achieve a simulated fixed base condition, boundary conditions were provided in the form of fixed supports on column foundations that restrained all the translational and rotational degrees of freedom. By removing the influence of soil-structure interaction, this assumption makes dynamic analysis easier and reflects idealized conditions typical in initial design and verification studies. To provide equal lateral displacement in each story, rigid diaphragms were also provided at every level to impose in-plane stiffness. For dynamic analysis, a 5 % 5 % 5%5 \%5% damping ratio was applied overall to account for inherent structural and material damping, but no other damping elements were specifically simulated.

Visualization of FEM setup

For enhanced transparency and reproducibility, schematic illustration of the finite element model has been provided (Fig.7), showing mesh distribution, boundary conditions, and areas of mesh refinement. Boundary conditions are indicated to mark fixed and restrained DOF at support nodes. All joints were assigned two degrees of freedom (DOF) per node, namely horizontal translation ( u u uuu ) and rotation ( θ θ theta\thetaθ ). The figure also indicates local mesh refinement where regions of interest such as beam-column joints and in the regions near the zones of support are higher in stress gradient. Such graphical inspection ensures that the model's constraints, degrees of freedom, and discretization scheme accurately portray the intended real-life phenomenon.
Figure 7: Finite element meshing strategy

Loading conditions

Dead load is the weight of the structural and non-structural elements themselves, i.e., beams, columns, and finishes. Dead load was automatically computed from cross-section dimensions and assigned material densities ( 25 kN / m 3 25 kN / m 3 25kN//m^(3)25 \mathrm{kN} / \mathrm{m}^{3}25kN/m3 for reinforced concrete). Live loads are temporary occupancy-related forces. A uniform live load was placed on every beam in accordance
with code requirements, using a value of 2.0 kN / m 2 2.0 kN / m 2 2.0kN//m^(2)2.0 \mathrm{kN} / \mathrm{m}^{2}2.0kN/m2. For seismic combination, a smaller portion of the live load (typically 25 % 25 % 25%25 \%25% or according to code) was included with the mass for dynamic analysis.
The seismic lateral forces were applied in the present study through the response spectrum analysis (RSA) method. A codebased design response spectrum was applied, taking into account parameters that include site class, seismic zone, importance factor, and structural damping ratio (assumed to be 5 % 5 % 5%5 \%5% for reinforced concrete frames).

Validation of the proposed equation

Overview of validation methodology

In this study, a number of examples are presented to assess the validity of the suggested equation to estimate buildings' fundamental period. Initially, two cases of moment-resisting frame (MRF) structures of reinforced concrete with uniform cross sections are investigated. Then, applying more models with varying section properties along their height. The fundamental periods derived from the suggested equation are compared with numerical model eigenvalue analysis solutions, which represent the benchmark solution. Comparisons are made using ETABS 2018 and MATLAB 2023 [38, 39] for reliability and accuracy.

Estimation of the fundamental period for a six-story reinforced concrete frame with uniform member sections

The following frame consists of six levels, each with a height of 3 meters and a span of 6 meters. The sections of both beams and columns are uniform throughout the building's entire height, as shown in Fig. 8 and Tab.1. This example was solved using the proposed method by first determining the equivalent lateral stiffness and then computing the mass transfer coefficient from one floor to the level of the next. This method is carried out in five iterations. The result of applying the proposed method is compared with the values from FEM as shown in Tab. 2.
Figure 8: Configuration of reinforced concrete MRF.
Floor Column Beam ( α ) = I b / I c ( α ) = I b / I c (alpha)=I_(b)//I_(c)(\alpha)=\mathrm{I}_{\mathrm{b}} / \mathrm{I}_{\mathrm{c}}(α)=Ib/Ic ( β 2 ) = L / h 2 β 2 = L / h 2 (beta_(2))=L//h_(2)\left(\beta_{2}\right)=\mathrm{L} / \mathrm{h}_{2}(β2)=L/h2
Section
( m ) ( m ) (m)(\mathrm{m})(m)
Section (m)| Section | | :---: | | $(\mathrm{m})$ |
I C I C I_(C)\mathrm{I}_{\mathrm{C}}IC
( m 4 ) m 4 (m^(4))\left(\mathrm{m}^{4}\right)(m4)
I_(C) (m^(4))| $\mathrm{I}_{\mathrm{C}}$ | | :---: | | $\left(\mathrm{m}^{4}\right)$ |
h 2 h 2 h_(2)\mathrm{h}_{2}h2
( m ) ( m ) (m)(\mathrm{~m})( m)
h_(2) (m)| $\mathrm{h}_{2}$ | | :---: | | $(\mathrm{~m})$ |
Section
( m ) ( m ) (m)(\mathrm{m})(m)
Section (m)| Section | | :---: | | $(\mathrm{m})$ |
I b I b I_(b)\mathrm{I}_{\mathrm{b}}Ib
( m 4 ) m 4 (m^(4))\left(\mathrm{m}^{4}\right)(m4)
I_(b) (m^(4))| $\mathrm{I}_{\mathrm{b}}$ | | :---: | | $\left(\mathrm{m}^{4}\right)$ |
L
( m ) ( m ) (m)(\mathrm{m})(m)
L (m)| L | | :---: | | $(\mathrm{m})$ |
1 st 1 st  1^("st ")1^{\text {st }}1st  to 6 th 6 th  6^("th ")6^{\text {th }}6th 
0.7 × 0.7 × 0.7 xx0.7 \times0.7×
0.7
0.7 xx 0.7| $0.7 \times$ | | :---: | | 0.7 |
0.02 3 0.7 × 0.7 0.7 × 0.7 0.7 xx0.70.7 \times 0.70.7×0.7 0.02 6 1.0 2
Floor Column Beam (alpha)=I_(b)//I_(c) (beta_(2))=L//h_(2) "Section (m)" "I_(C) (m^(4))" "h_(2) (m)" "Section (m)" "I_(b) (m^(4))" "L (m)" 1^("st ") to 6^("th ") "0.7 xx 0.7" 0.02 3 0.7 xx0.7 0.02 6 1.0 2| Floor | Column | | | Beam | | | $(\alpha)=\mathrm{I}_{\mathrm{b}} / \mathrm{I}_{\mathrm{c}}$ | $\left(\beta_{2}\right)=\mathrm{L} / \mathrm{h}_{2}$ | | :---: | :---: | :---: | :---: | :---: | :---: | :---: | :---: | :---: | | | Section <br> $(\mathrm{m})$ | $\mathrm{I}_{\mathrm{C}}$ <br> $\left(\mathrm{m}^{4}\right)$ | $\mathrm{h}_{2}$ <br> $(\mathrm{~m})$ | Section <br> $(\mathrm{m})$ | $\mathrm{I}_{\mathrm{b}}$ <br> $\left(\mathrm{m}^{4}\right)$ | L <br> $(\mathrm{m})$ | | | | $1^{\text {st }}$ to $6^{\text {th }}$ | $0.7 \times$ <br> 0.7 | 0.02 | 3 | $0.7 \times 0.7$ | 0.02 | 6 | 1.0 | 2 |
Table 1. Geometric and mechanical properties of frame elements adopted for analytical validation.
m eq m eq  m_("eq ")\mathrm{m}_{\text {eq }}meq 
( KN . s 2 / m ) KN . s 2 / m (KN.s^(2)//m)\left(\mathrm{KN} . \mathrm{s}^{2} / \mathrm{m}\right)(KN.s2/m)
m_("eq ") (KN.s^(2)//m)| $\mathrm{m}_{\text {eq }}$ | | :---: | | $\left(\mathrm{KN} . \mathrm{s}^{2} / \mathrm{m}\right)$ |
k eq k eq  k_("eq ")\mathrm{k}_{\text {eq }}keq 
( KN / m ) ( KN / m ) (KN//m)(\mathrm{KN} / \mathrm{m})(KN/m)
k_("eq ") (KN//m)| $\mathrm{k}_{\text {eq }}$ | | :---: | | $(\mathrm{KN} / \mathrm{m})$ |
ω ω omega\omegaω
( rad / sec ) ( rad / sec ) (rad//sec)(\mathrm{rad} / \mathrm{sec})(rad/sec)
omega (rad//sec)| $\omega$ | | :---: | | $(\mathrm{rad} / \mathrm{sec})$ |
T FEM T FEM  T_("FEM ")\mathrm{T}_{\text {FEM }}TFEM 
( sec ) ( sec ) (sec)(\mathrm{sec})(sec)
T_("FEM ") (sec)| $\mathrm{T}_{\text {FEM }}$ | | :---: | | $(\mathrm{sec})$ |
T calculated T calculated  T_("calculated ")\mathrm{T}_{\text {calculated }}Tcalculated 
( sec ) ( sec ) (sec)(\mathrm{sec})(sec)
( Eq .28 ) ( Eq .28 ) (Eq.28)(\mathrm{Eq} .28)(Eq.28)
T_("calculated ") (sec) (Eq.28)| $\mathrm{T}_{\text {calculated }}$ | | :---: | | $(\mathrm{sec})$ | | $(\mathrm{Eq} .28)$ |
Error%
709.8 243810 18.7 0.32 0.33 4 % 4 % 4%4 \%4%
"m_("eq ") (KN.s^(2)//m)" "k_("eq ") (KN//m)" "omega (rad//sec)" "T_("FEM ") (sec)" "T_("calculated ") (sec) (Eq.28)" Error% 709.8 243810 18.7 0.32 0.33 4%| $\mathrm{m}_{\text {eq }}$ <br> $\left(\mathrm{KN} . \mathrm{s}^{2} / \mathrm{m}\right)$ | $\mathrm{k}_{\text {eq }}$ <br> $(\mathrm{KN} / \mathrm{m})$ | $\omega$ <br> $(\mathrm{rad} / \mathrm{sec})$ | $\mathrm{T}_{\text {FEM }}$ <br> $(\mathrm{sec})$ | $\mathrm{T}_{\text {calculated }}$ <br> $(\mathrm{sec})$ <br> $(\mathrm{Eq} .28)$ | Error% | | :---: | :---: | :---: | :---: | :---: | :---: | | 709.8 | 243810 | 18.7 | 0.32 | 0.33 | $4 \%$ |
Table 2: The comparison between the fundamental period by the proposed equation with that calculated by FEM.
By comparing the results, it was found that the value of the fundamental period obtained from the suggested equation has an error margin of 4 % 4 % 4%4 \%4% compared to the results from the FEM analysis. This indicates the validity of the derived conclusions and the equations used.

Estimation of the fundamental period for a ten-story reinforced concrete frame with uniform member sections

In this example, the comparison is conducted on a 10-storey frame, where the floor height for the first five stories is 3.5 meters, and for the remaining five stories, it is 3 meters as shown in Fig.9. The geometric dimensions of the element of the structure are given in Tab. 3. The simplifying process involved nine iterations to convert the building into an equivalent single-story structure. As shown in Tab. 4, the error margin using the equation is approximately 3 % 3 % 3%3 \%3%, confirming the accuracy of the proposed equation. It can be stated that the proposed formula produces results with high accuracy.
Figure 9: Configuration of RC MRF building.
Floor Column Beam ( α ) = I b / I c ( α ) = I b / I c (alpha)=I_(b)//I_(c)(\alpha)=\mathrm{I}_{\mathrm{b}} / \mathrm{I}_{\mathrm{c}}(α)=Ib/Ic ( β 2 ) = L / h 2 β 2 = L / h 2 (beta_(2))=L//h_(2)\left(\beta_{2}\right)=\mathrm{L} / \mathrm{h}_{2}(β2)=L/h2
1 st 1 st  1^("st ")1^{\text {st }}1st  to 5 th 5 th  5^("th ")5^{\text {th }}5th  1 × 1 1 × 1 1xx11 \times 11×1 0.0833 3.5 1 × 1 1 × 1 1xx11 \times 11×1 0.0833 1.0 1.71
6 th 6 th  6^("th ")6{ }^{\text {th }}6th  to 10 th 10 th  10^("th ")10{ }^{\text {th }}10th  1 × 1 1 × 1 1xx11 \times 11×1 0.0833 3 1 × 1 1 × 1 1xx11 \times 11×1 0.0833 6 1.0 2
Floor Column Beam (alpha)=I_(b)//I_(c) (beta_(2))=L//h_(2) 1^("st ") to 5^("th ") 1xx1 0.0833 3.5 1xx1 0.0833 1.0 1.71 6^("th ") to 10^("th ") 1xx1 0.0833 3 1xx1 0.0833 6 1.0 2| Floor | Column | | | Beam | | | $(\alpha)=\mathrm{I}_{\mathrm{b}} / \mathrm{I}_{\mathrm{c}}$ | $\left(\beta_{2}\right)=\mathrm{L} / \mathrm{h}_{2}$ | | :--- | :--- | :--- | :--- | :--- | :--- | :--- | :--- | :--- | | $1^{\text {st }}$ to $5^{\text {th }}$ | $1 \times 1$ | 0.0833 | 3.5 | $1 \times 1$ | 0.0833 | | 1.0 | 1.71 | | $6{ }^{\text {th }}$ to $10{ }^{\text {th }}$ | $1 \times 1$ | 0.0833 | 3 | $1 \times 1$ | 0.0833 | 6 | 1.0 | 2 |
Table 3: Geometric and Section Properties of Columns and Beams for frame.
m eq m eq  m_("eq ")\mathrm{m}_{\text {eq }}meq 
( KN . s 2 / m ) KN . s 2 / m (KN.s^(2)//m)\left(\mathrm{KN} . \mathrm{s}^{2} / \mathrm{m}\right)(KN.s2/m)
m_("eq ") (KN.s^(2)//m)| $\mathrm{m}_{\text {eq }}$ | | :---: | | $\left(\mathrm{KN} . \mathrm{s}^{2} / \mathrm{m}\right)$ |
k eq k eq  k_("eq ")\mathrm{k}_{\text {eq }}keq 
( KN / m ) ( KN / m ) (KN//m)(\mathrm{KN} / \mathrm{m})(KN/m)
k_("eq ") (KN//m)| $\mathrm{k}_{\text {eq }}$ | | :---: | | $(\mathrm{KN} / \mathrm{m})$ |
ω ω omega\omegaω
( rad / sec ) ( rad / sec ) (rad//sec)(\mathrm{rad} / \mathrm{sec})(rad/sec)
omega (rad//sec)| $\omega$ | | :---: | | $(\mathrm{rad} / \mathrm{sec})$ |
T FEM T FEM  T_("FEM ")\mathrm{T}_{\text {FEM }}TFEM 
( sec ) ( sec ) (sec)(\mathrm{sec})(sec)
T_("FEM ") (sec)| $\mathrm{T}_{\text {FEM }}$ | | :---: | | $(\mathrm{sec})$ |
T calculated T calculated  T_("calculated ")\mathrm{T}_{\text {calculated }}Tcalculated 
( sec ) ( sec ) (sec)(\mathrm{sec})(sec)
( ( ((( Eq.28 ) ) )))
T_("calculated ") (sec) ( Eq.28 )| $\mathrm{T}_{\text {calculated }}$ | | :---: | | $(\mathrm{sec})$ | | $($ Eq.28 $)$ |
Error%
4019.4 596090 12.178 0.52 0.54 3 % 3 % 3%3 \%3%
"m_("eq ") (KN.s^(2)//m)" "k_("eq ") (KN//m)" "omega (rad//sec)" "T_("FEM ") (sec)" "T_("calculated ") (sec) ( Eq.28 )" Error% 4019.4 596090 12.178 0.52 0.54 3%| $\mathrm{m}_{\text {eq }}$ <br> $\left(\mathrm{KN} . \mathrm{s}^{2} / \mathrm{m}\right)$ | $\mathrm{k}_{\text {eq }}$ <br> $(\mathrm{KN} / \mathrm{m})$ | $\omega$ <br> $(\mathrm{rad} / \mathrm{sec})$ | $\mathrm{T}_{\text {FEM }}$ <br> $(\mathrm{sec})$ | $\mathrm{T}_{\text {calculated }}$ <br> $(\mathrm{sec})$ <br> $($ Eq.28 $)$ | Error% | | :---: | :---: | :---: | :---: | :---: | :---: | | 4019.4 | 596090 | 12.178 | 0.52 | 0.54 | $3 \%$ |
Table 4: Comparison between the fundamental period by the proposed equation with FEM.

Simplified equivalent mass calculations with validation examples

The equivalent mass calculation involves performing several iterations to reach the final equivalent mass for the SDOF (Single Degree of Freedom) system. This process represents a burden for the user and a waste of time. The proposed method has been simplified by avoiding multiple iterations for determining the equivalent mass. Instead, an equivalent curve has been derived that allows for a single-step calculation of the equivalent mass in terms of the storey numbers. This ultimately led to the curve shown in Fig. 10. This curve is used to determine the equivalent mass coefficient of the frame ( f 0 f 0 f_(0)\mathrm{f}_{0}f0 ) and convert the system into an SDOF (Single Degree of Freedom) based on the number of floors of the frame, as demonstrated in the following examples.
Tab. 5 presents the results of five models of frames with different story numbers, heights, and sections. The natural period of all the models was calculated one by one from the proposed equation and compared to the value taken from the ETABS software. Upon comparing the difference between the two values, we observe that the proposed equation has an error margin of no more than 6 % 6 % 6%6 \%6%, indicating its accuracy.
Figure 10: The cumulative seismic mass factor of the RC MRF system.

Numerical validation of the proposed period formula for frames with non-uniform crossSECTIONAL PROPERTIES

To validate the proposed Equation (28) for estimating the fundamental period of frame buildings with non-uniform cross-sectional properties along the height, a seven-storey frame was analyzed. Each storey had a height of 3 meters, and the span length was 6 meters, as shown in Figure 11.
A comparison of the equation's results is made using five different configurations in terms of beam-to-column crosssections. Each model has a different beam cross-section compared to the column cross-section; however, this difference remains consistent throughout the building's height, as shown in Tab. 6.
The analysis of models in Tab. 7 demonstrates that a reduction in the alpha coefficient value, which represents the beam-tocolumn inertia ratio, results in a greater difference between the value derived from the proposed equation and the value calculated by FEM. When the alpha value was 0.14 , the difference between the two values reached 10 % 10 % 10%10 \%10%. Therefore, a correction was made to the seismic mass coefficient to obtain the corrected ratio ( f / f 0 f / f 0 f//f_(0)f / f_{0}f/f0 ) as a function of alpha to achieve more accuracy. Fig. 12 shows the corrected seismic mass coefficient as a function of the alpha ratio used. To simplify the steps and obtain results quickly and with high accuracy, an equation was derived as a function of alpha to facilitate access to results faster and with higher accuracy, as stated in Eqn. (29).
Model number 1 2 3 4 5
Number of storeys 10 15 20 25 30
I b ( m ) 4 I b ( m ) 4 I_(b(m))^(4)\mathrm{I}_{\mathrm{b}(\mathrm{m})}{ }^{4}Ib(m)4 0.0341 0.0833 0.02 1.33 1.33
I c ( m ) I c ( m ) I_(c(m))\mathrm{I}_{\mathrm{c}(\mathrm{m})}Ic(m) 0.0341 0.0833 0.02 1.33 1.33
L ( m ) L ( m ) L(m)\mathrm{L}(\mathrm{m})L(m) 8.0 7.0 6.0 8.0 8.0
m eq ( KN . s 2 / m ) m eq  KN . s 2 / m m_("eq ")(KN.s^(2)//m)\mathrm{m}_{\text {eq }}\left(\mathrm{KN} . \mathrm{s}^{2} / \mathrm{m}\right)meq (KN.s2/m) 3420 16140 32701.5 526309 1309360
E ( KN / m 2 ) E KN / m 2 E(KN//m^(2))\mathrm{E}\left(\mathrm{KN} / \mathrm{m}^{2}\right)E(KN/m2) 24099790 24099790 24099790 24099790 24099790
K eq ( KN / m ) K eq  ( KN / m ) K_("eq ")(KN//m)\mathrm{K}_{\text {eq }}(\mathrm{KN} / \mathrm{m})Keq (KN/m) (Eq.16) 371262 937408 684556 12107600 14946270
ω ( rad / s ) ω ( rad / s ) omega(rad//s)\omega(\mathrm{rad} / \mathrm{s})ω(rad/s) 10.417 7.62 4.57 5.24 3.37
T ( sec ) T ( sec ) T(sec)\mathrm{T}(\mathrm{sec})T(sec) (Eqn.28) 0.6 0.82 1.373 1.2 1.85
T Fem ( sec ) T Fem  ( sec ) T_("Fem ")(sec)\mathrm{T}_{\text {Fem }}(\mathrm{sec})TFem (sec) 0.57 0.78 1.39 1.3 1.95
Error (%) 5% 5% 2% 6% 5%
Model number 1 2 3 4 5 Number of storeys 10 15 20 25 30 I_(b(m))^(4) 0.0341 0.0833 0.02 1.33 1.33 I_(c(m)) 0.0341 0.0833 0.02 1.33 1.33 L(m) 8.0 7.0 6.0 8.0 8.0 m_("eq ")(KN.s^(2)//m) 3420 16140 32701.5 526309 1309360 E(KN//m^(2)) 24099790 24099790 24099790 24099790 24099790 K_("eq ")(KN//m) (Eq.16) 371262 937408 684556 12107600 14946270 omega(rad//s) 10.417 7.62 4.57 5.24 3.37 T(sec) (Eqn.28) 0.6 0.82 1.373 1.2 1.85 T_("Fem ")(sec) 0.57 0.78 1.39 1.3 1.95 Error (%) 5% 5% 2% 6% 5%| Model number | 1 | 2 | 3 | 4 | 5 | | :--- | :--- | :--- | :--- | :--- | :--- | | Number of storeys | 10 | 15 | 20 | 25 | 30 | | $\mathrm{I}_{\mathrm{b}(\mathrm{m})}{ }^{4}$ | 0.0341 | 0.0833 | 0.02 | 1.33 | 1.33 | | $\mathrm{I}_{\mathrm{c}(\mathrm{m})}$ | 0.0341 | 0.0833 | 0.02 | 1.33 | 1.33 | | $\mathrm{L}(\mathrm{m})$ | 8.0 | 7.0 | 6.0 | 8.0 | 8.0 | | $\mathrm{m}_{\text {eq }}\left(\mathrm{KN} . \mathrm{s}^{2} / \mathrm{m}\right)$ | 3420 | 16140 | 32701.5 | 526309 | 1309360 | | $\mathrm{E}\left(\mathrm{KN} / \mathrm{m}^{2}\right)$ | 24099790 | 24099790 | 24099790 | 24099790 | 24099790 | | $\mathrm{K}_{\text {eq }}(\mathrm{KN} / \mathrm{m})$ (Eq.16) | 371262 | 937408 | 684556 | 12107600 | 14946270 | | $\omega(\mathrm{rad} / \mathrm{s})$ | 10.417 | 7.62 | 4.57 | 5.24 | 3.37 | | $\mathrm{T}(\mathrm{sec})$ (Eqn.28) | 0.6 | 0.82 | 1.373 | 1.2 | 1.85 | | $\mathrm{T}_{\text {Fem }}(\mathrm{sec})$ | 0.57 | 0.78 | 1.39 | 1.3 | 1.95 | | Error (%) | 5% | 5% | 2% | 6% | 5% |
Table 5: Accuracy assessment of the proposed period estimation formula for frames with varying numbers of storeys.
Figure 11: Configuration of MRF concrete structures.
Model number Column cross section(m) Beam cross section(m) Alpha ( I b / I c I b / I c I_(b)//I_(c)\mathrm{I}_{\mathrm{b}} / \mathrm{I}_{\mathrm{c}}Ib/Ic ) E(KN/m)
Model 1 0.8 × 0.8 0.8 × 0.8 0.8 xx0.80.8 \times 0.80.8×0.8 0.4096
Model 2 0.7 × 0.7 0.7 × 0.7 0.7 xx0.70.7 \times 0.70.7×0.7 0.24
Model 3 1 × 1 1 × 1 1xx11 \times 11×1 0.68 × 0.68 0.68 × 0.68 0.68 xx0.680.68 \times 0.680.68×0.68 0.178 24099790
Model 4 0.65 × 0.65 0.65 × 0.65 0.65 xx0.650.65 \times 0.650.65×0.65 0.21
Model 5 0.62 × 0.62 0.62 × 0.62 0.62 xx0.620.62 \times 0.620.62×0.62 0.14
Model number Column cross section(m) Beam cross section(m) Alpha ( I_(b)//I_(c) ) E(KN/m) Model 1 0.8 xx0.8 0.4096 Model 2 0.7 xx0.7 0.24 Model 3 1xx1 0.68 xx0.68 0.178 24099790 Model 4 0.65 xx0.65 0.21 Model 5 0.62 xx0.62 0.14 | Model number | Column cross section(m) | Beam cross section(m) | Alpha ( $\mathrm{I}_{\mathrm{b}} / \mathrm{I}_{\mathrm{c}}$ ) | E(KN/m) | | :--- | :--- | :--- | :--- | :--- | | Model 1 | | $0.8 \times 0.8$ | 0.4096 | | | Model 2 | | $0.7 \times 0.7$ | 0.24 | | | Model 3 | $1 \times 1$ | $0.68 \times 0.68$ | 0.178 | 24099790 | | Model 4 | | $0.65 \times 0.65$ | 0.21 | | | Model 5 | | $0.62 \times 0.62$ | 0.14 | |
Table 6: Geometric properties and flexural stiffness ratio ( I b / I c ) I b / I c (I_(b)//I_(c))\left(\mathrm{I}_{\mathrm{b}} / \mathrm{I}_{\mathrm{c}}\right)(Ib/Ic) for five analytical models with varying beam cross sections and constant column dimensions.
Column crosssection(m) Beam crosssection(m) α α alpha\alphaα ( I b / I c I b / I c I_(b)//I_(c)\mathrm{I}_{\mathrm{b}} / \mathrm{I}_{\mathrm{c}}Ib/Ic ) m eq m eq m_(eq)\mathrm{m}_{\mathrm{eq}}meq (KN.sec²/m) k eq k eq k_(eq)\mathrm{k}_{\mathrm{eq}}keq (KN/m) (Eq.15) ω ω omega\omegaω (rad/sec) T proposed equation (sec) T Fem T Fem  T_("Fem ")\mathrm{T}_{\text {Fem }}TFem  (sec) Error%
Model 1 1 × 1 1 × 1 1xx11 \times 11×1 0.8 × 0.8 0.8 × 0.8 0.8 xx0.80.8 \times 0.80.8×0.8 0.4096 35604 778590.5 4.67 1.34 1.282 5%
Model 2 1 × 1 1 × 1 1xx11 \times 11×1 0.7 × 0.7 0.7 × 0.7 0.7 xx0.70.7 \times 0.70.7×0.7 0.24 32701.5 684556.8 4.57 1.37 1.399 2%
Model 3 1 × 1 1 × 1 1xx11 \times 11×1 0.65 × 0.65 0.65 × 0.65 0.65 xx0.650.65 \times 0.650.65×0.65 0.178 31393.44 640356.4 4.51 1.39 1.492 7%
Model 4 1 × 1 1 × 1 1xx11 \times 11×1 0.68 × 0.68 0.68 × 0.68 0.68 xx0.680.68 \times 0.680.68×0.68 0.21 32167.44 666569.1 4.55 1.38 1.433 4%
Model 5 1 × 1 1 × 1 1xx11 \times 11×1 0.62 × 0.62 0.62 × 0.62 0.62 xx0.620.62 \times 0.620.62×0.62 0.14 30650.4 615260.5 4.48 1.40 1.56 10%
Column crosssection(m) Beam crosssection(m) alpha ( I_(b)//I_(c) ) m_(eq) (KN.sec²/m) k_(eq) (KN/m) (Eq.15) omega (rad/sec) T proposed equation (sec) T_("Fem ") (sec) Error% Model 1 1xx1 0.8 xx0.8 0.4096 35604 778590.5 4.67 1.34 1.282 5% Model 2 1xx1 0.7 xx0.7 0.24 32701.5 684556.8 4.57 1.37 1.399 2% Model 3 1xx1 0.65 xx0.65 0.178 31393.44 640356.4 4.51 1.39 1.492 7% Model 4 1xx1 0.68 xx0.68 0.21 32167.44 666569.1 4.55 1.38 1.433 4% Model 5 1xx1 0.62 xx0.62 0.14 30650.4 615260.5 4.48 1.40 1.56 10%| | Column crosssection(m) | Beam crosssection(m) | $\alpha$ ( $\mathrm{I}_{\mathrm{b}} / \mathrm{I}_{\mathrm{c}}$ ) | $\mathrm{m}_{\mathrm{eq}}$ (KN.sec²/m) | $\mathrm{k}_{\mathrm{eq}}$ (KN/m) (Eq.15) | $\omega$ (rad/sec) | T proposed equation (sec) | $\mathrm{T}_{\text {Fem }}$ (sec) | Error% | | :--- | :--- | :--- | :--- | :--- | :--- | :--- | :--- | :--- | :--- | | Model 1 | $1 \times 1$ | $0.8 \times 0.8$ | 0.4096 | 35604 | 778590.5 | 4.67 | 1.34 | 1.282 | 5% | | Model 2 | $1 \times 1$ | $0.7 \times 0.7$ | 0.24 | 32701.5 | 684556.8 | 4.57 | 1.37 | 1.399 | 2% | | Model 3 | $1 \times 1$ | $0.65 \times 0.65$ | 0.178 | 31393.44 | 640356.4 | 4.51 | 1.39 | 1.492 | 7% | | Model 4 | $1 \times 1$ | $0.68 \times 0.68$ | 0.21 | 32167.44 | 666569.1 | 4.55 | 1.38 | 1.433 | 4% | | Model 5 | $1 \times 1$ | $0.62 \times 0.62$ | 0.14 | 30650.4 | 615260.5 | 4.48 | 1.40 | 1.56 | 10% |
Table 7: Comparative evaluation of fundamental time periods derived from the proposed analytical equation and FEM results for five structural models with varying beam to column flexural stiffness.
Figure 12: Correction factor of a cumulative seismic mass factor.
(29) f 0 f = 2.3679 α 0.549 (29) f 0 f = 2.3679 α 0.549 {:(29)(f_(0))/(f)=2.3679alpha^(0.549):}\begin{equation*} \frac{f_{0}}{f}=2.3679 \alpha^{0.549} \tag{29} \end{equation*}(29)f0f=2.3679α0.549
where α α alpha\alphaα is the beam-to-column inertia ratio ( I b / I c ) , f 0 I b / I c , f 0 (I_(b)//I_(c)),f_(0)\left(\mathrm{I}_{\mathrm{b}} / \mathrm{I}_{\mathrm{c}}\right), \mathrm{f}_{0}(Ib/Ic),f0 is the cumulative factor for shifting the story masses when α α alpha\alphaα equals 1 , and f is the cumulative mass factor adjusted after correction for differences between the cross-sections of beams and columns. The comparison was conducted following the correction of the seismic mass coefficient since the alpha value was less than 1.0 for five different models, as indicated in Tab. 8. These comparisons were executed as shown in Tab. 9 using the revised f value. Upon doing the evaluation, favorable results were achieved in relation to the outputs of the FEM.
Model number Column crosssection (cm) Beam cross-section (cm) α α alpha\alphaα ( I b / I c I b / I c I_(b)//I_(c)\mathrm{I}_{\mathrm{b}} / \mathrm{I}_{\mathrm{c}}Ib/Ic ) M n M n M_(n)\mathrm{M}_{\mathrm{n}}Mn (KN.sec²/m)
Model 1 100 × 100 100 × 100 100 xx100100 \times 100100×100 80 × 80 80 × 80 80 xx8080 \times 8080×80 0.4096 23
Model 2 100 × 100 100 × 100 100 xx100100 \times 100100×100 70 × 70 70 × 70 70 xx7070 \times 7070×70 0.24 21.125
Model 3 100 × 100 100 × 100 100 xx100100 \times 100100×100 65 × 65 65 × 65 65 xx6565 \times 6565×65 0.178 20.28
Model 4 100 × 100 100 × 100 100 xx100100 \times 100100×100 68 × 68 68 × 68 68 xx6868 \times 6868×68 0.21 20.78
Model 5 100 × 100 100 × 100 100 xx100100 \times 100100×100 62 × 62 62 × 62 62 xx6262 \times 6262×62 0.14 19.8
Model number Column crosssection (cm) Beam cross-section (cm) alpha ( I_(b)//I_(c) ) M_(n) (KN.sec²/m) Model 1 100 xx100 80 xx80 0.4096 23 Model 2 100 xx100 70 xx70 0.24 21.125 Model 3 100 xx100 65 xx65 0.178 20.28 Model 4 100 xx100 68 xx68 0.21 20.78 Model 5 100 xx100 62 xx62 0.14 19.8| Model number | Column crosssection (cm) | Beam cross-section (cm) | $\alpha$ ( $\mathrm{I}_{\mathrm{b}} / \mathrm{I}_{\mathrm{c}}$ ) | $\mathrm{M}_{\mathrm{n}}$ (KN.sec²/m) | | :--- | :--- | :--- | :--- | :--- | | Model 1 | $100 \times 100$ | $80 \times 80$ | 0.4096 | 23 | | Model 2 | $100 \times 100$ | $70 \times 70$ | 0.24 | 21.125 | | Model 3 | $100 \times 100$ | $65 \times 65$ | 0.178 | 20.28 | | Model 4 | $100 \times 100$ | $68 \times 68$ | 0.21 | 20.78 | | Model 5 | $100 \times 100$ | $62 \times 62$ | 0.14 | 19.8 |
Table 8: Geometric properties and mass values (Mn) for structural models with varying beam-to-column flexural stiffness ratios.
Model number f 0 f 0 f_(0)\mathrm{f}_{0}f0 (Fig.9) f (Fig.11) m eq m eq m_(eq)\mathrm{m}_{\mathrm{eq}}meq (KN.sec 2 / m 2 / m ^(2)//m{ }^{2} / \mathrm{m}2/m ) k eq k eq  k_("eq ")\mathrm{k}_{\text {eq }}keq  (KN/m) (Eq.15) ω ω omega\omegaω (rad/sec) T proposed T proposed  T_("proposed ")\mathrm{T}_{\text {proposed }}Tproposed  equation(sec) (Eq.28) T Fem T Fem  T_("Fem ")\mathrm{T}_{\text {Fem }}TFem  (sec) Error (%)
Model 1 1548 1583.69 35604 778590.5 4.67 1.34 1.282 5%
Model 2 1548 1880.789 32701.5 684556.8 4.52 1.389 1.399 1%
Model 3 1548 1693.86 31393.44 640356.4 4.09 1.53 1.492 3%
Model 4 1548 2098.70 32167.44 666569.1 4.35 1.44 1.433 1%
Model 5 1548 2098.70 30650.4 615260.5 3.84 1.63 1.56 5%
Model number f_(0) (Fig.9) f (Fig.11) m_(eq) (KN.sec ^(2)//m ) k_("eq ") (KN/m) (Eq.15) omega (rad/sec) T_("proposed ") equation(sec) (Eq.28) T_("Fem ") (sec) Error (%) Model 1 1548 1583.69 35604 778590.5 4.67 1.34 1.282 5% Model 2 1548 1880.789 32701.5 684556.8 4.52 1.389 1.399 1% Model 3 1548 1693.86 31393.44 640356.4 4.09 1.53 1.492 3% Model 4 1548 2098.70 32167.44 666569.1 4.35 1.44 1.433 1% Model 5 1548 2098.70 30650.4 615260.5 3.84 1.63 1.56 5%| Model number | $\mathrm{f}_{0}$ (Fig.9) | f (Fig.11) | $\mathrm{m}_{\mathrm{eq}}$ (KN.sec ${ }^{2} / \mathrm{m}$ ) | $\mathrm{k}_{\text {eq }}$ (KN/m) (Eq.15) | $\omega$ (rad/sec) | $\mathrm{T}_{\text {proposed }}$ equation(sec) (Eq.28) | $\mathrm{T}_{\text {Fem }}$ (sec) | Error (%) | | :--- | :--- | :--- | :--- | :--- | :--- | :--- | :--- | :--- | | Model 1 | 1548 | 1583.69 | 35604 | 778590.5 | 4.67 | 1.34 | 1.282 | 5% | | Model 2 | 1548 | 1880.789 | 32701.5 | 684556.8 | 4.52 | 1.389 | 1.399 | 1% | | Model 3 | 1548 | 1693.86 | 31393.44 | 640356.4 | 4.09 | 1.53 | 1.492 | 3% | | Model 4 | 1548 | 2098.70 | 32167.44 | 666569.1 | 4.35 | 1.44 | 1.433 | 1% | | Model 5 | 1548 | 2098.70 | 30650.4 | 615260.5 | 3.84 | 1.63 | 1.56 | 5% |
Table 9: Comparative analysis of fundamental frequencies and time periods computed using the proposed analytical model and FEM for five structural models with corresponding mass and stiffness parameters.
In order to check the validity of the suggested formula for estimation of the fundamental time period of buildings, comparisons were also made with FEM results for different frame configurations, as presented in Figs. 13a and 13b. Fig. 13a shows the comparison for a height of 10 with different spans of 6 m , 7 m , 8 m , 9 m , 10 m , 12 m 6 m , 7 m , 8 m , 9 m , 10 m , 12 m 6m,7m,8m,9m,10m,12m6 \mathrm{~m}, 7 \mathrm{~m}, 8 \mathrm{~m}, 9 \mathrm{~m}, 10 \mathrm{~m}, 12 \mathrm{~m}6 m,7 m,8 m,9 m,10 m,12 m and 15 m . Here, the maximum error that was found was 5 % 5 % 5%5 \%5%. The proposed equation shows good concordance with FEM predictions, as most data points are in proximity to the 45 45 45^(@)45^{\circ}45 line. Statistical analysis also confirms this observation, as a high coefficient of determination ( R 2 R 2 R^(2)\mathrm{R}^{2}R2 ) of 0.86 and a low root mean square error (RMSE) of 0.0165 seconds indicate strong predictive accuracy. Fig. 13b shows the evaluation of a 20-story frame of constant 6-meter span, in which the ratio of alpha changes continuously from 0.24 to 0.21 , 0.17 0.21 , 0.17 0.21,0.170.21,0.170.21,0.17, and 0.14 . As in the former case, data points are very close to the 45 45 45^(@)45^{\circ}45 line, indicating a high degree of agreement between FEM results and the provided equation. This case yields even more impressive statistical results, with R 2 = 0.962 R 2 = 0.962 R^(2)=0.962\mathrm{R}^{2}=0.962R2=0.962 and RMSE = 0.0308 = 0.0308 =0.0308=0.0308=0.0308 seconds, demonstrating the powerful ability of the equation to represent the dynamic response of taller buildings. These results overall validate that the equation provided herein provides a very accurate estimate of the fundamental time period for a range of building configurations.
Figure 13: Scatter plots of the proposed analytical and estimated fundamental periods, a) 10 -storey building, and b) 20 -storey building.

Comprehensive validation of the proposed equation against FEM results and established FORMULAS FROM LITERATURE

The accuracy of the suggested formula for estimating the fundamental vibration period of reinforced concrete moment-resisting frame (RC MRF) structures is confirmed in Fig. 14 by comparing it with existing formulations from Goel and Chopra [5], Salama [6], Aninthaneni and Dhakal [23], as well as seismic design codes such as ASCE 7-05 [32] and UBC [1], using Finite Element Method (FEM) results as reference. Two structural layouts are considered: frames with a 5 m span (Fig. 14a) and frames with a 7 m span (Fig. 14b). The comparison is based on two important statistical indicators: the coefficient of determination ( R 2 ) R 2 (R^(2))\left(\mathrm{R}^{2}\right)(R2), which measures the degree of correlation between predicted and reference values, and the root mean square error (RMSE), which quantifies the average prediction error.
For the 5 m span illustration case (Fig. 14a), the new equation is well correlated with the FEM solution with R 2 = 0.998 R 2 = 0.998 R^(2)=0.998\mathrm{R}^{2}=0.998R2=0.998 and RMSE = 0.039 s = 0.039 s =0.039s=0.039 \mathrm{~s}=0.039 s. Both indicators confirm the high accuracy and reliability of the new model in estimating RC MRF frames' dynamic behavior. For comparison, Aninthaneni and Dhakal's equation [23] gives fair agreement R 2 = 0.81 , RMSE = 0.50 R 2 = 0.81 , RMSE = 0.50 R^(2)=0.81,RMSE=0.50\mathrm{R}^{2}=0.81, \mathrm{RMSE}=0.50R2=0.81,RMSE=0.50 s ), whereas Goel and Chopra's equation [5] is poor ( R 2 = 0.14 R 2 = 0.14 R^(2)=0.14\mathrm{R}^{2}=0.14R2=0.14, RMSE = 1.09 s = 1.09 s =1.09s=1.09 \mathrm{~s}=1.09 s ) with considerable departure from FEMbased periods. Likewise, Salama's equation [6] exhibits poor predictive strength ( R 2 = 0.62 R 2 = 0.62 R^(2)=0.62\mathrm{R}^{2}=0.62R2=0.62, RMSE = 0.71 s = 0.71 s =0.71s=0.71 \mathrm{~s}=0.71 s ) due to overestimation of stiffness. Of the seismic code equations, ASCE 7-05 [32] is moderately consistent with FEM solutions ( R 2 = 0.80 ) R 2 = 0.80 (R^(2)=0.80)\left(\mathrm{R}^{2}=0.80\right)(R2=0.80) but tends to underestimate the fundamental period for higher structures, which may result in overly conservative seismic force estimates. The most conservative estimates overall are given by UBC [1], with lower predicted periods at every height and relatively high R 2 = 0.86 R 2 = 0.86 R^(2)=0.86\mathrm{R}^{2}=0.86R2=0.86, indicating trend conformity more than actual accuracy.
In the case of 7 m span (Fig. 14b), the proposed equation maintains the predicting proficiency with R 2 = 0.999 R 2 = 0.999 R^(2)=0.999\mathrm{R}^{2}=0.999R2=0.999 and RMSE = 0.027 s = 0.027 s =0.027s=0.027 \mathrm{~s}=0.027 s, and makes it robust and span-insensitive. Aninthaneni and Dhakal's model [23] in the present case shows enhanced performance ( R 2 = 0.93 R 2 = 0.93 R^(2)=0.93\mathrm{R}^{2}=0.93R2=0.93, RMSE = 0.29 s = 0.29 s =0.29s=0.29 \mathrm{~s}=0.29 s ), depicting higher applicability for wide spans. Goel and Chopra's equation [5] still performs poorly ( R 2 = 0.016 R 2 = 0.016 R^(2)=0.016\mathrm{R}^{2}=0.016R2=0.016, RMSE = 1.12 s = 1.12 s =1.12s=1.12 \mathrm{~s}=1.12 s ), whereas Salama's equation [6] still has moderate accuracy ( R 2 = 0.56 R 2 = 0.56 (R^(2)=0.56:}\left(\mathrm{R}^{2}=0.56\right.(R2=0.56, RMSE = 0.75 s ) = 0.75 s {:=0.75(s))\left.=0.75 \mathrm{~s}\right)=0.75 s). The ASCE 7-05 [32] model once more has moderate correlation ( R 2 = 0.80 R 2 = 0.80 R^(2)=0.80\mathrm{R}^{2}=0.80R2=0.80 ) but still underestimates periods for high-rise structures. The UBC [1] formula has R 2 = 0.84 R 2 = 0.84 R^(2)=0.84\mathrm{R}^{2}=0.84R2=0.84 but is still too conservative, severely underestimating the fundamental period for all heights.
Overall, for both span cases, the new equation outperforms any other empirical and code-based formula in terms of prediction accuracy (lowest RMSE) and magnitude of correlation (highest R 2 R 2 R^(2)\mathrm{R}^{2}R2 ). Its excellent agreement with FEM solutions, regardless of span length, demonstrates its utility in the seismic analysis and design of RC MRF buildings in practice. These outcomes validate its usage as a rational alternative to conventional empirical period estimation formulas.
Figure 14 : Comparison of fundamental vibration periods estimated using the proposed equation with those predicted by FEM analysis and existing formulas for reinforced concrete moment-resisting frame (RCMRF) buildings : (a) frames with a 5 m span ; (b) frames with a 7 m span.

Summary of recent studies on fundamental period estimation and their methodological CHARACTERISTICS

In addition to the comparison with empirical code-based formulas and FEM results, it is important to situate the proposed analytical model within the context of recent research. Tab. 10 is an overview of recent advances in fundamental period prediction methods, comparing structural scope, methodologies, parameters treated, and main findings across selected studies. The development of research provides a clear trend toward the integration of data-driven approaches, namely machine learning (ML), with traditional analytical models both to enhance accuracy and enable interpretability. Kumar et al. [25] explored RC moment-resisting frames via comparative analysis of ML algorithms (ANN, GP, RT), where ANN was most accurate. Their results emphasize the requirement for geometric parameters such as number of bays and bay length to influence period predictions. Rahman et al. [26] extended their scope to steel-braced, RC structures and employed ML models that were enhanced with SHAP-based interpretability. Shan et al. [27] addressed the issue of prediction of existing high-rise RC buildings by integrating probabilistic ML with explicit accounting for material deterioration and structural aging. Karampinis et al. [28] bridged analytical and ML approaches, explaining the influence of plan irregularities and stiffness ratios through SHAP values, while also producing analytical equations with comparable performance to ML models.
Relative to recent advances in fundamental period predictive techniques (Tab.10), the deterministic analytical model developed herein offers a novel compromise between predictability, computational efficiency, and transparency. While several recent studies have applied machine learning (ML) techniques either alone or in combinations to enhance predictability, these approaches are often require large amount of data (with high-quality data sets), model-intensive (with sophisticated training), and may be infested with interpretability challenges in the lack of specific explainability tools (e.g., SHAP) [25, 26, 27, 28]. Compared to the previous, the new model is an analytical model based on merely the basic structural parameters of mass, height, and stiffness distribution and therefore can be applied without the computational overhead or data dependence of ML-based methods. Moreover, unlike probabilistic or data-driven methods, the present approach gives closed-form solutions that are inherently interpretable and therefore can be applied directly in early design and code development.
Study & year Structural Scope Methodology Key Parameters Considered Main Findings
Kumar et al. [25] RC momentresisting frames Algorithms for comparative machine learning (ANN, GP, RT) Height, length of the bays, number of bays, and material The influence of geometric irregularities was highlighted by ANN, which achieved the highest accuracy
Rahman et al. [26] RC buildings with steel bracing Interpretable machine learning (SHAP values) Height, bracing arrangement, and stiffness distribution Period was greatly shortened by bracing; important stiffness variables were identified by ML interpretability
Shan et al. [27] RC high-rise buildings Probabilistic machine learning building height, and material deterioration For older structures, probabilistic machine learning captured uncertainty and enhanced prediction.
Karampinis et al. [28] Structures with frames ML + Analytical hybrid with SHAP Height, stiffness ratios, and plan irregularities Created comprehensible analytical formulas with precision on par with the ML
Proposed Model RC momentresisting frames Deterministic Analytical Model Height, mass and stiffness distribution, Easy to understand, transparent, and computationally efficient; similar accuracy for a standard configuration
Study & year Structural Scope Methodology Key Parameters Considered Main Findings Kumar et al. [25] RC momentresisting frames Algorithms for comparative machine learning (ANN, GP, RT) Height, length of the bays, number of bays, and material The influence of geometric irregularities was highlighted by ANN, which achieved the highest accuracy Rahman et al. [26] RC buildings with steel bracing Interpretable machine learning (SHAP values) Height, bracing arrangement, and stiffness distribution Period was greatly shortened by bracing; important stiffness variables were identified by ML interpretability Shan et al. [27] RC high-rise buildings Probabilistic machine learning building height, and material deterioration For older structures, probabilistic machine learning captured uncertainty and enhanced prediction. Karampinis et al. [28] Structures with frames ML + Analytical hybrid with SHAP Height, stiffness ratios, and plan irregularities Created comprehensible analytical formulas with precision on par with the ML Proposed Model RC momentresisting frames Deterministic Analytical Model Height, mass and stiffness distribution, Easy to understand, transparent, and computationally efficient; similar accuracy for a standard configuration| Study & year | Structural Scope | Methodology | Key Parameters Considered | Main Findings | | :--- | :--- | :--- | :--- | :--- | | Kumar et al. [25] | RC momentresisting frames | Algorithms for comparative machine learning (ANN, GP, RT) | Height, length of the bays, number of bays, and material | The influence of geometric irregularities was highlighted by ANN, which achieved the highest accuracy | | Rahman et al. [26] | RC buildings with steel bracing | Interpretable machine learning (SHAP values) | Height, bracing arrangement, and stiffness distribution | Period was greatly shortened by bracing; important stiffness variables were identified by ML interpretability | | Shan et al. [27] | RC high-rise buildings | Probabilistic machine learning | building height, and material deterioration | For older structures, probabilistic machine learning captured uncertainty and enhanced prediction. | | Karampinis et al. [28] | Structures with frames | ML + Analytical hybrid with SHAP | Height, stiffness ratios, and plan irregularities | Created comprehensible analytical formulas with precision on par with the ML | | Proposed Model | RC momentresisting frames | Deterministic Analytical Model | Height, mass and stiffness distribution, | Easy to understand, transparent, and computationally efficient; similar accuracy for a standard configuration |
Table 10. Summary of recent studies on fundamental period prediction, highlighting structural scope, methodology, key parameters, and main findings, with emphasis on the positioning of the proposed deterministic analytical model relative to contemporary machine learning and hybrid approaches.

Sensitivity investigations

Sensitivity analysis was conducted to evaluate the influence of the main design parameters on the fundamental vibration period of RC moment-resisting frame (MRF) buildings. This investigation is based on numerical results of 125 building models, where the primary parameters are shown in Tab.11. They include the number of stories, total building height, typical storey height, span length, cross-sectional beam and column sizes, model mass, and concrete modulus of elasticity. The results of the sensitivity analysis are outlined in the following subsections.

Geometric parameters

The building's geometrical configuration is essential in determining its dynamic features, specifically its fundamental time period. The building's height and span length are among the most significant characteristics, as they directly affect the structure's overall stiffness and mass distribution. This subsection provides a comprehensive examination of how various elements affect the fundamental time period.

Effect of building height

Relative comparison provided in Fig. 15 (a-c) contrasts accuracy of different formulations in estimating the fundamental time period of moment-resisting frames having different spans ( 5 m , 6 m 5 m , 6 m 5m,6m5 \mathrm{~m}, 6 \mathrm{~m}5 m,6 m, and 7 m ), with Finite Element Method (FEM) results. The equations considered are those given by [ 5 , 6 , 23 ] [ 5 , 6 , 23 ] [5,6,23][5,6,23][5,6,23] and a new expression developed in the present work.
In Fig. 15a, for 5 m span frames, the proposed equation is in very good agreement with FEM-based calculations with R 2 = R 2 = R^(2)=\mathrm{R}^{2}=R2= 0.999 and RMSE = 0.02 s = 0.02 s =0.02s=0.02 \mathrm{~s}=0.02 s. Aninthaneni and Dhakal's equation [23], however, has R 2 = 0.60 R 2 = 0.60 R^(2)=0.60\mathrm{R}^{2}=0.60R2=0.60 and RMSE = 0.67 s = 0.67 s =0.67s=0.67 \mathrm{~s}=0.67 s, indicating a moderate relationship with some scatter. Goel and Chopra's formula [5] exhibits suboptimal performance, yielding an R 2 R 2 R^(2)\mathrm{R}^{2}R2 of 0.36 and a root mean square error (RMSE) of 0.84 s , whereas Salama's formula [6] is more accurate ( R 2 = 0.78 , RMSE = R 2 = 0.78 , RMSE = R^(2)=0.78,RMSE=\mathrm{R}^{2}=0.78, \mathrm{RMSE}=R2=0.78,RMSE= 0.49 s ) but less than the proposed model. For Fig. 15b, for a frame with a span of 6 m , the proposed model maintains high predictive fidelity ( R 2 = 0.999 R 2 = 0.999 R^(2)=0.999\mathrm{R}^{2}=0.999R2=0.999, RMSE = 0.013 s = 0.013 s =0.013s=0.013 \mathrm{~s}=0.013 s ), underscoring its robustness across structural configurations. Aninthaneni and Dhakal's equation [23] shows some improvement ( R 2 = 0.66 , RMSE = 0.55 s R 2 = 0.66 , RMSE = 0.55 s R^(2)=0.66,RMSE=0.55s\mathrm{R}^{2}=0.66, \mathrm{RMSE}=0.55 \mathrm{~s}R2=0.66,RMSE=0.55 s ), but yet maintains a very high margin of error. Goel and Chopra [5] approach is still imprecise ( R 2 = 0.38 R 2 = 0.38 R^(2)=0.38\mathrm{R}^{2}=0.38R2=0.38, RMSE = 0.76 s = 0.76 s =0.76s=0.76 \mathrm{~s}=0.76 s ), while Salama's formula shows quite good agreement ( R 2 = 0.83 R 2 = 0.83 (R^(2)=0.83:}\left(\mathrm{R}^{2}=0.83\right.(R2=0.83, RMSE = 0.40 s ) = 0.40 s {:=0.40(s))\left.=0.40 \mathrm{~s}\right)=0.40 s), though not as good as the suggested expression.
Figure 15: Comparison of fundamental periods for buildings with varying models at different spans: (a) 5 m , (b) 6 m , and (c) 7 m .
Model number Number of stories Building height (m) Storey height (m) Frame span (m) Column Dimension (m) Beam dimension (m) T equation (sec) T fem T fem  T_("fem ")\mathrm{T}_{\text {fem }}Tfem  (sec) Error %
1 3 9 3 5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5*0.5 0.169 0.17 1%
2 5 15 3 5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5*0.5 0.31 0.304 2%
3 6 18 3 5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.36 0.373 3%
4 7 21 3 5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.45 0.445 1%
5 8 24 3 5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.52 0.519 0%
6 9 27 3 5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.61 0.596 2%
7 10 30 3 5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.68 0.675 1%
8 12 36 3 5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.84 0.843 0%
9 15 45 3 5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 1.1 1.121 2%
10 20 60 3 5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 1.67 1.666 0%
11 22 66 3 5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 1.98 1.915 3%
12 25 75 3 5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 2.3 2.325 1%
13 27 81 3 5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 2.62 2.623 0%
14 29 87 3 5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5*0.5 2.95 2.942 0%
15 30 90 3 5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 3.05 3.109 2%
16 3 9 3 6 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.18 0.188 4%
17 5 15 3 6 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.32 0.335 4%
18 6 18 3 6 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.4 0.411 3%
19 7 21 3 6 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.48 0.489 2%
20 8 24 3 6 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.58 0.568 2%
21 9 27 3 6 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5*0.5 0.66 0.649 2%
22 10 30 3 6 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.73 0.731 0%
23 12 36 3 6 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.9 0.903 0%
24 15 45 3 6 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 1.2 1.18 2%
25 20 60 3 6 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 1.75 1.702 3%
26 22 66 3 6 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 1.98 1.934 2%
27 25 75 3 6 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 2.3 2.311 0%
28 27 81 3 6 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 2.6 2.582 1%
29 29 87 3 6 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 2.98 2.87 4%
30 30 90 3 6 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 3 3.02 1%
31 3 9 3 7 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.2 0.207 3%
32 5 15 3 7 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.35 0.368 5%
33 6 18 3 7 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.46 0.45 2%
34 7 21 3 7 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.53 0.534 1%
35 8 24 3 7 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.61 0.62 2%
36 9 27 3 7 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.7 0.706 1%
37 10 30 3 7 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.78 0.794 2%
38 12 36 3 7 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.98 0.974 1%
39 15 45 3 7 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 1.3 1.259 3%
40 20 60 3 7 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 1.7 1.78 4%
41 22 66 3 7 0.5*0.5 0.5*0.5 2 2.006 0%
Model number Number of stories Building height (m) Storey height (m) Frame span (m) Column Dimension (m) Beam dimension (m) T equation (sec) T_("fem ") (sec) Error % 1 3 9 3 5 0.5**0.5 0.5*0.5 0.169 0.17 1% 2 5 15 3 5 0.5**0.5 0.5*0.5 0.31 0.304 2% 3 6 18 3 5 0.5**0.5 0.5**0.5 0.36 0.373 3% 4 7 21 3 5 0.5**0.5 0.5**0.5 0.45 0.445 1% 5 8 24 3 5 0.5**0.5 0.5**0.5 0.52 0.519 0% 6 9 27 3 5 0.5**0.5 0.5**0.5 0.61 0.596 2% 7 10 30 3 5 0.5**0.5 0.5**0.5 0.68 0.675 1% 8 12 36 3 5 0.5**0.5 0.5**0.5 0.84 0.843 0% 9 15 45 3 5 0.5**0.5 0.5**0.5 1.1 1.121 2% 10 20 60 3 5 0.5**0.5 0.5**0.5 1.67 1.666 0% 11 22 66 3 5 0.5**0.5 0.5**0.5 1.98 1.915 3% 12 25 75 3 5 0.5**0.5 0.5**0.5 2.3 2.325 1% 13 27 81 3 5 0.5**0.5 0.5**0.5 2.62 2.623 0% 14 29 87 3 5 0.5**0.5 0.5*0.5 2.95 2.942 0% 15 30 90 3 5 0.5**0.5 0.5**0.5 3.05 3.109 2% 16 3 9 3 6 0.5**0.5 0.5**0.5 0.18 0.188 4% 17 5 15 3 6 0.5**0.5 0.5**0.5 0.32 0.335 4% 18 6 18 3 6 0.5**0.5 0.5**0.5 0.4 0.411 3% 19 7 21 3 6 0.5**0.5 0.5**0.5 0.48 0.489 2% 20 8 24 3 6 0.5**0.5 0.5**0.5 0.58 0.568 2% 21 9 27 3 6 0.5**0.5 0.5*0.5 0.66 0.649 2% 22 10 30 3 6 0.5**0.5 0.5**0.5 0.73 0.731 0% 23 12 36 3 6 0.5**0.5 0.5**0.5 0.9 0.903 0% 24 15 45 3 6 0.5**0.5 0.5**0.5 1.2 1.18 2% 25 20 60 3 6 0.5**0.5 0.5**0.5 1.75 1.702 3% 26 22 66 3 6 0.5**0.5 0.5**0.5 1.98 1.934 2% 27 25 75 3 6 0.5**0.5 0.5**0.5 2.3 2.311 0% 28 27 81 3 6 0.5**0.5 0.5**0.5 2.6 2.582 1% 29 29 87 3 6 0.5**0.5 0.5**0.5 2.98 2.87 4% 30 30 90 3 6 0.5**0.5 0.5**0.5 3 3.02 1% 31 3 9 3 7 0.5**0.5 0.5**0.5 0.2 0.207 3% 32 5 15 3 7 0.5**0.5 0.5**0.5 0.35 0.368 5% 33 6 18 3 7 0.5**0.5 0.5**0.5 0.46 0.45 2% 34 7 21 3 7 0.5**0.5 0.5**0.5 0.53 0.534 1% 35 8 24 3 7 0.5**0.5 0.5**0.5 0.61 0.62 2% 36 9 27 3 7 0.5**0.5 0.5**0.5 0.7 0.706 1% 37 10 30 3 7 0.5**0.5 0.5**0.5 0.78 0.794 2% 38 12 36 3 7 0.5**0.5 0.5**0.5 0.98 0.974 1% 39 15 45 3 7 0.5**0.5 0.5**0.5 1.3 1.259 3% 40 20 60 3 7 0.5**0.5 0.5**0.5 1.7 1.78 4% 41 22 66 3 7 0.5*0.5 0.5*0.5 2 2.006 0%| Model number | Number of stories | Building height (m) | Storey height (m) | Frame span (m) | Column Dimension (m) | Beam dimension (m) | T equation (sec) | $\mathrm{T}_{\text {fem }}$ (sec) | Error % | | :--- | :--- | :--- | :--- | :--- | :--- | :--- | :--- | :--- | :--- | | 1 | 3 | 9 | 3 | 5 | $0.5 * 0.5$ | 0.5*0.5 | 0.169 | 0.17 | 1% | | 2 | 5 | 15 | 3 | 5 | $0.5 * 0.5$ | 0.5*0.5 | 0.31 | 0.304 | 2% | | 3 | 6 | 18 | 3 | 5 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.36 | 0.373 | 3% | | 4 | 7 | 21 | 3 | 5 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.45 | 0.445 | 1% | | 5 | 8 | 24 | 3 | 5 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.52 | 0.519 | 0% | | 6 | 9 | 27 | 3 | 5 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.61 | 0.596 | 2% | | 7 | 10 | 30 | 3 | 5 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.68 | 0.675 | 1% | | 8 | 12 | 36 | 3 | 5 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.84 | 0.843 | 0% | | 9 | 15 | 45 | 3 | 5 | $0.5 * 0.5$ | $0.5 * 0.5$ | 1.1 | 1.121 | 2% | | 10 | 20 | 60 | 3 | 5 | $0.5 * 0.5$ | $0.5 * 0.5$ | 1.67 | 1.666 | 0% | | 11 | 22 | 66 | 3 | 5 | $0.5 * 0.5$ | $0.5 * 0.5$ | 1.98 | 1.915 | 3% | | 12 | 25 | 75 | 3 | 5 | $0.5 * 0.5$ | $0.5 * 0.5$ | 2.3 | 2.325 | 1% | | 13 | 27 | 81 | 3 | 5 | $0.5 * 0.5$ | $0.5 * 0.5$ | 2.62 | 2.623 | 0% | | 14 | 29 | 87 | 3 | 5 | $0.5 * 0.5$ | 0.5*0.5 | 2.95 | 2.942 | 0% | | 15 | 30 | 90 | 3 | 5 | $0.5 * 0.5$ | $0.5 * 0.5$ | 3.05 | 3.109 | 2% | | 16 | 3 | 9 | 3 | 6 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.18 | 0.188 | 4% | | 17 | 5 | 15 | 3 | 6 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.32 | 0.335 | 4% | | 18 | 6 | 18 | 3 | 6 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.4 | 0.411 | 3% | | 19 | 7 | 21 | 3 | 6 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.48 | 0.489 | 2% | | 20 | 8 | 24 | 3 | 6 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.58 | 0.568 | 2% | | 21 | 9 | 27 | 3 | 6 | $0.5 * 0.5$ | 0.5*0.5 | 0.66 | 0.649 | 2% | | 22 | 10 | 30 | 3 | 6 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.73 | 0.731 | 0% | | 23 | 12 | 36 | 3 | 6 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.9 | 0.903 | 0% | | 24 | 15 | 45 | 3 | 6 | $0.5 * 0.5$ | $0.5 * 0.5$ | 1.2 | 1.18 | 2% | | 25 | 20 | 60 | 3 | 6 | $0.5 * 0.5$ | $0.5 * 0.5$ | 1.75 | 1.702 | 3% | | 26 | 22 | 66 | 3 | 6 | $0.5 * 0.5$ | $0.5 * 0.5$ | 1.98 | 1.934 | 2% | | 27 | 25 | 75 | 3 | 6 | $0.5 * 0.5$ | $0.5 * 0.5$ | 2.3 | 2.311 | 0% | | 28 | 27 | 81 | 3 | 6 | $0.5 * 0.5$ | $0.5 * 0.5$ | 2.6 | 2.582 | 1% | | 29 | 29 | 87 | 3 | 6 | $0.5 * 0.5$ | $0.5 * 0.5$ | 2.98 | 2.87 | 4% | | 30 | 30 | 90 | 3 | 6 | $0.5 * 0.5$ | $0.5 * 0.5$ | 3 | 3.02 | 1% | | 31 | 3 | 9 | 3 | 7 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.2 | 0.207 | 3% | | 32 | 5 | 15 | 3 | 7 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.35 | 0.368 | 5% | | 33 | 6 | 18 | 3 | 7 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.46 | 0.45 | 2% | | 34 | 7 | 21 | 3 | 7 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.53 | 0.534 | 1% | | 35 | 8 | 24 | 3 | 7 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.61 | 0.62 | 2% | | 36 | 9 | 27 | 3 | 7 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.7 | 0.706 | 1% | | 37 | 10 | 30 | 3 | 7 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.78 | 0.794 | 2% | | 38 | 12 | 36 | 3 | 7 | $0.5 * 0.5$ | $0.5 * 0.5$ | 0.98 | 0.974 | 1% | | 39 | 15 | 45 | 3 | 7 | $0.5 * 0.5$ | $0.5 * 0.5$ | 1.3 | 1.259 | 3% | | 40 | 20 | 60 | 3 | 7 | $0.5 * 0.5$ | $0.5 * 0.5$ | 1.7 | 1.78 | 4% | | 41 | 22 | 66 | 3 | 7 | 0.5*0.5 | 0.5*0.5 | 2 | 2.006 | 0% |
Table 11: Comprehensive Sensitivity Analysis of Fundamental Time Periods Considering Variations in Building Height, Span Length, Beam-to-Column Flexural Stiffness Ratio ( I b / I c I b / I c I_(b)//I_(c)\mathrm{I}_{\mathrm{b}} / \mathrm{I}_{\mathrm{c}}Ib/Ic ), and Concrete Modulus of Elasticity (E)identify the most sensitive parameters affecting the fundamental period.
Model number Number of stories Building height (m) Storey height (m) Frame span (m) Column Dimension (m) Beam dimension (m) T equation (sec) T fem T fem  T_("fem ")\mathrm{T}_{\text {fem }}Tfem  (sec) Error %
42 25 75 3 7 0.5*0.5 0.5*0.5 2.35 2.368 1%
43 27 81 3 7 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 2.6 2.626 1%
44 29 87 3 7 0.5*0.5 0.5*0.5 2.89 2.896 0%
45 30 90 3 7 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 3 3.036 1%
46 20 60 3 5 0.6*0.6 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 1.63 1.639 1%
47 3 9 3 7 0.7*0.7 0.5*0.5 1.67 1.666 0%
48 5 15 3 7 0.8*0.8 0.5*0.5 1.7 1.719 1%
49 6 18 3 7 0.9*0.9 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 1.75 1.782 2%
50 7 21 3 7 1*1 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 1.83 1.848 1%
51 20 60 3 7 0.6*0.6 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 1.78 1.78 4%
52 20 60 3 7 0.7*0.7 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 1.85 1.782 0%
53 20 60 3 7 0.8*0.8 0.5*0.5 1.92 1.833 1%
54 20 60 3 7 0.9*0.9 0.5*0.5 1.94 1.905 1%
55 20 60 3 7 1*1 0.5*0.5 2 1.983 2%
56 20 60 3 4 1*1 1*1 1.35 1.382 2%
57 20 60 3 4.5 1*1 1*1 1.3 1.299 0%
58 20 60 3 5 1*1 1*1 1.25 1.239 1%
59 20 60 3 5.5 1*1 1*1 1.23 1.196 3%
60 20 60 3 6 1*1 1*1 1.2 1.167 3%
61 20 60 3 6.5 1*1 1*1 1.1 1.148 4%
62 20 60 3 7 1*1 1*1 1.13 1.138 1%
63 20 60 3 7.5 1*1 1*1 1.12 1.135 1%
64 20 60 3 8 1*1 1*1 1.13 1.137 1%
65 18 54 3 4 1*1 1*1 1.1 1.139 3%
66 18 54 3 4.5 1*1 1*1 1.09 1.076 1%
67 18 54 3 5 1*1 1*1 1.05 1.032 2%
68 18 54 3 5.5 0.5*0.5 0.5*0.5 1.01 1.001 1%
69 18 54 3 6 1*1 1*1 0.99 0.982 1%
70 18 54 3 6.5 1*1 1*1 0.97 0.972 0%
71 18 54 3 7 1*1 1*1 0.965 0.968 0%
72 18 54 3 7.5 1*1 1*1 0.96 0.97 1%
73 20 60 3 6 1*1 1*1 1.3 1.268 3%
74* 20 60 3 6 1*1 1*1 1.13 1.134 0%
75** 20 60 3 6 1*1 1*1 1.06 1.035 2%
76*** 20 60 3 6 1*1 1*1 0.95 0.959 1%
77 10 30 3 6 1*1 1*1 0.421 0.448 6%
78* 10 30 3 6 1*1 1*1 0.4 0.401 0%
79** 10 30 3 6 1*1 1*1 0.36 0.366 2%
80*** 10 30 3 6 0.34 0.342 1%
81 3 9 3 5 0.8*0.8 0.8*0.8 0.11 0.106 4%
82 5 15 3 5 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.8*0.8 0.19 0.195 3%
83 6 18 3 5 0.8*0.8 0.8*0.8 0.24 0.242 1%
Model number Number of stories Building height (m) Storey height (m) Frame span (m) Column Dimension (m) Beam dimension (m) T equation (sec) T_("fem ") (sec) Error % 42 25 75 3 7 0.5*0.5 0.5*0.5 2.35 2.368 1% 43 27 81 3 7 0.5**0.5 0.5**0.5 2.6 2.626 1% 44 29 87 3 7 0.5*0.5 0.5*0.5 2.89 2.896 0% 45 30 90 3 7 0.5**0.5 0.5**0.5 3 3.036 1% 46 20 60 3 5 0.6*0.6 0.5**0.5 1.63 1.639 1% 47 3 9 3 7 0.7*0.7 0.5*0.5 1.67 1.666 0% 48 5 15 3 7 0.8*0.8 0.5*0.5 1.7 1.719 1% 49 6 18 3 7 0.9*0.9 0.5**0.5 1.75 1.782 2% 50 7 21 3 7 1*1 0.5**0.5 1.83 1.848 1% 51 20 60 3 7 0.6*0.6 0.5**0.5 1.78 1.78 4% 52 20 60 3 7 0.7*0.7 0.5**0.5 1.85 1.782 0% 53 20 60 3 7 0.8*0.8 0.5*0.5 1.92 1.833 1% 54 20 60 3 7 0.9*0.9 0.5*0.5 1.94 1.905 1% 55 20 60 3 7 1*1 0.5*0.5 2 1.983 2% 56 20 60 3 4 1*1 1*1 1.35 1.382 2% 57 20 60 3 4.5 1*1 1*1 1.3 1.299 0% 58 20 60 3 5 1*1 1*1 1.25 1.239 1% 59 20 60 3 5.5 1*1 1*1 1.23 1.196 3% 60 20 60 3 6 1*1 1*1 1.2 1.167 3% 61 20 60 3 6.5 1*1 1*1 1.1 1.148 4% 62 20 60 3 7 1*1 1*1 1.13 1.138 1% 63 20 60 3 7.5 1*1 1*1 1.12 1.135 1% 64 20 60 3 8 1*1 1*1 1.13 1.137 1% 65 18 54 3 4 1*1 1*1 1.1 1.139 3% 66 18 54 3 4.5 1*1 1*1 1.09 1.076 1% 67 18 54 3 5 1*1 1*1 1.05 1.032 2% 68 18 54 3 5.5 0.5*0.5 0.5*0.5 1.01 1.001 1% 69 18 54 3 6 1*1 1*1 0.99 0.982 1% 70 18 54 3 6.5 1*1 1*1 0.97 0.972 0% 71 18 54 3 7 1*1 1*1 0.965 0.968 0% 72 18 54 3 7.5 1*1 1*1 0.96 0.97 1% 73 20 60 3 6 1*1 1*1 1.3 1.268 3% 74* 20 60 3 6 1*1 1*1 1.13 1.134 0% 75** 20 60 3 6 1*1 1*1 1.06 1.035 2% 76*** 20 60 3 6 1*1 1*1 0.95 0.959 1% 77 10 30 3 6 1*1 1*1 0.421 0.448 6% 78* 10 30 3 6 1*1 1*1 0.4 0.401 0% 79** 10 30 3 6 1*1 1*1 0.36 0.366 2% 80*** 10 30 3 6 0.34 0.342 1% 81 3 9 3 5 0.8*0.8 0.8*0.8 0.11 0.106 4% 82 5 15 3 5 0.8**0.8 0.8*0.8 0.19 0.195 3% 83 6 18 3 5 0.8*0.8 0.8*0.8 0.24 0.242 1%| Model number | Number of stories | Building height (m) | Storey height (m) | Frame span (m) | Column Dimension (m) | Beam dimension (m) | T equation (sec) | $\mathrm{T}_{\text {fem }}$ (sec) | Error % | | :--- | :--- | :--- | :--- | :--- | :--- | :--- | :--- | :--- | :--- | | 42 | 25 | 75 | 3 | 7 | 0.5*0.5 | 0.5*0.5 | 2.35 | 2.368 | 1% | | 43 | 27 | 81 | 3 | 7 | $0.5 * 0.5$ | $0.5 * 0.5$ | 2.6 | 2.626 | 1% | | 44 | 29 | 87 | 3 | 7 | 0.5*0.5 | 0.5*0.5 | 2.89 | 2.896 | 0% | | 45 | 30 | 90 | 3 | 7 | $0.5 * 0.5$ | $0.5 * 0.5$ | 3 | 3.036 | 1% | | 46 | 20 | 60 | 3 | 5 | 0.6*0.6 | $0.5 * 0.5$ | 1.63 | 1.639 | 1% | | 47 | 3 | 9 | 3 | 7 | 0.7*0.7 | 0.5*0.5 | 1.67 | 1.666 | 0% | | 48 | 5 | 15 | 3 | 7 | 0.8*0.8 | 0.5*0.5 | 1.7 | 1.719 | 1% | | 49 | 6 | 18 | 3 | 7 | 0.9*0.9 | $0.5 * 0.5$ | 1.75 | 1.782 | 2% | | 50 | 7 | 21 | 3 | 7 | 1*1 | $0.5 * 0.5$ | 1.83 | 1.848 | 1% | | 51 | 20 | 60 | 3 | 7 | 0.6*0.6 | $0.5 * 0.5$ | 1.78 | 1.78 | 4% | | 52 | 20 | 60 | 3 | 7 | 0.7*0.7 | $0.5 * 0.5$ | 1.85 | 1.782 | 0% | | 53 | 20 | 60 | 3 | 7 | 0.8*0.8 | 0.5*0.5 | 1.92 | 1.833 | 1% | | 54 | 20 | 60 | 3 | 7 | 0.9*0.9 | 0.5*0.5 | 1.94 | 1.905 | 1% | | 55 | 20 | 60 | 3 | 7 | 1*1 | 0.5*0.5 | 2 | 1.983 | 2% | | 56 | 20 | 60 | 3 | 4 | 1*1 | 1*1 | 1.35 | 1.382 | 2% | | 57 | 20 | 60 | 3 | 4.5 | 1*1 | 1*1 | 1.3 | 1.299 | 0% | | 58 | 20 | 60 | 3 | 5 | 1*1 | 1*1 | 1.25 | 1.239 | 1% | | 59 | 20 | 60 | 3 | 5.5 | 1*1 | 1*1 | 1.23 | 1.196 | 3% | | 60 | 20 | 60 | 3 | 6 | 1*1 | 1*1 | 1.2 | 1.167 | 3% | | 61 | 20 | 60 | 3 | 6.5 | 1*1 | 1*1 | 1.1 | 1.148 | 4% | | 62 | 20 | 60 | 3 | 7 | 1*1 | 1*1 | 1.13 | 1.138 | 1% | | 63 | 20 | 60 | 3 | 7.5 | 1*1 | 1*1 | 1.12 | 1.135 | 1% | | 64 | 20 | 60 | 3 | 8 | 1*1 | 1*1 | 1.13 | 1.137 | 1% | | 65 | 18 | 54 | 3 | 4 | 1*1 | 1*1 | 1.1 | 1.139 | 3% | | 66 | 18 | 54 | 3 | 4.5 | 1*1 | 1*1 | 1.09 | 1.076 | 1% | | 67 | 18 | 54 | 3 | 5 | 1*1 | 1*1 | 1.05 | 1.032 | 2% | | 68 | 18 | 54 | 3 | 5.5 | 0.5*0.5 | 0.5*0.5 | 1.01 | 1.001 | 1% | | 69 | 18 | 54 | 3 | 6 | 1*1 | 1*1 | 0.99 | 0.982 | 1% | | 70 | 18 | 54 | 3 | 6.5 | 1*1 | 1*1 | 0.97 | 0.972 | 0% | | 71 | 18 | 54 | 3 | 7 | 1*1 | 1*1 | 0.965 | 0.968 | 0% | | 72 | 18 | 54 | 3 | 7.5 | 1*1 | 1*1 | 0.96 | 0.97 | 1% | | 73 | 20 | 60 | 3 | 6 | 1*1 | 1*1 | 1.3 | 1.268 | 3% | | 74* | 20 | 60 | 3 | 6 | 1*1 | 1*1 | 1.13 | 1.134 | 0% | | 75** | 20 | 60 | 3 | 6 | 1*1 | 1*1 | 1.06 | 1.035 | 2% | | 76*** | 20 | 60 | 3 | 6 | 1*1 | 1*1 | 0.95 | 0.959 | 1% | | 77 | 10 | 30 | 3 | 6 | 1*1 | 1*1 | 0.421 | 0.448 | 6% | | 78* | 10 | 30 | 3 | 6 | 1*1 | 1*1 | 0.4 | 0.401 | 0% | | 79** | 10 | 30 | 3 | 6 | 1*1 | 1*1 | 0.36 | 0.366 | 2% | | 80*** | 10 | 30 | 3 | 6 | | | 0.34 | 0.342 | 1% | | 81 | 3 | 9 | 3 | 5 | 0.8*0.8 | 0.8*0.8 | 0.11 | 0.106 | 4% | | 82 | 5 | 15 | 3 | 5 | $0.8 * 0.8$ | 0.8*0.8 | 0.19 | 0.195 | 3% | | 83 | 6 | 18 | 3 | 5 | 0.8*0.8 | 0.8*0.8 | 0.24 | 0.242 | 1% |
Table 11:. continued .....
Model number Number of stories Building height (m) Storey height (m) Frame span (m) Column Dimension (m) Beam dimension (m) T equation (sec) T fem T fem  T_("fem ")\mathrm{T}_{\text {fem }}Tfem  (sec) Error %
84 7 21 3 5 0.8*0.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.292 0.293 0%
85 8 24 3 5 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.345 0.348 1%
86 9 27 3 5 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.4 0.405 1%
87 10 30 3 5 0.8*0.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.47 0.467 1%
88 12 36 3 5 0.8*0.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.6 0.603 0%
89 15 45 3 5 0.8*0.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.83 0.84 1%
90 20 60 3 5 0.8*0.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 1.34 1.333 1%
91 22 66 3 5 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 1.56 1.566 0%
92 25 75 3 5 0.8 0.8 0.8 0.8 0.8^(**)0.80.8{ }^{*} 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 2 1.955 2%
93 27 81 3 5 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 2.2 2.24 2%
94 29 87 3 5 0.8*0.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 2.5 2.547 2%
95 30 90 3 5 0.8*0.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 2.71 2.708 0%
96 3 9 3 6 0.8*0.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.11 0.117 6%
97 5 15 3 6 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.22 0.212 4%
98 6 18 3 6 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.26 0.262 1%
99 7 21 3 6 0.8*0.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.31 0.315 2%
100 8 24 3 6 0.8*0.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.37 0.37 0%
101 9 27 3 6 0.8*0.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.42 0.427 2%
102 10 30 3 6 0.8 0.8 0.8 0.8 0.8^(**)0.80.8{ }^{*} 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.488 0.488 0%
103 12 36 3 6 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.6 0.618 3%
104 15 45 3 6 0.8*0.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.83 0.839 1%
105 20 60 3 6 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 1.28 1.285 0%
106 22 66 3 6 0.8 0.8 0.8 0.8 0.8^(**)0.80.8{ }^{*} 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 1.49 1.493 0%
107 25 75 3 6 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 1.8 1.838 2%
108 27 81 3 6 0.8 0.8 0.8 0.8 0.8^(**)0.80.8{ }^{*} 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 2.1 2.09 0%
109 29 87 3 6 0.8 0.8 0.8 0.8 0.8^(**)0.80.8{ }^{*} 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 2.34 2.361 1%
110 30 90 3 6 0.8 0.8 0.8 0.8 0.8^(**)0.80.8{ }^{*} 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 2.5 2.503 0%
111 3 9 3 7 0.8*0.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.126 0.128 2%
112 5 15 3 7 0.8 0.8 0.8 0.8 0.8^(**)0.80.8{ }^{*} 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.23 0.231 0%
113 6 18 3 7 0.8 0.8 0.8 0.8 0.8^(**)0.80.8{ }^{*} 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.29 0.284 2%
114 7 21 3 7 0.8 0.8 0.8 0.8 0.8^(**)0.80.8{ }^{*} 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.34 0.34 0%
115 8 24 3 7 0.8 0.8 0.8 0.8 0.8^(**)0.80.8{ }^{*} 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.4 0.397 1%
116 9 27 3 7 0.8 0.8 0.8 0.8 0.8^(**)0.80.8{ }^{*} 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.45 0.456 1%
117 10 30 3 7 0.8 0.8 0.8 0.8 0.8^(**)0.80.8{ }^{*} 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.5 0.517 3%
118 12 36 3 7 0.8 0.8 0.8 0.8 0.8^(**)0.80.8{ }^{*} 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.64 0.647 1%
119 15 45 3 7 0.8 0.8 0.8 0.8 0.8^(**)0.80.8{ }^{*} 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.9 0.862 4%
120 20 60 3 7 0.8 0.8 0.8 0.8 0.8^(**)0.80.8{ }^{*} 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 1.29 1.283 1%
121 22 66 3 7 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 1.5 1.475 2%
122 25 75 3 7 0.8 0.8 0.8 0.8 0.8^(**)0.80.8{ }^{*} 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 1.785 1.792 0%
123 27 81 3 7 0.8 0.8 0.8 0.8 0.8^(**)0.80.8{ }^{*} 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 2 2.022 1%
124 29 87 3 7 0.8 0.8 0.8 0.8 0.8^(**)0.80.8{ }^{*} 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 2.26 2.268 0%
125 30 90 3 7 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 0.8 0.8 0.8 0.8 0.8**0.80.8 * 0.80.80.8 2.35 2.397 2%
Model number Number of stories Building height (m) Storey height (m) Frame span (m) Column Dimension (m) Beam dimension (m) T equation (sec) T_("fem ") (sec) Error % 84 7 21 3 5 0.8*0.8 0.8**0.8 0.292 0.293 0% 85 8 24 3 5 0.8**0.8 0.8**0.8 0.345 0.348 1% 86 9 27 3 5 0.8**0.8 0.8**0.8 0.4 0.405 1% 87 10 30 3 5 0.8*0.8 0.8**0.8 0.47 0.467 1% 88 12 36 3 5 0.8*0.8 0.8**0.8 0.6 0.603 0% 89 15 45 3 5 0.8*0.8 0.8**0.8 0.83 0.84 1% 90 20 60 3 5 0.8*0.8 0.8**0.8 1.34 1.333 1% 91 22 66 3 5 0.8**0.8 0.8**0.8 1.56 1.566 0% 92 25 75 3 5 0.8^(**)0.8 0.8**0.8 2 1.955 2% 93 27 81 3 5 0.8**0.8 0.8**0.8 2.2 2.24 2% 94 29 87 3 5 0.8*0.8 0.8**0.8 2.5 2.547 2% 95 30 90 3 5 0.8*0.8 0.8**0.8 2.71 2.708 0% 96 3 9 3 6 0.8*0.8 0.8**0.8 0.11 0.117 6% 97 5 15 3 6 0.8**0.8 0.8**0.8 0.22 0.212 4% 98 6 18 3 6 0.8**0.8 0.8**0.8 0.26 0.262 1% 99 7 21 3 6 0.8*0.8 0.8**0.8 0.31 0.315 2% 100 8 24 3 6 0.8*0.8 0.8**0.8 0.37 0.37 0% 101 9 27 3 6 0.8*0.8 0.8**0.8 0.42 0.427 2% 102 10 30 3 6 0.8^(**)0.8 0.8**0.8 0.488 0.488 0% 103 12 36 3 6 0.8**0.8 0.8**0.8 0.6 0.618 3% 104 15 45 3 6 0.8*0.8 0.8**0.8 0.83 0.839 1% 105 20 60 3 6 0.8**0.8 0.8**0.8 1.28 1.285 0% 106 22 66 3 6 0.8^(**)0.8 0.8**0.8 1.49 1.493 0% 107 25 75 3 6 0.8**0.8 0.8**0.8 1.8 1.838 2% 108 27 81 3 6 0.8^(**)0.8 0.8**0.8 2.1 2.09 0% 109 29 87 3 6 0.8^(**)0.8 0.8**0.8 2.34 2.361 1% 110 30 90 3 6 0.8^(**)0.8 0.8**0.8 2.5 2.503 0% 111 3 9 3 7 0.8*0.8 0.8**0.8 0.126 0.128 2% 112 5 15 3 7 0.8^(**)0.8 0.8**0.8 0.23 0.231 0% 113 6 18 3 7 0.8^(**)0.8 0.8**0.8 0.29 0.284 2% 114 7 21 3 7 0.8^(**)0.8 0.8**0.8 0.34 0.34 0% 115 8 24 3 7 0.8^(**)0.8 0.8**0.8 0.4 0.397 1% 116 9 27 3 7 0.8^(**)0.8 0.8**0.8 0.45 0.456 1% 117 10 30 3 7 0.8^(**)0.8 0.8**0.8 0.5 0.517 3% 118 12 36 3 7 0.8^(**)0.8 0.8**0.8 0.64 0.647 1% 119 15 45 3 7 0.8^(**)0.8 0.8**0.8 0.9 0.862 4% 120 20 60 3 7 0.8^(**)0.8 0.8**0.8 1.29 1.283 1% 121 22 66 3 7 0.8**0.8 0.8**0.8 1.5 1.475 2% 122 25 75 3 7 0.8^(**)0.8 0.8**0.8 1.785 1.792 0% 123 27 81 3 7 0.8^(**)0.8 0.8**0.8 2 2.022 1% 124 29 87 3 7 0.8^(**)0.8 0.8**0.8 2.26 2.268 0% 125 30 90 3 7 0.8**0.8 0.8**0.8 2.35 2.397 2%| Model number | Number of stories | Building height (m) | Storey height (m) | Frame span (m) | Column Dimension (m) | Beam dimension (m) | T equation (sec) | $\mathrm{T}_{\text {fem }}$ (sec) | Error % | | :--- | :--- | :--- | :--- | :--- | :--- | :--- | :--- | :--- | :--- | | 84 | 7 | 21 | 3 | 5 | 0.8*0.8 | $0.8 * 0.8$ | 0.292 | 0.293 | 0% | | 85 | 8 | 24 | 3 | 5 | $0.8 * 0.8$ | $0.8 * 0.8$ | 0.345 | 0.348 | 1% | | 86 | 9 | 27 | 3 | 5 | $0.8 * 0.8$ | $0.8 * 0.8$ | 0.4 | 0.405 | 1% | | 87 | 10 | 30 | 3 | 5 | 0.8*0.8 | $0.8 * 0.8$ | 0.47 | 0.467 | 1% | | 88 | 12 | 36 | 3 | 5 | 0.8*0.8 | $0.8 * 0.8$ | 0.6 | 0.603 | 0% | | 89 | 15 | 45 | 3 | 5 | 0.8*0.8 | $0.8 * 0.8$ | 0.83 | 0.84 | 1% | | 90 | 20 | 60 | 3 | 5 | 0.8*0.8 | $0.8 * 0.8$ | 1.34 | 1.333 | 1% | | 91 | 22 | 66 | 3 | 5 | $0.8 * 0.8$ | $0.8 * 0.8$ | 1.56 | 1.566 | 0% | | 92 | 25 | 75 | 3 | 5 | $0.8{ }^{*} 0.8$ | $0.8 * 0.8$ | 2 | 1.955 | 2% | | 93 | 27 | 81 | 3 | 5 | $0.8 * 0.8$ | $0.8 * 0.8$ | 2.2 | 2.24 | 2% | | 94 | 29 | 87 | 3 | 5 | 0.8*0.8 | $0.8 * 0.8$ | 2.5 | 2.547 | 2% | | 95 | 30 | 90 | 3 | 5 | 0.8*0.8 | $0.8 * 0.8$ | 2.71 | 2.708 | 0% | | 96 | 3 | 9 | 3 | 6 | 0.8*0.8 | $0.8 * 0.8$ | 0.11 | 0.117 | 6% | | 97 | 5 | 15 | 3 | 6 | $0.8 * 0.8$ | $0.8 * 0.8$ | 0.22 | 0.212 | 4% | | 98 | 6 | 18 | 3 | 6 | $0.8 * 0.8$ | $0.8 * 0.8$ | 0.26 | 0.262 | 1% | | 99 | 7 | 21 | 3 | 6 | 0.8*0.8 | $0.8 * 0.8$ | 0.31 | 0.315 | 2% | | 100 | 8 | 24 | 3 | 6 | 0.8*0.8 | $0.8 * 0.8$ | 0.37 | 0.37 | 0% | | 101 | 9 | 27 | 3 | 6 | 0.8*0.8 | $0.8 * 0.8$ | 0.42 | 0.427 | 2% | | 102 | 10 | 30 | 3 | 6 | $0.8{ }^{*} 0.8$ | $0.8 * 0.8$ | 0.488 | 0.488 | 0% | | 103 | 12 | 36 | 3 | 6 | $0.8 * 0.8$ | $0.8 * 0.8$ | 0.6 | 0.618 | 3% | | 104 | 15 | 45 | 3 | 6 | 0.8*0.8 | $0.8 * 0.8$ | 0.83 | 0.839 | 1% | | 105 | 20 | 60 | 3 | 6 | $0.8 * 0.8$ | $0.8 * 0.8$ | 1.28 | 1.285 | 0% | | 106 | 22 | 66 | 3 | 6 | $0.8{ }^{*} 0.8$ | $0.8 * 0.8$ | 1.49 | 1.493 | 0% | | 107 | 25 | 75 | 3 | 6 | $0.8 * 0.8$ | $0.8 * 0.8$ | 1.8 | 1.838 | 2% | | 108 | 27 | 81 | 3 | 6 | $0.8{ }^{*} 0.8$ | $0.8 * 0.8$ | 2.1 | 2.09 | 0% | | 109 | 29 | 87 | 3 | 6 | $0.8{ }^{*} 0.8$ | $0.8 * 0.8$ | 2.34 | 2.361 | 1% | | 110 | 30 | 90 | 3 | 6 | $0.8{ }^{*} 0.8$ | $0.8 * 0.8$ | 2.5 | 2.503 | 0% | | 111 | 3 | 9 | 3 | 7 | 0.8*0.8 | $0.8 * 0.8$ | 0.126 | 0.128 | 2% | | 112 | 5 | 15 | 3 | 7 | $0.8{ }^{*} 0.8$ | $0.8 * 0.8$ | 0.23 | 0.231 | 0% | | 113 | 6 | 18 | 3 | 7 | $0.8{ }^{*} 0.8$ | $0.8 * 0.8$ | 0.29 | 0.284 | 2% | | 114 | 7 | 21 | 3 | 7 | $0.8{ }^{*} 0.8$ | $0.8 * 0.8$ | 0.34 | 0.34 | 0% | | 115 | 8 | 24 | 3 | 7 | $0.8{ }^{*} 0.8$ | $0.8 * 0.8$ | 0.4 | 0.397 | 1% | | 116 | 9 | 27 | 3 | 7 | $0.8{ }^{*} 0.8$ | $0.8 * 0.8$ | 0.45 | 0.456 | 1% | | 117 | 10 | 30 | 3 | 7 | $0.8{ }^{*} 0.8$ | $0.8 * 0.8$ | 0.5 | 0.517 | 3% | | 118 | 12 | 36 | 3 | 7 | $0.8{ }^{*} 0.8$ | $0.8 * 0.8$ | 0.64 | 0.647 | 1% | | 119 | 15 | 45 | 3 | 7 | $0.8{ }^{*} 0.8$ | $0.8 * 0.8$ | 0.9 | 0.862 | 4% | | 120 | 20 | 60 | 3 | 7 | $0.8{ }^{*} 0.8$ | $0.8 * 0.8$ | 1.29 | 1.283 | 1% | | 121 | 22 | 66 | 3 | 7 | $0.8 * 0.8$ | $0.8 * 0.8$ | 1.5 | 1.475 | 2% | | 122 | 25 | 75 | 3 | 7 | $0.8{ }^{*} 0.8$ | $0.8 * 0.8$ | 1.785 | 1.792 | 0% | | 123 | 27 | 81 | 3 | 7 | $0.8{ }^{*} 0.8$ | $0.8 * 0.8$ | 2 | 2.022 | 1% | | 124 | 29 | 87 | 3 | 7 | $0.8{ }^{*} 0.8$ | $0.8 * 0.8$ | 2.26 | 2.268 | 0% | | 125 | 30 | 90 | 3 | 7 | $0.8 * 0.8$ | $0.8 * 0.8$ | 2.35 | 2.397 | 2% |
Table 11: continued .....
Note for Table 11: (), (), () denote models with different modulus of elasticity values, specifically 25 GPa, 30 GPa, and 35 GPa, respectively.
Fig. 15c, relative to a 7 m span, still confirms the goodness of the new model, which provides R 2 = 0.999 R 2 = 0.999 R^(2)=0.999\mathrm{R}^{2}=0.999R2=0.999 and RMSE = = === 0.018 s . Aninthaneni and Dhakal's formulation [23] has some improvement ( R 2 = 0.71 R 2 = 0.71 R^(2)=0.71\mathrm{R}^{2}=0.71R2=0.71, RMSE = 0.48 s = 0.48 s =0.48s=0.48 \mathrm{~s}=0.48 s ) but is still less accurate than that of the new model. Goel and Chopra's method [5] remains the least accurate ( R 2 = 0.37 R 2 = 0.37 R^(2)=0.37\mathrm{R}^{2}=0.37R2=0.37, RMSE = 0.72 s = 0.72 s =0.72s=0.72 \mathrm{~s}=0.72 s ), whereas Salama's equation [6] provides the most accurate estimate of the available models ( R 2 = 0.84 , RMSE = 0.36 s R 2 = 0.84 , RMSE = 0.36 s R^(2)=0.84,RMSE=0.36s\mathrm{R}^{2}=0.84, \mathrm{RMSE}=0.36 \mathrm{~s}R2=0.84,RMSE=0.36 s ).
In general, the statistical analysis for all three cases exhibits the continual outperformance of the new equation in predicting the fundamental time period of moment-resisting frames. Its R 2 R 2 R^(2)\mathrm{R}^{2}R2 values close to unity and substantially lower RMSE values than those of the existing models show both high precision and low error of prediction. This implies that the new model provides a more precise and generalized tool for the dynamic analysis of frame structures for a broad spectrum of changing span lengths.

Effect of span variation

The sensitivity study in Fig. 16a and Fig. 16b presents a comprehensive evaluation of the effect of building span length on the fundamental time period of reinforced concrete (RC) moment-resisting frames of two various building heights, 60 m and 54 m , respectively. The performance of the suggested analytical formula is compared to the established formulations of [ 5 , 6 , 23 ] [ 5 , 6 , 23 ] [5,6,23][5,6,23][5,6,23] with the benchmark being the finite element method (FEM).
As can be seen from the two figures, the FEM results show a reducing trend for the period over the increasing span length. This is due to the fact that there is more lateral stiffness in buildings with larger widths, and therefore, the natural period decreases. The equation we have developed predicts this trend accurately and agrees very well with the FEM results for the range of spans considered whereas the other models have large differences. To determine the predictive capability of both models, statistical measures in the form of the coefficient of determination ( R 2 R 2 R^(2)\mathrm{R}^{2}R2 ) and the Root Mean Square Error (RMSE) are provided. These measures criticize the correlation as well as predictive precision of both methods when compared to FEM results. Statistical comparison firmly confirms the higher performance of the new formula with the highest R 2 R 2 R^(2)\mathrm{R}^{2}R2 values ( 0.91 and 0.94 ) and the lowest RMSE values ( 0.025 and 0.016 ) for both building heights. These are indicators of high correlation with FEM results and low error in prediction and suggest the ability of the formula to capture the influence of variable span length on dynamic response.
In stark contrast, the other existing formulas all have negative R 2 R 2 R^(2)R^{2}R2 values, meaning that their predictions deviate so significantly from the FEM results. This is most evident for [5, 23], who have RMSE values above 0.4 and R 2 R 2 R^(2)\mathrm{R}^{2}R2 values below -40. Salama's formula [6] has a slightly better fit, but once more, it performs badly in comparison to the suggested method, having low correlation and moderate error. In conclusion, the comparison indicates the inability of the equations from the literature to effectively account for the influence of span length, especially in mid- to high-rise buildings. The suggested method is, however, invariably correct irrespective of height and span.
Figure 16: Variation of fundamental time period with span length for different building heights: (a) 60 m height; (b) 54 m height.
The member structural properties play an important role in the global lateral resistance and dynamic response of momentresisting frames. The beam to column inertia ratio and elastic modulus of members are parameters that have a direct influence on the global stiffness and thus the fundamental time period. Each of these parameters is treated separately in the following subsections to assess its own influence.

Effect of the beam to column inertia ratio

The sensitivity analysis shown in Fig. 17a and Fig. 17b evaluates the impact of the beam-to-column inertial ratio ( α α alpha\alphaα ) on the fundamental time period of reinforced concrete (RC) moment-resisting frames with spans of 5 m and 7 m , respectively. The work compares the performance of a new proposed equation with three others [ 5, 6, 23] based on results from finite element method (FEM) simulations as a benchmark reference.
In order to determine the validity of the proposed equation, two statistical measures are employed: the coefficient of determination ( R 2 R 2 R^(2)\mathrm{R}^{2}R2 ) and the Root Mean Square Error (RMSE). These statistical findings ratify the excellence of the new equation that in both ranges consistently provides the best approximation of FEM data. In Fig. 17a, the quasi-perfect R 2 = 0.95 R 2 = 0.95 R^(2)=0.95\mathrm{R}^{2}=0.95R2=0.95 and very low RMSE of 0.017 s reveal excellent agreement with numerical simulations. In Fig. 17b, although the correlation is slightly poorer ( R 2 = 0.85 R 2 = 0.85 R^(2)=0.85\mathrm{R}^{2}=0.85R2=0.85 ), the RMSE is still small at 0.045 s , which is still much better than all other models. Conversely, the existing equations, especially for the 7 m span, show extreme statistical deterioration, as reflected by large negative R 2 R 2 R^(2)\mathrm{R}^{2}R2 values and large RMSEs. Negative R 2 R 2 R^(2)\mathrm{R}^{2}R2 values of [ 5 , 6 , 23 ] [ 5 , 6 , 23 ] [5,6,23][5,6,23][5,6,23] in Fig. 17b show that these models cannot capture the trend in the FEM results. Their RMSE values ( 0.83 to 1.25 s ) also reflect their unsuitability for frames with larger spans and other stiffness properties. In conclusion, statistical analysis offers very strong evidence in favor of the adoption of the given method that possesses very high predictive power and good performance regarding the frame geometries. It takes into account the beam-to-column stiffness ratio, making it a better predictor of real structural behavior, a good and reliable tool for estimation of time periods in seismic frame design and analysis of RC frames.
Figure 17: Variation of the fundamental time period with respect to beam-to-column inertia ratio ( α = I b / Ic α = I b / Ic alpha=I_(b)//Ic\alpha=\mathrm{I}_{\mathrm{b}} / \mathrm{Ic}α=Ib/Ic ) for frames with spans of: (a) 5 m and (b) 7 m .

Effect of member elasticity

To investigate the influence of material stiffness, a sensitivity analysis is conducted by changing the value of Young's modulus over a reasonable range. For instance, the analysis of the frame model can be conducted using moduli of 20 GPa , 25 GPa , 30 GPa 25 GPa , 30 GPa 25GPa,30GPa25 \mathrm{GPa}, 30 \mathrm{GPa}25GPa,30GPa, and 35 GPa to model various material stiffnesses. A comparison between fundamental time periods from FEM analyses and the suggested equation in addition to the equations from literature [5, 6, 23] assists in determining the sensitivity of the equation to the variation of material properties.
Fig. 18 shows the variation in the fundamental time period of moment-resisting frames as a function of the modulus of elasticity (E), for 10-story and 20-story frames with a span of 6 m . The performance of the proposed equation is evaluated
against FEM results and benchmarked against existing equations by [5, 6, 23]. The statistical parameters coefficient of determination ( R 2 R 2 R^(2)\mathrm{R}^{2}R2 ) and root mean square error (RMSE) are used to assess the predictive accuracy of each equation.
In Fig. 18a, corresponding to the 10 -storey frame, the new formula demonstrates good consistency with FEM results with R 2 = 0.91 R 2 = 0.91 R^(2)=0.91\mathrm{R}^{2}=0.91R2=0.91 and low RMSE = 0.0146 s = 0.0146 s =0.0146s=0.0146 \mathrm{~s}=0.0146 s, which indicates high reliability and accuracy in capturing the role of material stiffness in dynamic response. The traditional formulas are of poor performance, like Aninthaneni and Dhakal's equation[23] provides a negative R 2 R 2 R^(2)\mathrm{R}^{2}R2 value of 1.71 and RMSE of 0.064 s , reflecting poor fit and inability to track the observed trend. Goel and Chopra's equation [5] gives R 2 = 1.3 R 2 = 1.3 R^(2)=-1.3\mathrm{R}^{2}=-1.3R2=1.3 and RMSE = 0.062 s = 0.062 s =0.062s=0.062 \mathrm{~s}=0.062 s, whereas Salama's equation [6] is the poorest with R 2 = R 2 = R^(2)=-\mathrm{R}^{2}=-R2= 14.4 and RMSE = 0.16 s = 0.16 s =0.16s=0.16 \mathrm{~s}=0.16 s. The negative R 2 R 2 R^(2)\mathrm{R}^{2}R2 for both of these models indicates that they are poorer in prediction, which further establishes their weakness.
In Fig. 18b, for the 20-storey frame, the superiority of the selected equation is clearer. It achieves an R 2 R 2 R^(2)\mathrm{R}^{2}R2 of 0.97 and an RMSE of 0.021 s , confirming its strength and versatility for taller structural configurations. Conversely, the remaining models show considerable disagreement with FEM results, with [5, 23] recording R 2 R 2 R^(2)\mathrm{R}^{2}R2 of -14.8 and -15.9 , and RMSE of 0.45 and 0.47 s , respectively. Although Salama's equation [6] shows comparatively better performance ( R 2 = 0.08 R 2 = 0.08 R^(2)=-0.08\mathrm{R}^{2}=-0.08R2=0.08, RMSE = 0.12 = 0.12 =0.12=0.12=0.12 s), it still fails to provide satisfactory accuracy. Generally, statistical comparison of 10 and 20 -story frames confirms that the new formula derived herein correctly outmatches existing empirical formulas in modeling the effect of modulus of elasticity on the fundamental time period. Its considerably higher R 2 R 2 R^(2)\mathrm{R}^{2}R2 values and significantly lower RMSEs confirm its adequacy for practical use in seismic analysis and design, especially where the change of material stiffness is of paramount concern to structural dynamics.
Figure 18: Variation of fundamental time period with respect to modulus of elasticity (E) for frames with 6 m span: (a) 10 -story frame, and (b) 20 -story frame, comparing FEM results with the proposed equation.
The total seismic mass of a building, comprising the weight of the structure and non-structural materials in addition to superimposed live loads, directly influences its fundamental period. An increase in mass is expected to result in a longer period, as it lowers the natural frequency of vibration for a given stiffness. To quantify this influence, a sensitivity analysis was performed for varying seismic mass values. Fig. 19 illustrates the chart of the fundamental period as a function of seismic mass comparing predictions from the suggested equation, the Aninthaneni & Dhakal model [23], Goel & Chopra model [5], and Salama model [6] with FEM results. As would be expected from fundamental vibration theory, the FEM computation does exhibit a clear-cut increasing trend for fundamental period with added seismic mass. The proposed equation accurately follows this trend and shows very good correlation with FEM computations of ( R 2 = 0.998 R 2 = 0.998 R^(2)=0.998\mathrm{R}^{2}=0.998R2=0.998, RMSE = = === 0.039), implying that it well-represents the dynamic mass-dependent behavior.
On the other hand, the Aninthaneni & Dhakal model [23], which is positively correlated ( R 2 = 0.81 R 2 = 0.81 R^(2)=0.81\mathrm{R}^{2}=0.81R2=0.81 ), substantially underestimates the period for heavy masses. This shows that the model is not entirely considering the influence of mass on the period for heavy structures. The Goel & Chopra [5] ( R 2 = 0.14 R 2 = 0.14 R^(2)=0.14\mathrm{R}^{2}=0.14R2=0.14 ) and Salama ( R 2 = 0.62 R 2 = 0.62 R^(2)=0.62\mathrm{R}^{2}=0.62R2=0.62 ) [6] models are considerably less responsive to mass change. Their predictions are relatively flat across the tested mass range, indicating that their models
implicitly perhaps are postulating a more rigid stiffness-dominated or geometry-governed period control with minimal explicit reference to seismic mass. As such, these models widely deviate from FEM projections for larger structures, leading to larger RMSE values ( 1.09 and 0.71 respectively). Overall, the discussion highlights the need for explicit and clear massdependence in empirical period equations, particularly in buildings whose seismic mass can change significantly by occupancy loads, non-structure elements, or large machinery. The proposed equation has this capability best, with performance prediction being accurate over a broad range of mass
Figure 19: Comparison of fundamental period variation with seismic mass as predicted by the proposed equation, Aninthaneni & Dhakal model, Goel & Chopra model, and Salama model against FEM results.

Key influencing parameters on the fundamental period: insights from sensitivity analysis

Sensitivity analysis highlighted that the fundamental period of moment-resisting frames is mostly governed by geometric parameters, in which height of the building is the most prominent parameter. Increase in height leads to a considerable increase in the period due to the significant reduction in the lateral stiffness and the associated increase in the lateral displacements. Seismic mass was the second most important parameter, since greater mass would elevate the system inertia and consequently the fundamental period. This consequence is always more definite in heavy structures, where code-based empirical relations simplifying or disregarding mass dependence do significantly deviate from finite element computations. Span length also had a measurable, though lesser, effect; longer spans reduce global stiffness, hence prolonging the period. Modulus of elasticity and beam-to-column inertia ratio, which are stiffness parameters, had fairly small effects in the typical ranges for reinforced concrete frames. The results generally indicate that seismic mass and height must be considered explicitly and accurately in predictive models and other parameters are second-order effective on the fundamental period.

\checkmark isual enhancement of sensitivity analysis results

To facilitate interpretation of the sensitivity analysis outcomes, the results are presented in the form of heat maps that allow for a direct visual comparison across different parameter combinations. Fig. 20 presents two heat maps that visually summarize the ranking of many parameter combinations according to their agreement with FEM results. To analyze the influence of different structural parameters on the accuracy of fundamental period prediction, heat maps were generated to separately evaluate the effect of the beam-to-column moment of inertia ratio ( α α alpha\alphaα ) and the span length, with variations examined across different numbers of stories. Fig. 20a illustrates the sensitivity of the fundamental period prediction error, classified into five levels ( 1 = 1 = 1=1=1= minimum error, 5 = 5 = 5=5=5= maximum error), as a function of the beam-to-column moment of inertia ratio ( α α alpha\alphaα ) and the number of stories. The results indicate that high values of α α alpha\alphaα (stiffer beams) give low errors for low- and mid-rise MRFs, but the error increases moderately with height. Very small values of α α alpha\alphaα (flexible beams) are always associated with higher errors, especially in tall MRFs, indicating that lower beam stiffness exaggerates discrepancies with FEM solutions. Intermediate α α alpha\alphaα values exhibit mixed behavior, performing well for mid- to high-rise cases but producing large errors in low-rise buildings. Overall, the trend reflects a non-linear, height-dependent relationship of α α alpha\alphaα with prediction accuracy, highlighting the requirement for proposed equations to explicitly account for stiffness ratios, particularly especially for low-rise models.
The heatmap shows in Fig. 20b how the level of prediction error for the fundamental period varies with frame span and building height (number of stories). Shorter spans ( 6 7 m 6 7 m 6-7m6-7 \mathrm{~m}67 m ) at lower error levels ( 1 2 1 2 1-21-212 ) are clustered, as well as at moderate heights ( 10 10 <= 10\leq 1010 stories), whereas taller buildings at moderate spans ( 9 12 m 9 12 m 9-12m9-12 \mathrm{~m}912 m ) perform erratically. The most significant error levels (4-5) occur mainly at long spans ( 12 m 12 m >= 12m\geq 12 \mathrm{~m}12 m ) with minimal stories or very tall structures, implying that both the extreme span-height ratios lead to greater deviation in period estimation.
Figure 20: Heat maps showing the sensitivity of fundamental period prediction error to variations in (a) α α alpha\alphaα versus number of stories and (b) span length versus number of stories.

Statistical indicators of agreement with FEM

To verify the efficacy of the provided analytical formula in estimating the fundamental time period, a comparison was also made using FEM results. Tab. 11 illustrates the comparison for a number of structural configurations that involve various span lengths, column and beam cross-sections, and building heights. The result shows a very good agreement between the values computed by the new approach and those computed using numerical analysis. FEM/proposed ranges from 0.96 to 1.04 , with an average value of approximately 1.00 . This reveals that the proposed model is not biased toward overestimation or underestimation of the fundamental period.
Besides, the coefficient of variation (COV) is just 2 % 2 % 2%2 \%2%, and the standard deviation of ratios is 0.02 , reflecting the uniformity and consistency of the proposed formulation. For the different structural configurations, the close agreement between numerical and analytical solutions verifies the capability of the proposed equation in simulating the dynamic properties of the framed buildings. This accuracy renders the proposed method an effective tool in early design and seismic analysis
Span (m) Col cross section (m) Beam cross section (m) Building height (m) Fundamental time period
FEM (sec) Proposed Equation (sec) FEM/proposed
4 0.4*0.4 0.3 0.3 0.3 0.3 0.3**0.30.3 * 0.30.30.3 30 0.99 0.95 1.04
4 0.4 0.4 0.4 0.4 0.4^(**)0.40.4{ }^{*} 0.40.40.4 0.4 0.4 0.4 0.4 0.4^(**)0.40.4{ }^{*} 0.40.40.4 30 0.775 0.77 1.006
4 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.4 0.4 0.4 0.4 0.4^(**)0.40.4{ }^{*} 0.40.40.4 30 0.756 0.78 0.969
5 0.4 0.4 0.4 0.4 0.4^(**)0.40.4{ }^{*} 0.40.40.4 0.4 0.4 0.4 0.4 0.4^(**)0.40.4{ }^{*} 0.40.40.4 30 0.822 0.82 1.002
5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 45 1.12 1.1 1.018
6 0.6 0.6 0.6 0.6 0.6**0.60.6 * 0.60.60.6 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 45 1.181 1.22 0.96
6 0.6 0.6 0.6 0.6 0.6**0.60.6 * 0.60.60.6 0.6 0.6 0.6 0.6 0.6**0.60.6 * 0.60.60.6 45 1.023 0.99 1.03
7 0.7 0.7 0.7 0.7 0.7**0.70.7 * 0.70.70.7 0.5 0.5 0.5 0.5 0.5**0.50.5 * 0.50.50.5 60 1.833 1.82 1.007
7 0.7*0.7 0.6*0.6 60 1.536 1.55 0.99
Mean St.deviation COV% 1
0.02
2
Span (m) Col cross section (m) Beam cross section (m) Building height (m) Fundamental time period FEM (sec) Proposed Equation (sec) FEM/proposed 4 0.4*0.4 0.3**0.3 30 0.99 0.95 1.04 4 0.4^(**)0.4 0.4^(**)0.4 30 0.775 0.77 1.006 4 0.5**0.5 0.4^(**)0.4 30 0.756 0.78 0.969 5 0.4^(**)0.4 0.4^(**)0.4 30 0.822 0.82 1.002 5 0.5**0.5 0.5**0.5 45 1.12 1.1 1.018 6 0.6**0.6 0.5**0.5 45 1.181 1.22 0.96 6 0.6**0.6 0.6**0.6 45 1.023 0.99 1.03 7 0.7**0.7 0.5**0.5 60 1.833 1.82 1.007 7 0.7*0.7 0.6*0.6 60 1.536 1.55 0.99 Mean St.deviation COV% 1 0.02 2 | Span (m) | Col cross section (m) | Beam cross section (m) | Building height (m) | Fundamental time period | | | | :--- | :--- | :--- | :--- | :--- | :--- | :--- | | | | | | FEM (sec) | Proposed Equation (sec) | FEM/proposed | | 4 | 0.4*0.4 | $0.3 * 0.3$ | 30 | 0.99 | 0.95 | 1.04 | | 4 | $0.4{ }^{*} 0.4$ | $0.4{ }^{*} 0.4$ | 30 | 0.775 | 0.77 | 1.006 | | 4 | $0.5 * 0.5$ | $0.4{ }^{*} 0.4$ | 30 | 0.756 | 0.78 | 0.969 | | 5 | $0.4{ }^{*} 0.4$ | $0.4{ }^{*} 0.4$ | 30 | 0.822 | 0.82 | 1.002 | | 5 | $0.5 * 0.5$ | $0.5 * 0.5$ | 45 | 1.12 | 1.1 | 1.018 | | 6 | $0.6 * 0.6$ | $0.5 * 0.5$ | 45 | 1.181 | 1.22 | 0.96 | | 6 | $0.6 * 0.6$ | $0.6 * 0.6$ | 45 | 1.023 | 0.99 | 1.03 | | 7 | $0.7 * 0.7$ | $0.5 * 0.5$ | 60 | 1.833 | 1.82 | 1.007 | | 7 | 0.7*0.7 | 0.6*0.6 | 60 | 1.536 | 1.55 | 0.99 | | | Mean St.deviation COV% | | | | 1 | | | | | | 0.02 | | | | | | 2 | |
Table 12: Comparison of the predicted fundamental periods using the proposed model and numerical modal analysis.

LIMITATIONS AND SCOPE OF APPLICABILITY

Although the proposed analytical model is of very high accuracy when used to estimate the fundamental time period of common RC moment-resisting frames, certain limitations of methodology should be taken into account in order to conclude its applicability range. First, the model has been calibrated and validated for regular configurations with evenly distributed stiffness and mass distribution. Its applicability to irregular buildings like those with setbacks, reentrant corners, torsional eccentricities, or mixed lateral systems (shear walls, braced frames) has not been proven and may vary from the results presented unless re-calibrated. This aligns with findings of Suthar and Purohit [40], who reported erratic deviations in the plan-irregular RC MRF durations from conventional height-based expected formulas, and proposed customized formulations with irregularity cases.
The FEM validation presumed fixed-base conditions, neglecting soil-structure interaction (SSI). Flexible foundation conditions, particularly in soft soils, can increase the period, alter lateral displacements, and modify damping characteristics. Camayang et al. [41] emphasized that SSI effects are typically underestimated in analytical equations, while they have a great influence on seismic response, and recommended site-specific assessments where foundation flexibility is expected.
Finally, the calibration process employed a linear-elastic FEM approach, which, although suitable for first-order approximation, cannot capture material nonlinearity, cracking, stiffness reduction, or energy dissipation processes that occur under strong seismic shaking. Mohamed et al. [17] demonstrated that cracking and nonlinear effects can significantly alter the fundamental period of RC moment-resisting frames compared to predictions in linear terms, which poses an imperative to consider these effects in model application. These limitations need to be overcome when the proposed model is applied in practice, and further study is encouraged to extend its application to irregular configurations, SSI conditions, and nonlinear structural responses.

Conclusions

It has been observed that many forms of fundamental vibration period equations with varying period predictions have been found in the literature, ranging from complex to simple. Because of the significant scatter and the structural differences among the available equations, most of which are based on regression analysis, this paper tries to develop a more reliable analytical formula incorporating both mass and stiffness parameters. The results show that the proposed equation provides good estimates of the fundamental period of MRFs. In the case of models analyzed herein, the maximum and average errors in the fundamental period obtained by using the proposed formula compared to those derived from eigenvalue analysis of the numerical model are no more than 5 % 5 % 5%5 \%5% and 6 % 6 % 6%6 \%6%, respectively. Thus, this equation gives an effective closed-form solution for estimating the natural vibration period. The proposed equation demonstrates very good agreement with the FEM results, as evidenced by high R 2 R 2 R^(2)\mathrm{R}^{2}R2 values and root mean square error (RMSE). These statistical metrics are important in assessing the predictive accuracy and general performance of the proposed model since they measure the extent of deviation and strength of correlation between predicted results and reference values.
The following significant conclusions are reached from the present study:
  • The current study proposes a simplified analytical approach to estimate the fundamental vibration period of reinforced concrete moment-resisting frame (RC-MRF) buildings. By idealizing the frame as an equivalent single-degree-of-freedom (SDOF) system, dynamic analysis is considerably simpler.
  • It was found that the first two storeys contribute predominantly to the overall lateral stiffness, which enabled the derivation of a simplified expression for equivalent frame rigidity. In addition, the study proposed a practical formula to compute the building's equivalent seismic mass, supported by a correction factor to account for differences in cross-sectional properties between beams and columns.
  • To enhance the usability of the mass estimation process, a curve was formulated to enable the usage of the proposed equation without sacrificing any accuracy. The final proposed model was validated through a comprehensive sensitivity study that demonstrated near correspondence with finite element method (FEM) outcomes. The model also demonstrated versatility across a range of structural configurations, including building height, span length, member stiffness ratio, and material elasticity.
  • Overall, the proposed formula is a robust and efficient method for seismic preliminary assessment and could be employed as a foundation for further detailed analytical studies or code applications. The proposed equation predicts the fundamental time period with good accuracy a mean FEM/proposed value of 1.0 and with small scatter (COV 2%).
The small variations due to span, section dimensions, and height have negligible influences, demonstrating the robustness of the equation. The proposed equation is of higher predictive accuracy compared to existing equations, and therefore, it is a reliable technique for the estimation of natural vibration periods. The technique has also been tested to be applicable in high-rise buildings up to 30 stories. The formula takes into account seismic mass and the strength of concrete in the considered direction. Vertical frame irregularities are also taken into account. In future work, refining the suggested formulas will be developed to incorporate the effect of the number of bays in high-rise buildings and soil-structure interaction parameter variations.

References

[1] UBC. (1997). Uniform building code: International Conference of Building Officials, Whittier, Calif.
[2] SEAOC. (1996). Structural Engineers Association of California , Recommended Lateral Force Requirements and Commentary, Seismological Committee, Structural Engineers Association of California, San Francisco, Calif.
[3] EGC. 2012. Egyptian Code for Computation of Loads and Forces in Structural and Building Work. Housing and Building Research Center, Cairo, Egypt.
[4] Goel, R. K. and Chopra, A. K. (1996). Evaluation of code formulas for fundamental period of buildings, 11th World Conf. on Earthquake Eng.
[5] Goel, R. K. and Chopra, A. K. (1997). Period Formulas for Moment-Resisting Frame Buildings, J. Struct. Eng., 123(11), pp. 1454-1461. DOI: https://doi.org/10.1061/(ASCE)0733-9445(1997)123:11(1454).
[6] Salama, M. (2015). Estimation of period of vibration for concrete moment-resisting frame buildings, HBRC Journal, 11(1), pp. 16-21. DOI: https://doi.org/10.1016/j.hbrcj.2014.01.006.
[7] Kaplan, O., Guney, Y., Dogangun, A. (2021). A period-height relationship for newly constructed mid-rise reinforced concrete buildings in Turkey," Engineering Structures, 232, 111807.
[8] TEC. (2018). Turkish Earthquake Code. Ministry of Environment and Urbanization, Ankara, Turkey.
[9] Alguhane, M., Fayed, N., Hussin, A., Ismail, M. (2016). Simplified equations for estimating the period of vibration of KSA existing building," Journal of Multidisciplinary Engineering Science and Technology, 3(3), pp. 4335-4343.
[10] Navarro, M., Vidal, F., Enomoto, T., Alcalá, F.J., García-Jerez, A., Sánchez, F.J. and Abeki, N. (2007). Analysis of the weightiness of site effects on reinforced concrete (RC) building seismic behaviour: The Adra town example (SE Spain). Earthquake engineering and structural dynamics, 36(10), pp.1363-1383. DOI: https://doi.org/10.1002/eqe. 685
[11] Gallipoli, M.R., Mucciarelli, M., Šket-Motnikar, B., Zupanćić, P., Gosar, A., Prevolnik, S., Herak, M., Stipčević, J., Herak, D., Milutinović, Z. and Olumćeva, T. (2010). Empirical estimates of dynamic parameters on a large set of European buildings. Bulletin of earthquake engineering, 8(3), pp.593-607.
[12] Guler, K., Yuksel, E. and Kocak, A. (2008). Estimation of the fundamental vibration period of existing RC buildings in Turkey utilizing ambient vibration records. Journal of Earthquake Engineering, 12(S2), pp.140-150.
DOI: https://doi.org/10.1080/13632460802013909.
[13] Michel, C., Guéguen, P., Lestuzzi, P., Bard, P.Y. (2010). Comparison between seismic vulnerability models and experimental dynamic properties of existing buildings in France. Bulletin of Earthquake Engineering, 8(6), pp.12951307.
[14] Oliveira, C.S. and Navarro, M. (2010). Fundamental periods of vibration of RC buildings in Portugal from in-situ experimental and numerical techniques. Bulletin of Earthquake Engineering, 8(3), pp.609-642. DOI: https://doi.org/10.1007/s10518-009-9162-1.
[15] Ditommaso, R., Vona, M., Gallipoli, M.R., Mucciarelli, M. (2013). Evaluation and considerations about fundamental periods of damaged reinforced concrete buildings. Natural Hazards and Earth System Sciences, 13(7), pp.1903-1912.
[16] Pan, T.C., Goh, K.S. and Megawati, K., 2014. Empirical relationships between natural vibration period and height of buildings in Singapore. Earthquake engineering and structural dynamics, 43(3), pp.449-465. DOI: https://doi.org/10.1002/eqe.2356.
[17] Mohamed, A.N., El Kashif, K.F., Salem, H.M. (2019). An investigation of the fundamental period of vibration for moment resisting concrete frames. Civil Engineering Journal, 5(12), pp.2626-2642.
[18] ASCE 7-16. (2016). American Society of Civil Engineers: Minimum Design Loads for Buildings and Other Structures, 7th ed. Reston, VA: American Society of Civil Engineers. DOI: https://doi.org/10.1061/9780784404454.
[19] Young, K. and Adeli, H. (2016). Fundamental period of irregular eccentrically braced tall steel frame structures. Journal of Constructional Steel Research, 120, pp.199-205. DOI: https://doi.org/10.1016/j.jcsr.2016.01.001.
[20] Young, K. and Adeli, H. (2014). Fundamental period of irregular moment-resisting steel frame structures," Structural Design Tall Build, 23(15), pp. 1141-1157. DOI: https://doi.org/10.1002/tal.1112.
[21] Sarkar, P., Prasad, A.M., Menon, D. (2010). Vertical geometric irregularity in stepped building frames. Engineering Structures, 32(8), pp. 2175-2182. DOI: https://doi.org/10.1016/j.engstruct.2010.03.020
[22] Ricci, P., Verderame, G.M., Manfredi, G. (2011). Analytical investigation of elastic period of infilled RC MRF buildings. Engineering structures, 33(2), pp.308-319. DOI: https://doi.org/ 10.1016/j.engstruct.2010.10.009.
[23] Aninthaneni, P.K. and Dhakal, R.P. (2016). Prediction of fundamental period of regular frame buildings. Bulletin of the New Zealand society for earthquake engineering, 49(2), pp.175-189.
[24] Alrudaini, T. (2023). Estimating vibration period of reinforced concrete moment resisting frame buildings, Res. Eng. Struct. Mater, 9, pp. 1417-1432.
[25] Kumar, P., Gogineni, A., Kumar, A., Modi, P. (2025). A comparative analysis of machine learning algorithms for predicting fundamental periods in reinforced concrete frame buildings. Iranian Journal of Science and Technology, Transactions of Civil Engineering, 49(3), 2257-2276. DOI: https://doi.org/ 10.1007/s40996-024-01560-0.
[26] Rahman, T., Momin, M. F., Provasha, A. A. (2025). Comprehensive analysis of structural parameters influencing the fundamental period of steel-braced RC buildings using machine learning interpretability. AI in Civil Engineering, 4(1), 7. DOI: https://doi.org/ 10.1007/s43503-025-00051-z.
[27] Shan, J., Huang, C., Wang, L., Loong, C. N. (2024). Data-driven prediction of natural period for existing RC high-rise buildings using probabilistic machine learning methods. Journal of Building Engineering, 90, p. 109394. DOI: https://doi.org/10.1016/j.jobe.2024.109394.
[28] Karampinis, I., Morfidis, K., Iliadis, L. (2024). Derivation of Analytical Equations for the Fundamental Period of Framed Structures Using Machine Learning and SHAP Values. Applied Sciences, 14(19), 9.9072.
DOI: https://doi.org/10.3390/app14199072.
[29] ATC. (1978). Applied Technological Council: Tentative provisions for the development of seismic regulation for buildings. Rep. No, ATC3-06, Applied Technological Council, Palo Alto.
[30] NEHRP. (1994). Recommended Provisions for the development of Seismic Regulations for New Buildings, Building Seismic Safety Council, Washington, D. C, 1. Federal Emergency Management Agency.
[31] NEHRP. (2000). Recommended Provisions for the development of Seismic Regulations for New Buildings, Building Seismic Safety Council, Washington, D. C, vol. 1. Federal Emergency Management Agency.
[32] ASCE 7-05. (2005). American Society of Civil Engineers: Minimum Design Loads for Buildings and Other Structures, 7th ed. Reston, VA: American Society of Civil Engineers. DOI: https://doi.org/10.1061/9780784404454.
[33] Chopra, (2001). Dynamics of Structures: Theory and Applications to Earthquake Engineering. Prentice-Hall Inc, Upper Saddle River, N.J.
[34] Asteris, P. G., Repapis, C. C., Repapi, E. V., Cavaleri, L. (2017). Fundamental period of infilled reinforced concrete frame structures. Structure and Infrastructure Engineering, 13(7), pp. 929-941.
DOI: https://doi.org/10.1080/15732479.2016.1227341.
[35] Saatcioglu M. and Humar, J. (2003). Dynamic analysis of buildings for earthquake-resistant design, Can. J. Civ. Eng., 30(2), pp. 338-359, DOI: https://doi.org/10.1139/102-108.
[36] BSLJ. (1987). Construction center of Japan .building standard law of Japan. Tokyo: construction center of Japan.
[37] NBCC. (1995). National Building Code of Canada. National Research Council of Canada, Canadian Commission on Building and Fire Codes, National Research Council of Canada. DOI: https://doi.org/10.4224/40001245.
[38] ETABS. (2018)., Structural Analysis Program (ETABS) 2018, Version 18.0.2, Computers and Structures Inc., USA.
[39] MATLAB. (2023). The MathWorks, Inc.: 2023[Online]
Avialable at: http://https://de. mathworks.com/de/help/ident/ref/goodnessoffit.html.
[40] Suthar, J. and Purohit, S. (2024). Fundamental period equations for plan irregular moment-resisting frame buildings. Advances in Civil and Architectural Engineering, 15(28), pp. 15-32. DOI: https://doi.org/10.13167/2024.28.2.
[41] Camayang, J. P. P., Padilla, J. A., De La Cruz, A. R., Bersamina, J. P. L. (2025). Soil-structure interaction in seismic response and structural performance: A systematic review. Advances in Structural Engineering. DOI: https://doi.org/10.1177/13694332251353607.