Next Article in Journal
Transparent Localized Haptics: Utilization of PVDF Actuators on Touch Displays
Next Article in Special Issue
Design and Experimental Assessment of a Vibration Control System Driven by Low Inertia Hydrostatic Magnetorheological Actuators for Heavy Equipment
Previous Article in Journal
Linear and Nonlinear Models for Drop Simulation of an Aircraft Landing Gear System with MR Dampers
Previous Article in Special Issue
Evolutionary Computation-Based Active Mass Damper Implementation for Vibration Mitigation in Slender Structures Using a Low-Cost Processor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Structural Optimization of Lightweight Composite Floors with Integrated Constrained Layer Damping for Vibration Control

by
Carlos M. C. Renedo
1,
Iván M. Díaz
1,*,
Jaime H. García-Palacios
2 and
Christian Gallegos-Calderón
1
1
Department of Continuum Mechanics and Theory of Structures, ETSI Caminos, Canales y Puertos, Universidad Politécnica de Madrid, Calle Profesor Aranguren 3, 28040 Madrid, Spain
2
Department of Hydraulics, Energy and Environmental Engineering, ETSI Caminos, Canales y Puertos, Universidad Politécnica de Madrid, Calle Profesor Aranguren 3, 28040 Madrid, Spain
*
Author to whom correspondence should be addressed.
Actuators 2023, 12(7), 288; https://doi.org/10.3390/act12070288
Submission received: 5 June 2023 / Revised: 3 July 2023 / Accepted: 10 July 2023 / Published: 13 July 2023

Abstract

:
Due to current architectural trends, contemporary public buildings are becoming open-plan spaces with much less weight and damping. Consequently, Vibration Serviceability Limit State (VSLS) due to human-induced vibrations has become an increasing concern for structural engineers, especially when designing offices, hospitals, or gymnasiums. When dealing with resonant vibrations, a slight increase in the floor-damping enables decreasing considerably the vibration level. The damping strategy studied in this work is usually known in the literature as Constrained Layer Damping (CLD) and consists of a viscoelastic layer constrained between the concrete slab and the steel beam of a lightweight composite floor. In this paper, a complete structural checking methodology has been developed for analyzing all the limit states that determine the final sizing of a steel–concrete composite floor treated with CLD, including a detailed analysis of the VSLS. The methodology has been used for setting a structural optimization problem for floors with and without CLD treatments. Thus, it has been demonstrated that the integration of CLD treatments at the design stage of the building allows the development of lighter floor structures with a smaller embodied carbon (EC) footprint, especially for long-span schemes with restrictive vibration limitations.

1. Introduction

Currently, composite floors are widely used in multistory commercial and office buildings due to their capability to span longer distances with stronger, stiffer, and lighter structural solutions, which are economical and offer a fast construction process. Their use often implies the reduction of the primary structure and the foundations of the building, thus reducing its overall embodied carbon (EC) [1]. Over recent decades, the increasing use of large open-plan spaces and the rise of the electronic office has led to the reduction of dead and live loads (1 kPa less) [2,3] due to the removal of full-height partitions and heavy furnishing elements. As a consequence, modern composite floors have lower damping ratios (around 1% to 3% [4,5,6,7]), lower mass, and lower fundamental natural frequencies (often below 10 Hz), making them more prone to vibrate under human-induced dynamic loading, such as walking in offices or rhythmic activities in gymnasiums [8]. Furthermore, restrictive vibration comfort levels (0.02–0.04 m/s2) are presently required for floors in calm environments such as residences, offices, and hospitals. In this context, the Vibration Serviceability Limit State (VSLS) has become a sizing criterion to consider when designing composite floors [9].
Excessive vibrations are often related to low-frequency floors (LFF) (with natural frequencies below 10 Hz) due to their resonance with the 2nd, 3rd, or even 4th harmonic of the human dynamic loading. In high-frequency floors (HFF) (above 10 Hz) it is considered that the vibration does not have a resonant nature and is made of consecutive transient decaying responses of the floor to each human footfall [5]. Knowing this, structural designers often try to include stiffer-enough steel members (that raise the natural frequencies of the floor to avoid resonance) or heavy concrete slabs that decrease the amplitude of the vibration [10]. This practice leads to a significant increase in the material used on the floor (especially relevant when the vibration limits to comply are very restrictive) affecting its lightness and increasing its EC as has been proved by Gonçalves and Pavic in [9]. This research studies the structural oversizing to be performed in a composite floor of 30 m × 44 m with bays of 6 m × 10 m when the VSLS to be met becomes more restrictive (changing from a limiting root mean square (rms) acceleration of 0.04 m/s2 to 0.02 m/s2). They concluded that the floor weight should be increased by 27% and its EC by 14% to accomplish the new comfort level. Alternatively, the integration of damping strategies into the floor’s design may be an effective way of improving their dynamic performance without increasing their structural mass, thus, minimizing their EC. This is a key point considering that floors are responsible for around 70% of the total EC of the superstructure of a building [11,12].
Two main types of damping strategies may be used to solve floor vibration problems. First, the use of inertial dampers, in their passive versions (most known as Tuned Mass Dampers or TMDs) [13], semi-active [14] and active versions [15] that apply counteract forces in real time to reduce the structural motion. Second, the use of dissipative strategies that enable increasing the inherent damping of the floor by wisely including specific elements that dissipate additional energy, such as viscous or Viscoelastic (VE) dampers. This paper studies one damping treatment related to this second group: Constrained Layer Damping (CLD). It consists of a thin VE layer embedded between two elastic bending members. The VE layer suffers shear deformation when the constraining elastic members vibrate in bending modes, hence additional energy is dissipated through a stress-strain shear hysteresis [16].
CLD treatments were developed and have been extensively applied in aerospace and mechanical engineering to mitigate broad-band frequency vibrations [17]. Their use in civil engineering is still limited, although, some studies have been carried out using them to improve the dynamic performance of composite floors. First, Nelson [18] (Figure 1a) and Farah et al. [19] (Figure 1b) applied a VE CLD treatment to the lower steel flanges of a lively floor, obtaining damping ratio increases from 3.5% to 5%. Ebrahimpour [20] (Figure 1c) confirmed these results with a similar treatment. Later, Ahmadi et al. [21] (Figure 1d) developed a new CLD retrofitting system applied to the concrete slab of the floor with similar damping ratio increases. In 2006, ARUP in collaboration with Richard Lee Steel Decking proposed the first commercial CLD solution to be integrated into a composite floor since the construction stage. This was called ‘Resotec’ (Figure 1e), and its applicability was described by Willford et al. [22] (Figure 2). This type of CLD treatment is the one studied in this paper. The major advance of this proposal over the previous ones was the location of the VE layer between the rib-deck concrete slab and the upper flange of the steel beam (as close as possible to the composite section centroid) where the shear strain to be achieved is maximum. Moreover, they recommended including the VE layer only for a percentage of the beam’s length near the supports, where the longitudinal shear strain is higher, whereas the central part of the composite beam remains connected to the shear using studs. Therefore, there is a trade-off between the achieved inherent damping and the loss of stiffness due to the partial shear connection in length. The CLD treatment itself used by ARUP has an overall thickness of 3 mm and is composed of two thin steel sheets of 1 mm that constrain a slim VE layer of 1 mm. Unlike previous applications, this CLD treatment was not conceived to be retrofitting solution but an integrated technology that needed to be considered when designing the composite floor. The optimal design of this CLD treatment and its integration into the design workflow of a composite floor were aspects not covered by Willford et al. [22].
Thus, this paper aims to provide a methodology for the structural checking and structural optimization of composite floors prone to vibrate with this integrated CLD treatment. The main objective is to prove that the use of this CLD treatment enables overcoming the VSLS without adding great amounts of additional mass to the floor and minimizing its EC.
The structural optimization of composite floors has been researched by many authors using different optimization algorithms, Objective Functions (OFs), and problem constraints (which are the limits states (LSs) to be met by the floor). Initially, the authors used the weight and the cost of the floor as OFs. Zahn [23], optimized the weight of a set of composite beams using two different codes. The VSLS was considered according to an early recommendation given by Murray in 1981 [24] based on a heel drop loading. He concluded that vibrations were the sizing criterion for spans between 3 m and 6 m. Later, Kim and Adeli [25,26] optimized the economic cost of composite floors with a simple OF considering the influences of concrete, steel, and studs. They affirmed that the VSLS was checked but without providing details. The cost optimization performed by Klanšek and Kravanja [27,28,29] used a complex OF accounting for material, energy consumption, and labor cost items required to manufacture a composite floor. They studied the use of different types of steel members for typical ranges of spans and loads, according to the Eurocodes but without considering the VSLS. The research developed by Kaveh and Ahangaran [30,31] followed a similar methodology as the one proposed by Zahn but they focused on testing the efficiency of novel meta-heuristic algorithms for structural optimization. Poitras et al. [32] was the first author to include a widely accepted methodology to check the VSLS (the AISC Steel Design Guide 11 Floor vibrations due to human activity [7]) for the optimization process of two composite floor examples. His work was used as a reference for another research performed by Kaveh and Ghafari [33] where it was concluded that the vibrations were one of the critical LSs when designing a floor bay of 8 m × 10 m. Yossef and Taher [34] performed another cost optimization of composite floors with castellated beams using a widely accepted methodology to check the VSLS (the SCI guideline Design of floors for vibration [5]). In recent years, with the rise of life-cycle assessment, many authors have started to include the EC (measured in kgCO2eq) of the floor as an OF to optimize. Roynon [35] has presented a manual optimization of framed buildings minimizing their EC and giving reference values of kgCO2eq/m2 for steel framed floors (between 80–250 kgCO2eq/m2 for regular spans and 150–300 kgCO2eq/m2 for long-span floors). Drevniok et al. [36] have developed an optimization methodology called ‘The Lightest Beam Method’ to study the amount of material that could have been reduced in a set of already designed and built steel floors without composite action. They consider the VSLS as a limitation in the natural frequency of 3 Hz, without assessing in detail the level of vibration. They concluded that the relaxation of serviceability requirements at the design stage could enable a reduction by more than 30% of the floors’ weight. Some authors as Whitworth and Tsavdaridis [37] and Kravanja et al. [38] have used the embodied energy (measured in J) to quantify the environmental impact of the floors. Finally, in the research performed by Gauch et al. [39] different structural types of floors (made of concrete, steel, and timber) are designed and compared in terms of cost and EC obtaining interesting conclusions especially applicable to a decision-making stage of the design process. They again consider the VSLS by imposing a limit on the floor fundamental frequency of 4 Hz.
As a conclusion, it is clear that the VSLS should be properly taken into account in terms of the floor’s acceleration response according to currently accepted guidelines as [5,7] or [6]. Thus, to the author’s knowledge, this paper is the first one addressing the multi-objective optimization (with the floor’s Weight and its EC as Ofs) of composite floors designed with and without integrated CLD treatments, and including a detailed verification of the VSLS.
The remainder of this paper is organized as follows: Section 2 describes the elements of the type of composite floors to be optimized. Section 3 outlines all the static limit states to be checked along the floor design. Section 4 explains the workflow used to check the VSLS and how the CLD treatment has been optimally designed for each analyzed case. Section 5 portrays the multi-objective structural optimization carried out describing the optimization procedure, the Ofs, and the optimization parameters used. Section 6 describes a parametric study in which different square floor bays with different spans have been optimized with and without the use of CLD treatments to illustrate the proposed methodology. Finally, Section 7 provides some conclusions.

