Joint Location–Allocation Model for Multi-Level Maintenance Service Network in Agriculture

: The maintenance service network is always designed as a multi-level service network to provide timely maintenance service for failed machinery, and is rarely studied in agriculture. Thus, this paper focuses on a three-level maintenance service network location–allocation problem in agriculture, which contains several spare part centres, service stations, and service units. This research aims to obtain the optimal location of spare part centres and service stations while determining service vehicle allocation results for service stations, and the problem can be called a multi-level facility location and allocation problem (MLFLAP). Considering contiguity constraints and hierarchical relationships, the proposed MLFLAP is formulated as a mixed-integer linear programming (MILP) model integrating with P-region and set covering location problems to minimize total service costs, including spare part centre construction costs, service vehicle usage costs, and service mileage costs of service stations. The Benders decomposition-based solution method with several improvements is then applied to decompose the original MLFLAP into master problem and subproblems to ﬁnd the optimal solutions effectively. Finally, a real-world case in China is proposed to evaluate the performance of the model and algorithm in agriculture, and sensitivity analysis is also conducted to demonstrate the impact of several parameters.


Introduction
Agricultural machinery, as key and general resources in agricultural production, are distributed over a wide geographical area in China [1].Agricultural machinery is prone to failure due to high temperatures, complicated environments, and long-term agricultural operations, especially in the busy season [2,3].To maintain stable production in agriculture, manufacturers always design a maintenance service network to supply maintenance for failed machinery [4,5].
In general, the maintenance service network is always designed as a multi-level service network in China, which contains several service facilities, e.g., a spare part centre and service station.This paper focuses on the design problem of the multi-level maintenance service network.In this research, the location of spare part centres and service stations is chosen from the potential location set while determining service vehicle allocation results, and the service relationship between two adjacent levels of service facilities is defined simultaneously, which can also be described to define these service stations' service area.In summary, the problem considered in this research can be described as a multi-level facility location-allocation problem in the maintenance service network (MLFLAP).
Numerous studies on service network design have been conducted and applied in different areas, such as supply chain networks [6,7], healthcare and emergency systems [8][9][10], maintenance service networks [11], and transportation and logistics systems [12,13].Service paradigms in agriculture bring new issues and motivate new research in the design of maintenance service networks.First, in agriculture, there is no cross-level service in the multi-level maintenance service network.That is, spare part centres are built to supply spare parts to service stations, and service stations are in charge of maintenance service for failed machinery.Thus, spare part centres cannot supply maintenance services for failed machinery directly, which is different from other service networks.In addition, different from discrete demand points in previous studies, maintenance demands in agriculture are spread across a wide geographical area.Thus, manufacturers always cluster some agricultural production areas into several service units, and maintenance demands in service units are calculated based on all failed machinery in these production areas.All service demands in the same service unit can be served by one service station.The region-based service region districting problem is a special problem in agricultural maintenance.Lastly, an adjacent production area always means a similar climate and crop-ripening time, so service units in connected areas allocated to the same service station can improve service efficiency.Contiguity constraints should be considered in this research to ensure that service units in a service region and service district are contiguous.
Thus, this paper considers the location-allocation problem in the multi-level maintenance service network, and this research aims to develop a practical solution approach to determine the location of spare part centres and service stations while determining the allocation results of service vehicles.
The contributions of this research are threefold.First, in agriculture, service stations and spare part centres need to be located considering hierarchical relationships while determining the service area of these service facilities.A defined number of service vehicles is expected to be allocated to these selected service stations to provide maintenance service for failed machinery simultaneously.Thus, a multi-level facility location-allocation problem (MLFLAP) in agriculture is proposed in this paper.Second, an MILP-based solution model is formulated for minimal service costs.Contiguity criteria are introduced to ensure that the service region served by one service station is contiguous, and multi-level constraints are applied to prohibit cross-level service.Lastly, a Benders decomposition approach with several improvements is developed to address the large-scale cases in reality efficiently, and sensitivity analysis is conducted for decision making in practical application.

