Multimodal processes prototyping subject to grid-like network and fuzzy operation time constraints

This paper describes a Grid-like Material Transportation Network (GMTN) in which several heterogeneous means of transportation (Automated Guided Vehicles (AGVs), hoists, lifts, etc.) interact with each other via common shared workstations to provide a variety of demand-responsive material handling operations. Different material handling transport modes provide movement of workpieces between workstations along their manufacturing routes in the GMTN and they can be seen as processes realized with synergic utilization of various local periodically acting unimodal processes. The main contribution of this research is the solution of a constraint satisfaction problem addressing AGVs fleet match-up scheduling subject to GMTN and fuzzy operation time constraints. In the presented case both production rate (production takt) and operations execution time are described by imprecise (fuzzy) data. In other words, the research’s objective concerns assessing grid-like networks of periodically acting local transportation modes from the perspective of possible mass-customized oriented requirements imposed on scheduling of multimodal flows of jobs assigned to certain technological routes passing through common shared workstations.


Introduction
The productivity of an AGV-served flow shop repetitively producing a set of different products depends on both the job flow sequencing and the material handling system's ability to achieve a pre-specified throughput.Normally implemented modelling frameworks of this type assume a deterministic character of such systems, where demand is known in advance and the processing times of each job on each machine are known constants.In order to improve systems ' [e.g. Flexible Manufacturing Systems (FMS)] productivity the AGVs fleet size and the scheduling/dispatching influence on the system throughput are usually examined.More formally, such problems belong to a class of deterministic cyclic scheduling problems for which throughput rate targets are established and non-identical, multistage job routings and machine assignments are fixed.A goal is to sequence and position job operations in a cyclic fashion so as to minimize a weighted sum of the job flow times subject to transportation, blocking and time window constraints.The general problem is known to be NP-hard (Levner et al. 2010;Groenemeyer 2012).
In contrast to the above mentioned, the problem considered in this research assumes a given transportation system encompassing a network of AGVs periodically circulating along cyclic routes while servicing workpieces load/unload operations.The considered performance measures such as a production flow takt time, a flow time, a flow rate and a job shop productivity depend on uncertainty requirements caused by availability of workers employed for workstations servicing, e.g. for jig/tool fastening and pick and place operation performance.For the sake of further consideration the following definitions are used: production flow is defined as the movement of material through segregated stages of the workshop, flow rate (throughput) is the number of flow units (e.g.workpieces) going through the workshop (seen as the job-shop) per unit time, flow time (total processing time) is the amount of time a flow unit spends in the job shop from beginning to end, job shop productivity is the output of the workshop in a specific period of time, and takt time (production rate) is the rate at which a finished product needs to be completed.
In that context the investigated question concerns AGVs fleet match-up scheduling subject to productivity constraints with reference to the schedule of the production flow imposed by the uncertainty implied by workers' availability for operation execution.Consequently it can be seen as an extension of the problems formulated within the framework of a grid-like multimodal transportation network in which an arbitrarily chosen Elementary Transport Networks (ETNs) supported by different transportation modes interact with each other via distinguished subsets of common shared workstations as to provide a variety of transportation/handling services (Bocewicz et al. 2014a, b), and then in (Bocewicz et al. 2015).Moreover, such commonly found regular structures, e.g.grid-like topology networks (Qiu et al. 2002;Balakrishnan and Chun 2007;Perera and Ratnayake 2015), are widely observed in different application domains (such as manufacturing, intercity freight transportation supply chains, multimodal passenger transport), provide additional motivation for this study.
The realistic examples supporting the adopted topological assumptions follow from numerous reports concerning mesh-like or grid-like as well as fractal-like structures of GMTNs.Such topologies have been considered for over twenty years see Qiu et al. (2002).At the same time fractal-like topologies are investigated as well see Silva (2010).Drawing upon the research listed above, our topological assumptions are motivated by the expected advantages gained from the GMTN layout flexibility being vital to improve manufacturing flexibility.This is especially the case when faced with shorter lead times, variable product mix and volumes.
Since, most manufacturing processes observed in their steady state are periodically recurring ones, cyclic scheduling methods can be considered.This kind of representation enables one to show which behaviours (consisting of many cyclic concurrent threads corresponding with routes of AGVs, routes of transported work pieces, etc.) can exist in any given structure (limited by the number of resources and their capacities).Many models and methods aimed at cyclic scheduling of concurrently flowing jobs have been considered to date.Among those, the mathematical programming approach (Von Kampmeyer 2006), max-plus algebra (Polak et al. 2004), Petri nets (Baruwa and Piera 2016), tabu search (Zheng et al. 2014) constraint logic programming (Bocewicz et al. 2014a, b) frameworks are the more frequently used.Most of these methods aim at finding a minimal cycle or maximal throughput while assuming deadlock-free processes flow.Very rarely, however, they concern production systems consisting of an AGV fleet arranged to transport workpieces between workstations (i.e.following a given production route) while moving in a carousel mode along cyclic path topologies.In that context the main question behind the study is: Can the AGV fleet assignment reach its goal subject to the constraints assumed on the concurrent multi-product manufacturing system?
The remainder of the paper is organized as follows: Sect. 2 provides a brief review of related works.Section 3 introduces a framework aimed at GMTN modeling in terms of System of Concurrently executed Cyclic Processes (SCCP) and provides the problem formulation concerning AGV fleet scheduling with a multi-product manufacturing flow schedule imposed by fuzzy operation time constraints.After discussion of some issues regarding multimodal processes prototyping, the solution to the considered constraint satisfaction problem is submitted in Sect. 4.An illustrative example providing results of computational experiments is given in Sect. 5. Concluding remarks and further work are presented in Sect.6.

