CHARACTERIZATION OF THERMAL AND MECHANICAL PROPERTIES OF POLYPROPYLENE–BASED COMPOSITES FOR FUEL CELL BIPOLAR PLATES AND DEVELOPMENT OF EDUCATIONAL TOOLS IN HYDROGEN AND FUEL CELL TECHNOLOGIES

In this project we developed conductive thermoplastic resins by adding varying amounts of three different carbon fillers: carbon black (CB), synthetic graphite (SG) and multi–walled carbon nanotubes (CNT) to a polypropylene matrix for application as fuel cell bipolar plates. This component of fuel cells provides mechanical support to the stack, circulates the gases that participate in the electrochemical reaction within the fuel cell and allows for removal of the excess heat from the system. The materials fabricated in this work were tested to determine their mechanical and thermal properties. These materials were produced by adding varying amounts of single carbon fillers to a polypropylene matrix (2.5 to 15 wt.% Ketjenblack EC-600 JD carbon black, 10 to 80 wt.% Asbury Carbons’ Thermocarb TC-300 synthetic graphite, and 2.5 to 15 wt.% of Hyperion Catalysis International’s FIBRILTM multi-walled carbon nanotubes) In addition, composite materials containing combinations of these three fillers were produced. The thermal conductivity results showed an increase in both through–plane and in–plane thermal conductivities, with the largest increase observed for synthetic graphite. The Department of Energy (DOE) had previously set a thermal conductivity goal of 20 W/m·K, which was surpassed by formulations containing 75 wt.% and 80 wt.% SG, yielding in–plane thermal conductivity values of 24.4 W/m·K and 33.6 W/m·K, respectively. In addition, composites containing 2.5 wt.% CB, 65 wt.% SG, and 6 wt.% CNT in PP had an in–plane thermal conductivity of 37 W/m·K. Flexural and tensile tests were conducted. All composite formulations exceeded the flexural strength target of 25 MPa set by DOE. The tensile and flexural modulus of the composites increased with higher concentration of carbon fillers. Carbon black and synthetic graphite caused a decrease in the tensile and flexural strengths of the composites. However, carbon nanotubes increased the composite tensile and flexural strengths. Mathematical models were applied to estimate through–plane and in–plane thermal conductivities of single and multiple filler formulations, and tensile modulus of single–filler formulations. For thermal conductivity, Nielsen’s model yielded accurate thermal conductivity values when compared to experimental results obtained through the Flash method. For prediction of tensile modulus Nielsen’s model yielded the smallest error between the predicted and experimental values. The second part of this project consisted of the development of a curriculum in Fuel Cell and Hydrogen Technologies to address different educational barriers identified by the Department of Energy. By the creation of new courses and enterprise programs in the areas of fuel cells and the use of hydrogen as an energy carrier, we introduced engineering students to the new technologies, policies and challenges present with this alternative energy. Feedback provided by students participating in these courses and enterprise programs indicate positive acceptance of the different educational tools. Results obtained from a survey applied to students after participating in these courses showed an increase in the knowledge and awareness of energy fundamentals, which indicates the modules developed in this project are effective in introducing students to alternative energy sources.

Improvement of termite resistance, dimensional stability and mechanical properties of pine wood by paraffin impregnation

Paraffin has been used as surface protection of wood throughout the ages but its use for impregnation to improve wood resistance to biodegradation is recent. This study determined the main improvements on wood properties with paraffin impregnation. Healthy Pinus pinaster Ait. wood was impregnated with paraffin at different levels using a hot–cold process. Weight gain, equilibrium moisture content and dimensional stability (ASE) at 35 and 65 % relative humidity, termite durability against Reticulitermes grassei (Clément), bending strength, bending stiffness (MOE) and Janka hardness were determined. Density increased from 0.57 to 0.99, ASE ranged between 38–96 % and 16–71 % for 35 and 65 % relative humidity, respectively. Equilibrium moisture content decreased from 9.9 and 12.0 % to 0.8 and 3.6 % for 35 and 65 % relative humidity. Termite durability improved from level 4 to level 3 of attack, and higher termite mortality was found in treated wood (52 % against 17 %). Bending strength (MOR) increased with paraffin weight gain, reaching a 39 % increase. MOE also increased by about 13 % for wood with a weight gain around 80 %. Janka hardness increased significantly reaching about 40 % for wood with 80 % weight gain. Paraffin impregnated wood has improved properties with regard to equilibrium moisture content, dimensional stability and density, bending strength and Janka hardness, and resistance against termites.