Literature Review
Facility location-allocation problems (FLAPs) have been a research focus for several years and are applied to determine the location of service facilities and the allocation of service resources [14,15].FLAPs have been applied in different areas, including emergency service systems [16], supply chain networks [17], transportation and logistics systems [18], urban health and services systems [19], and maintenance service networks [20].
In classical single-level FLAPs, the location of service facilities is selected and the allocation of service resources is also defined.However, there always exist multi-level facilities in service networks, and multi-level facility location-allocation problems (MLFAPs), as a subclass of FLAPs, have attracted more attention in the last two decades [15,21].For example, Doyen et al. [22] concentrated on the two-echelon humanitarian relief logistics system and developed a Lagrangian relaxation-based heuristic algorithm to determine the location of pre-and post-disaster rescue centres.Wu and Chu et al. [23] conducted research to determine the location of plants and warehouses in reality and combined the Lagrangian relaxation method and simulated annealing algorithm to obtain minimal transportation and production costs.Shu and Wu et al. [24] considered the location problem of warehouses and retailers in a two-echelon supply chain network and developed a novel cutting-plane approach to solve the real case with moderate size.Shavarani et al. [25] focused on the facility location problem in multi-level delivery systems with retailers, refuel stations, and warehouses, and developed a novel solution method to minimize transportation costs considering demands uncertainty and M/G/K queueing system.Abbassi et al. [26] considered a multi-objective optimization model to determine the location of intermediate healthcare facilities and two heuristic algorithms, including particle swarm optimization and a non-dominated sorting genetic algorithm, were designed to solve several case studies in healthcare supply chain logistics.
As for mathematical models for FLAPs, they are always formulated as mixed-integer program (MIP) models.For example, Zamani et al. [27] developed a mixed-integer programming model to solve congested facility location problems considering machine failures.Muffak et al. [28] considered the fixed and dynamic demand and developed a mixed-integer linear programming model to determine the location of public facilities.Govindan et al. [29] focused on facility location problems in food supply chain networks and formulated a multi-objective model to determine the number and location of service facilities.Talaei et al. [30] constructed an MIP-based optimization model to determine the location of collection/inspection centres and minimize the total costs in a closed-loop green supply chain system.Ghasemi et al. [16] developed a multi-objective optimization model to minimize total allocation costs while minimizing the total amount of relief suppliers.Tirkolaee et al. [31] considered the service facility location problem in urban service systems and formulated an MILP for minimal costs based on the robust optimization approach.
Solution approaches in FLAPs can be classified into two groups: exact approaches and heuristic algorithms [32,33].As for exact approaches, Sonmez et al. [34] designed a decomposition-based solution approach to obtain a near-optimal location plan and compared it with another exact method to demonstrate the performance of the developed approach.Fischetti et al. [35] combined the Benders decomposition method and commercial solver with the brand and cut method to obtain the best service facility location, and computational results were presented to show the method's efficiency compared with exact methods and heuristic algorithms.Irawan et al. [36] combined the integer linear programming method and aggregation approach to obtain the optimal location of distribution centres.Several experiments were conducted to demonstrate the efficiency of the developed method compared with other exact approaches.Han and Zhang et al. [37] developed a Benders decomposition-based exact approach to determine the location of the service facility and service region in each stage.Chandra et al. [38] developed an exact algorithm integrating with the convexification strategy and branch and cut method for solving larger cases in the real world.Christensen et al. [39] focused on the capacitated facility location problem.Lagrangean relaxation was applied for the nonlinear objective function, and the branch and bounded method was then introduced to obtain an optimal location plan with a lower bounding scheme.Qi et al. [40] proposed a brand-and-brunch-based solution method, and two valid inequalities were introduced with an approximation algorithm to obtain the optimal results.
Considering heuristic algorithms for FLAPs, Meng et al. [41] developed a hybrid genetic algorithm integrating with a logarithmic-quadratic proximal prediction-correction method to obtain the optimal location of retailers in a decentralized supply chain system.Jolia et al.
[42] designed an improved particle swarm optimization method to solve the multi-objective in dynamic facility location problem.Rahmaniani et al. [43] proposed a MIP-based optimization model, and an improved firefly algorithm-based heuristic approach was designed to solve a larger case with more than 100 nodes and 600 candidates.Halper et al. [44] concentrated on the mobile facility location problem and developed two local search neighbourhood algorithms to solve two subproblems in reality.Al-Rabiaah et al. [45] developed a heuristic solution approach integrating the greedy search and maximum coverage approach to determine the location and allocation of drone launching centres.Kang [46] considered these NP-hard problem and proposed a heuristic solution method based on maximum flow to find the optimal location.In addition, other heuristic algorithms were also applied to obtain the optimal location of a service facility in the service network, e.g., the simulated annealing algorithm [47], the kernel search heuristic algorithm [48], the tabu search method [49], and the greedy search heuristic approach [50].
In all, although several studies have been conducted on FLAPs, few studies exist on multi-level FLDPs in maintenance service networks, especially in agriculture.In addition, unlike in traditional studies, maintenance demands in connected areas are always clustered as a service unit in agriculture, and few studies consider facility location problems based on service district.Therefore, it is urgent to develop a more comprehensive solution framework considering contiguity criteria and multi-level constraints in maintenance service networks in agriculture.