2. Description of a Composite Floor with Integrated CLD Treatment

A composite floor with a given total length L and a total width B can be composed of different bays as represented in Figure 3. A bay is defined as a floor part between 4 columns and is composed of primary (1ry) beams spanning between columns and secondary (2ry) beams that span between these first ones. There are 3 types of bays depending on their location, (i) corner bays, (ii) edge bays, and (iii) internal bays. The first two groups are usually more prone to vibrate as they cannot mobilize as much mass when vibrating as the internal bays.
In this paper, the design of a composite floor with integrated CLD treatments is carried out according to the floor bay depicted in Figure 4. A given floor bay is geometrically defined by the length of its 1ry and 2ry beams, here known as  L 1  and  L 2 , respectively; the number of 2ry beams  N 2  and the distance between them  d 2 . The steel members of the floor are here defined through profile sections from the Universal Beam (UB) series (this catalog was chosen as it is large enough so that a meaningful optimization can be carried out). The profiles are represented with the integer values  P 1  and  P 2  which correspond with the profile positions in the UB catalog sorted according to the moment of inertia. S-275 steel has been considered. The concrete slab of the floor is defined by two parameters, (i) the commercial rib-deck used for the floor, which is defined with the integer value  R d  (according to its position on the list of Cofraplus 60 rib-decks, with 4 different gauge thicknesses of 0.7 mm; 0.9 mm; 1 mm and 1.2 mm), and, (ii) the concrete slab thickness over the ribs in centimeters represented with the integer value  h c . The concrete grade used is C-30. Steel studs have been included along the central region of the 1ry and 2ry beams to ensure a certain degree of shear connection between the steel members and the concrete slab. In the 1ry beams, the separation of studs depends on the required bending resistance of the composite section at mid-span. In the 2ry beams their separation is equal to the distance between the valleys of the rib-deck. Ductile shear studs made of steel with a diameter of 19 mm, an as-welded height of 95 mm, and an ultimate limit stress of  f u  = 450 Mpa, have been considered.
The CLD treatment has been included in all beams along a given length,  L v , i , near the supports (being “i” an index used to refer the 1ry or 2ry beam). Thus, the percentages of beam” length treated with CLD are obtained as  % C L D i = 100 λ i , where  λ i  is obtained as follows:
λ i = ( 2   L v , i ) / L i .
The thickness of the VE layer of these CLD treatments,  h v , i , needs to be small so three possible values of 0.5 mm; 1.0 mm and 2 mm have been considered for this parameter based on [40]. Moreover, the width of these VE layers,  b v , i , should be similar to the width of the top flange of their correspondent beams,  b f , i . For that reason,  b v , i  must take a value within the continuous range that goes from  0.5   b f , i  to  2   b f , i . The properties assigned to the VE material of the CLD treatment are in the range of commercial VE materials for vibration damping (for example, the HIP2, before cited [40]) i.e., a storage shear modulus  v  of 0.5 mPa, a loss factor  η v  with a value of 1 and a density of 1700 kg/m3. The loss factor of this type of material is usually measured based on the standard testing method ASTM:E756-05 [41].
In Table 1 the main material properties used along the paper are listed indicating, the name of the material, the element in which it is used, and its main mechanical properties.  ρ  is used to denote the density of the material,  E s t  for the Young modulus of the steel,  E c  for the Young modulus of the concrete,  f y  is the yielding strength of the steel, and  f u  is the ultimate tensile strength of the steel.
In this paper, for given a floor bay defined by a set of  L 1  and  L 2 , the rest of the floor parameters above described ( P 1 P 2 N 2 h c R d L v , 1 L v , 2 h v , 1 h v , 2 b v , 1  and  b v , 2 ) are optimized to minimize two oFs, while meeting all the LSs of the floor design.

3. Static LSs to Be Met by the Floor

3.1. Static LSs in the 1ry and 2ry Beams

All the beams of the floor are designed as simply supported and built without any propping system. Hence, two design stages have been verified: (i) the construction and (ii) the service life of the floor. The LSs checked in the 1ry and 2ry beams are summarized in Table 2. Each LS is defined through a safety factor  S F , which, if higher than 1, means that LS is met. To identify each  S F , a set of three subscripts separated by commas has been used:
  • Subscript 1 identifies the variable used to compute the safety factor. The following options are available: ‘ M + ’ means the sagging bending moment at mid-span, ‘ M + B ’ means the sagging bending moment at the location of the first shear stud in partially treated beams, ‘V’ denotes shear force at the support, and ‘ δ L L ’ means deflection under live load action.
  • Subscript 2 may be ‘c’ or ‘s’ depending on if the LS belongs to the construction stage or to the service-life stage, respectively.
  • Subscript 3 is the subscript ‘i’ used to denote the 1ry or 2ry beam.

3.1.1. Static LSs of 1ry and 2ry Beams: The Construction

The loads considered here are the self-weight,  S W c  and a construction live load  L L c  = 0.75 kN/m2. At this stage, the concrete is not a resisting element and it has a fresh density of 2500 kg/m3.
The verification has been done according to Eurocode 3 [42], as the only resisting element is the steel member of the floor. The Ultimate LS (ULS) of bending at mid-span has been checked with respect to the plastic bending moment of the steel section,  M p , s t , and the ULS of shear at the support section with respect to the plastic shear of the steel web,  V p , s t  (with the subscript ‘p’ indicating plastic and ‘st’ indicating steel).

3.1.2. Static LSs of 1ry and 2ry Beams: The Service Life

The loads considered at this stage are the self-weight  S W s  of the structural elements, (but now using a dried density for the concrete of 2400 kg/m3), a deal load  D L s  for flooring finishes of 1 kN/m2 and a live load  L L s  of 3 kN/m2.
The verification has been performed following Eurocode 4 [43], as the concrete slab is now a resisting element of the beam section. Four LSs have been checked at this stage, (i) the ULS of bending at the end of the CLD treatment, where the fist shear stud is located (ii) the ULS of bending at mid-span, (iii) the ULS of shear, verified only considering the resisting contribution of the web of the steel member, and (iv) the deflection serviceability LS (DSLS) at mid-span under the action of  L L s , which has been limited to  L i / 350 . Special attention must be paid to the verification of the ULS of bending and the DSLS of composite beams partially treated with CLD along their length.

3.1.3. ULS of Bending for Composite Beams Partially Treated with CLD

In a regular simply supported composite beam, the ULS of bending is verified by assuring that the resisting moment of the composite section at mid-span  M R d , c o m p  is higher than the service life acting design bending moment  M E d , s . This criterion is not sufficient for partially treated beams as the most critical section for bending may lie somewhere in between the first section connected to the shear and the mid-span section. This can be seen in Figure 5, where the bending resisting envelopes of a 15 m secondary beam treated with different  % C L D  have been compared with the laws of maximum design bending moment. It can be noted that along the CLD-treated length, the resisting bending moment is the one provided by the steel member,  M R d , s t , whereas, along the shear-connected central length its value increases due to the composite action. A simplified methodology to avoid checking the ULS of bending along the whole beam length has been developed based on checking two critical sections, (i) the mid-span and (ii) the first shear-connected section, here called Section B. A similar simplification was performed by Willford et al. [22]; However, in this contribution, it was assumed a degree of shear connection along the central connected region of 60% and it did not provide a detailed explanation of how this ULS should be computed. A more accurate description of this verification is provided below with adequate assumptions.
For checking the Section B the design bending moment at this section has been compared to the plastic bending moment of the steel member,  M p , s t .
To verify the mid-span section the following considerations have been made:
  • The effective breadth of the slab,  b e f f , i , contributing to the composite section has been computed as follows:
    b e f f , 1 = min ( L 1 2 L v 1 ) / 4 L 2 ,
    b e f f , 2 = min ( L 2 2 L v 2 ) / 4 d 2 .
    The expression used here is the same one as for regular simply supported composite beams, but replacing the beam length by the shear-connected length,  L c o n n , i = ( L i 2 L v , i ) . This decision is proposed by the authors as it is on the safe side.
  • Two ductile shear studs per section connected have been used to achieve a higher degree of shear connection in a shorter distance. Additionally, plastic redistribution of longitudinal shear forces between the different shear studs along the connected length has been assumed.
  • Partial shear connection theory has been assumed. Hence, the maximum longitudinal shear to be transferred within the steel–concrete interface,  N c  is limited by the plastic capacity of the shear connection. This means that the maximum longitudinal shear when considering full shear connection, named  N c f , needs to be reduced by a factor  α  called ’degree of shear connection’ which is computed as follows:
    α = P r d ( n / 2 ) N c f ,
    N c = α N c f ,
    where  P r d  is the shear resisting value of a single shear stud and n is the total number of shear studs along the shear-connected length of the beam.
    Finally, the resisting bending moment at mid-span has been computed following the methodology given in Eurocode 4 [43].

3.1.4. DSLS for Composite Beams Partially Treated with CLD

The authors propose the use of two bending stiffnesses along the beam length (see Figure 6): (i) the stiffness of the steel member for the CLD-treated region, and (ii) the bending stiffness of the composite section considering full-composite action (i.e., an infinitely rigid shear connection) and the cracking of the concrete slab. The slab” effective breadth used has been the same as the one described for the ULS given in Equations (2) and (3).

