Material Stiffness Research Articles

A safe highway driving surface using an open graded concrete

The use of open graded asphalt wearing courses for highway pavements is established technology but the opportunities to use a similar cementitious based surface type for concrete pavements has been limited by material technology and costs. Research work in Australia using new generation concrete admixtures and placing techniques has allowed the development of a stiff open graded concrete material which is expected to meet the same durability requirements as standard concrete used for roads. Open graded surface textures absorb vehicle noise emissions and minimise water film build up during rain events, thus resulting in quieter and safer driving conditions. In addition, experience has shown that high void contents at the surface assist in reducing water spray generation and glare reflection. This paper details the laboratory mix design test results and three field trials using hand placing techniques to define the limits of a wet on wet construction process. The concrete developed in the trials has high compressive strength with air voids in the range of 20 to 30%. The next stage in the development of the open graded concrete surfacing is in the use of mechanised placing techniques to make the process cost competitive to other asphalt alternatives and yet meet similar durability requirements to those currently expected in concrete road surfaces over a 40 year design period.

Theoretical exploration of structural, electronic, elastic, mechanical, and thermodynamic properties of MgCu4Sn intermetallic compound for engineering applications: first-principle calculations

The structural, electronic, elastic, mechanical, and thermodynamic properties of the ternary MgCu4Sn intermetallic compound are investigated by means of first-principle calculations within density functional theory (DFT), in combination with the quasi-harmonic Debye model. Local density approximation (LDA) and generalized gradient approximation (GGA) are made for electronic exchange–correlation potential energy. The lattice constant is in good agreement with experimental data. We determine the elastic constant tensor of the compound from the calculated stress–strain relation in both approximations. Once the elastic constants are obtained, the bulk modulus B, shear modulus G, Young’s modulus E, Poisson’s ratio ν, anisotropy factor A, and the ratio B/G for MgCu4Sn compound were deduced using Voigt-Reuss-Hill (VRH) approximation. The ground-state structure of MgCu4Sn is predicted to be thermodynamically and mechanically stable. The obtained band structure and density of states reveal metallic character of MgCu4Sn. The calculation results show also that this intermetallic crystal is a stiff, elastically anisotropic and ductile material. The Debye temperature is also determined from elastic constants. The temperature dependence of the constant volume heat capacity Cv, the entropy S, and the volumetric thermal expansion coefficient α in a quasi-harmonic approximation have been obtained from calculated energy E as a function of the volume V of a MgCu4Sn crystal and discussed for the first report.

Dynamic and quasi-static evaluation of stiffness properties of CLT: longitudinal MoE and effective rolling shear modulus

Cross-Laminated Timber (CLT) is an engineered wood product composed of solid layers of glued sawn timber. In this study, essential material stiffness parameters for CLT made from Norway spruce and Scots pine are evaluated. Specifically, the longitudinal modulus of elasticity (MoE) for longitudinally oriented layers and the effective rolling shear modulus for transversely oriented layers are the focus. By combining finite element (FE) analysis with four-point, out-of-plane bending tests using digital image correlation (DIC), a robust assessment of the effective rolling shear modulus of CLT layers is achieved. Additionally, eigenvalue analysis, applied to an FE model, along with resonance frequencies obtained from dynamic excitation of CLT, enables stable and simultaneous assessment of the dynamic longitudinal MoE and effective rolling shear modulus. Notably, while the dynamic MoE of longitudinal CLT layers is only 4% higher than the quasi-static local MoE, the dynamic effective rolling shear modulus of CLT layers is 40% higher than the quasi-static effective rolling shear modulus. This finding indicates a tangible viscoelastic behavior of wood concerning rolling shear.

Modulated Mechanical Properties of Epoxy-Based Hybrid Composites via Layer-by-Layer Assembly: An Experimental and Numerical Study.

In this study, epoxy-based composites were fabricated using a layer-by-layer assembly technique, and their mechanical properties were systematically evaluated. The inclusion of cellulose nanocrystals led to variations in the mechanical properties of the composites. These modified properties were assessed through tensile and flexural tests, with each layer cast to enhance strength. Due to the inherent characteristics of epoxy, a single specimen was fabricated through chemical bonding, even post-curing. This approach demonstrated that a three-layer structure, developed using the layer-by-layer method, exhibited improved elastic and flexural moduli compared to a single-layer composite. This improvement aligns with theoretical predictions, which suggest that stiffness increases when stiffer materials are positioned farther from the neutral axis in a layered structure. Furthermore, numerical analysis validated changes in stress distribution across each layer. Consequently, this method enables the production of composites with superior mechanical properties while minimizing the quantity of cellulose nanocrystals required.

