Search results

1 – 10 of 204
Article
Publication date: 2 May 2024

Gerasimos G. Rigatos

To provide high torques needed to move a robot’s links, electric actuators are followed by a transmission system with a high transmission rate. For instance, gear ratios of 100:1…

Abstract

Purpose

To provide high torques needed to move a robot’s links, electric actuators are followed by a transmission system with a high transmission rate. For instance, gear ratios of 100:1 are often used in the joints of a robotic manipulator. This results into an actuator with large mechanical impedance (also known as nonback-drivable actuator). This in turn generates high contact forces when collision of the robotic mechanism occur and can cause humans’ injury. Another disadvantage of electric actuators is that they can exhibit overheating when constant torques have to be provided. Comparing to electric actuators, pneumatic actuators have promising properties for robotic applications, due to their low weight, simple mechanical design, low cost and good power-to-weight ratio. Electropneumatically actuated robots usually have better friction properties. Moreover, because of low mechanical impedance, pneumatic robots can provide moderate interaction forces which is important for robotic surgery and rehabilitation tasks. Pneumatic actuators are also well suited for exoskeleton robots. Actuation in exoskeletons should have a fast and accurate response. While electric motors come against high mechanical impedance and the risk of causing injuries, pneumatic actuators exhibit forces and torques which stay within moderate variation ranges. Besides, unlike direct current electric motors, pneumatic actuators have an improved weight-to-power ratio and avoid overheating problems.

Design/methodology/approach

The aim of this paper is to analyze a nonlinear optimal control method for electropneumatically actuated robots. A two-link robotic exoskeleton with electropneumatic actuators is considered as a case study. The associated nonlinear and multivariable state-space model is formulated and its differential flatness properties are proven. The dynamic model of the electropneumatic robot is linearized at each sampling instance with the use of first-order Taylor series expansion and through the computation of the associated Jacobian matrices. Within each sampling period, the time-varying linearization point is defined by the present value of the robot’s state vector and by the last sampled value of the control inputs vector. An H-infinity controller is designed for the linearized model of the robot aiming at solving the related optimal control problem under model uncertainties and external perturbations. An algebraic Riccati equation is solved at each time-step of the control method to obtain the stabilizing feedback gains of the H-infinity controller. Through Lyapunov stability analysis, it is proven that the robot’s control scheme satisfies the H-infinity tracking performance conditions which indicate the robustness properties of the control method. Moreover, global asymptotic stability is proven for the control loop. The method achieves fast convergence of the robot’s state variables to the associated reference trajectories, and despite strong nonlinearities in the robot’s dynamics, it keeps moderate the variations of the control inputs.

Findings

In this paper, a novel solution has been proposed for the nonlinear optimal control problem of robotic exoskeletons with electropneumatic actuators. As a case study, the dynamic model of a two-link lower-limb robotic exoskeleton with electropneumatic actuators has been considered. The dynamic model of this robotic system undergoes first approximate linearization at each iteration of the control algorithm around a temporary operating point. Within each sampling period, this linearization point is defined by the present value of the robot’s state vector and by the last sampled value of the control inputs vector. The linearization process relies on first-order Taylor series expansion and on the computation of the associated Jacobian matrices. The modeling error which is due to the truncation of higher-order terms from the Taylor series is considered to be a perturbation which is asymptotically compensated by the robustness of the control algorithm. To stabilize the dynamics of the electropneumatically actuated robot and to achieve precise tracking of reference setpoints, an H-infinity (optimal) feedback controller is designed. Actually, the proposed H-infinity controller for the model of the two-link electropneumatically actuated exoskeleton achieves the solution of the associated optimal control problem under model uncertainty and external disturbances. This controller implements a min-max differential game taking place between: (i) the control inputs which try to minimize a cost function which comprises a quadratic term of the state vector’s tracking error and (ii) the model uncertainty and perturbation inputs which try to maximize this cost function. To select the stabilizing feedback gains of this H-infinity controller, an algebraic Riccati equation is being repetitively solved at each time-step of the control method. The global stability properties of the H-infinity control scheme are proven through Lyapunov analysis.

Research limitations/implications

