Real-Time Optimization of Precast Concrete Component Transportation and Storage

*e construction of prefabricated concrete structures involves many different types of precast concrete (PC) components. Improper arrangements of the transportation, storage, and hoisting of the PC components could lead to unnecessary relocation of PC components, causing a delay in construction progress. To address this problem, this paper presented a dynamic optimization method for PC component transportation and storage based on real-time scheduling and tracking. *e real-time schedule can be extracted from a 4D building information model (BIM), and the position of PC components during the whole transportation process is tracked by a proposed tracking system which integrates the global navigation satellite system (GNSS) and the radio frequency identification (RFID) technology. A transportation optimization model was built to obtain a reasonable transportation plan based on real-time construction progress. A cyclic operation network (CYCLONE) simulation model for component storing and hoisting was proposed to calculate the on-site transportation time. A storage optimization model was proposed by considering the real-time transportation information to optimize the storage mode and storage position in the yard. *e proposed models were solved in the context of a case study, which indicates that our method can reduce the on-site transportation time by 37% and effectively control the relocation times.


Introduction
Prefabricated concrete structure is recently gaining momentum in the architecture, engineering, and construction (AEC) industry. With this technique, a building product is constructed and delivered by assembling concrete components that have been precast in factories. Compared to traditional on-site concrete casting practice, the construction technique requires less on-site operation and labor resources and can yield higher construction quality within a shorter construction duration by causing less environmental pollution [1,2]. However, the construction of prefabricated concrete structures is a relatively complex process that involves a series of operations such as component transportation, storage, hoisting, and installation. During the process, the PC component transportation and storage are critically important operations [3,4], especially for the projects located in the city center, where space available for PC component storage is limited. e transportation plan of PC components determines whether there would be sufficient components and space on-site, respectively, for construction and for storage [3]. e storage layout affects the time of hoisting, relocation times, and then the construction efficiency. Here, the word "relocation" refers to the operation of removing the above components in a stack of PC components, in order to access the desired component to be installed according to the construction sequence of PC components. is operation is typical in projects where PC components are randomly piled up in a storage site and can cause unnecessary waste of time spent on finding and relocating components. Not being able to scientifically arrange the transportation plan and storage layout will result in problems such as storing difficulties and relocation, leading to undesirable delays in the construction schedule [4,5].
Despite the importance of transportation and storage of PC components, current industry practices are still confronted with many problems that hinder the PC component transportation and storage from being fully optimized. First, the transportation plan without considering the uncertainty of real-time progress of construction often brings difficulties in construction. When the as-built progress is ahead of the as-planned schedule or the delivery is delayed, it would lead to a shortage of available PC components on-site [3,6]. When the as-built progress is behind the as-planned schedule, excessive PC components would arrive at the construction site before the storage yard has sufficient space, which causes difficulties in PC component storage on-site [2]. Second, in current practice, PC components arriving at the construction site are usually stacked randomly in the storage yard without considering the real-time construction schedule and the arrival of subsequent components [4]. As a result, the chances of relocation of PC components increase. Research efforts have been made in component transportation [7][8][9][10] and storage layout optimization [11][12][13][14][15][16]. However, few of the existing researches have attempted to address the problem of transportation planning under the real-time construction schedule and optimize PC component storage layout by considering the impact of random storage and the sequence of arriving components. e research gap calls for a critical need to address the optimization problem of PC component transportation and storage in a context of the dynamically changing construction schedule and component position. Not meeting the need will continue to allow the transportation and storage plan inconsistent with the actual construction progress to misguide the construction practice, leading to extensive delays and economic loss. In light of the above limitations, this paper proposed and solved several real-time optimization models for the transportation plan and storage layout of PC components based on real-time scheduling and material position tracking. Recent advancements in computer graphics and sensing provide useful tools for incorporating dynamic construction information into our study: 4D building information modeling (BIM) allows for the collection of realtime schedule information, while radio frequency identification (RFID) and global navigation satellite system (GNSS) allow for the tracking of PC components during the whole transportation process [17][18][19][20][21][22]. e contribution of this research is threefold. First, a dynamic optimization model for the transportation plan was proposed based on real-time construction progress. e model dynamically uses the actual PC component demands retrieved from an as-built 4D BIM as one of its constraints and thus can yield a transportation plan that is consistent with the actual construction progress. Second, a systematic framework for PC component tracking was developed, which integrates GNSS and RFID techniques to enable the positioning of components in both the indoor and outdoor environments. e framework allows for the whole-process tracking of PC components and can provide real-time component arrival information for the following storage optimization. ird, a dynamic storage optimization model was proposed to guide the placement of each PC component by considering the real-time transportation tracking and construction schedule. e model can help reduce the relocation times and on-site transportation time and thus improve the construction efficiency.

