Internal Stress Distribution Research Articles

Tunable Mechanics and Micromechanism in Close-Knit Silicide-in-SiO2 Core-Shell Nanowires.

Bending/tension mechanics is one of the core issues for nanowires in flexible free-standing transport and sensor applications, but it remains a challenge to tailor the mechanical performance beyond the inherent properties. Herein, based on structure engineering, silicon-based Mn5Si3@SiO2 nanocables are proposed and demonstrated as versatile nanosystems. Except for outstanding toughness, large ultimate strain, and great strength, they display diverse mechanical behaviors such as simplex elasticity, plasticity, and viscoelasticity under different external conditions. The tunable performances originate from synergetic effects between the core and shell components, like the atomic bonding transitional interface and space confinement, which induce optimizing internal stress distribution and the dislocation evolution mechanism in the core. The related mechanical performance is revealed carefully. The bending and tension dynamic picture, quantitative force curve, stress-strain dependence, and the corresponding lattice evolution are acquired by in/ex situ characterizations and measurements. These results contribute to nanowire mechanical design and also expand to strain-regulated three-dimensional multifunctional nanosystems.

Mechanochemistry of Spiropyran under Internal Stresses of a Glassy Polymer.

Mechanophores are powerful molecular tools used to track bond rupture and characterize mechanical damage in polymers. The majority of mechanophores are known to respond to external stresses, and we report in this study the first precedent of a mechanochemical response to internal, residual stresses that accumulate during polymer vitrification. While internal stress is intrinsic to polymers that can form solids, we demonstrate that it can dramatically affect the mechanochemistry of spiropyran probes and alter their intramolecular isomerization barriers by up to 70 kJ mol-1. This new behavior of spiropyrans (SPs) enables their application for analysis of internal stresses distribution and their mechanochemical characterization on the molecular level. Spectroscopy and imaging based on SP mechanochemistry showed high topological sensitivity and allowed us to discern different levels of internal stress impacting various locations along the polymer chain. The nature of the developed technique allows for wide-field imaging of stress heterogeneities in polymer samples of irregular shapes and dimensions, making it feasible to directly observe molecular-level manifestations of mechanical stresses that accompany the formation of a vast number of solid polymers.

An investigation into the hammer toe effects on the lower extremity mechanics and plantar fascia tension: A case for a vicious cycle and progressive damage

Hammer toes are one of the common deformities of the forefoot that can lead to compensatory changes during walking in individuals with this condition. Predicting the adverse effects of tissue damage on the performance of other limbs is very important in the prevention of progressive damage. Finite element (FE) and musculoskeletal modeling can be helpful by allowing such effects to be studied in a way where the internal stresses in the tissue could be investigated. Hence, this study aims to investigate the effects of the hammer toe deformity on the lower extremity, especially on the plantar fascia functions. To compare the joint reactions of the hammer toe foot (HTF) and healthy foot (HF), two musculoskeletal models (MSM) of the feet of a healthy individual and that of a participant with hammer toe foot were developed based on gait analysis. A previously validated 3D finite element model which was constructed using Magnetic Resonance Imaging (MRI) of the diabetic participant with the hammer toe deformity was processed at five different events during the stance phase of gait.It was found that the hammer toe deformity makes dorsiflexion of the toes and the windlass mechanism less effective during walking. Specifically, the FE analysis results showed that plantar fascia (PF) in HTF compared to HF played a less dominant role in load bearing with both medial and lateral parts of PF loaded. Also, the results indicated that the stored elastic energy in PF was less in HTF than the HF, which can indicate a higher metabolic cost during walking. Internal stress distribution shows that the majority of ground reaction forces are transmitted through the lateral metatarsals in hammer toe foot, and the probability of fifth metatarsal fracture and also progressive deformity was subsequently increased. The MSM results showed that the joint reaction forces and moments in the hammer toe foot have deviated from normal, where the metatarsophalangeal joint reactions in the hammer toe were less than the values in the healthy foot. This can indicate a vicious cycle of foot deformity, leading to changes in body weight force transmission line, and deviation of joint reactions and plantar fascia function from normal. These in turn lead to increased internal stress concentration, which in turn lead to further foot deformities. This vicious cycle cause progressive damage and can lead to an increase in the risk of ulceration in the diabetic foot.