High-cycle constitutive model for traffic loading induced permanent deformation of coarse granular material incorporating particle breakage

Repetitive traffic loads lead to particle rearrangement and breakage in granular trackbeds and roadbeds, resulting in irreversible deformations that hinder normal operations. Current models for cyclic deformation primarily address rearrangement-induced densification but overlook breakage-induced degradation. Particle breakage disrupts inter-particle interlocking and creates finer debris, promoting volumetric contraction and weakening material stiffness. This study introduces a novel cyclic constitutive model for predicting plastic strain in coarse granular materials, emphasizing the role of particle breakage. The model extended from the existing cyclic densification model within the framework named “Shakedown Surface Model,” integrating breakage effects. In this model, shakedown thresholds, associated with material densification, are influenced by both plastic strain and breakage-induced loosening. Parameters for shearing contraction/dilation decrease with breakage, capturing breakage-induced contraction. Additionally, the model accounts for principal stress rotation from moving loads. Implemented via an implicit Euler-backward algorithm and Newton–Raphson iteration, the model was validated against cyclic triaxial and full-scale tests, demonstrating its accuracy in predicting the permanent deformation of granular materials under traffic loads.

Experimental and Simulation Study of Fiber-Reinforced and Foam-Filled Anti-Tetrachiral and Re-Entrant Structures

Negative Poisson's ratio (NPR) honeycomb structures exhibit excellent properties, including resistance to compression, impact resistance, and energy absorption. To better study NPR structures, this paper investigates two different methods to enhance the mechanical properties and energy absorption capabilities of re-entrant and anti-tetrachiral structures. Fiber reinforcement and polyurethane foam filling methods were used to address material enhancement issues, aiming to achieve structures with high specific energy absorption. Specimens were fabricated using selective laser sintering technology with three different material configurations, including nylon, carbon fiber-reinforced nylon, and glass fiber-reinforced nylon, for comparative analysis. Additionally, some specimens were filled with polyurethane foam. Quasi-static compression tests indicated that both glass fiber-reinforced nylon and foam filling significantly improved the material's stiffness, initial peak stress, and energy absorption, although the "auxetic" effect was weakened. Subsequently, finite element methods were used to simulate the deformation modes and mechanical properties of the two NPR structures, and the numerical results showed good agreement with the experimental results.

Theoretical study on cushioning isolation of cushioning materials under high impact environment

Abstract When a guided projectile is fired, the high pressure generated by the gunpowder gas in the chamber produces a high impact accelerationtion of tens of thousands of g in the projectile, which makes the components of the projectile withstand a huge instantaneous, high-energy, strong impact load. This type of load environment results in the deformation and failure of the internal mechanical components and damage to all type of electronic components, leading to the failure of the entire guidance system. Selecting the correct cushioning material to store or dissipate the impact energy so that the peak value of the impact is reduced to within the range of the component is an effective technical approach. In this study, a mechanical analysis of the high-impact environment inside a certain type of artillery was conducted, the axial impact load variation curve of the projectile was drawn, and a buffer isolation mathematical model of the buffer material was established. This study obtains a visualization interface using Python to solve the problem intelligently. By changing the mass of the mass block, the equivalent stiffness and equivalent damping value of the buffer material in the program script, researchers can use the model to draw the displacement, velocity and acceleration curves of the two mass blocks simultaneously, and obtain the buffering efficiency and maximum elastic compression of the buffer material. When the equivalent damping is the same, this study obtains the buffering efficiency curve and the maximum elastic compression curve changing with the stiffness through the model. The results show that in the range of elastic deformation, the buffering efficiency is proportional to the length of the buffering material but inversely proportional to its equivalent stiffness. The mathematical model provides technical support for the pre-design development and selection of cushioning materials in high-impact environments.

Optomagnetic Micromirror Arrays for Mapping Large Area Stiffness Distributions of Biomimetic Materials.

