Search results

1 – 10 of 723
Article
Publication date: 17 February 2023

Shengqian Li and Xiaofan Zhang

An active disturbance rejection controller (ADRC) based on model compensation is proposed in this paper. The method should first be taken a nominal model of the robot to…

Abstract

Purpose

An active disturbance rejection controller (ADRC) based on model compensation is proposed in this paper. The method should first be taken a nominal model of the robot to compensate. Subsequently, the uncertain external disturbance is estimated and compensated is used an expansion state observer (ESO) in real time, which can reduce the estimating range of observation for ESO. The purpose of this paper is to suggest a novel method to improve the system tracking performance, as well as the dynamic and static performance index.

Design/methodology/approach

A welding robot is a complicated system with uncertainty, time-varying, strong coupling and a nonlinear system; it is more complex as if it is used in an underwater environment, and it is difficult to establish an accurate dynamic model for an underwater welding robot. Aiming at the tracking control of an underwater welding robot, it is difficult to achieve the control performance requirements by the conventional proportional integral derivative method to realize automatic tracking of the seam.

Findings

The simulation experiment is carried out by MATLAB/Simulink, and the application experiment is recorded. The experimental results show that the control method is correct and effective, and the system’s tracking performance is stable, and the robustness and tracking accuracy of the system are also improved.

Originality/value

The seam gets plumper and smoother, with better continuity and no undercut phenomenon.

Details

Industrial Robot: the international journal of robotics research and application, vol. 50 no. 5
Type: Research Article
ISSN: 0143-991X

Keywords

Article
Publication date: 13 March 2024

Ziyuan Ma, Huajun Gong and Xinhua Wang

The purpose of this paper is to construct an event-triggered finite-time fault-tolerant formation tracking controller, which can achieve a time-varying formation control for…

Abstract

Purpose

The purpose of this paper is to construct an event-triggered finite-time fault-tolerant formation tracking controller, which can achieve a time-varying formation control for multiple unmanned aerial vehicles (UAVs) during actuator failures and external perturbations.

Design/methodology/approach

First, this study developed the formation tracking protocol for each follower using UAV formation members, defining the tracking inaccuracy of the UAV followers’ location. Subsequently, this study designed the multilayer event-triggered controller based on the backstepping method framework within finite time. Then, considering the actuator failures, and added self-adaptive thought for fault-tolerant control within finite time, the event-triggered closed-loop system is subsequently shown to be a finite-time stable system. Furthermore, the Zeno behavior is analyzed to prevent infinite triggering instances within a finite time. Finally, simulations are conducted with external disturbances and actuator failure conditions to demonstrate formation tracking controller performance.

Findings

It achieves improved performance in the presence of external disturbances and system failures. Combining limited-time adaptive control and event triggering improves system stability, increase robustness to disturbances and calculation efficiency. In addition, the designed formation tracking controller can effectively control the time-varying formation of the leader and followers to complete the task, and by adding a fixed-time observer, it can effectively compensate for external disturbances and improve formation control accuracy.

Originality/value

A formation-following controller is designed, which can handle both external disturbances and internal actuator failures during formation flight, and the proposed method can be applied to a variety of formation control scenarios and does not rely on a specific type of UAV or communication network.

Details

Aircraft Engineering and Aerospace Technology, vol. 96 no. 3
Type: Research Article
ISSN: 1748-8842

Keywords

Article
Publication date: 25 July 2023

Gerasimos G. Rigatos, Masoud Abbaszadeh, Bilal Sari and Jorge Pomares

A distinctive feature of tilt-rotor UAVs is that they can be fully actuated, whereas in fixed-angle rotor UAVs (e.g. common-type quadrotors, octorotors, etc.), the associated…

Abstract

Purpose

A distinctive feature of tilt-rotor UAVs is that they can be fully actuated, whereas in fixed-angle rotor UAVs (e.g. common-type quadrotors, octorotors, etc.), the associated dynamic model is characterized by underactuation. Because of the existence of more control inputs, in tilt-rotor UAVs, there is more flexibility in the solution of the associated nonlinear control problem. On the other side, the dynamic model of the tilt-rotor UAVs remains nonlinear and multivariable and this imposes difficulty in the drone's controller design. This paper aims to achieve simultaneously precise tracking of trajectories and minimization of energy dissipation by the UAV's rotors. To this end elaborated control methods have to be developed.