Stress distribution of the nonwoven fabric package during winding

During winding, the defects of wrinkles and buckling are prone to occur in the material package. This problem increases in larger packaging sizes, and it is related to the uneven distribution of internal stress. In order to study the internal stress distribution law of the wound package, it is necessary to analyze and discuss the mechanical behavior of the winding process. In this paper, a winding model based on thick-cylinder theory was established to deduce the internal stress distribution and change law of the wound package. Polypropylene melt-blown nonwoven fabric was used as the test material; it showed no deformation residue after the stretching test, which was the allowable stretching conditions of the nonwoven fabric from the low-strain stretching experiment. Based on the material mechanical properties, three tension control strategies (namely constant tension, constant torque, and taper tension winding) of the winding tension were analyzed and the variation law of winding stress was described. The internal stress of the ideal packaging was premised, which was determined as the winding tension control curve under even internal tension and can be called equal internal tension winding. The results showed that the properties of materials have a significant impact on packaging quality, and equal internal tension winding can effectively improve packaging quality and size. The research results can provide some reference for tension control in the winding process of nonwoven fabric or other textile materials.

Ultimate bearing capacity of UHPC-encased CFST medium-long columns under axial compression

To study the effect of slenderness ratio on the ultimate bearing capacity of UHPC (ultra-high performance concrete)-encased concrete-filled steel tubular (UC-CFST) columns, a total of nine specimens under axial compression were tested. A finite element model, considering the effects of material nonlinearity and interaction between concrete and steel tube, was developed. A full-range analysis was carried out on the typical load–deformation response of the UC-CFST medium-long columns. The experimental and numerical analysis results showed that with the increase in slenderness ratio, the deformation capacity of the UC-CFST columns increased, the ultimate bearing capacity decreased, and the failure mode changed from exceedance of compressive strength to loss of stability. The internal forces and stress distributions at the mid-height section of the columns at different stages were also investigated. Finally, the applicability of various methods for calculating the ultimate bearing capacity of ordinary concrete-encased CFST (OC-CFST) columns to UC-CFST columns was discussed, including the methods recommended by Chinese standard T/CECS 188-2019, European standard Eurocode 4, Australian/New Zealand joint standard AS/NZS 5100.6 : 2017 and American standard ANSI/AISC 360-16. It was found that the method recommended in AS/NZS 5100.6 : 2017 was suitable for calculating the ultimate bearing capacity of UC-CFST medium-long columns.

Thermal stresses in growing thermoviscoelastic cylinder and their evolution in the course of selective laser melting processing

AbstractThe evolution in time of stresses and strains in a solid during its manufacturing through selective laser melting (SLM) process is studied. The technological technique, which allows to reduce the intensity of internal and residual stresses is proposed and theoretically grounded. The basic concept of the technique is that the additional in‐homogeneous inductive heating in the vicinity of the growing boundary is applied in the processing. The term “growing boundary” refers to the part of the boundary, on which the additional melted material is added. The evolutionary problem for a growing thermoelastic body is considered in approach of thermal stresses (for slow additive processes) and in the approach of full coupled thermoelasticity (for fast processing). The time depend distribution of internal stresses and their residual counterparts has been evaluated in analytical form in the case of model cylindrical problem. The analysis of the influence of frequency and amplitude of electromagnetic induction as well as viscoelastic properties of heated material being added has been carried out. The results are analyzed numerically for Titanium and Copper materials.

Quasi-honeycomb grain morphologies strengthen passivating layers in Inconel-718 superalloy: Lessons learned from Additive Manufacturing

