Simulated Values Research Articles

Numerical and Experimental Analysis of Roller Hemming on Door Panel’s Curved and Straight-Edge of Flat Plane

Owing to its enhanced production efficiency, roller hemming has become the mainstream process for forming and joining metal sheets in the automotive industry. This study investigates the roller hemming process of a specific car door panel through a combination of experimental analysis and finite element analysis (FEA) on both straight-edge and curved-edge flat surfaces. Consequently, the mechanical properties of the door panel, including tensile strength, yield strength, modulus of elasticity, and Poisson’s ratio, were estimated through tensile testing and then underwent finite element modeling. The simulation results demonstrated the varying distribution of stress during the rolling hemming process, with the highest stress concentration observed in the bending area. Additionally, creepage and growing results were acquired from both simulation and experimental data to validate the precision of the numerical model. A comparison was made between the experimental and simulation results of the external forces exerted by the roller on the panel. In both straight- and curved-edge sections, the external force during final hemming exceeded that during pre-hemming, as revealed by experimental measurements of both normal and tangential external forces, surpassing their corresponding simulated values.

Comparative analysis of machine learning models for particle flow reconstruction

Abstract. Using the simulation data of different calorimeters, we investigate the use of machine learning and python algorithms for the simulation and reconstruction of particle flow energy in high-energy physics. We train models based on simulated pixel value image instead of common numerical data. We define two models that use different Regressor algorithms. Multi-layer Perceptron (MLP) Regressor offer stable and accurate prediction under less affection situation. Convolutional Neutral Networks (CNN) Regressor provided better and stronger modeling and reconstruction ability; however, it could cause overfitting results in certain situations. Furthermore, we optimize our models with the use of PyTorch. These models can serve as fast method for particle flow reconstruction in future studies and experiments.

Research on the viscosity characteristics of water glycol hydraulic fluids in deep-sea environment from molecular dynamics

The utilization of water glycol as a working medium represents a significant advancement in the realm of deep-sea hydraulic transmission systems. The viscosity of water glycol is a crucial parameter for hydrodynamic lubrication, yet it is profoundly influenced by the environmental conditions of the deep sea. This study employs the TIP4P/2005 water molecule model in conjunction with the optimized potentials for liquid simulations all-atom force field, and non-equilibrium molecular dynamics method to predict the viscosity of water glycol hydraulic fluids with varying compositions under deep-sea conditions. The simulation values agree well with the experimental results. Furthermore, this study introduces a fitting equation that accounts for the effects of composition, temperature, and pressure, enabling the prediction of the viscosity of water glycol hydraulic fluids within the 0–11 000 m sea water environment.

Design theory of self-centering column foot joint with dog bone weakened flange cover plate

Basing on the concept of post-earthquake repairability, plastic damage control, and self-centering, this study introduces a novel self-centering column foot joint with dog bone weakened flange cover plate (FCP) for prefabricated steel structure. The method and design theory for the target yield load and initial stiffness of the new joint were established. Parametric research was conducted using finite element simulation to assess the influence of the dog bone weakened FCP thickness, initial cable force, and axial compression ratio on the bearing performance of the new joint, thereby verifying the above method and theory. The results reveal the accuracy of the proposed method for determining the target yield load and initial stiffness, and provide practical design guidance for the new joint. The design theory values of the new joint and the simulation values demonstrate an error within 10%, which satisfies the engineering accuracy requirements. The corresponding procedures and suggestions for the design theory of this new joint were finally provided. Based on design theory, the overall plastic damage concentrates on the dog bone weakened FCP, successfully controlled by the new joint design, ensuring minimal damage to the main structure and demonstrating good bearing capacity and self-centering performance.

Ultraviolet Radiation Fields in Star-forming Disk Galaxies: Numerical Simulations with TIGRESS-NCR

