scholarly journals Viscoelastic Anisotropic Finite Element Mixture Model of Articular Cartilage using Viscoelastic Collagen Fibers and Validation with Stress Relaxation Data

2012 ◽  
Author(s):  
Matthew Alexander Griebel
2021 ◽  
Vol 22 (1) ◽  
Author(s):  
Xiaohui Zhang ◽  
Shuo Yuan ◽  
Jun Wang ◽  
Bagen Liao ◽  
De Liang

Abstract Background Recent studies have pointed out that arthroscopy, the commonly-used surgical procedure for meniscal tears, may lead to an elevated risk of knee osteoarthritis (KOA). The biomechanical factors of KOA can be clarified by the biomechanical analysis after arthroscopic partial meniscectomy (APM). This study aimed to elucidate the cartilage stress and meniscus displacement of the tibiofemoral joint under flexion and rotation loads after APM. Methods A detailed finite element model of the knee bone, cartilage, meniscus, and major ligaments was established by combining computed tomography and magnetic resonance images. Vertical load and front load were applied to simulate different knee buckling angles. At the same time, by simulating flexion of different degrees and internal and external rotations, the stresses on tibiofemoral articular cartilage and meniscus displacement were evaluated. Results Generally, the contact stress on both the femoral tibial articular cartilage and the meniscus increased with the increased flexion degree. Moreover, the maximum stress on the tibial plateau gradually moved backward. The maximum position shift value of the lateral meniscus was larger than that of the medial meniscus. Conclusion Our finite element model provides a realistic three-dimensional model to evaluate the influence of different joint range of motion and rotating tibiofemoral joint stress distribution. The decreased displacement of the medial meniscus may explain the higher pressure on the knee components. These characteristics of the medial tibiofemoral joint indicate the potential biomechanical risk of knee degeneration.


Cartilage ◽  
2021 ◽  
pp. 194760352098877
Author(s):  
Roy D. Bloebaum ◽  
Andrew S. Wilson ◽  
William N. Martin

Objective There has been a debate as to the alignment of the collagen fibers. Using a hand lens, Sir William Hunter demonstrated that the collagen fibers ran perpendicular and later aspects were supported by Benninghoff. Despite these 2 historical studies, modern technology has conflicting data on the collagen alignment. Design Ten mature New Zealand rabbits were used to obtain 40 condyle specimens. The specimens were passed through ascending grades of alcohol, subjected to critical point drying (CPD), and viewed in the scanning electron microscope. Specimens revealed splits from the dehydration process. When observing the fibers exposed within the opening of the splits, parallel fibers were observed to run in a radial direction, normal to the surface of the articular cartilage, radiating from the deep zone and arcading as they approach the surface layer. After these observations, the same samples were mechanically fractured and damaged by scalpel. Results The splits in the articular surface created deep fissures, exposing parallel bundles of collagen fibers, radiating from the deep zone and arcading as they approach the surface layer. On higher magnification, individual fibers were observed to run parallel to one another, traversing radially toward the surface of the articular cartilage and arcading. Mechanical fracturing and scalpel damage induced on the same specimens with the splits showed randomly oriented fibers. Conclusion Collagen fiber orientation corroborates aspects of Hunter’s findings and compliments Benninghoff. Investigators must be aware of the limits of their processing and imaging techniques in order to interpret collagen fiber orientation in cartilage.


2002 ◽  
Vol 124 (3) ◽  
pp. 273-280 ◽  
Author(s):  
Tammy L. Haut Donahue ◽  
M. L. Hull ◽  
Mark M. Rashid ◽  
Christopher R. Jacobs

