Minimize the Route Length Using Heuristic Method Aided with Simulated Annealing to Reinforce Lean Management Sustainability

Cost reduction is a cornerstone of the Lean administration’s sustainability through modify its algorithms scheme to become multi-useful. This paper focuses on control “movement” waste, to minimize pipeline, cabling and sewerage network deployments time, to avoid demurrages (i.e., constructor sectors) and quickens planning through two stages. The first belongs to the build constrained hybridization of published heuristic routing methods (e.g., S-Shape, Mid-point, Largest-Gap, Return, Ascending, FLA-5, FLA-6 [Flow Line Analysis], and Composite) to select the shortest path that serves many locations (i.e., Plan-A), while allowing for the modification of these locations during service (i.e., Plan-B). The new locations are grouped into two clusters, the first of which lay on the shortest preferred path, while the second cluster contains locations that do not lay on the preferred path and are therefore moved on the backlogs-list, then use Simulated Annealing when to serve them. Finally, the impact of the selected performance is investigated after studying its correlation with another published effective one under cost considerations. The computational results of proposed Minimize-Route-Length aided with simulated annealing (MRL-SA) significantly outperform others in terms of the performance of the routing heuristics and total costs and develop the Last Planner System, which has a good reputation in construction projects and approve the proposed algorithm to maintain its competitiveness sustainability.


Introduction
There are many flow-line analysis methods reported in the literature. Among these, there are eight famous methods that have many interdisciplinary citations. Two of the best methods have been selected in this paper via the correlation test. They are Flow Line Analysis (FLA)-5, S-shape, Mid-point, Ascending-which are effective in tackling bidirectional flow lines-and FLA-6, Return, Largest Gap and Composite-which are effective in tackling unidirectional flow lines, while combining FLA-6 and Composite methods with cost analysis. A hybrid search algorithm is proposed by combining the effective heuristic methods according to correlation outperformance. The optimum solution S may not serve all the required locations, entailing the creation of a backlogs-list that serves at the end of plan-A with a Simulated Annealing SA algorithm, to find out the optimum route length. Thirty-seven sets of data are created, which have 185 different examples and have been tested with the proposed algorithm Minimize-Route-Length aided with simulated annealing (MRL-SA) to enhance a planner system that has a good reputation in construction sectors, namely the Last Planner System, which has been developed to become the Lean Planner System, which can enhance this sector decision making.
The construction industry is overwhelmed by delays and has often suffered cost and time overrun. Alsehaimi and Koskela [1] reported that weak project management and competitiveness can be traced back to the absence of the implementation methodology and procedures, which was a dominant reason for the delay in the construction projects sector internationally, due to waste of cycle time and task variability. The Egypt Engineering syndicate indicated that up to 25% of construction is reworked due to poor planning, which loses between 40%-50% of potential competitiveness efficiency of total plan costs, and at least 10% of materials are wasted. Azevedo et. al., 2012 [2], depict how the idea mapping procedure, which joins subjective and quantitative strategies that communicate the capability of lean assembling, adds to an organization's natural and social supportability. The effect of incline toward profitability and procedure productivity was distinguished, yet the outcomes positively affect resource utilization as Lean recommended. conventional construction management ramparts had been stormed through construction academics and professionals to enhance plans output value to owners. Therefore, lean-intervene tools have emerged, and have been successfully implemented to enhance profits in Egypt. Remon Fayek Aziz et. al. (2013) [3] discussed the principles, methods, and implementation phases of lean in construction, showing waste type and how it could be minimized. Amongst the main waste which consumes time, cost and planning efforts are the drilling path and pipeline networks deployment. The study proposes using the Last Planner System (LPS) technique as a prevalent application of the lean construction concepts and methodologies. The proposed approach enhances the LPS by presenting suggestions to reduce the costs and efforts (pointed in Figure 1) of drilling or pipeline deployment in construction sector, and this is called the Lean Planner System (LPS), as illustrated in Figure 2. Effective utilization plans for various resources are essential to the competitiveness level. Mohammad Hamdy Elwany et al. (2013) [4] reported an integrated heuristics-based solution methodology that tackles assignment problems. A study by Dassonville L., et al., 2016 [5] was entitled "The challenges experienced in Lean implementation to maintain a sustainable management system". Laila M. Khodeir 2016 [6] looked at the interactions among lean and manageability (i.e., sustainability) standards on the administration procedures of planning and development ventures. To accomplish this, two methodologies were utilized, namely a literature review, and a correlation matrix to check the territory of association between both lean construction and sustainability standards. The survival of any organization depends upon its competitive edge. Even though Lean is one of the most powerful quality management improvement methodologies, according to Rodriguez D, et. al., 2016 [7], nearly two-thirds of the Lean implementations had been failures, and less than one-fifth of those implemented had sustained results. One of the most significant tasks of top management is to identify, understand and deploy the significant Lean practices through different sectors and applying its tools, such as quality circle, Kanban, just-in-time purchasing, etc., as cited by Cherrafi A, et. al., 2016 [8].
The objective of this paper is to find the routes that represent the shortest path cost, one that starts its trip and visits all locations only once and returns to the I/O point. This paper focuses on creating two sequential plans, Plan-A and Plan-B, by creating the random locations ψ i and ψ j using pairwise exchange, insertion and random insertion schemes. The objective values of the three schemes (i.e., minimum route length and its related execution in terms of time or cost) are compared by applying the proposed heuristic method that is supported with the SA (Simulated Annealing) technique. Plan-A means implementing published heuristic methods that have a minimum route length for specific locations, while Plan-B means modifying locations (i.e., adding or removing locations is allowed) and recreates another route that combines a new route with a previous one. If the new locations lay on the same direction of selected route, then no further actions are required, but if it is not on the tracking route, a backlogs list for non-serviced locations must be created, which is tackled after finishing the proposed route activities using the SA method. The researchers developed a new procedure by embedding a local search heuristic procedure in the SA technique to enhance its performance. Conversely, this work embedding SA to enhance the heuristic method, which found a better performance within the terms of the solution (i.e., route length) and the time taken for the execution of tasks. The SA technique was proposed to improve the cellular manufacturing (CM) systems [9]. The proposed algorithm was tested by 185 different and complex layout locations. The results show the superiority of the proposed SA-supported algorithm over the mathematical programming model in both route length and its related cost. Achieving these objectives supports the sustainability of the Lean management competitiveness [10].
Processes 2019, 7, x FOR PEER REVIEW 5 of 27 full-turnover assignment policy (no backtrack), which sorts all locations to be allocated.  There is another vision to switch the proposed plan from A to B. This divides the layout area into zones, which have similar treating proposed procedures, and will be assigned close to each other. However, this vision showed no significant impact on the performance. full-turnover assignment policy (no backtrack), which sorts all locations to be allocated.  There is another vision to switch the proposed plan from A to B. This divides the layout area into zones, which have similar treating proposed procedures, and will be assigned close to each other. However, this vision showed no significant impact on the performance. A hybrid algorithm was proposed to solve a single row layout design problem with unequal machines clearances [11]. The proposed algorithm was tested on many well-known data sets, and we compared the results. The proposed algorithm was found to have good effectiveness and showed improvements in solving layout problems. Genetic algorithms were proposed to determine a common linear sequence for multi-products with different operation sequences and facilities [12]. The proposed algorithm can be used to reduce the number of machines in the layout, the material handling distance and material handling cost. The Tabu-search algorithm was used to solve single row facility layout [13] problems and find the linear sequence in order to minimize the total material handling cost. The algorithm found a common linear machine layout for more than one product [14]. The selecting routing heuristics performance and efficiency were affected by the layout and the number of locations, especially in inventory material handling of the size requested, and batching methods affected by route length. In [15,16], the authors developed a Tabu-search heuristic for Vehicle Routing and Scheduling Problems. The authors considered vehicles with different capacities, and the access to some customers is restricted to particular vehicles [17,18]. Moreover, drivers must respect the maximum legal driving time per day. In their subsequent work [19], the authors adapt the algorithm to compare the results with those obtained by [20].
This paper aims to close this research gap by studying the published routing policies and how, by starting with Plan-A at specific time t, and another at t', the route could be modified into another strategy, Plan-B, according to the requested list L. Clearly, such deviations from the original route could significantly improve the performance of the initial routing policy, and perhaps even result in a situation where heuristic routing leads to better results in practice.