Literature review
Numerous papers have been published to address Vehicle Routing Problems (VRP) (Carić and Gold 2008), AGVs fleet assignment (Abara 1989;Hall et al. 2001), integrated machine scheduling and vehicle routing with time windows (Ullrich 2013) as well as multimodal network flow models (Haghani and Oh 1996).Because of the NP-hard nature of these problems, the proposed solution methods range from dynamic programming (Maneengam et al. 2012;Azi et al. 2012), through constraint programming (Kao and Chen 2013;Kilby and Urli 2016) to metaheuristic-driven evolutionary programming (Ullrich 2013;Reed et al. 2014).
Since VRPs assume searching for shortest job flows time subject to transportation, blocking and time window constraints it can be also seen as an extension of the Heterogeneous Fleet Vehicle Routing Problem with Time Windows (HFVRPTW), where workpieces transportation following production routes must occur in certain time windows, i.e. the interval is treated as part of the respective production takts, which varies from production order to production order.The goal of the HFVRPTW is to find transportation times for each AGV stream such that hard time window constraints on both kind of windows hold for a cyclic steady state of AGV flow.
Other perspectives provided by the integration of fleet assignment and maintenance planning (Moudani and Mora-Camino 2000;Lu et al. 2013) to AGV fleet scheduling is presented in Bzdyra et al. (2015) where a declarative framework enabling one to determine conditions employed in a decision support system aimed at small and medium size enterprises involved in a unique, multi project-like and mass customized oriented production is discussed.The unique production orders are grouped into a set of portfolio orders where each production order is treated as an activity network of common shared resources, known in advance, however imprecise operation times are allotted.The job routings may be arbitrary partial orders (i.e., they are not restricted to be serial), and the jobs may be re-entrant, meaning that the same machine may be visited multiple times before completion.The problem concerns scheduling a newly inserted project portfolio taking into account both: imprecise operations imposed by a multi-project environment and material handling constraints imposed by AGV availability.A declarative model is sought, enabling one to state a constraint satisfaction problem aimed at: Whether a given portfolio can be completed within the assumed time period in the given manufacturing system?
The material handling system of a flexible manufacturing system where vehicles serve workstations located along different modes of cyclic guided paths in a 3D environment composed of rails/lifts/hoists is considered in (Bocewicz et al. 2016).A production flow scheduling problem subject to fuzzy processing time constraints treated as a fuzzy constraint satisfaction problem is studied.The main objective is to provide sufficient conditions guaranteeing avoidance of time-consuming computer-simulation-based calculations of AGV fleet schedules, while keeping the throughput at the maximum achievable level.Nielsen et al. (2016) provides some generalization of the problem, stating the resourceconstrained multi-product scheduling problem for an AGV-served flow shop, where multiple material handling transport modes provide movement of workpieces between machining centers.
Consequently, the problem considered in this research can be seen as an extension and/or refining of the problems formulated in our former works (Bocewicz and Banaszak 2013;Bocewicz et al. 2014aBocewicz et al. , b, 2015Bocewicz et al. , 2016;;Nielsen et al. 2016) focused on a solution to a constraint satisfaction problem aimed at AGVs fleet match-up scheduling while taking into account fuzzy operation times constraints.In order to delimit the work addressed in this research we will first present contributions concerning some of the main issues.The schedulability analysis of the processes executed in the so-called systems of concurrent cyclic processes executed within a declarative modeling framework is provided in Bocewicz and Banaszak (2013).The considered representation provides a unified way for evaluating the performance of local cyclic processes as well as of the multimodal processes supported by these.We propose a recursive design approach to create sufficient conditions that guarantee the schedulability of both local and multimodal cyclic processes.In turn the solution to a constraint satisfaction problem aimed at AGV fleet match-up scheduling while taking into consideration the itineraries of concurrently manufactured product types is provided in Bocewicz et al. (2014b).An algebra-like driven approach to cyclic scheduling based on step-by-step composition of multimodal transportation network sub-structures is proposed.
The considered network of repetitively acting local transportation modes provides a framework for multimodal processes scheduling treated in terms of AGV fleet scheduling problems subject to fuzzy operation time constraints.The aim of the research is to show the universality of the employed declarative modeling driven approach.This approach makes it possible to treat the above mentioned scheduling problem in terms of the following two routine queries formulated either in: a direct way-What schedule of AGV fleet following a given set of operation times maximizes operation availability and productivity?, or a reverse way-What set of operation times guarantee that the resultant schedule of a given AGV fleet maximizes both operation availability and productivity?
3 Grid-like transportation network modelling Note, that a specific network topology of GMTN, e.g. as shown in Fig. 1a), assumes that it can be considered as composed of multiple Elementary Transport Networks (ETNs).Each ETN can be seen as an autonomous cell (or grid) encompassing cyclically circulating AGVs transporting workpieces among workstations.Workpieces moved along different ETNs while following different production routes can be viewed as multimodal processes passing along locally used transportation paths and workpiece processing operations.
The assumption behind these statements supposes a regular composition of ETNs, encompassing local cyclic steady states, implies a cyclic steady state of whole GMTN, i.e. guaranteeing congestion and deadlock-free flow of AGVs.In a mixed floor and overhead material handling transport system equipped with many unidirectional AGV and overhead hoist transport modes guaranteeing congestion and deadlocks-free operations plays a pivotal role in achieving a well-organized material flow.Our approach focuses on GMTN steady state behaviour following assumed steady state task requirements, enables deadlock-free scheduling and routing of a great number of multimodal transport means.

