# Near field effect on elasticity measurement for cartilage-bone structure using Lamb wave method

- Hao Xu
^{1}, - Shigao Chen
^{2}, - Kai-Nan An
^{3}and - Zong-Ping Luo
^{1}Email author

**16**:123

https://doi.org/10.1186/s12938-017-0417-9

© The Author(s) 2017

**Received: **29 June 2017

**Accepted: **26 October 2017

**Published: **30 October 2017

## Abstract

### Background

Cartilage elasticity changes with cartilage degeneration. Hence, cartilage elasticity detection might be an alternative to traditional imaging methods for the early diagnosis of osteoarthritis. Based on the wave propagation measurement, Shear wave elastography (SWE) become an emerging non-invasive elasticity detection method. The wave propagation model, which is affected by tissue shapes, is crucial for elasticity estimating in SWE. However, wave propagation model for cartilage was unclear.

### Methods

This study aimed to establish a wave propagation model for the cartilage-bone structure. We fabricated a cartilage-bone structure, and studied the elasticity measurement and wave propagation by experimental and numerical Lamb wave method (LWM).

### Results

Results indicated the wave propagation model satisfied the lamb wave theory for two-layered structure. Moreover, a near field region, which affects wave speed measurements and whose occurrence can be prevented if the wave frequency is larger than one critical frequency, was observed.

### Conclusion

Our findings would provide a theoretical foundation for further application of LWM in elasticity measurement of cartilage in vivo. It can help the application of LWM to the diagnosis of osteoarthritis.

## Keywords

## Background

Osteoarthritis is the most prevalent joint disease in senior people worldwide [1]. For end stage osteoarthritis, joint replacement surgery and pain management are the main clinical treatment options available due to the poor cartilage regeneration potential [2]. Therefore, early treatment, including preventing and even reversing the pathological changes of osteoarthritis, is desired. Fortunately, laboratory tests have partly verified the feasibility of early treatment [3]. Early treatment requires the early detection of osteoarthritis. However, early diagnosis of osteoarthritis cannot be achieved with traditional imaging diagnostic tools (i.e. X-ray and MRI) [4]. It is well established that the initial pathologic change observed in cartilage degeneration is focal glycosaminoglycan depletion. It changes the elastic properties of cartilage whereas the morphological appearance remains intact [5]. Previous studies have reported that a significant difference was observed in the surface elastic modulus of osteoarthritis (OA) grades 1 and 2 [6]. Therefore, compared to traditional imaging detection methods, elasticity detection may have the potential to early diagnose of osteoarthritis.

Shear wave elastography (SWE) is a newly developed non-invasive detection method based on elasticity detection [7, 8]. It has been widely used in the elasticity detection of various tissue lesions, such as liver fibrosis [7], breast lesions [7], and thyroid nodules [8]. However, commercial SWE is not acceptable for the elasticity diagnosis of cartilage owing to two major challenges. The first challenge is the difficulty of generation and measurement of mechanical waves in cartilage tissues. Commercial SWE machines can only detect tissues with Young’s modulus values less than 0.3 MPa [9, 10]. However, the Young’s modulus of cartilage exceeds the upper limit of commercial SWE machines [6]. In high stiffness tissues, mechanical wave motions generated by the ultrasound radiation force are too weak to be detected by the ultrasound probe. Therefore, a higher driving force and/or more sensitive measurement technology is required for cartilage elasticity measurement [7]. The second challenge is the appropriate wave propagation model in cartilage-bone structure. For tissue elasticity measurement with thin thicknesses (e.g. the myocardium [11] and arterial wall [12, 13]), the Lamb wave theory in single-layered structures is acceptable. Although the cartilage thickness is also thin, the Lamb wave propagation in a cartilage-bone structure is more complicated than that of myocardium because of the existence of the bone. Therefore, the Lamb wave propagation model in a cartilage-bone structure should be investigated.

The objective of this study is to develop the Lamb wave method (LWM) for elasticity measurement in cartilage-bone structure. Current study would provide a theoretical foundation for further application of LWM to cartilage in vivo. We fabricated a silicone-polypropylene plate to simulate the cartilage-bone structure. The lamb waves were generated by a mechanical shaker under different loading frequencies, and detected by a laser Doppler detector. The measured wave speed in the cartilage-bone structure was used to estimate silicone elasticity (simulated cartilage) using the Lamb wave propagation model. Numerical simulations were also carried out to deeply understand wave propagation in a cartilage-bone structure.

## Methods

