Evidence-Based Complementary and Alternative Medicine
Volume 2019 (2019), Article ID 5475426, 13 pages
https://doi.org/10.1155/2019/5475426

Research Article

A Mathematical Model to Study the Mechanical Information Induced by Lifting-Thrusting Needle

Yi Yu, Wei Yao, and Guanghong Ding

Department of Aeronautics and Astronautics, Fudan University, Shanghai Research Center of Acupuncture, 220 Handan Road, Shanghai, 200433, China

Correspondence should be addressed to Wei Yao; weiyao@fudan.edu.cn

Received 25 October 2018; Revised 7 January 2019; Accepted 11 March 2019; Published 8 April 2019

Academic Editor: Jeng-Ren Duann

Copyright © 2019 Yi Yu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Focusing on the mechanical effect of traditional Chinese acupuncture, this study builds a mathematical model that simulates the mechanical process of lifting-thrusting needle. Analytic and numerical solutions are obtained to explore the mechanical information (displacement, strain, stress, and energy) in the skin tissue. Our results show that needle manipulation leads to tissue displacement and mechanical stress field, but the needle should be inserted into the right position (about cm around the acupoint, where is the angular frequency) and enough depth (about 2 cm in lower limbs) to achieve effective mechanical stimuli; the tissue displacement decays with an increase of distance from the stimulus position, more rapidly at higher frequencies; there is an inverse relationship between the area of the ‘effective influence region’ (where shear strain > 0.2) and the stimulus frequency, which means larger needle movement is needed at higher frequencies to achieve a better curative effect; more energy is required to maintain high frequency manipulation. This study proposes a preliminary comprehension of the mechanical response around the needle during the acupuncture process.

1. Introduction

Acupuncture, a physiotherapy sourced from Traditional Chinese Medicine (TCM), now gains recognition all around the world. Its therapeutic effect has been widely accepted since the National Institutes of Health Hearing in 1997 [1]. Acupuncture is a biomechanical process, during which a needle is inserted into the skin at the acupoint, and mechanical stimulation to the tissue is generated when twisting or lifting the needle. Elucidating the mechanical and biological mechanism of acupuncture effect is one of the urgent problems in TCM as well as international researches on physiotherapy.

The premise to accomplish acupuncture effect is to get sensations (“De Qi” in Chinese), which means to maintain a certain intensity of stimulation by manipulating the needle. However, “De Qi” is usually described subjectively according to the experiences of skillful acupuncturists and the feelings of some sensitive patients; the underlying issues are lack of quantitative research. Therefore, it is important to make clear the mechanical sensation in acupuncture and to reveal the key factors related to it.

Numerous work and different explanations have been put forward to help understand the acupuncture process. Recent morphological studies have made clear that the acupoint is mainly composed of connective tissue. As shown in Figure 1, the dermic dense connective tissue and subcutaneous loose connective tissue form a three-dimensional collagen fiber network with numerous intertwining blood vessels, nerves, mast cells, and lymphatic vessels [2, 3]. Goh et al. hypothesized that there might be a relationship between the efficacy of acupuncture treatment and the depth of insertion needle [4]. Langevin et al. hypothesized that needle manipulation generates mechanical signals at acupoints [5], where longitudinal waves and transverse waves can propagate. Acupuncture needle manipulation results in sustained stretching by lifting-thrusting or rotating in the tissue of the acupoint [6]. Yu et al.’s experiment showed that the acupuncture effect and the force on the needle are very significant when the collagen structure is complete, and when the collagen is destroyed, the acupuncture effect is significantly weakened and the force on the needle body is very small [7]. Mvogo et al. constructed the analytic solutions of the tissue with collagen [8].

(a)

(b)

Figure 1: Acupoint morphology. (a) Basic pattern of organization of subcutaneous layers [9, 10]; (b) acupuncture location in lower limbs.