Systems of concurrent cyclic processes modelled in GTN framework
Consider a set of workstations deployed along paths in the transportation network, where fleet of AGVs service workpiece movements between workstations.Due to the organization of material handling operations AGVs circulate along transportation paths, loops in the network composed of multiple AGV loops connected through common shared workstations.In order to guarantee each workstation can be solely serviced by exactly one AGV at the moment, a dispatching rule choosing a transportation mean among currently accessing ones has to be assigned.In that contexts the considered transportation system can be seen as composed of a set of concurrently executed cyclic process synchronized by mutual exclusion protocol, i.e. a System of Concurrently executed Cyclic Processes (SCCP).In turn, SCCP can be directly represented in terms of Petri nets formalism being a powerful formal model for the analysis and modelling of concurrent systems, hence a rich theory has developed around them so far can be also easily used.A way in which the GMTN from Fig. 1a can be modelled as a SCCP is shown in Fig. 1b (Bocewicz et al. 2015).In other words, each SCCP consists of a set of periodically performed processes composed of material handling and processing operations.Processes interact with each other via common used resources in the form of either transportation or processing operations.Dispatching rules allotted to common shared resources provide a mutual exclusion protocol synchronizing processes interaction (Bocewicz and Banaszak 2013).
The considered GMTN consists of ETNs with the structure presented in Fig. 2a.In terms of SCCP representation, it will be called a SCCP-driven model of ETN, (SETN) for short.As the considered GMTN consists of a set of ETNs, therefore to each a-th ETN the relevant a-th SETN denoted as (a) SC belongs as well, see Fig. 2b.
The SETN shown in Fig. 2b, is composed of four local cyclic processes : P 1 , P 2 , P 3 , P 4 and their streams: P 1 1 , P 1 2 , P 1 3 , P 1 4 (distinguished in Fig. 2 by (a) P 1 , (a) P 2 , (a) P 3 , (a) P 4 , respectively-superscript (a) indicates the number of SETNs), associated to four AGVs following their routes in the ETN.Apart from local processes, two multimodal processes: m P 1 , m P 2 representing parts of manufacturing routes of two products W 1 and W 2 executed in the ETN, are considered.The processes follow given routes composed of workstations: R 1 − R 4 (resources (a) R c denoted by superscript (a) that refers to the number of the ETN) and resources representing transportation sectors (R 5 − R 13 ).Formally, local P and multimodal m P processes are defined in the following way: P = {P i |i = 1, . . ., n} − the set of local processes described by sequences of operations executed on resources R, where each P i is specified by the set of streams: , m P = {m P i |i = 1, . . ., w} − the set of multimodal processes described by sequences of some sub-sequences from local processes P, where each m P i is specified by the set of streams: (processes (a) m P i in Fig. 2 denoted by a superscript (a) which indicates the number of ETNs).
The operations of multimodal processes' streams require the execution of some local processes.For example, transport operations between resources in m P 1 (product W 1 , see Fig. 2) require the two streams of local processes: P 1 and P 4 .This means the routes of multimodal processes are also determined by the subsequences of routes of the local processes through which they have to be processed.It is assumed that processes are synchronized by a mutual exclusion protocol guaranteeing that only one process can be executed on a shared resource at a given time.
In order to describe the system shown in Fig. 1b and SETN in Fig. 2) let us introduce the following notations (Bocewicz et al. 2014a, b): lr(i) specifies the route of the stream of a local process P k i (the k-th stream of the i-th local process P i ), and its components define the resources used in course of process operations execution, where: p k i, j ∈ R, and R = {R 1 , R 2 , . . ., R c , . . ., R m }-denotes the set of resources used for the j-th operation in the k-th stream of the i-th local process; in the remainder of the paper, the j -th operation executed on the resource p k i, j in the stream P k i will be denoted by o k i, j ; lr (i)-is the length of the cyclic process route.For example, in the SETN seen in Fig. 2b) the routes of streams of processes P 1 , P 2 , P 3 and P 4 (e.g.executed by relevant AGVs), are defined by the following sequences: denotes the time of execution of the j-th operation o k i, j in the stream P k i .In the case of the SETN from Fig. 2b) t k i, j means either the transportation time assigned to the sector (i.e.resources R 5 − R 15 ) or the load/unload operation time assigned to the workstation (i.e.resources R 1 , R 2 , R 3 , R 4 , R 16 , R 17 ).).
i y β i y , γ i y − specifies the route of the stream m P k i from the multimodal process m P i (the k-th stream of the multimodal process m P i ), where: In the considered case of the SETN from Fig. 2b), streams of multimodal processes m P 1 , m P 2 are specified by the following sequences: ), being the concatenation of subsequences p 1 1 , p 1 4 and of p 1 3 , p 1 2 , respectively.In the remainder of the paper, the j -th operation of the stream m P k i will be denoted by mo k i, j .The operation times mt k i, j of multimodal processes are: the operation times of AGV moving between the workstations and operation times assigned to workstations R 1 − R 4 , R 16 , R 17 (in the model the setup times are included in the operation times).The operation time mt k i, j assigned to workstation R g is equal to: , where: x δ β,γ (l)-is the moment the process m P k i leaves R g (taken by local the stream P δ β ) and x ψ ϕ,φ (l)-is the moment the process m P k i enters R g (delivered by the stream is the set of the priority dispatching rules of local (v = 0) and multimodal processes (v = 1), where σ v c is the sequence components which determine the order in which streams of processes can be executed on resource R c .Dispatching rules that determine the order in which vehicles access (streams of local processes) common shared resources R 2 , R 8 , R 4 are the following: σ 0 2 = P 1 1 , P 1 2 , σ 0 8 = P 1 1 , P 1 2 , σ 0 4 = P 1 1 , P 1 4 , P 1 3 , P 1 2 , in turn, dispatching rules synchronizing multimodal processes execution are as follows: Therefore the structure of SCCP can be defined as the pair: where In simple terms the structure of the SCCP can be treated as a set of pipeline cyclic processes (P, m P) interacting via common shared resources R. Potential process conflicts are resolved in advance by the priority dispatching rules (Θ 0 , Θ 1 ).In that context, the pipeline-like processes typical for a straight-line production allows one to take into account both one-piece and variable batch/lot size flows.
In the case of GMTN (e.g.Fig. 1b) the structure SC (1) can be seen as a result of multiple compositions of ETNs (Fig. 2b).Consequently, the structure SC (1) consists of the set of SETNs SC * such that: (a) each substructure (a) SC ∈ SC * of the structure SC is defined similarly to (1): where: Rp a ⊂ R-the set of resources of the substructure (a) SC, SLp a = (Pp a , U p a , T p a , Θ 0 a -the structure of local processes of the substructure (a) SC including: the set of local processes Pp a ⊂ P, the set of processes routes U p a ⊂ U ; the set of execution times T p a ⊂ T of operations belonging to processes Pp a ; the set of dispatching rules Θ 0 a assigned to resources belonging to Rp a ; SMp a -the structure of multimodal processes of substructure (a) SC, including fragments of multimodal processes m P j (β, γ ) forming the set m Pp a (Bocewicz et al. 2015).(b) two sub-structures (a) SC, (b) SC ∈ SC * are called isomorphic if: Note that each ETN modelled in terms of a SCCP can be seen as an isomorphic substructure of SETN, i.e. (a) SC (see Fig. 2b) consisting of local processes (a) P 1 , (a) P 2 , (a) P 3 , (a) P 4 and parts belonging to multimodal processes (a) m P 1 , (a) m P 2 .

Uncertainty requirements
Besides the parameters describing the SC substructure (a) SC( 2), such as sets of routes, processes, dispatching rules, workers availability etc., variables reflecting its performance (i.e.job shop productivity, flow time and takt time) are also considered.In order to specify the constraints assumed on operation times the following variables (a) s i, j , (a) sT i have to be taken into account: • Operation time (a) mt i, j (the time the operation (a) mo k i, j executed in (a) SC) allocated to a workstation (resource) (a) R c from (a) SC is associated with the so called operation availability specified by the fuzzy set (a) s i, j = μ mt i, j , (a) s i, j ∈ [0, 1].The membership function μ (a) mt i, j determines workers' availability for operations handling on the resource (a) R c . (a) s i, j = 0 indicates that the corresponding value (a) mt i, j is out of worker's availability, while (a) s i, j = 1-indicates availability, with acceptable operation time value (a) mt i, j .The accepted production takt (a) mT c i of the i-th multimodal process has to follow the assumed Multimodal Process driven Awaited Productivity (a) sT i (MPAP (a) sT i -for short) measured along m P i in the a-th SETN, i.e. in the (a) SC.
Table 1 Resource operation availabilities for multimodal processes (a) m P 1 , (a) m P 2 of SETN from Fig. 2b Table 2 Awaited productivity measured along multimodal processes (a) m P 1 , (a) m P 2 in a production takt • The takt time of the process m P i , is the time between two consecutive moments of its subsequent streams completion: (a) mT c i = mx k i, j (l) − mx k−1 i, j (l).Similar to the operation availability (a) s i, j , the productivity (a) sT i is specified by the fuzzy set (a) sT i = μ (a) mT c i , (a) sT i ∈ [0, 1].The membership function μ (a) mT c i determines productivity following (a) mT c i in (a) SC. (a) sT i = 0 indicates that production performed with takt (a) mT c i is not profitable, while (a) sT i = 1 indicates that takt (a) mT c i ensures productivity.
The introduced variables can be illustrated on the following example limited to the SETN from Fig. 2b.Consider operation availability and productivity specified by the membership functions: μ (a) mt i, j , μ (a) mT c i as shown in Tables 1 and 2. Following from these, uncertainty requirements as well as constraints guaranteeing deadlock-free processes execution contribute to a declarative model of SCCP behavior is considered in the next section.

Problem formulation
The majority of problems stated within a SCCP framework can be seen as concerning the relationship between its structure SC and behavior (e.g. the schedule).In that context the main question concerns the SCCP's cyclic behavior, following requirements imposed by operations availability and awaited productivity subject to a given Θ? Since these kind of problems are NP-hard in nature, i.e. computationally intractable, hence the optimal solutions can be obtained only for small scale instances that are quite far from real-life cases (Bocewicz et al. 2014a, b;Levner et al. 2010).However, this disadvantage can be removed in case of a SCCP (modelling a GMTN) which can be seen as composed of a set of SETNs.In order to show this let us note that the considered problem can be decomposed into the following two sub-problems.The first one can be reduced to finding the SETN and corresponding set of dispatching rules Θ allocation of which guarantees cyclic behavior of the SCCP (consisting of set of isomorphic SETNs).The second sub-problem comes down to the question of whether the assumed uncertainty requirements regarding workers' availability and job shop productivity can be satisfied for the operation times T , mT specifying the considered SCCP?
The obtained pair: (a set of dispatching rules Θ, a set of operation times T , mT ) and determined by these, the system's behaviour can be seen as a distributed control procedure provided to the spatially distributed common shared resources while imposing the right initial process allocations.These kinds of control procedures can then be directly implemented in SCADA (Supervisory Control And Data Acquisition) systems that tie together decentralized facilities such as shared by technological processes AGVs and machine tools (under the given uncertainty requirements).
Above problem can be reduced to the question of whether it is possible to extract the SETN such that there exists a set of dispatching rules which guarantees cyclic behavior of the considered SCCP (representing the GMTN) while following fuzzy constraints imposed by awaited productivity and assumed operations' availability?Its graphical illustration is shown in Fig. 3.The research methodology used for this study consists of the following five stages (see Fig. 4).
Given SCCP (representing GMTN) and fuzzy constraints imposed on awaited productivity and operations availability SCCP Does there exist SETN and a set of dispatching rules allocation of which results in cyclic schedule encompassing cyclic behavior of considered SCCP and following assumed fuzzy constraints ?
Cyclic schedule following fuzzy constraints imposed by awaited productivity and assumed operations availability The main idea behind the proposed approach follows from the observation that the assessment of possible cyclic behaviors of GMTN (SCCP) can be replaced the by assessment of its ETN (SETN).Therefore, in order to consider possible deadlock-free/collision-free executions of processes executed in GMTN, let us focus on constraints of a Constraint Satisfaction Problem dedicated to its ETN.

Constraints
The SC (1) model emphasizes the structural characteristics of the modeled SCCP.Due to deadlock occurrence caused by the condition of mutual exclusion, not all behaviors of SCCP are allowed in these kinds of structures.The behavioral characteristics guaranteeing deadlockfree execution of processes can be specified in terms of the admissible states reachability concept, i.e. either driven by both cyclic steady states and leading to transient states (Bocewicz et al. 2012) or by the cyclic steady state space only (Bocewicz and Banaszak 2013).The second approach, which does not take into account the initial states leading through the transient states to the cyclic steady states, seems to be quite close to the cyclic scheduling methods widely used in many real-life cases (Von Kampmeyer 2006).In that context, the behavioral characteristics guaranteeing deadlock-free execution of processes while following the above mentioned uncertainty requirements can be specified by the following cyclic schedule: This kind of schedule X SC is defined as a sequence of ordered pairs describing the behavior of local (X, α) and multimodal (m X, mα) processes, where: -is the set of moments x k i, j (for l = 0-th cycle) of operations o k i, j beginning from the stream P k i and α-is the period of the cyclic schedule X .By analogy the set m X = mx 1 1,1 , . . ., mx k i, j . . .., mx lms(w) w,lm(w) consists of moments of multimodal processes operations beginning, and mα-is the period of the multimodal cyclic schedule m X.
Since each possible SCCP's behaviour can be specified by a schedule as a sequence of vectors' coordinates of which define resources currently allotted to operations over discrete time, hence any obtainable cyclic schedule X SC (3) can be seen as the representation of a cyclic steady state reachable in the state space generated by the SCCP structure (Bocewicz et al. 2012(Bocewicz et al. , 2015)).
Variables x k i, j /mx k i, j ∈ Z determine moments of operations beginning in the l-th cycle of the SCCP cyclic steady state:

Since values of x k
i, j /mx k i, j follow the system structure parameters, the cyclic behavior X SC is determined by SC (1).The constraints determining the admissible cyclic schedule X SC (3) are the following: Uncertainty requirements (Tables 1 and 2): The level operations' availability S and level of the productivity E (caused by X SC ) for the whole SCCP, are determined as a minimal value (a) s i, j among all the operations availabilities and a minimal MPAP productivity (a) sT i among all the multimodal processes, respectively: S = min i=1,2; j=1,...,4 (a) s i, j and E = min i=1,2; j=1,...,4 {(a)sT i } (4) Constraints describing the local processes execution: the moment of operation o k i, j beginning states for a maximum of both: the completion time of the operation o k i, j−1 preceding o k i, j , and the release time of the resource p k i, j awaiting for o k i, j execution: Moment of operationo k i, j begining = max moment of p k i, j release + lag time t , Moment of operation o k i, j−1 completion (5) Constraints describing multimodal processes execution: the moment of operation mo k i, j beginning is equal to the nearest admissible value (determined by set X k i, j of values mx k i, j ) being a maximum of both: the completion time of operation mo k i, j−1 preceding mo i, j , and the release time (lag time tm) of the resource mp k i, j awaiting mo i, j execution.
Moment of operation mo k i, j begining = max moment of mp k i, j release + lag time tm , moment of operation mo k i, j−1 completion where i, j determined by value of x δ β,γ , where: x δ β,γ -the moment of operation o β,γ beginning enabling execution of the operation mo k j , β B = min {ϕ ∈ B : ϕ ≥ β}.Constraints (4), ( 5) and ( 6) describe conditions guaranteeing execution the GMTN's processes subject to arbitrary given uncertainty requirements.A number of admissible schedules X SC (3), denoted by L, depends SC (1), and especially on the assumed operations times T of local processes.In that context different values of operation times t k i, j (i.e., representing transportation and layover times) lead to schedules with different levels of operations' availability S and productivity E (according to Tables 1 and 2).
An illustration of this fact can be seen in Fig. 5, where the number L of admissible schedules X SC is shown following from ( 5) and ( 6) for the SETN from Fig. 2b), while taking into account uncertainty requirements S and E (4).The results observed are obtained for operation times t k i, j ∈ {1 . . .6}. Space L from Fig. 5 consists of 104,946 (among 4 × 10 5 ) cyclic schedules X SC following constraints S > 0 and E > 0.

CSP-driven multimodal processes scheduling
Figure 6 presents details specifying each substructure, (a) SC, following the SCCP from Fig. 1b.Processes occurring in each substructure, e.g. (a) SC, (b) SC, (c) SC, …, ( f ) SC see Fig. 6b, are implemented in the same manner, this means the operations are performed along similar routes, the same dispatching rules are applied, etc.
In this context the introduced operator of substructures (a) SC composition can be used for the whole SC reconstruction: where SC-indicates a composition following ( 5) and (6), i.e. each substructure (a) SC is linked with the others by means of integrating the resources belonging to the same set of corresponding resources.For example, the structure (a) SC from Fig. 6b) is linked with (e) SC by the resource (a) R 16 , with (b) SC by the resource (a) R 17 and with (c) SC by the resources (a) R 14 and (a) R 15 .
Due to the same manner of process execution as well as the same manner of substructures composition, the cyclic schedule X SC (3) representing the behavior of the whole structure can be perceived as a composition of corresponding cyclic schedules: where: (a) X − the cyclic schedule of the substructure (a) SC: (a) X = (a) X, (a) α , (a) m X, (a) mα (9) (a) X/ (a) m X− the set of the moments of operations beginning of local/multimodal process operations occurring in the substructure (a) SC; (a) α / (a) mα-the period of local/multimodal processes executions; a) α, (b) α , (a) m X ∪ (b) m X, c • lcm (a) α, (b) α (10) where: c ∈ N + , lcm (a, b) -the least common multiple of a and b.In the considered system (i.e.composed of one kind of SETN-see Fig. 6) it is enough to know the schedule (a) X of substructure (a) SC to determine the schedule X .However, to make the composition (8) possible, it is necessary to ensure that the operations executed according to (a) X , does not lead to deadlocks.Moreover, the schedule (a) X is calculated for the "convoluted" structure (a) SC (see example shown in Fig. 6c).The convoluted structure (a) SC is obtained from (a) SC as a result of merging corresponding resources distinguished in neighbouring structures.In the example shown in Fig. 6c the structure (a) SC has been obtained by way of merging respectively: (a) R 1 and (a) R 17 ; (a) R 3 and (a) R 16 ; (a) R 5 and (a) R 14 ; (a) R 6 and (a) R 15 .It should be noted, that local processes executed in the considered resource-consolidated structures (a) SC are both collision-free and deadlock-free.This means the substructure (a) SC from which the structure (a) SC was born also interact with its neighboring structures (e.g. (a) SC, (b) SC, (c) SC, (e) SC) in a deadlock-free way.
Consider the structure SC (2) modeling a given GMTN, composed of ETNs specified by "convoluted" structures (a) SC.The periodicity of SETN specified by cyclic, i.e. specified by cyclic schedule X SC , behavior of the whole SC.
In that context the question corresponding with considered problem is reduced to the following: Which dispatching rules (a) Θ and operation times (a) T , (a) mT , guarantee that in the relevant structure (a) SC there exists the cyclic schedule (a) X such that following conditions hold: (a) S ≥ S T H and (a) E ≥ E T H ?
Consequently, the positive answer of this question allows one to determine Θ T and mT (based on (a) Θ, (a) T , (a) mT ) guaranteeing required values of S and E.
In order to determine (a) Θ, (a) T , (a) mT guaranteeing the attainability of the cyclic schedule (a) X within the structure of (a) SC, the following constraint satisfaction problem has to be considered (11): (a) X , (a) Θ, (a) T , (a) α , where (a) X , (a) Θ, (a) T , (a) α -decision variables, (a) X -the cyclic schedule (8) of substructure (a) SC, (a) Θdispatching rules of (a) SC, (a) T = (a) T, (a) mT -the sequence of local and multimodal operation times, (a) α = (a) α, (a) mα -the period of local and multimodal processes for (a) SC, D X , D , D T , D α − domains determining admissible values of decision variables: D X : mx k i, j , x k i, j ∈ Z; D T : mt k i, j , t k i, j ∈ {1 . . .10}, C-the set of constraints describing SCCP behavior, the set C includes the constraints describing the execution of local and multimodal processes (4), ( 5), ( 6) as well as uncertainty requirements constraints S ≥ S T H , E ≥ E T H .
The schedule (a) X that meets all the constraints from the given set C is the solution sought for the problem (11).It can be obtained using standard constraint programming driven software platforms such as ILOG, OzMozart, ECL i PSE (Sitek and Wikarek 2015).