Despite intensive studies on the effect of grain morphology, that results from Selective Laser Melting (SLM) fabrication, on the mechanical properties of Inconel-718 (IN718), there is a lack of insight into its effect on the alloy’s oxidation behavior. This work compared the oxidation of IN718 alloy manufactured by SLM (SLM-IN718) with that of a cast-rolled commercially available sample (Comm-IN718) when heating to 900 °C. Our results showed that the passivating chromia layer that forms on Comm-IN718 was resistant to rupture because it is a thin oxide layer with a favorable distribution of internal stress due to the alloy’s well-defined quasi-honeycomb primary grain morphology. In comparison, SLM-IN718 had a poorly defined grain morphology, with two distinct length-scales best described as dendritic structures. This morphology resulted in the formation of a thicker oxide scale without favorable distribution of the stress into the alloy grain boundaries. This led to the rupture of the passivation layer. These findings give us a new insight into how the primary microstructure of Inconel alloys can affect their oxidation behavior by affecting the stress distribution within surface oxide films. We are the first to show that morphology characteristic control plays a key role in improving Inconel oxidation.

Theoretical modelling of the effect of thermal delamination of an asphalt concrete pavement from a rigid foundation of a road or bridge

In the practice of road construction, one of the most common phenomena accompanied by delamination, subsequent cracking and destruction of the asphalt concrete pavement on a rigid (cement concrete or metal) base of a road or bridge is the effect of concentration of shear thermal stresses between the pavement and the base in the edge zones of the structure. They are caused by the fact that, as a rule, the coefficients of linear thermal expansion of the phases of the system have different values, which contributes to the occurrence of incompatible shrinkages and expansions in them. Under conditions of frequent temperature changes in heterogeneous asphalt concrete structures with thermomechanical incompatibility of their components, these effects can contribute to their accelerated aging. At the same time, with the thermomechanical compatibility of materials, a more favorable distribution of internal stresses of thermal and mechanical origin is achieved, which excludes premature degradation of the strength of the contacting phases and the entire system as a whole. Using the methods of strengthof materials and the finite element method, it has been established that under the conditions of a change in the temperature of the system during its seasonal and daily fluctuations, shear stresses are subjected to the highest concentration. They are localized in the edge zone of the plane of contact between the layers and increase with an increase in the thickness and modulus of elasticity of the upper layer. These stresses are the main reason for the occurrence of plastic deformations in these zones and subsequent delamination of the structure in them. It is proposed to reduce the concentration and level of generated high-gradient shear stresses by reducing the thickness of the asphalt concrete layer in these areas.

Photodimerization induced hierarchical and asymmetric iontronic micropatterns

Micropatterning various ion-based modality materials offers compelling advantages for functionality enhancement in iontronic pressure sensing, piezoionic mechanoreception, and skin-interfaced electrode adhesion. However, most existing patterning techniques for iontronic materials suffer from low flexibility and limited modulation capability. Herein, we propose a facile and robust method to fabricate hierarchical and asymmetrical iontronic micropatterns (denoted as HAIMs) through programmed regulation of the internal stress distribution and the local ionic migration among an iontronic host. The resultant HAIMs with arbitrarily regulated morphologies and region-dependent ionic electrical performance can be readily made via localized photodimerization of an anthracene-functionalized ionic liquid copolymer (denoted as An-PIL) and subsequent vapor oxidative polymerization of 3,4-ethylenedioxythiophene (EDOT). Based on the piezoionic effect within the resultant distinct doped PEDOT, HAIMs can serve as a scalable iontronic potential generator. Successful syntheses of these fascinating micropatterns may accelerate the development of patterned iontronic materials in a flexible, programmable, and functionally adaptive form.

The Influence of Internal Stress on the Nanocrystal Formation of Amorphous Fe73.8Si13B9.1Cu1Nb3.1 Microwires and Ribbons