### Simulated cartilage-bone structure preparation

To create an ideal model for the LWM test, a silicone-polypropylene plate was fabricated to simulate the actual cartilage-bone structure. The silicone plate was used to simulate the cartilage structure. A SYLGARD 184 silicone elastomer and a curing agent (Dow Corning, Midland, Michigan, USA) were mixed together in a mass ratio of 10:1 at room temperature. The mixture was poured into molds, degassed in a vacuum chamber for 20 min, and then maintained at 60 °C for 2 h to accelerate the curing process. The mold was removed when the silicone mixture was fully cured. From the mold, we obtained a silicone plate (0.8 cm × 20 cm × 20 cm). Then, the silicone plate was stuck to the polypropylene plate (3 cm × 20 cm × 20 cm) simulating the bone. Further, a cylinder silicone sample (diameter 19 mm and height 28 mm) and polypropylene sample (diameter 1.4 cm and height 3 cm) were fabricated for the uniaxial compression test. The measured density of the silicone and polypropylene samples were 1239 and 838 kg/m^{3}, respectively.

### Uniaxial compression

*E*= 2(1 +

*ν*)

*G*is the Young’s modulus,

*ν*is Poisson ratio and

*G*is shear modulus. H is the sample thickness, and S is the contact area. The Poisson ratio of SYLGARD 184 ranges from 0.45 to 0.5 [14].

### Lamb wave theory in single and two-layered structure

_{L}and \({\text{c}}_{\text{T}} = \sqrt {\mu /\rho }\) are longitudinal and transverse wave speeds, respectively, and 2 h is the thickness of the single layer structure. The frequency response of the bio-tissues could be modeled by Voigt model, so the shear modulus \(\mu = \mu_{ 1} + i\omega \mu_{ 2} ,\) where

*μ*

_{1}and

*μ*

_{2}are the shear elastic and viscous modulus. In current study, the viscosity coefficient

*μ*

_{2}was ignored for the silicone and polypropylene samples. To obtain cartilage elasticity coefficient in pure cartilage structure, Eq. (2) is fitted to the Lamb wave dispersion curves (wave speed versus frequency).

*u*

_{1}and

*u*

_{2}are displacement,

*σ*

_{12}and

*σ*

_{22}are stress, and the superscript m and n are the layer number. The field matrix g can be found in the literature [16]. The coefficients of matrix g are:

_{L}and \({\text{c}}_{\text{T}} = \sqrt {\mu /\rho }\) are longitudinal and transverse wave speeds, respectively. The frequency response of the bio-tissues could be modeled by Voigt model, so the shear modulus \(\mu = \mu_{ 1} + i\omega \mu_{ 2}\). In current study, the viscosity coefficient

*μ*

_{2}was ignored for the silicone and polypropylene samples.

To obtain cartilage elasticity coefficient in cartilage-bone structure, Eq. (8) is fitted to the Lamb wave dispersion curves (wave speed versus frequency).

### Experimental LWM tests

*x*

_{ i }as the distance between the shaking point and the Doppler detector point, where the subscript i means the ith measurement. In the present tests, the receiver distances varied from 2.5 to 7.5 cm (0.5 cm increments). The mean receiver distance could be defined as

Hence, the mean receiver distance in the experimental LWM tests was about 5 cm. Based on the particle velocity measured from different points, the wave speed could be estimated by the 2D fast Fourier transform (FFT) algorithm [11].

### Numerical LWM (finite element method)

A 2D axisymmetric finite element method (FEM) was used to investigate lamb wave propagation in the cartilage-bone plate. In the present study, the FEM package was the ABAQUS/CAE (SIMULIA, Providence, RI). The plate was represented with an axisymmetric element (CAX4R). The structure shape, material parameters, and the wave speed measurement method were established according to the experimental LWM tests described above. A fixed boundary condition was prescribed to the lateral surface of the cartilage-bone structure. A cyclic motion with a maximum displacement of 10 μm was applied to the middle of the cartilage surface as a point source. The cyclic motion had a loading frequency ranging from 200 to 3000 Hz.

## Results

Shear modulus estimated by uniaxial compression, lamb wave method in cartilage-bone and pure cartilage structure

Method | Shear modulus (MPa) |
---|---|

Uniaxial compression | 0.77 ± 0.18 |

Lamb wave method (cartilage-bone structure) | 0.86 ± 0.06 |

Lamb wave method (pure cartilage structure) | 0.80 ± 0.03 |