Computational experiments
Consider GMTN from Fig. 1a where four different products W 1 −W 4 , are manufactured along distinguished routes followed by multimodal processes m P 1 -m P 4 .The GMTN modeled in terms of SCCP (see Fig. 1b) reflects itself in the structure SC (1).
In that context, taking into account the first stage of the search methodology sketched in Fig. 4, let us consider the SETN and its structure (a) SC as shown in Fig. 6b.The SETN consists of four local cyclic processes (a) P 1 , (a) P 2 , (a) P 3 , (a) P 4 corresponding to the operations of four AGVs and two multimodal processes (a) m P 1 , (a) m P 2 encompassing appropriate parts of m P 1 -m P 4 .In turn, due to the second stage of the search methodology sketched in Fig. 4, the dispatching rules allotted to shared resources in the convoluted representation of (a) SC, i.e. (a) SC as shown in Fig. 6c, are sought.Consequently, the third stage of the search methodology sketched in Fig. 4, the response to the question: Does there exist the dispatching rules (a) Θ guaranteeing the cyclic schedule (a) X ? is sought.Assuming all operation times t k i, j occurring in (a) T are unit, i.e., t k i, j = 1 the PS REX i (11) has been implemented and then computed in the constraint programming environment OzMozart system (version 1.3.2,Dual Core 2.67, GHz, 2.0, GB RAM).The first solution obtained in less than 1 s is shown in the Table 3.
The obtained dispatching rules guarantee cyclic steady state execution of processes completing the convoluted representation (a) SC.  a) σ 1/17 = (a) P 1 4 , (a) P 1 1 (a) R 4 (a) σ 4 = (a) P 1 1 , (a) P 1 4 , (a) P 1 3 , (a) P 1 2 a) σ 8 = (a) P 1 2 , (a) P 1 1 a) σ 3/16 = (a) P 1 2 , (a) P 1 3 , (a) P 1 1 It should be noted, that obtained solution for operation times, assumed to be unique, is also valid for other values from N. Since the state evolution observed in SCCP depends entirely on the occurrence of discrete events over time, the SCCP as well as (a) SC can be treated in terms of discrete event systems (DES) (Cuninghame-Green and Butkovič 1995;Floudas and Lin 2005).Consequently, assuming an expected productivity measured by production takts (a) mT c 1 and (a) mT c 2 of multimodal processes (a) m P 1 , (a) m P 2 (see Table 2) as well as expected operation times, (see Table 1) the fourth stage of search methodology sketched in Fig. 4, seeks the response to the question below: Does there exist a plan of production flow (specified by (a) T , X SC ), guaranteeing the following requirements hold: S ≥ 0.7 (level of operations' availability is not less than S T H = 0.7), E ≥ 0.7, (level of job shop productivity, i.e. minimal level value among the all multimodal processes are not less than E T H =0.7)?.
The response to this question depends on the parameters (a) X and (a) T specifying (a) SC (from Fig. 6b), i.e., on the solution to the problem PS REX i (11).This problem has been implemented and then computed in the constraint programming environment OzMozart system (version 1.3.2,Dual Core 2.67, GHz, 2.0, GB RAM).According to Fig. 5 the search process took into account the space L of admissible solutions (schedules) containing over 4 × 10 5 instances.The final solution was obtained in less than one second.The space L does not contain any admissible production flow plan following the assumed degrees S and E, see Fig. 7.
Therefore, the following request with lower values of S T H and E T H has been considered: Does there exist a plan of production flow (specified by (a) T , X SC ) guaranteeing the following requirements hold S ≥ 0.6, E ≥ 0.6?
The final solution was obtained in less than one second as well.The set of admissible solutions L is shown in Fig. 8.The sample of schedule (a) X , obtained for the convoluted structure (a) SC (from Fig. 6c), and following (a) 2) while leading to (0.66; 0.6), is shown in Fig. 9.
In turn, the obtained schedule (a) X enables one to determine the final schedule X SC (8)due to the fifth stage of the methodology shown in Fig. 4. The Gantt chart shown in Fig. 11 provides a graphical illustration of the process m P 2 (distinguished in Fig. 10) execution guaranteeing workers availability S on the level 0.66 and productivity E on the level 0.6 (production takt 24 u.t.).Summing up the space L contains an admissible production flow plan following assumed degrees S and E, see Fig. 8.The obtained schedule X SC , (following S = 0.66 and E = 0.6) is shown in Fig. 11.
The above discussed issues can be supplemented by evaluating the performance indexes, such as: flow rate, flow time, takt time and job shop productivity.The example results collected   4 The values of the flow rate, the flow time, the takt time and the job shop productivity E achieved by (a) SC from Fig. 6b and SCCP from Fig. 1b for the schedule depicted by the Gantt chart in Fig. 11 (a) SC from Fig. 6b SCCP from Fig. 1b Flow rate (streams/unit time) 0,083(3) 0,16(6) job shop productivity E 0,6 0,6 a) mP 1 (a)   10 in Table 4 are obtained for (a) SC from Fig. 6b, and SCCP from Fig. 1b, and the schedule X SC following S = 0.66, E = 0.6 and depicted by the Gantt chart in Fig. 11.The flow rate equal to 0,083(3) is calculated due to the observation that 2 workpieces can be completed per 24 u.t.within the ETN modeled by (a) SC.This value is twice as small as the 0,16(6) observed in the SCCP from Fig. 1b where 4 workpieces can be completed in the 24 u.t.period.In turn, the flow time calculated due to the following formulae F T i = lz • (a) mT c i − T i depends on a number lz of SETNs passed by multimodal processes and an amount of u.t.determining the faster completion of process performance T i < (a) mT c i .Therefore, in case of SCCP from Fig. 1b the flow times for m P 1 and m P 2 are equal to F T 1 = F T 2 = 3•24−2 = 70 and in case of m P 3 and m P 4 are equal to F T 3 = F T 4 = 3•24−1 = 71.
To summarize let us point out that computational experiments have been conducted for SCCP composed of 8 SETNs, i.e. (a) SC, (b) SC, …, (h) SC-see Fig. 10, or in other words consisting 32 local and 4 multimodal processes as well as 118 resources.Due to proposed research methodology the set of priority dispatching rules Θ guaranteeing reachability of SCCP cyclic steady state have been obtained in less than 1 s.In order to evaluate benefits following from this GMTN decomposition driven approach many additional experiments have been conducted.Some results are collected in Table 5.It should be noted that in the case of SCCP from Fig. 10 (number of SETN = 8) the first solution, i.e.Θ, providing cyclic steady state has been obtained after 5400 s, while in the case of SETN from Fig. 6b, the first solution, i.e. (a) Θ, has been obtained in time less than 1 s.

