Advance Modeling of a Skid-Steering Mobile Robot for Remote ...
Advance Modeling of a Skid-Steering Mobile Robot for Remote ... Advance Modeling of a Skid-Steering Mobile Robot for Remote ...
3.2 Characterization of the Vibrations 40order Taylor expansion, i.e. cos x ≈ 1, sin x ≈ x, so that we can rewrite (3.6) as follows:a x = ˙v x + gψa y = ˙v y − p ax ˙ω z − gθ(3.7)a z = ˙v z + p ax ˙ω y − gThe data are acquired from the accelerometers at a sample rate of 1 kHz by using LabVIEWand then exported and saved in an Excel format, ready to be loaded and analyzed in Matlab.The Matlab functions that are used for the frequency analysis of the signals are defined inAppendix A. For every value of the control input ω ∗ z contained in the variable rangeW, themain program first parses the data by calling the function parseData(), which loads thedata from the Excel file, computes their Fast Fourier Transform (FFT), converts them fromVoltage values to acceleration values and eventually plots them and their FFT. To improvethe accuracy and reduce the computation time of the FFT function implemented in Matlab,the FFT should be computed on the closest lower power of 2 of samples [20]. In our case, asthe data are acquired for 20 seconds, only the central 2 14 over 20000 samples are consideredfor the FFT. Then, by calling the function filterData(), the program low-pass filters thedata at 50 Hz by using a 10 th order Butterworth filter in order to cut off the noise due to alternatecurrent which works at 60 Hz. Such a choice of the cut-off frequency is also motivatedby the fact that, by taking a first look to the FFT of all the data, there never appear significantfrequency components higher than 50 Hz, except for the alternate current frequency around60 Hz. After low-pass filtering, the amplitude of the signal is calculated as difference betweentheir maximum and minimum value, the FFT of the low-pass filtered data is computedand the three highest peaks in the FFT are found by calling the function findPeakFreq().The findPeakFreq() function also band-pass filter the data by a narrow band, generally2 Hz, centered on the maximum peak frequency. Finally, the original data are plotted togetherwith the low-pass and band-pass filtered data, in order to check the reliability of themaximum peak frequency obtained, and all the amplitudes and peak frequencies are plottedwith respect to the robot angular velocity.Let us consider the data acquired while swiveling the robot in place at different desired angularvelocity ω ∗ z, separately on concrete, tile and carpet. In order to check the repeatabilityof the data, the accelerations are also acquired three times for each angular velocity ω ∗ z. Becausethe wheels start lifting from the floor generally when ω ∗ z > 64 deg , making the robots
3.2 Characterization of the Vibrations 41continuously hopping on its wheels and eventually making the motors stalling, we will beconsidering only the data acquired for ω ∗ z ≤ 64 degs .Some example graphs of the low-pass filtered data with their FFT are provided in AppendixB. The graphs representing the average of the data amplitudes and the average of the threehighest peak frequencies with respect to the angular velocity ω ∗ z and the type of ground arealso provided in Appendix B.Although the data acquired during different trials present a sort of irregularity in their peakfrequencies, due to the irregularity of the ground, there are some features which characterizethem. In particular, by looking at the graphs depicted in Appendix B, the main features canbe listed as follows:• There are, usually, two main frequency components (generally the two highest peak inthe FFT), which can be clearly identified for every angular velocity ω ∗ z and for all thethree types of ground.• The two main components are, usually, one at ¨low¨ (< 10 Hz) and the other at ¨high¨frequency (> 20 Hz), when swiveling the robot with relatively low angular velocities(ω ∗ z < 30 degs ).• The two main components are both at ¨low¨ frequency (< 10 Hz), when swiveling therobot with relatively high angular velocities (ω ∗ z > 30 degs ).• When running the robot on carpet, the ¨low¨ frequency component is always the dominantcomponent.• When running the robot on concrete and tile, the ¨high¨ frequency component is, usually,the dominant component for relatively low angular velocities (ω ∗ z < 30 deg ), whilesthe ¨low¨ frequency component is, usually, the dominant component for relatively highangular velocities (ω ∗ z > 30 degs ).• The ¨high¨ frequency components seem to increase when increasing the robot angularvelocity until 30 degs .• When running the robot on carpet, the amplitude of the oscillations increases almostlinearly with respect to the angular velocity.• When running the robot on concrete and tile, the amplitude of the oscillations is relativelysmall and constant for low angular velocities (ω ∗ z< 30 deg ), while it becomess
- Page 1 and 2: Università degli Studi di GenovaRo
- Page 3 and 4: IINothing great in the World has be
- Page 5: CONTENTSIV6 Simulation 786.1 Simuli
- Page 8 and 9: LIST OF FIGURESVII6.4 (a) Block sch
- Page 10 and 11: List of Tables5.1 Table of the geom
- Page 12 and 13: Introduction 2boDynamics, VGo Commu
- Page 14 and 15: Introduction 4The robot employed, i
- Page 16 and 17: Chapter 1Hardware and SoftwareIn th
- Page 18 and 19: 1.1 The Pioneer 3-AT Robot 8four ho
- Page 20 and 21: 1.2 The sensors 10control resolutio
- Page 22 and 23: 2.1 State of the Art 12not negligib
- Page 24 and 25: 2.1 State of the Art 14velocity vec
- Page 26 and 27: 2.1 State of the Art 16Figure 2.3:
- Page 28 and 29: 2.1 State of the Art 18Finally, by
- Page 30 and 31: 2.1 State of the Art 20constraint (
- Page 32 and 33: 2.1 State of the Art 22approximates
- Page 34 and 35: 2.2 Generalized Modeling 24⎡1 0¯
- Page 36 and 37: 2.2 Generalized Modeling 26V = [v T
- Page 38 and 39: 2.2 Generalized Modeling 28where T(
- Page 40 and 41: 2.2 Generalized Modeling 30⎡g R(J
- Page 42 and 43: 2.2 Generalized Modeling 32where we
- Page 44 and 45: Chapter 3Experimental DataIn this c
- Page 46 and 47: 3.1 The accelerometer 36(a)(b)Figur
- Page 48 and 49: 3.2 Characterization of the Vibrati
- Page 52 and 53: 3.2 Characterization of the Vibrati
- Page 54 and 55: 3.3 Characterization of the Tire Fo
- Page 56 and 57: 3.3 Characterization of the Tire Fo
- Page 58 and 59: 3.3 Characterization of the Tire Fo
- Page 60 and 61: 3.3 Characterization of the Tire Fo
- Page 62 and 63: 3.3 Characterization of the Tire Fo
- Page 64 and 65: 3.3 Characterization of the Tire Fo
- Page 66 and 67: 4.1 Tire Lateral Force 56Figure 4.1
- Page 68 and 69: 4.1 Tire Lateral Force 58velocity b
- Page 70 and 71: 4.1 Tire Lateral Force 600 ≤ µ i
- Page 72 and 73: 4.3 Tire Vertical Force 62f ix (t)
- Page 74 and 75: Chapter 5Identification of the Robo
- Page 76 and 77: 5.1 Identification of the Geometric
- Page 78 and 79: 5.2 Identification of the Tire Dyna
- Page 80 and 81: 5.2 Identification of the Tire Dyna
- Page 82 and 83: 5.2 Identification of the Tire Dyna
- Page 84 and 85: 5.2 Identification of the Tire Dyna
- Page 86 and 87: 5.2 Identification of the Tire Dyna
- Page 88 and 89: Chapter 6SimulationIn this chapter,
- Page 90 and 91: 6.1 Simulink Model 80On the right h
- Page 92 and 93: 6.1 Simulink Model 82Figure 6.3: Bl
- Page 94 and 95: 6.1 Simulink Model 84(a)(b)Figure 6
- Page 96 and 97: 6.1 Simulink Model 86that v icx can
- Page 98 and 99: 6.2 Simulation Results 88(a)(b)Figu
3.2 Characterization <strong>of</strong> the Vibrations 41continuously hopping on its wheels and eventually making the motors stalling, we will beconsidering only the data acquired <strong>for</strong> ω ∗ z ≤ 64 degs .Some example graphs <strong>of</strong> the low-pass filtered data with their FFT are provided in AppendixB. The graphs representing the average <strong>of</strong> the data amplitudes and the average <strong>of</strong> the threehighest peak frequencies with respect to the angular velocity ω ∗ z and the type <strong>of</strong> ground arealso provided in Appendix B.Although the data acquired during different trials present a sort <strong>of</strong> irregularity in their peakfrequencies, due to the irregularity <strong>of</strong> the ground, there are some features which characterizethem. In particular, by looking at the graphs depicted in Appendix B, the main features canbe listed as follows:• There are, usually, two main frequency components (generally the two highest peak inthe FFT), which can be clearly identified <strong>for</strong> every angular velocity ω ∗ z and <strong>for</strong> all thethree types <strong>of</strong> ground.• The two main components are, usually, one at ¨low¨ (< 10 Hz) and the other at ¨high¨frequency (> 20 Hz), when swiveling the robot with relatively low angular velocities(ω ∗ z < 30 degs ).• The two main components are both at ¨low¨ frequency (< 10 Hz), when swiveling therobot with relatively high angular velocities (ω ∗ z > 30 degs ).• When running the robot on carpet, the ¨low¨ frequency component is always the dominantcomponent.• When running the robot on concrete and tile, the ¨high¨ frequency component is, usually,the dominant component <strong>for</strong> relatively low angular velocities (ω ∗ z < 30 deg ), whilesthe ¨low¨ frequency component is, usually, the dominant component <strong>for</strong> relatively highangular velocities (ω ∗ z > 30 degs ).• The ¨high¨ frequency components seem to increase when increasing the robot angularvelocity until 30 degs .• When running the robot on carpet, the amplitude <strong>of</strong> the oscillations increases almostlinearly with respect to the angular velocity.• When running the robot on concrete and tile, the amplitude <strong>of</strong> the oscillations is relativelysmall and constant <strong>for</strong> low angular velocities (ω ∗ z< 30 deg ), while it becomess