Tether pre-tension within vertebral body tethering reduces motion of the spine and influences coupled motion: a finite element analysis.

Anterior Vertebral Body Tethering (VBT) is a novel fusionless treatment option for selected adolescent idiopathic scoliosis patients which is gaining widespread interest. The primary objective of this study is to investigate the effects of tether pre-tension within VBT on the biomechanics of the spine including sagittal and transverse parameters as well as primary motion, coupled motion, and stresses acting on the L2 superior endplate. For that purpose, we used a calibrated and validated Finite Element model of the L1-L2 spine. The VBT instrumentation was inserted on the left side of the L1-L2 segment with different cord pre-tensions and submitted to an external pure moment of 6 Nm in different directions. The range of motion (ROM) for the instrumented spine was measured from the initial post-VBT position. The magnitudes of the ROM of the native spine and VBT-instrumented with pre-tensions of 100 N, 200 N, and 300 N were, respectively, 3.29°, 2.35°, 1.90° and 1.61° in extension, 3.30°, 3.46°, 2.79°, and 2.17° in flexion, 2.11°, 1.67°, 1.33° and 1.06° in right axial rotation, and 2.10°, 1.88°, 1.48° and 1.16° in left axial rotation. During flexion-extension, an insignificant coupled lateral bending motion was observed in the native spine. However, VBT instrumentation with pre-tensions of 100 N, 200 N, and 300 N generated coupled right lateral bending of 0.85°, 0.81°, and 0.71° during extension and coupled left lateral bending of 0.32°, 0.24°, and 0.19° during flexion, respectively. During lateral bending, a coupled extension motion of 0.33-0.40° is observed in the native spine, but VBT instrumentation with pre-tensions of 100 N, 200 N, and 300 N generates coupled flexion of 0.67°, 0.58°, and 0.42° during left (side of the implant) lateral bending and coupled extension of 1.28°, 1.07°, and 0.87° during right lateral bending, respectively. Therefore, vertebral body tethering generates coupled motion. Tether pre-tension within vertebral body tethering reduces the motion of the spine.

On the contribution of solid and fluid behavior to the modeling of the time-dependent mechanics of tendons under semi-confined compression.

The present work aims to investigate whether it is possible to identify and quantify the contributions of the interstitial fluid and the solid skeleton to the overall time-dependent behavior of tendons based on a single mechanical test. For this purpose, the capabilities of three different time-dependent models (a viscoelastic, a poroelastic and a poroviscoelastic) were investigated in the modeling of the experimental behavior obtained from semi-confined compression with stress relaxation tests transverse to collagen fibers. The main achieved result points out that the poroviscoelastic model was the only one capable to characterize both the experimental responses of the force and volume changes of the tissue samples. Moreover, further analysis of this model shows that while the kinematics of the sample are mainly governed by the fluid flow (pore pressure contribution of the model), the behavior intrinsically associated with the viscoelastic solid skeleton makes a significant contribution to the experimental force response. This study reinforces the importance of taking both the experimental kinematics and kinetics of tendon tissues into account during the constitutive characterization procedure.

Is the fluid volume fraction equal to the water content in tendons? Insights on biphasic modeling.

The mass density of highly hydrated soft tissues is generally assumed to be very close to that of the water, resulting that the fluid mass fraction (water content) being equal to the fluid volume fraction. Within this context, the present study aims to investigate whether such an assumption actually holds for tendon tissues and to what extent it may affect the constitutive characterizations based on biphasic (poroelastic) models. Once the water content was assessed by a classical drying assay, the fluid volume fraction was obtained based on an image segmentation approach. The main achieved results point out that the fluid volume fraction is ∼20% higher than the water content in the studied tendons (flexor digitorum profundus bovine tendons). Based on this, it is shown that the use of the water content instead of the fluid volume fraction may considerably bias the results drawn by biphasic modeling of tendons. Accordingly, a proper measurement of the fluid volume fraction is then required.

A numerical study of the contact geometry and pressure distribution along the glenoid track.

