Open access
Technical Papers
Aug 4, 2021

Lifting Sequence Optimization of Luffing Tower Cranes Considering Motion Paths with Dynamic Programming

Publication: Journal of Construction Engineering and Management
Volume 147, Issue 10

Abstract

The lifting sequence of luffing tower cranes is a key factor in the normal operation of construction projects. Accurate estimation of total lifting time is necessary for lifting sequence optimization. Existing formulations for lifting time estimation have two main deficiencies. One is that the relationship between the hoist motion path of the hook and the luffing motion path of the boom is often neglected. The other is that lifting delays resulting from the horizontal motion path of the boom caused by the relative locations of the crane and the lifting start and end points are rarely taken into account. To address those limitations, this paper proposes a lifting sequence optimization model (LSOM) considering motion paths with dynamic programming. The effectiveness of the proposed model is evaluated by comparing it with three conventional lifting strategies [first in–first serve (FIFS), shortest job first (SJF), and nearest neighbor first (NNF)]. The results show that LSOM achieves a shorter total lifting time and higher crane utilization when compared with FIFS, SJF, and NNF.

Introduction

With continuous urbanization around the world, more construction projects take place in dense urban areas, which limit the space of construction sites (Ibrahim et al. 2018). Construction sites with limited space pose a severe challenge to the layout of tower cranes (Tam et al. 2001; Abdelmegid et al. 2015; Ji and Leite 2018). According to their luffing modes, tower cranes can be mainly classified as T-type and luffing tower cranes (Wu et al. 2020). Compared with the T-type tower crane with a fixed boom, the luffing tower crane is easier to lay out because of its structural characteristics, which are ideal for construction sites with limited space (Li et al. 2012).
Construction sites with limited space typically have higher requirements for space utilization (Al-Hussein et al. 2006; Irizarry and Karan 2012; Marzouk and Abubakr 2016). As one of the main space users, the luffing tower crane plays an important role in space utilization, and thus the lifting sequence is a key factor in normal project operation (Lee et al. 2015; Younes and Marzouk 2018). Total lifting time is one of the most basic elements in evaluating the effectiveness of a lifting sequence. That is, the less the total lifting time in the sequence, the higher the space utilization (Younes and Marzouk 2018; Wu and García de Soto 2020).
Having an accurate estimate of total lifting time is beneficial to lifting sequence optimization. Similar to that of the T-type tower crane, the lifting cycle time of the luffing tower crane mainly depends on its potential no-load lifting from the last demand point to the current supply point, essential loaded lifting from the current supply point to the current demand point, and potential lifting delay to avoid space conflicts with other cranes. However, the lifting details are different due to the separate luffing modes of the two crane types (Zhang et al. 1999; Lee et al. 2015; Younes and Marzouk 2018). Compared with that of the T-type tower crane, luffing tower crane lifting time estimation has received little attention. Therefore, it is necessary to establish a lifting time estimation formulation for the luffing tower crane according to its motion characteristics (Lee et al. 2015; Younes and Marzouk 2018). Specifically, it is necessary to consider the relationship between the hoist motion path of the hook and the luffing motion path of the boom as well as the lifting delay resulting from the horizontal motion path of the boom caused by the relative locations of the crane and the lifting start and end points.
An effective optimization method can enhance lifting sequence optimization. First in–first serve (FIFS), shortest job first (SJF), and nearest neighbor first (NNF) are conventional lifting strategies that decrease total lifting time in some cases, but cannot guarantee optimal total lifting time (Huang and Wong 2018). Lifting sequence optimization is determining the optimal lifting task and crane at each stage in the lifting process. Dynamic programming can realize the global optimum of multistage decision-making which matches the lifting sequence problem (Lazarev and Werner 2009; Bürgy et al. 2020). Although dynamic programming is seldom applied to such optimization problems, it is expected to be effective.
To sum up, this paper proposes a lifting sequence optimization model (LSOM) to minimize total luffing tower crane lifting time. The estimation formulation is explicitly established by taking the crane’s motion characteristics into account, and dynamic programming is used to enhance the optimization solution. The remainder of this paper is structured as follows. In the section “Literature Review,” studies on lifting sequence optimization are reviewed and summarized. In the section “Lifting Sequence Optimization Model,” the LSOM is established. The section “Example,” evaluates LSOM by comparing it with FIFS, SJF, and NNF in a realistic situation. In the section “Discussion,” the significance and limitation of LSOM are discussed. The section “Conclusion,” draws conclusions and discusses future work.

Literature Review

In construction projects, determination of the lifting sequence of tower cranes means sequence planning of lifting tasks and corresponding cranes (Wu and García de Soto 2020). In order to improve lifting efficiency, the lifting sequence is usually optimized using conventional lifting strategies such as FIFS, SJF, and NNF and, in some cases, specific lifting sequence optimization models (Zavichi et al. 2014; Huang and Wong 2018). Common objective metrics are minimum total lifting time, total lifting cost, pending time, and so on; the exact choice depends on the needs of project teams (Monghasemi et al. 2016; Al Hattab et al. 2017).
Lifting time estimation formulation is the basis of lifting sequence optimization (Rodriguez-Ramos and Francis 1983). In the 1970s, the relationship between lifting time, crane location, and lifting start and end point locations was established, providing the opportunity for quantitative evaluation of lifting time (Warszawski 1973). The luffing motion of T-type tower cranes is via trolley motion on the boom track, while that of luffing tower cranes is via boom motion; therefore lifting time estimation formulations for the two are different (Younes and Marzouk 2018). The formulation for T-type cranes mainly considers radial motion time of the trolley and slewing motion time of the boom in the horizontal plane; hook hoist motion time in the vertical plane; and loading and unloading time (Zhang et al. 1996, 1999). According to lifting locations (of T-type tower cranes and lifting start and end points), crane specifications (e.g., trolley radial speed, boom slewing speed, hook hoist speed), motion coordination (radial and slewing motion coordination in the horizontal plane and motion coordination in the horizontal and vertical planes), and lifted object characteristics (e.g., weight, loading and unloading time), lifting time can be estimated (Huang et al. 2011; Lien and Cheng 2014; Wang et al. 2015; Moussavi Nadoushani et al. 2017).
Unlike the formulation for T-type tower cranes in the horizontal plane, that for luffing tower cranes mainly takes into account the duration of boom luffing and slewing motion. Based on the Cartesian coordinate system, a formulation for luffing tower crane motion duration in the horizontal and vertical planes was proposed by Lee et al. (2015). In this formulation, luffing and slewing motion in the horizontal plane and hook hoist motion in the vertical plane were considered consecutive. To estimate the simultaneous motion time for luffing tower cranes, another formulation added the slewing and luffing motion coordination in the horizontal plane and the motion coordination in the horizontal and vertical planes (Younes and Marzouk 2018).
A type of lifting time estimation formulation was created according to the historical data of construction projects, including multiple linear regression, multilayer feed-forward neural network, general regression neural network, and group method of data handling (Leung and Tam 1999; Leung et al. 2001; Tam et al. 2002; Tam and Tong 2003). The factors affecting lifting time can be classified as controllable and uncontrollable. The controllable factors are manipulated by planning, such as crane specifications and supply point locations; the uncontrollable factors include demand point locations and lifted object weights (Leung et al. 2001). The factors affecting loaded and no-load lifting time are partly different. For loaded lifting time, factors such as lifting height, lifted object areas, and supply point locations have a greater impact (Tam et al. 2002). For no-load lifting time, factors such as lifting height and motion coordination have a more significant impact (Tam et al. 2002). Therefore, formulations for loaded and no-load lifting time estimation were created separately (Leung and Tam 1999; Leung et al. 2001; Tam et al. 2002; Tam and Tong 2003).
Compared with the multiple linear regression model and the multilayer feed-forward neural network, the general regression neural network shows good prediction of loaded and no-load lifting time. However, it cannot describe the relationship between input and output variables (Leung et al. 2001; Tam et al. 2002). The performance of the group method of data handling is slightly lower than that of the general regression neural network, but still achieves significant and acceptable prediction; meanwhile, it can reveal the relationship between input and output variables in the form of a complicated set of mathematical expressions (Tam et al. 2002).
Consideration of lifting delay to avoid space conflicts among tower cranes cannot be neglected in lifting time estimation (Hwang 2012; Wu et al. 2020). Four main criteria are proposed to identify these conflicts. The first criterion regards the location of a crane and those of its lifting start and end points as three points of a triangle. If two triangles intersect, space conflict occurs and increases as the number of intersection among six sides of two triangles increases (Zhang et al. 1999). The second criterion regards crane coverage as a circle. If two circles intersect, space conflict occurs and increases as the intersection area increases (Shapira and Simcha 2009; Irizarry and Karan 2012). The third criterion considers hook position. If two crane hooks are within the intersection area, space conflict is likely (Al Hattab et al. 2018). The fourth criterion further formalizes the previous criteria according to category, height, and motion path. Space conflict occurs in the following situations: (1) the trolley motion path of a taller T-type tower crane intersects with the boom motion path of a shorter T-type tower crane; (2) the boom slewing motion path of a luffing tower crane intersects with that of another luffing tower crane irrespective of the height ranking between them; (3) the trolley motion path of the taller T-type tower crane intersects with the boom slewing motion path of the shorter luffing tower crane; and (4) the boom slewing path of a taller luffing tower crane intersects with the boom motion path of a shorter T-type tower crane (Wu et al. 2020).
In earlier studies, the possibility of space conflicts was incorporated into tower crane planning as the objective index or constraint (Zhang et al. 1999; Irizarry and Karan 2012). Later studies began to focus on quantitative evaluation of the lifting delay caused by potential space conflicts among tower cranes (Al Hattab et al. 2018; Younes and Marzouk 2018; Wu et al. 2020).
Due to the luffing mode of luffing tower cranes, the luffing motion of the boom is associated with the hoist motion of the hook. Therefore, the hoist motion distance of the hook depends not only on the height difference between the lifting start and end points but also on the boom vertical motion distance. However, the relationship between the hook hoist motion path and the boom luffing motion path is often neglected, which leads to overestimation or underestimation of the hoist motion time of the hook. The lifting delay of luffing tower cranes is related to the relative locations of the cranes and the lifting start and end points, and the boom horizontal motion path caused by different relative locations can result in various lifting delays. However, such factors are rarely considered, making delay time estimation inaccurate. In this paper, the lifting time estimation formulation for luffing tower cranes presented here is aimed at solving these two problems.
Optimization methods include mathematical programming and metaheuristic algorithms (Wu and García de Soto 2020). In mathematical programming, the lifting sequence is formulated as an integer programming problem to minimize total lifting time (Zavichi et al. 2014; Huang and Wong 2018). In metaheuristic algorithms, harmony search and tabu search are used to search the optimal lifting sequence in terms of the least deviation of pending time and minimum total lifting time, respectively (Monghasemi et al. 2016; Al Hattab et al. 2017; Wu and García de Soto 2020).
As a branch of mathematical programming, dynamic programming was developed by Bellman (1957). Dynamic programming transforms the multistage decision-making problem into a series of single-stage decision-making subproblems according to its characteristics, and then solves the subproblems one by one to reach the optimal decision-making sequence of the whole problem (Lazarev and Werner 2009; de Matos et al. 2015). It solves each subproblem only once, and if there is a common subproblem, the stored state is used directly without repeated decision making (Lazarev and Werner 2009; de Matos et al. 2015; Grevtsov et al. 2017). Dynamic programming has been widely applied in various sectors with significant results (Konak et al. 2011; Grevtsov et al. 2017; Bürgy et al. 2020). Since it is seldom applied to lifting sequence optimization for luffing tower cranes, this paper addresses how to integrate such optimization problems with dynamic programming to extend its application.