3.2. Static LSs in the Rib-Deck Slab

The rib-deck slab has been analyzed as a multi-span beam with as many supports as secondary beams are used within a floor bay. Again, two stages have been checked: (i) the construction and, (ii) the service life of the rib-deck. Table 3 summarizes the set of LSs checked in the rib-deck slab.
Again, each LS is related to a  S F  defined with a set of two subscripts:
  • Subscript 1 defines the variable used to compute the  S F : ‘ M + ’ refers to the sagging bending moment at the edge span of the slab, ‘ M ’ denotes the hogging bending moment at the internal support of the edge span, ‘V’ means the shear force at the internal support of the edge span, ‘ δ L L ’ means the maximum deflection under the live load action at the edge span, and ‘ δ S W ’ indicates the maximum deflection under the self-weight action at the edge span.
  • Subscript 2 again may be ‘c’ or ‘s’ depending on the stage analyzed.

3.2.1. Static LSs in the Rib-Deck Slab: The Construction

No propping system has been considered during the construction, hence, the only resisting element in this stage is the steel decking. The same loads considered for the beam” verification have been used here. Four LSs have been checked, (i) the ULS of the sagging bending moment in the edge span, (ii) the ULS of the hogging bending moment in the second support, (iii) the ULS of shear in the second support (for this case, the crippling of the steel decking webs has been considered to be the most restrictive shear failure mode) and, (iv) the SLS of deflection in the edge span during the concreting process limited to  d 2 / 180  to control the ponding effect. In each case, the construction live load was assumed to be present in the worst spans. The resisting parameters of the steel decking have been taken from the Cofraplus 60 technical sheets.

3.2.2. Static LSs in the Rib-Deck Slab: The Service Life

During the service life, the composite section of the rib-deck slab must be verified. The same four LSs as for the construction stage have been checked but considering the resisting parameters of the composite section: (i) the ULS of sagging bending moment. For computing the resisting moment no reinforcement has been used, and the contribution of the steel decking has been accounted through using a degree of shear connection with the concrete (as done before with the beam” shear connection). The limit shear tension within the steel–concrete interface has been assumed to be 0.1 mPa [44]. (ii) The ULS of hogging bending moment. For computing the resisting moment, the contribution of the steel decking has been neglected, thus, a regular reinforced concrete section has been analyzed. An appropriate hogging reinforcement has been sized for each case being extended over a length equal to  d 2 / 3  from the supports. To obtain both the acting sagging and hogging bending moments, a redistribution due to cracking of the bending law of a 15% has been assumed. (iii) The ULS of shear. The resisting shear of the composite slab has been obtained as a result of two contributions, the shear resistance of the steel decking and the shear resistance of the concrete ribs without any shear reinforcement. (iv) The DSLS. The deflection at the edge span under the action of the worst live load placing has been computed as the average value between those obtained with and without considering cracking in the concrete. This has been limited to  d 2 / 350 .

4. VSLS of Floors Partially Treated with CLD

The VSLS of the floor has been assessed according to the simplified methodology provided by the AISC Steel Design Guide 11 [7]. This guideline has been reported to be the most accurate one when predicting the VSLS of composite floors according to [8,45]. The AISC method quantifies the vibration in terms of peak acceleration; however, this paper provides equivalent root mean square (rms) acceleration values, as these enable obtaining the response factors widely used in other guidelines and in future versions of Eurocodes to establish vibration limits.

4.1. Modal Parameters of the Floor

4.1.1. Fundamental Natural Frequency of the Floor

The fundamental natural frequency of a floor bay,  f n , has been obtained by computing the fundamental natural frequencies of the primary and secondary beams from the maximum static deflections,  δ 1  and  δ 2 , under the floor’s weight considered for the VSLS, and then applying the Dunkerle”s approximation as follows:
f n = 18 1000 ( δ 1 + δ 2 ) .
These deflections have been computed considering a load per unit of length for each beam,  q V S L S , i , equal to:
q V S L S , i = q S W , i + q D L , i + 0.1   q L L , i ,
where  q S W , i q D L , i  and  q L L , i  are the self-weight, dead load, and live load per unit of length of the beam, respectively. Then, the following expression has been derived to compute the deflection of a partially CLD-treated beam. The authors propose to apply the Maxwell-Mohr method to the law of elastic curvatures of a simply supported beam, with two different bending stiffnesses and uniformly loaded:
δ i = L i 4 q V S L S , i 384 3 λ i 4 + 8 λ i 3 E s t   I s t , i + 3 λ i 4 8 λ i 3 + 5 E c   I c , i ,
where  E s t  and  E c  are the Young modulus of the steel and the concrete, respectively;  I s t , i  is the moment of inertia of the steel profile used for the beams, and  I c , i  is a concrete-homogenized moment of inertia of the composite section present along the shear-connected region of the beam.  λ i  is the proportion of CLD-treated beam in parts per unit computed in Equation (1). For 1ry beams just supporting one 2ry beam at mid-span,  δ 1  has been increased by a factor of 1.3.

4.1.2. Effective Weight of the Floo’s Fundamental Mode of Vibration

This magnitude has been computed following the AISC Design Guideline 11 as follows. This method is composed of 3 steps:
  • Computing the effective widths,  B e f f , 1  and  B e f f , 2 , respectively associated with the simply supported bending modes of the 1ry and 2ry beams of the floor, following an already calibrated formulation described below:
    B e f f , 1 = L 1 C 1 I e f f , 2 / d 2 I e f f , 1 / L 2 1 / 4   < L ,
    B e f f , 2 = L 2 C 2 I s l a b I e f f , 2 / d 2 1 / 4   < B ,
    where  I s l a b  is the concrete-homogenized moment of inertia of one meter of composite rib-deck slab,  C 1  and  C 2  are calibration coefficients depending on the type of composite floor beams and floor bay to be analyzed, and  I e f f , 1  and  I e f f , 2  are the effective concrete-homogenized moments of inertia of the 1ry and 2ry beams, respectively. They have been computed as follows:
    I e f f , i = L i 4   5   q V S L S , i 384   E c   δ i .
  • Computing the effective weights,  W e f f , 1  and  W e f f , 2 , respectively associated with the simply supported bending modes of the 1ry and 2ry beams:
    W e f f , 1 = B e f f , 1   L 1   q V S L S , 1 2   L 2 ,
    W e f f , 2 = K 2 B e f f , 2 L 2 q V S L S , 2 2 d 2 ,
    where  K 2  is a factor equal to 1.5 for floor bays with adjacent bays having secondary beams of length higher than 0.7  L 2 , and equal to 1 in other cases.
  • Computing the final effective weight  W e f f  of the combined mode of vibration:
    W e f f = W e f f , 1 δ 1 δ 1 + δ 2 + W e f f , 2 δ 2 δ 1 + δ 2 .

4.1.3. Intrinsic Damping Ratio of the Floo’s Fundamental Mode of Vibration

Two main sources of damping contribute to defining the final value of the fundamental modal damping ratio,  ξ n . These are: (i) the intrinsic damping of the structure,  ξ i n t  (the usual source of energy dissipation that depends on the type of floor structure, the floor finishes, and the furniture present on the floor) and, (ii) the additional damping ratio provided by the VE CLD treatment,  ξ C L D . Thus, the following equation has been used:
ξ n = ξ i n t + ξ C L D .
ξ i n t  has been computed assuming an electronic office fit-out with ceiling and ductwork. Hence, according to the AISC guideline,  ξ i n t = 0.01 + 0.01 + 0.005 = 0.025 .

4.1.4. Additional Damping Ratio Provided by the CLD Treatment

In this paper,  ξ C L D  has been computed using a simplified new methodology proposed by the author based on solving the problem of a simply supported ’sandwich beam’ partially treated with a VE core. This method is divided into 4 sub-steps.
  • Obtain  N 1  and  N 2 , the number of 1ry and 2ry beams involved in the fundamental mode of vibration, respectively. This is performed using the effective widths,  B e f f , i  computed before, assuming that they define the effective floor area contributing to the vibration when assessing the VSLS:
    N 1 = B e f f , 1 L 2 + 1 B e f f , 2 L 1 ,
    N 2 = B e f f , 2 d 2 + 1 B e f f , 1 L 2 .
  • Compute the additional damping ratio,  ξ C L D , i , provided by the applied CLD treatments to isolated 1ry and 2ry beams. This has been done by analyzing a simply supported beam partially treated with CLD and characterized by the following parameters  L i L v , i P i R d h c b e f f , i  and with a CLD treatment defined by  h v , i  and  b v , i .
    The value of  ξ C L D , i  depends on four dimensionless parameters:
    • The loss factor of the VE material  η v , the higher  η v  the better. In this paper, this parameter has been assumed to be constant and equal to 1.
    • The percentage of the beam length treated with CLD or  % C L D . The higher this parameter, the higher the damping enhancement.
    • The so-called ‘geometric parameter’ Y, computed as the ratio between the following two bending stiffness:
      Y = ( E c I c , i ) ( E s t I s t , i ) .
      where the numerator and denominator are the bending stiffnesses of a  100 %   C L D  beam where the VE core has been replaced by a shear connection with infinite or zero shear stiffness, respectively. This ratio oscillates between 1 and 3 for composite floor beams. The higher the Y, the higher the damping.
    • The so-called ’shear parameter’ of the beam, g, which represents the shear stiffness of the VE core (it should be noted that g is strongly dependent on the parameters of the VE treatment). This parameter is computed as follows:
      g = v b v , i   L i 2 h v , i E s t   A s t , i + E c   A s l a b , i E s t   A s t , i   E c   A s l a b , i ,
      where  A s t , i  and  A s l a b , i  are the areas of the profile and the slab that belong to the beam section to be analyzed, respectively. For each beam, there is an optimum  g o p t  that provides the maximum additional damping to the beam. The value of  g o p t  depends on Y η v  and  % C L D  of the beam as demonstrated in [46]. The dependency of  g o p t  with Y can be neglected in the range of Y values adopted by floor beams. In the same way, the VE materials used for this purpose have similar values of  η v  close to 1, so this dependency can be also neglected.
