Image-based Finite Element Research Articles

Intravascular Ultrasound Image-Based Finite Element Modeling Approach for Quantifying In Vivo Mechanical Properties of Human Coronary Artery.

Quantifying the mechanical properties of coronary arterial walls could provide meaningful information for the diagnosis, management, and treatment of coronary artery diseases. Since patient-specific coronary samples are not available for patients requiring continuous monitoring, direct experimental testing of vessel material properties becomes impossible. Current coronary models typically use material parameters from available literature, leading to significant mechanical stress/strain calculation errors. Here, we would introduce a finite element model-based updating approach (FEMBUA) to quantify patient-specific in vivo material properties of coronary arteries based on medical images. In vivo cine intravascular ultrasound (IVUS) and virtual histology (VH)-IVUS images of coronary arteries were acquired from a patient with coronary artery disease. Cine IVUS images showing the vascular movement over one cardiac cycle were segmented, and two IVUS frames with maximum and minimum lumen circumferences were selected to represent the coronary geometry under systolic and diastolic pressure conditions, respectively. VH-IVUS image was also segmented to obtain the vessel contours, and a layer thickness of 0.05 cm was added to the VH-IVUS contours to reconstruct the coronary geometry. A computational finite element model was created with an anisotropic Mooney-Rivlin material model used to describe the vessel's mechanical properties and pulsatile blood pressure conditions prescribed to the coronary luminal surface to make it contract and expand. Then, an iterative updating approach was employed to determine the material parameters of the anisotropic Mooney-Rivlin model by matching minimum and maximum lumen circumferences from the computational finite element model with those from cine IVUS images. This image-based finite element model-based updating approach could be successfully extended to determine the material properties of arterial walls in various vascular beds and holds the potential for risk assessment of cardiovascular diseases.

Image-based finite element model stiffness and vBMD by single and dual energy CT reconstruction kernel

Single-energy quantitative computed tomography (SEQCT) provides volumetric bone mineral density (vBMD) measures for bone analysis and input to image-based finite element models (FEMs). Dual-energy CT (DECT) improves vBMD by accounting for voxel-specific material variations utilizing scans at multiple x-ray energies. vBMD is also altered by reconstruction kernel that cannot be accounted for using calibration phantoms. This study compared vBMD and FEM stiffness derived from SEQCT and DECT images reconstructed with two common kernels. SEQCT and DECT images of cadaveric shoulders (n = 10) were collected using standard (STD) and boneplus (BONE) kernels. Hounsfield Units were converted to vBMD using specimen-specific calibrations. DECT STD and BONE images were generated using an established material decomposition method with 40 and 90 keV simulated monochromatic images. A proximal humerus bone section below the anatomic neck was used for vBMD analysis and FEM generation. FEMs were loaded to 1 % apparent strain for stiffness measurements.Between STD and BONE kernel images, average vBMD differed 0.9 mgK2HPO4/cc and 4.1 mg K2HPO4/cc, in SEQCT and DECT images, respectively. Significant differences occurred in DECT images (p = 0.001). BONE reconstructed images produced higher vBMD measures across both SEQCT and DECT images. The difference between STD and BONE in both SEQCT- and DECT-based FEMs persisted, with larger estimated stiffness in BONE models. For six of the models DECT-based had higher stiffness than SEQCT-based models using the same kernel, although these models differed between STD and BONE kernels. Differences in stiffness between STD and BONE derived models were similar across image types (DECT: 17.5 kN/mm; SEQCT: 19.0 kN/mm). Stiffness values were significantly different within SECT kernels and between SEQCT BONE and DECT STD models. This study shows important differences in vBMD and FEM stiffness that occur due to CT-based imaging parameters alone. These results indicate that consistent imaging parameters should be used for vBMD analysis and FEM input to avoid systematic measurement errors.

Patient-specific design and fabrication of a trapeziometacarpal joint orthosis using a computed tomography image-based finite element model

Patient-specific design and fabrication of a trapeziometacarpal joint orthosis using a computed tomography image-based finite element model