A new device termed "Optomagnetic Micromirror Arrays" (OMA) is demonstrated capable of mapping the stiffness distribution of biomimetic materials across a 5.1 mm × 7.2 mm field of view with cellular resolution. The OMA device comprises an array of 50000 magnetic micromirrors with optical grating structures embedded beneath an elastic PDMS film, with biomimetic materials affixed on top. Illumination of a broadband white light beam onto these micromirrors results in the reflection of microscale rainbow light rays on each micromirror. When a magnetic field is applied, it causes each micromirror to tilt differently depending on the local stiffness of the biomimetic materials. Through imaging these micromirrors with low N.A. optics, a specific narrow band of reflection light rays from each micromirror is captured. Changing a micromirror's tilt angle also alters the color spectrum it reflects back to the imaging system and the color of the micromirror image it represents. As a result, OMA can infer the local stiffness of the biomimetic materials through the color change detected on each micromirror. OMA offers the potential for high-throughput stiffness mapping at the tissue-level while maintaining spatial resolution at the cellular level.

Full-scale testing of the bending behaviour of UHPC gravity–grouted sleeve–prestressed anchor joints in assembled frame tunnels

Joints are weak parts of a prefabricated tunnel frame structure, and the joint bending resistance is a crucial factor controlling the deformation and bearing capacity of the tunnel structure. This study proposes a joint with dual anchoring of ultrahigh-performance concrete (UHPC) gravity-grouted sleeves and prestressed steel reinforcement. The material strain, opening amount, deflection, stiffness, rotation angle, and ductility of the grouted sleeve-prestressed anchor joint and grouted sleeve joint types were tested through full-scale bending tests, and with the help of noncontact measurement technology, the dynamic process of joint bending failure was captured. In addition, the finite element method was used to evaluate the influences of the prestressing magnitude and the position of the prestressing steel bars on the bending performance of the joints. The results indicate that the damage to the concrete joint is concentrated at the top of the joint and that the damage to the top surface of the joint concrete is distributed in an "m" shape, indicating that the damage to the concrete in the area between the sleeves is relatively extensive. Owing to the turning points of grouting material fracture and joint yielding, the joint bending process clearly shows a three-stage change. The grouted sleeve-prestressed anchor joint exhibits greater bending stiffness than the grouted sleeve joint does, and the slip of the grouting material in the sleeve is smaller. Prestressing plays a key role in enhancing the bearing capacity and deformation resistance of joints. Considering that excessive prestress may increase the brittleness of the structure, the prestress of the grouted sleeve-prestressed anchor joint is recommended to be 200 kN. This study provides theoretical support for the joint design and construction of prefabricated tunnel frame structures.

Mechanical Characterization and Modeling of Large-Format Lithium-Ion Battery Cell Electrodes and Separators for Real Operating Scenarios

This study presents a novel application-oriented approach to the mechanical characterization and subsequent modeling of porous electrodes and separators in lithium-ion cells to gain a better understanding of their real mechanical operating behavior. An experimental study was conducted on the non-linear stiffness of LiNi0.8Co0.15Al0.05O2 and graphite electrodes as well as PE separators, harvested from large-format lithium-ion cells, using compression tests. The mechanical response of the components was determined for different operating conditions, including nominal stress levels, mechanical loading rates, and mechanical cycles. The presented work describes the test procedure, the experimental setup, and an objective evaluation method, allowing for a detailed summary of the observed mechanical behavior. A distinct nominal stress level and mechanical cycle dependency of the non-linear stiffnesses of the porous materials were found. However, no clear dependency on compression rate was observed. Based on the experimental data, a poroelastic mechanical model was utilized to predict the non-linear behavior of these porous materials under real mechanical operating scenarios with a normalized root-mean-squared error less than 5.5%. The results provide essential new insights into the mechanical behavior of porous electrodes and separators in lithium-ion cells under real operating conditions, enabling the accelerated development of high-performing and safe batteries for various applications.

Modeling the uncoupled damage-healing behavior of self-healing cementitious material with phase-field method

Microbially induced calcium carbonate precipitation (MICCP), which is a biochemical process that utilizes microbial metabolic activities to precipitate calcium carbonate (CaCO3), can be used for damage-healing by crack filling. In this study, we develop a computational uncoupled damage-healing framework to predict the fracture-healing responses of a multifunctional vascular cementitious composite. The framework integrates a numerical model for self-healing cementitious materials into a damage model developed using the phase-field method. The numerical self-healing model is utilized to compute the MICCP healing based on the chemical and enzyme kinetics of microbial activities, whereas the phase-field method is employed to capture the fracture response of the material. This framework comprises three stages: Stage 1 — damage, Stage 2 — healing, and Stage 3 — re-damage of the healed structure. Using the load–displacement curve of the healed cementitious material, the proposed uncoupled damage-healing models can predict the recovery in material stiffness based on the MICCP healing ratio. The proposed framework can be applied to a wide range of fracture-healing problems in self-healing cementitious materials utilizing MICCP as the healing methodology.