Pneumatic actuators are characterized by high nonlinearities which are due to air compressibility, thermodynamics and valves behavior and thus pneumatic robots require elaborated nonlinear control schemes to ensure their fast and precise positioning. Among the control methods which have been applied to pneumatic robots, one can distinguish differential geometric approaches (Lie algebra-based control, differential flatness theory-based control, nonlinear model predictive control [NMPC], sliding-mode control, backstepping control and multiple models-based fuzzy control). Treating nonlinearities and fault tolerance issues in the control problem of robotic manipulators with electropneumatic actuators has been a nontrivial task.

Practical implications

The novelty of the proposed control method is outlined as follows: preceding results on the use of H-infinity control to nonlinear dynamical systems were limited to the case of affine-in-the-input systems with drift-only dynamics. These results considered that the control inputs gain matrix is not dependent on the values of the system’s state vector. Moreover, in these approaches the linearization was performed around points of the desirable trajectory, whereas in the present paper’s control method the linearization points are related with the value of the state vector at each sampling instance as well as with the last sampled value of the control inputs vector. The Riccati equation which has been proposed for computing the feedback gains of the controller is novel, so is the presented global stability proof through Lyapunov analysis. This paper’s scientific contribution is summarized as follows: (i) the presented nonlinear optimal control method has improved or equally satisfactory performance when compared against other nonlinear control schemes that one can consider for the dynamic model of robots with electropneumatic actuators (such as Lie algebra-based control, differential flatness theory-based control, nonlinear model-based predictive control, sliding-mode control and backstepping control), (ii) it achieves fast and accurate tracking of all reference setpoints, (iii) despite strong nonlinearities in the dynamic model of the robot, it keeps moderate the variations of the control inputs and (iv) unlike the aforementioned alternative control approaches, this paper’s method is the only one that achieves solution of the optimal control problem for electropneumatic robots.

Social implications

The use of electropneumatic actuation in robots exhibits certain advantages. These can be the improved weight-to-power ratio, the lower mechanical impedance and the avoidance of overheating. At the same time, precise positioning and accurate execution of tasks by electropneumatic robots requires the application of elaborated nonlinear control methods. In this paper, a new nonlinear optimal control method has been developed for electropneumatically actuated robots and has been specifically applied to the dynamic model of a two-link robotic exoskeleton. The benefit from using this paper’s results in industrial and biomedical applications is apparent.

Originality/value

A comparison of the proposed nonlinear optimal (H-infinity) control method against other linear and nonlinear control schemes for electropneumatically actuated robots shows the following: (1) Unlike global linearization-based control approaches, such as Lie algebra-based control and differential flatness theory-based control, the optimal control approach does not rely on complicated transformations (diffeomorphisms) of the system’s state variables. Besides, the computed control inputs are applied directly on the initial nonlinear model of the electropneumatic robot and not on its linearized equivalent. The inverse transformations which are met in global linearization-based control are avoided and consequently one does not come against the related singularity problems. (2) Unlike model predictive control (MPC) and NMPC, the proposed control method is of proven global stability. It is known that MPC is a linear control approach that if applied to the nonlinear dynamics of the electropneumatic robot, the stability of the control loop will be lost. Besides, in NMPC the convergence of its iterative search for an optimum depends on initialization and parameter values selection and consequently the global stability of this control method cannot be always assured. (3) Unlike sliding-mode control and backstepping control, the proposed optimal control method does not require the state-space description of the system to be found in a specific form. About sliding-mode control, it is known that when the controlled system is not found in the input-output linearized form the definition of the sliding surface can be an intuitive procedure. About backstepping control, it is known that it cannot be directly applied to a dynamical system if the related state-space model is not found in the triangular (backstepping integral) form. (4) Unlike PID control, the proposed nonlinear optimal control method is of proven global stability, the selection of the controller’s parameters does not rely on a heuristic tuning procedure, and the stability of the control loop is assured in the case of changes of operating points. (5) Unlike multiple local models-based control, the nonlinear optimal control method uses only one linearization point and needs the solution of only one Riccati equation so as to compute the stabilizing feedback gains of the controller. Consequently, in terms of computation load the proposed control method for the electropneumatic actuator’s dynamics is much more efficient.

Details

Robotic Intelligence and Automation, vol. ahead-of-print no. ahead-of-print
Type: Research Article
ISSN: 2754-6969

Keywords

Article
Publication date: 7 May 2024

