Computer simulating a trial of a load-bearing implant: example of an intramedullary prosthesis

Computational modelling is becoming ever more important for obtaining regulatory approval for new medical devices. An accepted approach is to infer performance in a population from an analysis conducted for an idealised or ‘average’ patient; we present here a method for predicting the performance of an orthopaedic implant when released into a population—effectively simulating a clinical trial. Specifically we hypothesise that an analysis based on a method for predicting the performance in a population will lead to different conclusions than an analysis based on an idealised or ‘average’ patient. To test this hypothesis we use a finite element model of an intramedullary implant in a bone whose size and remodelling activity is different for each individual in the population. We compare the performance of a low Young’s modulus implant (View the MathML source) to one with a higher Young’s modulus (200 GPa). Cyclic loading is applied and failure is assumed when the migration of the implant relative to the bone exceeds a threshold magnitude. The analysis for an idealised of ‘average’ patient predicts that the lower modulus device survives longer whereas the analysis simulating a clinical trial predicts no statistically-significant tendency (p=0.77) for the low modulus device to perform better. It is concluded that population-based simulations of implant performance–simulating a clinical trial–present a very valuable opportunity for more realistic computational pre-clinical testing of medical devices.

Non-linear finite-element analysis of punching capacities of steel-concrete bridge deck slabs

Punching failure is the common failure mode in concrete bridge deck slabs when these structural components are subjected to local patch loads, such as tyre loads. Past research has shown that reinforced concrete slabs in girder–slab type bridges have a load-carrying capacity far greater than the ultimate static loads predicted by traditional design methods, because of the presence of compressive membrane action. However, due to the instability problems from punching failure, it is difficult to predict ultimate capacities accurately in numerical analyses. In order to overcome the instability problems, this paper establishes an efficient non-linear finite-element analysis using the commercial finite-element package Abaqus. In the non-linear finite-element analysis, stabilisation methods were adopted and failure criteria were established to predict the ultimate punching behaviour of deck slabs in composite steel–concrete bridges. The proposed non-linear finite-element analysis predictions showed a good correlation on punching capacities with experimental tests.

PLASTID DNA RESTRICTION ANALYSIS LINKS THE HETEROMORPHIC PHASES OF AN APOMICTIC RED ALGAL LIFE-HISTORY

Gymnogongrus sp. (Phyllophoraceae) from Nova Scotia, Canada, identified tentatively as G. devoniensis (Greville) Schotter, grows in association with an Erythrodermis-like crust that forms chains of tetrasporangia or bisporangia. The crust resembles tetrasporophytic phases of other Gymnogongrus species, but in culture both it and the G. ?devoniensis gametophytes cycle independently by apomictic reproduction.

Cold-formed steel sections with web openings subjected to web crippling under two-flange loading conditions - part I: Tests and finite element analysis

The results of 82 web crippling tests are presented, with 20 tests conducted on channel sections without web openings and 62 tests conducted on channel sections with web openings. The tests consider both end-two-flange and interior-two-flange loading conditions. In the case of the tests with web openings, the hole was located directly under the concentrated load. The concentrated load was applied through bearing plates; the effect of different bearing lengths is investigated. In addition, the cases of both flanges fastened and unfastened to the support is considered. A non-linear elasto-plastic finite element model is described, and the results compared against the laboratory test results; a good agreement was obtained in terms of both strength and failure modes.

Strong-field photodetachment of H- by few-cycle laser pulses

