diff --git a/notes.html b/notes.html index 7b850d2..3433b73 100644 --- a/notes.html +++ b/notes.html @@ -3,7 +3,7 @@ "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd">
- +This report is also available as a pdf.
++
Figure 1: Accuracy and Repeatability
Figure 2: Position Resolution
Figure 3: Position Stability
@@ -216,11 +189,11 @@ The stability is always related to the time frame taken into account. (Figure-Limited stiffness, play and friction will induce an hysteresis for a positioning system as shown in Figure 4. +Limited stiffness, play and friction will induce an hysteresis for a positioning system as shown in Figure 4.
@@ -228,7 +201,7 @@ The hysteresis can actually help estimating the play and friction present in the
-
Figure 4: Stiffness, play and Friction
@@ -244,7 +217,7 @@ Ways to make the hysteresis smaller:-The position uncertainty of a system can be estimated as follow (Figure 5): +The position uncertainty of a system can be estimated as follow (Figure 5):
\begin{equation} \text{Position Uncertainty} = \text{play} + 2 \times \text{Virtual Play} @@ -257,7 +230,7 @@ where the virtual play can be estimated as follow: \end{equation} -
Figure 5: Hysterestis, play and virtual play
@@ -281,11 +254,11 @@ Note that it is very difficult to make a system with constant friction in practi-Estimate the virtual play of the system in Figure 6 with following characteristics: +Estimate the virtual play of the system in Figure 6 with following characteristics:
Figure 6: Studied system for “Case 1”
@@ -333,43 +306,43 @@ And finally:There exist many conventional elements for constraining DoFs. Some of them are:
Figure 7: Ball Joint
Figure 8: Ball Bearing
Figure 9: Roller Bearing
Figure 10: Roller Rail Guide
@@ -377,19 +350,19 @@ Some of them are:-An example of a complaint element is shown in Figure 11. +An example of a complaint element is shown in Figure 11.
-
Figure 11: Example of 1dof constrained compliant element
@@ -399,28 +372,28 @@ An example of a complaint element is shown in Figure 11
Figure 12: Leaf springs
Figure 13: Folded Leaf springs
Figure 14: Flexure Pivots (5dof constrained)
@@ -428,11 +401,11 @@ These are generally used in combination with other folded leaf springs.-Parallel guiding can be made using two leaf springs (Figure 15): +Parallel guiding can be made using two leaf springs (Figure 15):
Figure 15: Parallel guiding
Figure 16: Example of bucklink
Figure 17: Reinforced leaf springs
@@ -472,11 +445,11 @@ This sag is predictible and reproducible:-Figure 18 shows a rotation compliant mechanism: +Figure 18 shows a rotation compliant mechanism:
Figure 18: Example of rotation stage using leaf springs
@@ -492,11 +465,11 @@ Figure 18 shows a rotation compliant mechanism:-Figure 19 shows a Z translation mechanism: +Figure 19 shows a Z translation mechanism:
Figure 19: Z translation using 5 struts
-An alternative is to use folder leaf springs (Figure 20), and this avoid the parasitic rotation. +An alternative is to use folder leaf springs (Figure 20), and this avoid the parasitic rotation.
-
Figure 20: Z translation using 5 folded leaf springs
@@ -528,22 +501,22 @@ An alternative is to use folder leaf springs (Figure 20-An X-Y-Rz stage can be done either using 3 struts (Figure 21) or using 3 folded leaf springs (Figure 22). +An X-Y-Rz stage can be done either using 3 struts (Figure 21) or using 3 folded leaf springs (Figure 22).
-
Figure 21: X,Y,Rz using 3 struts
-The compliant mechanism shown in Figure 23 only constrain the rotation about the y-axis. +The compliant mechanism shown in Figure 23 only constrain the rotation about the y-axis.
-
Figure 23: 5dof motion, only the Ry is constrained
@@ -567,8 +540,8 @@ The compliant mechanism shown in Figure 23 only constr-An example of a complex compliant mechanism is shown in Figure 24. +An example of a complex compliant mechanism is shown in Figure 24.
-
Figure 24: Design concept
-Figure 25 shown a reinforced part to avoid buckling and improve vertical stiffness. +Figure 25 shown a reinforced part to avoid buckling and improve vertical stiffness.
-
Figure 25: Use leaf springs instead of linear roller bearings
@@ -606,11 +579,11 @@ Figure 25 shown a reinforced part to avoid buckling an-A X-Y-Rz stage is shown in Figure 26. +A X-Y-Rz stage is shown in Figure 26. To make this stage usable for nano-metric positioning, the following ideas where used:
Figure 26: Example of X-Y-Rz positioning stage
@@ -653,11 +626,11 @@ To make this stage usable for nano-metric positioning, the following ideas where-Figure 27 shows a parallel mechanism that should be converted to a compliant mechanism. +Figure 27 shows a parallel mechanism that should be converted to a compliant mechanism. Its characteristics are:
Figure 27: Example of a parallel stage that should be converting to a compliant mechanism
@@ -686,11 +659,11 @@ The goals are to:-The solution is shown in Figure 28. +The solution is shown in Figure 28.
-Thin plates are very important for compliant mechanisms. @@ -737,7 +710,7 @@ where \(A\) is the area of the cross section.
-
Figure 29: A plate under torsion
@@ -745,8 +718,8 @@ where \(A\) is the area of the cross section.The close profile has much more torsional stiffness than the open profile. @@ -757,7 +730,7 @@ Just by opening the tube, we have a much smaller torsional stiffness (but almost
-
Figure 30: Stiffness comparison open and closed tube (torsion)
@@ -766,22 +739,22 @@ Just by opening the tube, we have a much smaller torsional stiffness (but almostWe have similar behavior with an open/closed box. -If we remove one side of the cube shown in Figure 31, we would have much smaller torsional stiffness along the axis perpendicular to the removed side. +If we remove one side of the cube shown in Figure 31, we would have much smaller torsional stiffness along the axis perpendicular to the removed side.
-
Figure 31: Closed box.
-If we use triangles, we obtain high torsional stiffness as shown in Figure 32. +If we use triangles, we obtain high torsional stiffness as shown in Figure 32.
-
Figure 32: Open box (double triangle)
@@ -793,11 +766,11 @@ On way to reinforce it is using triangles.-A nice way to have a 1dof flexure guiding with stiff frame is shown in Figure 33. +A nice way to have a 1dof flexure guiding with stiff frame is shown in Figure 33.
-
Figure 33: Box with integrated flexure guiding
@@ -807,12 +780,12 @@ A nice way to have a 1dof flexure guiding with stiff frame is shown in FigureQuestion: in chip manufacturing, how do developments in optical lithography impact the mechatronic design? @@ -830,15 +803,15 @@ Main developments:
In this presentation, only the exposure step is discussed (lithography).
-
Figure 34: Chip manufacturing loop
@@ -846,19 +819,19 @@ In this presentation, only the exposure step is discussed (lithography).
Figure 35: Imaging process - basics
@@ -866,8 +839,8 @@ This will induce a sinusoidal wave on the wafer as shown in Figure -Before, one chip was illumating at a time, but people wanted to make bigger chips. @@ -880,7 +853,7 @@ This implied many requirements in dynamics and accuracy!
-
Figure 36: From stepper to scanner
@@ -888,8 +861,8 @@ This implied many requirements in dynamics and accuracy!Both the reticle stage and wafer stage are moving. @@ -914,7 +887,7 @@ Which are solved by: -
Figure 37: Machine based on the dual stage scanners
@@ -922,20 +895,20 @@ Which are solved by:Water is used between the lens and the wafer to increase the “NA” and thus decreasing the “critical dimension”.
-The “hood” is there to prevent any bubble to enter the illumination area (Figure 38). +The “hood” is there to prevent any bubble to enter the illumination area (Figure 38). The position of the “hood” is actively control to follow the wafer stage (that can move in z direction and tilt).
-Three solutions are used for the positioning control of the “hood” system (Figure 39): +Three solutions are used for the positioning control of the “hood” system (Figure 39):
Figure 38: Hood System
Figure 39: Control system for the “hood”
@@ -959,8 +932,8 @@ Three solutions are used for the positioning control of the “hood” sThe multiple patterning approach adds few mechatronics challenges: @@ -980,16 +953,16 @@ This was solved by:
-Each stage is controlled with 6dof lorentz short stroke actuators (Figure 40). +Each stage is controlled with 6dof lorentz short stroke actuators (Figure 40). The magnet stage can move horizontally (due to reaction forces of the wafer stages): it asks as a balance mass.
-
Figure 40: Machine layout
@@ -997,8 +970,8 @@ The magnet stage can move horizontally (due to reaction forces of the wafer stagVacuum is required which implies: @@ -1026,7 +999,7 @@ Wafer stage: -
Figure 41: Schematic of the ASML EUV machine
@@ -1034,22 +1007,22 @@ Wafer stage:
Figure 42: The CD will be 8nm
-In order to do so, high “opening” of the optics is required which is very challenges because the reflectiveness of mirror is decreasing as high angle of incidence (Figure 43). +In order to do so, high “opening” of the optics is required which is very challenges because the reflectiveness of mirror is decreasing as high angle of incidence (Figure 43).
-
Figure 43: Change of reflection of a mirror as a function of the angle of indicence
@@ -1057,8 +1030,8 @@ In order to do so, high “opening” of the optics is required which isChallenges: @@ -1080,8 +1053,8 @@ In order to do so, high “opening” of the optics is required which is
The conclusions are: @@ -1102,12 +1075,12 @@ The conclusions are:
Weakly damped flexible modes of the mechanism can limit the performance of motion control systems. @@ -1118,36 +1091,36 @@ For discrete time controlled systems, there can be an additional limitation: ali
-
Figure 44: Example of high frequency lighlty damped resonances
-The aliasing of signals is well known (Figure 45). +The aliasing of signals is well known (Figure 45).
-However, aliasing in systems can also happens and is schematically shown in Figure 46. +However, aliasing in systems can also happens and is schematically shown in Figure 46.
-
Figure 45: Aliasing of Signals
Figure 46: Aliasing of Systems
-The poles of the system will be aliased and their location will change in the complex plane as shown in Figure 47. +The poles of the system will be aliased and their location will change in the complex plane as shown in Figure 47.
@@ -1163,7 +1136,7 @@ Therefore, the damping of the aliased resonances are foreseen to have larger dam
-
Figure 47: Aliasing of poles in the complex plane
@@ -1178,7 +1151,7 @@ Let’s consider two systems with a resonance:-Then looking at the same systems in the digital domain, one can see thathen the resonance is above the Nyquist frequency (Figure 48): +Then looking at the same systems in the digital domain, one can see thathen the resonance is above the Nyquist frequency (Figure 48):
Figure 48: Aliazed resonance shown on the Bode Diagram
Figure 49: Higher resonance frequency
@@ -1205,26 +1178,26 @@ Therefore, when identifying a low damped resonance, it could be that it comes fo-The aliased modes can for instance comes from local modes in the actuators that are lightly damped and at high frequency (Figure 50) +The aliased modes can for instance comes from local modes in the actuators that are lightly damped and at high frequency (Figure 50)
-
Figure 50: Local vibration mode that will be alized
-The proposed idea to better model aliasing resonances is to include more modes in the FEM software as shown in Figure 51 and then perform an order reduction in matlab. +The proposed idea to better model aliasing resonances is to include more modes in the FEM software as shown in Figure 51 and then perform an order reduction in matlab.
-
Figure 51: Common procedure and proposed procedure to include aliazed resonances
@@ -1232,32 +1205,32 @@ The proposed idea to better model aliasing resonances is to include more modes i
Figure 52: Example of the effect of aliased resonance on the open-loop
Figure 53: Example of the effect of aliased resonance on sensitivity function
@@ -1265,8 +1238,8 @@ The proposed idea to better model aliasing resonances is to include more modes iConcept: @@ -1290,7 +1263,7 @@ Similarly, \(\omega_{0zi}\) is the natural frequency \(\xi_{zi}\) is the damping
-Examples (Figure 54): +Examples (Figure 54):
Figure 54: Magnitude, Phase and Phase delay of 3 filters
@@ -1310,26 +1283,26 @@ Similarly, \(\omega_{0zi}\) is the natural frequency \(\xi_{zi}\) is the damping-The budgeting of the phase lag is done by expressing the phase lag of each element by a time delay (Figure 55) +The budgeting of the phase lag is done by expressing the phase lag of each element by a time delay (Figure 55)
-
Figure 55: Typical control loop with several phase lag / time delays
-The equivalent delay of each element are listed in Figure 56. +The equivalent delay of each element are listed in Figure 56.
-The filter order can be chosen depending on the frequency of the resonance. -Some example of Butterworth filters are shown in Figure 57 and summarized in Figure 58. +Some example of Butterworth filters are shown in Figure 57 and summarized in Figure 58.
-
Figure 57: Example of few Butterworth filters
Figure 58: Butterworth filters
@@ -1361,8 +1334,8 @@ Some example of Butterworth filters are shown in Figure 57The equivalent delay of a low pass (here second order) depends on its damping, since: @@ -1370,7 +1343,7 @@ The equivalent delay of a low pass (here second order) depends on its damping, s
-
Figure 59: Change of the phase delay with the damping of the filter
@@ -1379,8 +1352,8 @@ The equivalent delay of a low pass (here second order) depends on its damping, sThe phenomenon of aliasing of resonances: @@ -1415,12 +1388,12 @@ Anti-aliasing filter design:
-Technical choice: flexure based delta robot (Figure 60). +Technical choice: flexure based delta robot (Figure 60).
Figure 60: Picture of the Delta Robot
Figure 61: x1, x2 x3 are the motor positions. f1,f2 f3 are the force motors. x,y,z are the position of the final point in cartesian coordinates
@@ -1464,8 +1437,8 @@ The control architecture should be as simple as possible.Lagrange equations are used to model the dynamics of the delta robot. @@ -1473,33 +1446,33 @@ The motor positions are used as the general coordinate system.
-The system is then linearized around the working point (Figure 62). +The system is then linearized around the working point (Figure 62).
-
Figure 62: Linearized equations of the Delta Robot
-Then the parameters are identified from experiment (Figure 63). +Then the parameters are identified from experiment (Figure 63).
-
Figure 63: Identification fo the transfer function from \(F_1\) to \(x_1\)
-The measurement of the coupling is move complicated as shown in Figure 64. +The measurement of the coupling is move complicated as shown in Figure 64.
-
Figure 64: Problem of identifying the coupling between F1 and x2 at low frequency
@@ -1507,8 +1480,8 @@ The measurement of the coupling is move complicated as shown in Figure -Control requirements: @@ -1521,7 +1494,7 @@ Control requirements: -
Figure 65: Control concept used for the Delta robot
@@ -1529,11 +1502,11 @@ Control requirements:-A 3 axis servo control board as been developed (Figure 66) which includes: +A 3 axis servo control board as been developed (Figure 66) which includes:
-Step response of the current control loop is shown in Figure 66. +Step response of the current control loop is shown in Figure 66.
-
Figure 66: Step response for the current control loop
@@ -1567,40 +1540,40 @@ Step response of the current control loop is shown in Figure --XY renishaw interferometers used to verify the performance of the system (Figure 67). +XY renishaw interferometers used to verify the performance of the system (Figure 67).
-
Figure 67: Experimental setup to verify the performances of the system
-Some results are shown in Figures 68, 69 and 70. +Some results are shown in Figures 68, 69 and 70.
-
Figure 68: Circuit motion results and point to point motion results
Figure 69: Step response of the system
Figure 70: Measured dynamical errors
@@ -1609,8 +1582,8 @@ Some results are shown in Figures 68, -As a conclusion, here are the identified conditions for precise and high dynamic positioning: @@ -1629,19 +1602,19 @@ Resonances at mid frequencies are an issue for further improvements.
Flexible eigenmodes are present in every system component and leads to::
@@ -1649,25 +1622,25 @@ Flexible eigenmodes are present in every system component and leads to::
-
Figure 71: Limitation of the control bandwidth due to flexible eigenmodes
Figure 72: Coupling due to flexible eigenmodes
-In order to estimate the performances of a system, the sensitivity function can be used (Figure 73). +In order to estimate the performances of a system, the sensitivity function can be used (Figure 73).
-
Figure 73: Bode plot of a typical Sensitivity function
@@ -1675,11 +1648,11 @@ In order to estimate the performances of a system, the sensitivity function can-There are different way to analyse the sensitivity function base on different plants (Figure 74): +There are different way to analyse the sensitivity function base on different plants (Figure 74):
Figure 74: Visual representation of the three systems
@@ -1704,11 +1677,11 @@ One loop is closed at a time, and the coupling effects are taken into account.-In order to compare the use of the three systems to estimate the performances of a MIMO system, the system shown in Figure 75 is used. +In order to compare the use of the three systems to estimate the performances of a MIMO system, the system shown in Figure 75 is used. The 4 top masses are used to represent a payload that will add coupling in the system due to its resonances.
@@ -1717,7 +1690,7 @@ A diagonal PID controller is used. -
Figure 75: Schematic representation of the example system
@@ -1725,24 +1698,24 @@ A diagonal PID controller is used.-The bode plot of the MIMO system is shown in Figure 76 where we can see the resonances in the off-diagonal elements. +The bode plot of the MIMO system is shown in Figure 76 where we can see the resonances in the off-diagonal elements.
-
Figure 76: Bode plot of the full MIMO system
-In Figure 77 is shown that the sensitivity function computed from the SISO system is not correct. +In Figure 77 is shown that the sensitivity function computed from the SISO system is not correct. Whereas for the “interaction method” system, it is correct and almost match the full system sensibility. However, as expected, the off-diagonal sensibilities are not modelled.
-
Figure 77: Bode plots of sensitivity functions
@@ -1750,11 +1723,11 @@ However, as expected, the off-diagonal sensibilities are not modelled.-The conclusion are the following and summarized in Figure 78: +The conclusion are the following and summarized in Figure 78:
The goal of this project is to perform a topology optimization of a 6dof magnetic levitated stage suitable for vacuum.
-For the current system (Figure 79), the bandwidth is limited by the short-stroke dynamics (eigenfrequencies). +For the current system (Figure 79), the bandwidth is limited by the short-stroke dynamics (eigenfrequencies).
@@ -1797,7 +1770,7 @@ The goal here is to make the eigen-frequency higher as this will allow more band
-
Figure 79: Schematic of the 6dof levitating stage
@@ -1805,15 +1778,15 @@ The goal here is to make the eigen-frequency higher as this will allow more band-More precisely, the goal is to automatically maximize the three eigen-frequencies of the system shown in Figure 80. +More precisely, the goal is to automatically maximize the three eigen-frequencies of the system shown in Figure 80.
-
Figure 80: System to be optimized
@@ -1821,16 +1794,16 @@ More precisely, the goal is to automatically maximize the three eigen-frequencieThe manufacturing process must be embedded in the optimization such that the obtained design is producible. -The process is shown in Figure 81. +The process is shown in Figure 81.
-Problem: for a given volume, maximize the eigen-frequencies of the system.
-To do so, the system is discretized into small elements (Figure 82). +To do so, the system is discretized into small elements (Figure 82). Then, a Finite Element Analysis is performed to compute the eigen-frequencies of the system. Finally, for each element, the “gradient is computed” and we determine if material should be added or removed.
-This is done in 3D. The individual 1mm x 1mm x 1mm elements are shown in Figure 82. +This is done in 3D. The individual 1mm x 1mm x 1mm elements are shown in Figure 82. The number of elements is 1 million (=> 15 minutes per iteration to compute the 3 eigen-frequencies).
-
Figure 82: Results of the topology optimization and zoom to see individual elements
@@ -1865,26 +1838,26 @@ The number of elements is 1 million (=> 15 minutes per iteration to compute t-The obtained mass and eigen-frequencies of the optimized system and the solid equivalents are compared in Figure 83. +The obtained mass and eigen-frequencies of the optimized system and the solid equivalents are compared in Figure 83.
-
Figure 83: Comparison of the obtained performances
-Identification on the realized system shown that the obtained eigen-frequencies are very closed to the estimated ones (Figure 84). +Identification on the realized system shown that the obtained eigen-frequencies are very closed to the estimated ones (Figure 84).
-
Figure 84: Results very close to simulation (~1% for the eigen frequencies)
@@ -1892,8 +1865,8 @@ Identification on the realized system shown that the obtained eigen-frequenciesGoal: Need for higher quality FRF models that are used to: @@ -1927,7 +1900,7 @@ High quality FRFs requires careful design of excitation \(w\).
-Typical experimental identification of the FRFs is shown in Figure 85. +Typical experimental identification of the FRFs is shown in Figure 85.
@@ -1939,7 +1912,7 @@ The design trade-off is: -
Figure 85: schematic of the identification of the FRF
@@ -1968,14 +1941,14 @@ For MIMO systems:The classical way to estimate MIMO FRFs is the following:
Figure 86: Example of a SISO approach to identify MIMO FRFs
-When having a MIMO approach and choosing both the direction and gain of the excitation inputs, we can obtained much better FRFs uncertainty while still fulfilling the constraints (Figure 87). +When having a MIMO approach and choosing both the direction and gain of the excitation inputs, we can obtained much better FRFs uncertainty while still fulfilling the constraints (Figure 87).
-
Figure 87: Example of the MIMO approach that gives much better FRFs
@@ -2003,15 +1976,15 @@ When having a MIMO approach and choosing both the direction and gain of the exciThe optimization problem is to minimize the model uncertainty by choosing the design variables which are the magnitude and direction of the inputs \(w\).
-The optimization is a two step process as shown in Figure 88: +The optimization is a two step process as shown in Figure 88:
Figure 88: Two step optimization process
@@ -2036,8 +2009,8 @@ In this work, two algorithms are proposed and not further detailed here.Experimental identification of a 7x8 MIMO plant was performed in for different cases: @@ -2050,24 +2023,24 @@ Experimental identification of a 7x8 MIMO plant was performed in for different c
-The obtained FRFs are shown in Figure 89. +The obtained FRFs are shown in Figure 89.
-
Figure 89: Obtained MIMO FRFs
-A comparison of one of the obtained FRFs is shown in Figure 90. +A comparison of one of the obtained FRFs is shown in Figure 90. It is quite clear that the MIMO approach can give much lower FRF uncertainty. The RR proposed algorithm is giving the best results
-
Figure 90: Example of one of the obtained FRF
@@ -2075,8 +2048,8 @@ The RR proposed algorithm is giving the best results