Andong Liu, Yawen Zhang, Jiayun Fu, Yuankun Yan and Wen-An Zhang

In response to the issue of traditional algorithms often falling into local minima or failing to find feasible solutions in manipulator path planning. The purpose of this paper is…

Abstract

Purpose

In response to the issue of traditional algorithms often falling into local minima or failing to find feasible solutions in manipulator path planning. The purpose of this paper is to propose a 3D artificial moment method (3D-AMM) for obstacle avoidance for the robotic arm's end-effector.

Design/methodology/approach

A new method for constructing temporary attractive points in 3D has been introduced using the vector triple product approach, which generates the attractive moments that attract the end-effector to move toward it. Second, distance weight factorization and spatial projection methods are introduced to improve the solution of repulsive moments in multiobstacle scenarios. Third, a novel motion vector-solving mechanism is proposed to provide nonzero velocity for the end-effector to solve the problem of limiting the solution of the motion vector to a fixed coordinate plane due to dimensionality constraints.

Findings

A comparative analysis was conducted between the proposed algorithm and the existing methods, the improved artificial potential field method and the rapidly-random tree method under identical simulation conditions. The results indicate that the 3D-AMM method successfully plans paths with smoother trajectories and reduces the path length by 20.03% to 36.9%. Additionally, the experimental comparison outcomes affirm the feasibility and effectiveness of this method for obstacle avoidance in industrial scenarios.

Originality/value

This paper proposes a 3D-AMM algorithm for manipulator path planning in Cartesian space with multiple obstacles. This method effectively solves the problem of the artificial potential field method easily falling into local minimum points and the low path planning success rate of the rapidly-exploring random tree method.

Details

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

Keywords

Article
Publication date: 6 May 2024

Ahmed Taibi, Said Touati, Lyes Aomar and Nabil Ikhlef

Bearings play a critical role in the reliable operation of induction machines, and their failure can lead to significant operational challenges and downtime. Detecting and…

Abstract

Purpose

Bearings play a critical role in the reliable operation of induction machines, and their failure can lead to significant operational challenges and downtime. Detecting and diagnosing these defects is imperative to ensure the longevity of induction machines and preventing costly downtime. The purpose of this paper is to develop a novel approach for diagnosis of bearing faults in induction machine.

Design/methodology/approach

To identify the different fault states of the bearing with accurately and efficiently in this paper, the original bearing vibration signal is first decomposed into several intrinsic mode functions (IMFs) using variational mode decomposition (VMD). The IMFs that contain more noise information are selected using the Pearson correlation coefficient. Subsequently, discrete wavelet transform (DWT) is used to filter the noisy IMFs. Second, the composite multiscale weighted permutation entropy (CMWPE) of each component is calculated to form the features vector. Finally, the features vector is reduced using the locality-sensitive discriminant analysis algorithm, to be fed into the support vector machine model for training and classification.

Findings

The obtained results showed the ability of the VMD_DWT algorithm to reduce the noise of raw vibration signals. It also demonstrated that the proposed method can effectively extract different fault features from vibration signals.

Originality/value

This study suggested a new VMD_DWT method to reduce the noise of the bearing vibration signal. The proposed approach for bearing fault diagnosis of induction machine based on VMD-DWT and CMWPE is highly effective. Its effectiveness has been verified using experimental data.

Details

COMPEL - The international journal for computation and mathematics in electrical and electronic engineering , vol. ahead-of-print no. ahead-of-print
Type: Research Article
ISSN: 0332-1649

Keywords

Article
Publication date: 8 May 2024

Hongze Wang

Many practical control problems require achieving multiple objectives, and these objectives often conflict with each other. The existing multi-objective evolutionary reinforcement…

Abstract

Purpose

Many practical control problems require achieving multiple objectives, and these objectives often conflict with each other. The existing multi-objective evolutionary reinforcement learning algorithms cannot achieve good search results when solving such problems. It is necessary to design a new multi-objective evolutionary reinforcement learning algorithm with a stronger searchability.

Design/methodology/approach

The multi-objective reinforcement learning algorithm proposed in this paper is based on the evolutionary computation framework. In each generation, this study uses the long-short-term selection method to select parent policies. The long-term selection is based on the improvement of policy along the predefined optimization direction in the previous generation. The short-term selection uses a prediction model to predict the optimization direction that may have the greatest improvement on overall population performance. In the evolutionary stage, the penalty-based nonlinear scalarization method is used to scalarize the multi-dimensional advantage functions, and the nonlinear multi-objective policy gradient is designed to optimize the parent policies along the predefined directions.

