Next Article in Journal
An Impedimetric Biosensor Based on Ionic Liquid-Modified Graphite Electrodes Developed for microRNA-34a Detection
Previous Article in Journal
Fully Bayesian Prediction Algorithms for Mobile Robotic Sensors under Uncertain Localization Using Gaussian Markov Random Fields
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Sensor Based on a Spherical Parallel Mechanism for the Measurement of Fluid Velocity: Physical Modelling and Computational Analysis

by
Roque Saltarén
1,2,*,
Gerardo Portilla
1,2,*,
Alejandro R. Barroso
1,2,* and
Juan Cely
1,2,*
1
Universidad Politécnica de Madrid, Escuela Técnica Superior De Ingenieros Industriales, 28006 Madrid, Spain
2
Centro de Automática y Robótica, 28006 Madrid, Spain
*
Authors to whom correspondence should be addressed.
Sensors 2018, 18(9), 2867; https://doi.org/10.3390/s18092867
Submission received: 24 July 2018 / Revised: 27 August 2018 / Accepted: 27 August 2018 / Published: 30 August 2018
(This article belongs to the Section Physical Sensors)

Abstract

:
In this article, a new method was developed to measure the velocity of a fluid using a sensor, based on the use of a spherical parallel mechanism with three degrees-of-freedom (DOF). This sensor transforms the kinetic energy of the fluid into potential energy by deforming the parallel mechanism. This deformation is due to the impact of the fluid on a sphere attached to the platform of the parallel mechanism. Through the acquisition of data from a sensor using an inertial measurement unit (IMU) in the sphere, an algorithm calculates the velocity and direction of the fluid. In this article, a mathematical model of the mechanism and an algorithm for correctly measuring the velocity and direction of the fluid is developed; this algorithm is tested through a simulation in the Adams software, and the MATLAB software is used to execute the algorithm. The results show that the algorithm calculates the velocity and the direction of the fluid correctly, demonstrating the technical feasibility of the sensor.

1. Introduction

Fluid velocity sensors can be classified by their measuring principles, whether they use propeller meters, ball meters, electromagnetic meters, or acoustic meters [1]. Unlike these types of sensors, a novel sensor for measuring the velocity and direction of the fluid is proposed in this article. This is possible using a spherical parallel mechanism with three degrees-of-freedom (DOF). The parallel mechanism (PM) can be defined as a multi-degree-of-freedom mechanism, composed of a moving platform and a base connected by at least two serial kinematic chains in parallel [2]. One of the first applications of this type of mechanism is believed to be the tire testing machine introduced by Gough and Whitewall [3], followed by the motion simulation platform built by Stewart [4]. Another example of this class of mechanisms is the Delta robot [5].
A feature of the spherical parallel mechanism is that the degrees-of-freedom of its end effector is rotational. For example, Gosselin and Angeles [6] proposed a spherical three-degree-of-freedom parallel manipulator. This manipulator uses three actuated spherical chains with rotary actuators with axes converging to a point that is the center of rotation. Another parallel mechanism is that of Lande [7]; this has linear actuators ensure a motion of the joint centers, and the rotation around the plate normal is obtained by an articulated stem. Lu and Zhang [8] propose a new parallel mechanism kind: a three-leg, five-DoF 2 SPS + 1 UPU. The researchers Saltaren and Sabater [9] simplified this parallel mechanism to a three-DOF rotational effector with the configuration 2 UPS + 1 RU. We chose to use this last spherical parallel mechanism as a sensor to measure the velocity of a fluid because of its simple mechanical features and its few singularities in the workspace.
There are other devices proposed in the literature that can potentially be used as flow sensors. In these types of devices, the measurement of forces and torques are based on strain gages, for example, such as that implemented in the “compliant-parallel-mechanisms” that appear in [10,11]. The flow measurement for to be implemented in devices, based on stress–strain concepts, should be different from the kinematic models of the parallel mechanisms (as proposed in this article). For this reason, the compliant devices would require specific adaptations and different algorithm formulation in order to measure flows.
Regarding the use of a parallel mechanism as a sensor, the research by Eileen and Ying [12] proposed a novel 3D force sensor based on a parallel mechanism, which displayed good sensitivity and high resolution. Yao and Hou [13] conducted a theoretical and experimental investigation into the isotropy performance of a pre-stressed force sensor based on a Stewart platform. Lu and Wang [14] designed and analyzed a novel force/torque sensor for a hybrid hand with three fingers.
With these properties of parallel mechanisms, the contribution of this article is a novel kind of sensor for measuring the velocity of a fluid based on these mechanisms, as evidenced by the patent number ES2525394B2 (see the section on patents). This sensor measures the velocity and direction of the fluid through a parallel mechanism that has passive actuators (using springs). The platform of the mechanism is connected to a spherical body; within this body, there is an inertial measurement unit (IMU). When the fluid collides with the sphere, deforming the parallel mechanism, it is possible to know the angles of the platform and the acceleration of the sphere using the IMU. These angles and acceleration are input into an algorithm that is developed in this article, and the velocity and direction of the fluid are calculated. This is possible with the model’s mathematics and the dynamics of the parallel mechanism.

1.1. Motivation

The development of the proposed sensor has arisen from the need to have a real-time velocity measurement system in order to control the navigation of the underwater robot [15]. Usually, the velocity measurement of the underwater robot is made by acoustic signals and sophisticated algorithms. Kalman filters and particle filters are used for the treatment of these signals [16,17]. Normally, this method uses the combination of signals from an inertial navigation system (INS) and a speed log doppler (DVL) system to provide 3D velocity to the navigation control of the underwater vehicle. Theoretically, the INS/DVL signal combination requires sufficient beam measurements (at least three) to calculate the 3D velocity. However, in cases where the DVL only has limited beam measurements (less than three), the slightly coupled navigation system can no longer function, and only the INS works. Therefore, the navigation error accumulates over time [18]. The proposed sensor, unlike the usual acoustic methods, works on the conversion of kinetic energy into pressure. Due to the features of the velocity sensor, it is possible to believe that it would have a very significant impact for underwater robot uses.