Design/methodology/approach

A solution of the nonlinear control problem of tilt-rotor UAVs is attempted using a novel nonlinear optimal control method. This method is characterized by computational simplicity, clear implementation stages and proven global stability properties. At the first stage, approximate linearization is performed on the dynamic model of the tilt-rotor UAV with the use of first-order Taylor series expansion and through the computation of the system's Jacobian matrices. This linearization process is carried out at each sampling instance, around a temporary operating point which is defined by the present value of the tilt-rotor UAV's state vector and by the last sampled value of the control inputs vector. At the second stage, an H-infinity stabilizing controller is designed for the approximately linearized model of the tilt-rotor UAV. To find the feedback gains of the controller, an algebraic Riccati equation is repetitively solved, at each time-step of the control method. Lyapunov stability analysis is used to prove the global stability properties of the control scheme. Moreover, the H-infinity Kalman filter is used as a robust observer so as to enable state estimation-based control. The paper's nonlinear optimal control approach achieves fast and accurate tracking of reference setpoints under moderate variations of the control inputs. Finally, the nonlinear optimal control approach for UAVs with tilting rotors is compared against flatness-based control in successive loops, with the latter method to be also exhibiting satisfactory performance.

Findings

So far, nonlinear model predictive control (NMPC) methods have been of questionable performance in treating the nonlinear optimal control problem for tilt-rotor UAVs because NMPC's convergence to optimum depends often on the empirical selection of parameters while also lacking a global stability proof. In the present paper, a novel nonlinear optimal control method is proposed for solving the nonlinear optimal control problem of tilt rotor UAVs. Firstly, by following the assumption of small tilting angles, the state-space model of the UAV is formulated and conditions of differential flatness are given about it. Next, to implement the nonlinear optimal control method, the dynamic model of the tilt-rotor UAV undergoes approximate linearization at each sampling instance around a temporary operating point which is defined by the present value of the system's state vector and by the last sampled value of the control inputs vector. The linearization process is based on first-order Taylor series expansion and on the computation of the associated Jacobian matrices. The modelling error, which is due to the truncation of higher-order terms from the Taylor series, is considered to be a perturbation that is asymptotically compensated by the robustness of the control scheme. For the linearized model of the UAV, an H-infinity stabilizing feedback controller is designed. To select the feedback gains of the H-infinity controller, an algebraic Riccati equation has to be repetitively solved at each time-step of the control method. The stability properties of the control scheme are analysed with the Lyapunov method.

Research limitations/implications

There are no research limitations in the nonlinear optimal control method for tilt-rotor UAVs. The proposed nonlinear optimal control method achieves fast and accurate tracking of setpoints by all state variables of the tilt-rotor UAV under moderate variations of the control inputs. Compared to past approaches for treating the nonlinear optimal (H-infinity) control problem, the paper's approach is applicable also to dynamical systems which have a non-constant control inputs gain matrix. Furthermore, it uses a new Riccati equation to compute the controller's gains and follows a novel Lyapunov analysis to prove global stability for the control loop.

Practical implications

There are no practical implications in the application of the nonlinear optimal control method for tilt-rotor UAVs. On the contrary, the nonlinear optimal control method is applicable to a wider class of dynamical systems than approaches based on the solution of state-dependent Riccati equations (SDRE). The SDRE approaches can be applied only to dynamical systems which can be transformed to the linear parameter varying (LPV) form. Besides, the nonlinear optimal control method performs better than nonlinear optimal control schemes which use approximation of the solution of the Hamilton–Jacobi–Bellman equation by Galerkin series expansions. The stability properties of the Galerkin series expansion-based optimal control approaches are still unproven.

Social implications

The proposed nonlinear optimal control method is suitable for using in various types of robots, including robotic manipulators and autonomous vehicles. By treating nonlinear control problems for complicated robotic systems, the proposed nonlinear optimal control method can have a positive impact towards economic development. So far the method has been used successfully in (1) industrial robotics: robotic manipulators and networked robotic systems. One can note applications to fully actuated robotic manipulators, redundant manipulators, underactuated manipulators, cranes and load handling systems, time-delayed robotic systems, closed kinematic chain manipulators, flexible-link manipulators and micromanipulators and (2) transportation systems: autonomous vehicles and mobile robots. Besides, one can note applications to two-wheel and unicycle-type vehicles, four-wheel drive vehicles, four-wheel steering vehicles, articulated vehicles, truck and trailer systems, unmanned aerial vehicles, unmanned surface vessels, autonomous underwater vessels and underactuated vessels.

