Early changes in biochemical markers of bone formation during teriparatide therapy correlate with improvements in vertebral strength in men with glucocorticoid-induced osteoporosis

Summary Changes of the bone formation marker PINP correlated positively with improvements in vertebral strength in men with glucocorticoid-induced osteoporosis (GIO) who received 18-month treatment with teriparatide, but not with risedronate. These results support the use of PINP as a surrogate marker of bone strength in GIO patients treated with teriparatide. Introduction To investigate the correlations between biochemical markers of bone turnover and vertebral strength estimated by finite element analysis (FEA) in men with GIO. Methods A total of 92 men with GIO were included in an 18-month, randomized, open-label trial of teriparatide (20 μg/day, n = 45) and risedronate (35 mg/week, n = 47). High-resolution quantitative computed tomography images of the 12th thoracic vertebra obtained at baseline, 6 and 18 months were converted into digital nonlinear FE models and subjected to anterior bending, axial compression and torsion. Stiffness and strength were computed for each model and loading mode. Serum biochemical markers of bone formation (amino-terminal-propeptide of type I collagen [PINP]) and bone resorption (type I collagen cross-linked C-telopeptide degradation fragments [CTx]) were measured at baseline, 3 months, 6 months and 18 months. A mixed-model of repeated measures analysed changes from baseline and between-group differences. Spearman correlations assessed the relationship between changes from baseline of bone markers with FEA variables. Results PINP and CTx levels increased in the teriparatide group and decreased in the risedronate group. FEA-derived parameters increased in both groups, but were significantly higher at 18 months in the teriparatide group. Significant positive correlations were found between changes from baseline of PINP at 3, 6 and 18 months with changes in FE strength in the teriparatide-treated group, but not in the risedronate group. Conclusions Positive correlations between changes in a biochemical marker of bone formation and improvement of biomechanical properties support the use of PINP as a surrogate marker of bone strength in teriparatide-treated GIO patients.

High resolution quantitative computed tomography-based assessment of trabecular microstructure and strength estimates by finite-element analysis of the spine, but not DXA, reflects vertebral fracture status in men with glucocorticoid-induced osteoporosis

High-resolution quantitative computed tomography (HRQCT)-based analysis of spinal bone density and microstructure, finite element analysis (FEA), and DXA were used to investigate the vertebral bone status of men with glucocorticoid-induced osteoporosis (GIO). DXA of L1–L3 and total hip, QCT of L1–L3, and HRQCT of T12 were available for 73 men (54.6±14.0years) with GIO. Prevalent vertebral fracture status was evaluated on radiographs using a semi-quantitative (SQ) score (normal=0 to severe fracture=3), and the spinal deformity index (SDI) score (sum of SQ scores of T4 to L4 vertebrae). Thirty-one (42.4%) subjects had prevalent vertebral fractures. Cortical BMD (Ct.BMD) and thickness (Ct.Th), trabecular BMD (Tb.BMD), apparent trabecular bone volume fraction (app.BV/TV), and apparent trabecular separation (app.Tb.Sp) were analyzed by HRQCT. Stiffness and strength of T12 were computed by HRQCT-based nonlinear FEA for axial compression, anterior bending and axial torsion. In logistic regressions adjusted for age, glucocorticoid dose and osteoporosis treatment, Tb.BMD was most closely associated with vertebral fracture status (standardized odds ratio [sOR]: Tb.BMD T12: 4.05 [95% CI: 1.8–9.0], Tb.BMD L1–L3: 3.95 [1.8–8.9]). Strength divided by cross-sectional area for axial compression showed the most significant association with spine fracture status among FEA variables (2.56 [1.29–5.07]). SDI was best predicted by a microstructural model using Ct.Th and app.Tb.Sp (r2=0.57, p<0.001). Spinal or hip DXA measurements did not show significant associations with fracture status or severity. In this cross-sectional study of males with GIO, QCT, HRQCT-based measurements and FEA variables were superior to DXA in discriminating between patients of differing prevalent vertebral fracture status. A microstructural model combining aspects of cortical and trabecular bone reflected fracture severity most accurately.

Finite element analyses of human vertebral bodies embedded in polymethylmethalcrylate or loaded via the hyperelastic intervertebral disc models provide equivalent predictions of experimental strength