The recent adiabatic saddle-point method of Shearer et al. [ Phys. Rev. A 84 033409 (2011)] is applied to study strong-field photodetachment of H- by few-cycle linearly polarized laser pulses of frequencies near the two-photon detachment threshold. The behavior of the saddle points in the complex-time plane for a range of laser parameters is explored. A detailed analysis of the influence of laser intensities [(2×1011)–(6.5 × 1011) W/cm2], midinfrared laser wavelengths (1800–2700 nm), and various values of the carrier envelope phase (CEP) on (i) three-dimensional probability detachment distributions, (ii) photoangular distributions (PADs), (iii) energy spectra, and (iv) momentum distributions are presented. Examination of the probability distributions and PADs reveal main lobes and jetlike structures. Bifurcation phenomena in the probability distributions and PADs are also observed as the wavelength and intensity increase. Our simulations show that the (i) probability distributions, (ii) PADs, and (iii) energy spectra are extremely sensitive to the CEP and thus measuring such distributions provides a useful tool for determining this phase. The symmetrical properties of the electron momentum distributions are also found to be strongly correlated with the CEP and this provides an additional robust method for measuring the CEP of a laser pulse. Our calculations further show that for a three-cycle pulse inclusion of all eight saddle points is required in the evaluation of the transition amplitude to yield an accurate description of the photodetachment process. This is in contrast to recent results for a five-cycle pulse.

Stress correlation between instrumentation and simulation analysis of the die for high pressure die casting

A strain gauge instrumentation trial on a high pressure die casting ‘HPDC’ die was compared to a corresponding simulation model using Magmasoft® casting simulation software at two strain gauge rosette locations. The strains were measured during the casting cycle, from which the von Mises stress was determined and then compared to the simulation model. The von Mises stress from the simulation model correlated well with the findings from the instrumentation trial, showing a difference of 5.5%, ~ 10 MPa for one strain gauge rosette located in an area of low stress gradient. The second rosette was in a region of steep stress gradient, which resulted in a difference of up to 40%, ~40 MPa between the simulation and instrumentation results. Factors such as additional loading from die closure force or metal injection pressure which are not modelled by Magmasoft® were seen to have very little influence on the stress in the die, less than 7%.

Allometric analysis of Sitka spruce branches:Mechanical versus hydraulic design principles

The geometry of tree branches can have considerable effect on their efficiency in terms of carbon export per unit carbon investment in structure. The purpose of this study was to evaluate different design criteria using data describing the form of Picea sitchensis branches. Allometric analysis of the data suggests that resources are distributed to favour shoots with the greatest opportunity for extension into new space, with priority to the extension of the leader. The distribution of allometric relations of links (branch elements) was tested against two models: the pipe model, based on hydraulic transport requirements, and a static load model based on the requirement of shoots to provide mechanical resistance to static loads. Static load resistance required the load parameter to be proportional to the link radius raised to the power of 4. This was shown to be true within a 95% statistical confidence limit. The pipe model would require total distal length to be proportional to link radius squared but the measured branches did not conform well to this model. The comparison suggests that the diameters of branch elements were more related to the requirements for mechanical load. The cost of following a hydraulic design principle (the pipe model) in terms of mechanical efficiency was estimated and suggested that the pipe model branch would not be mechanically compromised but would use structural resources inefficiently. Resource allocation among branch elements was found to be consistent with mechanical stability criteria but also indicated the possibility of allocation based on other criteria, such as potential light interception by shoots. The evidence suggests that whilst branch topology increments by reiteration of units of morphogenesis, the geometry follows a functional design pattern.

A streamlined life-cycle assessment and decision tool for used tyres recycling

In order to achieve progress towards sustainable resource management, it is essential to evaluate options for the reuse and recycling of secondary raw materials, in order to provide a robust evidence base for decision makers. This paper presents the research undertaken in the development of a web-based decision-support tool (the used tyres resource efficiency tool) to compare three processing routes for used tyres compared to their existing primary alternatives. Primary data on the energy and material flows for the three routes, and their alternatives were collected and analysed. The methodology used was a streamlined life-cycle assessment (sLCA) approach. Processes included were: car tyre baling against aggregate gabions; car tyre retreading against new car tyres; and car tyre shred used in landfill engineering against primary aggregates. The outputs of the assessment, and web-based tool, were estimates of raw materials used, carbon dioxide emissions and costs. The paper discusses the benefits of carrying out a streamlined LCA and using the outputs of this analysis to develop a decision-support tool. The strengths and weakness of this approach are discussed and future research priorities identified which could facilitate the use of life cycle approaches by designers and practitioners.