Figure 4 shows two significant difference between the pure cartilage and cartilage-bone structure. The wave dispersion curve is the first difference between the pure cartilage and cartilage-bone structure (Fig. 4). In the pure cartilage structure, the wave speed increased slowly and approached a steady value under the higher wave frequency (about 24.2 m/s). In the cartilage-bone structure, with a higher wave frequency, the wave speed increased sharply to a peak speed (about 165 m/s), then dropped to a steady value (about 25.1 m/s). Results demonstrate that the existence of bone had a significant effect on the wave speed in cartilage. Nevertheless, both the steady values in pure cartilage and cartilage-bone structure seemed to be approximately the same as the value of the estimated Rayleigh wave speed (\({\text{c}} = \frac{1}{1.05}\sqrt {\mu /\rho } = 23.74 {\text{m}}/{\text{s }}\) [17]). Actually, the lamb wave speed converges to the Rayleigh wave speed as the frequency increases [18]. The frequency range of wave speed measurement is the another difference between the pure cartilage and cartilage-bone structure (Fig. 4). For the pure cartilage structure, the wave speed in the entire frequency region could be measured; however, for the cartilage-bone structure, the wave speeds could only be estimated when the frequency was higher than the critical frequency *f*
_{
critical
} 1000 Hz. Two questions emerged from these observations: in cartilage-bone structure, why could the wave speed be estimated only when the frequency was higher than the critical frequency, and what may affect the critical frequency? It is crucial for the elasticity measurement in actual cartilage-bone structure.

*f*

_{ critical }1000 Hz (Fig. 5c, d). The wave speeds could only be estimated when the loading frequency was higher than the critical frequency \(f_{critical}\) of 1000 Hz.

*x*

_{ mean }of 25 cm; and (2) mean receiver distance

*x*

_{ mean }of 1.5 cm. Figure 6a compares the relationship between the wave speed and the frequency under different mean receiver distances. It can be seen that the higher the mean receiver distance

*x*

_{ mean }, the lower the critical frequency. We speculate that this phenomenon may be due to the near field effect. It is well known that a near field region that is associated with the wavelength λ exists in wave propagation. When a receiver is located too close to the driving force point, wave propagation is impacted by near field effects. The task was finding a way to estimate the critical frequency theoretically. Based on the theoretical model in Fig. 6a, we estimated the frequency–wavelength curve in Fig. 6b. In addition, the numerical LWM values in Fig. 6a were also plotted in Fig. 6b as critical frequency vs. 2

*x*

_{ mean }. A good agreement was observed. Therefore, the receiver is located in the near field region when the wavelength \(\left( {\lambda = c/f} \right)\) is larger than 2

*x*

_{ mean }, i.e.

## Discussion

Elasticity measurement of simulated cartilage in cartilage-bone structure was developed using experimental LWM tests. The measured shear modulus were in good agreement with uniaxial compression tests. In addition, a near field region was found to exist in the LWM tests. When the receiver was located in the near field region, the wave speed estimated by the current method was invalid. To prevent occurrence of the near field effect and measure a valid wave speed, the wave frequency should be higher than the critical frequency. The critical frequency was about the ratio of wave speed and 2*x*
_{
mean
} in current study. That is, the shorter mean receiver distance *x*
_{
mean
} required a higher critical wave frequency.

The maximum Young’s modulus for the commercial SWE machine is about 0.3 MPa, and the simulated cartilage in the present study had a Young’s modulus of about 2.27 MPa. Hence, the simulated cartilage elasticity exceeded the maximum value of the commercial SWE machine, i.e., the commercial SWE machines (such as the SSI [9] or VTIQ [10]) could not be used in the present study. To study wave propagation in a cartilage-bone structure with higher elasticity, a mechanical shaker and laser Doppler detector were adopted to enhance the driving force and measurement sensitivity [19], respectively. However, this method cannot be translated to elasticity measurement of cartilage in vivo. In spite of this limitation, this method could be used to study the wave propagation model [17] and the near field effect, and could establish the theoretical foundation for the LWM in actual cartilage. Nevertheless, it remains challenging to set up the measurement machine for a cartilage-bone structure. In the future, it would be preferable to use ultrasound [7] or MRE [20] to measure wave motion. Certainly, the measurement’s accuracy and driving force should be increased [7].