Quantitative computer tomography (QCT)-based finite element (FE) models of vertebral body provide better prediction of vertebral strength than dual energy X-ray absorptiometry. However, most models were validated against compression of vertebral bodies with endplates embedded in polymethylmethalcrylate (PMMA). Yet, loading being as important as bone density, the absence of intervertebral disc (IVD) affects the strength. Accordingly, the aim was to assess the strength predictions of the classic FE models (vertebral body embedded) against the in vitro and in silico strengths of vertebral bodies loaded via IVDs. High resolution peripheral QCT (HR-pQCT) were performed on 13 segments (T11/T12/L1). T11 and L1 were augmented with PMMA and the samples were tested under a 4° wedge compression until failure of T12. Specimen-specific model was generated for each T12 from the HR-pQCT data. Two FE sets were created: FE-PMMA refers to the classical vertebral body embedded model under axial compression; FE-IVD to their loading via hyperelastic IVD model under the wedge compression as conducted experimentally. Results showed that FE-PMMA models overestimated the experimental strength and their strength prediction was satisfactory considering the different experimental set-up. On the other hand, the FE-IVD models did not prove significantly better (Exp/FE-PMMA: R²=0.68; Exp/FE-IVD: R²=0.71, p=0.84). In conclusion, FE-PMMA correlates well with in vitro strength of human vertebral bodies loaded via real IVDs and FE-IVD with hyperelastic IVDs do not significantly improve this correlation. Therefore, it seems not worth adding the IVDs to vertebral body models until fully validated patient-specific IVD models become available.

Comparative effects of Teriparatide and Risedronate in glucocorticoid‐induced osteoporosis in men: 18‐month results of the EuroGIOPs trial

Data on treatment of glucocorticoid-induced osteoporosis (GIO) in men are scarce. We performed a randomized, open-label trial in men who have taken glucocorticoids (GC) for ≥3 months, and had an areal bone mineral density (aBMD) T-score ≤ –1.5 standard deviations. Subjects received 20 μg/d teriparatide (n = 45) or 35 mg/week risedronate (n = 47) for 18 months. Primary objective was to compare lumbar spine (L1–L3) BMD measured by quantitative computed tomography (QCT). Secondary outcomes included BMD and microstructure measured by high-resolution QCT (HRQCT) at the 12th thoracic vertebra, biomechanical effects for axial compression, anterior bending, and axial torsion evaluated by finite element (FE) analysis from HRQCT data, aBMD by dual X-ray absorptiometry, biochemical markers, and safety. Computed tomography scans were performed at 0, 6, and 18 months. A mixed model repeated measures analysis was performed to compare changes from baseline between groups. Mean age was 56.3 years. Median GC dose and duration were 8.8 mg/d and 6.4 years, respectively; 39.1% of subjects had a prevalent fracture, and 32.6% received prior bisphosphonate treatment. At 18 months, trabecular BMD had significantly increased for both treatments, with significantly greater increases with teriparatide (16.3% versus 3.8%; p = 0.004). HRQCT trabecular and cortical variables significantly increased for both treatments with significantly larger improvements for teriparatide for integral and trabecular BMD and bone surface to volume ratio (BS/BV) as a microstructural measure. Vertebral strength increases at 18 months were significant in both groups (teriparatide: 26.0% to 34.0%; risedronate: 4.2% to 6.7%), with significantly higher increases in the teriparatide group for all loading modes (0.005 < p < 0.015). Adverse events were similar between groups. None of the patients on teriparatide but five (10.6%) on risedronate developed new clinical fractures (p = 0.056). In conclusion, in this 18-month trial in men with GIO, teriparatide showed larger improvements in spinal BMD, microstructure, and FE-derived strength than risedronate.

Experimental Validation of Finite Element Analysis of Human Vertebral Collapse Under Large Compressive Strains

Osteoporosis-related vertebral fractures represent a major health problem in elderly populations. Such fractures can often only be diagnosed after a substantial deformation history of the vertebral body. Therefore, it remains a challenge for clinicians to distinguish between stable and progressive potentially harmful fractures. Accordingly, novel criteria for selection of the appropriate conservative or surgical treatment are urgently needed. Computer tomography-based finite element analysis is an increasingly accepted method to predict the quasi-static vertebral strength and to follow up this small strain property longitudinally in time. A recent development in constitutive modeling allows us to simulate strain localization and densification in trabecular bone under large compressive strains without mesh dependence. The aim of this work was to validate this recently developed constitutive model of trabecular bone for the prediction of strain localization and densification in the human vertebral body subjected to large compressive deformation. A custom-made stepwise loading device mounted in a high resolution peripheral computer tomography system was used to describe the progressive collapse of 13 human vertebrae under axial compression. Continuum finite element analyses of the 13 compression tests were realized and the zones of high volumetric strain were compared with the experiments. A fair qualitative correspondence of the strain localization zone between the experiment and finite element analysis was achieved in 9 out of 13 tests and significant correlations of the volumetric strains were obtained throughout the range of applied axial compression. Interestingly, the stepwise propagating localization zones in trabecular bone converged to the buckling locations in the cortical shell. While the adopted continuum finite element approach still suffers from several limitations, these encouraging preliminary results towardsthe prediction of extended vertebral collapse may help in assessing fracture stability in future work.