Literature Review
2.1. Review of Component Transportation. Component transportation plays an important role in PC component construction. Not being able to deliver the PC components on time can bring difficulties to storage or construction progress [2,7]. To avoid the difficulties, the transportation plan of PC components needs to be arranged reasonably. Some efforts have been made on transportation plan arrangement. Sometimes, conflict exists between precast concrete component factories and the contractor because the contractor prefers earlier delivery which forces the factories to compress their production process [8]. To resolve this conflict, Zhai et al. [8] proposed three models to optimize the transportation plan and achieve a balance between the factories and the contractor. e information inconsistence for different stakeholders is another problem for transportation plan arrangement. Xu et al. [7] proposed a cloudbased fleet management platform and share transportation information among the stakeholders through it, which helps to arrange the transportation plan more reasonably. Besides, Chang et al. [9] optimized the transportation plan aiming to maximize the transport vehicle utilization rate on the basis of analyzing the geometrical characteristics and installation sequence. In order to reduce the transportation cost as well as meet the requirements, Tang et al. [10] proposed a bilevel programming model to optimize the transportation plan.

Real-Time Tracking Techniques for PC Components.
Many scholars have applied RFID technology for PC component or construction material tracking and information collection of the production supply chain and construction site [17][18][19][20]23]. e BIM is a repository of architectural engineering that integrates the information (geometry, state, location, etc.) of a building in a 3D visualization manner. By combining the BIM with RFID technology, the superiority of visualizing and tracking PC components in a BIM was demonstrated [24]. A BIM platform-integrated RFID enabled the visibility and traceability of PC components in manufacturing, logistics, and on-site construction [21]. However, the transportation information in the platform is manually updated by the driver, and real-time transportation tracking is not realized. As a widely used positioning technology, the GNSS has been applied to PC component tracking. GPS technology was used to track the transportation process of PC components, based on which the time when the components will arrive at the construction site was predicted [25]. e RFID, GPS, and GMS technologies were also used comprehensively to collect the real-time location and other information of construction materials on-site [22,26].
A 4D BIM combines the construction progress with a 3D BIM, which can dynamically simulate the construction process and display the real-time construction progress. A model which integrated the BIM and capital path method (CPM) was proposed to guide the construction site layout at different stages of construction [27]. e 4D BIM was also applied on the ETO (engineering-to-order) prefabricated project to simulate the construction process and analyze the working time of each process [28]. e method was used to arrange the transportation plan of PC components, which helped improve the construction efficiency and reduce inventory, and also played a key role in the logistics management of materials. To provide an accurate design of the transportation order, quantity, and type of the PC component, a logistics planning and control model for PC components based on the construction simulation results obtained from the 4D BIM was proposed [29].

Review of Component Storage Optimization.
e study of storage optimization for PC components was first carried out in the field of prefabricated factories [30,31]. Based on the consideration of component size, storage method, and safety distance, the concept of storage zoning was proposed, and a component storage and transportation optimization model with the goal of minimizing storage cost and transportation cost was developed [30]. e storage of PC components is an important task for the layout design of the construction site. At present, extensive studies have been focused on the optimization of materials and equipment layout. By dividing the construction site into several segments with equal size to store materials or facilities in, the construction site layout was optimized using simulation and the genetic algorism (GA) [11]. To yield the most appropriate size and location coordinates to store the equipment and materials on-site, the actual travel path was used to calculate the transportation distance on-site rather than the linear path in a BIM-based dynamic optimization model [12]. For a congested construction site where interior building spaces are also used for material storage, construction logistics optimization models were still established [13,14]. Additionally, different construction stages of construction [15] and the uncertainty of transportation distance and frequency [16] were also considered in the optimization of the layout of the construction site.