Development of a continuum-based, meshless, finite element modeling approach for representation of trabecular bone indentation

Implant subsidence into the underlying trabecular bone is a common problem in orthopaedic surgeries; however, the ability to pre-operatively predict implant subsidence remains limited. Current state-of-the-art computational models for predicting subsidence have issues addressing this clinical problem, often resulting from the size and complexity of existing subject-specific, image-based finite element (FE) models. The current study aimed to develop a simplified approach to FE modeling of subject-specific trabecular bone indentation resulting from implant penetration.Confined indentation experiments of human trabecular bone with flat- and sharp-tip indenters were simulated using FE analysis. A generalized continuum-level approach using a meshless smoothed particle hydrodynamics (SPH) approach and an isotropic crushable foam (CF) material model was developed for the trabecular bone specimens. Five FE models were generated with CF material parameters calibrated to cadaveric specimens spanning a range of bone mineral densities (BMD). Additionally, an alternative model configuration was developed that included consideration of bone marrow, with bone and marrow material parameters assigned to elements randomly according to bone volume (BV%) measurements of experimental specimens, owing to the non-uniform nature of trabecular bone tissue microstructure.Statistical analysis found significant correlation between the shapes of the numerical and experimental force-displacement curves. FE models accurately captured the bone densification patterns observed experimentally. Inclusion of marrow elements offered improved response prediction of the flat-tip indenter tests.Ultimately, the developed approach demonstrates the ability of a generalizable continuum-level SPH approach to capture bone variability using clinical bone imaging metrics without needing detailed image-based geometries, a significant step towards simplified subject-specific modeling of implant subsidence.

Heterogeneous Mechanical Stress and Interstitial Fluid Flow Predictions Derived from DCE-MRI for Rat U251N Orthotopic Gliomas.

Mechanical stress and fluid flow influence glioma cell phenotype in vitro, but measuring these quantities in vivo continues to be challenging. The purpose of this study was to predict these quantities in vivo, thus providing insight into glioma physiology and potential mechanical biomarkers that may improve glioma detection, diagnosis, and treatment. Image-based finite element models of human U251N orthotopic glioma in athymic rats were developed to predict structural stress and interstitial flow in and around each animal's tumor. In addition to accounting for structural stress caused by tumor growth, our approach has the advantage of capturing fluid pressure-induced structural stress, which was informed by in vivo interstitial fluid pressure (IFP) measurements. Because gliomas and the brain are soft, elevated IFP contributed substantially to tumor structural stress, even inverting this stress from compressive to tensile in the most compliant cases. The combination of tumor growth and elevated IFP resulted in a concentration of structural stress near the tumor boundary where it has the greatest potential to influence cell proliferation and invasion. MRI-derived anatomical geometries and tissue property distributions resulted in heterogeneous interstitial fluid flow with local maxima near cerebrospinal fluid spaces, which may promote tumor invasion and hinder drug delivery. In addition, predicted structural stress and interstitial flow varied markedly between irradiated and radiation-naïve animals. Our modeling suggests that relative to tumors in stiffer tissues, gliomas experience unusual mechanical conditions with potentially important biological (e.g., proliferation and invasion) and clinical consequences (e.g., drug delivery and treatment monitoring).

Comparison of Left Ventricular Function Derived from Subject-Specific Inverse Finite Element Modeling Based on 3D ECHO and Magnetic Resonance Images.