Problem Description
To supply timely maintenance for failed machinery in the busy farming seasons, manufacturers always construct a multi-layer maintenance service network.There exist three layers in the service network, which include spare part centres, service stations, and service units.Service units are a set of connected production areas including agricultural machinery, just like a county in China, and a connected service area always means a similar climate and ripening time, which represents maintenance demands.Service stations with several service vehicles are built by manufacturers to supply maintenance for the demand for service units.
Thus, the multi-layer facility location-allocation problem (MLFLAP) can be described as follows.There are n service units with defined maintenance demands in agriculture.In order to supply maintenance, l service stations are selected from the potential location set of service stations to supply maintenance service, and m service vehicles with defined service capacity are allocated to the selected service stations.In addition, several spare part centres are built by the manufacturer to supply spare parts for service stations.
Assuming that there are 20 service units, four service stations are selected to supply service for these service units, and nine service vehicles are allocated to the selected service stations.In Figure 1, service stations are located in service units 4, 8, 14, and 17, and the service units served by each service station are distinguished by bold lines.As for service station 4, there exists one service vehicle that can provide maintenance service for four service units, e.g., 3, 4, 5, and 10.As for service station 17, two service vehicles are allocated to this service station, and five service units, such as 11, 12, 16, 17, and 18, are served by service station 17.In addition, an indefinite number of spare part centres are built by the manufacturer to supply spare parts for service stations, which are served by the nearest spare part centres.As shown in Figure 1, in order to supply spare parts for service stations, two spare part centres are selected in service units 12 and 9, and the maximal distribution distance for these spare part centres is shown with black dotted lines.Thus, the spare part centre in 9 can supply spare parts for service stations 4, 8, and 14, while the spare part centre in 12 can supply spare parts for service stations 8 and 17.
Thus, the proposed MLFLAP can be divided into several subproblems: (1) Service station location and allocation problem in agriculture, which is to decide the location of service stations while allocating service vehicles to these stations.(2) Spare part centre location problem, which is to determine the location of spare part centres.(3) Service unit allocation problem, which is to assign these service units to selected service stations.The objective of the proposed MLFLAP is to minimize total costs, which include total service mileage costs from service stations to service units, total usage costs of service vehicles, and total construction costs for spare part centres.
Several assumptions are presented in this research: (1) The maintenance requirement in service units is known.
(2) One service unit can only be served by one service station.
(3) Several service units served by one service station are continuous.
(4) There is a defined service capacity for service vehicles, and there exists a maximal distribution distance for spare part centres to service stations.Thus, the proposed MLFLAP can be divided into several subproblems: (1) Service station location and allocation problem in agriculture, which is to decide the location of service stations while allocating service vehicles to these stations.(2) Spare part centre location problem, which is to determine the location of spare part centres.(3) Service unit allocation problem, which is to assign these service units to selected service stations.The objective of the proposed MLFLAP is to minimize total costs, which include total service mileage costs from service stations to service units, total usage costs of service vehicles, and total construction costs for spare part centres.
Several assumptions are presented in this research: (1) The maintenance requirement in service units is known.
(2) One service unit can only be served by one service station.
(3) Several service units served by one service station are continuous.
(4) There is a defined service capacity for service vehicles, and there exists a maximal distribution distance for spare part centres to service stations.

Nomenclature
As for the proposed MLFLAP, two subproblems, such as the service station location and allocation problem and the spare part centre location problem, are considered in this paper, which can be seen as a joint optimization problem integrating with the P-region problem and set covering location problems.Thus, this paper formulates an MILP-based solution model.Several parameters and variables are defined first: (1) Sets: The number of service units; m The number of service vehicles; l The number of chosen service stations; Distance from the potential location of service station j to service units i;