Knowing this, the authors have proposed a simplified methodology to obtain  ξ C L D , i , assuming that it depends only on g and  % C L D ξ C L D , i  has been obtained using the set of curves depicted in Figure 7.
These curves have been obtained for a set of 15 m composite beams with a Cofraplus 60 rib-deck slab with a width of 3 m and a height of 0.12 m, made of lightweight concrete. The steel profiles used varied depending on the  % C L D  analyzed (UB 533X210X82 for 10  % C L D , UB 610X178X82 for 2—40  % C L D , UB 610X178X100 for 5—60  % C L D , and UB 610X210X133 for 7—100  % C L D ). A detailed FE model of the beams was built using SHELL181 elements for the steel profile and SOLID185 elements for the VE treatment and the concrete slab. The applied CLD treatment had a thickness of 1 mm and a width equal to the width of the upper steel flange. The value of  η v  used was 1. Different complex modal analyses of each beam were performed varying the shear stiffness of the VE core to obtain the set of  ξ C L D , i  curves as a function of g.
Knowing this and given a partially treated floor beam, the geometry of the CLD treatment is designed by selecting the right values of  h v , i  and  b v , i  within the available ones, to provide a shear parameter as close as possible to  g o p t , as shown in Figure 7 for a 50  % C L D  beam.
3.
Compute the bending modal strain energies,  U b e n d , i  of isolated 1ry and 2ry beams, as follows:
U b e n d , i = 1 2 0 L i M V S L S , i ( s ) 2 E ( s ) I ( s ) d s ,
where  M V S L S , i ( x ) E ( x )  and  I ( x )  are the laws of bending moments (under the loads considered for the VSLS), Young modulus, and Moment of inertia along the beam length. Assuming the scheme of a simply supported beam uniformly loaded and divided into two regions of different inertia, the following expression is derived:
U b e n d , i = q V S L S , i 2   L i 5 1920 3   λ i 5 15   λ i 4 + 20   λ i 3 E s t I s t , i 3   λ i 5 15   λ i 4 + 20   λ i 3 8 E c I c , i .
4.
Calculate the final  ξ C L D  of the floor’s fundamental mode of vibration as a weighted value of  ξ C L D , 1  and  ξ C L D , 2  depending on their contribution to the total modal strain energy of the floor:
ξ C L D = ξ C L D , 1   N 1   U b e n d , 1 + ξ C L D , 2   N 2   U b e n d , 2 N 1   U e e n d , 1 + N 2   U b e n d , 2 .

4.2. Response of the Floor

For checking the VSLS, two types of dynamic floor responses have been computed under walking excitation according to the methodology given by [7]: (i) the resonant response of vibration modes with frequencies lower than 9 Hz, and (ii) the impulsive transient produced by a footfall excitation.

4.2.1. Low-Frequency Resonant Response

A harmonic punctual load given by the following expression has been used to model the walking excitation acting at the anti-node of the vibration mode:
F h ( t ) = Q   α h   s i n ( 2 π f n t ) ,
where t represents the time, Q is the weight of an average human assumed to be 700 N and  α h  is a dynamic loading factor depending on the loading frequency (assumed equal to  f n ) and given by:
α h = 0.83   e 0.35 f n .
The resonant rms acceleration of the floor is then computed according to:
a r m s , r e s = K c r e s   Q   0.83 e 0.35 f n 2   2   ξ n   W e f f / a g ,
where  K r e s  is a factor equal to 0.5 that accounts for incomplete resonant build-up from walking, and, it considers the fact that the walker and the potentially annoyed person are not simultaneously at the same location of maximum modal acceleration. Finally, a resonant response factor,  R r e s  has been computed as follows (with  a r m s , r e s  in m/s2):
R r e s = a r m s , r e s 0.005 .

4.2.2. Impulsive Response of the Floor

The consecutive impulsive transient responses produced by a floor when excited by human footfalls may also be excessive in comfort terms; thus, this vibration has also been checked. For this the effective impulse function,  I e f f  (with units of [Ns], derived by Willford and Young [6]) has been used as excitation:
I e f f = 42 f s t e p 1.43 f n 1.30 ,
where  f s t e p  is the human pacing frequency. The use of higher pacing frequencies results in higher dynamic impulsive responses, thus, a  f s t e p  = 2.6 Hz has been selected, as a conservative value. The initial peak acceleration just after the impulse is given by:
a p , i m p = K c   2 π f n I e f f W e f f / a g ,
where  K c  is a factor introduced by the AISC guideline equal to 1.3. This code also amplifies the impulsive response with a higher mode factor equal to 2, which considers the contribution of higher modes to the impulsive response. However, to the authors’ knowledge, this factor is over-conservative and, in fact, it is not used in other guidelines as [6]. From here,  a r m s , i m p , the rms acceleration of the decaying response during the time lapse from one footfall to the next one, is computed according to:
a r m s , i m p = f s t e p   0 ( 1 / f s t e p ) a p , i m p   e 2 π f n ξ n t   s i n ( 2 π f n t ) 2   d t ,
a r m s , i m p a p , i m p 2 1 e 4 π f n ξ n f s t e p 4 π f n ξ n f s t e p f o r ξ n < 0.1 .
Finally, an impulsive response factor,  R i m p  is computed as indicated below (with  a r m s , i m p  in m/s2):
R i m p = a r m s , i m p 0.005 f n < 8 Hz , R i m p = a r m s , i m p 2 π f n   0.0001 f n > 8 Hz .

4.2.3. Comparison with VSLS Limit

A limiting response factor,  R l i m , must be set depending on the vibration comfort level required by the floor use. Hence, two possible safety factors for checking the VSLS can be built: (i)  S F R m a x , s  when considering checking both the resonant and impulsive responses, and (ii)  S F R r e s , s  when only checking the resonant response of the floor:
S F R m a x , s = m a x R r e s , R i m p R l i m ,
S F R r e s , s = R r e s R l i m .

5. Optimization Problem Definition

In this section, the authors propose the definition of an optimization problem for the structural design of composite floors partially treated with CLD. The floor optimization addressed is a discrete multi-objective optimization problem with two OFs ( f 1 ( x ̲ )  and  f 2 ( x ̲ ) ) depending on a vector  x ̲  of seven design variables defined by integer numbers (from  x ̲ 1  to  x ̲ 7 ), and with six design constraints defined by the functions  g 1 ( x ̲ )  to  g 6 ( x ̲ ) . The optimization problem can be formulated as follows:
min f 1 ( x ̲ ) = M ( x ̲ ) min f 2 ( x ̲ ) = E C ( x ̲ ) s . t . g 1 ( x ̲ ) = U L S 1 ( x ̲ ) 0 g 2 ( x ̲ ) = U L S 2 ( x ̲ ) 0 g 3 ( x ̲ ) = D S L S 1 ( x ̲ ) 0 g 4 ( x ̲ ) = D S L S 2 ( x ̲ ) 0 g 5 ( x ̲ ) = L S s l a b ( x ̲ ) 0 g 6 ( x ̲ ) = V S L S ( x ̲ ) 0 6 x ̲ 1 = h c 20 1 x ̲ 2 = R d 4 1 x ̲ 3 = P 1 94 1 x ̲ 4 = P 2 94 0 x ̲ 5 = 10 λ 1 10 0 x ̲ 6 = 10 λ 2 10 1 x ̲ 7 = ( N 2 1 ) 10 x ̲ Z ,
where  M ( x ̲ )  and  E C ( x ̲ )  are the functions to compute the mass of a floor bay and its embodied upfront carbon per unit of floor area, respectively. Additionally,  U L S 1 ( x ̲ ) U L S 2 ( x ̲ ) D S L S 1 ( x ̲ ) D S L S 2 ( x ̲ ) L S s l a b ( x ̲ )  and  V S L S ( x ̲ )  are aggregate functions of the LSs checked for the different floor elements, and are described in detail in the following subsection devoted to the constraints of the problem. They provide a value of 0 if the concerning LSs have been met and positive values if not.

5.1. Design Variables

The seven design variables are integers since the structural optimization of a floor is subjected to the adoption of discrete values on its geometrical parameters to ease the construction process of the floor.
The floor optimization performed assumes a given floor bay defined by fixed distances between columns  L 1  and  L 2  that belongs to a given floor with fixed total length (L) and fixed total width (B). This bay is optimized by finding the optimal values of the following variables:  h c R d P 1 P 2 10 λ 1 10 λ 2  and  ( N 2 1 ) . The remaining floor-defining parameters are obtained when checking the different LSs as follows:
  • n i : the number of shear studs used in each beam is calculated to meet the bending ULS of the beam at mid-span.
  • b v , i  and  h v , i : the dimensions of the CLD treatment used in each floor beam are obtained to maximize the  ξ C L D , i  of each beam.
  • ρ h , s l a b : the amount of hogging reinforcement of the slab is computed to meet the ULS of the hogging bending moment of the slab.

5.2. Objective Functions

Two objective functions have been minimized: (i) The mass of the floor bay per unit of area,  M ( x ̲ ) , and (ii) the embodied upfront carbon of the floor bay per unit of area,  E C ( x ̲ ) . The choice of these two OFs is based on the fact that to tackle the vibration problem from the design perspective, two main strategies can be adopted: (i) to increase the mass of the floor, by thickening the concrete slab, which has a great impact on the floor’s  M ( x ̲ )  but not much in the  E C ( x ̲ ) , or (ii) to stiffen the steel members of the floor, which increases much more the  E C ( x ̲ )  than the  M ( x ̲ ) . This stays true when using low-strength concrete (with low cement content) and non-recycled steel, as in this case the carbon factor of the steel is much higher than the concrete one.
The computation of  E C ( x ̲ )  has been carried out according to the guideline ’How to Compute Embodied Carbon’ by Gibbons et al. [47]. This study has only considered the influence of the embodied upfront carbon of the floor (i.e., from modules A1 to A5 of the floor’s life-cycle). For the computation of  E C ( x ̲ ) , the following floor elements have been considered: concrete slab, slab reinforcement, rib-deck sheet, steel beams, shear studs, and CLD treatment. Modules A1 to A3 (those belonging to the Product Stage of the life-cycle) are the main contributors to the final upfront EC of the floor. Thus, a list of the A1–A3 Carbon Factors used (denoted as  E C F A 1 A 3 , n ) for the nth materials of the floor, is provided in Table 4.
The floor studied is supposed to be built in a UK city. The steel and concrete have been assumed to be provided by a national and a local supplier, through road transport of 300 km and 50 km, respectively. An ‘A4 Carbon Factor’ for road transport of 0.10749 gCO2eq/kg/km has been used.

