Research Article  Open Access
Free and Forced Vibration Analysis of an Infilled Steel Frame: Experimental, Numerical, and Analytical Methods
Abstract
Structural frames with masonry infill panels make up a significant portion of the buildings constructed in earthquakeprone areas prior to the developing of the seismic design standards. In this paper, the effects of masonry panels on the vibration response of an infilled steelframe building are investigated. Various ambient and steady state forced vibration tests are carried out to realize the dynamic characteristics of the system. 3D finite element models of the building with and without infill panels are provided based on marcomodeling theorem. A set of analytical approximate formulas are also derived to estimate the vibrational period. The natural frequencies of the building are computed using numerical, analytical, and experimental methods. The results show that neglecting the effect of infill panels leads to considerable error. Moreover, it is shown that there is good agreement among the results obtained by the three methods considering the effect of infill panels.
1. Introduction
Steelframe buildings with masonry infill walls are a type of construction widely used throughout many developing countries. Hollow clay tile blocks, hollow concrete blocks, and normal bricks are used in infill walls (Figure 1(a)). The infill walls, being traditionally nonengineered, have asbuilt properties which at the design stage are almost impossible to estimate reliably and/or to specify and at the construction stage are hard to control [1]. However, the effects of infill walls on the structural properties of a building have been recognized by engineers and studied for a long time [2]. The infill walls are supposed to increase the building stiffness. They also may have some undesirable effects on the building performance, such as enhancing the soft storey mechanism or causing short column effects. Chaker and Cherifati [3] suggested that plane stress finite elements provide a better representation of the inplane initial stiffness of the infill panels under the small strain condition. In order to simulate the behavior of the infill wall, different types of models can be used [4]:(i)micromodeling in which the effect of mortar joints are considered as a discrete element. In this approach, the brick and mortar are modeled as continuum elements and interface between the brick and mortar is modeled by an interface joint element (Figure 1(b)).(ii)mesomodeling in which the bricks are modeled by continuum elements, but the mortar joint and its interface with bricks are modeled together as an interface element (Figure 1(c)).(iii)macromodeling in which a single numerical model represents the infill panel effect [2]. Macromodeling is considered in two levels, that is, (1) homogenized model in which the brick, mortar, and the interface are modeled as one continuum element (Figure 1(d)), and (2) strut model in which the infill panel is modeled by one (Figure 1(e)) or more struts (Figure 1(f)) in each direction.
(a)
(b)
(c)
(d)
(e)
(f)
Micromodeling is relatively timeconsuming especially for analysis of large structures [5]. Different mechanisms have been used for the micro and mesomodeling, for example, nonlinear orthotropic model [6], smeared crack model [7], and dilatant interface constitutive model [8]. The infill panel is macromechanically modeled as homogeneous material and the average properties are used for the whole panel [9]. Modeling the infill panel as an equivalent strut is the most conventional technique. Table 1 summarizes different strut modeling approaches [4]. In addition, all the macro, meso and micromodels for simulating the presence of the infill panels are affected by a strong sensitivity to the parameters of the model. Sensitivity and (aleatory and/or epistemic) uncertainty analysis of infilled frames were investigated by Panagiotakos and Fardis [10], Uva et al. [11], Celarec et al. [12], and Celarec and DolĹˇek [13].