Limitations of Existing Researches.
Despite the advancements made by the above researches, existing studies on PC component transportation and storage still have the following limitations: (1) e existing researches on transportation plan arrangement overlook the influence of uncertainty of the construction schedule. e difference between the real-time progress and the designed schedule should be considered.
(2) Transportation tracking information is not taken into account in the existing researches, leading to an overlook of the impact of upcoming components on storage optimization.
(3) Few efforts have been made to optimize the storage of PC components in the construction site. e impact of random storage and repeated relocation on construction efficiency has not been considered. In addition, different types of PC components are simplified as one kind of material [12], which is not accurate enough to decide the location to store and hoist each type of component.

Methodologies
In order to address the limitations of existing researches, this paper proposed two real-time optimization models for the transportation plan and storage layout, respectively, and a CYCLONE simulation model of storing and hoisting. As shown in Figure 1, the models are based on the information retrieved from a 4D BIM and a real-time tracking system of PC components. e 4D BIM is used to obtain the real-time demand of PC components. e real-time tracking system integrated with RFID and GNSS is used to obtain the location of PC components, through which the arriving time, type, and quantity of components can be predicted. e transportation optimization model is solved by a linear programming solver, while the storage layout optimization model is solved by CYCLONE simulation and genetic algorithms.

Transportation Optimization Based on Actual Demands.
A model is proposed to optimize the transportation plan of PC components by considering the real-time construction progress and demand. A 4D BIM integrates the construction schedule with the BIM, from which the real-time construction progress can be obtained and the actual demand of components under the current progress can be further analyzed, to assist the formulation of the transportation plan.

Objectives.
e model assumes that a prefabricated project requires I different types of PC components supplied by J different factories. e on-site staff releases transportation orders of PC components to the factories according to the current construction schedule and a reasonable estimation. e required PC components will then be delivered to the construction site within T days after the order is released. e variable for the proposed model can be expressed as X(i, j, t), which stands for the quantity of type i components arriving at the site on the tth day from the jth factory (t = 1, 2, . . ., T, j = 1, 2, . . ., J, i = 1, 2, . . ., I, and X(i,j,t) is an integer). e model has the following two optimization objectives: (1) Minimizing the transportation cost Z 1 : the model assumes that the same type of vehicle is used to transport the PC components, and the transportation cost of unit distance is fixed. e transportation cost is calculated by the following equation:

Advances in Civil Engineering
where b(j, t) is the number of transport vehicles sent from the jth factory on the tth day and d(j) is the distance from the jth factory to the construction site. (2) Minimizing the difference between the actual amount and the rational amount of the PC component inventory on-site Z 2 : although the model is based on real-time progress, there still exist uncertainties in the subsequent construction. Increasing the component inventory on-site can meet the needs when the construction speed exceeds the expected speed. However, once the construction speed is lower than the expected speed, there may be no space on-site to store the excessive components. Considering the uncertainty of the demand and schedule, the on-site inventory buffer coefficient u is introduced to determine a rational amount of the PC component inventory, which is u times the demand. e objective is to minimize the difference between the inventory amount and this rational amount, as follows: where Q(i, t) is the inventory amount of the ith component on the tth day and D(i, t + 1) is the demand of the ith component on the t + 1th day.

Constrains
(1) Production Capacity Constraints. e quantity of PC components which are available in each factory each day is limited, so the transportation plan is subject to the production capacity of the factories, as shown by where C(i, j) is the quantity of the ith type of component that the j component factory is able to provide. (2) Demand Constraints. In order to meet the construction demand, the number of components stored on-site should be no less than the demand for the next day, as shown by the following equations: where Q 0 (i) is the initial stocking amount of the ith component.

(3) Storage Yard Area Constraints. e transportation plan
should be formulated within the scope of on-site storage yard capacity, as illustrated by the following equation: where A(i) is the maximum quantity of the ith component allowed to be stored in the storage yard and Q(i, t)/A(i) represents the percentage of the storage yard area that the ith component is occupying.