The mechanical effect of needle insertion and manipulation is the key to acupuncture; however, researches related with the quantitative mechanical effect of acupuncture are limited. Deleuze et al. has numerically simulated the features of subcutaneous interstitial flow induced by acupuncture [11]. Thiriet et al.’s mathematical model analyzed the different effects between acupuncture positioned at and out of acupoints [12]. We have also constructed a mathematical model to analyze the needle rotation at acupoints [13]. In this study, the previous work is extended, and the acupuncture of lifting-thrusting manipulation is studied through the establishment of a mechanical model based on the theory of viscoelastic mechanics. The influences of needle factors, such as stimulus frequency, movement amplitude, radius, and insertion depth, have been discussed in order to explore the acupuncture mechanical effect.

2. Methods

2.1. Model

Soft tissues are usually assumed as viscoelastic materials; some are more “elastic” and some are more “viscous”; they exhibit obvious hysteresis, stress relaxation, and other viscoelastic phenomena [14]. In our model, the skin tissue is treated as a semi-infinite viscoelastic body and the needle is treated as a rigid body. As shown in Figure 2, the tissue is divided into two parts. Part represents the tissue from the plane of the tip location to infinity ( ≥0) where tissue displacement is caused by the movement of the needle tip (point A). Part is the needle manipulation region (– ≤ < 0), where is the needle insertion depth and the standard is 2 cm in lower limbs.

Figure 2: Diagrammatic sketch of lifting-thrusting manipulation.

2.2. Equations

The soft tissue is supposed as a homogeneous isotropic incompressible viscoelastic body and fit for Feng's quasilinear viscoelastic theory of small amplitude displacement [15]. For acupuncture manipulation, the needle movement causes small amplitude oscillatory displacements in the low stress range (physiological state) of the tissue. The boundary conditions are specified as simple harmonic functions of time; therefore, the variables can be expressed as

where is the displacement component which is the distance from the original (equilibrium) position, is the coordinate component, is the time, is the angular frequency, is the component of stress tensor, and is the component of strain tensor.

Based on the strain-displacement relation and excluding the effect of volume force, we get the governing equation of motion in linear viscoelastic medium

where and are the viscoelastic material parameters and ρ is the density. Postulate the tissue is incompressible, that is, , and putting (1) into (2), we get

Rewrite it in tensor form

Suppose the tissue around the needle is a semi-infinite space. This is an axisymmetric problem (,). Replace with and with ; then we have (5) in cylindrical coordinate

where and are the displacement amplitude in r and z direction, respectively. Defining the dimensionless parameters: , where is the amplitude of needle movement. Then (5) is changed to

Based on the strain-displacement relation

where is calculated by explicit difference scheme

And the corresponding stress is

where , , and is the phase difference between stress and strain. The shear stress amplitude .

In a cycle, the energy dissipation of one-unit volume caused by lifting-thrusting is

The energy dissipation of the whole domain is

Define as the max energy stored

Therefore, the mechanical quality factor Q which reflect the ratio of energy storage to energy dissipation is

2.3. Boundary Conditions

The motion of the boundary is

where r₀ is the radius of the needle. Then

2.4. Analytic Solutions of Part (1)

Let ; (6a) is expressed as

Let ; then we get

Equation (17b) is a 1-order Bessel equation, and the solutions are

where J₁, Y₁ are 1-order Bessel function and Neumann function, respectively.

Inputting , we get . Because and , based on the boundary condition , we get ; therefore

Because , then .

Put into (6b),

Let ; then we get

Equation (21b) is a 0-order Bessel equation; therefore the solutions are

Inputting , we get . Because and , based on the boundary condition , we get ; therefore

Suppose the tissue is incompressible, that is, ; therefore

Putting (19) and (23) into (24), we get

Based on the recurrence relation of Bessel’s function

we get

Putting (27) into (25), we get

and then . Therefore

2.5. Calculation Solutions of Part (2)

Wang et al. have simplified the tissue as one-dimensional question [15]; therefore (6a) and (6b) is changed to

Based on , the solutions are

Both (29) and (31) do not fit very much in this part; therefore, we suppose (31) as the upper boundary condition (), (29) works as the lower boundary condition ( = 0), linear interpolation is adopted to calculate part , that is,

2.6. Parameters