Originality/value

The proposed nonlinear optimal control method is a novel and genuine result and is used for the first time in the dynamic model of tilt-rotor UAVs. The nonlinear optimal control approach exhibits advantages against other control schemes one could have considered for the tilt-rotor UAV dynamics. For instance, (1) compared to the global linearization-based control schemes (such as Lie algebra-based control or flatness-based control), it does not require complicated changes of state variables (diffeomorphisms) and transformation of the system's state-space description. Consequently, it also avoids inverse transformations which may come against singularity problems, (2) compared to NMPC, the proposed nonlinear optimal control method is of proven global stability and the convergence of its iterative search for an optimum does not depend on initialization and controller's parametrization, (3) compared to sliding-mode control and backstepping control the application of the nonlinear optimal control method is not constrained into dynamical systems of a specific state-space form. It is known that unless the controlled system is found in the input–output linearized form, the definition of the associated sliding surfaces is an empirical procedure. Besides, unless the controlled system is found in the backstepping integral (triangular) form, the application of backstepping control is not possible, (4) compared to PID control, the nonlinear optimal control method is of proven global stability and its performance is not dependent on heuristics-based selection of parameters of the controller and (5) compared to multiple-model-based optimal control, the nonlinear optimal control method requires the computation of only one linearization point and the solution of only one Riccati equation.

Details

International Journal of Intelligent Unmanned Systems, vol. 12 no. 1
Type: Research Article
ISSN: 2049-6427

Keywords

Article
Publication date: 12 September 2023

Gerasimos G. Rigatos, Masoud Abbaszadeh, Pierluigi Siano and Jorge Pomares

Permanent magnet synchronous spherical motors can have wide use in robotics and industrial automation. They enable three-DOF omnidirectional motion of their rotor. They are…

Abstract

Purpose

Permanent magnet synchronous spherical motors can have wide use in robotics and industrial automation. They enable three-DOF omnidirectional motion of their rotor. They are suitable for several applications, such as actuation in robotics, traction in electric vehicles and use in several automation systems. Unlike conventional synchronous motors, permanent magnet synchronous spherical motors consist of a fixed inner shell, which is the stator, and a rotating outer shell, which is the rotor. Their dynamic model is multivariable and strongly nonlinear. The treatment of the associated control problem is important.

Design/methodology/approach

In this paper, the multivariable dynamic model of permanent magnet synchronous spherical motors is analysed, and a nonlinear optimal (H-infinity) control method is developed for it. Differential flatness properties are proven for the spherical motors’ state-space model. Next, the motors’ state-space description undergoes approximate linearization with the use of first-order Taylor series expansion and through the computation of the associated Jacobian matrices. The linearization process takes place at each sampling instance around a time-varying operating point, which is defined by the present value of the motors’ state vector and by the last sampled value of the control input vector. For the approximately linearized model of the permanent magnet synchronous spherical motors, a stabilizing H-infinity feedback controller is designed. To compute the controller’s gains, an algebraic Riccati equation has to be repetitively solved at each time-step of the control algorithm. The global stability properties of the control scheme are proven through Lyapunov analysis. Finally, the performance of the nonlinear optimal control method is compared against a flatness-based control approach implemented in successive loops.

Findings

Due to the nonlinear and multivariable structure of the state-space model of spherical motors, the solution of the associated nonlinear control problem is a nontrivial task. In this paper, a novel nonlinear optimal (H-infinity) control approach is proposed for the dynamic model of permanent magnet synchronous spherical motors. The method is based on approximate linearization of the motor’s state-space model with the use of first-order Taylor series expansion and the computation of the associated Jacobian matrices. Furthermore, the paper has introduced a different solution to the nonlinear control problem of the permanent magnet synchronous spherical motor, which is based on flatness-based control implemented in successive loops.

Research limitations/implications

The presented control approaches do not exhibit any limitations, but on the contrary, they have specific advantages. In comparison to global linearization-based control schemes (such as Lie-algebra-based control), they do not make use of complicated changes of state variables (diffeomorphisms) and transformations of the system's state-space description. The computed control inputs are applied directly to the initial nonlinear state-space model of the permanent magnet spherical motor without the intervention of inverse transformations and thus without coming against the risk of singularities.

