Simulation Of Deformation Research Articles

Real-time simulation for multi-component biomechanical analysis using localized tissue constraint progressive transfer learning

In virtual surgical training, it is crucial to achieve real-time, high-fidelity simulation of the tissue deformation. The anisotropic and nonlinear characteristics of the organ with multi-component make accurate real-time deformation simulation difficult. A localized tissue constraint progressive transfer learning method is proposed in this paper, where the base-compensated dual-output transfer learning strategy and the localized tissue constraint progressive learning architecture are developed. The proposed strategy enriches the multi-component biomechanical dataset to fully represent complex force-displacement with minimal high-quality data. Meanwhile, the proposed architecture adopts focused and progressive model to accurately describe tissues with varied biomechanical properties rather than singular homogeneous model. We made comparison with 4 state-of-the-art (SOTA) methods in simulating multi-component biomechanical deformations of organs with 100 pairs of testing data. Results show that the accuracy of our method is 50% higher than other methods in different validation matrix. And our method can stably simulate the deformations in 0.005 s per frame, which largely improves the computing efficiency.

A yarn-to-fabric multiscale modeling strategy for auxetic behavior analysis of negative Poisson’s yarn-based fabric

Auxetic behavior is critical to determining the deformation and mechanical performance of materials. Auxetic mechanisms of the woven fabric have previously not been well understood, especially for those containing negative Poisson’s ratio yarn (NPRY). Most theoretical models were developed from the macro geometric structure which lacks analysis of yarn levels. This paper aims to use numeral analysis combined with finite element modeling to investigate the auxetic behavior of the NPRY-based fabric. Firstly, the theoretical model that reveals the relationship between the auxetic behaviors of yarn and NPRY-based fabric was established based on component yarn interweaving rules, deformation, and stress. Following, a strain-dominated deformation simulation was further performed to clarify the micro stress and auxetic mechanism of fabric by calculating the stress–strain response and NPR value. Moreover, the experiments of NPRY-based fabric were also developed to validate the theoretical and simulation results, Results indicate that the model can well predict the auxetic behavior of fabric with an error under 7%. The micro–macro analysis of the auxetic behavior from yarn to fabric provides a fundamental understanding of the auxetic mechanism of NPRY-based fabric. Additionally, the yarn-to-fabric model can used to study the auxetic behavior of NPR fabric with variable parameters thus supplying theoretical guidance for engineering.

Computer simulation help predict the frame deformation following a Venus-A transcatheter aortic valve implantation in patients with pure aortic regurgitation: a retrospective study.

Patient-specific computer simulation of transcatheter aortic valve implantation (TAVI) predicts the interaction between an implanted device and the surrounding anatomy. In this study, we validated the predictive value of computer simulation for the frame deformation following a Venus-A TAVI implant in patients with pure aortic regurgitation (AR). Furthermore, we used the validated computational model to evaluate the anchoring mechanism within the same cohort. This was a retrospective study. FEops HEARTguide technology was used to simulate the virtual implantation of a Venus-A valve model in a patient-specific geometry. The predicted frame deformation was quantitatively compared to the postoperative device deformation at multiple levels. The outward forces acting on the frame were extracted for each patient and the total outward force acting around the aortic annular (AA) and sinotubular junction (STJ) planes were recorded. Thirty patients were enrolled in the study with 10 in the migration group and 20 in the non-migration group. The dimensions of the simulated and observed frames had good correlations at Dmax (R2=0.88), Dmin (R2=0.91), perimeter (R2=0.92), and area (R2=0.92). The predicted outward force acting on the frame at the AA level was comparable between the migration and no-migration groups. The predicted outward force acting on the frame at the STJ level was always significantly higher in the migration group than the no migration group at different bandwidths: 3 mm (P=0.002), 5 mm (P=0.005), 10 mm (P=0.002). Patient-specific computer simulation of TAVI accurately predicted frame deformation in Chinese patients with pure AR. The forces at the STJ facilitated stabilization of the device within the aortic root, which might be used as a discriminator to identify patients at risk of device migration prior to intervention.

Effectiveness of thermoviscoplastic material models in predicting the thermomechanical behavior of rolled homogenous armor steel

