Nothing Special   »   [go: up one dir, main page]

You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (220)

Search Parameters:
Keywords = rough terrain

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
23 pages, 3670 KiB  
Article
Vegetation Succession Patterns at Sperry Glacier’s Foreland, Glacier National Park, MT, USA
by Ami Bryant, Lynn M. Resler, Dianna Gielstra and Thomas Pingel
Land 2025, 14(2), 306; https://doi.org/10.3390/land14020306 - 2 Feb 2025
Viewed by 601
Abstract
Plant colonization patterns on deglaciated terrain give insight into the factors influencing alpine ecosystem development. Our objectives were to use a chronosequence, extending from the Little Ice Age (~1850) terminal moraine to the present glacier terminus, and biophysical predictors to characterize vegetation across [...] Read more.
Plant colonization patterns on deglaciated terrain give insight into the factors influencing alpine ecosystem development. Our objectives were to use a chronosequence, extending from the Little Ice Age (~1850) terminal moraine to the present glacier terminus, and biophysical predictors to characterize vegetation across Sperry Glacier’s foreland—a mid-latitude cirque glacier in Glacier National Park, Montana, USA. We measured diversity metrics (i.e., richness, evenness, and Shannon’s diversity index), percent cover, and community composition in 61 plots. Field observations characterized drainage, concavity, landform features, rock fragments, and geomorphic process domains in each plot. GIS-derived variables contextualized the plots’ aspect, terrain roughness, topographic position, solar radiation, and curvature. Overall, vegetation cover and species richness increased with terrain age, but with colonization gaps compared to other forelands, likely due to extensive bedrock and slow soil development, potentially putting this community at risk of being outpaced by climate change. Generalized linear models revealed the importance of local site factors (e.g., drainage, concavity, and process domain) in explaining species richness and Shannon’s diversity patterns. The relevance of field-measured variables over GIS-derived variables demonstrated the importance of fieldwork in understanding alpine successional patterns and the need for higher-resolution remote sensing analyses to expand these landscape-scale studies. Full article
Show Figures

Figure 1