Evaluation of mechanical properties of composite of pinhão with polyrethane derived from castor oil.

2016

Investigation on mechanical behaviour of AM60 magnesium alloys

Purpose: In this work, tension, impact, bend and fatigue tests were conducted in an AM60 magnesium alloy. The effects of environmental temperature and loading rates on impact and tension behavior of the alloy were also investigated. Design/methodology/approach: The tests were conducted using an Instron universal testing machine. The loading speed was changed from 1 mm/min to 300 mm/min to gain a better understanding of the effect of strain rate. To understand the failure behavior of this alloy at different environmental temperatures, Charpy impact test was conducted in a range of temperatures (-40~35°C). Plane strain fracture toughness (KIC) was evaluated using compact tension (CT) specimen. To gain a better understanding of the failure mechanisms, all fracture surfaces were observed using scanning electron microscopy (SEM). In addition, fatigue behavior of this alloy was estimated using tension test under tension-tension condition at 30 Hz. The stress amplitude was selected in the range of 20~50 MPa to obtain the S-N curve. Findings: The tensile test indicated that the mechanical properties were not sensitive to the strain rates applied (3.3x10-4~0.1) and the plastic deformation was dominated by twining mediated slip. The impact energy is not sensitive to the environmental temperature. The plane strain fracture toughness and fatigue limit were evaluated and the average values were 7.6 MPa.m1/2 and 25 MPa, respectively. Practical implications: Tested materials AM60 Mg alloy can be applied among others in automotive industry aerospace, communication and computer industry. Originality/value: Many investigations have been conducted to develop new Mg alloys with improved stiffness and ductility. On the other hand, relatively less attention has been paid to the failure mechanisms of Mg alloys, such as brittle fracture and fatigue, subjected to different environmental or loading conditions. In this work, tension, impact, bend and fatigue tests were conducted in an AM60 magnesium alloy.

Mechanics of deformable glacier beds

My focus in this thesis is to contribute to a more thorough understanding of the mechanics of ice and deformable glacier beds. Glaciers flow under their own weight through a combination of deformation within the ice column and basal slip, which involves both sliding along and deformation within the bed. Deformable beds, which are made up of unfrozen sediment, are prevalent in nature and are often the primary contributors to ice flow wherever they are found. Their granular nature imbues them with unique mechanical properties that depend on the granular structure and hydrological properties of the bed. Despite their importance for understanding glacier flow and the response of glaciers to changing climate, the mechanics of deformable glacier beds are not well understood.

Our general approach to understanding the mechanics of bed deformation and their effect on glacier flow is to acquire synoptic observations of ice surface velocities and their changes over time and to use those observations to infer the mechanical properties of the bed. We focus on areas where changes in ice flow over time are due to known environmental forcings and where the processes of interest are largely isolated from other effects. To make this approach viable, we further develop observational methods that involve the use of mapping radar systems. Chapters 2 and 5 focus largely on the development of these methods and analysis of results from ice caps in central Iceland and an ice stream in West Antarctica. In Chapter 3, we use these observations to constrain numerical ice flow models in order to study the mechanics of the bed and the ice itself. We show that the bed in an Iceland ice cap deforms plastically and we derive an original mechanistic model of ice flow over plastically deforming beds that incorporates changes in bed strength caused by meltwater flux from the surface. Expanding on this work in Chapter 4, we develop a more detailed mechanistic model for till-covered beds that helps explain the mechanisms that cause some glaciers to surge quasi-periodically. In Antarctica, we observe and analyze the mechanisms that allow ocean tidal variations to modulate ice stream flow tens of kilometers inland. We find that the ice stream margins are significantly weakened immediately upstream of the area where ice begins to float and that this weakening likely allows changes in stress over the floating ice to propagate through the ice column.

Anisotropic behavior in plasticity and ductile fracture of an aluminum alloy

