Evolutionary design of bone scaffolds with reference to material selection

The favourable scaffold for bone tissue engineering should have desired characteristic features, such as adequate mechanical strength and three-dimensional open porosity, which guarantee a suitable environment for tissue regeneration. In fact, the design of such complex structures like bone scaffolds is a challenge for investigators. One of the aims is to achieve the best possible mechanical strength-degradation rate ratio. In this paper we attempt to use numerical modelling to evaluate material properties for designing bone tissue engineering scaffold fabricated via the fused deposition modelling technique. For our studies the standard genetic algorithm was used, which is an efficient method of discrete optimization. For the fused deposition modelling scaffold, each individual strut is scrutinized for its role in the architecture and structural support it provides for the scaffold, and its contribution to the overall scaffold was studied. The goal of the study was to create a numerical tool that could help to acquire the desired behaviour of tissue engineered scaffolds and our results showed that this could be achieved efficiently by using different materials for individual struts. To represent a great number of ways in which scaffold mechanical function loss could proceed, the exemplary set of different desirable scaffold stiffness loss function was chosen. © 2012 John Wiley & Sons, Ltd.

An implicit RBF meshless method for a fractal mobile/immobile transport model

Fractional differential equation is used to describe a fractal model of mobile/immobile transport with a power law memory function. This equation is the limiting equation that governs continuous time random walks with heavy tailed random waiting times. In this paper, we firstly propose a finite difference method to discretize the time variable and obtain a semi-discrete scheme. Then we discuss its stability and convergence. Secondly we consider a meshless method based on radial basis functions (RBF) to discretize the space variable. By contrast to conventional FDM and FEM, the meshless method is demonstrated to have distinct advantages: calculations can be performed independent of a mesh, it is more accurate and it can be used to solve complex problems. Finally the convergence order is verified from a numerical example is presented to describe the fractal model of mobile/immobile transport process with different problem domains. The numerical results indicate that the present meshless approach is very effective for modeling and simulating of fractional differential equations, and it has good potential in development of a robust simulation tool for problems in engineering and science that are governed by various types of fractional differential equations.

Lateral movements in composite high-rise buildings under seismic action

Daring human nature has already led to the construction of high-rise buildings in naturally challenging geological regions and in worse environments of the world. However; literature review divulges that there is a lag in research of certain generic principles and rules for the prediction of lateral movement in multistorey construction. The present competitive trend orders the best possible used of available construction material and resources. Hence; the mixed used of reinforced concrete with structural steel is gaining prevalence day by day. This paper investigates the effects of Seismic load on composite multistorey building provided with core wall and trusses through FEM modelling. The results showed that increased rigidity corresponds to lower period of vibration and hence higher seismic forces. Since Seismic action is a function of mass and response acceleration, therefore; mass increment generate higher earthquake load and thus cause higher impact base shear and overturning movement. Whereas; wind force depends on building exposed, larger the plan dimension greater is the wind impact. Nonetheless; outriggers trusses noticeably contribute, in improving the serviceability of structure subjected to wind and earthquake forces.

Experimental study on mechanical properties of pumpkin tissue

Purpose: The purpose of this study was to calculate mechanical properties of tough skinned vegetables as a part of Finite Element Modelling (FEM) and simulation of tissue damage during mechanical peeling of tough skinned vegetables. Design/methodology: There are some previous studies on mechanical properties of fruits and vegetables however, behaviour of tissue under different processing operations will be different. In this study indentation test was performed on Peel, Flesh and Unpeeled samples of pumpkin as a tough skinned vegetable. Additionally, the test performed in three different loading rates for peel: 1.25, 10, 20 mm/min and 20 mm/min for flesh and unpeeled samples respectively. The spherical end indenter with 8mm diameter used for the experimental tests. Samples prepare from defect free and ripped pumpkin purchased from local shops in Brisbane, Australia. Humidity and temperature were 20-55% and 20-250C respectively. Findings: Consequently, force deformation and stress and strain of samples were calculated and shown in presented figures. Relative contribution (%) of skin to different mechanical properties is computed and compared with data available from literature. According the results, peel samples had the highest value of rupture force (291N) and as well as highest value of firmness (1411Nm-1). Research limitations/implications: The proposed study focused on one type of tough skinned vegetables and one variety of pumpkin however, more tests will give better understandings of behaviours of tissue. Additionally, the behaviours of peel, unpeeled and flesh samples in different speed of loading will provide more details of tissue damages during mechanical loading. Originality/value: Mechanical properties of pumpkin tissue calculated using the results of indentation test, specifically the behaviours of peel, flesh and unpeeled samples were explored which is a new approach in Finite Element Modelling (FEM) of food processes. Keywords: Finite Element Modelling (FEM), relative contribution, firmness, toughness and rupture force.

Studies on bending limitations for the optimal fit of orthopaedic bone plates

Distal tibial fractures are now commonly treated via intermedullary plate fixation due to higher rates of union and lower rates of postoperative complications. However, patient specific bone morphology demands manual deformation of the plate to ensure appropriate fit along the bone Distal tibial fractures are now commonly treated via intermedullary plate fixation due to higher rates of union and lower rates of postoperative complications. However, patient specific bone morphology demands manual deformation of the plate to ensure appropriate fit along the bone contours, and depending on the material of the plate, different outcomes have been reported along with postoperative complications. A comparative analysis of Stainless Steel 316L and Ti-6Al-4V alloys was carried to estimate the safe bending limit for appropriate fits. The results from the ANSYS FEA simulations were validated with experiments based on ASTM F382-99. It is found that SS316L is better suited for large deformations (up to 16˚ in proximal tip and 7.5˚ in distal end) and Ti for smaller deformation contours (up to 3˚ in proximal tip and 1.8˚ in distal end). The results of this study have profound implications for the choice of plates based on preliminary radiographical fracture examinations to ensure better fixation and higher rates of union of distal tibial fractures.

Experimental investigation of wheel/rail rolling contact at railhead edge

Wheel-rail rolling contact at railhead edge, such as a gap in an insulated rail joint, is a complex problem; there are only limited analytical, numerical and experimental studies available on this problem in the academic literature. This paper describes experimental and numerical investigations of railhead strains in the vicinity of the edge under the contact of a loaded wheel. A full-scale test rig was developed to cyclically apply wheel/rail rolling contact load to the edge zone of the railhead. An image analysis technique was employed to determine the railhead vertical, lateral and shear strain components. The vertical strains determined using the image analysis method have been validated with the strain gauge measurements and used for the calibration of a 3D nonlinear Finite Element Model (FEM) that simulates the wheel/rail contact at the railhead edge and use suitable boundary conditions commensurate to the experimental setup. The FEM was then used to determine other states of strains.

Numerical study of impact forces on railway sleepers under wheel flat

Wheel-rail interaction is one of the most important research topics in railway engineering. It includes track vibration, track impact response and safety of the track. Track structure failures caused by impact forces can lead to significant economic loss for track owners through damage to rails and to the sleepers beneath. The wheel-rail impact forces occur because of imperfections on the wheels or rails such as wheel flats, irregular wheel profile, rail corrugation and differences in the height of rails connected at a welded joint. The vehicle speed and static wheel load are important factors of the track design, because they are related to the impact forces under wheel-rail defects. In this paper, a 3-Dimensional finite element model for the study of wheel flat impact is developed by use of the FEA software package ANSYS. The effects of the wheel flat to impact force on sleepers with various speeds and static wheel loads under a critical wheel flat size are investigated. It has found that both wheel-rail impact force and impact force on sleeper induced by wheel flat are varying nonlinearly by increasing the vehicle speed; both impact forces are nonlinearly and monotonically increasing by increasing the static wheel load. The relationships between both of impact forces induced by wheel flat and vehicles speed or static load are important to the track engineers to improve the design and maintenance methods in railway industry.

A non-linear dynamic analysis of heavy haul track line sleepers under wheel flat

Significant wheel-rail dynamic forces occur because of imperfections in the wheels and/or rail. One of the key responses to the transmission of these forces down through the track is impact force on the sleepers. Dynamic analysis of nonlinear systems is very complicated and does not lend itself easily to a classical solution of multiple equations. Trying to deduce the behaviour of track components from experimental data is very difficult because such data is hard to obtain and applies to only the particular conditions of the track being tested. The finite element method can be the best solution to this dilemma. This paper describes a finite element model using the software package ANSYS for various sized flat defects in the tread of a wheel rolling at a typical speed on heavy haul track. The paper explores the dynamic response of a prestressed concrete sleeper to these defects.

A dynamic wheel-rail impact analysis of railway track under wheel flat by Finite Element Analysis

Wheel–rail interaction is one of the most important research topics in railway engineering. It involves track impact response, track vibration and track safety. Track structure failures caused by wheel–rail impact forces can lead to significant economic loss for track owners through damage to rails and to the sleepers beneath. Wheel–rail impact forces occur because of imperfections in the wheels or rails such as wheel flats, irregular wheel profiles, rail corrugations and differences in the heights of rails connected at a welded joint. A wheel flat can cause a large dynamic impact force as well as a forced vibration with a high frequency, which can cause damage to the track structure. In the present work, a three-dimensional (3-D) finite element (FE) model for the impact analysis induced by the wheel flat is developed by use of the finite element analysis (FEA) software package ANSYS and validated by another validated simulation. The effect of wheel flats on impact forces is thoroughly investigated. It is found that the presence of a wheel flat will significantly increase the dynamic impact force on both rail and sleeper. The impact force will monotonically increase with the size of wheel flats. The relationships between the impact force and the wheel flat size are explored from this finite element analysis and they are important for track engineers to improve their understanding of the design and maintenance of the track system.

Magnetic flux analysis in superconducting Fault Current Limiters under short-circuit conditions

In order to obtain a more compact Superconducting Fault Current limiter (SFCL), a special geometry of core and AC coil is required. This results in a unique magnetic flux pattern which differs from those associated with conventional round core arrangements. In this paper the magnetic flux density within a Fault Current Limiter (FCL) is described. Both experimental and analytical approaches are considered. A small scale prototype of an FCL was constructed in order to conduct the experiments. This prototype comprises a single phase. The analysis covers both the steady state and the short-circuit condition. Simulation results were obtained using commercial software based on the Finite Element Method (FEM). The magnetic flux saturating the cores, leakage magnetic flux giving rise to electromagnetic forces and leakage magnetic flux flowing in the enclosing tank are computed.

A hybrid smoothed finite element method (H-SFEM) to solid mechanics problems

In this paper, a hybrid smoothed finite element method (H-SFEM) is developed for solid mechanics problems by combining techniques of finite element method (FEM) and Node-based smoothed finite element method (NS-FEM) using a triangular mesh. A parameter is equipped into H-SFEM, and the strain field is further assumed to be the weighted average between compatible stains from FEM and smoothed strains from NS-FEM. We prove theoretically that the strain energy obtained from the H-SFEM solution lies in between those from the compatible FEM solution and the NS-FEM solution, which guarantees the convergence of H-SFEM. Intensive numerical studies are conducted to verify these theoretical results and show that (1) the upper and lower bound solutions can always be obtained by adjusting ; (2) there exists a preferable at which the H-SFEM can produce the ultrasonic accurate solution.

Porohyperelastic analysis to explore mechanical properties of chondrocytes using numerical modeling and experiments : a finite element study

There are many continuum mechanical models have been developed such as liquid drop models, solid models, and so on for single living cell biomechanics studies. However, these models do not give a fully approach to exhibit a clear understanding of the behaviour of single living cells such as swelling behaviour, drag effect, etc. Hence, the porohyperelastic (PHE) model which can capture those aspects would be a good candidature to study cells behaviour (e.g. chondrocytes in this study). In this research, an FEM model of single chondrocyte cell will be developed by using this PHE model to simulate Atomic Force Microscopy (AFM) experimental results with the variation of strain rate. This material model will be compared with viscoelastic model to demonstrate the advantages of PHE model. The results have shown that the maximum value of force applied of PHE model is lower at lower strain rates. This is because the mobile fluid does not have enough time to exude in case of very high strain rate and also due to the lower permeability of the membrane than that of the protoplasm of chondrocyte. This behavior is barely observed in viscoelastic model. Thus, PHE model is the better model for cell biomechanics studies.

Visualization of thermal characteristics around the absorber tube of a standard parabolic trough thermal collector by 3D simulation

Three dimensional conjugate heat transfer simulation of a standard parabolic trough thermal collector receiver is performed numerically in order to visualize and analyze the surface thermal characteristics. The computational model is developed in Ansys Fluent environment based on some simplified assumptions. Three test conditions are selected from the existing literature to verify the numerical model directly, and reasonably good agreement between the model and the test results confirms the reliability of the simulation. Solar radiation flux profile around the tube is also approximated from the literature. An in house macro is written to read the input solar flux as a heat flux wall boundary condition for the tube wall. The numerical results show that there is an abrupt variation in the resultant heat flux along the circumference of the receiver. Consequently, the temperature varies throughout the tube surface. The lower half of the horizontal receiver enjoys the maximum solar flux, and therefore, experiences the maximum temperature rise compared to the upper part with almost leveled temperature. Reasonable attributions and suggestions are made on this particular type of conjugate thermal system. The knowledge that gained so far from this study will be used to further the analysis and to design an efficient concentrator photovoltaic collector in near future.

A finite element model of seat cushion indentation with a soft tissue human occupant model

Effective digital human model (DHM) simulation of automotive driver packaging ergonomics, safety and comfort depends on accurate modelling of occupant posture, which is strongly related to the mechanical interaction between human body soft tissue and flexible seat components. This paper presents a finite-element study simulating the deflection of seat cushion foam and supportive seat structures, as well as human buttock and thigh soft tissue when seated. The three-dimensional data used for modelling thigh and buttock geometry were taken on one 95th percentile male subject, representing the bivariate percentiles of the combined hip breadth (seated) and buttock-to-knee length distributions of a selected Australian and US population. A thigh-buttock surface shell based on this data was generated for the analytic model. A 6mm neoprene layer was offset from the shell to account for the compression of body tissue expected through sitting in a seat. The thigh-buttock model is therefore made of two layers, covering thin to moderate thigh and buttock proportions, but not more fleshy sizes. To replicate the effects of skin and fat, the neoprene rubber layer was modelled as a hyperelastic material with viscoelastic behaviour in a Neo-Hookean material model. Finite element (FE) analysis was performed in ANSYS V13 WB (Canonsburg, USA). It is hypothesized that the presented FE simulation delivers a valid result, compared to a standard SAE physical test and the real phenomenon of human-seat indentation. The analytical model is based on the CAD assembly of a Ford Territory seat. The optimized seat frame, suspension and foam pad CAD data were transformed and meshed into FE models and indented by the two layer, soft surface human FE model. Converging results with the least computational effort were achieved for a bonded connection between cushion and seat base as well as cushion and suspension, no separation between neoprene and indenter shell and a frictional connection between cushion pad and neoprene. The result is compared to a previous simulation of an indentation with a hard shell human finite-element model of equal geometry, and to the physical indentation result, which is approached with very high fidelity. We conclude that (a) SAE composite buttock form indentation of a suspended seat cushion can be validly simulated in a FE model of merely similar geometry, but using a two-layer hard/soft structure. (b) Human-seat indentation of a suspended seat cushion can be validly simulated with a simplified human buttock-thigh model for a selected anthropomorphism.

Safety enhancement of water-filled road safety barrier using interaction of composite materials

Road safety barriers are used to redirect traffic at roadside work-zones. When filled with water, these barriers are able to withstand low to moderate impact speeds up to 50kmh-1. Despite this feature, Portable Water-filled barriers (PWFB) face challenges such as large lateral displacements, tearing and breakage during impact; especially at higher speeds. This study explores the use of composite action to enhance the crashworthiness of PWFBs and enable their usage at higher speeds. Initially, energy absorption capability of water in PWFB is investigated. Then, composite action of the PWFB with the introduction of steel frame is considered to evaluate its enhanced impact performance. Findings of the study show that the initial height of the impact must be lower than the free surface level of water in a PWFB in order for the water to provide significant crash energy absorption. In general, an impact of a road barrier with 80% filled is a good estimation. Furthermore, the addition of a composite structure greatly reduces the probability of tearing by decreasing the strain and impact energy transferred to the shell container. This allows the water to remain longer in the barrier to absorb energy via inertial displacements and sloshing response. Information from this research will aid in the design of new generation roadside safety structures aimed to increase safety in modern roadways.