The early stages of nanocrystallization in amorphous Fe73.8Si13B9.1Cu1Nb3.1 ribbons and microwires were compared in terms of their internal stress effects. The microstructure was investigated by the X-ray diffraction method. Classical expressions of crystal nucleation and growth were modified for microwires while accounting for the internal stress distribution, in order to justify the XRD data. It was assumed that, due to the strong compressive stresses on the surface part and tensile stresses on the central part, crystallization on the surface part of the microwire proceeded faster than in the central part. The results revealed more rapid nanocrystallization in microwires compared to that in ribbons. During the initial period of annealing, the compressive surface stress of a microwire caused the formation of a predominantly crystallized surface layer. The results obtained open up new possibilities for varying the high-frequency properties of microwires and their application in modern sensorics.

Thermo-mechanical coupling simulation of PBX explosive charge in launch environment

The distribution of internal stress and temperature of explosive charge in launch environment is significant to determine the launch safety of explosive charge. PBX is composed of energetic particles, binder and metal powder, hence it can perform the viscoelastic mechanical properties under the external force loading condition. The generalized Maxwell constitutive equation is embedded into LS-DYNA program by secondary development technology, and verified rational after the simulation of launch safety experiment of explosive charge. Utilizing the constitutive equation and considering the influence of viscoelastic work and chemical reaction heat, the distribution of internal stress and temperature of explosive charge in launch environment is acquired. The results indicate that the viscoelastic work and chemical reaction heat produce a 4K highest temperature rise of explosive charge, the explosive charge without defects cannot ignite in the launch process. The results provide a certain guiding significance for the launch safety evaluation of explosive charge.

Research on an Ultrasonic Longitudinal Critically Refracted Wave Detection Method for the Depth Distribution of Stress

Aluminum alloy components typically have structural characteristics such as large size and complex shape, making the in situ non-destructive detection of internal residual stress in these structures a challenge that the manufacturing sector has tried to solve. Ultrasonic longitudinal critically refracted (LCR) waves have shown good sensitivity to normal stress in the horizontal direction and could be used to detect the distribution of internal residual stress in components, offering an advantage not shared by other detection methods. In this study, we investigated the propagation mode of LCR waves in a 2A14 aluminum alloy component and established the characterization model of the average normal stress of LCR waves in different depth ranges. The blocking effect of LCR waves by a groove with a depth equal to twice the wavelength was analyzed and experimentally verified using a machined aluminum alloy test specimen. Then, the propagation depths of LCR waves in the aluminum alloy at different frequencies were determined. A load test on a cantilever beam based on the stress depth distribution model was designed, and the stress characterization model and LCR waves’ propagation depth were further verified by the self-developed LCR wave stress detection system. The test results showed that the LCR wave could accurately detect the depth distribution of stress and could serve as a useful tool for evaluating the depth distribution of normal stress inside aluminum alloy components.

Numerical analysis of deformation and failure mechanism of metal truss structure of spacecraft under re-entry aerothermodynamic environment

