
Citation: Esther Menendez, Paula Garcia-Fraile, Raul Rivas. Biotechnological applications of bacterial cellulases[J]. AIMS Bioengineering, 2015, 2(3): 163-182. doi: 10.3934/bioeng.2015.3.163
[1] | Guohui Zhang, Jinghe Sun, Xing Liu, Guodong Wang, Yangyang Yang . Solving flexible job shop scheduling problems with transportation time based on improved genetic algorithm. Mathematical Biosciences and Engineering, 2019, 16(3): 1334-1347. doi: 10.3934/mbe.2019065 |
[2] | Kongfu Hu, Lei Wang, Jingcao Cai, Long Cheng . An improved genetic algorithm with dynamic neighborhood search for job shop scheduling problem. Mathematical Biosciences and Engineering, 2023, 20(9): 17407-17427. doi: 10.3934/mbe.2023774 |
[3] | Jianguo Duan, Mengting Wang, Qinglei Zhang, Jiyun Qin . Distributed shop scheduling: A comprehensive review on classifications, models and algorithms. Mathematical Biosciences and Engineering, 2023, 20(8): 15265-15308. doi: 10.3934/mbe.2023683 |
[4] | Cong Zhao, Na Deng . An actor-critic framework based on deep reinforcement learning for addressing flexible job shop scheduling problems. Mathematical Biosciences and Engineering, 2024, 21(1): 1445-1471. doi: 10.3934/mbe.2024062 |
[5] | Zilong Zhuang, Zhiyao Lu, Zizhao Huang, Chengliang Liu, Wei Qin . A novel complex network based dynamic rule selection approach for open shop scheduling problem with release dates. Mathematical Biosciences and Engineering, 2019, 16(5): 4491-4505. doi: 10.3934/mbe.2019224 |
[6] | Yifan Gu, Hua Xu, Jinfeng Yang, Rui Li . An improved memetic algorithm to solve the energy-efficient distributed flexible job shop scheduling problem with transportation and start-stop constraints. Mathematical Biosciences and Engineering, 2023, 20(12): 21467-21498. doi: 10.3934/mbe.2023950 |
[7] | Shaofeng Yan, Guohui Zhang, Jinghe Sun, Wenqiang Zhang . An improved ant colony optimization for solving the flexible job shop scheduling problem with multiple time constraints. Mathematical Biosciences and Engineering, 2023, 20(4): 7519-7547. doi: 10.3934/mbe.2023325 |
[8] | Wenqiang Zhang, Chen Li, Mitsuo Gen, Weidong Yang, Zhongwei Zhang, Guohui Zhang . Multiobjective particle swarm optimization with direction search and differential evolution for distributed flow-shop scheduling problem. Mathematical Biosciences and Engineering, 2022, 19(9): 8833-8865. doi: 10.3934/mbe.2022410 |
[9] | Xuyin Wang, Weiguo Liu, Lu Li, Peizhen Zhao, Ruifeng Zhang . Resource dependent scheduling with truncated learning effects. Mathematical Biosciences and Engineering, 2022, 19(6): 5957-5967. doi: 10.3934/mbe.2022278 |
[10] | Yufan Zheng, Rafiq Ahmad . Feature extraction and process planning of integrated hybrid additive-subtractive system for remanufacturing. Mathematical Biosciences and Engineering, 2020, 17(6): 7274-7301. doi: 10.3934/mbe.2020373 |
Remanufacturing, as an industrial process of restoring discarded products/components back to their useful lives, is of growing importance due to the emerging pressure of legislation and increasing awareness of environmental conservation. As an ultimate form of recycling, remanufacturing maintains much of the value added from original material and reprocessing conservation, leading to lower production costs and improved profits [1]. A recent wall street journal revealed that large scale remanufacturing in the United States employs more than 500, 000 people and contributes to approximately $100 billion of goods sold each year [2]. More successful industry applications can be seen in automobile remanufacturing, aerospace remanufacturing and electronic remanufacturing [3, 4, 5].
Compared to manufacturing, remanufacturing is more complicated in the way that (1) the supply of returned products is unpredictable in timing and quantities; (2) the quality and composition of returned products varies; and (3) the process routings are not necessarily fixed but rather adapt to the actual conditions of products/components [6]. Such characteristics have led many new methodologies to deal with various operation management issues in remanufacturing as summarized in the survey literature [7]. One of the challenges, which is the focus of this paper, is remanufacturing job shop scheduling. With the high level of uncertainty and resource limitation, it becomes extremely difficult to carry out various remanufacturing tasks for different products and still achieve maximal system performance.
Guide and his colleagues are the first group of researchers that looked into the remanufacturing scheduling. Through simulation, they examined the performance of several static dispatching rules in a repair shop. They assumed the system has stochastic arrivals, exponentially distributed process times and probabilistic process routings [8]. Extended from this work, an analytical model was further developed to explore the impact of those rules on system performance in terms of weighted average sojourn time [9]. While the pursued simulation and model analysis provided insights to practicing managers in a repair shop, much of work is needed for complex job shop settings where various products/parts with uncertain remanufacturing routings compete for constrained resources.
A second line of work has focused on economic lot scheduling with random returns. For instance, Tang and Teunter considered a hybrid system where (re)manufacturing operations of multiple products are performed on the same production line [10]. The work formulated the problem as a mixed integer linear program to determine the optimal cycle time and production starting times for (re)manufacturing. The research was further extended to relax the constraint of common cycle time and a single (re)manufacturing lot per item per cycle [11]. By dividing the returns into different quality grades, Sun et al. investigated a lot scheduling problem where the remanufacturing rate increases as the quality grade increases and holding costs for serviceable products are higher than returns [12]. The objective is to minimize the average total cost by optimizing the acquisition lot size and scheduling the remanufacturing sequence. However, the deterministic nature of the assumptions (i.e., constant demand and returns, and constant remanufacturing rate) made in those works failed to address the uncertainty inherent in remanufacturing. Moreover, the authors considered remanufacturing as a single-stage process. Given that refurbishing used products up to like-new condition tends to involve a sequence of recovery operations, such single-stage scheduling methods present practical limitation.
A third line of researches have tackled the multi-stage and multi-product remanufacturing scheduling problem. For instance, Kim et al. considered a remanufacturing system with a single disassembly workstation, parallel flow-shop-type reprocessing lines, and a single reassembly workstation [13]. A hybrid scheduling heuristic is proposed to determine the sequences of products to be disassembled, to be reprocessed at each reprocessing lines and to be reassembled. The authors, however, treated the problem deterministic and static. Luh et al. presented another work that considered a similar remanufacturing system but tried to address the underlying uncertainty [14]. The main novelty of the research is that it considered stochastic asset arrival time and part repair time to model the dynamics of rotable inventory for optimization. Nevertheless, the authors presumed that each asset goes through a fixed series of overhaul operations. Giglio et al. studied a production system which produces multi-class single-level products through both manufacturing of raw materials and remanufacturing of return products [15]. A mixed-integer programming formulation is proposed to obtain the optimal lot size and job shop schedules with minimum total cost. This work presented the same limitation that the recovery process routings of returns are deterministic. By considering the alternatives of remanufacturing process routes for worn cores, Zhang et al. proposed a simulation based genetic algorithm approach for remanufacturing process planning and scheduling, where the processing time of each alternative routes are assume to be known a priori [16].
Summarizing the findings of the above discussion, no study has comprehensively dealt with remanufacturing job shop scheduling in subjecting to alternative recovery operations, random operation times and limited resource conflicts. In a practical remanufacturing system, the quality of returned products/cores varies, ranging from slightly used with minor blemishes to significantly damaged and requiring extensive repair. Such differences in the condition of products/cores strongly affect the set of recovery operations necessary to bring them up to a quality standard. In addition, the time that each recovery operation takes varies as well with respect to different damages of the cores. Those stochastic natures in remanufacturing triggers the potential conflicts of limited machine resources to a noticeable extend. Such characteristics have led a number of studies to investigate the operation management issues in remanufacturing environment by using different modeling methods, i.e., through analytic model [8, 9, 11, 12, 17], through Mixed Integer Programming model [10, 13, 15, 16], through Lagrangian Relaxation [14] and through Petri nets model [18, 19]. While our previous work [17] and [20] investigated the remanufacturing processes of worn cores with uncertain failure conditions and proposed an analytical model to study such dynamics in recovery process routings through probabilistic measures, how to incorporate those dynamic measures into remanufacturing scheduling has yet to be undertaken.
The colored timed Petri nets (CTPN) is an extension of ordinary Petri nets where the time and color are introduced to describe the logic structure and dynamic behavior of complex systems [21]. By taking advantage of both the well-formed formalism of CTPN and the robustness of Simulated Annealing (SA) for global optimization, the paper tackles this challenge and makes the following contributions. First, a colored timed Petri net is designed to explicitly model the dynamic of remanufacturing processes, such as various process routings, different uncertain operation times for cores, and real-time machine resource conflicts. With the color attributes in Petri nets, two types of decision points, recovery routing selection and resource dispatching, are introduced and embedded in places of CTPN model. With time attributes in Petri nets, the temporal aspect of recovery operations for cores as well as the evolution dynamics in cores' operational stages is mathematically analyzed. Second, a hybrid meta-heuristic algorithm embedded scheduling strategy over CTPN is proposed to search for the optimal recovery routings for worn cores and their recovery operation sequences on workstations, in minimizing the total production cost.
The rest of the paper is organized as follows. Section 2 presents the CTPN model. Section 3 describes the hybrid meta-heuristics based on SA and minimum slack time (MST) dispatching rule. A case study is conducted in Section 4, followed by the conclusion and future research in Section 5.
This paper considers a typical remanufacturing job shop where its clients drop their used products with an expectation to get them fully recovered within a fixed time frame. The remanufacturing processes are traditionally complicated including disassembly that disassembles a worn product into components/cores, cores' recovery operations and assembly that assembles recovered cores back into the product. Figure 1 gives an organizational structure of the remanufacturing job shop.
In this case, each used product (e.g., the lth used product) arrives at the job shop with an arrival time AT(l) and a due date Due(l). The used product is first disassembled into what so called "cores" that are disassembled subassemblies or components. The cores are then classified into "reusable", "repairable", and "disposal" through the inspection process as depicted in Figure 1. This paper focuses on the recovery activities only. The challenge arises when dealing with the repairable cores with different levels of damage that need to go through a set of recovery operations to restore their size and property.
In a remanufacturing job shop, the quality of repairable cores commonly ranges from slightly used with minor blemishes to significantly damaged and requiring extensive recovery. Such quality variations strongly affect the set of recovery operations necessary to bring them up to a quality standard. As exemplified in Figure 2, there are several process routings for the remanufacturing of used machine tool cores of different types and various damage conditions. For instance, the lathe spindles can be inspected and classified into either "abrasion", "corrosion", or "micro-cracks", on which a set of specific recovery operations are designated for their remanufacturing.
With a wide variety of recovery technologies being applied into remanufacturing industry, some recovery process routings for cores with a certain type of damage can be alternative. For instance, if a spindle is detected with severe abrasion, two typical process routings would be "grinding → chromium electroplating" or "grinding → cold welding → fine grind" to regain the surface accuracy. Each alternative process routing consists of a set of predesigned recovery operations, with each operation being performed by a certain kind of machine equipment in a workstation, resulting in a specified recovery quality and unit processing cost [22, 23]. As an upper level of process planning, different combinations of recovery process routings for cores would lead to much variation in remanufacturing cost and real-time machine workload in remanufacturing job shop.
With the recovery routings determined in the upper level, the core then travels to and is processed by a set of recovery workstations sequentially according to the designed process flow. However, the recovery operation time for a particular core in a given workstation is not fixed but rather dependent on the actual condition of the core. While different routings might demand operation times on the same workstation, various cores would compete for constrained machine resource to be recovered.
Our objective is to consider such variation for efficient resource utilization in fulfilling the demand with the minimum cost (i.e., operation cost and tardiness penalty). The remanufacturing job shop scheduling problem is strongly NP-hard due to: a) assignment decisions of reparable cores to recovery process routings and b) sequencing decisions of recovery operations of cores in each workstation. For that purpose, the paper makes the following assumptions:
(1) Each workstation has only one machine and each machine can process at most one task at a time.
(2) Each task is non-preemptive, requiring one and only one machine at a time.
(3) The transportation time between buffers is negligible.
(4) Penalty is charged if the system does not meet the due date of a product.
The ordinary of Petri nets (PN) does not include time and color, which limits its capability to describe the logic structure and behavior of modeled systems, but not its evolution over time and color [24, 25]. In view of the specific characteristics of remanufacturing system to be modeled, this section introduces time and color attributes to a Petri net and proposes a CTPN. The model considers not only recovery process routings and stochastic operation times, but also integrates recovery routing assignments and efficient resource dispatching.
Definition 1: A colored timed Petri net (CTPN) is defined as 7-tuple: CTPN = (P, T, I, O, M, U, C).
(1) P = S ∪ W, where the set of places in S represents buffers, the set of places in W stands for workstations.
(2) T = Ta ∪ Te, where the finite set of transitions in Ta represent recovery operations, and the ones in Te stand for transportation.
(3) I: P×T→{0, 1} is the input function that defines the set of ordered pairs (pi, tk), where I(pi, tk) = 1, if pi is an input place for tk (i.e., ∀pi $ \in $·tk), otherwise 0.
(4) O: P×T→{0, 1} is the output function that defines the set of ordered pairs (pi, tk), where O(pi, tk) = 1, if pi is an output place for tk (i.e., ∀pi $ \in $·tk·), otherwise 0.
(5) M: P→{0, 1, 2, …} is a marking vector, where M(pi) represents the number of tokens in pi.
(6) U: Ta→Θ is a non-zero time function that assigns to each recovery operation t∈Ta, while U: Tb→Θ is a zero time function such that the transportation time between buffers is negligible.
(7) C: P→Σ is a color function that assigns to each place p∈ S ∪W with a color set C(p).
Figure 3 gives an example of the proposed CTPN that corresponds to recovery process routings for used machine tools in Figure 2. The description for the places and transitions in CTPN are given Table 1. For the modeling and optimization purpose, some basic notations are introduced. Table 2 shows the index used in the CTPN model. Table 3 and Table 4 gives the model parameters and decision variables in the CTPN, respectively.
Places and transitions in CTPN | Description |
s0 | The output buffer of the inspection center |
si (i=1, 2, ..., 12) | The input buffer of each recovery workstation |
sz (z=13, 14, 15, 16) | The output buffer of each recovery workstation |
s17 | The input buffer of the reassembly center |
wk (k=1, 2, ..., 12) | Workstation places |
ti, j (i, j=1, 2, ..., 12, i≠j) | Recovery operations |
t0, i, ti, z, ti,17 (i=1, 2, ..., 12, z=13, 14, 15, 16) | Transportation transitions |
Index | Description |
L= {1, 2, ..., L} | index set of used products |
D= {1, 2, ..., D} | index set of component types |
K= {1, 2, ..., K} | index set of workstations |
H= {1, 2, ..., H} | index set of recovery process routings |
Parameters | Description |
ϑld | the dth core disassembled from the lth product |
wk | the kth workstation that performs a specific recovery operation |
Rh | the hth recovery process routing designed for worn cores and is represented as a set of recovery operations |
pre(tij, Rh) | the operational stages of recovery operation tij in Rh |
Last(Rh) | the last recovery operation in Rh |
aldh | a binary variable with a value 1 indicates the hth routing for the core ϑld |
u(tij, ϑld) | the time for a core ϑld to be processed in transition tij |
AT(si, ϑld) | the arrival time of the core ϑld at buffer place si |
ST(wk, ϑld) | the start time of the core ϑld being processed by wk |
Com(wk, ϑld) | the completion time of the core ϑld processed by wk |
ρ(ϑld) | the makespan of the component ϑld |
χ(ϑld) | the tardiness of the component ϑld |
$\zeta (l)$ | the tardiness of the lth product |
Probld | the probability of a core ϑld being selected for routing reassignment |
$\Pi _{ld}^{h, x}$ | the difference of the total average waiting time between the hth routing and the xth one for the core ϑld |
WTk | the average waiting time per core consumed in the kth workstation |
h | the total average waiting time per core consumed in the hth routing |
$v_{ld}^{h, x}$ | the probability of the xth routing being selected to replace the current hth routing for the core ϑld |
Slackxy | slack time of a core ϑxy in scheduling horizon |
Decision variables | Description |
rldh | a binary variable with a value 1 indicates that the hth process routing is assigned to a core ϑld |
Pri(wk, ϑld) | the priority for the dth component of the lthused product to be processed by the kth workstation |
(1) Modeling of recovery process routings
As stated earlier, the recovery process routing for a particular core is not deterministic but rather contingent on its damage level. When a core goes through inspection center and enter in the buffer place s0, it will be assigned with a specified process routing (i.e., Rh), which is represented as a set of recovery operations that a core should go through sequentially for its recovery.
For the cores with alternative process routings to be recovered, decisions should be made on which "reasonable" process routings to be assigned to them. In the proposed CTPN model, the buffer place s0 operates as a decision agent of choosing optimal recovery process routings for cores, by utilizing a hybrid scheduling algorithm elaborated in Section 3.
When the process routings for cores are determined, cores of different types and stochastic damage conditions require going through various recovery operations. In the CTPN model of Figure 3, different types of core tokens in s0 might go to transitions among t0, 12, t0, 5, t0, 1, t0, 3, t0, 2, t0, 6, t0, 8, t0, 11; tokens in s13 may either go to t13, 3, t13, 10 or t13, 17; tokens in s16 might enters into either t16, 17 or t16, 5. In addition, cores of the same type and the same damage condition might go through different recovery processes. For example, when a spindle is detected with abrasion, its process routing can either be R1 = {t03, t3, 14, t14, 9, t9, 17} or R2 = {t03, t3, 14, t14, 4, t4, 15, t15, 5, t5, 17}. Then in Figure 3, a token in s14 representing a spindle with abrasion might fire transition t14, 9 or t14, 4.
Thus, to make sure where the tokens in a buffer place should go, a unique color is introduced and associated with buffer places to distinguish the cores going through different recovery operations.
Definition 2: The color of a token ϑld in a buffer place si is defined as C(sj, ϑld) = cij if the token goes to tij after being released from si.
Let M(si) denote the number of tokens in place si regardless of token color, and M(si, cij) the number of tokens in place si at M that have the color cij. Therefore, $M({s_i}) = \sum\limits_j {M({s_i}, {c_{ij}})} $.
Based on Definition 2, the color attached to a token in a buffer place changes through transition firings. Such color evolution represents the complex operational stages of a core according to its recovery process routing. For instance, if the recovery process routing of the core ϑld is determined (i.e., rldh is known), its color evolution can then be derived as follows.
$C({{s}_{j}},{{\vartheta }_{ld}})={{c}_{jp}},\ \text{if}\ \left\{ C(si,ϑld)=cijpre(tjp,Rh)=pre(tij,Rh)+1, ∀rldh=1si∈∙tij, sj∈∙tjp \right.$
|
(1) |
(2) Modeling of shared resources
In a remanufacturing job shop, various cores would compete for constrained resource to be recovered. An input buffer is then designated for each workstation in our CTPN model to accommodate all the cores that are waiting for processing. When multiple tokens occur in the input buffer place sk of the kth workstation, decisions should be made on which core to be processed next as soon as the machine resource is ready. In the proposed CTPN model, each input buffer place si (i = 1, 2, ..., 12) is considered as a decision agent for resource dispatching and determines the operation sequences (priority) for the cores to be processed by the workstation.
Definition 3: The priority for the core ϑld to be processed by the kth workstation is defined as Pri(wk, ϑld). The priority is a decision variable for scheduling optimization. The smaller the value is, the higher priority for the core to be processed by workstation wk.
Definition 4: A transition tij in the CTPN developed above is enabled at a marking M to process the core ϑld with a color cij, if and only if:
$M({s_i},{c_{ij}}) \geqslant 1{,^{}}{s_i} \in {}^ \bullet {t_{ij}}$ | (2) |
$M({w_k}) = 1{,^{}}{w_k} \in {}^ \bullet {t_{ij}} \cap t_{ij}^ \bullet $ | (3) |
$Pri({w_k},{\vartheta _{ld}}) = \min (Pri({w_k},{\vartheta _{xy}})) \;{\rm{and}} \;C({s_i},{\vartheta _{xy}}) = {c_{ij}},\forall x{,^{}}y$ | (4) |
Definition 5: An enabled transition tij can fire in a marking M with respect to a color cij yielding a new marking M' denoted as $M[({t_{ij}}, {c_{ij}}) > M'$:
$M'({s_i},{c_{ij}}) = M({s_i},{c_{ij}}) - 1{,^{}}{s_i} \in {}^ \bullet {t_{ij}}$ | (5) |
$M'({w_k}) = M'({w_k}){,^{}}{w_k} \in {}^ \bullet {t_{ij}} \cap t_{ij}^ \bullet $ | (6) |
Definition 6: When a transition tij is fired, the marking M' would transfer to M" according to:
$M''({s_i},{c_{ij}}) = M'({s_i},{c_{ij}}){,^{}}{s_i} \in {}^ \bullet {t_{ij}}$ | (7) |
$M''({w_k}) = M'({w_k}) + 1{,^{}}{w_k} \in {}^ \bullet {t_{ij}} \cap t_{ij}^ \bullet $ | (8) |
$M''({{s}_{j}},{{c}_{jp}})=M'({{s}_{j}},{{c}_{jp}})+1\ \text{and}\ C({{s}_{j}},{{\vartheta }_{ld}})={{c}_{jp}}\ \text{if}\left\{ sj∈t∙ij, sj∈∙tjppre(tjp,Rh)=pre(tij,Rh)+1 \right.$
|
(9) |
According to the above transition enabling and firing rules, the system always dispatch an available workstation wk to the core with the highest priority (i.e., the lowest Pri(wk, ϑld) value) in the queue to maximize the system performance.
(3) Modeling of recovery operation times
When associating time with operation transitions, the CTPN model is able to describe the temporal aspect of recovery operations.
Definition 7: The time needed for a token ϑld to be processed in transition tij is defined as u(tij, ϑld).
The processing times required to perform a necessary recovery operation is highly variable since the quality and composition of returned products/parts varies. It is reported that the processing time data obtained from a remanufacturing facility follows an exponential distribution in order to simulate the large range of possible values [8]. The amount of processing time for a remanufacturing operation is widely modeled as an exponential function of the wear associated with a part [19]. In this section, we models the recovery operation time as the sum of average operation time and extended time, where the latter is a exponential function of the quality condition of the core being processed. To that end, quality testing of a core is conducted before recovery, on which an inspection score, Score(ϑld)∈[0, 1], is assigned to the core based. It is assumed that the better the quality, the higher the score and the shorter the operation time. So the relationship between the quality score and the processing time is modeled as follows.
$u({t_{ij}},{\vartheta _{ld}}) = \frac{{ - \ln (Score({\vartheta _{ld}}))}}{{\beta ({w_k})}} + \lambda ({w_k}){,_{}}{t_{ij}} \in w_k^ \bullet { \cap ^ \bullet }{w_k}$ | (10) |
where β(wk) is a control factor for workstation wk whose value is determined based on statistics of historical data, and λ(wk) is the average operation time of cores in wk. Whenever a component goes into a transition for processing, it must stay in transition tj for more thanu(tj, ϑld) time units and then can be released into another buffer place.
(4) Temporal aspect of recovery operations
In the CTPN model, with the core tokens dynamically entering the buffer places, a set of transitions are then enabled and fired according to the enabling and firing rules. With the dynamic behavior of the CTPN, the remanufacturing operations are subject to the following constraints.
Operation start time requirements: Recovery operations cannot be started until the core has arrived at the buffer sk of a corresponding workstation wk. When the core arrived in the buffer sk, it need waiting until all other cores with higher priorities are processed.
$ ST({w_k},{\vartheta _{ld}}) = {\text{max\{ }}A{\text{T}}({s_i},{\vartheta _{ld}}),Com({w_k},{\vartheta _{ij}}){\text{\} }}{{\text{,}}^{}}\forall Pri({w_k},{\vartheta _{ld}}) = Pri({w_k},{\vartheta _{ij}}) + 1,{s_i} \in {}^ \bullet {t_{ij}},{t_{ij}} \in {}^ \bullet {w_k} \cap w_k^ \bullet $ | (11) |
Processing completion time requirements: When a core goes into a workstation for processing, it must stay in a specific machine for more thanu(ϑld, wk) time units and then can be released into the successor buffer.
$Com({w_k},{\vartheta _{ld}}) = ST({w_k},{\vartheta _{ld}}) + u({t_{ij}},{\vartheta _{ld}}){,^{}}\forall {t_{ij}} \in w_k^ \bullet { \cap ^ \bullet }{w_k}$ | (12) |
Arrival time constraints: When a core is finished processing by a workstation, it can be transported immediately into another buffer of the successor workstation by ignoring the transportation time.
$A{\rm{T}}({s_i},{\vartheta _{ld}}) = Com({w_k},{\vartheta _{ld}}),\;\forall \left\{ rldh=1pre(wv,Rh)=pre(wk,Rh)+1si∈∙tij,tij∈∙wv∩w∙v \right.$
|
(13) |
Using Eqs 11, 12 and 13 recursively, the makespan and tardiness of the component ϑld, i.e., ρ(ϑld) and $\chi ({\vartheta _{ld}})$, can be derived as shown in Eqs 14 and 15. The tardiness of the lth product is then calculated in Eq 16.
$\rho ({\vartheta _{ld}}) = Com({w_k},{\vartheta _{ld}}) - AT(l){,^{}}\forall Last({R_h}) = {t_{ij}}{,^{}}{t_{ij}} \in w_k^ \bullet { \cap ^ \bullet }{w_k}$ | (14) |
$\chi ({\vartheta _{ld}}) = \max \{ 0{,^{}}\rho ({\vartheta _{ld}}) - Due(l)\} $ | (15) |
$\zeta (l) = \mathop {\max }\limits_d \chi ({\vartheta _{ld}})$ | (16) |
Generally speaking, there are two ways to solve a scheduling problem by PN. One is to construct their reachability graph. However, it suffers from the state explosion problem [26]. Since we considered two types of decision variables in the proposed CTPN, using reachability graph for scheduling optimization might not be efficient. The other is to introduce scheduling algorithms into a PN model, which is adopted in this paper. The embedding style is optionally through places or transitions in a PN. This work uses the former mechanism and adopts two types of decision places in the proposed CTPN model. The decisions of choosing optimal recovery process routings for cores are embedded in buffer place s0, while the input buffer place si (I = 1, 2, ..., 12) makes decisions on dispatching machine resource in workstation to perform recovery operations for the cores that are waiting in queue. Considering the two types of scheduling decisions, a hybrid meta-heuristic algorithm using simulated annealing and minimum slack time rule is proposed and embedded in these decision places to search for global optimal process plans and schedules with a quick convergence speed.
The scheduling objective is to minimize both operation cost and tardiness penalty as shown in Eq 17. There are two types of decision variables: (1) recovery process routing for a given core (i.e., rldh), and (2) operation sequences for a given core in a workstation (i.e., Pri(wk, ϑld)).
Minimize:
$
TC = \sum\limits_l {\sum\limits_d {\sum\limits_{} {u({t_{ij}},{\vartheta _{ld}}) \cdot \phi ({w_k})} + \sum\limits_l {\zeta (l) \times \varpi } } } ,\;\forall \left\{ rldh=1tij∈Rhandtij∈w∙k \right.
$
|
(17) |
where φ(wk) is unit operation cost of the kth workstation and ϖ is unit penalty cost for product.
SA is a stochastic gradient method for global optimization that has been applied to a wide variety of sophisticated combinatorial problems [27, 28, 29]. While it conducts local searches, it is capable of exploring the solution space stochastically to prevent from being trapped in a local optimum. Starting from an initial solution g0, SA uses a certain mechanism to generate a neighborhood solution g' with an acceptance probability AP.
$AP = \left\{ 1,ifΔ<0exp(−Δ/Ω),otherwise \right.$
|
(18) |
where Δ = f(g')-f(g0), Ω is called temperature, a global time-varying parameter, and f(g) the optimization objective function. In this paper, SA is used to determine the optimal process routings for cores. At each iteration when a neighbor solution is generated, MST to be elaborated in section 3.2 is then applied to determine the operation sequences of the cores in each workstation. The parameter Ω is reduced by a factor η at each iteration and the chance of choosing an inferior solution decreases as well. The nested SA-MST search process continues until the stopping criterion is met.
In this section, the CTPN model is integrated and embedded with SA algorithm and MST dispatching rules to obtain the optimal remanufacturing routings and schedules. At each iteration of SA algorithm, the feasible recovery routing is first generated by the Initial/Neighborhood Solution Generation process. Then the CTPN model is established based on the arrival time of parts, the predetermined recovery routings and the processing time of each operations. Whenever multiple parts arrived in a same buffer place in the CTPN, MST rules is used for resource dispatching. The CTPN model evolves to a final version until the resource conflicts in all the buffer places are tackled by the SMT. Finally, the operation cost and tardiness penalty of a SA solution at each iteration can be calculated according to the time aspects marked in the CTPN, and the next iteration of SA algorithm goes on until the stopping criterion is met. The flowchart of the CTPN-based hybrid meta-heuristics scheduler is shown in Figure 4.
(1) SA solution representation
In this paper, a 3-dimension matrix R = [rldh]L×D×H is proposed as representation of SA solutions. Each element rldh is a binary variable, one indicating the hth recovery process routing is assigned to the core ϑld and zero otherwise. Given that not all routings are feasible to a specific core, the viability of a randomly generated R is controlled via a logic AND operation in Eq 19:
$ R = \text{ }R\wedge A = \text{ }[{{r}_{ldh}}\wedge {{a}_{ldh}}] $ | (19) |
where A = [aldh] is a L-by-D-by-H mask whose element aldh = 1 indicates the feasibility of the hth routing for the core ϑld.
(2) Initial solution generation
The choice of an initial solution is vital to the performance of the algorithm [30]. To that end, a cost function of a routing (i.e., the hth routing) for a core ϑld is defined in Eq 20. The smaller the value of qldh, the greater chance the hth routing is selected for the core, i.e., rldh = 1.
${q_{ldh}} = \frac{{\sum\limits_j {u({t_{ij}},{\vartheta _{ld}}) \times \phi ({w_k})} }}{{\sum\limits_h {\sum\limits_j {u({t_{ij}},{\vartheta _{ld}}) \times \phi ({w_k})} } }},{t_{ij}} \in {R_h},{r_{ldh}} = 1,{t_{ij}} \in {}^ \bullet {w_k} \;{\rm{I}}\; w_k^ \bullet $ | (20) |
(3) Neighborhood solution generation
Traditional SA randomly generates a neighborhood solution g' from the current solution g by using a predefined move mechanism [31, 32]. In this section, a multiple moves mechanism is used for neighborhood solution generation [33, 34]. That is, for each current solution, randomly selects an element in the matrix R and uses one type of moves (i.e., ζ-opt, ζ∈{1, 2, …, H-h}) to exchange rldh with rld[h+ζ] for a new neighborhood solution.
To improve the convergence of SA, this work opts out the random selection process. Instead, at each annealing step, the proposed SA identifies cores with the potential to improve their routing assignment and to minimize their tardiness. To that end, χ(ϑld) is used to evaluate the urgency of routing reassignment for cores. The probability of a core being selected for routing reassignment is defined as follows.
$ pro{b_{ld}} = \chi ({\vartheta _{ld}})/\sum\limits_l {\sum\limits_d {\chi ({\vartheta _{ld}})} } $ | (21) |
For any identified core, the selection of an alternative yet better routing is determined based on an efficiency index, $\Pi _{ld}^{h, x}$, which calculates the difference of the total average waiting time between the current process routing (i.e., the hth routing) and the alternative one (i.e., xth routing):
$\Pi _{ld}^{h,x} = {{\mathit{\Psi}} _h} - {{\mathit{\Psi}} _x},\forall \left\{ rldh=1x∈{1,2,⋯,H}−{h}rldx=0,aldx=1 \right.$
|
(22) |
${\Psi _h} = \sum\limits_{{t_{ij}} \in {R_h},{w_k} \in {}^ \bullet {t_{ij}}} {W{T_k}} {\Psi _x} = \sum\limits_{{t_{ij}} \in {R_x},{w_k} \in {}^ \bullet {t_{ij}}} {W{T_k}} $
|
(23) |
$W{T_k} = \frac{{\sum\limits_l {\sum\limits_d {(ST({w_k},{\vartheta _{ld}}) - AT({s_i},{\vartheta _{ld}}))} } }}{{{N_k}}}$ | (24) |
where WTk is the average waiting time per core consumed in the kth workstation, Ψh and Ψx are the total average waiting time per core consumed in the hth and xth process routing, respectively, and Nk is the total number of cores being processed by workstation wk. When the core ϑld is switched from the current routing to the other one with a larger $\Pi _{ld}^{h, x}$, it is more likely to relieve resource conflict and accomplish the task with a greatly reduced tardiness. Therefore, the probability of the xth routing being selected to replace the current hth routing for the core is defined as
$v_{ld}^{h,x} = \Pi _{ld}^{h,x}/\sum\limits_{x∈{1,2,⋯,H}−{h}aldx=1 } {\Pi _{ld}^{h,x}} $
|
(25) |
(4) Stopping criterion
In this proposed algorithm, two thresholds are pre-defined: the maximum number of inefficient iterations where the global optimum is not updated Imax, and the maximum number of outer loop iterations olmax. Our SA will stop whenever either of the thresholds is reached.
With a SA solution, each core follows the derived process routing to be recovered step by step. When multiple cores enter a workstation and compete for the same machine resource, a secondary optimization is conducted to determine the operation sequences of cores to be processed by the workstation, aiming to minimize the total tardiness penalty. The dispatching rules have been widely used for operation sequence optimization to address job shop scheduling problems [35]. Minimum Slack Time (MST), as one of typical heuristic sequencing rules, is widely used for the scheduling problems. MST measures the "urgency" of a job by its slack time, which is defined as the difference between its due date and production time. That is, the shorter the slack time, the higher priority to be processed. In the case of MST dispatching rules, the jobs might be processed as soon as possible in order to be finished right before its due date. As a result, the tardiness of all the jobs can be minimized. Recently, MST rule have been widely used for operation sequence optimization and are embedded with some meta-heuristic approaches to address various scheduling problems [36, 37]. In this section, MST is chosen for resource dispatching and is embedded with SA algorithm to search for the optimal schedules in minimizing both the operation cost and the penalty cost with related to tardiness.
For a specific core, the slack time of its recovery operation can be calculated as the temporal difference between the due date, the current time, and the remaining operation times, as formulated in Eq 26. MST sequences jobs in the descending order of slack time. Whenever a workstationwk releases a component ϑuv(∀Pri(wk, ϑuv) = q > 0) and turns back into idle state, it will pick up a core ϑld with the least slack time ($\exists $ slackld = min{slackxy}) from the waiting queue (∀ϑxy, AT(si, ϑxy)-Com(wk, ϑld) ≤ 0), and then immediately process it in the (q + 1)th sequence, as shown in Eq 27.
$slac{k_{xy}} = Due(l) - Com({w_k},{\vartheta _{uv}}) - \sum\limits_{pre({t_{cd}},{R_h}) \geqslant pre({t_{ij}},{R_h})} {u({t_{cd}},{\vartheta _{xy}})} \mathop ,\limits^{} \mathop {}\limits^{} \\\forall {r_{xyh}} = 1{,^{}}{t_{ij}} \in {R_h}{,^{}}{t_{cd}} \in {R_h}{,^{}}{t_{ij}} \in {}^ \bullet {w_k} \cap w_k^ \bullet {,^{}}Pri({w_k},{\vartheta _{uv}}) = q$ | (26) |
$Pri({w_k},{\vartheta _{ld}}) = q + 1,\;if\;\\\left\{ AT(si,ϑxy)−Com(wk,ϑuv)≤0,∀ruvh=1,tij∈Rh,si∈∙tij,tij∈∙wk∩w∙k,Pri(wk,ϑuv)=qslackld=min{slackxy},∀rxyh=1,tij∈Rh \right.$
|
(27) |
The solution obtained by MST is a 3-dimension matrix G = [gldk]L×D×K, where gldk, a non-negative integer, represents the operation sequence of the core ϑld to be performed by workstation wk. gldk = 0 means that the core ϑld is not processed by wk. Once the matrix G is derived, Pri(wk, ϑld) in CTPN are determined and thereby the total cost can be calculated.
With the recovery process routing R and operation sequence matrix G being derived through the proposed hybrid meta-heuristics, a set of transitions are then be enabled and fired as core tokens dynamically enter or release from the buffer places.
To fully understand the above concepts and algorithms, the remanufacturing of a batch of obsolete machine tools in an example system is used to demonstrate the remanufacturing scheduling. The workflows of cores in this remanufacturing shop and its corresponding CTPN are shown in Figs. 2 and 3. In order to verify the efficiency of the proposed method, three cases are considered for simulation. In the baseline case, cores are recovered through a set of fixed process routings. An available workstation takes MST rule for resource dispatching. In the second case, the traditional SA uses a standard neighborhood generation method and is embodied with MST rule for routing assignment and resource dispatching. In the third case, the recovery routing for a core and the resource dispatching for a workstation are guided through the proposed SA/MST by using the proposed neighborhood solution generation method.
All the algorithms are implemented in Matlab 2009. The simulation running duration is set to be 30 days. Everyday's work hours are 24 hours. Data of the first 3 days are thrown off as warm-up consideration. For the simulation, the input data is populated as follows.
(1) Since that the machine resources in the floor shop are fixed and limited, the machine workloads and the intension of resource conflicts increases with the rise of arrival rate, and vice versa. In order to simulated a set of resource conflict scenarios, the random arrival of used machine tools follows a Poisson distribution with an arrival rate per hour π $ \in ${5, 6, 7, 8, 9, 10, 11}. The designed arrival rate can limit the machine workloads (i.e., the average length of waiting queue on the machines) in the floor shop within a satisfactory range [0.75, 3.6].
(2) The due date (in days) for each used machine tool is randomly generated by a uniform distribution, Due(l) ~ U[5, 10].
(3) Each core disassembled from used machine tools is inspected before its recovery and is assigned a quality score Score(ϑld) through an Exponential distribution, Score(ϑld))~ Γ(1, τ) where τ ∈ [0.08, 0.10].
(4) The control factor and average operation time for each workstation is set to β(wk) = 0.2 and λ(wk) = 15.
(5) The operation cost per hour in workstations 7 and 9 (i.e., chromium electroplating and laser cladding) is $100, while that in other workstations is $50.
(6) The penalty cost per day per machine tool is set to be ϖ ∈ {20, 60, 100} in dollars.
(7) The parameters in SA are set toolmax = 5000, ilmax = 100, Imax = 30, Ω = 169340, η = 0.998.
(8) Considering that the total number (L) of machine tools generated in each run of the simulation varies, our comparison focuses on normalized objective values. In particular, the comparison evaluates the total cost per product (tc), operation cost per product (pc), tardiness penalty per product (dc), and the average waiting time per core consumed in each workstation (wt).
Our first set of experiments evaluates how the three methods respond to system demand and various penalty cost rates. In each trial, we run the three methods independently. It should be noted that the baseline case only employs MST rule for resource scheduling, while the proposed SA/MST and the standard one run continuously until they meet the stopping criteria. As the arrival rate of used machine tools changes from 5 to 11, the average total cost tc is obtained and compared among the three methods. This process is repeated with different penalty cost rates (i.e., ϖ = 20, 60, and 100) as shown in Tables 5, 6, and 7, respectively.
Arrival rate | The average total cost per machine tool | Improvement of the standard method over the baseline one tc1–tc2 | Improvement of the proposed method over the standard one tc2–tc3 | ||
Baseline casetc1 | Standard SA/MSTtc2 | Proposed SA/MSTtc3 | |||
π=5 | 186.66 | 117.02 | 117.02 | 69.64 | 0 |
π=6 | 228.37 | 133.48 | 118.21 | 94.89 | 15.27 |
π=7 | 293.05 | 172.45 | 146.89 | 120.6 | 25.56 |
π=8 | 335.87 | 208.06 | 168.69 | 127.81 | 39.37 |
π=9 | 447.52 | 266.47 | 216.25 | 181.05 | 50.22 |
π=10 | 581.63 | 352.28 | 270.91 | 229.35 | 81.37 |
π=11 | 714.29 | 461.86 | 335.77 | 252.43 | 126.09 |
Sample mean | 398.20 | 244.52 | 196.25 | 153.68 | 48.27 |
Sample standard deviation | 193.44 | 125.46 | 82.67 | 68.93 | 43.15 |
Arrival rate | The average total cost per machine tool | Improvement of the standard method over the baseline onetc1–tc2 | Improvement of the proposed method over the standard onetc2–tc3 | ||
Baseline casetc1 | Standard SA/MSTtc2 | Proposed SA/MSTtc3 | |||
π=5 | 215.76 | 117.83 | 117.83 | 97.93 | 0 |
π=6 | 327.61 | 160.31 | 140.12 | 167.3 | 20.19 |
π=7 | 455.9 | 226.47 | 185.02 | 229.43 | 41.45 |
π=8 | 658.09 | 342.57 | 273.22 | 315.52 | 69.35 |
π=9 | 891.23 | 515.96 | 406.38 | 375.27 | 109.58 |
π=10 | 1137.96 | 738.41 | 559.67 | 399.55 | 178.74 |
π=11 | 1530.54 | 986.65 | 775.72 | 543.89 | 210.93 |
Sample mean | 745.30 | 441.17 | 351.14 | 304.13 | 90.03 |
Sample standard deviation | 472.28 | 324.23 | 244.76 | 151.99 | 80.25 |
Arrival rate | The average total cost per machine tool | Improvement of the standard method over the baseline onetc1–tc2 | Improvement of the proposed method over the standard onetc2–tc3 | ||
Baseline casetc1 | Standard SA/MSTtc2 | Proposed SA/MSTtc3 | |||
π=5 | 270.84 | 117.32 | 117.32 | 153.52 | 0 |
π=6 | 463.08 | 224.8 | 186.58 | 238.28 | 38.22 |
π=7 | 687.42 | 416.46 | 275.13 | 270.96 | 141.33 |
π=8 | 1063.8 | 655.68 | 453.47 | 408.12 | 202.21 |
π=9 | 1495.67 | 957.36 | 706.49 | 538.31 | 250.87 |
π=10 | 1917.52 | 1231.56 | 918.76 | 685.96 | 312.8 |
π=11 | 2403.39 | 1483.89 | 1120.77 | 919.5 | 363.12 |
Sample mean | 1185.96 | 726.72 | 539.79 | 459.24 | 186.94 |
Sample standard deviation | 789.17 | 517.36 | 385.07 | 273.81 | 135.61 |
The performance of the standard SA/MST and the baseline one is first evaluated and compared by calculating the cost difference under different arrival rates, as shown in Figure 5. It is clear that the proposed SA/MST outperforms the baseline method in terms of generating less cost, particularly in reacting to the larger arrival rates. This is reasonable since that when a substantial amount of cores are fighting for limited resources for to be recovered, the standard SA/MST operates on routings reassignment of cores to reduce the ever intensified resource conflicts and thus minimize a relatively large amount of tardiness penalty to lower the total cost, while the baseline method is not capable of addressing so by fixing process routings of cores. The results in Figure 5 can be also interpreted as evidence that the hybrid SA/MST is promising for remanufacturing scheduling optimization.
The differences of total cost per product obtained by the two hybrid methods are also calculated and plotted in Figure 6. The following remarks can be made regarding the performance of the proposed SA/MST over the standard one.
(1) When the arrival rate is low (e.g., π = 5), the system has sufficient enough resources to handle remanufacturing demand. Thus the benefit of our proposed SA/MST in comparison to the standard one is minor. As the system demand increases with the rising arrival rate, more and more used machine tools are fighting for the limited resource to be remanufactured. The proposed method then outperforms the standard method through a more targeted annealing process to properly reassign recovery routings of cores and resolve ever intensified resource conflict. Therefore, the advantage of the proposed method over the standard one becomes significant.
(2) In a similar fashion, by comparing the corresponding row in Tables 5, 6 and 7, the relative improvement of the proposed method over the standard one becomes more and more significant with the increasing penalty cost. Using π = 8 as an example, the improvement increases from 11% when ϖ = 20, to 17% when ϖ = 60, and to 30% when ϖ = 100.
To further show the significant difference between the performances of the two hybrid algorithms, a sample mean, a sample standard deviation and confident interval estimation are introduced. A relative value σ = –(tc3–tc2)/tc2 is defined to characterize the improvement of the best solution obtained by the proposed SA/MST over the one obtained by the standard one. A 95% confidence interval of the improvement rate σ is conducted using the data in Tables 5–7. The approximate 100(1–α)% confidence interval for σ is defined as
$ ⌢σ±t1−α/2(n−1)SD√nor ⌢σ−t1−α/2(n−1)SD√n≤σ≤⌢σ+t1−α/2(n−1)SD√n $
|
(28) |
where ${\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{\sigma }}$ is the sample mean of σ based on a sample of size n; SD is the sample standard deviation; ${{t}_{{1-\alpha }/{2}\; }}(n-1)$is the 100(1–α)% percentage point of a t-distributed with n–1 degree of freedom. The 95% confidence intervals for the improvement σ under three scenarios (i.e., ϖ = 20, 60 and 100) are given in Eqs 29, 30 and 31, respectively. It is obvious that the 95% confidence interval for σ lies completely above zero (i.e., $\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{\sigma } > 0$), which provides strong evidence that the proposed SA/MST is better than the standard one in terms of generating less total cost.
$0.0814 \leqslant {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{\sigma } _1} \leqslant 0.2455$ | (29) |
$0.1042 \leqslant {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{\sigma } _2} \leqslant 0.3118$ | (30) |
$0.1414 \leqslant {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{\sigma } _3} \leqslant 0.3911$ | (31) |
The results in Tables 5–7 represents significant enhancement of the proposed methodology in terms of overall total cost per product, especially when the system is over utilized and its responsiveness to fluctuating market demand is important.
In the second set of experiments, the performance of the two hybrid algorithms is further evaluated by calculating the normalized objective values achieved after pre-specified time breaks (Table 8), where the arrival rate and the penalty cost are set as π = 8 and ϖ = 60. The following remarks can be made regarding algorithms' performance.
Computation time breaks | Standard SA/MST | Proposed SA/MST | ||||||
tc1 | pc1 | dc1 | wt1 | tc2 | pc2 | dc2 | wt2 | |
After 0sec | 399.25 | 112.16 | 287.09 | 76.80 | 399.25 | 112.16 | 287.09 | 76.80 |
After 30sec | 380.38 | 116.78 | 263.60 | 72.37 | 371.14 | 114.62 | 256.51 | 55.21 |
After 60sec | 393.70 | 113.71 | 279.99 | 61.76 | 341.42 | 118.56 | 222.86 | 48.34 |
After 90sec | 393.70 | 113.71 | 279.99 | 58.90 | 311.93 | 120.51 | 191.43 | 46.49 |
After 120sec | 386.11 | 115.62 | 270.49 | 58.75 | 292.79 | 121.75 | 171.04 | 44.54 |
After 150sec | 385.56 | 119.96 | 265.60 | 58.42 | 280.06 | 123.28 | 156.78 | 41.52 |
After 180sec | 384.17 | 121.18 | 262.99 | 57.90 | 274.14 | 124.64 | 149.51 | 41.47 |
After 210sec | 383.35 | 122.86 | 260.49 | 56.51 | 273.33 | 126.31 | 147.03 | 38.40 |
After 240sec | 378.22 | 124.65 | 252.57 | 54.73 | 273.22 | 132.79 | 140.42 | 35.67 |
No time limits* | 342.57 | 129.41 | 213.16 | 48.07 | 273.22 | 132.79 | 140.42 | 35.67 |
* "No time limit" means the simulation keeps the algorithm running until it meets the stopping criteria elaborated above. |
(1) As the computation time increases, both algorithms strive to search for a better scheduling solution that maintains a good trade-off between operation cost and tardiness penalty. It should be noted that both algorithms start at the same initial solution (i.e., zero time break). Since the initial solution is generated with the goal of minimizing operation cost only, it presents a relatively larger penalty cost. The annealing step in both methods then begins to rearrange routings for cores that ultimately balances the overall workload of workstations and thus reduces waiting time to lower the tardiness penalty. As shown in Table 8, at the 210-second time break, the optimal solutions generated by both methods have a bit larger pc and a quite smaller dc, resulting in a much lower total cost per product.
(2) The proposed SA/MST method converges much rapidly than the standard one, providing a better solution within a given short time frame. This is intuitively understandable since the proposed SA/MST uses resource conflict and tardiness penalty to guide neighborhood solution generation while the random mechanism in the standard SA fails to address so, resulting in a relatively larger tardiness penalty.
All simulation output analysis verifies the effectiveness of the proposed hybrid SA/MST algorithm for the remanufacturing scheduling.
Remanufacturing systems are faced with a greater degree of uncertainty and complexity than traditional manufacturing systems, leading to the need for planning and control systems designed to deal with the added uncertainty and complexity. Although many new methodologies have been led up to deal with various operation management issues in remanufacturing environments, no study has comprehensively dealt with remanufacturing job shop scheduling in subjecting to alternative recovery operations, random operation times and limited resource conflicts. A CTPN is introduced to model the dynamics of remanufacturing process, such as various process routings, uncertain operations times and resource conflicts. With time and color attributes in PN, two types of decision variables are linked with places in CTPN and the evolution of system dynamics in recovery operations are mathematically analyzed. With the support of the computation model via CTPN, a hybrid meta-heuristic using SA and MST rule is proposed and embedded with CTPN to sample large search space efficiently and search for a good scheduling solution in minimizing total production cost. The performance of the proposed SA/MST algorithm is compared against another two cases: baseline case with fixed recovery process routings and case 2 using standard SA/MST. The comparison results provide strong evidence that the proposed scheduling method is of significant importance in achieving minimum production cost especially when the system is overloaded and its responsiveness to shared resource conflicts is important.
Our research can be also extended in several directions. For instance, the integration of alternative recovery routing selection and resource dispatching is a NP-hard problem. Therefore, it is worth of investigation to take the structure advantage of CTPN for a more efficient optimization heuristics. Another challenge is to take into consideration of unexpected system disruptions (i.e., random machine breakdown) in remanufacturing scheduling. The future work also includes further validation of our methodology using more factory data.
This work was supported in part by the Fundamental Research Funds for the Central Universities (SWU22300503), the National Natural Science Foundation of China (51875480) and the National Natural Science Foundation of China (51475059).
All authors declare that there is no conflict of interests in this paper.
[1] | Keshk SMAS (2014) Bacterial cellulose production and its industrial applications. J Bioproces Biotechniq 4: 2. |
[2] |
Rajwade JM, Paknikar KM, Kumbhar JV (2015) Applications of bacterial cellulose and its composites in biomedicine. Appl Microbiol Biotechnol 99: 2491-2511. doi: 10.1007/s00253-015-6426-3
![]() |
[3] | Araujo S, Moreira da Silva F, Gouveia IC (2015) The role of technology towards a new bacterial-cellulose-based material for fashion design. J Ind Intell Inf 3: 168-172. |
[4] |
Castellane TCL, Persona MR, Campanharo JC, et al. (2015) Production of exopolysaccharide from rhizobia with potential biotechnological and bioremediation applications. Int J Biol Macromol 74: 515-522. doi: 10.1016/j.ijbiomac.2015.01.007
![]() |
[5] |
Ross P, Weinhouse H, Aloni Y, et al. (1987) Regulation of cellulose synthesis in Acetobacter xylinum by cyclic diguanylic acid. Nature 325: 279-281. doi: 10.1038/325279a0
![]() |
[6] | Ude S, Arnolod DL, Moon CD, et al. (2006) Biofilm formation and cellulose expression among diverse environmental Pseudomonas isolates. Environ Microbiol 8: 1997-2011. |
[7] |
Zogaj X, Nimtz M, Rohde M, et al. (2001) The multicellular morphotypes of Salmonella typhimurium and Escherichia coli produce cellulose as the second component of the extracellular matrix. Mol Microbiol 39: 1452-1463. doi: 10.1046/j.1365-2958.2001.02337.x
![]() |
[8] | Romling U (2002) Molecular biology of cellulose production in bacteria. ResMicrobiol 153: 205-212. |
[9] |
Solano C, Garcia B, Valle J, et al. (2002) Genetic analysis of Salmonella enteritidis biofilm formation: critical role of cellulose. Mol Microbiol 43: 793-808. doi: 10.1046/j.1365-2958.2002.02802.x
![]() |
[10] | Mohite BV, Patil SV (2014) Physical, structural, mechanical and thermal characterization of bacterial cellulose by G. hansenii NCIM 2529. Carbohyd Polym106: 132-141. |
[11] | Morgan JL, Strumillo J, Zimmer J (2013) Crystallographic snapshot of cellulose synthesis and membrane translocation. Nature 493: 181-186. |
[12] | Amikam D, Benziman M (1989) Cyclic diguanylic acid and cellulose synthesis in Agrobacterium tumefaciens. J Bacteriol 171: 6649-6655. |
[13] | Romling U (2012) Cyclic di-GMP, an established secondary messenger still speeding up. Environ Microbiol 14: 1817-1829. |
[14] | Whiteley CG, Lee DJ (2014) Bacterial diguanylate cyclases: Structure, function and mechanism in exopolysaccharide biofilm development. Biotech Adv 33: 124-41. |
[15] | Morgan JL, Mcnamara JT, Zimmer J (2014) Mechanism of activation of bacterial cellulose synthase by cyclic di-GMP. Nature Struc Biol 21: 489-496. |
[16] | Dazzo FB, Truchet GL, Sherwood JE, et al. (1984) Specific phases of root hair attachment in the Rhizobium trifolii-clover symbiosis. Appl Environ Microbiol 48: 1140-1150. |
[17] |
Mateos PF, Baker DL, Philip-Hollingsworth S, et al. (1995) Direct in situ identification of cellulose microfibrils associated with Rhizobium leguminosarum biovar trifolii attached to the root epidermis of white clover. Can J Microb 41: 202-207. doi: 10.1139/m95-028
![]() |
[18] | Flemming HC, Wingender J (2010) The biofilm matrix. Nat Rev 8: 624-633. |
[19] |
McCrate OA, Zhou X, Reichhardt C, et al. (2013) Sum of the Parts: Composition and Architecture of the Bacterial Extracellular Matrix. J Mol Biol 425: 4286-4294. doi: 10.1016/j.jmb.2013.06.022
![]() |
[20] |
Koeck DE, Pechtl A, Zerlov VV, et al., (2014) Genomics of cellulolytic bacteria Curr Opin Biotech 29: 171-183. doi: 10.1016/j.copbio.2014.07.002
![]() |
[21] |
Lombard V, Ramulu HG, Drula E, et al. (2014) The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res 42: D490-D495. doi: 10.1093/nar/gkt1178
![]() |
[22] |
Davies G, Henrissat B (1995) Structures and mechanisms of glycosyl hydrolases. Structure 3: 853-859. doi: 10.1016/S0969-2126(01)00220-9
![]() |
[23] |
Bourne Y, Henrissat B (2001) Glycoside hydrolases and glycosyltransferases: families and functional modules. Curr Opin Struct Biol 11: 593-600. doi: 10.1016/S0959-440X(00)00253-0
![]() |
[24] |
Sadhu S, Maiti TK (2013) Cellulase production by bacteria: a review. British Microbiology Research Journal 3: 235-258. doi: 10.9734/BMRJ/2013/2367
![]() |
[25] | Fuji K, Satomi M, Fukui Y, et al. (2013) Streptomyces abietis sp. nov., a cellulolytic bacterium isolated from soil of a pine forest. Int J Syst Evol Microbiol 63: 4754-4759. |
[26] | Cole JK, Gieler BA, Heisler DL, et al. (2013) Kallotenue papyrolyticum gen. nov., sp. nov., a cellulolytic and filamentous thermophile that represents a novel lineage (Kallotenuales ord. nov., Kallotenuaceae fam. nov.) within the class Chloroflexia. Int J Syst Evol Microbiol 63: 4675-4682. |
[27] | Podosokorskaya OA, Bonch-Osmolovskaya EA, Novikov AA, et al. (2013) Ornatilinea apprima gen. nov., sp. nov., a cellulolytic representative of the class Anaerolineae. Int J Syst Evol Microbiol 63: 86-92. |
[28] | Hatamoto M, Kaneshige M, Nakamura A, et al. (2014) Bacteroides luti sp. nov., an anaerobic, cellulolytic and xylanolytic bacterium isolated from methanogenic sludge. Int J Syst Evol Microbiol 64: 1770-1774. |
[29] | Kusube M, Sugihara A, Moriwaki Y, et al. (2014) Alicyclobacillus cellulosilyticus sp. nov., a thermophilic, cellulolytic bacterium isolated from steamed Japanese cedar chips from a lumbermill. Int J Syst Evol Microbiol 64: 2257-2263. |
[30] | Horino H, Fujita T, Tonouchi A (2014) Description of Anaerobacterium chartisolvens gen. nov., sp. nov., an obligately anaerobic bacterium from Clostridium rRNA cluster III isolated from soil of a Japanese rice field, and reclassification of Bacteroides cellulosolvens Murray et al. 1984 as Pseudobacteroides cellulosolvens gen. nov., comb. nov. Int J Syst Evol Microbiol 64: 1296-1303. |
[31] | Bing W, Wang H, Zheng B, et al. (2015) Caldicellulosiruptor changbaiensis sp. nov., a cellulolytic and hydrogen-producing bacterium from a hot spring. Int J Syst Evol Microbiol 65: 293-297. |
[32] | Koeck DE, Ludwig W, Wanner G, et al. (2015) Herbinix hemicellulosilytica, gen. nov., sp. nov., a thermophilic cellulose-degrading bacterium isolated from a thermophilic biogas reactor. Int J Syst Evol Microbiol [in press]. |
[33] | Menéndez E, Ramírez-Bahena MH, Fabryová A, et al. (2015) Pseudomonas coleopterorum sp. nov., a cellulase producing bacterium isolated from the bark beetle Hylesinus fraxini. Int J Syst Evol Microbiol [in press]. |
[34] | Huang S, Sheng P, Zhang H (2012) Isolation and Identification of Cellulolytic Bacteria from the Gut of Holotrichia parallela Larvae (Coleoptera: Scarabaeidae) Int J Mol Sci 13: 2563-2577. |
[35] | Hameed A, Shahina M, Lai WA, et al. (2015) Oricola cellulosilytica gen. nov., sp. nov., a cellulose-degrading bacterium of the family Phyllobacteriaceae isolated from surface seashore water, and emended descriptions of Mesorhizobium loti and Phyllobacterium myrsinacearum. A van Leeuw J Microb 107: 759-771. |
[36] | Sethi S, Datta A, Gupta L, et al. (2013) Optimization of Cellulase Production from Bacteria Isolated from Soil. ISRN Biotechnology 2013: 985685. |
[37] | Howard JA, Hood E (2005) Bioindustrial and biopharmaceutical products produced in plants. Adv Agron 85: 91-124. |
[38] | Yu LX, Gray BN, Rutzke CJ, et al. (2007) Expression of thermostable microbial cellulases in the chloroplasts of nicotine-free tobacco. J Biotechnol1 31: 362-369. |
[39] |
Rastogi G, Muppidi GL, Gurram RN, et al. (2009) Isolation and characterization of cellulose-degrading bacteria from the deep subsurface of the Homestake gold mine, Lead, South Dakota, USA. J Ind Microbiol Biotechnol 36: 585-598. doi: 10.1007/s10295-009-0528-9
![]() |
[40] |
Trivedi N, Gupta V, Kumar M, et al. (2011) An alkali-halotolerant cellulase from Bacillus flexus isolated from green seaweed Ulva lactuca. Carbohyd Polym 83: 891-897. doi: 10.1016/j.carbpol.2010.08.069
![]() |
[41] |
DiPasquale L, Romano I, Picardello G, et al., (2014) Characterization of a native cellulase activity from an anaerobic thermophilic hydrogen-producing bacterium Thermosipho sp. strain 3. Ann Microbiol 64: 1493-1503. doi: 10.1007/s13213-013-0792-9
![]() |
[42] |
Dalmaso GZL, Ferreira D, Vermelho AB (2015) Marine extremophiles: a source of hydrolases for biotechnological applications. Mar Drugs 13: 1925-1965. doi: 10.3390/md13041925
![]() |
[43] |
Kinet R, Destain J, Hiligsmann S, et al. (2015) Thermophilic and cellulolytic consortium isolated from composting plants improves anaerobic digestion of cellulosic biomass: Toward a microbial resource management approach. Bioresource Technol 189: 138-144. doi: 10.1016/j.biortech.2015.04.010
![]() |
[44] |
Voget S, Steele HL, Streit WR (2006) Characterization of a metagenome-derived halotolerant cellulase. J Biotechnol 126: 26-36. doi: 10.1016/j.jbiotec.2006.02.011
![]() |
[45] | Duan CJ, Feng JX (2010) Mining metagenomes for novel cellulase genes. Biotechnol Lett 32: 1765-1775. |
[46] |
Liu J, Liu W, Zhao X, et al. (2011) Cloning and functional characterization of a novel endo-β-1,4-glucanase gene from a soil-derived metagenomic library. Appl Microbiol Biotechnol 89: 1083-1092. doi: 10.1007/s00253-010-2828-4
![]() |
[47] | Adrio JL, Demain AL (2014) Microbial enzymes: tools for biotechnological processes. Biomolecules 4: 117-139. |
[48] | Mori T, Kamei I, Hirofumi H, et al. (2014) Identification of novel glycosyl hydrolases with cellulolytic activity against crystalline cellulose from metagenomic libraries constructed from bacterial enrichment cultures. Springer Plus 3: 365. |
[49] | Attri S, Garg G (2014) Isolation of microorganisms simultaneously producing xylanase, pectinase and cellulase enzymes using cost effective substrates. J Innovative Biol 1: 45-50. |
[50] | Singh G, Singh AK (2014) Alternative substrates for the amylase and cellulase production with rhizobial isolates. Int J Avd Res Sci Technol 3: 79-85. |
[51] | Teather RM, Wood PJ (1982) Use of Congo red-polysaccharide interactions in enumeration and characterization of cellulolytic bacteria from the bovine rumen. Environ Microbiol 43: 777-780. |
[52] |
Kasana RC, Salwan R, Dhar H, et al. (2008) A rapid and easy method for the detection of microbial cellulases on agar plates using gram’s iodine. Curr Microbiol 57: 503-507. doi: 10.1007/s00284-008-9276-8
![]() |
[53] | Taha M, Kadali KK, Al-Hothaly K, et al. (2015) An effective microplate method (Biolog MT2) for screening native lignocellulosic-straw-degrading bacteria. Ann Microbiol [in press]. |
[54] |
Zhang YHP, Himmel ME, Mielenz JR (2006) Outlook of cellulase improvement: screening and selection strategies. Biotechnol Adv 24: 452-481. doi: 10.1016/j.biotechadv.2006.03.003
![]() |
[55] | Sadhu S, Ghosh PK, Aditya G, et al. (2014) Optimization and strain improvement by mutation for enhanced cellulase production by Bacillus sp. (MTCC10046) isolated from cow dung. J King Saud Univ 26: 323-332. |
[56] | Sangkharak K, Vangsirikul P, Janthachat S (2012) Strain improvement and optimization for enhanced production of cellulase in Cellulomonas sp. TSU-03. African J Microbiol Res 6: 1079-1084. |
[57] | Maki M, Leung KT, Qin W (2009) The prospects of cellulose-producing bacteria for the bioconversion of lignocellulosic biomass. Int J Biol Sci 5: 500-516. |
[58] |
Zhao H, Chockalingam K, Chen Z (2002) Directed evolution of enzymes and pathways for industrial biocatalysis. Curr Opin Biotechnol 13: 104-110. doi: 10.1016/S0958-1669(02)00291-4
![]() |
[59] | Cherry JR, Fidantsef AL (2003) Directed evolution of industrial enzymes: an update. Curr Opin Biotechnol 14: 438-443. |
[60] |
Hasunuma T, Okazaki F, Okai N, et al. (2013) A review of enzymes and microbes for lignocellulosic biorefinery and the possibility of their application to consolidated bioprocessing technology. Bioresource Technol 135: 513-522. doi: 10.1016/j.biortech.2012.10.047
![]() |
[61] | Garvey M, Klose H, Fischer R, et al. (2014) Cellulases for biomass degradation: comparing recombinant cellulase expression platforms. Trends Biotechnol 31: 581-593. |
[62] |
Juturu V, Wu JC (2014) Microbial cellulases: Engineering, production and applications. Renew Sust Energ Rev 33: 188-203. doi: 10.1016/j.rser.2014.01.077
![]() |
[63] |
Lambertz C, Garvey M, Klinger J, et al. (2014) Challenges and advances in the heterologous expression of cellulolytic enzymes: a review. Biotechnol Biofuels 7: 135. doi: 10.1186/s13068-014-0135-5
![]() |
[64] |
Munjal N, Jaewd K, Wajid S, et al. (2015) A Constitutive Expression System for Cellulase Secretion in Escherichia coli and Its Use in Bioethanol Production. PLoS ONE 10: e0119917. doi: 10.1371/journal.pone.0119917
![]() |
[65] | Chung D, Cha M, Guss AM, et al. (2015) Direct conversion of plant biomass to ethanol by engineered Caldicellulosiruptor bescii. P Natl Acad Sci U S A 111: 8931-8936. |
[66] |
Robledo M, Jimenez-Zurdo JI, Soto MJ, et al. (2011) Development of functional symbiotic white clover root hairs and nodules requires tightly regulated production of rhizobial cellulase CelC2. Mol Plant Microbe Interact 24: 798-807. doi: 10.1094/MPMI-10-10-0249
![]() |
[67] | Batthacharya AS, Batthacharya A, Pletschke BI (2015) Synergism of fungal and bacterial cellulases and hemicellulases: a novel perspective for enhanced bio-ethanol production. Biotechnol Lett [in press]. |
[68] |
Demain AL, Vaishnav P (2009) Production of recombinant proteins by microbes and higher organisms. Biotechnol Adv 27: 297-306. doi: 10.1016/j.biotechadv.2009.01.008
![]() |
[69] |
Streatfield AJ (2007) Approaches to achieve high-level heterologous protein production in plants. Plant Biotechnol J 5: 2-15. doi: 10.1111/j.1467-7652.2006.00216.x
![]() |
[70] | Ziegelhoffer T, Will J, Austin-Phillips S (1998) Expression of bacterial cellulase genes in transgenic alfalfa (Medicago sativa L.), potato (Solanum tuberosum L.) and tobacco (Nicotiana tabacum L.) Mol Breeding 5: 309-318. |
[71] | Brunecki R, Selig MJ, Vinzant TB, et al. (2011) In planta expression of A. cellulolyticus Cel5A endocellulase reduces cell wall recalcitrance in tobacco and maize. Biotechnol Biofuels 4: 1. |
[72] |
Petersen K, Bock R (2011) High-level expression of a suite of thermostable cell wall-degrading enzymes from the chloroplast genome. Plant Mol Biol 76: 311-321. doi: 10.1007/s11103-011-9742-8
![]() |
[73] |
Lilly M, Fierobe HP, van Zyl WH, et al. (2009) Heterologous expression of a Clostridium minicellulosome in Saccharomyces cerevisiae. FEMS Yeast Res 9: 1236-1249. doi: 10.1111/j.1567-1364.2009.00564.x
![]() |
[74] |
Bayer EA, Shimon LJW, Lamed R (1998) Cellulosomes: structure and ultrastructure. J Struct Biol 124: 221-234. doi: 10.1006/jsbi.1998.4065
![]() |
[75] |
Bayer EA, Belaich JP, Sholam Y, et al. (2004) The cellulosomes: multienzyme machines for degradation of plant cell wall polysaccharides. Ann Rev Microbiol 58: 521-554. doi: 10.1146/annurev.micro.57.030502.091022
![]() |
[76] |
Fontes CMGA, Gilbert HJ (2010) Cellulosomes: highly efficient nanomachines designed to deconstruct plant cell wall complex carbohydrates. Ann Rev Biochem 79: 655-681. doi: 10.1146/annurev-biochem-091208-085603
![]() |
[77] |
Shinoda S, Kanamasa S, Arai M (2012) Cloning of an endoglycanase gene from Paenibacillus cookii and characterization of the recombinant enzyme. Biotechnol Lett 34: 281-286. doi: 10.1007/s10529-011-0759-5
![]() |
[78] |
Shi R, Li Z, Ye Q, et al. (2013) Heterologous expression and characterization of a novel thermo-halotolerant endoglucanase Cel5H from Dictyoglomus thermophilum. Bioresource Technol 142: 338-344. doi: 10.1016/j.biortech.2013.05.037
![]() |
[79] |
Wei KSC, Teoh TC, Koshy P, et al. (2015) Cloning, expression and characterization of the endoglucanase gene from Bacillus subtilis UMC7 isolated from the gut of the indigenous termite Macrotermes malaccensis in Escherichia coli. Electronic J Biotechnol 18: 103-109. doi: 10.1016/j.ejbt.2014.12.007
![]() |
[80] |
Bhat MK (2000) Cellulases and related enzymes in biotechnology. Biotechnology Adv 18: 355-383. doi: 10.1016/S0734-9750(00)00041-0
![]() |
[81] | Wilson DB (2009) Cellulases and biofuels. Curr Opin Biotechnol20: 295-299. |
[82] | Kuhad RC, Gupta R, Singh A. (2011) Microbial cellulases and their industrial applications. Enzyme Res 2011: 280696. |
[83] |
McMullan G, Meehan C, Connely M, et al. (2001) Microbial decolourisation and degradation of textile dyes. Appl Microbiol Biotechnol 56: 81-87. doi: 10.1007/s002530000587
![]() |
[84] |
Anish R, Rahman MS, Rao M (2007) Application of cellulases from an alkalothermophilic Thermomonospora sp. in biopolishing of denims. Biotechnol Bioeng 96: 48-56. doi: 10.1002/bit.21175
![]() |
[85] | Ladeira SA, Cruz E, Delatorre AB, et al. (2015) Cellulase production by thermophilic Bacillus sp. SMIA-2 and its detergent compatibility. Electronic J Biotechnol 18: 110-115. |
[86] | Yu M, Qiu Y, Chen W, et al., (2015) Action modes of recombinant endocellulase, EGA, and its domains on cotton fabrics. Biotechnol Lett [in press]. |
[87] | Singh K (2015) Role of Enzymes in Fruit juices Clarification during Processing: A review. Int J Biol Technology 6: 114-124. |
[88] | Cinar I (2004) Effects of cellulase and pectinase concentrations on the colour yield of enzyme extracted plant carotenoids. Process Biochem 40: 945-949. |
[89] |
Wilkins MR, Widmer WW, Grohmann K, et al. (2007) Hydrolysis of grapefruit peel waste with cellulase and pectinase enzymes. Bioresource Technol 98: 1596-1601. doi: 10.1016/j.biortech.2006.06.022
![]() |
[90] | Meyer AS, Jepsen SM, Sorensen NS (1998) Enzymatic Release of Antioxidants for Human Low-Density Lipoprotein from Grape Pomace. J Agr Food Chem 46: 2399-2446. |
[91] |
Bamforth CW (2009) Current perspectives on the role of enzymes in brewing. J Cereal Sci 50: 353-357 doi: 10.1016/j.jcs.2009.03.001
![]() |
[92] |
Himmel ME, Ruth MF, Wyman CE (1999) Cellulase for commodity products from cellulosic biomass. Curr Opin Biotechnol 10: 358-364. doi: 10.1016/S0958-1669(99)80065-2
![]() |
[93] |
Dhiman TR, Zaman MS, MacQueen IS, et al. (2002) Influence of corn processing and frequency of feeding on cow performance. J Dairy Sci 85: 217-226. doi: 10.3168/jds.S0022-0302(02)74070-8
![]() |
[94] | Beauchemin KA, Colombatto D, Morgavi DP, et al. (2003) Use of exogenous fibrolytic enzymes to improve feed utilization by ruminants. J Anim Sci 81: E37-E47. |
[95] |
Wongputtissin P, Khanongnuch C, Kongbuntad W, et al. (2014) Use of Bacillus subtilis isolates from Tua-nao towards nutritional improvement of soya bean hull for monogastric feed application. Lett App Microbiol 59: 328-333. doi: 10.1111/lam.12279
![]() |
[96] | Shatya TA, Khan M (2014) Diversity of glycosyl hydrolase enzymes from metagenome and their application in food industry. J Food Sci 79: R2149-R 2156. |
[97] |
Pottkamper J, Barthen P, Ilmberger N, et al. (2009) Applying metagenomics for the identification of bacterial cellulases that are stable in ionic liquids. Green Chem 11: 957-965. doi: 10.1039/b820157a
![]() |
[98] | Eriksson KEL (1990) Biotechnology in the pulp and paper industry.Wood Sci Technol 24: 79-101. |
[99] |
Viesturs U, Leite M, Eisimonte M, et al. (1999) Biological deinking, technology for the recycling of office waste papers. Bioresource Technol 67: 255-265. doi: 10.1016/S0960-8524(98)00119-9
![]() |
[100] |
Bajpai P (1999) Application of Enzymes in the Pulp and Paper Industry. Biotechnol Prog 15: 147-157. doi: 10.1021/bp990013k
![]() |
[101] |
Garcia O, Torres AL, Colom JF, et al. (2002) Effect of cellulase-assisted refining on the properties of dried and never-dried eucalyptus pulp. Cellulose 9: 115-125. doi: 10.1023/A:1020191622764
![]() |
[102] |
Lynd LR, Laser MS, Bransby D, et al. (2008) How biotech can transform biofuels. Nat Biotechnol 26: 169-172. doi: 10.1038/nbt0208-169
![]() |
[103] |
Yang SJ, Kataeva I, Wiegel J, et al. (2010) Classification of ‘Anaerocellum thermophilum’ strain DSM 6725 as Caldicellulosiruptor bescii sp. nov. Int J Syst Evol Microbiol 60: 2011-2015 doi: 10.1099/ijs.0.017731-0
![]() |
[104] |
Kanafusa-Shinkai S, Wakayama J, Tsukamoto K, et al. (2013) Degradation of microcrystalline cellulose and non-pretreated plant biomass by a cell-free extracellular cellulase/hemicellulase system from the extreme thermophilic bacterium Caldicellulosiruptor bescii. J Biosci Bioeng 115: 64-70. doi: 10.1016/j.jbiosc.2012.07.019
![]() |
[105] | Scully SM, Orlygsson J (2015) Recent advances in second generation ethanol production by thermophilic bacteria. Energies 8: 1-30. |
[106] | Assareh R, Zahiri HS, Noghabi KA, et al. (2012) Characterization of the newly isolated Geobacillus sp. T1, the efficient cellulose-producer on untreated barley and wheat straws. Bioresource Technol 120: 99-105. |
[107] |
Gaur R, Tiwati S (2015) Isolation, production, purification and characterization of an organic-solvent-thermostable alkalophilic cellulase from Bacillus vallismortis RG-07. BMC Biotechnol 15: 19. doi: 10.1186/s12896-015-0129-9
![]() |
[108] | Li Y, Horsman M, Wu N, et al. (2008) Biofuels from microalgae. Biotechnol Prog 24: 815-820. |
[109] |
Muñoz C, Hidalgo C, Zapata M, et al. (2014) Use of cellulolytic marine bacteria for enzymatic pretreatment in microalgal biogas production. App Environ Microbiol 80: 4199-4206. doi: 10.1128/AEM.00827-14
![]() |
[110] |
Rinaudo M (2006) Chitin and chitosan: Properties and applications. Prog Polym Sci 31: 603-632. doi: 10.1016/j.progpolymsci.2006.06.001
![]() |
[111] |
Pillai CKS, Paul W, Sharma CP (2009) Chitin and chitosan polymers: Chemistry, solubility and fiber formation. Prog Polym Sci 34: 641-678. doi: 10.1016/j.progpolymsci.2009.04.001
![]() |
[112] |
Zhang J, Xia W, Liu P, et al. (2010) Chitosan Modification and Pharmaceutical/Biomedical Applications. Mar Drugs 8: 1962-1987. doi: 10.3390/md8071962
![]() |
[113] |
Qin C, Zhou B, Zeng L, et al. (2004) The physicochemical properties and antitumor activity of cellulase-treated chitosan. Food Chem 84: 107-115. doi: 10.1016/S0308-8146(03)00181-X
![]() |
[114] |
Lin SB, Lin YC, Chen HH (2009) Low molecular weight chitosan prepared with the aid of cellulase, lysozyme and chitinase: Characterisation and antibacterial activity. Food Chem 116: 47-53. doi: 10.1016/j.foodchem.2009.02.002
![]() |
[115] |
Liu J, Xia W (2006) Purification and characterization of a bifunctional enzyme with chitosanase and cellulase activity from commercial cellulase. Biochem Eng J 30: 82-87. doi: 10.1016/j.bej.2006.02.005
![]() |
[116] |
Xie Y, Wei Y, Hu J (2010) Depolymerization of Chitosan with a Crude Cellulase Preparation from Aspergillus niger. Appl Biochem Biotechnol 160: 1074-1083. doi: 10.1007/s12010-009-8559-2
![]() |
[117] | Pedraza-Reyes M, Gutierrez-Corona F (1997) The bifunctional enzyme chitosanase-cellulase produced by the gram-negative microorganism Myxobacter sp. AL-1 is highly similar to Bacillus subtilis endoglucanases. Arch Microbiol 168: 321-327. |
[118] |
Tanabe T, Morinaga K, Fukamizo T, et al. (2003) Novel Chitosanase from Streptomyces griseus HUT 6037 with Transglycosylation Activity. Biosci Biotechnol Biochem 67: 354-364. doi: 10.1271/bbb.67.354
![]() |
[119] | Sinha S, Tripathi P, Chand S (2012) A New Bifunctional Chitosanase Enzyme from Streptomyces sp. and Its Application in Production of Antioxidant Chitooligosaccharides. Appl Biochem Biotechnol 167: 1029-1039. |
[120] |
Zhang Z, Wang Y, Ruan J (1998) Reclassification of Thermomonospora and Microtetraspora. J Syst Bacteriol 48: 411-422. doi: 10.1099/00207713-48-2-411
![]() |
[121] | Lao G, Changas GS, Jung ED, et al. (1991) DNA Sequences of Three 3-1,4-Endoglucanase Genes from Thermomonospora fusca. J Bacteriol 173: 3397-3407. |
[122] | Ali WAA, Gondal ZI, Yammahi AAAK, et al. (2013) A case of small bowel obstruction due to phytobezoars. J Surg Case Report 7: 1. |
[123] | Kramer SJ, Pochapin MB (2012) Gastric Phytobezoar Dissolution with Ingestion of Diet Coke and Cellulase. Gastroen Hepatol 8: 770-772. |
[124] |
Ladas SD, Kamberoglou D, Karamanolis G, et al. (2013) Systematic review: Coca-Cola can effectively dissolve gastric phytobezoars as a first-line treatment. Aliment Pharmacol Ther 37: 169-173. doi: 10.1111/apt.12141
![]() |
[125] |
Pinos N, Moreno-Merino S, Congregado M (2015) Phytobezoar by aloe vera as long term complication after oesophagectomy resolved using cellulase. Int J Surg Case Report 13: 37-39. doi: 10.1016/j.ijscr.2015.05.008
![]() |
[126] |
Robert C, Bernalier-Donadille A (2003) The cellulolytic micro£ora of the human colon: evidence of microcrystalline cellulose-degrading bacteria in methane-excreting subjects. FEMS Microbiol Ecol 46: 81-89. doi: 10.1016/S0168-6496(03)00207-1
![]() |
[127] | Wedekind KJ, Mansfield HR, Montgomery L (1988) Enumeration and isolation of cellulolytic and hemicellulolytic bacteria from human feces. App Environ Microbiol 54: 1530-1535 |
[128] | Robert C, Chassard C, Lawson PA, et al. (2007) Bacteroides cellulosilyticus sp. nov., a cellulolytic bacterium from the human gut microbial community. Int J Syst Evol Microbiol 57: 1516-1520 |
[129] | Chassard C, Delmas E, Robert C, et al. (2012) Ruminococcus champanellensis sp. nov., a cellulose-degrading bacterium from human gut microbiota. Int J Syst Evol Microbiol 62: 138-143. |
[130] | Martinez AJ, Visvesvara GS (1997) Free-living, amphizoic and opportunistic amebas. Brain Pathol 97: 583-98. |
[131] |
Lakhundi S, Siddiqui R, Khan NA (2015) Cellulose degradation: a therapeutic strategy in the improved treatment of Acanthamoeba infections. Parasite Vector 8: 23. doi: 10.1186/s13071-015-0642-7
![]() |
[132] |
Loiselle M, Anderson KW (2003) The use of cellulase in inhibiting biofilm formation from organisms commonly found on medical implants. Biofouling 19: 77-85. doi: 10.1080/0892701021000030142
![]() |
[133] |
Ma L, Conover M, Lu H, et al. (2009) Assembly and development of the Pseudomonas aeruginosa biofilm matrix. PLoS Pathog 5: e1000354. doi: 10.1371/journal.ppat.1000354
![]() |
[134] | Rajasekharan SK, Ramesh S (2013) Cellulase inhibits Burkholderia cepacia biofilms on diverse prosthetic materials. Polish J Microbiol 62: 327-330. |
[135] |
Huertas MG, Zarate L, Acosta IC, et al. (2014) Klebsiella pneumoniae yfiRNB operon affects biofilm formation, polysaccharide production and drug susceptibility. Microbiology 160: 2595-2606. doi: 10.1099/mic.0.081992-0
![]() |
[136] | Glick BR (2012) Plant growth-promoting bacteria: mechanisms and applications. Scientifica, 2012. |
[137] | Haldar S, Sengupta S (2015) Plant-microbe cross-talk in the rhizosphere: insight and biotechnological potential. Open Microbiol J 9: 1-7. |
[138] |
El-Tarabily KA, SykesML, Kurtböke ID, et al. (1996) Synergistic effects of a cellulose-producing Micromonospora carbonacea and an antibiotic-producing Streptomyces violascens on the suppression of Phytophthora cinnamomi root rot of Banksia grandis. Can J Botany 74: 618-624. doi: 10.1139/b96-078
![]() |
[139] |
Naing KW, Anees M, Nguyen XH, et al. (2014) Biocontrol of Late Blight Disease (Phytophthora capsici) of Pepper and the Plant Growth Promotion by Paenibacillus ehimensis KWN38. J Phytopathol 162: 367-376. doi: 10.1111/jph.12198
![]() |
[140] | Han W, He M (2010) The application of exogenous cellulase to improve soil fertility and plant growth due to acceleration of straw decomposition. Bioresource Technol 101: 3724-3731. |
[141] | Compant S, Reiter B, Sessitsch A, et al. (2005) Endophytic colonization of Vitis vinifera l. by plant growth-promoting bacterium Burkholderia sp. strain PsJN App Environ Microbiol 71: 1685-1693. |
1. | Rudrajeet Pal, Yasaman Samie, Armaghan Chizaryfard, Demystifying process-level scalability challenges in fashion remanufacturing: An interdependence perspective, 2021, 286, 09596526, 125498, 10.1016/j.jclepro.2020.125498 | |
2. | Shuo Zhu, Hua Zhang, Zhigang Jiang, Bernard Hon, A carbon efficiency upgrading method for mechanical machining based on scheduling optimization strategy, 2020, 15, 2095-0233, 338, 10.1007/s11465-019-0572-8 | |
3. | Wenjie Wang, Guangdong Tian, Honghao Zhang, Kangkang Xu, Zheng Miao, Modeling and scheduling for remanufacturing systems with disassembly, reprocessing, and reassembly considering total energy consumption, 2021, 0944-1344, 10.1007/s11356-021-17292-x | |
4. | Guangyuan Zhang, Xiaonan Gao, Zhenfang Zhu, Fengyv Zhou, Dexin Yu, Determination of the location of the needle entry point based on an improved pruning algorithm, 2022, 19, 1551-0018, 7952, 10.3934/mbe.2022372 | |
5. | Xingwang Shen, Shimin Liu, Can Zhang, Jinsong Bao, Intelligent material distribution and optimization in the assembly process of large offshore crane lifting equipment, 2021, 159, 03608352, 107496, 10.1016/j.cie.2021.107496 | |
6. | Yaping Ren, Xinyu Lu, Hongfei Guo, Zhaokang Xie, Haoyang Zhang, Chaoyong Zhang, A Review of Combinatorial Optimization Problems in Reverse Logistics and Remanufacturing for End-of-Life Products, 2023, 11, 2227-7390, 298, 10.3390/math11020298 | |
7. | Wenyu Zhang, Jun Wang, Xiangqi Liu, Shuai Zhang, A new uncertain remanufacturing scheduling model with rework risk using hybrid optimization algorithm, 2023, 1614-7499, 10.1007/s11356-023-26219-7 | |
8. | Jun Guo, Junfeng Zou, Baigang Du, Kaipu Wang, Integrated scheduling for remanufacturing system considering component commonality using improved multi-objective genetic algorithm, 2023, 182, 03608352, 109419, 10.1016/j.cie.2023.109419 | |
9. | Wenyu Zhang, Jiaxuan Shi, Shuai Zhang, Mengjiao Chen, Scenario-Based Robust Remanufacturing Scheduling Problem Using Improved Biogeography-Based Optimization Algorithm, 2023, 53, 2168-2216, 3414, 10.1109/TSMC.2022.3225443 |
Places and transitions in CTPN | Description |
s0 | The output buffer of the inspection center |
si (i=1, 2, ..., 12) | The input buffer of each recovery workstation |
sz (z=13, 14, 15, 16) | The output buffer of each recovery workstation |
s17 | The input buffer of the reassembly center |
wk (k=1, 2, ..., 12) | Workstation places |
ti, j (i, j=1, 2, ..., 12, i≠j) | Recovery operations |
t0, i, ti, z, ti,17 (i=1, 2, ..., 12, z=13, 14, 15, 16) | Transportation transitions |
Index | Description |
L= {1, 2, ..., L} | index set of used products |
D= {1, 2, ..., D} | index set of component types |
K= {1, 2, ..., K} | index set of workstations |
H= {1, 2, ..., H} | index set of recovery process routings |
Parameters | Description |
ϑld | the dth core disassembled from the lth product |
wk | the kth workstation that performs a specific recovery operation |
Rh | the hth recovery process routing designed for worn cores and is represented as a set of recovery operations |
pre(tij, Rh) | the operational stages of recovery operation tij in Rh |
Last(Rh) | the last recovery operation in Rh |
aldh | a binary variable with a value 1 indicates the hth routing for the core ϑld |
u(tij, ϑld) | the time for a core ϑld to be processed in transition tij |
AT(si, ϑld) | the arrival time of the core ϑld at buffer place si |
ST(wk, ϑld) | the start time of the core ϑld being processed by wk |
Com(wk, ϑld) | the completion time of the core ϑld processed by wk |
ρ(ϑld) | the makespan of the component ϑld |
χ(ϑld) | the tardiness of the component ϑld |
$\zeta (l)$ | the tardiness of the lth product |
Probld | the probability of a core ϑld being selected for routing reassignment |
$\Pi _{ld}^{h, x}$ | the difference of the total average waiting time between the hth routing and the xth one for the core ϑld |
WTk | the average waiting time per core consumed in the kth workstation |
h | the total average waiting time per core consumed in the hth routing |
$v_{ld}^{h, x}$ | the probability of the xth routing being selected to replace the current hth routing for the core ϑld |
Slackxy | slack time of a core ϑxy in scheduling horizon |
Decision variables | Description |
rldh | a binary variable with a value 1 indicates that the hth process routing is assigned to a core ϑld |
Pri(wk, ϑld) | the priority for the dth component of the lthused product to be processed by the kth workstation |
Arrival rate | The average total cost per machine tool | Improvement of the standard method over the baseline one tc1–tc2 | Improvement of the proposed method over the standard one tc2–tc3 | ||
Baseline casetc1 | Standard SA/MSTtc2 | Proposed SA/MSTtc3 | |||
π=5 | 186.66 | 117.02 | 117.02 | 69.64 | 0 |
π=6 | 228.37 | 133.48 | 118.21 | 94.89 | 15.27 |
π=7 | 293.05 | 172.45 | 146.89 | 120.6 | 25.56 |
π=8 | 335.87 | 208.06 | 168.69 | 127.81 | 39.37 |
π=9 | 447.52 | 266.47 | 216.25 | 181.05 | 50.22 |
π=10 | 581.63 | 352.28 | 270.91 | 229.35 | 81.37 |
π=11 | 714.29 | 461.86 | 335.77 | 252.43 | 126.09 |
Sample mean | 398.20 | 244.52 | 196.25 | 153.68 | 48.27 |
Sample standard deviation | 193.44 | 125.46 | 82.67 | 68.93 | 43.15 |
Arrival rate | The average total cost per machine tool | Improvement of the standard method over the baseline onetc1–tc2 | Improvement of the proposed method over the standard onetc2–tc3 | ||
Baseline casetc1 | Standard SA/MSTtc2 | Proposed SA/MSTtc3 | |||
π=5 | 215.76 | 117.83 | 117.83 | 97.93 | 0 |
π=6 | 327.61 | 160.31 | 140.12 | 167.3 | 20.19 |
π=7 | 455.9 | 226.47 | 185.02 | 229.43 | 41.45 |
π=8 | 658.09 | 342.57 | 273.22 | 315.52 | 69.35 |
π=9 | 891.23 | 515.96 | 406.38 | 375.27 | 109.58 |
π=10 | 1137.96 | 738.41 | 559.67 | 399.55 | 178.74 |
π=11 | 1530.54 | 986.65 | 775.72 | 543.89 | 210.93 |
Sample mean | 745.30 | 441.17 | 351.14 | 304.13 | 90.03 |
Sample standard deviation | 472.28 | 324.23 | 244.76 | 151.99 | 80.25 |
Arrival rate | The average total cost per machine tool | Improvement of the standard method over the baseline onetc1–tc2 | Improvement of the proposed method over the standard onetc2–tc3 | ||
Baseline casetc1 | Standard SA/MSTtc2 | Proposed SA/MSTtc3 | |||
π=5 | 270.84 | 117.32 | 117.32 | 153.52 | 0 |
π=6 | 463.08 | 224.8 | 186.58 | 238.28 | 38.22 |
π=7 | 687.42 | 416.46 | 275.13 | 270.96 | 141.33 |
π=8 | 1063.8 | 655.68 | 453.47 | 408.12 | 202.21 |
π=9 | 1495.67 | 957.36 | 706.49 | 538.31 | 250.87 |
π=10 | 1917.52 | 1231.56 | 918.76 | 685.96 | 312.8 |
π=11 | 2403.39 | 1483.89 | 1120.77 | 919.5 | 363.12 |
Sample mean | 1185.96 | 726.72 | 539.79 | 459.24 | 186.94 |
Sample standard deviation | 789.17 | 517.36 | 385.07 | 273.81 | 135.61 |
Computation time breaks | Standard SA/MST | Proposed SA/MST | ||||||
tc1 | pc1 | dc1 | wt1 | tc2 | pc2 | dc2 | wt2 | |
After 0sec | 399.25 | 112.16 | 287.09 | 76.80 | 399.25 | 112.16 | 287.09 | 76.80 |
After 30sec | 380.38 | 116.78 | 263.60 | 72.37 | 371.14 | 114.62 | 256.51 | 55.21 |
After 60sec | 393.70 | 113.71 | 279.99 | 61.76 | 341.42 | 118.56 | 222.86 | 48.34 |
After 90sec | 393.70 | 113.71 | 279.99 | 58.90 | 311.93 | 120.51 | 191.43 | 46.49 |
After 120sec | 386.11 | 115.62 | 270.49 | 58.75 | 292.79 | 121.75 | 171.04 | 44.54 |
After 150sec | 385.56 | 119.96 | 265.60 | 58.42 | 280.06 | 123.28 | 156.78 | 41.52 |
After 180sec | 384.17 | 121.18 | 262.99 | 57.90 | 274.14 | 124.64 | 149.51 | 41.47 |
After 210sec | 383.35 | 122.86 | 260.49 | 56.51 | 273.33 | 126.31 | 147.03 | 38.40 |
After 240sec | 378.22 | 124.65 | 252.57 | 54.73 | 273.22 | 132.79 | 140.42 | 35.67 |
No time limits* | 342.57 | 129.41 | 213.16 | 48.07 | 273.22 | 132.79 | 140.42 | 35.67 |
* "No time limit" means the simulation keeps the algorithm running until it meets the stopping criteria elaborated above. |
Places and transitions in CTPN | Description |
s0 | The output buffer of the inspection center |
si (i=1, 2, ..., 12) | The input buffer of each recovery workstation |
sz (z=13, 14, 15, 16) | The output buffer of each recovery workstation |
s17 | The input buffer of the reassembly center |
wk (k=1, 2, ..., 12) | Workstation places |
ti, j (i, j=1, 2, ..., 12, i≠j) | Recovery operations |
t0, i, ti, z, ti,17 (i=1, 2, ..., 12, z=13, 14, 15, 16) | Transportation transitions |
Index | Description |
L= {1, 2, ..., L} | index set of used products |
D= {1, 2, ..., D} | index set of component types |
K= {1, 2, ..., K} | index set of workstations |
H= {1, 2, ..., H} | index set of recovery process routings |
Parameters | Description |
ϑld | the dth core disassembled from the lth product |
wk | the kth workstation that performs a specific recovery operation |
Rh | the hth recovery process routing designed for worn cores and is represented as a set of recovery operations |
pre(tij, Rh) | the operational stages of recovery operation tij in Rh |
Last(Rh) | the last recovery operation in Rh |
aldh | a binary variable with a value 1 indicates the hth routing for the core ϑld |
u(tij, ϑld) | the time for a core ϑld to be processed in transition tij |
AT(si, ϑld) | the arrival time of the core ϑld at buffer place si |
ST(wk, ϑld) | the start time of the core ϑld being processed by wk |
Com(wk, ϑld) | the completion time of the core ϑld processed by wk |
ρ(ϑld) | the makespan of the component ϑld |
χ(ϑld) | the tardiness of the component ϑld |
$\zeta (l)$ | the tardiness of the lth product |
Probld | the probability of a core ϑld being selected for routing reassignment |
$\Pi _{ld}^{h, x}$ | the difference of the total average waiting time between the hth routing and the xth one for the core ϑld |
WTk | the average waiting time per core consumed in the kth workstation |
h | the total average waiting time per core consumed in the hth routing |
$v_{ld}^{h, x}$ | the probability of the xth routing being selected to replace the current hth routing for the core ϑld |
Slackxy | slack time of a core ϑxy in scheduling horizon |
Decision variables | Description |
rldh | a binary variable with a value 1 indicates that the hth process routing is assigned to a core ϑld |
Pri(wk, ϑld) | the priority for the dth component of the lthused product to be processed by the kth workstation |
Arrival rate | The average total cost per machine tool | Improvement of the standard method over the baseline one tc1–tc2 | Improvement of the proposed method over the standard one tc2–tc3 | ||
Baseline casetc1 | Standard SA/MSTtc2 | Proposed SA/MSTtc3 | |||
π=5 | 186.66 | 117.02 | 117.02 | 69.64 | 0 |
π=6 | 228.37 | 133.48 | 118.21 | 94.89 | 15.27 |
π=7 | 293.05 | 172.45 | 146.89 | 120.6 | 25.56 |
π=8 | 335.87 | 208.06 | 168.69 | 127.81 | 39.37 |
π=9 | 447.52 | 266.47 | 216.25 | 181.05 | 50.22 |
π=10 | 581.63 | 352.28 | 270.91 | 229.35 | 81.37 |
π=11 | 714.29 | 461.86 | 335.77 | 252.43 | 126.09 |
Sample mean | 398.20 | 244.52 | 196.25 | 153.68 | 48.27 |
Sample standard deviation | 193.44 | 125.46 | 82.67 | 68.93 | 43.15 |
Arrival rate | The average total cost per machine tool | Improvement of the standard method over the baseline onetc1–tc2 | Improvement of the proposed method over the standard onetc2–tc3 | ||
Baseline casetc1 | Standard SA/MSTtc2 | Proposed SA/MSTtc3 | |||
π=5 | 215.76 | 117.83 | 117.83 | 97.93 | 0 |
π=6 | 327.61 | 160.31 | 140.12 | 167.3 | 20.19 |
π=7 | 455.9 | 226.47 | 185.02 | 229.43 | 41.45 |
π=8 | 658.09 | 342.57 | 273.22 | 315.52 | 69.35 |
π=9 | 891.23 | 515.96 | 406.38 | 375.27 | 109.58 |
π=10 | 1137.96 | 738.41 | 559.67 | 399.55 | 178.74 |
π=11 | 1530.54 | 986.65 | 775.72 | 543.89 | 210.93 |
Sample mean | 745.30 | 441.17 | 351.14 | 304.13 | 90.03 |
Sample standard deviation | 472.28 | 324.23 | 244.76 | 151.99 | 80.25 |
Arrival rate | The average total cost per machine tool | Improvement of the standard method over the baseline onetc1–tc2 | Improvement of the proposed method over the standard onetc2–tc3 | ||
Baseline casetc1 | Standard SA/MSTtc2 | Proposed SA/MSTtc3 | |||
π=5 | 270.84 | 117.32 | 117.32 | 153.52 | 0 |
π=6 | 463.08 | 224.8 | 186.58 | 238.28 | 38.22 |
π=7 | 687.42 | 416.46 | 275.13 | 270.96 | 141.33 |
π=8 | 1063.8 | 655.68 | 453.47 | 408.12 | 202.21 |
π=9 | 1495.67 | 957.36 | 706.49 | 538.31 | 250.87 |
π=10 | 1917.52 | 1231.56 | 918.76 | 685.96 | 312.8 |
π=11 | 2403.39 | 1483.89 | 1120.77 | 919.5 | 363.12 |
Sample mean | 1185.96 | 726.72 | 539.79 | 459.24 | 186.94 |
Sample standard deviation | 789.17 | 517.36 | 385.07 | 273.81 | 135.61 |
Computation time breaks | Standard SA/MST | Proposed SA/MST | ||||||
tc1 | pc1 | dc1 | wt1 | tc2 | pc2 | dc2 | wt2 | |
After 0sec | 399.25 | 112.16 | 287.09 | 76.80 | 399.25 | 112.16 | 287.09 | 76.80 |
After 30sec | 380.38 | 116.78 | 263.60 | 72.37 | 371.14 | 114.62 | 256.51 | 55.21 |
After 60sec | 393.70 | 113.71 | 279.99 | 61.76 | 341.42 | 118.56 | 222.86 | 48.34 |
After 90sec | 393.70 | 113.71 | 279.99 | 58.90 | 311.93 | 120.51 | 191.43 | 46.49 |
After 120sec | 386.11 | 115.62 | 270.49 | 58.75 | 292.79 | 121.75 | 171.04 | 44.54 |
After 150sec | 385.56 | 119.96 | 265.60 | 58.42 | 280.06 | 123.28 | 156.78 | 41.52 |
After 180sec | 384.17 | 121.18 | 262.99 | 57.90 | 274.14 | 124.64 | 149.51 | 41.47 |
After 210sec | 383.35 | 122.86 | 260.49 | 56.51 | 273.33 | 126.31 | 147.03 | 38.40 |
After 240sec | 378.22 | 124.65 | 252.57 | 54.73 | 273.22 | 132.79 | 140.42 | 35.67 |
No time limits* | 342.57 | 129.41 | 213.16 | 48.07 | 273.22 | 132.79 | 140.42 | 35.67 |
* "No time limit" means the simulation keeps the algorithm running until it meets the stopping criteria elaborated above. |