On the other hand the vibrational behavior of the frame structures was investigated by different researchers. Foutch studied the vibrational characteristics of a ninestory reinforcement concrete frame [30]. Jain studied the Imperial County Services Building analytically and compared the response with those recorded by strong motion accelerographs during the 1979 Imperial Valley earthquake [31]. However, the fundamental frequencies from the numerical results and the field measurements were not in close agreement because of the foundation flexibility. Jennings et al. [32] tested a forced vibration experiment on a twentytwo story steelframed structure and extracted first six natural frequencies. Foutch [33] compared the natural frequencies and 3D mode shapes of a twenty story steel structure using forced vibration test and finite element method. He concluded that the properties of the primarily translational modes are predicted reasonably well, but adequate prediction of torsional motions is not obtained. Trifunac [34] compared ambient and forced vibration techniques for testing fullscale structures under small level excitations. He concluded that both methods provide mutually consistent results.
Beside the different experimental and numerical techniques that are used for free and forced vibration analysis of structures, the dynamic behavior of the buildings and towers can be estimated by analytical methods. Free vibration of beamlike structures was studied by several researchers [35â€“37]. Most of the previous studies have been oriented towards the closedform solutions of the free flexural vibration of a beam. Li et al. [38] derived the exact solution for free flexural vibration of the stepped cantilever beams subjected to axial loads. They used the exact solution of a prismatic beam accompanied by the transfer matrix method. However, effects of the soilstructure interaction (SSI), shear deformation, and the joint flexibility were not considered. Lin and Chang [39] investigated free vibration analysis of a multispan beam with arbitrary number of flexible constraints. They assumed that each span of the continuous beam obeys Timoshenko beam theory. FirouzAbadi et al. [40] studied the transverse free vibrations of typical truncated nonuniform EulerBernoulli beams by using a WKB global approximation. Sina et al. [41] proposed a new beam theory for free vibration analysis of functionally graded beams. The beam properties were assumed to be varied through the thickness following a power law distribution. Also, it was assumed that the lateral normal stress of the beam is zero. Carrera et al. [42] proposed a hierarchical finite element method based on the Carrera Unified Formulation for free vibration analysis of beam with arbitrary section geometries. Ghasemzadeh et al. [43] proposed a set of approximate formulas to determine the natural frequencies of the structures considering the panel zone flexibility and SSI effects. AkgĂ¶z and Civalek [44] studied vibration response of the nonhomogenous and nonuniform microbeams in conjunction with BernoulliEuler beam and modified couple stress theorem. Li et al. [45] derived a closed form solution for free vibration analysis of axially inhomogeneous beams with different end conditions.
In this paper, an old steel building in which the outer frames have been filled by masonry materials is selected as case study. Various steady state forced vibration tests as well as ambient vibration test are carried out for realize the frameâ€™s dynamic characteristics. Also, a set of experiments are used to find out the compressive and shear strength of the masonry infill panels. Some approximate formulas are proposed for free vibration analysis of the steel frame, in which the structure is idealized as prismatic cantilever flexuralshear beam. In addition, 3D finite element model of the building is developed and the macromodeling approach is used for simulation of the infill panels. Finally, the natural frequencies resulting from different methods are compared.
2. Description of the Case Study Building
The selected structure is a 36yearold fivestory steelframe building located in Tehran, Iran. The outer frames of the building are filled by brick and masonry materials (Figure 2(a)). The dimensions are 40.70â€‰m Ă— 14.70â€‰m. The building is designed according to the Iranian National Building Code (INBC) [46] for a very high seismic zone. The building is located on soil type III (where the average shear wave velocity to a depth of 30â€‰m would be 180â€‰m/s to 360â€‰m/s). The foundation is strip footing with the depth of 80â€‰cm and the width of 150â€‰cm. The structural system is moment resisting frame in the EW direction and concentrically braced frame in the NS direction. A rigid diaphragm can be assumed according to the roof system used in this building. Youngâ€™s modulus, E, and the Poissonâ€™s ratio, , of the material for the main steel frame are 2.1E5â€‰MPa and 0.3, respectively. Dead load equal to DL = 650â€‰kg/m^{2} for all stories except the roof level (DL = 600â€‰kg/m^{2}) and live load equal to LL = 350â€‰kg/m^{2} for all stories except the second floor (LL = 750â€‰kg/m^{2}) are considered based on the asbuilt drawings. The typical plan view of the building for all the stories is shown in Figure 2(b).
(a)
(b)
3. Experimental Tests
3.1. Instruments
The generated vibrations in structural system should be recorded using the appropriate instruments for further analyses. Instruments used for measurement in the present study are data logger and accelerometers (Figure 3). EDX1500A is a sixteencanal data collector which is able to measure variable frequencies up to 10â€‰KHz (Figure 3(a)). Furthermore, it includes a powerful monitoring system and can be used for revision and process of collected data. EDX1500A includes four signal conditioner cards in which one of them is VAQ60A and the others are CDV60A. Each card is used for connecting four accelerometers. CDV60A card is used just for miniature accelerometers and VAQ60A card is used just for connection of the server accelerometers. It is noteworthy that all the accelerometers should be calibrated before using in the final measurements. For this purpose VAQtest is required for ensuring the proper performance of the accelerometers.
(a)
(b)
There are various methods for shaking of the structure based on purpose of the research. One of the simplest and inexpensive methods is based on the rotation of out of center masses using a motor. For this purpose a shaker was designed by the authors (Figure 4(b)). The shaker basically includes two same masspoints rotating in opposite the directions with the same velocity. Layout and theory of this method are shown in Figure 4(a). As it is clear, the resultant force in this method has a sinusoidal nature just in one direction.
(a)
(b)
3.2. Steady State Forced Vibration Test
In this test, shaker is installed on the roof of the building approximately near the center of stiffness. In addition, uniaxial accelerometers are installed in all the stories near the center of the stiffness. By shaking the structure in NS direction and analyzing the recorded data, the natural frequency of building is measured as â€‰Hz. By the same method for EW direction, â€‰Hz.
The shape of the entire frequencyresponse wave is controlled by the amount of damping in the system; so it is possible to derive the damping ratio from many different properties of the curve. The most conventional method is bandwidth or halfpower method, in which the damping ratio is determined from the frequencies at which the response is 1/ times the peak response. With this method, the damping ratio is computed as about 6% for this structure.
3.3. Ambient Vibration Test
Considering that all sides of the selected building are free, it is possible to apply ambient vibration test without using any external shaker. For this purpose, data loggers are installed near the center of the stiffness in each story in NS and EW directions and results are recorded on a stormy day. Resulting frequencies are the same as those obtained from forced vibration test except that the wind is used as a natural shaker in the ambient vibration test.
4. Numerical Model
In this section finite element model of the structure is prepared and the natural frequencies are extracted. For this purpose all the properties of the infill panels should be extracted based on the standard tests.
4.1. Material Properties of Infill Panels
All the infill panels should be in the well condition without any cracking. Then, based on Iranian rehabilitation guideline (2007), chapter 76112, shear strength of the mortar should be evaluated as follows [46]: â€śAt least one test is required per 300â€‰m^{2} of the infill panels. At least 8 tests are required for each building to obtain the reliable results. Moreover, the dispersion of the test locations should be selected in a way that covers all parts of the buildingâ€ť. Implementation of this experiment is shown in Figure 5. The mortar average bedjoint shear strength is calculated as where is the load magnitude due to first movement of the masonry unit, is the stress due to gravity loads at the test location, and is total net area of the bed joints above and below the test unit. The results of the tests and corresponding calculations are summarized in Table 2. Based on this table, the average shear stress in the mortar is 2.25â€‰kg/cm^{2}; however, the acceptable value for the mortar shear stress should be calculated using the probabilistic method, chapter 76112, in a way that 80% of the test results should have greater value of the final shear stress [47]. Thus, the final value of the mortar shear stress is 1.93â€‰kg/cm^{2} in the present case. Considering that the final value is less than 2.0â€‰kg/cm^{2}, the minimum requirements for the mortar strength are not satisfied and there is need for rehabilitation and strengthens.