A coupled aerodynamic-thermoelastic analysis framework is presented to simulate the deformation and failure mechanism of metal truss structure for spacecraft under re-entry aerothermodynamic environment. The three-dimensional dynamic thermo-mechanical coupling finite element algorithm is developed by computable modeling of thermoelasticity and heat conduction equations to solve the thermo-mechanical coupling response of spacecraft structure. The aerodynamic boundary conditions such as temperature, heat flux and pressure distribution are computed and provided by the gas-kinetic unified algorithm (GKUA) based on the Boltzmann model equation. The linear interpolation method is designed and employed to couple the boundary conditions between the external aerodynamic flow field and the solid structure interface. The present thermo-mechanical coupling algorithm is verified on one-dimensional and two-dimensional transient heat flow problems with analytical solutions, and good consistency is achieved. The reliability of the GKUA for solving the boundary conditions of the external flow field and the dynamic thermo-mechanical coupling algorithm for the thermal response behavior of the internal structure material is verified by the consistency calculation of the internal and external temperature distribution. Then, the aerodynamic-thermoelastic analysis framework is further applied to study the deformation behavior of structural thermo-mechanical responses and failure mechanism of a hollow sphere under re-entry aerothermodynamic environment with emphasis placed on the influence of flight speeds, rarefied effects, algorithm schemes and damping effects. The numerical results obtained by the presented algorithm including the internal temperature, thermal stress and deformation distribution law are compared and analyzed. Finally, the deformation and failure mechanism of spacecraft structure during re-entry process are solved and analyzed by the present algorithm framework. The displacement, temperature and stress distribution contours of the structure are presented to make a preliminary understanding and evaluation on the process of spacecraft deformation and failure disintegration. It is indicated that the present coupled aerodynamic-thermoelastic analysis framework can provide a reliable, applicable and efficient platform for the thermo-mechanical coupling responses and failure mechanism of re-entry spacecraft in the end of life.

Time-resolved in-situ dislocation density evolution during martensitic transformation by high-energy-XRD experiments: A study of C content and cooling rate effects

Recent studies on the mechanical behavior of martensitic steels show the importance of considering lath martensite as a “polycrystalline aggregate” with a large distribution of the local yield strength. The latter depends on local microstructural features, such as lath sizes, carbon distribution (in solid solution, segregated to defects, in carbides) and density of defects Internal stresses have to be considered as an additional contribution to the yield strength participating to enlarge the distribution of the local flow stresses at the microstructure scale.This study focuses on one of these microstructural features, the dislocation density. Its evolution is followed in situ upon martensitic transformation by High Energy X-Ray Diffraction experiments on a synchrotron beamline. A previously introduced method, based on modified Williamson-Hall (mWH) analysis, is improved in order to take account of the martensite lattice tetragonality when applying the mWH method.The influence of the steel carbon content (0.11, 0.21 and 0.31 wt.% C steels) and of the cooling rate (-10, -50, -100 °C/s and water-quench) on the dislocation density evolution in martensite during the martensitic transformation is established. The effect of the cooling rate on the dislocation density is low between -10 and -100 °C/s, but becomes visible after the water-quench. The higher the C content, the higher the dislocation densities at the end of the transformation. The steels with higher C content also show a wider distribution of the local dislocation density and, therefore, of the local yield strength.

Methodology for Mapping the Residual Stress Field in Serviced Rails Using LCR Waves.

Non-destructive stress measurement by ultrasonic testing is based on calculating the acoustoelastic modulus obtained from the relationship between material stress and sound wave velocity. A critically refracted longitudinal (LCR) wave, which is a bulk longitudinal wave penetrating below and parallel to the surface below an effective depth, is most suitable for ultrasonic stress measurement tests because it exhibits a relatively large change in travel time in response to a change in stress. In particular, the residual stress distribution through the thickness of the subject can be calculated if transducers of different frequencies are applied because of the characteristic of propagation to different depths of penetration depending on the frequency. The main purpose of this study was to visualize the internal or residual stress distribution through the thickness of rails using LCR waves. To this end, LCR probes with different center frequencies were designed and manufactured, and the residual stress values ​​of an unused railroad rail and two used railroad rails operated under different conditions were calculated. This was done using the ultrasonic signals received from each probe, of which the distributions were mapped. Through these mapping results, different residual stress values could be calculated according to the depth. The differences in residual stress generation and distribution according to the conditions surrounding the contact between train wheels and rails, and their characteristics, were visualized and analyzed. As a result, it could be concluded that the non-destructive evaluation technique using LCR waves could detect differences in the residual stress of a rail, and thus can be used to measure the residual stress of the rail accurately.

Application of a New Alloy and Post Processing Procedures for Laser Cladding Repairs on Hypereutectoid Rail Components.