Figure 1
<p>Sperry Glacier’s foreland facing north toward the sparsely vegetated Little Ice Age terminal moraines.</p>
Full article ">Figure 2
<p>Examples of vegetation establishment in a variety of microsites, including ledges (<b>a</b>), fissures (<b>b</b>), and rocks and boulders (<b>c</b>,<b>d</b>). Photos: A. Bryant 2022 and 2023.</p>
Full article ">Figure 3
<p>Map of the study area, including Sperry Glacier and Sperry Glacier’s foreland with glacier margin dates and 61 sample plot locations spanning from the 1850 LIA moraine to the present glacier margin.</p>
Full article ">Figure 4
<p>A poorly drained, concave plot (<b>a</b>); poorly drained, straight plot (<b>b</b>); and moderately drained, convex plot (<b>c</b>). Photos: A. Bryant 2022.</p>
Full article ">Figure 5
<p>Box plots depicting the distribution of vegetation cover (VC) (%) and species richness (SR) within terrain age ranges.</p>
Full article ">Figure 6
<p>Box plots showing the differences in species richness (SR) for each variable included in the best GLM model (for terrain age, see <a href="#land-14-00306-f005" class="html-fig">Figure 5</a>). The letters depict post hoc test groups. Shared letters indicate ranges that are not significantly different, while different letters indicate a significant difference (<span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">Figure 7
<p>Box plots showing the differences in Shannon’s diversity (SD) for each variable included in the best GLM model. The letters depict post hoc test groups. Shared letters indicate ranges that are not significantly different, while different letters indicate a significant difference (<span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">Figure 8
<p>A simplified interpretation of species richness (SR) trends over years since deglaciation from reported data in published glacial studies, including Sperry Glacier (this study), Jamtalferer Glacier in the Australian Eastern Alps [<a href="#B4-land-14-00306" class="html-bibr">4</a>], Skaftafellsjӧkull Glacier in Southern Iceland sampled in 2007 and 2014 [<a href="#B20-land-14-00306" class="html-bibr">20</a>], and Humboldt Glacier in Venezuela [<a href="#B80-land-14-00306" class="html-bibr">80</a>].</p>
Full article ">
19 pages, 9749 KiB  
Article
Numerical Simulation of Debris Flow Behavior over a Series of Groundsills
by Chyan-Deng Jan, Yi-Chao Zeng and Litan Dey
Water 2025, 17(3), 293; https://doi.org/10.3390/w17030293 - 21 Jan 2025
Viewed by 489
Abstract
Debris flows propagating in natural environments often encounter irregular terrain features, such as bottom roughness and man-made structures like groundsills, which significantly influence their behavior and dynamics. In practice, groundsills are commonly used as debris flow mitigation structures. This study examines the effects [...] Read more.
Debris flows propagating in natural environments often encounter irregular terrain features, such as bottom roughness and man-made structures like groundsills, which significantly influence their behavior and dynamics. In practice, groundsills are commonly used as debris flow mitigation structures. This study examines the effects of a beam-type groundsill array on the flow behavior of sediment mixtures in an inclined channel using numerical simulations. The sediment mixtures, modeled as Bingham fluids, were tested as they flowed over groundsill arrays with varying densities, characterized by the spacing-to-height ratio (d/h) ranging from 2 to 10. The findings indicate that interaction with the groundsills produces a hydraulic jump-like flow, reaching a height approximately 2.2 times the approach flow depth across different array densities. High-density arrays (d/h4) substantially hindered flow propagation, reducing front velocities but leading to sediment buildup upstream of the groundsills. Conversely, low-density arrays (d/h>4) facilitated smoother flow with higher velocities. These insights into the relationship between array density, flow behavior, and sediment trapping provide valuable guidance for optimizing groundsill array designs to effectively reduce the mobility of gravity-driven flows of non-Newtonian fluids (such as snow avalanches, debris, lava, or mudflows) and mitigate the associated risks. Full article
(This article belongs to the Section Hydraulics and Hydrodynamics)
Show Figures

Figure 1

Figure 1
<p>A series of groundsills constructed along the debris flow gully in Hualien County, Taiwan.</p>
Full article ">Figure 2
<p>Particle size distribution curve of the clay–silt materials used as sediment in the rheometer tests.</p>
Full article ">Figure 3
<p>(<b>a</b>) Schematic diagram illustrating the rheological measurement setup using a rotational rheometer [<a href="#B21-water-17-00293" class="html-bibr">21</a>]. (<b>b</b>) Rheological characteristics of the sediment mixture with <math display="inline"><semantics> <mrow> <mi>ϕ</mi> <mo>=</mo> <mn>0.30</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>Schematic diagrams of the experimental setup for simulating sediment mixture flow over a beam-type array of groundsills in an inclined channel: (<b>a</b>) side view and (<b>b</b>) top view at a fixed channel slope, <math display="inline"><semantics> <mrow> <mi>θ</mi> <mo>=</mo> <mn>15</mn> <mo>°</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 5
<p>Illustrations of beam-type groundsill arrays with (<b>a</b>) high array density (closely spaced groundsills), (<b>b</b>) low array density (widely spaced groundsills) in 2D and 3D, and (<b>c</b>) schematic diagram for the relation between <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> </mrow> </semantics></math> and channel slope, <math display="inline"><semantics> <mrow> <mi>θ</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 6
<p>Configuration of the numerical model domain for simulating flow over an array of groundsills at a fixed slope: (<b>a</b>) boundary conditions (S: symmetry and W: wall conditions) and (<b>b</b>) computational mesh.</p>
Full article ">Figure 7
<p>Visualization of the instantaneous velocity magnitude (m/s) of sediment mixtures propagating over a low array density (specifically, <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math>) of groundsills at different simulation times.</p>
Full article ">Figure 8
<p>Temporal evolution of the front position of sediment mixtures for varying groundsill array densities.</p>
Full article ">Figure 9
<p>Contour plots for the fraction of sediment mixture at simulation times of 0.7, 1.0, and 1.5 s for flow over groundsill arrays with densities of (<b>a</b>) <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math>, (<b>b</b>) <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>4</mn> </mrow> </semantics></math>, (<b>c</b>) <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>6</mn> </mrow> </semantics></math>, (<b>d</b>) <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>8</mn> </mrow> </semantics></math>, and (<b>e</b>) <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 9 Cont.
<p>Contour plots for the fraction of sediment mixture at simulation times of 0.7, 1.0, and 1.5 s for flow over groundsill arrays with densities of (<b>a</b>) <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math>, (<b>b</b>) <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>4</mn> </mrow> </semantics></math>, (<b>c</b>) <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>6</mn> </mrow> </semantics></math>, (<b>d</b>) <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>8</mn> </mrow> </semantics></math>, and (<b>e</b>) <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 10
<p>Contour plots for fraction of sediment mixture at simulation time <math display="inline"><semantics> <mrow> <mi>t</mi> <mo>=</mo> <mn>30</mn> <mo> </mo> <mi mathvariant="normal">s</mi> </mrow> </semantics></math> for flow over groundsill arrays with densities of <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mn>6</mn> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <mn>10</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 11
<p>Streamwise profiles of (<b>a</b>) flow depth and (<b>b</b>) depth-averaged velocity of sediment mixtures propagating over groundsill arrays at <math display="inline"><semantics> <mrow> <mi>t</mi> <mo>=</mo> <mn>1.5</mn> <mo> </mo> <mi mathvariant="normal">s</mi> </mrow> </semantics></math>, comparing cases with and without groundsills for different array densities.</p>
Full article ">Figure 12
<p>Schematic representation of the formation of hydraulic jump-like flows at the first groundsill for array density <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>4</mn> </mrow> </semantics></math>.</p>
Full article ">
16 pages, 6767 KiB  
Article
Terrain Irregularity Sensing by Evaluating Feet Coordinate Standard Deviation
by Tomas Luneckas, Mindaugas Luneckas and Dainius Udris
Appl. Sci. 2025, 15(1), 411; https://doi.org/10.3390/app15010411 - 4 Jan 2025
Viewed by 473
Abstract
Locomotion over rough terrain is still a problem yet to be solved for legged robots. One of the problems arises from the inability to identify terrain roughness during locomotion, which could be crucial for decision-making and successful task completion. Our proposed terrain roughness [...] Read more.
Locomotion over rough terrain is still a problem yet to be solved for legged robots. One of the problems arises from the inability to identify terrain roughness during locomotion, which could be crucial for decision-making and successful task completion. Our proposed terrain roughness method is inspired by the observation that humans can sense their limb position in space without looking at them, which allows us to estimate obstacle heights. This method is based on robot feet coordinate standard deviation (further referred to as SD) parameter evaluation. SD values could be categorized to represent different terrain roughness, and such categories could be useful for selecting different gaits for different terrains. In this paper, we investigate the possibility of using already known feet coordinates to evaluate terrain roughness by calculating their standard deviation (SD). We present simulation results that show that the SD value only depends on terrain roughness and is not influenced by large terrain slopes. Experiments were conducted with real robots while walking over obstacles with different gaits to validate the method. This research mainly aims to test how robot gaits influence SD parameters for terrain roughness evaluation. The experimental results showed that the SD parameter calculated from the robot’s foot coordinates can be used to evaluate terrain roughness. The robot’s gaits have little to no influence on the SD parameter. Full article
(This article belongs to the Section Robotics and Automation)
Show Figures

Figure 1

Figure 1
<p>Hexapod robot 3D CAD model with shown body and feet coordinate frames.</p>
Full article ">Figure 2
<p>Robot stances and ground clearance.</p>
Full article ">Figure 3
<p>Robot feet positions on low-roughness terrain.</p>
Full article ">Figure 4
<p>Simulation results for low- and high-roughness terrain. The arrows show when the robot starts to climb up the slope.</p>
Full article ">Figure 5
<p>Hexapod robot model.</p>
Full article ">Figure 6
<p>Function diagram of robot‘s control system.</p>
Full article ">Figure 7
<p>Robot’s foot with limit switch for indicating when the ground is touched.</p>
Full article ">Figure 8
<p>Hexapod gait diagram shows leg transfer sequences for (<b>a</b>) tripod, (<b>b</b>) tetrapod, (<b>c</b>) wave, and (<b>d</b>) ripple. Here: LF—left front, LM—left middle, LH—left hind, RF—right front, RM—right middle, and RH—right hind legs.</p>
Full article ">Figure 9
<p>Feet trajectory for adaptive locomotion projection in <math display="inline"><semantics> <mrow> <mi>y</mi> <mi>z</mi> </mrow> </semantics></math> plane. <span class="html-italic">h</span> is step height; <span class="html-italic">l</span> is step length.</p>
Full article ">Figure 10
<p>Feet <span class="html-italic">z</span> coordinate average values during locomotion over a rectangular obstacle using different gaits with deviation compensation using robot body inclinations; uncertainty ±4 mm. (<b>a</b>) Left front, (<b>b</b>) right front, (<b>c</b>) left middle, (<b>d</b>) right middle, (<b>e</b>) left hind, and (<b>f</b>) right hind legs.</p>
Full article ">Figure 11
<p>Standard deviation variation during robot locomotion over a rectangular obstacle using different gaits with deviation compensation using robot body inclinations.</p>
Full article ">
21 pages, 19119 KiB  
Article
Caterpillar-Inspired Multi-Gait Generation Method for Series-Parallel Hybrid Segmented Robot
by Mingyuan Dou, Ning He, Jianhua Yang, Lile He, Jiaxuan Chen and Yaojiumin Zhang
Biomimetics 2024, 9(12), 754; https://doi.org/10.3390/biomimetics9120754 - 11 Dec 2024
Viewed by 788
Abstract
The body structures and motion stability of worm-like and snake-like robots have garnered significant research interest. Recently, innovative serial–parallel hybrid segmented robots have emerged as a fundamental platform for a wide range of motion modes. To address the hyper-redundancy characteristics of these hybrid [...] Read more.
The body structures and motion stability of worm-like and snake-like robots have garnered significant research interest. Recently, innovative serial–parallel hybrid segmented robots have emerged as a fundamental platform for a wide range of motion modes. To address the hyper-redundancy characteristics of these hybrid structures, we propose a novel caterpillar-inspired Stable Segment Update (SSU) gait generation approach, establishing a unified framework for multi-segment robot gait generation. Drawing inspiration from the locomotion of natural caterpillars, the segments are modeled as rigid bodies with six degrees of freedom (DOF). The SSU gait generation method is specifically designed to parameterize caterpillar-like gaits. An inverse kinematics solution is derived by analyzing the forward kinematics and identifying the minimum lifting segment, framing the problem as a single-segment end-effector tracking task. Three distinct parameter sets are introduced within the SSU method to account for the stability of robot motion. These parameters, represented as discrete hump waves, are intended to improve motion efficiency during locomotion. Furthermore, the trajectories for each swinging segment are determined through kinematic analysis. Experimental results validate the effectiveness of the proposed SSU multi-gait generation method, demonstrating the successful traversal of gaps and rough terrain. Full article
(This article belongs to the Section Locomotion and Bioinspired Robotics)
Show Figures

Figure 1

Figure 1
<p>Natural caterpillar locomotion pattern. (<b>a</b>) Natural caterpillar locomotion sequence (the red dashed line represents stable segment; the yellow dashed line represents swinging segment). (<b>b</b>) Schematic diagram of natural caterpillar segments.</p>
Full article ">Figure 2
<p>Nine-state of one segment motion trajectory.</p>
Full article ">Figure 3
<p>The hump formed on the natural caterpillar locomotion in a single segment 9-state trajectory. The illustration of the hump formed in the SSU method (red segment ((<b>3</b>)–(<b>6</b>) left) is the segment that is about to enter the swinging phase during the stance phase; red segment (right) is the segment that has ended the swinging phase during the stance phase. The yellow segment is the swinging segment in the swinging phase).</p>
Full article ">Figure 4
<p>Footfall-pattern diagram of nature caterpillar gait.</p>
Full article ">Figure 5
<p>Robot mechanism and variables. (<b>a</b>) 3-RSR. (<b>b</b>) 4-3-RSR.</p>
Full article ">Figure 6
<p>The kinematics analysis of 4-3-RSR robot SSU parameters <math display="inline"><semantics> <mrow> <msub> <mi>N</mi> <mrow> <mi>min</mi> </mrow> </msub> </mrow> </semantics></math>. In the 3-RSR parallel mechanism, (<b>a</b>) the relationship of the distal plate center in axis <math display="inline"><semantics> <mi>X</mi> </semantics></math> coordinate component and base angle <math display="inline"><semantics> <mrow> <msub> <mi>θ</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>θ</mi> <mn>2</mn> </msub> <mo>,</mo> <msub> <mi>θ</mi> <mn>3</mn> </msub> </mrow> </semantics></math>; (<b>b</b>) the relationship of the pitch angle of the distal plate and base angle <math display="inline"><semantics> <mrow> <msub> <mi>θ</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>θ</mi> <mn>2</mn> </msub> <mo>,</mo> <msub> <mi>θ</mi> <mn>3</mn> </msub> </mrow> </semantics></math>; (<b>c</b>) the relationship of the pitch angle of the distal plate and distal plate center in axis <math display="inline"><semantics> <mi>X</mi> </semantics></math> coordinate component. (<b>d</b>) The 2-3-RSR mechanism and variables.</p>
Full article ">Figure 7
<p>The robot posture when <math display="inline"><semantics> <mrow> <mi>φ</mi> <mo>=</mo> <mrow> <mn>1</mn> <mo>/</mo> <mn>2</mn> </mrow> </mrow> </semantics></math>.</p>
Full article ">Figure 8
<p>SSU gait generation flowchart.</p>
Full article ">Figure 9
<p>The SSU swinging segment trajectory. (<b>a</b>) The gaits sequence for <math display="inline"><semantics> <mrow> <mfenced> <mrow> <mi>m</mi> <mo>−</mo> <mn>1</mn> </mrow> </mfenced> <mi>th</mi> </mrow> </semantics></math> segment trajectory <math display="inline"><semantics> <mrow> <msub> <mover accent="true"> <mi>s</mi> <mo stretchy="false">^</mo> </mover> <mrow> <msub> <mi>c</mi> <mn>1</mn> </msub> </mrow> </msub> <mfenced> <mi>t</mi> </mfenced> </mrow> </semantics></math>. (<b>b</b>) The gaits sequence for <math display="inline"><semantics> <mrow> <mfenced> <mrow> <mi>m</mi> <mo>−</mo> <mn>1</mn> </mrow> </mfenced> <mi>th</mi> </mrow> </semantics></math> segment trajectory <math display="inline"><semantics> <mrow> <msub> <mover accent="true"> <mi>s</mi> <mo stretchy="false">^</mo> </mover> <mrow> <mi>c</mi> <mn>2</mn> </mrow> </msub> <mfenced> <mi>t</mi> </mfenced> </mrow> </semantics></math>. (<b>c</b>) The trajectory of <math display="inline"><semantics> <mrow> <mfenced> <mrow> <mi>m</mi> <mo>−</mo> <mn>2</mn> </mrow> </mfenced> <mi>th</mi> </mrow> </semantics></math> segment progressive. (<b>d</b>) The compensate trajectory <math display="inline"><semantics> <mrow> <msub> <mover accent="true"> <mi>s</mi> <mo stretchy="false">^</mo> </mover> <mrow> <mi>c</mi> <mn>1</mn> </mrow> </msub> <mfenced> <mi>t</mi> </mfenced> </mrow> </semantics></math> of <math display="inline"><semantics> <mrow> <mfenced> <mrow> <mi>m</mi> <mo>−</mo> <mn>1</mn> </mrow> </mfenced> <mi>th</mi> </mrow> </semantics></math> segment. (<b>e</b>) The compensate trajectory <math display="inline"><semantics> <mrow> <msub> <mover accent="true"> <mi>s</mi> <mo stretchy="false">^</mo> </mover> <mrow> <mi>c</mi> <mn>2</mn> </mrow> </msub> <mfenced> <mi>t</mi> </mfenced> </mrow> </semantics></math> of <math display="inline"><semantics> <mrow> <mfenced> <mrow> <mi>m</mi> <mo>−</mo> <mn>1</mn> </mrow> </mfenced> <mi>th</mi> </mrow> </semantics></math> segment progressive.</p>
Full article ">Figure 10
<p>Three gaits pattern of 4-3-RSR robot. (<b>a</b>) The 1-1-1-1-1 gait, (<b>b</b>) 1-1-2-1 gait, and (<b>c</b>) 1-2-2 gait. Footfall-pattern diagram of the (<b>d</b>) 1-1-1-1-1 gait, (<b>e</b>) 1-1-2-1 gait, and (<b>f</b>) 1-2-2 gait.</p>
Full article ">Figure 11
<p>The 4-3-RSR robot. (<b>a</b>) Three rotary joints replace the sphere joint. (<b>b</b>) The 4-3-RSR robot press plate (left) and main view (right).</p>
Full article ">Figure 12
<p>Joint trajectories of 4-3-RSR robot. (<b>a</b>) The 1-1-1-1-1 gait, (<b>b</b>) 1-1-2-1 gait, and (<b>c</b>) 1-2-2 gait, where (1) (2) (3) (4) illuminate the 1st, 2nd, 3rd, and 4th 3-RSR parallel mechanism joint trajectories.</p>
Full article ">Figure 13
<p>Three gaits experiment of the 4-3-RSR robot. (<b>a</b>) The 1-1-1-1-1 gait, (<b>b</b>) 1-1-2-1 gait, and (<b>c</b>) 1-2-2 gait. (The red dotted line represents the stable segment, and the yellow dotted line represents the swinging segment).</p>
Full article ">Figure 14
<p>Locomotion of the 4-3-RSR robot rectilinear gait.</p>
Full article ">Figure 15
<p>The 1-1-1-1-1-1 gait crossing gaps.</p>
Full article ">Figure 16
<p>The 1-1-1-1-1-1 gait on roughness terrain.</p>
Full article ">
15 pages, 10958 KiB  
Article
ARS: AI-Driven Recovery Controller for Quadruped Robot Using Single-Network Model
by Han Sol Kang, Hyun Yong Lee, Ji Man Park, Seong Won Nam, Yeong Woo Son, Bum Su Yi, Jae Young Oh, Jun Ha Song, Soo Yeon Choi, Bo Geun Kim, Hyun Seok Kim and Hyouk Ryeol Choi
Biomimetics 2024, 9(12), 749; https://doi.org/10.3390/biomimetics9120749 - 10 Dec 2024
Viewed by 881
Abstract
Legged robots, especially quadruped robots, are widely used in various environments due to their advantage in overcoming rough terrains. However, falling is inevitable. Therefore, the ability to overcome a falling state is an essential ability for legged robots. In this paper, we propose [...] Read more.
Legged robots, especially quadruped robots, are widely used in various environments due to their advantage in overcoming rough terrains. However, falling is inevitable. Therefore, the ability to overcome a falling state is an essential ability for legged robots. In this paper, we propose a method to fully recover a quadruped robot from a fall using a single-neural network model. The neural network model is trained in two steps in simulations using reinforcement learning, and then directly applied to AiDIN-VIII, a quadruped robot with 12 degrees of freedom. Experimental results using the proposed method show that the robot can successfully recover from a fall within 5 s in various postures, even when the robot is completely turned over. In addition, we can see that the robot successfully recovers from a fall caused by a disturbance. Full article
(This article belongs to the Special Issue Bio-Inspired and Biomimetic Intelligence in Robotics: 2nd Edition)
Show Figures

Figure 1

Figure 1
<p>The proposed method, ARS, consists of two phases. In phase 1, fall recovery is executed by learning a self-righting and standing-up policy and training a behavior selector that chooses between them based on the state. In phase 2, a single neural network is trained to recover from a fall using reinforcement learning, incorporating imitation loss to replicate the actions of the hierarchical behavior-based controller learned in phase 1.</p>
Full article ">Figure 2
<p>Comparison of standard deviation of action rate and joint acceleration during fall recovery in a simulation.</p>
Full article ">Figure 3
<p>Comparison of standard deviation of action rate and joint acceleration during fall recovery in real robot experiment.</p>
Full article ">Figure 4
<p>Fall recovery from different postures. (<b>a</b>,<b>b</b>) A snapshot of the experiment showcasing the robot initialized in random poses and recovering from a fall.</p>
Full article ">Figure 5
<p>Fall recovery from a completely inverted posture.</p>
Full article ">Figure 6
<p>Snapshot of AiDIN-VIII recovering after being kicked using a single neural network trained with the proposed framework.</p>
Full article ">Figure 7
<p>Experiment in a 0.1 m high stair environment that was not encountered during the training of ARS.</p>
Full article ">Figure 8
<p>This graph depicts the learning performance with and without imitation loss. The blue line is when trained using imitation loss, and the pink line is when trained with a single neural network without imitation loss. The line represents the moving average of the reward, while the lighter color in the background represents the raw data of the reward.</p>
Full article ">Figure 9
<p>The red solid line represents the desired torque calculated through the PD controller, while the black dashed line represents the actual torque calculated using the current measured from the robot’s motor driver. Graph (<b>a</b>) displays the results of applying the policy learned without considering a hip torque limit reward, whereas graph (<b>b</b>) illustrates the results of applying the policy learned while considering a hip torque limit reward. Snapshots (<b>c</b>,<b>d</b>) depict experimental snapshots using a real robot without and with the hip torque limit reward, respectively.</p>
Full article ">
33 pages, 17311 KiB  
Article
Development of a Virtual Telehandler Model Using a Bond Graph
by Beatriz Puras, Gustavo Raush, Javier Freire, Germán Filippini, Pedro Roquet, Manel Tirado, Oriol Casadesús and Esteve Codina
Machines 2024, 12(12), 878; https://doi.org/10.3390/machines12120878 - 4 Dec 2024
Viewed by 892
Abstract
Recent technological advancements and evolving regulatory frameworks are catalysing the integration of renewable energy sources in construction equipment, with the objective of significantly reducing greenhouse gas emissions. The electrification of non-road mobile machinery (NRMM), particularly self-propelled Rough-Terrain Variable Reach Trucks (RTVRT) equipped with [...] Read more.
Recent technological advancements and evolving regulatory frameworks are catalysing the integration of renewable energy sources in construction equipment, with the objective of significantly reducing greenhouse gas emissions. The electrification of non-road mobile machinery (NRMM), particularly self-propelled Rough-Terrain Variable Reach Trucks (RTVRT) equipped with telescopic booms, presents notable stability challenges. The transition from diesel to electric propulsion systems alters, among other factors, the centre of gravity and the inertial matrix, necessitating precise load capacity determinations through detailed load charts to ensure operational safety. This paper introduces a virtual model constructed through multiphysics modelling utilising the bond graph methodology, incorporating both scalar and vector bonds to facilitate detailed interconnections between mechanical and hydraulic domains. The model encompasses critical components, including the chassis, rear axle, telescopic boom, attachment fork, and wheels, each requiring a comprehensive three-dimensional treatment to accurately resolve spatial dynamics. An illustrative case study, supported by empirical data, demonstrates the model’s capabilities, particularly in calculating ground wheel reaction forces and analysing the hydraulic self-levelling behaviour of the attachment fork. Notably, discrepancies within a 10% range are deemed acceptable, reflecting the inherent variability of field operating conditions. Experimental analyses validate the BG-3D simulation model of the telehandler implemented in 20-SIM establishing it as an effective tool for estimating stability limits with satisfactory precision and for predicting dynamic behaviour across diverse operating conditions. Additionally, the paper discusses prospective enhancements to the model, such as the integration of the virtual vehicle model with a variable inclination platform in future research phases, aimed at evaluating both longitudinal and lateral stability in accordance with ISO 22915 standards, promoting operator safety. Full article
(This article belongs to the Section Vehicle Engineering)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Telescopic machine. Source: AUSA; (<b>b</b>) virtual model (20-SIM animation tool) of telescopic machine.</p>
Full article ">Figure 2
<p>Main causes of overturn accident. Source: Health &amp; Safety Executive, HSE (UK) [<a href="#B7-machines-12-00878" class="html-bibr">7</a>].</p>
Full article ">Figure 3
<p>Modelling dynamic behaviour of the telehandler. Blue line: modelisation; red line: experimentation.</p>
Full article ">Figure 4
<p>Icons representing 3D prismatic joint; 3D rotation joint; spherical joint; rigid body; 3D rotation (R); and 3D transformation between point A and B (T) and bond graph 3D dynamics (PJ).</p>
Full article ">Figure 5
<p>Telehandler model (3D bond graph scheme).</p>
Full article ">Figure 6
<p>Bond graph representation of platform submodel (3D bond graph scheme).</p>
Full article ">Figure 7
<p>Bond graph representation of rear axle mechanism model and Steering System Model (3D bond graph scheme).</p>
Full article ">Figure 8
<p>Bond graph representation of tyre/soil interaction model (3D bond graph scheme).</p>
Full article ">Figure 9
<p>Bond graph representation of Telescopic Arm System Model (3D bond graph scheme): 1—boom, 6—telescopic arm, 8—attachment unit (fork), 10—load, 4 and 5—lift cylinders, 9—extension cylinder, 6 and 7—tilt cylinders, 2 and 3—slave cylinders.</p>
Full article ">Figure 10
<p>(<b>a</b>) Real cylinder, (<b>b</b>) Bond graph representation of hydraulic cylinder submodel (3D bond graph scheme), (<b>c</b>) Prismatic Joint 3D Bond Graph, (<b>d</b>) Hydraulic cylinder 1D Bond Graph.</p>
Full article ">Figure 11
<p>Hydraulic circuit corresponding to the actuation of the telescopic arm and its attachment.</p>
Full article ">Figure 12
<p>Hydraulic actuator system of boom arm (1D-BG submodel scheme).</p>
Full article ">Figure 13
<p>Load-holding valve, also called overcentre valve (1D-BG submodel scheme).</p>
Full article ">Figure 14
<p>Hydraulic block of directional control valves (1D-BG submodel scheme).</p>
Full article ">Figure 15
<p>Directional control valve, DCV (1D-BG submodel scheme).</p>
Full article ">Figure 16
<p>Experimental values of the tyre stiffness.</p>
Full article ">Figure 17
<p>Numerical simulation results: vertical motion of tyre’s centre of mass. (<b>a</b>) For different values of the tyre stiffness; (<b>b</b>) for damping coefficient = 1 kN s/m.</p>
Full article ">Figure 18
<p>(<b>a</b>) Left frontal view; (<b>b</b>) right frontal view of instrumented T164 prototype.</p>
Full article ">Figure 19
<p>Experimental ground reaction forces during test of lift movement (upward and downward) with the loads on the fork at 0 kg and 1600 kg, when the telescopic arm is fully retracted.</p>
Full article ">Figure 20
<p>Experimental ground reaction forces as a function of mass on the fork attachment. Solid line: maximum values; dashed line: minimum values).</p>
Full article ">Figure 21
<p>Experimental evolution of pressures in the chambers of the lift cylinder during the upward and downward movement of the lift arm for four load conditions.</p>
Full article ">Figure 22
<p>Experimental operating values of the overcentre valve during the upward and downward movement of the lift arm for four load conditions.</p>
Full article ">Figure 23
<p>Numerical ground reaction forces on the wheels (N) and time (s) due to lifting and lowering 1020 kg load, Pos. E.</p>
Full article ">Figure 24
<p>Experimental ground reaction forces and numerical for test: lifting and lowering a 640 kg mass: (<b>a</b>) extension in Pos. A and boom up; (<b>b</b>) Pos. A and boom down; (<b>c</b>) Pos. E and boom up; (<b>d</b>) Pos. E and boom down.</p>
Full article ">Figure 25
<p>% of difference between numerical and experimental results: (<b>a</b>) extended arm Pos. A; (<b>b</b>) Pos. E.</p>
Full article ">Figure 26
<p>Ground reaction forces and extension telescopic position for 640 kg mass on fork. Solid line: experimental values; dashed line: numerical values.</p>
Full article ">Figure 27
<p>Numerical values of fork self-levelling.</p>
Full article ">
32 pages, 11087 KiB  
Article
Path Planning and Motion Control of Robot Dog Through Rough Terrain Based on Vision Navigation
by Tianxiang Chen, Yipeng Huangfu, Sutthiphong Srigrarom and Boo Cheong Khoo
Sensors 2024, 24(22), 7306; https://doi.org/10.3390/s24227306 - 15 Nov 2024
Viewed by 2078
Abstract
This article delineates the enhancement of an autonomous navigation and obstacle avoidance system for a quadruped robot dog. Part one of this paper presents the integration of a sophisticated multi-level dynamic control framework, utilizing Model Predictive Control (MPC) and Whole-Body Control (WBC) from [...] Read more.
This article delineates the enhancement of an autonomous navigation and obstacle avoidance system for a quadruped robot dog. Part one of this paper presents the integration of a sophisticated multi-level dynamic control framework, utilizing Model Predictive Control (MPC) and Whole-Body Control (WBC) from MIT Cheetah. The system employs an Intel RealSense D435i depth camera for depth vision-based navigation, which enables high-fidelity 3D environmental mapping and real-time path planning. A significant innovation is the customization of the EGO-Planner to optimize trajectory planning in dynamically changing terrains, coupled with the implementation of a multi-body dynamics model that significantly improves the robot’s stability and maneuverability across various surfaces. The experimental results show that the RGB-D system exhibits superior velocity stability and trajectory accuracy to the SLAM system, with a 20% reduction in the cumulative velocity error and a 10% improvement in path tracking precision. The experimental results also show that the RGB-D system achieves smoother navigation, requiring 15% fewer iterations for path planning, and a 30% faster success rate recovery in challenging environments. The successful application of these technologies in simulated urban disaster scenarios suggests promising future applications in emergency response and complex urban environments. Part two of this paper presents the development of a robust path planning algorithm for a robot dog on a rough terrain based on attached binocular vision navigation. We use a commercial-of-the-shelf (COTS) robot dog. An optical CCD binocular vision dynamic tracking system is used to provide environment information. Likewise, the pose and posture of the robot dog are obtained from the robot’s own sensors, and a kinematics model is established. Then, a binocular vision tracking method is developed to determine the optimal path, provide a proposal (commands to actuators) of the position and posture of the bionic robot, and achieve stable motion on tough terrains. The terrain is assumed to be a gentle uneven terrain to begin with and subsequently proceeds to a more rough surface. This work consists of four steps: (1) pose and position data are acquired from the robot dog’s own inertial sensors, (2) terrain and environment information is input from onboard cameras, (3) information is fused (integrated), and (4) path planning and motion control proposals are made. Ultimately, this work provides a robust framework for future developments in the vision-based navigation and control of quadruped robots, offering potential solutions for navigating complex and dynamic terrains. Full article
Show Figures

Figure 1

Figure 1
<p>Simplified box model of the Lite3P quadruped robotic dog.</p>
Full article ">Figure 2
<p>Internal sensor arrangement of the quadruped robotic dog.</p>
Full article ">Figure 3
<p>Dynamic control flowchart.</p>
Full article ">Figure 4
<p>MPC flowchart.</p>
Full article ">Figure 5
<p>WBC flowchart [<a href="#B30-sensors-24-07306" class="html-bibr">30</a>].</p>
Full article ">Figure 6
<p>Robot coordinates and joint point settings [<a href="#B30-sensors-24-07306" class="html-bibr">30</a>].</p>
Full article ">Figure 7
<p>Intel D435i and velodyne LIDAR.</p>
Full article ">Figure 8
<p>ICP diagram.</p>
Full article ">Figure 9
<p>Comparison of before and after modifying the perception region.</p>
Full article ">Figure 10
<p>Point cloud processing flowchart.</p>
Full article ">Figure 11
<p>{p, v} generation: (<b>a</b>) the creation of {p, v} pairs for collision points; (<b>b</b>) the process of generating anchor points and repulsive vectors for dynamic obstacle avoidance [<a href="#B41-sensors-24-07306" class="html-bibr">41</a>].</p>
Full article ">Figure 12
<p>Overall framework of 2D EGO-Planner.</p>
Full article ">Figure 13
<p>Robot initialization and control process in Gazebo simulation: (<b>a</b>) Gazebo environment creation, (<b>b</b>) robot model import, (<b>c</b>) torque balance mode activation, and (<b>d</b>) robot stepping and rotation in simulation.</p>
Full article ">Figure 14
<p>Joint rotational angles of FL and RL legs.</p>
Full article ">Figure 15
<p>Joint angular velocities of FL and RL legs.</p>
Full article ">Figure 16
<p>Torque applied to FL and RL joints during the gait cycle.</p>
Full article ">Figure 17
<p>The robot navigating in a simple environment using a camera.</p>
Full article ">Figure 18
<p>The robot navigating in a complex environment using a camera.</p>
Full article ">Figure 19
<p>A 2D trajectory showing start and goal positions, obstacles, and rough path.</p>
Full article ">Figure 20
<p>Initial environment setup.</p>
Full article ">Figure 21
<p>The robot starts navigating in a simple environment with a static obstacle (brown box).</p>
Full article ">Figure 22
<p>Dynamic Obstacle 1 introduced: the robot detects a new obstacle and recalculates its path.</p>
Full article ">Figure 23
<p>Dynamic Obstacle 2 introduced: after avoiding the first obstacle, a second obstacle is introduced and detected by the planner.</p>
Full article ">Figure 24
<p>Approaching the target: the robot adjusts its path to approach the target point as the distance shortens.</p>
Full article ">Figure 25
<p>Reaching the target: the robot completes its path and reaches the designated target point.</p>
Full article ">Figure 26
<p>Real-time B-spline trajectory updates in response to dynamic obstacles. Set 1 (orange) shows the initial path avoiding static obstacles. When the first dynamic obstacle is detected, the EGO-Planner updates the path (Set 2, blue) using local optimization. A second obstacle prompts another adjustment (Set 3, green), guiding the robot smoothly towards the target as trajectory updates become more frequent.</p>
Full article ">Figure 27
<p>The robot navigating a simple environment using SLAM.</p>
Full article ">Figure 28
<p>The robot navigating a complex environment using SLAM.</p>
Full article ">Figure 29
<p>A 2D trajectory showing start and goal positions, obstacles, and the planned path in a complex environment using SLAM.</p>
Full article ">Figure 30
<p>Navigation based on RGB-D camera.</p>
Full article ">Figure 31
<p>Navigation based on SLAM.</p>
Full article ">Figure 32
<p>Velocity deviation based on RGB-D camera.</p>
Full article ">Figure 33
<p>Velocity deviation based on SLAM.</p>
Full article ">Figure 34
<p>Cumulative average iterations.</p>
Full article ">Figure 35
<p>Cumulative success rate.</p>
Full article ">
13 pages, 3510 KiB  
Article
Laboratory Validation of 3D Model and Investigating Its Application to Wind Turbine Noise Propagation over Rough Ground
by James Naylor and Qin Qin
Wind 2024, 4(4), 363-375; https://doi.org/10.3390/wind4040018 - 7 Nov 2024
Viewed by 642
Abstract
In an investigation into how wind turbine noise interacts with the surrounding terrain, its propagation over rough ground is simulated using a parabolic equation code using a modified effective impedance model, which characterizes the effects of a three-dimensional, rigid roughness within a relatively [...] Read more.
In an investigation into how wind turbine noise interacts with the surrounding terrain, its propagation over rough ground is simulated using a parabolic equation code using a modified effective impedance model, which characterizes the effects of a three-dimensional, rigid roughness within a relatively long wavelength limit (ka1). The model is validated by comparison to experiments conducted within an anechoic chamber wherein different source–receiver geometries are considered. The relative sound pressure level spectra from the parabolic equation code using the modified effective impedance model highlight a sensitivity to the roughness parameters. At a low frequency and far distance, the relative sound pressure level decreased as the roughness coverage increased. A difference of 4.9 dB has been reported. The simulations highlight how the roughness shifts the ground effect dips, resulting in the sound level at the distance of 2 km being altered. However, only the monochromatic wave has been discussed. Further work on broadband noise is desirable. Furthermore, due to the long wavelength limit, only a portion of audible wind turbine noise can be investigated. Full article
Show Figures

Figure 1

Figure 1
<p>Roughness types. (<b>a</b>) Closely packed; (<b>b</b>) Medium packed; (<b>c</b>) Sparsely packed.</p>
Full article ">Figure 2
<p>Imaginary part of the normalized effective impedance calculated using the Attenborough–Tolstoy impedance model for various roughness types.</p>
Full article ">Figure 3
<p>Diagram of experimental setup.</p>
Full article ">Figure 4
<p>Comparison of the measured excess attenuation to the analytical and parabolic equation theory using the Attenborough–Tolstoy effective impedance model. (<b>a</b>) Closely packed roughness; (<b>b</b>) Medium packed roughness; (<b>c</b>) Sparsely packed roughness. For all figures (<b>a</b>–<b>c</b>), the sound source height, receiver height and horizontal distance between the two were 100 mm, 100 mm and 400 mm, respectively, for situation 1.</p>
Full article ">Figure 4 Cont.
<p>Comparison of the measured excess attenuation to the analytical and parabolic equation theory using the Attenborough–Tolstoy effective impedance model. (<b>a</b>) Closely packed roughness; (<b>b</b>) Medium packed roughness; (<b>c</b>) Sparsely packed roughness. For all figures (<b>a</b>–<b>c</b>), the sound source height, receiver height and horizontal distance between the two were 100 mm, 100 mm and 400 mm, respectively, for situation 1.</p>
Full article ">Figure 5
<p>Comparison of the measured relative sound pressure level, Rel. SPL, to the analytical and parabolic equation theory using the Attenborough–Tolstoy effective impedance model for a closely packed roughness. The sound source height, receiver height and horizontal distance between the two were 200 mm, 20 mm, and 100 mm, respectively, for situation 2.</p>
Full article ">Figure 6
<p>Simplified diagram of the scenario in mind for the parabolic equation method simulations.</p>
Full article ">Figure 7
<p>Relative sound pressure level, Rel. SPL, against propagation range, measured at 2 m height. Calculated for various roughness types using the parabolic equation method.</p>
Full article ">Figure 8
<p>Relative sound pressure level, Rel. SPL, map close to ground level. Calculated using the parabolic equation method. (<b>a</b>) Sparsely packed roughness; (<b>b</b>) Medium packed roughness; (<b>c</b>) Closely packed roughness.</p>
Full article ">Figure 8 Cont.
<p>Relative sound pressure level, Rel. SPL, map close to ground level. Calculated using the parabolic equation method. (<b>a</b>) Sparsely packed roughness; (<b>b</b>) Medium packed roughness; (<b>c</b>) Closely packed roughness.</p>
Full article ">
23 pages, 21957 KiB  
Article
Terrain Analysis According to Multiscale Surface Roughness in the Taklimakan Desert
by Sebastiano Trevisani and Peter L. Guth
Land 2024, 13(11), 1843; https://doi.org/10.3390/land13111843 - 5 Nov 2024
Viewed by 906
Abstract
Surface roughness, interpreted in the wide sense of surface texture, is a generic term referring to a variety of aspects and scales of spatial variability of surfaces. The analysis of solid earth surface roughness is useful for understanding, characterizing, and monitoring geomorphic factors [...] Read more.
Surface roughness, interpreted in the wide sense of surface texture, is a generic term referring to a variety of aspects and scales of spatial variability of surfaces. The analysis of solid earth surface roughness is useful for understanding, characterizing, and monitoring geomorphic factors at multiple spatiotemporal scales. The different geomorphic features characterizing a landscape exhibit specific characteristics and scales of surface texture. The capability to selectively analyze specific roughness metrics at multiple spatial scales represents a key tool in geomorphometric analysis. This research presents a simplified geostatistical approach for the multiscale analysis of surface roughness, or of image texture in the case of images, that is highly informative and interpretable. The implemented approach is able to describe two main aspects of short-range surface roughness: omnidirectional roughness and roughness anisotropy. Adopting simple upscaling approaches, it is possible to perform a multiscale analysis of roughness. An overview of the information extraction potential of the approach is shown for the analysis of a portion of the Taklimakan desert (China) using a 30 m resolution DEM derived from the Copernicus Glo-30 DSM. The multiscale roughness indexes are used as input features for unsupervised and supervised learning tasks. The approach can be refined both from the perspective of the multiscale analysis as well as in relation to the surface roughness indexes considered. However, even in its present, simplified form, it can find direct applications in relation to multiple contexts and research topics. Full article
(This article belongs to the Section Land, Soil and Water)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Reprojected COP DEM (30 m resolution, UTM F44) of the area of interest overlaid on the hillshade.</p>
Full article ">Figure 2
<p>Sentinel-2 true color RGB image (bands 4, 3, and 2) of the study area, with the main dune morphologies labeled.</p>
Full article ">Figure 3
<p>Main dune morphologies in the study area, visualized using Sentinel-2 imagery (<b>a</b>), hillshade (<b>b</b>), and residual DEM (<b>c</b>). From top to bottom: network/transverse dunes, longitudinal and transverse dunes, and dome-shaped dunes.</p>
Full article ">Figure 4
<p>Mixed morphologies in the area of interest, visualized using Sentinel-2 imagery (<b>a</b>), hillshade (<b>b</b>), and residual DEM (<b>c</b>). From top to bottom: outcropping bedrock with shadow and linear dunes, fluvial morphology, and a flat area with minor dune morphologies.</p>
Full article ">Figure 5
<p>RA direction, where the RA strength is higher than 0.3, overlaid on the hillshade (<b>a</b>) and the residual DEM (<b>b</b>) calculated for level L2.</p>
Full article ">Figure 6
<p>Omnidirectional short-range roughness (m) for the different resolutions. Different color scales for each diagram.</p>
Full article ">Figure 7
<p>Roughness anisotropy strength at different resolutions. Different color scales for each diagram.</p>
Full article ">Figure 8
<p>RGB image (each band normalized) of 3 omnidirectional roughness indexes computed at different resolutions (B = L1; G = L2; R = L4). Despite the high correlation of the three indexes, they differentiate very well the morphological features of the area. For example, they markedly highlight the characteristic smoothness of interdune areas of the longitudinal dunes south of the mountain ridge.</p>
Full article ">Figure 9
<p>RGB image (each band normalized) of 3 anisotropy strength roughness indexes computed at different resolutions (B = L1; G = L4; R = L16). In the dune fields north of the mountains, long-wavelength anisotropic features prevail; in contrast, for the southern longitudinal dunes, shorter anisotropic features (L4) are highlighted.</p>
Full article ">Figure 10
<p>Landscape clustered according to multiscale surface roughness indexes. The cluster centers in terms of OR and RA are described in <a href="#land-13-01843-f011" class="html-fig">Figure 11</a>.</p>
Full article ">Figure 11
<p>Cluster centers of the 7 classes resulting from K-means clustering for OR and RA at the different levels.</p>
Full article ">Figure 12
<p>MRI clustering results in the area of the northern dune field, characterized by network and transverse dunes. Clustering results (<b>d</b>), Sentinel-2 imagery (<b>a</b>), hillshade (<b>b</b>), and residual DEM (<b>c</b>).</p>
Full article ">Figure 13
<p>MRI clustering results in the area of the southern longitudinal dune fields. Clustering results (<b>d</b>), Sentinel-2 imagery (<b>a</b>), hillshade (<b>b</b>), and residual DEM (<b>c</b>).</p>
Full article ">Figure 14
<p>MRI clustering results in the area with fluvial morphology, outcropping bedrock, and dome dune fields. Clustering results (<b>d</b>), Sentinel-2 imagery (<b>a</b>), hillshade (<b>b</b>), and residual DEM (<b>c</b>).</p>
Full article ">Figure 15
<p>Manual classification of crest lines (<b>a</b>) for large dunes using visual analysis of slope (<b>b</b>), profile curvature (<b>c</b>), and residual DEM (<b>d</b>). Crest lines are associated with high positive profile curvature, strongly positive residual DEM, and low slope. These locations are then located in areas in which the neighborhood is characterized by an abrupt variation in the selected geomorphometric derivatives.</p>
Full article ">Figure 16
<p>Probability of observing a crest obtained by means of RF considering the GDs integrated with the MRIs (<b>a</b>) and only the five GDs (<b>b</b>) to obtain details of the study area, which is located on the western mountain ridge. The RF model integrating the MRIs provides a more focused prediction of crest lines of large dunes. In (<b>c</b>), the prediction of the crest lines of the two RF models is compared. Pixels with a probability higher than 0.8 have been classified as crests. The transparent color is where both models predicted a not-crest pixel, green is where both models predicted a crest, and red and blue are where, respectively, only RF GDs and RF GDs + MRIs predicted a crest.</p>
Full article ">Figure 17
<p>Variables’ importance in the two RF models according to the mean decrease in the Gini index ((<b>a</b>), RF based on GDs; (<b>b</b>), RF based on GDs integrated with MRIs).</p>
Full article ">Figure 18
<p>Prediction of crest lines with the RF model based on GDs and MRIs of an unseen area ((<b>c</b>), green box) external to the one with reference data used for training and testing ((<b>c</b>), red box). The reference crest lines (<b>a</b>) have been manually digitized by means of visual analysis of the profile curvature, the residual DEM, and the slope; the predicted crest lines have been derived as crests of all of the pixels with a probability above 0.8. The predicted crest lines are compared with the reference data (<b>b</b>). Green pixels are correctly classified as crests; red and blue pixels are incorrectly classified, respectively, as crests and not crests.</p>
Full article ">
22 pages, 6213 KiB  
Article
Simulation of the Neutral Atmospheric Flow Using Multiscale Modeling: Comparative Studies for SimpleFoam and Fluent Solver
by Zihan Zhao, Lingxiao Tang and Yiqing Xiao
Atmosphere 2024, 15(10), 1259; https://doi.org/10.3390/atmos15101259 - 21 Oct 2024
Viewed by 744
Abstract
The reproduced planetary boundary layer (PBL) wind is commonly applied in downscaled simulations using commercial CFD codes with Reynolds-averaged Navier–Stokes (RANS) turbulence modeling. When using the turbulent inlets calculated by numerical weather prediction models (NWP), adjustments of the turbulence eddy viscosity closures and [...] Read more.
The reproduced planetary boundary layer (PBL) wind is commonly applied in downscaled simulations using commercial CFD codes with Reynolds-averaged Navier–Stokes (RANS) turbulence modeling. When using the turbulent inlets calculated by numerical weather prediction models (NWP), adjustments of the turbulence eddy viscosity closures and wall function formulations are concerned with maintaining the fully developed wind profiles specified at the inlet of CFD domains. The impact of these related configurations is worth discussing in engineering applications, especially when commercial codes restrict the internal modifications. This study evaluates the numerical performances of open-source OpenFOAM 2.3.0 and commercial Fluent 17.2 codes as supplementary scientific comparisons. This contribution focuses on the modified turbulence closures to incorporate turbulent profiles produced from mesoscale PBL parameterizations and the modified wall treatments relating to the roughness length. The near-ground flow features are evaluated by selecting the flat terrains and the classical Askervein benchmark case. The improvement in near-ground wind flow under the downscaled framework shows satisfactory performance in the open-source CFD platform. This contributes to engineers realizing the micro-siting of wind turbines and quantifying terrain-induced speed-up phenomena under the scope of wind-resistant design. Full article
(This article belongs to the Section Atmospheric Techniques, Instruments, and Modeling)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Schematic of the full-scale CFD domain with flat terrain. Five positions of wind profiles (D0~D2500) are extracted along the streamwise direction, where the number represents the straight-line distance from the inlet. The purple gridlines represent the mesh grid. (<b>b</b>) Details of the surface-grid resolution in the streamwise direction.</p>
Full article ">Figure 2
<p>(<b>a</b>,<b>b</b>) are the structured grids on Askervein Hill case with central refinements close to measurements (Line A-A, AA-AA, RS and HT). Note that the positive X and Y axes represent the east and north directions, respectively. The green color gridlines represent the surface mesh grid. (<b>c</b>) Details of the XZ plane grid along with the hilltop HT. The red gridlines represent the spatial mesh grid, and blue background color is the working space during grid generation.</p>
Full article ">Figure 3
<p>The concerned parameters in the <math display="inline"><semantics> <mrow> <mi>k</mi> <mo>−</mo> <mi>ε</mi> </mrow> </semantics></math> turbulence solver using the multiscale modeling.</p>
Full article ">Figure 4
<p>Summary of the inlet conditions used in the flat terrain case. (<b>a</b>) Wind magnitude profile, (<b>b</b>) turbulent kinetic energy profile, (<b>c</b>) dissipation profile.</p>
Full article ">Figure 5
<p>Summary of the three-hour averaged wind profiles extracted from WRF solutions on the Askervein TU-03B field campaign. Twenty-eight positions are specified on the west (W-WRF) and south (S-WRF) of CFD lateral boundaries. (<b>a</b>–<b>c</b>) Wind components and directions. (<b>d</b>,<b>e</b>) Turbulent kinetic energy and dissipation profiles from WRF PBL parameterizations.</p>
Full article ">Figure 6
<p>Illustration of offline data transfer between mesoscale coarse profiles data (<b>left</b>) and CFD fine-grained boundaries (<b>right</b>). The different colors represent the interpolated variables.</p>
Full article ">Figure 7
<p>Comparison of inlet and outlet wind velocity profiles (<b>a</b>,<b>d</b>), turbulent kinetic profiles (<b>b</b>,<b>e</b>), and turbulent kinetic energy dissipation (<b>c</b>,<b>f</b>). The model constants and wall functions are configured with the standard treatments listed in <a href="#atmosphere-15-01259-t003" class="html-table">Table 3</a>. The different color lines represent the position extracted in <a href="#atmosphere-15-01259-f001" class="html-fig">Figure 1</a>. (FL-Fluent code, OF-OpenFOAM code).</p>
Full article ">Figure 7 Cont.
<p>Comparison of inlet and outlet wind velocity profiles (<b>a</b>,<b>d</b>), turbulent kinetic profiles (<b>b</b>,<b>e</b>), and turbulent kinetic energy dissipation (<b>c</b>,<b>f</b>). The model constants and wall functions are configured with the standard treatments listed in <a href="#atmosphere-15-01259-t003" class="html-table">Table 3</a>. The different color lines represent the position extracted in <a href="#atmosphere-15-01259-f001" class="html-fig">Figure 1</a>. (FL-Fluent code, OF-OpenFOAM code).</p>
Full article ">Figure 8
<p>Comparison of inlet and outlet wind velocity profiles (<b>a</b>,<b>d</b>), turbulent kinetic profiles (<b>b</b>,<b>e</b>), and turbulent kinetic energy dissipation (<b>c</b>,<b>f</b>). The numerical cases are configured with the standard model constants and modified wall functions as listed in <a href="#atmosphere-15-01259-t003" class="html-table">Table 3</a>. The different color lines represent the position extracted in <a href="#atmosphere-15-01259-f001" class="html-fig">Figure 1</a>.</p>
Full article ">Figure 9
<p>Comparison of inlet and outlet wind velocity profiles (<b>a</b>,<b>b</b>), turbulent kinetic profiles (<b>c</b>,<b>d</b>), and turbulent kinetic energy dissipation (<b>e</b>,<b>f</b>). The numerical cases are configured with the modified model constants and the standard wall functions, as listed in <a href="#atmosphere-15-01259-t003" class="html-table">Table 3</a>. The different color lines represent the position extracted in <a href="#atmosphere-15-01259-f001" class="html-fig">Figure 1</a>.</p>
Full article ">Figure 9 Cont.
<p>Comparison of inlet and outlet wind velocity profiles (<b>a</b>,<b>b</b>), turbulent kinetic profiles (<b>c</b>,<b>d</b>), and turbulent kinetic energy dissipation (<b>e</b>,<b>f</b>). The numerical cases are configured with the modified model constants and the standard wall functions, as listed in <a href="#atmosphere-15-01259-t003" class="html-table">Table 3</a>. The different color lines represent the position extracted in <a href="#atmosphere-15-01259-f001" class="html-fig">Figure 1</a>.</p>
Full article ">Figure 10
<p>Profiles of terrain-induced wind speed-up <math display="inline"><semantics> <mrow> <mo>Δ</mo> <mi>S</mi> </mrow> </semantics></math> and turbulent kinetic energy <math display="inline"><semantics> <mi>k</mi> </semantics></math>. Numerical cases are configured with the standard model constants (SC) and the wall function (SW). (<b>a</b>,<b>b</b>) <math display="inline"><semantics> <mrow> <mo>Δ</mo> <mi>S</mi> </mrow> </semantics></math> along line A-A and line AA-AA. (<b>c</b>) <math display="inline"><semantics> <mi>k</mi> </semantics></math> along lines A-A. The horizontally distributed profiles are 10 m AGL in height. (<b>d</b>) The vertical <math display="inline"><semantics> <mi>k</mi> </semantics></math> profile in the hilltop HT site. D-HT and D-CP represent the distance from HT and CP positions.</p>
Full article ">Figure 11
<p>Contours of the horizontal wind velocities (<b>a</b>,<b>b</b>) and turbulent kinetic energy (<b>c</b>,<b>d</b>) distributed with 10 m AGL height. The black and white circles in (<b>a</b>,<b>b</b>) denote the maximum and minimum wind speeds in the hilltop and leeward regions, respectively. The numerical cases are configured with the standard model constants (SC) and near-wall treatments (SW).</p>
Full article ">Figure 12
<p>The modified wall functions (MW) are configured in the numerical cases. (<b>a</b>,<b>b</b>) are <math display="inline"><semantics> <mrow> <mo>Δ</mo> <mi>S</mi> </mrow> </semantics></math> along line A-A and line AA-AA. (<b>c</b>) <math display="inline"><semantics> <mi>k</mi> </semantics></math> along lines A-A. The horizontally distributed profiles are 10 m AGL in height. (<b>d</b>) The vertical <math display="inline"><semantics> <mi>k</mi> </semantics></math> profile in the hilltop HT site.</p>
Full article ">Figure 13
<p>The horizontally distributed contours of wind velocities (<b>a</b>,<b>b</b>) and turbulent kinetic energy (<b>c</b>,<b>d</b>) are configured using the modified near-wall functions.</p>
Full article ">Figure 14
<p>The modified turbulence model constants (MC) are configured in the numerical case. (<b>a</b>,<b>b</b>) are <math display="inline"><semantics> <mrow> <mo>Δ</mo> <mi>S</mi> </mrow> </semantics></math> along line A-A and line AA-AA. (<b>c</b>) <math display="inline"><semantics> <mi>k</mi> </semantics></math> along lines A-A. The horizontally distributed profiles are 10 m AGL in height. (<b>d</b>) The vertical <math display="inline"><semantics> <mi>k</mi> </semantics></math> profile in the hilltop HT site.</p>
Full article ">Figure 15
<p>As for <a href="#atmosphere-15-01259-f011" class="html-fig">Figure 11</a>, the horizontally distributed contours of wind velocities (<b>a</b>,<b>b</b>) and turbulent kinetic energy (<b>c</b>,<b>d</b>) are configured using the modified turbulence closures.</p>
Full article ">Figure 15 Cont.
<p>As for <a href="#atmosphere-15-01259-f011" class="html-fig">Figure 11</a>, the horizontally distributed contours of wind velocities (<b>a</b>,<b>b</b>) and turbulent kinetic energy (<b>c</b>,<b>d</b>) are configured using the modified turbulence closures.</p>
Full article ">
31 pages, 33586 KiB  
Article
A Fuzzy Pure Pursuit for Autonomous UGVs Based on Model Predictive Control and Whole-Body Motion Control
by Yaoyu Sui, Zhong Yang, Haoze Zhuo, Yulong You, Wenqiang Que and Naifeng He
Drones 2024, 8(10), 554; https://doi.org/10.3390/drones8100554 - 6 Oct 2024
Cited by 1 | Viewed by 1129
Abstract
In this paper, we propose an adaptive fuzzy pure pursuit trajectory tracking algorithm for autonomous unmanned ground vehicles (UGVs), addressing the challenges of accurate and stable navigation in complex environments. Traditional pure pursuit methods with fixed look-ahead distances struggle to maintain precision in [...] Read more.
In this paper, we propose an adaptive fuzzy pure pursuit trajectory tracking algorithm for autonomous unmanned ground vehicles (UGVs), addressing the challenges of accurate and stable navigation in complex environments. Traditional pure pursuit methods with fixed look-ahead distances struggle to maintain precision in dynamic and uneven terrains. Our approach uniquely integrates a fuzzy control algorithm that allows for real-time adjustments of the look-ahead distance based on environmental feedback, thereby enhancing tracking accuracy and smoothness. Additionally, we combine this with model predictive control (MPC) and whole-body motion control (WBC), where MPC forecasts future states and optimally adjusts control actions, while WBC ensures coordinated motion of the UGV, maintaining balance and stability, especially in rough terrains. This integration not only improves responsiveness to changing conditions but also enables dynamic balance adjustments during movement. The proposed algorithm was validated through simulations in Gazebo and real-world experiments on physical platforms. In real-world tests, our algorithm reduced the average trajectory tracking error by 45% and the standard deviation by nearly 50%, significantly improving stability and accuracy compared to traditional methods. Full article
(This article belongs to the Special Issue Advances in Guidance, Navigation, and Control)
Show Figures

Figure 1

Figure 1
<p>The quadruped robot used in the physical experiments is about to take on the challenge of tracking an S-shaped curve in complex outdoor terrain.</p>
Full article ">Figure 2
<p>Simplified model of quadruped robot.</p>
Full article ">Figure 3
<p>Pure pursuit for quadruped robots.</p>
Full article ">Figure 4
<p>Fuzzy logic diagram.</p>
Full article ">Figure 5
<p>Fuzzy inputs and outputs.</p>
Full article ">Figure 6
<p>The control quantity surfaces for <math display="inline"><semantics> <mrow> <mo>Δ</mo> <msub> <mi>k</mi> <mi>v</mi> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>Δ</mo> <msub> <mi>k</mi> <mi>w</mi> </msub> </mrow> </semantics></math> are listed as (<b>a</b>) the control quantity surfaces for <math display="inline"><semantics> <mrow> <mo>Δ</mo> <msub> <mi>k</mi> <mi>v</mi> </msub> </mrow> </semantics></math> and (<b>b</b>) the control quantity surfaces for <math display="inline"><semantics> <mrow> <mo>Δ</mo> <msub> <mi>k</mi> <mi>w</mi> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 7
<p>MPC–WBC framework diagram.</p>
Full article ">Figure 8
<p>Simulation in Gazebo.</p>
Full article ">Figure 9
<p>Quadruped robot used for simulation.</p>
Full article ">Figure 10
<p>Ground undulation 3 cm test (walk in line).</p>
Full article ">Figure 11
<p>Robot desired trajectory vs. actual trajectory (3 cm, line).</p>
Full article ">Figure 12
<p>Position errors, angular error, and velocities during robot walking (3 cm, line).</p>
Full article ">Figure 13
<p>Three joint torques of one leg (3 cm, line).</p>
Full article ">Figure 14
<p>(<b>a</b>) Value of <math display="inline"><semantics> <mrow> <msub> <mi>k</mi> <mi>v</mi> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>k</mi> <mi>w</mi> </msub> </mrow> </semantics></math>; (<b>b</b>) value of <math display="inline"><semantics> <mi>l</mi> </semantics></math>.</p>
Full article ">Figure 15
<p>Ground undulation 6 cm test (work in line).</p>
Full article ">Figure 16
<p>Robot desired trajectory vs. actual trajectory (6 cm, line).</p>
Full article ">Figure 17
<p>Position errors, angular error, and velocities during robot walking (6 cm, line).</p>
Full article ">Figure 18
<p>Three joint torques of one leg (6 cm, line).</p>
Full article ">Figure 19
<p>(<b>a</b>) Value of <math display="inline"><semantics> <mrow> <msub> <mi>k</mi> <mi>v</mi> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>k</mi> <mi>w</mi> </msub> </mrow> </semantics></math>; (<b>b</b>) value of <math display="inline"><semantics> <mi>l</mi> </semantics></math>.</p>
Full article ">Figure 20
<p>Ground undulation 6 cm test (walk in circle).</p>
Full article ">Figure 21
<p>Robot desired trajectory vs. actual trajectory (6 cm, circle).</p>
Full article ">Figure 22
<p>Position errors, angular error, and velocities during robot walking (6 cm, circle).</p>
Full article ">Figure 23
<p>Three joint torques of one leg (6 cm, circle).</p>
Full article ">Figure 24
<p>(<b>a</b>) Value of <math display="inline"><semantics> <mrow> <msub> <mi>k</mi> <mi>v</mi> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>k</mi> <mi>w</mi> </msub> </mrow> </semantics></math>; (<b>b</b>) value of <math display="inline"><semantics> <mi>l</mi> </semantics></math>.</p>
Full article ">Figure 25
<p>USLGO1.</p>
Full article ">Figure 26
<p>Demonstration of experimental robot walking routes in complex outdoor environments. (<b>a</b>–<b>d</b>) represent the 4 parts of the S-curve respectively.</p>
Full article ">Figure 27
<p>Comparison of real-world robot trajectories with different algorithms.</p>
Full article ">Figure 28
<p>Position errors, angular error, and velocities during robot walking (in real situation).</p>
Full article ">Figure 29
<p>(<b>a</b>) Three joint torques of one leg (in real situation); (<b>b</b>) value of <math display="inline"><semantics> <mi>l</mi> </semantics></math>.</p>
Full article ">
16 pages, 34354 KiB  
Article
Autonomous Vehicles Traversability Mapping Fusing Semantic–Geometric in Off-Road Navigation
by Bo Zhang, Weili Chen, Chaoming Xu, Jinshi Qiu and Shiyu Chen
Drones 2024, 8(9), 496; https://doi.org/10.3390/drones8090496 - 18 Sep 2024
Viewed by 1606
Abstract
This paper proposes an evaluating and mapping methodology of terrain traversability for off-road navigation of autonomous vehicles in unstructured environments. Terrain features are extracted from RGB images and 3D point clouds to create a traversal cost map. The cost map is then employed [...] Read more.
This paper proposes an evaluating and mapping methodology of terrain traversability for off-road navigation of autonomous vehicles in unstructured environments. Terrain features are extracted from RGB images and 3D point clouds to create a traversal cost map. The cost map is then employed to plan safe trajectories. Bayesian generalized kernel inference is employed to assess unknown grid attributes due to the sparse raw point cloud data. A Kalman filter also creates density local elevation maps in real time by fusing multiframe information. Consequently, the terrain semantic mapping procedure considers the uncertainty of semantic segmentation and the impact of sensor noise. A Bayesian filter is used to update the surface semantic information in a probabilistic manner. Ultimately, the elevation map is utilized to extract geometric characteristics, which are then integrated with the probabilistic semantic map. This combined map is then used in conjunction with the extended motion primitive planner to plan the most effective trajectory. The experimental results demonstrate that the autonomous vehicles obtain a success rate enhancement ranging from 4.4% to 13.6% and a decrease in trajectory roughness ranging from 5.1% to 35.8% when compared with the most developed outdoor navigation algorithms. Additionally, the autonomous vehicles maintain a terrain surface selection accuracy of over 85% during the navigation process. Full article
Show Figures

Figure 1

Figure 1
<p>Navigation of an autonomous vehicle in unstructured wild environments with different terrains and obstacles.</p>
Full article ">Figure 2
<p>Autonomous navigation framework overview. The algorithm evaluates terrain traversability using RGB images and 3D point clouds. The traversability map guarantees safe and smooth path planning.</p>
Full article ">Figure 3
<p>Semantic cost projection. Some examples of projection: green represents easily traversable areas such as asphalt and concrete, yellow represents rough surfaces such as grass and sand, orange represents rocks, and blue represents obstacles.</p>
Full article ">Figure 4
<p>Comparison of obstacle identification results of the three methods in two real-world scenarios. In the scenario, Case A (first row) has two hidden obstacle areas (two puddle areas 1, 2); Case B (second row) has three false obstacles (three individual tall grassy areas 3, 4, and 5). (<b>a</b>) RGB image inputs (corresponding to the two images in the first column) for both cases; (<b>b</b>) semantic segmentation for both cases; (<b>c</b>) geometric gridmaps for both cases; (<b>d</b>) fused passivity maps for both cases.</p>
Full article ">Figure 5
<p>Results of evaluating the terrain for geometry only (top) and geometric-semantic fusion (bottom). (<b>a</b>) Scenario 1, (<b>b</b>) Scenario 2, and (<b>c</b>) Scenario 3. The areas that the geometric method gets wrong are shown in blue, and the areas where the geometry is neglected are displayed in red.</p>
Full article ">Figure 6
<p>Autonomous vehicle trajectories when navigating in four various unstructured scenarios using proposed method (black), RSPMP (blue), PUTN (red), BGK (orange), and DWA (yellow). (<b>a</b>) Scenario 1; (<b>b</b>) Scenario 2; (<b>c</b>) Scenario 3; (<b>d</b>) Scenario 4. It can be observed that proposed method allows the autonomous vehicle to navigate on smooth, low-cost surfaces and maintains a short trajectory.</p>
Full article ">Figure 7
<p>Quantitative comparison of navigation performance between various methods. (<b>a</b>) Scenario 1; (<b>b</b>) Scenario 2; (<b>c</b>) Scenario 3; (<b>d</b>) Scenario 4.</p>
Full article ">
19 pages, 9948 KiB  
Article
Traversability Analysis and Path Planning for Autonomous Wheeled Vehicles on Rigid Terrains
by Nan Wang, Xiang Li, Zhe Suo, Jiuchen Fan, Jixin Wang and Dongxuan Xie
Drones 2024, 8(9), 419; https://doi.org/10.3390/drones8090419 - 23 Aug 2024
Viewed by 1219
Abstract
Autonomous vehicles play a crucial role in three-dimensional transportation systems and have been extensively investigated and implemented in mining and other fields. However, the diverse and intricate terrain characteristics present challenges to vehicle traversability, including complex geometric features such as slope, harsh physical [...] Read more.
Autonomous vehicles play a crucial role in three-dimensional transportation systems and have been extensively investigated and implemented in mining and other fields. However, the diverse and intricate terrain characteristics present challenges to vehicle traversability, including complex geometric features such as slope, harsh physical parameters such as friction and roughness, and irregular obstacles. The current research on traversability analysis primarily emphasizes the processing of perceptual information, with limited consideration for vehicle performance and state parameters, thereby restricting their applicability in path planning. A framework of traversability analysis and path planning methods for autonomous wheeled vehicles on rigid terrains is proposed in this paper for better traversability costs and less redundancy in path planning. The traversability boundary conditions are established first based on terrain and vehicle characteristics using theoretical methods to determine the traversable areas. Then, the traversability cost map for the traversable areas is obtained through simulation and segmented linear regression analysis. Afterward, the TV-Hybrid A* algorithm is proposed by redefining the path cost functions of the Hybrid A* algorithm through the simulation data and neural network method to generate a more cost-effective path. Finally, the path generated by the TV-Hybrid A* algorithm is validated and compared with that of the A* and Hybrid A* algorithms in simulations, demonstrating a slightly better traversability cost for the former. Full article
Show Figures

Figure 1

Figure 1
<p>An example of rigid terrain composed of compacted soil.</p>
Full article ">Figure 2
<p>The proposed framework of traversability analysis and path planning methods for autonomous wheeled vehicles on rigid terrains.</p>
Full article ">Figure 3
<p>The relationship among vehicle speed, vehicle load, and terrain parameters: (<b>a</b>) The maximum average driving speed of vehicles under different loads on varying terrain slopes with favorable friction and roughness conditions; (<b>b</b>) The maximum average driving speed of vehicles under different loads on varying friction conditions of flat terrain with favorable roughness conditions; (<b>c</b>) The maximum average driving speed of vehicles under different loads on varying roughness conditions of flat terrain with favorable friction conditions. The <span class="html-italic">x</span>-axis values 1–8 correspond to the roughness categories A–H.</p>
Full article ">Figure 3 Cont.
<p>The relationship among vehicle speed, vehicle load, and terrain parameters: (<b>a</b>) The maximum average driving speed of vehicles under different loads on varying terrain slopes with favorable friction and roughness conditions; (<b>b</b>) The maximum average driving speed of vehicles under different loads on varying friction conditions of flat terrain with favorable roughness conditions; (<b>c</b>) The maximum average driving speed of vehicles under different loads on varying roughness conditions of flat terrain with favorable friction conditions. The <span class="html-italic">x</span>-axis values 1–8 correspond to the roughness categories A–H.</p>
Full article ">Figure 4
<p>A three-layer neural network to analyze the impact of terrain parameters on path costs.</p>
Full article ">Figure 5
<p>The simulation validation framework for traversability analysis and path planning methods on rigid terrains.</p>
Full article ">Figure 6
<p>The initial environmental information for the simulation validation: (<b>a</b>) The terrain DEM model; (<b>b</b>) The terrain height map obtained from the DEM model; (<b>c</b>) The terrain gradient map obtained from the DEM model; (<b>d</b>) The terrain physical distribution information of friction marked in gray with roughness marked with filling; (<b>e</b>) The terrain obstacles with height and type information.</p>
Full article ">Figure 6 Cont.
<p>The initial environmental information for the simulation validation: (<b>a</b>) The terrain DEM model; (<b>b</b>) The terrain height map obtained from the DEM model; (<b>c</b>) The terrain gradient map obtained from the DEM model; (<b>d</b>) The terrain physical distribution information of friction marked in gray with roughness marked with filling; (<b>e</b>) The terrain obstacles with height and type information.</p>
Full article ">Figure 6 Cont.
<p>The initial environmental information for the simulation validation: (<b>a</b>) The terrain DEM model; (<b>b</b>) The terrain height map obtained from the DEM model; (<b>c</b>) The terrain gradient map obtained from the DEM model; (<b>d</b>) The terrain physical distribution information of friction marked in gray with roughness marked with filling; (<b>e</b>) The terrain obstacles with height and type information.</p>
Full article ">Figure 7
<p>The simulation model of the articulated chassis vehicle on the designed terrain.</p>
Full article ">Figure 8
<p>Traversability boundary map with black areas indicating impassable areas.</p>
Full article ">Figure 9
<p>The traversability cost maps with the impassable areas marked in black: (<b>a</b>) Height cost map; (<b>b</b>) Slope cost map along <span class="html-italic">x</span>-direction gradient; (<b>c</b>) Slope cost map along <span class="html-italic">y</span>-direction gradient; (<b>d</b>) Friction cost map; (<b>e</b>) Roughness cost map; (<b>f</b>) Obstacles cost map.</p>
Full article ">Figure 9 Cont.
<p>The traversability cost maps with the impassable areas marked in black: (<b>a</b>) Height cost map; (<b>b</b>) Slope cost map along <span class="html-italic">x</span>-direction gradient; (<b>c</b>) Slope cost map along <span class="html-italic">y</span>-direction gradient; (<b>d</b>) Friction cost map; (<b>e</b>) Roughness cost map; (<b>f</b>) Obstacles cost map.</p>
Full article ">Figure 10
<p>The paths generated by the Hybrid A* and TV-Hybrid A* planners.</p>
Full article ">
10 pages, 1022 KiB  
Technical Note
A Simple Path to the Small Perturbation Method for Scattering from Slightly Rough Dielectric Surfaces
by Antonio Iodice and Pasquale Imperatore
Remote Sens. 2024, 16(16), 3035; https://doi.org/10.3390/rs16163035 - 18 Aug 2024
Viewed by 1173
Abstract
We propose a perturbative method to compute electromagnetic scattering from slightly rough dielectric surfaces, which leads to the same result as the usual Small Perturbation Method (SPM) in a surprisingly simple way. The proposed method is based on three pillars: the volumetric perturbative [...] Read more.
We propose a perturbative method to compute electromagnetic scattering from slightly rough dielectric surfaces, which leads to the same result as the usual Small Perturbation Method (SPM) in a surprisingly simple way. The proposed method is based on three pillars: the volumetric perturbative approach, the reciprocity theorem, and a proper approximation of the electric field within the perturbation volume, that we name Internal Field Approximation (IFA). The proposed new mathematical derivation of the SPM turns out to be much simpler and more concise than the classical one. In addition, being based on a volumetric perturbation approach, it has the potential of dealing in future with surface and volume scattering within a unitary framework, which is useful in modelling scattering from, e.g., vegetated soil, snow-covered terrain, and inhomogeneous soils. Therefore, although the presented result is mainly theoretical, it can have important applications in remote sensing. Full article
(This article belongs to the Section Engineering Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Geometry of the problem.</p>
Full article ">Figure 2
<p>Decomposition scheme for the perturbation volume.</p>
Full article ">
24 pages, 24217 KiB  
Article
Evaluating the Impact of DEM Spatial Resolution on 3D Rockfall Simulation in GIS Environment
by Maria P. Kakavas, Paolo Frattini, Alberto Previati and Konstantinos G. Nikolakopoulos
Geosciences 2024, 14(8), 200; https://doi.org/10.3390/geosciences14080200 - 29 Jul 2024
Cited by 1 | Viewed by 1206
Abstract
Rockfalls are natural geological phenomena characterized by the abrupt detachment and freefall descent of rock fragments from steep slopes. These events exhibit considerable variability in scale, velocity, and trajectory, influenced by the geological composition of the slope, the topography, and other environmental conditions. [...] Read more.
Rockfalls are natural geological phenomena characterized by the abrupt detachment and freefall descent of rock fragments from steep slopes. These events exhibit considerable variability in scale, velocity, and trajectory, influenced by the geological composition of the slope, the topography, and other environmental conditions. By employing advanced modeling techniques and terrain analysis, researchers aim to predict and control rockfall hazards to prevent casualties and protect properties in areas at risk. In this study, two rockfall events in the villages of Myloi and Platiana of Ilia prefecture were examined. The research was conducted by means of HY-STONE software, which performs 3D numerical modeling of the motion of non-interacting blocks. To perform this modeling, input files require the processing of base maps and datasets in a GIS environment. Stochastic modeling and 3D descriptions of slope topography, based on Digital Elevation Models (DEMs) without spatial resolution limitations, ensure multiscale analysis capabilities. Considering this capability, seven freely available DEMs, derived from various sources, were applied in HY-STONE with the scope of performing a large number of multiparametric analyses and selecting the most appropriate and efficient DEM for the software requirements. All the necessary data for the multiparametric analyses were generated within a GIS environment, utilizing either the same restitution coefficients and rolling friction coefficient or varying ones. The results indicate that finer-resolution DEMs capture detailed terrain features, enabling the precise identification of rockfall source areas and an accurate depiction of the kinetic energy distribution. Further, the results show that a correct application of the model to different DEMs requires a specific parametrization to account for the different roughness of the models. Full article
(This article belongs to the Special Issue Earth Observation by GNSS and GIS Techniques)
Show Figures

Figure 1

Figure 1
<p>Pictures of the rockfall events in the Myloi and Platiana region in relation to the Hellenic region.</p>
Full article ">Figure 2
<p>Pictures of the rockfall in the Myloi region illustrate the final position of rock blocks. Red arrows indicate the final positions of the rock masses at the conclusion of the rockfall event.</p>
Full article ">Figure 3
<p>Picture taken from the Platiana region after the rockfall event. The slope picture is taken from the Greek Cadastral (2008).</p>
Full article ">Figure 4
<p>Reforestation procedures in the slope next to the Platiana village. The slope pictures are taken from the Greek Cadastral.</p>
Full article ">Figure 5
<p>Myloi area before (<b>left</b> image) and after (<b>right</b> image) removing the vegetation through Cloud Compare Software.</p>
Full article ">Figure 6
<p>HY-STONE results for the site of Myloi using parameters (E<sub>N</sub>, E<sub>T</sub>, and A<sub>T</sub>) calibrated on the UAV DEM. The transit frequencies of blocks are shown for (<b>a</b>) UAV DEM, (<b>b</b>) Greek Cadastral DEM, (<b>c</b>) ALOS AW3D30 DEM, (<b>d</b>) ASTER GDEM, and (<b>e</b>) SRTM30 DEM.</p>
Full article ">Figure 7
<p>HY-STONE results for the site of Myloi using parameters (E<sub>N</sub>, E<sub>T</sub>, and A<sub>T</sub>) calibrated on the UAV DEM. The maximum translation kinetic energies are shown for (<b>a</b>) UAV DEM, (<b>b</b>) Greek Cadastral DEM, (<b>c</b>) ALOS AW3D30 DEM, (<b>d</b>) ASTER GDEM, and (<b>e</b>) SRTM30 DEM.</p>
Full article ">Figure 8
<p>Slope profiles from (<b>a</b>) UAV DSE, (<b>b</b>) Greek Cadastral DEM, (<b>c</b>) ALOS AW3D30 DEM, (<b>d</b>) ASTER GDEM, and (<b>e</b>) SRTM30 DEM.</p>
Full article ">Figure 9
<p>HY-STONE results, in terms of transit frequency for the site of Platiana. (<b>a</b>) Greek Cadastral DEM, (<b>b</b>) ALOS AW3D30 DEM, (<b>c</b>) ASTER GDEM, (<b>d</b>) SRTM30 DEM, (<b>e</b>) SRTM90 DEM, and (<b>f</b>) TanDEM_X.</p>
Full article ">Figure 10
<p>HY-STONE results in terms of the maximum translation kinetic energy for the site of Platiana. (<b>a</b>) Greek Cadastral DEM, (<b>b</b>) ALOS AW3D30 DEM, (<b>c</b>) ASTER GDEM, (<b>d</b>) SRTM30 DEM, (<b>e</b>) SRTM90 DEM, and (<b>f</b>) TanDEM_X.</p>
Full article ">Figure 11
<p>Slope profiles from (<b>a</b>) Greek Cadastral DEM, (<b>b</b>) ALOS AW3D30 DEM, (<b>c</b>) ASTER GDEM, (<b>d</b>) SRTM30 DEM, (<b>e</b>) SRTM90 DEM, and (<b>f</b>) TanDEM_X.</p>
Full article ">Figure 12
<p>HY-STONE results by using the Greek Cadastral DEM with the optimal coefficients shown in <a href="#geosciences-14-00200-t003" class="html-table">Table 3</a> for the site of Myloi: (<b>a</b>) the transit frequency and (<b>b</b>) the maximum translation kinetic energy.</p>
Full article ">Figure 13
<p>HY-STONE results by using the ALOS AW3D30 DEM with optimal coefficients shown in <a href="#geosciences-14-00200-t004" class="html-table">Table 4</a> for the site of Platiana: (<b>a</b>) the transit frequency and (<b>b</b>) the maximum translation kinetic energy.</p>
Full article ">Figure 14
<p>HY-STONE results by using the ASTER GDEM with the optimal coefficients shown in <a href="#geosciences-14-00200-t005" class="html-table">Table 5</a> for the site of Platiana: (<b>a</b>) the transit frequency and (<b>b</b>) the maximum translation kinetic energy.</p>
Full article ">Figure 15
<p>HY-STONE results by using the SRTM30 DEM with the optimal coefficients shown in <a href="#geosciences-14-00200-t006" class="html-table">Table 6</a> for the site of Platiana: (<b>a</b>) the transit frequency and (<b>b</b>) the maximum translation kinetic energy.</p>
Full article ">Figure 16
<p>HY-STONE results by using the SRTM90 DEM with the optimal coefficients shown in <a href="#geosciences-14-00200-t007" class="html-table">Table 7</a> for the site of Platiana: (<b>a</b>) the transit frequency and (<b>b</b>) the maximum translation kinetic energy.</p>
Full article ">Figure 17
<p>HY-STONE results by using the TanDEM_X with the optimal coefficients shown in <a href="#geosciences-14-00200-t008" class="html-table">Table 8</a> for the site of Platiana: (<b>a</b>) the transit frequency and (<b>b</b>) the maximum translation kinetic energy.</p>
Full article ">Figure 18
<p>The effect of spatial resolution in the normal (E<sub>N</sub>) and tangential (E<sub>T</sub>) restitutions and the rolling friction (A<sub>T</sub>) coefficients. For 30 m and 90 m, three and two DEMS are available, respectively.</p>
Full article ">
Back to TopTop