Presentation SBO
Presentation SBO
Presentation SBO
Full Professor
Simulation-based Industrial Engineering Department
National Engineering School of Tunis (ENIT)
Optimization for
Smart Industry Head of the Industry 4.0 Laboratory
Email: amel.jaoua@enit.utm.tn
amel.jaoua@polymtl.ca
Professional Experience
2013 – Today
Full Professor - Industrial Engineering Department
@ National Engineering School of Tunis
Co-Founder and Head of the Industry 4.0 Laboratory
2009 – Today
Consultant and Expert in Simulation and Optimization for Decision-Making and
Control in Manufacturing and Service Systems
2009 – 2013
Assistant Professor – Mathematic and Industrial Engineering Department
Speciality Simulation of Industrial Systems.
@ Ecole Polytechnique Montreal-Canada
How the MITACS Canada Research Program Offers Funding for International Students ?
5
Agenda
1. Place of Simulation-based Optimization in Smart Industry
• 1.1 Definition
• 1.2 An example of Simulation Project at Pratt & Whitney Canada
2. Discrete Event Simulation Modelling (DES)
• 2.1 Low-Immersion in a Warehouse System
• 2.2 Simulation by Hand of Queueing behavior
• 2.3 DES model development
3. Coupling Simulator with Optimizer
• 3.1 Simulation-based Optimization (SBO) theory
• 3.2 UML sequencing and searching Algorithms
4. Application for Flexible Manufacturing System (FMS)
• 4.1 Complexity of the FMS
• 4.2 Formulation of the SBO problem
• 4.3 Implementation of the SBO for Makespan Minimization
5. Application in an Industry 4.0 Laboratory
• 5.1 Low-Immersion in the Laboratory 4.0
• 5.2 Definition of the Optimization Problem
• 5.3 Developing the Simulation model with ARENA
• 5.4 Solving the Optimization Problem by coupling ARENA with OptQuest
6
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
1.Place of Simulation-based optimization in Digital Twin concept
8
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
1.2 An example of Simulation Project
10
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
Agenda
13
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
2.2 Simulation by Hand of Queueing behavior
To do:
Service
Delay
L(t) Exit of Pr.3 &
Exit of Pr.1 & Exit of Pr.2 & Exit of Pr.4 & Exit of Pr.5 &
Entry of Pr.2 in Server Entry of Pr.3 in Server Entry of Pr.4 in Server Entry of Pr.5 in Server Entry of Pr.6 in Server
0
5 6 10 12 15 18 20 24 25 30
Entry of Pr.6 in Queue
Time (minutes)
Entry of Pr.2 in Queue Entry of Pr.3 in Queue Entry of Pr.4 in Queue Entry of Pr.5 in Queue
Entry of Pr.1 in Server
15
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
2.2 Simulation by Hand of Queueing behavior
∞
1 1 30
𝐿 = 𝑖𝑇𝑖 = න 𝐿 𝑡 𝑑𝑡 = 1.5 𝑝𝑟𝑜𝑑𝑢𝑐𝑡𝑠
𝑇 30 0
𝑖=0
𝑁 𝑄 𝑄 𝑄 𝑄 𝑄 𝑄
1 𝑄 𝑊1 + 𝑊2 + 𝑊3 + 𝑊4 + 𝑊5 + 𝑊6
WQ = 𝑊𝑖 = = 2.5 𝑚𝑖𝑛𝑢𝑡𝑒𝑠
N 6
𝑖=1
As we are interested by the steady-state performances of 𝐿 and 𝑊𝑄 we let 𝑇 → ∞
16
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
Limitation of analytical
modeling?
2.3 DES model development
Interconnected conditional
Non-Markovian Processes
flow behavior
Estimate 𝑳 and 𝑾𝑸
for the Warehouse
system?
Warehouse process
Returns Arrive: Quality control No Disposal
Functional?
Yes
Average outgoing
quality AOQ
Manufacturing process
Yes
Manufacturing at
Row materiels Check the stock
demand
arrive level xF(t)
d / ( 1 -AOQ)
No
Sufficient?
Manufacturing at
max
17
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
Steps in Simulation Project?
2.3 DES model development
Problem Definition
Model Translation
Verified ?
Implementation
Non
Validated ?
Non Non
18
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
Steps in Simulation project
2.3 DES model development
To do :
Draw an analogy to the previously exposed Pratt & Whitney Simulation Project for engine maintenance ?
•Conceptual Model : Define the boundaries and assumptions, then use a well-established formalism such
as UML activity diagrams or BPMN to represent the flow of processes.
•Data collection: Collect historical data from existing systems. This could include operational logs,
maintenance records, production data, etc. Then convert the collected data into a form that can be directly
input into the simulation model. This often includes standardizing formats, cleaning up datasets, fitting data
to appropriate probability distributions.
19
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
Steps in Simulation project
2.3 DES model development
•Model translation: coding the model using general purpose ( C#, Phython, etc.) or specific
simulation languages (Arena, Anylogic, Simulink, Simio, etc.)
• Verified?: Debug the computer program by using tracing to detect potential bugs.
•Validated?: Compare the model's performance indicators with the actual system performance data
gathered. Utilize statistical comparison methods, such as the Mann-Whitney U Test, for validation. If
inconsistencies are revealed by the statistical tests, the conceptual model and/or the data collection phases
must be revised.
•Design of Experiments to Optimize decision: Formulate the simulation-based optimization problem
to find efficient decision i.e. that leads to better system performances. Coupling the simulator with an
Optimization tools for the execution of the designed experiments.
20
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
Design of Experiments
2.3 DES model development
• Etc.
21
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
Agenda
23
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
3. Coupling Simulation with Optimizer 3.1 SBO Theory
• Decision Variables
• Objective Function
• Constraints
An estimator of the expected value can be obtained by the sample mean over N
replications:
𝑁
1
𝑓ҧ ≡ 𝐿(𝑥, 𝑤𝑖 )
𝑁
𝑖=1
Where the number of replication N is defined according to the desired width of the
confidence interval, generally fixed at significance level (α) of 5%.
OptQuest: SNDModel.p:
User
Optimizer Simulator
Optimize (Variance)
[[TimeUp]or[AutomaticStop]]
alt Simulate(Path)
[Variance == 0]
PerformanceIndicators (Cost)
Destroy
[Else]
[ConfLevel Met]
Simulate (Path, Stream) CMRG:
Disk
RNG
Request RN (Stream)
PerformanceIndicators (Cost)
Evaluate(ConfLevel)
Compute Half-width
Destroy
• Finite and Small Feasible Region: Ranking and selection methods can be effectively
utilized.
• Combinatorially Large Feasible Region: Metaheuristic algorithms are more suitable
for optimization.
Optquest also handle multi-objective optimization using the weighted sum method.
31
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
4. Application for FMS 4.1 Complexity of the FMS
Number of AGV Speed Dispatching rule Makespan Half-width Lower bound Makesapn Upper bound Makespan
34
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
4. Application for FMS 4.2 Formulation of the SBO problem
The Optimization problem consist on finding the optimal/suboptimal combination of the number of AGV, their
speed and the dispatching rule that leads to minimize the Makespan.
35
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
4. Application for FMS 4.2 Implementation of the SBO
To implement the Simulated Annealing algorithm, we need to define the problem-specific parameters and to
determine the annealing parameters.
Problem-specific parameters are the evaluation function and the neighbourhood structure
Select an initial state X0;
Select an initial temperature T0; 1 2 3
Select the epoch length E; X at searching step (i-1)
Repeat 3 8 3
Set i = 0
Repeat Random choice of three range
Generate a neighbourhood state Xi;
Evaluate by simulation mechanism fi = f(Xi);
Calculate δ := fi - fi-1; 1 2 3
New X to test at
If δ<0 then Xi-1 := Xi 9 30 1
Step (i)
else if random(0,1) <exp(-δ /T) then Xi-1 := Xi ;
i = i+1;
until i=E;
T := α T; Generation of neighborhood move
until T= Tf
Jaoua, A., Gamache, M., & Riopel, D. (2012). Specification of an
intelligent simulation-based real time control architecture:
Application to truck control system. Computers in Industry, 63(9),
882-894.
36
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
4. Application for FMS 4.2 Formulation of the SBO problem
42
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
5. Application in an Industry 4.0 Laboratory 5.1 Low-immersion
43
5. Application in an Industry 4.0 Laboratory 5.2 Formulation of the SBO problem
The optimization problem consists of finding the optimal number of operators to assign to
“Assembly_workstation_1” and to the “Packaging workstation”, along with the AMR
allocation rate for each assembly workstation, in order to Minimize the Makespan.
To do:
45
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
5. Application in an Industry 4.0 Laboratory 5.3 Developing the Simulation Model
46
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
5. Application in an Industry 4.0 Laboratory 5.4 Solving with OptQuest
48
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
A glimpse into the future of SBO
[https://github.com/neprev].
49
A glimpse into the future of SBO
50
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca
References
• Jaoua, A., Masmoudi, S., & Negri, E. (2024). Digital twin-based reinforcement learning framework: application
to autonomous mobile robot dispatching. International Journal of Computer Integrated Manufacturing, 1-24.
• Layeb, S. B., A. Jaoua, A. Jbira, and Y. Makhlouf. 2018. “A Simulation-Optimization Approach for Scheduling in
Stochastic Freight Transportation.” Computers & Industrial Engineering 126:99–110.
https://doi.org/10.1016/j.cie.2018. 09.021 .
• Glover, F., M. Laguna, and R. Marti. 2000. “Scatter Search and Path Relinking: Advances and Applications.”
Control and Cybernetics 29 (3): 653–684.
• Gosavi, A. 2015. Simulation-Based Optimization, 47–55. Berlin: Springer.
• S.B. Saoud, A. Jaoua, N.B. Saoud, Efficient Simulator Based on Meta-Heuristic for FMS and AGV systems Design
and Control, Journal of Computer Science, vol. 4, n°1, pp. 8-14, 2005
• Fu, M. C. 2002. “Feature Article: Optimization for Simulation: Theory Vs. Practice.” INFORMS Journal on
Computing 14 (3): 192–215. https://doi.org/10.1287/ijoc.14.3.192.113
51
Professor Amel Jaoua, Ph.D. amel.jaoua@polymtl.ca