The development of a laser cladding repair strategy is critical for the continued growth of heavy-haul railway networks. Premium hypereutectoid rails have undergone laser cladding using a new martensitic stainless-steel alloy, 415SS, developed for high carbon rails after standard cladding metals were found to be incompatible. Non-destructive neutron diffraction techniques were used to measure the residual stress in different layers generated across a dissimilar metal joint during laser cladding. The internal stress distribution across the cladding, heat-affected zone (HAZ), and substrate was measured in the untempered rail, after 350 °C and 540 °C heat treatment procedures and two surface grinding operations. The martensitic 415SS depositions produce compressive stress in the cladding, regardless of tempering procedures, which may inhibit fatigue crack propagation whilst grinding operations locally relive surface stress. Balancing tensile stresses were recorded below the fusion boundary in the HAZ due to thermal gradients altering the microstructure. The combination of 540 °C tempering and 0.5 mm surface layer removal produced a desirable combination of compression in the cladding deposition with significantly reduced tensile stresses in the HAZ. A comparison with the current literature shows that this alloy achieves a unique combination of desirable hardness, low tensile stress, and compression in the cladding layer. Data obtained during strain scanning has been used to determine the location of microstructural changes at the fusion boundary and HAZ through correlation of the stress, strain, full width at half maximum (FWHM), and intensity profiles. Therefore, neutron diffraction can be used for both the accurate measurement of internal residual stress and to obtain microstructural information of a metallurgical join non-destructively.

Numerical Analysis of Microcrack Propagation Characteristics and Influencing Factors of Serrated Structural Plane.

The serrated structural plane is the basic unit of structural plane morphology. However, the understanding of its internal stress distribution, failure mode and crack evolution law was not clear enough in previous studies. In this paper, the shear mechanical properties of the serrated structural planes were studied by numerical simulation, and the crack evolution law of the serrated structural planes and the effects of four microscopic parameters on the shear properties were analyzed. The results show that: (1) the number of microcracks increases with the increase in normal stress; the crack expansion rate is slow before the shear stress reaches the peak. After the shear stress reaches the peak, the crack expansion rate continues to increase, and the microcracks keep sprouting and expanding, and the number of microcracks tends to stabilize when the shear stress reaches the residual shear strength. (2) The particle contact stiffness ratio and parallel bond stiffness ratio were negatively correlated with the shear strength; and the particle contact modulus and parallel bond modulus were positively correlated with the shear strength. As the particle contact modulus and parallel bond modulus increase, the peak shear displacement gradually decreases. The parallel bond stiffness ratio has a negative correlation with the peak shear displacement. This study is expected to provide theoretical guidance for the microscopic parameter calibration and shear mechanical analysis of serrated structural planes. (3) Several XGBoost, WOA-XGBoost, and PSO-XGBoost algorithms are introduced to construct the quantitative prediction model, and the comparative analysis found that WOA-XGBoost has the best fitting effect and can be used for the prediction of shear strength. When using this model to calculate the weight shares of micro-parameters, it was found that has the greatest influence on shear strength, followed by ; and had the least influence.

Finite Element Analysis of Proximal Femur Bionic Nail (PFBN) Compared with Proximal Femoral Nail Antirotation and InterTan in Treatment of Intertrochanteric Fractures.

