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.
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.
Single-cell analysis has become the interest of a wide range of biological and biomedical engineering research. It could provide precise information on individual cells, leading to important knowledge regarding human diseases. To perform single-cell analysis, it is crucial to isolate the individual cells before further manipulation is carried out. Recently, microfluidic biochips have been widely used for cell trapping and single cell analysis, such as mechanical and electrical detection. This work focuses on developing a finite element simulation model of single-cell trapping system for any types of cells or particles based on the hydrodynamic flow resistance (Rh) manipulations in the main channel and trap channel to achieve successful trapping. Analysis is carried out using finite element ABAQUS-FEA™ software. A guideline to design and optimize single-cell trapping model is proposed and the example of a thorough optimization analysis is carried out using a yeast cell model. The results show the finite element model is able to trap a single cell inside the fluidic environment. Fluid's velocity profile and streamline plots for successful and unsuccessful single yeast cell trapping are presented according to the hydrodynamic concept. The single-cell trapping model can be a significant important guideline in designing a new chip for biomedical applications.
The present study was conducted to compare the stability of four commercially available implants by investigating the focal stress distributions and relative micromotion using finite element analysis. Variations in the numbers of pegs between the implant designs were tested. A load of 750 N was applied at three different glenoid positions (SA: superior-anterior; SP: superior-posterior; C: central) to mimic off-center and central loadings during activities of daily living. Focal stress distributions and relative micromotion were measured using Marc Mentat software. The results demonstrated that by increasing the number of pegs from two to five, the total focal stress volumes exceeding 5 MPa, reflecting the stress critical volume (SCV) as the threshold for occurrence of cement microfractures, decreased from 8.41 to 5.21 % in the SA position and from 9.59 to 6.69 % in the SP position. However, in the C position, this change in peg number increased the SCV from 1.37 to 5.86 %. Meanwhile, micromotion appeared to remain within 19-25 µm irrespective of the number of pegs used. In conclusion, four-peg glenoid implants provide the best configuration because they had lower SCV values compared with lesser-peg implants, preserved more bone stock, and reduced PMMA cement usage compared with five-peg implants.
The functions of dams are for water supply, flood control and hydroelectric power generation. A concrete faced rockfill dam (CFRD) is preferred by dam consultants due to its many advantages. It is designed to withstand all applied loads, namely gravity load due to its massive weight and hydrostatic load due to water thrust from the reservoir. The structural response of the Bakun CFRD which is the second highest dam in the world after the Shuibuya Dam, was analyzed using finite element method. A two-dimensional plane strain finite element analysis of the non-linear Duncan-Chang hyperbolic model was used in the study of the dam in respect to the deformation and stresses of the main dam of Bakun CFRD project. The Dead-Birth-Ghost element technique was used to simulate sequences of the construction of the dam. The comparison of rigid and flexible foundation on the behaviour of the dam is discussed in this study. In the finite element modeling, the concrete slab on the upstream was represented through a six-noded element, while the interface characteristic between the dam body and the concrete slab was modeled using an interface element. The maximum settlement and stresses of the cross-section were found and their distribution was discussed and tabulated in the form of contours. The effect of reservoir filling loading had a gradual effect on the dam response behaviour. Comparisons with no water impoundment in the dam were also discussed.
In electromagnetic dosimetry, anatomical human models are commonly obtained by segmentation of magnetic resonance imaging or computed tomography scans. In this paper, a human head model extracted from thermal infrared images is examined in terms of its applicability to specific absorption rate (SAR) calculations. Since thermal scans are two-dimensional (2D) representation of surface temperature, this allows researchers to overcome the extensive computational demand associated with 3D simulation. The numerical calculations are performed using the finite-difference time-domain method with mesh sizes of 2 mm at 900 MHz plane wave irradiation. The power density of the incident plane wave is assumed to be 10 W/m(2). Computations were compared with a realistic anatomical head model. The results show that although there were marked differences in the local SAR distribution in the various tissues in the two models, the 1 g peak SAR values are approximately similar in the two models.
In this study, a hybridized neuro-genetic optimization methodology realized by embedding finite element analysis (FEA) trained artificial neural networks (ANN) into genetic algorithms (GA), is used to optimize temperature control in a ceramic based continuous flow polymerase chain reaction (CPCR) device. The CPCR device requires three thermally isolated reaction zones of 94 degrees C, 65 degrees C, and 72 degrees C for the denaturing, annealing, and extension processes, respectively, to complete a cycle of polymerase chain reaction. The most important aspect of temperature control in the CPCR is to maintain temperature distribution at each reaction zone with a precision of +/-1 degree C or better, irrespective of changing ambient conditions. Results obtained from the FEA simulation shows good comparison with published experimental work for the temperature control in each reaction zone of the microfluidic channels. The simulation data are then used to train the ANN to predict the temperature distribution of the microfluidic channel for various heater input power and fluid flow rate. Once trained, the ANN analysis is able to predict the temperature distribution in the microchannel in less than 20 min, whereas the FEA simulation takes approximately 7 h to do so. The final optimization of temperature control in the CPCR device is achieved by embedding the trained ANN results as a fitness function into GA. Finally, the GA optimized results are used to build a new FEA model for numerical simulation analysis. The simulation results for the neuro-genetic optimized CPCR model and the initial CPCR model are then compared. The neuro-genetic optimized model shows a significant improvement from the initial model, establishing the optimization method's superiority.
The mechanical behavior of graphene/polymer interfaces in the graphene-reinforced epoxy nanocomposite is one of the factors that dictates the deformation and damage response of the nanocomposites. In this study, hybrid molecular dynamic (MD) and finite element (FE) simulations of a graphene/polymer nanocomposite are developed to characterize the elastic-damage behavior of graphene/polymer interfaces under a tensile separation condition. The MD results show that the graphene/epoxy interface behaves in the form of elastic-softening exponential regressive law. The FE results verify the adequacy of the cohesive zone model in accurate prediction of the interface damage behavior. The graphene/epoxy cohesive interface is characterized by normal stiffness, tensile strength, and fracture energy of 5 × 10-8 (aPa·nm-1), 9.75 × 10-10 (nm), 2.1 × 10-10 (N·nm-1) respectively, that is followed by an exponential regressive law with the exponent, α = 7.74. It is shown that the commonly assumed bilinear softening law of the cohesive interface could lead up to 55% error in the predicted separation of the interface.
This study presents an efficient finite element analysis technique which shows great versatility in
modelling of precast composite flooring system subjected to static loadings. The method incorporates sliding and opening in the analysis of composite structures using the interface element which was specifically designed to simulate the actual behaviour at the interfaces between contacting materials. A three-dimensional finite element model of the precast composite slab which exhibits discontinuous behaviour was performed to demonstrate the potential and applicability of the proposed method of analysis. The results of the analysis demonstrate that the overall response of a discontinuous system to external loading is significantly affected by the bonding condition at the interfaces between the contacting materials.
Conventional seismic rehabilitation methods may not be suitable for some buildings owing to their high cost and time-consuming foundation work. In recent years, viscoelastic dampers (VEDs) have been widely used in many mid- and high-rise buildings. This study introduces a viscoelastic passive control system called rotary rubber braced damper (RRBD). The RRBD is an economical, lightweight, and easy-to-assemble device. A finite element model considering nonlinearity, large deformation, and material damage is developed to conduct a parametric study on different damper sizes under pushover cyclic loading. The fundamental characteristics of this VED system are clarified by analyzing building structures under cyclic loading. The result show excellent energy absorption and stable hysteresis loops in all specimens. Additionally, by using a sinusoidal shaking table test, the effectiveness of the RRBD to manage the response displacement and acceleration of steel frames is considered. The RRBD functioned at early stages of lateral displacement, indicating that the system is effective for all levels of vibration. Moreover, the proposed damper shows significantly better performance in terms of the column compression force resulting from the brace action compared to chevron bracing (CB).
This paper presents the design and evaluation of a high force density fishbone shaped electrostatic comb drive actuator. This comb drive actuator has a branched structure similar to a fishbone, which is intended to increase the capacitance of the electrodes and hence increase the electrostatic actuation force. Two-dimensional finite element analysis was used to simulate the motion of the fishbone shaped electrostatic comb drive actuator and compared against the performance of a straight sided electrostatic comb drive actuator. Performances of both designs are evaluated by comparison of displacement and electrostatic force. For both cases, the active area and the minimum gap distance between the two electrodes were constant. An active area of 800 × 300 μm, which contained 16 fingers of fishbone shaped actuators and 40 fingers of straight sided actuators, respectively, was used. Through simulation, improvement of drive force of the fishbone shaped electrostatic comb driver is approximately 485% higher than conventional electrostatic comb driver. These results indicate that the fishbone actuator design provides good potential for applications as high force density electrostatic microactuator in MEMS systems.
This article uses finite volume and finite element methods for optimization of the artificial hair cell sensor. The performance of the sensor was investigated for different materials such as sicon and polysilicon and by varying hair cell dimensions including width and length. The silicon material which has low young modulus was proposed based on the simulation performance. The performance of the hair cell sensor was achieved by increasing the hair cell length while increasing the width did not significantly influence the performance. The
performance of the sensor was studied for its viscous force, deflection, von mises stress and sensitivity. From the simulation, the hair cell with a length of 1600 µm and 80 µm width was suggested for the subsequent analysis. Another way to improve the performance was by modifying the hair cell geometry and it was proved that the modified hair cell was more sensitive, based on the deflection. The angle of flow that hit the hair cell also affected the deflection of the sensor where the zero angle flow which was parallel to the substrate was the most effective angle. The limitations of the performance of hair cell for various fluid velocity were also discussed in this paper.
The mechanical properties of tissue scaffolds are essential in providing stability for tissue repair and growth. Thus, the ability of scaffolds to withstand specific loads is crucial for scaffold design. Most research on scaffold pores focuses on grids with pore size and gradient structure, and many research models are based on scaffolding with vertically arranged holes. However, little attention is paid to the influence of the distribution of holes on the mechanical properties of the scaffold. To address this gap, this research investigates the effect of pore distribution on the mechanical properties of tissue scaffolds. The study involves four types of scaffold designs with regular and staggered pore arrangements and porosity ranging from 30% to 80%. Finite element analysis (FEA) was used to compare the mechanical properties of different scaffold designs, with von-Mises stress distribution maps generated for each scaffold. The results show that scaffolds with regular vertical holes exhibit a more uniform stress distribution and better mechanical performance than those with irregular holes. In contrast, the scaffold with a staggered arrangement of holes had a higher probability of stress concentration. The study emphasized the importance of balancing porosity and strength in scaffold design.
Experimental findings indicate that in-situ chondrocytes die readily following impact loading, but remain essentially unaffected at low (non-impact) strain rates. This study was aimed at identifying possible causes for cell death in impact loading by quantifying chondrocyte mechanics when cartilage was subjected to a 5% nominal tissue strain at different strain rates. Multi-scale modelling techniques were used to simulate cartilage tissue and the corresponding chondrocytes residing in the tissue. Chondrocytes were modelled by accounting for the cell membrane, pericellular matrix and pericellular capsule. The results suggest that cell deformations, cell fluid pressures and fluid flow velocity through cells are highest at the highest (impact) strain rate, but they do not reach damaging levels. Tangential strain rates of the cell membrane were highest at the highest strain rate and were observed primarily in superficial tissue cells. Since cell death following impact loading occurs primarily in superficial zone cells, we speculate that cell death in impact loading is caused by the high tangential strain rates in the membrane of superficial zone cells causing membrane rupture and loss of cell content and integrity.
The relationship between microarchitecture to the failure mechanism and mechanical properties can be assessed through experimental and computational methods. In this study, both methods were utilised using bovine cadavers. Twenty four samples of cancellous bone were extracted from fresh bovine and the samples were cleaned from excessive marrow. Uniaxial compression testing was performed with displacement control. After mechanical testing, each specimen was ashed in a furnace. Four of the samples were exemplarily scanned using micro-computed tomography (μCT) and three dimensional models of the cancellous bones were reconstructed for finite element simulation. The mechanical properties and the failure modes obtained from numerical simulations were then compared to the experiments. Correlations between microarchitectural parameters to the mechanical properties and failure modes were then made. The Young's modulus correlates well with the bone volume fraction with R² = 0.615 and P value 0.013. Three different types of failure modes of cancellous bone were observed: oblique fracture (21.7%), perpendicular global fracture (47.8%), and scattered localised fracture (30.4%). However, no correlations were found between the failure modes to the morphological parameters. The percentage of error between computer predictions and the actual experimental test was from 6 to 12%. These mechanical properties and information on failure modes can be used for the development of synthetic cancellous bone.
Understanding the pathomechanics involved in rheumatoid arthritis (RA) of the wrist provides valuable information, which will invariably allow various therapeutic possibilities to be explored. The computational modelling of this disease permits the appropriate simulation to be conducted seamlessly. A study that underpins the fundamental concept that produces the biomechanical changes in a rheumatoid wrist was thus conducted through the use of finite element method. The RA model was constructed from computed tomography datasets, taking into account three major characteristics: synovial proliferation, cartilage destruction and ligamentous laxity. As control, a healthy wrist joint model was developed in parallel and compared. Cartilage was modelled based on the shape of the articulation while the ligaments were modelled with linear spring elements. A load-controlled analysis was performed simulating physiological hand grip loading conditions. The results demonstrated that the diseased model produced abnormal wrist extension and stress distribution as compared to the healthy wrist model. Due to the weakening of the ligaments, destruction of the cartilage and lower bone density, the altered biomechanical stresses were particularly evident at the radioscaphoid and capitolunate articulations which correlate to clinical findings. These results demonstrate the robust finding of the developed RA wrist model, which accurately predicted the pathological process.
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.
Purpose: Some wheelchair users continue to struggle in maneuvering a wheelchair and navigating through manual doors. Several smart wheelchairs and robotic manipulators were developed to minimize such challenges facing disabled people. Disappointingly, a majority of these high-tech solutions are restricted to laboratories and are not extensively available as commercial products. Previously, a low-tech wheelchair accessory (arc-shaped with many wheels) for pushing doors was modelled and simulated. This work demonstrates the fabrication and testing of the first-generation prototype of the accessory.Materials and methods: The accessory has side portions with a straight arrangement of wheels and a front portion with a straight-arc-straight arrangement of wheels. The accessory was fabricated using conventional manufacturing, off-the-shelf components, and 3D printed ABS fasteners. Stress analysis simulations were done for the fasteners that attach the front accessory to the wheelchair frame. The proof-of-concept of the prototype installed onto a powered wheelchair was tested with a door and an obstacle, each with ∼50 N resistance force.Results: Prototype tests demonstrate the ability of the accessory along with the mechanical robustness of the 3D printed fasteners to push open doors allowing easy navigation through doors and to push/glide against obstacles. The accessory is foldable and detachable.Conclusion: The low-cost of the accessory makes it affordable to many users intending to improve their quality of life. The current study provides an engineering perspective of the accessory, and a clinical perspective is crucial. Other potential applications of the wheelchair accessory include use with scooters, walkers and stretchers.Implications for rehabilitationLow-cost, low-tech accessory is foldable and detachable.Accessory is effective for pushing doors and pushing/gliding against obstacles.Protective nature of the front accessory could prove highly beneficial to some wheelchair users.
The stiffness response or load-deformation/displacement behavior is the most important mechanical behavior that frequently being utilized for validation of the mathematical-physical models representing the mechanical behavior of solid objects in numerical method, compared to actual experimental data. This numerical study aims to investigate the linear-nonlinear stiffness behavior of carbon fiber-reinforced polymer (CFRP) composites at material and structural levels, and its dependency to the sets of individual/group elastic and damage model parameters. In this regard, a validated constitutive damage model, elastic-damage properties as reference data, and simulation process, that account for elastic, yielding, and damage evolution, are considered in the finite element model development process. The linear-nonlinear stiffness responses of four cases are examined, including a unidirectional CFRP composite laminate (material level) under tensile load, and also three multidirectional composite structures under flexural loads. The result indicated a direct dependency of the stiffness response at the material level to the elastic properties. However, the stiffness behavior of the composite structures depends both on the structural configuration, geometry, lay-ups as well as the mechanical properties of the CFRP composite. The value of maximum reaction force and displacement of the composite structures, as well as the nonlinear response of the structures are highly dependent not only to the mechanical properties, but also to the geometry and the configuration of the structures.
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.