This paper studies the dynamic and elevated temperature behavior of rolled homogeneous armor steel, which is crucial for applications such as blast, impact, crash, and ballistic resistance. The deformation behavior is analyzed through tensile tests. The tensile tests of RHA steel are conducted at quasi-static 10−4s−1≤εṗ≤10−1s−1 to dynamic 10−1s−1≤εṗ≤36.439s−1 strain-rates and at high temperatures ranging from room temperature (RT,27°C) to 500°C. Phenomenological Johnson-Cook (JC) models and semi-physical Rusinek-Klepaczko (RK) material models are employed in finite element (FE) simulations of tensile tests. The JC failure model, with a small modification, is used along with the material models to simulate the loss of load-bearing capacity of the component. All necessary material parameters are extracted from the test data. The material models and the damage model are integrated into ABAQUS CAE FE software through a user-defined material subroutine (UMAT). The simulated outcomes acquired from diverse material models are validated with relevant experimental findings, and the effectiveness of each material model is evaluated through qualitative and quantitative assessments. Observations reveal limitations in the existing models to accurately replicate material behavior under tension, particularly in the dynamic strain rate, and elevated temperature range. Consequently, a modified version of the JC model (MJC) is explored to better account for the impact of temperature on strain hardening and strain rate hardening. The proposed modification in the material model significantly enhances the predictive accuracy of FE simulations of tensile deformation across the entire spectrum of studied strain rates and temperatures compared to the original JC and RK models.

Seismicity changes and numerical simulation of coseismic deformation following the 2022 Ms6.8 Luding earthquake in Sichuan, China

Seismicity changes and numerical simulation of coseismic deformation following the 2022 Ms6.8 Luding earthquake in Sichuan, China

Modeling time-dependent deformation in concrete: A fractional calculus method with finite element implementation

This study presents a novel numerical method to simulate the time-dependent behavior of concrete materials. The deformation response of both traditional and fractional calculus (FC) viscoelastic models is analyzed by the proposed method, in which models’ one-dimensional constitutive relationships are extended to the three-dimensional form. By applying numerical and finite difference methods based on Caputo-type fractional integration, the stress-strain behavior of the FC viscoelastic model is discretized over the time scale. This method exhibits broad application prospects, and can be employed to represent various forms of viscoelastic models. For the concrete material, suitable FC viscoelastic models are developed. To facilitate practical engineering applications of the proposed model, a numerical solution algorithm is implemented in the finite element (FE) analysis through the User-defined Material Mechanical Behavior (UMAT) interface of the commercial FE software ABAQUS. Finally, FE results are compared with the experimental results from several concrete creep tests. The consistency between the FE results and experimental data confirms the effectiveness of the proposed model in describing concrete behavior. The proposed method and model provide theoretical and numerical support for a more profound understanding and simulation of the time-dependent deformation of RC specimens in practical engineering applications.

Design analysis and development of a gantry robot for multi-layer 3-D concrete printing with simulation and experimental validation

The present study investigates the design analysis and development of a framed gantry robot (FGR) for 3-D concrete printing. Finite element simulations are performed for design validation and parameter estimation to conduct detailed analysis of the concrete printer. The kinematic and dynamic models of FGR are used for workspace analysis, motion control study, and motor sizing for printer development. The motion trajectory of the printer nozzle tip is evaluated using MATLAB simulations. An integrated programmable logic controller (PLC) with an in-built human-machine interface (HMI) through control algorithms has been implemented for interactive, user-friendly, and precise motion control. A separate control system is developed to monitor the pumping power and feed rate of mortar mix during printing. The simulated control scheme is validated experimentally using trajectory tracking for various input printing paths through HMI, with a maximum tracking error of 1 mm. Successful concrete printing trials have been conducted using various mix-design, following extensive evaluation of different printer and material parameters at the laboratory scale. The proposed design concept, mechanisms, interactive control interface, and tracking accuracy have been validated through simulation and experimental results, for ensuring quality concrete prints by the developed printer. Material testing parameters have been evaluated experimentally for modeling various mix-designs. To further evaluate printing performance, simulations of layer deformation in printed structures are conducted using finite element analysis, incorporating non-linear material models and various contact layer interactions. The simulated results of printed layers of different shapes of structures show a close resemblance with the experimental concrete printing trials.