The analytic solution (29) shows that λ affects the amplitude attenuation along . Langevin et al.’s research on acupuncture demonstrates the displacement is small when is over 2 cm [10]. Therefore, we suppose the displacement amplitude at = 2 cm is a small amount compared to the amplitude at = 0 cm (the plane of the needle tip), that is ; we can determine λ ≈ 0.25. The standard is 0.5 cm and is 0.02 cm.

The biological viscoelasticity is determined by [16]

where .

and will diverge at if have an imaginary part. Obviously, λ>0; notice that the imaginary part of is small compared to the real part; for example, = 1.2 − 0.03 when ω=2π and = 76 − 0.38 when ω=16π; therefore we replace with its module to void divergence.

3. Results

3.1. Boundary Solutions

Obviously, the upper () and lower ( = 0) boundary solutions determine displacement distributions in the skin tissue. The dimensionless tissue displacements ( at the upper and lower boundaries and the at the lower boundary) at two stimuli frequencies (ω = 2π and ω = 16π) are plotted in Figure 3. The distributions of are very similar at upper and lower boundaries.

(a)

(b)

Figure 3: Amplitudes of tissue displacements at the upper () and lower ( = 0) boundaries. (a) ω = 2π; (b) ω = 16π.

Figure 4 shows the corresponding absolute values of displacements ( and ) at lower boundary ( = 0). Additionally, the combined displacement is plotted. It shows that U is nearly the same as and is relatively small, especially at high stimulus frequency (ω=16π, in Figure 4(b)). Therefore, we will only study U afterwards. The influence of stimulus frequency on the displacement (U) is shown in Figure 5. It is clear that U decreases with the increase of R, and at higher frequencies, it decreases more rapidly.

(a)

(b)

Figure 4: Amplitudes of tissue displacements at = 0. (a) ω = 2π; (b) ω = 16π.

Figure 5: Amplitudes of tissue displacements under different frequencies.

The influence of needle movement amplitude () is shown in Figure 6. Figure 6(a) plots the absolute values of tissue displacement () at the lower boundary with different s (0.25, 0.5, 0.75, and 1.0cm). The range of is largely dependent on , with a nearly linear relationship at . Figure 6(b) shows the dimensionless displacement , which attenuates a little more slowly along R direction with an increasing . Figure 7 demonstrates in the cases of different needle radius (r₀). It shows that r₀ has a minimal effect on U.

(a)

(b)

Figure 6: Amplitudes of tissue displacements under different u₀s at ω = 2π; (a) . (b) .

(a)

(b)

Figure 7: Amplitudes of tissue displacements under different r₀s. (a) ω = 2π; (b) ω = 16π.

3.2. Displacement Amplitude (U) Distribution in the Skin Tissue

Figure 8 shows the displacement amplitude U distribution at ω = 2π and ω = 16π. Comparing Figures 8(a) and 8(b) shows that the tissue displacement around the needle () decays along the radial direction R under both low and high frequency stimulus, but U decays less rapidly under low frequency (U ≤ 50% where r > 0.6 cm) than that under high frequency (U ≤ 50% where r > 0.07 cm). Moreover, U also decays quickly under the needle tip ( > 0).

(a)

(b)

Figure 8: Amplitudes of tissue displacements; the black lines represent the displacement at , = 0 and . (a) ω = 2π; (b) ω = 16π.

3.3. Shear Strains () and Stresses () under Different Frequencies

Figure 9 shows the contour lines of shear strain amplitude in absolute values in the tissue. The needle is located at R = 0, (shown as white lines in Figure 9). The maximum shear strain occurs at the needle tip, and again, strains decay along the radial direction, with more rapid attenuation under high frequency stimulus. Thus, the low frequency stimulus affects larger area than high frequency stimulus. Take the area where > 0.2 as the effective influence region , and define Z and R at as and , respectively; then =×≈1×0.94 at ω = 2π and ≈ 1×0.08 at ω = 16π, shown as the black rectangles in Figure 9. It turns out that the stimulus frequency has little effect on but has significant influence on . Figure 10 shows the contour lines of dimensionless shear stress amplitude in absolute values (); it shows the same tendency as .

