An Optimization-Based Motion Planner for Car-like Logistics Robots on Narrow Roads

Thanks to their strong maneuverability and high load capacity, car-like robots with non-holonomic constraints are often used in logistics to improve efficiency. However, it is difficult to plan a safe and smooth optimal path in real time on the restricted narrow roads of the logistics park. To solve this problem, an optimization-based motion planning method inspired by the Timed-Elastic-Band algorithm is proposed, called Narrow-Roads-Timed-Elastic-Band (NRTEB). Three optimization modules are added to the inner and outer workflow of the Timed-Elastic-Band framework. The simulation results show that the proposed method achieves safe reversing planning on narrow roads while the jerk of the trajectory is reduced by 72.11% compared to the original method. Real-world experiments reveal that the proposed method safely and smoothly avoids dynamic obstacles in real time when navigating forward and backward. The motion planner provides a safer and smoother trajectory for car-like robots on narrow roads in real time, which greatly enhances the safety, robustness and reliability of the Timed-Elastic-Band planner in logistics parks.


Introduction
With the continuous development of service robot technology, the logistics industry is gradually transforming and upgrading. The traditional manual handling and logistics management model is being eliminated. The automation mode with logistics robots is widely used at the core of factories and production enterprises [1,2]. The intelligent operation of logistics handling has become the trend of industry development. Among them, motion planning is the key technology of autonomous navigation, which determines whether the robot can safely and efficiently reach the target. Motion planning is the search for a collision-free optimal trajectory from the current position to a local target in the dynamic environment around the robot [3]. Assuming that the localization [4], mapping [5], perception [6], and routing [7] modules perform well while the robot navigates on the wide-structured roads, this type of general motion planning approach has been well studied [8][9][10][11]. However, in logistics parks, robots travel on diverse and complex roads. Due to the size limitation of the park, the environment is mostly unstructured roads and presents a network or ring shape [12]. Unlike structured roads, unstructured roads generally do not have lane lines. Therefore, the sensors cannot recognize lane information. The drivable area cannot be determined based on the lane lines [13]. A safe path can only be reasonably planned based on the location of obstacles. To make matters worse, the park inevitably contains some narrow curves and straights. Due to the narrow exploration space, higher demands are placed on the robot motion planning capability. First, the robot needs to track the global path smoothly. To ensure the integrity of the goods, the robot's speed during transport needs to be as smooth as possible. Second, the transportation process has the potential to involve reversing conditions. Non-holonomic constrained car-like robots also present a challenge in reversing on narrow roads. Finally, the robot often needs to interact with dynamic obstacles such as pedestrians and vehicles in the park during navigation.
Ensuring collision-free transportation is the basis for realizing motion planning in the park. Various types of algorithms have been made available with poor applicability in this type of scenario [14]. Although the Timed-Elastic-Band algorithm is widely used in various mobile robot navigation systems for its advantages such as rapid solution and high applicability, it is still not sufficient to fully cope with these problems. Based on the above motivation, this paper proposes the Narrow-Roads-Timed-Elastic-Band (NRTEB) method for narrow space planning. This aims to improve the safety, stability and robustness of car-like logistics robots for motion planning on narrow roads. The main contributions of this paper are: (1) We reconstruct the inner and outer workflow of the Timed-Elastic-Band algorithm to improve its motion planning performance on narrow roads.
(2) The jerk is added to the objective function in the form of a soft constraint, which improves the smoothness of the speed profile without significantly increasing the computation time.
(3) A strategy is designed to enhance safety when navigating in reverse, thereby reducing the risk of collision with narrow-road edges.
(4) A sampling-based obstacle avoidance strategy is proposed for the presence of obstacles on local targets to improve the safety and stability of narrow-road obstacle avoidance.
A comparison of the method in this paper with the framework of the TEB is shown in Figure 1.