Three-dimensional echocardiography (3D ECHO) and magnetic resonance (MR) imaging are frequently used in patients and animals to evaluate heart functions. Inverse finite element (FE) modeling is increasingly applied to MR images to quantify left ventricular (LV) function and estimate myocardial contractility and other cardiac biomarkers. It remains unclear, however, as to whether myocardial contractility derived from the inverse FE model based on 3D ECHO images is comparable to that derived from MR images. To address this issue, we developed a subject-specific inverse FE model based on 3D ECHO and MR images acquired from seven healthy swine models to investigate if there are differences in myocardial contractility and LV geometrical features derived using these two imaging modalities. We showed that end-systolic and end-diastolic volumes derived from 3D ECHO images are comparable to those derived from MR images (R2=0.805 and 0.969, respectively). As a result, ejection fraction from 3D ECHO and MR images are linearly correlated (R2=0.977) with the limit of agreement (LOA) ranging from -17.95% to 45.89%. Using an inverse FE modeling to fit pressure and volume waveforms in subject-specific LV geometry reconstructed from 3D ECHO and MR images, we found that myocardial contractility derived from these two imaging modalities are linearly correlated with an R2 value of 0.989, a gradient of 0.895, and LOA ranging from -6.11% to 36.66%. This finding supports using 3D ECHO images in image-based inverse FE modeling to estimate myocardial contractility.

Biomechanical impact of different isthmus positions in mandibular first molar root canals: a finite element analysis.

This study used image-based finite element analysis (FEA) to assess the biomechanical changes in mandibular first molars resulting from alterations in the position of the root canal isthmus. A healthy mandibular first molar, characterized by two intact root canals and a cavity-free surface, was selected as the subject. A three-dimensional model for the molar was established using scanned images of the patient's mandibular teeth. Subsequently, four distinct finite element models were created, each representing varied root canal morphologies: non-isthmus (Group A), isthmus located at the upper 1/3 of the root (Group B), middle 1/3 of the root (Group C), and lower 1/3 of the root (Group D). A static load of 200 N was applied along the tooth's longitudinal axis on the occlusal surface to simulate regular chewing forces. The biomechanical assessment was conducted regarding the mechanical stress profile within the root dentin. The equivalent stress (Von Mises stress) was used to assess the biomechanical features of mandibular teeth under mechanical loading. In Group A (without an isthmus), the maximum stress was 22.2MPa, while experimental groups with an isthmus exhibited higher stresses, reaching up to 29.4MPa. All maximum stresses were concentrated near the apical foramen. The presence of the isthmus modified the stress distribution in the dentin wall of the tooth canal. Notably, dentin stresses at specific locations demonstrated differences: at 8mm from the root tip, Group B: 13.6MPa vs. Group A: 11.4MPa; at 3mm from the root tip, Group C: 14.2MPa vs. Group A: 4.5MPa; at 1mm from the root tip, Group D: 25.1MPa vs. Group A: 10.3MPa. The maximum stress in the root canal dentin within the isthmus region was located either at the top or bottom of the isthmus. A root canal isthmus modifies the stress profile within the dentin. The maximum stress occurs near the apical foramen and significantly increases when the isthmus is located closer to the apical foramina.

Image-based finite element modeling of air flow and thermal transport in Al-fiber-sintered porous materials

In this study, the pressure drop and heat transfer in a heat-transfer tube filled with a sintered porous medium comprising Al fibers were investigated using computational fluid dynamics (CFD) simulations. We reconstructed the sintered fibrous porous structure and generated a computational mesh using X-ray computed tomography data, and the simulated pressure drop and heat transfer agreed well with those reported by Enoki et al. [2021]. The differences in both the form coefficient and the permeability between two different samples can be explained from the differences in both the specific solid surface area and the porosity. Further, the CFD simulations indicated that thermal conduction of the Al solid phase enhanced the heat exchange between the air and Al fibers. To examine the effect of solid-phase thermal conduction, we generated a model of a heat-transfer tube having an ideal wire mesh structure, in which we introduced the interfacial thermal conductivity between Al fibers and the inner tube wall as additional factors. We also performed CFD simulations with four different shell-region thicknesses and found that even a very narrow gap of 5 to 10 μm heavily affected the heat-exchange performance because of the low thermal conductivity of the shell region.

Myocardial biomechanical effects of fetal aortic valvuloplasty

