Hip fracture has been identified as one of the main health problems in the elderly. To improve the accuracy in assessing subject-specific hip fracture risk, this study proposed normalized cortical bone thickness (NCBT) estimated from patient’s hip DXA as an alternative predictor of hip fracture risk. Hip fracture risk index (HFRI) derived from DXA-based finite element model was utilized as a baseline for evaluating the effectiveness of NCBT in predicting hip fracture risk. It was found that NCBT at the lateral side of the narrowest femoral neck had the strongest correlation with femoral neck HFRI among the six locations of the proximal femur. This study suggests that it is possible to use NCBT as a surrogate for a quick evaluation of hip fracture risk. Yet its clinical performances such as sensitivity to therapy effectiveness and the ability to discriminate clinical fracture cases will be investigated in a future study.

In this thesis, a load path caused by arching action in reinforced concrete slabs is described and validated using a three-dimensional model. Currently, the CHBDC enforces a 4 meter girder spacing requirement in the design of deck slabs. The aim of this thesis is to investigate the load path induced by arching action in deck slabs with a wide range of girder spacing. To do this, a two-dimensional model was developed to examine the path of horizontal stress and was validated using a FEM. A parametric study showed that girder spacing does not affect the development of restraining stress while cantilever width does. As well, cracking of the slab is necessary for arching action to occur. These results help with future development of a rational model to be used by bridge designers.

Dowel-type steel to timber connections are commonly used in timber structure. The load carrying capacity and the stress distribution within the connection area are complicated and the failure behavior of a connection depends on many parameters. The main purpose of this thesis was to verify, using the data obtained from previous experiments, the conventional design method of European Code 5(EC5) (hand calculation) for dowel type joints subjected to pure bending moment and other alternative design methods based on the finite element method (FEM) including the use of the mean stress approach and the extended finite element method (XFEM). Finite element models were created in the software ABAQUS. The models were then used to predict the load bearing capacity and compare this to the experimental results. In addition parametric studies were performed with modifications of material properties and other parameters. The closest prediction in relation to the test results was obtained using XFEM where the predicted capacity was 3.82% larger than the experimental result. An extension of the mean stress method going from a 2D-formulation to a 3D-formulation was verified as well. A general conclusion drawn from this work is that the numerical modelling approaches used should also be suitable for application to complex connections and situations involving other loading situations than pure tension.

The main aim of this research is directed towards the study of thin-walled rectangular planform silos with a view to maximising their structural efficiency. In thin plates of the type making up the wall, membrane action may increase the load carrying capability and current design guides make no account of this. Designing rectangular silos with this in mind can lead to significant structural savings. The core of the research involves using the finite element method to study the patterns of pressure exerted by the weight of a granular bulk solid on the walls of the silo structure. The stored granular solid must use an elastic-plastic material law in order to account for large deformations that can occur in a thin-walled structure. The need for this type of constitutive law led to the investigation of bulk solid properties and shows that parameters that have previously been used to categorise bulk solids may not be sufficient to describe all aspects of their behaviour. The finite element model created uses material constitutive laws that can be found in a number of packages. The required granular material parameters can be determined from a number of simple tests. This approach aims to enable engineers to routinely use similar models when designing silos. The results obtained from the finite element model exhibited some anomalies that had been observed in previous work. These were mainly apparent in the form of localised pressure peaks near the base of the model. These effects were investigated and possible mechanisms that lead to them were proposed. The results from the finite element model were compared to previous experimental work and existing theories. The model was then used to conduct parametric surveys on square and rectangular planform silos and the distribution of pressure across the wall compared to previous predictive models. Finally, a scale thin-walled metal silo was constructed and pressure measurements on filling with pea gravel made. These are compared to predictions made by the finite element model.

Concrete-filled tubular (CFT) columns have been widely adopted in building structures owing to their superior structural performance, such as enhanced load bearing capacity, compared to hollow tubes. Circular, square and rectangular hollow sections are most commonly used in the past few decades. Elliptical hollow section (EHS) available recently is regarded as a new cross-section for the CFT columns due to its attractive appearance, optional orientation either on major axis or minor axis and improved structural efficiency. The state of the research in terms of elliptical columns, tubular joints between EHSs and connections with CFT columns, etc., are reviewed in this thesis, showing a lack of investigations on EHSs, especially on beam to elliptical column connections which are essential in framed structures. The structural behaviour of elliptical column to I-beam connections under bending is studied in this thesis to fill the research gap. Overall ten specimens with various joint assemblies were tested to failure to highlight the benefits of adopting concrete infill and stiffeners in the columns. A three-dimensional finite element model developed by using ABAQUS software is presented and verified against obtained experimental results, which shows acceptable accuracy and reliability in predicting failure modes of the connections and their moment capacities. Parametric studies were performed to access the main parameters that affecting the bending behaviour of the connections. A simple hand calculation method in terms of ultimate moment capacity is proposed according to experiments conducted for connections with concrete-filled columns.