Researches on Point-to-Point Motion Planning Methods
For autonomous mobile robots that need to accomplish specific tasks, their operation process can be divided into segments of point-to-point motion planning [15]. For example, logistics robots continuously move from the current target to the next target in a set of freight handling processes. For point-to-point motion planning methods, they can be broadly classified into graph search-based ones, sampling-based ones, curve-interpolationbased ones, optimization-based ones and deep reinforcement learning-based ones [16]. One of the earliest point-to-point path planning methods is the graph search method. A well-known method is the shortest path algorithm proposed by Dijkstra [17]. In order to improve the search efficiency, the A* algorithm was proposed later [18]. This algorithm significantly reduces the search time by estimating the cost from the current location to the target through a heuristic function. The graph search algorithm is the most mature global path planning method, of which the A* algorithm has been widely used in the fields of autonomous driving and robotics. However, it only considers static obstacles and the path does not contain time information, so it is not applicable to dynamic local trajectory planning. Sampling-based schemes were initially studied to deal with pathfinding problems in high-dimensional spaces. One of the most representative is the rapid-exploration random tree (RRT) [19] method. A series of sampling-based studies were then enriched, such as RRT* [20] and informed-RRT* [21]. These both aim to improve the sampling efficiency and thus speed up the pathfinding process, but the paths are usually not smooth enough and optimality is not guaranteed. Curve interpolation methods are often used to generate smooth, comfortable trajectories. Common interpolated curves are polynomial [22], Bezier [23], and spline curves [24]. Since they are interpolation operations based on waypoints, the quality of local trajectory planning results depends on the quality of global path planning. In recent years, numerical optimization-based methods have attracted a great deal of attention [25]. This class of methods models the search problem as an optimization problem. Various constraints such as velocity, acceleration, minimum steering radius, etc. are incorporated into a unified model for solving the problem. Because of its ability to explicitly handle dynamic obstacles and various types of constraints, it is widely used in various autonomous driving and robotic systems. For example, Baidu's motion planning module for unmanned vehicles extensively uses this type of approach to generate a short trajectory that is safe, comfortable, and collision-free [26]. In the field of robotics, the most popular one belonging to this category is the TEB local trajectory planner [27]. However, the TEB method ignores the comfort of speed change [28], which is detrimental to logistics handling robots. Although deep reinforcement learning-based motion planning methods are becoming a hot research topic, most of the work is still in the simulation stage. Such methods have not been adopted in the industry for the direct control of robots because safety and stability cannot be guaranteed [29].

Studies on Narrow Road Motion Planning Methods
A small amount of literature exists in the area of intelligent driving in narrow-road scenarios. The "narrow road assistant" proposed in [30] provided drivers with collision warning and steering support for driving on narrow roads. The driver-adapted ADAS proposed in [31] matched the driver's driving style to provide appropriate steering assistance. Sato, T. et al. [32] presented two driving assistance system models based on drivers' cognitive processes that are useful for achieving smooth passage on narrow roads and improving driver performance. In addition to advanced driver assistance systems applied on structured roads, some researchers have explored methods for the motion planning of mobile robots in narrow, unstructured road environments. Tian et al. [33] proposed a local path planning method for turning around in a narrow environment which combines Bezier curve and the minimum turning radius arc. The safety of planning in narrow roads is improved, but the real-time performance is poor. Li et al. [34] improved the planning performance of RRT in narrow roads or cluttered environments by fusing bidirectional RRT and rapidly exploring random vines. However, these paths are not smooth. Therefore, they cannot be provided to car-like robots for execution.