Validation of stress magnetic resonance imaging of the canine stifle joint with and without an intact cranial cruciate ligament.

OBJECTIVE To validate use of stress MRI for evaluation of stifle joints of dogs with an intact or deficient cranial cruciate ligament (CrCL). SAMPLE 10 cadaveric stifle joints from 10 dogs. PROCEDURES A custom-made limb-holding device and a pulley system linked to a paw plate were used to apply axial compression across the stifle joint and induce cranial tibial translation with the joint in various degrees of flexion. By use of sagittal proton density-weighted MRI, CrCL-intact and deficient stifle joints were evaluated under conditions of loading stress simulating the tibial compression test or the cranial drawer test. Medial and lateral femorotibial subluxation following CrCL transection measured under a simulated tibial compression test and a cranial drawer test were compared. RESULTS By use of tibial compression test MRI, the mean ± SD cranial tibial translations in the medial and lateral compartments were 9.6 ± 3.7 mm and 10 ± 4.1 mm, respectively. By use of cranial drawer test MRI, the mean ± SD cranial tibial translations in the medial and lateral compartments were 8.3 ± 3.3 mm and 9.5 ± 3.5 mm, respectively. No significant difference in femorotibial subluxation was found between stress MRI techniques. Femorotibial subluxation elicited by use of the cranial drawer test was greater in the lateral than in the medial compartment. CONCLUSIONS AND CLINICAL RELEVANCE Both stress techniques induced stifle joint subluxation following CrCL transection that was measurable by use of MRI, suggesting that both methods may be further evaluated for clinical use.

End caps prevent nail migration in elastic stable intramedullary nailing in paediatric femoral fractures: a biomechanical study using synthetic and cadaveric bones.

End caps are intended to prevent nail migration (push-out) in elastic stable intramedullary nailing. The aim of this study was to investigate the force at failure with and without end caps, and whether different insertion angles of nails and end caps would alter that force at failure. Simulated oblique fractures of the diaphysis were created in 15 artificial paediatric femurs. Titanium Elastic Nails with end caps were inserted at angles of 45°, 55° and 65° in five specimens for each angle to create three study groups. Biomechanical testing was performed with axial compression until failure. An identical fracture was created in four small adult cadaveric femurs harvested from two donors (both female, aged 81 and 85 years, height 149 cm and 156 cm, respectively). All femurs were tested without and subsequently with end caps inserted at 45°. In the artificial femurs, maximum force was not significantly different between the three groups (p = 0.613). Push-out force was significantly higher in the cadaveric specimens with the use of end caps by an up to sixfold load increase (830 N, standard deviation (SD) 280 vs 150 N, SD 120, respectively; p = 0.007). These results indicate that the nail and end cap insertion angle can be varied within 20° without altering construct stability and that the risk of elastic stable intramedullary nailing push-out can be effectively reduced by the use of end caps.

Comparison of proximal femur and vertebral body strength improvements in the FREEDOM trial using an alternative finite element methodology

Denosumab reduced the incidence of new fractures in postmenopausal women with osteoporosis by 68% at the spine and 40% at the hip over 36 months compared with placebo in the FREEDOM study. This efficacy was supported by improvements from baseline in vertebral (18.2%) strength in axial compression and femoral (8.6%) strength in sideways fall configuration at 36 months, estimated in Newtons by an established voxel-based finite element (FE) methodology. Since FE analyses rely on the choice of meshes, material properties, and boundary conditions, the aim of this study was to independently confirm and compare the effects of denosumab on vertebral and femoral strength during the FREEDOM trial using an alternative smooth FE methodology. Unlike the previous FE study, effects on femoral strength in physiological stance configuration were also examined. QCT data for the proximal femur and two lumbar vertebrae were analyzed by smooth FE methodology at baseline, 12, 24, and 36 months for 51 treated (denosumab) and 47 control (placebo) subjects. QCT images were segmented and converted into smooth FE models to compute bone strength. L1 and L2 vertebral bodies were virtually loaded in axial compression and the proximal femora in both fall and stance configurations. Denosumab increased vertebral body strength by 10.8%, 14.0%, and 17.4% from baseline at 12, 24, and 36 months, respectively (p < 0.0001). Denosumab also increased femoral strength in the fall configuration by 4.3%, 5.1%, and 7.2% from baseline at 12, 24, and 36 months, respectively (p < 0.0001). Similar improvements were observed in the stance configuration with increases of 4.2%, 5.2%, and 5.2% from baseline (p ≤ 0.0007). Differences between the increasing strengths with denosumab and the decreasing strengths with placebo were significant starting at 12 months (vertebral and femoral fall) or 24 months (femoral stance). Using an alternative smooth FE methodology, we confirmed the significant improvements in vertebral body and proximal femur strength previously observed with denosumab. Estimated increases in strength with denosumab and decreases with placebo were highly consistent between both FE techniques.