Model
Solving. e model proposed above is a multiobjective integer linear programming model with T × I × J integer variables and is solved by a linear programming  solver. e linear weighting method is used to construct the evaluation function. Normalization has been performed to ensure the consistency of dimensions of the two objectives, as is shown by the following equation: where a and b are the weights of the two objectives, respectively, and a + b � 1; min Z m and max Z m (m � 1, 2), respectively, represent the minimum and maximum values of the single objective model under the same constraints. Gurobi solver is used to solve the model. As a large-scale mathematical programming solver, the Gurobi solver can provide quick and precise solutions for multivariable and multiconstraint optimization models, such as integer programming, linear programming, and quadratic programming. A MATLAB interface provided by the Gurobi solver was used to derive the solutions of the model.

Whole-Process Tracking of PC Components.
A systematic framework is developed to enable the whole-process tracking of PC components, which is mainly composed of three parts, i.e., component factory management, transportation tracking, and on-site storage management.

Factory Management. RFID registration techniques
are used in the stage of factory management. Before PC components are transported out of the factory, an RFID tag is embedded in each component. e embedded tag is scanned by a handheld RFID reader to upload the information, such as the PC component ID, type, manufacturer, and current state, to the system. On being loaded onto transport vehicles, the PC components need to be scanned again to bind the components with the transport vehicles. In the meantime, the state of the PC components in the system database is changed on the handheld reader. e process is shown in Figure 2(a).

Transportation Tracking.
e transport vehicles are equipped GNSS-positioning devices to provide real-time geographical information of the vehicles. Since each vehicle has been associated with the PC components it carries at the factory management stage, the transportation route and the geographic position information of the PC components can also be obtained in real time. With the position information and an estimation of transportation speed, the arriving time of the PC components can be predicted to guide the planning of the storage layout in advance. Figure 2(b) shows an example of the transportation route and the position of a vehicle.

On-Site Storage Management.
On arrival, the RFID tags embedded in the PC components are scanned to update the location information of the components to the system. According to the storage location determined by the storage optimization model, the PC components are placed into the storage yard. During construction, the component is hoisted from the location which is determined by the simulation model. A BIM of the storage yard is created and dynamically updated to reflect the actual PC component layout, as is shown in Figure 2(c).

Premises and Assumptions.
e proposed model is based on the following premises and assumptions: (1) e real-time schedule obtained from the 4D BIM and transportation tracking information is readily available so that the arrival time, type, and quantity of in-transit components in a certain time period (T days) can be predicted. (2) e time period (T days) of the storage optimization model is set equal to that of the transportation optimization model proposed above to ensure that the transportation plan and tracking information are available.
(3) T is divided into K shorter time periods with the equal interval A, on which the simulation of component storing and hoisting is based. At the beginning of each shorter time period, the components that have arrived during the previous shorter period are uniformly stacked. e proposed model only aims to optimize the storage position of PC components at the beginning of the first shorter period of T days to ensure that the newest schedule and transportation tracking information are used. Besides, component storing and hoisting in a certain time period (T days) are considered in the storage optimization model to analyze the impact of component storage in the first period on subsequent construction and components' storage. (4) In order to ensure that the quantity of PC component storage on-site at least satisfies the construction demand in one day, the model assumes that the components demanded in one day have been already stored on-site.  [11,30], this paper divides a storage yard into E equal-sized storage zones. Only components of the same storage mode can be stacked in the same storage zone before the storage zone is empty again, and the number of storage zones is determined by the site manager. Each storage zone is further divided into several bays, where a piece or a pile of components is stored. Because the components of the same storage mode tend to have similar sizes, the quantity of bays in each zone is determined Advances in Civil Engineering 5 according to the storage mode in the storage zone. Figure 3 shows a schematic diagram of a storage yard (E = 6), where e and m represent the storage zone number and the bay number, respectively.