Structural Evaluation of a Three-Wheeler VehicleUsing Simulation Tools: A Case Study in Mexico

The mototaxi, a three-wheeled vehicle equipped with a roof, is a widely utilized mode of transportation in Mexico. Typically, it is employed for short-distance journeys in exchange for payment, similar to the operation of a conventional taxi. This study conducts astructural analysis of a mototaxi-type vehicle utilized in Mexico to assess its performance and safety. It underscores the significance of this mode of transportation, widely relied upon by numerous individuals. A product design and development methodology was employed, utilizing torsional deformation simulations to validate the new geometry. The objective was to minimize torsions as much as possible, thereby enhancing the motorcycle taxi’s safety and ensuring the vehicle’s correct positioning. Through computeraided design, the prevailing torsions within the casing were assessed, establishing the operating conditions to which the system is commonlysubjected. The findings from the chosen vehicular structure reveal a flexural rigidity of 6,508.15 N/mm, torsional rigidity of 27.35 KNm/°, and a range of natural frequencies between 8-21 Hz. These values indicate favorable resistance against bending forces and operational frequency. However, the torsional results exhibit deficiencies, suggesting an unsafe structure for all motorcycle taxi occupants. Consequently, technology developers and national legislators should prioritize enhancing the structural integrity of such vehicles.

A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening

A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. The method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.

Simulation of Automobile Structural Member Deformation and Crash via Isogeometric Analysis

We conducted an isogeometric analysis (IGA) to evaluate the performance of automobile structural member deformation and crash. In automobile crash analysis, ensuring the accuracy of the acceleration, velocity, and load in time series, as well as the structural deformation behavior, is important. To maintain the aforementioned consistency, accurately reproducing the bending and buckling of structural members is indispensable. In this study, we firstly computed the bending and buckling of structural members using IGA and validated its performance by comparing the results with those of a conventional finite element analysis and experiments. In addition, we utilized IGA for the crash analysis of an automobile body.

Study on the influence of rainfall infiltration on the stability of dump slope

The present strength reduction finite element methods usually only reduce strength parameters, which result in the deformation state of earth-rock aggregate slope are different from the practical ones. The article holds that under rainfall infiltration, not only the strength parameters of earth-rock aggregate slope will reduce, but also the deformation parameters will also change with the changing of soils stress levels and strength parameters. Using the developed FDM (Finite Difference Method), the numerical simulation for earth-rock aggregate slope deformation under different rainfall infiltration depths are conducted, and the results of numerical calculation are compared with those of the laboratory model tests. Then the influence law of rainfall infiltration depths on the road-embankment slope stability and sliding arc changing characteristics are analyzed by strength reduction finite element methods.

A hybrid thermo-mechanical phase-field model for anisotropic brittle fracture

The phase-field modeling approach is integrated with the isogeometric-meshfree approach to investigate the coupled thermo-mechanical fracture behaviors of composite materials, accounting for material anisotropy. The present approach is formulated within a thermodynamic framework, exploring the effects of both thermal and mechanical factors on crack evolution. The crack driving force in composites is determined by the collective contributions from both fibers and the matrix. Moreover, a hybrid computational framework, enhanced with an adaptive mesh refinement technique, is proposed for the efficient implementation of thermo-mechanical fracture. Finally, simulations of the thermal deformation and fracture processes in both homogeneous and composite materials are implemented to verify the present approach. The impact of the weak material anisotropy, such as anisotropic thermal conductivity and material strength, along with the strong anisotropy coefficient on fracture patterns in orthotropic composites is also investigated. The simulation results demonstrate that the developed approach can predict complex failure patterns in composites under thermo-mechanical loading, thus offering insights for the design of composites.

Atomic simulation of phosphorus segregation in nickel grain boundaries and its impact on deformation mechanisms