With numerical simulations that employ adaptive ray-tracing (ART) for radiative transfer at the same time as evolving gas magnetohydrodynamics, thermodynamics, and photochemistry, it is possible to obtain a high-resolution view of ultraviolet (UV) fields and their effects in realistic models of the multiphase interstellar medium. Here, we analyze results from TIGRESS-NCR simulations, which follow both far-UV (FUV) wavelengths, important for photoelectric heating and polycyclic aromatic hydrocarbon excitation, and the Lyman continuum (LyC), which photoionizes hydrogen. Considering two models, representing solar neighborhood and inner-galaxy conditions, we characterize the spatial distribution and time variation of UV radiation fields, and quantify their correlations with gas. We compare four approximate models for the FUV to simulated values to evaluate alternatives when full ART is infeasible. By convolving FUV radiation with density, we produce mock maps of dust emission. We introduce a method to calibrate mid-IR observations, for example from JWST, to obtain high-resolution gas surface density maps. We then consider the LyC radiation field, finding most of the gas exposed to this radiation to be in ionization–recombination equilibrium and to have a low neutral fraction. Additionally, we characterize the ionization parameter as a function of the environment. Using a simplified model of the LyC radiation field, we produce synthetic maps of emission measure (EM). We show that the simplified model can be used to extract an estimate of the neutral fraction of the photoionized gas and mean free path of ionizing radiation from observed EM maps in galaxies.

An improved representation of aerosol in the ECMWF IFS-COMPO 49R1 through the integration of EQSAM4Climv12 – a first attempt at simulating aerosol acidity

Abstract. The atmospheric composition forecasting system used to produce the Copernicus Atmosphere Monitoring Service (CAMS) forecasts of global aerosol and trace gas distributions, the Integrated Forecasting System (IFS-COMPO), undergoes periodic upgrades. In this study we describe the development of the future operational cycle 49R1 and focus on the implementation of the thermodynamical model EQSAM4Clim version 12, which represents gas–aerosol partitioning processes for the nitric acid–nitrate and ammonia–ammonium couples and computes diagnostic aerosol, cloud, and precipitation pH values at the global scale. This information on aerosol acidity influences the simulated tropospheric chemistry processes associated with aqueous-phase chemistry and wet deposition. The other updates of cycle 49R1 concern wet deposition, sea-salt aerosol emissions, dust optics, and size distribution used for the calculation of sulfate aerosol optics. The implementation of EQSAM4Clim significantly improves the partitioning of reactive nitrogen compounds, decreasing surface concentrations of both nitrate and ammonium in the particulate phase, which reduces PM2.5 biases for Europe, the US, and China, especially during summertime. For aerosol optical depth there is generally a decrease in the simulated wintertime biases and for some regions an increase in the summertime bias. Improvements in the simulated Ångström exponent are noted for almost all regions, resulting in generally good agreement with observations. The diagnostic aerosol and precipitation pH calculated by EQSAM4Clim have been compared to ground observations and published simulation results. For precipitation pH, the annual mean values show relatively good agreement with the regional observational datasets, while for aerosol pH the simulated values over continents are quite close to those simulated by ISORROPIA II. The use of aerosol acidity has a relatively smaller impact on the aqueous-phase production of sulfate compared to the changes in gas-to-particle partitioning induced by the use of EQSAM4Clim.

Simulation Analysis of 3-D Airflow and Temperature Uniformity of Paddy in a Laboratory Drying Oven

This study analyzed the effects of airflow characteristics on the temperature distribution and drying uniformity of paddy during convective drying. Simulations of the drying process with varying airflow inlet and outlet positions were conducted using COMSOL Multiphysics 6.1 software. The determination coefficient (R2) between the simulated data and experimental values of Sample1 (S1), Sample2 (S2), and Sample3 (S3) was calculated, and its average values were 0.964, 0.963, 0.963, and 0.967, respectively. This study demonstrates that the airflow direction and outlet location have a significant impact on the temperature uniformity of the paddy. The vortex structure generated by the obstruction of the sidewalls and paddy influences both the airflow and temperature distribution within the drying chamber. When the outlet was on the left side and the inlet airflow was in a vertical orientation (VO), the temperature distribution of the paddy exhibited higher temperatures in the edge regions and lower temperatures in the center, with a maximum temperature difference of around 16 °C. The time required for the temperature to reach equilibrium with the outlet positioned on the left was 28.6% shorter than with the outlets positioned in the center or on both sides. Moreover, the temperature uniformity of the three paddy samples was better under this condition. The developed model accurately reflected the paddy drying process. It could also be used to analyze the optimal heating uniformity, providing a technical basis for the design of grain dryers.

Molecular Dynamics Simulation of the Viscosity Enhancement Mechanism of P-n Series Vinyl Acetate Polymer–CO2