Based on the conducted tests and considering Iranian rehabilitation guideline [46] and FEMA356 [48] recommendations for calculation of the compressive strength of the mortar, the basic parameters are as follows: brick compressive and tensile strength, and , are 38.5â€‰kg/cm^{2} and 5.5â€‰kg/cm^{2}, respectively, and the mortar compressive strength, , is 10.0â€‰kg/cm^{2}. In addition, the thickness of the brick and the mortar are â€‰cm and â€‰cm, respectively. Thus, the compressive strength can be calculated as The resulting is 26.66â€‰kg/cm^{2}. Based on the FEMA356 [48], the factor to translate lowerbound masonry properties to the expected strength masonry properties is 1.3. Also, it recommends a factor to translate compressive strength to the modules of elasticity. The final value for the compressive strength and the modulus of elasticity are 34.658â€‰kg/cm^{2} and 19061â€‰kg/cm^{2}, respectively.
4.2. Developing Macromodels
As mentioned before, there are various methods in order to simulate the infill panelsâ€™ effects. In the present study, both the macromodeling based on continuum elements (Figure 1(d)), and the compressivestrut element (Figure 1(e)) are used. All the beams and columns are modeled using â€śBeamâ€ť elements while all the beamcolumn connections are assumed to be rigid. SSI effect is considered by modeling the foundation and its surrounding environment as massless medium using â€śSolidâ€ť elements. The massless foundation is extended about three times of the building height in all directions. It should be noticed that there is no need to apply the nonreflecting boundary conditions on the farend boundaries of the foundation model, because no seismic analysis is performed in this paper. In the macromodeling based on the continuum elements, the infill panel is directly modeled by 3D â€śShellâ€ť elements considering appropriate properties of the main panel. In the macromodeling based on equivalent strut approach, the infill panel is replaced with a diagonal compressivestrut element of width, , and the same thickness and modulus of elasticity of infill panel using â€śLinkâ€ť elements. This is a twonode pinnedend element which connects two opposite corners of the frame. The element is capable of transferring only compressive pressure in the axial direction. In this method, the equivalent width, , for the elastic inplane behavior prior to cracking can be calculated based on Figure 6 [46]: where in which is the column height between centerlines of beams, is the height of infill panel, is the expected modulus of elasticity of frame material, is the expected modulus of elasticity of infill material, is the moment of inertia of column, is the diagonal length of infill panel, is the thickness of infill panel and equivalent strut, is the angle whose tangent is the infill heighttolength aspect ratio, and is the coefficient used to determine the equivalent width of the strut.
Consequently, three different models are provided in this study for simulation of the building.(i)The structure is modeled without infill panels effects (reference case).(ii)The structure is modeled considering the infill panels as a compressive struts.(iii)The structure is modeled considering the infill panels as 3D shell elements.
5. Analytical Approximate Solution
5.1. Fundamentals and Assumptions
In order to derive a set of appropriate approximate formulas for analytical vibration analysis of the structure, a doubly symmetric structure in plan is selected as a sample which is subjected to uniformly distributed gravity loads at story levels as shown in Figure 7(a). It is also assumed that both the beams and the columns have uniform sections throughout the height of the building. Furthermore, floor slabs are considered as rigid diaphragms in their own plane so that the relative displacements between frames are restricted.
(a)
(b)
The structure is idealized as a prismatic cantilever beam with flexural rigidity , shear rigidity , axial distributed compression force , and mass per unit length , as shown in Figure 7(b). Floor masses and gravity loads at story levels are also replaced by concentrated masses and concentrated forces , respectively. SSI is modeled using the axial and torsional springs. The general algorithm for the analytical solution is summarized in Figure 8 and will be explained in detail in the next subsections.
5.2. Concept of Energy Method
5.2.1. Kinetic and Potential Energy
Assuming that the frame is vibrating freely, the kinetic energy is given by [43]: in which is transverse displacement of the frame, dots indicate the differentiation with respect to time , is the total number of stories, and is the height of the structure. The transverse displacement of the structure is expressed as where and are contribution of beam (frame) and springs (soil) in the lateral displacement, respectively, is portion of the displacement induced by bending deformation and is portion of the displacement induced by the shear deformation of the beam. is transverse displacement due to foundation sliding and is transverse displacement due to foundation rotation.
The potential energy of the system, , consists of the strain energy, , and the work done by external loads, : where where is the shear shape factor of the crosssection, considering the effects of uneven distribution of shear deformation over the crosssection, is the stiffness coefficient against sliding of the foundation, is the stiffness coefficient against rotation of the foundation, and is rotation of the torsional spring. Consider in which represents the axial shortening due to lateral deformation and is expressed as
5.2.2. Energy Equation
Assuming that the structure is vibrating freely in simple harmonic motion, the transverse displacement and the velocity of vibration are expressed as where is the natural vibration frequency. The maximum potential energy of the system over a vibration cycle is associated with the maximum displacement; that is, and also the maximum kinetic energy of the system over a vibration cycle is associated with maximum velocity: The principle of energy conservation states that the total energy in a freely vibrating system without damping is constant; thus, must be equal to which yields where As seen, (14) is of the same form as the frequency expression for a single degree of freedom (SDOF) system. In other words, by restricting the lateral displacement of the structure to a single shape function, , which defines the mode shape, the structure behaves as a generalized SDOF system [47]. Parameters and are generalized mass and generalized stiffness of the system, respectively, and is generalized geometric stiffness of the system which considers the effect of axial loads on the natural frequency. Assuming that all the stories have equal height of , all floor slabs have equal mass of , and gravity loads at story levels are equal to . Thus, (15) can be simplified as
5.2.3. Proper Shape Function
The accuracy of frequency of vibration depends on the assumed shape function, . The more accurate the shape function is, the more accurate the result would be. It is worth mentioning that is not a particular function, and any shape function satisfying the geometric boundary conditions can be used. By using the deflected shape of the equivalent beam due to a selected set of forces, the geometric boundary conditions are automatically satisfied. Several methods are available to determine the deflected shape of the structure, but arguably, the most promising one is the principle of virtual work: where and are the real shear and moment functions along the length of the structure, respectively, and and are the virtual force functions in the structure arising from the application of unit load. and are the real shear and moment at the base of the structure, respectively, and finally and are the virtual forces at the base of the structure from the application of unit load, respectively.
5.3. Calculation of the Natural Frequency
Considering the SSI effects, the deflected shape of the structure, , under lateral uniform distributed load can be written as where The subscript, , represents the fixity of the lowermost elevation and subscript, , represents the effects of soil flexibility. Using (14) for extracting the natural frequencies in the present case, new parameters can be derived as follows: Different terms are defined as follows in which the subscripts and represent the portions associated with distributed axial load and concentrated axial forces, respectively:
5.4. Properties of Equivalent Beam
5.4.1. Flexural Rigidity
The modulus of elasticity, , for the equivalent beam should be the same as the modulus of elasticity of the real structure. The moment of inertia, , of the equivalent beam is in which is the crosssectional area of the th column, is the distance of the th column from axis, and is the distance of the th column from axis.
5.4.2. Shear Rigidity
(1) Moment Resisting Frame. In order to calculate the shear rigidity per unit web frame, , a subassemblage is extracted from the frame assuming that inflection points occur at midspan of the beams on either side of the joint and at midheight of the columns above and below the joint. This subassemblage is shown in Figure 9. It is further assumed that the beams on either side of the joint are of the same section and length and that the columns above and below the joint are the same section and length. is obtained by equalizing the displacement of the subassemblage to displacement of a shear element. Considering completely rigid joints, the shear rigidity is computed as follows: where in which is the height of column (story height), is the length of the beam (bay width), and are depth of the beam and depth of the column, respectively, and are moment of inertia of the beam and the column, respectively, and and are crosssectional area of the beam and the column, respectively.
(2) Braced Frame. The equivalent shear rigidity of double bracings is where and are the crosssectional area and Youngâ€™s modulus of the brace, respectively. For a single and tensiononly bracing, the coefficient 2 in (25) must be replaced with 1.
(3) Infilled Frame. In order to consider the effect of infill panel, shear rigidity of the infill panel should be added to the equivalent shear rigidity of the system:
(4) Dual System. The equivalent shear rigidity of the dual system is
6. Results
In this section, the results obtained from numerical simulations, experimental tests, analytical approximate formulas, and the codified proposed values for natural frequencies are compared. Table 3 summarizes all the methods used in this study and also the corresponding abbreviations. Based on this table, two experimental methods, six numerical models with and without infill panel effect on rigid/flexible foundation, four analytical approximate formulations, and finally three codifiedbased methods are used in order to estimate the natural frequencies.