Studies on Timed-Elastic-Band Methods
In the scope of the TEB method, a variety of improved algorithms have successively been proposed. TEB adopts a strategy of predictive control to guide the robot to accomplish online trajectory planning tasks in a dynamic environment [35]. The first control action is finally adopted as the robot's command [36]. The TEB planning algorithm proposed by C. Rösmann [37] is an evolution of the elastic band (EB) [38], which explicitly adds the time intervals to its optimization. The original TEB planner might cause a trajectory across obstacles which are owed to the optimization being trapped in local minima. Subsequently, C. Rösmann proposed an improved method for generating alternative suboptimal trajectory clusters based on distinctive topologies [39]. After the TEB algorithm was adopted by the ROS navigation stack [40], it has gradually become a hot topic of study among scholars. Keller et al. [41] studied the optimal obstacle avoidance trajectory generation method with the TEB framework. Their approach improves the obstacle avoidance capability of the TEB, but does not talk about the real-time planning problem. Ulbrich et al. [42] proposed a variant of TEB that uses loose ends in planning to produce more human-like paths. However, the authors admit that the method still has the probability of falling into local minima and failing to provide safe trajectories. S. Smith et al. [43] presented an egocentric map representation for TEB in an uncertain environment to solve the problem of mismatch between optimization graph and grid-based map. This reduces the planning time but is only validated in a simulated environment, so it is not sufficiently convincing to provide guidance for a real robot. A car-like sweeping robot system based on TEB local trajectory planning is proposed in [15]. The method improves the sweeping coverage, but the speed change during navigation is not smooth. The previous works have more or less improved the trajectory smoothness, solution stability, and obstacle avoidance safety of TEB. However, none of them can perfectly solve all the problems of TEB in narrow-road planning, especially the reverse planning and reverse obstacle avoidance scenarios that logistics robots often encounter. In addition, most of the studies only show the optimization effect of a certain aspect. There is a lack of real-time results presentation, which is very important for motion planning. Therefore, this paper intends to comprehensively address the problems arising from TEB in narrow-road planning. The effectiveness and superiority of the method is verified by simulation experiments and real-world experiments. This will provide good guidance for car-like robots driving on the narrow roads of logistics parks.
There are several reasons why the TEB algorithm can be widely used in robot motion planning: (1) Strong solving ability: To reduce costs, logistics robots usually have limited computational resources. TEB is able to solve optimization problems very quickly. In the objective function of TEB optimization, many error terms are partial [37]. As such, the Hessian matrix is sparse so that the solution is rapid. (2) Wide applicability: TEB is not only applicable to differential wheels but also the chassis of the Ackermann structure [44]. This makes a lot of sense for ROS [45] because, in the early stage of development, ROS studied mostly round universal wheel robots and rarely considered car-like robots [46]. In addition, TEB considered different vehicle chassis contour models, different obstacle types, and whether the obstacles are dynamic, so it is more adaptable to obstacle avoidance in dynamic environments. (3) Insensitive to the smoothness of global reference path: TEB selects the farthest point of the global path within the local sensing range as the local planning terminal. The shape of the trajectory is adjusted by various constraints, so it does not require the global path to conform to the kinematic constraints of the robot. Unlike Baidu Apollo's local planner [26], it heavily relies on global path smoothing, so it reduces the overhead time of path smoothing. However, there are still the following problems that limit its application to car-like robots on narrow roads. (1) Speed planning results are not smooth: Although the non-holonomic constraint ensures smooth paths, both angular and linear speed profiles are not smooth. This is due to the fact that the TEB does not constrain the jerk of the entire trajectory. The unsmoothed speed planning results mean that there is a possibility of unstable control, which is contrary to the need for the robot motion and steering to be as smooth as possible in narrow-road scenarios. (2) Narrow road reversing planning is flawed: At the beginning of a new planning cycle, the TEB method selects a waypoint from the global path as a local target based on the forward-looking distance. When the global path is behind the robot, the heading of the local target is opposite the robot heading, because the global planning starts from the robot and points to the target. Based on this problem, TEB generates a parking trajectory, which is unacceptable for reversing behavior in confined environments. The robot is highly susceptible to collisions with road edges or cannot turn around due to the minimum turning radius. (3) Planning fails when obstacles fall on the local target: A local target is a waypoint obtained from the global path. Its position is fixed at the beginning of the planning cycle. If at this point the sensors detect a dynamic obstacle moving to the position of the local target, it will cause the planning to fail because the target is unreachable. This can be solved by updating the global path, however, global planning is unreliable as the module with the slowest update frequency and computation takes a long time. Additionally, the introduction of dynamic obstacles makes the search space non-convex and may lead the optimization into local minima.