Carbon dioxide (CO2) drive is one of the effective methods to develop old oil fields with high water content for tertiary oil recovery and to improve the recovery rate. However, due to the low viscosity of pure CO2, it is not conducive to expanding the wave volume of the mixed phase, which leads to difficulty utilizing the residual oil in vertical distribution and a low degree of recovery in the reservoir. By introducing viscosity enhancers, it is possible to reduce the two-phase fluidity ratio, expanding the degree of longitudinal rippling and oil recovery efficiency. It has been proven that the acetate scCO2 tackifier PVE can effectively tackify CO2 systems. However, little research has been reported on the microscopic viscosity enhancement mechanism of scCO2 viscosity enhancers. To investigate the influence of a vinyl acetate (VAc) functional unit on the viscosity enhancement effect of the CO2 system, PVE (Polymer–Viscosity–Enhance, P-3) was used as the parent, the proportion of VAc was changed, and the molecules P-1 and P-2 were designed to establish a molecular dynamics simulation model for the P-n-CO2 system. The molecules in the system under the conditions of 70 °C-10 MPa, 80 °C-10 MPa, and 70 °C-20 MPa were simulated; the viscosity of the system was calculated; and the error between the theoretical and simulated values of the viscosity in the CO2 system was relatively small. The difference between P-n molecular structure and system viscosity was analyzed at multiple scales through polymer molecular dynamics simulations and used the molecular radial distribution function, system density, accessible surface area, radius of gyration, minimum intermolecular distance, and minimum number of intermolecular contacts as indicators. This study aimed to elucidate the viscosity enhancement mechanism, and the results showed that the higher the proportion of VAc introduced into the molecules of P-n-scCO2 viscosities, the larger the molecular amplitude, the larger the effective contact area, and the greater the viscosity of the system. Improvement in the contact efficiency between the ester group on the P-n molecule and CO2 promotes the onset of solvation behavior. This study on the microscopic mechanism of scCO2 tackifiers provides a theoretical approach for the design of new CO2 tackifiers.

Robust malicious software detection and classification using global whale optimization algorithm with deep learning approach

Software malware detection and classification leverage sophisticated procedures and methods from the cybersecurity domain for identifying and categorizing malicious software, generally called malware. This procedure analyses code behaviour, file structures, and other features to distinguish between benign and malicious programs. Machine learning (ML) and artificial intelligence (AI) are vital in this domain, allowing the progress of dynamic and adaptive systems that identify novel and developing malware attacks. By training on massive datasets of benign and malicious instances, these systems learn patterns and signatures indicative of malware. This lets them correctly categorize and respond to potential attacks in real-time. This study presents a Global Whale Optimization Algorithm with Neutrosophic Logic for Software Malware Detection and Classification (GWOANL-SMDC) technique. The GWOANL-SMDC technique secures the software via the Android malware recognition process. Primarily, the GWOANL-SMDC technique employs the Neutrosophic Cognitive Maps (NCM) model for the feature selection process. The GWOANL-SMDC technique uses a convolutional long short-term memory (ConvLSTM) model for software malware detection. At last, the GWOA-based parameter tuning is performed to improve the performance of the ConvLSTM model. The simulation values of the GWOANL-SMDC technique are examined on the malware dataset. The obtained results ensured that the GWOANL-SMDC technique improved capability in detecting software malware.

Experimental study and simulation of a solid desiccant cooling system installed in northern Algeria

An experimental investigation of a solid desiccant cooling system for a research laboratory on the Algerian coast has been carried out. The system utilized flat-plate solar collectors for regeneration and was evaluated through a combination of individual component characterization, system modeling (using TRNSYS), and experimental validation. The results obtained have been validated and show that the efficiency of the rotary heat exchanger reached 65 %, the evaporative humidifier efficiency reached 94 %, and the desiccant wheel achieved efficiencies of 73 % on the process side and 66 % on the regeneration side. This translates to a significant dehumidification effect, reducing absolute humidity in the treated air by 9 g/kg, for the supply air temperature effectively maintained at a comfortable 24 °C. The close agreement between the simulated and experimental COP values (1.13 and 1.08, respectively) validates the chosen approach and confirms the potential of this technology for energy-efficient building cooling in the hot and humid Mediterranean climate of northern Algeria.

A Passenger Ship Emergency Evacuation Line Efficiency Evaluation Operator Considering Safety Capacity and its Application