Figure 10 shows the natural frequencies of the considered building in NS and EW directions using various methods. As mentioned before, the results of the natural frequencies based on forced vibration and ambient vibration tests have the same values and are 2.70â€‰Hz in NS and 2.50â€‰Hz in EW directions. Based on the numerical models, it can be concluded that neglecting infill panel effect leads to decreasing the frequency meaningfully. Also, modeling the foundation as a rigid medium leads to increasing the frequency a little in all models. Modeling infill panel as a compressive strut increases frequencies in both directions. The differences between N3 and T1 are about 11% for NS direction and 10% in EW direction. Macromodeling based on shell elements increases the natural frequencies than to macromodeling based on the compressivestrut theorem. In this condition, the differences between N5 and the experimental tests are 5.6% and 6% for NS and EW directions, respectively. As seen, modeling the infill panel based on the continuum element theorem leads to realistic behavior than to modeling it as a compressive strut.
Table 4 summarizes the percentage of the error between the numerical, approximate formulas, and the codified based methods with the exact values obtained from the experimental tests. Using the proposed approximate formulas neglecting the infill panels leads to almost the same results as those obtained from numerical simulations. On the other hand, using the proposed technique considering both infill panel and SSI effects leads to very close results to the experimental test. In this condition, the percentage of the errors between F4 and T1 are only 1.5% and 2.4% for NS and EW directions, respectively. Once again, it should be mentioned that the proposed formulation does not account for the torsional deformation effect. That means that the term is not incorporated in the formulas. Therefore, when a building is unsymmetrical, the error of the formula is meaningful. The value of the error depends to the distance between the center of mass and the center of rigidity.