The glenoid track geometry and the contact forces acting on the glenohumeral joint at static positions of 30°, 60°, 90° and 120° of abduction with 90° of external rotation were evaluated using a finite element model of the shoulder that, differently from most usual approximations, accounts the humeral head translations and the deformable-to-deformable non-spherical joint contact. The model was based on data acquired from clinical exams of a single subject, including the proximal humerus, scapula, their respective cartilages concerning the glenohumeral joint, and the rotator cuff and deltoid muscles. The forces acting on the glenohumeral joint were estimated using a simulation framework consisting of an optimization procedure allied with finite element analysis that seeks the minimum muscle forces that stabilize the joint. The joint reaction force magnitude increases up to 680.25 N at 90° of abduction and decreases at further positions. From 60° onward the articular contact remains at the anterior region of the glenoid cartilage and follows an inferior to superior path at the posterior region of the humeral head cartilage. The maximum contact pressure of 3.104 MPa occurs at 90° abduction. Although translating inferiorly throughout the movement, the projection of the humeral head center at the glenoid plane remains at the central region of the glenoid surface. The model results qualitatively matched the trends observed in the literature and supports the consideration of the translational degrees of freedom to evaluate the joint contact mechanics.

An experimental-numerical method for the calibration of finite element models of the lumbar spine.

We present a systematic and automated stepwise method to calibrate computational models of the spine. For that purpose, a sequential resection study on one lumbar specimen (L2-L5) was performed to obtain the individual contribution of the IVD, the facet joints and the ligaments to the kinematics of the spine. The experimental data was prepared for the calibration procedure in such manner that the FE model could reproduce the average motion of the 10 native spines from former cadaveric studies as well as replicate the proportional change in ROM after removal of a spinal structure obtained in this resection study. A Genetic Algorithm was developed to calibrate the properties of the intervertebral discs and facet joints. The calibration of each ligament was performed by a simple and novel technique that requires only one simulation to obtain its mechanical property. After calibration, the model was capable of reproducing the experimental results in all loading directions and resections.

Numerical study on the effect of stent shape on suture forces in stent-grafts.

Second-generation stent-grafts (SGs) have addressed many of the mechanical problems reported for first-generation endoprostheses, such as graft tear and stent rupture; however, suture wear and detachment due to pulsatile fatigue remains an issue. Numerical studies on the mechanical behavior of these endoprostheses usually model the attachment between stents and graft as a continuous ''tie'' constraint, which does not provide information on the mechanical loads acting on individual sutures. This paper presents a suitable approach for Finite Element (FE) simulations of SGs which allows for a qualitative evaluation of the loads acting on sutures. Attachment between stents and graft is modeled as rigid beams at discrete locations of the endoprostheses, and the reaction forces on the beams are analyzed. This modeling strategy is employed for four different SG models (two Z-stented commercial models and two circular-stented models) subjected to a severe 180° U-bend, followed by intraluminal pressurization. Results show that, for all models, the majority of sutures is experiencing fluctuating forces within a cardiac cycle (between 80 and 120 mmHg), which points to pulsatile fatigue as potential failure mode. In addition, the highest loads are concentrated in kinks and, for Z-stented models, at the apexes of stents. Moreover, suture loads for circular-stented models are lower than for Z-stented models, indicating better resistance to suture detachment. All these observations are in line with experimental results published in the literature, and, therefore, the procedure herein proposed may serve as a valuable tool in the development of new SG models with better suture resistance to pulsatile wear and fatigue.

Glenoid track evaluation by a validated finite-element shoulder numerical model.