Optimization Variables.
e variable X of the storage optimization model is set as the state (type and quantity of components) of each bay in each storage zone of the storage yard. In this way, the size of variables is determined by the storage yard only and easy to control, which can help improve the solving efficiency. e state of the yard is updated every time the component storage is completed in a shorter time period. erefore, the variable X consisted of K storage states, and each state X k represents the storage state of period k after storage, k = 1, 2, . . ., K.
e storage state at the end of each period consists of the states of each bay. A column vector is used to represent the storage state of a bay. e column vector can be divided into two parts, as shown in Figure 4. e first part is the first element of the vector, which indicates the component storage mode of the bay-0 for no component stored, 1 for the vertical storage, 2 for the tiled storage, 3 for the sametype stacking, and 4 for the mixed-type stacking. e remaining elements make up the second part, which indicate the type and quantity of components. e column vector length U is determined by the maximum allowable amount of the mixed-type storage components h max , U � h max + 1, which will be explained below.
When the storage mode of a zone is vertical storage, tiled storage, or same-type storage, the stacking order in a bay of the zone is not considered. Note that each bay can only contain a vertical storage component or a tiled storage component. Hence, for the above three storage modes, only the first three elements of the column vector are needed, which represent the storage condition of the bay, the component type, and the quantity, respectively. e remaining elements are set as 0. An example of such vectors is shown in Figure 4(a). When the storage mode of the zone is mixed-type stacking, the second to the last elements, respectively, indicate the type of the components stored in the bay from top to bottom. As shown in Figure 4(b), 4 components are piled up in the bay, and the types of them from the top layer to the bottom layer are i 2 , i 3 , i 4 , and i 4 , respectively. By combining the column vectors of all the M bays in E zones at all the K shorter periods, the storage state variable, which is a matrix of size U × (KEM), can be obtained, and M is determined by the largest number of bays among the E storage zones.

Objectives
(1) Minimizing the weighted transportation distance P in the construction site: the transportation distance is represented by the linear distance from the center point of the storage zone to the tower crane.
Considering that components with larger size are more difficult to transport than smaller ones, the transportation difficulty coefficient w i is introduced to indicate the transportation difficulty of the type i component. e objective function is formulated as where y i,k,e is the quantity of type i components hoisted from the storage zone e, during period k, and d e is the linear distance from the storage zone e to the tower crane.   Advances in Civil Engineering (2) Minimizing relocation times R: the objective function is formulated as the following equation: where r i,k is the number of relocation times when the type i component is hoisted in period k.

Constraints
(1) Space Constraints. e number of bays must not exceed the maximum number of bays allowed in a storage zone, as illustrated by where z is the storage mode of a bay; A z is the maximum allowed number of bays under the current storage mode; k � 1, 2, . . ., K; e � 1, 2, . . ., E; and m � 1, 2, . . ., M.

Advances in Civil Engineering
(2) Stacking Height Constraints. For components that can be piled up, the number of layers must not exceed the maximum number allowed by safety considerations: where H 3 max and H 4 max are the maximum allowable number of piled-up layers of the same-type stacking components and the mixed-type stacking components, respectively, and m � 1, 2, . . ., M.
(3) Total Amount Constraint. e constraint of the total component amount can be formulated as where Q k,i is the quantity of the components of type i at the end of period k; D i,k represents the demand of the component of type i during the k period; and q i,k is the components of type i stored in period k; k = 1, 2, ..., K. Second, when the storage mode of a zone is vertical storage, tiled storage, or same-type storage, only the first 3 rows of the column vector are occupied with the storage mode, the component type, and the quantity and the elements of the remaining rows must be 0.

GA-Simulation-Based Solution for the Storage Optimization Model.
e proposed storage optimization model is a multiobjective integer programming model with U × (KEM) integer variables and is solved by a simulation-based genetic algorithm (GA). e flowchart of the solving process is shown in Figure 5. e genetic algorithm is used to generate variables that satisfy the constraints. e generated variables and parameters are then input to the simulation model of the component storing and hoisting process for transportation time calculation. e GA optimization operation persists until the fitness value of the results in certain generation meets the terminating condition. After getting the optimized storage state by the GA-simulation method, the specific storage location of each component is to be calculated by the difference between the optimized state and the state just before the storage happens, which is obtained from the simulation model. e storage location of each PC component is output as a result of the model at the end.