The segregation behavior of phosphorus on the grain boundary (GB) structure and mechanical response has been studied in the Ni–P system. Fourteen twist and tilt bi-crystals have been selected from the Olmsted dataset. Then, hybrid Monte Carlo/Molecular dynamics simulations (MC/MD) are used to obtain the equilibrium GB at 300 K and 1000 K, respectively. The present study indicates that solute phosphorus dissolved in the nickel matrix is inhomogeneous segregated to all of the GBs. GB characters including GB free volume and energy are found to correlate with segregation. The participated phase with Ni3P composition was observed only at 1000 K and for the lower energy GBs. Furthermore, we performed simulations of tensile deformation for the GBs with and without segregation. The twist and tilt GBs initially exhibit linear and non-linear elasticity behavior, respectively. The twist GBs demonstrate an increased yield strength with the help of phosphorus segregation. Besides, the thermal annealing process at 1000 K promotes uniform segregation at GB, leading to further strength enhancement. Conversely, the strength is weakened in segregated tilt GBs. This work illustrates the relationships among composition, GB structure, and mechanical properties.

A coupled FEM-FFT concurrent multiscale method for the deformation simulation of CFRPs laminate

Carbon fiber reinforced polymer composites (CFRPs) laminate is increasingly used in aircraft structures and its assembly deformation has great influence on the aircraft aerodynamic profile and mechanical performance. Due to the inherent multiscale nature of CFRPs, the traditional macroscopic models are unable to incorporate their true microstructural details, resulting in them being unable to provide both the macro and micro deformations simultaneously. To address this issue, this work proposed a coupled FEM-FFT concurrent multiscale method for the deformation simulation of CFRPs laminate, where the mechanical response of each macro point of CFRPs laminate was from the corresponding unidirectional representative volume element (UD RVE). The FEM (Finite Element Method) was used for the analysis at the macroscopic scale, while the FFT (Fast Fourier Transform) based method was employed in the UD RVE to speed up the calculation. The elastoplastic behavior of resin matrix in the UD RVE was modelled with the parabolic yield criterion. The proposed concurrent multiscale method was validated respectively by the tension, compression and shearing experiments of CFRPs laminate. The results show that there is a good agreement both for the full strain field and load–displacement curve, demonstrating the effectiveness of the proposed concurrent multiscale method.

Density dependent constitutive model for Bi-2212 powder compression deformation

Bi-2212 HTS materials are fabricated into multi-filamentary wires via powder-in-tube (PIT) method followed by proper heat treatment to obtain superconductivity, but how to predict the large compression deformation behaviors of the Bi-2212 powder is critical to design the processing of the Bi-2212 HTS wire. Drucker Prager/Cap (DPC) model was the most commonly used model for powders including Bi-2212 with soil-like mechanical behavior to consider its shear failure as well as hydrostatic compression. However, the parameters for DPC Cap evolve with densities change and the original model is inadequate to precisely describe the densification process of Bi-2212 powder with large strain. In this study, the modified DPC model with density dependent parameters was introduced for Bi-2212 powder compressions by measuring the failure strength and hydrostatic compressive behavior under different density states. The DPC yield surface was plotted with an evolution trend of non-linear outward expansion with density increased. FEM model of uniaxial compression based on the as-introduced model was built with subroutine VUSDFLD applied. The distribution of Mises stress and relative density were analyzed. The axial stress-density curve for FEM and experimental results were normalized and quantitively evaluated by Mean Square Error (MSE). The introduced model shows good convergence and could match the experimental results well with normalized MSE of 0.000207 and Root Mean Square Error (RMSE) of 0.0144, indicating the mean error percentage of 1.44%. The model introduced in this article provides supports toward large strain deformation simulation of Bi-2212 powder.

Efficient boundary conditions identification in thermal simulation of the spindle system with reduced order model and differential evolution algorithm

The boundary conditions in thermal simulation of the spindle system are very complicated and the empirically calculated ones will usually result in a significant discrepancy between experiment and simulation. Thus, the optimal boundary conditions were identified by treating it as an inverse optimization problem in this paper. Moreover, to advance the optimization efficiency, a reduced order model is established to replace the time-consuming thermal simulation of the spindle system. By superposition of truncated eigenmodes and accurately predicted modal coefficients, the reduced order model can reconstruct the field very similar to thermal simulation. Experimental reconstruction of both the temperature field and thermal deformation field by the reduced order model demonstrated its at least 360 times speedup with 99.8 % accuracy compared to the actual thermal simulation. With the accurate reduced order model as a lightweight digital twin and using the differential evolution algorithm, three types of boundary conditions, i.e., the heat generation rates, the convective heat transfer coefficients and the thermal contact resistances, under the shaft rotation speed of 4,000r/min were identified within 11s. The maximum temperature simulation error was reduced from 85.6 % to 6.6 % and the thermal deformation simulation error was reduced from 60.9 % to 10.8 %.