Genetic variants influencing human aging from late-onset Alzheimer's disease (LOAD) genome-wide association studies (GWAS)

Genetics plays a crucial role in human aging with up to 30% of those living to the mid-80s being determined by genetic variation. Survival to older ages likely entails an even greater genetic contribution. There is increasing evidence that genes implicated in age-related diseases, such as cancer and neuronal disease, play a role in affecting human life span. We have selected the 10 most promising late-onset Alzheimer's disease (LOAD) susceptibility genes identified through several recent large genome-wide association studies (GWAS). These 10 LOAD genes (APOE, CLU, PICALM, CR1, BIN1, ABCA7, MS4A6A, CD33, CD2AP, and EPHA1) have been tested for association with human aging in our dataset (1385 samples with documented age at death [AAD], age range: 58-108 years; mean age at death: 80.2) using the most significant single nucleotide polymorphisms (SNPs) found in the previous studies. Apart from the APOE locus (rs2075650) which showed compelling evidence of association with risk on human life span (p = 5.27 × 10(-4)), none of the other LOAD gene loci demonstrated significant evidence of association. In addition to examining the known LOAD genes, we carried out analyses using age at death as a quantitative trait. No genome-wide significant SNPs were discovered. Increasing sample size and statistical power will be imperative to detect genuine aging-associated variants in the future. In this report, we also discuss issues relating to the analysis of genome-wide association studies data from different centers and the bioinformatic approach required to distinguish spurious genome-wide significant signals from real SNP associations.

A whole-genome massively parallel sequencing analysis of BRCA1 mutant oestrogen receptor-negative and -positive breast cancers

BRCA1 encodes a tumour suppressor protein that plays pivotal roles in homologous recombination (HR) DNA repair, cell-cycle checkpoints, and transcriptional regulation. BRCA1 germline mutations confer a high risk of early-onset breast and ovarian cancer. In more than 80% of cases, tumours arising in BRCA1 germline mutation carriers are oestrogen receptor (ER)-negative; however, up to 15% are ER-positive. It has been suggested that BRCA1 ER-positive breast cancers constitute sporadic cancers arising in the context of a BRCA1 germline mutation rather than being causally related to BRCA1 loss-of-function. Whole-genome massively parallel sequencing of ER-positive and ER-negative BRCA1 breast cancers, and their respective germline DNAs, was used to characterize the genetic landscape of BRCA1 cancers at base-pair resolution. Only BRCA1 germline mutations, somatic loss of the wild-type allele, and TP53 somatic mutations were recurrently found in the index cases. BRCA1 breast cancers displayed a mutational signature consistent with that caused by lack of HR DNA repair in both ER-positive and ER-negative cases. Sequencing analysis of independent cohorts of hereditary BRCA1 and sporadic non-BRCA1 breast cancers for the presence of recurrent pathogenic mutations and/or homozygous deletions found in the index cases revealed that DAPK3, TMEM135, KIAA1797, PDE4D, and GATA4 are potential additional drivers of breast cancers. This study demonstrates that BRCA1 pathogenic germline mutations coupled with somatic loss of the wild-type allele are not sufficient for hereditary breast cancers to display an ER-negative phenotype, and has led to the identification of three potential novel breast cancer genes (ie DAPK3, TMEM135, and GATA4).

Cloning and analysis of a Trichinella pseudospiralis muscle larva secreted serine protease gene

Nematode parasites of the genus Trichinella are intracellular and distinct life cycle stages invade intestinal epithelial and skeletal muscle cells. Within the genus, Trichinella spiralis and Trichinella pseudospiralis exhibit species-specific differences with respect to host-parasite complex formation and host immune modulation. Parasite excretory-secretory (ES) proteins play important roles at the host-parasite interface and are thought to underpin these differences in biology. Serine proteases are among the most abundant group of T. spiralis ES proteins and multiple isoforms of the muscle larvae-specific TspSP-1 serine protease have been identified. Recently, a similar protein (TppSP-1) in T. pseudospiralis muscle larvae was identified. Here we report the cloning and characterisation of the full-length transcript of TppSP-1 and present comparative data between TspSP-1 and TppSP-1.