5.3. Design Constrains

Six different design constraint functions have been used:  U L S 1 ( x ̲ ) U L S 2 ( x ̲ ) D S L S 1 ( x ̲ )  and  D S L S 2 ( x ̲ )  for the ULS and DSLS of 1ry and 2ry beams, respectively;  L S s l a b ( x ̲ )  for all the LSs of the rib-deck slab, and  V S L S ( x ̲ )  for the VSLS of the whole floor. These functions are built using safety factor functions  S F ( x ̲ )  (that compute for a given floor configuration  x ̲ , all the  S F  listed in previous sections) in combination with the following Modified Heaviside function  H ( S F )  (that gives a value of 0 if the  S F  value is lower than 1, and provides 1 if the input is greater or equal than 1):
H : R 0 , 1 S F H ( S F ) H ( S F ) = 1 i f S F 1 0 i f S F < 1 .
Thus, the constraint equations used are defined as follows:
U L S i ( x ̲ ) = 1 [ H ( S F M + , c , i ( x ̲ ) ) + H ( S F V , c , i ( x ̲ ) ) + H ( S F M + , s , i ( x ̲ ) ) + + H ( S F M + B , s , i ( x ̲ ) ) + H ( S F V , s , i ( x ̲ ) ) ] / 5 ,
D S L S i ( x ̲ ) = 1 H ( S F δ L L , s , i ( x ̲ ) ) ,
L S s l a b ( x ̲ ) = 1 [ H ( S F M + , c ( x ̲ ) ) + H ( S F M , c ( x ̲ ) ) + H ( S F V , c ( x ̲ ) ) + H ( S F δ S W , c ( x ̲ ) ) + + H ( S F M + , s ( x ̲ ) ) + H ( S F M , s ( x ̲ ) ) + H ( S F V , s ( x ̲ ) ) + H ( S F δ L L , s ( x ̲ ) ) ) ] / 8 ,
V S L S ( S F ( x ̲ ) ) = 1 H ( S F R m a x , s ( x ̲ ) ) 1 H ( S F R r e s , s ( x ̲ ) ) .
These functions provide a value greater than zero if not all the LSs involving them are fulfilled, and equal to zero if all the LSs are met. The function  V S L S ( x ̲ )  has two possible definitions depending on whether the impulsive acceleration wants to be limited.

5.4. Optimization Algorithm

The evolutionary optimization algorithm ‘Non-Dominated Sorting Genetic Algorithm II’, best known as NSGA-II, for constrained and discrete multi-objective optimization, has been used in this paper. The algorithm version used is available in the python library ‘pymoo’ and has been implemented according to [49]. A random sampling of integer values has been performed to generate the initial solutions. A constraint handling method based on the principle ‘Feasibility First’ has been used to avoid the definition of any penalty function. This methodology was proposed by [50] and uses a fitness function (applied to each solution) that depends on the current population. When tournament selection is applied to these fitness values, feasible solutions are always emphasized over infeasible ones. A simulated binary crossover based on [50] with a crossover index  η c  = 3, and a probability of crossover for each variable of  p c  = 0.5, has been used. A polynomial mutation has been employed with a mutation parameter  η m  = 3 and a probability of mutation  p m  = 0.5. The histograms used for crossover and mutation have been rounded to deal with integer variables. Finally, for each optimization performed 100 generations with a population of 100, have been used.

6. Parametric Study

The optimization problem proposed has been used to develop a parametric study in which many different floor bays are optimized with and without CLD. This has been built by combining four different parameters which are needed to optimize a given floor bay:
  • L i : The length of the 1ry and 2ry beams of the floor. The parametric study has been focused on analyzing square floor bays, where  L 1 = L 2 L i  has been varied from 4.5 m to 19.5 m each 1.5 m, analyzing a whole set of 11 possibilities for this parameter.
  • R l i m : The limiting response factor of the floor bay to define the VSLS. Four different possibilities have been considered for  R l i m . First, a  R l i m = , which corresponds to floor designs in which the VSLS has not been checked, here denoted as ‘Statically designed floors’.  R l i m = 8 , which is the limitation used for regular electronic offices.  R l i m = 4  is the limit used for quiet spaces like silent offices or libraries. Finally,  R l i m = 2  would apply to hospital floors.
  • S F V S L S : The  S F  used to define the VSLS. Two possibilities are contemplated,  S F R m a x  which limits the impulsive and resonant floor responses, and  S F R r e s  which only limits the resonant response.
  • C L D i n t : This indicates if the CLD treatment has been integrated into the design. Two possibilities are studied,  C L D i n t  = CLD and  C L D i n t  = NO CLD.
The following vectors provide a summary of the four parameters used and the values they can adopt:
L i = [ 4.5 ,   6.0 ,   7.5 ,   9.0 ,   10.5 ,   12.0 ,   13.5 ,   15.0 ,   16.5 ,   18.0 ,   19.5 ] , R l i m = [ 2 ,   4 ,   8 ,   ] , S F V S L S = [ S F R m a x ,   S F R r e s ] , C L D i n t = [ N O   C L D ,   C L D ] .
A total of 176 cases (resulting from combining all the different possibilities 11 × 4 × 2 × 2 = 176) have been analyzed. One particular case can be, for example, the one defined by the vector [9, 4,  S F R m a x , NO CLD], which corresponds to a 9 m × 9 m floor bay, designed to comply with a maximum response factor of 4, including the impulsive and resonant responses, and without any integrated CLD treatment.
For each optimized floor bay seven design parameters result from the optimization problem if the CLD is integrated:  h c R d P 1 P 2 λ 1 λ 2  and  ( N 2 1 ) . If the CLD is not used, this number reduces to five as  λ 1  and  λ 2  are set to 0.
The geometry configuration of the whole floor to which the floor bays optimized belong is a row of bays sharing a secondary beam between them, as depicted in Figure 8.

6.1. Results

For each case, a multi-objective optimization has been performed. Figure 9 and Figure 10 provide the Pareto fronts of the achieved global optima for different floor spans, different vibration limitations and, with and without CLD treatment. The first one has been obtained using  S F R m a x  and the second one making use of  S F R r e s . Shaded areas in green, blue, and red represent the difference in the Pareto space between the optimal designs obtained without and with CLD. The bigger these areas, the higher the effectiveness of the CLD treatment.
For each studied case, the lightest designs (those with an optimal mass M and located at the left edge of each Pareto front) have been represented in terms of their M and  E C  depending on the span. Figure 11 and Figure 12 correspond to floors designed using  S F R m a x  and  S F R r e s , respectively. Polynomial trend curves have been included in these charts for a better interpretation of the results. The shaded areas in these figures have the same function as in the previous ones.
Finally, within the previous figures, the black circles and lines represent the static designs in which the VSLS has not been checked. When imposing a more restrictive vibration limitation, these floor designs must be changed in terms of mass, stiffness, or damping to meet the new VSLS. To provide a meaningful interpretation of these, the  R r e s  and the  R i m p  of these statically designed floors have been represented together with their  f n  and  W e f f / g  in Figure 13a,b, respectively. Shaded areas in Figure 13a represent the vibration limitation used in the parametric study. Shaded areas in Figure 13b represent the frequency restriction of 3 or 4 Hz that some design codes impose on long-span floors to avoid excessive low-frequency vibrations.

6.2. Discussion

The first finding to note in Figure 9, Figure 10, Figure 11 and Figure 12 is that the statically designed floors are the ones with less M and  E C . Also, the Pareto front in these cases converges to one point, as the concrete slab is kept as light as possible.
It is also clear from the results obtained that the more restrictive the vibration limitation, the higher the oversizing of the floor to meet the VSLS.
In floors designed without CLD, when imposing more restrictive  R l i m  values (especially 4 and 2), the Pareto front begins to open. This widening effect is more evident in two cases:
  • In long-span floors with span > 12 m (as can be seen in Figure 9d,e and Figure 10d,e). This can be explained by observing in Figure 13a,b that long-span floors have a dynamic response dominated by  R r e s . They also have low values of  f n . On the one hand, the excessive resonant vibration may be tackled by stiffening the steel profiles of the floor (which implies a high  E C  cost, but has a low repercussion in M) to increase its  f n  and thus, reduce the amplitude of the exciting human harmonic  α h . On the other hand, a solution with less  E C  cost but more impact on the final floor’s M, is to thicken the concrete slab to reduce the floor resonant response. These two possibilities are the extreme solutions of the Pareto fronts.
  • In short-span floors designed using  S F R m a x  (compare, for example, Figure 9a,b with respect to Figure 10a,b). Again, looking at Figure 13a,b, it is noticeable that short-span floors have higher values of  R r e s  and  R i m p , higher values of  f n , and lower values of  W e f f . In these floors,  R r e s  can be optimally reduced with a minimum stiffening of the floor that effectively rises  f n  the right amount to meet the VSLS (this explains why in Figure 10a,b, the Pareto front converges to 1 point). Nevertheless, reducing the floors’  R i m p  requires, either a major stiffening of the floor, (which reduces the value of the effective impulse loading  I e f f )  or a substantial increase of its  W e f f  (to reduce the impulsive response). Hence, controlling the impulsive vibration  R i m p  produces a significant oversizing of short-span floors compared to when it is not controlled (see results of Figure 11 and Figure 12 for span values lower than 12 m).