1.2. Outline

The rest of this article is organized as follows. In Section 2, we describe the sensor as in the patent, as well as the model’s mathematics, the algorithm for the sensor, and the experimental testing of the algorithm. Section 3 shows the results of the experimentation, and Section 4 presents the discussion about the results obtained. Finally, Section 5 presents the patent for this sensor.

2. Materials and Methods

2.1. Sensor Description

Figure 1a shows the sensor as described in the patent (see the section on the patent). This sensor has three important parts: the sphere, neck, and the parallel mechanism. The sphere is a hollow part, and inside this there is an inertial measurement unit (IMU). When the fluid drags this sphere, the IMU calculates the inclination angle and the acceleration of the sphere, as seen in Figure 1a. This inclination is transmitted through the neck to the parallel mechanism. This mechanism has three DOF (see Figure 1b) with the configuration 2 UPS + 1 RU. It has three passive actuators, two prismatic (L1 and L2) and one rotational (A1). These actuators are passive because they only use springs. The joints of the parallel mechanism include three universal joints (B1, A2, and A3): one in the rotational actuator and two in the base of the prismatic actuator. It also has two spherical joints (B2 and B3) above the prismatic actuators.
The working method of this sensor is as follows. When the fluid drags the sphere, this deforms the parallel mechanism; by using the IMU, it is possible to know the inclination angle of the platform of the parallel mechanism and, by using inverse kinematics, the deformation of the springs is calculated. Knowing the deformation of the springs, we can ascertain the force of each actuator. Then, by using dynamics, the drag force in the sphere is calculated. Finally, the velocity of the fluid by the hydrodynamic equation of the drag force is calculated because it is proportional to the velocity of the fluid.

2.2. Kinematics of the Parallel Mechanism

Figure 2 shows a simplified diagram of the parallel mechanism, where the axes of the coordinates in Joint 1 can be seen. It should also be noted that the triangle of the platform is isosceles. The line between 1 and 2 represents the rotation joint (M), and the lines between 3 and 4 (L1) and between 5 and 6 (L2) show passive linear actuators.

2.2.1. Calculation of the Inverse Kinematics

To develop the inverse kinematics, the position of each joint is calculated with the rotation matrices [19,20], all in the frame of reference of X, Y, and Z located in Joint 1.
After calculating at each point, we obtain
[ x 3 y 3 z 3 ] = R ( q x , q y , q z ) [ b e 0 ] ,   [ x 5 y 5 z 5 ] = R ( q x , q y , q z ) [ b e 0 ] ,
where qx, qy, and qz are the angles of orientation of the platform in Joint 1. Therefore, the actuator would be determined by the equation
L 1 = ( x 3 x 4 ) 2 + ( y 3 y 4 ) 2 + ( z 3 z 4 ) 2 ,
L 2 = ( x 5 x 6 ) 2 + ( y 5 y 6 ) 2 + ( z 5 z 6 ) 2 ,
M = q z ,
In this way, the displacements of the actuators L1 and L2 and the rotation of M are calculated, depending on the orientation angles.

2.2.2. Calculation of the Jacobian Matrix

The calculation of the Jacobian matrix was performed by using the screw theory [21,22]. Therefore, after calculating the twist in each joint, we obtain
$ ^ 1 = [ u 1 i b i   x   u 1 i ] ,   $ ^ 2 = [ u 2 i b i   x   u 2 i ] ,   $ ^ 3 = [ 0 u 3 i ] ,   $ ^ 4 = [ u 4 i a i   x   u 4 i ] , $ ^ 5 = [ u 5 i a i   x   u 5 i ] , $ ^ 6 = [ u 6 i a i   x   u 6 i ] ,
$ ^ 7 = [ u 7 0 ] ,   $ ^ 8 = [ u 8 0 ] ,   $ ^ 9 = [ u 9 0 ] ,
Now, we identify a screw ($r) that is reciprocal [23,24] to all the joints, which is given by
$ ^ r = [ ( a 2 x u 3 , i ) T u 3 , i T u 9 T 0 T ] ,
After applying the product between reciprocals [25,26], we see that
$ ^ r , i T   $ p = v i $ ^ r , i T   $ i ,
In this way, the Jacobian matrix is obtained
  J = [ ( a 2 x u 3 , 1   ) T ( a 2 x u 3 , 2 ) T u 9 T ] ,

2.2.3. Analysis of Workspace and Singularities of the Parallel Mechanism

An important limitation of the parallel mechanism is that singular configurations may exist within the workspace, where the mechanism can gain one or more degrees-of-freedom and become uncontrollable. Furthermore, the actuator forces may become very large and may result in a breakdown of the mechanism. Therefore, it is of primary importance to avoid singularities in a given workspace [27].
For this reason, we analyze the workspace that is calculated through the inverse kinematic explained above, and the singularity is determined using Equation (9), where the determinant of the Jacobian matrix is zero. At that point, there is a singularity.
  det ( J )   = 0 ,
Figure 3 shows the workspace (red surface) and the singularities (yellow surface) of the parallel mechanism. We can see that this mechanism has very few singularities, because the singularities are almost entirely placed outside the workspace.

2.3. Dynamics

The dynamics of this parallel mechanism were solved by means of the virtual work principle [28,29] in which the equation that describes the dynamics of a parallel mechanism is the following [30]
  J T τ + F p + i = 1   2 ( J i 1 T   F i 1 + J i 2 T   F i 2 ) + J 3 T   F 3   = τ a   ,