Experiment-free multiscale simulation of residual deformation in non-crimp fabric composites

The effectiveness of multiscale modeling, consisting of quantum chemical reaction path calculations, molecular dynamics simulations, and microscopic and macroscopic finite element analysis, developed for the process-induced residual deformation of carbon fiber-reinforced plastic was experimentally evaluated. In this study, non-crimp fabric composites with an arbitrary thermoset resin were fabricated, and process-induced deformation was measured. Process-induced residual deformations due to cure shrinkage and thermal shrinkage of cross-ply laminates under conditions similar to those in the experiment were predicted using multiscale modeling. The predicted deformation shape and maximum deformation value are in good agreement with the fabrication experiments. Multiscale modeling can consider the combined factors that affect the deformation (i.e. size, elastic moduli, coefficient of thermal expansion, and fiber volume fraction) and can provide important knowledge regarding the development of high-performance composite structures and for stable manufacturing.

Electric-driven shape memory deformation control and verification of composite films based on PVA/WPU/SCFs

Based on the research of Polyvinyl alcohol (PVA)/Waterborne Polyurethane (WPU)/Short Carbon Fibers (SCFs) novel electrochemically shape memory composites, this study delved into an exploration of factors influencing the shape memory deformation of shape memory film samples when subjected to voltage stimulation. These main affecting factors encompassed sample thickness, activation voltage, and wire layout. Subsequently, based on the above conclusions, the viscoelastic constitutive model of PVA/WPU/SCFs composite was established based on Maxwell model and the deformation deflection of the composite film was predicted by shape memory recovery, and the electrothermal coupling simulation analysis and shape memory recovery and deformation simulation verification were carried out. Finally, based on the shape memory recovery performance of the composite film, a new reciprocating moving deformation structure was conceptualized, characterized by distinct cycles of motion, each analyzed individually.

Multi-phase-field lattice Boltzmann simulations of semi-solid simple shear deformation in thin film

Semi-solid deformation during casting often results in significant solidification defects, such as segregation bands. Consequently, the development of a numerical simulation tool is crucial for accurately replicating semi-solid deformation. In our previous study, we applied a multi-phase-field lattice Boltzmann (MPF-LB) model to semi-solid deformation, facilitating seamless simulation from polycrystalline solidification to semi-solid deformation in a two-dimensional (2D) problem. This study extends the 2D MPF-LB model to a three-dimensional (3D) problem and develops a simulation method for semi-solid simple shear deformation in thin films. To enhance the efficiency of the 3D semi-solid simulation, we implemented parallel computations using multiple graphics processing units. Through a discussion of the relationships among the stress–strain curve, grain rearrangement behavior, and fluid flow, we confirmed that the developed 3D MPF-LB model successfully reproduced the characteristic phenomena of semi-solid deformation, and has high potential to investigate the nuanced mechanisms of semi-solid deformation.

Inrushing failure characteristics of deep shaft structures in soft soils based on material point method simulation

In areas with soft soil, inrushing may occur during the construction of deep shaft structures under high water pressures, causing destructive accidents. Most studies on inrushing have used the finite element method. However, it is difficult to achieve a large deformation simulation because of mesh limitations. The material point method is not constrained by the mesh, and the material properties and motion information of the material are stored at the material points. Therefore, it can be used to effectively simulate large soil deformations. In this study, a numerical model for a deep shaft structure in Shanghai was formulated using the material point method to simulate the evolution of inrushing failure. The likelihood of inrushing damage at different excavation depths was analyzed. Soil heave at the bottom of the pit in the absence of inrushing was examined. The structural and environmental factors affecting the inrushing process are discussed. Strategies have been proposed from different perspectives, such as structure, environment, and construction methods, to reduce the damage caused by inrushing, thereby providing a reference for the safer construction of deep shaft structures.