Biomechanical Prediction of Crack Angle Effect on Local Tensile Stress of Articular Cartilage

Jia Yu, Meng Zhang, Wen Zhang and Zong-Ping Luo*

Department of Orthopaedics, Medical College, Soochow University, China

*Corresponding Author:
Zong-Ping Luo
Department of Orthopedics, The First Affiliated Hospital of Soochow University
Orthopedic Institute, Medical College, Soochow University, 708 Renmin Rd, Suzhou, Jiangsu, 215007 China
E-mail: [email protected]

Received date: October 10, 2016; Accepted date: October 11, 2016; Published date: October 15, 2016

Citation: Luo ZP. Biomechanical Prediction of Crack Angle Effect on Local Tensile Stress of Articular Cartilage. J Bone Rep Recommendations. 2016, 2:4.Doi: 10.4172/2469-6684.100030

Copyright: © 2016 Luo ZP, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Visit for more related articles at Journal of Bone Reports & Recommendations


Purpose: Articular cartilage is an important load-bearing structure in synovial joint. Posttraumatic osteoarthritis may occur after cracking insult to cartilage, while not all insult to cartilage could progress into osteoarthritis. The initiating mechanism of osteoarthritis is still not well understood. Tensile stress is considered to have strong relationships with cartilage degradation. We hypothesize that cracks on articular cartilage under compression will induce tensile stress around the crack region, which may lead to chondrocytes change and cartilage degeneration.

Methods: In this study, using a 2D finite element model and animal validation study, the relationship between crack angle and local tensile stress was investigated.

Results: Under pure compressive load, 45 degrees oblique crack on cartilage could induce the highest peak tensile stress (3.9 times to compressive pressure) while no tensile stress was found in perpendicular crack cartilage. Animal study also confirmed finite element prediction that oblique crack may produce type I collagen around crack tip region.

Conclusions: This finding implies that crack angle could be an important factor contributing to initiate osteoarthritis after a traumatic injury. New strategies for prevention and early treatment of osteoarthritis might benefit from this study.


Synovial joint; electrolytes; Osteoarthritis; Chondrocytes; Traumatic injury; Type I collagen; Ligaments; Type II collagen


The normal mature articular cartilage is smooth and lubricated, functioning as almost frictionless load-bearing surface in synovial joints, mainly cushioning compressive load. Posttraumatic osteoarthritis (PTOA) may occur after traumatic injury to the joint and around 12% of the overall symptomatic OA burden can be attributed to PTOA of main joints [1]. However, cartilage is an anural tissue where events such as fibrillation and cracking of its structure cannot be detected by an individual. Besides, a traumatic injury such as intra-articular fracture might reduce its integrity and lead to osteoarthritis [2]. The risk of OA following a significant joint trauma that disrupts the articular surface has been reported to be much higher than the risk of OA in those who did not have a traumatic injury. However, not all the insult to articular cartilage lead to osteoarthritis; less than half of traumatic insult to cartilage would cause osteoarthritis-like changes [3].

An insult such as impact on articular cartilage or intra-articular trauma may produce crack lesions which may lead to a different local stress environment. Alteration of the local mechanical environment around the injury site might determine early OA progress since mechanical factors have been strongly implicated in the development of OA, although exact mechanisms of the initial events in cartilage degradation lead to OA are still not elucidated. PTOA degeneration is associated with loss of type II collagen, which is the principle component of the extracellular matrix of articular cartilage, being replaced by type I collagen. Type I collagen is typically found in ligaments and tendons since type I collagen fibrils have enormous tensile strength needed in these structures. Change in tensile stress has been related to osteoarthritis progress but without conclusive evidence.

Crack lesions are the typical result of insults to articular cartilage [4]. In-depth investigation of early tensile stress of crack lesions could offer insight into OA initiation. We hypothesize that crack on articular cartilage under compression will induce tensile stress around the crack region, which may lead to chondrocytes change and cartilage degeneration. In this study, we will use finite element analysis to explore local tensile stress around crack region and validate simulation predictions with animal model observation.


Finite element model

A simplified two-dimensional rectangular model was generated with Abaqus 6.11 finite element package (Simulia, USA), representing coronal cross-section of articular cartilage of a rabbit distal femur. The 2D rectangular model size was 4.5 mm × 0.50 mm according to the normal sample size.

A 4-node bilinear plane element (CPE4R) which has quadrilateral strain and reduced integration with hourglass control was used in the model because its good compressive and bending behaviour according to ABAQUS user’s manual. The element size was set at 10 μm × 10 μm after model convergence verification. The total mesh number is about 18500. The rectangular model was assumed as homogenous and isotropic. Elastic Young’s modulus (E=8MPa) and Poisson ratio (v=0.42) was assigned according to literature [5]. The bottom edge of rectangle was constrained in vertical direction and the top surface was uniformly compressed with 0.15 MPa pressure to simulate normal loading to knee joint during walking [6].