BACKGROUND: The limits of the glenoid track have been defined through methods that do not take properly into account the physiological articular forces involved in the articular contact, which may interfere with its size. Finite elements numerical models can simulate joint forces more realistically. OBJECTIVE: To evaluate the glenoid track in a finite element numerical model of the shoulder. METHODS: We developed a finite element numerical model of the shoulder, based on imaging exams of a volunteer, including the proximal humerus, scapula, their respective articular cartilages, and the rotator cuff muscles. An algorithm to balance the weight of the arm calculated muscle, wrapping, and articular reaction forces. The model has freedom of translation in three axes. The articular contact characteristics and glenoid track's dimensions according to the literature references were evaluated in 60°, 90° and 120° of abduction, all at the 90° external rotation. RESULTS: The model's anatomy and physiology were validated. The value of the glenoid track (according to Yamamoto's parameters) was 86% of glenoid length at 90° abduction before loading of forces, and 79% afterwards. The glenoid track at 60°, 90° and 120° of abduction (Omori's parameters) corresponded, respectively, to 71%, 88% and 104% of glenoid length before loading of forces, and 76%, 84% and 103% afterward. CONCLUSION: The numerical model is suitable for the shoulder articular contact evaluation. The articular contact analysis ratifies the glenoid track concept and contributes to its evolution. This value is influenced by glenohumeral joint forces, which should be considered for the analysis. LEVEL OF EVIDENCE: Basic Science Study; Computer Modelling.

On multiscale boundary conditions in the computational homogenization of an RVE of tendon fascicles.

Present study provides a numerical investigation on multiscale boundary conditions in the computational homogenization of a representative volume element (RVE) of tendon fascicles. A three-dimensional hexagonal-helicoidal finite element RVE composed of two material phases (collagen fibers and cells) and three finite strain viscoelastic models (collagen fibrils, matrix of fibers and cells) compose the multiscale model. Due to the unusual helical geometry of the RVE, the performance of four multiscale boundary conditions is evaluated: the linear boundary displacements model, the minimally constrained model and two mixed boundary conditions allying characteristics of both, linear and minimal models. Numerical results concerning microscopic kinematic fields and macroscopic stress-strain curves point out that one of the mixed models is able to predict the expected multiscale mechanics of the RVE, presenting sound agreement with experimental facts reported in literature, for example: characteristic non-linear shape of the stress-strain curves; macroscopic energy loss by hysteresis; axial rotation of fascicles observed in tensile tests; collagen fibrils are the main load-bearing components of tendons; cells contribute neither to the stiffness nor to the macroscopic energy loss. Moreover, the multiscale model provides important insights on the micromechanics of tendon fascicles, predicting a non-homogeneous and relevant strain localization on cells, even under physiological macroscopic strain amplitudes.

An experimental and numerical study on the transverse deformations in tensile test of tendons.

The transverse deformations of tendons assessed in tensile tests seems to constitute a controversial issue in literature. On the one hand, large positive variations of the Poisson's ratio have been reported, indicating volume reduction under tensile states. On the other hand, negative values were also observed, pointing out an auxetic material response. Based on these experimental observations, the following question is raised: Are these large and discrepant transverse deformations intrinsically related to the constitutive response of tendons or they result from artifacts of the mechanical test setup? In order to provide further insights to this question, an experimental and numerical study on the transverse kinematics of tendons was carried out. Tensile experiments were performed in branches of deep digital flexor tendons of domestic porcine, where the transverse displacements were measured by two high-speed, high-accuracy optical digital micrometers placed transversely to one another. Aiming at a better understanding of the effects of the mechanical test setup in the transverse measurements, a three-dimensional finite element model is proposed to resemble the tensile experiment. The main achieved results strongly support the following hypotheses regarding tensile tests of tendons: the clamping region considerably affects the kinematics of the specimen even at a large distance from the clamps; the transverse deformations are mainly ruled by stiff fibers embedded in a soft matrix; the generalization of the Poisson's ratio to draw conclusions about changes in volume of tendons may lead to misinterpretations.

A transversely isotropic coupled hyperelastic model for the mechanical behavior of tendons.

Several constitutive models for fibrous soft tissues used in literature provide a completely isotropic response when fibers are compressed. However, recent experimental investigations confirm the expectation that tendons behave anisotropically during compression tests. Motivated by these facts, the present manuscript presents an appropriate choice of hyperelastic potentials able to predict the coupled mechanical behaviors of tendons under both tensile and compressive loads with a relatively small number of material parameters. The high stiffness of tendons under tensile tests is handled by a transversely isotropic model while the coupled compressive response is modeled by means of a Fung-type potential in terms of Seth-Hill's generalized strain tensors. In present study the logarithm strain measure is used instead of the usually employed Green-Lagrange strain. After a parameter identification procedure, the resulting model showed ability to satisfactorily reproduce the experimental data. Details on the analytical material tangent modulus are provided. Present results will then enhance further researches related to tendon dissipative effects and numerical multiscale investigations.