Finally, the results of the experimental tests are compared with codified proposed formulas. All three design codes are proposed an empirical formulation in the form of , in which is the height of the structure and and are the constant values which differ by the type of the lateral resisting system. Considering that the lateral resisting system in the current case is a combination of moment resisting and braced and infill panel systems, the value of is 0.05 based on INBC and 0.0488 based on UBC and ASCE [48] codes. The value of the is 0.75 based on all three codes (all values are presented in SI system). Using the codifiedbased method leads to the same frequency for both directions of the building because this method does not directly account for the direction. Also, it seems that using UBC and ASCE methods leads to a little better result in comparison with INBC.
Figure 11 shows the first mode shape of the considered building in two orthogonal directions using experimental test and numerical models. These mode shapes are extracted based on the relative lateral displacements of the structure in numerical models and also values of the resonance test which describe average peak of Fourier spectrum in different stories. As seen, in both directions N5 model has the closest mode shape to the actual mode shape of building and N1 model has the least similarity to the real mode shape.
(a)
(b)
7. Conclusions
This paper presents vibration analysis of a steelframe building considering the infill panelsâ€™ effect through experimental, numerical, and analytical approaches. In order to extract the actual natural frequencies of the building various steady state forced and ambient vibration experimental tests were carried out. It was found that the results of two methods meaning the forced and ambient vibration are the same. In addition, a set of appropriate tests were carried out to realize the compressive and shear strength of the masonry materials of the infill panels. It was observed that the materials are not in good condition and this old building needs to be rehabilitated.
A set of approximate formulas were proposed for free vibration analysis of steel structures with various structural systems, where the frame is idealized as prismatic cantilever flexuralshear beam. In the proposed technique, soilstructure interaction is considered as axial and torsional springs whose potential energy is formulated and incorporated into overall potential energy of the structure.
Moreover, a set of threedimensional finite element models of the building were provided. In these models main structural systems are modeled by â€śbeamâ€ť elements while the infill panels are simulated macromechanically by either continuum elements or equivalent struts.
It was found that modeling infill panels has significant effect on vibration characteristics of the structures and neglecting their effects can leads to errors about 35%â€“45% in natural frequencies. Considering soilstructure interaction leads to more flexibility in the system and decreases frequency of system. Generally, modeling infill panels by continuum elements lead to higher frequencies and decrease the percentage of the errors between the finite element model and experimental tests. It was concluded that using detailed finite element model of the structure, modeling infill panels using shell elements, and also considering the soilstructureinteraction may leads to more accurate result.
It was shown that the proposed approximate formulas have good capability in estimation of the natural frequency of steel structures. They can be used for computing the frequency of various structural systems by less computational efforts. In the present study, there is only 1.5%â€“2.5% error when using the approximate formulas compared to the forced vibration test. Also, it was shown that the codified proposed formulas have a close estimation to the test results and lead to errors of about 2%â€“12% in considered building.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 M. N. Fardis, S. N. Bousias, G. Franchioni, and T. B. Panagiotakos, â€śSeismic response and design of RC structures with planeccentric masonry infills,â€ť Earthquake Engineering and Structural Dynamics, vol. 28, no. 2, pp. 173â€“191, 1999. View at: Google Scholar
 K. A. Ghassan, Nonductile behavior of reinforced concrete frames with masonry infill panels subjected to inplane loading [Ph.D. thesis], University of Illinois at Chicago, USA, Chicago, Ill, USA, 1998.
 A. A. Chaker and A. Cherifati, â€śInfluence of masonry infill panels on the vibration and stiffness characteristics of RC frame buildings,â€ť Earthquake Engineering & Structural Dynamics, vol. 28, no. 9, pp. 1061â€“1065, 1999. View at: Google Scholar
 S. Sattar, Influence of masonry infill walls and other building characteristics on seismic collapse of concrete frame buildings [Ph.D. thesis], University of Colorado, Boulder, Colo, USA, 2013.
 P. B. Shing, H. R. Lofti, A. Barzegarmehrabi, and J. Bunner, â€śFinite element analysis of shear resistance of masonry wall panels with and without confining frames,â€ť in Proceedings of the 10th World Conference on Earthquake Engineering, pp. 2581â€“2586, A. A. Balkema, Rotterdam, The Netherlands, 1992. View at: Google Scholar
 M. Dhanasekar and A. W. Page, â€śInfluence of brick masonry infill properties on the behaviour of infilled frames,â€ť Proceedings of the Institution of Civil Engineers, vol. 81, no. 2, pp. 593â€“605, 1986. View at: Publisher Site  Google Scholar
 T. C. Liauw and C. Q. Lo, â€śMultibay infilled frames without shear connectors,â€ť ACI Structural Journal, vol. 85, no. 4, pp. 423â€“428, 1988. View at: Google Scholar
 A. B. Mehrabi and P. B. Shing, â€śFinite element modeling of masonryinfilled RC frames,â€ť Journal of Structural Engineering, vol. 123, no. 5, pp. 604â€“613, 1997. View at: Publisher Site  Google Scholar
 H. R. Lotfi and P. B. Shing, â€śAn appraisal of smeared crack models for masonry shear wall analysis,â€ť Computers and Structures, vol. 41, no. 3, pp. 413â€“425, 1991. View at: Publisher Site  Google Scholar
 T. B. Panagiotakos and M. N. Fardis, â€śSeismic response of infilled RC frames structures,â€ť in Proceedings of the 11th World Conference on Earthquake Engineering, Paper No. 225, Acapulco, Mexico, 1996. View at: Google Scholar
 G. Uva, F. Porco, and A. Fiore, â€śAppraisal of masonry infill walls effect in the seismic response of RC framed buildings: a case study,â€ť Engineering Structures, vol. 34, no. 1, pp. 514â€“526, 2012. View at: Publisher Site  Google Scholar
 D. Celarec, P. Ricci, and M. Dolšek, â€śThe sensitivity of seismic response parameters to the uncertain modelling variables of masonryinfilled reinforced concrete frames,â€ť Engineering Structures, vol. 35, pp. 165â€“177, 2012. View at: Publisher Site  Google Scholar
 D. Celarec and M. Dolšek, â€śThe impact of modelling uncertainties on the seismic performance assessment of reinforced concrete frame buildings,â€ť Engineering Structures, vol. 52, pp. 340â€“354, 2013. View at: Publisher Site  Google Scholar
 M. Holmes, â€śSteel frame with brickwork and concrete infilling,â€ť ICE Proceedings, vol. 19, pp. 473â€“478, 1961. View at: Google Scholar
 B. StaffordSmith, â€śLateral stiffness of infilled frames,â€ť Journal of Structural Division, vol. 88, pp. 183â€“199, 1962. View at: Google Scholar
 B. StaffordSmith and C. Carter, â€śA method of analysis for infilled frames,â€ť ICE Proceedings, vol. 44, pp. 31â€“48, 1969. View at: Google Scholar
 R. J. Mainstone and G. A. Weeks, â€śThe influence of a bounding frame on the racking stiffness and strengths of brick walls,â€ť in Proceedings of the 2nd International Brick Masonry Conference (SIBMAC '70), Building Research Station, England, UK, 1970. View at: Google Scholar
 R. J. Mainstone, â€śOn the stiffness and strengths of infilled frames,â€ť ICE Proceedings, vol. 49, no. 2, p. 230, 1971. View at: Publisher Site  Google Scholar
 R. E. Klingner and V. V. Bertero, â€śEarthquake resistance of infilled frames,â€ť Journal of the Structural Division, vol. 104, no. 6, pp. 973â€“989, 1978. View at: Google Scholar
 L. TeChang and K. KwokHung, â€śNonlinear behaviour of nonintegral infilled frames,â€ť Computers and Structures, vol. 18, no. 3, pp. 551â€“560, 1984. View at: Publisher Site  Google Scholar
 C. A. Syrmakezis and V. Y. Vratsanou, â€śInfluence of infill walls to RC frames Response,â€ť in Proceedings of the 8 th European Conference on Earthquake Engineering (EAEE '86), pp. 47â€“53, European Association for Earthquake Engineering, Istanbul, Turkey, 1986. View at: Google Scholar
 R. Zarnic and M. Tomazevic, â€śAn experimentally obtained method for evaluation of the behavior of masonry infilled RC frames,â€ť in Proceedings of the 9th World Conference on Earthquake Engineering, pp. 163â€“168, 1988. View at: Google Scholar
 T. Schmidt, â€śAn approach of modelling masonry infilled frames by the F.E. method and a modified equivalent strut method, Darmstadt Concrete,â€ť Annual Journal on Concrete and Concrete Structures, pp. 185â€“194, 1989. View at: Google Scholar
 C. Z. Chrysostomou, Effects of degrading infill walls on the nonlinear seismic response of twodimensional steel frames [Ph.D. thesis], Cornell University Press, Ithaca, NY, USA, 1991.
 A. Saneinejad and B. Hobbs, â€śInelastic design of infilled frames,â€ť Journal of Structural Engineering, vol. 121, no. 4, pp. 634â€“650, 1995. View at: Publisher Site  Google Scholar
 F. J. Crisafulli, Seismic behavior of reinforced concrete structures with masonry infills [Ph.D. thesis], University of Canterbury, Christchurch, New Zealand, 1997.
 R. D. Flanagan and R. M. Bennett, â€śArching of masonry infilled frames: comparison of analytical methods,â€ť Practice Periodical on Structural Design and Construction, vol. 4, no. 3, pp. 105â€“110, 1999. View at: Publisher Site  Google Scholar
 W. W. ElDakhakhni, M. Elgaaly, and A. A. Hamid, â€śThreestrut model for concrete masonryinfilled steel frames,â€ť Journal of Structural Engineering, vol. 129, no. 2, pp. 177â€“185, 2003. View at: Publisher Site  Google Scholar
 M. Dolšek and P. Fajfar, â€śThe effect of masonry infills on the seismic response of a fourstorey reinforced concrete frame—a deterministic assessment,â€ť Engineering Structures, vol. 30, no. 7, pp. 1991â€“2001, 2008. View at: Publisher Site  Google Scholar
 D. A. Foutch, Study of the vibration characteristics of two multistory building [Ph.D. thesis], California Institute of Technology, Pasadena, Calif, USA, 1977.
 S. K. Jain, â€śContinuum models for dynamics of buildings,â€ť Journal of Engineering Mechanics, vol. 110, no. 12, pp. 1713â€“1730, 1984. View at: Publisher Site  Google Scholar
 P. C. Jennings, R. B. Matthiesen, and J. Brent Hoerner, â€śForced vibration of a tall steelframe building,â€ť Earthquake Engineering and Structural Dynamics, vol. 1, pp. 107â€“132, 1972. View at: Publisher Site  Google Scholar
 D. S. A. Foutch, â€śThe vibrational characteristics of a twelve storey steel frame building,â€ť Earthquake Engineering and Structural Dynamics, vol. 6, no. 3, pp. 265â€“294, 1978. View at: Google Scholar
 M. D. Trifunac, â€śComparisons between ambient and forced vibration experiments,â€ť Earthquake Engineering and Structural Dynamics, vol. 1, no. 2, pp. 133â€“150, 1972. View at: Publisher Site  Google Scholar
 N. M. Auciello and G. Nolè, â€śVibrations of a cantilever tapered beam with varying section properties and carrying a mass at the free end,â€ť Journal of Sound and Vibration, vol. 214, no. 1, pp. 105â€“118, 1998. View at: Publisher Site  Google Scholar
 J. Wu and C. Chen, â€śAn exact solution for the natural frequencies and mode shapes of an immersed elastically restrained wedge beam carrying an eccentric tip mass with mass moment of inertia,â€ť Journal of Sound and Vibration, vol. 286, no. 3, pp. 549â€“568, 2005. View at: Publisher Site  Google Scholar
 D.W. Chen and J.S. Wu, â€śThe exact solutions for the natural frequencies and mode shapes of nonuniform beams with multiple springmass systems,â€ť Journal of Sound and Vibration, vol. 255, no. 2, pp. 299â€“322, 2003. View at: Publisher Site  Google Scholar
 Q. S. Li, H. Cao, and G. Li, â€śAnalysis of free vibrations of tall buildings,â€ť Journal of Engineering Mechanics, vol. 120, no. 9, pp. 1861â€“1876, 1994. View at: Publisher Site  Google Scholar
 H. Lin and S. C. Chang, â€śFree vibration analysis of multispan beams with intermediate flexible constraints,â€ť Journal of Sound and Vibration, vol. 281, no. 12, pp. 155â€“169, 2005. View at: Publisher Site  Google Scholar
 R. D. FirouzAbadi, H. Haddadpour, and A. B. Novinzadeh, â€śAn asymptotic solution to transverse free vibrations of variablesection beams,â€ť Journal of Sound and Vibration, vol. 304, no. 3–5, pp. 530â€“540, 2007. View at: Publisher Site  Google Scholar
 S. A. Sina, H. M. Navazi, and H. Haddadpour, â€śAn analytical method for free vibration analysis of functionally graded beams,â€ť Materials and Design, vol. 30, no. 3, pp. 741â€“747, 2009. View at: Publisher Site  Google Scholar
 E. Carrera, M. Petrolo, and P. Nali, â€śUnified formulation applied to free vibrations finite element analysis of beams with arbitrary section,â€ť Shock and Vibration, vol. 18, no. 3, pp. 485â€“502, 2011. View at: Publisher Site  Google Scholar
 H. Ghasemzadeh, H. RahmaniSamani, and M. Mirtaheri, â€śVibration analysis of steel structures including the effect of panel zone flexibility based on the energy method,â€ť Earthquake Engineering and Engineering Vibration, vol. 12, pp. 587â€“598, 2013. View at: Publisher Site  Google Scholar
 B. Akgöz and Ö. Civalek, â€śFree vibration analysis of axially functionally graded tapered BernoulliEuler microbeams based on the modified couple stress theory,â€ť Composite Structures, vol. 98, pp. 314â€“322, 2013. View at: Publisher Site  Google Scholar
 X.F. Li, Y.A. Kang, and J.X. Wu, â€śExact frequency equations of free vibration of exponentially functionally graded beams,â€ť Applied Acoustics, vol. 74, no. 3, pp. 413â€“420, 2013. View at: Publisher Site  Google Scholar
 Instruction for Seismic Rehabilitation of Existing Buildings, Code No. 360, Management and Planning Organization (Office of Deputy for Technical Affairs), 2007.
 A. K. Chopra, Dynamics of Structures, Prentice Hall, Englewood Cliffs, NJ, USA, 1995.
 ASCEFEMA, â€śPrestandard and commentary for the seismic rehabilitation of buildings,â€ť FEMA 356, 2000. View at: Google Scholar
Copyright
Copyright © 2014 Mohammad Amin HaririArdebili et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.