TDCS is a form of non-invasive neurostimulation that is comprised of injection of current via strategically placed scalp electrodes into targeted areas of the brain. TDCS has shown therapeutic benefit for numerous clinical applications. This technique has not however been widely adopted due to high variability of response to the stimulation. Current state of the art methods for optimizing tDCS are based on FEM models that generally model tissue as isotropic and homogeneous and do not take into account inter subject variability of head tissue electrical properties. We therefore develop an in-vivo method of measuring and analyzing bioimpedance spectroscopy measurements of the head to estimate change to tDCS dose in neural tissues for different subjects. Finite element simulations are implemented on a realistic MRI derived head model. 5\% random Gaussian noise is added. Experimental bioimpedance measurements are taken of the heads of 8 subjects. We simulate sensitivity distribution and impedance for a variety of 2 and 4 electrode configurations over a wide frequency range. We also extract Cole parameters and implement PCA on simulated and experimental impedance. We demonstrate that the Cole model of the head can be accurately approximated by the sum in series of Cole systems of each tissue. Comparison of Cole parameters from various simulated electrode configurations show statistical differences (paired t test $p<.05$). PCA shows that close to 100\% of the variance between two impedance spectra is described along a single principal component. Variation described by the second principal component increases as a function of increasing inter electrode gap which may be related to changes in dose. FEM and experimentally derived Cole parameters show different trends for various electrode configurations, good agreement is however shown for the PCA results. The outcome of this research may lead to a higher tDCS efficacy by improving standardization and control of stimulation by relation of dose and bioimpedance spectra characteristics.

The human body has a unique mechanism for perceiving surrounding temperatures. When an object is in contact with the skin, we do not feel its temperature. Instead, we feel the temperature change that is caused on our skin by that object. The faster the heat is transferred, the more intense the thermal sensation is. In this dissertation, a new dynamic thermal display method, where different rates of warm and cold are applied on the skin to generate a unique sensation, is presented. The new method can be related in a wide range of applications including thermal haptics and virtual reality. To understand the perception of temperature and the general thermal state of the human body, the first aspect of this dissertation focuses on investigating the interaction between temperature change and perception on a large scale. Three field surveys were conducted inside airconditioned buildings to investigate the change in the thermal state and temperature perception of occupants when the room temperature changes. The results showed that the participants’ prediction of constant operating temperature was poor, however, their prediction was significantly improved when temperature changes were presented. In order to more accurately investigate the perception of temperature on the skin, a new thermal display method using multiple-channel thermal actuators was developed. The principle of this method is to apply slow and fast rates of temperature change simultaneously on the skin. The slowly changing temperatures are below the perceptual threshold of the thermal receptors, therefore will not be detected whereas the quickly changing temperatures are above the perceptual threshold, hence, will be detected. The idea here is to keep the average surface temperature of the skin constant, however a person will perceive a sensation of continuous cooling. This method was tested through a series of experiments, and the results showed that it is capable of generating a continuous cooling sensation without changing the average temperature of the stimulation area. Multiple variations of this method were tested including different heating and cooling rates of change, different skin locations and patterns of stimuli. Also, a continuous warming was generated using similar concept. To further investigate the temperature distribution that is caused by this method and its effect on the skin, a computational simulation was conducted. An approximate model of the skin was used to monitor its surface temperature and record the temperatures in the stimulation area when the continuous cooling method is applied. The results of the simulation showed that the temperature under the surface of the stimulation area was affected by the continuous cooling method that was applied on the skin model, however this method did not affect the average surface temperature of the skin. These findings may later determine the efficiency and intensity of the method of continuous cooling, and allow us to investigate different technically challenging variations of this method.

