Mechanical Behavior Of Materials Research Articles

Constitutive Model for Thermal-Oxygen-Aged EPDM Rubber Based on the Arrhenius Law.

Ethylene-propylene-diene monomer (EPDM) is a key engineering material; its mechanical characterization is important for the safe use of the material. In this paper, the coupled effects of thermal degradation temperature and time on the tensile mechanical behavior of EPDM rubber were investigated. The tensile stress-strain curves of the aged and unaged EPDM rubber show strong nonlinearity, demonstrating especially rapid stiffening as the strain increases under small deformation. The popular Mooney-Rivlin and Ogden (N = 3) models were chosen to fit the test data, and the results indicate that neither of the classical models can accurately describe the tensile mechanical behavior of this rubber. Six hyperelastic constitutive models, which are excellent for rubber with highly nonlinearity, were employed, and their abilities to reproduce the stress-strain curve of the unaged EPDM were assessed. Finally, the Davis-De-Thomas model was found to be an appropriate hyperelastic model for EPDM rubber. A Dakin-type kinetic relationship was employed to describe the relationships between the model parameters and aging temperature and time, and, combined with the Arrhenius law, a thermal aging constitutive model for EPDM rubber was established. The ability of the proposed model was checked by independent testing data. In the moderate strain range of 200%, the errors remained below 10%. The maximum errors of the prediction results at 85 °C for 4 days and 100 °C for 2 and 4 days were computed to be 17.06%, 17.51% and 19.77%, respectively. This work develops a theoretical approach to predicting the mechanical behavior of rubber material that has suffered thermal aging; this approach is helpful in determining the safe long-term use of the material.

Impact of Regular and Irregular Pore Distributions on the Elasticity of Porous Materials: A Microstructure-Free Finite Element Study.

Conventional analytical formulas for predicting the effective Young's modulus of porous materials often rely on simplifying assumptions and do not explicitly incorporate microstructural information. This study investigates the impact of regular versus irregular pore distributions on the stiffness of porous materials using microstructure-free finite element modeling (MF-FEM). After conducting a convergence study, MF-FEM predictions were validated against experimental data and used to assess the accuracy of commonly employed analytical models. The results demonstrate that materials with irregular microstructures exhibit a rapid decrease in Young's modulus, approaching zero at porosities slightly greater than 50%. In contrast, regular microstructures show a more gradual decline, maintaining significant stiffness until the porosity exceeds 90%. Additionally, the study reveals that some analytical formulas align better with irregular microstructures while others are more suited to regular ones, attributable to the underlying assumptions of these models. These findings underscore the necessity of considering pore distribution patterns in modeling to accurately predict the mechanical behavior of porous materials.

Effect of helium bubbles on the mobility of edge dislocations in copper

Abstract Helium bubbles can form in materials upon exposure to irradiation. It is well known that the presence of helium bubbles can cause changes in the mechanical behavior of materials. To improve the lifetime of nuclear components, it is important to understand deformation mechanisms in helium-containing materials. In this work, we investigate the interactions between edge dislocations and helium bubbles in copper using molecular dynamics (MD) simulations. We focus on the effect of helium bubble pressure (equivalently, the helium-to-vacancy ratio) on the obstacle strength of helium bubbles and their interaction with dislocations. Our simulations predict significant differences in the interaction mechanisms as a function of helium bubble pressure. Specifically, bubbles with high internal pressure are found to exhibit weaker obstacle strength as compared to low-pressure bubbles of the same size due to the formation of super-jogs in the dislocation. Activation energies and rate constants extracted from the MD data confirm this transition in mechanism and enable upscaling of these phenomena to higher length-scale models.

Multiscale modelling and characterisation of fused filament fabricated neat and graphene nanoplatelet reinforced G-polymers