Materials and Methods
For the above-mentioned problems arising in narrow road environments, we reconstruct the inner and outer workflow of the TEB algorithm. The improvements in the internal optimization loop are presented first, corresponding to the orange module in Figure 1. Then, the two working modules in the external working loop are introduced, corresponding to the blue modules in Figure 1.

Motion Planning Modeling
Car-like robots with an Ackermann structure traveling on the configuration space, which is shown in Figure 2. This paper adopts the kinematic model [47] for the Ackermann chassis.
where (x, y, θ) is the configuration space state, (x, y) is a two-dimensional Cartesian coordinate, θ denotes the heading angle of the robot, v is the longitudinal speed, δ f is the front steering angle, and L is the wheelbase of the vehicle. The time intervals between poses are added to the optimization variables. Therefore, TEB not only modifies the geometry of the trajectory but also optimizes the velocity and acceleration, thus solving the high-dimensional search problem of motion planning. The sequence of poses in configuration space is expressed as: This denotes a robot pose in the configuration space. The time interval between each two poses is denoted by ∆T, and the set of time interval vertices is expressed as: where ∆T i is the time interval between q i and q i+1 . Thus, the TEB trajectory is represented as a combination of a sequence of poses and a sequence of time intervals, which are the optimization variables of the subsequent solution: To facilitate the transformation of multi-objective optimization problems into graph optimization problems, the TEB algorithm transforms the hard constraint into a soft one. Therefore, this paper adopts the piecewise continuous, differentiable uniform cost function described in [27] to penalize the violation of constraints: where x r presents the bound value. n is the polynomial order, R is the scaling, and expresses the safety margin. In this paper, we choose x r = 0.5, = 0.05, R = 0.25, n = 6. The multi-objective function of TEB can be expressed as: where f k (S) represents the deviation of each constraint from the desired value, and γ k denotes the corresponding weights. f v = e(v i , v max , , R, n) represents the velocity constraint of the trajectory, where v i is obtained by differentiating the adjacent pose q i ,q i+1 from their time intervals ∆T i . f a = e(a i , a max , , R, n) represents the acceleration limitation of the trajectory, where a i is obtained from the three adjacent poses q i , q i+1 , q i+2 with their two time intervals ∆T i , ∆T i+1 . f path = e(q i , p i , , R, n) is the target-following constraint. f obs = e(−d min,j , −r min , , R, n) indicates that the obstacle constraint for those poses that are less than the minimum safe distance r min . f time = (∑ n i=1 ∆T i ) 2 is the shortest time constraint. f nh represents a nonholonomic constraint between two adjacent poses which guarantees the generation of a trajectory that conforms to the vehicle kinematics. We refer to [39] for more details on constraint modeling.
Finally, Equation (6) is transformed into a graph optimization problem, as shown in Figure 3b. It is solved with the optimization framework of G2O [48] to obtain the optimal motion: where S * is the optimal TEB trajectory obtained from the optimal solution.