Anisotropic mechanical behavior is investigated for an aluminum alloy of 6K21-IH T4 both in plastic deformation and ductile fracture. Anisotropic plastic deformation is characterized by uniaxial tensile tests of dog-bone specimens, while anisotropy in ductile fracture is illustrated with specimens with a central hole, notched specimens and shear specimens. All these specimens are cut off at every 15º from the rolling direction. The r-values and uniaxial tensile yield stresses are measured from the tensile tests of dog-bone specimens. Then the anisotropic plasticity is modeled by a newly proposed J2-J3 criterion under non-associate flow rule (non-AFR). The testing processes of specimens for ductile fracture analysis are simulated to extract the maximum plastic strain at fracture strokes as well as the evolution of the stress triaxiality and the Lode parameter in different testing directions. The measured fracture behavior is described by a shear-controlled ductile fracture criterion proposed by Lou et al. (2014. Modeling of shear ductile fracture considering a changeable cut-off value for stress triaxiality. Int. J. Plasticity 54, 56-80) for different loading directions. It is demonstrated that the anisotropic plastic deformation is described by the J2-J3 criterion with high accuracy in various loading conditions including shear, uniaxial tension and plane strain tension. Moreover, the anisotropy in ductile fracture is not negligible and cannot be modeled by isotropic ductile fracture criteria. Thus, an anisotropic model must be proposed to accurately illustrate the directionality in ductile fracture.

Compósitos poliméricos a base de fibras de licuri: efeitos da hibridização e do envelhecimento ambiental acelerado

The gradual replacement of conventional materials by the ones called composite materials is becoming a concern about the response of these composites against adverse environmental conditions, such as ultraviolet radiation, high temperature and moist. Also the search for new composite using natural fibers or a blend of it with synthetic fibers as reinforcement has been studied. In this sense, this research begins with a thorough study of microstructural characterization of licuri fiber, as a proposal of alternative reinforcement to polymeric composites. Thus, a study about the development of two composite laminates was done. The first one, involving only the fiber of licuri and the second comprising a hybrid composite based of fiber glass E and the fiber of licuri, in order to know the performance of the fiber when of fiber across the hybridization process. The laminates were made in the form of plates using the tereftálica ortho-polyester resin as matrix. The composite laminate made only by licuri fiber had two reinforcing fabric layers of unidirectional licuri and the hybrid composite had two reinforcing layers of unidirectional licuri fabric and three layers of fiber short glass-E mat. Finally, both laminates was exposed to aging acceleration in order to study the influence of environmental degradation involving the mechanical properties and fracture characteristics thereof. Regarding the mechanical properties of composites, these were determined through uniaxial tensile tests, uniaxial compression and three bending points for both laminates in original state, and uniaxial tensile tests and three bending points after accelerated aging. As regards the study of structural degradation due to aging of the laminates, it was carried out based on microscopic analysis and microstructure, as well as measuring weight loss. The characteristics of the fracture was performed by macroscopic and microscopic (optical and SEM) analysis. In general, the laminated composites based on fiber licuri showed some advantages in their responses to environmental aging. These advantages are observed in the behavior related to stiffness as well as the microstructural degradation and photo-oxidation processes. However, the structural integrity of this laminate was more affected in case the action of uniaxial tensile loads, where it was noted a lower rate of withholding his last resistance property

Efeitos do envelhecimento ambiental acelerado em compósitos poliméricos