Nomenclature
As for the proposed MLFLAP, two subproblems, such as the service station location and allocation problem and the spare part centre location problem, are considered in this paper, which can be seen as a joint optimization problem integrating with the P-region problem and set covering location problems.Thus, this paper formulates an MILP-based solution model.Several parameters and variables are defined first: (1) Sets: (2) Parameters: The number of service units; m The number of service vehicles; l The number of chosen service stations; Distance from the potential location of service station j to service units i; D k j Distance between the potential location of service station j and potential location of spare part centre k; c 1 The construction costs per spare part centre; c 2 The usage costs per service vehicle; c 3 The service mileage costs for service stations to service units per kilometre; r i Maintenance demand requirements in each service unit i; D1 max Maximal service distance for service stations to service units; D2 max Maximal service distance for spare part centres; C The maximal service capacity per service vehicle; M A constant more than the number of service units.
(3) Variables: x k If the spare part centre is set in location k, the value is 1; otherwise, it is 0; y j If the service station is set in location j, the value is 1; otherwise, it is 0; If service unit i is served by the potential location of service station j, the value is 1; otherwise, it is 0; v j The number of service vehicles allocated to potential locations j; u k j If the potential location of service station j can be served by the potential location of spare part centre k, the value is 1; otherwise, it is 0; The contiguity evaluation for optimal location-allocation results.The variable for contiguity evaluation to indicate the imaginary flow from service unit i to service unit i , which is served by service station j.

Mathematical Model
The MILP-based optimization model P 0 is used to obtain the minimal service costs, which is shown as follows: The objective function (1) aims to obtain minimal service costs, which include total construction costs for spare part centres, total usage costs for service vehicles, and total service mileage costs for service stations: Constraints ( 2) and ( 3) ensure that all these chosen service stations can be served by the spare part centres.Constraint (4) implies the maximal service distance for spare part centres.Constraint (5) restricts the number of service vehicles allocated to service stations, and constraint (6) ensures the service capacity of service stations.Constraint (7) defines the number of chosen service stations.Constraints ( 8) and ( 9) specify that one service unit can only be served by one service station.Constraints (10) and ( 11) introduce a limitation that all service units must served by the chosen service stations.Constraint (12) enforces the maximal service distance for the service station.Constraints ( 13)-( 15) refer to Han and Zhang [37] and are extended for contiguity constraints in service regions served by service stations.Constraint ( 13) is applied to calculate the imaginary flow from the service unit to the chosen service stations.Constraints ( 14) and ( 15) are introduced to prohibit the imaginary flow from other service stations.Constraints ( 16)-( 18) specify the range of decision variables.

MP and SP in Benders Decomposition
Considering two integer variables w and v in the proposed MLFLAPs, a variant Benders decomposition, which is called compound Benders decomposition in this paper, is introduced to decompose the original MLFLAPs into the master problem (MP) and two subproblems (SP 1 and SP 2 ).In MP, the location of service stations and spare part centres are determined in this paper using the decision variables x, y, z, and u.As for SP 1 , decision variables v are introduced to determine the service vehicle allocation, and in SP 2 , decision variables w are introduced for contiguity evaluation of service units served by one service station.

MP and SP in Benders Decomposition
As for decision variables v in SP 1 , they are applied to obtain minimal service vehicle usage costs and are subject to constraints ( 5)-( 6) and (17).Thus, the dual SP 1 for variable v is shown as DSP 1 .
Thus, during the Benders decomposition approach, the obtained facility location results and service assignment plan will be introduced to DSP 1 to solve the DSP 1 .If there is no feasible solution for DSP 1 , there is no feasible solution for the original MLFLAPs, which means the number of service vehicles cannot meet all maintenance demands in the service network.If the DSP 1 is unbounded, there exist optimal solutions for original MLFLAPs, and the generated Benders feasibility cuts in (26) are added in MP.If the DSP 1 is feasible and bounded, the generated Benders optimality cuts in (27) are added in MP.
As for SP 2 , the facility location results and service assignment results obtained in MP are introduced for contiguity evaluation.If the SP 2 returns a feasible solution, the service units served by service stations are contiguous.Otherwise, it demonstrates that there exists at least one service unit which is not contiguous with other service units served by one service station.Selected locations of service stations that are not contiguous are denoted as J*.Thus, as for these service units served by these service stations in J*, at least one service unit should be removed or one new service unit should be added, which can be described as shown in constraint (28).Thus, the combinatorial Benders cuts in (28) are generated and added to MP.