Fetal critical aortic stenosis with evolving hypoplastic left heart syndrome (CAS-eHLHS) can progress to a univentricular (UV) birth malformation. Catheter-based fetal aortic valvuloplasty (FAV) can resolve stenosis and reduce the likelihood of malformation progression. However, we have limited understanding of the biomechanical impact of FAV and subsequent LV responses. Therefore, we performed image-based finite element (FE) modeling of 4 CAS-eHLHS fetal hearts, by performing iterative simulations to match image-based characteristics and then back-computing physiological parameters. We used pre-FAV simulations to conduct virtual FAV (vFAV) and compared pre-FAV and post-FAV simulations. vFAV simulations generally enabled partial restoration of several physiological features toward healthy levels, including increased stroke volume and myocardial strains, reduced aortic valve (AV) and mitral valve regurgitation (MVr) velocities, reduced LV and LA pressures, and reduced peak myofiber stress. FAV often leads to aortic valve regurgitation (AVr). Our simulations showed that AVr could compromise LV and LA depressurization but it could also significantly increase stroke volume and myocardial deformational stimuli. Post-FAV scans and simulations showed FAV enabled only partial reduction of the AV dissipative coefficient. Furthermore, LV contractility and peripheral vascular resistance could change in response to FAV, preventing decreases in AV velocity and LV pressure, compared with what would be anticipated from stenosis relief. This suggested that case-specific post-FAV modeling is required to fully capture cardiac functionality. Overall, image-based FE modeling could provide mechanistic details of the effects of FAV, but computational prediction of acute outcomes was difficult due to a patient-dependent physiological response to FAV.

In situ X-ray imaging and numerical modeling of damage accumulation in C/SiC composites at temperatures up to 1200 °C

Carbon fiber reinforced silicon carbide matrix composites (C/SiC) have emerged as key materials for thermal protection systems owing to their high strength-to-weight ratio, high-temperature durability, resistance to oxidation, and outstanding reliability. However, manufacturing defects deteriorate the mechanical response of these composites under extreme thermal-force coupling conditions, prompting significant research attention. This study demonstrates a customized in situ loading device compatible with synchrotron radiation facilities, enabling high spatial and temporal resolution recording of internal material damage evolution and failure behavior under thermal-force coupling conditions. Infrared thermal radiation units in a confocal configuration were used to create ultra-high-temperature environments, offering advantages of compactness, rapid heating, and versatility. In situ tensile tests were conducted on C/SiC samples in a nitrogen atmosphere at both room temperature and 1200 °C. The high-resolution image data demonstrate various failure phenomena, such as matrix cracking and pore linkage. Image-based finite element simulations indicate that the temperature-dependent variation of the failure mechanism is attributable to thermal residual stresses and defect-induced stress concentrations. This work seamlessly integrates extreme mechanical testing methods with in situ observation techniques, providing a comprehensive solution for accurately quantifying crack initiation, pore connection, and failure behavior of C/SiC composites.

A Finite Element Model to Investigate the Stability of Osteochondral Grafts Within a Human Tibiofemoral Joint

Osteochondral grafting has demonstrated positive outcomes for treating articular cartilage defects by replacing the damaged region with a cylindrical graft consisting of bone with a layer of cartilage. However, factors that cause graft subsidence are not well understood. The aim of this study was to develop finite element (FE) models of osteochondral grafts within a tibiofemoral joint, suitable for an investigation of parameters affecting graft stability. Cadaveric femurs were used to experimentally calibrate the bone properties and graft-bone frictional forces for use in corresponding image-based FE models, generated from µCT scan data. Effects of cartilage defects and osteochondral graft repair were measured by examining contact pressure changes using further in vitro tests. Here, six defects were created in the femoral condyles, which were subsequently treated with osteochondral autografts or metal pins. Matching image-based FE models were created, and the contact patches were compared. The bone material properties and graft-bone frictional forces were successfully calibrated from the initial tests with good resulting levels of agreement (CCC = 0.87). The tibiofemoral joint experiment provided a range of cases that were accurately described in the resultant pressure maps and were well represented in the FE models. Cartilage defects and repair quality were experimentally measurable with good agreement in the FE model pressure maps. Model confidence was built through extensive validation and sensitivity testing. It was found that specimen-specific properties were required to accurately represent graft behaviour. The final models produced are suitable for a range of parametric testing to investigate immediate graft stability.