where J is the Jacobian matrix of the parallel mechanism, τ is a vector of forces and torques of the actuators, τ a is a vector of torques of Joint 1, F P is the force of platform, J i 1 T is the Jacobian matrix of the link i of Actuator 1, F i 1 is the actuator cylinder force i, J i 2 T is the Jacobian matrix of the link i of Actuator 2, F i 2 is the actuator piston force I, J 3 T is the Jacobian of the rotational actuator in Joint 2, and F 3 is the torque of the rotational actuator in Joint 2.

2.3.1. Resulting Force of Each Body Due to Buoyancy

The fluids generate a pushing force. The resultant force between the weight and that force [31,32] would be resolved using
  f m i   = m i   g ρ g   v i   ,
where f m i is the resulting force due to buoyancy, m i is the mass of the body i, g is the gravity vector, ρ is the density of the fluid, and v i is the volume of the body i. Clearing mass and gravity, we get
  f m i   = m i g ( 1 ( ρ   ρ i   ) )     ,
where ρ i is the density of the body i. We can see that the relationship of densities is a constant. Hence, to simplify the equation, it can be replaced by a constant λ i in such a way that the resulting force would be
  f m i   = m i g   λ i ,

2.3.2. Matrix of Forces and Torques of Each Actuator

The force of actuators is calculated with the spring displacement of each actuator.
  τ = | f 1   f 2 τ m | ,
where f 1 is the force of actuator L1, f 2 is the force of actuator L2, and τ m is the torque of actuator M.
These forces are calculated as
f 1 = Δ L 1   k L Δ v 3 , 1   B L ,
  f 2 = Δ L 2   k L Δ v 3 , 2   B L ,
  f 2 = Δ L 2   k L Δ v 3 , 2   B L ,
where k is the stiffness coefficient and B is the damping coefficient of the springs.

2.3.3. Platform of the Parallel Mechanism

Gravity and buoyancy are the forces of action, and the degrees-of-freedom of the platform are angular. The wrench will only have torque and not force, so this is resolved by
  F p = [ F n ] = [ 0 I c p w ˙ p w p x I c p w p + C   x   m p x   g λ p ] ,
where C is the centroid of the platform, m p is the mass of the platform, I c p the inertia of the platform, w p is the angular velocity of the platform, g is the gravity vector, and λ p is the relationship between the densities of the platform and fluid (see Equation (13)).

2.3.4. Cylinder of Linear Actuator i

The cylinder is affected by gravity and its linear acceleration.
  F i 1   = [ F n ] = [ m i 1 ( g λ i 1 a i 1 ) I c i 1 w ˙ i w i x I c i 1 w i ] ,
where m i 1 is the mass of the cylinder, a i 1 is the linear acceleration of the cylinder, w i is the angular velocity of the actuator Li in the lower joint (4 or 6), I c i 1 is its inertia, and λ i 1 is the relationship between the densities of the cylinder and fluid (see Equation (13)).

2.3.5. Piston of the Linear Actuator i

In the same manner, the following is solved by
F i 2 = [ F n ] = [ m i 2 ( g λ i 2 a i 2 ) I c i 2 w ˙ i w i x I c i 2 w i ] ,
where m i 2 is the mass of the piston, a i 2 is the linear acceleration of the piston, I c i 2 is its inertia, and λ i 2 is the relationship between the densities of the piston and fluid (see Equation (13)).

2.3.6. Rotating Actuator

This joint is only affected by gravity, buoyancy, and angular velocity,
F 3 = [ F n   ] = [ m 3 g λ 3 I c 2 w ˙ 3 w 3 x I c 2 w 3 ] ,
where   m 3 is the mass of the rotation actuator, w ˙ 3 is its angular acceleration, w 3 is its angular velocity, I c 2 is its inertia, and the λ i 2 is the relationship between the density of the piston and fluid (see Equation (13)).

2.3.7. Velocities and Angular Accelerations

The angular velocities and accelerations of the cylinder are calculated, as well as the piston and its linear accelerations. The vectors of the parallel mechanism will be defined first.
The unitary vectors of each actuator are
  S 1   = [ ( p 3 x p 4 x )   , ( p 3 y p 4 y )   ,   ( p 3 z p 4 z ) ] T | [ ( p 3 x p 4 x ) , ( p 3 y p 4 y )   , ( p 3 z p 4 z ) ] T |   ,
S 2   = [ ( p 5 x p 6 x )   ,   ( p 5 y p 6 y )   ,     ( p 5 z p 6 z ) ] T | [ ( p 5 x p 6 x )     ,     ( p 5 y p 6 y )   ,     ( p 5 z p 6 z ) ] T | ,
  S 3   = [ 0 , 0 , 1 ] T   ,
where S 1 is the unit vector of the actuator L1, S 2 of the unit vector actuator L2, and S 3 is the unit vector of actuator M in the center of the frame of reference in Joint 1.
Then, the vectors of the platform are
a 1 = [ ( p 3 x p 1 x )   , ( p 3 y p 1 y )   ,   ( p 3 z p 1 z ) ] ,
  a 2   = [ ( p 5 x p 1 x )   , ( p 5 y p 1 y )   ,   ( p 5 z p 1 z ) ]   ,
h = [ ( p 2 x p 1 x )   , ( p 2 y p 1 y )   ,   ( p 2 z p 1 z ) ] ,
Now, the angular speeds of the actuators are calculated.
In the case of actuator L1, the tangential speed [33] of actuator L1 and that of the platform in Joint 3 are equalized. Therefore, we have Equation (28):
w 1 x   L 1   = w p   x   a 1   ,
Solving this, we get
  w 1 x   L 1 x   S 1   = w p   x   a 1   ,
  w 1 = ( 1 L 1   )   S 1 x   w p   x   a 1   ,
Thus, Equation (30) describes the angular velocity of actuator L1 in Joint 4. In the same manner, actuator L2 is solved and obtained
  w 2 = ( 1 L 2   )   S 2 x   w p   x     a 2   ,