Conclusions
The proposed declarative approach, aimed at AGV fleet scheduling stated in terms of constraint satisfaction problem representation, provides a unified framework for performance evaluation of local as well as multimodal processes supported by these.The considered problem is reduced to searching for a multimodal processes schedule that guarantee the requirements of workers' availability S and productivity E hold while taking into account a set of priority dispatching rules synchronizing local processes execution.
The proposed approach leads to solutions based on sufficient conditions that allow the network designer to compose ETNs in such a way as to obtain the final AGV scheduling system with required quantitative and qualitative behavior features.This provides a step towards a method to replace the exhaustive search for the admissible control by an iterative structural feature transportation system design guaranteeing its required behavior.
In our model setup times are included in the processing times, i.e. similar as in most of classical machine scheduling research where setup times are considered either as negligible or as parts of the processing times of jobs.Of course, since interest about scheduling problems with setup times has grown significantly during the last four decades, see: (Yang and Liao 1999), this issue will be taken into account in our further work as well.
Future research will be aimed at grid communication networks, while devoted to scalability extensions of the proposed declarative-model-driven approach.Since the regularity of grid transportation networks provide a modelling framework enabling significant simplification of the evaluation of the overall network performance the ordered fuzzy numbers concept (Kosiński et al. 2013) and hybrid techniques (Sitek and Wikarek 2015) joining MLP and CP may dramatically decrease computation time requirements.Coupling both the grid-like structure based methodology of transportation networks and a new problem modelling concepts may result in a new generation of DSS providing online solutions to real-life-sized communication transportation problems.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/),which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Fig. 1
Fig. 1 An example of GMTN (a) and its SCCP model (b)