Serviceable Sequence Actions for Route Direction
The proposed plan aims to select the minimum route length between main successive locations ψ i of extensions (such as electrical transformers, water traps, etc.) that are related in terms of the time spent performing each activity, as illustrated in Figure 1. There are five sequences of activities to minimize the route length and time (moving to specific location (i.e., travel), search about priority direction, cost locations, serviceable actions, setup if necessary, and other sudden reactions).
Blair 2013, states that the performance and efficiency of any sequential activities' plan depends on four tactical decisions [21]: layout design and its selected locations, assignment policy, priority facilities movement and routing policy-these four are a seed for creating the proposed algorithm shown in Figure 3.
The Classification of Proposed systems emphasizes that layout design is one important factor affecting the processes performance and routing distance. It takes n times to finish n (locations); the numbers of trips to serve all locations are equal to the number of locations requested. The factors to be considered in internal layout include: number of locations, length and width of avenue, number and shape of cross avenue and position of main extension source (I/O point). Is there a cross avenue, yes or no? If yes, how many (I/O) locations? Rear and front cross avenue. There are three types of layouts. The first is a parallel avenues case with an I/O-Point located in the middle, the second and third are Vertical Avenues, but the I/O-point is located in the middle and lower left corner to sketch a Flying-V or fishbone. There is another movement direction in ascending order, according to the full-turnover assignment policy (no backtrack), which sorts all locations to be allocated.
There is another vision to switch the proposed plan from A to B. This divides the layout area into zones, which have similar treating proposed procedures, and will be assigned close to each other. However, this vision showed no significant impact on the performance.

