Computational Particle Mechanics最新文献

Ejector deep hole drilling achieves high-quality boreholes in production processes. High feed rates are applied to ensure a high productivity level, requiring reliable chip removal from the cutting zone for a stable process. Therefore, a constant metalworking fluid flow under high volume flow rates or high pressure is required. Experimental results show a vortex formation at the outer cutting edge. This vortex can lead to delayed chip removal from the cutting zone, and ultimately, it can lead to chip clogging and result in drill breakage due to increased torque. This paper investigates modified drill head designs using the smoothed particle hydrodynamics method. The investigated modifications include various designs of the chip mouth covering. Besides graphical analysis based on flow visualizations, flow meters are placed at the tool’s head to evaluate the impact of the modifications on the flow rate and possible increased resistance and relocation of the fluid flow from the outer cutting edge to other parts of the tool. The simulation results for the reference design show the experimentally observed vortex formation, validating the simulation model. By adding the tool’s rotation in the SPH simulation, which is not included in the experiments for observation reasons, the vortex formation is positively influenced. In addition, some designs show promising results to further mitigate the vortex formation while maintaining a sufficient fluid flow around the cutting edges.

Virtual crime scene investigation using numerical models has the potential to assist in the forensic investigation of firearm-related fatalities, where ethical concerns and expensive resources limit the scope of physical experiments to comprehend the post-impact biomechanics comprehensively. The human cranial numerical model developed in this study incorporates three main components (skin, skull, and brain) with dynamic biomaterial properties. The virtual model provides valuable insights into the post-impact biomechanics of cranial ballistic injuries, particularly in high-speed events beyond conventional investigative capabilities, including the velocity of ejected blood backspatter, cavitation collapsing, and pressure waves. The validation of the numerical model, both quantitatively and qualitatively, demonstrates its ability to replicate similar bone fractures, entrance wound shapes, and backward skin ballooning observed in physical experiments of the human cranial geometry. The model also yields similar temporary cavity sizes, wound sizes, and blood backspatter time against the physical cranial model, aiding in bloodstain pattern analysis. Additionally, the numerical model enables exploration of ballistic factors that vary in each crime scene environment and influence cranial injuries, such as projectile type, velocity, impact location, and impact angle. These established injury patterns contribute to crime scene reconstruction by providing essential information on projectile trajectory, discharge distance, and firearm type, assisting in the resolution of court cases. In conclusion, the developed human cranial geometry in this study offers a reliable tool for investigating firearm-related cranial injuries, serving as a statistical reference in forensic science. Virtual crime scene investigationsusing these modelshave the potential to enhance the accuracy and efficiency of forensic analyses.

When metals and alloys are exposed to ultrasonic vibrations (UV), a softening behavior occurs, caused by the phenomenon of acousto-plasticity. To obtain accurate results in a deformation analysis, this phenomenon must be included in the formulation of the constitutive material model. In this work, an acoustic-plastic model is proposed to capture the effects of ultrasonic vibrations during machining. The desired effect is to modify the chip morphology to reduce the magnitude of the cutting forces and thus reduce the energy consumption of the process. The study focuses on the modeling of ultrasonic vibration-assisted micromachining (VAMM). The particle finite element method is used and extended to perform a thermo-mechanical analysis capable of capturing the responses of conventional micromachining (CMM) and VAMM operations of 32 HRC stainless steel. The cutting speed and UV parameters, including amplitude and frequency, are integrated into the Johnson–Cook constitutive model to account for the effects of acoustic softening on the machining characteristics. The results show that the influence of UV on microcutting leads to thinner chips and lower cutting force. In the VAMM operations, an average reduction in cutting forces of 20% is achieved at five different cutting speeds. In addition, the contact length between the tool and chip decreases at different cutting speeds from 29% to a maximum of 44%. Furthermore, the thermal analysis results show that there is a negligible temperature change during the CMM and VAMM simulations, indicating that the study of the machining process can focus exclusively on its mechanical aspects when performed at the microscale. The predicted average chip thickness and effective shear angle of the workpiece material are in strong agreement with the experimental results, emphasizing the importance of considering acoustic softening in VAMM studies.

