Analysis of the compressive behaviour of the three-dimensional printed porous titanium for dental implants using a modified cellular solid model

A set of cylindrical porous titanium test samples were produced using the three-dimensional printing and sintering method with samples sintered at 900 °C, 1000 °C, 1100 °C, 1200 °C or 1300 °C. Following compression testing, it was apparent that the stress-strain curves were similar in shape to the curves that represent cellular solids. This is despite a relative density twice as high as what is considered the threshold for defining a cellular solid. As final sintering temperature increased, the compressive behaviour developed from being elastic-brittle to elastic-plastic and while Young's modulus remained fairly constant in the region of 1.5 GPa, there was a corresponding increase in 0.2% proof stress of approximately 40-80 MPa. The cellular solid model consists of two equations that predict Young's modulus and yield or proof stress. By fitting to experimental data and consideration of porous morphology, appropriate changes to the geometry constants allow modification of the current models to predict with better accuracy the behaviour of porous materials with higher relative densities (lower porosity).

Automated mixed dimensional modelling from 2D and 3D CAD models

The motivation for this paper is to present procedures for automatically creating idealised finite element models from the 3D CAD solid geometry of a component. The procedures produce an accurate and efficient analysis model with little effort on the part of the user. The technique is applicable to thin walled components with local complex features and automatically creates analysis models where 3D elements representing the complex regions in the component are embedded in an efficient shell mesh representing the mid-faces of the thin sheet regions. As the resulting models contain elements of more than one dimension, they are referred to as mixed dimensional models. Although these models are computationally more expensive than some of the idealisation techniques currently employed in industry, they do allow the structural behaviour of the model to be analysed more accurately, which is essential if appropriate design decisions are to be made. Also, using these procedures, analysis models can be created automatically whereas the current idealisation techniques are mostly manual, have long preparation times, and are based on engineering judgement. In the paper the idealisation approach is first applied to 2D models that are used to approximate axisymmetric components for analysis. For these models 2D elements representing the complex regions are embedded in a 1D mesh representing the midline of the cross section of the thin sheet regions. Also discussed is the coupling, which is necessary to link the elements of different dimensionality together. Analysis results from a 3D mixed dimensional model created using the techniques in this paper are compared to those from a stiffened shell model and a 3D solid model to demonstrate the improved accuracy of the new approach. At the end of the paper a quantitative analysis of the reduction in computational cost due to shell meshing thin sheet regions demonstrates that the reduction in degrees of freedom is proportional to the square of the aspect ratio of the region, and for long slender solids, the reduction can be proportional to the aspect ratio of the region if appropriate meshing algorithms are used.

Application of a mechanobiological simulation technique to stents used clinically

Many cardiovascular diseases are characterised by the restriction of blood flow through arteries. Stents can be expanded within arteries to remove such restrictions; however, tissue in-growth into the stent can lead to restenosis. In order to predict the long-term efficacy of stenting, a mechanobiological model of the arterial tissue reaction to stress is required. In this study, a computational model of arterial tissue response to stenting is applied to three clinically relevant stent designs. We ask the question whether such a mechanobiological model can differentiate between stents used clinically, and we compare these predictions to a purely mechanical analysis. In doing so, we are testing the hypothesis that a mechanobiological model of arterial tissue response to injury could predict the long-term outcomes of stent design. Finite element analysis of the expansion of three different stent types was performed in an idealised, 3D artery. Injury was calculated in the arterial tissue using a remaining-life damage mechanics approach. The inflammatory response to this initial injury was modelled using equations governing variables which represented tissue-degrading species and growth factors. Three levels of inflammation response were modelled to account for inter-patient variability. A lattice-based model of smooth muscle cell behaviour was implemented, treating cells as discrete agents governed by local rules. The simulations predicted differences between stent designs similar to those found in vivo. It showed that the volume of neointima produced could be quantified, providing a quantitative comparison of stents. In contrast, the differences between stents based on stress alone were highly dependent on the choice of comparison criteria. These results show that the choice of stress criteria for stent comparisons is critical. This study shows that mechanobiological modelling may provide a valuable tool in stent design, allowing predictions of their long-term efficacy. The level of inflammation was shown to affect the sensitivity of the model to stent design. If this finding was verified in patients, this could suggest that high-inflammation patients may require alternative treatments to stenting.

Automatic dimensional reduction and meshing of stiffened thin-wall structures