Introduction: Aiming at evaluating the efficiency of passenger ship emergency evacuation routes, this study proposes an evaluation operator for passenger ship emergency evacuation routes, and then conducts simulation research on the passenger ship evacuation process. Method: The main innovations of this study are as follows. Firstly, this study proposes an evaluation operator for passenger ship emergency evacuation lines based on consideration of the passage capacity, ship heeling effect, traffic flow obstruction, personnel evacuation mood, accident location, and other factors. Secondly, this study discretizes the emergency evacuation path of passenger ships into a three-dimensional topological network structure and uses simulation data to calculate the traffic capacity between nodes. By comparing the calculated value with the simulated value, a risk assessment method for emergency evacuation is given. Thirdly, this study takes the three-story ro-ro passenger ship in the European "SAFEGUARD" project as an example and gives three passenger ship evacuation simulation processes under different passenger flow densities. Result: The results of the calculation example show that the risk-prone areas are mainly concentrated in the passenger ship stairs, and the risk value increases with the passenger flow. Conclusion: The calculation results verify the effectiveness of the emergency evacuation line efficiency evaluation operator proposed in this study.

Elucidating Protein Structures in the Gas Phase: Traversing Configuration Space with Biasing Methods.

Achieving accurate characterization of protein structures in the gas phase continues to be a formidable challenge. To tackle this issue, the present study employs Molecular Dynamics (MD) simulations in tandem with enhanced sampling techniques (methods designed to efficiently explore protein conformations). The objective is to identify suitable structures of proteins by contrasting their calculated Collision Cross-Section (CCS) with those observed experimentally. Significant discrepancies were observed between the initial MD-simulated and experimentally measured CCS values through Ion Mobility-Mass Spectrometry (IMS-MS). To bridge this gap, we employed two distinct enhanced sampling methods, Harmonic Biasing Potential and Adaptive Biasing Force, which help the proteins overcome energy barriers to adopt more compact configurations. These techniques leverage the radius of gyration as a reaction coordinate (guiding parameter), guiding the system toward compressed states that potentially match experimental configurations more closely. The guiding forces are only employed to overcome existing barriers and are removed to allow the protein to naturally arrive at a potential gas phase configuration. The results demonstrated close alignment (within ∼4%) between simulated and experimental CCS values despite using different strengths and/or methods, validating their efficacy. This work lays the groundwork for future studies aimed at optimizing biasing methods and expanding the collective variables used for more accurate gas-phase structural predictions.

Simulation of Summer Groundnut for Yield and Yield Attributing Traits Using CROPGRO-Peanut Model

Field experiments were carried out at College farm, B. A. College of Agriculture, Anand Agricultural University, Anand. The DSSAT v 4.6 CROPGRO-Peanut model was used to predict the phenology of groundnut crop under combinations of three sowing dates and four groundnut cultivars. The model was calibrated with a 2015 dataset of growth and phenological parameters for estimating the genetic coefficients of all four cultivar and was validated with a 2016 dataset of the same parameters. Result found that the model was able to reasonably simulate the pod yield, kernel yield and haulm yield with per cent error ranging (± 10.06) between observed and simulated value for all cultivars under different sowing dates Simulations of yield and yield attributing characters using the calibrated model were found to be quite accurate.

Development of a chemical mechanism for ammonia/hydrogen blends for engines

ABSTRACT Ammonia/hydrogen blend fuels have been extensively used as alternative fuel for engines. The combustion and emissions characteristics of ammonia/hydrogen engines have been investigated utilising a CFD model. Developing a robust and accurate chemical reaction mechanism for ammonia/hydrogen blend fuels is of utmost importance. In this paper, a chemical mechanism has been developed and validated for the combustion characteristics of ammonia/hydrogen, which includes 31 species and 223 reactions. Extensive validation has been conducted using experimental data as the basis. The obtained data indicate that the simulated values for the ignition delay times of the shock tube and the concentration of the main components in the stirred jet reactor align well with the measured values. The simulated laminar flame speed under turbulent flow conditions shows a remarkable agreement with the corresponding experimental data. The study compared simulated flame evolution values with experimental data from a constant volume combustion bomb. The results indicate that the proposed mechanism is able to produce favourable flame evolution and combustion characteristics. The presented detailed mechanism demonstrates credible overall performance in terms of combustion behaviours, indicating its suitability for effectively modelling ammonia/hydrogen combustion in real engines.

Design and study of crashworthiness for square vicsek fractal hierarchical multicellular tubes under axial compression