The Framework of Benders Decomposition
Thus, the Benders decomposition approach is applied to obtain the location-allocation results in a multi-level maintenance service network.As shown in Figure 2, the detailed information of the developed solution algorithm is described as follows.
Step 2. MP is solved to obtain a feasible solution (x, y, z, u), and LB (Lower Bounded) is updated.
Step 3. The obtained service station location results y and service assignments plan z are introduced to DSP 1 .If the DSP 1 is infeasible, there are no results for the original problem.If the DSP 1 is unbounded, the generated Benders feasibility cuts (26) are added in MP, and a feasible solution is found to update UB (Upper Bounded).And if the DSP 1 is feasible and bounded, the generated optimality cuts (27) are added in MP, and UB is updated.
Step 4. The obtained service station location results y and service assignments plan z are introduced to SP 2 .If it returns a feasible solution, service units served by service stations are contiguous.Otherwise, the generated Benders cuts (28) are added in MP.
Step 5. Calculate the difference between UB and LB, and determine if SP 2 is feasible.If SP 2 is feasible and the values of UB and LB are equal, the optimal results of MLFLAPs are obtained.Otherwise, repeat Step 2 to Step 4. Step 1. Input initial data.
Step 2. MP is solved to obtain a feasible solution ( , , , x y z u ), and LB (Lower Bounded) is updated.
Step 3. The obtained service station location results y and service assignments plan z are introduced to DSP1.If the DSP1 is infeasible, there are no results for the original problem.If the DSP1 is unbounded, the generated Benders feasibility cuts (26) are added in MP, and a feasible solution is found to update UB (Upper Bounded).And if the DSP1 is feasible and bounded, the generated optimality cuts ( 27) are added in MP, and UB is updated.
Step 4. The obtained service station location results y and service assignments plan z are introduced to SP2.If it returns a feasible solution, service units served by service stations are contiguous.Otherwise, the generated Benders cuts (28) are added in MP.
Step 5. Calculate the difference between UB and LB, and determine if SP2 is feasible.If SP2 is feasible and the values of UB and LB are equal, the optimal results of MLFLAPs are obtained.Otherwise, repeat Step 2 to Step 4.

Case Study
In this section, a real case in East China is conducted to demonstrate the efficiency of the proposed mathematical model and solution algorithm.The real maintenance data in three provinces, e.g., Shandong, Anhui, and Jiangsu, is shown first and the optimal location-allocation plan is then obtained in the next part.Finally, several analyses are conducted to verify the application and performance of the developed mathematical model and solution approach.

Case Study
In this section, a real case in East China is conducted to demonstrate the efficiency of the proposed mathematical model and solution algorithm.The real maintenance data in three provinces, e.g., Shandong, Anhui, and Jiangsu, is shown first and the optimal location-allocation plan is then obtained in the next part.Finally, several analyses are conducted to verify the application and performance of the developed mathematical model and solution approach.

Case Description
In this part, the maintenance service network is built, including three provinces in East China, namely Shandong, Anhui, and Jiangsu, and the maintenance data are obtained from a leading agricultural machinery manufacturer in China.As shown in Figure 3, there exist 251 service units in agriculture, and maintenance demands in these service units are shown in Figure 4.The manufacturer wants to construct 25 service stations in these 251 service units.In addition, the manufacturer can provide 75 service vehicles to be allocated to these service stations.The maximal service capacity per service vehicle is 300, and the usage costs per service vehicle are CNY 800.
In order to satisfy the maintenance demands efficiently, the maximal service distance for service stations is set as 200 km.As for spare part centres, the construction costs for per spare part centre are set as CNY 10,000, and the maximal distribution distance for spare part centres is set as 360 km.In addition, the travelling distance between two service units is obtained from Baidu Map, and the service mileage costs (from service stations to service units) are set as CNY 0.5.
East China, namely Shandong, Anhui, and Jiangsu, and the maintenance data are obtained from a leading agricultural machinery manufacturer in China.As shown in Figure 3, there exist 251 service units in agriculture, and maintenance demands in these service units are shown in Figure 4.The manufacturer wants to construct 25 service stations in these 251 service units.In addition, the manufacturer can provide 75 service vehicles to be allocated to these service stations.The maximal service capacity per service vehicle is 300, and the usage costs per service vehicle are CNY 800.In this part, the maintenance service network is built, including three provinces in East China, namely Shandong, Anhui, and Jiangsu, and the maintenance data are obtained from a leading agricultural machinery manufacturer in China.As shown in Figure 3, there exist 251 service units in agriculture, and maintenance demands in these service units are shown in Figure 4.The manufacturer wants to construct 25 service stations in these 251 service units.In addition, the manufacturer can provide 75 service vehicles to be allocated to these service stations.The maximal service capacity per service vehicle is 300, and the usage costs per service vehicle are CNY 800.

