Discussion of : constructing summary statistics for approximate Bayesian computation: semi-automatic approximate Bayesian computation

This is a discussion of the journal article: "Construcing summary statistics for approximate Bayesian computation: semi-automatic approximate Bayesian computation". The article and discussion have appeared in the Journal of the Royal Statistical Society: Series B (Statistical Methodology).

Sub-modelling for the ratchetting failure of insulated rail joints

Insulated rail joints are critical for train safety as they control electrical signalling systems; unfortunately they exhibit excessive ratchetting of the railhead near the endpost insulators. This paper reports a three-dimensional global model of these joints under wheel–rail contact pressure loading and a sub-model examining the ratchetting failures of the railhead. The sub-model employs a non-linear isotropic–kinematic elastic–plastic material model and predicts stress/strain levels in the localised railhead zone adjacent to the endpost which is placed in the air gap between the two rail ends at the insulated rail joint. The equivalent plastic strain plot is utilised to capture the progressive railhead damage adequately. Associated field and laboratory testing results of damage to the railhead material suggest that the simulation results are reasonable.

Rolling contact fatigue in rail - insulated rail joints (IRJ)

The insulated rail joint (IRJ) is an essential component in a track circuit that controls the signaling system. Failure of IRJs leads to improper functioning of the signals,with potential for catastrophic results. Therefore, IRJs are regarded as safety-critical sections of rail network; hence, all of their components must be maintained in pristine design condition.

Structural and material characterisation of insulated rail joints

Insulated rail joints are designed in a similar way to butt jointed steel structural systems, the difference being a purpose made gap between the main rail members to maintain electrical insulation for the proper functioning of the track circuitry at all times of train operation. When loaded wheels pass the gap, they induce an impact loading with the corresponding strains in the railhead edges exceeding the plastic limit significantly, which lead to metal flow across the gap thereby increasing the risk of short circuiting and impeding the proper functioning of the signalling and broken rail identification circuitries, of which the joints are a critical part. The performance of insulated rail joints under the passage of the wheel loading is complex due to the presence of a number of interacting components and hence is not well understood. This paper presents a dynamic wheel-rail contact-impact modelling method for the determination of the impact loading; a brief description of a field experiment to capture strain signatures for validating the predicted impact loading is also presented. The process and the results of the characterisation of the materials from virgin, in-service and damaged insulated rail joints using neutron diffraction method are also discussed.

Residual stresses in rail-ends from the in-service insulated rail joints using neutron diffraction

Insulated rail joints (IRJs) are an integral part of the rail track signaling system and pose significant maintenance and replacement costs due to their low and fluctuating service lives. Failure occurs mainly in rail head region, bolt- holes of fishplates and web-holes of the rails. Propagation of cracks is influenced by the evolution of internal residual stresses in rails during rail manufacturing (hot-rolling, roller-straightening, and head-hardening process), and during service, particularly in heavy rail haul freight systems where loads are high. In this investigation, rail head accumulated residual stresses were analysed using neutron diffraction at the Australian Nuclear Science and Technology Organisation (ANSTO). Two ex-service two head-hardened rail joints damaged under different loading were examined and results were compared with those obtained from an unused rail joint reference sample in order to differentiate the stresses developed during rail manufacturing and stresses accumulated during rail service. Neutron diffraction analyses were carried out on the samples in longitudinal, transverse and vertical directions, and on 5mm thick sliceed samples cut by Electric Discharge Machining (EDM). For the rail joints from the service line, irrespective of loading conditions and in-service times, results revealed similar depth profiles of stress distribution. Evolution of residual stress fields in rails due to service was also accompanied by evidence of larger material flow based on reflected light and scanning electron microscopy studies. Stress evolution in the vicinity of rail ends was characterised by a compressive layer, approximately 5 mm deep, and a tension zone located approximately 5- 15mm below the surfaces. A significant variation of d0 with depth near the top surface was detected and was attributed to decarburization in the top layer induced by cold work. Stress distributions observed in longitudinal slices of the two different deformed rail samples were found to be similar. For the undeformed rail, the stress distributions obtained could be attributed to variations associated with thermo-mechanical history of the rail.