Comparative Study of Core Material's Stiffness on Sandwich Panel with Composite Face Sheets beyond the Yield Point

The study aims to explore the performance of sandwich panels exceeding the yield point stiffness of the core material. Sandwich panels have gained growing attention among designers owing to their excellent corrosion properties, and lightweight, and speedy installation process. They have been applied in numerous industrial sectors, including aerospace, architectural, marine, and transportation. Typically, sandwich panels are composed of a single central core sandwiched by a pair of outer face sheets, where the core is normally developed using softer materials compared to the face sheets. Given that past studies have primarily focused on sandwich panels in the elastic range, this present study explored the performance of sandwich panels exceeding the yield point stiffness of the core material. The univariate search optimization method was utilized to assess the elastic modulus ratio of the core (typically foam) to the face sheet (composite material). The load was elevated in a quasi-static order until the face sheets reached their yield point. Subsequently, the panel was simulated using the finite element analysis commercial package ANSYS APDL, with simply supported boundary conditions used on all sides of the panel. The proposed model was verified by comparing the numerical and experimental data from recent literature. Based on the results, the panel's increased load-carrying capacity corresponded as the core material stiffness exceeded its yield limit. Moreover, the transmission of load to the face sheets increased as the core stiffness decreased. In summary, stiffer core materials caused the sandwich panel to behave more as isotopic face sheets. Thus, the face sheets yielded ahead of the core material.

3D modeling of geological and geotechnical soil characterization using the MASW method: a case study in Southern Ijaw LGA, Bayelsa State, Nigeria

The Multichannel Analysis of Surface Waves (MASW) method is pivotal for non-invasive subsurface shear-wave velocity profiling, essential in geotechnical and seismic investigations. This study aims to model 3D geological and geotechnical soil properties in Southern Ijaw LGA, Bayelsa State, Nigeria, using MASW. Data collection involved a 12-channel ABEM Terraloc Mark 6 geophone system across nine survey points, complemented by nine boreholes and two Standard Penetration Tests (SPT). The results revealed a shallow water table at 0.30 m and stratified soil layers ranging from soft clay to dense sands. SPT N-values increased from 2 at the surface to 34 in deeper layers, reflecting significant soil densification and strength enhancement with depth. Shear-wave velocities (Vs) ranged from 207.11 m/s at 0.87 m to 502.88 m/s at 12.9 m, while compressional-wave velocities (Vp) ranged from 388.71 m/s to 948.98 m/s. Poisson’s ratio was consistent at 0.3, indicating uniform properties across depths. Elastic moduli, including shear modulus (µ), Young’s modulus (E), and bulk modulus (K), increased with depth, indicating greater subsurface material stiffness. The 3D geological model delineated distinct layers: soft clay (0–3.61 m), silty clay (3.61–5.19 m), fine sand (5.19–6.55 m), medium sand (6.55–7.83 m), and medium to coarse sand (7.83–12.90 m). These layers exhibited specific engineering properties, with ultimate and allowable bearing capacities ranging from 154.38 to 543.71 kPa and 51.46 to 181.24 kPa, respectively. MASW-derived N-values showed a strong correlation with traditional SPT N-values (R2 = 0.9401) and shear wave velocities and SPT N-values (R2 = 0.961), confirming MASW's effectiveness for geotechnical characterization. This enhances the precision of 3D soil models and supports more informed engineering decisions. The findings significantly advance the understanding of soil properties, particularly in the Niger Delta.

Deformation and fracture of lithosphere-inspired polymeric multi-layer composites

Inspired by the diversity of structures and patterns inherent in the earth's lithosphere, this study endeavors to enhance the interplay between stiffness and toughness through the introduction of a new class of polymeric multi-layer composite materials termed by the authors as "lithomers". Structured single-edge notched bending specimens were fabricated using a combination of additive manufacturing and casting, employing two different methacrylate-thiol resins. The outer layers exhibit a stiff and brittle characteristic, while the layer in between is compliant in nature. Three types of lithomers with wave-like structures and one with a rectilinear structure were investigated regarding their stiffness and toughness in a 3-point bending setup. The results were compared with those of a pure stiff matrix material. The findings revealed that fracture toughness increased regardless of the interlayer's shape compared to the pure matrix material. Correspondingly, this enhancement in fracture toughness correlated with a reduction in stiffness. The most balanced results in terms of stiffness and fracture toughness were achieved, with the lithomer having a wave-like structure in its initial stage. It exhibited a roughly 27 times improvement in fracture toughness with a moderate decrease in stiffness of approx. 1/5 compared to the pure matrix material.