The advantages of the use of vegetable fibers on the synthetic fibers, such as glass fibers, in the reinforcements in composites are: low cost, low density, good tenacity, good thermal properties and reduced use of instruments for their treatment or processing. However, problems related to poor performance of some mechanical natural fibers, have hindered its direct use in structural elements. In this sense, the emergence of alternative materials such as hybrids composites, involving natural and synthetic fibers, has been encouraged by seeking to improve the performance of structural composites based only on natural fibers. The differences between the physical, chemical and mechanical properties of these fibers, especially facing the adverse environmental conditions such as the presence of moisture and ultraviolet radiation, is also becoming a concern in the final response of these composites. This piece of research presents a comparative study of the strength and stiffness between two composite, both of ortoftalic polyester matrix, one reinforced with fibers of glass-E (CV) and other hybrid reinforced with natural fibers of curauá and fiberglass-E (CH). All the comparative study is based on the influence of exposure to UV rays and steam heated water in composites, simulating the aging environment. The conditions for the tests are accelerated through the use of the aging chamber. The composites will be evaluated through tests of uniaxial static mechanical traction and bending on three points. The composite of glass fiber and hybrid manufacturing industry are using the rolling manual (hand lay-up) and have been developed in the form of composites. All were designed to meet possible structural applications such as tanks and pipes. The reinforcements used in composites were in the forms of short fiber glass-E quilts (450g/m2 - 5cm) of continuous wires and fuses (whose title was of 0.9 dtex) for the curauá fibers. The results clearly show the influence of aging on the environmental mechanical performance of the composite CV and CH. The issues concerning the final characteristics of the fracture for all types of cargoes studied were also analyzed

Free radical entrapment and crystallinity of resin composites after accelerated aging as a function of the expiration date

This study evaluated the spin concentration and the crystallinity in different classifications of dental composites as a function of the material condition (new, aged and expired). Specimens were obtained according to the factors: composites: Filtek P60, Filtek Z250, Filtek Z350XT, and Filtek Silorane; and material conditions: new, aged, and expired. The syringe composites underwent an accelerated aging protocol (Arrhenius model). The magnetic properties of the composites were characterized using Electron Paramagnetic Resonance (EPR) and the concentration of spins (number of spins/mass) was calculated. The crystallinity of the composites tested was characterized with X-ray diffraction (XRD). Filtek P60 and Filtek Z250 presented similarities in terms of spin concentration and crystallinity, irrespective of the material condition. The aging protocol influenced the composite Filmic Z350XT that exhibited a significant increase in the spin concentration. Besides, lower intensity peaks of the organic matrix and amorphous silica were also observed for both aged and expired Filtek Z350XT. Although a significant lower spin concentration was observed for the silorane composite in comparison to that of the methacrylates, a decrease in the relative intensity of peaks of the amorphous region related to the organic components in the diffractograms was observed. The material conditions tested influence the crystallinity and the magnetic properties of the composites evaluated. (C) 2014 Elsevier Ltd. All rights reserved.

Resin composite characterizations following a simplified protocol of accelerated aging as a function of the expiration date

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Environmental effects on viscoelastic behavior of carbon fiber/PEKK thermoplastic composites

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Mixed Mode Fracture Behaviour of Piezoelectric Materials

Piezoelectrics present an interactive electromechanical behaviour that, especially in recent years, has generated much interest since it renders these materials adapt for use in a variety of electronic and industrial applications like sensors, actuators, transducers, smart structures. Both mechanical and electric loads are generally applied on these devices and can cause high concentrations of stress, particularly in proximity of defects or inhomogeneities, such as flaws, cavities or included particles. A thorough understanding of their fracture behaviour is crucial in order to improve their performances and avoid unexpected failures. Therefore, a considerable number of research works have addressed this topic in the last decades. Most of the theoretical studies on this subject find their analytical background in the complex variable formulation of plane anisotropic elasticity. This theoretical approach bases its main origins in the pioneering works of Muskelishvili and Lekhnitskii who obtained the solution of the elastic problem in terms of independent analytic functions of complex variables. In the present work, the expressions of stresses and elastic and electric displacements are obtained as functions of complex potentials through an analytical formulation which is the application to the piezoelectric static case of an approach introduced for orthotropic materials to solve elastodynamics problems. This method can be considered an alternative to other formalisms currently used, like the Stroh’s formalism. The equilibrium equations are reduced to a first order system involving a six-dimensional vector field. After that, a similarity transformation is induced to reach three independent Cauchy-Riemann systems, so justifying the introduction of the complex variable notation. Closed form expressions of near tip stress and displacement fields are therefore obtained. In the theoretical study of cracked piezoelectric bodies, the issue of assigning consistent electric boundary conditions on the crack faces is of central importance and has been addressed by many researchers. Three different boundary conditions are commonly accepted in literature: the permeable, the impermeable and the semipermeable (“exact”) crack model. This thesis takes into considerations all the three models, comparing the results obtained and analysing the effects of the boundary condition choice on the solution. The influence of load biaxiality and of the application of a remote electric field has been studied, pointing out that both can affect to a various extent the stress fields and the angle of initial crack extension, especially when non-singular terms are retained in the expressions of the electro-elastic solution. Furthermore, two different fracture criteria are applied to the piezoelectric case, and their outcomes are compared and discussed. The work is organized as follows: Chapter 1 briefly introduces the fundamental concepts of Fracture Mechanics. Chapter 2 describes plane elasticity formalisms for an anisotropic continuum (Eshelby-Read-Shockley and Stroh) and introduces for the simplified orthotropic case the alternative formalism we want to propose. Chapter 3 outlines the Linear Theory of Piezoelectricity, its basic relations and electro-elastic equations. Chapter 4 introduces the proposed method for obtaining the expressions of stresses and elastic and electric displacements, given as functions of complex potentials. The solution is obtained in close form and non-singular terms are retained as well. Chapter 5 presents several numerical applications aimed at estimating the effect of load biaxiality, electric field, considered permittivity of the crack. Through the application of fracture criteria the influence of the above listed conditions on the response of the system and in particular on the direction of crack branching is thoroughly discussed.