Nonlinear analysis of plastic ball grid array solder joints

A nonlinear finite element analysis was carried out to investigate the viscoplastic deformation of solder joints in a ball grid array (BGA) package under temperature cycle. The effects of constraint on print circuit board (PCB) and stiffness of substrate on the deformation behaviour of the solder joints were also studied. A relative damage stress was adopted to analyze the potential failure sites in the solder joints. The results indicated that high inelastic strain and strain energy density were developed in the joints close to the package center. On the other hand, high constraint and high relative damage stress were associated with the joint closest to the edge of the silicon chip. The joint closest to the edge of the silicon chip was regarded as the most susceptible failure site if cavitation instability is the dominant failure mechanism. Increase the external constraint on the print circuit board (PCB) causes a slight increase in stress triaxiality (m/eq) and relative damage stress in the joint closest to the edge of silicon die. The relative damage stress is not sensitive to the Young’s modulus of the substrate.

Second-order elastic finite element analysis of structures using a single element per member

Finite element frame analysis programs targeted for design office application necessitate algorithms which can deliver reliable numerical convergence in a practical timeframe with comparable degrees of accuracy, and a highly desirable attribute is the use of a single element per member to reduce computational storage, as well as data preparation and the interpretation of the results. To this end, a higher-order finite element method including geometric non-linearity is addressed in the paper for the analysis of elastic frames for which a single element is used to model each member. The geometric non-linearity in the structure is handled using an updated Lagrangian formulation, which takes the effects of the large translations and rotations that occur at the joints into consideration by accumulating their nodal coordinates. Rigid body movements are eliminated from the local member load-displacement relationship for which the total secant stiffness is formulated for evaluating the large member deformations of an element. The influences of the axial force on the member stiffness and the changes in the member chord length are taken into account using a modified bowing function which is formulated in the total secant stiffness relationship, for which the coupling of the axial strain and flexural bowing is included.

Visible and near infrared spectroscopy of hayabusa reentry using semi-autonomous tracking

Aground-based tracking camera and coaligned slitless spectrograph were used to measure the spectral signature of visible radiation emitted from the Hayabusa capsule as it entered into the Earth’s atmosphere in June 2010. Good quality spectra were obtained, which showed the presence of radiation from the heat shield of the vehicle and the shock-heated air in front of the vehicle. An analysis of the blackbody nature of the radiation concluded that the peak average temperature of the surface was about (3100± 100)K. Line spectra from oxygen and nitrogen atoms were used to infer a peak average shock-heated gas temperature of around((7000±400))K.

A biomechanical analysis of growing rods used in the management of early onset scoliosis

INTRODUCTION Managing spinal deformities in young children is challenging, particularly early onset scoliosis (EOS). Surgical intervention is often required if EOS has been unresponsive to conservative treatment particularly with rapidly progressive curves. An emerging treatment option for EOS is fusionless scoliosis surgery. Similar to bracing, this surgical option potentially harnesses growth, motion and function of the spine along with correcting spinal deformity. Dual growing rods are one such fusionless treatment, which aims to modulate growth of the vertebrae. The aim of this study was to ascertain the extent to which semiconstrained growing rods (Medtronic, Sofamor, Danek, Memphis, TN) with a telescopic sleeve component, reduce rotational constraint on the spine compared with standard "constrained / rigid" rods and hence potentially provide a more physiological mechanical environment for the growing spine. METHODS Six 40-60kg English Large White porcine spines served as a model for the paediatric human spine. Each spine was dissected into a 7 level thoracolumbar multi-segment unit (MSU), removing all non-ligamentous soft tissues and leaving 3cm of ribs either side. Pure nondestructive axial rotation moments of ±4Nm at a constant rotation rate of 8deg.s-1 were applied to the mounted MSU spines using a biaxial Instron testing machine. Displacement of each vertebral level was captured using a 3D motion tracking system (Optotrak 3020, Northern Digital Inc, Waterloo, ON). Each spine was tested in an un-instrumented state first and then with appropriately sized semi-constrained growing rods and rigid rods in alternating sequence. The rods were secured by multi-axial pedicle screws (Medtronic CD Horizon) at levels 2 and 6 of the construct. The range of motion (ROM), neutral zone (NZ) size and stiffness (Nm.deg-1) were calculated from the Instron load-displacement data and intervertebral ROM was calculated through a MATLAB algorithm from Optotrak data. RESULTS Irrespective of the order of testing, rigid rods significantly reduced the total ROM compared with semi-constrained rods (p<0.05) with in a significantly stiffer spine for both left and right axial rotation (p<0.05). Analysing the intervertebral motion within the instrumented levels 2-6, rigid rods showed reduced ROM compared with semi-constrained growing rods and compared with un-instrumented motion segments. CONCLUSION Semi-constrained growing rods maintain similar stiffness in axial rotation to un-instrumented spines, while dual rigid rods significantly reduce axial rotation. Clinically the effect of semi-constrained growing rods as observed in this study is that they would be expected to allow growth via the telescopic rod components while maintaining the axial flexibility of the spine, which may reduce occurrence of the crankshaft phenomenon.

