motion planning of circular water tank
TRANSCRIPT
![Page 1: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/1.jpg)
Motion Planning of Water Tank
Zhi Li
5/8/2014
![Page 2: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/2.jpg)
Outline
Problem Definition
Mathematical Derivation
MATLAB Simulation
Conclusion
Future Perfection
![Page 3: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/3.jpg)
Problem Definition
Math Derivation
MATLAB Simulation
ConclusionFuture
Perfection
ProblemSloshing when moving tanks containing liquids from one place to another
GoalSimulate tank dynamicsMinimize sloshing
ApproachesCalculate flatnessSelect proper output function
![Page 4: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/4.jpg)
Problem Definition
Math Derivation
MATLAB Simulation
ConclusionFuture
Perfection
Assumptions
Irrotational potential flow
Shallow water
Flat Bottom
1-D motion of tank
Petit, 2002
![Page 5: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/5.jpg)
Problem Definition
Math Derivation
MATLAB Simulation
ConclusionFuture
Perfection
Notationsh(x,y,t) = Actual water depthH(x,y,t) = Actual depth minus average depthD = Location of the tank’s bottom centerv(x,y,t) = Velocity of water
Basic EquationsContinuity
Momentum
Petit, 2002
![Page 6: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/6.jpg)
Problem Definition
Math Derivation
MATLAB Simulation
ConclusionFuture
Perfection
Define Input & OutputInput D(t) – Position of tank center on the bottom
Output y(t) – Center of gravity of tank’s first half mass and second half mass
Dubois, 1999
![Page 7: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/7.jpg)
Problem Definition
Math Derivation
MATLAB Simulation
ConclusionFuture
Perfection
System dynamics for 1D model
Recall d’Alembert’s formula
Plug into the BC’s to solve for D(t)
![Page 8: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/8.jpg)
Problem Definition
Math Derivation
MATLAB Simulation
ConclusionFuture
Perfection
System dynamics for 1D model
This approach can be extended to 2D case
![Page 9: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/9.jpg)
Problem Definition
Math Derivation
MATLAB Simulation
ConclusionFuture
Perfection
System dynamics for 2D rectangular tank
Recall:H = Fluctuation of surfaceD = Location of tank bottom centerK, I = Global coordinate systemk, i = Local systemΩ = Tank boundary
2D model can be decoupled to get 1D model
![Page 10: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/10.jpg)
Problem Definition
Math Derivation
MATLAB Simulation
ConclusionFuture
Perfection
System dynamics for circular tank
Laplace transform
SOV in cylindrical coordinates
Can solve for H
![Page 11: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/11.jpg)
Problem Definition
Math Derivation
MATLAB Simulation
ConclusionFuture
Perfection
System dynamics for circular tank
![Page 12: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/12.jpg)
Problem Definition
Math Derivation
MATLAB Simulation
ConclusionFuture
Perfection
Selecting y(t) function
Steady state at start and end points – ControllableTank moves from 0 to d in time T – D(t) from 0 to dGot four BC’s to restrict y(t):
y(0) = 0, y(T) = d, y’(0) = 0, y’(T) = 0
y(t) should roughly look like this:
y(t)
t
![Page 13: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/13.jpg)
Problem Definition
Math Derivation
MATLAB Simulation
ConclusionFuture
Perfection
Four y(t) functions are tested
Cosine Function 3rd Order Polynomial
Error Function Arc + Line + Arc
![Page 14: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/14.jpg)
Problem Definition
Math Derivation
MATLAB Simulation
ConclusionFuture
Perfection
D(t) – y(t) Plots2 x 2 Rectangular tankWater depth h = 1Travel time T = 5Time step dt = 0.001Travel distance d = 3.19
![Page 15: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/15.jpg)
Problem Definition
Math Derivation
MATLAB Simulation
ConclusionFuture
Perfection
D(t) – y(t) Plots2 x 2 Rectangular tankWater depth h = 1Travel time T = 5Time step dt = 0.001Travel distance d = 3.19
![Page 16: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/16.jpg)
Problem Definition
Math Derivation
MATLAB Simulation
ConclusionFuture
Perfection
Mid-presentation Review
Two systems: Rectangular Tank and Circular Tank
Three output functions: Cosine, Polynomial and Erf
Goal: Minimize H
![Page 17: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/17.jpg)
Problem Definition
Math Derivation
MATLAB Simulation
ConclusionFuture
Perfection
Simulation Video
![Page 18: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/18.jpg)
Problem Definition
Math Derivation
MATLAB Simulation
ConclusionFuture
Perfection
Analysis
The standard deviation of the water surface can be used to evaluate sloshing
Rectangular Tank Circular Tank
Cosine Function 0.0528 0.0323
Polynomial 0.0584 0.0310
Error Function 0.0830 0.0464
![Page 19: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/19.jpg)
Problem Definition
Math Derivation
MATLAB Simulation
ConclusionFuture
Perfection
Conclusion
Cosine and 3rd order polynomial outputs have
better performance
Circular tank produces less sloshing than
rectangular tank
![Page 20: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/20.jpg)
Problem Definition
Math Derivation
MATLAB Simulation
ConclusionFuture
Perfection
Future Perfection
Simulate non-linear moving trajectories
Extend to more complicated geometries
Find proper initial and final values of output
Test more output functions
![Page 21: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/21.jpg)
Problem Definition
Math Derivation
MATLAB Simulation
ConclusionFuture
Perfection
References
[1] N. Petit, P. Rouchon,”Dynamics and solutions to some control problems for water-tank systems”, IEEE Trans. Contr. Syst. Technol., vol. 47, pp.594 -609 2002
[2] F. Dubois, N. Petit and P. Rouchon, “Motion planning and nonlinear simulation for a tank containing a fluid”, European Control Conference, 1999
[3] D. Ho, B. Grunloh and N. Jenson, “Motion planning of a water tank with differential flatness and optimal control”, CE291 Final Project, 2012
![Page 22: Motion Planning of Circular Water Tank](https://reader031.vdocument.in/reader031/viewer/2022012507/6182c27d68a768589561dd3d/html5/thumbnails/22.jpg)
Problem Definition
Math Derivation
MATLAB Simulation
ConclusionFuture
Perfection
Thanks!