Behaviour of FRP confined concrete in square columns

A significant amount of research has been conducted on FRP-confined circular columns, but much less is known about rectangular/square columns in which the effectiveness of confinement is much reduced. This paper presents the results of experimental investigations on low strength square concrete columns confined with FRP. Axial compression tests were performed on ten intermediate size columns. The tests results indicate that FRP composites can significantly improve the bearing capacity and ductility of square section reinforced concrete columns with rounded corners. The strength enhancement ratio is greater the lower the concrete strength and also increases with the stiffness of the jacket. The confined concrete behaviour was predicted according to the more accepted theoretical models and compared with experimental results. There are two key parameters which critically influence the fitting of the models: the strain efficiency factor and the effect of confinement in non-circular sections.

Influência do confinamento na resistência e ductilidade de pilares curtos de concreto de ultra alta resistência submetidos à compressão centrada

Neste estudo foram analisados experimentalmente o comportamento de 24 pilares curtos de Concreto de Ultra Alta Resistência - CUAR, confinados por armaduras helicoidais, avaliando especificamente os acréscimos de resistência e ductilidade obtidos com diferentes níveis de pressão lateral de confinamento. Na etapa experimental foram realizados ensaios de pilares curtos de CUAR com as seguintes características: - seção circular de 7,2 cm de diâmetro e comprimento de 23 cm, e quatro níveis de resistência à compressão do concreto sendo eles, 165, 175, 200 e 229 MPa, dosados sem e com adição de fibras metálicas; - diferentes espaçamentos das armaduras helicoidais, de modo que fossem obtidas situações com baixo, médio e alto índice de confinamento e taxa de armadura longitudinal fixa. Os ensaios de compressão centrada foram realizados com controle de deslocamento, de modo que foram obtidas as curvas força x deslocamento completas. Constatou-se que a seção resistente dos pilares de CUAR é a formada pelo núcleo de concreto confinado, área delimitada pelo eixo da armadura transversal. Observou-se que o CUAR com fibras metálicas apresenta maior deformação do núcleo de concreto confinado em relação ao núcleo de concreto confinado de CUAR sem adição de fibras metálicas, indicando dessa forma, que os pilares de CUAR com fibras metálicas apresentam comportamento mais dúctil. Para as situações de alto confinamento foram gerados ao concreto do núcleo confinado significativos acréscimos de resistência e deformação axial, aumentando a resistência do concreto confinado em relação a resistência do concreto não confinado em: 82,26%, 75,34%, 90,46% e 70,51%, respectivamente, e as deformações axiais do concreto confinado em relação a deformação axial do concreto não confinado em: 433%, 474%, 647% e 550%. Finalmente, acredita-se que os resultados obtidos poderão trazer subsídios para aplicações futuras desta técnica de confinamento na construção de novos elementos estruturais e no reforço de pilares submetidos a elevados níveis de solicitação axial.

Statistical tests of load-unload response ratio signals by lattice solid model: Implication to tidal triggering and earthquake prediction

Statistical tests of Load-Unload Response Ratio (LURR) signals are carried in order to verify statistical robustness of the previous studies using the Lattice Solid Model (MORA et al., 2002b). In each case 24 groups of samples with the same macroscopic parameters (tidal perturbation amplitude A, period T and tectonic loading rate k) but different particle arrangements are employed. Results of uni-axial compression experiments show that before the normalized time of catastrophic failure, the ensemble average LURR value rises significantly, in agreement with the observations of high LURR prior to the large earthquakes. In shearing tests, two parameters are found to control the correlation between earthquake occurrence and tidal stress. One is, A/(kT) controlling the phase shift between the peak seismicity rate and the peak amplitude of the perturbation stress. With an increase of this parameter, the phase shift is found to decrease. Another parameter, AT/k, controls the height of the probability density function (Pdf) of modeled seismicity. As this parameter increases, the Pdf becomes sharper and narrower, indicating a strong triggering. Statistical studies of LURR signals in shearing tests also suggest that except in strong triggering cases, where LURR cannot be calculated due to poor data in unloading cycles, the larger events are more likely to occur in higher LURR periods than the smaller ones, supporting the LURR hypothesis.