Regarding the effectiveness of the CLD integration in the final design of the floor, the following conclusions may be extracted:
  • The CLD integration enables increasing the damping ratio of the floor  ξ n . Hence, it is mainly effective when implemented in floors in which the resonant response  R r e s  is the dominant one, i.e., long-span floors with values of  R l i m  around 4 or 2. This can be appreciated in Figure 9d,e and Figure 10d,e, and also a bit in Figure 9c and Figure 10c. Table 5 provides 15 m and 18 m floors, the percentage of oversizing with respect to the statically designed cases (those not meeting the VSLS), for solutions located at the middle of the Pareto front. When CLD is not used, this oversizing in terms of M and  E C  is around 100% for floors with  R l i m  = 2, and around 50% for floors with  R l i m  = 4. When the CLD is used the oversizing decreases to around 20% for the M, and an average value of 34% for the  E C , in floors with  R l i m  = 2. When  R l i m  = 4 the oversizing decreases to 2% for the M, and 5% for the  E C . This decrease is even more evident when checking the VSLS using  S F R r e s .
  • The CLD does not result as effective in short-span floors (those with a span < 10 m). On the one hand, When assessing the VSLS according to  S F R m a x , a substantial big increase of  ξ n  does not have a significant impact on reducing the impulsive response of the floor. On the other hand, when using  S F R r e s , the resonant response of these floors seems to be better tackled by minimally stiffening the floors rather than increasing their damping. Both conclusions are clear when looking at Figure 11 and Figure 12 for spans lower than 10 m.
  • Floors with a span of around 12 m represent a transition zone between short-span and long-span floors. In this intermediate range of spans the CLD efficacy is perceptible but not as high as on long-span floors.
  • For floor designs in which M is the minimum possible, (as those depicted in Figure 11 and Figure 12) the CLD effect can be summarized as follows: for floors with a  R l i m  = 4, the CLD enables a reduction by around 50 kg/m2 and 100 kgCO2eq/m2 with respect to the cases when it is not used. In floors with  R l i m  = 2, this reduction increases with the span, with an average of 200 kg/m2 and 250 kgCO2eq/m2.

7. Conclusions

This paper studies the integration of CLD treatments into the design workflow of composite floors prone to vibrate. This treatment consists of a thin VE layer included between the steel member and the concrete slab of composite floor beams for a proportion of their length near the supports. This technology enables increasing the floor-damping ratio when vibrating in vertical bending modes, which allows for the reduction in the amount of additional mass or stiffness typically increased to overcome the VSLS.
A constrained discrete multi-objective optimization problem has been proposed to design a floor bay with different CLD treatments applied on their 1ry and 2ry beams. Seven design variables of the floor have been considered The design constraints of the problem are the different LSs of the floor, and two OFs have been used: the embodied carbon and the mass of the floor per unit of area.
Finally, a parametric study for different square floor bays, with spans varying from 4.5 to 19.5 m, has been carried out to compare the optimal structural solutions obtained with and without making use of the CLD treatment. Four different limits of vibration have been imposed from less to more restrictive. The results obtained indicate that for long-span floors (>12 m) the reduction in terms of mass and EC is substantial. The CLD enables the reduction in the structural oversizing in terms of mass from values of 100% or 50% to around 20% or 2% for floors meeting a VSLS limited by R factors of 2 or 4, respectively. Moreover, when the impulsive vibration of the floor is not checked, this enhancement is even higher.
Future work will be focused on studying the efficacy of this CLD treatment when implemented on lightweight concrete and timber floors.

Author Contributions

Conceptualization, C.M.C.R., I.M.D. and C.G.-C.; methodology, C.M.C.R., I.M.D., J.H.G.-P. and C.G.-C.; software, C.M.C.R., J.H.G.-P. and C.G.-C.; validation, C.M.C.R. and C.G.-C.; formal analysis, C.M.C.R.; investigation, C.M.C.R.; resources, C.M.C.R. and C.G.-C.; data curation, C.M.C.R.; writing—original draft preparation, C.M.C.R.; writing—review and editing, C.M.C.R., I.M.D., J.H.G.-P. and C.G.-C.; visualization, C.M.C.R.; supervision, I.M.D. and J.H.G.-P.; project administration, I.M.D. and J.H.G.-P.; funding acquisition, I.M.D. and J.H.G.-P. All authors have read and agreed to the published version of the manuscript.

Funding

The authors acknowledge Grant PID2021-127627OB-I00 (Transport Infrastructures subjected to dynamic loading: assessment techniques for the sustainability, intelligent maintenance and comfort) funded by Ministerio de Ciencia e Innovación, Agencia Estatal de Investigación and 10.13039/501100011033 FEDER, European Union. Christian Gallegos-Calderón expresses his gratitude to the Secretariat of Higher Education, Science, Technology and Innovation of Ecuador (SENESCYT) for the PhD scholarship CZ02-000167-2018.

Acknowledgments

The authors would like to thank the help provided by Pablo Vidal Fernández on the previous developments related to this work.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AISCAmerican Institute of Steel Construction
CLDConstrained layer damping
DSLSDeflection serviceability limit state
DLDead load
ECEmbodied carbon
HFFHigh-frequency floors
HIPHeathcote Industrial Polymers
LLLive load
LFFLow-frequency floors
LSLimit state
NSGA-IINon-Dominated Sorting Genetic Algorithm
OFObjective function
SCISteel Construction Institute
SFSafety factor
SWSelf-weight
ULSUltimate limit state
VEVisco-elastic
VSLSVibration serviceability limit state
1ryPrimary
2rySecondary
A s t , i Area of the steel profile of beam `i
A s l a b , i Area of the slab section belonging to the section of the composite beam `i
a g Acceleration of gravity
a p , i m p Impulsive peak acceleration
a r m s , r e s Root mean square resonant acceleration
a r m s , i m p Root mean square impulsive acceleration
BTotal width of the floor parallel to the primary beams
B e f f , i Effective width of the area involved in the fundamental mode of vibration of beams `i
b e f f , i Effective breadth of the slab of the composite beam `i
b f , i Width of the top steel flange of the beam `i
b v , i Width of the VE core of the CLD treatment applied to the beam `i
C 1 Calibration coefficient
C 2 Calibration coefficient
D S L S Functions to compute the DSLS of primary and secondary floor beams
d 2 Separation between secondary beams
E c Young modulus of concrete
ELaw of Young modulus along the beam `i
E s t Young modulus of steel
E C Embodied carbon per unit of area of the floor
F h Human-induced dynamic punctual force
f n Fundamental natural frequency of the floor
f s t e p Pacing frequency of the human walking load
f u Ultimate tensile strength
f y Steel yielding stress
f 1 , f 2 Objective function
G v Storage modulus of the VE material
gShear parameter
g o p t Optimum shear parameter
g 1 g 6 Design constraint functions
HModified Heaviside function
h c Height of the concrete slab over the ribs in cm
h v , i Height of the VE core of the CLD treatment applied to the beam `i
I c , i Concrete-Homogenized Moment of the composite section of beam `i
I e f f Effective impulse due to human footfall
I e f f , i Effective concrete-homogenized moment of inertia of beam `i
I s l a b Concrete-homogenized moment of inertia of 1 m of slab
I s t , i Moment of inertia of the steel profile of beam `i
iSubscript to indicate the type of beam: 1 for primary and 2 for secondary
K c Calibration factor
K c r e Resonant build-up factor
K 2 Calibration factor
LTotal length of the floor parallel to the secondary beams
L i Length of the beam `i
L v , i Half of the CLD-treated length of a floor beam `i
L S s l a b Functions to compute the LSs of concrete slab
MMass per unit of area of the floor
M E d Design sagging bending moment
M E d , s Service-life design sagging bending moment of a beam
M p , s t Plastic bending moment of a steel profile
M R d Resisting sagging bending moment
M R d , c o m p Resisting sagging bending moment of a composite steel–concrete section
M R d , s t Resisting sagging bending moment of a steel section
M V S L S , i Bending moments law under the loads considered for the VSLS in beam `i
N c Axial force in the concrete when a partial degree of shear connection is considered.
N c f Axial force in the concrete when the full degree of shear connection is considered.
N i Number of beams of type `i’ involved in the fundamental vibration mode of the floor
nTotal number of shear studs
P i Profile number associated with the beam `i
P r d Shear resisting force of a shear stud
p c Crossover probability
p m Probability of mutation
QAverage human weight
q D L , i Load per unit of length in beam `i’ due to deal load
q L L , i Load per unit of length in beam `i’ due to live load
q S W , i Load per unit of length in beam `i’ due to self-Weight
q V S L S , i Load per unit of length in beam `i’ considered for the VSLS
R d Number to designate the rib-deck
R i m p Impulsive response factor
R r e s Resonant response factor
S F M + , c Safety factor for sagging bending moment of the slab in construction
S F M , c Safety factor for hogging bending moment of the slab in construction
S F M + , s Safety factor for sagging bending moment of the slab in service life
S F M , s Safety factor for hogging bending moment of the slab in service life
S F M + , c , i Safety factor for sagging bending moment in the construction of beam `i
S F M + , s , i Safety factor for sagging bending moment at mid-span in the service life of beam `i
S F M + B , s , i Safety factor for sagging bending moment at section B in the service life of beam ‘i
S F R m a x , s Safety factor of VSLS considering impulsive and resonant floor response.
S F R r e s , s Safety factor of VSLS considering only resonant floor response.
S F V , c Safety factor for shear of the slab in construction
S F V , s Safety factor for shear of the slab in service life
S F V , c , i Safety factor for shear in construction of beam `i
S F δ L L , s Safety factor for deflection of the slab in service life
S F δ L L , s , i Safety factor for deflection at mid-span of beam `i
S F δ S W , c Safety factor for deflection of the slab in construction
sLongitudinal spatial coordinate along a beam element
tTime variable
U b e n d , i Modal strain energy of bending of a beam `i
U L S Functions to compute the ULS of primary and secondary floor beams
V p , s t Plastic shear force of the web of a steel profile
V S L S Function to compute the VSLS of the floor
W e f f Effective weight associated with the fundamental mode of vibration of a floor bay
W e f f , i Effective weight associated with the fundamental mode of vibration of beam `i
x ̲ Vector of design variables
YGeometric parameter of a CLD-treated beam
α h Dynamic loading factor
Δ ξ Additional damping ratio
δ i Maximum static deflection of beam ‘i’ under the loads considered for the VSLS
η c Parameter of crossover
e t a m Parameter of mutation
η v Loss factor or the VE material
λ i Parts per unit of CLD-treated length in a beam with subscript ‘i
ξ n Modal damping ratio of the fundamental mode of vibration
ξ i n t Damping ratio of the floor due to its intrinsic energy dissipation capacity
ξ C L D Additional modal damping ratio of the floor provided by the CLD treatment
ξ C L D , i Additional modal damping ratio provided by the CLD treatment to a beam ‘i
ρ Density