Optimal Location-Allocation Plan
In this section, the Benders decomposition solution approach is developed to solve the proposed MLFLAPs, and all computing programs are coded on a personal computer with an AMD 4.00 GHz processor.
First, the commercial solver GUROBI is introduced to solve the proposed MLFLAPs directly, but the optimal results can only be obtained after about 20 h (71,933 s) due to the massive decision variables and constraints, which may not be applicable in reality.The developed Benders decomposition approach is then applied to this problem.The optimal solutions are obtained with three iterations in 8727 s (2.4 h), and the minimal service costs in this case is 5.836 × CNY 10 5 .The location-allocation results of service stations and spare part centres are shown in Table 1 and Figure 5.As shown in Figure 5, 25 service stations and two spare part centres are built in this service network.First, service stations are usually selected in these service units with a higher number of service demands, such as service units 1, 12, 79, 85, 106, 138, 153, 202, and 238.Service units served by one service station are distinguished by different colours, and all 251 of these service units are always allocated to the nearest service stations.For example, the service station in 73 provides maintenance service for three service units, e.g., 73, 76, and 78.Otherwise, the service station in 220 provides maintenance service for seventeen service units, e.g., 3, 178, 206, 215, 216, 218, 220, 221, 222, 223, 239, 240, 241, 242, 243, 244, and 246.As shown in Figure 5, 25 service stations and two spare part centres are built in this service network.First, service stations are usually selected in these service units with a higher number of service demands, such as service units 1, 12, 79, 85, 106, 138, 153, 202, and 238.Service units served by one service station are distinguished by different colours, and all 251 of these service units are always allocated to the nearest service stations.For example, the service station in 73 provides maintenance service for three service units, e.g., 73, 76, and 78.Otherwise, the service station in 220 provides maintenance service for seventeen service units, e.g.In addition, there are 60 service vehicles are allocated to these selected service stations based on the service demands in their service area.For example, four service vehicles are allocated to the service station in 12 with a total maintenance demand of 1168, while only one service vehicle is allocated to the service station in 94 with a total maintenance demand of 246.
As for spare part centres, two spare part centres are set in service units 118 and 207, and the service area of these two spare part centres is shown in Figure 5 with the black dotted lines.They are responsible for supplying spare parts for service stations in the top and bottom areas in Figure 5, respectively.In this case, service stations in 8, 10, 12, 43, and 138 can be served by these two spare part centres.In addition, there are 60 service vehicles are allocated to these selected service stations based on the service demands in their service area.For example, four service vehicles are allocated to the service station in 12 with a total maintenance demand of 1168, while only one service vehicle is allocated to the service station in 94 with a total maintenance demand of 246.
As for spare part centres, two spare part centres are set in service units 118 and 207, and the service area of these two spare part centres is shown in Figure 5 with the black dotted lines.They are responsible for supplying spare parts for service stations in the top and bottom areas in Figure 5, respectively.In this case, service stations in 8, 10, 12, 43, and 138 can be served by these two spare part centres.