Pre-intervention myocardial stress is a good predictor of aortic valvoluplasty outcome for fetal critical aortic stenosis and evolving HLHS.

Fetal critical aortic stenosis with evolving hypoplastic left heart syndrome (CAS-eHLHS) causes biomechanical and functional aberrations, leading to a high risk of progression to hypoplastic left heart syndrome (HLHS) at birth. Fetal aortic valvuloplasty (FAV) can resolve outflow obstruction and may reduce progression risk. However, it is currently difficult to accurately predict which patients will respond to the intervention and become functionally biventricular (BV) at birth, as opposed to becoming functionally univentricular (UV). This prediction is important for patient selection, parental counselling, and surgical planning. Therefore, we investigated whether biomechanics parameters from pre-FAV image-based computations could robustly distinguish between CAS-eHLHS cases with BV or UV outcomes in a retrospective cohort. To do so we performed image-based finite element biomechanics modelling of nine CAS-eHLHS cases undergoing intervention and six healthy fetal control hearts, and found that a biomechanical parameter, peak systolic myofibre stress, showed a uniquely large difference between BV and UV cases, which had a larger magnitude effect than echocardiography parameters. A simplified equation was derived for quick and easy estimation of myofibre stress from echo measurements via principal component analysis. When tested on a retrospective cohort of 37 CAS-eHLHS cases, the parameter outperformed other parameters in predicting UV versus BV outcomes, and thus has a high potential of improving outcome predictions, if incorporated into patient selection procedures. Physiologically, high myocardial stresses likely indicate a healthier myocardium that can withstand high stresses and resist pathological remodelling, which can explain why it is a good predictor of BV outcomes. KEY POINTS: Predicting the morphological birth outcomes (univentricular versus biventricular) of fetal aortic valvuloplasty for fetal aortic stenosis with evolving HLHS is important for accurate patient selection, parental counselling and management decisions. Computational simulations show that a biomechanics parameter, pre-intervention peak systolic myofibre stress, is uniquely robust in distinguishing between such outcomes, outperforming all echo parameters. An empirical equation was developed to quickly compute peak systolic myofibre stress from routine echo measurements and was the best predictor of outcomes among a wide range of parameters tested.

Geometric determinants of the mechanical behavior of image-based finite element models of the intervertebral disc.

The intervertebral disc is an important structure for load transfer through the spine. Its injury and degeneration have been linked to pain and spinal fractures. Disc injury and spine fractures are associated with high stresses; however, these stresses cannot be measured, necessitating the use of finite element (FE) models. These models should include the disc's complex structure, as changes in disc geometry have been linked to altered mechanical behavior. However, image-based models using disc-specific structures have yet to be established. This study describes a multiphasic FE modeling approach for noninvasive estimates of subject-specific intervertebral disc mechanical behavior based on medical imaging. The models (n = 22) were used to study the influence of disc geometry on the predicted global mechanical response (moments and forces), internal local disc stresses, and tractions at the interface between the disc and the bone. Disc geometry was found to have a strong influence on the predicted moments and forces on the disc (R2 = 0.69-0.93), while assumptions regarding the side curvature (bulge) of the disc had only a minor effect. Strong variability in the predicted internal disc stresses and tractions was observed between the models (mean absolute differences of 5.1%-27.7%). Disc height had a systematic influence on the internal disc stresses and tractions at the disc-to-bone interface. The influence of disc geometry on mechanics highlights the importance of disc-specific modeling to estimate disc injury risk, loading on the adjacent vertebral bodies, and the mechanical environment present in disc tissues.

Comparative Biomechanical Evaluation of Unilateral and Bilateral Cages in Posterior Lumbar Interbody Fusion: Endplates Subsidence, Pedicle Screw Loosening and Implant Stability

