1. the control of elementary siphons
2. finding the relevant approach to finding conditions for control of dependent siphons;
3. finding the way of making the ES3PR deadlock-free.
3.2. Procedure of Setting the Supervisor for ES3PR
In [
6,
7] it can be found the following definition:
Defnition 15 [
6,
7]: Let
S be a strict minimal siphon in an ES
3PR net model of a plant (
Nµ0, Mµ0) where
Nµ0 = ◯
i∈{1, …, n} Ni = (
P0 ∪
PS ∪
PR, T, Fµ0,
Wµ0). Let {
α, β, …, γ} ⊆ {1, 2, …,
n} such that
i ∈ {
α, β, …, γ}, |
ThS| ∩
p ∈
PSi ≠ ∅ and ∀
j ∈ {1, 2, …,
n} \ {
α, β, …, γ}, |
ThS| ∩ p ∈
PSj ≠ ∅. For
S, a nonnegative
P-vector
kS is constructed as follows:
Step 1: ∀p ∈ P0 ∪ PS ∪ PR, kS(p): = 0;
Step 2: ∀p ∈ |Th(S)|, kS(p): = hS(p), where Th(S) = Σp∈|Th(S)| hS(p)p;
Step 3: ∀i ∈ {α, β, …, γ}, let ps ∈ ||ThS|| ∩ PSi be such a place that ∀pt ∈ SP (pu, p0i), pu ∈ ps••, pt ∉ |T h(S)|. Suppose that there are m such places p1s, p2s, …, pms. Assuredly, we have {pis |i = 1, 2, …, m} ⊆ |Th(S)| ∩ PSi. ∀pis let piν ∈ SP (p0i, pis) be such a place that hS(piν) > hS(pw), ∀pw ∈ SP (p0i, pis). ∀px ∈ SP (p0i, pis), kS(px): = hS(piν).
Step 4: ∀Im’ ⊆ {1, 2, …, m}. ∀py ∈ ∩ i∈Im’ SP (p0i, pis), kS(py): = hS(piz), where piz ∈ |Th(S)| ∩ PSi, and ∄p ∈ Th(S) ∩ PSi, s.t. hS(p) > hS(piz).
In general, such procedure makes possible to finalize the supervisor synthesis, especially to set markings of the supervisor monitors VSi, i = 1, 2, …
3.3. Illustrative Example of ES3PR—Siphons and P-Invariants
Let us apply the theory of siphon-based control, introduced above, on the example of ES
3PR. Consider the GPN in
Figure 1 modelling an AMS. Here, the set of operation places is
PS = {
p1,
p2,
p3,
p4, p5,
p7,
p8,
p9}, the set of resource places is
PR = {
p11,
p12,
p13,
p14, p15, p16}, and the set of idle places is
P0 = {
p6,
p10}.
As we can see in
Figure 1, there exists one arc with the weight equal to 2
—notice the arc from
t6 to
p13. Consequently, in
Figure 1 is GPN. The state space of the net
N expressed by the reachability tree (RT) has 308 nodes (including the initial state). Because of such large amount of nodes, the RT cannot be displayed here. The minimal siphons are introduced in
Table 1.
As we can see in
Table 1, there exist 11 minimal siphons in this GPN model of AMS. They can be enumerated e.g., by the tool GPenSIM [
13,
14,
15] or by another PN tools. Because the traps corresponding to the crossed siphons are marked, crossed siphons may be omitted from the list of siphons. Hence, only 4 of them are strict minimal siphons—
S3, S7, S8 and
S10. Other 7 siphons are not relevant because they are equal to corresponding traps.
Let us renumber siphons S3, S7, S8 and S10 to siphons S3, S1, S2 and S4, respectively. Consequently, S1 = {p5, p8, p15, p16}, S2 = {p4, p9, p14, p15}, S3 = {p4, p13} and S4 = {p5, p9, p14, p15, p16}. In the form of row P-vectors they are as follows:
S1 = (0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1)
S2 = (0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0)
S3 = (0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0)
S4 = (0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 1)
Corresponding matrix [
λ] = [
S1T, S2T, S3T, S4T]
T can be utilized (analogically to [
1] in the case of S
3PR paradigm) for finding potential interconnections with the supervisor eliminating the deadlocks. Namely, when we multiply matrix [
λ] by [
N], being (16 × 11) dimensional incidence matrix of GPN in
Figure 1, we obtain the matrix [
η]
= [
λ]. [
N]
= [
η1T, η2T, η3T, η4T]
Tconsisting of the following
T-vectors
η1 = (0, 0, 0, 0, −1, 1, 0, −1, 1, 0, 0)
η2 = (0, 0, 0, −1, 1, 0, 0, 0, −1, 1, 0)
η3 = (0, 0, 0, −1, 0, 1, 0, 0, 0, 0, 0)
η4 = (0, 0, 0, −1, 0, 1, 0, −1, 0, 1, 0).
Hence, we have λS1 = p5 + p8 + p15 + p16, λS2 = p4 + p9 + p14 + p15, λS3 = p4 + p13, λS4 = p5 + p9 + p14 + p15 + p16 and ηS1 = − t5 + t6 − t8 + t9, ηS2 = − t4 + t5 − t9 + t10, ηS3 = − t4 + t6, and ηS4 = − t4 + t6 − t8 + t10. It is easy to verify that ηS4 = ηS1 + ηS2. However, this set of vectors ηSi , i = 1, …, 4 is not linearly independent, because η4 = η1 + η2, i.e., rank([η]) = 3, and rank([η]) = NES = 3 < rank([N]) = 8. It means that there are three elementary siphons and one strongly dependent siphon. Thus, ΠE = {S1, S2, S3} and ΠD= {S4}.
In
Figure 1 we have
p01 = p6,
p02 = p10,
PS1 = {
p1, …, p5},
PS2 = {
p7, …, p9},
PR1 = {
p11, …, p16}, and
PR2 = {
p14, …, p16}. Unfortunately, the system is deadlocked.
GPN in
Figure 1 has eight minimal
P-invariants (see
Table 2). They may also be computed e.g., by PN tools [
13,
14,
15]. Namely,
I1 = p2 + p11,
I2 = p1 + p12,
I3 = p3 + 2p4 + p13,
I4 = p3 + p9 + p14,
I5 = p4 + p8 + p15,
I6 = p5 + p7 + p16,
I7 = p1 + p2 + p3 + p4 + p5 + p6 and
I8 = p7 + p8 + p9 + p10.
P-invariants in the form of row vectors we can see in
Table 2. Negative elements in the vectors
ηi in S
3PR paradigm (see [
1]) mean that directed arcs are emerging from the GPN model and enter through corresponding transition
ti (index
i depends on the position of the transition in a vector
ηj) to the supervisor, while positive elements represent directed arcs in opposite direction, i.e., from the supervisor to GPN model through corresponding transitions. However in case of ES
3PR as well as S
4PR paradigms the supervisor synthesis is not so simple like in S
3PR paradigm. The structural properties presented here will be applied in the continuation of this Example introduced in the
Section 4.3.
If we applied these parameters (which are suitable for S3PR) to ES3PR, we would find that the responding supervisor does not comply for such paradigm of the AMS model, because it does not prevent deadlocks.