Problem Statement in Lean Planner System
The aim of this work is to develop a heuristic algorithm for discovering the shortest route length to serve some locations that lie on a specific layout. The proposed heuristic method aided with SA (MRL-SA) is combined to maximize the number of exchanges in the sequence for the locations in the backlogs-list, as described above. The locations in the backlogs-list have two classifications, unidirectional and bidirectional, which could not support Plan-A/Plan-B when the added locations were not on the proposed routes. Some 185 examples of different locations distribution were generated for unidirectional (i.e., clockwise) and bidirectional (i.e., counterclockwise) movement by choosing interaction values from 1:9 using uniform distribution.
A residential compound contains a large number of building and avenues that need to deploy infrastructure cables and hydrant piping via the sequential procedures list, illustrated in Figure 1,

Problem Statement in Lean Planner System
The aim of this work is to develop a heuristic algorithm for discovering the shortest route length to serve some locations that lie on a specific layout. The proposed heuristic method aided with SA (MRL-SA) is combined to maximize the number of exchanges in the sequence for the locations in the backlogs-list, as described above. The locations in the backlogs-list have two classifications, unidirectional and bidirectional, which could not support Plan-A/Plan-B when the added locations were not on the proposed routes. Some 185 examples of different locations distribution were generated for unidirectional (i.e., clockwise) and bidirectional (i.e., counterclockwise) movement by choosing interaction values from 1:9 using uniform distribution.
A residential compound contains a large number of building and avenues that need to deploy infrastructure cables and hydrant piping via the sequential procedures list, illustrated in Figure 1, taking into account, the proposed plan created and updated in two aspects, plan-A created according to MRL. In the proposed planning considered in this paper, a contractor faces dynamic stochastic [9] extensions for serviceable items (locations) ψ i through Plan-A (electrical transformers/hydrant pipe or other) and the dynamic stochastic modification of Plan-A to become plan-B ψ j in each period over a finite time horizon. The probability distributions of locations under service, whether at Plan-A or at Plan-B, may be disposed of the specific layout if they are off the route. Only the mean values and the maximum deviations from the mean values can be estimated in each period. The contractor must decide on the number of serviceable items ψ i , and if they are serviceable at the beginning of each period t. Once the serviceable items ψ i in period t have been observed, the distance to the serviceable items at the end of period t, which is equal to the distance to the serviceable items at the end of period t-1, plus the new distance that belongs to the locations at Plan-A (lpA: location on Plan-A route ψ i ), and are serviceable in period t less realized ψ i in period t, is obtained. If there are two choices to select direction when modifying Plan-A ψ i to switched to Plan-B ψ j , the route length yield to cost analysis, which should be carried out for the proposed change of direction [16], and the route processing cost per item is incurred. If this is neglected and continued to implement Plan-A (i.e., backlogging), costs for items will be added at the end of Plan-A. The contractor can install or dispose the items inserted in Plan-B (lpB ψ j ) and must decide on the number of items that are (re)serviceable or disposed of at the beginning of each period t + 1. A variable disposal cost per item disposed is incurred when the contractor disposes the items of lpB ψ i . Once the lpB items in period t have been observed, the route length to the lpB items at the end of period t-which is equal to the route length of the lpB at the end of period t-1 plus the realized serviceable lpB items in period t, minus the (re)serviceable or disposed items in period t, according to the cost analysis, is obtained. If the ending route length of the (re)serviceable items are in two different directions, a route processing cost per item is incurred (operational cost-saving analysis) to make a decision. Note that the ending route of the (re)serviceable items cannot be negative. Without loss of generality, we assume that serviceable, (re)serviceable and disposable items from Plan-A or Plan-B occur instantaneously through executing Plan-A. In the following, important parameters and decision variables that are used in selecting the route direction after starting at t are given. In order to change the listed locations during the execution of the selected proposed plan, a push to create Plan-B is necessary (i.e., some of the serviceable locations may be added while others may be disposed of through executing Plan-A). The problem assumptions are as follow:

Structure and Layout
Shape: Rectangular Number of blocks: 2 < X < 2000 Main source point location: Lower left corner Avenue width: As OSHA instructed Service assignment: Random/proposed Proposed policy: Single or groups Zoning: None Proposed system: Stochastic and dynamic ψ i, plan A or ψ j, plan B Electrical transformers, hydrants valves, etc.

The Layout Assumption
The maximum capability for the study is based on these assumptions. The understudy consists of 2000 locations, J = 40 avenues running between the back and front and N = 100 locations per avenue, each avenue is two-sided and has 100 ψ i and can be serviceable by being published on each side. For simplicity, the published random assignment policy in this paper deals with locations, not the number of ψ i in every location vs. the potential bias, which is likely to occur in the location choices cited by Yilun Shang 2019 [22].