The creation of idealised, dimensionally reduced meshes for preliminary design and optimisation remains a time-consuming, manual task. A dimensionally reduced model is ideal for assessing design changes through modification of element properties without the need to create a new geometry or mesh. In this paper, a novel approach for automating the creation of mixed dimensional meshes is presented. The input to the process is a solid model which has been decomposed into a non-manifold assembly of smaller volumes with different meshing significance. Associativity between the original solid model and the dimensionally reduced equivalent is maintained. The approach is validated by means of a free-free modal analysis on an output mesh of a gas turbine engine component of industrial complexity. Extensions and enhancements to this work are also discussed.

Issues in fe modelling of FRP strengthened masonry arch bridges

Externally bonding of FRP composites is an effective technique for retrofitting historical masonry arch structures. A major failure mode in such strengthened structures is the debonding of FRP from the masonry. The bond behaviour between FRP and masonry thus plays a crucial role in these structures. Major challenges exist in the finite element modelling of such structures, such as modelling of mixed Mode-I and Mode-II bond behaviour between the FRP and the curved masonry substrate, modelling of existing damages in the masonry arches, consideration of loading history in the unstrengthened and strengthened structure etc. This paper presents a rigorous FE model for simulating FRP strengthened masonry arch structures. A detailed solid model was developed for simulating the masonry and a mixed-mode interface model was used for simulating the FRP-to-masonry bond behaviour. The model produces results in very close agreement with test results.

Investigation of various growth mechanisms of solid tumour growth within the linear-quadratic model for radiotherapy

The standard linear-quadratic survival model for radiotherapy is used to investigate different schedules of radiation treatment planning to study how these may be affected by different tumour repopulation kinetics between treatments.

Case study of the uniaxial anisotropic spin-1 bilinear-biquadratic Heisenberg model on a triangular lattice

We study the spin-1 model on a triangular lattice in the presence of a uniaxial anisotropy field using a cluster mean-field (CMF) approach. The interplay among antiferromagnetic exchange, lattice geometry, and anisotropy forces Gutzwiller mean-field approaches to fail in a certain region of the phase diagram. There, the CMF method yields two supersolid phases compatible with those present in the spin-1/2 XXZ model onto which the spin-1 system maps. Between these two supersolid phases, the three-sublattice order is broken and the results of the CMF approach depend heavily on the geometry and size of the cluster. We discuss the possible presence of a spin liquid in this region.

Solubility parameter based screening methods for early stage formulation development of itraconazole amorphous solid dispersions

Objectives: This article uses conventional and newly extended solubility parameter (δ) methods to identify polymeric materials capable of forming amorphous dispersions with itraconazole (itz). Methods: Combinations of itz and Soluplus, Eudragit E PO (EPO), Kollidon 17PF (17PF) or Kollidon VA64 (VA64) were prepared as amorphous solid dispersions using quench cooling and hot melt extrusion. Storage stability was evaluated under a range of conditions using differential scanning calorimetry and powder X-ray diffraction. Key findings: The rank order of itz miscibility with polymers using both conventional and novel δ-based approaches was 17PF > VA64 > Soluplus > EPO, and the application of the Flory–Huggins lattice model to itz–excipient binary systems corroborated the findings. The solid-state characterisation analyses of the formulations manufactured by melt extrusion correlated well with pre-formulation screening. Long-term storage studies showed that the physical stability of 17PF/vitamin E TPGS–itz was poor compared with Soluplus and VA64 formulations, and for EPO/itz systems variation in stability may be observed depending on the preparation method. Conclusion: Results have demonstrated that although δ-based screening may be useful in predicting the initial state of amorphous solid dispersions, assessment of the physical behaviour of the formulations at relevant temperatures may be more appropriate for the successful development of commercially acceptable amorphous drug products.

High speed DSC (hyper-DSC) as a tool to measure the solubility of a drug within a solid or semi-solid matrix

Conventional differential scanning calorimetry (DSC) techniques are commonly used to quantify the solubility of drugs within polymeric-controlled delivery systems. However, the nature of the DSC experiment, and in particular the relatively slow heating rates employed, limit its use to the measurement of drug solubility at the drug's melting temperature. Here, we describe the application of hyper-DSC (HDSC), a variant of DSC involving extremely rapid heating rates, to the calculation of the solubility of a model drug, metronidazole, in silicone elastomer, and demonstrate that the faster heating rates permit the solubility to be calculated under non-equilibrium conditions such that the solubility better approximates that at the temperature of use. At a heating rate of 400 degrees C/min (HDSC), metronidazole solubility was calculated to be 2.16 mg/g compared with 6.16 mg/g at 20 degrees C/min. (C) 2005 Elsevier B.V. All rights reserved.