The angular velocity of the motor would be given simply by the following equation.
  w 3 = S 3   x   w p   ,  
Next, the angular accelerations of each actuator in its lower joints (4 or 6) will be calculated in the same way as Equation (28). However, matching the tangential accelerations, we see that the angular acceleration in actuator L1 would be
  w 1   ˙ = ( 1 L 1 ) x   S 1 x   a c 3 2   L 1 ˙ w 1   ,
where a c 3 is the acceleration in Joint 3 which is determined by
  a c 3 = w ˙ p   x   a 1   + w p x ( w p x a 1 ) ,
In the same way, the acceleration is solved for actuator L2 as
  w 2   ˙ = ( 1 L 2 ) x   S 2 x   a c 5 2   L 1 ˙ w 2   ,
where a c 3 is the acceleration in Joint 3 which is determined by
  a c 5 = w ˙ p   x   a 2   + w p x ( w p x a 2 )   ,
and the angular acceleration of the motor is only determined by
  w 3   ˙ = d w 3 d t     ,

2.3.8. Jacobian Matrix of Each Actuator

The velocity of the center of mass of the cylinder would be given by
  v i 1   = c i 1 ( w i   x   S i )   ,
  v i 1   = C i L i [ S i 2 S i 2 a i ] X ˙ ,
where C i is the center of mass of the cylinder and X ˙ is the velocity of the final effector. Now, we solve for the piston using
  v i 2   = c i 2 ( w i   x   S i ) + L i ˙ S i   ,  
  v i 2   = [ c i 2 L i S i 2   + S i   S i   T c i 2 L i   x   S i 2 x   a i S i   S i   T   x   a i   ] ,
Finally, the angular velocity of the actuators with respect to that of the effector would be given by
  w i = 1 L i   ( S i   x   V p + S i x ( w p x   a i ) ) ,
  w i = 1 L i   [     S i       S i x     a i   ]   [ V p W p ]   ,
where V p is the speed of the end effector of the parallel mechanism. In this case, as it is angular, this speed is zero, and W p is the angular velocity of the end effector. Therefore, the Jacobian matrix of the cylinder is
  J i 1   = 1 L i   [ C i S i 2 C i S i 2 a i S i   S i x   a i ] ,
and the Jacobian matrix of the piston is
  J i 2   = 1 L i   [ c i 2 S i 2   + S i   S i   T   c i 2 x   S i 2 x   a i S i   S i   T   x   a i S i   S i x   a i ] ,
and the Jacobian axis of the rotational actuator, as it is directly attached to the effector of the parallel mechanism, would be
  J 3 = [ S 3   ]   ,

2.3.9. Torque Generated by Fluid Dynamic Forces in the Sphere

The fluid dynamic forces in a body are calculated as follows [34]
  M   v ˙ + C ( v )   v + D ( v ) v + G = F s ,
where F s is the force resulting from the fluid dynamic forces in a body, M is the is the inertia matrix (including added mass), C is the vector of the Coriolis and centripetal force, D is the damping effect, and G is the buoyancy effect.
As this case is a sphere, the equation is resolved as
( m s + m s   ) a s + f a + m s g λ s = F S ,
where m s is the mass of the sphere, m s is the added mass in the sphere by the fluid, a s is the acceleration vector of the sphere obtained by the IMU, f a is the vector of the drag force of the fluid in the sphere, and λ s is the relationship of the densities between the sphere and fluid (see Equation (13)).
In this case, the force is not considered to be Coriolis because it is assumed that the sensor is not in an external referential system that rotates, and the centripetal force is very small with respect to the drag force. This is because the distance of the axis of rotation of the parallel mechanism and the sphere is small. This allows us to obtain the drag force as the only dependent factor on the velocity of the fluid.
The added mass in the sphere by the fluid can be calculated as [35]
  m s = 0.5 ρ v s     ,
where ρ is the density of the fluid and v s is the volume of the sphere.
Finally, the torque generated in Joint 1 (see Figure 2) of the parallel mechanism in the frame of reference is determined by
  τ a = r   x   F S     ,  
where r is the vector of the center of the frame of reference to the sensor’s IMU (see Figure 2).
Since we are interested in finding the drag force, it will be necessary to invert r . For this reason, we arrange the equation as an antisymmetric matrix [36]. Finally, replacing Equation (10), the drag force is calculated by
  [ r ]   1 x     (   J T τ + F p + i = 1 2 ( J i 1 T   F i 1 + J i 2 T   F i 2 ) + J 3 T   F 2   ) m s g λ s ( m s + m s )   a s = f a   ,

2.3.10. Velocity and Direction of the Fluid

The velocity and direction of the fluid is possible to calculate with the equation of the drag force of the fluids [34]
  v = f a   C D 0.5 ρ A ,
where C D is the drag coefficient of the sphere, ρ is the fluid density, and A is the cross area of the sphere.
The drag coefficient of the sphere varies with the Reynolds number, as can be seen in Figure 4. This coefficient is constant when the Reynolds number is 103 Re, with a value of 0.47 . The researchers Milkailov and Silvia Freire [37] show different equations for calculating the drag coefficient. According to this criterion, the sphere is designed to keep the drag coefficient constant.
Finally, the direction of the fluid is calculated with the unit vector of the drag force ( f a )
  u ^ = f a   f a   ,

2.3.11. Dynamic Performance

Since the sensor is a parallel mechanism and has mass and springs, this will have a dynamic behavior of the second order, in such a way that the transference function of the second-order system is as [38]
  Y ( s )   = k p   w n 2 s 2 + 2 ζ w n s + w n 2   ,
where ζ is the damping ratio, w n is the natural frequency of the system, and k p is the gain. For good design criteria ζ is set to 0.7 because, in this condition, the system is more robust to the variations in the parameters of the plant or actuator. Thus, with this value, it is possible to calculate the damping of the parallel mechanism with the following equation [38]
  B = ζ   2   K m     ,