Findings

The penalty-based nonlinear scalarization method can force policies to improve along the predefined optimization directions. The long-short-term optimization method can alleviate the exploration-exploitation problem, enabling the algorithm to explore unknown regions while ensuring that potential policies are fully optimized. The combination of these designs can effectively improve the performance of the final population.

Originality/value

A multi-objective evolutionary reinforcement learning algorithm with stronger searchability has been proposed. This algorithm can find a Pareto policy set with better convergence, diversity and density.

Details

Robotic Intelligence and Automation, vol. ahead-of-print no. ahead-of-print
Type: Research Article
ISSN: 2754-6969

Keywords

Article
Publication date: 18 April 2024

Stefano Costa, Eugenio Costamagna and Paolo Di Barba

A novel method for modelling permanent magnets is investigated based on numerical approximations with rational functions. This study aims to introduce the AAA algorithm and other…

Abstract

Purpose

A novel method for modelling permanent magnets is investigated based on numerical approximations with rational functions. This study aims to introduce the AAA algorithm and other recently developed, cutting-edge mathematical tools, which provide outstandingly fast and accurate numerical computation of potentials and vector fields.

Design/methodology/approach

First, the AAA algorithm is briefly introduced along with its main variants and other advanced mathematical tools involved in the modelling. Then, the analysis of a circular Halbach array with a one-pole pair is carried out by means of the AAA-least squares method, focusing on vector potential and flux density in the bore and validating results by means of classic finite element software. Finally, the investigation is completed by a finite difference analysis.

Findings

AAA methods for field analysis prove to be strikingly fast and accurate. Results are in excellent agreement with those provided by the finite element model, and the very good agreement with those from finite differences suggests future improvements. They are also easy programming; the MATLAB code is less than 200 lines. This indicates they can provide an effective tool for rapid analysis.

Research limitations/implications

AAA methods in magnetostatics are novel, but their extension to analogous physical problems seems straightforward. Being a meshless method, it is unlikely that local non-linearities can be considered. An aspect of particular interest, left for future research, is the capability of handling inhomogeneous domains, i.e. solving general interface problems.

Originality/value

The authors use cutting-edge mathematical tools for the modelling of complex physical objects in magnetostatics.

Details

COMPEL - The international journal for computation and mathematics in electrical and electronic engineering , vol. ahead-of-print no. ahead-of-print
Type: Research Article
ISSN: 0332-1649

Keywords

Article
Publication date: 2 May 2024

Xin Fan, Yongshou Liu, Zongyi Gu and Qin Yao

Ensuring the safety of structures is important. However, when a structure possesses both an implicit performance function and an extremely small failure probability, traditional…

Abstract

Purpose

Ensuring the safety of structures is important. However, when a structure possesses both an implicit performance function and an extremely small failure probability, traditional methods struggle to conduct a reliability analysis. Therefore, this paper proposes a reliability analysis method aimed at enhancing the efficiency of rare event analysis, using the widely recognized Relevant Vector Machine (RVM).

Design/methodology/approach

Drawing from the principles of importance sampling (IS), this paper employs Harris Hawks Optimization (HHO) to ascertain the optimal design point. This approach not only guarantees precision but also facilitates the RVM in approximating the limit state surface. When the U learning function, designed for Kriging, is applied to RVM, it results in sample clustering in the design of experiment (DoE). Therefore, this paper proposes a FU learning function, which is more suitable for RVM.

Findings

Three numerical examples and two engineering problem demonstrate the effectiveness of the proposed method.

Originality/value

By employing the HHO algorithm, this paper innovatively applies RVM in IS reliability analysis, proposing a novel method termed RVM-HIS. The RVM-HIS demonstrates exceptional computational efficiency, making it eminently suitable for rare events reliability analysis with implicit performance function. Moreover, the computational efficiency of RVM-HIS has been significantly enhanced through the improvement of the U learning function.

Details

Engineering Computations, vol. ahead-of-print no. ahead-of-print
Type: Research Article
ISSN: 0264-4401

Keywords

Article
Publication date: 19 April 2024