To determine the desirable bonding parameters of the soybean-bonded particle model for accidentally simulating the working process of a pneumatic soybean seed-metering device. Taking the compressive destructive force (F_c,p) derived from the uniaxial compression test of soybean seeds as the evaluation index for the compression simulation tests. The Plackett–Burman and the steepest ascent tests were executed to identify the centroids of the influential factors that substantially affect the bonding force of the soybean-bonded particle model. The optimal values of the significance influencing variables were determined based on the Box–Behnken response surface test. The results indicated that the effect of bonded disk radius (R_B,p) between fraction particles on the F_c,p was extremely significant, and the effects of the restitution coefficient (e_p-steel) and static friction coefficient (μ_p-steel) of soybean-steel, normal stiffness per unit area (k_n,p) and critical normal stress (σ_max,p) were found to be statistically significant. The preferred values identified by Box–Behnken response surface test were 0.520 for e_p-steel, 0.274 for μ_p-steel, 4.082 × 10⁷ N/m³ for k_n,p, 3.517 × 10⁵ Pa for σ_max,p, and 0.982 mm for R_B,p, respectively. The compressive destructive force of soybean seeds was 211.32 N at this point, which was 0.2% less than the measured value of 211.74 N. The results of comparing the grain morphologies during the actual and simulated compressions indicated that the compression states had a superior consistency. It was determined that the DEM simulation input parameters for the soybean-bonded particle model calibrated were proven to be effective and dependable. The investigation presented in this paper can be utilized to effectively analyze the working process of the pneumatic soybean seed-metering devices through coupled simulation. It can also serve as a reference for other researchers to construct a particle model for DEM simulation using the BPM approach.

This paper presents a novel method for smoothed particle hydrodynamics (SPH) with thin-walled rigid structures. Inspired by the direct forcing immersed boundary method, this method employs a moving least square method for the velocity interpolation instead of the linear interpolation. It reduces oscillations due to changing relative positions between fluid grids and structures. It also simplifies thin-walled rigid structure simulations by eliminating the need for multiple layers of boundary particles, and improves computational accuracy and stability in three-dimensional scenarios. Results of the impulsively started plate test demonstrate that the proposed method obtains smooth velocity and pressure, as well as a good match to the references results of the vortex wake development. Results of the flow past cylinder test show that the proposed method avoids mutual interference on both side of the boundary, while accurately calculating the forces acting on structure. By comparing to linear least square direct forcing scheme and the diffusive direction scheme, advantages of lower oscillation and higher accuracy are proven. Results of flow past a sphere further indicate the stability of the proposed method for three-dimensional simulations.

Existing literature on true triaxial and torsional shear tests indicate that the mechanical response of a granular assembly is significantly influenced by the magnitude of the intermediate principal stress ratio. The present study aims to explore the mechanism behind such effects in reference to the particle-level interaction using 3D DEM simulations. In this regard, true triaxial numerical simulations have been carried out with constant minor principal stress and varying (b) values employing rolling resistance-type contact model to mimic particle shape. The numerical simulations have been validated against the true triaxial experiments reported in the literature for dense Santa Monica beach sand. The macro-level shearing response of the granular assembly has been examined in terms of the evolution of stress ratio and volumetric strain for different rolling resistance coefficients. Further, such macro-level response has been assessed in reference to the micro-scale attributes, e.g. average contact force, number of interparticle contacts, mechanical coordination number, contact normal orientation, and fabric tensor as well as meso-scale attribute like strong contact force network. Lade’s failure surface has been adopted to represent the stress and fabric at peak state in the octahedral plane, and mathematical expressions have been proposed relating the failure surface parameters to the rolling resistance coefficient.

Graphical abstract

Discrete element method (DEM) is a powerful tool for the dynamic simulation of irregular non-spherical particle systems. The efficient integration of the rotational motions of numerous particles in DEM poses a big challenge. This paper presents six explicit time integration algorithms, comprising three first-order algorithms and three second-order algorithms, for the rotational motions of non-spherical particles based on the theory of unit quaternion group S(3). The proposed algorithms based on Cayley map do not contain any transcendental function and have high efficiency. The numerical examples underscore the superiority of the first-order symplectic Euler Cayley algorithm (SECay) and the second-order central difference Cayley algorithm (CDCay) in terms of both efficiency and accuracy. In the testing cases of granular systems, SECay and CDCay demonstrate approximately 80% reduction in computational time for the time integration part, compared to the improved predictor–corrector direct multiplication method (IPCDM). Therefore, SECay and CDCay emerge as promising tools for non-spherical DEM simulations.

Motivated by the spatial variability observed in geological profiles, this study explored the feasibility of using discrete element method (DEM) to capture the effect of layered spatial variability into overall soil performance. The spatial variability of packing densities, particle Young’s modulus (E), and frictional properties (μ) within specimens was studied. It was observed that samples with similar overall void ratios exhibited comparable small-strain stiffness and shearing behaviors. In contrast, the coordination number and particle stress transmission demonstrated significant sensitivity to the layer-wise spatial variability in packing densities. Regarding the spatial variability effect of particle-scale E values, this study illustrates that spatial variability strongly affects the stiffness contributions of individual layers. Specifically, layers with higher E values are capable of transferring much stress and stiffness. For the spatial variability effect of frictional property, a degree of consistency in shearing behaviors was observed among specimens with similar average frictional characteristics, while layers with lower frictional property were identified as potential initial failure junctures. Overall, this study validates the utility of employing a DEM code for analyzing both the macroscopic behavior and localized vulnerabilities within complex granular systems, presenting profound implications for engineering practices.