As a step towards developing a finite element model of the knee that can be used to study how the variables associated with a meniscal replacement affect tibio-femoral contact, the goals of this study were 1) to develop a geometrically accurate three-dimensional solid model of the knee joint with special attention given to the menisci and articular cartilage, 2) to determine to what extent bony deformations affect contact behavior, and 3) to determine whether constraining rotations other than flexion/extension affects the contact behavior of the joint during compressive loading. The model included both the cortical and trabecular bone of the femur and tibia, articular cartilage of the femoral condyles and tibial plateau, both the medial and lateral menisci with their horn attachments, the transverse ligament, the anterior cruciate ligament, and the medial collateral ligament. The solid models for the menisci and articular cartilage were created from surface scans provided by a noncontacting, laser-based, three-dimensional coordinate digitizing system with an root mean squared error (RMSE) of less than 8 microns. Solid models of both the tibia and femur were created from CT images, except for the most proximal surface of the tibia and most distal surface of the femur which were created with the three-dimensional coordinate digitizing system. The constitutive relation of the menisci treated the tissue as transversely isotropic and linearly elastic. Under the application of an 800 N compressive load at 0 degrees of flexion, six contact variables in each compartment (i.e., medial and lateral) were computed including maximum pressure, mean pressure, contact area, total contact force, and coordinates of the center of pressure. Convergence of the finite element solution was studied using three mesh sizes ranging from an average element size of 5 mm by 5 mm to 1 mm by 1 mm. The solution was considered converged for an average element size of 2 mm by 2 mm. Using this mesh size, finite element solutions for rigid versus deformable bones indicated that none of the contact variables changed by more than 2% when the femur and tibia were treated as rigid. However, differences in contact variables as large as 19% occurred when rotations other than flexion/extension were constrained. The largest difference was in the maximum pressure. Among the principal conclusions of the study are that accurate finite element solutions of tibio-femoral contact behavior can be obtained by treating the bones as rigid. However, unrealistic constraints on rotations other than flexion/extension can result in relatively large errors in contact variables.


1986 ◽  
Vol 29 (4) ◽  
pp. 544-548 ◽  
Author(s):  
Joel C. Kahane ◽  
Alice R. Kahn

Collagen fiber organization in the articular surfaces of the cricoarytenoid joint (CAJ) was studied using a pinpricking technique used in biomechanical research in orthopedics. Four male human formalin preserved specimens (3 months to 20 years) and 6 male freshly autopsied specimens (19 to 30 yrs) were studied. Specimens were dissected using the stereomicroscope. Distinctive patterns of articular cartilage slits reflect the orientation of collagen fibers in the cricoid and arytenoid articular surfaces. The orientation of the collagen fibers reinforces the articular surfaces along the principle path of CAJ motion. No age related differences were found. This suggests that the orientation of collagen fibers in the CAJ articular surfaces is prenatally determined rather than significantly influenced by postnatal mechanical factors.


1996 ◽  
Vol 118 (3) ◽  
pp. 399-404 ◽  
Author(s):  
T. L Norman ◽  
V. C. Saligrama ◽  
K. T. Hustosky ◽  
T. A. Gruen ◽  
J. D. Blaha

A tapered femoral total hip stem with a debonded stem-cement interface and an unsupported distal tip subjected to constant axial load was evaluated using two-dimensional (2D) axisymmetric finite element analysis. The analysis was performed to test if the mechanical condition suggest that a “taper-lock” with a debonded viscoelastic bone cement might be an alternative approach to cement fixation of stem type cemented hip prosthesis. Effect of stem-cement interface conditions (bonded, debonded with and without friction) and viscoelastic response (creep and relaxation) of acrylic bone cement on cement mantle stresses and axial displacement of the stem was also investigated. Stem debonding with friction increased maximum cement von Mises stress by approximately 50 percent when compared to the bonded stem. Of the stress components in the cement mantle, radial stresses were compressive and hoop stresses were tensile and were indicative of mechanical taper-lock. Cement mantle stress, creep and stress relaxation and stem displacement increased with increasing load level and with decreasing stem-cement interface friction. Stress relaxation occur predominately in tensile hoop stress and decreased from 1 to 46 percent over the conditions considered. Stem displacement due to cement mantle creep ranged from 614 μm to 1.3 μm in 24 hours depending upon interface conditions and load level.


2019 ◽  
Vol 2 ◽  
pp. 52-60 ◽  
Author(s):  
Oleg Markov ◽  
Oleksiy Gerasimenko ◽  
Leila Aliieva ◽  
Alexander Shapoval