AbstractIn this study, the fused filament fabrication of novel biodegradable G-polymer is investigated, including also filaments reinforced by graphene nanoplatelets. A combined framework of computational multiscale modelling and experimental studies is applied to characterise and model the mechanical behaviour of composite and nanocomposite materials, fabricated using a 3D-printer by depositing G-polymer filaments in +45$$^\circ$$ ∘ and $$-45^\circ$$ - 45 ∘ angles. The investigation is performed on both hot-extruded single filaments and 3D-printed specimens. Experimental results of hot-extruded single filaments are utilized on the microscale to extract the macroscale constitutive models. At the microscale level, representative volume elements with different percentages of penetration are analyzed to compute the effective orthotropic elastic material properties. A comprehensive comparison study is conducted using different micromechanical models, multiscale simulation outputs, and the mechanical properties resulting from experiments. The accuracy of the results obtained from the homogenization technique is validated against realistic finite element microstructural models and experimental measurements. The results demonstrate that computational homogenization, when coupled with accurate property characterisation, serves as a reliable tool for predicting the elastic response of 3D-printed parts. Furthermore, the proposed framework reduces the need for extensive experimental replication and lowers manufacturing costs.

Revamp of the eco-friendly grouting materials by mixing basalt fiber for earthen sites conservation: Compressive performance and constitutive models

Reinforcement techniques, including anchorage and sealing with grout, are frequently employed to protect earthen sites and minimize the impact of fissures, the effectiveness of these approaches relies heavily on the mechanical properties of grout. To enhance the strength and ductility of the conventional grout, an eco-friendly composite grouting material mixture basalt fibers of varying lengths and contents is prepared. The DIC-based UCS tests were conducted to evaluate its compressive performance. Constitutive models of the basalt fiber-reinforced grout, including elastic perfectly-plastic (EPP) and median nonlinear (MNL) constitutive models, were established respectively based on the equal energy rule and cross-sectional data analysis for simplification application and precision evaluation purposes. The results show that the mixture of basalt fibers resulted in a distinct semi-brittle failure mode for the grout, with a significant increase in its compressive strength, ductility, compressive fracture energy, and residual bearing capacity. However, the elastic modulus will be lower when the fibers are relatively longer. Considering the differences in the relations between performance parameters and fiber content, and length, the highest compressive strength was reached at the fiber length of 6 mm and the content of 0.2 wt%, while maximum ductility was achieved at 18 mm and 0.4 wt%. The fracture energy exhibits an approximate positive correlation with fiber content. In addition, the EPP and MNL constitutive models can effectively characterize the compressive mechanical behavior of semi-brittle materials to satisfy different requirements on the precision of engineering analysis. This study contributes valuable insights into the development of high-ductility grouting materials of earthen sites, and the establishment of constitutive models for semi-brittle materials.

Atomistic insight into laser shock peening response of α-titanium: an experimentally verified simulation study

Laser shock peening (LSP) can induce severe plastic deformation in the surface zone of metallic materials, which potentially enhances their mechanical properties. Although it is widely recognized that the improvement of mechanical properties can be attributed to microstructure evolution, the microstructural characteristics and response under different LSP parameters are complex, leading to distinct mechanical behaviors of materials at the macroscale. Besides, the exceedingly brief duration of LSP also poses considerable challenges to capture the microstructure evolution within the surface layer of materials. Therefore, to grasp a first in-situ atomistic insight into LSP response of α-titanium, molecular dynamics method is applied to simulate laser incidence at the [101¯0] crystallographic orientation with various energy densities. The results show that the lattice reorientation and dislocation slip dominate the plastic deformation at low laser energy density, while stacking faults and crystal disordering are the primary microstructural features at high laser energy density. The above shock-induced microstructures are verified by transmission electron microscopy (TEM) observations. Further, varied microscale plastic deformation modes under varying shock energies are proposed based on Schmid factors and energy barriers. The findings bridge laser energy to post-processing microstructure of α-titanium, which paves the way for improving mechanical properties of materials by optimizing LSP treatment parameters.