Practical implications

The motion control problem of spherical motors is nontrivial because of the complicated nonlinear and multivariable dynamics of these electric machines. So far, there have been several attempts to apply nonlinear feedback control to permanent magnet-synchronous spherical motors. However, due to the model’s complexity, few results exist about the associated nonlinear optimal control problem. The proposed nonlinear control methods for permanent magnet synchronous spherical motors make more efficient, precise and reliable the use of such motors in robotics, electric traction and several automation systems.

Social implications

The treated research topic is central for robotic and industrial automation. Permanent magnet synchronous spherical motors are suitable for several applications, such as actuation in robotics, traction in electric vehicles and use in several automation systems. The solution of the control problem for the nonlinear dynamic model of permanent magnet synchronous spherical motors has many industrial applications and therefore contributes to economic growth and development.

Originality/value

The proposed nonlinear optimal control method is novel compared to past attempts to solve the optimal control problem for nonlinear dynamical systems. Unlike past approaches, in the new nonlinear optimal control method, linearization is performed around a temporary operating point, which is defined by the present value of the system's state vector and by the last sampled value of the control inputs vector and not at points that belong to the desirable trajectory (setpoints). Besides, the Riccati equation which is used for computing the feedback gains of the controller is new, and so is the global stability proof for this control method. Compared to nonlinear model predictive control, which is a popular approach for treating the optimal control problem in industry, the new nonlinear optimal (H-infinity) control scheme is of proven global stability, and the convergence of its iterative search for the optimum does not depend on initial conditions and trials with multiple sets of controller parameters. It is also noteworthy that the nonlinear optimal control method is applicable to a wider class of dynamical systems than approaches based on the solution of state dependent Riccati equations (SDRE). The SDRE approaches can be applied only to dynamical systems which can be transformed into the linear parameter varying form. Besides, the nonlinear optimal control method performs better than nonlinear optimal control schemes, which use approximation of the solution of the Hamilton–Jacobi–Bellman equation by Galerkin series expansions. Furthermore, the second control method proposed in this paper, which is flatness-based control in successive loops, is also novel and demonstrates substantial contribution to nonlinear control for robotics and industrial automation.

Article
Publication date: 5 December 2023

Zhirui Zhao, Lina Hao, Guanghong Tao, Hongjun Liu and Lihua Shen

This study discusses the tracking trajectory issue of the exoskeleton under the bounded disturbance and designs an useful tracking trajectory control method to solve it. By using…

124

Abstract

Purpose

This study discusses the tracking trajectory issue of the exoskeleton under the bounded disturbance and designs an useful tracking trajectory control method to solve it. By using the proposed control method, the tracking error can be successfully convergence to the assigned boundary. Meanwhile, the chattering effect caused by the actuators is already reduced, and the tracking performance of the pneumatic artificial muscles (PAMs) elbow exoskeleton is improved effectively.

Design/methodology/approach

A prescribed performance sliding mode control method was developed in this study to fulfill the joint position tracking trajectory task on the elbow exoskeleton driven by two PAMs. In terms of the control structure, a dynamic model was built by conforming to the adaptive law to compensate for the time variety and uncertainty exhibited by the system. Subsequently, a super-twisting algorithm-based second-order sliding mode control method was subjected to the exoskeleton under the boundedness of external disturbance. Moreover, the prescribed performance control method exhibits a smooth prescribed function with an error transformation function to ensure the tracking error can be finally convergent to the pre-designed requirement.

Findings

From the theoretical perspective, the stability of the control method was verified through Lyapunov synthesis. On that basis, the tracking performance of the proposed control method was confirmed through the simulation and the manikin model experiment.

Originality/value

As revealed by the results of this study, the proposed control method sufficiently applies to the PAMs elbow exoskeleton for tracking trajectory, which means it has potential application in the actual robot-assisted passive rehabilitation tasks.

Details

Industrial Robot: the international journal of robotics research and application, vol. 51 no. 1
Type: Research Article
ISSN: 0143-991X

Keywords

Article
Publication date: 23 January 2024

Li Li, Hui Ye and Xiaohua Meng