where B is the damping of the parallel mechanism, K is the stiffness of the parallel mechanism, and m is the mass of the mechanism. The relationship of the stiffness of the parallel mechanism with the stiffness (ki) of the actuators is related to the Jacobian matrix, as [27]
  K = J T k i   J   ,  
where k i = [ k L , k L , k m ] T .
The stiffness of the parallel mechanism is determined by the response time we want for the sensor. This is the settling time ( t s ) and, for a second-order system, the response is usually 2% and can be obtained using the equation [38]
  t s = 4 ζ w n ,
where [39]
  w n   = K I   ,
I is the rotational mass moment of inertia since the parallel mechanism has rotational grades of freedom. Hence, depending on the settling time that we want the sensor to work with, we will determine the natural frequency ( w n ) of the parallel mechanism and, with Equation (58), determine the stiffness of the parallel mechanism. Finally, with Equation (56), we calculate the stiffness (ki) of each actuator.

2.4. Practical Implementation and Testing of the Sensor

2.4.1. Algorithm for the Sensor

Figure 5 shows the algorithm for calculating the velocity and direction of a fluid through a parallel mechanism, using the model mathematics described in this article. The algorithm begins by entering the constants, such as inertias, masses, spring proprieties, geometry, etc. Next, the values of the angles and acceleration measured by the IMU are obtained. Then, we calculate the displacement of the spring of each passive actuator with the equation of the inverse kinematic (see Equations (2)–(4)). We also calculate the Jacobian matrix (see Equation (8)) and, with the Jacobian matrix, we calculate the velocity of each actuator. Finally, with this data, it is possible to calculate the force and torque of the actuators (see Equation (14)).
Next, the algorithm calculates the drag force with the equation dynamic (see Equation (51)). Note that a part of the equation has been shown in blue; this part represents the dynamic forces due to the inertia of the parallel mechanism and, if the inertia of the mechanism is very low with respect to the drag force, it is possible not to consider this part of the equation because the forces that will generate the inertia of the mechanism will be very low. However, this should be done only if it is necessary to reduce the computational cost.
After obtaining the vector of forces of the drag force, the velocity is calculated with the hydrodynamic equation (see Equation (52)) and the direction of the fluid (see Equation (53)). Finally, with this result, we obtain the velocity and the direction of the fluid only with the acquisition of data on the angles and accelerations of the parallel mechanism through the IMU sensor.

2.4.2. Experiment

Testing the algorithm employed a co-simulation of MATLAB and Adams software. The algorithm was written in MATLAB, and the sensor of the parallel mechanism was designed in Adams software. The fluid considered in this simulation is water with a density ρ = 1000   kg / m 3 . Figure 6 shows the simulation process. Firstly, the velocity of the water is input. In this case, a velocity of 1 m/s in the x-axis was input. With sinusoidal variation in the model, Adams transformed this velocity into the drag force applied in the sphere.
This force deforms the structure of the parallel mechanism, and a virtual inertial sensor within the sphere measures the angles and accelerations, simulating the IMU sensor. These angles and accelerations were input into the algorithm in MATLAB, and this calculated the velocity and direction of the water. These results were compared with the velocity of the input. In this experiment, a test was also performed to understand the dynamic behavior of the sensor in terms of frequency response and measuring range.
Table 1 shows the physical characteristics of this sensor for the test. The springs’ constants have been calculated for a settling time of 0.15 s of the response of the sensor, as explained in the section on dynamic performance.

3. Results

The simulation time was 70 s. Figure 7 shows a comparison between the input velocity of the water (red line) and the velocity of the water calculated with the algorithm (blue line). This figure shows the results for 70 s (Figure 7a) and in the first 0.2 s (Figure 7b). Figure 8 shows the direction of the water calculated by the algorithm. The triangle is a reference to the geometry of the sensor. Figure 9 shows the response of the sensor to the frequency to determine the bandwidth of this sensor. Finally, Figure 10 shows the response of the sensor to several velocities. In this case, the velocity of the water that increases in linear form (red line) was input. The blue line represents the velocity of the water calculated by the algorithm. The determinant of the Jacobian matrix is shown in Figure 9b to better understand the singularity of the parallel mechanism.

4. Discussion

The results of Figure 7 show that the algorithm can calculate the velocity of the fluid through the angles and accelerations obtained by the sensor’s IMU. The settling time of the sensor is 0.15 s (Figure 7b). This value can increase or decrease, depending of the uses of the sensor. The stiffness of the parallel mechanism can be designed in such a way that the natural frequency of the mechanism increases or decreases the settling time (see Equation (57)).
On the other hand, Figure 8 shows the direction of water calculated by the algorithm. This is calculated to match the direction of the water that was input, as has been previously explained. The water has a direction in the x-axis. This is possible to calculate because the drag force has the same direction as the water, so the unitary vector of the drag force and unitary vector of the velocity water are equal.
Figure 9 shows the response of the frequency of the sensor. In this figure, we can see that the sensor response is good until a frequency of 40.6 rad/s. Above this value, the response of the sensor increases and the phase change to 180° degrees. This is because the parallel mechanism comes into resonance. Since the natural frequency of this mechanism is 55.6 rad/s, it is possible calculate the natural frequency with Equation (58). To graph the response of the frequency of the sensor, the transfer function was determined as
  Y ( s )   = 1.0592   ( 55.6 ) 2 s 2 + 2 ( 0.7 ) ( 55.6 ) s + ( 55.6 ) 2   ,
The range of measurement of the sensor can be seen in Figure 10a. Due to inertia and the stiffness of the parallel mechanism, the sensor measures correctly until a velocity of 3.41 m/s. This is because, at higher velocities, the drag force increases until it destroys the sensor. As can be seen in Figure 10b, the determinant of the Jacobian matrix (see Equation (8)) is zero at 61 s. This means that the parallel mechanism entered into a singularity, failing the sensor, as shown in Figure 10a, at 61 s. It is possible to increase this breaking point if we increase the stiffness of the parallel mechanism and also increase the natural frequency, hence decreasing the settling time. The main idea is to maintain the mechanism in its workspace because, as we have seen in Figure 3, there are no singularities that have an effect within the workspace.