Oguzhan Ozcelebi, Jose Perez-Montiel and Carles Manera

Might the impact of the financial stress on exchange markets be asymmetric and exposed to regime changes? Departing from the existing literature, highlighting that the domestic…

15

Abstract

Purpose

Might the impact of the financial stress on exchange markets be asymmetric and exposed to regime changes? Departing from the existing literature, highlighting that the domestic and foreign financial stress in terms of money market have substantial effects on exchange market, this paper aims to investigate the impacts of the bond yield spreads of three emerging countries (Mexico, Russia, and South Korea) on their exchange market pressure indices using monthly observations for the period 2010:01–2019:12. Additionally, the paper analyses the impact of bond yield spread of the US on the exchange market pressure indices of the three mentioned emerging countries. The authors hypothesized whether the negative and positive changes in the bond yield spreads have varying effects on exchange market pressure indices.

Design/methodology/approach

To address the research question, we measure the bond yield spread of the selected countries by using the interest rate spread between 10-year and 3-month treasury bills. At the same time, the exchange market pressure index is proxied by the index introduced by Desai et al. (2017). We base the empirical analysis on nonlinear vector autoregression (VAR) models and an asymmetric quantile-based approach.

Findings

The results of the impulse response functions indicate that increases/decreases in the bond yield spreads of Mexico, Russia and South Korea raise/lower their exchange market pressure, and the effects of shocks in the bond yield spreads of the US also lead to depreciation/appreciation pressures in the local currencies of the emerging countries. The quantile connectedness analysis, which allows for the role of regimes, reveals that the weights of the domestic and foreign bond yield spread in explaining variations of exchange market pressure indices are higher when exchange market pressure indices are not in a normal regime, indicating the role of extreme development conditions in the exchange market. The quantile regression model underlines that an increase in the domestic bond yield spread leads to a rise in its exchange market pressure index during all exchange market pressure periods in Mexico, and the relevant effects are valid during periods of high exchange market pressure in Russia. Our results also show that Russia differs from Mexico and South Korea in terms of the factors influencing the demand for domestic currency, and we have demonstrated the role of domestic macroeconomic and financial conditions in surpassing the effects of US financial stress. More specifically, the impacts of the domestic and foreign financial stress vary across regimes and are asymmetric.

Originality/value

This study enriches the literature on factors affecting the exchange market pressure of emerging countries. The results have significant economic implications for policymakers, indicating that the exchange market pressure index may trigger a financial crisis and economic recession.

Details

International Journal of Emerging Markets, vol. ahead-of-print no. ahead-of-print
Type: Research Article
ISSN: 1746-8809

Keywords

Article
Publication date: 7 May 2024

Xinzhe Li, Qinglong Li, Dasom Jeong and Jaekyeong Kim

Most previous studies predicting review helpfulness ignored the significance of deep features embedded in review text and instead relied on hand-crafted features. Hand-crafted and…

Abstract

Purpose

Most previous studies predicting review helpfulness ignored the significance of deep features embedded in review text and instead relied on hand-crafted features. Hand-crafted and deep features have the advantages of high interpretability and predictive accuracy. This study aims to propose a novel review helpfulness prediction model that uses deep learning (DL) techniques to consider the complementarity between hand-crafted and deep features.

Design/methodology/approach

First, an advanced convolutional neural network was applied to extract deep features from unstructured review text. Second, this study used previous studies to extract hand-crafted features that impact the helpfulness of reviews and enhance their interpretability. Third, this study incorporated deep and hand-crafted features into a review helpfulness prediction model and evaluated its performance using the Yelp.com data set. To measure the performance of the proposed model, this study used 2,417,796 restaurant reviews.

Findings

Extensive experiments confirmed that the proposed methodology performs better than traditional machine learning methods. Moreover, this study confirms through an empirical analysis that combining hand-crafted and deep features demonstrates better prediction performance.

Originality/value

To the best of the authors’ knowledge, this is one of the first studies to apply DL techniques and use structured and unstructured data to predict review helpfulness in the restaurant context. In addition, an advanced feature-fusion method was adopted to better use the extracted feature information and identify the complementarity between features.

研究目的