CYCLONE Simulation Model for PC Component
Storing and Hoisting. In order to calculate the on-site transportation time, the process of the PC component storing and hoisting is simulated. CYCLONE has been applied in construction simulation of dams, underground powerhouses, tunnels [32,33], etc. In this paper, the CYCLONE simulation technology is adopted, as shown in Figure 6. e simulation runs on a period basis. After the components are hoisted for N times in period k, the components arriving during the period are uniformly stacked at the end of the period. us, the functional node is set in the simulation model. e number of cycles N is the number of components being hoisted during the period, which is determined by the average construction speed. Since the storage of PC components does not occupy the tower crane and has no effect on the construction progress, the time for component storage is ignored in calculating the transportation time on-site. e calculation process of on-site transportation time is shown in Figure 7. e data input at the beginning of the simulation is as follows: the initial state of the storage yard, the hoisting order I k (k = 1, 2, . . ., K), the total quantity of components to be hoisted N sum , the quantity of components in each period N k (k = 1, 2, . . ., K), the time cost for each relocation, the time to connect a component and the crane, the vertical and horizontal speed of the tower crane, the distance between storage zones and the tower crane, the number of optimization periods K, and the state of the storage yard X. e same type of components can be stored in different places in the storage yard. Hence when a component is needed, the place from where the component is hoisted should be decided. Two principles are generally used to determine the place, i.e., minimizing relocation times and minimizing distance. ese two principles might conflict with each other sometimes. To compare the influence of the two principles on the construction speed, two simulation rules are proposed: one considers the distance first and the other considers relocation times first. e on-site transportation time t(n) of a component can be calculated as follows: 8 Advances in Civil Engineering where t d is the time to connect the unit component and crane;   erefore, the on-site transportation time obtained by the simulation model can reflect the fitness value of the solution. Besides, the equality constraint (total amount constraint) is difficult to satisfy in generating the initial population and the genetic operation, so a penalty function is used to reflect whether the equality constraint is satisfied. e fitness function is as follows:

Genetic Algorithm for Solving Storage Optimization
where T T is the on-site transportation time obtained by simulation; M is the penalty coefficient; and n p is the deviation between the two sides of the equality constraint.
(2) Coding Method. e variable of the model is a matrix, which cannot be simply and directly represented by a conventional one-bit coding method (binary coding, real coding, etc.). erefore, the matrix coding is used to represent the variable, as shown in Figure 8.
(3) Initial Population. e initial population of the model is difficult to generate with a traditional random method for its large dimensions and complicated constraints. erefore, the following two principles are proposed to guide the generation of initial population: First, only after one storage zone is fully filled, the next storage zone can be used to store components of the same storage mode. Second, the generated initial population and genetic operation should meet the area constraint, the height constraint, and the logical constraint, and the quantity constraint is included as a penalty function.
(4) Genetic Operation. e algorithm adopts the roulette selection strategy, and the probability of being selected as parents increases with the fitness value. An elitism strategy is used to ensure that good genes can be inherited to offsprings. e offspring population is generated through the following genetic operations: values than an average value will have a lower crossover rate to allow their good genes to be kept in offspring population. On the contrary, individuals with a lower fitness rate will have higher crossover rates. e operator can ensure the diversity of the population and prevent it from falling into a local optimum. e formula for calculating the crossover rate is shown as where P c is the crossover rate; f max is the maximum fitness value of the population; f avg is the average fitness value of the population; f is the larger fitness value of the two individuals to be crossed; and k 1 and k 2 are parameters and k 1 < k 2 . e single-point crossover is used as the crossover mode, and the crossover node n is generated randomly. e node n represents period n in the storage model. As shown in Figure 9(a), the storage state from period 1 to period n of the parent i is combined with the storage states from period n + 1 to period N (the last period) of the parent j to constitute a new offspring individual.
(ii) Mutation operator: the mutation node n is randomly generated. e variables before the node are reserved and the variables after the node are regenerated, which means that the storage state after the nth period is reestablished, as shown in Figure 9(b). (iii) In order to prevent the solution falling into a local optimum, some of the offspring individuals will be regenerated to introduce new genes into the population when the optimal solution does not change within the certain number of iterations.

Project Background.
e D project is located in Z Province, China. It consists of 13 townhouses, 9-storey bungalows, and 4 high-rise buildings with a total construction area of 100595.6 m 2 . e structure type is the frame shear wall structure combined with PC prefabrication. PC composite beams, PC composite slabs, PC stairs, and other PC components are used in the project. RFID and BeiDou (Chinese Global Navigation Satellite System) positioning systems are used to track PC components in real time.
e project uses the 4D BIM to reflect the construction progress and to carry out dynamic project management. Figure 10 shows interfaces of the 4D BIM. e component quantity demanded in next two days is calculated through the BIM based on the current progress, and a list of PC component requirements is shown in the table at the bottom right of Figure 10. In order to implement real-time tracking and visualization of PC components on the yard, the storage yard is also separately modeled. e storage yard BIM is modeled based on the information obtained from the wholeprocess tracking system of PC components and displays the state of the storage yard in real time. Connecting the storage yard BIM with RFID technology supports to realize the whole-process management of PC components on-site.