Through extensive examination of thin-walled structures, the efficacy of enhancing the crashworthiness of thin-walled tubes via layered structures has been confirmed. In this study, we introduce a novel design called the square Vicsek fractal hierarchical multicellular tube (SVFHMT), which integrates principles from fractal science, specifically Vicsek fractals. Unlike conventional layered structures that align along a single direction, our design incorporates layers evenly distributed at five positions throughout the structure, resulting in a more uniform material distribution compared to traditional layering methods. In this paper, the verified finite element model is used to analyze the crashworthiness of the structure, and the effects of levels, wall thickness, wall thickness ratio of different levels and height to diameter ratio on the crashworthiness of the structure are discussed. The mean crushing force (MCF) of SVFHMT is predicted by Simplified Super Folding Element theory. Our findings highlight the significant influence of layering, wall thickness, and ratios of wall thickness between layers on the crashworthiness of the structure. Specifically, for structures of equal mass, the energy absorption of second-order SVFHMTs and third-order SVFHMTs demonstrates substantial improvements compared to square nine-cell tubes, with increases of 51.05% and 91.70%, respectively. Moreover, when maintaining equal mass for second-order SVFHMTs with differing ratios of wall thickness between layers, specific energy absorption ranges from 15.69 kJ/kg to 21.29 kJ/kg. Notably, the maximum error between theoretical and simulated values of the mean crushing force for each order of SVFHMT is only 7.34%. This underscores the effectiveness of our novel layer design strategy, which diverges significantly from conventional approaches and offers valuable insights for advancing layered structures further.

Clinician assessment of kidney function from plasma creatinine values during critical illness: A scenario-based international multi-professional survey

PurposeDuring critical illness interpretation of serum creatinine is affected by non-steady state conditions, reduced creatinine generation, and altered distribution. We evaluated healthcare professionals' ability to adjudicate underlying kidney function, based on simulated creatinine values. MethodsWe developed an online survey, incorporating 12 scenarios with simulated trajectories of creatinine based on profiles of muscle mass, GFR and fluid balance using bespoke kinetic modelling. Participants predicted true underlying GFR (<5, 5–14, 15–29, 30–44, 45–59, 60–90, >90 ml.min−1.1.73 m−2) and AKI stage (stages 1–3, defined as 33 %, 50 %, 66 % decrease in GFR from baseline) during the first 7-days and at ICU discharge. Results100 of 103 respondents from 16 countries, 94 completed 1 or more scenarios. 43(43 %) were senior physicians, 74(74 %) critical care and 31(31 %) nephrology physicians. Over the first 7-days, true GFR was correctly estimated 43 % of the time and underlying AKI stage in 57 % of patient days. At ICU discharge GFR was predicted 35 % of the time. At all timepoints, over and under-estimation of GFR was observed. ConclusionParticipants displayed marked variation in estimation of kidney function, suggesting difficulty in accounting for multiple confounders. There is need for alternative, unbiased measures of kidney function in critical illness to avoid misclassifying kidney disease.

Estimation of Regional Electricity Consumption Using National Polar-Orbiting Partnership’s Visible Infrared Imaging Radiometer Suite Night-Time Light Data with Gradient Boosting Regression Trees

With the rapid development of society and economy, the growth of electricity consumption has become one of the important indicators to measure the level of regional economic development. This paper utilizes NPP-VIIRS nighttime light remote sensing data to model electricity consumption in parts of southern China. Four predictive models were initially selected for evaluation: LR, SVR, MLP, and GBRT. The accuracy of each model was assessed by comparing real power consumption with simulated values. Based on this evaluation, the GBRT model was identified as the most effective and was selected to establish a comprehensive model of electricity consumption. Using the GBRT model, this paper analyzes electricity consumption in the study area across different spatial scales from 2013 to 2022, demonstrating the distribution characteristics of electricity consumption from the pixel level to the city scale and revealing the close relationship between electricity consumption and regional economic development. Additionally, this paper examines trends in electricity consumption across various temporal scales, providing a scientific basis for the optimal allocation of energy and the effective distribution of power resources in the study area. This analysis is of great significance for promoting balanced economic development between regions and enhancing energy efficiency.

Improving sustainability of rice-canola rotation through water and nitrogen management in a humid region