Effect of Injection Molding Melt Temperatures on PLGA Craniofacial Plate Properties during In Vitro Degradation.

The purpose of this article is to present mechanical and physicochemical properties during in vitro degradation of PLGA material as craniofacial plates based on different values of injection molded temperatures. Injection molded plates were submitted to in vitro degradation in a thermostat bath at 37 ± 1°C by 16 weeks. The material was removed after 15, 30, 60, and 120 days; then bending stiffness, crystallinity, molecular weights, and viscoelasticity were studied. A significant decrease of molecular weight and mechanical properties over time and a difference in FT-IR after 60 days showed faster degradation of the material in the geometry studied. DSC analysis confirmed that the crystallization occurred, especially in higher melt temperature condition. DMA analysis suggests a greater contribution of the viscous component of higher temperature than lower temperature in thermomechanical behavior. The results suggest that physical-mechanical properties of PLGA plates among degradation differ per injection molding temperatures.

Influence of Processing Conditions on the Mechanical Behavior and Morphology of Injection Molded Poly(lactic-co-glycolic acid) 85:15.

Two groups of PLGA specimens with different geometries (notched and unnotched) were injection molded under two melting temperatures and flow rates. The mechanical properties, morphology at the fracture surface, and residual stresses were evaluated for both processing conditions. The morphology of the fractured surfaces for both specimens showed brittle and smooth fracture features for the majority of the specimens. Fracture images of the notched specimens suggest that the surface failure mechanisms are different from the core failure. Polarized light techniques indicated birefringence in all specimens, especially those molded with lower temperature, which suggests residual stress due to rapid solidification. DSC analysis confirmed the existence of residual stress in all PLGA specimens. The specimens molded using the lower injection temperature and the low flow rate presented lower loss tangent values according to the DMA and higher residual stress as shown by DSC, and the photoelastic analysis showed extensive birefringence.

Biomechanical performance of Bio Cross-Pin and EndoButton for ACL reconstruction at femoral side: a porcine model

Abstract Introduction The method of graft fixation is critical in anterior cruciate ligament (ACL) reconstruction surgery. Success of surgery is totally dependent on the ability of the implant to secure the graft inside the bone tunnel until complete graft integration. The principle of EndoButton is based on the cortical suspension of the graft. The Cross-Pin is based on graft expansion. The aim of this study was to evaluate the biomechanical performance of EndoButton and Bio Cross-Pin to fix the hamstring graft at femoral side of porcine knee joints and evaluate whether they are able to support of loading applied on graft during immediate post-operative tasks. Methods Fourteen ACL reconstructions were carried out in porcine femurs fixing superficial flexor tendons with Titanium EndoButton (n = 7) and with 6 × 50 mm HA/PLLA Bio Cross-Pin (n = 7). A cyclic loading test was applied with 50-250 N of tensile force at 1 Hz for 1000 cycles. The displacement was measured at 20, 100, 500 and 1000 load cycles to quantify the slippage of the graft during the test. Single-cycle load-to-failure test was performed at 50 N/mm to measure fixation strength. Results The laxity during cyclic loading and the displacement to failure during single-cycle test were lower for the Bio Cross-Pin fixation (8.21 ± 1.72 mm) than the EndoButton (11.20 ± 2.00 mm). The Bio Cross-Pin (112.22 ± 21.20 N.mm–1) was significantly stiffer than the EndoButton fixation (60.50 ±10.38 N.mm–1). There was no significant difference between Bio Cross-Pin (failure loading: 758.29 ± 188.05 N; yield loading: 713.67 ± 192.56 N) and EndoButton strength (failure loading: 672.52 ± 66.56 N; yield loading: 599.91 ± 59.64 N). Both are able to support the immediate post-operative loading applied (445 N). Conclusion The results obtained in this experiment indicate that the Bio Cross-Pin technique promote stiffer fixation during cyclic loading as compared with EndoButton. Both techniques are able to support the immediate post-operative loading applied.