Model Establishment and Solution.
e construction of building #21 in the project was taken as an example. e transportation plan was released each two days in the project. Consequently, the optimization period T of both the transportation optimization model and the storage optimization model was set as 2 days. According to the current progress, the construction of a standard floor should be completed in the period, and a total of 68 PC composite slabs of 42 types are going to be used in the construction. Assuming that the construction speed was uniformly distributed, 34 components would be constructed per day. Table 1 shows the demand and supply of the PC components. e "De." and "Seq." columns show the demand and construction sequence of each type of component in the optimization period (2 days), respectively. e "Fac." column shows the code number of the factory that is able to supply the components, and the corresponding quantity that can be supplied per day by the factories is shown in the "Supply per day" column. e quantity of PC components initially stored in the storage yard was exactly equal to the quantity demanded on the first day. e storage yard had a size of 10 × 21 m 2 , which was divided into six 5 × 7 m 2 storage zones. e horizontal linear distance from each storage zone to the tower crane was 15.3 m, 16.9 m, 19.3 m, 20.5 m, 23.3 m, and 24.3 m, respectively. e distance between each two bays in a storage zone was at least 1 m. e PC composite slabs were stored in the mixed-type storage mode, and two storage bays were set in each storage zone considering the maximum size of components. e components cannot pile up with more than 6 layers in a bay for the safety consideration.

Solution of Transportation Optimization Model.
ree factories, i.e., factories I, II, and III, supplied PC components for the project, which were 200 km, 160 km, and 300 km away from the construction site, respectively. e  Table 1. PC components were transported using 10-meter long flatbed trailers, and each trailer can be loaded with 12 PC components.
According to the site situation, the value of the on-site inventory buffer coefficient k was set as 1.5 [25]. e weights of the two objectives (i.e., minimum transportation cost and minimum difference between actual and rational amount of the inventory) were set as 0.6 and 0.4, respectively. e transportation optimization model was solved by the Gurobi solver, and the optimal fitness value is 0.0949. e optimized transportation plan for two days is shown in Table 2, and the component quantity of each type is shown in the brackets. Component requirements in the next two days

12
Advances in Civil Engineering

Solution of Storage Optimization Model.
e optimization period was two days, and the working hours of each day was 8 hours.
erefore, the optimization time period T was 16 hours. e PC components were stored each 2 hours (A = 2), and the number of periods K = T/A = 8. e number of storage zones E equaled 6, and the number of bays in each storage zone M equaled 2. Hence, the size of the matrix variable was 7 × 96. In this case study, the optimization model was used to obtain the storage location of components arriving during the first one of the 8 periods. erefore, the model was solved according to the real-time transportation tracking information at the beginning of the first period. At this point, the expected delivery time and the sequence of the PC components in the next 8 time periods can be predicted, as shown in Table 3.
(1) Results of the GA-Simulated Solution.
e maximum weight of the PC components in this project was 1.7 t, and the nominal speed of the tower crane was constant when the hoisted components are less than 3 t. erefore, the transportation difficulty coefficient w i of each component is 1.
e nominal lifting speed of the tower crane is 40 m/min, the horizontal speed is 20 m/min, and the lifting height of the components is set to 12 m. e time to connect a component with the crane is 4 minutes, and the time cost per relocation is 5 minutes. e number of functional node cycles N in the simulation model is the number of components hoisted during each time period. According to the construction speed, N of the eight periods is 9, 8, 9, 8, 9, 8, 9, and 8, respectively. e GA-simulation method was used to solve the model. e population size was set as 1000, and the Advances in Civil Engineering 13 number of iterations was 2000. Two parameters of the adaptive crossover rate k 1 and k 2 were 0.75 and 0.8, respectively, and the mutation rate was 0.1. In the simulation model, the two principles are separately used to calculate onsite transportation time. Figure 11 shows the convergence curves of the two principles.
In terms of minimum relocation times, the fitness value of the optimal solution is 2.094. e number of relocation times in the optimization period (2 days) is 8. e horizontal transportation distance is 1247 meters, and the on-site transportation time is 478 minutes. In terms of minimum distance, the fitness value of the optimal solution is 2.017. e number of relocation times is 12. e horizontal transportation distance is 1229 meters, and the transportation time is 496 minutes. e results demonstrate that the optimal solution obtained by the principle of minimum relocation time has a larger fitness value and a shorter on-site transportation time.
In the case study, the difference of horizontal transportation distance of each storage zone was relatively small, and the tower crane had a relatively high transportation speed. erefore, the impact of difference between the two principles on on-site transportation time was no more than 1 minute on each component. Consequently, the principle of minimum relocation time was more applicable in this case. Figure 12 shows the storage state of period 1 optimized by the principle of minimum relocation time.
e storage position of each component can be obtained by subtracting the initial storage state from the state after storage in period 1, as is in the table "Components to be stored" of Figure 12. e number of layers was counted from top to bottom. e components in a darker colour were already stacked before optimization, while the components in a lighter colour were to be stacked after optimization.
(2) Result Simulation. e state after storage in each time period can be obtained according to the optimal solution. e storage positions of PC components in period k determine the position where a component should be hoisted from during period k. erefore, the simulation model was called again to calculate the position and transportation time of each component in period k. e results of this simulation operation can be used to direct the component hoisting in construction. As an example, Table 4 shows the position and on-site transportation time of each component during period 3.