Fig. 2
Fig. 2 The a-th ETN belonging to the GMTN from Fig. 1 (a), and its SETN's representation specified by (a) SC (b) the set of resources, SL = P, U, T, Θ 0 − -characterizes the structure of local processes of SCCP, i.e.P-the set of local processes, U -the set of local process routes, T -the set of local process operations times, Θ 0the set of dispatching priority rules.SM = m P, mΠ, mT, Θ 1 -characterizes the structure of the multimodal processes of the SCCP, i.e. m P-the set of multimodal processes, M-the set of multimodal process routes, mT -the set of multimodal process operation times, Θ 1 -the set of dispatching priority rules.
(a) SC corresponds to exactly one resource R ω ∈ Rp b of the structure (b) SC,• each process P i ∈ Pp a /m P i ∈ m Pp a (local as well as multimodal) of the substructure(a) SC corresponds to exactly one processP μ ∈ Pp b /m P μ ∈ m Pp b of the structure (b) SC, • routes p i ∈ U p a /mp i ∈ mPp a and p μ ∈ U p b /mp μ ∈ m Pp b of the corresponding processes are sequences consisting of corresponding resources, • each operation o k i, j /mo k i, j executed within the substructure (a) SC corresponds to exactly one operation o ϕ μ,ν /mo ϕ μ,υ executed within the substructure (b) SC, • dispatching rules σ v c of the resources are sequences consisting of elements s v c,d indicating the corresponding streams of processes.

Fig. 5
Fig. 5 Space of admissible schedules L including L instances following the SETN from Fig. 2b parameterized by S and E (a), "top-down view"-elevation of L onto plane S × E (b)

Fig. 6
Fig. 6 Part of the SCCP from Fig. 1b (a), the a-th SETN specified by (i) SC (b), the convoluted representation of (a) SC (c)

Fig. 11
Fig. 11Cyclic schedule for multimodal process m P 2 from Fig.10

Table 3
Dispatching rules allotted to shared resources in(a)SC from Fig.6c

Table 5
Resolution time of the cyclic steady state reachability problem for different scale SCCPs * SCCP from Fig.