Jerk Limitation
The linear jerk is the derivative of linear acceleration while the angular jerk is the derivative of angular acceleration. Constraining the jerk smooths out the speed profile of the trajectory, thus increasing the stability of the motion on narrow roads. A linear jerk is derived by the difference method from the four adjacent poses of S and the three time intervals between them, as shown in Figure 3a. The linear velocity is first represented by the sequence of poses q: The linear acceleration is obtained from the three adjacent poses with their two time intervals, and thus the linear acceleration is represented by the linear velocity v i , v i+1 and the time interval ∆T, ∆T i+1 : Finally, the linear jerk is represented by the acceleration a vi , a vi+1 with their three adjacent time intervals ∆T, ∆T i+1 , ∆T i+2 as follows: In order to obtain the exact difference result, it is necessary to ensure that α + β + γ = 1. ∆T i+1 is correlated with both a vi and a vi+1 , as seen from Figure 3. ∆T i and ∆T i+2 are only correlated with a vi , a vi+1 , respectively. Thus, ∆T i+1 is given greater weight. This paper let Similarly, the angular jerk can be obtained as follows: where a ω indicates the angular acceleration of the robot.
The constraint function of the linear jerk and angular jerk is expressed as: where J vmax is the threshold of the linear jerk and J ω max is the threshold of the angular jerk. Therefore, the multi-objective function of NRTEB can be expressed as: where γ J v is the optimization weight of the linear jerk constraint and γ J ω is the optimization weight of the angular jerk constraint.
With the above definition, the multi-objective problem is transformed into a graph optimization problem, where four poses and three time intervals constitute the minimum hyper-graph considering jerk, as shown in Figure 3a. A simplified example structure of the NRTEB hyper-graph is shown in Figure 3b.

Reverse Planning Enhancement
This paper has already explained the shortcomings of the original TEB in reverse planning in Section 2.3. Due to the problem of local target reference heading, TEB always tries to generate a parking trajectory when the global path is behind the vehicle, as shown in Figure 4a. This reversing strategy is extremely unsafe in narrow roads. To enhance the safety and robustness of reversing planning, this paper intends to solve the problem from local trajectory planning cycle instead of the global planning module. When the global path is behind the vehicle and the robot is far from the global target, NRTEB is able to generate a direct reversing trajectory, as shown in Figure 4b. When the robot approaches the global target, it takes this goal as the end point of optimization. As such, the robot's terminating heading problem is also solved. A vector product is used to represent the position of a point in relation to the robot. If the vector product of the body direction vector and vector formed by the selected point and the robot center of mass is negative, the picked point is considered to be behind the robot, i.e., Point b in Figure 5. The reference heading of the local target is corrected by judging this condition, so as to keep the heading of the planning start and the planning end consistent. The workflow for handling reversing is shown in Algorithm A1 in Appendix A, which indicates the corresponding processBackward module in Figure 1.

Obstacle Avoidance Strategy
The essence of TEB planning is to periodically solve an optimization problem. The search space becomes non-convex when dynamic obstacles appear on the local targets, which leads to an optimization that may fall into local minima. In order to make the local target reachable and convexize the search space, an obstacle avoidance strategy that samples first and searches later to determine the local target location is proposed as shown in Figure 6. x gtl = x lt − l sin θ y gtl = y lt + l cos θ θ gtl = θ + ∆θ (14) where ∆θ and l are the parameters determined by the working environment. From the symmetry relationship, the pose of the right-guided turning point is: Secondly, on the basis of guided turning points, the set of obstacle avoidance points is sampled according to the preset parameters and geometric relationships, as shown in Figure 6b. Then, the set of sampling points is traversed until a point that meets the collision constraint is found as the new local target for the current planning cycle.
To obtain the sampling point, the assist point is first calculated according to Equation (16): where (x gt , y gt ) is the position of guided turning point and ∆l denotes the determined lateral variation.
With the help of an assist point, the sampling point is then calculated by Equation (17): where (x sp , y sp , θ sp ) is the pose of sampling point, ∆s denotes the determined longitudinal variation, and ∆φ denotes the determined angular variation. The decision process for the local target is shown in Algorithm A2 in Appendix A, which indicates the corresponding processObstacleOnLocalTarget function in Figure 1. Finally, the algorithm introduces two hyperparameters, l and ∆θ, which represent the amount of lateral distance and angle change, respectively. There is a relationship between the choice of these two parameters and the minimum turning radius because the car-like robot is subject to non-holonomic constraints. A simple proof is as follows: Assume that the distance between the robot and the local target is s, according to the Pythagorean theorem: e = l 2 + s 2 The robot kinematics are geometrically related as shown in Figure 7, so the relationship between ∆θ and the vehicle turning radius R is given by: Bringing Equation (18) into Equation (19), when the value of l is determined, the range of values of ∆θ is: where R min represents the minimum turning radius of the robot. Therefore, the value of ∆θ needs to be determined according to the constraint range of Equation (20), while l needs to be selected as appropriate according to the road width. s depends on the lookahead distance which can be set as a hyperparameter in NRTEB.