Discussion.
In order to evaluate the effectiveness of the dynamic storage optimization model, a control group model which does not consider the real-time transportation tracking information was developed, and the corresponding on-site transportation time was calculated. In the control group model, PC components which were on the way to the construction site or planned to be transported were not considered. e optimization time period T of the control group model was 2 hours. In order to keep the optimization time T of the control group model consistent with that of the original model, the storage of PC components in 8 periods was optimized separately, by assuming that the actual arrival time was the same as the results predicted through the transportation tracking. e simulation principle in the control group model was minimum relocation times. e on-site transportation time, the horizontal transportation distance, and the number of relocation times in each period were obtained. Figure 13 shows the comparison of the results of the control group model and the original model.
As shown in Figure 13, the differences of the number of relocation times and on-site transportation time between the two models are not significant until the later periods of the simulation. In period 5, the on-site transportation time even exceeds the duration of a period (i.e., 120 minutes), which means that too many times of relocation of PC components lead to delay in construction. is is because the components hoisted in the earlier periods had been initially stored in the storage yard before optimization. Nevertheless, the components arriving during optimization periods would be used in subsequent periods. Without transportation tracking to obtain the real-time information of the components' arrival, the control group model ignores the impact of oncoming components, resulting in a significant increase of relocation times and the on-site transportation time in later periods.
In the case study, the storage yard is small, while the quantity of the component types is relatively huge. erefore, even though the transportation tracking information has been considered, there still exist relocations in the later periods of the original model. However, compared with the situation where the transportation tracking information is not considered, the number of relocation times is reduced considerably.

Conclusions
Given the uncertainty of the construction schedule and transportation condition, this study developed the real-time optimization models for the transportation plan and storage layout of PC components based on real-time scheduling and material position tracking. A 4D BIM was adopted to retrieve the real-time schedule. A systematic framework was developed to enable the whole-process tracking of PC components. In order to solve the storage layout optimization model, a CYCLONE simulation model was developed  to simulate the process of component storing and hoisting, and genetic algorithms were used to optimize the on-site transportation time and the relocation times. A case study was implemented in a project located in Z Province in China. e results demonstrate the effectiveness of the proposed models from the following three aspects: (1) Compared with the situation without transportation tracking, the results obtained by the storage optimization model can effectively reduce the number of relocation times and on-site transportation time in the PC component construction. rough the real-time tracking system, the whole-process tracking and systematic management of PC components were realized, which can provide the information from leaving factory, transportation, storage on-site, to hoisting and construction.
Further research is needed to address the following limitation: In the proposed models, the transportation time on-site is calculated by considering the linear distance from the storage position to the tower crane position instead of the specific storage position to the installation position, which is not accurate enough. In addition, the installation time of a component is also ignored. In the future study, the BIM can be used to calculate the transportation time of PC components from the yard to the specific installation position, and the coordinated operation of multiple tower cranes needs to be studied to establish a more accurate simulation model to analyze the on-site transportation time.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.