The stabilising role of the costovertebral joints and spinous processes on thoracic spine stability

Introduction The importance of in vitro biomechanical testing in today’s understanding of spinal pathology and treatment modalities cannot be stressed enough. Different studies have used differing levels of dissection of their spinal segments for their testing protocols[1, 2]. The aim of this study was to assess the impact of removing the costovertebral joints and partial resection of the spinous process sequentially, on the stiffness of the immature thoracic bovine spinal segment. Materials and Methods Thoracic spines from 6-8 week old calves were used. Each spine was dissected and divided into motion segments with 5cm of attached rib on each side and full spinous processes including levels T4-T11 (n=28). They were potted in polymethylemethacrylate. An Instron Biaxial materials testing machine with a custom made jig was used for testing. The segments were tested in flexion/extension, lateral bending and axial rotation at 37⁰C and 100% humidity, using moment control to a maximum 1.75 Nm with a loading rate of 0.3 Nm per second. They were first tested intact for ten load cycles with data collected from the tenth cycle. Progressive dissection was performed by removing first the attached ribs, followed by the spinous process at its base. Biomechanical testing was carried out after each level of dissection using the same protocol. Statistical analysis of the data was performed using repeated measures ANOVA. Results In combined flexion/extension there was a significant reduction in stiffness of 16% (p=0.002). This was mainly after resection of the ribs (14%, p=0.024) and mainly occurred in flexion where stiffness reduced by 22% (p=0.021). In extension, stiffness dropped by 13% (p=0.133). However there was no further significant change in stiffness on resection of the spinous process (<1%) (p=1.00). In lateral bending there was a significant decrease in stiffness of 13% (p<0.001). This comprised a drop of 11% on resection of the ribs (p=0.009) and a further 8% on resection of the spinous process (p=0.014). There was no difference between left and right bending. In axial rotation there was no significant change in stiffness after each stage of dissection (p=0.253). There was no difference between left and right rotation. Conclusion The costovertebral joints play a significant role in providing stability to the bovine thoracic spine in both flexion/extension and lateral bending, whereas the spinous processes play a minor role. Both elements have little effect on axial rotation stability.

Demonstration of nonlinear absorption in Au semi-continuous film by electrical measurement

We have demonstrated the nonlinear absorption at 532 nm wavelength in an Au semi-continuous film (SF) resulting from smearing of the Fermi distribution and diffusion of conduction electrons into the substrate. The Au SF was irradiated by a pulsed laser with 8 ns pulse width at 532 nm in near resonance with the interband transition of the Au. We determined the temperature increase in the SF for different intensities by electrical measurement. We calculated the temperature increase by using a 1D heat transport equation; comparing the results of the calculation with measured values for the temperature increase, revealed the nonlinear absorption in the Au SF. We employed this deviation from linear behaviour to determine the nonlinear absorption coefficient.

Effects of plasterboard joints on the fire performance of cold-formed steel walls