Results and Discussion
In this section, firstly, the trajectory of our method is compared with the advanced motion planning method in a constrained environment. Then, a simulation environment containing narrow corridors is built in Gazebo (https://github.com/osrf/gazebo, accessed on 15 January 2010) to compare the motion performance, reversing performance, and realtime performance of each method. Finally, we evaluate the reversing and obstacle avoidance performance of the proposed method in the narrow-road scenario on a real car-like robot platform. The experimental results of the article are reproducible with the same robot structure and the same experimental parameters. However, note that it needs to be assumed that the localization, perception, and global path planning modules work properly during the experiment.

Simulation and Analysis
The method of this paper is compared with the original TEB and MPC applied to the car-like robot model. The simulation experiments were conducted on a laptop with an Intel Core i5-6300HQ CPU 2.30 GHz × 4. The influence of algorithm parameters on the trajectory is also discussed.

Trajectory Motion Performance Comparison
The motion planner usually updates the trajectory according to a certain frequency, so the trajectory has been dynamically varying as the robot moves and the environment changes. In order to accurately compare the motion performance of trajectories generated by different motion planners in the same environment, the start and end point are fixed to simulate the robot position and local target.  Table 1. The trajectory generated by NRTEB is shown in Figure 8. Table 1. Kino-dynamic parameters and optimization parameters.

Constraint Parameters Values Application Scope
Maximum The motion performance of the trajectories generated by the three planners is shown in Figure 9. All three methods generate the correct obstacle avoidance path where MPC has the closest distance to the obstacle and the highest collision risk, but can reach the end point faster. While the paths generated by TEB and NRTEB are relatively safer. The speed profile of TEB and MPC is steeper and the yaw change is more rapid, which poses a more stringent challenge to the tracking stability. In contrast, the NRTEB has the smoothest speed profile, which significantly reduces jitter and oscillation during travel. This is important for logistics robots to avoid shaking and thus ensure the integrity of the products. This paper also explores the effect of planning parameters on the performance of the NRTEB planner for trajectories. The different jerk bounds or optimization weights are set to compare the planning performance in the same environment, respectively. The results are shown in Figure 10. The tighter the jerk constraint, the greater the fluctuation of the planning path. It can be seen from the heading angle, linear velocity, and angular velocity curves that the tighter the jerk limit, the smoother the curve, and the longer the motion time at the same time, where the motion curve is smoothest for a jerk threshold of 0.1 and steepest for a jerk threshold of 4. An increase in the optimization weight of jerk further improves the motion smoothness and motion time. In this paper, we suggest that the constraint range of jerk is between 0.1 and 5, and the optimization weight of jerk is between 1 and 10. If there is a higher requirement for the motion time, the optimization weight of the trajectory motion duration can be further adjusted upward, thus making a balanced trade-off between motion smoothness and motion time. The global path is provided and recorded using the A* algorithm, and the same global path is kept published during the experiments. The simulation environment containing narrow corridors is shown in Figure 11. The key parameters for the simulation experiments are shown in Table 2. All parameters do not exceed the kino-dynamic constraints of the car-like robot model. Figure 11. The simulation environment. The green curve represents the global path recorded in advance. The red curve with arrows represents the motion trajectory and the blue arrow indicates the local target.

Constraint Parameters Values
Maximum linear velocity (m/s) 0. Firstly, the pre-recorded global path is published, and different motion planning algorithms are used to generate local motion trajectories. Then, the motion state of the robot is recorded in real-time from the starting point to the global target with a frequency of 2 Hz. The traces of the robot's odometry are shown in Figure 12. The results of the actual robot path length and average lateral tracking error are shown in Table 3. Among them, MPC generates the shortest tracking path but has the largest tracking error, which makes it significantly riskier to collide with walls at the corners of narrow roads. TEB and NRTEB have a smaller tracking error, which allows the robot to stay in the middle of the narrow road. NRTEB requires a larger turning radius because it limits the jerk, but does not significantly affect the tracking performance.  The robot status data during navigation is shown in Table 4. It should be noted that, to reflect the degree of change in motion, the average values in Table 4 need to be averaged after taking the absolute values of the jerk data first. The motion of MPC took 117.495 s, the motion of TEB took 124.516 s, and the motion of NRTEB took 137.017 s. The results again verified that NRTEB can improve the speed smoothness by constraining the jerk of the trajectory. However, the side effect is to increase the steering radius and motion duration. It can be solved by adjusting the constraint range of jerk, the magnitude of jerk optimization weight and the weight of the trajectory motion duration, as described in the previous section. The average planning time consumption of MPC is 154 ms, which cannot guarantee real-time planning when the CPU computing power is limited, while the planning efficiency of TEB and NRTEB is significantly improved. It is further demonstrated that adding the jerk to the objective function for optimization does not significantly increase the planning time consumption since the jerk term is partial with respect to the whole trajectory, thus ensuring real-time planning.

Comparison of Reversing Trajectory
When the global path is behind the robot, TEB generates the parking trajectory for turning around, but this does not work on narrow roads, as shown in Figure 13a,b. Due to the minimum steering radius constraint of the car-like robot, the trajectory is very likely to cause the robot to collide with the boundary or unable to turn around, thus leading to planning failure. MPC is even unable to generate a trajectory. NRTEB generates a straight reversing trajectory by modifying the reference heading of the local target, as shown in Figure 13c, which solves the problem of robot oscillation due to a minimum turning radius constraint when reversing on narrow roads. The initial state of the robot in this part of the experiment is the terminal state of the previous part, i.e., the robot pose is [1, 9, π] T and the global target is set to [1.5, 7, − π 4 ] T . The global reference path generated by the A* algorithm located behind the robot is shown in Figure 13d. The diagram shows that NRTEB guides the robot out of the narrow dead zone and eventually reaches a new global target, which is not possible with the other two methods.

Realistic Obstacle Avoidance Scenarios
Given that the simulation environment ignores the data transmission time delay and does not reflect the obstacle information of the robot in real narrow-road environments, the obstacle avoidance performance of NRTEB and TEB is compared through our selfdeveloped car-like logistics robot platform, as shown in Figure 14. The components of the automated guided vehicle (AGV) are shown in Table 5. Vehicle chassis parameters are shown in Table 6.  The AGV is tested on a 2.75 m-wide road with traffic guardrails. The planning parameters of NRTEB are shown in Table 7. When the vehicle is driving in the middle of the road, its distance from the edge of the road is 1.375 m, so set l to 1.0. According to Equation (20), the maximum value of ∆θ is 67.59 • , so choose ∆θ as 60. Our system integrates with the global road network planning module in Autoware (https://www.autoware.org/, accessed on 1 August 2015) to provide a global reference path. Figure 15a shows the real environment captured by the camera, while Figure 15b shows the visualization of the AGV in Rviz. Table 7. NRTEB parameters.

Planner Parameters
Values Lookahead_dist (  To validate our planner, it is tested repeatedly on the loop road as shown in Figure 15. During navigation, pedestrians pass through the global reference path with irregular motion. A global target behind the robot is also set to test its reversal planning capabilities. The comparison experiments are conducted on the same narrow roads while keeping the same parameters. The results are shown in Figure 16. When an obstacle appears on the global reference path of the AGV, as shown in Figure 16a, the TEB planner suffers from planning confusion which is caused by the fact that the search space becomes non-convex due to the intervention of dynamic obstacles, so that the TEB planning falls into local minima, as shown in Figure 16e, while NRTEB updates the optimization endpoint by sampling search to make the non-convex problem convex so it can avoid falling into local minima and generate optimal trajectories, as shown in Figure 16i. When the AGV has enough lateral displacement after following the first few points, the trajectory following the global reference path is generated again after the local target leaves obstacles at the next moment, as shown in Figure 16b,f,j, and at the current time, the TEB planner also generates the correct trajectory but closer to the obstacle, causing tracking to be more difficult. When the global target is behind the AGV and in the same direction as the AGV, the TEB planner generates a parking trajectory, which leads the vehicle to becoming trapped at the edge of the road by not being able to turn around on narrow roads, as shown in Figure 16g, and the NRTEB planner generates a reasonable reversing trajectory, as shown in Figure 16k. In addition, our improvement in obstacle avoidance is also applicable to the reversing scenario. If there are obstacles on the local target in the reversing process, as shown in Figure 16d, the TEB planner also does not generate a reasonable trajectory, as shown in Figure 16h, and the NRTEB planner looks for a new local target without obstacles as the endpoint of the current planning cycle to generate a correct trajectory, as shown in Figure 16l.
The planning performance of the real-world navigation is shown in Table 8. The results show that NRTEB significantly improves the smoothness of speed changes, although it increases some navigation elapsed time. When collision detection is performed, the planning time consumption is approximately 100 milliseconds and no more than 50 milliseconds during normal driving. If the computational power of IPC is insufficient, the size of the local cost map can be reduced appropriately and the maximum frequency of motion planning can be decreased at the same time. It should be noted that the NRTEB solution can be performed under the current parameters, but when the planning frequency and the local cost map update frequency are too high-causing the computational capacity of the IPC to be exceeded-the motion planning solution no longer converges, which can be solved by replacing the CPU with a more computationally powerful one.

Conclusions
In this paper, a motion planning algorithm for logistics car-like robots was proposed. It greatly improves the planning capability of a TEB algorithm in a narrow-road environment. It makes up for the shortcomings of TEB in terms of an unsmooth speed profile, unsafe reverse planning, and unstable dynamic obstacle avoidance. Compared with TEB, NRTEB plans a smoother speed profile. The linear jerk of the trajectory is 55.48% smaller and the angular jerk is 72.11% smaller than that of TEB. This improves the comfort of robot motion and greatly reduces the probability of good being damaged due to sudden speed changes. When navigating without obstacles, NRTEB tracks the global path in real time and stably. NRTEB's tracking error is 24.46% smaller than MPC, and the planning time is reduced by 98.66%. In addition, NRTEB enables safe reverse navigation on narrow roads, which is not possible with the other two comparison methods. More significantly, real-world experiments demonstrate that NRTEB safely avoids dynamic obstacles in real time. It even achieves narrow-road obstacle avoidance whilst reversing.
NRTEB's planning efficiency can be further improved. Future work will deploy local map update modules on GPU for parallel processing. We recommend applying this method to car-like robots that need to work frequently in narrow-road environments, including but not limited to logistics, inspection, search, and rescue, etc.  Get the local_target(x, y, θ) from the last point of local_path → local_target 5: D local ← Euclidean_distance(P local_target , P robot ) 6: Calculate the direction vector between the robot and the local target → V lt 7: Calculate the direction vector of the robot → V robot 8: dot product ← V lt · V robot 9: if dot product < 0 then 10: θ local_target ← θ local_target + π 11: else 12: θ local_target ← θ local_target