Coupled model of initiation and propagation of corrosion in reinforced concrete

Corrosion of reinforcing steel in concrete due to chloride ingress is one of the main causes of the deterioration of reinforced concrete structures. Structures most affected by such a corrosion are marine zone buildings and structures exposed to de-icing salts like highways and bridges. Such process is accompanied by an increase in volume of the corrosión products on the rebarsconcrete interface. Depending on the level of oxidation, iron can expand as much as six times its original volume. This increase in volume exerts tensile stresses in the surrounding concrete which result in cracking and spalling of the concrete cover if the concrete tensile strength is exceeded. The mechanism by which steel embedded in concrete corrodes in presence of chloride is the local breakdown of the passive layer formed in the highly alkaline condition of the concrete. It is assumed that corrosion initiates when a critical chloride content reaches the rebar surface. The mathematical formulation idealized the corrosion sequence as a two-stage process: an initiation stage, during which chloride ions penetrate to the reinforcing steel surface and depassivate it, and a propagation stage, in which active corrosion takes place until cracking of the concrete cover has occurred. The aim of this research is to develop computer tools to evaluate the duration of the service life of reinforced concrete structures, considering both the initiation and propagation periods. Such tools must offer a friendly interface to facilitate its use by the researchers even though their background is not in numerical simulation. For the evaluation of the initiation period different tools have been developed: Program TavProbabilidade: provides means to carry out a probability analysis of a chloride ingress model. Such a tool is necessary due to the lack of data and general uncertainties associated with the phenomenon of the chloride diffusion. It differs from the deterministic approach because it computes not just a chloride profile at a certain age, but a range of chloride profiles for each probability or occurrence. Program TavProbabilidade_Fiabilidade: carries out reliability analyses of the initiation period. It takes into account the critical value of the chloride concentration on the steel that causes breakdown of the passive layer and the beginning of the propagation stage. It differs from the deterministic analysis in that it does not predict if the corrosion is going to begin or not, but to quantifies the probability of corrosion initiation. Program TavDif_1D: was created to do a one dimension deterministic analysis of the chloride diffusion process by the finite element method (FEM) which numerically solves Fick’second Law. Despite of the different FEM solver already developed in one dimension, the decision to create a new code (TavDif_1D) was taken because of the need to have a solver with friendly interface for pre- and post-process according to the need of IETCC. An innovative tool was also developed with a systematic method devised to compare the ability of the different 1D models to predict the actual evolution of chloride ingress based on experimental measurements, and also to quantify the degree of agreement of the models with each others. For the evaluation of the entire service life of the structure: a computer program has been developed using finite elements method to do the coupling of both service life periods: initiation and propagation. The program for 2D (TavDif_2D) allows the complementary use of two external programs in a unique friendly interface: • GMSH - an finite element mesh generator and post-processing viewer • OOFEM – a finite element solver. This program (TavDif_2D) is responsible to decide in each time step when and where to start applying the boundary conditions of fracture mechanics module in function of the amount of chloride concentration and corrosion parameters (Icorr, etc). This program is also responsible to verify the presence and the degree of fracture in each element to send the Information of diffusion coefficient variation with the crack width. • GMSH - an finite element mesh generator and post-processing viewer • OOFEM – a finite element solver. The advantages of the FEM with the interface provided by the tool are: • the flexibility to input the data such as material property and boundary conditions as time dependent function. • the flexibility to predict the chloride concentration profile for different geometries. • the possibility to couple chloride diffusion (initiation stage) with chemical and mechanical behavior (propagation stage). The OOFEM code had to be modified to accept temperature, humidity and the time dependent values for the material properties, which is necessary to adequately describe the environmental variations. A 3-D simulation has been performed to simulate the behavior of the beam on both, action of the external load and the internal load caused by the corrosion products, using elements of imbedded fracture in order to plot the curve of the deflection of the central region of the beam versus the external load to compare with the experimental data.