A new thermo-optical system with a fractional Caputo operator for a rotating spherical semiconductor medium immersed in a magnetic field

PurposeUnderstanding the mechanical and thermal behavior of materials is the goal of the branch of study known as fractional thermoelasticity, which blends fractional calculus with thermoelasticity. It accounts for the fact that heat transfer and deformation are non-local processes that depend on long-term memory. The sphere is free of external stresses and rotates around one of its radial axes at a constant rate. The coupled system equations are solved using the Laplace transform. The outcomes showed that the viscoelastic deformation and thermal stresses increased with the value of the fractional order coefficients.Design/methodology/approachThe results obtained are considered good because they indicate that the approach or model under examination shows robust performance and produces accurate or reliable results that are consistent with the corresponding literature.FindingsThis study introduces a proposed viscoelastic photoelastic heat transfer model based on the Moore-Gibson-Thompson framework, accompanied by the incorporation of a new fractional derivative operator. In deriving this model, the recently proposed Caputo proportional fractional derivative was considered. This work also sheds light on how thermoelastic materials transfer light energy and how plasmas interact with viscoelasticity. The derived model was used to consider the behavior of a solid semiconductor sphere immersed in a magnetic field and subjected to a sudden change in temperature.Originality/valueThis study introduces a proposed viscoelastic photoelastic heat transfer model based on the Moore-Gibson-Thompson framework, accompanied by the incorporation of a new fractional derivative operator. In deriving this model, the recently proposed Caputo proportional fractional derivative was considered. This work also sheds light on how thermoelastic materials transfer light energy and how plasmas interact with viscoelasticity. The derived model was used to consider the behavior of a solid semiconductor sphere immersed in a magnetic field and subjected to a sudden change in temperature.

Exploring grain size influence on tensile behavior of 316 H austenitic stainless steel at high temperature: A phenomenological dislocation model

The high temperature mechanical behavior of 316 H austenitic stainless steel with different grain size was investigated in this study through the use of a novel phenomenological dislocation model finite element method. The study was conducted by producing steel specimens with varying grain sizes through adjustments in annealing temperature and time, followed by high temperature tensile tests at 550°C. The finite element method, when integrated with a phenomenological dislocation model that considers grain size parameters, effectively captured the mechanical response under different grain size conditions. Additionally, the two-variable Kocks-Mecking(KM) model accurately captured the grain size effect in 316 H austenitic steel and effectively depicted its tensile flow behavior. The model’s predictions were based on the evolution of forest dislocation density and mobile dislocation density, proving to be a reliable tool for analyzing the microstructural evolution of 316 H austenitic steel specimens with varying grain sizes. This study provides insight into the effect of grain size on the high temperature strength of austenitic stainless steels and demonstrates the utility of a novel phenomenological model finite element method for predicting the mechanical behavior of polycrystalline materials.

Evaluation of the Mechanical Properties of PLA Material Used For 3D Printing Solar E-Hub Component

The advent of additive manufacturing (AM) technologies revolutionized design and production processes across various industries. In renewable energy, AM enabled new possibilities for optimizing solar e-hub (solar energy harvesting module) configurations, enhancing efficiency and performance. This study examined critical considerations such as material selection, durability, and cost-effectiveness in solar hub development. This fast-prototyping technique was controlled by computer-aided design (CAD) software like CREO Parametric 7.0 and Creality Slicer 4.8. Experimental results indicated that PLA (Polylactic acid) materials exhibited superior strength, with an impact energy of 4.8 Joules. The deformation study revealed that the maximum load of 22 MPa aligned with the ultimate tensile strength of PLA, and a hardness test result of 83.1 HRF featured its exemplary hardness properties. These findings advanced the understanding of using AM to investigate mechanical behaviours of PLA materials and optimize solar e-hub configurations for portable device applications, promoting sustainable energy solutions and the adoption of renewable energy technologies. In addition, the successful implementation of this approach will enable the renewable energy sectors to minimize the carbon foot-prints towards helping the global net-zero emissions by aligning the circular economy approach.