Lifting Sequence Optimization Model

In this section, LSOM is elaborated. The model comprises two parts: (1) lifting time estimation considering motion paths and (2) integration of lifting time estimation and dynamic programming.

Loaded Lifting Time

Loaded lifting is essential in a lifting task. Fig. 1 shows loaded lifting in the horizontal and vertical planes; (Cxki,Cyki,Czki), (Sxmi,Symi,Szmi), (Exni,Eyni,Ezni), and (Sxm*i,Sym*i,Szm*i) are the locations of crane k, lifting start point m (current supply point), lifting end point n (current demand point), and transformed lifting start point m* (current transformed supply point) in lifting task i, respectively.
Fig. 1. Loaded lifting in (a) horizontal; and (b) vertical planes.
Loaded lifting time (LTk,m,ni) depends on loading time (Tlk,m,ni), loaded motion time (Tmk,m,ni), and unloading time (Tuk,m,ni), as expressed in Eq. (1)
LTk,m,ni=Tlk,m,ni+Tmk,m,ni+Tuk,m,ni
(1)
Loaded motion time (Tmk,m,ni) is determined by horizontal motion time (Thk,m,ni), vertical motion time (Tvk,m,ni), and motion coordination in the horizontal and vertical planes (β) between 0 and 1, as expressed in Eq. (2); β being equal to 0 or 1 means simultaneous or consecutive motion and typically defaults to 1 (Zhang et al. 1999)
Tmk,m,ni=max(Thk,m,ni,Tvk,m,ni)+β×min(Thk,m,ni,Tvk,m,ni)
(2)
Horizontal motion time (Thk,m,ni) is determined by boom slewing motion time (Tωk,m,ni), boom luffing motion time (Tβk,m,ni), and slewing and luffing motion coordination in the horizontal plane (α) between 0 and 1, as expressed in Eq. (3); α being equal to 0 or 1 means simultaneous or consecutive motion and typically defaults to 0.25 (Zhang et al. 1999)
Thk,m,ni=max(Tωk,m,ni,Tβk,m,ni)+α×min(Tωk,m,ni,Tβk,m,ni)
(3)
Boom slewing motion time (Tωk,m,ni) is determined by the slewing angle of the boom from lifting start point m to lifting end point n (θk,m,ni) and by boom slewing speed (Vωki), as expressed in Eq. (4). Slewing angle (θk,m,ni) is calculated using Eqs. (5)–(8)
Tωki=θk,m,niVωki
(4)
θk,m,ni=arccos(ρ(Eni)2+ρ(Smi)2(lm,ni)22×ρ(Eni)×ρ(Smi))0arccos(ρ(Eni)2+ρ(Smi)2(lm,ni)22×ρ(Eni)×ρ(Smi))π
(5)
ρ(Smi)=(SxmiCxki)2+(SymiCyki)2
(6)
ρ(Eni)=(ExniCxki)2+(EyniCyki)2
(7)
lm,ni=(ExniSxmi)2+(EyniSymi)2
(8)
Boom luffing motion time (Tβk,m,ni) is determined by the luffing angle from the transformed lifting start point m* to lifting end point n (θk,m*,ni) and by luffing speed (Vβki), as expressed in Eq. (9). Luffing angle (θk,m*,ni) is calculated using Eqs. (10)–(13)
Tβki=θk,m*,niVβki
(9)
θk,m*,ni=|arccos(ρ(Sm*i)BLki)arccos(ρ(Eni)BLki)|0arccos(ρ(Sm*i)BLki)π,0arccos(ρ(Eni)BLki)π
(10)
ρ(Sm*i)=(Sxm*iCxki)2+(Sym*iCyki)2
(11)
Sxm*i=Cxki+ρ(Smi)×ExniCxkiρ(Eni)
(12)
Sym*i=Cyki+ρ(Smi)×EyniCykiρ(Eni)
(13)
Vertical motion time (Tvk,m,ni) is determined by hoist motion distance of the hook from lifting start point m to lifting end point n (HDk,m,ni) and by hoist speed of the hook (Vvki), as expressed in Eq. (14). Hoist motion distance of the hook (HDk,m,ni) is determined by the height difference between lifting end point n (Ezni) and lifting start point m (Szmi) and by boom vertical motion distance from the transformed lifting start point m* to lifting end point n (BDk,m*,ni). Boom vertical motion distance (BDk,m*,ni) is calculated using Eq. (15)
Tvki=HDk,m,niVvki
(14)
BDk,m*,ni=|(BLki)2(ρ(Sm*i))2(BLki)2(ρ(Eni))2|
(15)
If the horizontal distance between crane k and lifting start point [ρ(Smi)] is equal to the horizontal distance between crane k and lifting end point n [ρ(Eni)], as shown in Fig. 2, hoist motion distance of the hook (HDk,m,ni) is calculated using Eq. (16)
HDk,m,ni=|EzniSzmi|ρ(Smi)=ρ(Eni)
(16)
Fig. 2. Loaded lifting in (a) horizontal; and (b) vertical planes when ρ(Smi)=ρ(Eni).
If the horizontal distance between crane k and lifting start point m [ρ(Smi)] is greater than the horizontal distance between crane k and lifting end point n [ρ(Eni)], it is classified as one of three situations. In the first situation, the height of lifting end point n (Ezni) is less than or equal to that of lifting start point m (Szmi), and the hoist motion distance of the hook (HDk,m,ni) is calculated using the first fraction of Eq. (17), as shown in Figs. 3(a and b). In the second situation, the height of lifting end point n (Ezni) is greater than that of lifting start point m (Szmi), and this height difference is less than or equal to the boom vertical motion distance from transformed lifting start point m* to lifting end point n (BDk,m*,ni); the hoist motion distance of the hook (HDk,m,ni) is calculated using the second fraction of Eq. (17), as shown in Figs. 3(a and c). In the third situation, the height of lifting end point n (Ezni) is greater than that of lifting start m (Szmi), and this difference is greater than the boom vertical motion distance from transformed lifting start point m* to lifting end point n (BDk,m*,ni); the hoist motion distance of the hook (HDk,m,ni) is calculated using the third fraction of Eq. (17), as shown in Figs. 3(a and d)
HDk,m,ni={|EzniSzmi|+BDk,m*,niifρ(Smi)>ρ(Eni)andEzniSzmi0BDk,m*,ni|EzniSzmi|ifρ(Smi)>ρ(Eni)  and0<EzniSzmiBDk,m*,ni|EzniSzmi|BDk,m*,niifρ(Smi)>ρ(Eni)  andBDk,m*,ni<EzniSzmi
(17)
Fig. 3. Loaded lifting in (a) horizontal; and (b–d) vertical planes when ρ(Smi)>ρ(Eni).
If the horizontal distance between crane k and lifting start point [ρ(Smi)] is less than the horizontal distance between crane k and lifting end point n [ρ(Eni)], it is also classified as one of three situations. In the first situation, the height of lifting start point m (Szmi) is less than or equal to that of the lifting end point n (Ezni), and the hoist motion distance of the hook (HDk,m,ni) is calculated using the first fraction of Eq. (18), as shown in Figs. 4(a and b). In the second situation, the height of lifting start point m (Szmi) is greater than that of lifting end point n (Ezni) and this height difference is less than or equal to the boom vertical motion distance from transformed lifting start point m* to lifting end point n (BDk,m*,ni); the hoist motion distance of the hook (HDk,m,ni) is calculated using the second fraction of Eq. (18), as shown in Figs. 4(a and c). In the third situation, the height of lifting start point m (Szmi) is greater than that of lifting end point n (Ezni) and this height difference is greater than the boom vertical motion distance from transformed lifting start point m* to lifting end point n (BDk,m*,ni); the hoist motion distance of the hook (HDk,m,ni) is calculated using the third fraction of Eq. (18), as shown in Figs. 4(a and d)
HDk,m,ni={|EzniSzmi|+BDk,m*,niif  ρ(Smi)<ρ(Eni)  andSzmiEzni0BDk,m*,ni|EzniSzmi|ifρ(Smi)<ρ(Eni)   and0<SzmiEzniBDk,m*,ni|EzniSzmi|BDk,m*,niif  ρ(Smi)<ρ(Eni)   andBDk,m*,ni<SzmiEzni
(18)
Fig. 4. Loaded lifting in (a) horizontal; and (b–d) vertical planes when ρ(Smi)<ρ(Eni).

No-Load Lifting Time

No-load lifting in a lifting task depends on the current lifting number of crane k (Cki) and the locations of lifting start point m (Smi; last demand point) and lifting end point n (Eni; current supply point) in lifting task i. This is classified as one of three situations.

Situation 1

In Situation 1, the current lifting number of crane k (Nki) is 0, so crane k starts moving from the current supply point. Lifting task i does not have a no-load lifting process, and no-load lifting time (NTk,m,ni) is equal to 0, as expressed in the first fraction of Eq. (19).

Situation 2

In Situation 2, the current lifting number of crane k (Nki) is not 0, and the locations of lifting start point m (Smi) and lifting end point n (Eni) are the same. Therefore, crane k starts moving from the current supply point, lifting task i does not have a no-load lifting process, and no-load lifting time (NTk,m,ni) is equal to 0, as expressed in the second fraction of Eq. (19).

Situation 3

In Situation 3, the current lifting number of crane k (Nki) is not 0, and the locations of lifting start point m (Smi) and lifting end point n (Eni) are not the same. Therefore, crane k starts moving from the last demand point, and lifting task i has a no-load lifting process. Estimation of no-load lifting time (NTk,m,ni) is the same as that of loaded lifting time (LTk,m,ni) except for loading time (Tlk,m,ni) and unloading time (Tuk,m,ni). It only needs to replace the current supply point and the current demand point with the last demand point and the current supply point, respectively, as expressed in the third fraction of Eq. (19)
NTk,m,ni={0ifNki=00ifNki0  andSmi=Enimax(Thk,m,ni,Tvk,m,ni)+β×min(Thk,m,ni,Tvk,m,ni)ifNki0  andSmiEni
(19)

Lifting Delay Time

Lifting delay in a lifting task depends on the location of crane k (Cki), lifting start point m (Smi), and lifting end point n (Eni) in current lifting task i and relative to the location of crane k (Cki), lifting start point m (Smi), and lifting end point n (Eni) in previous lifting task i. This is represented by one of the following four situations.

Situation 1

In Situation 1, none of the three stages of crane k (Cki)—boom at lifting start point m (Smi), the motion path, and lifting end point n (Eni)—intersects with the boom of crane k (Cki), and lifting task i does not have a lifting delay. There is only one case in this situation (Case 1 in Fig. 5):
Case 1: the boom of crane k (Cki) at lifting start point m (Smi), the motion path, and lifting end point n (Eni) do not intersect with the boom of k (Cki) at lifting start point m (Smi), the motion path, or lifting end point n (Eni).
Fig. 5. Sixteen cases of relative locations of two lifting tasks.

Situation 2

In Situation 2, the only stage of crane k (Cki)—the boom at lifting start point m (Smi)—intersects with the boom of crane k (Cki), and lifting task i has a lifting delay. There are six cases in this situation (Cases 2–7 in Fig. 5):
Case 2: the boom of crane k (Cki) at lifting start point m (Smi) intersects with the boom of crane k (Cki) at lifting start point m (Smi) and does not intersect with the boom of k (Cki) at the motion path or lifting end point n (Eni).
Case 3: the boom of crane k (Cki) at lifting start point m (Smi) intersects with the boom of crane k (Cki) at lifting start point m (Smi) and the motion path, and does not intersect with the boom of crane k (Cki) at lifting end point n (Eni).
Case 4: the boom of crane k (Cki) at lifting start point m (Smi) intersects with the boom of crane k (Cki) at lifting start point m (Smi), the motion path, and lifting end point n (Eni).
Case 5: the boom of crane k (Cki) at lifting start point m (Smi) intersects with the boom of crane k (Cki) at the motion path and does not intersect with the boom of crane k (Cki) at lifting start point m (Smi) or lifting end point n (Eni).
Case 6: the boom of crane k (Cki) at lifting start point m (Smi) intersects with the boom of crane k (Cki) at the motion path and lifting end point n (Eni) and does not intersect with the boom of crane k (Cki) at lifting start point m (Smi).
Case 7: the boom of crane k (Cki) at lifting start point m (Smi) intersects with the boom of crane k (Cki) at lifting end point n (Eni) and does not intersect with the boom of crane k (Cki) at lifting start point m (Smi) or the motion path.

Situation 3

In Situation 3, at least two stages of crane k (Cki), including the boom at lifting start point m (Smi), intersect with the boom of crane k (Cki), and lifting task i has a lifting delay. There are six cases in this situation (Cases 8–13 in Fig. 5):
Case 8: the boom of crane k (Cki) at lifting start point m (Smi) intersects with the boom of crane k (Cki) at lifting start point m (Smi) and does not intersect with the boom of crane k (Cki) at the motion path or lifting end point n (Eni); the boom of crane k (Cki) at the motion path and/or lifting end point n (Eni) intersect(s) with the boom of crane k (Cki) at lifting start point m (Smi); neither of them intersects with the boom of crane k (Cki) at the motion path or lifting end point n (Eni).
Case 9: the boom of crane k (Cki) at lifting start point m (Smi) intersects with the boom of crane k (Cki) at lifting start point m (Smi) and does not intersect with the boom of crane k (Cki) at the motion path or lifting end point n (Eni); the boom of crane k (Cki) at the motion path and/or lifting end point n (Eni) intersect(s) with the boom of crane k (Cki) at the motion path; neither of them intersects with the boom of crane k (Cki) at lifting end point n (Eni).
Case 10: the boom of crane k (Cki) at lifting start point m (Smi) intersects with the boom of crane k (Cki) at lifting start point m (Smi) and does not intersect with the boom of crane k (Cki) at the motion path or lifting end point n (Eni); the boom of crane k (Cki) at the motion path and/or lifting end point n (Eni) intersect(s) with the boom of crane k (Cki) at lifting end point n (Eni).
Case 11: the boom of crane k (Cki) at lifting start point m (Smi) intersects with the boom of crane k (Cki) at the motion path and does not intersect with the boom of crane k (Cki) at lifting end point n (Eni); the boom of crane k (Cki) at the motion path and/or lifting end point n (Eni) intersect(s) with the boom of crane k (Cki) at lifting start point m (Smi) and/or the motion path; neither of them intersects with the boom of crane k (Cki) at lifting end point n (Eni).
Case 12: the boom of crane k (Cki) at lifting start point m (Smi) intersects with the boom of crane k (Cki) at the motion path and does not intersect with the boom of crane k (Cki) at lifting end point n (Eni); the boom of crane k (Cki) at the motion path and/or lifting end point n (Eni) intersect(s) with the boom of crane k (Cki) at lifting end point n (Eni).
Case 13: the boom of crane k (Cki) at lifting start point m (Smi) intersects with the boom of crane k (Cki) at lifting end point n (Eni); the boom of crane k (Cki) at the motion path and/or lifting end point n (Eni) intersect(s) with the boom of crane k (Cki) at lifting start point m (Smi), the motion path, and/or lifting end point n (Eni),

Situation 4

In Situation 4, at least one stage of crane k (Cki), excluding the boom at lifting start point m (Smi), intersects with the boom of crane k (Cki), and lifting task i has a lifting delay. There are three cases in this situation (Cases 14–16 in Fig. 5):
Case 14: the boom of crane k (Cki) at the motion path and/or lifting end point n (Eni) intersects with the boom of crane k (Cki) at lifting start point m (Smi); neither of them intersects with the boom of k (Cki) at the motion path or lifting end point n (Eni).
Case 15: the boom of crane k (Cki) at the motion path and/or lifting end point n (Eni) intersect(s) with the boom of crane k (Cki) at the motion path; neither of them intersects with the boom of crane k (Cki) at lifting end point n (Eni).
Case 16: the boom of crane k (Cki) at the motion path and/or lifting end point n (Eni) intersect(s) with the boom of crane k (Cki) at lifting end point n (Eni).
According to the distinction between loaded and no-load lifting of current lifting task i and previous lifting task i, estimation of lifting delay time (DTk,m,ni) in the 16 cases is divided into four situations, as shown in Fig. 6, which assumes that the two lifting tasks start at the same time and that i is completed before i.
Fig. 6. Estimation of lifting delay time in 16 cases.
In Situation 1, lifting tasks i and i are both in no-load lifting, and the estimation of lifting delay time (DTk,m,ni) is divided into two fractions.
In Case 1, 2, 5, 6, 7, 8, or 14, crane k (Cki) starts moving once crane k (Cki) starts moving; lifting delay time (DTk,m,ni) is 0, as expressed in the first fraction of Eq. (20).
In Case 3, 4, 9, 10, 11, 12, 13, 15, or 16, crane k (Cki) starts moving once crane k (Cki) ends moving; lifting delay time (DTk,m,ni) is the difference between the end time of no-load lifting in lifting task i (ENTk,m,ni) and the start time of no-load lifting in lifting task i (SNTk,m,ni), as expressed in the second fraction of Eq. (20)
DTk,m,ni={0if  Case1,2,5,6,7,8,or14ENTk,m,niSNTk,m,niifCase3,4,9,10,11,12,13,15,or16
(20)
In Situation 2, lifting tasks i and i are in no-load and loaded-lifting, respectively, and the estimation of lifting delay time (DTk,m,ni) is divided into four fractions:
In Case 1, 5, 6, or 7, crane k (Cki) starts moving once crane k (Cki) starts loading; lifting delay time (DTk,m,ni) is 0, as expressed in the first fraction of Eq. (21).
In Case 2, 8, or 14, crane k (Cki) starts moving once crane k (Cki) starts moving; lifting delay time (DTk,m,ni) is the difference between the start time of the loaded motion stage of loaded lifting in lifting task i (STmk,m,ni) and the start time of the loading stage of loaded lifting in lifting task i (STlk,m,ni), as expressed in the second fraction of Eq. (21).
In Case 3, 9, 11, or 15, crane k (Cki) starts moving once crane k (Cki) ends moving; lifting delay time (DTk,m,ni) is the difference between the end time of the loaded motion stage of loaded lifting in lifting task i (ETmk,m,ni) and the start time of the loading stage of loaded lifting in lifting task i (STlk,m,ni), as expressed in the third fraction of Eq. (21).
In Case 4, 10, 12, 13, or 16, crane k (Cki) starts moving once crane k (Cki) ends unloading; lifting delay time (DTk,m,ni) is the difference between the end time of the unloading stage of loaded lifting in lifting task i (ETuk,m,ni) and the start time of the loading stage of loaded lifting in lifting task i (STlk,m,ni), as expressed in the fourth fraction of Eq. (21)
DTk,m,ni={0ifCase1,5,6,or7STmk,m,niSTlk,m,niifCase2,8,or14ETmk,m,niSTlk,m,niifCase3,9,11,or15ETuk,m,niSTlk,m,niifCase4,10,12,13,or16
(21)
In Situation 3, lifting tasks i and i are in loaded and no-load lifting, respectively, and the estimation of lifting delay time (DTk,m,ni) is divided into three fractions:
In Case 1, 2, 8, or 14, crane k (Cki) starts loading once crane k (Cki) starts moving; lifting delay time (DTk,m,ni) is 0, as expressed in the first fraction of Eq. (22).
In Case 3, 4, 5, 6, 7, 11, 12, or 13, crane k (Cki) starts loading once crane k (Cki) ends moving; lifting delay time (DTk,m,ni) is the difference between the end time of no-load lifting in lifting task i (ENTk,m,ni) and the start time of no-load lifting in lifting task i (SNTk,m,ni), as expressed in the second fraction of Eq. (22).
In Case 9, 10, 15, or 16, crane k (Cki) starts loading once crane k (Cki) starts moving and waits to start moving until k (Cki) ends moving; lifting delay time (DTk,m,ni) is 0 if the end time of the loading stage of loaded lifting in lifting task i (ETlk,m,ni) is not earlier than the end time of no-load lifting in lifting task i (ENTk,m,ni); otherwise, lifting delay time (DTk,m,ni) is the difference between the end time of no-load lifting in lifting task i (ENTk,m,ni) and the end time of the loading stage of loaded lifting in lifting task i (ETlk,m,ni), as expressed in the third fraction of Eq. (22)
DTk,m,ni={0ifCase1,2,8,or14ENTk,m,niSNTk,m,niifCase3,4,5,6,7,11,12,or130  or  ENTk,m,niETlk,m,niifCase9,10,15,or16
(22)
In Situation 4, lifting tasks i and i are both in loaded lifting, and the estimation of lifting delay time (DTk,m,ni) is divided into ten fractions:
In Case 1, crane k (Cki) starts loading once crane k (Cki) starts loading. Lifting delay time (DTk,m,ni) is 0, as expressed in the first fraction of Eq. (23).
In Case 2 or 8, crane k (Cki) starts loading once crane k (Cki) starts moving. Lifting delay time (DTk,m,ni) is the difference between the start time of the loaded motion stage of loaded lifting in lifting task i (STmk,m,ni) and the start time of the loading stage of loaded lifting in the lifting task i (STlk,m,ni), as expressed in the second fraction of Eq. (23).
In Case 3, 5, or 11, crane k (Cki) starts loading once crane k (Cki) ends moving. Lifting delay time (DTk,m,ni) is the difference between the end time of the loaded motion stage of loaded lifting in lifting task i (ETmk,m,ni) and the start time of the loading stage of loaded lifting in lifting task i (STlk,m,ni), as expressed in the third fraction of Eq. (23).
In Case 4, 6, 7, or 13, crane k (Cki) starts loading once crane k (Cki) ends unloading. Lifting delay time (DTk,m,ni) is the difference between the end time of the unloading stage of loaded lifting in lifting task i (ETuk,m,ni) and the start time of the loading stage of loaded lifting in lifting task i (STlk,m,ni), as expressed in the fourth fraction of Eq. (23).
In Case 9, crane k (Cki) starts loading once crane k (Cki) starts moving and waits to start moving until crane k (Cki) ends moving. Lifting delay time (DTk,m,ni) is the difference between the start time of the loaded motion stage of loaded lifting in lifting task i (STmk,m,ni) and the start time of the loading stage of loaded lifting in lifting task i (STlk,m,ni) if the end time of the loading stage of loaded lifting in lifting task i (ETlk,m,ni) is not earlier than the end time of the loaded motion stage of loaded lifting in lifting task i (ETmk,m,ni); otherwise, lifting delay time (DTk,m,ni) is their difference plus the difference between the end time of the loaded motion stage of loaded lifting in lifting task i (ETmk,m,ni) and the end time of the loading stage of loaded lifting in lifting task i (ETlk,m,ni), as expressed in the fifth fraction of Eq. (23).
In Case 10, crane k (Cki) starts loading once crane k (Cki) starts moving and waits to start moving until crane k (Cki) ends unloading. Lifting delay time (DTk,m,ni) is the difference between the start time of the loaded motion stage of loaded lifting in lifting task i (STmk,m,ni) and the start time of the loading stage of loaded lifting in lifting task i (STlk,m,ni) if the end time of the loading stage of loaded lifting in lifting task i (ETlk,m,ni) is not earlier than the end time of the unloading stage of loaded lifting in lifting task i (ETuk,m,ni); otherwise, lifting delay time (DTk,m,ni) is their difference plus the difference between the end time of the unloading stage of loaded lifting in lifting task i (ETuk,m,ni) and the end time of the loading stage of loaded lifting in lifting task i (ETlk,m,ni), as expressed in the sixth fraction of Eq. (23).
In Case 12, crane k (Cki) starts loading once crane k (Cki) ends moving and waits to start moving until crane k (Cki) ends unloading. Lifting delay time (DTk,m,ni) is the difference between the end time of the loaded motion stage of loaded lifting in lifting task i (ETmk,m,ni) and the start time of the loading stage of loaded lifting in lifting task i (STlk,m,ni) if the end time of the loading stage of loaded lifting in lifting task i (ETlk,m,ni) is not earlier than the end time of the unloading stage of loaded lifting in lifting task i (ETuk,m,ni); otherwise, lifting delay time (DTk,m,ni) is their difference plus the difference between the end time of the unloading stage of loaded lifting in lifting task i (ETuk,m,ni) and the end time of the loading stage of loaded lifting in lifting task i (ETlk,m,ni), as expressed in the seventh fraction of Eq. (23).
In Case 14, crane k (Cki) starts loading once crane k (Cki) starts loading and waits to start moving until crane k (Cki) starts moving. Lifting delay time (DTk,m,ni) is 0 if the end time of the loading stage of loaded lifting in lifting task i (ETlk,m,ni) is not earlier than the start time of the loaded motion stage of loaded lifting in lifting task i (STmk,m,ni); otherwise, lifting delay time (DTk,m,ni) is the difference between the start time of the loaded motion stage of loaded lifting in lifting task i (STmk,m,ni) and the end time of the loading stage of loaded lifting in lifting task i (ETlk,m,ni), as expressed in the eighth fraction of Eq. (23).
In Case 15, crane k (Cki) starts loading once crane k (Cki) starts loading and waits to start moving until crane k (Cki) ends moving. Lifting delay time (DTk,m,ni) is 0 if the end time of the loading stage of loaded lifting in lifting task i (ETlk,m,ni) is not earlier than the end time of the loaded motion stage of loaded lifting in lifting task i (ETmk,m,ni); otherwise, lifting delay time (DTk,m,ni) is the difference between the end time of the loaded motion stage of loaded lifting in lifting task i (ETmk,m,ni) and the end time of the loading stage of loaded lifting in lifting task i (ETlk,m,ni), as expressed in the ninth fraction of Eq. (23).
In Case 16, crane k (Cki) starts loading once crane k (Cki) starts loading and waits to start moving until crane k (Cki) ends unloading. Lifting delay time (DTk,m,ni) is 0 if the end time of the loading stage of loaded lifting in lifting task i (ETlk,m,ni) is not earlier than the end time of the unloading stage of loaded lifting in lifting task i (ETuk,m,ni); otherwise, lifting delay time (DTk,m,ni) is the difference between the end time of the unloading stage of loaded lifting in lifting task i (ETuk,m,ni) and the end time of the loading stage of loaded lifting in lifting task i (ETlk,m,ni), as expressed in the tenth fraction of Eq. (23)
DTk,m,ni={0ifCase1STmk,m,niSTlk,m,niifCase2  or  8ETmk,m,niSTlk,m,niifCase3,5,or11ETuk,m,niSTlk,m,niifCase4,6,7,or13STmk,m,niSTlk,m,ni  or  STmk,m,niSTlk,m,ni+ETmk,m,niETlk,m,niifCase9STmk,m,niSTlk,m,ni  or  STmk,m,niSTlk,m,ni+ETuk,m,niETlk,m,niifCase10ETmk,m,niSTlk,m,ni  or  ETmk,m,niSTlk,m,ni+ETuk,m,niETlk,m,niifCase120  or  STmk,m,niETlk,m,niifCase140  or  ETmk,m,niETlk,m,niifCase150  or  ETuk,m,niETlk,m,niifCase16
(23)

Dynamic Programming

The lifting process can be divided into interrelated stages according to the lifting sequence, and a lifting task and a crane should be determined at each stage to minimize total lifting time. In this way, lifting sequence optimization is formulated into dynamic programming.
The decision-making process for lifting tasks and cranes at different stages is shown in Fig. 7, where n, i, and k refer to stage n, lifting task i, and crane k, respectively. The state variables s0,s1,,sn represent lifting state information such as lifting tasks performed before or after a given stage. For example, s1 is lifting state information after Stage 1 and before Stage 2. The decision variables x1(i,k),x2(i,k),,xn(i,k) represent the performed lifting task and the performing crane at a given stage. For example, decision x1(2,3) indicates Task 2 performed by Crane 3 at Stage 1. The effect variables (i.e., stage index functions)—e1(s0,x1(i,k)),e2(s1,x2(i,k)),,en(sn1,xn(i,k))—represent current total lifting time at a given stage when a previous state and a decision are determined. For example, e1{s0,x1(2,3)} indicates current total lifting time at Stage 1 if lifting state information before Stage 1 is given and Task 2 is performed by Crane 3 at Stage 1
sn=Tn(sn1,xn(i,k))
(24)
E1,n=max{e1(s0,x1(i,k)),e2(s1,x2(i,k)),,en(sn1,xn(i,k))}
(25)
Tt1,n=minE1,n
(26)
Fig. 7. Decision-making process for lifting tasks and cranes at different stages.
The state transfer function is expressed in Eq. (24). State sn after stage n is determined according to state sn1 before stage n and the decision xn(i,k) at stage n. The relationship between the process index function and the stage index functions is expressed in Eq. (25), where process index E1,n in Stages 1n is the maximum of the stage index functions e1(s0,x1(i,k)),e2(s1,x2(i,k)),,en(sn1,xn(i,k))—that is, the maximum current total lifting time at all stages. The optimal process index function is expressed in Eq. (26), and the optimal process index Tt1,n in Stages 1n is the minimum of existing process indexes E1,n—that is, the minimum total lifting time in existing solutions.
Fig. 8 shows the dynamic programming pseudocode for lifting sequence optimization. The input is the set of lifting tasks {1,2,,i} and the set of cranes {1,2,,k}. The output is the minimum total lifting time Tt1,n and its corresponding lifting sequence. Dynamic programming adopts the bottom-up recursive method in Stages 1n, and the state after a stage depends on its previous state. The set of tasks {1,2,,i} and the set of cranes {1,2,,k} form a two-level nested loop from outside to inside. If Stage n exceeds the last stage of the lifting process, this loop is skipped and the next loop starts directly; else, the following three conditions are evaluated successively.
If the number of existing states is greater than or equal to the current stage value, remove the last state until it is one less than the current stage value.
If the criteria for avoiding infeasible solutions are satisfied, determine state sn according to Eq. (24); else, skip this loop and start the next loop. The three criteria to be satisfied are as follows: (1) the current lifting task has not yet been performed; (2) the locations of supply and demand points as well as the weight in the current task are covered by the capacity of the performing crane; and (3) the masts of other cranes are not covered by the motion range of the performing crane.
If the current stage value is equal to the number of lifting tasks, determine process index E1,n and optimal process index Tt1,n according to Eqs. (25) and (26), respectively.
Fig. 8. Dynamic programming pseudocode for lifting sequence optimization.

Example

To evaluate the effectiveness of LSOM, an example with 12 lifting tasks and two luffing tower cranes in a 75-story building project with a total height of 348.5 m was tested. The example is a typical arrangement in construction projects, involving lifting tasks with different weights and locations as well as potential space conflicts between luffing tower cranes. FIFS, SJF, and NNF are compared with LSOM in this example, and lifting tasks are performed according to number, loaded lifting time, and no-load lifting time, respectively.
The structure is based on a concrete-filled steel tube frame with a reinforced concrete core tube and outrigger trusses. The installation of Cranes C1 and C2 is by external climbing. There are seven supply points and nine demand points. Supply points S1, S2, and S3, and demand point D9 are on the ground; Demand points D5, D6, D7, and D8 are on the 23rd floor of the outer frame; common supply/demand points S4 (D1), S5 (D2), S6 (D3), and S7 (D4) are on the 29th floor of the core tube, as shown in Fig. 9. The coordinates of the cranes and supply and demand points are summarized in Table 1.
Fig. 9. Layout of the example elements.
Table 1. Coordinates of cranes and supply and demand points
NameCoordinate (X)Coordinate (Y)Coordinate (Z)
C114.0500.000
C214.0500.750
S132.65019.6000.200
S240.3500.1000.200
S340.2207.6000.200
S4 (D1)4.7758.200124.750
S5 (D2)4.7758.200124.750
S6 (D3)4.7758.200124.750
S7 (D4)4.7758.200124.750
D513.33017.65098.650
D613.33017.65098.650
D713.33017.65098.650
D813.33017.65098.650
D926.65029.2500.200
Lifting weights and supply and demand points are provided in Table 2. The lifting weight of Tasks 1–4 is 11 tons, and the supply and demand points are on the ground floor and the 29th floor of the core tube, respectively. The lifting weight of Tasks 5–8 is 15 t, and the supply and demand points are on the ground floor and the 23rd floor of the outer frame, respectively. The lifting weight of Tasks 9–12 is 12 t, and the supply and demand points are on the 29th floor of the core tube and the ground floor, respectively.
Table 2. Lifting task characteristics
TaskWeight (t)Supply pointDemand point
111S1D1
211S1D2
311S1D3
411S1D4
515S2D5
615S2D7
715S3D6
815S3D8
912S4D9
1012S5D9
1112S6D9
1212S7D9
Crane specifications are listed in Table 3; they consist of boom length, maximum working radius, minimum working radius, maximum lifting capacity, whole-range luffing time, boom slewing speed, and hook hoist speed. These are 55.0 m, 52.5 m, 4.4 m, 50.0 t, 1.5 min, 0.8 rpm, and 75  m/min, respectively.
Table 3. Crane specifications
Boom length (m)Maximum working radius (m)Minimum working radius (m)Maximum lifting capacity (t)Whole-range luffing time (min)Slewing speed of boom (rpm)Hoist speed of hook (m/min)
55.052.54.450.01.50.875
All tests are done in Visual Studio version 2015 with C# on a laptop with a Windows 7 64-bit operating system, Intel(R) Corei7-6700HQ CPU @ 2.60GHz, and 16 GB of RAM. The FIFS, SJF, and NNF computation time is less than 1 min, and that of LSOM about 1 h.
Lifting sequences for FIFS, SJF, NNF, and LSOM are provided in Table 4. In FIFS, the sequence is Tasks 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, and 12. Tasks 1, 2, 3, 4, 5, and 6 are performed by Crane C1 successively; Tasks 7, 8, 9, 10, 11, and 12 are performed by Crane C2 successively; total lifting time is 154.27 min. In SJF, the lifting task sequence is 1, 2, 4, 3, 11, 12, 9, 10, 8, 6, 5, and 7; Tasks 1, 2, 4, 3, 6, and 5 are performed by Crane C1 successively; Tasks 11, 12, 9, 10, 8, and 7 are performed by Crane C2 successively; total lifting time is 126.83 min. In NNF, the lifting sequence is 12, 7, 10, 8, 11, 9, 5, 1, 2, 3, 6, and 4; Tasks 5, 1, 2, 3, 6, and 4 are performed by Crane C1 successively; Tasks 12, 7, 10, 8, 11, and 9 are performed by Crane C2 successively; total lifting time is 166.61 min. In LSOM, the lifting task sequence is 3, 9, 1, 8, 2, 4, 5, 12, 7, 6, 10, and 11; Tasks 3, 1, 2, 4, 5, and 6 are performed by Crane C1 successively; Tasks 9, 8, 12, 7, 10, and 11 are performed by Crane C2 successively; the total lifting time is 89.22 min.
Table 4. FIFS, SJF, NNF, and LSOM lifting sequences
SequenceFIFSSJFNNFLSOM
TaskCraneSchedule (min)TaskCraneSchedule (min)TaskCraneSchedule (min)TaskCraneSchedule (min)
StartEndDurationStartEndDurationStartEndDurationStartEndDuration
1110.006.036.03110.006.036.031220.008.108.10310.006.296.29
2216.0314.108.07216.0314.108.07728.1035.0326.93920.008.148.14
33114.1022.438.334114.1022.368.2610235.0343.768.73116.2914.608.31
44122.4330.938.503122.3630.858.498243.7670.5926.83828.1434.9726.83
55130.9359.7828.8511222.3638.8316.4711270.5979.188.592114.6022.678.07
66159.7888.0628.2812238.8349.0310.209279.1889.4610.284122.6730.938.26
77259.7886.4726.699249.0359.3110.285179.18105.8226.645130.9359.7828.85
88286.47114.7528.2810259.3169.7510.4411105.82113.387.5612234.9743.608.63
992114.75123.688.938269.7596.5826.8321113.38121.468.087243.6070.5326.93
10102123.68134.1210.446169.7598.5528.8031121.46129.798.336159.7888.0628.28
11112134.12144.079.955198.55126.8328.2861129.79158.5928.8010270.5379.268.73
12122144.07154.2710.207298.55126.8328.2841158.59166.618.0211279.2689.229.96
Compared with FIFS, SJF, and NNF, LSOM decreases total lifting time by 65.05, 37.61, and 77.39 min, respectively, which corresponds to respective decrease rates of 42.17%, 29.65%, and 46.45%.
Crane utilization in FIFS, SJF, NNF, and LSOM is listed in Table 5. It signifies the ratio of each crane’s lifting time, excluding delay time, to total lifting time and indicates that the greater the crane utilization, the greater the crane’s space utilization. In FIFS, utilization of Cranes 1 and 2 is 57.09% and 61.25%, respectively. In SJF, utilization of Cranes 1 and 2 is 69.33% and 74.12%, respectively. In NNF, utilization of Cranes 1 and 2 is 52.47% and 53.70%, respectively. In LSOM, utilization of Cranes 1 and 2 is 98.71% and 100.00%, respectively. Compared with FIFS, SJF, and NNF, LSOM increases utilization of Cranes 1 and 2 by 41.62%, 38.75%, 29.38%, 25.88%, 46.24%, and 46.30%, respectively, which corresponds to increase rates of 72.90%, 63.27%, 42.38%, 34.92%, 88.13%, and 86.22%.
Table 5. Crane utilization using FIFS, SJF, NNF, and LSOM
CraneUtilization (%)
FIFSSJFNNFLSOM
157.0969.3352.4798.71
261.2574.1253.70100.00

Discussion

The proposed LSOM with dynamic programming based on motion characteristics of luffing tower cranes has demonstrated its effectiveness in an example of 12 lifting tasks and two cranes. The example involves only lifting sequence optimization, which is the research boundary of this study. LSOM can be integrated in planning, such as layout optimization, as well.
The luffing mode of luffing tower cranes makes the hoist motion of the hook dependent on the luffing motion of the boom, and the hoist motion distance of the hook depends not only on the height difference between the lifting start and end points but also on the vertical motion distance of the boom. In construction projects, lifting task start and end points typically differ in horizontal distance and height. This is particularly evident in large-scale and high-rise construction projects, enhancing the impact of the hoist motion of the hook on total lifting time. Thus, the connection between hoist motion distance and boom luffing motion has been taken into account in the proposed formulation.
Estimation of hoist motion distance of the hook is classified in two ways: one way is based on the respective horizontal distance between the crane and the lifting start point and between the crane and the lifting end point; the other way is based on the height difference between the lifting start and end points and the boom vertical motion distance. This classification can estimate the hoist motion distance of the hook more effectively and reduces the risk of overestimating and underestimating the actual situation.
For construction sites with limited space, lifting delay caused by avoiding space conflicts among luffing tower cranes is common and often has a great impact on total lifting time. The horizontal motion path of the boom created by the relative locations of the crane and the lifting start and end points can result in delays, so it is essential to evaluate the impact of delays accurately. In the proposed formulation, the boom’s horizontal motion path is divided into three stages: lifting start point, motion path, and lifting end point. The intersection of these stages determines the relative locations of two lifting tasks, which have been classified as four situations consisting of 16 cases. According to the distinction between loaded and no-load lifting, the specific lifting delays for 16 cases have been detailed and quantified to refine estimation of lifting delay time and make them more reliable.
In construction projects with loose time, conventional lifting strategies such as FIFS, SJF, and NNF provide adequate lifting-sequence solutions. However, in projects with limited time, these strategies are not guaranteed. Therefore, dynamic programming has been used here to achieve the global optimum of multistage decision-making in lifting sequence optimization. By analyzing the lifting sequence characteristics of luffing tower cranes, the lifting process is divided into stages according to lifting sequence, and a lifting task and a crane are determined at each stage in terms of minimum total lifting time.
The example results have shown that LSOM is significantly better than FIFS, SJF, and NNF in reducing total lifting time and improving crane utilization. Although its computation takes more time when compared with the other strategies in the example, this problem can be effectively solved using compiled languages, such as C++, and parallel computing. The pros and cons of each strategy could be weighed to select the most suitable one when specific total lifting time, crane utilization, and computation time are required.
The proposed formulation for loaded and no-load lifting time is based on ideal environmental conditions. In practice, bad weather and complex site characteristics may result in increased lifting difficulty and longer lifting time. The method of multiplying lifting time by numerical parameters to account for impact factors may be useful. The determination of impact factors, although beyond the scope of this study, can be made according to data from similar construction projects.
The proposed formulation for lifting delay time is presented in the basic form of two lifting tasks. In practice, however, no matter how many lifting tasks there are, a space conflict always results from the intersection of two luffing tower cranes. Thus, when there are multiple potential space conflicts among more than two lifting tasks, the situation can be dealt with according to the basic form of two lifting tasks. However, when a lifting task faces more than one space conflict, estimation of lifting delay time becomes more complex, and extensions need to be investigated.

Conclusion

Considering the importance of the lifting sequence of luffing tower cranes in construction projects, a lifting sequence optimization model based on motion paths has been proposed to minimize total lifting time. The model consists of loaded lifting time, no-load lifting time, and lifting delay time, integrating these parameters into dynamic programming. Specifically, it (1) quantitatively evaluates the influence of the boom’s luffing motion on the hoist motion distance of the hook; (2) refines the solution to various lifting delays resulting from the boom’s horizontal motion path based on different relative locations of the crane and the lifting start and end points; and (3) integrates the lifting sequence optimization problem and dynamic programming. The effectiveness of LSOM is evaluated by comparing it with FIFS, SJF, and NNF. Total lifting time using LSOM decreases by 65.05, 37.61, and 77.39 min, respectively, which corresponds to decrease rates of 42.17%, 29.65%, and 46.45%. Utilization of Cranes 1 and 2 in the example increases by 41.62%, 38.75%, 29.38%, 25.88%, 46.24%, and 46.30%, respectively, which corresponds to increase rates of 72.90%, 63.27%, 42.38%, 34.92%, 88.13%, and 86.22%. The results show that LSOM effectively decreases total lifting time and increases crane utilization.
In some cases, lifting tasks have fixed lifting sequences or special time limits, which affects the motion path of luffing tower cranes, so corresponding constraints should be added to lifting sequence optimization. For example, the shoring system needs to be lifted before the structural steel is lifted, and concrete needs to be lifted within a specific period of time. The proposed model has to be adjusted when considering such conditions. This issue will be addressed in future work.

Data Availability Statement

Some or all data, models, or code that support the findings of this study are available from the corresponding author upon reasonable request.

References

Abdelmegid, M. A., K. M. Shawki, and H. Abdel-Khalek. 2015. “GA optimization model for solving tower crane location problem in construction sites.” Alexandria Eng. J. 54 (3): 519–526. https://doi.org/10.1016/j.aej.2015.05.011.
Al Hattab, M., E. Zankoul, and F. R. Hamzeh. 2017. “Near-real-time optimization of overlapping tower crane operations: A model case study.” J. Comput. Civ. Eng. 31 (4): 05017001. https://doi.org/10.1061/(ASCE)CP.1943-5487.0000666.
Al Hattab, M., M. Zankoul, M. Barakat, and F. Hamzeh. 2018. “Crane overlap and operational flexibility: Balancing utilization, duration, and safety.” Constr. Innovation 18 (1): 43–63. https://doi.org/10.1108/CI-11-2016-0062.
Al-Hussein, M., M. A. Niaz, H. Yu, and H. Kim. 2006. “Integrating 3D visualization simulation for tower crane operations on construction sites.” Autom. Constr. 15 (5): 554–562. https://doi.org/10.1016/j.autcon.2005.07.007.
Bellman, R. 1957. Dynamic programming. Princeton, NJ: Princeton University Press.
Bürgy, R., A. Hertz, and P. Baptiste. 2020. “An exact dynamic programming algorithm for the precedence-constrained class sequencing problem.” Comput. Oper. Res. 124 (Dec): 105063. https://doi.org/10.1016/j.cor.2020.105063.
de Matos, V. L., A. B. Philpott, and E. C. Finardi. 2015. “Improving the performance of stochastic dual dynamic programming.” J. Comput. Appl. Math. 290 (Dec): 196–208. https://doi.org/10.1016/j.cam.2015.04.048.
Grevtsov, N. M., S. A. Kumakshev, and A. M. Shmatkov. 2017. “Optimization of the flight trajectory of a non-manoeuvrable aircraft to minimize fuel consumption by the dynamic programming method.” J. Appl. Math. Mech. 81 (5): 368–374. https://doi.org/10.1016/j.jappmathmech.2018.03.004.
Huang, C., and C. K. Wong. 2018. “Optimization of crane setup location and servicing schedule for urgent material requests with non-homogeneous and non-fixed material supply.” Autom. Constr. 89 (May): 183–198. https://doi.org/10.1016/j.autcon.2018.01.015.
Huang, C., C. K. Wong, and C. M. Tam. 2011. “Optimization of tower crane and material supply locations in a high-rise building site by mixed-integer linear programming.” Autom. Constr. 20 (5): 571–580. https://doi.org/10.1016/j.autcon.2010.11.023.
Hwang, S. 2012. “Ultra-wide band technology experiments for real-time prevention of tower crane collisions.” Autom. Constr. 22 (Mar): 545–553. https://doi.org/10.1016/j.autcon.2011.11.015.
Ibrahim, A., O. El-Anwar, and M. Marzouk. 2018. “Socioeconomic impact assessment of highly dense-urban construction projects.” Autom. Constr. 92 (Aug): 230–241. https://doi.org/10.1016/j.autcon.2018.04.001.
Irizarry, J., and E. P. Karan. 2012. “Optimizing location of tower cranes on construction sites through GIS and BIM integration.” J. Inf. Technol. Constr. 17 (23): 351–366.
Ji, Y., and F. Leite. 2018. “Automated tower crane planning: Leveraging 4-dimensional BIM rule-based checking.” Autom. Constr. 93 (Sep): 78–90. https://doi.org/10.1016/j.autcon.2018.05.003.
Konak, A., M. R. Bartolacci, and B. Gavish. 2011. “A dynamic programming approach for batch sizing in a multi-stage production process with random yields.” Appl. Math. Comput. 218 (4): 1399–1406. https://doi.org/10.1016/j.amc.2011.06.022.
Lazarev, A. A., and F. Werner. 2009. “A graphical realization of the dynamic programming method for solving NP-hard combinatorial problems.” Comput. Math. Appl. 58 (4): 619–631. https://doi.org/10.1016/j.camwa.2009.06.008.
Lee, D., H. Lim, H. Cho, and K.-I. Kang. 2015. “Optimization of luffing-tower crane location in tall building construction.” KICEM J. Constr. Eng. Project Manage. 5 (4): 7–11. https://doi.org/10.6106/JCEPM.2015.5.4.007.
Leung, A. W. T., and C. M. Tam. 1999. “Prediction of hoisting time for tower cranes for public housing construction in Hong Kong.” Constr. Manage. Econ. 17 (3): 305–314. https://doi.org/10.1080/014461999371510.
Leung, A. W. T., C. M. Tam, and D. K. Liu. 2001. “Comparative study of artificial neural networks and multiple regression analysis for predicting hoisting times of tower cranes.” Build. Environ. 36 (4): 457–467. https://doi.org/10.1016/S0360-1323(00)00029-9.
Li, B., B. Cai, H. Liu, F. Wang, and J. Zhang. 2012. “Research on anti-backward-tipping device of luffing jib tower crane.” In Vol. 229–231 of Proc., Applied Mechanics and Materials, 521–525. Bäch, Switzerland: Trans Tech Publications. https://doi.org/10.4028/www.scientific.net/AMM.229-231.521.
Lien, L.-C., and M.-Y. Cheng. 2014. “Particle bee algorithm for tower crane layout with material quantity supply and demand optimization.” Autom. Constr. 45 (Sep): 25–32. https://doi.org/10.1016/j.autcon.2014.05.002.
Marzouk, M., and A. Abubakr. 2016. “Decision support for tower crane selection with building information models and genetic algorithms.” Autom. Constr. 61 (Jan): 1–15. https://doi.org/10.1016/j.autcon.2015.09.008.
Monghasemi, S., M. R. Nikoo, and J. Adamowski. 2016. “Sequential ordering of crane service requests considering the pending times of the requests: An approach based on game theory and optimization techniques.” Autom. Constr. 70 (Oct): 62–76. https://doi.org/10.1016/j.autcon.2016.06.006.
Moussavi Nadoushani, Z. S., A. W. A. Hammad, and A. Akbarnezhad. 2017. “Location optimization of tower crane and allocation of material supply points in a construction site considering operating and rental costs.” J. Constr. Eng. Manage. 143 (1): 04016089. https://doi.org/10.1061/(ASCE)CO.1943-7862.0001215.
Rodriguez-Ramos, W. E., and R. L. Francis. 1983. “Single crane location optimization.” J. Constr. Eng. Manage. 109 (4): 387–397. https://doi.org/10.1061/(ASCE)0733-9364(1983)109:4(387).
Shapira, A., and M. Simcha. 2009. “Measurement and risk scales of crane-related safety factors on construction sites.” J. Constr. Eng. Manage. 135 (10): 979–989. https://doi.org/10.1061/(ASCE)CO.1943-7862.0000066.
Tam, C. M., A. W. T. Leung, and D. K. Liu. 2002. “Nonlinear models for predicting hoisting times of tower cranes.” J. Comput. Civ. Eng. 16 (1): 76–81. https://doi.org/10.1061/(ASCE)0887-3801(2002)16:1(76).
Tam, C. M., and K. L. T. Tong. 2003. “GA-ANN model for optimizing the locations of tower crane and supply points for high-rise public housing construction.” Constr. Manage. Econ. 21 (3): 257–266. https://doi.org/10.1080/0144619032000049665.
Tam, C. M., K. L. T. Tong, and W. K. W. Chan. 2001. “Genetic algorithm for optimizing supply locations around tower crane.” J. Constr. Eng. Manage. 127 (4): 315–321. https://doi.org/10.1061/(ASCE)0733-9364(2001)127:4(315).
Wang, J., X. Zhang, W. Shou, X. Wang, B. Xu, M. J. Kim, and P. Wu. 2015. “A BIM-based approach for automated tower crane layout planning.” Autom. Constr. 59 (Nov): 168–178. https://doi.org/10.1016/j.autcon.2015.05.006.
Warszawski, A. 1973. “Analysis of transportation methods in construction.” J. Constr. Div. 99 (1): 191–202. https://doi.org/10.1061/JCCEAZ.0000373.
Wu, K., and B. García de Soto. 2020. “Spatiotemporal modeling of lifting task scheduling for tower cranes with a tabu search and 4D simulation.” Front. Built Environ. 6: 79. https://doi.org/10.3389/fbuil.2020.00079.
Wu, K., B. García de Soto, and F. Zhang. 2020. “Spatio-temporal planning for tower cranes in construction projects with simulated annealing.” Autom. Constr. 111 (Mar): 103060. https://doi.org/10.1016/j.autcon.2019.103060.
Younes, A., and M. Marzouk. 2018. “Tower cranes layout planning using agent-based simulation considering activity conflicts.” Autom. Constr. 93 (Sep): 348–360. https://doi.org/10.1016/j.autcon.2018.05.030.
Zavichi, A., K. Madani, P. Xanthopoulos, and A. A. Oloufa. 2014. “Enhanced crane operations in construction using service request optimization.” Autom. Constr. 47 (Nov): 69–77. https://doi.org/10.1016/j.autcon.2014.07.011.
Zhang, P., F. C. Harris, and P. O. Olomolaiye. 1996. “A computer-based model for optimizing the location of a single tower crane.” Build. Res. Inf. 24 (2): 113–123. https://doi.org/10.1080/09613219608727511.
Zhang, P., F. C. Harris, P. O. Olomolaiye, and G. D. Holt. 1999. “Location optimization for a group of tower cranes.” J. Constr. Eng. Manage. 125 (2): 115–122. https://doi.org/10.1061/(ASCE)0733-9364(1999)125:2(115).

Information & Authors

Information

Published In

Go to Journal of Construction Engineering and Management
Journal of Construction Engineering and Management
Volume 147Issue 10October 2021

History

Received: Dec 9, 2020
Accepted: Apr 6, 2021
Published online: Aug 4, 2021
Published in print: Oct 1, 2021
Discussion open until: Jan 4, 2022

Authors

Affiliations

Keyi Wu, Ph.D. [email protected]
Postdoctoral Associate, S.M.A.R.T. Construction Research Group, Div. of Engineering, New York Univ. Abu Dhabi, Abu Dhabi 129188, United Arab Emirates (corresponding author). Email: [email protected]
Borja García de Soto, Ph.D., M.ASCE https://orcid.org/0000-0002-9613-8105
P.E.
Assistant Professor, S.M.A.R.T. Construction Research Group, Div. of Engineering, New York Univ. Abu Dhabi, Abu Dhabi 129188, United Arab Emirates. ORCID: https://orcid.org/0000-0002-9613-8105

Metrics & Citations

Metrics

Citations

Download citation

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

Cited by

  • Multiuser Virtual Reality-Enabled Collaborative Heavy Lift Planning in Construction, Journal of Construction Engineering and Management, 10.1061/JCEMD4.COENG-14102, 150, 4, (2024).

View Options

Media

Figures

Other

Tables

Share

Share

Copy the content Link

Share with email

Email a colleague

Share