ObjectiveTo compare the biomechanical properties of proximal femur bionic nail (PFBN), proximal femoral nail antirotation (PFNA) and InterTan in the treatment of elderly intertrochanteric fractures AO/OTA 31‐A1.3 by finite element analysis.MethodsWe used Mimics, Unigraphics and other software to establish normal femur and AO/OTA 31‐A1.3 fracture models, and reconstructed PFBN, PFNA and InterTan intramedullary nail models, and assembled them on the fracture model. The ANSYS software was used to compare the femoral von Mises stress distribution, deformation distribution, and internal fixation stress distribution of each group under a load of 2100 N.ResultsIt could be seen that the femoral maximum stress, femoral maximum displacement, and maximum stress of internal fixation of the PFBN group were lower than those in the PFNA group and the InterTan group. The maximum femoral stress of the PFBN was 190.25 MPa, while the maximum stress of the femur of the PFNA and InterTan groups were 238.41 Mpa and 226.97 Mpa. The maximum femoral displacement of each group were located at the top of the femoral head, and the maximum displacement of the PFBN group was 14.373 mm, and the maximum displacement values of the PFNA and InterTan groups were 19.49 and 15.225 mm. For the stress distribution of intramedullary nail, the maximum stress of the three kinds of internal fixation was located on the main nail. The maximum stress of PFBN was 1191.8 MPa, compared with 2142.8 MPa for PFNA and 1702.3 MPa for InterTan. And the maximum stress on the PFBN pressure nail was 345.35 MPa, compared with 868.6 MPa for the PFNA spiral blade and 545.5 MPa for InterTan interlocking twin nails.ConclusionCompared with PFNA and InterTan, PFBN has better mechanical properties. The biomechanical characteristics of PFBN are more advantageous than PFNA and InterTan internal fixation system in the treatment of femoral intertrochanteric fractures.

An investigation of the ankle contact forces in a foot with hammer toe deformity. A comparison of patient-specific approaches using finite element modeling and musculoskeletal simulation.

The internal forces and stresses in the tissue are important as they are linked to the risk of mechanical trauma and injuries. Despite their value, the internal stresses and forces cannot be directly measured in-vivo. A previously validated 3D finite element model (FEM) was constructed using Magnetic Resonance Imaging (MRI) of a person with diabetes and hammer toe deformity. The foot model simulated at five different instances during the stance phase of gait. The internal stress distribution on the talus that was obtained using the FEM simulation, was used to calculate the joint reaction force at the ankle joint. In addition, the musculoskeletal model (MSM) of the participant with hammer toe foot was developed based on the gait analysis and was used to determine the muscle forces and joint reactions. The result showed that the vertical reaction forces obtained from the FEM and MSM follow a similar trend through the stance phase of gait cycle and are significantly correlated ( R=0.99). The joint reaction forces obtained through the two methods do not differ for the first 25% of the gait cycle, while the maximum difference was ∼0.7 Body weight that was observed at 50% of the stance phase. Clinical Relevance: Finite element modeling and musculoskeletal simulation can shed light on the internal forces at the ankle in pathological conditions such as hammer toe. The similarities and differences observed in the joint reaction forces calculated from the two methods can have implications in assessing the effect of clinical interventions.

Atomistic simulations of plasticity heterogeneity in gradient nano-grained FCC metals

Gradient nano-grained (GNG) metals are an emerging class of heterogeneous structured materials with extraordinary mechanistic performances and unnatural plasticity mechanisms, which are unachievable in traditional homogeneous nanocrystalline metals. Here, the metal-type dependence of plasticity heterogeneity in three different GNG structured face-centered cubic metals (Cu, Al, and Ni) was studied by large-scale molecular dynamic simulation, from the plasticity perspectives of dislocation, deformation twinning, and grain boundary. The dislocation and deformation twinning-based plasticity show clear metal-type dependence in GNG structure of different selected metals. For grain boundary activities, GNG Cu, Al, and Ni are equivalent in the overall grain boundary sliding and micro-cracking process, but GNG Ni shows markedly higher grain boundary-free volume and stress due to the grain size gradient-induced twinning process. We find that GNG Al shows the weakest plasticity-induced strengthening due to its weak strain hardening and ununiform plasticity distribution, while GNG Ni exhibits prominent plasticity-induced strengthening because of the obvious extra strain hardening process. Besides, we further discuss the reasonable gradient grain size design of GNG structure, considering the critical grain size of the Hall-Petch and inverse Hall-Petch relationship, which will demonstrate a more rational plasticity heterogeneity and internal atomic stress distribution.