Enhancing shear strength predictions of rocks using a hierarchical ensemble model

Shear strength (SS) parameters are essential for understanding the mechanical behavior of materials, particularly in geotechnical engineering and rock mechanics. This study proposes a novel hierarchical ensemble model (HEM) to predict SS parameters: cohesion (C\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$C$$\\end{document}) and angle of internal friction (φ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varphi$$\\end{document}). The HEM addresses the limitations of traditional machine learning models. Its performance was validated using leave-one-out cross-validation (LOOCV) and out-of-bag (OOB) evaluation methods. The model's accuracy was assessed with R-squared correlation (R2), absolute average relative error percentage (AAREP), Taylor diagrams, and quantile–quantile plots. The computational results demonstrated that the proposed HEM outperforms previous studies using the same database. The model predicted φ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varphi$$\\end{document} and C\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$C$$\\end{document} with R2 values of 0.93 and 0.979, respectively. The AAREP values were 1.96% for φ and 4.7% for C\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$C$$\\end{document}. These results indicate that the HEM significantly improves the prediction quality of φ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varphi$$\\end{document} and C\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$C$$\\end{document}, and exhibits strong generalization capability. Sensitivity analysis revealed that σ_3maxσ3max (maximum principal stress) had the greatest impact on modeling both φ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varphi$$\\end{document} and C\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$C$$\\end{document}. According to uncertainty analysis, the LOOCV and OOB had the widest uncertainty bands for the φ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varphi$$\\end{document} and C\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$C$$\\end{document} parameters, respectively.

Experimental and Numerical Study on the Blast Resistant Performance of Geopolymer Concrete

Geopolymer, with its notable benefit of low carbon dioxide emissions, holds the potential to substantially curtail environmental pollution. According to the existing related research on geopolymer materials, it is obvious that it has great development potential in many engineering application fields, and it is a new generation of green and environmentally friendly recycled materials. Nowadays, there is a growing concern regarding explosion protection. Explosions near buildings can cause catastrophic damages on the building external and internal structure, and the most important thing is that can cause injuries and loss of life to the occupants of these buildings. This study investigates the mechanical performance of the fiber-reinforced geopolymer concrete under explosive testing. Furthermore, the finite element analysis models have been established through LS-DYNA software to simulate the explosive testing using Structure-Arbitrary Lagrangian Eulerian Method (S-ALE). The model is used to assess the dynamic mechanical behavior of geopolymer materials.

Influence of Substrate Location on Mechanical Behaviour of Glass Fibre Composite Materials with Embedded Printed Electronics

Influence of Substrate Location on Mechanical Behaviour of Glass Fibre Composite Materials with Embedded Printed Electronics

Prediction of Mechanical Properties of Nano-Clay-Based Biopolymeric Composites.

An understanding of the mechanical behavior of polymeric materials is crucial for making advancements in the applications and efficiency of nanocomposites, and encompasses their service life, load resistance, and overall reliability. The present study focused on the prediction of the mechanical behavior of biopolymeric nanocomposites with nano-clays as the nanoadditives, using a new modeling and simulation method based on Comsol Multiphysics software 6.1. This modeling considered the complex case of flake-shaped nano-clay additives that could form aggregates along the polymeric matrix, varying the nanoadditive thickness, and consequently affecting the resulting mechanical properties of the polymeric nanocomposite. The polymeric matrix investigated was biopolyamide 11 (BIOPA11). Several BIOPA11 samples reinforced with three different contents of nano-clays (0, 3, and 10 wt%), and with three different nano-clay dispersion grades (employing three different extrusion screw configurations) were obtained by the compounding extrusion process. The mechanical behavior of these samples was studied by the experimental tensile test. The experimental results indicate an enhancement of Young's modulus as the nano-clay content was increased from 0 to 10 wt% for the same dispersion grades. In addition, the Young's modulus value increased when the dispersion rate of the nano-clays was improved, showing the highest increase of around 93% for the nanocomposite with 10 wt% nano-clay. A comparison of the modeled mechanical properties and the experimental measurements values was performed to validate the modeling results. The simulated results fit well with the experimental values of Young's modulus.