5. Patents

Gonzalo Emmanuel Ejarque Rinaldini, Roque Jacinto Saltaren Pazmiño, Gabriel Armando Poletti Ruiz, and Rafael Aracil Santonja. “Dispositivo y método para la medición de corrientes de fluidos mediante mecanismo paralelo esférico actuado por fuerzas de arrastre.” ES. Patent no. ES2525394B2. 22 December 2014.

Author Contributions

R.S. contributed to the invention of the sensor, supervision, project administration, and review of the article. G.P. contributed to the development of the mathematical model, algorithm, simulation using the software, and writing of the article. A.R.B. and J.C. contributed to data collection, analysis, and review of the article.

Funding

The research leading to these results has received funding from the Spanish Government CICYT project ref. DPI2014-57220-C2-1-P, DPI2013-49527-EXP, the Universidad Politécnica de Madrid project Ref. AL14-PID-15, and the RoboCity2030-III-CM project (Robótica aplicada a la mejora de la calidad de vida de los ciudadanos. Fase III; S2013/MIT-2748). It was funded by Programas de Actividades I+D en la Comunidad de Madrid and co-funded by structural funds of the EU.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Neauman, G. Ocean Currents; Elsevier: New York, NY, USA, 1968. [Google Scholar]
  2. Merlet, J.P. Parallel Robots, 2nd ed.; Springer: Dordrecht, The Nertherlands, 2006. [Google Scholar]
  3. Gough, V.E.; Whitehall, S.G. Universal Tyre Testing Machine; International Technical Congress FISITA; Institution of Mechanical Engineers: London, UK, 1961. [Google Scholar]
  4. Stewart, D. Platform with six degrees of freedom. Proc. Inst. Mech. Eng. 1965, 180, 371–386. [Google Scholar] [CrossRef]
  5. Clavel, R. A fast robot with parallel geometry. In Proceedings of the 18th International Symposium on Industrial Robots, Lausanne, France, 26–28 April 1988. [Google Scholar]
  6. Gosselin, C.M.; Angeles, J. The optimum kinematic design of a spherical three-degree-of-freedom parallel manipulator. J. Mech. Transm. Autom. Des. 1989, 111, 202–207. [Google Scholar] [CrossRef]
  7. Lande, M.A.; David, R.J. Articulation for manipulator arm. U.S. Patent 4300362, 17 November 1981. [Google Scholar]
  8. Lu, Y.; Zhang, X.; Sui, C.; Han, J.; Hu, B. Kinematics/statics and workspace analysis of a 3-leg 5-DoF parallel manipulator with a UPU-type composite active constrained leg. Robotica 2013, 31, 183–191. [Google Scholar] [CrossRef]
  9. Saltaren, R.J.; Sabater, J.M.; Yime, E.; Azorin, J.M.; Aracil, R.; Garcia, N. Performance evaluation of spherical parallel platforms for humanoid robots. Robotica 2006, 25, 257–267. [Google Scholar] [CrossRef]
  10. Hao, G.; Murphy, M.; Luo, X. Development of a compliant-mechanism-based compact three-axis force sensor for high-precision manufacturing. In Proceedings of the International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, Boston, MA, USA, 2–5 August 2015; American Society of Mechanical Engineers: New York, NY, USA, 2015. [Google Scholar]
  11. Hao, G.; Kong, X. A normalization-based approach to the mobility analysis of spatial compliant multi-beam modules. Mech. Mach. Theory 2013, 59, 1–9. [Google Scholar] [CrossRef]
  12. Yang, E. Design and Sensitivity Analysis Simulation of a Novel 3D Force Sensor Based on a Parallel Mechanism. Sensors 2016, 16, 2147. [Google Scholar] [CrossRef] [PubMed]
  13. Yao, J.T.; Hou, Y.L.; Chen, J.; Lu, L.; Zhao, Y.S. Theoretical analysis and experiment research of a statically indeterminate pre-stressed six-axis force sensor. Sens. Actuators A Phys. 2009, 150, 1–11. [Google Scholar] [CrossRef]
  14. Lu, Y.; Wang, Y.; Ye, N.; Chen, L. Development of a novel sensor for hybrid hand with three fingers and analysis of its measure performances. Mech. Syst. Signal Process. 2017, 83, 116–129. [Google Scholar] [CrossRef]
  15. Álvarez, C.; Saltaren, R.; Aracil, R.; García, C. Concepción, Desarrollo y Avances en el Control de Navegación de Robots Submarinos Paralelos: El Robot Remo-I. Rev. Iberoam. Autom. Inf. Ind. 2009, 6, 92–100. [Google Scholar] [CrossRef]
  16. Rundtop, R.; Frank, K. Experimental evaluation of hydroacoustic instruments for ROV navigation along aquaculture net pens. Aquac. Eng. 2016, 74, 143–156. [Google Scholar] [CrossRef]
  17. Zhao, B.; Blanke, M.; Skjetne, R. Particle filter ROV navigation using hydroacoustic position and speed log measurements. In Proceedings of the American Control Conference (ACC), Montreal, QC, Canada, 27–29 June 2012. [Google Scholar]
  18. Liu, P.; Wang, B.; Deng, Z.; Fu, M. INS/DVL/PS Tightly Coupled Underwater Navigation Method with Limited DVL Measurements. IEEE Sens. J. 2018, 18, 2994–3002. [Google Scholar] [CrossRef]
  19. Fu, K.; Gonzalez, R.; Lee, C. ROBOTICS: Control, Snsing, Visión, and Intelligence, 1st ed.; McGraw-Hill Education: New York, NY, USA, 1989. [Google Scholar]
  20. Wen, L. Robot Analisis, the Mechanics of Serial and Parallel Manipulators; John Wiley & Sons: New York, NY, USA, 1999. [Google Scholar]
  21. Davidson, J.; Hunt, K. Robots and Screw Theory; Oxford University Press: Oxford, UK, 2002. [Google Scholar]
  22. Zhao, J.; Feng, Z.; Chu, F. Advanced Theory of Constraint and Motion Analysis for Robots Mechanisms; Academic Press: Cambridge, MA, USA, 2014. [Google Scholar]
  23. An, H.S.; Lee, J.H.; Lee, C.; Seo, T.W.; Lee, J.W. Geometrical kinematic solution of serial spatial manipulators using screw theory. Mech. Mach. Theory 2017, 116, 404–418. [Google Scholar] [CrossRef]
  24. Saltarén, R.; Puglisi, L.J.; Sabater, J.M.; Yime, E. Thomás. In Robótica Aplicada; Dextra Editorial: Madrid, Spain, 2017; pp. 314–315. [Google Scholar]
  25. Alvarado, J.G.; Agundis, A.R.; Garduño, H.R.; Ramírez, B.A. Kinematics of an asymmetrical three-legged parallel manipulator by means of the screw theory. Mech. Mach. Theory 2010, 45, 1013–1023. [Google Scholar] [CrossRef]
  26. Alvarado, J.G.; Martínez, M.R.; Alici, G. Kinematics and singularity analyses of a 4-dof parallel manipulator using screw theory. Mech. Mach. Theory 2006, 41, 1048–1061. [Google Scholar] [CrossRef]
  27. Taghirad, H.D. Parallel Robots: Mechanics and Control; CRC Press: New York, NY, USA, 2013. [Google Scholar]
  28. Franz, Z. Mechanics of Solids and Fluids; Springer: New York, NY, USA, 2012; pp. 225–256. [Google Scholar]
  29. Antman, S.S.; Osborn, J.E. The principle of virtual work and integral laws of motion. Arch. Ration. Mech. Anal. 1979, 69, 231–262. [Google Scholar] [CrossRef]
  30. Khalil, W.; Ibrahim, O. General solution for the dynamic modeling of parallel robots. J. Intell. Robot. Syst. 2007, 49, 19–37. [Google Scholar] [CrossRef]
  31. Mott, R.L. Mecánica de Fluidos Aplicada; Pearson Educación: Naucalpan de Juárez, Mexico, 1996. [Google Scholar]
  32. Agrawal, A.; Prasad, B.; Viswanathan, V.; Panda, S.K. Dynamic modeling of variable ballast tank for spherical underwater robot. In Proceedings of the International Conference on Industrial Technology, Cape Town, South Africa, 25–28 February 2013. [Google Scholar]
  33. Feastherstone, R. Rigid Body Dynamics Algorithms; Springer: Sydney, Autralia, 2008. [Google Scholar]
  34. Antonelli, G. Underwater Robots; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  35. Sumer, B.M.; Fredsoe, J. Hydrodynamics around Cylindrical Structures; World Scientific Publishing: Singapore, 1997. [Google Scholar]
  36. Kanatani, K. Understanding Geometric Algebra, Japan; CRC Press: New York, NY, USA, 2015. [Google Scholar]
  37. Mikhailov, M.D.; Silva Freire, A.P. The drag coefficient of a sphere: An approximation using Shanks transform. Powder Technol. 2013, 237, 432–435. [Google Scholar] [CrossRef]
  38. Ogata, K. Ingeniería de Control Moderna, 5th ed.; Pearson Educación: Madrid, Spain, 2010. [Google Scholar]
  39. Thomson, W.T. Theory of Vibration with Applications; CRC Press: New York, NY, USA, 1993. [Google Scholar]