Compósitos poliméricos a base de fibras de licuri: efeitos da hibridização e do envelhecimento ambiental acelerado

The gradual replacement of conventional materials by the ones called composite materials is becoming a concern about the response of these composites against adverse environmental conditions, such as ultraviolet radiation, high temperature and moist. Also the search for new composite using natural fibers or a blend of it with synthetic fibers as reinforcement has been studied. In this sense, this research begins with a thorough study of microstructural characterization of licuri fiber, as a proposal of alternative reinforcement to polymeric composites. Thus, a study about the development of two composite laminates was done. The first one, involving only the fiber of licuri and the second comprising a hybrid composite based of fiber glass E and the fiber of licuri, in order to know the performance of the fiber when of fiber across the hybridization process. The laminates were made in the form of plates using the tereftálica ortho-polyester resin as matrix. The composite laminate made only by licuri fiber had two reinforcing fabric layers of unidirectional licuri and the hybrid composite had two reinforcing layers of unidirectional licuri fabric and three layers of fiber short glass-E mat. Finally, both laminates was exposed to aging acceleration in order to study the influence of environmental degradation involving the mechanical properties and fracture characteristics thereof. Regarding the mechanical properties of composites, these were determined through uniaxial tensile tests, uniaxial compression and three bending points for both laminates in original state, and uniaxial tensile tests and three bending points after accelerated aging. As regards the study of structural degradation due to aging of the laminates, it was carried out based on microscopic analysis and microstructure, as well as measuring weight loss. The characteristics of the fracture was performed by macroscopic and microscopic (optical and SEM) analysis. In general, the laminated composites based on fiber licuri showed some advantages in their responses to environmental aging. These advantages are observed in the behavior related to stiffness as well as the microstructural degradation and photo-oxidation processes. However, the structural integrity of this laminate was more affected in case the action of uniaxial tensile loads, where it was noted a lower rate of withholding his last resistance property

An experimental and finite element investigation of the biomechanics of vertebral compression fractures