A coupled crystal inelasticity-phase field model for crack growth in polycrystalline nitinol microstructures

Abstract This paper introduces a comprehensive computational framework, comprising a finite deformation crystal inelasticity constitutive model and phase field model, for modeling crack growth in superelastic nitinol polycrystalline microstructures. The crystal inelasticity model represents crystal stretching and lattice rotation from elastic mechanisms, as well as local inelastic deformation due to austenite-martensite phase transformation. The phase field formulation decomposes the Helmholtz free energy density into stored elastic energy, phase transformation energy, and crack surface energy components. The elastic energy accounts for tension-compression asymmetry with the formation of the crack through a spectral decomposition. Kinetic Monte Carlo simulations generate equilibrium area fractions of different surface orientations, which serve as weights for the surface energy. An adaptive wavelet-enhanced hierarchical finite element (FE) model is introduced to alleviate high computational overhead in phase field crack simulations. Simulations with the coupled inelasticity phase field model are conducted under various loading conditions including Mode-I tension, a quasi-static Kalthoff experiment, and cyclic loading of polycrystalline microstructures. Crack propagation is effectively predicted by this model, providing valuable insights into the material mechanical behavior with growing cracks.

Experimental and Constitutive Modeling Investigations of the Mechanical Behaviors of a Gravelly Soil Material Under Large-Size Triaxial Cyclic Tests

Experimental and Constitutive Modeling Investigations of the Mechanical Behaviors of a Gravelly Soil Material Under Large-Size Triaxial Cyclic Tests

Data-driven integration of synthetic representative volume elements and machine learning for improved microstructure-property linkage and material performance in ceramics

Ceramic materials, characterized by their heat resistance, dielectric properties, and mechanical performance, are often compromised by brittleness due to microstructural inclusions. These inclusions, such as defects, secondary phases, and pores, are critically analyzed for their impact on material strength and fracture properties. In this study, the linkage between microstructure and properties in ceramic materials is explored through a methodological approach that combines experimental observations with physics-based and machine learning models. A data-driven approach has been employed, utilizing synthetic Representative Volume Elements (RVEs) derived from X-ray computed tomography scans of ceramics. The methodology involves an automated finite element (FE) simulation process for progressive failure analysis under uniaxial compression and tension. The analysis is conducted using statistical, data-driven, and machine learning techniques, including principal component analysis and k-means clustering, to assess the microstructural features' impact on material performance. The study focuses on the use of RVEs for accurately capturing essential microstructural characteristics and addresses the challenges in developing synthetic models and the limitations of simulation capabilities. The study's findings reveal that less uniform inclusion distribution and a higher standard deviation in inclusion size correlate to lower mechanical performance. By applying data-driven methods, this research contributes to the optimization of material performance and the establishment of structure-property relationships, with a particular emphasis on the influence of inclusions and defects on the mechanical behavior of ceramic materials.

Uncovering all possible dislocation locks in face-centered cubic materials

Dislocation reactions and locks play an important role in the plastic deformation and mechanical behavior of crystalline materials. Various types of dislocation locks in face-centered cubic (FCC) materials have been reported in the literature pertaining to material-specific molecular-dynamic simulations and high-resolution transmission electron microscopy observations. However, it is unknown how many dislocation locks are possible, and how immobile all the dislocation locks are, with respect to each other. Here we present a discrete mathematics-based approach to reveal all possible dislocation locks in the FCC crystal structure. Totally eight types of dislocation locks are uncovered, resulting from all possible reactions of mobile/glissile (namely, perfect and Shockley partial) dislocations with (a) non-coplanar Burgers vectors which reside on two slip planes intersecting at both obtuse and acute angles and (b) coplanar Burgers vectors. We redefine the degree of dislocation lock immobility based on misorientations between non-close-packed lock planes and close-packed {111} slip planes. The subsequently derived sequences for the dislocation lock immobility and formation tendency are rationalized by the reported experimental and dislocation-dynamics simulation results.