Sensitivity Analysis of Selected Parameters
In this section, sensitivity analysis is conducted to evaluate the effect of several parameters on the final results, which include the following: (a) the number of service stations, l; (b) the maximal service distance for service stations, D1 max .
First, this paper conducts a sensitivity analysis of the number of service stations, l.Several examples are conducted to find the minimum value of l, and the computational results show that there is no optimal solution when the value of l is less than 9 due to the restrictions on the number and the maximal service distance of the service vehicles.After determining the lower bounded of l, sensitivity analysis of the number of service stations l is conducted after increasing the value of l from 9. Computational results with different values of l are shown in Figure 6.Total service costs experience a downward trend as l increases from 9 to 35, and the reduction becomes smaller with the increase in the value of l.For example, total service costs reduce by 12.31% when l changes from 9 to 10, whereas total service costs reduce by 3.37% when l increases from 34 to 35.Thus, in this case, the number of service stations may be set as 20 to 30 from a comprehensive perspective.tions, l; (b) the maximal service distance for service stations, D1max.
First, this paper conducts a sensitivity analysis of the number of service stations, l.Several examples are conducted to find the minimum value of l, and the computational results show that there is no optimal solution when the value of l is less than 9 due to the restrictions on the number and the maximal service distance of the service vehicles.After determining the lower bounded of l, sensitivity analysis of the number of service stations l is conducted after increasing the value of l from 9. Computational results with different values of l are shown in Figure 6.Total service costs experience a downward trend as l increases from 9 to 35, and the reduction becomes smaller with the increase in the value of l.For example, total service costs reduce by 12.31% when l changes from 9 to 10, whereas total service costs reduce by 3.37% when l increases from 34 to 35.Thus, in this case, the number of service stations may be set as 20 to 30 from a comprehensive perspective.This paper then changes the value of D1max to examine the effect of maximal service distance of service stations on total service costs.Multi-level facility location-allocation problems without the constraints of maximal service distance are first obtained, and the total service costs are 5.678 × CNY 10 5 with the maximal service distance of 237 km.After that, the maximal service distance for service stations is reduced from 240 km by 10 km each time, and optimal results are obtained using the proposed solution algorithm, which is shown in Figure 7.The optimal results show that the total service costs increase with the decrease in the value of Dmax, and there is no optimal solution until the maximal service distance Dmax is less than 120 km.Thus, manufacturers usually choose a suitable maximum service distance to ensure that they can obtain a small service cost while improving customer satisfaction.Therefore, the maximal service distance may be set as 130 km to 160 km from a comprehensive perspective.This paper then changes the value of D1 max to examine the effect of maximal service distance of service stations on total service costs.Multi-level facility location-allocation problems without the constraints of maximal service distance are first obtained, and the total service costs are 5.678 × CNY 10 5 with the maximal service distance of 237 km.After that, the maximal service distance for service stations is reduced from 240 km by 10 km each time, and optimal results are obtained using the proposed solution algorithm, which is shown in Figure 7.The optimal results show that the total service costs increase with the decrease in the value of D max , and there is no optimal solution until the maximal service distance D max is less than 120 km.Thus, manufacturers usually choose a suitable maximum service distance to ensure that they can obtain a small service cost while improving customer satisfaction.Therefore, the maximal service distance may be set as 130 km to 160 km from a comprehensive perspective.

Analysis of Contiguity Constraints
In the MIP-based optimization model, contiguity constraints ( 13)-( 16) are introduced to ensure all service units served by one service station are geographically contiguous.Thus, to analyse the effect of the contiguity constraints in MLFLAPs, optimal solutions

Analysis of Contiguity Constraints
In the MIP-based optimization model, contiguity constraints ( 13)-( 16) are introduced to ensure all service units served by one service station are geographically contiguous.Thus, to analyse the effect of the contiguity constraints in MLFLAPs, optimal solutions without constraints ( 13)-( 16) are obtained in this section.The total service costs without contiguity constraints are 5.760 × CNY 10 5 , a 1.30% reduction compared with the costs in Section 5.2.
The location of service stations and districting results without contiguity constraints are shown in Figure 8.The service districts served by service stations 1, 50, 12, and 122, are not contiguous, and service units 39, 40, 137, and 171 are served by the corresponding service stations across other service districts.In reality, the cross-regional maintenance service shown in Figure 8 is not generally a feasible solution in agriculture, which results in a long service latency and resource waste [5].In agriculture, the decision-maker may prefer to assign these service units to the nearest service stations, such as allocating service unit 39 to service station 50, service unit 40 to service station 43, service unit 137 to service station 138, and service unit 171 to service station 173.In all, contiguity constraints should be considered to prohibit cross-region services and ensure a contiguous service area for managers in practical application.In all, contiguity constraints should be considered to prohibit cross-region services and ensure a contiguous service area for managers in practical application.

The Performance of Benders Decomposition Approach
The Benders decomposition approach is analysed in this section.First, the cases in Figures 6 and 7 are solved by the Gurobi solver, which cannot be obtained in less than 12 h.Especially for some cases in Figures 6 and 7, the optimal solutions cannot be obtained in less than 24 h.But for the Benders decomposition approach, all optimal solutions can be obtained in 3 h.
In addition, the Gurobi and Benders decomposition solvers are applied to solve several cases with different sizes.The calculation time using different solution methods is shown in Table 2.As shown in Table 2, the calculation time of the Gurobi solver is shorter than that of the Benders decomposition algorithm only when the number of service units is less than 50.When the number of service units is greater than 50, the solution time using the commercial Gurobi solver will be much greater than that of the Benders decomposition algorithm.Even when the number of service units exceeds 200, the Benders decomposition algorithm can still obtain the optimal solution in 3 h, while the optimal solution cannot be obtained within 24 h using the commercial Gurobi solver.