References

  1. Hawkins, W.; Orr, J.; Ibell, T.; Shepherd, P. A design methodology to reduce the embodied carbon of concrete buildings using thin-shell floors. Eng. Struct. 2020, 207, 110195. [Google Scholar] [CrossRef]
  2. Setareh, M. Office floor vibrations: Evaluation and assessment. Proc. Inst. Civ. Eng. Struct. Build. 2014, 167, 187–199. [Google Scholar] [CrossRef]
  3. Hawkins, W.; Peters, A.; Mander, T. A weight off your mind: Floor loadings and the climate emergency. Struct. Eng. 2021, 99, 18–20. [Google Scholar] [CrossRef]
  4. Joint Research Center—European Convention for Constructional Steelworks. Design of Floor Structures for Human Induced Vibration; Joint Research Center: Brussels, Belgium, 2009. [Google Scholar]
  5. Smith, A.L.; Hicks, S.J.; Devine, P. Design of Floors for Vibration, 2nd ed.; SCI: Singapore, 2009; pp. 1–113. [Google Scholar]
  6. Willford, M.R.; Young, P. A Design Guide for Footfall Induced Vibration of Structures, 1st ed.; The Concrete Society: Cardiff, UK, 2006; p. 82. [Google Scholar]
  7. Murray, T.M.; Allen, D.E.; Ungar, E.E.; Davis, D.B. AISC Steel Design Guide 11. Vibrations of Steel-Framed Structural Systems Due to Human Activity, 2nd ed.; AISC: Chicago, IL, USA, 2016; p. 130. [Google Scholar]
  8. Royvaran, M.; Avci, O.; Davis, B. Analysis of floor vibration evaluation methods using a large database of floors framed with W-Shaped members subjected to walking excitation. J. Constr. Steel Res. 2020, 164, 105764. [Google Scholar] [CrossRef]
  9. Gonçalves, M.S.; Pavic, A. Environmental impact of structural modifications in office floors to satisfy vibration service ability. Eur. Assoc. Struct. Dyn. 2020, 1, 1924–1931. [Google Scholar] [CrossRef]
  10. Murray, T.M.; Ungar, E.E.; Davis, D.B. Facts for Steel Buildings, Number 5: Vibration; AISC: Chicago, IL, USA, 2018; pp. 1–50. [Google Scholar]
  11. Jayasinghe, A.; Orr, J.; Hawkins, W.; Ibell, T.; Boshoff, W.P. Comparing different strategies of minimising embodied carbon in concrete floors. J. Clean. Prod. 2022, 345, 131177. [Google Scholar] [CrossRef]
  12. Regúlez, B.; Faria, D.M.; Todisco, L.; Fernández Ruiz, M.; Corres, H. Sustainability in construction: The urgent need for a new ethics. Struct. Concr. 2023, 24, 1893–1913. [Google Scholar] [CrossRef]
  13. Gallegos-Calderón, C.; Renedo, C.M.; Pulido, M.D.G.; Díaz, I.M. A frequency-domain procedure to design TMDs for lively pedestrian structures considering Human–Structure Interaction. Structures 2022, 43, 1187–1199. [Google Scholar] [CrossRef]
  14. Soria, J.M.; Díaz, I.M.; García-Palacios, J.H. Further steps towards the tuning of inertial controllers for broadband-frequency-varying structures. Struct. Control Health Monit. 2020, 27, e2461. [Google Scholar] [CrossRef]
  15. Ramírez-Senent, J.; Gallegos-Calderón, C.; García-Palacios, J.H.; Díaz, I.M. Active control of human-induced vibrations on lightweight structures via electrodynamic actuator dynamics inversion. J. Vib. Control 2022. [Google Scholar] [CrossRef]
  16. Jones, D.I. Handbook of Viscoelastic Vibration Damping; John Wiley & Sons: Hoboken, NJ, USA, 2001; p. 391. [Google Scholar]
  17. Baz, A.M. Active and Passive Vibration Damping; John Wiley & Sons: Hoboken, NJ, USA, 2019. [Google Scholar] [CrossRef]
  18. Nelson, F.C. The use of visco-elastic material to dampen vibrations in buildings and large structures. Am. Inst. Steel Constr. Eng. J. 1968, 5, 72–78. [Google Scholar]
  19. Farah, A.; Ibrahim, I.M.; Green, R. Canadian Journal of Revue canadienne Civil bngineering de genie civil. Can. J. Civ. Eng. 1977, 4, 405–4011. [Google Scholar] [CrossRef]
  20. Ebrahimpour, A. Chapter 18, Externally Bonded FRP Composite and Viscoelastic Materials for Mitigating Vibrations of Floor Systems; CRC Press: Boca Raton, FL, USA, 2013; pp. 313–329. [Google Scholar] [CrossRef]
  21. Ahmadi, H.; Goodchild, I.; Fuller, K.; Canisius, G. Modelling Dynamic Behaviour of Constrained Layer Damped Floors Using Finite Element Analysis; IRCO: Singapore, 2002; pp. 583–595. [Google Scholar]
  22. Willford, M.; Young, P.; Algaard, W.H. A constrained layer damping system for composite floors. Struct. Eng. 2006, 84, 31–38. [Google Scholar]
  23. Zahn, M.C. The Economies of LRFD in Composite Floor Beams. Eng. J. 1987, 24, 87–92. [Google Scholar]
  24. Murray, T.M. Acceptability Criterion for Occupant-Induced Floor Vibrations. Eng. J. 1981, 18, 62–70. [Google Scholar]
  25. Kim, H.; Adeli, H. Discrete cost optimization of composite floors using a floating-point genetic algorithm. Eng. Optim. 2001, 33, 485–501. [Google Scholar] [CrossRef]
  26. Adeli, H.; Kim, H. Cost optimization of composite floors using neural dynamics model. Commun. Numer. Methods Eng. 2001, 17, 771–787. [Google Scholar] [CrossRef]
  27. Klanšek, U.; Kravanja, S. Cost estimation, optimization and competitiveness of different composite floor systems-Part 2: Optimization based competitiveness between the composite i beams, channel-section and hollow-section trusses. J. Constr. Steel Res. 2006, 62, 449–462. [Google Scholar] [CrossRef]
  28. Klanšek, U.; Kravanja, S. Cost estimation, optimization and competitiveness of different composite floor systems-Part 1: Self-manufacturing cost estimation of composite and steel structures. J. Constr. Steel Res. 2006, 62, 434–448. [Google Scholar] [CrossRef]
  29. Kravanja, S.; Žula, T.; Klanšek, U. Multi-parametric MINLP optimization study of a composite I beam floor system. Eng. Struct. 2017, 130, 316–335. [Google Scholar] [CrossRef]
  30. Kaveh, A.; Ahangaran, M. Discrete cost optimization of composite floor system using social harmony search model. Appl. Soft Comput. J. 2012, 12, 372–381. [Google Scholar] [CrossRef]
  31. Kaveh, A.; Abadi, A.S.M. Cost optimization of a composite floor system using an improved harmony search algorithm. J. Constr. Steel Res. 2010, 66, 664–669. [Google Scholar] [CrossRef]
  32. Poitras, G.; Lefrançois, G.; Cormier, G. Optimization of steel floor systems using particle swarm optimization. J. Constr. Steel Res. 2011, 67, 1225–1231. [Google Scholar] [CrossRef]
  33. Kaveh, A.; Ghafari, M.H. Optimum design of steel floor system: Effect of floor division number, deck thickness and castellated beams. Struct. Eng. Mech. 2016, 59, 933–950. [Google Scholar] [CrossRef]
  34. Yossef, N.M.; Taher, S. Cost Optimization of Composite Floor Systems with Castellated Steel Beams. Pract. Period. Struct. Des. Constr. 2019, 24, 04018035. [Google Scholar] [CrossRef]
  35. Roynon, J. Emboided Carbon: Structural Sensitivity Study. Struct. I Struct. 2020, 1–7. [Google Scholar]
  36. Drewniok, M.P.; Campbell, J.; Orr, J. The Lightest Beam Method—A methodology to find ultimate steel savings and reduce embodied carbon in steel framed buildings. Structures 2020, 27, 687–701. [Google Scholar] [CrossRef]
  37. Whitworth, A.H.; Tsavdaridis, K.D. Embodied Energy Optimization of Steel-Concrete Composite Beams using a Genetic Algorithm. Procedia Manuf. 2020, 44, 417–424. [Google Scholar] [CrossRef]
  38. Kravanja, S.; Klanšek, U.; Žula, T. Mass, direct cost and energy life-cycle cost optimization of steel-concrete composite floor structures. Appl. Sci. 2021, 11, 316. [Google Scholar] [CrossRef]
  39. Gauch, H.L.; Hawkins, W.; Ibell, T.; Allwood, J.M.; Dunant, C.F. Carbon vs. cost option mapping: A tool for improving early-stage design decisions. Autom. Constr. 2022, 136, 104178. [Google Scholar] [CrossRef]
  40. Heathcote Industrial Polymers. Desing Properties of HIP2 Polymer; Heathcote Industrial Polymers: Staffordshire, UK, 2015; pp. 1–3. [Google Scholar]
  41. ASTM:E756-05; Standard Test Method for Measuring Vibration-Damping Properties of Materials. American Society for Testing and Materials. ASTM International: West Conshohocken, PA, USA, 2010.
  42. EN 1993-1-1:2005; Eurocode 3: Design of Steel Structures—Part 1-1: General Rules and Rules for Buildings. European Committee of Standarization: Brussels, Belgium, 2005.
  43. EN 1994-1-1:2004; Eurocode 4: Design of Composite steel and Concrete Structures—Part 1-1: General Rules and Rules for Buildings. European Committee of Standarization: Brussels, Belgium, 2004.
  44. AcelorMittal. GlobalFloor. Cofraplus 6. Safe Load Tables. pp. 1–12. Available online: http://www.arcelormittal-construction.com.pl/admin/user_files/page_attachments/154/arval_cofraplus-60_en_view.pdf (accessed on 4 June 2023).
  45. Muhammad, Z.; Reynolds, P.; Avci, O.; Hussein, M. Review of Pedestrian Load Models for Vibration Serviceability Assessment of Floor Structures. Vibration 2018, 2, 1–24. [Google Scholar] [CrossRef] [Green Version]
  46. Rao, D.K. Frequency and Loss Factors of Sandwich Beams Under Various Boundary Conditions. J. Mech. Eng. Sci. 1978, 20, 271–282. [Google Scholar] [CrossRef]
  47. Gibbons, O.P.; Orr, J.J.; Archer-Jones, C.; Arnold, W.; Green, D. How to Calculate Embodied Carbon, 2nd ed.; The Institution of Structural Engineers: London, UK, 2022; p. 58. [Google Scholar]
  48. Flannery, B.P.; Mares, J.W. Greenhouse Gas Index for Products in 39 Industrial Sectors: Synthetic Rubber; Resources for the Future: Washington, DC, USA, 2022; pp. 1–16. [Google Scholar]
  49. Deb, K.; Pratap, A.; Agarwal, S.; Meyarivan, T. A Fast and Elitist Multiobjective Genetic Algorithm: NSGA-II. IEEE Trans. Evol. Comput. 2002, 6, 182–197. [Google Scholar] [CrossRef] [Green Version]
  50. Deb, K.; Agrawal, S. A Niched-Penalty Approach for Constraint Handling in Genetic Algorithms; Springer: Berlin/Heidelberg, Germany, 1999; pp. 235–243. [Google Scholar]