Biomechanical testing of virtual meniscus implants made from a bi-phasic silk fibroin-based hydrogel and polyurethane via finite element analysis

ObjectiveTo investigate the suitability of different material compositions and structural designs for 3D-printed meniscus implants using finite element analysis (FEA) to improve joint function after meniscal injury and guide future implant development. DesignThis experimental study involved in-silico testing of a meniscus model developed from two materials: a specially formulated hydrogel composed of silk fibroin (SF), gelatine, and decellularized meniscus-derived extracellular matrix (MD-dECM), and polyurethane (PU) with stiffness levels of 54 and 205 MPa. Both single-material implants and a two-volumetric meniscus model with an SF/gelatine/MD-dECM core and a PU shell were analysed using FEA to simulate the biomechanical performance under physiological conditions. ResultsThe hydrogel alone was found to be unsuitable for long-term use due to instability in material properties beyond two weeks. PU 54 closely replicated the biomechanical properties of an intact meniscus, particularly in terms of contact pressure and stress distribution. However, hybrid implants combining PU 54 with hydrogel showed potential but required further optimization to reduce stress peaks. In contrast, implants with a PU 205 shell generated higher induced stresses, increasing the risk of material failure. ConclusionsFEA proves to be a valuable tool in the design and optimization of meniscal implants. The findings suggest that softer PU 54 is a promising material for mimicking natural meniscus properties, while stiffer materials may require design modifications to mitigate stress concentrations. These insights are crucial for refining implant designs and selecting appropriate material combinations before physical prototype production, potentially reducing costs, time, and the risk of implant failure.

Trading Spaces: Adaptive Subspace Time Integration for Contacting Elastodynamics

We construct a subspace simulator that adaptively balances solution improvement against system size. The core components of our simulator are an adaptive subspace oracle, model, and parallel time-step solver algorithm. Our in-time-step adaptivity oracle continually assesses subspace solution quality and candidate update proposals while accounting for temporal variations in deformation and spatial variations in material. In turn our adaptivity model is subspace agnostic. It allows application across subspace representations and expresses unrestricted deformations independent of subspace choice. We couple our oracle and model with a custom-constructed parallel time-step solver for our enriched systems that exposes a pair of user tolerances which provide controllable simulation quality. As tolerances are tightened our model converges to full-space solutions (with expected cost increases). On the other hand, as tolerances are relaxed we obtain output-bound simulation costs. We demonstrate the efficacy of our approach across a wide range of challenging nonlinear materials models, material stiffnesses, heterogeneities, dynamic behaviors, and frictionally contacting conditions, obtaining scalable and efficient simulations of complex elastodynamic scenarios.

Technical note: Effects of chemical, radiational, and thermal treatment of bone tissue on material stiffness and anisotropy.

We evaluated the effects of four treatments for pathogen inactivation in bone tissue (ethanol storage, formalin fixation, gamma irradiation, and heat treatment via autoclave) on stiffness and anisotropy values in bone samples. Cortical bone samples from the humerus of 14 bovine specimens were subjected to Knoop microindentation analysis in longitudinal and transverse planes of section and four indenter orientations within each section. From each specimen, individual samples were assigned to one of five treatment conditions: 50% ethanol-saline solution, formalin immersion, gamma irradiation, autoclave, and buffered saline (controls). Each sample was microindented 100 times and the resultant stiffness data were analyzed by a resampled factorial analysis of variance. First- and second-order interactions as well as main effects of indenter orientation, treatment, and specimen were significant for both transverse and longitudinal sections, with the sole exception that indenter orientation-treatment interaction was nonsignificant for longitudinal sections. Interaction plots reveal that thermal (autoclave) and irradiation treatments depress stiffness values the most, while patterns of indentation stiffness at different orientations are relatively unaffected. Patterns of anisotropy are relatively unaffected by preservative treatments, but elastic modulus changes are consistent and unambiguous. Formalin and ethanol treatments are most comparable to controls and represent the best preservative media for mechanical testing. These options, however, are likely to be the least effective for ensuring the inactivation or sterilization of potentially contaminated samples.

An Alternative Finite Element Formulation to Predict Ductile Fracture in Highly Deformable Materials