Practical Application in Agriculture
In practical application, agricultural machinery always requires the construction of a multi-level maintenance service network including spare part centres and service stations to supply spare parts and provide maintenance service in a broad geographic agricultural operation area.For example, a famous agricultural machinery manufacturer in China has constructed several spare part centres in Harbin, Tianjin, Fuyang, and so on, and the manufacturer has built more than 300 service stations in China.
Thus, the proposed MIP-based optimization model and solution approach can be applied to solve the multi-level facilities location and service resource allocation problem, which can reduce service costs effectively and improve maintenance service.
Furthermore, agricultural maintenance service network software is developed to determine the location of service stations and spare part centres in China.As shown in Figure 9, the practical application in Henan Province is described based on the proposed solution approach.
Thus, the proposed MIP-based optimization model and solution approach can be applied to solve the multi-level facilities location and service resource allocation problem, which can reduce service costs effectively and improve maintenance service.
Furthermore, agricultural maintenance service network software is developed to determine the location of service stations and spare part centres in China.As shown in Figure 9, the practical application in Henan Province is described based on the proposed solution approach.

Conclusions and Future Work
This paper considered a multi-level facility location-allocation problem in agricultural maintenance service networks to find the optimal location of service stations and spare part centres while allocating service units to service stations.Considering multilevel constraints and contiguity criteria, a MILP-based optimal model was formulated to minimize total service costs.A Benders decomposition-based solution algorithm with several improvements was developed to address the MILP model.Finally, a real case in East China was introduced to demonstrate the performance of the developed mathematical

Conclusions and Future Work
This paper considered a multi-level facility location-allocation problem in agricultural maintenance service networks to find the optimal location of service stations and spare part centres while allocating service units to service stations.Considering multi-level constraints and contiguity criteria, a MILP-based optimal model was formulated to minimize total service costs.A Benders decomposition-based solution algorithm with several improvements was developed to address the MILP model.Finally, a real case in East China was introduced to demonstrate the performance of the developed mathematical model and solution approach, and sensitivity analysis was conducted to evaluate the effect of several parameters.
Computational results demonstrated that total service costs experience a downward trend with the increase in the number of service stations and the reduction becomes smaller with the increase in the number of service stations.In addition, total service costs also increase with the decrease in the value of the maximum service distance.In this case, the number of service stations may be set as 20 to 30, and the maximal service distance may be set as 130 km to 160 km from a comprehensive perspective.In addition, contiguity constraints should be considered to prohibit cross-region services and ensure a contiguous service area for managers in practical application.Lastly, considering the solution time, the proposed Benders decomposition can be applied in different size cases especially in the large size case.
There exist several limitations in this work.First, the maintenance demands are calculated based on the number of agricultural machines and agricultural production areas; thus, the maintenance demands in one service unit can be determined based on a more scientific method, combining historical data, the number of agricultural machines, and the machine learning method.In addition, agricultural machinery may move dynamically to other service units during the busy farming seasons, and the influence of the dynamic movement of agricultural machinery on maintenance demands is not considered in this paper.
Thus, this research can be improved in several ways in the future.On the one hand, considering the uncertain maintenance demand in agricultural service networks, robust optimization should be considered to construct a more stable maintenance service network.On the other hand, there are several service facilities participating in maintenance service in the busy farming season, and the location of these dynamic service facility locations should be considered in further studies.

Figure 1 .
Figure 1.A multi-level service network case.

KA
Set of service units indexed by i or i'; J Set of the potential location of service stations indexed by j, J I  ; Set of the potential location of spare part centres indexed by k, K I  ; Set of adjacent service unit pairs.(2) Parameters: n

Figure 1 .
Figure 1.A multi-level service network case.
Set of service units indexed by i or i ; J Set of the potential location of service stations indexed by j, J ⊆ I; K Set of the potential location of spare part centres indexed by k, K ⊆ I; A Set of adjacent service unit pairs.

Figure 5 .
Figure 5.The location-allocation results in multi-level service network.

Figure 5 .
Figure 5.The location-allocation results in multi-level service network.

Figure 6 .
Figure 6.Computational results with different values of l.

Figure 6 .
Figure 6.Computational results with different values of l.

Figure 7 .
Figure 7. Computational results with different values of D max .

Figure 9 .
Figure 9.The practical application in Henan Province.

Figure 9 .
Figure 9.The practical application in Henan Province.

Table 1 .
The location allocation results in the service network.

Table 2 .
The performance of the Bender decomposition approach.