Figure 1. Evolution in the application of CLD treatments to composite floors. (a) CLD Application of Nelson [18] over the lower flange. (b) CLD Application of Farah et al. [19] under the lower flange. (c) CLD Application of Ebrahimpour [20] with a CFRP constraining layer. (d) CLD Application of Ahmadi et al. [21] with a constraining concrete slab. (e) Resotec CLD Application developed by Willford et al. [22].
Figure 1. Evolution in the application of CLD treatments to composite floors. (a) CLD Application of Nelson [18] over the lower flange. (b) CLD Application of Farah et al. [19] under the lower flange. (c) CLD Application of Ebrahimpour [20] with a CFRP constraining layer. (d) CLD Application of Ahmadi et al. [21] with a constraining concrete slab. (e) Resotec CLD Application developed by Willford et al. [22].
Actuators 12 00288 g001
Figure 2. ‘Resotec’ CLD treatment installed along 50% of the beam’s length (50% CLD) Willford et al. [22].
Figure 2. ‘Resotec’ CLD treatment installed along 50% of the beam’s length (50% CLD) Willford et al. [22].
Actuators 12 00288 g002
Figure 3. Scheme of a floor layout composed of bays.
Figure 3. Scheme of a floor layout composed of bays.
Actuators 12 00288 g003
Figure 4. Elements of a floor bay with an integrated CLD treatment.
Figure 4. Elements of a floor bay with an integrated CLD treatment.
Actuators 12 00288 g004
Figure 5. Envelopes of resisting bending moment,  M R d , compared to the laws of maximum design bending moment,  M E d , for different  % C L D  treatment.
Figure 5. Envelopes of resisting bending moment,  M R d , compared to the laws of maximum design bending moment,  M E d , for different  % C L D  treatment.
Actuators 12 00288 g005
Figure 6. Equivalent moments of inertia homogenized to the concrete used to compute the deflection of beams with different  % C L D  treatment.
Figure 6. Equivalent moments of inertia homogenized to the concrete used to compute the deflection of beams with different  % C L D  treatment.
Actuators 12 00288 g006
Figure 7. Set of design curves used to compute  ξ C L D , i  in a simply supported isolated beam of the floor.
Figure 7. Set of design curves used to compute  ξ C L D , i  in a simply supported isolated beam of the floor.
Actuators 12 00288 g007
Figure 8. Floor configuration composed of edge bays used for the parametric study.
Figure 8. Floor configuration composed of edge bays used for the parametric study.
Actuators 12 00288 g008
Figure 9. Pareto fronts of the floors in terms of M and  E C  designed with  S F R m a x  and for different spans: (a) 6 m. (b) 9 m. (c) 12 m. (d) 15 m (e) 18 m.
Figure 9. Pareto fronts of the floors in terms of M and  E C  designed with  S F R m a x  and for different spans: (a) 6 m. (b) 9 m. (c) 12 m. (d) 15 m (e) 18 m.
Actuators 12 00288 g009
Figure 10. Pareto fronts of the floors in terms of M and  E C  designed with  S F R r e s  and for different spans: (a) 6 m. (b) 9 m. (c) 12 m. (d) 15m (e) 18 m.
Figure 10. Pareto fronts of the floors in terms of M and  E C  designed with  S F R r e s  and for different spans: (a) 6 m. (b) 9 m. (c) 12 m. (d) 15m (e) 18 m.
Actuators 12 00288 g010
Figure 11. Floors designed with  S F R m a x , optimal designs in terms of mass per square meter. (a) Mass per square meter M. (b) Embodied carbon per square meter  E C .
Figure 11. Floors designed with  S F R m a x , optimal designs in terms of mass per square meter. (a) Mass per square meter M. (b) Embodied carbon per square meter  E C .
Actuators 12 00288 g011
Figure 12. Floors designed with  S F R r e s , optimal designs in terms of mass per square meter. (a) Mass per square meter M. (b) Embodied carbon per square meter  E C .
Figure 12. Floors designed with  S F R r e s , optimal designs in terms of mass per square meter. (a) Mass per square meter M. (b) Embodied carbon per square meter  E C .
Actuators 12 00288 g012
Figure 13. Statically designed floors without CLD and with  R l i m = . (a) Dynamic response of the floors in terms of  R r e s  and  R i m p . (b) Dynamic parameters of the floors  f n  and  W e f f .
Figure 13. Statically designed floors without CLD and with  R l i m = . (a) Dynamic response of the floors in terms of  R r e s  and  R i m p . (b) Dynamic parameters of the floors  f n  and  W e f f .
Actuators 12 00288 g013
Table 1. Material properties of the floor.
Table 1. Material properties of the floor.
MaterialElement ρ  [kg/m3] E st  [gPa] f y  [mPa]
Steel S275Profiles7850210275
Steel B500SReinforcement7850210500
Steel S350GDRib-deck7850210350
ρ  [kg/m3 ] E st  [gPa] f u  [mPa]
Steel S235J2+Shear studs7850190450
ρ  [kg/m3] E c  [gPa] f ck  [mPa]
Concrete C30Slab25003030
ρ  [kg/m3] v  [mPa] η v  [-]
HIP2CLD treatment17000.71
Table 2. Static LSs checked in the 1ry and 2ry beams of the floor with integrated CLD treatment.
Table 2. Static LSs checked in the 1ry and 2ry beams of the floor with integrated CLD treatment.
StageLimit StateTypeLocationLoadsSectionSafety Factor
ConstructionULSBendingActuators 12 00288 i001 S W + L L c Steel S F M + , c , i
ShearActuators 12 00288 i002 S W + L L c Steel S F V , c , i
Service-LifeULSBendingActuators 12 00288 i003 S W + D L + L L s Steel S F M + B , s , i
BendingActuators 12 00288 i004 S W + D L + L L s Composite S F M + , s , i
ShearActuators 12 00288 i005 S W + D L + L L s Steel S F V , s , i
SLSDeflectionActuators 12 00288 i006 L L s Composite S F δ L L , s , i
Table 3. Static LSs checked in the rib-deck slab of the floor.
Table 3. Static LSs checked in the rib-deck slab of the floor.
StageLimit StateTypeLocationLoadsSectionSafety Factor
ConstructionULSBendingActuators 12 00288 i007 S W + L L c Steel S F M + , c S F M , c
ShearActuators 12 00288 i008 S W + L L c Steel S F V , c
SLSDeflection pondingActuators 12 00288 i009 S W + L L c Steel S F δ L L , c
Service-LifeULSBendingActuators 12 00288 i010 S W + D L + L L s Composite S F M + , s S F M , s
ShearActuators 12 00288 i011 S W + D L + L L s Composite S F V , s
SLSDeflectionActuators 12 00288 i012 L L s Composite S F δ L L , s
Table 4. Carbon Factors  E C F A 1 A 3 , n  used for the different materials of the floor.
Table 4. Carbon Factors  E C F A 1 A 3 , n  used for the different materials of the floor.
MaterialFloor Element ECF A 1 A 3 , n  [kgCO2eq/kg]
ConcreteSlab0.10 1
SteelReinforcement0.76 1
SteelBeam profiles1.74 1
SteelStuds1.74 1
Galvanized SteelRib-deck sheet2.87 1
Galvanized SteelCLD Sheets2.87 1
VE materialCLD VE core6.00 2
1 Recommended default values for projects in the UK given by Table 2.3 of [47].2 Assuming the use of a VE material similar to the natural or butyl rubber. Taken from [48].
Table 5. Oversizing in terms of M and  E C  of floor designs with 15 m and 18 m span, that meet different VSLSs, with respect to the statically designed floors ( R l i m = ).
Table 5. Oversizing in terms of M and  E C  of floor designs with 15 m and 18 m span, that meet different VSLSs, with respect to the statically designed floors ( R l i m = ).
VSLS Type L i [ m ] CLD TreatmentM Oversizing [%]EC Oversizing [%]
S F R m a x  − R l i m  = 215NO CLD84139
CLD4570
18NO CLD134134
CLD1828
S F R r e s  −  R l i m  = 215NO CLD84139
CLD935
18NO CLD134155
CLD516
S F R m a x  −  R l i m  = 415NO CLD5156
CLD27
18NO CLD4046
CLD25
S F R r e s  −  R l i m  = 415NO CLD6058
CLD14
18NO CLD4043
CLD24
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Renedo, C.M.C.; Díaz, I.M.; García-Palacios, J.H.; Gallegos-Calderón, C. Structural Optimization of Lightweight Composite Floors with Integrated Constrained Layer Damping for Vibration Control. Actuators 2023, 12, 288. https://doi.org/10.3390/act12070288

AMA Style

Renedo CMC, Díaz IM, García-Palacios JH, Gallegos-Calderón C. Structural Optimization of Lightweight Composite Floors with Integrated Constrained Layer Damping for Vibration Control. Actuators. 2023; 12(7):288. https://doi.org/10.3390/act12070288

Chicago/Turabian Style

Renedo, Carlos M. C., Iván M. Díaz, Jaime H. García-Palacios, and Christian Gallegos-Calderón. 2023. "Structural Optimization of Lightweight Composite Floors with Integrated Constrained Layer Damping for Vibration Control" Actuators 12, no. 7: 288. https://doi.org/10.3390/act12070288

APA Style

Renedo, C. M. C., Díaz, I. M., García-Palacios, J. H., & Gallegos-Calderón, C. (2023). Structural Optimization of Lightweight Composite Floors with Integrated Constrained Layer Damping for Vibration Control. Actuators, 12(7), 288. https://doi.org/10.3390/act12070288

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop