Next Article in Journal
On Modulus Statistical Convergence in Partial Metric Spaces
Next Article in Special Issue
On a Neumann Problem with an Intrinsic Operator
Previous Article in Journal
AHD-SLE: Anomalous Hyperedge Detection on Hypergraph Symmetric Line Expansion
Previous Article in Special Issue
Exponential Stability of Dynamical Systems on Time Scales with Application to Multi-Agent Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Precise Obstacle Avoidance Movement for Three-Wheeled Mobile Robots: A Modified Curvature Tracking Method

1
Faculty of Science, Yibin University, Yibin 644000, China
2
School of Mathematics and Statistics, Guizhou University, Guiyang 550025, China
*
Author to whom correspondence should be addressed.
Axioms 2024, 13(6), 389; https://doi.org/10.3390/axioms13060389
Submission received: 24 April 2024 / Revised: 27 May 2024 / Accepted: 31 May 2024 / Published: 8 June 2024
(This article belongs to the Special Issue Recent Developments in Stability and Control of Dynamical Systems)

Abstract

:
This paper proposes a precise motion control strategy for a three-wheeled mobile robot with two driven rear wheels and one steered front wheel so that an obstacle avoidance motion task is able to be well implemented. Initially, the motion laws under nonholonomic constraints are expounded for the three-wheeled mobile robot in order to facilitate the derivation of its dynamic model. Subsequently, a prescribed target curve is converted into a speed target through the nonholonomic constraint of zero lateral speed. A modified dynamical tracking target that is aligned with the dynamic model is then developed based on the relative curvature of the prescribed curve. By applying this dynamical tracking target, path tracking precision is enhanced through appropriate selection of a yaw motion speed target, thus preventing speed errors from accumulating during relative curvature tracking. On this basis, integral sliding mode control and feedback linearization methods are adopted for designing robust controllers, enabling the accurate movement of the three-wheeled mobile robot along a given path. A theoretical analysis and simulation results corroborate the effectiveness of the proposed trajectory tracking control strategy in preventing off-target deviations, even with significant speed errors.

1. Introduction

Wheeled mobile robots, as a kind of frequently used mechanical model, are driven by actuators acting on wheels. Due to their simple structure, high maneuverability, and low energy consumption, wheeled mobile robots have garnered increasing attention in the fields of human transporters and humanoid robots [1,2]. They can be categorized based on the number of wheels, ranging from single-wheeled to multi-wheeled configurations. Among these, the three-wheeled mobile robot (TWMR for short) stands out due to its stability in the vertical direction—offering better stability compared to single or two-wheeled counterparts—and greater flexibility and smaller footprint compared to four-wheeled models [3]. These advantages make TWMRs widely applicable across various domains. For instance, the Canguro three-wheeled robot developed by Chiba University of Japan in 2008 exemplifies the practical use of TWMRs for following users and carrying items. Additionally, the advancement of autonomous driving technology for wheeled mobile robots holds promise for enhancing traffic safety, road utilization, and travel efficiency in the future. Accordingly, the most important feature for automated driving for TWMRs is the ability to avoid obstacles for safety. To achieve satisfactory obstacle avoidance performance, two steps need to be taken into account [4,5]. One is to provide a reasonable trajectory curve from the current position to a safe location, and the other is to develop an effective trajectory tracking strategy so that the TWMR can accurately track the desired path. In general, according to the types of wheels, TWMRs are mainly categorized into three groups, including two driven wheels with one steered wheel, two driven wheels with one castor wheel, and three omnidirectional wheels [6,7]. Among these, the case of two driven wheels with one steered wheel is characterized by a simpler structure, higher load capacity, and lower control torque requirements. Despite these advantages, the complex motion laws under nonholonomic constraints have limited the extent of research on this type, making it the focus of investigation in this paper.
Trajectory tracking control is a fundamental concern for TWMRs. Current research predominantly centers on kinematic models for designing speed controllers [8,9]. Nevertheless, in practical implementations, the required speeds are usually achieved by either forces or torques. Consequently, dynamic models should be considered to devise torque controllers [10,11], which are typically implemented through current or voltage controllers. For instance, a novel three-level dynamic tracking controller has been proposed that takes into account actuators and power stage subsystems; this is more accordant with practical situations [12,13]. Moreover, actual motion tasks often necessitate TWMRs to follow a specified path, whereas dynamic models only describe the relationship between speed and control torque. To address this obstruction, it is necessary to convert the required path curve to a speed target format to ensure compatibility with the dynamic model. In this sense, the original motion issue can be transformed to a conventional tracking control problem. Taking the two-wheeled inverted pendulum as an example, a prescribed trajectory curve is transformed into the tangential motion speed and yaw motion speed targets, thereby streamlining the initial motion control design problem [14].
Furthermore, in order to design an accurate and effective controller for achieving satisfactory trajectory tracking performance, it is imperative to thoroughly elucidate the motion laws of the TWMR. In recent years, many scholars have paid extensive attention to the motion laws of two-wheeled mobile robots [15,16]. However, to the best of our knowledge, limited work has explored the holonomic and nonholonomic constraints specific to TWMRs. In particular, in some existing studies, nonholonomic constraints are often overlooked. For example, the lateral speed of a three-wheeled vehicle is sometimes treated as a degree of freedom rather than a nonholonomic constraint [17,18]. Equipping TWMRs with three omnidirectional wheels enhances their motion flexibility and significantly reduces the number of nonholonomic constraints [19,20]. Moreover, in TWMRs featuring two driven wheels and one castor wheel, the nonholonomic constraints are akin to those of two-wheeled mobile robots, as the castor wheel serves a supporting role [21,22]. However, for the case of two differentially driven wheels and one steered front wheel, the nonholonomic constraints become more intricate, thus posing challenges in deriving dynamic models and motion laws [23]. Specifically, when the steering angle of the front wheel is considered as a control variable, the final dynamic model is different from the conventional form. As a consequence, the dynamical tracking target proposed in our prior work [24] does not align with this dynamic model. This significantly heightens the complexity of control design, particularly for achieving satisfactory trajectory tracking performance.
Inspired by the above analysis, this paper proposes a modified curvature tracking method to achieve satisfactory obstacle avoidance performance for the TWMR. The main contributions lie in the following.
(i)
To be more realistic, the motion laws of TWMRs are elucidated by specifically considering the steering angle of the front wheel as a control variable. Subsequently, a standardized method is proposed to derive the dynamic model by sufficiently using the obtained motion laws.
(ii)
A modified dynamical tracking target that fits the dynamic model very well is proposed to resolve the accurate trajectory tracking problem of the TWMR. It is composed of two parts: a yaw motion speed target and a tangential motion speed target. The former employs a suitable smooth function to design the yaw motion speed target, reducing cumulative location errors. The latter is reliant on actual yaw motion speed and allows for continuous adjustments. Consequently, the original motion issue transforms into a standard path tracking control problem, enabling TWMRs to adeptly follow desired paths by tracking the relative curvature.
The remainder of the paper is organized as following. First, a system description, kinematic model, and dynamic model are given in Section 2. Then, in Section 3, the desired trajectory is converted to a modified dynamical tracking target that matches well with the dynamic model. After that, control design for an obstacle avoidance problem is discussed in Section 4. As an application in Section 5, two numerical examples are developed to demonstrate the robustness and effectiveness of the proposed method. Finally, Section 6 draws the conclusions.

2. Modeling of the TWMR

In this section, kinematic and dynamic models are established for the TWMR. The corresponding parameters and variables of the TWMR are defined in Table 1.

2.1. System Description

Consider a TWMR comprising a steered front wheel and two differentially driven rear wheels, as depicted in Figure 1, where the two driven wheels are mounted on each side of the wheel axis. The driven torques are provided by the actuator acting on two rear wheels. In addition, in terms of the center rod P O , one end is anchored to the center of the wheel axis while the other end is flexibly connected to the steered wheel. Here, we assume that the rear wheels roll purely forward with no lateral slip for modeling the TWMR conveniently. As discussed in the introduction, the primary objective is to ensure accurate movement of the midpoint O along an identical path r ^ ( t ) = ( x ^ ( t ) , y ^ ( t ) ) while avoiding collisions with obstacles.

2.2. Kinematic and Dynamic Models

To facilitate the establishment of the dynamic model, we introduce a generalized coordinate denoted as q = ( x , y , θ , θ l , θ r ) T to represent the state of the TWMR.
At first, we examine the constraints imposed on the two driven rear wheels. Considering that the wheels roll purely forward without lateral slip, point O is subjected to the following two nonholonomic constraints:
x ˙ sin θ + y ˙ cos θ = 0 , v = x ˙ cos θ + y ˙ sin θ .
Here, the first equation indicates that the lateral speed is zero, and the second equation represents that the tangential speed at point O is the length of the tangential vector of the trajectory curve. Due to the nonholonomic constraint of zero lateral speed, point O moves along a smooth curve. We see from Equation (1) that the tangential speed at point O can be expressed as v = x ˙ 2 + y ˙ 2 . Thus, Equation (1) is equivalent to the following equation:
x ˙ = v cos θ , y ˙ = v sin θ .
As we know, the two driven wheels are subjected to two holonomic constraints, which are expressed as
v = r 2 ( θ ˙ l + θ ˙ r ) , θ ˙ = r d ( θ ˙ l θ ˙ r ) .
Equation (3) describes the relationship between the motion speed of point O and the rotation speed of two driven wheels. Solving for θ ˙ r and θ ˙ l from Equation (3), in terms of v and θ ˙ , one has
θ ˙ l = 1 r v + d 2 r θ ˙ , θ ˙ r = 1 r v d 2 r θ ˙ .
Now, we would like to consider the nonholonomic constraint of the steered front wheel. From Figure 1, we observe that point P represents both the front endpoint of the intermediate link and the center point of the front wheel simultaneously. For the former, it is subjected to the following nonholonomic constraint due to zero lateral velocity:
x ˙ P sin φ + y ˙ P cos φ = 0 .
From a geometric perspective, Equation (5) actually indicates that the speed of point P in the vertical direction of O P must be zero. Otherwise, if the nonholonomic constraint (5) is not satisfied, the lateral speed of point O must be nonzero, which would contradict the first nonholonomic constraint in Equation (1). For the latter, the velocity schematic of point P is shown in Figure 1. Along with Equation (5), one has
y ˙ P x ˙ P = tan φ = θ ˙ a v ,
which yields
θ ˙ = v tan φ a .
Equation (6) represents a significant nonholonomic constraint that describes the relation between the steering angle φ of the front wheel and the yaw motion of the TWMR. Substituting the first equation of (2) into Equation (6), one has
x ˙ tan φ + a θ ˙ cos θ = 0 .
So far, two holonomic and three nonholonomic constraints of the TWMR have been clarified. Substituting Equation (3) into Equation (1), and together with Equation (7), all nonholonomic constraints are expressed as
x ˙ sin θ + y ˙ cos θ = 0 , x ˙ cos θ + y ˙ sin θ r 2 ( θ ˙ l + θ ˙ r ) = 0 , x ˙ tan φ + a θ ˙ cos θ = 0 .
Equation (8) can be rewritten in a matrix form as
F ( q ) 3 × 5 q ˙ 5 × 1 = 0 3 × 1 ,
where
F ( q ) = sin θ cos θ 0 0 0 cos θ sin θ 0 r 2 r 2 tan φ 0 a cos θ 0 0 .
On the other hand, from Equations (2), (4) and (6), all constraints of the TWMR are summarized by
x ˙ = v cos θ , y ˙ = v sin θ , θ ˙ = v tan φ a , θ ˙ l = 1 r v + d 2 r θ ˙ , θ ˙ r = 1 r v d 2 r θ ˙ .
Similarly, the kinematic model (10) can be further rewritten in matrix form as
q ˙ 5 × 1 = S ( q ) 5 × 2 V 2 × 1 ,
where
S ( q ) = cos θ 0 sin θ 0 tan φ a 0 1 r d 2 r 1 r d 2 r , V = ( v , θ ˙ ) T .
In fact, by multiplying F ( q ) by S ( q ) , it is easy to obtain
F ( q ) 3 × 5 · S ( q ) 5 × 2 = 0 3 × 2 .
Now, we would like to derive the dynamic model of the TWMR. In accordance with the Euler–Lagrange equation of nonholonomic mechanical systems, the dynamic equation of the TWMR is formulated as
d d t ( L q ˙ ) L q = E ( q ) τ + F T ( q ) λ ,
where τ is the input motor-torque vector, E ( q ) represents the control input matrix, and λ denotes the so-called Lagrange multiplier. The terms τ and E ( q ) are respectively given by
τ = τ l τ r , E ( q ) = 0 0 0 0 1 0 0 0 0 0 0 1 T .
After being further calculated in Appendix A, the dynamic Equation (14) is rewritten in a state equation form as
v ˙ = I a r 2 φ ˙ sin φ Δ cos φ v + r a 2 cos 2 φ Δ u 2 , θ ¨ = ω ˙ = c u 1 ,
where Δ = M 1 a 2 r 2 cos 2 φ + 2 I w a 2 cos 2 φ + I r 2 sin 2 φ , c = r d I w , u 1 = T l T r , u 2 = T l + T r .
In terms of the dynamic Equation (16), u 1 is an intermediate control variable, and u 2 , φ are the final control variables. Since system (16) is decoupled, u 1 can be firstly designed by the second equation. Then, according to Equation (6), φ can be calculated by
θ ˙ = v tan φ a = 0 t c u 1 ( τ ) d τ .
Furthermore, substituting the obtained φ into the first equation of (16), u 2 can be finally designed by common control methods. In this way, when u 1 , u 2 are obtained, the desired control torques T l and T r for the TWMS can be derived from T l = 1 2 ( u 1 + u 2 ) and T r = 1 2 ( u 2 u 1 ) , respectively.
Remark 1. 
In practical applications, φ, as the central control variable, determines the yaw motion of the TWMS, which is achievable through steering wheel operation. U 2 = T l + T r provides the driving force for the TWMS by controlling the throttle. Once these two control variables are determined, the differential automatically achieves the different speeds of the left and right wheels under nonholonomic constraints. On this basis, precise motion control of the trailer is finally realized.

3. Dynamical Trajectory Planning for an Obstacle Avoidance Problem

As shown in Figure 2, the control task of the TWMR is to follow an ideal target trajectory curve while avoiding collisions with obstacles. Assume that the shape of the planning trajectory is approximated as a cycloid curve, as can be seen by the blue line. The parametric representation of the target trajectory curve is given by
r ^ ( t ) = ( x ^ ( t ) , y ^ ( t ) ) = ( α [ μ t + β sin ( μ t + β ) ] , α [ 1 cos ( μ t + β ) ] ) ,
where t [ 0 , 2 π 2 β μ ] , and α , β , μ are undetermined parameters of the non-standard cycloid curve that are determined by actual needs. In fact, β is a truncation parameter that affects the bending degree of the target curve. When β = 0 , r ^ ( t ) represents the whole cycloid curve. The proportion of r ^ ( t ) in the cycloid curve decreases with the increment of β . The term α affects the pulling up and down of the target curve. Since the height of h is fixed, h ^ increases with α . The term μ affects the pulling left and right of the target curve and increases with D. Note that dynamic Equation (16) only describes the relationship between the driven torques and the motion speeds. However, the actual motion task is to make the TWMR track a desired curve precisely. Thus, to match the dynamic equation, it is reasonable to transform the required path curve to a special speed target form.
For this purpose, assume that the smooth curve r ^ ( t ) = ( x ^ ( t ) , y ^ ( t ) ) is an actual motion trajectory of the TWMR. It follows from Equation (2) that
x ^ ˙ = v ^ cos θ ^ = x ^ ˙ 2 + y ^ ˙ 2 cos θ ^ , y ^ ˙ = v ^ sin θ ^ = x ^ ˙ 2 + y ^ ˙ 2 sin θ ^ ,
where θ ^ represents an angle between the tangent vector of the smooth curve and the horizontal axis. Taking the derivative of both sides of Equation (19) with respect to time reads
x ^ ¨ = 1 v ^ ( x ^ ˙ x ^ ¨ + y ^ ˙ y ^ ¨ ) cos θ ^ v ^ sin θ ^ θ ^ ˙ , y ^ ¨ = 1 v ^ ( x ^ ˙ x ^ ¨ + y ^ ˙ y ^ ¨ ) sin θ ^ + v ^ cos θ ^ θ ^ ˙ .
From Equations (19) and (20), we obtain
x ^ ˙ y ^ ¨ x ^ ¨ y ^ ˙ = ( x ^ ˙ 2 + y ^ ˙ 2 ) θ ^ ˙ .
Thus, the tangential speed and yaw motion speed of the TWMR are formulated as
v ^ = x ^ ˙ 2 + y ^ ˙ 2 , θ ^ ˙ = x ^ ˙ y ^ ¨ x ^ ¨ y ^ ˙ x ^ ˙ 2 + y ^ ˙ 2 = k ( s ( t ) ) v ^ , s ( t ) = 0 t v ^ ( τ ) d τ ,
where k ( s ( t ) ) = x ^ ˙ y ^ ¨ x ^ ¨ y ^ ˙ ( x ^ ˙ 2 + y ^ ˙ 2 ) 3 2 denotes the relative curvature of the actual motion path. From Equations (19) and (21), we see that the speed ( v ^ , θ ^ ˙ ) uniquely determines a motion trajectory curve r ^ ( t ) . However, a prescribed motion trajectory is able to be presented by different speed forms. Note that the inherent characteristic of a trajectory curve is the relative curvature. As a result, the motion trajectory r ^ ( t ) is presented in a general speed form as
v ^ = ϕ ˙ ( t ) , θ ^ ˙ = k ( ϕ ( t ) ) v ^ ( t ) .
Here, ϕ ˙ ( t ) is a one-by-one arbitrary smooth function that can be designed with regard to actual needs, and k ( ϕ ( t ) ) is the relative curvature function related to ϕ ˙ ( t ) . In an actual problem, both the initial and terminal speeds of the TWMR are zero. In this case, the tangential motion speed target is designed as
ϕ ˙ ( t ) = l γ 2 t e γ t ,
where l is the required motion distance, and γ is a proper parameter for the actual needs. Furthermore, the target curve (22) is described in a dynamical tracking target form as
v ^ = ϕ ˙ ( t ) , θ ^ ˙ = ω ^ = k ( s ) v ( t ) , s = 0 t v ( τ ) d τ ,
where v ( t ) is the actual tangential motion speed [24]. Actually, path tracking control based on Equation (24) is a curvature tracking method in essence. Hence, if the relative curvature is well tracked, the TWMR follows the ideal path accurately through dynamical tracking target (24).
In general, if Equation (24) is applied to track a smooth curve, the tangential motion speed target v ^ should be firstly tracked. Then, the output of the actual tangential motion speed is used to construct the yaw motion speed target ω ^ . Therefore, the control design of the tangential speed subsystem should be firstly performed; then, the yaw rotation subsystem is considered. However, according to the dynamic Equation (16) analyzed in Section 2, φ should be calculated from the intermediate control variable u 1 . This implies that the yaw rotation subsystem needs to be addressed first. Thus, the designed dynamic tracking target (24) cannot align with dynamic Equation (16).
To address this issue, the dynamical tracking target (24) is further converted to a modified form as
θ ^ ˙ = x ^ ˙ y ^ ¨ x ^ ¨ y ^ ˙ x ^ ˙ 2 + y ^ ˙ 2 ϕ ˙ ( t ) , v ^ = θ ˙ k ( s ( t ) ) , s ( t ) = 0 t v ( τ ) d τ ,
where ϕ ˙ ( t ) can be chosen flexibly for better tracking performance, and θ ˙ is the actual yaw motion speed. Obviously, the modified dynamical tracking target (25) is able to match the dynamic Equation (16) very well, while retaining the merits of tracking the relative curvature. To this end, the relative curvature is calculated by
k ( t ) = 1 4 α sin μ t + β 2 .
Moreover, the arc-length function of the trajectory curve is given by
s ( t ) = 0 t α μ 2 ( 1 cos ( μ τ + β ) d τ = 4 α cos β 2 4 α cos μ t + β 2 ,
and the whole length of the required path curve is expressed by
l = s ( 2 π 2 β μ ) = 8 α cos β 2 .
Solving sin μ t + β 2 from (27) and substituting it into (26) gives the desired relative curvature:
k ( s ) = 1 16 α 2 16 α 2 cos 2 β 2 s 2 + 8 α s cos β 2 .
In this way, the modified dynamical tracking target (25) for the target trajectory curve (18) is eventually designed as
θ ^ ˙ = x ^ ˙ y ^ ¨ x ^ ¨ y ^ ˙ x ^ ˙ 2 + y ^ ˙ 2 ϕ ˙ ( t ) = μ 2 ϕ ˙ ( t ) , v ^ = 16 α 2 16 α 2 cos 2 β 2 s 2 + 8 α cos β 2 s θ ˙ ,
where s [ 0 , 8 α cos β 2 ] . The maximum value of s represents the length of the entire trajectory curve, which is calculated by Equation (28).

4. Trajectory Tracking Control Design

In this section, we design two controllers ( u 2 , φ ) for dynamic Equation (16) to ensure precise tracking of the target curve by the TWMR through the dynamic tracking target (30). According to the analysis in Section 3, we initially address the yaw motion subsystem to derive the intermediate control u 1 and subsequently obtain the terminal control φ through Equation (17). Next, we consider the tangential speed subsystem to design the terminal control u 2 .
With regard to the yaw motion speed subsystem in the second equation of (16), we set X = x 1 x 2 = θ θ ˙ . Then, the yaw motion speed subsystem can be rewritten as
X ˙ = A X + B u 1 ,
here A = 0 1 0 0 and B = 0 c . Moreover, denoting X ^ = x ^ 1 x ^ 2 = θ ^ θ ^ ˙ , Y = y 1 y 2 = X X ^ and η : = A X ^ X ^ ˙ with θ ^ = 0 t θ ^ ˙ ( τ ) d τ , subsystem (31) is transformed into an error system as
Y ˙ = A Y + B u 1 + η .
When ϕ ˙ ( t ) = l γ 2 t e γ t is used to construct the yaw motion speed target, one can see that θ ^ ˙ approaches zero as t 0 , and its infinite integral is convergent. Accordingly, we reasonably introduce the following quadratic performance index with a large weight with regard to the yaw motion speed error, described by
J = 1 2 0 + Y T Q Y + u 1 T R u 1 d t .
In Equation (33), R and Q represent weight matrices. Due to the use of a dynamic target tracking concept, it is critical to minimize the speed tracking error (i.e., the corresponding curvature tracking error) to achieve precise tracking of the target curve. Therefore, for the parameter selection of the weight matrices discussed in Section 5, we prioritize assigning a relatively larger weight to the speed error to ensure precise tracking of the target curve.
Clearly, by employing an optimal controller u 1 , the yaw motion speed error would be small via minimizing the performance index J. On the basis of LQR optimal control theory, the intermediate control variable u 1 is designed by
u 1 = R 1 B T [ P Y + b ] .
Here, P and b are given by
PA A T P + P B R 1 B T P Q = 0 , b ˙ = [ A B R 1 B T P ] T b P η , b ( + ) = 0 .
The optimal controller (34) consists of two components. The first component, R 1 B T P Y , is the fundamental solution to the optimal control problem, where P is determined by the first equation of Equation (35) (i.e., the Riccati equation). The second component, R 1 B T b , is primarily used to account for the effect of the known disturbance η in system (32), where b satisfies the second equation of Equation (35).
Together with Equation (17), the terminal control variable φ is formulated by
φ = arctan ( a 0 t c u 1 ( τ ) d τ v ) .
Since the terminal control variable φ is obtained, the tangential motion speed subsystem in the first equation of (16) is taken into consideration. To address linearization errors and uncertainties, an integral sliding mode control method is introduced as follows [25,26]. We set u ¯ = I a r 2 φ ˙ sin φ Δ cos φ v + r a 2 cos 2 φ Δ u 2 and z = v v ^ . Then, the tangential speed subsystem is transformed into an error subsystem as
z ˙ = u ^ 2 ,
where u ^ 2 = u ¯ v ^ ˙ , and v ^ ˙ = k ( s ) ( 4 α cos β 2 s ) θ ˙ v 2 + 1 k ( s ) c u 1 . Since u ^ 2 is designed as u ^ 2 = k z , it is easy to obtain u ¯ = k z + v ^ ˙ . Then, the fundamental part of the integral sliding mode controller u 2 is obtained as
u 20 = Δ r a 2 cos 2 φ ( k z + v ^ ˙ + I a r 2 φ ˙ sin φ Δ cos φ ( z + v ^ ) ) .
Furthermore, to achieve satisfactory obstacle avoidance performance, the tracking target v ^ is chosen as the integral sliding mode manifold ϝ ( v ^ ( t ) ) = 0 , where the sliding mode function is expressed by
ϝ ( v ( t ) ) : = G [ v ( t ) v ( 0 ) ] G 0 t v ^ ˙ ( τ ) d τ ,
with G > 0 being an appropriate parameter. In this case, the switching control part of u 2 can be constructed by
u 21 = G 1 ( ε + G E v ^ v ) · sgn ( ϝ ( v ( t ) ) ) ,
where ε is a switching control parameter, and E is the maximum energy limit of the whole of the uncertainties. As a consequence, the robust controller u 2 is finally designed by
u 2 = u 20 + u 21 .
In summary, from Equations (36) and (41), the desired controllers ( φ , u 2 ) are eventually designted as
φ = arctan ( a c v 0 t R 1 B T [ P ( τ ) Y ( τ ) + b ( τ ) ] d τ ) , u 2 = Δ r a 2 cos 2 φ ( k z + v ^ ˙ + I a r 2 φ ˙ sin φ Δ cos φ ( z + v ^ ) ) G 1 ( ε + G E v ^ v ) · sgn ( ϝ ( v ) ) .
In light of the above analyses, the detailed control design process is as follows in Table 2:

5. Simulation Results

In this section, we present two groups of simulations to demonstrate the robustness and effectiveness of the proposed tracking control strategy. First, we perform a comparison of the trajectory tracking performance between dynamical and statical tracking targets. After that, we report the trajectory tracking precision of the proposed method when subjected to unstable speed error subsystems.
For numerical simulations, the relevant parameters and initial values of the TWMR are configured as M B = 3 kg , M w = 2 kg , I w = 0.01 kg m 2 , I B = 0.16 kg m 2 , I w d = 0.005 kg m 2 , a = 0.4 m , r = 0.1 m , d = 0.2 m , μ = 1 , g = 9.8 m / s 2 , D = 10 m , h = 5 m , h ^ = 1 m , R = 1 , Q = diag ( 1 , 100 ) . Thus, the dynamic Equation (16) is expressed as
v ˙ = 10.8 φ ˙ tan φ 37 cos 2 φ + 27 v + 53.33 cos 2 φ 37 cos 2 φ + 27 u 2 , θ ¨ = 50 u 1 .
Since the initial speed of the dynamic tracking target is zero, it would be better to stop the three-wheeled mobile robot before performing the motion control task. Therefore, we can take the initial conditions as θ ( 0 ) = 1.3147 , v ( 0 ) = 0 , and θ ˙ ( 0 ) = 0 .
According to the geometric relationship in Figure 2, the parameters α and β of the target curve (18) are determined by
α π α ( β sin β ) = D , 2 α α ( 1 cos β ) = h + h ^ ,
where h ^ is the safety margin. Solving the above equations gives
α = 3.2056 , β = 0.5121 .
Thus, the target trajectory curve (18) is approximately given by
r ^ ( t ) = ( 3.2056 · [ t + 0.5121 sin ( t + 0.5121 ) ] , 3.2056 · [ 1 cos ( t + 0.5121 ) ] ) .

5.1. A Simulation Comparison between the Dynamical and the Statical Tracking Targets

We first suppose that the yaw motion speed error subsystem is affected by an uncertain factor d 1 = ( 5 y 1 , 5 y 2 ) T related to the error state, and (32) becomes
Y ˙ = A Y + B u 1 + η + d 1 .
Meanwhile, the tangential motion speed error subsystem is assumed to be affected by an uncertain factor d 2 = 0.1 sin t + 2 z ( t ) . After that, the corresponding tangential motion speed error subsystem can be expressed as
z ˙ = u ^ 2 + d 2 .
For the switching control part u 21 , let G = 1 , ε = 0.1 , E = 0.2 and k = 600 . With these preparations, all the quantities required for φ and u 2 are available on hand.
According to Equation (30), the modified dynamical tracking target of the target trajectory curve (46) is designed by
θ ^ ˙ = 3 4 t e 1 2 t , v ^ = 10.5476 s 2 + 24.8089 s θ ˙ , s [ 0 , 24.8089 ] .
where s = 0 t v ( τ ) d τ . Meanwhile, the statical tracking target (21) is expressed by
v ^ = 2 α sin t + β 2 , θ ^ ˙ = k ( t ) v ^ ( t ) = 1 2 , t [ 0 , 2 π 2 β ] .
If (50) is adopted to design the controller (42), the initial speed error is non-zero and is described as
v ( 0 ) v ^ ( 0 ) = 2 α sin β 2 = 1.6239 , θ ˙ ( 0 ) θ ^ ˙ ( 0 ) = 1 2 . .
As shown in Figure 3a, the tangential motion and yaw motion speed errors are sufficiently small when employing such a dynamic tracking target. This results in minimal accumulative location errors for the TWMR. Consequently, the actual trajectory nearly coincides with the prescribed path, as depicted in Figure 3b. However, as observed in Figure 4a, the initial speed error is significant, and the accumulative location errors increase when using the static tracking target. Correspondingly, as shown in Figure 4b, the actual trajectory deviates significantly from the target trajectory curve.

5.2. Accurate Trajectory Tracking with Unstable Speed Error Subsystems

In fact, although the speed error subsystems are unstable, actual motion trajectory can closely follow the prescribed path via dynamical tracking target. In what follows, two types of uncertain factors acting on the yaw motion speed error subsystem are explored to further verify the influence of uncertainties on path tracking.
On the one hand, we denote the uncertainty as d 1 = d 11 = ( 12 y 1 , 11.2 y 2 ) T , which is an internal system uncertainty related to the error state. In this case, the yaw motion speed error subsystem becomes unstable when Equation (42) is applied. As shown in Figure 5, despite increasing speed errors, the actual motion trajectory can achieve satisfactory tracking performance if the tangential motion speed error remains sufficiently small. This implies that throughout the motion process, the actual location error remains manageable, and the speed error does not accumulate significantly due to the modified dynamic tracking target.
On the other hand, we choose the uncertainty as d 1 = d 12 = ( sin 4 t , cos 4 t ) T , which represents an external disturbance independent of the error state. Figure 6a shows that both tangential speed error and yaw motion speed error are not convergent to zero. However, as shown in Figure 6b, the TWMR can still effectively follow the prescribed path if the tangential motion speed error remains sufficiently small. Moreover, the corresponding control inputs of the TWMR is decribed yin Figure 7a,b.
In the following, we calculate the square error between the actual motion trajectory and the target trajectory to validate the effects of uncertainties on trajectory tracking error.
Let ( x a ( k ) , y a ( k ) ) denote a known point on the actual motion trajectory. For the same horizontal coordinate x a ( k ) , we define ( x a ( k ) , y ^ a ( k ) ) as the corresponding point on the target trajectory curve. Since x a ( k ) is known, according to Equation (46), we have
x a ( k ) = 3.2056 · [ t k + 0.5121 sin ( t k + 0.5121 ) ] , y ^ a ( k ) = 3.2056 · [ 1 cos ( t k + 0.5121 ) ] ) ,
in which case, y ^ a ( k ) can be calculated. Then, the square error between ( x a ( k ) , y a ( k ) ) and ( x a ( k ) , y ^ a ( k ) ) can be expressed by ( y a ( k ) y ^ a ( k ) ) 2 . As a consequence, the sum of these errors is denoted as J a = k = 1 N ( y a ( k ) y ^ a ( k ) ) 2 with time step h = 0.001 and N = 20 h . By calculating the trajectory tracking error in Figure 6b, we obtain J a = 6.3217 . Accordingly, the formula for calculating the average error per point is given by J a N = 6.3217 20000 < 0.02 . Therefore, the quantitative results show that by applying the proposed control strategy the TWMR can track the target path accurately despite various uncertainties.

6. Conclusions

The nonholonomic constraints and motion laws of a TWMR are comprehensively elucidated in this paper. Compared to existing related research, our study exhibits two distinctive features. (i) In order to agree well with engineering practices, the steering angle of the front wheel is considered as a control variable rather than a state variable. In this sense, the dynamic model becomes more complicated, and an intermediate control variable has to be introduced to determine the steering angle. (ii) By incorporating the relative curvature of a prescribed curve, a modified dynamical tracking target that matches well the dynamic model is constructed. On this basis, the original obstacle avoidance motion issue is transformed to a conventional tracking control problem. The core essence of the modified dynamic tracking target lies in the relative curvature, which significantly enhances trajectory tracking accuracy and robustness. Moreover, the yaw motion speed target can be chosen freely for better tracking performance in practical scenarios. Finally, simulation results indicate that trajectory deviation can be greatly reduced, even in the presence of unstable speed error subsystems.
The modified dynamical tracking target can be widely utilized to precisely track a given smooth trajectory curve for various nonholonomic mechanical structures such as robots, aircraft, and steamers. Moreover, the dynamical tracking target method can be further extended to a space trajectory curve by incorporating curve torsion, which is our future research.

Author Contributions

The authors X.W. and Y.Z. contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

The work was supported by the National Natural Science Foundation under grant (12162006). Also, it was supported by Cultivation Project of Yibin University (2022PY31), and the Guizhou Provincial Science and Technology Projects (ZK[2023] 037, [2023]119).

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Appendix A.1. The Detailed Derivation Process of the Dynamic Model

Since the system energy is mainly provided by the driven motors on the right and left rear wheels, all posture variables ( x , y , θ ) T are achieved by the coupling action of nonholonomic constraints. In what follows, we will calculate the system kinetic energy. As depicted in Figure 1, position relations between point O and the center points of the two rear wheels are described by
x w l = x + d 2 sin θ , y w l = y d 2 cos θ , x w r = x d 2 sin θ , y w r = y + d 2 cos θ .
Here, ( x w l , y w l ) and ( x w r , y w r ) represent the coordinates of the center points of the left and right rear wheels, respectively. Similarly, position relations among points C, P, and O are given by
x C = x + a 2 cos θ , y C = y + a 2 sin θ , x P = x + a cos θ , y P = y + a sin θ .
In the first place, the sum of the translational energy and rotational kinetic energy of the two rear wheels can be formulated as
T w = 1 2 M w ( v w l 2 + v w r 2 ) + 1 2 I w ( θ ˙ r 2 + θ ˙ l 2 ) + I w d θ ˙ 2 = M w ( x ˙ 2 + y ˙ 2 + d 2 4 θ ˙ 2 ) + 1 2 I w ( θ ˙ r 2 + θ ˙ l 2 ) + I w d θ ˙ 2 ,
where v w r , v w l are the tangential motion speeds of the center points of the right and left rear wheels, respectively. Moreover, the sum of the translational energy and rotational kinetic energy of the center rod can be presented as
T b = 1 2 M b [ x ˙ 2 + y ˙ 2 + a 2 4 θ ˙ 2 a θ ˙ ( x ˙ sin θ y ˙ cos θ ) ] + 1 2 I b θ ˙ 2 .
Furthermore, the sum of the translational energy and rotational kinetic energy of the front wheel can be expressed as
T f = 1 2 ( M w + I w r 2 ) [ x ˙ 2 + y ˙ 2 + a 2 θ ˙ 2 2 a θ ˙ ( x ˙ sin θ y ˙ cos θ ) ] .
Therefore, the Lagrangian function is obtained by
L = T w + T b + T f .
Substituting L into Equation (14), the dynamic equation is derived as
A ( q ) q ¨ + U ( q , q ˙ ) q ˙ = E ( q ) τ + F T ( q ) λ
where
A ( q ) = M 1 0 M 2 a sin θ 0 0 0 M 1 M 2 a cos θ 0 0 M 2 a sin θ M 2 a cos θ I 0 0 0 0 0 I w 0 0 0 0 0 I w , U ( q , q ˙ ) = 0 0 M 2 a θ ˙ cos θ 0 0 0 0 M 2 a θ ˙ sin θ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ,
with M 1 = M b + 3 M w + 1 r 2 I w , M 2 = 1 2 M b + M w + 1 r 2 I w , and I = I b + 2 I w d + a 2 4 M b + d 2 2 M w + a 2 M w + a 2 r 2 I w . Taking the derivative of both sides of Equation (11) with respect to time and substituting it into Equation (A6) results in
A ( q ) S ( q ) V ˙ + A ( q ) S ˙ ( q ) V + U ( q , q ˙ ) S ( q ) V = E ( q ) τ + F T ( q ) λ .
Multiplying Equation (A8) by S T ( q ) from the left, and together with Equation (13), yields
( S T ( q ) A ( q ) S ˙ ( q ) + S T ( q ) U ( q , q ˙ ) S ( q ) ) V + S T ( q ) A ( q ) S ( q ) V ˙ = S T ( q ) E ( q ) τ .
After being further calculated, the dynamic Equation (A9) is rewritten in a state equation form as
v ˙ = I a r 2 φ ˙ sin φ Δ cos φ v + r a 2 cos 2 φ Δ u 2 , θ ¨ = ω ˙ = c u 1 ,
where Δ = M 1 a 2 r 2 cos 2 φ + 2 I w a 2 cos 2 φ + I r 2 sin 2 φ , c = r d I w , u 1 = T l T r , u 2 = T l + T r .

References

  1. Guo, F.; Yang, H.; Wu, X. Model-based deep learning for low-cost IMU dead reckoning of wheeled mobile robot. IEEE Trans. Ind. Electron. 2023, 1, 7531–7541. [Google Scholar] [CrossRef]
  2. Zhao, L.; Li, J.; Li, H. Double-loop tracking control for a wheeled mobile robot with unmodeled dynamics along right angle roads. ISA Trans. 2023, 136, 525–534. [Google Scholar] [CrossRef] [PubMed]
  3. Parhi, D.R.; Deepak, B. Kinematic model of three wheeled mobile robot Journal of Mechanical Engineering Research. J. Mech. Eng. Res. 2011, 3, 307–318. [Google Scholar]
  4. He, W.; Li, Z.J.; Chen, C.L.P. A survey of human-centered intelligent robots: Issues and challenges. IEEE-CAA J. Autom. Sin. 2017, 4, 602–609. [Google Scholar] [CrossRef]
  5. Lin, H.Y.; Tsai, S.H.; Chen, K.Y. Design and implementation of the trajectory tracking and dynamic obstacle avoidance of wheeled mobile robot based on T-S fuzzy model. Int. J. Fuzzy Syst. 2023, 25, 2423–2438. [Google Scholar] [CrossRef]
  6. Saenz, A.; Santiban˜ez, V.; Bugarin, E. Velocity control of an omnidirectional wheeled mobile robot using computed voltage control with visual feedback: Experimental results. Int. J. Control. Autom. Syst. 2021, 19, 1089–1102. [Google Scholar] [CrossRef]
  7. Savaee, E.; Rahmani Hanzaki, A. A new algorithm for calibration of an omni-directional wheeled mobile robot based on effective kinematic parameters estimation. J. Intell. Robot. Syst. 2021, 101, 28. [Google Scholar] [CrossRef]
  8. Huang, J.; Chen, J.; Fang, H. An overview of recent progress in high-order nonholonomic chained system control and distributed coordination. J. Control. Decis. 2015, 2, 64–85. [Google Scholar] [CrossRef]
  9. Li, H.P.; Yan, W.S.; Shi, Y. A receding horizon stabilization approach to constrained nonholonomic systems in power form. Syst. Control. Lett. 2017, 99, 47–56. [Google Scholar] [CrossRef]
  10. Korayem, M.H.; Shafei, A.M. Motion equation of nonholonomic wheeled mobile robotic manipulator with revolute-prismatic joints using recursive gibbs-appell formulation. Appl. Math. Model. 2015, 39, 1701–1716. [Google Scholar] [CrossRef]
  11. Kalin´ski, K.J.; Buchholz, C. Mechatronic design of strongly nonlinear systems on a basis of three wheeled mobile platform. Mech. Syst. Signal Process. 2015, 52-53, 700–721. [Google Scholar] [CrossRef]
  12. Garcia-Sa´nchez, J.R.; Tavera-Mosqueda, S.; Silva-Ortigoza, R.; Herna´ndez-Guzma´n, V.M. A novel dynamic three-level tracking controller for mobile robots considering actuators and power stage subsystems: Experimental assessment. Sensors 2020, 20, 4959. [Google Scholar] [CrossRef]
  13. Silva-Ortigoza, J.R. Garcia-Sa´nchez, R.; Herna´ndez-Guzma´n, V.M.; Marquez-Sanchez, C. Trajectory tracking control for a differential drive wheeled mobile robot considering the dynamics related to the actuators and power stag. IEEE Lat. Am. Trans. 2016, 14, 657–664. [Google Scholar] [CrossRef]
  14. Zheng, Y.J.; Zheng, J.C.; Shao, K.; Zhao, H.; Man, Z.H.; Sun, Z. Adaptive fuzzy sliding mode control of uncertain nonholonomic wheeled mobile robot with external disturbance and actuator saturation. Inf. Sci. 2024, 663, 120303. [Google Scholar] [CrossRef]
  15. Airton, D.S.; Frank, S. A robotic walker based on a two-wheeled inverted pendulum. J. Intell. Robot. Syst. 2017, 86, 17–34. [Google Scholar]
  16. Guo, Z.Q.; Xu, J.X.; Lee, T.H. Design and implementation of a new sliding mode controller on an underactuated wheeled inverted pendulum. J. Frankl. Inst. 2014, 351, 2261–2282. [Google Scholar] [CrossRef]
  17. Furuichi, H.; Huang, J.; Fukuda, T.; Matsuno, T. Switching dynamic modeling and driving stability analysis of three-wheeled narrow tilting vehicle. IEEE-ASME Trans. Mechatronics 2014, 19, 1309–1322. [Google Scholar] [CrossRef]
  18. Chatterjee, M.; Kale, M.; Chaudhari, B. A dynamic stability control for electric narrow tilting three wheeled vehicle using integrated multivariable controller. Transp. Res. Part Transp. Environ. 2019, 66, 58–75. [Google Scholar] [CrossRef]
  19. Hacene, N.; Mendil, B. Fuzzy behavior-based control of three wheeled omnidirectional mobile robot. Int. J. Autom. Comput. 2019, 16, 163–185. [Google Scholar] [CrossRef]
  20. Viet, T.D.; Doan, P.T.; Hung, N. Tracking control of a three-wheeled omnidirectional mobile manipulator system with disturbance and friction. J. Mech. Sci. Technol. 2012, 26, 2197–2211. [Google Scholar] [CrossRef]
  21. Guo, Y.J.; Yu, L.; Xu, J.M. Robust finite-time trajectory tracking control of wheeled mobile robots with parametric uncertainties and disturbances. J. Syst. Sci. Complex. 2019, 32, 1358–1374. [Google Scholar] [CrossRef]
  22. Yang, H.J.; Wang, S.Z.; Zuo, Z.Q.; Li, P. Trajectory tracking for a wheeled mobile robot with an omnidirectional wheel on uneven ground. IET Control. Theory Appl. 2020, 14, 921–929. [Google Scholar] [CrossRef]
  23. Zhou, Y.S.; Wang, Z.H. Motion controller design of wheeled inverted pendulum with an input delay via optimal control theory. J. Optim. Theory Appl. 2016, 168, 625–645. [Google Scholar] [CrossRef]
  24. Zhou, Y.S.; Wang, Z.H.; Chung, K.K. Turning motion control design of a two-wheeled inverted pendulum using curvature tracking and optimal control. J. Optim. Theory Appl. 2019, 181, 634–652. [Google Scholar] [CrossRef]
  25. Zhou, Y.S.; Wang, Z.H. Robust motion control of a two-wheeled inverted pendulum with an input delay based on optimal integral sliding mode manifold. Nonlinear Dyn. 2016, 85, 2065–2074. [Google Scholar] [CrossRef]
  26. Mobayen, S. Chaos synchronization of uncertain chaotic systems using composite nonlinear feedback based integral sliding mode control. ISA Trans. 2018, 77, 100–111. [Google Scholar] [CrossRef]
Figure 1. The model of the investigated TWMR.
Figure 1. The model of the investigated TWMR.
Axioms 13 00389 g001
Figure 2. Schematic diagram of the obstacle avoidance motion of the TWMR.
Figure 2. Schematic diagram of the obstacle avoidance motion of the TWMR.
Axioms 13 00389 g002
Figure 3. Actual motion speed errors and motion trajectory of the TWMR via the dynamical tracking target. (a) Actual motion speed errors of the TWMR; (b) actual motion trajectory curve of the TWMR.
Figure 3. Actual motion speed errors and motion trajectory of the TWMR via the dynamical tracking target. (a) Actual motion speed errors of the TWMR; (b) actual motion trajectory curve of the TWMR.
Axioms 13 00389 g003
Figure 4. Actual motion speed errors and motion trajectory of the TWMR via the statical target. (a) Actual motion speed errors of the TWMR; (b) actual motion trajectory curve of the TWMR.
Figure 4. Actual motion speed errors and motion trajectory of the TWMR via the statical target. (a) Actual motion speed errors of the TWMR; (b) actual motion trajectory curve of the TWMR.
Axioms 13 00389 g004
Figure 5. Actual motion speed errors and motion trajectory of the TWMR under the internal uncertainty d 11 . (a) Actual motion speed errors of the TWMR; (b) Actual motion trajectory curve of the TWMR.
Figure 5. Actual motion speed errors and motion trajectory of the TWMR under the internal uncertainty d 11 . (a) Actual motion speed errors of the TWMR; (b) Actual motion trajectory curve of the TWMR.
Axioms 13 00389 g005
Figure 6. Actual motion speed errors and motion trajectory of the TWMR under the external disturbance d 12 . (a) Actual motion speed errors of the TWMR; (b) Actual motion trajectory curve of the TWMR.
Figure 6. Actual motion speed errors and motion trajectory of the TWMR under the external disturbance d 12 . (a) Actual motion speed errors of the TWMR; (b) Actual motion trajectory curve of the TWMR.
Axioms 13 00389 g006
Figure 7. Angle control input and torque control input of the TWMR under external disturbance d 12 . (a) Angle control input of the TWMR; (b) torque control input of the TWMR.
Figure 7. Angle control input and torque control input of the TWMR under external disturbance d 12 . (a) Angle control input of the TWMR; (b) torque control input of the TWMR.
Axioms 13 00389 g007
Table 1. Parameters and variables of the TWMR.
Table 1. Parameters and variables of the TWMR.
NotationDefinition
T l ,   T r Torques provided by actuators acting on the left and right rear wheels
θ l ,   θ r Rotation angles of the left and right wheels
x,  yThe coordinates of point O in the X Y plane
x C ,   y C The coordinates of the center point C in the X Y plane
x P ,   y P The coordinates of point P in the X Y plane
φ Steering angle of the front wheel
θ Heading angle of the TWMR
vTangential speed of the TWMR, with v : = x ˙ 2 + y ˙ 2
M w Mass of the wheel
rRadius of the wheel
M b Mass of the center rod
gGravity acceleration
aLength of the center rod
dDistance between the two rear wheels along the axle center
I w Moment of inertia of the wheel along the wheel axis direction
I w d Moment of inertia of the wheel about Z-axis through the center of the wheel
I b Moment of inertia of the center rod about the Z-axis through point O
Table 2. The control Design Process of the TWMR.
Table 2. The control Design Process of the TWMR.
StepProcess
Step 1.Design an appropriate target trajectory curve r ^ ( t ) based on the obstacle avoidance motion task.
Step 2.Convert the target trajectory curve r ^ ( t ) into the dynamic speed target form (30).
Step 3.Consider the yaw rotation subsystem in the dynamic equation, design an intermediate variable controller u 1 based on the designed yaw rotation speed target, and obtain the steering angle controller φ from Equation (17).
Step 4.Substitute φ into the forward speed subsystem and design the torque controller u 2 in conjunction with the forward speed target, which includes the relative curvature of the desired trajectory.
Step 5.Substitute Equation (42) into the dynamic equation for numerical simulation to verify the effectiveness and robustness of the designed controllers.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wen, X.; Zhou, Y. Precise Obstacle Avoidance Movement for Three-Wheeled Mobile Robots: A Modified Curvature Tracking Method. Axioms 2024, 13, 389. https://doi.org/10.3390/axioms13060389

AMA Style

Wen X, Zhou Y. Precise Obstacle Avoidance Movement for Three-Wheeled Mobile Robots: A Modified Curvature Tracking Method. Axioms. 2024; 13(6):389. https://doi.org/10.3390/axioms13060389

Chicago/Turabian Style

Wen, Xiangrong, and Yusheng Zhou. 2024. "Precise Obstacle Avoidance Movement for Three-Wheeled Mobile Robots: A Modified Curvature Tracking Method" Axioms 13, no. 6: 389. https://doi.org/10.3390/axioms13060389

APA Style

Wen, X., & Zhou, Y. (2024). Precise Obstacle Avoidance Movement for Three-Wheeled Mobile Robots: A Modified Curvature Tracking Method. Axioms, 13(6), 389. https://doi.org/10.3390/axioms13060389

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