Figure 1. This figure shows the sensor parts and the working method. (a) The sensor and its parts; in this figure, we see the sphere, the sensor’s inertial measurement unit (IMU), and the parallel mechanism. (b) The parts of the parallel mechanism. L1 and L2 are the prismatic passive actuators, and A1 is the rotational passive actuator. A2, A3, and B1 are universal joints and B2 and B3 are spherical joints. Finally point P is the center of the spherical motion.
Figure 1. This figure shows the sensor parts and the working method. (a) The sensor and its parts; in this figure, we see the sphere, the sensor’s inertial measurement unit (IMU), and the parallel mechanism. (b) The parts of the parallel mechanism. L1 and L2 are the prismatic passive actuators, and A1 is the rotational passive actuator. A2, A3, and B1 are universal joints and B2 and B3 are spherical joints. Finally point P is the center of the spherical motion.
Sensors 18 02867 g001
Figure 2. Simplified diagram of the parallel mechanism, where the plane triangles are isosceles and the line of length ‘e’ is the median of the triangle. The sphere is also represented, and within it is the IMU. F s is the force resulting from the fluid dynamic forces, and r is the vector of the center of the coordinates to the IMU.
Figure 2. Simplified diagram of the parallel mechanism, where the plane triangles are isosceles and the line of length ‘e’ is the median of the triangle. The sphere is also represented, and within it is the IMU. F s is the force resulting from the fluid dynamic forces, and r is the vector of the center of the coordinates to the IMU.
Sensors 18 02867 g002
Figure 3. Workspace (red surface) and singularities (yellow surface) of the parallel mechanism. In this figure, almost all the of the singularity surface is placed outside the workspace. There are very few singularities in the possible rotation of the parallel mechanism.
Figure 3. Workspace (red surface) and singularities (yellow surface) of the parallel mechanism. In this figure, almost all the of the singularity surface is placed outside the workspace. There are very few singularities in the possible rotation of the parallel mechanism.
Sensors 18 02867 g003
Figure 4. Drag coefficient (CD) of a sphere versus the Reynolds number (Re). We can see that the drag coefficient is constant, starting from 103 Re; with an approximate value of 0.47 [37].
Figure 4. Drag coefficient (CD) of a sphere versus the Reynolds number (Re). We can see that the drag coefficient is constant, starting from 103 Re; with an approximate value of 0.47 [37].
Sensors 18 02867 g004
Figure 5. Algorithm for calculating the velocity and direction of the fluid with the parallel mechanism.
Figure 5. Algorithm for calculating the velocity and direction of the fluid with the parallel mechanism.
Sensors 18 02867 g005
Figure 6. This figure shows the method of simulation for testing the algorithm. The sensor is designed in the Adams software, where all the effects of the dynamics and hydrodynamics of the mechanism are simulated. In this Adams model, the water velocity is input, and the model in Adams software transforms this velocity into drag forces applied in the sphere. The model has a sensor in the sphere that calculates the angles and accelerations as an IMU sensor. These angles and accelerations are input into the algorithm, and the velocity and direction of the water are calculated.
Figure 6. This figure shows the method of simulation for testing the algorithm. The sensor is designed in the Adams software, where all the effects of the dynamics and hydrodynamics of the mechanism are simulated. In this Adams model, the water velocity is input, and the model in Adams software transforms this velocity into drag forces applied in the sphere. The model has a sensor in the sphere that calculates the angles and accelerations as an IMU sensor. These angles and accelerations are input into the algorithm, and the velocity and direction of the water are calculated.
Sensors 18 02867 g006
Figure 7. The results of the algorithm. The red line is the input velocity of the water and the blue line is the velocity of the water calculated with the algorithm. (a) The results for 70 s. In the first seconds, the sensor measures a little late; it can be seen in (b) that the settling time is 0.15 s. Then, as seen in (a), the sensor calculates the velocity correctly.
Figure 7. The results of the algorithm. The red line is the input velocity of the water and the blue line is the velocity of the water calculated with the algorithm. (a) The results for 70 s. In the first seconds, the sensor measures a little late; it can be seen in (b) that the settling time is 0.15 s. Then, as seen in (a), the sensor calculates the velocity correctly.
Sensors 18 02867 g007
Figure 8. The results of the direction of the water calculated by the algorithm. This vector matches the water velocity vector. The triangle represents the parallel mechanism.
Figure 8. The results of the direction of the water calculated by the algorithm. This vector matches the water velocity vector. The triangle represents the parallel mechanism.
Sensors 18 02867 g008
Figure 9. This figure shows the response dynamics of the sensor in terms of frequency. We can see that the sensor can measure correctly up to a lower frequency of 40.6 rad/s. Up to a frequency of 55.6 rad/s the system enters into resonance, as this value is the natural frequency of the parallel mechanism.
Figure 9. This figure shows the response dynamics of the sensor in terms of frequency. We can see that the sensor can measure correctly up to a lower frequency of 40.6 rad/s. Up to a frequency of 55.6 rad/s the system enters into resonance, as this value is the natural frequency of the parallel mechanism.
Sensors 18 02867 g009
Figure 10. (a) The measure range of the sensor. The red line is the input velocity of the water. In this case, the velocity is increasing in a linear form. The blue line is the result of the algorithm. We can see that the sensor response is up to 3.41 m/s. (b) The determinant of the Jacobian matrix (see Equation (8)). We can see that thr is zero at 61 s. This means that the parallel mechanism enters into a singularity, failing the sensor, as shown in (a), at 61 s.
Figure 10. (a) The measure range of the sensor. The red line is the input velocity of the water. In this case, the velocity is increasing in a linear form. The blue line is the result of the algorithm. We can see that the sensor response is up to 3.41 m/s. (b) The determinant of the Jacobian matrix (see Equation (8)). We can see that thr is zero at 61 s. This means that the parallel mechanism enters into a singularity, failing the sensor, as shown in (a), at 61 s.
Sensors 18 02867 g010
Table 1. Physical characteristics of the sensor flow.
Table 1. Physical characteristics of the sensor flow.
Geometry (m)Masses (kg)Spring Constants
L1 = 0.250
L2 = 0.250
h = 0.25
h = 0.250
b = 0.20
a = 0.40
Diameter of the sphere = 0.1 m
Triangular platforms = 0.1
L1 = 0.2
L2 = 0.2
h = 0.01
h = 0.01
Sphere = 0.001
k L = 320 N/m
B L = 8 (N.s)/m
k m = 18.33 (N.m)/°
B m = 0.458 (N.m.s)/°