To simulate a cut-lesion on the articular cartilage, a crack effect was added. The crack length was 300 ± 25.0 μm as it was measured from cross-sectional film of animal knee study. To examine the effect of crack angle, the length was set as 300 μm and various crack angles from 0 to 90 degrees, at interval of 15 degrees, were inflicted to the middle node of top surface, as depicted in (Figure 1). The crack itself was set as frictionless contact.


Figure 1: Loading condition for 2D cartilage finite element model.

Animal Study

Animal model

Four months old New Zealand white rabbits with a weight between 2.5 and 3.0 kg were obtained from the central animal department breeding facility of the Soochow University. The study was approved by the Institutional Animal Care and Use Committee. They were kept in pairs in cages with sawdust bedding and had free access to a standard diet and tap water.

Fifteen rabbits were randomly divided into 3 groups. Rabbits were anesthetized by chloral hydrate (0.1 g/mL) by vein injection with the animals lying supine on the operating table. After shaving the knee joint, the skin of their hind legs was disinfected with entoiodine for 3 min, and then a Para patellar skin incision was made on the medial side of the joint. The patella was then dislocated laterally to expose the fossa intercondyloidea of the distal femur. The crack lesion insult was carefully made in the middle of the femoral trochlear groove using the tip of a shave knife, and the crack lesion was angled in two directions: one perpendicular to the articular surface (90°) and the other oblique (45°). Each rabbit only received one crack lesion. The knee joint was irrigated with 0.9% sodium solution after surgical operation. The patella was repositioned and the capsule and skin were sutured carefully. After operation procedures, all the animals received antibiotics (penicillin) to prevent infection.

Histology and immunohistochemistry

For routine histology, animals (n=5) were killed via overdose of chloral hydrate after 2, 6 and 20 weeks. Whole knee joints were dissected and fixed in 10% neutral buffered formalin for two days. The samples were then dehydrated and embedded in paraffin after decalcification in buffered 10% EDTA of pH7.4 for 4 weeks. Afterwards, horizontal plane sections of femoral trochlear groove (5 μm) were cut into slides by a microtome (RM2165, Leica, Nussloch, Germany).

Type I and type II collagen were then evaluated immunehistochemically. Cartilage slides were deparaffinised in xylene and incubated with testicular hyaluronidase (2 mg/mL) (H3506-1G, Sigma-Aldrich Co, St. Louis, USA) for 60 min at 37°C. Afterwards, the slides were treated with 1% H2O2 for 30 min at room temperature (RT) to eliminate endogenous peroxidase activity, and 1.5% horse serum for 30 min at RT to eliminate nonspecific protein absorption. The slides were put exposure to the primary antibody at 37°C. The primary antibodies were anti-collagen type I or type II, 1 mg/mL (Ab6308, Abcam, San Francisco, USA), diluted with phosphate-buffered saline. The samples were washed with PBS and used by the secondary antibodies (PK6200, Vector Laboratories, Burlingame, USA) for 60 min and to avidin at RT. Color was developed with 0.3% diaminobenzidine tetrahydrochloride (Sk4100, Vector Laboratories, Burlingame, USA).


Different cracking angles will produce a different tensile stress around the crack tip (Figure 2). The obtained results showed that pure compression load of 0.15 MPa pressure would result into 1.5% strain at the superior surface of articular cartilage and there was almost uniform strain of 1.5% of whole intact cartilage. At pure compression, a 90 degrees crack did not induce tensile stress, while the highest peak tensile stress (0.59 MPa, or 3.9 times of compressive pressure loading) was reached with a 45 degrees cracking (Figure 3). With a cracking of 45 degrees, the highest peak tensile strain was predicted, which is 6.81% or 4.54 times of compressive load strain on the surface of articular cartilage (Figures 4 and 5).


Figure 2: Tensile stress distributions at different crack angles. a. 15 degrees crack. b. 30 degrees crack. c. 45 degrees crack. d. 60 degrees. e. 75 degrees crack. f. 90 degrees crack.


Figure 3: a) Peak tensile stress at different crack angles. b) Ratio of peak tensile stress to compressive pressure varied with different crack angles.


Figure 4: Tensile strain distributions with different crack angles. a) 15 degrees crack. b) 30 degrees crack. c) 45 degrees crack. d) 60 degrees. e) 75 degrees crack. f) 90 degrees crack.


Figure 5: Peak tensile stress varied with different crack angles.

Animal validation