Implant sinking, pedicle screw loosening, and instability are commonly observed complications in posterior lumbar interbody fusion (PLIF) surgeries, attributed to a range of mechanical, biomechanical, and environmental factors. To address these challenges, the utilization of unilateral cages positioned obliquely, along with bilateral posterior instrumentation, and employing a material with a comparable modulus to cortical bone, has shown promising efficacy. The present study employed image-based finite element analysis (FEA) to evaluate the influences on cage sinking, screw loosening, and construct stability. The outcomes revealed that obliquely positioned unilateral fusion cage constructs with posterior instrumentation yielded superior biomechanical outcomes, characterized by reduced range of motion and stress distortion at the cage-endplate and screw-bone junctions. Furthermore, these findings indicated a lower risk of fractures and diminished deformations in normal and traumatic events. Hence, the utilization of biocompatible cage materials and structural symmetry holds critical significance for achieving biomechanical success in PLIF surgery.

Multiscale morphological analysis of bone microarchitecture around Mg-10Gd implants.

The utilization of biodegradable magnesium (Mg)-based implants for restoration of bone function following trauma represents a transformative approach in orthopaedic application. One such alloy, magnesium-10weight percent gadolinium (Mg-10Gd), has been specifically developed to address the rapid degradation of Mg while enhancing its mechanical properties to promote bone healing. Previous studies have demonstrated that Mg-10Gd exhibits favorable osseointegration; however, it exhibits distinct ultrastructural adaptation in comparison to conventional implants like titanium (Ti). A crucial aspect that remains unexplored is the impact of Mg-10Gd degradation on the bone microarchitecture. To address this, we employed hierarchical three-dimensional imaging using synchrotron radiation in conjunction with image-based finite element modelling. By using the methods outlined, the vascular porosity, lacunar porosity and the lacunar-canaliculi network (LCN) morphology of bone around Mg-10Gd in comparison to Ti in a rat model from 4 weeks to 20 weeks post-implantation was investigated. Our investigation revealed that within our observation period, the degradation of Mg-10Gd implants was associated with significantly lower (p<0.05) lacunar density in the surrounding bone, compared to Ti. Remarkably, the LCN morphology and the fluid flow analysis did not significantly differ for both implant types. In summary, a more pronounced lower lacunae distribution rather than their morphological changes was detected in the surrounding bone upon the degradation of Mg-10Gd implants. This implies potential disparities in bone remodelling rates when compared to Ti implants. Our findings shed light on the intricate relationship between Mg-10Gd degradation and bone microarchitecture, contributing to a deeper understanding of the implications for successful osseointegration.

Microstructure Image-Based Finite Element Methodology to Design Abradable Coatings for Aero Engines

Upgrading abradable or wearable coatings in the high-temperature zone of aero engines is advised to increase the efficiency and high-density power in gas turbine engines for military or commercial fixed-wing and rotary-wing aircraft. The development of these coated materials is also motivated by minimization of the number of failures in the blade, as well as increasing their resistance to wear and erosion. It is suggested that abradable coatings or seals be used to accomplish this goal. The space between the rotor and the shroud is minimized thanks to an abradable seal at the blade’s tip. Coatings that can withstand abrasion are often multiphase materials sprayed through thermal spray methods, and which consist of a metal matAzmeerix, oxide particles, and void space. The maintenance of an ideal blend of qualities, such as erosion resistance and hardness, during production determines a seal’s effectiveness. The objective of this research is to develop microstructure-based modelling methodology which will mimic the coating wear process and subsequently help in designing the abradable coating composition. Microstructure modelling, meshing, and wear analysis using many tools such as Fusion360, Hyper Mesh, and LS-Dyna, have been employed to develop an abradable coating model and perform wear analysis using a simulated rub rig test. The relation between percentage composition and morphology variations of metal, oxide, and voids to the output parameters such as hardness, abradability, and other mechanical properties is explored using simulated finite analysis models of real micrographic images of abradable coatings.

Optimization of mechanical properties of bio-inspired Voronoi structures by genetic algorithm

Numerous patterns in nature comprise aesthetics structures and intriguing geometries. The geometric factors behind these patterns include area, aspect ratio, circularity, orientation, and regularity. In this study, the quasi-3D porous structures with varying regularities were systematically described and generated by the Voronoi segmentation. Then, genetic algorithms were applied to optimize the mechanical properties of these structures under tension condition. The image-based finite element simulation could calculate each structure's fitness value (strength). Afterward, structures with higher fitness values were selected and underwent the replication process, including crossover and mutation. New structures were generated based on the input structures. The above steps driven by the genetic algorithms were repeated 30 iterations until the tensile strength of the structures was improved and converged. Both simulation and tensile testing results as well as digital image correlation of specimens fabricated by stereolithography 3D printing showed that the Voronoi structures optimized by the genetic algorithm could enhance their stiffness, strength, and toughness values by ∼30%. This research has the potential to be applied in the fields of structural materials and biomimetic micro-aerial vehicles, which require lightweight, strength, and toughness simultaneously.

Tensile behaviors of nanoporous phenolic composites reinforced by 3D needle-punched preforms with different weave patterns

Needle-punched preform reinforced nanoporous phenolic composite (NPC) is a cost-effective lightweight ablative material for thermal protection systems. It consists of a needle-punched preform composed of alternating stacked fiber felt and stacked woven fabric, which plays a crucial role in the thermal insulation efficiency and mechanical properties. Herein, needle-punched preforms with different types of fabrics, including plain, twill, and satin, are deliberately designed for the reinforcement of NPC, and their tensile properties and fracture mechanism are investigated through experimentation and the CT image-based finite element method. At the same preform density, the composite with plain fabric exhibits the highest tensile strength of 146 ± 6.3 MPa, a tensile modulus of 8.2 ± 0.2 GPa, and an elongation of 2.1 ± 0.1%, due to having the highest number of crossover points and the largest undulation of fiber yarns. Furthermore, the CT image-based simulation reveals the failure mechanism of NPC. The von mises stress of woven fabric is higher than that of fiber felt, which demonstrates the woven fabric is the main load-bearing structure and determines the mechanical of composite. These results will provide useful guidance for the structural design and optimization of needle-punched preform reinforced composites.

Finite element analysis (FEA) of femur to predict biomechanical properties and its validation

The aim of this preliminary study is to predict the biomechanical properties as strain of the femur using computed tomography medical image-based finite element analysis. Further, the finite element results were compared with the experimental results for validation purpose. The main objective of the study was to establish and validate a finite element model used to analyse biomechanical properties. For this, four sheep femur specimens were used to predict normal strain in x-x direction. Computed tomography scans were taken from each femur sample. Subsequently, an accurate image-based 3D model of each femur specimen was segmented and reconstructed with its material properties. Furthermore, the experimental test was performed at 8° bone inclination in the coronal plane and a compressive load was applied at displacement rate of 0.2 mm/s to the femur head in a vertical downward direction until fracture. Consequently, the numerical results of predicted normal strain were compared with the experimental measured normal strain. The average deviation between finite element and experimental results was found to be −11.11%–20.10% in strains. Hence, it is concluded that the results obtained using finite element analysis are in good agreement with the experimental results.

Tensile failure analysis of laser stake-welded T-joints: X-ray tomography based conformal finite element and fractal theoretical approaches

To promote the application of complex structural materials in large size, the interfacial joint is the key component for sandwich structure to wipe off the sample size limitation and fabricate multi-material structures. However, the laser welding process inevitably induced void defects in the welded joint, influencing the mechanical properties of the whole structure. In this paper, the aluminum alloy laser stake-welded T-joints were fabricated in a laser welding system. The geometric deviations and internal void defects of T-joints were visualized and quantified by micro X-ray computed tomography (μCT). Two image-based finite element models were reconstructed from CT images to discuss the effect of different manufacturing-induced defects on the tensile behaviour of T-joints. The geometric deviations of T-joints make a response for the decrease of strength and the void defects result in the decrease of Young’s modulus of T-joints. The fractal dimension theory, based on the value of void volume and the distance of the void to the interface obtained from CT images, was developed to predict the strength of the T-joint. The results matched well with the experimental results.