Osteoporosis is a disease characterized by low bone mass and micro-architectural deterioration of bone tissue, with a consequent increase in bone fragility and susceptibility to fracture. Osteoporosis affects over 200 million people worldwide, with an estimated 1.5 million fractures annually in the United States alone, and with attendant costs exceeding $10 billion dollars per annum. Osteoporosis reduces bone density through a series of structural changes to the honeycomb-like trabecular bone structure (micro-structure). The reduced bone density, coupled with the microstructural changes, results in significant loss of bone strength and increased fracture risk. Vertebral compression fractures are the most common type of osteoporotic fracture and are associated with pain, increased thoracic curvature, reduced mobility, and difficulty with self care. Surgical interventions, such as kyphoplasty or vertebroplasty, are used to treat osteoporotic vertebral fractures by restoring vertebral stability and alleviating pain. These minimally invasive procedures involve injecting bone cement into the fractured vertebrae. The techniques are still relatively new and while initial results are promising, with the procedures relieving pain in 70-95% of cases, medium-term investigations are now indicating an increased risk of adjacent level fracture following the procedure. With the aging population, understanding and treatment of osteoporosis is an increasingly important public health issue in developed Western countries. The aim of this study was to investigate the biomechanics of spinal osteoporosis and osteoporotic vertebral compression fractures by developing multi-scale computational, Finite Element (FE) models of both healthy and osteoporotic vertebral bodies. The multi-scale approach included the overall vertebral body anatomy, as well as a detailed representation of the internal trabecular microstructure. This novel, multi-scale approach overcame limitations of previous investigations by allowing simultaneous investigation of the mechanics of the trabecular micro-structure as well as overall vertebral body mechanics. The models were used to simulate the progression of osteoporosis, the effect of different loading conditions on vertebral strength and stiffness, and the effects of vertebroplasty on vertebral and trabecular mechanics. The model development process began with the development of an individual trabecular strut model using 3D beam elements, which was used as the building block for lattice-type, structural trabecular bone models, which were in turn incorporated into the vertebral body models. At each stage of model development, model predictions were compared to analytical solutions and in-vitro data from existing literature. The incremental process provided confidence in the predictions of each model before incorporation into the overall vertebral body model. The trabecular bone model, vertebral body model and vertebroplasty models were validated against in-vitro data from a series of compression tests performed using human cadaveric vertebral bodies. Firstly, trabecular bone samples were acquired and morphological parameters for each sample were measured using high resolution micro-computed tomography (CT). Apparent mechanical properties for each sample were then determined using uni-axial compression tests. Bone tissue properties were inversely determined using voxel-based FE models based on the micro-CT data. Specimen specific trabecular bone models were developed and the predicted apparent stiffness and strength were compared to the experimentally measured apparent stiffness and strength of the corresponding specimen. Following the trabecular specimen tests, a series of 12 whole cadaveric vertebrae were then divided into treated and non-treated groups and vertebroplasty performed on the specimens of the treated group. The vertebrae in both groups underwent clinical-CT scanning and destructive uniaxial compression testing. Specimen specific FE vertebral body models were developed and the predicted mechanical response compared to the experimentally measured responses. The validation process demonstrated that the multi-scale FE models comprising a lattice network of beam elements were able to accurately capture the failure mechanics of trabecular bone; and a trabecular core represented with beam elements enclosed in a layer of shell elements to represent the cortical shell was able to adequately represent the failure mechanics of intact vertebral bodies with varying degrees of osteoporosis. Following model development and validation, the models were used to investigate the effects of progressive osteoporosis on vertebral body mechanics and trabecular bone mechanics. These simulations showed that overall failure of the osteoporotic vertebral body is initiated by failure of the trabecular core, and the failure mechanism of the trabeculae varies with the progression of osteoporosis; from tissue yield in healthy trabecular bone, to failure due to instability (buckling) in osteoporotic bone with its thinner trabecular struts. The mechanical response of the vertebral body under load is highly dependent on the ability of the endplates to deform to transmit the load to the underlying trabecular bone. The ability of the endplate to evenly transfer the load through the core diminishes with osteoporosis. Investigation into the effect of different loading conditions on the vertebral body found that, because the trabecular bone structural changes which occur in osteoporosis result in a structure that is highly aligned with the loading direction, the vertebral body is consequently less able to withstand non-uniform loading states such as occurs in forward flexion. Changes in vertebral body loading due to disc degeneration were simulated, but proved to have little effect on osteoporotic vertebra mechanics. Conversely, differences in vertebral body loading between simulated invivo (uniform endplate pressure) and in-vitro conditions (where the vertebral endplates are rigidly cemented) had a dramatic effect on the predicted vertebral mechanics. This investigation suggested that in-vitro loading using bone cement potting of both endplates has major limitations in its ability to represent vertebral body mechanics in-vivo. And lastly, FE investigation into the biomechanical effect of vertebroplasty was performed. The results of this investigation demonstrated that the effect of vertebroplasty on overall vertebra mechanics is strongly governed by the cement distribution achieved within the trabecular core. In agreement with a recent study, the models predicted that vertebroplasty cement distributions which do not form one continuous mass which contacts both endplates have little effect on vertebral body stiffness or strength. In summary, this work presents the development of a novel, multi-scale Finite Element model of the osteoporotic vertebral body, which provides a powerful new tool for investigating the mechanics of osteoporotic vertebral compression fractures at the trabecular bone micro-structural level, and at the vertebral body level.