Thermal sensation is one of the most dynamic stimulus-response systems in the human body. It is relied upon for safety, comfort and general equilibrium of the human body. Thermal sensation is dependent upon many variables such as area of effected skin, rate of temperature change and location of stimulation. It has been shown that certain rates of change can intensify the sensation of heating or cooling. Conversely, sufficiently low rates of change can go undetected by the skin. As such, the thermal response system can be manipulated by the proper combination of applied hot and cold stimuli. Previous research has shown that through precise application of an asymmetrically heated and cooled thermal display, a sensation of constant cooling can be perceived. This thesis seeks to (1) explore the heat flux characteristics of the thermal display through the use of computer simulations, (2) test a hypothesis about the relationship between thermal sensation and heat flux and (3) examine modifications of the thermal display patterns with the intention of producing more intense thermal sensations. To characterize the heat flux patterns produced by the thermal display, finite element simulations, performed using commercially available software ANSYS©. Simulations are conducted on individual heating and cooling rates to examine the expected values of heat flux as temperatures approach and diverge from skin temperature. Evaluated in the cylindrical coordinate system (axial, angular and radial), the simulations showed a slight nonlinear heat flux generation at the beginning of heating and cooling, but after the initial transient period, this gave way to a strong linear generation of increasing or decreasing heat flux. Simulations were performed that represent the physical experiments implemented in pre- vious research. These simulations were done in two parts: the first examines a small subregion with finer detail on the area between heating and cooling stimuli, the second is a larger scale examination of the heat flux profile of the thermal display. Initially it was observed that directly under the thermal stimulus, in the radial direction, the heat flux was almost perfectly in-phase with the oscillation of temperature whereas between the stimuli, it was nearly 180 degrees out of phase. The heat flux in the axial and angular directions under the thermal stimulus were negligible. Additionally, between stimuli, the values were nearly 180 degrees out of phase with temperature. Additionally, it was observed that the heat flux profiles for all patterns used in the thermal display were approximately identical. From the data gathered by the simulations in conjunction with the thermal sensation data from previous research, a linear relationship is hypothesized that relates these two quantities. This relationship was then used to determine the theoretical thermal sensations of newly developed thermal display patterns in order to determine which are best suited for future physical experimentation.

A finite element model of fibre metal laminate (FML) fuselage was developed in order to evaluate its impact response under survivable crash event. To create a reliable crash finite element (FE) model of FML fuselage, a ‘building block approach’ is adapted. It involves a series of validation and verification tasks in order to establish reliable material and damage models, verified impact model with structural instability and large displacement and verified individual fuselage structure under crash event. This novel development methodology successfully produced an FE model to simulate crash of both aluminium alloy and FML fuselage under survivable crash event using ABAQUS/Explicit. On the other hand, this allows the author to have privilege to evaluate crashworthiness of fuselage that implements FML fuselage skin for the whole fuselage section for the first time in aircraft research field and industry. The FE models consist of a two station fuselage section with one meter longitudinal length which is based on commercial Boeing 737 aircraft. For FML fuselage, the classical aluminium alloy skin was replaced by GLARE grade 5-2/1. The impact response of both fuselages was compared to each other and the results were discussed in terms of energy dissipation, crushing distance, failure modes, failure mechanisms and acceleration response at floor-level. Overall, it was observed that FML fuselage responded similarly to aluminium alloy fuselage with some minor differences which conclusively gives great confidence to aircraft designer to use FML as fuselage skin for the whole fuselage section. In terms of crushing distance, FML fuselage skin contributed to the failure mechanisms of the fuselage section that lead to higher crushing distance than in aluminium alloy fuselage. The existence of various failure modes within FML caused slight differences from the aluminium fuselage in terms of deformation process and energy dissipation. These complex failure modes could potentially be manipulated to produce future aircraft structure with better crashworthiness performance.

With massive infrastructures built in US, timely condition assessment of these infrastructures becomes critical to daily traffic and economics. Due to high cost, long time consumption of direct condition assessment methods, such as closing traffic for sensor installation and monitoring, indirect bridge monitoring has become a promising method. However, the technology is in its initial stage and needs substantial refinement. In this research, virtual simulation approaches, both in 2D and 3D, are used to model the bridge and vehicle interaction through ABAQUS. Artificial Damages were embedded to the model according to different locations and different levels of intensities. With the modelled outcomes, the hypothesis of identifying damages through the responses of the vehicle will be tested. From the simulated vehicle responses, bridge frequencies and damage locations and sizes could be identified accurately through short time flourier transformation and mode shape difference.