Considering the unmeasurable states of the systems and the previewed reference signal, a novel fuzzy observer-based preview controller, which is a mixed controller of the fuzzy…

Abstract

Purpose

Considering the unmeasurable states of the systems and the previewed reference signal, a novel fuzzy observer-based preview controller, which is a mixed controller of the fuzzy observer-based controller, fuzzy integrator and preview controller, is considered to address the tracking control problem.

Design/methodology/approach

The authors employ an augmentation technique to construct an augmented error system for uncertain T-S fuzzy discrete-time systems with time-varying uncertainties. Additionally, the authors obtain the corresponding linear matrix inequality (LMI) conditions for designing the preview controller.

Findings

This paper discusses the preview tracking problem for nonlinear systems. First, considering the unmeasurable states of the systems and the previewed reference signal, a novel fuzzy observer-based preview controller, which is a mixed controller of the fuzzy observer-based controller, fuzzy integrator, and preview controller, is considered to address the tracking control problem. Then, using the fuzzy Lyapunov functional with the linear matrix inequality (LMI) technique, new sufficient conditions for the asymptotic stability of the augmented system are derived by applying the LMI technique. The preview controller and fuzzy observer can be designed in one step. Finally, a numerical example is used to illustrate the effectiveness of the results.

Originality/value

An augmented error system is successfully constructed by the state augmentation approach. A novel preview controller is designed to address the tracking control problem. The preview controller and fuzzy observer can be designed in one step.

Details

Engineering Computations, vol. 41 no. 1
Type: Research Article
ISSN: 0264-4401

Keywords

Article
Publication date: 4 October 2022

Mohammad Bajelani, Morteza Tayefi and Man Zhu

This study aims to minimize the risk of costly failures of flight tests during the path tracking control design, and a noble approach has been proposed in this study to put the…

Abstract

Purpose

This study aims to minimize the risk of costly failures of flight tests during the path tracking control design, and a noble approach has been proposed in this study to put the whole vehicle-in-the simulation loop. Working with the real system is essential for developing intelligent and data-driven controllers for multirotor drones which needs learning the drones' nonlinear complicated dynamics. The vehicle-in-the-loop (VIL) platform developed in this paper is a safe and effective solution to deal with this problem.

Design/methodology/approach

To avoid risky flight test during controller design, the multirotor is hinged to a shaft that allows the multirotor's angular motion but restricts translational motion. The test-bed includes the real system attitude dynamics and the simulation of the position dynamics to model the complete flight based on real-time reactions of the vehicle. For the authors' case study, a hexacopter angular motion provides the real-time attitude data in translational motion simulation loop. To test the set-up, a proportional-integral-derivative (PID) and a brain emotional learning-based intelligent controller (BELBIC) is implemented for tracking of circle and 8-shape flight trajectories.

Findings

The results show that the platform helps the intelligent controller to learn the system dynamics without worrying about the failure in the early stages of the design and in the real-world flight test. Although the hexacopter translational dynamics is modeled in simulation, the authors still have highly accurate attitude dynamics matching the requirement of the control loop design. The comparison of the two controllers also shows that the performance of BELBIC is better than PID in this test.

Originality/value

The research background is reviewed in the introduction section. The other sections are originally developed in this paper.

Details

International Journal of Intelligent Unmanned Systems, vol. 11 no. 4
Type: Research Article
ISSN: 2049-6427

Keywords

Article
Publication date: 7 June 2023

Bingwei Gao, Wei Zhang, Lintao Zheng and Hongjian Zhao

The purpose of this paper is to design a third-order linear active disturbance rejection controller (LADRC) to improve the response characteristics and robustness of the…

Abstract

Purpose

The purpose of this paper is to design a third-order linear active disturbance rejection controller (LADRC) to improve the response characteristics and robustness of the electrohydraulic servo system.

Design/methodology/approach

The LADRC was designed by replacing the nonlinear functions in each part of ADRC with linear functions or linear combinations, and the parameters of each part of the LADRC were connected with their bandwidth through the pole configuration method to reduce the required tuning parameters, and used an improved grey wolf optimizer to tune the LADRC parameters.

Findings

The anti-interference control simulation and experiment on the LADRC, ADRC and proportion integration differentiation (PID) were carried out to test the robustness, anti-interference ability and superiority of the designed LADRC. The simulation and experiment results showed that the LADRC control and anti-interference control had excellent performance, and because of its simple structure and fewer parameters, LADRC was easier to implement and had a better control effect and anti-interference.

