[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
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

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (1,860)

Search Parameters:
Keywords = Discrete Element Method

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
17 pages, 3511 KiB  
Article
Numerical Simulation Analysis of Different Excavation Parameters for TBM 3D Disc Cutters Based on the Discrete Element Method
by Feng Liang, Chenyuan Pei, Weibang Luo, Minglong You and Fei Tan
Appl. Sci. 2025, 15(1), 38; https://doi.org/10.3390/app15010038 - 24 Dec 2024
Abstract
This study provides a theoretical foundation for optimizing tunnel boring machine (TBM) excavation parameters under diverse geological conditions, offering significant engineering value by enhancing construction efficiency and reducing costs. As the development of underground spaces advances, TBMs play a pivotal role in tunnel [...] Read more.
This study provides a theoretical foundation for optimizing tunnel boring machine (TBM) excavation parameters under diverse geological conditions, offering significant engineering value by enhancing construction efficiency and reducing costs. As the development of underground spaces advances, TBMs play a pivotal role in tunnel excavation. TBMs enhance safety in excavation by mechanically breaking rock, reducing the reliance on explosives, and the associated risks of blasts. The shield support minimizes surrounding rock collapse, advanced geological forecasting mitigates risks posed by complex geologies, and intelligent monitoring systems improve operational safety. To enhance TBM efficiency and safety, this study developed a 3D simulation model of rock breaking by disc cutters using the discrete element method. This study systematically examined the effects of excavation parameters, including disc-cutter diameter, cutter spacing, and penetration, on rock-breaking performance. The findings reveal, that as the disc-cutter diameter increases, the rolling force also increases, while the rock-breaking specific energy initially rises and then declines. The 19-inch disc cutter demonstrated a superior rock-breaking efficiency in conventional excavation operations. At a cutter spacing of 60 mm, the rock-breaking specific energy reached its lowest value, representing optimal efficiency. Furthermore, as the penetration increased, both the rolling force and rock fragmentation volume grew, whereas the specific energy decreased, further improving the rock-breaking efficiency. Full article
30 pages, 1856 KiB  
Article
Calibration of Discrete Element Simulation Parameters and Model Construction for the Interaction Between Coastal Saline Alkali Soil and Soil-Engaging Components
by Nan Xu, Zhenbo Xin, Jin Yuan, Zenghui Gao, Yu Tian, Chao Xia, Xuemei Liu and Dongwei Wang
Agriculture 2025, 15(1), 7; https://doi.org/10.3390/agriculture15010007 - 24 Dec 2024
Abstract
There are approximately 36.7 million hectares of saline alkali land available in China. To enhance the comprehensive utilization value of coastal saline alkali land and boost crop yields in such areas, it is essential to conduct research on optimizing the operational performance of [...] Read more.
There are approximately 36.7 million hectares of saline alkali land available in China. To enhance the comprehensive utilization value of coastal saline alkali land and boost crop yields in such areas, it is essential to conduct research on optimizing the operational performance of high-performance soil contact components in light of the soil characteristics of coastal saline alkali land. Discrete element simulation can be employed to investigate the operational mechanisms of various key components. Nevertheless, at present, there is a dearth of discrete element models for the key physical parameters and soil structure of coastal saline alkali land soil. In this article, typical coastal saline alkali field soil was sampled, and the physical properties of the saline alkali soil, including salt content, moisture content, particle size distribution, and particle size, as well as intrinsic parameters such as soil compaction, density, Poisson’s ratio, and shear modulus, were measured. The Hertz Mindlin with Bonding contact model was employed. Physical experiments on soil accumulation angles at different depths were carried out using the cylindrical lifting method. Subsequently, by means of the discrete element method and the BBD experimental design method, a response surface model was established, and an optimization analysis was performed on the optimal parameters for the soil–soil collision recovery coefficient, static friction coefficient, and dynamic friction coefficient at each depth. Test benches for measuring the collision recovery coefficient, static friction coefficient, and rolling friction coefficient of saline alkali soil at -65Mn were set up, calculation formulas for each parameter were derived, and the contact parameters between soil at different depths and 65Mn were obtained. The results of the sliding friction angle test on different depths of saline alkali soil at -65Mn were further verified using the discrete element method, with a maximum error of 3.11%, which falls within the allowable range. This suggests that the calibration results of the discrete element simulation parameters for the interaction between soil and contact components are reliable, providing data and model support for future research on enhancing the operational performance of high-performance contact components. Full article
(This article belongs to the Special Issue Intelligent Agricultural Equipment in Saline Alkali Land)
11 pages, 16045 KiB  
Article
Study of Ventilation Strategies in a Passenger Aircraft Cabin Using Numerical Simulation
by S. M. Abdul Khader, John Valerian Corda, Kevin Amith Mathias, Gowrava Shenoy, Kamarul Arifin bin Ahmad, Augustine V. Barboza, Sevagur Ganesh Kamath and Mohammad Zuber
Computation 2025, 13(1), 1; https://doi.org/10.3390/computation13010001 - 24 Dec 2024
Abstract
Aircraft cabins have high occupant densities and may introduce the risk of COVID-19 contamination. In this study, a segment of a Boeing 767 aircraft cabin with a mixing type of air distribution system was investigated for COVID-19 deposition. A section of a Boeing [...] Read more.
Aircraft cabins have high occupant densities and may introduce the risk of COVID-19 contamination. In this study, a segment of a Boeing 767 aircraft cabin with a mixing type of air distribution system was investigated for COVID-19 deposition. A section of a Boeing 737-300 cabin, featuring four rows with 28 box-shaped mannequins, was used for simulation. Conditioned air entered through ceiling inlets and exited near the floor, simulating a mixed air distribution system. Cough droplets were modeled using the Discrete Phase Model from two locations: the centre seat in the second row and the window seat in the fourth row. These droplets had a mean diameter of 90 µm, an exhalation velocity of 11.5 m/s and a flow rate of 8.5 L/s. A high-quality polyhedral mesh of about 7.5 million elements was created, with a skewness of 0.65 and an orthogonality of 0.3. The SIMPLE algorithm and a second-order upwind finite volume method were used to model airflow and droplet dynamics. It was found that the ceiling accounted for the maximum concentration followed by the seats. The concentration of deposits was almost 50% more when the source was at window as compared to the centre seat. The Covid particles resided for longer duration when the source was at the centre of the cabin than when it was located near the widow. Full article
(This article belongs to the Special Issue Advances in Computational Methods for Fluid Flow)
Show Figures

Figure 1

Figure 1
<p>3D CAD model of a section of a Boeing 767-300 cabin.</p>
Full article ">Figure 2
<p>Airflow distribution for a mixed ventilation system: (<b>a</b>) vector plot from the literature review [<a href="#B15-computation-13-00001" class="html-bibr">15</a>] and (<b>b</b>) vector plot from the current study.</p>
Full article ">Figure 3
<p>Contaminant transmission with the source at the second-row centre seat: (<b>a</b>) particle concentration plot and (<b>b</b>) cough droplet distribution.</p>
Full article ">Figure 4
<p>Contaminant transmission with the source at the fourth-row window seat: (<b>a</b>) particle concentration plot and (<b>b</b>) cough droplet distribution.</p>
Full article ">Figure 5
<p>Comparison of the particle concentration distribution with the source at different locations: (<b>a</b>) second-row centre seat and (<b>b</b>) fourth-row window seat (SV—side view; TV—top view).</p>
Full article ">Figure 6
<p>Comparison of cough droplet distribution in different zones of the aircraft cabin (isometric view): (<b>a</b>) source—-second-row centre seat and (<b>b</b>) source—fourth-row window seat.</p>
Full article ">Figure 7
<p>Comparison of the vector plot of the airflow and cough droplet streamlines: (<b>a</b>) source–second-row centre seat and (<b>b</b>) source—fourth-row window seat.</p>
Full article ">Figure 8
<p>Cough droplet deposits on various surfaces of the aircraft cabin from different contaminant sources.</p>
Full article ">Figure 9
<p>Cough droplet residence time in various zones of the aircraft cabin from different contaminant sources.</p>
Full article ">
25 pages, 23734 KiB  
Article
Automated Mulched Transplanting of Angelica Seedlings Using a Pneumatic Sowing Device
by Hengtai Wang, Wei Sun, Hucun Wang and Petru A. Simionescu
Agronomy 2024, 14(12), 3076; https://doi.org/10.3390/agronomy14123076 - 23 Dec 2024
Abstract
To address the challenges of labor-intensive, inefficient, and inconsistent manual hole sowing and transplanting of Angelica sinensis in rain-fed hilly regions of Northwest China, a pneumatic hole-sowing device was designed based on the principle of electromagnetically controlled, high-speed reciprocating cylinder motion. Considering the [...] Read more.
To address the challenges of labor-intensive, inefficient, and inconsistent manual hole sowing and transplanting of Angelica sinensis in rain-fed hilly regions of Northwest China, a pneumatic hole-sowing device was designed based on the principle of electromagnetically controlled, high-speed reciprocating cylinder motion. Considering the agronomic requirements for transplanting mulched Angelica sinensis, the device’s structure and operational parameters were optimized. The key mechanisms involved in hole sowing and seedling placement were analyzed. A pneumatic circuit system, controlled by a relay circuit, was established, and a hole-sowing mechanism with a delayed closure effect was designed. Using the Discrete Element Method (DEM) and Multi-Body Dynamics (MBD) coupling technology, a simulation of the hole-sowing process was conducted to evaluate the device’s performance and its impact on soil disturbance and hole reformation in the seedbed. Prototype device performance tests were conducted, using qualified seeding depth under mulch and hole spacing as indicators. When the theoretical hole spacing was 30 cm and the hole-sowing frequency was 60 plants/(min·row), the soil bin test results indicated a seeding depth qualification rate of 93%, a misalignment rate of 3%, and a spacing qualification rate of 83%; the field test results showed a qualified seeding depth rate under mulch of 96%, the hole misalignment rate was 5%, and the spacing qualified rate was 86%. The pneumatic hole-sowing device’s performance meets the agronomic requirements for vertical transplanting of Angelica sinensis seedlings. This research can serve as a reference for designing planting machinery for rhizomatous medicinal plants. Full article
(This article belongs to the Special Issue Advances in Data, Models, and Their Applications in Agriculture)
Show Figures

Figure 1

Figure 1
<p>Morphological parameters of Angelica Seedlings, with <span class="html-italic">L<sub>R</sub></span> root length, <span class="html-italic">L<sub>Z</sub></span> seedling length, <span class="html-italic">B</span> seedling width, <span class="html-italic">D<sub>r</sub></span> root diameter.</p>
Full article ">Figure 2
<p>Scatter plot of morphological parameters of 50 <span class="html-italic">Angelica</span> seedlings.</p>
Full article ">Figure 3
<p>Agronomy of transplanting <span class="html-italic">Angelica sinensis</span> with plastic film.</p>
Full article ">Figure 4
<p>Structure diagram of the entire machine.</p>
Full article ">Figure 5
<p>Detailed view of the pneumatic hole-sowing device components.</p>
Full article ">Figure 6
<p>Split-cone hole-sowing device.</p>
Full article ">Figure 7
<p>Delayed closing mechanism.</p>
Full article ">Figure 8
<p>Force analysis of the hole maker during soil entry.</p>
Full article ">Figure 9
<p>Stress analysis of seedlings. 1. Hole maker. 2. Angelica seedling. Note: G: Weight of the Angelica seedlings, N; <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>F</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> </mrow> </semantics></math>: Support force from the left inner wall, N; <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>F</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math>: Support force from the right inner wall, N; <span class="html-italic">F</span>: Resultant force of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>F</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>F</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math>, N.</p>
Full article ">Figure 10
<p>Schematic diagram of pneumatic system. 1. Air compressor. 2. Air storage tank. 3. Stainless steel ball valve. 4. Air treatment unit. 5. Speed control valve. 6. Solenoid valve. 7. Quick exhaust valve. 8. Muffler. 9. High-speed cylinder. 10. Drain valve.</p>
Full article ">Figure 11
<p>Pneumatic and electric circuit design in FluidSim-P.</p>
Full article ">Figure 12
<p>Control system circuit flow chart.</p>
Full article ">Figure 13
<p>Schematic circuit diagram of control system.</p>
Full article ">Figure 14
<p>Dynamic model constraint addition.</p>
Full article ">Figure 15
<p>Dynamic model of the hole-sowing device.</p>
Full article ">Figure 16
<p>EDEM–RecurDyn coupling process.</p>
Full article ">Figure 17
<p>Soil disturbance states at different times.</p>
Full article ">Figure 18
<p>Effect of hole formation after soil reflux.</p>
Full article ">Figure 19
<p>TCC electric four-drive soil bin test trolley structure. 1. Rubber tires. 2. Hydraulic pump. 3. Oil tank. 4. Hydraulic power output drive motor. 5. Brake air pump. 6. Walking drive motor. 7. Water tank. 8. Lifting device oil pump. 9. Sprinkling rod. 10. Rear suspension. 11. Semi-mounted lifting device.</p>
Full article ">Figure 20
<p>The division of the test area.</p>
Full article ">Figure 21
<p>Field test.</p>
Full article ">Figure 22
<p>Field experiment results.</p>
Full article ">
16 pages, 4871 KiB  
Article
Research on the Method for Solving the Safety Factor of Rock Slope Based on Deformation Monitoring Warning Threshold
by Xiaoyan Wei and Xiuli Zhang
Appl. Sci. 2024, 14(24), 12061; https://doi.org/10.3390/app142412061 - 23 Dec 2024
Abstract
In view of the fact that field monitoring information can more intuitively and accurately reflect the stability state of slopes, this paper takes the warning threshold of slope deformation rate monitoring as the slope stability evaluation standard, and puts forward a method for [...] Read more.
In view of the fact that field monitoring information can more intuitively and accurately reflect the stability state of slopes, this paper takes the warning threshold of slope deformation rate monitoring as the slope stability evaluation standard, and puts forward a method for solving the safety coefficient of rocky slopes. The discrete element method (3DEC), which is suitable for rocky slopes, is selected as the numerical calculation tool, the convergence criterion of its strength reduction method is modified to the slope deformation rate threshold, and the method is realized by the bifurcation method through its built-in FISH programming language. The results of the classical case show that, by selecting the slope deformation rate threshold in the appropriate interval, the results of this paper’s method are very close to those of the finite unit stress method and the limit equilibrium method, verifying the reliability of this paper’s method. Further, the method of this paper is applied to an open-pit mine slope project, based on the slope deformation on-site monitoring data and through the time series prediction method to determine the slope deformation rate warning threshold, using the threshold as an evaluation criterion to solve the slope’s coefficient of safety. The calculation results show that the slope’s coefficient of safety in natural working conditions is 1.086, being basically stable. However, with continuous rainfall, the slope’s body gradually becomes saturated, the slope’s coefficient of safety is reduced to 0.987, and the slope’s safety is reduced to 0.987. After continuous rainfall and gradual saturation of the slope, the coefficient of safety decreases to 0.987, resulting in destabilization and destruction, which is consistent with the site conditions. Full article
(This article belongs to the Section Civil Engineering)
21 pages, 13380 KiB  
Article
Macro-Mesoscopic Failure Mechanism Based on a Direct Shear Test of a Cemented Sand and Gravel Layer
by Long Qian, Xingwen Guo, Qinghui Liu, Xin Cai and Xiaochuan Zhang
Buildings 2024, 14(12), 4078; https://doi.org/10.3390/buildings14124078 - 23 Dec 2024
Abstract
In order to explore the influence of different layer treatment methods on the macro- and meso-mechanical properties of cemented sand and gravel (CSG), in this paper, the shear behavior of CSG material was simulated by a three-dimensional particle flow program (PFC3D) based on [...] Read more.
In order to explore the influence of different layer treatment methods on the macro- and meso-mechanical properties of cemented sand and gravel (CSG), in this paper, the shear behavior of CSG material was simulated by a three-dimensional particle flow program (PFC3D) based on the results of direct shear test in the laboratory. In shear tests, untreated CSG samples with interface coating mortar and chiseling were used, and granular discrete element software (PDC3D 7.0) was used to establish mesoscopic numerical models of CSG samples with the above three interface treatment methods, in order to reveal the effects of interface treatment methods on the interface strength and damage mechanism of CSG samples. The results show that, with the increase in normal stress, the amount of aggregate falling off the shear failure surface increases, the bump and undulation are more obvious, and the failure mode of the test block is inferred to be extrusion friction failure. The shear strength of the mortar interface is 40% higher than that of the untreated interface, and the failure surface is smooth and flat under different normal stresses. The shear strength of the chiseled interface is 10% higher than that of the untreated interface, and the failure surface fluctuates significantly under different normal stresses. Through the analysis of the fracture evolution process in the numerical simulation, it is found that the fracture of the sample at the mortar interface mainly expands along the mortar–aggregate interface and the damage mode is shear slip. However, the cracks of the samples at the gouged interface are concentrated on the upper and lower sides of the interface, and the damage mode is tension–shear. The failure mode of the samples without surface treatment is mainly tensile and shear failure, and the failure mode gradually changes to extrusion friction failure. Full article
Show Figures

Figure 1

Figure 1
<p>Specimen material.</p>
Full article ">Figure 2
<p>Experimental procedure flowchart.</p>
Full article ">Figure 3
<p>Fabrication process of CSG specimens with no treatment forms.</p>
Full article ">Figure 4
<p>Fabrication process of CSG specimens with different interface treatment forms.</p>
Full article ">Figure 5
<p>Direct shear test device.</p>
Full article ">Figure 6
<p>Shear stress–shear displacement curve.</p>
Full article ">Figure 7
<p>Results of interface strength fitting.</p>
Full article ">Figure 8
<p>Interface with no treatment.</p>
Full article ">Figure 9
<p>Damage pattern of spreading mortar specimens.</p>
Full article ">Figure 10
<p>Damage pattern of chiseling specimens.</p>
Full article ">Figure 11
<p>Numerical model of the specimen.</p>
Full article ">Figure 12
<p>Comparison of numerical simulation results with experimental results for different interfaces of treatment.</p>
Full article ">Figure 13
<p>Number of fracture development in the specimen–shear displacement curve.</p>
Full article ">Figure 13 Cont.
<p>Number of fracture development in the specimen–shear displacement curve.</p>
Full article ">Figure 14
<p>Internal fracture distribution during shearing of the specimen.</p>
Full article ">Figure 14 Cont.
<p>Internal fracture distribution during shearing of the specimen.</p>
Full article ">Figure 15
<p>Thickness of fracture distribution at the end of specimen shear test.</p>
Full article ">
15 pages, 10783 KiB  
Article
Evaluation of the Effects of Rainfall Infiltration Boundaries on the Stability of Unsaturated Soil Slopes Using the Particle Flow Code
by Jian Zhang, Fangrui Hu, Qi Zhang, Jun Wang, Wenting Deng, Li Zhang and Xiaoquan Shao
Water 2024, 16(24), 3704; https://doi.org/10.3390/w16243704 - 22 Dec 2024
Viewed by 235
Abstract
Rainfall infiltration is the primary triggering factor for the instability of unsaturated slopes. At present, rainfall-induced landslides are mainly considered to be influenced by the overall infiltration conditions, while few investigations have been conducted on the influence of infiltration boundaries on slope instability. [...] Read more.
Rainfall infiltration is the primary triggering factor for the instability of unsaturated slopes. At present, rainfall-induced landslides are mainly considered to be influenced by the overall infiltration conditions, while few investigations have been conducted on the influence of infiltration boundaries on slope instability. This study proposes a rainfall infiltration method using a discrete element model (DEM), which is based on saturated–unsaturated seepage theory. The influence of three infiltration boundaries on the instability of homogeneous unsaturated soil slopes was studied. The results showed that the infiltration rate of a rainfall-covered slope crest was faster than that of rainfall-covered slope surfaces. A transient saturated zone was formed on the slope surface after a certain duration of rainfall. Rain infiltration boundary conditions significantly impact the saturation distribution, seepage field, failure mode, and failure period. The safety and stability factors for the rainfall-covered slope crest and full area decreased monotonically with the increase in rainfall duration, while there was a brief increase at the initial stage of rainfall before a quick decline for rainfall-covered slope surfaces. This research provides a preliminary exploration of the impact of rainfall boundary conditions on the instability of slopes, offering a reference basis for DEM simulations that consider slope stability under the influence of rainfall infiltration. Full article
Show Figures

Figure 1

Figure 1
<p>Flowchart of the analysis of seepage behavior with the DEM.</p>
Full article ">Figure 2
<p>The geometry and boundary conditions of a simple slope.</p>
Full article ">Figure 3
<p>The balls in the DEM models and monitoring points.</p>
Full article ">Figure 4
<p>The stress–strain curves of soil under different confining pressures.</p>
Full article ">Figure 5
<p>The variation in the degree of saturation during the rainfall process at different monitoring points: (<b>a</b>) rainfall-covered slope crest; (<b>b</b>) rainfall-covered slope surface; (<b>c</b>) full rainfall-covered slope area.</p>
Full article ">Figure 6
<p>Vector field of the seepage velocity of the slope in different rainfall infiltration scenarios: (<b>a</b>) rainfall-covered slope crest; (<b>b</b>) rainfall-covered slope surface; (<b>c</b>) full rainfall-covered slope area.</p>
Full article ">Figure 7
<p>Crack distribution, landslide debris, and slope displacement during slope failure: (<b>a</b>) rainfall-covered slope crest; (<b>b</b>) rainfall-covered slope surface; (<b>c</b>) full rainfall-covered slope area.</p>
Full article ">Figure 8
<p>Failure pattern obtained from the model test of a homogeneous sandy slope in [<a href="#B29-water-16-03704" class="html-bibr">29</a>]: (<b>a</b>) rainfall-covered slope crest; (<b>b</b>) rainfall-covered slope surface; (<b>c</b>) full rainfall-covered slope area.</p>
Full article ">Figure 9
<p>Evolution law of crack numbers in different rainfall infiltration scenarios.</p>
Full article ">Figure 10
<p>Evolution law of the slope safety factor during rainfall.</p>
Full article ">
20 pages, 4823 KiB  
Article
Design and Preliminary Evaluation of a Precision Cylindrical Air-Assisted Drill Sowing Device for Rapeseed, Wheat, and Rice
by Alfarog H. Albasheer, Qingxi Liao, Lei Wang, Anas Dafaallah Abdallah and Jianxin Lin
Agriculture 2024, 14(12), 2355; https://doi.org/10.3390/agriculture14122355 - 21 Dec 2024
Viewed by 316
Abstract
To address challenges in seed feeding stability and seeding uniformity in agricultural practices, this study aimed to introduce a cylindrical air-assisted drill sowing device (CADSD) designed for rapeseed, wheat, and rice (RWR). The device features a prototype hill-feeding mechanism that addresses problems related [...] Read more.
To address challenges in seed feeding stability and seeding uniformity in agricultural practices, this study aimed to introduce a cylindrical air-assisted drill sowing device (CADSD) designed for rapeseed, wheat, and rice (RWR). The device features a prototype hill-feeding mechanism that addresses problems related to seed feeding, airflow disruptions, and seed–wall collisions. Comprehensive bench tests, Discrete Element Method (DEM) simulations, and preliminary field experiments were conducted to evaluate the seed-feeding stability characteristics and optimize the structural parameters of the air-assisted drill sowing system, enhancing seeding uniformity and operational efficiency. The optimal operating speed range is between 4 and 5 km/h. When the seed feeding speed is 30 to 38 r/min, the coefficient of variation of the seed supply rate stability is less than 0.55%, and the relative error between the theoretical and the experimental actual values of the RWR supply rate regression model is less than 2%, further supporting the effectiveness of the device. A preliminary field test revealed a seeding uniformity coefficient of variation (CV) of 3.44% and an emergence rate of 88%, closely aligning with the desired metrics. The CADSD effectively sows multiple crop types with improved precision and uniformity, handling diverse seed types and sizes without requiring equipment modifications, highlighting its innovative impact on agricultural technology in the precise seeding of RWR. Full article
(This article belongs to the Section Agricultural Technology)
Show Figures

Figure 1

Figure 1
<p>Schematic diagram of an air-assisted centralized metering system for RWR. 1. Double disc diches. 2. Seed delivery pipe. 3. Distributor 4. Booster pipe. 5. Venturi tube 6. Centralize seed feeding device. 7. Seed box. 8. High-pressure blower.</p>
Full article ">Figure 2
<p>Sketch of the centralized seed-feeding device; (<b>a</b>) RWR shell sections and barrier Structure of centralized seed-feeding device; (<b>b</b>) Seed feeding mechanism. 1. Shell section barrier; 2. Wheat and rice section; 3. Rapeseed section; 4. Seed shell; 5. Seed layer; 6. Seed feeding mechanism; 7. Prototype hole wheel; 8. Seed fall mouth; 9. Blank wheel; 10. Wheat and rice hole wheel; 11. Division plate; 12. Rapeseed hole wheel; 13. Transmission shaft.</p>
Full article ">Figure 3
<p>Mechanics analysis of seeds population in filling zone. <span class="html-italic">ω</span> is the angular velocity of the cone hole wheel; <span class="html-italic">r</span> is the radius of the cone hole wheel; <span class="html-italic">β</span> is the seed’s initial filling position; <span class="html-italic">F<sub>f</sub></span> is the friction between microsegments and holes of the population; <span class="html-italic">F<sub>n</sub></span> is the lateral pressure of the population on the microsegments of the population; <span class="html-italic">F<sub>N</sub></span> is the type pore supports the population microsegment; <span class="html-italic">x</span> is <span class="html-italic">x</span>-axis; <span class="html-italic">G</span> is the gravity of population segments; <span class="html-italic">F<sub>c</sub></span> is the inertial centrifugal force of the population microsegment; <span class="html-italic">y</span> is <span class="html-italic">y</span>-axis.</p>
Full article ">Figure 4
<p>Mechanics analysis of seeds in the feeding zone. <span class="html-italic">ω</span> is the angular velocity; <span class="html-italic">ε</span> is the angle between the axis and the horizontal plane; <span class="html-italic">F<sub>f</sub></span> is the friction between microsegments and holes of the population; <span class="html-italic">θ</span> is the hole inclination angle; <span class="html-italic">F<sub>N1</sub></span> is the type pore supports the population microsegment; <span class="html-italic">x</span> is the <span class="html-italic">x</span>-axis; <span class="html-italic">F<sub>c</sub></span> is the inertial centrifugal force of the population microsegment; <span class="html-italic">y</span> is the <span class="html-italic">y</span>-axis.</p>
Full article ">Figure 5
<p>Sectional structure of model-hole. 1. Type hole; 2. Skew gear tooth; (<span class="html-italic">A</span>) the type hole width, mm; (<span class="html-italic">H</span>) the type hole depth, mm.</p>
Full article ">Figure 6
<p>DEM simulation model. (<b>a</b>) Rapeseed seed; (<b>b</b>) Wheat seed; (<b>c</b>) Rice seed; (<b>d</b>) Seed supply metering device.</p>
Full article ">Figure 7
<p>Seed metering performance test bench. 1. Distribution device; 2. Seed guide pipe; 3. Seed box; 4. Transmission pipeline; 5. RWR seed supply device; 6. Venturi tube.</p>
Full article ">Figure 8
<p>Field seeding test of rice.</p>
Full article ">Figure 9
<p>DEM simulation model stability test result for RWR supply at a different rotational speed. (<b>A</b>) Simulation test result of rapeseed supply at different rotational speeds. (<b>B</b>) Simulation test result of wheat supply at different rotational speeds. (<b>C</b>) Simulation test result of rice supply at different rotational speeds. (<b>a</b>) Seed supply at rotational speed 30 r/min. (<b>b</b>) Seed supply at rotational speed 38 r/min.</p>
Full article ">Figure 10
<p>Seeding performance of CADSD for RWR. (<b>A</b>) Rapeseed seeding performance. (<b>B</b>) Rice seeding performance. (<b>C</b>) Wheat seeding performance. (<b>a</b>) Total seed feed rates and coefficient of variation with respect to the working speed of the hole wheel; (<b>b</b>) Comparison of seed feed rates of each row; (<b>c</b>) Variations of seeding feed rates of each row with respect to the working speed of the hole wheel.</p>
Full article ">Figure 10 Cont.
<p>Seeding performance of CADSD for RWR. (<b>A</b>) Rapeseed seeding performance. (<b>B</b>) Rice seeding performance. (<b>C</b>) Wheat seeding performance. (<b>a</b>) Total seed feed rates and coefficient of variation with respect to the working speed of the hole wheel; (<b>b</b>) Comparison of seed feed rates of each row; (<b>c</b>) Variations of seeding feed rates of each row with respect to the working speed of the hole wheel.</p>
Full article ">Figure 10 Cont.
<p>Seeding performance of CADSD for RWR. (<b>A</b>) Rapeseed seeding performance. (<b>B</b>) Rice seeding performance. (<b>C</b>) Wheat seeding performance. (<b>a</b>) Total seed feed rates and coefficient of variation with respect to the working speed of the hole wheel; (<b>b</b>) Comparison of seed feed rates of each row; (<b>c</b>) Variations of seeding feed rates of each row with respect to the working speed of the hole wheel.</p>
Full article ">Figure 11
<p>Regression surface for the correlation Z = f (X, Y) between the quantity planted, seed density, and hole wheel rotation speed.</p>
Full article ">Figure 12
<p>Rice field experiment.</p>
Full article ">
28 pages, 1662 KiB  
Review
Numerical Simulation of Earthquake Impacts on Marine Structures: A Comprehensive Review
by Adel Kabi, Jersson X. Leon-Medina and Francesc Pozo
Buildings 2024, 14(12), 4039; https://doi.org/10.3390/buildings14124039 - 19 Dec 2024
Viewed by 322
Abstract
Marine and underwater structures, such as seawalls, piers, breakwaters, and pipelines, are particularly susceptible to seismic events. These events can directly damage the structures or destabilize their supporting soil through phenomena like liquefaction. This review examines advanced numerical modeling approaches, including CFD, FEM, [...] Read more.
Marine and underwater structures, such as seawalls, piers, breakwaters, and pipelines, are particularly susceptible to seismic events. These events can directly damage the structures or destabilize their supporting soil through phenomena like liquefaction. This review examines advanced numerical modeling approaches, including CFD, FEM, DEM, FVM, and BEM, to assess the impacts of earthquakes on these structures. These methods provide cost-effective and reliable simulations, demonstrating strong alignment with experimental and theoretical data. However, challenges persist in areas such as computational efficiency and algorithmic limitations. Key findings highlight the ability of these models to accurately simulate primary forces during seismic events and secondary effects, such as wave-induced loads. Nonetheless, discrepancies remain, particularly in capturing energy dissipation processes in existing models. Future advancements in computational capabilities and techniques, such as high-resolution DNS for wave–structure interactions and improved near-field seismoacoustic modeling show potential for enhancing simulation accuracy. Furthermore, integrating laboratory and field data into unified frameworks will significantly improve the precision and practicality of these models, offering robust tools for predicting earthquake and wave impacts on marine environments. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Hull with transverse stiffeners CAD detail, (<b>b</b>) Preparation of mesh and for a hull with transverse stiffeners and (<b>c</b>) result of the vibration mode of the hull transversely stiffened at frequency 11.209 Hz [<a href="#B17-buildings-14-04039" class="html-bibr">17</a>].</p>
Full article ">Figure 2
<p>CFD-DEM simulation of particle ejection test. (<b>a</b>) Setup and (<b>b</b>) particle motion trajectory with and without Magnus force [<a href="#B21-buildings-14-04039" class="html-bibr">21</a>].</p>
Full article ">Figure 3
<p>Marine current turbine. Wake geometry of IBEM model at different operating conditions. From left to right, TSR = 3, 6, 9. The diameter of the turbine was 700 mm [<a href="#B27-buildings-14-04039" class="html-bibr">27</a>].</p>
Full article ">Figure 4
<p>Schematic of numerical wave tank: (<b>a</b>) cross-section and (<b>b</b>) plan view [<a href="#B31-buildings-14-04039" class="html-bibr">31</a>].</p>
Full article ">Figure 5
<p>(<b>a</b>) STL files for the bottom geometry and cylinder, and (<b>b</b>) Computational domain with bottom slope and vertical cylinder [<a href="#B31-buildings-14-04039" class="html-bibr">31</a>]. The dimensions correspond to those described in <a href="#buildings-14-04039-f004" class="html-fig">Figure 4</a>.</p>
Full article ">Figure 6
<p>Results of waves2Foam simulations in four time steps from 31.10 s, until 31.90 s [<a href="#B31-buildings-14-04039" class="html-bibr">31</a>].</p>
Full article ">Figure 7
<p>Computational domain and coordinate system for DNS of wind over steep and breaking waves [<a href="#B36-buildings-14-04039" class="html-bibr">36</a>]. (<b>a</b>) 3D View of the waves, (<b>b</b>) 2D view dash line shows the level <math display="inline"><semantics> <mrow> <mi>z</mi> <mo>=</mo> <mn>0</mn> </mrow> </semantics></math> and red line indicate the wave.</p>
Full article ">Figure 8
<p>Beach profiles at Pont del Petroli. The original beach profile from the design report is indicated by a blue line. In red, the two profiles surveyed by LIM/UPC before and after storm Gloria [<a href="#B39-buildings-14-04039" class="html-bibr">39</a>].</p>
Full article ">Figure 9
<p>2D view of a pipe inside the lattice Boltzmann grid points [<a href="#B43-buildings-14-04039" class="html-bibr">43</a>].</p>
Full article ">
21 pages, 10204 KiB  
Article
Numerical Modeling of Engineering-Scale Jointed Coal Mass and Confining Pressure Effect
by Feiteng Zhang, Xiangyu Wang, Qingcang Wang, Jianbiao Bai, Dingchao Chen, Menglong Li and Shiqi Sun
Appl. Sci. 2024, 14(24), 11844; https://doi.org/10.3390/app142411844 - 18 Dec 2024
Viewed by 368
Abstract
In coal masses, joints serve as crucial components influencing mechanical responses and failure mechanisms. The joint distribution complicates the acquisition of physical parameters for engineering-scale coal masses, and traditional laboratory or in situ tests often suffer from inaccuracies and high costs. Therefore, examination [...] Read more.
In coal masses, joints serve as crucial components influencing mechanical responses and failure mechanisms. The joint distribution complicates the acquisition of physical parameters for engineering-scale coal masses, and traditional laboratory or in situ tests often suffer from inaccuracies and high costs. Therefore, examination of the dynamic properties of engineering-scale coal mass is generally performed by numerical simulations. This paper proposes a model construction approach using two-dimension Particle Flow Code (PFC2D) software to study the mechanical properties of engineering-scale jointed coal mass, addressing the limitations of conventional models by integrating scale effects, accuracy, and computational efficiency. Firstly, the distribution characteristics and mechanical parameters of the joints in the coal mass were obtained based on field statistics and laboratory experiments. The parameters of the laboratory-scale model were calibrated by the numerical matching method. The discrete element model for the engineering-scale coal mass was constructed by the step-by-step matching method. The confining pressure effect on the coal mass under a biaxial loading condition was studied, while the strength change, fissure evolution, and failure mechanism under different confining pressures and fissure degrees were investigated. Based on the simulation results, a quantitative relationship was established between the mechanical parameters, fissure degree, and confining pressure under compression conditions. Ultimately, the failure zone ahead of the working face and the distribution of the abutment pressure were assessed using the mechanics parameters of coal masses with diverse joint distributions. Full article
Show Figures

Figure 1

Figure 1
<p>Construction of numerical model under the uncalibrated parameters.</p>
Full article ">Figure 2
<p>Numerical uniaxial compression of laboratory-scale coal sample (NT—numerical test, LT—laboratory test, 1—compression-shear failure, 2—inclined splitting failure, 3—nearly vertical splitting failure, 4—regional failure).</p>
Full article ">Figure 3
<p>Loading curves of numerical model of coal samples of different sizes.</p>
Full article ">Figure 4
<p>Numerical matching of coal direct shear test.</p>
Full article ">Figure 5
<p>Calculation of fissure degrees in different directions.</p>
Full article ">Figure 6
<p>Coal model sampling lines and description of fissure degree.</p>
Full article ">Figure 7
<p>Schematic diagram of the numerical model of uniaxial compression of coal under different fissure degrees.</p>
Full article ">Figure 8
<p>Biaxial loading stress curve of models with different fissure degrees.</p>
Full article ">Figure 9
<p>Change rules of deviatoric stress under different fissure degrees and confining pressures.</p>
Full article ">Figure 10
<p>Evolution rules of fissures under different fissure degrees and confining pressures (for the proportion of joint fissures, the solid line with corresponding color represents joint fissures, while the dashed line represents matrix fissures).</p>
Full article ">Figure 10 Cont.
<p>Evolution rules of fissures under different fissure degrees and confining pressures (for the proportion of joint fissures, the solid line with corresponding color represents joint fissures, while the dashed line represents matrix fissures).</p>
Full article ">Figure 11
<p>The failure mode under different fissure degrees and confining pressures.</p>
Full article ">Figure 11 Cont.
<p>The failure mode under different fissure degrees and confining pressures.</p>
Full article ">Figure 12
<p>The angle of main fissure under different confining pressures.</p>
Full article ">Figure 13
<p>Stress distribution of longwall top-coal caving.</p>
Full article ">Figure 14
<p>The distribution of failure tendency coefficient and stress in coal seam under different fissure degrees.</p>
Full article ">Figure 15
<p>Relationship between top-coal damage range and fissure degree.</p>
Full article ">
17 pages, 19002 KiB  
Article
Study on the Failure Mechanism and Movement Characteristics Prediction of Gongdang Landslide in Linzhi, China
by Yuezu Huang, Yuanzhong Li, Yubin Zhao, Faming Zhang, Xiaokai Li, Huaqing Zhang and Xiaolong Zhang
Water 2024, 16(24), 3649; https://doi.org/10.3390/w16243649 - 18 Dec 2024
Viewed by 338
Abstract
Instability of landslide accumulation bodies is one of the common geological hazards under the influence of rainfall and water impoundment, especially under the transformation of rainfall patterns caused by global climate changes. Owing to the fact that determining the landslide potential failure mode [...] Read more.
Instability of landslide accumulation bodies is one of the common geological hazards under the influence of rainfall and water impoundment, especially under the transformation of rainfall patterns caused by global climate changes. Owing to the fact that determining the landslide potential failure mode is vital for preventing landslide disasters, this paper takes the Gongdang landslide as the research object to study the landslide deformation mechanism and predict movement characteristics. Firstly, the geological conditions of the study area and landslide were determined according to the field investigations; secondly, the physical and mechanical parameters of the sliding mass were clarified through laboratory tests. Moreover, the particle flow code (PFC) method was utilized to simulate the potential failure process of the landslide based on the three-dimensional numerical model according to the geological features and the micro-parameters. The results showed that the landslide deformation process lasted approximately 640 s with the stage characteristics of displacement and velocity and presented the evolutionary process with the local instability deformation. The simulation results are of practical significance and application value by effectively illustrating the potential deformation and failure process of the Gongdang landslide, which provides a reference for predicting and preventing the potential failure process of geological hazards in similar engineering through field investigations, laboratory tests, and numerical simulation. Full article
Show Figures

Figure 1

Figure 1
<p>Location and image of the Gongdang landslide.</p>
Full article ">Figure 2
<p>The geological map of the Gongdang landslide.</p>
Full article ">Figure 3
<p>The typical section map of the Gongdang landslide.</p>
Full article ">Figure 4
<p>The geological structure background of the study area and the distribution of the geological hazards in the study area [<a href="#B41-water-16-03649" class="html-bibr">41</a>].</p>
Full article ">Figure 5
<p>The theory and calculation flow of the PFC.</p>
Full article ">Figure 6
<p>The PFC numerical model of the Gongdang landslide.</p>
Full article ">Figure 7
<p>The acquisition and calibration of the laboratory tests: (<b>a</b>) the specific gravity test; (<b>b</b>) the soil sample after direct shear test; (<b>c</b>) X-ray diffraction test.</p>
Full article ">Figure 8
<p>The schematic diagram of micro-parameters: (<b>a</b>) the numerical model used for the virtual direct shear test; (<b>b</b>) the shear strength results; (<b>c</b>) the stress–strain curves.</p>
Full article ">Figure 9
<p>The diagram of the average velocity and average displacement during the movement process: (<b>a</b>) the average velocity; (<b>b</b>) the average displacement.</p>
Full article ">Figure 10
<p>The diagram of the velocity field of different stages: (<b>a</b>) stage 1 (t = 20 s); (<b>b</b>) stage 2 (t = 60 s); (<b>c</b>) stage 3 (t = 300 s); (<b>d</b>) stage 4 (t = 600 s).</p>
Full article ">Figure 11
<p>The diagram of the displacement field of different stages: (<b>a</b>) stage 1 (t = 20 s); (<b>b</b>) stage 2 (t = 60 s); (<b>c</b>) stage 3 (t = 300 s); (<b>d</b>) stage 4 (t = 600 s).</p>
Full article ">Figure 12
<p>The locations of the monitoring points of the landslide.</p>
Full article ">Figure 13
<p>The monitoring result of the velocity: (<b>a</b>) all monitoring points; (<b>b</b>) rear edge; (<b>c</b>) middle part; (<b>d</b>) front edge.</p>
Full article ">Figure 14
<p>The monitoring results of the displacement: (<b>a</b>) all monitoring points; (<b>b</b>) rear edge; (<b>c</b>) middle part; (<b>d</b>) front edge.</p>
Full article ">Figure 15
<p>The locations of the cross and longitudinal lines.</p>
Full article ">Figure 16
<p>The diagram of the cross and longitudinal section maps: (<b>a</b>) the cross-section map; (<b>b</b>) the longitudinal section map.</p>
Full article ">Figure 17
<p>The simulation results with different parameters of the particles: (<b>a</b>) the maximum values of average displacement; (<b>b</b>) the maximum values of average velocity.</p>
Full article ">
26 pages, 19228 KiB  
Article
The Coupled Thermal Response Analysis of Green Roofs Based on the Discrete Element Method
by Chang Liu, Xiaoyong Zhang, Mingjie Jiang, Shengnan Zhu, Zhuan Wang and Jianxu Long
Buildings 2024, 14(12), 3984; https://doi.org/10.3390/buildings14123984 - 15 Dec 2024
Viewed by 680
Abstract
As an effective energy-saving measure, green roofs significantly improve the thermal environment of buildings by covering the roof with vegetation and soil. This paper compares the thermal transfer performance of concrete roofs and green roofs under different temperature conditions. First, a uniaxial compression [...] Read more.
As an effective energy-saving measure, green roofs significantly improve the thermal environment of buildings by covering the roof with vegetation and soil. This paper compares the thermal transfer performance of concrete roofs and green roofs under different temperature conditions. First, a uniaxial compression discrete element method (DEM) was used to calibrate the mesoscopic parameters of concrete, ensuring an accurate representation of concrete properties. The results indicate that green roofs have significant insulation effects under high-temperature conditions in summer. After being exposed to high temperatures for 5 h, the temperature of the green roof was 23.4 degrees Celsius lower than that of the ordinary concrete roof. In addition, different initial temperatures of the model also have a certain impact on heat transfer. The higher the initial temperature, the slower the temperature increase under high-temperature conditions. In winter, the green roof significantly delays the cooling at the top of the building, demonstrating excellent thermal insulation performance. The maximum temperature difference compared with the concrete roof is 8 °C. Finally, there is an exponential relationship between the thermal resistivity of the green roof and the temperature. In conclusion, green roofs have significant energy-saving and environmental protection value. Full article
Show Figures

Figure 1

Figure 1
<p>Diagram of green roof field test.</p>
Full article ">Figure 2
<p>Schematic diagram of particle heat transfer.</p>
Full article ">Figure 3
<p>The thermomechanical coupling model for the DEM.</p>
Full article ">Figure 4
<p>Simulation of uniaxial compression using discrete element method: (<b>a</b>) concrete sample, (<b>b</b>) schematic diagram of uniaxial compression discrete element.</p>
Full article ">Figure 5
<p>Stress–strain curves of uniaxial compression tests for different particle sizes.</p>
Full article ">Figure 6
<p>Simulation of stress–strain and crack number curves for concrete under uniaxial compression.</p>
Full article ">Figure 7
<p>Discrete element simulation results of concrete uniaxial compression test: (<b>a</b>) particle displacement cloud map, (<b>b</b>) particle displacement vector diagram, (<b>c</b>) crack propagation diagram.</p>
Full article ">Figure 8
<p>Concrete roof and green roof models: (<b>a</b>) structural schematic diagram, (<b>b</b>) discrete element model schematic diagram (unit: cm).</p>
Full article ">Figure 9
<p>Distribution of heat conduction pipes between the green roof and concrete.</p>
Full article ">Figure 10
<p>Temperature variation curve in summer.</p>
Full article ">Figure 11
<p>Discrete element model of high-temperature heat transfer between green roof and concrete roof.</p>
Full article ">Figure 12
<p>Cloud map of temperature changes at different times between green roof and concrete roof under high-temperature conditions.</p>
Full article ">Figure 13
<p>Temperature changes of heat pipes at different times for green roofs and concrete roofs under high-temperature conditions.</p>
Full article ">Figure 14
<p>The particle temperature changes with time at different locations under high-temperature conditions: (<b>a</b>) concrete roof, (<b>b</b>) green roof.</p>
Full article ">Figure 15
<p>Particle temperature curves of concrete and green roof at different positions at time t = 5 h under high-temperature conditions.</p>
Full article ">Figure 16
<p>The change curve of porosity and coordination number of concrete roof and green roof with time under high-temperature conditions.</p>
Full article ">Figure 17
<p>Cloud chart of particle temperature distribution with an initial temperature of 20 °C under high-temperature conditions: (<b>a</b>) concrete roof, (<b>b</b>) green roof.</p>
Full article ">Figure 18
<p>The variation of particle temperatures with different initial temperatures under high-temperature conditions over time: (<b>a</b>) concrete roof, (<b>b</b>) green roof.</p>
Full article ">Figure 19
<p>Discrete element model of low-temperature heat transfer between green roof and concrete roof.</p>
Full article ">Figure 20
<p>Cloud map of temperature changes at different times between green roof and concrete roof under low-temperature conditions.</p>
Full article ">Figure 21
<p>Temperature changes of heat pipes at different times for green roofs and concrete roofs under low-temperature conditions.</p>
Full article ">Figure 22
<p>The particle temperature changes with time at different locations under high-temperature conditions: (<b>a</b>) concrete roof, (<b>b</b>) green roof.</p>
Full article ">Figure 23
<p>Particle temperature curves of concrete and green roof at different positions at time t = 5 under low-temperature conditions.</p>
Full article ">Figure 24
<p>The change curve of porosity and coordination number of concrete roof and green roof with time under low-temperature conditions.</p>
Full article ">Figure 25
<p>Comparison of particle displacement cloud map of different green roof thermal resistance at t = 5 h.</p>
Full article ">Figure 26
<p>Variation curve of particle temperature at <span class="html-italic">Z</span> = 0.3 m of green roof with time under different thermal resistivity.</p>
Full article ">Figure 27
<p>Particle temperature curve of green roof at <span class="html-italic">Z</span> = 0.3 m under different thermal resistivity.</p>
Full article ">
19 pages, 7335 KiB  
Article
Mechanical Behavior of Marine Soft Soil with Different Water Contents Under Cyclic Loading
by Yajun Liu, Heng Zhang, Yindong Sun, Ke Wu and Wenbin Xiao
J. Mar. Sci. Eng. 2024, 12(12), 2307; https://doi.org/10.3390/jmse12122307 - 15 Dec 2024
Viewed by 415
Abstract
This study integrates macroscopic dynamic triaxial tests with microscopic discrete element simulations to comprehensively examine the dynamic deformation characteristics of marine soft soils under cyclic loading. Unlike previous research that typically focuses solely on experimental or numerical methods, this approach combines both techniques [...] Read more.
This study integrates macroscopic dynamic triaxial tests with microscopic discrete element simulations to comprehensively examine the dynamic deformation characteristics of marine soft soils under cyclic loading. Unlike previous research that typically focuses solely on experimental or numerical methods, this approach combines both techniques to enable a holistic analysis of soil behavior. The dynamic triaxial tests assessed macroscopic responses, including strain evolution and energy dissipation, under varying dynamic stress ratios, confining pressures, and water contents. Concurrently, discrete element simulations uncovered the microscopic mechanisms driving these behaviors, such as particle rearrangement, porosity variations, and shear zone development. The results show that (1) The strain range of marine soft soils increases significantly with higher dynamic stress ratios, confining pressures, and water contents; (2) Cumulative dynamic strain and particle displacement intensify at water contents of 50% and 55%. However, at a water content of 60%, the samples exhibit significant damage characterized by the formation of shear bands throughout the entire specimen; (3) As water content increases, energy dissipation in marine soft soils accelerates under lower confining pressures but increases more gradually under higher confining pressures. This behavior is attributed to enhanced particle packing and reduced pore space at elevated confining pressures. This integrated methodology not only enhances analytical capabilities but also provides valuable engineering insights into the dynamic response of marine soft soils. The findings offer essential guidance for the design and stabilization of marine soft soil infrastructure in coastal urban areas. Full article
Show Figures

Figure 1

Figure 1
<p>Triaxial apparatus.</p>
Full article ">Figure 2
<p>Numerical simulation model sample.</p>
Full article ">Figure 3
<p>The curve between strain and number of cycles obtained from triaxial specimen and discrete element analysis. (<b>a</b>) Test results. (<b>b</b>) Discrete element analysis results.</p>
Full article ">Figure 4
<p>Stress–strain curves of the first 10 cycles under different moisture content conditions. (<b>a</b>) σ<sub>c</sub> = 50 kPa, <span class="html-italic">η</span><sub>d</sub> = 0.15. (<b>b</b>) σ<sub>c</sub> = 100 kPa, <span class="html-italic">η</span><sub>d</sub> = 0.15. (<b>c</b>) σ<sub>c</sub> = 100 kPa, <span class="html-italic">η</span><sub>d</sub> = 0.25. (<b>d</b>) σ<sub>c</sub> = 200 kPa, <span class="html-italic">η</span><sub>d</sub> = 0.15.</p>
Full article ">Figure 5
<p>Stress–strain curves of the first 10 cycles under different dynamic stress ratios. (<b>a</b>) σ<sub>c</sub> = 50 kPa, <span class="html-italic">w</span> = 40%. (<b>b</b>) σ<sub>c</sub> = 100 kPa, <span class="html-italic">w</span> = 45%. (<b>c</b>) σ<sub>c</sub> = 100 kPa, <span class="html-italic">w</span> = 50%. (<b>d</b>) σ<sub>c</sub> = 200 kPa, <span class="html-italic">w</span> = 60%.</p>
Full article ">Figure 6
<p>Stress–strain curves of the first 10 cycles under different confining pressures. (<b>a</b>) <span class="html-italic">w</span> = 40%, <span class="html-italic">η</span><sub>d</sub> = 0.15. (<b>b</b>) <span class="html-italic">w</span> = 45%, <span class="html-italic">η</span><sub>d</sub> = 0.15. (<b>c</b>) <span class="html-italic">w</span> = 50%, <span class="html-italic">η</span><sub>d</sub> = 0.25. (<b>d</b>) <span class="html-italic">w</span> = 45%, <span class="html-italic">η</span><sub>d</sub> = 0.15.</p>
Full article ">Figure 7
<p>Schematic diagram of energy dissipation in marine soft soil under a dynamic stress ratio of 0.15.</p>
Full article ">Figure 8
<p>Schematic diagram of energy dissipation in marine soft soil at 100 kPa perimeter pressure condition.</p>
Full article ">Figure 9
<p>Particle displacement distribution of samples with a moisture content of 40.</p>
Full article ">Figure 10
<p>Particle displacement distribution of samples with a moisture content of 50.</p>
Full article ">Figure 11
<p>Particle displacement distribution of samples with a moisture content of 60.</p>
Full article ">Figure 12
<p>Deformation curve of porosity with cycle times under different water content conditions. (<b>a</b>) σ<sub>c</sub> = 50 kPa, <span class="html-italic">η</span><sub>d</sub> = 0.15. (<b>b</b>) σ<sub>c</sub> = 100 kPa, <span class="html-italic">η</span><sub>d</sub> = 0.15. (<b>c</b>) σ<sub>c</sub> = 100 kPa, <span class="html-italic">η</span><sub>d</sub> = 0.25. (<b>d</b>) σ<sub>c</sub> = 200 kPa, <span class="html-italic">η</span><sub>d</sub> = 0.15.</p>
Full article ">Figure 13
<p>Deformation curve of porosity with number of cycles under different dynamic stress ratios. (<b>a</b>) <span class="html-italic">w</span> = 40%, σ<sub>c</sub> = 100 kPa. (<b>b</b>) <span class="html-italic">w</span> = 45%, σ<sub>c</sub> = 100 kPa. (<b>c</b>) <span class="html-italic">w</span> = 50%, σ<sub>c</sub> = 100 kPa. (<b>d</b>) <span class="html-italic">w</span> = 60%, σ<sub>c</sub> = 100 kPa.</p>
Full article ">Figure 14
<p>Deformation curve of porosity with number of cycles. (<b>a</b>) <span class="html-italic">w</span> = 40%, <span class="html-italic">η</span><sub>d</sub> = 0.15. (<b>b</b>) <span class="html-italic">w</span> = 45%, <span class="html-italic">η</span><sub>d</sub> = 0.15. (<b>c</b>) <span class="html-italic">w</span> = 50%, <span class="html-italic">η</span><sub>d</sub> = 0.15. (<b>d</b>) <span class="html-italic">w</span> = 60%, <span class="html-italic">η</span><sub>d</sub> = 0.15.</p>
Full article ">
20 pages, 5882 KiB  
Article
Contact Parameter Calibration for Discrete Element Potato Minituber Seed Simulation
by Kai Chen, Xiang Yin, Wenpeng Ma, Chengqian Jin and Yangyang Liao
Agriculture 2024, 14(12), 2298; https://doi.org/10.3390/agriculture14122298 - 14 Dec 2024
Viewed by 580
Abstract
The discrete element method (DEM) has been widely applied as a vital auxiliary technique in the design and optimization processes of agricultural equipment, especially for simulating the behavior of granular materials. In this study, the focus is placed on accurately calibrating the simulation [...] Read more.
The discrete element method (DEM) has been widely applied as a vital auxiliary technique in the design and optimization processes of agricultural equipment, especially for simulating the behavior of granular materials. In this study, the focus is placed on accurately calibrating the simulation contact parameters necessary for the V7 potato minituber seed DEM simulation. Firstly, three mechanical tests are conducted, and through a combination of actual tests and simulation tests, the collision recovery coefficient between the seed and rubber material is determined to be 0.469, the static friction coefficient is 0.474, and the rolling friction coefficient is 0.0062. Subsequently, two repose angle tests are carried out by employing the box side plates lifting method and the cylinder lifting method. With the application of the response surface method and a search algorithm based on Matlab 2019, the optimal combination of seed-to-seed contact parameters, namely, the collision recovery coefficient, static friction coefficient, and rolling friction coefficient, is obtained, which are 0.500, 0.476, and 0.043, respectively. Finally, the calibration results are verified by a seed-falling device that combines collisions and accumulation, and it is shown that the relative error between the simulation result and the actual result in the verification test is small. Thus, the calibration results can provide assistance for the design and optimization of the potato minituber seed planter. Full article
(This article belongs to the Section Agricultural Technology)
Show Figures

Figure 1

Figure 1
<p>Potato minituber seed triaxial size.</p>
Full article ">Figure 2
<p>Miniature potato compression test.</p>
Full article ">Figure 3
<p>Potato minituber seed simulation model. (<b>a</b>) ellipsoidal model, (<b>b</b>) spherical model.</p>
Full article ">Figure 4
<p>Potato minituber seed—rubber plate free fall test. (<b>a</b>) actual test, (<b>b</b>) simulation test.</p>
Full article ">Figure 5
<p>Collision recovery coefficient fits the curve.</p>
Full article ">Figure 6
<p>Slope sliding test. (<b>a</b>) actual test, (<b>b</b>) simulation test.</p>
Full article ">Figure 7
<p>Static friction coefficient fitting curve.</p>
Full article ">Figure 8
<p>Inclined plane rolling test. (<b>a</b>) actual test, (<b>b</b>) simulation test.</p>
Full article ">Figure 9
<p>Rolling friction coefficient fitting curve.</p>
Full article ">Figure 10
<p>Repose angle test device. (<b>a</b>) box side plates lifting method, (<b>b</b>) cylinder lifting method.</p>
Full article ">Figure 11
<p>Simulation test of box side plates lifting method. (<b>a</b>) initial state, (<b>b</b>) side plates lifting, (<b>c</b>) test completed.</p>
Full article ">Figure 12
<p>Simulation test of cylinder lifting method. (<b>a</b>) initial state, (<b>b</b>) cylinder lifting, (<b>c</b>) test completed.</p>
Full article ">Figure 13
<p>Repose angle image processing and angle acquisition. (<b>a</b>) initial image, (<b>b</b>) binarization, (<b>c</b>) boundary extraction, (<b>d</b>) linear fitting.</p>
Full article ">Figure 14
<p>3D response surface diagram of interaction term BC. (<b>a</b>) interaction terms AB, (<b>b</b>) interaction terms AC, (<b>c</b>) interaction terms BC.</p>
Full article ">Figure 15
<p>Test verification device physical drawing.</p>
Full article ">Figure 16
<p>Simulation seed drop test. (<b>a</b>) seeds began to pile up, (<b>b</b>) seed pile is almost complete.</p>
Full article ">
22 pages, 4715 KiB  
Article
Design and Testing of a Low-Speed, High-Frequency Straw Chopping and Returning Machine Using a Constant Breath Cam Mechanism
by Han Lin, Jin He, Guangyuan Zhong and Hanyu Yang
Agriculture 2024, 14(12), 2293; https://doi.org/10.3390/agriculture14122293 - 14 Dec 2024
Viewed by 314
Abstract
Straw incorporation offers significant advantages in agricultural crop cultivation systems. Mechanized methods constitute the predominant approach, potentially reducing yield costs and enhancing operational efficiency. The imperative to enhance the quality of straw chopping within the field is of particular significance, as suboptimal chopping [...] Read more.
Straw incorporation offers significant advantages in agricultural crop cultivation systems. Mechanized methods constitute the predominant approach, potentially reducing yield costs and enhancing operational efficiency. The imperative to enhance the quality of straw chopping within the field is of particular significance, as suboptimal chopping quality can engender a cascade of issues, particularly seeding blockages. The straw chopping pass rate (CPR) is a pivotal metric for assessing the quality of straw chopping. Therefore, enhancing the CPR during the straw chopping process is necessary. This study introduces a novel maize-straw-chopping device with the ground as its supporting base. This device facilitates the rapid vertical chopping of maize straw through a constant breath cam transmission mechanism. Critical parameters were determined to optimize the performance of the chopping device by establishing mathematical models and kinematic simulation analysis methods. With the help of Rocky 2022.R2 software, the influence of the rotational velocity of the draft, tractor velocity, and blade edge angles on the CPR during the operation of the device was analyzed. The Box–Behnken test methodology was used to carry out a three-factor, three-level orthogonal rotation test to obtain the optimal working parameter combination. The results indicated that the maximum CPR value was achieved with a draft rotational velocity of 245 rpm, a tractor velocity of 3.8 km/h, and a blade edge angle of 20.75°. Finally, field validation experiments were conducted under these optimized conditions, with the average CPR of maize straw reaching an impressive 91.45%. These findings have significant implications for enhancing crop production practices. Full article
(This article belongs to the Section Agricultural Technology)
Show Figures

Figure 1

Figure 1
<p>Overall structure of straw-chopping and -returning machine: (1) sliders; (2) blade; (3) plates; (4) cam; (5) drive wheels; (6) shaft; (7) sprocket.</p>
Full article ">Figure 2
<p>Low-speed and high-frequency working principle. The yellow represents maize straw; The blue represents plates; The green represents cam; The red represents blade.</p>
Full article ">Figure 3
<p>Theoretical contour lines of cams with different parameters “<span class="html-italic">n</span>”.</p>
Full article ">Figure 4
<p>Schematic diagram of cam mechanism movement. <span class="html-italic">Oxy</span> is a stationary coordinate system; <span class="html-italic">Ox</span><sub>1</sub><span class="html-italic">y</span><sub>1</sub> is a dynamic coordinate system; AB is the straight line where <span class="html-italic">oy</span> is located; A<sub>1</sub>B<sub>1</sub> is the straight line where <span class="html-italic">oy</span><sub>1</sub> is located; <span class="html-italic">O</span><sub>1</sub><span class="html-italic">x</span><sub>2</sub><span class="html-italic">y</span><sub>2</sub> is a moving coordinate system; <span class="html-italic">r</span><sub>0</sub> is the radius of the roller; <span class="html-italic">δ</span> is the arbitrary angular displacement through which the roller rotates; <span class="html-italic">θ</span> is the arbitrary angular displacement through which the roller rotates; <span class="html-italic">ω</span> is the angular velocity of cam rotation; <span class="html-italic">T</span> represents the coordinates of the contact point between the roller and the cam.</p>
Full article ">Figure 5
<p>Variation law of the pressure angle under different factors.</p>
Full article ">Figure 6
<p>Mechanical model of the straw chopping process.</p>
Full article ">Figure 7
<p>Description of a bonded sphero-cylinder model.</p>
Full article ">Figure 8
<p>The simulation process during operation (<b>a</b>) Front view; (<b>b</b>) top view. (1) Stubble; (2) straw; (3) chopping device; (4) blade; (5) broken straw.</p>
Full article ">Figure 8 Cont.
<p>The simulation process during operation (<b>a</b>) Front view; (<b>b</b>) top view. (1) Stubble; (2) straw; (3) chopping device; (4) blade; (5) broken straw.</p>
Full article ">Figure 9
<p>Analysis of the interaction between two factors on CPR. (<b>a</b>) Influence of factor A and factor B on CPR; (<b>b</b>) influence of factor B and factor C on CPR. Different colors represent different straw CPR.</p>
Full article ">Figure 10
<p>Field validation experiment.</p>
Full article ">
Back to TopTop