大多数先前预测评论有用性的研究忽视了嵌入在评论文本中的深层特征的重要性, 而主要依赖手工制作的特征。手工制作和深层特征具有高解释性和预测准确性的优势。本研究提出了一种新颖的评论有用性预测模型, 利用深度学习技术来考虑手工制作特征和深层特征之间的互补性。

研究方法

首先, 采用先进的卷积神经网络从非结构化的评论文本中提取深层特征。其次, 本研究利用先前研究中提取的手工制作特征, 这些特征影响了评论的有用性并增强了其解释性。第三, 本研究将深层特征和手工制作特征结合到一个评论有用性预测模型中, 并使用Yelp.com数据集对其性能进行评估。为了衡量所提出模型的性能, 本研究使用了2,417,796条餐厅评论。

研究发现

广泛的实验验证了所提出的方法优于传统的机器学习方法。此外, 通过实证分析, 本研究证实了结合手工制作和深层特征可以展现出更好的预测性能。

研究创新

据我们所知, 这是首个在餐厅评论预测中应用深度学习技术, 并结合了结构化和非结构化数据来预测评论有用性的研究之一。此外, 本研究采用了先进的特征融合方法, 更好地利用了提取的特征信息, 并识别了特征之间的互补性。

Article
Publication date: 25 April 2024

Muhammad Tariq, Muhammad Azam Khan and Niaz Ali

This study aims to investigate the effect of monetary policy on housing prices for US economy. It specifically examines whether nominal or real interest rates are the key drivers…

Abstract

Purpose

This study aims to investigate the effect of monetary policy on housing prices for US economy. It specifically examines whether nominal or real interest rates are the key drivers behind fluctuations in housing prices in US.

Design/methodology/approach

Monthly data from January 1991 to July 2023 and various appropriate analytical tools such as unit root tests, Johansen’s cointegration test, vector error correction model (VECM), impulse response function and Granger causality test were applied for the data analysis.

Findings

The Johansen cointegration findings reveal the presence of a long-term relationship among the variables. VECM results indicate a negative correlation between nominal and real interest rates and housing prices in both the short and long terms, suggesting that a strict monetary policy can help in controlling the housing price increase in the USA. However, housing prices are more responsive to changes in nominal interest rates than to real interest rates. Additionally, the study reveals that the COVID-19 pandemic contributed to the upsurge in housing prices in the USA.

Originality/value

This study contributes by examining the role that nominal or real interest rates play in shaping housing prices in the USA. Moreover, given the recent significant upsurge in housing prices, this study presents a unique opportunity to investigate whether these price increases are influenced by the Federal Reserve's monetary policy decisions regarding nominal or real interest rates. Additionally, using monthly data, this study provides a deeper understanding of the fluctuations in housing prices and their connection to monetary policy tools.

Details

International Journal of Housing Markets and Analysis, vol. ahead-of-print no. ahead-of-print
Type: Research Article
ISSN: 1753-8270

Keywords

Article
Publication date: 18 April 2024

Anton Salov

The purpose of this study is to reveal the dynamics of house prices and sales in spatial and temporal dimensions across British regions.

Abstract

Purpose

The purpose of this study is to reveal the dynamics of house prices and sales in spatial and temporal dimensions across British regions.

Design/methodology/approach

This paper incorporates two empirical approaches to describe the behaviour of property prices across British regions. The models are applied to two different data sets. The first empirical approach is to apply the price diffusion model proposed by Holly et al. (2011) to the UK house price index data set. The second empirical approach is to apply a bivariate global vector autoregression model without a time trend to house prices and transaction volumes retrieved from the nationwide building society.

Findings

Identifying shocks to London house prices in the GVAR model, based on the generalized impulse response functions framework, I find some heterogeneity in responses to house price changes; for example, South East England responds stronger than the remaining provincial regions. The main pattern detected in responses and characteristic for each region is the fairly rapid fading of the shock. The spatial-temporal diffusion model demonstrates the presence of a ripple effect: a shock emanating from London is dispersed contemporaneously and spatially to other regions, affecting prices in nondominant regions with a delay.

Originality/value

The main contribution of this work is the betterment in understanding how house price changes move across regions and time within a UK context.

Details

International Journal of Housing Markets and Analysis, vol. ahead-of-print no. ahead-of-print
Type: Research Article
ISSN: 1753-8270

Keywords

Access

Year

Last month (204)

Content type

Earlycite article (204)
1 – 10 of 204