This paper presents the effect of plasterboard joints on the fire performance of cold-formed steel walls. Plasterboard joints are unavoidable. However, they can be arranged in a way that they do not significantly influence the fire performance of cold-formed steel walls. Hence a research study into the effects of plasterboard joints on the fire performance of plasterboard lined cold-formed steel walls was undertaken using both full-scale fire tests and numerical studies. In this study a back-blocking technique was used to eliminate the plasterboard joints being located over the studs. Instead plasterboard joints were used between studs with 150 mm wide plasterboards as back-blocks. Both experimental and numerical results from this study show that the fire resistance rating of single plasterboard lined cold-formed steel walls can be increased by 25% through the use of a back-blocking joint arrangement in comparison to the traditional plasterboard joint arrangement over the studs.

A biomechanical analysis of growing rods used in the management of early onset scoliosis

INTRODUCTION Managing spinal deformities in young children is challenging, particularly early-onset scoliosis (EOS). Any progressive spinal deformity particularly in early life presents significant health risks for the child and a challenge for the treating surgeon. Surgical intervention is often required if EOS has been unresponsive to conservative treatment particularly with rapidly progressive curves. An emerging treatment option particularly for EOS is fusionless scoliosis surgery. Similar to bracing this surgical option potentially harnesses growth, motion and function of the spine along with correcting spinal deformity. Dual growing rods is one such fusionless treatment, which aims to modulate growth of the vertebrae. The aim of this study was to ascertain the extent to which semi-constrained growing rods (Medtronic, Memphis, TN) with a telescopic sleeve component, reduce rotational constraint on the spine compared with standard rigid rods and hence potentially provide a more physiological mechanical environment for the growing spine. METHODS Six 40-60kg English Large White porcine spines served as a model for the paediatric human spine. Each spine was dissected into 7 level thoracolumbar multi-segment unit (MSU) spines, removing all non-ligamentous soft tissues. Appropriately sized semi-constrained growing rods and rigid rods were secured by multi-axial screws (Medtronic) prior to testing in alternating sequences for each spine. Pure nondestructive moments of +/4Nm at a constant rotation rate of 8deg/s was applied to the mounted MSU spines. Displacement of each level was captured using an Optotrak (Northern Digital Inc, Waterloo, ON). The range of motion (ROM), neutral zone (NZ) size and stiffness (Nm/deg) were calculated from the Instron load-displacement data and intervertebral ROM was calculated through a MATLAB algorithm from Optotrak data. RESULTS Irrespective of sequence order rigid rods significantly reduced the total ROM (deg) than compared to semi-constrained rods (p<0.05) and resulted in a significantly stiffer (Nm/deg) spine for both left and right axial rotation testing (p<0.05). Analysing the intervertebral motion within the instrumented levels, rigid rods showed reduced ROM (Deg) than compared to semi-constrained growing rods and the un-instrumented (UN-IN) test sequences. CONCLUSION The semi-constrained growing rods maintained rotation similar to UN-IN spines while the rigid rods showed significantly reduced axial rotation across all instrumented levels. Clinically the effect of semi-constrained growing rods evaluated in this study is that they will allow growth via the telescopic rod components while maintaining the axial rotation ability of the spine, which may also reduce the occurrence of the crankshaft phenomenon.

Reversible and semi-blind relational database watermarking

Database watermarking has received significant research attention in the current decade. Although, almost all watermarking models have been either irreversible (the original relation cannot be restored from the watermarked relation) and/or non-blind (requiring original relation to detect the watermark in watermarked relation). This model has several disadvantages over reversible and blind watermarking (requiring only watermarked relation and secret key from which the watermark is detected and original relation is restored) including inability to identify rightful owner in case of successful secondary watermarking, inability to revert the relation to original data set (required in high precision industries) and requirement to store unmarked relation at a secure secondary storage. To overcome these problems, we propose a watermarking scheme that is reversible as well as blind. We utilize difference expansion on integers to achieve reversibility. The major advantages provided by our scheme are reversibility to high quality original data set, rightful owner identification, resistance against secondary watermarking attacks, and no need to store original database at a secure secondary storage.