Coupled Temperature-Displacement Modelling of Injection Stretch Blow Moulding of PET bottles using Buckley Model

This study presents a fully coupled temperature–displacement finite element modelling of the injection stretch-blow moulding (ISBM) process of polyethylene terephthalate (PET) bottles using ABAQUS with a view to optimising the process conditions. A physically-based material model (Buckley model) was used to predict the mechanical behaviour of PET at temperatures slightly above its glass transition temperature. A model incorporating heat transfer between the stretch rod, the preform and the mould was built using axisymmetric solid elements. Extensive finite element analyses were carried out to predict the deformation, the distribution and history of strain and temperature during ISBM of a 20 g–330 ml bottle, which was made in an in situ test on a Sidel SB06 machine. Comparisons of numerical results with the measurements demonstrate that the model can satisfactorily model the sidewall thickness and material distributions. It is also shown that significant non-linear differentials exist in temperature and strain in both bottle thickness and length directions during the process. This justifies the employment of a volume approach to accurately predict the final mechanical properties of the bottles governed by the orientation and crystallinity which are highly temperature and strain dependent.

A Phase-field Model for Computer Simulation of Lamellar Structure Formation in y-TiAI

In the present paper, a phase-field model is developed to simulate the formation and evolution of lamellar microstructure in γ-TiAl alloys. The mechanism of formation of TiAl lamellae proposed by Denquin and Naka is incorporated into the model. The model describes the formation and evolution of the face-centered cubic (fcc) stacking lamellar zone followed by the subsequent appearance and growth of the γ-phase, involving both the chemical composition change by atom transfer and the ordering of the fcc lattice. The thermodynamics of the model system and the interaction between the displacive and diffusional transformations are described by a non-equilibrium free energy formulated as a function of concentration and structural order parameter fields. The long-range elastic interactions, arising from the lattice misfit between the α, fcc (A1) and the various orientation variants of the γ-phase are taken into account by incorporating of the elastic strain energy into the total free energy. Simulation studies based on the model successfully predicted some essential features of the lamellar structure. It is found that the formation and evolution of the lamellar structure are predominantly controlled by the minimization of the elastic energy of the interfaces between the different fcc stacking groups, low-symmetry product phase γ and the high-symmetry α-phase, as well as between the various orientation variants of the product phase.

Towards An Acoustic Model-Based Poroelastic Imaging Method: I. Theoretical Foundation

The ultrasonic measurement and imaging of tissue elasticity is currently under wide investigation and development as a clinical tool for the assessment of a broad range of diseases, but little account in this field has yet been taken of the fact that soft tissue is porous and contains mobile fluid. The ability to squeeze fluid out of tissue may have implications for conventional elasticity imaging, and may present opportunities for new investigative tools. When a homogeneous, isotropic, fluid-saturated poroelastic material with a linearly elastic solid phase and incompressible solid and fluid constituents is subjected to stress, the behaviour of the induced internal strain field is influenced by three material constants: the Young's modulus (E(s)) and Poisson's ratio (nu(s)) of the solid matrix and the permeability (k) of the solid matrix to the pore fluid. New analytical expressions were derived and used to model the time-dependent behaviour of the strain field inside simulated homogeneous cylindrical samples of such a poroelastic material undergoing sustained unconfined compression. A model-based reconstruction technique was developed to produce images of parameters related to the poroelastic material constants (E(s), nu(s), k) from a comparison of the measured and predicted time-dependent spatially varying radial strain. Tests of the method using simulated noisy strain data showed that it is capable of producing three unique parametric images: an image of the Poisson's ratio of the solid matrix, an image of the axial strain (which was not time-dependent subsequent to the application of the compression) and an image representing the product of the aggregate modulus E(s)(1-nu(s))/(1+nu(s))(1-2nu(s)) of the solid matrix and the permeability of the solid matrix to the pore fluid. The analytical expressions were further used to numerically validate a finite element model and to clarify previous work on poroelastography.