In the animal study, immunohistochemical results of type I collagen showed the 45 degrees crack had 0.13 ± 0.07 of the mean optical density around the crack at 2 weeks increasing progressively without statistically difference to 0.18 ± 0.05 at 6 weeks, and to 0.19 ± 0.09 at 20 weeks (Figure 6). In comparison, the type II collagen and proteoglycan did not change significantly. At perpendicular crack condition, Type I collagen was not found, and type II collagen and proteoglycan were unchanged at 2, 6 and 20 weeks.


Figure 6: Immunohistochemical staining of type I collagen showed increase from 2 weeks to 20 weeks around the 45 degrees crack but no presentation around the 90o crack. a) 45 degrees crack at 2 weeks. b) 45 degrees crack at 6 weeks. c) 45 degrees crack at 20 weeks. d) 90 degrees crack at 2 weeks. e) 90 degrees crack at 6 weeks. f) 90 degrees crack at 20 weeks.


A characteristic feature of OA initiation is cracking, fibrillation, and wear of cartilage; those could be treated as a mechanically driven process. Tensile failure of cartilage has been of particular interest, since it was generally believed that vertical cracks in cartilage were initiated by relatively high tensile stresses on the articular surface. In this study was found that tensile stress has strong relationship with crack angle. The oblique crack insult would induce remarkable tensile stress and increased strain around the crack tip which could influence feedback loops of cartilage degradation process.

Cartilage is composed of fluid, electrolytes, chondrocytes, collagen fibers, proteoglycans and other glycoproteins [5,7,8]. The dense collagen fiber network provides load-bearing support and tensile strength to the diarthrodial joints. OA is the progress of degradation and loss of articular cartilage, subchondral bone remodeling, and chronic inflammation of the synovial membrane.

Its etiology, although is not well understood. OA is considered to result from a complex system of interacting mechanical, biological, biochemical, molecular and enzymatic feedback loops [9].

While studies via computational models of joint contact show that the tensile stress on the surface is quite low, still some tensile stress exists within the cartilage. Under normal physiologic condition, articular cartilage would have 1-3% compressive strain [10] and normal synovial joints have a relatively low coefficient of friction of about 0.001. In this study, with a compressive strain set at 1.5% (which is within the normal range), we found that under pure compression load without any shear stress, the oblique crack could induce tensile stress around crack tip as well. This finding implies that crack angle could be a contributing factor to initiate OA after a traumatic injury.

Previous studies have showed that 2D FE models can simultaneously predict stress distribution and reaction forces during swelling, confined compression, indentation and unconfined compression of articular cartilage. In this study, 2D model was considered to be accurate enough to predict stress under pure compression. Other studies related to torsion and shear loading could benefit of implementation of 3D models.

This study present some disadvantages. The mechanical response of cartilage has relationship with the flow of fluid through the tissue. When deformed, fluid flows through the cartilage and across the articular surface. Recognizing that fluid flow and deformation are interdependent and has led to the modeling of cartilage as a mixture of fluid and solid components, which is a biphasic model [11,12]. Under pure compression static load, single-phase elastic solid model is relative accurate to describe the behavior of cartilage with crack insult. Thus, considering the applied loading condition, we did not consider the inclusion of fluid into this quasi-static model.

The cartilage matrix is occupied by proteoglycan aggregates, which is limited by the entangling collagen framework. When cartilage is compressed by loads, the negatively charged sites on aggrecan are pushed closer together, which increases their mutual repulsive force and adds to the compressive stiffness of the cartilage [13]. Damage to the collagen framework also reduces the compressive stiffness of the tissue, since the aggregated proteoglycans are contained less efficiently. Type X collagen has been shown to be up-regulated in osteoarthritic cartilage, and also by mechanical loading in normal cartilage [14]. In this study, detailed micro structural of cartilage was not modeled, therefore the effect of aggregated proteoglycans and constraint of fiber on tensile stress were not considered, which may underestimate maximum tensile stress in cartilage.

In this study animal models of cut insult on articular cartilage were used. Animal models are used to provide a controlled environment for studying the progression of OA. Although OA may be induced by numerous means, models based on disruption of the mechanical environment of the joint, either by surgical alteration of periarticular structures or by abnormal joint load, are commonly used and benefit from its implementation [15].


Oblique cracks may remarkably induce local tensile stress to amplitudes of several times of loading and could lead to formation or type I collagen. This finding suggests that tensile stress is sensitive to crack angle on articular cartilage, which might be one of the mechanism to explain the initiation of PTOA.

Conflicts of Interest Statement

The authors hereby certify that the work which is reported in the manuscript has not received financial support from any company or other commercial source, and the authors do not have any special financial interest with other organization.


This work was supported by the National Natural Science Foundation of China (81320108018, 31270995 and 11572211), the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).


Select your language of interest to view the total content in your interested language

Viewing options

Post your comment

Share This Article

Flyer image

Post your comment

captcha   Reload  Can't read the image? click here to refresh