Degenerative disc disease (DDD) is a common condition in elderly population that can be painful and can significantly affect individual's quality of life. Diagnosis of DDD allows prompt corrective actions but it is challenging due to the absence of any symptoms at early stages. In studying disc degeneration, measurement of the range of motion (RoM) and loads acting on the spine are crucial factors. However, direct measurement of RoM involves increased instrumentation and risk. In this paper, an innovative method is proposed for calculating RoM, emphasizing repeatability and reliability by considering the posterior thickness of the spine. This is achieved by offsetting the position of markers in relation to the actual vertebral loci. Three geometrically identical finite element models of L3-L4 are developed from a CT scan with different types of elements, and thereafter, mesh element-related metrics are provided for the assessment of the quality of models. The model with the best mesh quality is used for further analysis, where RoM are within ranges as reported in literature and in vivo experiment results. Various kinds of stresses acting on individual components including facet joints are analysed for normal and abnormal loading conditions. The results showed that the stresses in abnormal load conditions for all components including cortical (76.67 MPa), cancellous (69.18 MPa), annulus (6.30 MPa) and nucleus (0.343 MPa) are significantly greater as compared to normal loads (49.96 MPa, 44.2 MPa, 4.28 MPa and 0.23 MPa respectively). However, stress levels for both conditions are within safe limits (167-215 MPa for cortical, 46 MPa for the annulus and 3 MPa for facets). The results obtained could be used as a baseline motion and stresses of healthy subjects based on their respective lifestyles, which could benefit clinicians to suggest corrective actions for those affected by DDD.
Complication rates of anterior cruciate ligament reconstruction (ACL-R) were reported to be around 15% although it is a common arthroscopic procedure with good outcomes. Breakage and migration of fixators are still possible even months after surgery. A fixator with optimum stability can minimise those two complications. Factors that affect the stability of a fixator are its configuration, material, and design. Thus, this paper aims to analyse the biomechanical effects of different types of fixators (cross-pin, interference screw, and cortical button) towards the stability of the knee joint after ACL-R. In this study, finite element modelling and analyses of a knee joint attached with double semitendinosus graft and fixators were carried out. Mimics and 3-Matic softwares were used in the development of the knee joint models. Meanwhile, the graft and fixators were designed by using SolidWorks software. Once the meshes of all models were finished in 3-Matic, simulation of the configurations was done using MSC Marc Mentat software. A 100-N anterior tibial load was applied onto the tibia to simulate the anterior drawer test. Based on the findings, cross-pin was found to have optimum stability in terms of stress and strain at the femoral fixation site for better treatment of ACL-R.
Pilon fractures can be caused by high-energy vertical forces which may result in long-term patient immobilization. Many experts in orthopedic surgery recommend the use of a Delta external fixator for type III Pilon fracture treatment. This device can promote immediate healing of fractured bone, minimizing the rate of complications as well as allowing early mobilization. The characteristics of different types of the Delta frame have not been demonstrated yet. By using the finite element method, this study was conducted to determine the biomechanical characteristics of six different configurations (Model 1 until Model 6). CT images from the lower limb of a healthy human were used to reconstruct three-dimensional models of foot and ankle bones. All bones were assigned with isotropic material properties and the cartilages were assigned to exhibit hyperelasticity. A linear link was used to simulate 37 ligaments at the ankle joint. Axial loads of 70 and 350 N were applied at the proximal tibia to simulate the stance and swing phase. The metatarsals and calcaneus were fixed distally in order to prevent rigid body motion. A synthetic ankle bone was used to validate the finite element model. The simulated results showed that Delta3 produced the highest relative micromovement (0.09 mm, 7 μm) during the stance and swing phase, respectively. The highest equivalent von Mises stress was found at the calcaneus pin of the Delta4 (423.2 MPa) as compared to others. In conclusion, Delta1 external fixator was the most favorable option for type III Pilon fracture treatment. Graphical abstract ᅟ.
The provision of the most suitable rehabilitation treatment for stroke patient remains an ongoing challenge for clinicians. Fully understanding the pathomechanics of the upper limb will allow doctors to assist patients with physiotherapy treatment that will aid in full arm recovery. A biomechanical study was therefore conducted using the finite element (FE) method. A three-dimensional (3D) model of the human wrist was reconstructed using computed tomography (CT)-scanned images. A stroke model was constructed based on pathological problems, i.e. bone density reductions, cartilage wane, and spasticity. The cartilages were reconstructed as per the articulation shapes in the joint, while the ligaments were modelled using linear links. The hand grip condition was mimicked, and the resulting biomechanical characteristics of the stroke and healthy models were compared. Due to the lower thickness of the cartilages, the stroke model reported a higher contact pressure (305 MPa), specifically at the MC1-trapezium. Contrarily, a healthy model reported a contact pressure of 228 MPa. In the context of wrist extension and displacement, the stroke model (0.68° and 5.54 mm, respectively) reported a lower magnitude than the healthy model (0.98° and 9.43 mm, respectively), which agrees with previously reported works. It was therefore concluded that clinicians should take extra care in rehabilitation treatment of wrist movement in order to prevent the occurrence of other complications. Graphical abstract ᅟ.
Carbon nanotubes (CNTs) are potentially ideal tips for atomic force microscopy (AFM) due to the robust mechanical properties, nanoscale diameter and also their ability to be functionalized by chemical and biological components at the tip ends. This contribution develops the idea of using CNTs as an AFM tip in computational analysis of the biological cells. The proposed software was ABAQUS 6.13 CAE/CEL provided by Dassault Systems, which is a powerful finite element (FE) tool to perform the numerical analysis and visualize the interactions between proposed tip and membrane of the cell. Finite element analysis employed for each section and displacement of the nodes located in the contact area was monitored by using an output database (ODB). Mooney-Rivlin hyperelastic model of the cell allows the simulation to obtain a new method for estimating the stiffness and spring constant of the cell. Stress and strain curve indicates the yield stress point which defines as a vertical stress and plan stress. Spring constant of the cell and the local stiffness was measured as well as the applied force of CNT-AFM tip on the contact area of the cell. This reliable integration of CNT-AFM tip process provides a new class of high performance nanoprobes for single biological cell analysis.
Stainless steels are increasingly used in construction today, especially in harsh environments, in which steel corrosion commonly occurs. Cold-formed stainless steel structures are currently increasing in popularity because of its efficiency in load-bearing capacity and its appealing architectural appearance. Cold-rolling and press-braking are the cold-working processes used in the forming of stainless steel sections. Press braking can produce large cross-sections from thin to thick-walled sections compared to cold-rolling. Cold-forming in press-braked sections significantly affect member behaviour and joints; therefore, they have attained great attention from many researchers to initiate investigations on those effects. This paper examines the behaviour of residual stress distribution of stainless steel press-braked sections by implementing three-dimensional finite element (3D-FE) technique. The study proposed a full finite element procedure to predict the residual stresses starting from coiling-uncoiling to press-braking. This work considered material anisotropy to examine its effect on the residual stress distribution. The technique adopted was compared with different finite element techniques in the literature. This study also provided a parametric study for three corner radius-to-thickness ratios looking at the through-thickness residual stress distribution of four stainless steels (i.e., ferritic, austenitic, duplex, lean duplex) in which have their own chemical composition. In conclusion, the comparison showed that the adopted technique provides a detailed prediction of residual stress distribution. The influence of geometrical aspects is more pronounced than the material properties. Neglecting the material anisotropy shows higher shifting in the neutral axis. The parametric study showed that all stainless steel types have the same stress through-thickness distribution. Moreover, R/t ratios' effect is insignificant in all transverse residual stress distributions, but a slight change to R/t ratios can affect the longitudinal residual stress distribution.
The curvature correction factor is an important parameter in the stress calculation formulation of a helical extension spring, which describes the effect of spring wire curvature on the stress increase towards its inner radius. In this study, the parameters affecting the curvature correction factor were investigated through theoretical and numerical methods. Several finite element (FE) models of an extension spring were generated to obtain the distribution of the tensile stress in the spring. In this investigation, the hook orientation and the number of coils of the extension spring showed significant effects on the curvature correction factor. These parameters were not considered in the theoretical model for the calculation of the curvature correction factor, causing a deviation between the results of the FE model and the theoretical approach. A set of equations is proposed for the curvature correction factor, which relates both the spring index and the number of coils. These equations can be applied directly to the design of extension springs with a higher safety factor.
Triple-layered composite plates are created by joining three composite layers using shear connectors. These layers, which are assumed to be always in contact and able to move relatively to each other during deformation, could be the same or different in geometric dimensions and material. They are applied in various engineering fields such as ship-building, aircraft wing manufacturing, etc. However, there are only a few publications regarding the calculation of this kind of plate. This paper proposes novel equations, which utilize Mindlin's theory and finite element modelling to simulate the forced vibration of triple-layered composite plates with layers connected by shear connectors subjected to a moving load. Moreover, a Matlab computation program is introduced to verify the reliability of the proposed equations, as well as the influence of some parameters, such as boundary conditions, the rigidity of the shear connector, thickness-to-length ratio, and the moving load velocity on the dynamic response of the composite plate.
The environmental condition in which the Royal Malaysian Airforce is currently operating its aircraft is prone to corrosion. This is due to the high relative humidity and temperature. With most of its aircraft being in the legacy aircraft era, the aircraft's main construction consists of the aluminium 2024 material. However, this material is prone to corrosion, thus reducing fatigue life and leading to fatigue failure. Using the concept of either Safe Life or Damage Tolerance as its fatigue design philosophy, the RMAF adopts the Aircraft Structure Integrity Program (ASIP) to monitor its structural integrity. With the current problem of not having the structural limitation on corrosion-damaged structure, the RMAF has embarked on its fatigue testing method. Finite Element (FE) studies and flight tests were conducted, and the outcome is summarized. The conclusion is that the longeron tested on the aircraft can withstand the operational load, and its yield strength is below the ultimate yield strength of the material. These research outcomes will also enhance the ASIP for other aircraft platforms in the RMAF fleet for its structure life assessment or service life extension program.
Wire ropes undergo a fretting fatigue condition when subjected to axial and bending loads. The fretting behavior of wires are classified as line contact and trellis point of contact. The experimental study on the fatigue of wire ropes indicates that most of the failure occurs due to high localized stresses at trellis point of contact. A continuum damage mechanics approach was previously proposed to estimate the fatigue life estimation of wire ropes. The approach majorly depends on the high value of localized stresses as well as the micro-slippage occurs at the contact region. Finite element approach has been used to study radial and axial distribution of stresses and displacement in order to clearly understand the evolution of stresses and existence of relative displacements between neighboring wires under various loading and frictional conditions. The relative movements of contacting wires are more when friction is not considered. In the presence of friction, the relative movement occurs at the boundaries of the contact region. The location of microslip in the presence of friction is backed by the experimental observation stating the crack is initiated at or the outer boundary of the contact spot. The existence of slip is due to different displacement of outer and central wires.
The purpose of this study was to synthesize high-quality recycled α-Fe2O3 to improve its complex permittivity properties by reducing the particles to nanosize through high energy ball milling. Complex permittivity and permeability characterizations of the particles were performed using open-ended coaxial and rectangular waveguide techniques and a vector network analyzer. The attenuation characteristics of the particles were analyzed with finite element method (FEM) simulations of the transmission coefficients and electric field distributions using microstrip model geometry. All measurements and simulations were conducted in the 8-12 GHz range. The average nanoparticle sizes obtained after 8, 10 and 12 h of milling were 21.5, 18, and 16.2 nm, respectively, from an initial particle size of 1.73 µm. The real and imaginary parts of permittivity increased with reduced particle size and reached maximum values of 12.111 and 0.467 at 8 GHz, from initial values of 7.617 and 0.175, respectively, when the particle sizes were reduced from 1.73 µm to 16.2 nm. Complex permeability increased with reduced particle size while the enhanced absorption properties exhibited by the nanoparticles in the simulations confirmed their ability to attenuate microwaves in the X-band frequency range.
The compressive behaviour of column members can be considerably affected by local buckling, material yielding and local end conditions. In this paper, the effects of the loading conditions at the ends of plain channel section columns subjected to uniformly compressed loading, and fixed conditions at the column ends with respect to global rotations, was examined. Finite element simulation was employed to look at the post-buckled response of thin-walled, plain channel section columns that covered the complete loading history of the compression columns from the onset of elastic local buckling through the nonlinear elastic and elastoplastic post-buckling phases of behaviour to final collapse and unloading. Two types of loading conditions were considered: the first was one that has been used practically in tests whereby one end is loaded with a moving top platen while the other end is fixed at the lower platen, but, for the second loading condition, both ends were loaded with equally moving top and lower platens. These two conditions were shown to lead to quite different characteristic interactive responses of the columns due to mode jumping in the buckling mode for the locally rotationally constrained case.
Most orthopaedic cases that involved with bone fracture are normally treated with medical implants. To be noticed that some precautions in terms of biomechanical and biomaterial properties are necessary for a successful post-sur- gery process. The biomechanical evaluation of implants could be carried out using computing and engineering technologies. However, in the computer simulation, some assumptions are needed as the limitations on computer resources and data input. This review focuses on the current method of developing the finite element model for patients with specific values of material properties for lower limb part such as hip, knee and ankle joint. Previous literature was reviewed from which keywords and search engines were identified. In this review, inclusion and exclusion criteria were used to limit the literature search. We reviewed the state-of-the-art in this area and provide recommendations for future research. In conclusion, the previous published reports illustrated different methods to develop numerical models.
Introduction: Detecting breast cancer at earlier stage is crucial to increase the survival rate. Mammography as the golden screening tool has shown to be less effective for younger women due to denser breast tissue. Infrared Ther- mography has been touted as an adjunct modality to mammography. Further investigation of thermal distribution in breast cancer patient is important prior to its clinical interpretation. Therefore, thermal profiling using 3D compu- tational simulation was carried out to understand the effect of changes in size and location of tumour embedded in breast to the surface temperature distribution at different breast densities. Methods: Extremely dense (ED) and pre- dominantly fatty dense (PF) breast models were developed and simulated using finite element analysis (FEA). Pennes’ bioheat equation was adapted to show the heat transfer mechanism by providing appropriate thermophysical prop- erties in each tissue layer. 20 case studies with various tumour size embedded at two asymmetrical positions in the breast models were analysed. Quantitative and qualitative analyses were performed by recording the temperature values along the arc of breast, calculating of temperature difference at the peaks and comparing multiple thermal images. Results: Bigger size of tumour demands a larger increase in breast surface temperatures. As tumour is located far from the centre of the breast or near to the edge, there was a greater shift of temperature peak. Conclusion: Size and location of tumour in various levels of breast density should be considered as a notable factor to thermal profile on breast when using thermography for early breast cancer detection.
Blood flow through a bifurcated artery with the presence of an overlapping stenosis located at parent’s arterial lumen under the action of a uniform external magnetic field is studied in this paper. Blood is treated as an electrically conducting fluid which exhibits the Magnetohydrodynamics principle and it is characterized by a Newtonian fluid model. The governing equations are discretized using a stabilization technique of finite element known as Galerkin least-squares. The maximum velocity and pressure drop evaluated in this present study are compared with the results found in previous literature and COMSOL Multiphysics. The solutions found in a satisfactory agreement, thus verify the source code is working properly. The effects of dimensionless parameters of Hartmann and Reynolds numbers in the fluid’s velocity and pressure are examined in details with further scientific discussions.
Numerical simulation of the behaviour of blood flow through a stenosed bifurcated
artery with the presence of single mild stenosis at parent artery is investigated. The
flow analysis applies the incompressible, steady, three-dimensional Navier-Stokes equations
for non-Newtonian generalized power law fluids. Behaviour of blood flow is simulated
numerically using COMSOL Multiphysicsthat based on finite element method.The
results showthe effect of severity of stenosis on flow characteristics such as axial velocity
and its exhibit flow recirculation zone for analysis on streamlines pattern.
Mechanical failure of zirconia-based full-arch implant-supported fixed dental prostheses (FAFDPs) remains a critical issue in prosthetic dentistry. The option of full-arch implant treatment and the biomechanical behaviour within a sophisticated screw-retained prosthetic structure have stimulated considerable interest in fundamental and clinical research. This study aimed to analyse the biomechanical responses of zirconia-based FAFDPs with different implant configurations (numbers and distributions), thereby predicting the possible failure sites and the optimum configuration from biomechanical aspect by using finite element method (FEM). Five 3D finite element (FE) models were constructed with patient-specific heterogeneous material properties of mandibular bone. The results were reported using volume-averaged von-Mises stresses (σVMVA) to eliminate numerical singularities. It was found that wider placement of multi-unit copings was preferred as it reduces the cantilever effect on denture. Within the limited areas of implant insertion, the adoption of angled multi-unit abutments allowed the insertion of oblique implants in the bone and wider distribution of the multi-unit copings in the prosthesis, leading to lower stress concentration on both mandibular bone and prosthetic components. Increasing the number of supporting implants in a FAFDPs reduced loading on each implant, although it may not necessarily reduce the stress concentration in the most posterior locations significantly. Overall, the 6-implant configuration was a preferable configuration as it provided the most balanced mechanical performance in this patient-specific case.
Denture fracture is a common clinical complication caused by improper material selection, design, or fabrication technique. This study aimed to investigate the effect of two attachment systems on fracture risk of the implant-overdentures (IOD) via finite element analysis (FEA), using the force distributions obtained from patients' occlusal analyses and to compare the obtained results with the clinical complications associated with these attachments. A three-dimensional jaw model comprised of the edentulous bones was constructed. Three types of mandibular prostheses including complete denture (CD) (model LCD), IOD with Locator attachment (model LID-L), and IOD with telescopic attachment (model LID-T), as well as a maxillary CD (model UCD) were assembled. The vertical occlusal forces at anterior and posterior quadrants were obtained from the patients wearing mandibular CDs or IODs. The FEA results were further compared with the mechanical failures of different prostheses observed at patient recalls. In overall, the fracture risk of mandibular prostheses was lower than the maxillary compartments. The UCD opposing LCD underwent higher strains than that opposing LID-L and LID-T, which was mostly concentrated at the anterior mid-palatal polished surface. On the other hand, LID-L showed the lowest strain, followed by LID-T, and LCD. The obtained results were consistent with the clinical complications observed in the patient recalls.
Electrospinning is a simple and efficient process for producing sub-micron fibres. However, the process has many variables, and their effects on the non-woven mesh of fibres is complex. In particular, the effects on the mechanical properties of the fibre meshes are poorly understood. This paper conducts a parametric study, where the concentration and bloom strength of the gelatin solutions are varied, while all electrospinning process parameters are held constant. The effects on the fibrous meshes are monitored using scanning electron microscopy and mechanical testing under uniaxial tension. Mesh mechanical properties are relatively consistent, despite changes to the solutions, demonstrating the robustness of electrospinning. The gel strength of the solution is shown to have a statistically significant effect on the morphology, stiffness and strength of the meshes, while the fibre diameter has surprisingly little influence on the stiffness of the meshes. This experimental finding is supported by finite element analysis, demonstrating that the stiffness of the meshes is controlled by the volume fraction, rather than fibre diameter. Our results demonstrate the importance of understanding how electrospinning parameters influence the pore size of the meshes, as controlling fibre diameter alone is insufficient for consistent mechanical properties.
NiTi arch wires are used widely in orthodontic treatment due to its superelastic and biocompatibility properties. In brackets configuration, the force released from the arch wire is influenced by the sliding resistances developed on the arch wire-bracket contact. This study investigated the evolution of the forces released by a rectangular NiTi arch wire towards possible intraoral temperature and deflection changes. A three dimensional finite element model was developed to measure the force-deflection behavior of superelastic arch wire. Finite element analysis was used to distinguish the martensite fraction and phase state of arch wire microstructure in relation to the magnitude of wire deflection. The predicted tensile and bending results from the numerical model showed a good agreement with the experimental results. As contact developed between the wire and bracket, binding influenced the force-deflection curve by changing the martensitic transformation plateau into a slope. The arch wire recovered from greater magnitude of deflection released lower force than one recovered from smaller deflection. In contrast, it was observed that the plateau slope increased from 0.66N/mm to 1.1N/mm when the temperature was increased from 26°C to 46°C.