(a)

(b)

Figure 9: Contour lines of in the tissue; the white lines at represent the needle; the black rectangles represent the approximate area where > 0.2. (a) ω = 2π; (b) ω = 16π.

(a)

(b)

Figure 10: Contour lines of in the tissue, the white lines at represent the needle. (a) ω = 2π; (b) ω = 16π.

Table 1 lists s and s under different frequencies. Figure 11 shows that there is an approximately linear relationship between s and the reciprocal of the stimulus frequency 1/ωs. The same relationship also exists between s and 1/ωs.

Table 1: s and s under different angular frequencies.

ω
2π
4π
6π
8π
10π
12π
14π
16π

0.94
0.56
0.36
0.26
0.20
0.16
0.12
0.08

0.5×0.94
0.5×0.56
0.5×0.36
0.5×0.26
0.5×0.20
0.5×0.16
0.5×0.12
0.5×0.08

(a)

(b)

Figure 11: Relationship between s and s on one side and 1/ωs on the other. Diamonds represent calculating results and blue lines represent the fitting curve.

3.4. Energy Dissipation (ΔE) and Mechanical Quality Factor (Q) under Different Frequencies

Figure 12 shows ΔE and Q under different frequencies. Figure 12(a) shows that ΔE is larger and more intensified around the needle at high frequencies than that at low frequencies. For the mechanical quality factor Q, a nearly linear relationship between s and the stimulus frequencies can be found (shown in Figure 12(b)). This means that energy storage ratio increases with the increase of stimulus frequency, and the tissue demonstrates more “elastic” at high frequencies.

(a)

(b)

Figure 12: ΔE and Q under different frequencies. (a) Δ at = 0; (b) .

4. Discussion

Langevin et al. claimed that acupuncture manipulation would lead to tissue displacement and generate a mechanical stress field [10]. This phenomenon is well demonstrated in our simulation. Moreover, the local displacement, strain, stress, and energy near the acupoint (or the needle) have been quantitively analyzed. The characteristics of the tissue displacement are concluded as it is larger at low frequencies (see Figure 8) given other conditions the same (the needle movement u₀ and the material properties, etc.), and , the displacement will decay more slowly and extend relatively further from the needle at low frequencies (see Figure 5).

Main characteristics in the mechanical process of acupuncture are the stimulus amplitude (u₀) and the angular frequency (ω). However, there is a lack of consistency of these characteristics. Wang XM pointed out that mechanical information can only influence nearer region of the acupoint [17]. According to our simulation, the resultant tissue displacement is nearly proportional to u₀ (see Figure 6(a)). It is also shown that the displacement, strain, and stress under low frequency stimulus decay more slowly than those under high frequency (see Figures 8–10). These results indicate that the mechanical information (displacement, strain, and stress) can propagate further under low frequency stimulus. Figure 13 shows the characteristic curves of different acupuncturists, where significant differences in and ω can be found, but is larger at high frequencies; thus tends to be constant [18]. Figure 11(a) shows that there is nearly an inverse relationship between s and ωs, or , where is a constant. Given , then the radius of the effective influence area . Thus a constant is equivalent to a constant ; therefore, to get good acupuncture effect, a sufficient area of tissue should be deformed, which may be an alternative interpolation of acupuncture sensations in the mechanical viewpoint. However, our model is not able to figure out what “sufficient area” should be, which requires further experimental research. But suggestions on achieving mechanical sensations can be made. Since needle movement is a damage process to the tissue and slip between the needle movement and tissue movement may occur at large (we suppose there is no slip because of small ), thus is limited. Therefore, to achieve better effect (larger effective influence region), low frequency stimulus is a good choice. It was reported that the frequency of acupuncture sensations is concentrated at low frequency [18]. Ding et al. reported the mean frequency to achieve acupuncture sensations by experienced acupuncturists is 1.2 Hz [19].

Figure 13: Characteristic curves of five professional acupuncturists.

Figure 7 shows that needle diameter has no effect on displacement. Li et al. reported needling of different diameters induced the consistent change of fascial connective tissues and collagenous fiber arrangement [20]. Figures 9 and 10 show that strain and stress decay quickly under the needle tip. That means the needle must be inserted into the right depth (about 2 cm in lower limbs) and position ( cm around the acupoint in our simulation) to induce enough strain and stress. A survey showed that acupuncturists think deep insertion was easier to induce acupuncture sensations than shallow insertion [21].

Because of the limitation of experimental means, it is impossible to measure tissue strain and stress in vivo, tissue displacement, or the stress on the needle are exploited to infer the tissue strain [22]. But there is a complicated relationship between strain and displacement; it is controversial to infer stress field without considering needle frequency. Our model predicted the relationship between the main deformation areas and frequencies (see Figure 11). But there is no agreement on the threshold of to achieve acupuncture effect and tissue is more complicated (anisotropy, nonlinearity) than our model [23]; further experiments are required to testify the results.

Figure 12(a) shows that energy dissipation is relatively small under low frequency compared with that under high frequency; Figure 10(b) shows that there is a linear relationship between mechanical quality factors and frequencies. Obviously, a high energy is required to maintain a high frequency manipulation.

5. Conclusion

In this study, we have illustrated that lifting-thrusting needle will lead to tissue displacements and generate a mechanical stress field. The displacement correlated with needle frequency and decayed more slowly under low frequency than that under high frequency. Strain and stress correlated with frequency too, and there exists an inverse relationship between the main deformation areas and frequencies. Energy dissipation is also affected by frequency. This is a preliminary study of mechanical effects of acupuncture manipulation, isotropy, and quasi-linear viscoelastic theory adopted to simplify soft tissue, while Fox et al. observed anisotropic tissue motion induced by acupuncture [24]. Our results may reflect local effects around the needle; anisotropy must be considered if we study the mechanical information propagation in the tissue.

Data Availability

The data used to support the findings of this study are included within the article. We declared that materials described in the manuscript, including all relevant raw data, will be freely available to any scientist wishing to use them for noncommercial purposes, without breaching participant confidentiality.

Disclosure

This article does not contain any studies with human participants or animals performed by any of the authors.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was funded by National Natural Science Foundation of China (81574053 and 81473750) and Shanghai Key Laboratory of Acupuncture Mechanism and Acupoint Function (14DZ2260500).

References

NIH. Acupuncture. in Consensus Development Conference Statement. 1997, https://consensus.nih.gov/1997/1997acupuncture107html.htm:.

L. Fei, G. H. Ding, and E. Y. Chen, “Research progress in physical basis and functional characterization of meridians,” in WHO Report consultation meeting on traditional and modern medicine: harmonizing the two approaches, Manila, Philippines, 2000.

M. F. Luo, X. T. Dong, X. J. Somg, J. Jin, J. L. Zhann, and Y. Han, “Study on the dynamic compound structure composed of mast cells, blood vessels, and nerves in rat acupoint,” Evidence-Based Complementary and Alternative Medicine, vol. 2013, Article ID 160651, 4 pages, 2013.

Y. L. Goh, C. E. Ho, and B. Zhao, “Acupuncture and depth: future direction for acupuncture research,” Evidence-based Complementary and Alternative Medicine, vol. 2014, Article ID 871217, 5 pages, 2014.

H. M. Langevin, D. L. Churchill, and M. J. Cipolla, “Mechanical signaling through connective tissue: a mechanism for the therapeutic effect of acupuncture,” The FASEB Journal, vol. 15, no. 12, pp. 2275–2282, 2001.

H. M. Langevin, “The science of stretch,” The Scientist, vol. 27, no. 5, p. 32, 2013.

X. Yu, G. Ding, H. Huang, J. Lin, W. Yao, and R. Zhan, “Role of collagen fibers in acupuncture analgesia therapy on rats,” Connective Tissue Research, vol. 50, no. 2, pp. 110–120, 2009.

A. Mvogo, G. H. Ben-Bolie, and T. C. Kofané, “Energy transport in the three coupled α−polypeptide chains of collagen molecule with long-range interactions effect,” Chaos, Article ID 063115, pp. 1–9, 2015.

C. Stecco, C. Tiengo, and A. Stecco, “Fascia redefined: anatomical features and technical relevance in fascial flap surgery,” Surgical and Radiologic Anatomy, vol. 35, no. 5, pp. 369–376, 2013.

H. M. Langevin, E. E. Konofagou, G. J. Badger et al., “Tissue displacements during acupuncture using ultrasound elastography techniques,” Ultrasound in Medicine & Biology, vol. 30, no. 9, pp. 1173–1183, 2004.

Y. Deleuze, M. Thiriet, and T. W. Sheu, “Modeling and simulation of the interstitial medium deformation induced by the needle manipulation during acupuncture,” Communications in Computational Physics, vol. 18, no. 04, pp. 850–867, 2015.

M. Thiriet, Y. Deleuze, and T. W. Sheu, “A biological model of acupuncture and its derived mathematical modeling and simulations,” Communications in Computational Physics, vol. 18, no. 04, pp. 831–849, 2015.

W. Yao, Y. Yu, and G. Ding, “A hybrid method to study the mechanical information induced by needle rotating,” Mathematical Methods in the Applied Sciences, vol. 41, no. 15, pp. 5939–5950, 2018.

S. Nicolle, J. Decorps, B. Fromy, and J. Palierne, “New regime in the mechanical behavior of skin: strain-softening occurring before strain-hardening,” Journal of the Mechanical Behavior of Biomedical Materials, vol. 69, pp. 98–106, 2017.

L. Wang, G. H. Tao, M. D. Tao, and G. H. Ding, “The calculation of energy distribution based on two respective techniques of acupuncture therapy,” Journal of Medical Biomechanics, vol. 18, no. 4, pp. 195–201, 2003 (Chinese).

Y. C. Fung, Mechanical Properties of Living Tissues, Springer, Berlin, Germany, 1993.

X.-M. Wang, “Lifting-thrusting and rotating manipulations :a comparison on energy input,” Chinese Acupuncture & Moxibustion, vol. 31, no. 1, pp. 71–74, 2011 (Chinese).

T.-Y. Liu, L. Kuai, H.-Y. Yang, and M. Gao, “Preliminary research on the standardization of acupuncture manipulation,” Chinese Acupuncture & Moxibustion, vol. 28, no. 5, pp. 356–358, 2008.

G. H. Ding, X. Y. Shen, J. H. Dai et al., “The difference of acupuncture manipulation frequencies of lifting-thrusting and twirling-rotating at deqi and non-deqi,” Chinese Acupuncture & Moxibustion, vol. 22, no. 10, pp. 679–681, 2002 (Chinese).

C. Li, B. Chen, T. H. Hu, and L. Chen, “Effects of acupuncture with different filiform needles on tissues,” Chinese Acupuncture & Moxibustion, vol. 35, no. 8, pp. 801–805, 2015 (Chinese).

Y. L. Ren, T. P. Guo, H. B. Du et al., “A survey of the practice and perspectives of chinese acupuncturists on deqi,” Evidence-Based Complementary and Alternative Medicine, vol. 2015, Article ID 684708, 8 pages, 2015.

Y. S. Lee, W. M. Jung, I. S. Lee et al., “Visualizing motion patterns in aupuncture manipulation,” Journal of Visualized Experiments, vol. 113, Article ID e54213, pp. 1–7, 2016.

L. Falland-Cheung, M. Scholze, P. F. Lozano et al., “Mechanical properties of the human scalp in tension,” Journal of the Mechanical Behavior of Biomedical Materials, vol. 84, pp. 188–197, 2018.

J. R. Fox, W. Gray, C. Koptiuch, G. J. Badger, and H. M. Langevin, “Anisotropic tissue motion induced by acupuncture needling along Intermuscular connective tissue planes,” The Journal of Alternative and Complementary Medicine, vol. 20, no. 4, pp. 290–294, 2014.