Implementation of particle-scale rotation in the 3-d Lattice Solid Model

In this study, 3-D Lattice Solid Model (LSMearth or LSM) was extended by introducing particle-scale rotation. In the new model, for each 3-D particle, we introduce six degrees of freedom: Three for translational motion, and three for orientation. Six kinds of relative motions are permitted between two neighboring particles, and six interactions are transferred, i.e., radial, two shearing forces, twisting and two bending torques. By using quaternion algebra, relative rotation between two particles is decomposed into two sequence-independent rotations such that all interactions due to the relative motions between interactive rigid bodies can be uniquely decided. After incorporating this mechanism and introducing bond breaking under torsion and bending into the LSM, several tests on 2-D and 3-D rock failure under uni-axial compression are carried out. Compared with the simulations without the single particle rotational mechanism, the new simulation results match more closely experimental results of rock fracture and hence, are encouraging. Since more parameters are introduced, an approach for choosing the new parameters is presented.

Relativistic laser pulse compression in plasmas with a linear axial density gradient

The self-compression of a relativistic Gaussian laser pulse propagating in a non-uniform plasma is investigated. A linear density inhomogeneity (density ramp) is assumed in the axial direction. The nonlinear Schrodinger equation is first solved within a one-dimensional geometry by using the paraxial formalism to demonstrate the occurrence of longitudinal pulse compression and the associated increase in intensity. Both longitudinal and transverse self-compression in plasma is examined for a finite extent Gaussian laser pulse. A pair of appropriate trial functions, for the beam width parameter (in space) and the pulse width parameter (in time) are defined and the corresponding equations of space and time evolution are derived. A numerical investigation shows that inhomogeneity in the plasma can further boost the compression mechanism and localize the pulse intensity, in comparison with a homogeneous plasma. A 100 fs pulse is compressed in an inhomogeneous plasma medium by more than ten times. Our findings indicate the possibility for the generation of particularly intense and short pulses, with relevance to the future development of tabletop high-power ultrashort laser pulse based particle acceleration devices and associated high harmonic generation. An extension of the model is proposed to investigate relativistic laser pulse compression in magnetized plasmas.

A numerical method to quantify the axial shortening of vertical elements in concrete buildings

Differential axial shortening, distortion and deformation in high rise buildings is a serious concern. They are caused by three time dependent modes of volume change; “shrinkage”, “creep” and “elastic shortening” that takes place in every concrete element during and after construction. Vertical concrete components in a high rise building are sized and designed based on their strength demand to carry gravity and lateral loads. Therefore, columns and walls are sized, shaped and reinforced differently with varying concrete grades and volume to surface area ratios. These structural components may be subjected to the detrimental effects of differential axial shortening that escalates with increasing the height of buildings. This can have an adverse impact on other structural and non-structural elements. Limited procedures are available to quantify axial shortening, and the results obtained from them differ because each procedure is based on various assumptions and limited to few parameters. All these prompt to a need to develop an accurate numerical procedure to quantify the axial shortening of concrete buildings taking into account the important time varying functions of (i) construction sequence (ii) Young’s Modulus and (iii) creep and shrinkage models associated with reinforced concrete. General assumptions are refined to minimize variability of creep and shrinkage parameters to improve accuracy of the results. Finite element techniques are used in the procedure that employs time history analysis along with compression only elements to simulate staged construction behaviour. This paper presents such a procedure and illustrates it through an example. Keywords: Differential Axial Shortening, Concrete Buildings, Creep and Shrinkage, Construction Sequence, Finite Element Method.

A STUDY INTO THE BEHAVIOR OF AN ALUMINUM FOAM UNDER COMPRESSION

The characterization of a closed-cell aluminum foam with the trade name Alporas is carried out here under compression loading for a nominal cross-head speed of 1 mm/min. Foam samples in the form of cubes are tested in a UTM and the average stress-strain behavior is obtained which clearly displays a plateau strength of approximately 2 MPa. It is noted that the specific energy absorption capacity of the foam can be high despite its low strength which makes it attractive as a material for certain energy-absorbing countermeasures. The mechanical behavior of the present Alporas foam is simulated using cellular (i.e. so-called microstructure-based) and solid element-based finite element models. The efficacy of the cellular approach is shown, perhaps for the first time in published literature, in terms of prediction of both stress-strain response and inclined fold formation during axial crush under compression loading. Keeping in mind future applications under impact loads, limited results are presented when foam samples are subjected to low velocity impact in a drop-weight test set-up.