Originality/value

For the problems of parameter perturbation, unknown interference and inaccurate model in the electrohydraulic position servo system, the designed third-order LADRC has good tracking accuracy and anti-interference, has few parameters and is conducive to promotion.

Details

Robotic Intelligence and Automation, vol. 43 no. 3
Type: Research Article
ISSN: 2754-6969

Keywords

Article
Publication date: 23 October 2023

Yerui Fan, Yaxiong Wu and Jianbo Yuan

This study aims to improve the muscle model control performance of a tendon-driven musculoskeletal system (TDMS) to overcome disadvantages such as multisegmentation and strong…

Abstract

Purpose

This study aims to improve the muscle model control performance of a tendon-driven musculoskeletal system (TDMS) to overcome disadvantages such as multisegmentation and strong coupling. An adaptive network controller (ANC) with a disturbance observer is established to reduce the modeling error of the musculoskeletal model and improve its antidisturbance ability.

Design/methodology/approach

In contrast to other control technologies adopted for musculoskeletal humanoids, which use geometric relationships and antagonist inhibition control, this study develops a method comprising of three parts. (1) First, a simplified musculoskeletal model is constructed based on the Taylor expansion, mean value theorem and Lagrange–d’Alembert principle to complete the decoupling of the muscle model. (2) Next, for this simplified musculoskeletal model, an adaptive neuromuscular controller is designed to acquire the muscle-activation signal and realize stable tracking of the endpoint of the muscle-driven robot relative to the desired trajectory in the TDMS. For the ANC, an adaptive neural network controller with a disturbance observer is used to approximate dynamical uncertainties. (3) Using the Lyapunov method, uniform boundedness of the signals in the closed-loop system is proved. In addition, a tracking experiment is performed to validate the effectiveness of the adaptive neuromuscular controller.

Findings

The experimental results reveal that compared with other control technologies, the proposed design techniques can effectively improve control accuracy. Moreover, the proposed controller does not require extensive considerations of the geometric and antagonistic inhibition relationships, and it demonstrates anti-interference ability.

Originality/value

Musculoskeletal robots with humanoid structures have attracted considerable attention from numerous researchers owing to their potential to avoid danger for humans and the environment. The controller based on bio-muscle models has shown great performance in coordinating the redundant internal forces of TDMS. Therefore, adaptive controllers with disturbance observers are designed to improve the immunity of the system and thus directly regulate the internal forces between the bio-muscle models.

Details

Robotic Intelligence and Automation, vol. 43 no. 6
Type: Research Article
ISSN: 2754-6969

Keywords

Article
Publication date: 25 August 2023

Shuai Yue, Ben Niu, Huanqing Wang, Liang Zhang and Adil M. Ahmad

This paper aims to study the issues of adaptive fuzzy control for a category of switched under-actuated systems with input nonlinearities and external disturbances.

Abstract

Purpose

This paper aims to study the issues of adaptive fuzzy control for a category of switched under-actuated systems with input nonlinearities and external disturbances.

Design/methodology/approach

A control scheme based on sliding mode surface with a hierarchical structure is introduced to enhance the responsiveness and robustness of the studied systems. An equivalent control and switching control rules are co-designed in a hierarchical sliding mode control (HSMC) framework to ensure that the system state reaches a given sliding surface and remains sliding on the surface, finally stabilizing at the equilibrium point. Besides, the input nonlinearities consist of non-symmetric saturation and dead-zone, which are estimated by an unknown bounded function and a known affine function.

Findings

Based on fuzzy logic systems and the hierarchical sliding mode control method, an adaptive fuzzy control method for uncertain switched under-actuated systems is put forward.

Originality/value

The “cause and effect” problems often existing in conventional backstepping designs can be prevented. Furthermore, the presented adaptive laws can eliminate the influence of external disturbances and approximation errors. Besides, in contrast to arbitrary switching strategies, the authors consider a switching rule with average dwell time, which resolves control problems that cannot be resolved with arbitrary switching signals and reduces conservatism.

Details

Robotic Intelligence and Automation, vol. 43 no. 5
Type: Research Article
ISSN: 2754-6969

Keywords

1 – 10 of 723