ORIGINALARBEITEN/ORIGINALS https://doi.org/10.1007/s10010-022-00605-y Forsch Ingenieurwes (2023) :617–626 Navigation method utilizing floor-integrated inductive power supply modules for omnidirectional AGVs Carolin Carmen Brenner1 · Alexander Enssle2 · Robert Schulz1 · Nejila Parspour2 Received: 20 January 2022 / Accepted: 14 November 2022 / Published online: 8 December 2022 © The Author(s) 2022 Abstract This paper answers the research question: Can the contactless induced energy supply from a novel inductive floor be used to navigate omnidirectional automated guided vehicles (AGVs)? In contrast to existing systems a novel inductive floor enables AGVs traveling through production without charging breaks. This floor consists of tiles with inductive modules, which supply the AGV with energy. In addition to supplying power to the AGV, the inductive modules are also intended to guide the vehicle through production. To enable such a guidance sensors placed in the AGV measure the induced voltages of the floor. To answer the research question these voltages are calculated with the help of an electromagnetic simulation of the AGV’s travel on the inductive tiles. To estimate the position as well as rotation of the AGV depending on the simulated voltages as inputs a novel algorithm is presented. During the travel the AGV is able to move in arbitrary directions independently of its orientation. To control the omnidrectional AGV consistently without singularities, a transformation in Omni-Curve-Parameters (OCP) is proposed. As simulation case study a four wheeled steering- and velocity controlled AGV is introduced. For the evaluation a novel motion model depending on the input OCP is presented. This model is compared to the estimation of the position to verify the accuracy and the reproducibility of the algorithm. Methodik zur Navigation anhand bodenintegrierter induktiver Stromversorgungsmodule für omnidirektionale FTFs Zusammenfassung Dieser Beitrag beantwortet die Forschungsfrage: Kann anhand der kontaktlosen Energieversorgung eines neuartigen in- duktiven Bodens eine Navigation von fahrerlosen Transportfahrzeugen (FTFs) erfolgen? � Carolin Carmen Brenner carolin.brenner@ift.uni-stuttgart.de Alexander Enssle alexander.enssle@iew.uni-stuttgart.de Robert Schulz robert.schulz@ift.uni-stuttgart.de Nejila Parspour nejila.parspour@iew.uni-stuttgart.de 1 Institute of Mechanical Handling and Logistics, University of Stuttgart, Holzgartenstraße 15, 70174 Stuttgart, Germany 2 Institute of Electrical Energy Conversion, University of Stuttgart, Pfaffenwaldring 47, 70569 Stuttgart, Germany K https://doi.org/10.1007/s10010-022-00605-y http://crossmark.crossref.org/dialog/?doi=10.1007/s10010-022-00605-y&domain=pdf http://orcid.org/0000-0003-3972-5466 618 Forsch Ingenieurwes (2023) 87:617–626 Im Gegensatz zu bestehenden Systemen ermöglicht ein neuartiger induktiver Boden die Fahrt von FTF durch die Pro- duktion ohne Ladepausen. Dieser Boden besteht aus Platten mit induktiven Modulen, die das FTF mit Energie versorgen. Neben der Versorgung des FTFs sollen die induktiven Module das Fahrzeug auch durch die Produktion führen. Um ein Leiten des FTF durch den Boden umzusetzen, messen im FTF platzierte Sensoren die induzierten Spannungen. Zur Be- antwortung der Forschungsfrage werden diese Spannungen mit Hilfe einer elektromagnetischen Simulation der Fahrt des FTFs auf den induktiven Platten berechnet. Zur Schätzung der Position, so wie der Rotation des FTF wird ein neuartiger Algorithmus anhand dieser simulierten Spannungen als Eingangsgrößen vorgestellt. Während der Fahrt kann sich das FTF unabhängig von seiner Orientierung in beliebige Richtungen bewegen. Um das FTF konsistent ohne Singularitäten zu führen, wird eine Transformation in Omni-Kurven-Parametern (OCP) dargelegt. Als Fallbeispiel wird ein vierrädriges lenk- und geschwindigkeitsgesteuertes FTF vorgestellt. Für die Bewertung der Kursbestimmung und der Spurführung wird ein neuartiges Bewegungsmodell in Abhängigkeit von den OCP vorgeschlagen. Dieses Modell wird mit der Schätzung der Position verglichen, um die Genauigkeit und Reproduzierbarkeit des Algorithmus zu verifizieren. 1 Introduction Automated guided vehicles (AGVs) are an essential com- ponent of flexible and adaptable production of the future. To enable flexible and modular production systems, AGVs offer great advantages due to their maneuverability: In con- trast to conveyor technologies, they are able to serve differ- ent stations adaptably, which results in a flexible, versatile and freely configurable production [7, 11, 21]. To increase the operating time of AGVs a novel inductive floor enables AGVs to be dynamically supplied with energy during their travel. This paper addresses the research question: Can the contactless induced energy supply from a novel inductive floor be used to navigate omnidirectional AGVs? With such a navigation no extra tracking sensors are needed. The motion characteristics of AGVs can be distinguished in vehicles with two degrees of freedom (DOF) and omni- directional vehicles with three degrees of freedom. Due to their steering kinematics, AGVs with two DOF can travel along any target trajectory with a restriction regarding their chassis orientation. In contrast, vehicles with three DOF are able to move in every direction, independent of the chassis orientation [19]. Using three degrees of freedom enables for example smaller radii, which optimizes use of space [19]. Therefore maneuvering time is shorter, because omnidirec- tional AGVs can travel in the desired direction without first changing their orientation. This facilitates faster cycle times in production. Due to the high cost of warehouse facilities, optimizing space saves money [6, 12]. To utilize the maneuverability of the omnidirectional AGV a navigation which leads the AGV through production is needed. In general, navigation of an AGV solves three important challenges [5]: � Tracking the position of the AGV; � predicting where the current path will lead it without cor- rective measures; � deciding on necessary adjustments to safely reach a tar- get, possibly along a specified path. While localization describes the process of determining an object in space, tracking has the goal to identify the position of an object over time [9]. To navigate the AGV through a production layout, tracking can be realized using different kinds of markers in the environment. A grid of discrete optical or inductive floor marks can be formed (distance between 1m and 10m), providing ab- solute positions along a free-form track [20]. Simultaneous localization and mapping (SLAM) algorithms and Monte Carlo localization as introduced in [8], are additionally used for the navigation of the AGV between the marks, including position and orientation. Laser, ultrasonic, radar or camera systems are used, to detect the position and direc- tion of the vehicle. These systems use off-the-floor fixing marks as reflectors on walls or pillars to deduce the posi- tion. Instead of these fixing marks, contours are also used as natural marks [18]. The marks have to be unambigu- ous and sufficiently numerous, which limits the flexibility of the production using these systems [20]. A millimeter- scale accuracy is possible [22]. Active transmitting marks include radio transmitters as well as optical transmitters. This method includes global positioning system (GPS) out- side buildings and stationary miniature transmitters inside buildings. Navigation using active transmitting marks al- lows extended ranges [20]. The system is accurate in the order of few centimetres [13] and facilitates as the other introduced systems omnidirectional movement. Physical guidance systems often use optical, inductive or mechanical lanes. In these systems a primary circuit con- sists of one coil along the track of the AGV [18]. One or more coils are attached to the underside of the vehicle to detect this track. A deviation from this track is detected by the voltage difference of the coils in the vehicle [4]. Thereby the path is corrected accordingly in the vehicle’s control system. This physical lane, which is installed in the floor, restricts the path of the AGV. Consequently it is only able to move along a fixed line. Hence transport routes must be defined before the tracks are installed in the floor. Subsequent changes to the routes need structural changes K Forsch Ingenieurwes (2023) 87:617–626 619 to the track in the floor, which makes a production layout inflexible. In order to fulfill a successful navigation along a desired path, it is also important to decide on necessary adjustments. Odometry is one of the most important techniques to cope with such problems. Encoder data is used to track the mo- tion progress from a certain starting position. Due to wheel slip, ground roughness, etc., errors occur [10]. Therefore, to use odometry for localization, a fusion is made with, for example, the sensor data of the navigation techniques mentioned above. All of these introduced navigation systems separate the navigation and energy supply: The AGV additionally has to be equipped with a battery. Such a battery can be charged wirelessly, however, the wireless charging module is fixed in place and must be reached by the AGV before energy can be supplied. A novel inductive raised floor system supplies AGVs with energy at any given point, independent of fixed charg- ing spots. It is intended to make production equipment posi- tioning and supply efficient [14] and thus enables flexible, adaptable production. For this purpose, a concept is de- scribed in [15], which enables the supply of energy by an inductive system for mobile production and transport sys- tems. The individual modules, each installed in a raised floor panel, are able to provide their energy selectively. The floor consists of individual tiles with up to five modules in each tile. By moving the energy supply point, in which individual modules provide the energy selectively, the con- tinuous transfer of energy to an AGV is ensured. Thus, the AGV is guided omnidirectional through the production and logistics areas. Simultaneously it is supplied with energy by individual modules under the AGV. The prototype of such an omnidirectional AGV on the inductive floor is shown in Fig. 1. The intelligent floor was conceptualized and individual tiles were manufactured. Fig. 1 Study case: Scooty Inductive. Scooty is an omnidirectional AGV equipped with four wheels controllable in steering and velocity. It is guided with the help of an inductive floor. With the floor tiles, which include inductive modules, the AGV is able to reach every arbitrary position independent of its orientation. During its travel it is supplied with energy A whole floor full of inductive tiles will be manufactured in the near future. In this paper the induced voltages from this inductive floor are used to navigate omnidirectional automated guided vehicles (AGVs). More specifically, the main contributions are: � The proposal of the Omni-Curve-Parameters (OCP) for the control of an omnidirectional AGV. This is the first work, to the knowledge of the authors, that uses this Transformation to control an AGV to avoid singularities. � Depending on the induced voltages of the inductive floor tiles a finite element method (FEM) is used to simulate the movement of the AGV. A novel algorithm, which processes the measure data, combined with an odome- try model in a Kalman filter, that uses these simulated voltages is presented to estimate the position on the floor tiles. � A novel motion model depending on the input OCP is proposed. This is compared to the estimation of the posi- tion to verify the accuracy and the reproducibility of the algorithm. 2 Structure of Localization and Motion Control of the AGV with the Inductive Floor AGVs usually travel with their fronts centered along the track. Omnidirectional vehicles offer more opportunities to move through production. The omnidirectional AGV in Fig. 2 is capable of reaching every goal independent of its orientation. An arbitrary path is represented by the dashed red line. The vehicle is supposed to follow this illustrated path. For this purpose it starts, equivalent to the already existing inductive lanes, with the front side perpendicularly aligned to the track. In the second position, the vehicle con- Fig. 2 Exemplary movement of an omnidirectional AGV with an in- ductive floor. The given path of the AGV is marked dashed in red. During its travel along the path the omnidirectional AGV is able to ro- tate. In the first position it starts with the front centered along the track (rotation and slip angle eˇ = 0ı). In the second position the AGV ro- tated along the way and drives with abouteˇ = 35ı. In the third position it is following the path witheˇ = 90ı rotation K 620 Forsch Ingenieurwes (2023) 87:617–626 Fig. 3 Structure of the novel navigation method for the omnidirectional AGV with the inductive floor. A desired position of the vehicle is specified in OCP (ˇ; �n; vn). The angle ˇ0 is introduced as angular offset to the global coordinate system. With the evaluated speed utarget and orientation of the wheels ˛target, the specification of the vehicle model is simulated. The vehicle reference movement evaluates the movement of the AGV on the inductive floor. To determine the measurements xm and ym for the Kalman filter the voltages induced by the inductive floor and measured in the AGV (us ) are used as inputs for the in Sect. 5 introduced algorithm. The Kalman filter estimates xpos, ypos and rotation (rot) of the vehicle. The introduced concept can be integrated into the control loop outlined. In this case vn is not controlled for stability reasons tinues to follow the path, but it has turned by a desired angle during its travel. It follows the track with its front right cor- ner (eˇ = 35ı). As shown in this example, omnidirectional vehicles are able to smoothly position themselves in every desired orientation to the lane during their travel. The AGV used in the study case and shown in Fig. 1 consists of four traction and four steering drives. With these four wheels, controllable in speed u = .u1; u2; u3; u4/ T and orientation ˛= .˛1; ˛2; ˛3; ˛4/ T it is positioned on the in- ductive tiles in an omnidirectional manner. Four additional sensor coils are mounted in the vehicle, which Fig. 4 shows schematically, marked S1 to S4. These sensors are sym- metrically arranged around the center of the vehicle and determine the induced voltages of the floor’s magnetic field. The coordinate system on an inductive module used in the following is also mapped in Fig. 4. Various components are needed to track the vehicle with the inductive floor. An overview of the most important quantities are shown in Fig. 3. A desired position of the Fig. 4 Schematic representation of the vehicle on an inductive floor tile. The sensors for detecting the induced voltages integrated in the AGV are named S1–S4. The indicated coordinate system applies throughout all following calculations vehicle is specified in OCP and described in Sect. 3. With the evaluated speed and orientation of the wheels, the spec- ification of the vehicle model is simulated. To determine the measurements xm and ym for the Kalman filter a novel algorithm is introduced in Sect. 5. 3 Omni-Curve-Parameters as Universal Transformation for Omnidirectional AGVs Omni-Curve-Parameters presented in [2] are used to specify the omnidirectional movement of the AGV. The OCP describe the coordinated motion of a vehi- cle with an arbitrary number and types of wheels, using three physically independent parameters: nominal veloc- ity vn, nominal curvature �n and slip angle ˇ. With these OCP, a vehicle with arbitrary variably arranged wheels can be controlled intuitively in an omnidirectional manner. Two of these parameters (�n and ˇ) describe the motion modes and only influence the steering configuration of the vehicle wheels [2]. Therefore, the desired path and the orientation of the vehicle towards the path is described with the configuration parameters ˇ and �n. To rotate the AGV during travel, as described in Fig. 2, the angle ˇ0 is introduced. ˇ0 rotates the fixed vehicle co- ordinate system xAGV, yAGV counterclockwise, as shown in Fig. 5. The result is coordinate system x0 AGV, y0 AGV. The slip angle ˇ is specified according to the vehicle orientation in x0 AGV, y0 AGV plane towards the desired path. If ˇ0 is zero, the slip angle ˇ is specified according to the vehicle ori- entation in xAGV, yAGV plane towards the desired path. ˇ0 and ˇ adds up to eˇ, as shown in Fig. 2. K Forsch Ingenieurwes (2023) 87:617–626 621 Fig. 5 Specification of the OCP for the omnidirectional AGV. Nomi- nal curvature �n is used to describe the curvature of the path. Slip angle ˇ describes the orientation of the vehicle towards the path. A change of the rotation of the velocity vn is realized with a rotation of the xAGV, yAGV plane through an angle ˇ0. The OCP are then applied in the rotated coordinate system x0 AGV, y0 AGV In general, a path in a topological space (W) is a curve f W Œa; b� ! W with domain [a,b] (with a > b and a,b as real numbers). This path can be parameterized by a curve. E.g. a curve can be defined as image of the path [1]. There- fore the nominal curvature is used, to describe the path of the AGV without singularities. For example, in Fig. 5, the nominal curvature is chosen according to the curvature of the desired path marked in red. The last parameter, the nominal velocity is the desired velocity adapted to the AGV. Nominal velocity, is scaled by a maximum fixed, vehicle- specific velocity [2]. The major advantage of using OCP in the case of omnidi- rectional inductive guidance is the unambiguous specifica- tion of the steering target angles and target velocities. Con- sequently the AGV is controlled consistently and smoothly. Otherwise, in vehicles with many steered wheels, singu- larities occur due to the over-determinism, which leads to a high control effort. An example is given in [3], where chassis, which provide non-holonomic and omnidirectional motion capabilities, are kinematically modelled and con- trolled. Using OCP, each independently steered and indepen- dently driven motor is specified. The calculated target ve- locities and target steering angles are chosen to satisfy the condition that the pole rays of all wheels intersect in a particular point. This allows controlled cornering with- out singularities. Consequently, the overdetermined system is controllable in an omnidirectional manner without any computationally expensive conditions. The presented approach of steering the vehicle with OCP is already implemented, manually tested and verified on the presented vehicle. 4 Reference Vehicle Movement Generated reference values of the OCP are the input val- ues for the vehicle model (see Fig. 3). The AGV model is required for the simulation of the navigation of the AGV. In the simulation, it is assumed that the motors directly ap- proach the specified reference values. Therefore, the OCP are used for the vehicle reference model. Using the OCP as generalized model, it is applicable for any omnidirectional AGV. The transformation of the motor steering angles and the motor velocities back into the OCP is not unique. The following vehicle reference is used as test case with constant input �n. The relationship of the presented parameters is shown in Fig. 6. Driving a curve with curvature-radius R and angle � and slip angle ˇ = 0, as shown in Fig. 6 the velocity vAGV in the xpos, ypos plane (as shown in Fig. 4) on the given tile- coordinates can be determined: vAGV;ˇ=0 = 0 @ vxpos;ˇ=0 vypos;ˇ=0 1 1 A = 0 @ cos.�/vn sin.�/vn 1 1 A ; (1) with time t1: � = Z t1 0 ! dt: (2) For driving a curve with an arbitrary slip angle ˇ the vehi- cles position is calculated as described in (1) with a rotation around the instantaneous center of rotation (ICR): TAGV!ICR = 0 @ cos.ˇ/ − sin.ˇ/ 0 sin.ˇ/ cos.ˇ/ 0 0 0 1 1 A : (3) Fig. 6 Important Parameters to mathematically describe the motion of the AGV. R describes the radius of the path, the AGV is driving. This path describing a circle section is characterized with the angle � . The side-slip angle rotates the instantaneous center of rotation (ICR) coun- terclockwise. The velocity of the vehicle is named vAGV K 622 Forsch Ingenieurwes (2023) 87:617–626 The movement depending on ˇ is consequently expressed as: vAGV = TAGV!ICR � vAGV;ˇ=0: (4) The path rAGV with: rAGV = 0 @ xpos ypos 1 1 A (5) is calculated: rAGV = Z t1 0 vAGV dt: (6) In addition to the vehicle reference model, the vehicle sensor signals of the sensors detecting the voltages, are required (see Fig. 3). For this purpose data is gathered by an electromagnetic simulation. Based on the xpos- and ypos of the vehicle calculated in the vehicle movement, the sensor coil values are evaluated at the corresponding position and output from the vehicle model. Accordingly, the measured induced voltage values of the four sensors and the individual motor speeds and motor angles are the output parameters of the vehicle model. The mathematical approach of the OCP has already been manually tested and verified on the presented vehicle. 5 Estimation of the Vehicle Position with the Inductive Floor In Fig.7 an overview of the FEM-model is illustrated. The inductive floor consists of a cross-shaped arrangement of square coils. That way, the concept of the inductive floor can be modularly extended in x- and y-direction. The square- shaped secondary coil is located underneath the vehicle, Fig. 7 FEM-model of the inductive floor and the vehicle model. The secondary coil is included in the vehicle, the primary coil in the floor. The sense coils on the bottom side of the vehicle are used to calculate the position-dependent induced voltages. [16] which itself is modeled with an aluminum panel of the vehicle’s bottom extents. Ferrite cores on primary and sec- ondary side improve the flux guidance of the system. The four sense coils are symmetrically arranged in the corners of the secondary coil system. In principle, the excitation of primary and secondary power coils leads to a geometry- and thus position-depen- dent spatial magnetic field distribution. The sense coils are arranged in such way that they form a weak magnetic cou- pling with the secondary power coil. This is due to the fact that the secondary coil, is square shaped in a first approxi- mation. Following Ampere’s Law, the magnetic field caused by the coil’s current mainly occurs around its sides’ lengths. The sense coils are only slightly affeced by this particular magnetic field contribution. Consequently, the position of the vehicle in relation to the primary power coils dominates the effect on the sense coils. With the parametrized FEM-model, the coupling factors (or mutual inductances M ) of each sense coil with each power coil can be simulated for any arbitrary position of the vehicle on the inductive floor. The knowledge of the coupling factors (or mutual inductances M ) can be used to extract the position-dependent induced voltage in each sense coil using Faraday’s law: usn = −M di dt : (7) As shown in Fig.7 each tile module of the inductive floor has a symmetrical structure, as there is a square core in the centre. The voltage maps of the four sensors mounted in the vehicle, depending on the x- and y-position (as de- fined in Fig. 7) are shown in Fig. 8. The 3D plots indicate that the largest offset voltage is measured in the center of the inductive module. Toward the edge of the module the voltages are decreasing. The coordinate origin of the ve- hicle, as shown in Fig. 4, is on the bottom module in the simulation, shown in Fig. 8. Accordingly, the following ap- plies: If the vehicle center is exactly above the inductive modules center, the sensor signals of S1 and S3, as well as the measured voltages of S2 and S4 cancel each other out. On the other hand, if the AGV drives with a positive y-offset, the detection of a higher voltage at sensor coils S3 and S4 occurs. With a y-offset in the negative direction of the vehicle, a correspondingly higher voltage value occurs at sensor coils S1 and S2. For an AGV positioned as shown in Fig. 4 rotated with ˇ in positive direction, as in the first position shown in Fig. 2, the rotation can also be evalu- ated. In this case the sensor signal S4 is high because it is near the center of the module and the voltage of S2 is low, because it is far away from the center of the module. The sensors installed in the back show the opposite behavior as the front: S1 detects a high voltage and S3 a low voltage. Based on this symmetrical behaviour of the coils mounted K Forsch Ingenieurwes (2023) 87:617–626 623 Fig. 8 Simulation of the signals of the four sensors (S1–S4) installed in the AGV. The x and and y-position relates to the center of the AGV. Since the sensors are not centered in the origin of the AGV coordinate system, the voltages are not zero at this point. But in the origin the sense-signals cancel each other out. Mesh resolution depends on the resolution of the electromagnetic simulation in the vehicle, it is consequently possible to determine the position of the vehicle on the inductive floor modules. For the evaluation of the vehicle position on the basis of the coil signals on an inductive module, an algorithm is required to evaluate the xpos-position and the ypos-position of the AGV. Therefore every sensor coil k has its own evaluation matrix Uk s : Uk s = � uk s.i;j / � i=1;:::;gIj=1;2;3 (8) with g as number of simulation data points. The first column includes the induced voltages, the second column the inher- ent measured xm-position and the third column the inherent ym-position. xm and ym correspond to xpos and ym coordi- nate system, shown in Fig. 4. They serve as measurement variables for the Kalman filter, as illustrated in Fig. 3. In practice the coil signals of the S1–S4 sensors can diverge due to production accuracies. Therefore the intro- duced matrix Uk s of the particular sensors is replaced with experiment data. This ensures the practical feasibility of the introduced algorithm. With the measured voltage uk m at sensor k the four near- est neighbours with the smallest euclidean distance: d � uk s.i;1/; uk m � = v u u t g X i=1 � uk s.i;1/ − uk m �2 (9) are evaluated. In this presented implementation, the four nearest neighbours to the measured um;k value are searched for. However because of measurement inaccuracy and man- ufacturing tolerances in practice, the search for more neigh- bours may turn out to be useful and can be trivially ex- tended. Since the y-offset of the vehicle is equivalent at the front (S2 and S4), these two are used for evaluation. The same applies to the rear sensors (S1 and S3). Using the deter- mined neighbours nk1:::k4 of the measured uk m value, the er ror is evaluated. For this purpose, the error efront of the ne ighbours of S2 and S4 is calculated respectively. Based on the minimum of the error efront the position of the vehi- cle is determined. Subsequently, the rear position of the ve hicle is calculated. To evaluate the position, if the resolution of the Us;k matrix is high, more nearest neighbours can be searched and the expected value can be calculated. If more measuring points are used than in this simu- lation, searching the matrix can be very computationally expensive. Therefore, the described evaluation is supple- mented. A prediction of the next position value is possible with an adaption of the search space. Around the last esti- mated value uev the evaluation matrix Uk s is adapted: Uk s = .uk s.i;j //i=i.uev/;:::;i.uev+Ra/Ij=1;2;3 (10) The radius Ra of this downsizing process increases propor- tionally to the driven speed and can be adjusted accordingly. Initially, however, an initialization on the floor tile is needed to find the appropriate starting point. The evaluated xm- and ym-position is used as measure- ment combined with an odometry model of the AGV in a Kalman filter. The advantage of a Kalman filter is that noisy inaccurate measurement values and inaccuracies due to wheel slip and ground roughness etc. can be suppressed. K 624 Forsch Ingenieurwes (2023) 87:617–626 In this case, a trivial AGV model, depending on the mo- tor velocities u and on the steering angles ˛ is used for prediction. Estimate the x-position xpos as: Px = � 0 1 0 0 � x + � 0 bx ; � (11) where: Px = � xpos Pxpos � (12) with: bx = cos.˛1/ Pu1 + cos.˛2/ Pu2 + cos.˛3/ Pu3 + cos.˛4/ Pu4: (13) and xpos = � 1 1 � x: (14) For the y-position ypos equivalently: Py = � 0 1 0 0 � y + � 0 by � (15) where: Py = � ypos Pypos � (16) with: by = sin.˛1/ Pu1 + sin.˛2/ Pu2 + sin.˛3/ Pu3 + sin.˛4/ Pu4 (17) and ypos = � 1 1 � y: (18) This xpos and ypos are also used to estimate the orientation o of the vehicle within the coordinate system of the inductive modules: o = arctan.ypos; xpos/: (19) Alternatively, this can also be determined using the y-offsets at the front and rear of the vehicle, respectively. 6 Simulative Evaluation of Inductive Omnidirectional Navigation The presented algorithm is implemented inMatlab Simulink. For this purpose, two different vehicle movements of the introduced AGV with a vehicle specific curvature limit of �g = 0.4845m−1 on the inductive floor are simulated and evaluated. The two paths of the AGV are shown in red and green in Fig. 9. The vehicle travels over two floor tiles. For the first evaluation (red path), a constant slip angle ˇ = −3ı is used as input for the OCP transformation and a constant nominal velocity vn = 0.0041. vn in this study case results in a velocity of 0.01986ms−1. This velocity is chosen that low, because it shows the inaccuracies dur- ing the travel. For a higher velocity the algorithm is also working. The nominal curvature is set to zero, �n = 0, to evaluate the translatory movement of the vehicle. The result of the reconstruction of the position of the vehicle with the inductive floor is shown in Figs. 10, 11 and 12. Fig. 10 shows the evaluation of the x-position of the AGV. After the vehicle center leaves one inductive floor module the x-position is set to zero. The overtraveled modules of the vehicle, calculated in solid blue agree with the simulation Fig. 9 Simulative scenario on the inductive floor. The virtual target path for the AGV is marked in red. The vehicle drives with �n = 0, vn = 0.0041 and ˇ = −3ı. The second virtual path with �n = −0.2, vn = 0.0041 and ˇ = −2ı is marked in dotted green Fig. 10 Simulative evaluation of the x-position of the vehicle driving over two inductive modules. The calculation of the vehicle model is shown in dashed orange. The estimation of the position of the Kalman filter is shown in solid blue K Forsch Ingenieurwes (2023) 87:617–626 625 Fig. 11 Simulative evaluation of the y-position of the vehicle driv- ing over two inductive modules. The calculation of the AGV model is shown in dashed orange, the estimation of the position after the Kalman filter is shown in solid blue Fig. 12 Simulative evaluation of the rotation of the vehicle on the in- ductive floor. The calculation of the vehicle model is shown in dashed orange, the estimation of the position after the Kalman filter is shown in solid blue Fig. 13 Simulative evaluation of the AGV travel on the inductive floor. In this scenario the AGV drives a curvature with �n = −0.2; vn = 0.0041 and a slip angle ˇ = −2ı results shown in dashed orange. Since the simulation of the inductive coils is performed with steps of 0.0255m, small jumps occur in the estimated values. These can be reduced with more accurate electromagnetic simulation or meas- urement values. With the constant slip angle the constant y-deviation of the vehicle from the centre of the vehicle (shown in Fig. 11) is correctly estimated by the Kalman filter (solid blue) compared to the created vehicle model (dashed orange). Because of the simulation steps in y- direction of 0.00425m there are inaccuracies in this order. The rotation of ˇ = −3ı of the AGV with respect to the floor coordinate system is also evaluated. This show case is presented in Fig. 12. The rotation is inaccurate on the first inductive module because of the initial conditions and the inaccuracies of the y-position. On the second module the rotation is evaluated successfully. In the next simulation experiment (green path in Fig. 9) the AGV is simulated to drive a curvature with �n = −0.2, a side-slip angle of ˇ = −2ı and nominal velocity vn = 0.0041. The result of the simulation is shown in Fig. 13. xpos-position and rotation are almost similar to the first ex- periment and therefore not shown. The evaluation of the travel with a curvature, shown in Fig. 13 is also success- fully evaluated by the Kalman filter (solid blue) compared to the vehicle model. The inaccuracy is as in the transla- tional drive 0.00425m. In practice the simulation values can be replaced with measurements on one inductive module to increase the prac- tical functionality of the algorithm. 7 Conclusion This paper answered the research question: Can the con- tactless induced energy supply from a novel inductive floor be used to navigate omnidirectional AGVs? To ensure the AGV follows the requested path trough production a control algorithm is necessary. For controlling the omnidirectional AGV consistently without singularities, a transformation in Omni-Curve-Parameters (OCP) is pro- posed. To use the floor for tracking, sensors placed in the AGV measure induced voltages. Depending on these in- duced voltages of the inductive floor a finite element method is used to simulate the movement. Depending on the simu- lated voltages a novel algorithm is presented which evalu- ates the voltages to positions on the floor. These positions are combined with an odometry model and a Kalman filter is used to estimate the tracking on the floor. As case study a four wheeled steering- and velocity controlled AGV is introduced. For the evaluation a novel motion model de- pending on the input OCP is proposed. Simulation results of this vehicle with the inductive floor demonstrate, that the K 626 Forsch Ingenieurwes (2023) 87:617–626 tracking of a vehicle on the inductive floor are accurate to 0.425cm order and the results are reproducible. In practice, the accuracy of the method needs to be inves- tigated in the future. It will probably not work as accurately as in the simulation, due to manufacturing inaccuracies, for example. In addition, an investigation of the possible speed in practice must be carried out. If the plates are passed over too quickly, clear position detection can no longer be guar- anteed. Since the presented simulated model is based on physical principles and parts of it could already be veri- fied, the research question is answered theoretically. In the near future, the implementation into reality still has to take place. The presented algorithm can also be used for wireless charging as extension e.g. introduced in [17]. Based on the introduced estimation method of the po- sition with the inductive floor, path planning can be per- formed. With the position estimated in the Kalman filter this can be considered equivalent to the existing navigation methods. The presented floor is conceptually created and individual plates are already being tested as demonstrators. The introduced vehicle model and control algorithm de- pending on the OCP have already been tested on the vehicle prototype and the functionality could already be verified. In further research work, a methodology will be devel- oped, to implement a control algorithm for an omnidirec- tional vehicle in OCP for different guidance control types and a validation is expected in the near future. Funding Open Access funding enabled and organized by Projekt DEAL. Open Access This article is licensed under a Creative Commons At- tribution 4.0 International License, which permits use, sharing, adapta- tion, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, pro- vide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4. 0/. References 1. Brown R (2006) Topology and groupoids: a geometric account of general topology, homotopy types and the fundamental groupoid. BookSurge, Deganwy. ISBN 978-1-4196-2722-4 (REV Updtd and Exp Ver ed) 2. Colomb A, Brenner C (2020) Konzept zur intuitiven Steuerung om- nidirektionaler Flurfoerderzeuge mit beliebiger Radkonfiguration. Logist J Proc. https://doi.org/10.2195/LJ_PROC_COLOMB_DE_ 202012_01 3. Connette C (2013) Kinematische Modellierung und Regelung om- nidirektionaler, nicht-holonomer Fahrwerke. Fraunhofer-Verlag, Stuttgart. ISBN 978-3-8396-0564-6. 4. Demuth R (2021) Induktiver Lenksensor. https://www.goetting.de/ komponenten/19370-19380. Accessed 22 Nov 2021 5. DIN 13312 (2005) Navigation – Concepts, abbreviations, letter symbols, graphical symbols. Berlin, Germany: Deutsches Institut für Normung e.V 6. Foith-Förster P, Bauernhansl T (2015) Changeable and reconfig- urable assembly systems – A structure planning approach in auto- motive manufacturing. 15. Internationales Stuttgarter Symposium, Bargende, pp 2198–7440 7. Hofmann M (2018) Intralogistikkomponenten für die Automo- bilproduktion ohne Band und Takt – erste Prototypen. Logist J. https://doi.org/10.2195/LJ_PROC_HOFMANN_DE_201811_01 (Proceedings) 8. Li L, Schulze L (2018) Comparison and Evaluation of SLAM Allgorithms for AGV Navigation. Production engineering and ma- nagement. Proceedings 8th international conference, Lemgo. ISBN 978-3-946856-03-0. 9. Manley E, Al Nahas H, Deogun J (2006) Localization and tracking in sensor systems https://doi.org/10.1109/SUTC.2006.83 10. Marín L, Vallés M, Soriano N et al (2022) Multi sensor fusion framework for indoor-outdoor localization of limited resource mo- bile robots https://doi.org/10.3390/s131014133 11. Peter N, Reinhart G, Abele E (2008) Wandlungsfähige Produktion- ssysteme: Heute die Industrie von morgen gestalten. Verlag PZH Produktionstechnisches Zentrum, Garbsen 12. Schumacher P (2021) Technische Logistik Flurförderzeuge 2021/2022. Aufbrechen der Perlenkette 13. Shi D, Mi H, Collins EG et al (2020) An indoor low-cost and high- accuracy localization approach for agvs. IEEE Access 8:50,085– 50,090. https://doi.org/10.1109/ACCESS.2020.2980364 14. Stillig J, Parspour N (2020a) Der intelligente Boden: Eine In- frastrukturplattform für die wandelbare Produktion. 1 Stuttgarter Tagung zur Zukunft der Automobilproduktion, Stuttgart, pp 38–39 15. Stillig J, Parspour N (2020b) Novel concept for wireless power transfer modules. 2nd International Conference on Broadband Communications, Yogyakarta. https://doi.org/10.1109/BCWSP 50066.2020.9249465 16. Stillig J, Enssle A, Parspour N (2021) Modular wireless power transfer system for thesupply of mobile industrial production equipment. IEEE PELS Workshop on Emerging Technologies: Wireless Power Transfer (WoW), San Diego. https://doi.org/10. 1109/WoW51332.2021.9462864 17. Tan L, Li C, Li J et al (2021) Mesh-based accurate positioning strategy of ev wireless charging coil with detection coils. IEEE Trans Ind Informatics 17(5):3176–3185. https://doi.org/10.1109/ TII.2020.3007870 18. Ullrich G (2014) Fahrerlose Transportsysteme https://doi.org/10. 1007/978-3-8348-2592-6 19. VDI 2510 (2005) Automated Guided Vehicle Systems (AGVS). Düsseldorf, Germany: Verein Deutscher Ingenieure e.V. 20. VDI 4451 (2005) Compatibility of automated guidedvehicle sys- tems (AGVS) Sensor systems for navigation and control. Düssel- dorf, Germany: Verein Deutscher Ingenieure e.V. 21. Wiendahl HP, Reichardt J, Nyhuis P (2009) Handbuch Fabrikpla- nung – Konzept, Gestaltung und Umsetzung wandlungsfähiger Pro- duktionsstätten. Carl Hanser, München, Wien. ISBN 978-3-446- 43702-9. 22. Xiong J, Liu Y, Ye X et al (2016) A hybrid lidar-based indoor nav- igation system enhanced by ceiling visual codes for mobile robots. International Conference on Robotics and Biomimetics, Qingdao. https://doi.org/10.1109/ROBIO.2016.7866575 K http://creativecommons.org/licenses/by/4.0/ http://creativecommons.org/licenses/by/4.0/ https://doi.org/10.2195/LJ_PROC_COLOMB_DE_202012_01 https://doi.org/10.2195/LJ_PROC_COLOMB_DE_202012_01 https://www.goetting.de/komponenten/19370-19380 https://www.goetting.de/komponenten/19370-19380 https://doi.org/10.2195/LJ_PROC_HOFMANN_DE_201811_01 https://doi.org/10.1109/SUTC.2006.83 https://doi.org/10.3390/s131014133 https://doi.org/10.1109/ACCESS.2020.2980364 https://doi.org/10.1109/BCWSP50066.2020.9249465 https://doi.org/10.1109/BCWSP50066.2020.9249465 https://doi.org/10.1109/WoW51332.2021.9462864 https://doi.org/10.1109/WoW51332.2021.9462864 https://doi.org/10.1109/TII.2020.3007870 https://doi.org/10.1109/TII.2020.3007870 https://doi.org/10.1007/978-3-8348-2592-6 https://doi.org/10.1007/978-3-8348-2592-6 https://doi.org/10.1109/ROBIO.2016.7866575 Navigation method utilizing floor-integrated inductive power supply modules for omnidirectional AGVs Abstract Zusammenfassung Introduction Structure of Localization and Motion Control of the AGV with the Inductive Floor Omni-Curve-Parameters as Universal Transformation for Omnidirectional AGVs Reference Vehicle Movement Estimation of the Vehicle Position with the Inductive Floor Simulative Evaluation of Inductive Omnidirectional Navigation Conclusion References