Abstract An alternative finite element formulation to predict ductile damage and fracture in highly deformable materials is presented. For this purpose, a finite-strain elastoplastic model based on the Gurson–Tvergaard–Needleman (GTN) formulation is employed, in which the level of damage is described by the void volume fraction (or porosity). The model accounts for large strains, associative plasticity, and isotropic hardening, as well as void nucleation, coalescence, and material failure. To avoid severe damage localization, a nonlocal enrichment is adopted, resulting in a mixed finite element whose degrees-of-freedom are the current positions and nonlocal porosity at the nodes. In this work, 2D triangular elements of linear-order and plane-stress conditions are used. Two systems of equations have to be solved: the global variables system, involving the degrees-of-freedom; and the internal variables system, including the damage and plastic variables. To this end, a new numerical strategy has been developed, in which the change in material stiffness due to the evolution of internal variables is embedded in the consistent tangent operator regarding the global system. The performance of the proposed formulation is assessed by three numerical examples involving large elastoplastic strains and ductile fracture. Results confirm that the present formulation is capable of reproducing fracture initiation and evolution, as well as necking instability. Convergence analysis is also performed to evaluate the effect of mesh refinement on the mechanical response. In addition, it is demonstrated that the nonlocal parameter alleviates damage localization, providing smoother porosity fields.

Effects of TiB2 particles on deformation behavior, softening mechanisms and recrystallization texture of hot-compressed Fe-TiB2 composites

Fe-TiB2 composites, also termed as high modulus steel, offer a promising solution to the challenge of achieving both lightweight and high stiffness materials. However, the presence of TiB2 particles in Fe-TiB2 composites results in poor hot-workability. Therefore, understanding the effects of TiB2 particles on the hot deformation behavior, dynamic recrystallization (DRX), and microstructural evolution of Fe-TiB2 composites is crucial for optimizing their hot-working process. In this study, we elucidated the effects of TiB2 particles on deformation behavior and dynamic softening behavior by conducting a series of isothermal compression tests on as-cast Fe-TiB2 composites and as-cast base alloys (control group) at temperatures of 800–1200 °C, strains of 0.36–1.2, and strain rates of 0.01–1 s⁻1. Using electron backscatter diffraction, we characterized the microstructures of composites and base alloys, showing that TiB2 particles induce a DRX process through particle stimulated nucleation (PSN) at low temperatures and promote continuous dynamic recrystallization (CDRX) at high temperatures. The presence of TiB2 particles have significantly affected the dislocation movement and distribution, which changes the deformation energy distribution and thus facilitates different DRX behaviors under various thermal deformation conditions. Additionally, the microstructure resulting from DRX through PSN exhibits significant texture weakening and grain refinement, presenting a promising method for fabricating ultrafine-grained Fe-TiB2 composites.

Stress transfer, stress shielding, and micro-movement analysis in the implanted femoral bone region with various hip implant stem materials: A finite element approach

This study analyzed stress distribution, stress shielding, and bone-implant micro-movement in different zones of implanted femur bone with hip implant stem composed of Inconel 718, 316 Stainless Steel, Ti6Al7Nb, Ti6Al4V, Ti-12Mo-6Zr-2Fe (TMZF), and Ti-33.6Nb-4Sn (TNS). A 3D model of a well-known hip implant “Summit” was placed in an adult femur bone, and walking, standing, and sitting load instances were simulated as loading conditions using a finite element model (FEM). Increasing hip implant material hardness reduces stress transmission in Gruen zones 1, 2, 5, 6, and 7 of the femur bone but increases it in Gruen zone 3. It also enhances stress holding in the implant body and elevates the Stress Shielding index (SSI) at seven Gruen zones for each load case. TMZF and TNS hip implants in Gruen regions 5 and 7 have negative SSI values under specific load situations. At Gruen zone 3, the SSI value for implanted femur bone increases with hip implant stem material hardness and SSI values variation in implanted femur bone for hip implants of different materials at each loading scenario is within 10 %. However, SSI values do not follow any pattern at Gruen zone 4. Less stiff implants demonstrate increased micro-movement in bone-implant contact zones 1 and 2 due to higher stress transfer. Greater stiffness (Inconel 718) caused increased micro-movement in zone 3. The overall micro-movement (Average and Highest) found an incremental trend for hip implants with higher stiffer to lower materials. The lowest stiff material (TNS) is still in the range of appropriate micro-movement for promoting osteointegration. This study indicates TNS as a superior hip implant stem material due to its enhanced stress distribution, reduced SSI in most Gruen zones on implanted femurs, and acceptable micro-movement over the implant-bone area of contact. Further research is required to minimize SSI, particularly in Gruen zone 3, when using TNS hip implants in implanted femurs.