It is shown that when modeling the processes of forging and stamping, it is necessary to take into account not only the hardening of the material, but also softening, which occurs during hot processing. Otherwise, the power parameters of the deformation processes are precisely determined, which leads to the choice of more powerful equipment. Softening accounting (processes of stress relaxation) will allow to accurately determine the stress and strain state (SSS) of the workpiece, as well as the power parameters of the processes of deformation. This will expand the technological capabilities of these processes. Existing commercial software systems for modeling hot plastic deformations based on the finite element method (FEM) do not allow this. This is due to the absence in these software products of the communication model of the component deformation rates and stresses, which would take into account stress relaxation. As a result, on the basis of the Maxwell visco-elastic model, a relationship is established between deformation rates and stresses. The developed model allows to take into account the metal softening during a pause after hot deformation. The resulting mathematical model is tested by experiment on different steels at different temperatures of deformation. The process of steels softening is determined using plastometers. It is established experimentally that the model developed by 89 ... 93 % describes the rheology of the metal during hot deformation. The relationship between the components of the deformation rates and stresses is established, which allows to obtain a direct numerical solution of plastic deformation problems without FED iterative procedures, taking into account the real properties of the metal during deformation. As a result, the number of iterations and calculations has significantly decreased.


Author(s):  
Sadegh Naserkhaki ◽  
Jacob L. Jaremko ◽  
Greg Kawchuk ◽  
Samer Adeeb ◽  
Marwan El-Rich

The spinal load sharing and mechanical stresses developed in the spine segments due to mechanical loads are dependent on the unique spinal anatomy (geometry and posture). Variation in spinal curvature alters the load sharing of the lumbar spine as well as the stiffness and stability of the passive tissues. In this paper, effects of lumbar spine curvature variation on spinal load sharing under compressive Follower Load (FL) are investigated numerically. 3D nonlinear Finite Element (FE) models of three ligamentous lumbosacral spines are developed based on personalized geometries; hypo-lordotic (Hypo-L), normal (Normal-L) and hyper-lordotic (Hyper-L) cases. Analysis of each model is performed under Follower Load and developed stress in the discs and forces in the collagen fibers are investigated. Stresses on the discs vary in magnitude and distribution depending on the degree of lordosis. A straight hypo-lordotic spine shows stresses more equally distributed among discs while a highly curved hyper-lordotic spine has stresses concentrated at lower discs. Stresses are uniformly distributed in each disc for Hypo-L case while they are concentrated posteriorly for Hyper-L case. Also, the maximum force in collagen fibers is developed in the Hyper-L case. These differences might be clinically significant related to back pain.


1985 ◽  
Vol 107 (1) ◽  
pp. 231-237 ◽  
Author(s):  
A. Kaufman

A simplified inelastic analysis computer program (ANSYMP) was developed for predicting the stress-strain history at the critical location of a thermomechanically cycled structure from an elastic solution. The program uses an iterative and incremental procedure to estimate the plastic strains from the material stress-strain properties and a plasticity hardening model. Creep effects can be calculated on the basis of stress relaxation at constant strain, creep at constant stress or a combination of stress relaxation and creep accumulation. The simplified method was exercised on a number of problems involving uniaxial and multiaxial loading, isothermal and nonisothermal conditions, dwell times at various points in the cycles, different materials, and kinematic hardening. Good agreement was found between these analytical results and nonlinear finite element solutions for these problems. The simplified analysis program used less than 1 percent of the CPU time required for a nonlinear finite element analysis.


2021 ◽  
Vol 7 (2) ◽  
pp. 363-366
Author(s):  
Thomas Reuter ◽  
Christof Hurschler

Abstract Mechanical parameters of hard and soft tissues are explicit markers for quantitative tissue characterization. In this study, we present a comparison of biphasic material properties of equine articular cartilage estimated from stress relaxation (ε = 6 %, t = 1000 s) and creep indentation tests (F = 0.1 N, t = 1000 s). A biphasic 3D-FE-based method is used to determine the biomechanical properties of equine articular cartilage. The FE-model computation was optimized by exploiting the axial symmetry and mesh resolution. Parameter identification was executed with the Levenberg- Marquardt-algorithm. Additionally, sensitivity analyses of the calculated biomechanical parameters were performed. Results show that the Young’s modulus E has the largest influence and the Poisson’s ratio of ν ≤ 0.1 is rather insensitive. The R² of the fit results varies between 0.882 and 0.974 (creep model) and between 0.695 and 0.930 (relaxation model). The averaged parameters E and k determined from the creep model yield higher values in comparison to the relaxation model. The differences can be traced back to the experimental settings and to the biphasic material model.


Sign in / Sign up

Export Citation Format

Share Document