Numerical analysis of stiffener runout sections

The recent trend of incorporating more composite material in primary aircraft structures has highlighted the vulnerability of stiffened aerostructures to through-thickness stresses, which may lead to delamination and debonding at the skin-stiffener interface, leading to collapse. Stiffener runout regions are particularly susceptible to this problem and cannot be avoided due to the necessity to terminate stiffeners at rib intersections or at cutouts, interrupting the stiffener load path. In this paper, experimental tests relating to two different stiffener runout specimens are presented and the failure modes of both specimens are discussed in detail. A thinner-skinned specimen showed sudden and unstable crack propagation, while a thicker-skinned specimen showed initially unstable but subsequent stable crack growth. Detailed finite element models of the two specimens are developed, and it is shown how such models can explain and predict the behaviour and failure mode of stiffener runouts. The models contain continuum shell elements to model the skin and stiffener, while cohesive elements using a traction-separation law are placed at the skin-stiffener interface to effectively model the debonding which promotes structural failure.

The Behavior of compressively loaded stiffener runout specimens - Part II: Finite element analysis

The termination of stiffeners in composite aircraft structures give rise to regions of high interlaminar shear and peel stresses as the load in the stiffener is diffused into the skin. This is of particular concern in co-cured composite stiffened structures where there is a relatively low resistance to through-thickness stress components at the skin-stiffener interface. In Part I, experimental results of tested specimens highlighted the influence of local design parameters on their structural response. Indeed some of the observed behavior was unexpected. There is a need to be able to analyse a range of changes in geometry rapidly to allow the analysis to form an integral part of the structural design process.

This work presents the development of a finite element methodology for modelling the failure process of these critical regions. An efficient thick shell element formulation is presented and this element is used in conjuction with the Virtual Crack Closure Technique (VCCT) to predict the crack growth characteristics of the modelled specimens. Three specimens were modelled and the qualitative aspects of crack growth were captured successfully. The shortcomings in the quantitative correlation between the predicted and observed failure loads are discussed. There was evidence to suggest that high through-thickness compressive stresses enhanced the fracture toughness in these critical regions.

Design and analysis of dual-rail circuits for security applications

Dual-rail encoding, return-to-spacer protocol, and hazard-free logic can be used to resist power analysis attacks by making energy consumed per clock cycle independent of processed data. Standard dual-rail logic uses a protocol with a single spacer, e.g., all-zeros, which gives rise to energy balancing problems. We address these problems by incorporating two spacers; the spacers alternate between adjacent clock cycles. This guarantees that all gates switch in every clock cycle regardless of the transmitted data values. To generate these dual-rail circuits, an automated tool has been developed. It is capable of converting synchronous netlists into dual-rail circuits and it is interfaced to industry CAD tools. Dual-rail and single-rail benchmarks based upon the advanced encryption standard (AES) have been simulated and compared in order to evaluate the method and the tool.

Design technique for mm-wave IC realization of the load network of switched-mode class-EF power amplifier

Analysis and synthesis of the new Class-EF power amplifier (PA) are presented in this paper. The proposed circuit offers means to alleviate some of the major issues faced by existing Class-EF and Class-EF PAs, such as (1) substantial power losses due to parasitic resistance of the large inductor in the Class-EF load network, (2) unpredictable behaviour of practical lumped inductors and capacitors at harmonic frequencies, and (3) deviation from ideal Class-EF operation mode due to detrimental effects of device output inductance at high frequencies. The transmission-line load network of the Class-EF PA topology elaborated in this paper simultaneously satisfies the Class-EF optimum impedance requirements at fundamental frequency, second, and third harmonics as well as simultaneously providing matching to the circuit optimum load resistance for any prescribed system load resistance. Furthermore, an elegant solution using an open and short-circuit stub arrangement is suggested to overcome the problem encountered in the mm-wave IC realizations of the Class-EF PA load network due to lossy quarter-wave line. © 2010 IEICE Institute of Electronics Informati.