Intensive cropping systems face significant environmental and economic challenges due to current water and nutrient management practices. This study addresses the limited research on optimizing these practices in rice-canola rotations, focusing on enhancing sustainability through improved water and nitrogen (N) management. Utilizing the DSSAT crop growth model, we evaluated three irrigation scenarios for rice- I1 (100 % irrigation requirement: IR), I2 (80 % IR), and I3 (60 % IR)-alongside three N application levels: N1 (30 % below RAN: regional average N), N2 (RAN), and N3 (30 % above RAN). Key sustainability indicators, including water productivity (WP), nitrogen use efficiency (NUE), and economic productivity, were analyzed. The DSSAT model demonstrated strong predictive accuracy, with yield differences of only 236 kg ha⁻¹ for rice and 121 kg ha⁻¹ for canola between observed and simulated values. A 40 % reduction in water use increased rice yield to a maximum of 5654 kg ha⁻¹ under I3N2 conditions. Additionally, this reduction in water consumption resulted in a 27 % decrease in production costs, significantly enhancing WP, net benefits, and gross benefits of rice by 46.7 %, 43.1 %, and 47.9 %, respectively. Conversely, a 30 % increase in N application for canola correlated with a 9 % yield increase, underscoring the importance of tailored N strategies. Furthermore, N3 raised WP, NUE, net benefits, and gross benefits of canola by 6.7 %, 1.4 %, 9.1 %, and 9.2 %, respectively, compared with current N application level. Overall, this study illustrates the potential for substantial water savings in rice cultivation while recommending stable or reduced N application levels. These management strategies enhance farm performance, lower production costs, and contribute to the sustainability of rice-canola cropping systems, promoting both agricultural productivity and environmental stewardship.

Simulated resections and responsive neurostimulator placement can optimize postoperative seizure outcomes when guided by fast ripple networks.

In medication-resistant epilepsy, the goal of epilepsy surgery is to make a patient seizure free with a resection/ablation that is as small as possible to minimize morbidity. The standard of care in planning the margins of epilepsy surgery involves electroclinical delineation of the seizure-onset zone and incorporation of neuroimaging findings from MRI, PET, single-photon emission CT and magnetoencephalography modalities. Resecting cortical tissue generating high-frequency oscillations has been investigated as a more efficacious alternative to targeting the seizure-onset zone. In this study, we used a support vector machine (SVM), with four distinct fast ripple (FR: 350-600 Hz on oscillations, 200-600 Hz on spikes) metrics as factors. These metrics included the FR resection ratio, a spatial FR network measure and two temporal FR network measures. The SVM was trained by the value of these four factors with respect to the actual resection boundaries and actual seizure-free labels of 18 patients with medically refractory focal epilepsy. Leave-one-out cross-validation of the trained SVM in this training set had an accuracy of 0.78. We next used a simulated iterative virtual resection targeting the FR sites that were of highest rate and showed most temporal autonomy. The trained SVM utilized the four virtual FR metrics to predict virtual seizure freedom. In all but one of the nine patients who were seizure free after surgery, we found that the virtual resections sufficient for virtual seizure freedom were larger in volume (P < 0.05). In nine patients who were not seizure free, a larger virtual resection made five virtually seizure free. We also examined 10 medically refractory focal epilepsy patients implanted with the responsive neurostimulator system and virtually targeted the responsive neurostimulator system stimulation contacts proximal to sites generating FR at highest rates to determine if the simulated value of the stimulated seizure-onset zone and stimulated FR metrics would trend towards those patients with a better seizure outcome. Our results suggest the following: (i) FR measures can accurately predict whether a resection, defined by the standard of care, will result in seizure freedom; (ii) utilizing FR alone for planning an efficacious surgery can be associated with larger resections; (iii) when FR metrics predict the standard-of-care resection will fail, amending the boundaries of the planned resection with certain FR-generating sites may improve outcome and (iv) more work is required to determine whether targeting responsive neurostimulator system stimulation contact proximal to FR generating sites will improve seizure outcome.

Modeling and chewing parameters calibration of Euryale ferox based on discrete element method.

Euryale ferox was chosen for this study to examine its mechanical properties during chewing. Experiments and the discrete element method were used to conduct the study. Initially, the intrinsic and contact parameters of E. ferox were established through physical tests. The maximum compressive force of breakage and chewing mechanical properties (hardness, springiness, and chewiness) were measured using a texture profile analyzer (TPA). A Hertz-Mindlin with bonding model of the E. ferox was constructed. Optimal values for significant factors, including the normal stiffness per unit area, shear stiffness per unit area, and bonding radius, were obtained through single-factor, Plackett-Burman, steepest ascent, and Box-Behnken response surface tests, with the maximum compressive force as the index. Under optimal parameter combination conditions, the relative error between the simulated and experimental values of the maximum compressive force was 0.79%, and the relative errors between the simulated and TPA test values of all indicators were less than 5.65%. This study provides a valuable reference for simulating the textural characteristics of agricultural products.