Analysis of end-restraint effects in DEM triaxial tests using linear and hertz contact models
Abstract
This study investigates how end restraint influences the relative peak-response performance of Linear and Hertz contact models in discrete element simulations of sand triaxial tests. Consolidated-drained tests were conducted on Xiamen medium sand at confining pressures of 100, 200, 300, and 400 kPa. A flexible-boundary model was developed through PFC–FLAC3D coupling. Both contact models were calibrated against the same experimental deviatoric stress–axial strain curve at 100 kPa; therefore, this case represents calibration rather than independent validation. With all micromechanical parameters fixed, the 200–400 kPa tests were used for independent validation. Based on peak stress and peak strain, the Linear model produced smaller deviations at 100 and 200 kPa, whereas the Hertz model performed better at 300 and 400 kPa. At 400 kPa, the Hertz model yielded a lower relative peak-stress error (1.19% versus 2.37%) and a smaller peak-strain deviation (0.4 versus 0.6% points) than the Linear model. For the modeled interface contrast µ = 0.2 for frictional ends and µ = 0 for ideal smooth ends—frictional ends increased simulated peak strength by approximately 0.9%–8.0%. This difference decreased with increasing confinement, from 5.57% to 0.90% for the Linear model and from 7.41% to 1.59% for the Hertz model between 100 and 400 kPa. End friction also restricted radial deformation near the loading platens and altered the specimen deformation pattern. These results show that contact-model agreement depends on both stress level and end boundary condition.
Related
- Child restraint system ownership and use in selected urban districts of Shenzhen and Shantou: a cross-sectional comparison
- A micro-learning-based education program about the principles of physical restraint using a mobile application and nursing students’ knowledge, attitudes, and performance regarding physical restraint for the elderly
- Host–Microbe Crosstalk in Sjögren’s Disease: Epithelial Type I Interferon Amplification, Loss of Short-Chain Fatty Acid Restraint, and the Effector-Phase Focus of Validated Therapy
- Behavioral and molecular characterization of fatigue-like manifestations triggered by acute restraint stress and LPS or saline
- From De-escalation to Restraint: A Qualitative Study of Emergency Clinician Decision-Making in Patient Agitation.
BibTeX
@article{liu2026analysis,
title = {{Analysis of end-restraint effects in DEM triaxial tests using linear and hertz contact models}},
author = {Jungang Liu and Lingxin Ge and Yi Zhang and Yifei Li and Jiawei Yu},
journal = {Scientific Reports},
year = {2026},
month = oct,
doi = {10.1038/s41598-026-75056-4},
url = {https://www.nature.com/articles/s41598-026-75056-4},
note = {CC BY-NC-ND 4.0},
}
RIS
TY - JOUR
TI - Analysis of end-restraint effects in DEM triaxial tests using linear and hertz contact models
AU - Jungang Liu
AU - Lingxin Ge
AU - Yi Zhang
AU - Yifei Li
AU - Jiawei Yu
PY - 2026
DA - 2026/10/09
JO - Scientific Reports
DO - 10.1038/s41598-026-75056-4
UR - https://www.nature.com/articles/s41598-026-75056-4
AB - This study investigates how end restraint influences the relative peak-response performance of Linear and Hertz contact models in discrete element simulations of sand triaxial tests. Consolidated-drained tests were conducted on Xiamen medium sand at confining pressures of 100, 200, 300, and 400 kPa. A flexible-boundary model was developed through PFC–FLAC3D coupling. Both contact models were calibrated against the same experimental deviatoric stress–axial strain curve at 100 kPa; therefore, this case represents calibration rather than independent validation. With all micromechanical parameters fixed, the 200–400 kPa tests were used for independent validation. Based on peak stress and peak strain, the Linear model produced smaller deviations at 100 and 200 kPa, whereas the Hertz model performed better at 300 and 400 kPa. At 400 kPa, the Hertz model yielded a lower relative peak-stress error (1.19% versus 2.37%) and a smaller peak-strain deviation (0.4 versus 0.6% points) than the Linear model. For the modeled interface contrast µ = 0.2 for frictional ends and µ = 0 for ideal smooth ends—frictional ends increased simulated peak strength by approximately 0.9%–8.0%. This difference decreased with increasing confinement, from 5.57% to 0.90% for the Linear model and from 7.41% to 1.59% for the Hertz model between 100 and 400 kPa. End friction also restricted radial deformation near the loading platens and altered the specimen deformation pattern. These results show that contact-model agreement depends on both stress level and end boundary condition.
SN - 2045-2322
ER -