Graft tendon slippage with metallic and bioabsorbable interference screws under cyclic load: a biomechanical study in a porcine model

Introduction The rupture of the anterior cruciate ligament (ACL) is the most common type of knee injury. Reconstructive surgery is the ‘gold standard’ treatment. During the immediate post-operative period, the fixation of the graft is entirely dependent on the ability of the grafted implant to be secured inside the bone tunnel under the cyclical loads associated with daily tasks. Poor fixation can lead to graft slippage, thus impairing the healing and integration of the graft. The aim of this study was to evaluate the biomechanical performance of tendon graft fixation devices with metallic and bioabsorbable interference screws. Methods Twenty ACL reconstructions were carried out in porcine tibias using deep flexor tendons to fix 9 × 20 mm metallic (n=10) and PLLA 70/30 bioabsorbable screws (n=10). To verify the ability of a construct to resist immediate postoperative (PO) rehabilitation protocols for immediate load bearing, a cyclic loading test was applied with 50-250 N of tensile force at 1 Hz for 1000 cycles, and the displacement was measured at 10, 50, 100, 500 and 1000 load cycles to quantify the slippage of the graft during the test. After the cyclic loading test, a single-cycle load-to-failure test was applied. Results The slippage of the graft using metallic screws did not differ (P = 0.616) from that observed when using bioabsorbable screws. Conclusion The results obtained in this experiment indicate that metallic screws may promote a similar amount of graft slippage during low cyclic loading as bioabsorbable screws. Additionally, there was no difference in the biomechanical performance of these two types of screws during high failure loads.

Simulação numérica tridimensional da mecânica do joelho humano / Three-dimensional numerical simulation of human knee joint mechanics

OBJETIVO: Por ser a articulação mecanicamente mais solicitada de nossa estrutura e pelo grande número de lesões associadas, motivaram a construção de um modelo tridimensional da articulação do joelho humano para simular a cinemática da articulação e obter as solicitações mecânicas nos principais ligamentos durante o movimento de flexão do joelho. Essas informações podem futuramente ser empregada como ferramenta de apoio à decisão médica em ortopedia, fornecendo subsídios na escolha do procedimento cirúrgico. MÉTODOS: Método dos Elementos Finitos foi utilizado para construir um modelo biomecânico, tridimensional, da articulação do joelho. Nesse modelo com seis graus de liberdade é aplicado movimento de flexão/extensão sendo os demais cinco graus de liberdade governados pelas interações entre os componentes da articulares. RESULTADOS: Foram obtidas informações dos movimentos, das rotações interna/externa e adução/abdução, das translações anterior/posterior, lateral/medial e superior/inferior e dos esforços nos quatro principais ligamentos articulares, no decorrer de um amplo movimento de flexão/extensão. Estes valores foram comparados, de forma qualitativa, com valores equivalentes obtidos na literatura. CONCLUSÃO: A análise de resultados permitiu observar que vários aspectos cinemáticos são satisfatoriamente reproduzidos. A pré-carga inicial dos ligamentos e o posicionamento das inserções ligamentares no modelo mostraram-se variáveis relevantes nos resultados.

OBJECTIVE: The knee joint is the part of our structure upon which most mechanical demands are placed and a large number of lesions are associated to it. These factors motivated the construction of a three-dimensional model of the human knee joint in order to simulate joint kinematics and obtain the mechanical demands on the main ligaments during knee flexion movements. METHODS: The finite elements method was used to build a three-dimensional, biomechanical model of the knee joint. In this model with six degrees of freedom, the flexion/extension movement is applied, while the other five degrees of freedom are governed by the interactions between joint components. RESULTS: Data was collected on the movements, on the internal/external and adduction/ abduction rotations, on the anterior/posterior, lateral/medial and upper/lower translations, and on the forces acting upon the four main joint ligaments, during a wide flexion/extension movement. These values were qualitatively compared with comparable values available in the literature. CONCLUSIONS: It was observed through an analysis of the results that several kinematic aspects are satisfactorily reproduced. The initial pre-load of the ligaments and the positioning of the ligament insertions in the model were shown to be relevant variables in the results.