Effects of Rubber Core on the Mechanical Behaviour of the Carbon-Aramid Composite Materials Subjected to Low-Velocity Impact Loading Considering Water Absorption.

The large-scale use of composite materials reinforced with carbon-aramid hybrid fabric in various outdoor applications, which ensures increased mechanical resistance including in impact loadings, led to the need to investigate the effects of aggressive environmental factors (moisture absorption, temperature, thermal cycles, ultra-violet rays) on the variation of their mechanical properties. Since the literature is still lacking in research on this topic, this article aims to compare the low-velocity impact behaviour of two carbon-aramid hybrid composite materials (with and without rubber core) and to investigate the effects of water absorption on impact properties. The main objectives of this research were as follows: (i) the investigation of the mechanical behavior in tests for two impact energies of 25 J and 50 J; (ii) comparison of the results obtained in terms of the force, displacement, velocity, and energy related to the time; (iii) analysis of the water absorption data; (iii) low-velocity impact testing of wet specimens after saturation; (iv) comparison between the impact behaviour of the wet specimens with that of the dried ones. One of the main findings was that for the wet specimens without rubber core, absorbed impact energy was 16% less than the one recorded for dried specimens at an impact energy of 50 J. The failure modes of the dried specimens without rubber core are breakage for both carbon and aramid fibres, matrix cracks, and delamination at matrix-fibre interfaces. The degradation for the wet specimens with rubber core is much more pronounced because the decrease in the absorbed impact energy was 53.26% after 10,513 h of immersion in water and all the layers were broken.

Strength and fabric anisotropy of granular materials under true triaxial configurations using DEM

This paper investigates the strength and fabric anisotropy of granular materials under true triaxial configurations using the discrete element method. A clump logic based on the multi-sphere approach was utilized to replicate realistic particle shapes. The evolutions of deviatoric stress and volumetric strains were found to be independent of mean effective stress, and intermediate principal stress parameter (b). From a macroscopic viewpoint, all samples exhibited initial strain hardening followed by softening behaviour, stabilizing at large deviatoric strains due to the dense state of assemblies. The peak state friction angles showed dependency on the b-value, aligning closely with experimental findings. Microscopic parameters such as mechanical coordination number and sliding fraction were found to be approximately independent of b-values. Additionally, non-coaxiality was observed for various b-values except for specific shear modes, aligning with previous findings using spherical particles. These results highlight the capability of clumped particles to simulate the true mechanical behaviour of granular materials and contribute to our understanding of their complex response under different loading conditions.

Analysis of the mechanical behavior of porous materials containing two populations of voids under dynamic spherical loading

A computational homogenization analysis is performed on three-dimensional representative volume elements (RVE) that contain two distinct populations of voids, to demonstrate the influence of the interaction between cavities. RVEs are constructed as cubic elastic perfectly-plastic matrices embedding two families of spherical voids, and subjected to dynamic loading under homogeneous kinematic boundary conditions. Multiple microstructure models are considered by varying the number, position, and size of voids to evaluate the micro-inertia contribution to the overall macroscopic stress. Velocity fields within numerical RVEs are investigated to reveal the interaction of voids and their key role in the dynamic macroscopic response of the porous material. Based on numerical simulation results, a homogenization analytical model considering void interaction is proposed to describe the mechanical behavior under dynamic loadings. This model relies on adjusting the strain rate level at the boundary of unit cells composing the porous material. Comparison between our numerical and analytical results with those obtained using the classical Taylor homogenization scheme highlights the limitations of the Taylor model in the case of porous materials under dynamic loading.