The cartilage-bone structure was simulated by silicon-polypropylene plate in the present study. In terms of material properties, the actual cartilage is the porous hydrous material, and its elasticity (Young’s modulus 5.98–26.51 MPa [6]) is slightly higher than the elasticity of simulated cartilage (Young’s modulus 2.27 MPa). In terms of the structure’s shape, the actual cartilage has a curved surface and an uneven thickness. In spite some distinctions, current simulated silicon-polypropylene plate captured the main features of actual cartilage-bone structure, i.e., the actual cartilage-bone structure has a two-layered structure, and the elastic modulus of bone is much higher than cartilage. Therefore, the wave propagation model and the near field effect established in present study are reliable. The theoretical model developed in this study can also be used in LWM of actual cartilage.

Like the arterial wall and myocardium, cartilage has a “thin thickness”. Even so, there exist two important differences. First, the wave propagation models in single and two-layered structures were significant different because of the existence of bone. In a pure cartilage structure, which is similar to the arterial wall and myocardium (single-layered structure), the wave speed increased and approached the Rayleigh wave speed at a higher wave frequency [11]; in a cartilage-bone structure (two-layered structure), at a higher wave frequency, the wave speed increased sharply to a peak speed and then dropped to the Rayleigh wave speed. Secondly, the near field effect significantly affect the wave speed measurement in the cartilage-bone structure. The measured wave speed was invalid when the wave frequency was lower than the critical frequency. In contrast, this phenomenon was not observed in the pure cartilage structure (single-layered structure). The study of the critical frequency in single-layered structure is beyond the scope of this study.

The present study reveals the near field effect of LWM in a cartilage-bone structure. However, this study neglected certain factors that might influence the near field effect, and these factors should be explored in the future. First, as documented in various publications, soft tissues are inherently viscoelastic [21]. A large viscosity is linked to a large wave speed. To maintain the mean receiver distance, the critical frequency must be increased according to Eq. (10). Hence, neglecting the cartilage viscosity may underestimate the critical frequency. Second, the Poisson’s ratio could affect the longitudinal and transverse wave speed, which would change the near field region and affect the critical frequency. In future LWM studies of cartilage, the actual Poisson’s ratio of cartilage should be studied.

When the receiver was located in the near field region, a near field effect emerged in the LWM tests. It caused a deviation between the measured wave speed and the theoretical model. However, current wave inversion methods (2D FFT method) take no account of the near field effect. Hence, the wave frequency must be higher than the critical frequency to avoid the near field effect. If we have to measure the wave speed in the near field in the future, a more advanced wave inversion method that taking account of the near field effect will be required.

## Conclusion

Lamb wave method shows a potential for elasticity detection of cartilage. In this study, the elasticity measurement and wave propagation model in a cartilage-bone structure was investigated by experimental and numerical LWM tests. Results showed that the wave propagation model satisfied the lamb wave theory for a two-layered structure. Moreover, there existed the near field effect in the cartilage-bone structure. To prevent the occurrence of the near field effect and improve the validity of the wave speed measurement, the wave frequency should be higher than the critical frequency. The present study established a theoretical foundation for experimental LWM tests in cartilage, and it is the first step, to the best of the author’s knowledge, towards the ultimate goal of using LWM to diagnose osteoarthritis.

## Declarations

### Authors’ contributions

HX carried out protocol development, data collection, and data analysis. SGC, KNA and ZPL handled project development. All authors read and approved the final manuscript.

### Acknowledgements

The authors thank Jian-Qiang Mo for preparing the simulated cartilage structure.

### Competing interests

The authors declare that they have no competing interests.

### Availability of data and materials

All data analysed during this study are included in this published article.

### Funding

This project is funded by the National Natural Science Foundation of China (81320108018, 31570943 and 31270995), Natural Science Foundation of Jiangsu Province (CN) (BK20160305), Innovation and Entrepreneurship Program of Jiangsu Province, and Priority Academic Program Development of Jiangsu Higher Education Institutions.

### Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access**This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

## Authors’ Affiliations

## References

- Zhang Y, Jordan JM. Epidemiology of osteoarthritis. Clin Geriatr Med. 2011;26:355–69.View ArticleGoogle Scholar
- Sarzi-Puttini P, Cimmino MA, Scarpa R, Caporali R, Parazzini F, Zaninelli A, et al. Osteoarthritis: an overview of the disease and its treatment strategies. Semin Arthritis Rheum. 2005;35:1–10.View ArticleGoogle Scholar
- Anderson DD, Chubinskaya S, Guilak F, Martin JA, Oegema TR, Olson SA, et al. Post-traumatic osteoarthritis: improved understanding and opportunities for early intervention. J Orthop Res. 2011;29:802–9.View ArticleGoogle Scholar
- Palmer AJR, Brown CP, McNally EG, Price AJ, Tracey I, Jezzard P, et al. Non-invasive imaging of cartilage in early osteoarthritis. Bone Joint J. 2013;95(6):738–46.View ArticleGoogle Scholar
- Knecht S, Vanwanseele B, Stüssi E. A review on the mechanical quality of articular cartilage—implications for the diagnosis of osteoarthritis. Clin Biomech. 2006;21:999.View ArticleGoogle Scholar
- Niu H, Liu C, Li A, Wang Q, Wang Y, Li D, et al. Relationship between triphasic mechanical properties of articular cartilage and osteoarthritic grade. Sci China Life Sci. 2012;55:444–51.View ArticleGoogle Scholar
- Sarvazyan AP, Urban MW, Greenleaf JF. Acoustic waves in medical imaging and diagnostics. Ultrasound Med Biol. 2013;39(7):1133–46.View ArticleGoogle Scholar
- Zaleska-Dorobisz U, Kaczorowski K, Pawluś A, Puchalska A, Inglot M. Ultrasound elastography—review of techniques and its clinical applications. Adv Clin Exp Med. 2014;23:645–55.View ArticleGoogle Scholar
- Zhang ZJ, Fu SN. Shear elastic modulus on patellar tendon captured from supersonic shear imaging: correlation with tangent traction modulus computed from material testing system and test-retest reliability. PLoS ONE. 2013;8:e68216.View ArticleGoogle Scholar
- Liu H, Zhao LX, Xu G, Yao MH, Zhang AH, Xu HX, et al. Diagnostic value of virtual touch tissue imaging quantification for benign and malignant breast lesions with different sizes. Int J Clin Exp Med. 2015;8:13118.Google Scholar
- Nenadic IZ, Urban MW, Mitchell SA, Greenleaf JF. Lamb wave dispersion ultrasound vibrometry (LDUV) method for quantifying mechanical properties of viscoelastic solids. Phys Med Biol. 2011;56(7):2245–64.View ArticleGoogle Scholar
- Maksuti E, Widman E, Larsson D, Urban MW, Larsson M, Bjällmark A. Arterial stiffness estimation by shear wave elastography: validation in phantoms with mechanical testing. Ultrasound Med Biol. 2016;42:308–21.View ArticleGoogle Scholar
- Couade M, Pernot M, Prada C, Messas E, Emmerich J, Bruneval P, et al. Quantitative assessment of arterial wall biomechanical properties using shear wave imaging. Ultrasound Med Biol. 2010;36:1662–76.View ArticleGoogle Scholar
- Du P, Lin IK, Lu H, Zhang X. Extension of the beam theory for polymer bio-transducers with low aspect ratios and viscoelastic characteristics. J Micromech Microeng. 2010;20:95013–6.View ArticleGoogle Scholar
- Lamb H. On waves in an elastic plate. In: Proceedings of the Royal Society of London. 1917. p. 114–28.Google Scholar
- Lowe MJ. Matrix techniques for modeling ultrasonic waves in multilayered media. IEEE Trans Ultrason Ferroelectr Freq Control. 1995;42:525–42.View ArticleGoogle Scholar
- Zhang X, Qiang B, Hubmayr RD, Urban MW, Kinnick R, Greenleaf JF. Noninvasive ultrasound image guided surface wave method for measuring the wave speed and estimating the elasticity of lungs: a feasibility study. Ultrasonics. 2011;51:289–95.View ArticleGoogle Scholar
- Nenadic IZ, Urban MW, Aristizabal S, Mitchell SA, Humphrey TC, Greenleaf JF. On Lamb and Rayleigh wave convergence in viscoelastic tissues. Phys Med Biol. 2011;56:6723–38.View ArticleGoogle Scholar
- Zhang X, Osborn TG, Pittelkow MR, Qiang B, Kinnick RR, Greenleaf JF. Quantitative assessment of scleroderma by surface wave technique. Med Eng Phys. 2011;33:31–7.View ArticleGoogle Scholar
- Neu CP, Arastu HF, Curtiss S, Reddi AH. Characterization of engineered tissue construct mechanical function by magnetic resonance imaging. J Tissue Eng Regen Med. 2009;3:477–85.View ArticleGoogle Scholar
- Urban W, Chen S, Fatemi M. A review of shearwave dispersion ultrasound vibrometry (sduv) and its applications. Curr Med Imaging Rev. 2012;8:27–36.View ArticleGoogle Scholar