Share and Cite

MDPI and ACS Style

Saltarén, R.; Portilla, G.; Barroso, A.R.; Cely, J. A Sensor Based on a Spherical Parallel Mechanism for the Measurement of Fluid Velocity: Physical Modelling and Computational Analysis. Sensors 2018, 18, 2867. https://doi.org/10.3390/s18092867

AMA Style

Saltarén R, Portilla G, Barroso AR, Cely J. A Sensor Based on a Spherical Parallel Mechanism for the Measurement of Fluid Velocity: Physical Modelling and Computational Analysis. Sensors. 2018; 18(9):2867. https://doi.org/10.3390/s18092867

Chicago/Turabian Style

Saltarén, Roque, Gerardo Portilla, Alejandro R. Barroso, and Juan Cely. 2018. "A Sensor Based on a Spherical Parallel Mechanism for the Measurement of Fluid Velocity: Physical Modelling and Computational Analysis" Sensors 18, no. 9: 2867. https://doi.org/10.3390/s18092867

APA Style

Saltarén, R., Portilla, G., Barroso, A. R., & Cely, J. (2018). A Sensor Based on a Spherical Parallel Mechanism for the Measurement of Fluid Velocity: Physical Modelling and Computational Analysis. Sensors, 18(9), 2867. https://doi.org/10.3390/s18092867

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop