Finite element modelling of load bearing cold-formed steel wall systems under fire conditions

Light Gauge Steel Framing (LSF) walls are made of cold-formed, thin-walled steel lipped channel studs with plasterboard linings on both sides. However, these thin-walled steel sections heat up quickly and lose their strength under fire conditions despite the protection provided by plasterboards. A new composite wall panel was recently proposed to improve the fire resistance rating of LSF walls, where an insulation layer was used externally between the plasterboards on both sides of the wall frame instead of using it in the cavity. A research study using both fire tests and numerical studies was undertaken to investigate the structural and thermal behaviour of load bearing LSF walls made of both conventional and the new composite panels under standard fire conditions and to determine their fire resistance rating. This paper presents the details of finite element models of LSF wall studs developed to simulate the structural performance of LSF wall panels under standard fire conditions. Finite element analyses were conducted under both steady and transient state conditions using the time-temperature profiles measured during the fire tests. The developed models were validated using the fire test results of 11 LSF wall panels with various plasterboard/insulation configurations and load ratios. They were able to predict the fire resistance rating within five minutes. The use of accurate numerical models allowed the inclusion of various complex structural and thermal effects such as local buckling, thermal bowing and neutral axis shift that occurred in thin-walled steel studs under non-uniform elevated temperature conditions. Finite element analyses also demonstrated the improvements offered by the new composite panel system over the conventional cavity insulated system.

Ultimate flexural limit states analysis of prestressed concrete sleeper

Railway is one of the most important, reliable and widely used means of transportation, carrying freight, passengers, minerals, grains, etc. Thus, research on railway tracks is extremely important for the development of railway engineering and technologies. The safe operation of a railway track is based on the railway track structure that includes rails, fasteners, pads, sleepers, ballast, subballast and formation. Sleepers are very important components of the entire structure and may be made of timber, concrete, steel or synthetic materials. Concrete sleepers were first installed around the middle of last century and currently are installed in great numbers around the world. Consequently, the design of concrete sleepers has a direct impact on the safe operation of railways. The "permissible stress" method is currently most commonly used to design sleepers. However, the permissible stress principle does not consider the ultimate strength of materials, probabilities of actual loads, and the risks associated with failure, all of which could lead to the conclusion of cost-ineffectiveness and over design of current prestressed concrete sleepers. Recently the limit states design method, which appeared in the last century and has been already applied in the design of buildings, bridges, etc, is proposed as a better method for the design of prestressed concrete sleepers. The limit states design has significant advantages compared to the permissible stress design, such as the utilisation of the full strength of the member, and a rational analysis of the probabilities related to sleeper strength and applied loads. This research aims to apply the ultimate limit states design to the prestressed concrete sleeper, namely to obtain the load factors of both static and dynamic loads for the ultimate limit states design equations. However, the sleepers in rail tracks require different safety levels for different types of tracks, which mean the different types of tracks have different load factors of limit states design equations. Therefore, the core tasks of this research are to find the load factors of the static component and dynamic component of loads on track and the strength reduction factor of the sleeper bending strength for the ultimate limit states design equations for four main types of tracks, i.e., heavy haul, freight, medium speed passenger and high speed passenger tracks. To find those factors, the multiple samples of static loads, dynamic loads and their distributions are needed. In the four types of tracks, the heavy haul track has the measured data from Braeside Line (A heavy haul line in Central Queensland), and the distributions of both static and dynamic loads can be found from these data. The other three types of tracks have no measured data from sites and the experimental data are hardly available. In order to generate the data samples and obtain their distributions, the computer based simulations were employed and assumed the wheel-track impacts as induced by different sizes of wheel flats. A valid simulation package named DTrack was firstly employed to generate the dynamic loads for the freight and medium speed passenger tracks. However, DTrack is only valid for the tracks which carry low or medium speed vehicles. Therefore, a 3-D finite element (FE) model was then established for the wheel-track impact analysis of the high speed track. This FE model has been validated by comparing its simulation results with the DTrack simulation results, and with the results from traditional theoretical calculations based on the case of heavy haul track. Furthermore, the dynamic load data of the high speed track were obtained from the FE model and the distributions of both static and dynamic loads were extracted accordingly. All derived distributions of loads were fitted by appropriate functions. Through extrapolating those distributions, the important parameters of distributions for the static load induced sleeper bending moment and the extreme wheel-rail impact force induced sleeper dynamic bending moments and finally, the load factors, were obtained. Eventually, the load factors were obtained by the limit states design calibration based on reliability analyses with the derived distributions. After that, a sensitivity analysis was performed and the reliability of the achieved limit states design equations was confirmed. It has been found that the limit states design can be effectively applied to railway concrete sleepers. This research significantly contributes to railway engineering and the track safety area. It helps to decrease the failure and risks of track structure and accidents; better determines the load range for existing sleepers in track; better rates the strength of concrete sleepers to support bigger impact and loads on railway track; increases the reliability of the concrete sleepers and hugely saves investments on railway industries. Based on this research, many other bodies of research can be promoted in the future. Firstly, it has been found that the 3-D FE model is suitable for the study of track loadings and track structure vibrations. Secondly, the equations for serviceability and damageability limit states can be developed based on the concepts of limit states design equations of concrete sleepers obtained in this research, which are for the ultimate limit states.

Porohyperelastic analysis of single osteocyte using AFM and inverse FEA

Osteocytes are the mature cells and perform as mechanosensors within the bone. The mechanical property of osteocytes plays an important role to fulfill these functions. However, little researches have been done to investigate the mechanical deformation properties of single osteocytes. Atomic Force Microscopy (AFM) is a state-of-art experimental facility for high resolution imaging of tissues, cells and any surfaces as well as for probing mechanical properties of the samples both qualitatively and quantitatively. In this paper, the experimental study based on AFM is firstly used to obtain forceindentation curves of single round osteocytes. The porohyperelastic (PHE) model of a single osteocyte is then developed by using the inverse finite element analysis (FEA) to identify and extract mechanical properties from the experiment results. It has been found that the PHE model is a good candidature for biomechanics studies of osteocytes.

Porohyperelastic analysis of single chondrocyte using AFM and inverse FEA

The aim of this paper is to determine the creep and relaxation responses of single chondrocytes in vitro. Firstly, Atomic Force Microscopy (AFM) was used to obtain the force-indentation curves of single chondrocytes at the strain-rate of 7.05 s-1. This result was then employed in inverse finite element analysis (FEA) using porohyperelastic (PHE) idealization of the cells to determine their mechanical properties. The PHE model results agreed well with AFM experimental data. This PHE model was then utilized to study chondrocyte’s creep and relaxation behaviors. The results revealed that the effect of fluid was predominant for cell’s mechanical behaviors and that the PHE is a good model for biomechanics studies of chondrocytes.

Using patient specific finite element models to understand the biomechanics of paediatric spinal deformity surgery

Adolescent idiopathic scoliosis (AIS) is a spinal deformity, which may require surgical correction by attaching rods to the patient’s spine using screws inserted into the vertebrae. Complication rates for deformity correction surgery are unacceptably high. Determining an achievable correction without overloading the adjacent spinal tissues or implants requires an understanding of the mechanical interaction between these components. We have developed novel patient specific modelling software to create individualized finite element models (FEM) representing the thoracolumbar spine and ribcage of scoliosis patients. We are using these models to better understand the biomechanics of spinal deformity correction.

Stress relaxation analysis of single chondrocytes using porohyperelastic model based on AFM experiments

The aim of this study is to investigate the stress relaxation behavior of single chondrocytes using the Porohyperelastic (PHE) model and inverse Finite Element Analysis (FEA). Firstly, based on Atomic Force Microscopy (AFM) technique, we have found that the chondrocytes exhibited stress relaxation behavior. We explored the mechanism of this stress relaxation behavior and concluded that the intracellular fluid exuding out from the cells during deformation plays the most important role in the stress relaxation. Next, we have applied the inverse FEA technique to determine necessary material parameters for PHE model to simulate this stress relaxation behavior as this model is proven capable of capturing the non-linear behavior and the fluid-solid interaction during the stress relaxation of the single chondrocytes. It is observed that this PHE model can precisely capture the stress relaxation behavior of single chondrocytes and would be a suitable model for cell biomechanics.

Computational strategies for surface fitting using thin plate spline finite element methods

Thin plate spline finite element methods are used to fit a surface to an irregularly scattered dataset [S. Roberts, M. Hegland, and I. Altas. Approximation of a Thin Plate Spline Smoother using Continuous Piecewise Polynomial Functions. SIAM, 1:208--234, 2003]. The computational bottleneck for this algorithm is the solution of large, ill-conditioned systems of linear equations at each step of a generalised cross validation algorithm. Preconditioning techniques are investigated to accelerate the convergence of the solution of these systems using Krylov subspace methods. The preconditioners under consideration are block diagonal, block triangular and constraint preconditioners [M. Benzi, G. H. Golub, and J. Liesen. Numerical solution of saddle point problems. Acta Numer., 14:1--137, 2005]. The effectiveness of each of these preconditioners is examined on a sample dataset taken from a known surface. From our numerical investigation, constraint preconditioners appear to provide improved convergence for this surface fitting problem compared to block preconditioners.

An integrated modeling and risk analysis tool for mine water interactions

Poor mine water management can lead to corporate, environmental and social risks. These risks become more pronounced as mining operations move into areas of water scarcity and/or increase climatic variability while also managing increased demand, lower ore grades and increased strip ratios. Therefore, it is vital that mine sites better understand these risks in order to implement management practices to address them. Systems models provide an effective approach to understand complex networks, particularly across multiple scales. Previous work has represented mine water interactions using systems model on a mine site scale. Here, we expand on that work by present an integrated tool that uses a systems modeling approach to represent mine water interactions on a site and regional scale and then analyses the risks associated with events stemming from those interactions. A case study is presented to represent three indicative corporate, environmental and social risks associated with a mine site that exists in a water scarce region. The tool is generic and flexible, and can be used in many scenarios to provide significant potential utility to the mining industry.

Comparison of eight published static finite element models of the intact lumbar spine : predictive power of models improves when combined together

Finite element (FE) model studies have made important contributions to our understanding of functional biomechanics of the lumbar spine. However, if a model is used to answer clinical and biomechanical questions over a certain population, their inherently large inter-subject variability has to be considered. Current FE model studies, however, generally account only for a single distinct spinal geometry with one set of material properties. This raises questions concerning their predictive power, their range of results and on their agreement with in vitro and in vivo values. Eight well-established FE models of the lumbar spine (L1-5) of different research centres around the globe were subjected to pure and combined loading modes and compared to in vitro and in vivo measurements for intervertebral rotations, disc pressures and facet joint forces. Under pure moment loading, the predicted L1-5 rotations of almost all models fell within the reported in vitro ranges, and their median values differed on average by only 2° for flexion-extension, 1° for lateral bending and 5° for axial rotation. Predicted median facet joint forces and disc pressures were also in good agreement with published median in vitro values. However, the ranges of predictions were larger and exceeded those reported in vitro, especially for the facet joint forces. For all combined loading modes, except for flexion, predicted median segmental intervertebral rotations and disc pressures were in good agreement with measured in vivo values. In light of high inter-subject variability, the generalization of results of a single model to a population remains a concern. This study demonstrated that the pooled median of individual model results, similar to a probabilistic approach, can be used as an improved predictive tool in order to estimate the response of the lumbar spine.

Porohyperelastic finite element model for the kangaroo humeral head cartilage based on experimental study and the consolidation theory

Solid-extracellular fluid interaction is believed to play an important role in the strain-rate dependent mechanical behaviors of shoulder articular cartilages. It is believed that the kangaroo shoulder joint is anatomically and biomechanically similar to human shoulder joint and it is easy to get in Australia. Therefore, the kangaroo humeral head cartilage was used as the suitable tissue for the study in this paper. Indentation tests from quasi-static (10-4/sec) to moderately high strain-rate (10-2/sec) on kangaroo humeral head cartilage tissues were conduced to investigate the strain-rate dependent behaviors. A finite element (FE) model was then developed, in which cartilage was conceptualized as a porous solid matrix filled with incompressible fluids. In this model, the solid matrix was modeled as an isotropic hyperelastic material and the percolating fluid follows Darcy’s law. Using inverse FE procedure, the constitutive parameters related to stiffness, compressibility of the solid matrix and permeability were obtained from the experimental results. The effect of solid-extracellular fluid interaction and drag force (the resistance to fluid movement) on strain-rate dependent behavior was investigated by comparing the influence of constant, strain dependent and strain-rate dependent permeability on FE model prediction. The newly developed porohyperelastic cartilage model with the inclusion of strain-rate dependent permeability was found to be able to predict the strain-rate dependent behaviors of cartilages.

A poroviscohyperelastic model for numerical analysis of mechanical behavior of single chondrocyte

The aim of this paper is to utilize a poroviscohyperelastic (PVHE) model which is developed based on the porohyperelastic (PHE) model to explore the mechanical deformation properties of single chondrocytes. Both creep and relaxation responses are investigated by using FEM models of micropipette aspiration and AFM experiments, respectively. The newly developed PVHE model is compared thoroughly with the SnHS and PHE models. It has been found that the PVHE can accurately capture both creep and stress relaxation behaviors of chondrocytes better than other two models. Hence, the PVHE is a promising model to investigate mechanical properties of single chondrocytes.

Explicit finite element modelling of bridge girder bearing pedestals

Bridge girder bearings rest on pedestals to transfer the loading safely to the pier headstock. In spite of the existence of industry guidelines, due to construction complexities, such guidelines are often overlooked. Further, there is paucity of research on the performance of pedestals, although their failure could cause exorbitant maintenance costs. Although reinforced concrete pedestals are recommended in the industry design guidelines, unreinforced concrete and/ or epoxy glue pedestals are provided due to construction issues; such pedestals fail within a very short period of service. With a view to understanding the response of pedestals subject to monotonic loading, a three-dimensional nonlinear explicit finite element micro-model of unreinforced and reinforced concrete pedestals has been developed. Contact and material nonlinearity have been accounted for in the model. It is shown that the unreinforced concrete pedestals suffer from localised edge stress singularities, the failure of which was comparable to those in the field. The reinforced concrete pedestals, on the other hand, distribute the loading without edge stress singularity, again conforming to the field experience.