Parameters of Proposed Cost Analysis Plan for Branching Case
If the locations are added for Plan-A to create Plan-B and recalculate the route direction, the first step assumes the method for processing cost analysis to select the minimum cost direction. If the cost=N cost=0 Plan B , selects Plan-A items ψ i branches and neglects Plan-B items ψ j to backlogging after finishing all locations in Plan-A proposals ψ i , or applies hybrid modification of items into Plan-B ψ j and adopts its execution according Figure 3. C PA : Cost of distance drilling (route length) to item ψ i locations in Plan-A requested. C PB : Cost of distance drilling (route length) to item ψ j locations in Plan-B that executes through Plan-A and in its route. C dpB : Cost of distance drilling for disposing locations ψ j that serviced plan-A or were canceled during Plan-B. h PA : Item processing cost to serviceable ψ i at locations in plan-A for requested list L at period t. h PB : Item processing cost to serviceable ψ i locations appeared at plan-B, if in the same proposed route. b : Item backlogging cost to serviceable ψ j locations at plan-B after finishing Plan-A.
Number of locations ψ i of plan-A at the beginning of the proposed planning horizon at time (t). I R 0 : Number of locations ψ j of plan-B at the beginning of the proposed planning horizon at time (t').
r t : Uncertain number of ψ j locations in plan-B and serviceable or hybrid at plan-A in period t + ∆t.
Estimate of the mean number of ψ i + ψ j locations updated to serviceable in period t. r t : Maximum deviation from the mean published number of ψ i locations updated to serviceable in period t. d t : Uncertain serviceable locations' requested for the list L in period t. d t : Estimate of the mean requested list L in period t. d t : Maximum deviation from the mean requested list L in period t. τ d t , τ r t : Maximum number of uncertain requested and uncertain (re)request # of locations (route length) that can simultaneously deviate from their mean (i.e., d t and r t ) until the end of period t.

Decision Variables
S B t : # of serviceable ψ j locations at Plan-B that is updated of Plan-A in period t'. X D t : # of disposed items from the list L in period t' due to cancelation or against cost analysis step and backlogged to end of the route activities. I S t : # of backlogging locations ψ i+ j for requested all locations at the end of period t of Plan-A and Plan-B.

The Contribution of This Study
The question is, how to make a hybrid between different methods to achieve the maximum benefits of the research objective, when implemented through Lean Planner System. 1-. Proposed new routing heuristic method [23] to avoid demurrage and control those costs. 2-. Investigating the impact of proposed regions on the performance of routing heuristics.
The proposed layout consists of: a fixed two dimensional rectangular inside has only one location (J avenue x N two-sided locations per avenue), with one source feeding (I/O-point) located at the extreme lower left corner, as shown in Figure 4. The two end avenues are located near and parallel to the back and front walls of the boundaries. The definition of each symbol on the layout is shown at each location and illustrated in Table 1.
The locations of the ψ is indicated as (the first location in any avenue j) where j {1,2, . . . ,J} is (j − 1)N + 1, the second location in any avenue j is (j − 1)N + 3, . . . , (j − 1)N + i where i {1,2, . . . ,N} and so on, to reach the last location at any avenue j is j N , but the routing is restarted when ∆S becomes a straight distance back and returns to the same point to start next trip of connection. The distance between any two locations x in j x and y in j y is the shortest distance which are calculated as: The traversed distance (backtrack), the current avenue to the back avenue, then going to the next serviceable location ψ i . 2.
The distance if the returned to the front avenue, then going the next serviceable location.
Then follow the minimum distance.
The distance of serviceable all locations in Plan-A is calculated as shown in the next equation [24]: where:   Table 1. Symbols used in the proposed plan model illustrated in Figure 3.

x, y
Coordinates of serviceable items such as electrical transformers or hydrants or any similar items.

ΔB
The distance between the last ψ , ψ in the vertical avenue and the back avenue.
j Avenue number on suggest layout plan ΔA The distance between two adjacent avenues.

N
Number of serviceable targeted locations at each avenue whether in Plan A/plan B

J
Total number of avenues, whether horizontal, vertical, or diagonal ΔF The distance between the first location and the front avenue. on, to reach the last location at any avenue j is jN, but the routing is restarted when ΔS becomes a straight distance back and returns to the same point to start next trip of connection. The distance between any two locations x in jx and y in jy is the shortest distance which are calculated as:  Table 1. Symbols used in the proposed plan model illustrated in Figure 3.
x, y Coordinates of serviceable items such as electrical transformers or hydrants or any similar items.

∆B
The distance between the last ψ i , ψ j in the vertical avenue and the back avenue.
j Avenue number on suggest layout plan ∆A The distance between two adjacent avenues.

Number of serviceable targeted locations at each avenue whether in Plan A/plan B ∆l
The distance between two adjacent locations: ψ i , ψ j .

J
Total number of avenues, whether horizontal, vertical, or diagonal ∆F The distance between the first location and the front avenue.

Cost Analysis Formulation
Subject to: where The objective function (2) is the total test cost in the branch case for a new location added to Plan-A but not on the suggested route, to a serviceable specific location or if disposed. Constraints (7) and (4) keep track of the distance at the end of period t or back locations at the end of period t for any possible realization of ψ i+ j . Note that is the capability level of the service requested ψ i+ j at the end of period t. Constraints (9) ensure that the capability level is nonnegative at the end of period t. Constraints (10) are for the non-negativity of variables. The variables q t , ϑ it , w t , ε it , β T and ∅ tT for 1 ≤ i ≤ t ≤ T are used to transform the uncertain model into a traceable model. Next, we show that there are flaws in constraints (9) and (13) of the model, given the above.
To reformulate the model illustrated in (2)-(10) to a traceable model and to control the uncertainty behavior, the right-hand sides of variables S A t , S B t , X D t , I S t shown in the proposed cost analysis are in step with constraints (15)- (16). The proposed step inferred, as deduced from Ben-Tal et al. (2004), the following inequality: The constraint (11) represents an uncertain parameter, and α t for 1 ≤ t ≤ T is an expression involving some decision variables, which can equivalently be reformulated as The proposed step will yield a component-disposal policy that adapts the ψ i+ j disposal location, based on the realization of requested and updated requests. The ψ i+ j disposal policy, which will be feasible regardless of the realization of requested and updated requests, is determined by the realized ψ i+ j as requested and updated, and by the values of

Lean Planner System-Phase I [Proposed MRL Heuristic Method]
The minimum route length plan is the main objective of this study. This method has one main drawback, which has a negative impact on its performance. This drawback is that when the route must move from the current avenue to the next one, it just follows the ascending high frequency locations of the ψ i+ j list. To eliminate this drawback, the proposed solution will decide whether to traverse its avenue or return back through the front avenue, based on the shortest routing with minimum cost analysis. The proposed plan has been updated during the connection mission, where some of locations were removed the list and others were returned to the requested list during implementation. The updated ψ j request will be classified to get into Plan-A implementation or disposed, to avoid obstruction or downtime during the implementation mission. This explanation reveals that the problem must be processed as an uncertainty case related to time. The detailed steps of proposed heuristic are as shown in next algorithm.

Step1:
Arrange the requested list L of (Plan-A) according to the serviceable locations ψ i in an ascending according to specific layout.

Step2:
The route starts from the feeding source point to the first location ψ i in L

Step3:
The serviceable all ψ i in (Plan-A) list is subject to one of previous published methods (such as: S-Shape, Mid-point, Return, L-Gap, Return, Composite or Ascending) and uses the same method that have minimum route to service ψ j in (Plan-B) that id in the same proposed direction and stops at the last one ψ i+ j in the last avenue, then goes to step 10.

Step4:
Resolve serviceable ψ j (which are on the proposed route) items in (Plan-B) and considers last ψ i is the starting point, and keeps the previous calculated distance route in variable δ b = λ i via the same published method stated in step 3.

Step5:
If there are ψ j not on the proposed route, a backlogs list is created for them Step6: λ im is the distance between the current location I and m, where m is the first serviceable location in the next avenue, according to the list L of (Plan-A). λ jn is the distance between the current location j and n, where n is the last serviceable location in the next avenue, according to the list L of (Plan-B), which is based on resolving the new distribution of serviceable ψ i+ j on the same route.

Step7:
If the ψ j is not on the (Plan-A) route, and needs to reverse its direction, it must be subject to the step 7 analysis.

Step8:
If cost=N cost=0 Plan A ≤ cost=N cost=0 Plan B , select (Plan-A) items ψ i branch and neglect (Plan B) items ψ j to backlogging after completing all (Plan-A) proposal items ψ i , or carry out hybrid modification of items into (Plan-B) ψ j and merge the suggested routes.

Step9:
If min [(λ im > λ jn + ∆F)] < y, reverse the sequence of the next avenue in list L.
Step10: Continue serviceable the next item in L, then go to step 3 again.
Step11: Go back to the front avenue, then go to the feeding source point.

Performance Measure
The performance is a measurable response that leads to an understanding of efficiency and effectiveness (checking and handling the location right), which is important to determine the total route distance λ d i , which is required for the service of all ψ and make a comparison with the minimum route, according to the published methods, and uses the deviation %Dv ∝ = where it represents the percentage of deviation for the suggested route distance ∝ from the minimum routing distance, and λ d i ∝ is the total distance of the routing heuristic ∝ .

Illustrative Example for implementing MRL-SA (ψ i Distribution)
This example explains how the route distance is calculated based on the MRL-SA heuristic, assuming the compound have five avenues, each of which have 16 two-sided building locations, and the feeding source point is located at lower extreme left corner, and the black cells are the locations of ψ i to be serviced, as shown in Figure 5. The requested list L is required for serviceable locations # as follows (80, 67, 16, 29, 54, 84, 75, 23, 19, 8, 50, and 12) at period t for Plan-A and 1, 2, 4, 5, 6, 7, 41, 63, at period t' for Plan-B. Where, (∆F = ∆B = 1 unit, ∆A = 3 units, ∆D = 0 unit, ∆l = 1 unit), by applying the distance equations, which are described above. During the implementation according to time, the ψ 41, 63 is located on the proposed route, and therefore will serviceable according to the cost and distance analysis, but other added locations shall be saved in the backlogs list L b . The proposed MRL-SA tackles L b via the simulated annealing method after arriving at the end location in the proposed route. Step2: The plan aims to go from the feeding source point to the location ψ 3 (the distance λ i = 3 units) Step3: The plan aims to service all ψ i in the current avenue, and stop at the last ψ 16 in the same avenue in this illustrative example (the distance λ i = 3 + 5 + 4 + 4 = 16 units). Step4: Calculates the distance between the location (ψ 16 , ψ 19 ), which is 20 units, and (ψ 16 , ψ 29 ), which is 10 units to arrive at the same location (Plan-A). [Take into account costs analysis equations] Step5: The direction will reverse locations in this avenue because 10 < 20, and go to location ψ 29 (the total distance λ i = 16 + 10 = 26 units). Step6: The proposed plan will collect ψ i from locations ψ 29 , ψ 23 , ψ 19 and the total travelling distance at end of this avenue is λ i = 16 + 10 = 26 units. Step7: At t', the location ψ 41 is added to the layout and placed on the track route. Therefore, there is no need to use cost analysis equations and continue with the MRL method. Step8: Calculates the distance between locations (ψ 19 , ψ 41 ), which is 11 units, while the distance between (ψ 19 , ψ 50 ) is 20 units. [Take into account costs analysis equations] Step9: The track route will follow the ascending arrangement and go to location 4, because 11< 20.

only six units. [Take into account costs analysis equations]
Step12: The direction will reverse the ascending arrangement, because 6<21, and go to location ψ 67 (the route distance λ i = 57 + 6 = 63 units). Step13: At t', location ψ 63 is added to layout (Plan B), but not effects on direction because it between two terminal locations (ψ 54 , ψ 67 ). Therefore, this location will be serviced on the track route. Step14: The proposed route will serve locations ψ 67 , ψ 63 , ψ 54 , and the route distance will reach λ i = 65 + 4 + 9 = 78 units.
Step18: The proposed route direction will return to the front avenue and back to the main upstream source distance λ i = 100 + 16 = 116 units. Step19: The proposed route will return to the feeding source point and the route distance λ i = 116 + 12 = 128 units.
This solution is done via implementing MRL (i.e., Plan-A), and create a backlog-list for evolving mandates for new locations (Plan-B), which will tackle by Simulated Annealing SA to enhance the solution.

Lean Planner System Phase-II: Backlogs List Tackling Via SA
The optimum sequence for the backlogs list was obtained for both the bidirectional and unidirectional route. In this work, a seed sequence layout route was obtained by MRL (i.e., based on hybridizing superiority heuristic methods). The SA method increases the performance of implementation by finding out the best sequence of backlogs list locations by computing the objective value, as given below. Plan-B, modifies Plan-A at t' by inserting new locations ψ i , before and after the proposed method starts at time t, which the backlogs list tackles via SA and rearranges the locations ψ j randomly, then measures the route length and related cost, then adds that to the previous data obtained by MRL. In this work, the initial temperature of 500 • C was chosen via expert. The solution can be accepted with an associated probability of 0.92, when it was inferior by 60%, relative to the original solution (expected by proposed heuristic method). Therefore, the initial setting was followed as given below: P = exp (−λ ∆ i /T), i.e., 0.92 = exp(−60/T) or T = 500 • C, Then the SA algorithm starts at a temperature of 500 • C, with temperature reduction factor of 0.92. In this work, the numbers of temperature stages were taken as 42. Then the final temperature was, T f = 500 * (0.92) 42 = 20 • C.
The initial seed sequence appeared during implementing Plan-A or a serviceable Plan-B on route were 2-1-6-3-4-7-5, according to t'. If the initial temperature = 500 • C, ACCEPT = 0 and TOTAL = 0 are set as 12 different sequences from the initial backlogs list sequence. Objective values for the 12 sequences were computed, and ∆ length/unit (i.e., double the distance between the objective and real route length) were calculated as shown in Table 2. Among the twelve ∆ route length/unit, −20 was the least value, its corresponding sequence was taken as preferred sequence 2-7-1-6-3-4-5, its objective value was 36, and the TOTAL and ACCEPT increased by one.

MRL-SA Verification
This section focuses on analyzing the results of the 185 different implemented distributions of different locations for all the published methods that are compared with proposed heuristic, to find a relationship between them and a similarity level to create verification level for the proposed method. Also, the objectives of the enterprises were collected according to a generic framework [25] using five-point Likert scales, which revealed four independent variables as Losses cost, route length, traveling time and processing time. Table 3 Table 3 (i.e., the means and standard deviations of different heuristics), and was found to be 4.54 (S.D = 0.205), 4.49 (0.402), 4.21 (0.0518) and 3.69 (0.562), respectively, for enterprises objectives. The frequencies were calculated to test the normality of data and it was found that the data used in present research were normally distributed. Select the value above average in the minimum or maximum columns, and the minimum standard deviation is less than the average and maximum mean. The impact factors of measurement objectives illustrated in Table 4 which prelude to correlation which emerge the hybridization level via correlation analysis which discussed in Table 5. Table 5 illustrates the effective of Composite method. Therefore, this work interests with compare the results of many examples between MRL-SA and Composite as illustrate in Table 6.   In this work, two heuristic search algorithms, one combining MRL with SA, and the other combining composite with SA, which was selected from a previous correlation test, were proposed and compared together, and the result of the comparison is illustrated in Table 6.
If the relationship between Plan-A and Plan-B is positive according to impacted deployment shown in Table 4 (as A improved or was implemented in minimum time, B also improved in a vice versa negative relationship), the magnitude of Cor(A, B) measures the strength of the relationship. The closer Cor(A, B) is to 1 or -1 and the sign indicates the direction of the relationship. If the Cor(A, B) > 0, this implies that B and A are positively related. Conversely, Cor(B, A) < 0 implies they are negatively related. If p-value (p < 0.0001) is much less than α = 0.05, this indicates a very high significance. In this table, we analyzed the four types of enterprise needs separately to clarify their affect independently, as illustrated in Table 5.
Thirty-seven sets of input data were generated from the uniform distribution from one to nine for the unidirectional flow path, as well as the bidirectional flow paths [8]. The above combination of the selected heuristic method and the simulated annealing method was tested in the 185 sets of generated data.

Computational Verification Results
The computational results are presented as follows, based on ANOVA and statistical analysis shown in Tables 7 and 8 for measuring the different impacts on the performance of the routing technique if the layout to the parallel sectors such as upper/lower or right/middle/left, is divided. Figure 6 illustrates the average percentage deviation of each routing. It is noticeable that the best routing is MRL-SA without considering the requested size for every location, while Figure 7 shows the interaction between the routing policies by changing the requested locations and controlling the number of ψ i (i.e., until 1300 locations in the layout) at a specific time related to the other locations, according to the published plan. Figure 7 illustrates that the largest gap is best at low duty ψ i until a 4% modification of Plan-A, whether MRL is better than others if it has more than 160 locations from 1300. This value push to improve the proposed heuristic published plan divides into two areas, the upper half, which consists of regions (C, G, K, D, H, L), and the lower half, which consists of regions (A, E, I, B, F, J) and three vertical areas, as shown in Table 9 and illustrated in Figure 8, while the ANOVA results are illustrated in Table 10.     Table 9.
Computational results of the proposed model were solved by a solver with its default settings on a hp-core i3 Precision Workstation running under Windows 10. The I/O place selection in the left, right or middle will not be significant for the proposed heuristic, but S-Shape is preferable if it is in  Table 9.
Computational results of the proposed model were solved by a solver with its default settings on a hp-core i3 Precision Workstation running under Windows 10. The I/O place selection in the left, right or middle will not be significant for the proposed heuristic, but S-Shape is preferable if it is in the right place and the composite method is preferable if it is in the middle, as illustrated in Figures 9-11. Also, the layout is divided to Qi zones. The initial deviation of the request locations in Plan-A, and the modified locations in Plan-B are set to zero (i.e., I S 0 = I R 0 = 0). the location requested, updated and disposing are taken as S A = 7, S B = 4, X D = 2, respectively. h PB = 4 ≤ h PA = 5 The locations backlogging is set to b ∈ {3, 4, 5, 6}, as shown in Table 8, to ensure all possible relations (i.e., less than, equal, or greater than) between h PB and b, h PA and b. The uncertain request has a mean d t ∈ {18, 20, 22} and a standard deviation of σ ∈ {2, 4}. In the proposed framework, the uncertain request d t can take any value from the interval [d t − 2σ, d t + 2σ], i.e.,d t = 2σ. Similarly; the uncertain updated requests have a mean r t ∈ {14, 16, 18} and a standard deviation σ ∈ {2, 4}. As a result, the uncertain updated locations r t can take any value from the interval [r t − 2σ, r t + 2σ], i.e.,r t = 2σ. Thus, combining all the parameter settings that have generated 72 instances for the average-case performance assessment, in total. While the costs of uncertainty τ d t are generated such that the probability of infeasibility at constraints (4) and (5) for t is less than 5%, the proposed costs of uncertainty τ r t are generated such that the probability of infeasibility at constraint (6) for t is less than 0.1%.   Figure 11. Also, the layout is divided to Qi zones. The initial deviation of the request locations in Plan-A, and the modified locations in Plan-B are set to zero (i.e., = = 0). the location requested, updated and disposing are taken as = 7, = 4, = 2, respectively. ℎ = 4 ≤ ℎ = 5 The locations backlogging is set to b ∈ {3, 4, 5, 6}, as shown in Table 8, to ensure all possible relations (i.e., less than, equal, or greater than) between ℎ , ℎ . The uncertain request has a mean ̅ ∈ {18,20,22} and a standard deviation of ∈ {2,4}. In the proposed framework, the uncertain request can take any value from the interval [ − 2 , ̅ + 2 ], . . , = 2 . Similarly; the uncertain updated requests have a mean ̅ ∈ {14,16,18} and a standard deviation ∈ {2,4}. As a result, the uncertain updated locations can take any value from the interval [ − 2 , ̅ + 2 ], . . , ̂ = 2 . Thus, combining all the parameter settings that have generated 72 instances for the average-case performance assessment, in total. While the costs of uncertainty are generated such that the probability of infeasibility at constraints (4) and (5) for t is less than 5%, the proposed costs of uncertainty are generated such that the probability of infeasibility at constraint (6) for t is less than 0.1%.     . Comparison under middle area classification. Figure 11. Comparison under middle area classification.

Conclusions
In this work, the Last Planner System used in the construction projects sector has evolved by a hybrid search algorithm that was used by combining MRL and SA through two sequential phases, and the sum of the route length value (e.g., illustrated in Figure 12) of the main route that results from MRL and backlogs-list of the unidirectional or bidirectional flow path to become the Lean Planner System. The performance of the proposed algorithm was tested using a wide set of data, and the results showed that the proposed MRL-SA algorithm was performing well, as illustrated in Figure 13. The main conclusions for the paper are summarized in the following points:

Conclusions
In this work, the Last Planner System used in the construction projects sector has evolved by a hybrid search algorithm that was used by combining MRL and SA through two sequential phases, and the sum of the route length value (e.g., illustrated in Figure 12) of the main route that results from MRL and backlogs-list of the unidirectional or bidirectional flow path to become the Lean Planner System. The performance of the proposed algorithm was tested using a wide set of data, and the results showed that the proposed MRL-SA algorithm was performing well, as illustrated in Figure  13. The main conclusions for the paper are summarized in the following points:

Conclusions
In this work, the Last Planner System used in the construction projects sector has evolved by a hybrid search algorithm that was used by combining MRL and SA through two sequential phases, and the sum of the route length value (e.g., illustrated in Figure 12) of the main route that results from MRL and backlogs-list of the unidirectional or bidirectional flow path to become the Lean Planner System. The performance of the proposed algorithm was tested using a wide set of data, and the results showed that the proposed MRL-SA algorithm was performing well, as illustrated in Figure  13. The main conclusions for the paper are summarized in the following points: When the requests reach up to 60% of all locations that must be serviceable, the proposed heuristic MRL-SA presents the optimum route length, which has not proved costly when compared to the reviewed heuristic methods, but if the request size is considered, the Composite heuristic When the requests reach up to 60% of all locations that must be serviceable, the proposed heuristic MRL-SA presents the optimum route length, which has not proved costly when compared to the reviewed heuristic methods, but if the request size is considered, the Composite heuristic method is preferred to achieve the objective of the proposed Plan-A. Therefore, Table 6 illustrates the results of the analyses and a comparison between the proposed MRL-SA and the composite heuristic method. The published three policies (S-Shape, Largest Gap and Mid-point) have nearly the same performance. There are two cases affecting MRL-SA performance, first when the serviceable locations fall to less than 4.2%, while the second emerged when the requested locations for serviceable becomes more than 14% of the total requested using S-shape or MRL-SA, and these will be integrated together to achieve the objective.
The paper revealed that when the layout was divided into two halves (upper and lower) and the farthest points about the I/O position, the integration of the MRL-SA and Composite is preferred as illustrated in Figure 13 and can be deduced from Table 8. While, if the requested locations lay in the lower half, the integration between MRL-SA and the Return or Mid-point are preferred to execute Plan-A and Plan-B, as illustrated in Table 10.
When the requested serviceable locations are distributed vertically in three parallel regions, the Plan's implementation of MRL-SA or another method will result in a short route, but with differences in cost, according to the requested distribution case.

Future Work
The same proposed algorithm that targets a reduction in the movement waste of Lean management, can be used for other applications. For example, improving the processes layout that has minimum backtrack direction and minimum transportation, whether in shop-floor or in inventory handling, which has been widely implemented as cited by Joël, R [26], but nowadays, this algorithm can be programmed to be a Mobile application to guide patients in hospitals to reduce their walks between examining rooms. This policy targets a reduction in the likelihood of patients meeting in aisles, which is especially relevant in the context of combat COVID19, as suggested by the WHO recommendations.