Research Article  Open Access
A Mathematical Model to Study the Mechanical Information Induced by LiftingThrusting Needle
Abstract
Focusing on the mechanical effect of traditional Chinese acupuncture, this study builds a mathematical model that simulates the mechanical process of liftingthrusting 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 threedimensional 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 liftingthrusting 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)
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 liftingthrusting 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 semiinfinite 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.
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 aswhere 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 straindisplacement relation and excluding the effect of volume force, we get the governing equation of motion in linear viscoelastic mediumwhere and are the viscoelastic material parameters and ρ is the density. Postulate the tissue is incompressible, that is, , and putting (1) into (2), we getRewrite it in tensor formSuppose the tissue around the needle is a semiinfinite space. This is an axisymmetric problem (,). Replace with and with ; then we have (5) in cylindrical coordinatewhere 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 toBased on the straindisplacement relationwhere is calculated by explicit difference schemeAnd the corresponding stress iswhere , , and is the phase difference between stress and strain. The shear stress amplitude .
In a cycle, the energy dissipation of oneunit volume caused by liftingthrusting is
The energy dissipation of the whole domain isDefine as the max energy storedTherefore, 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 iswhere r_{0} is the radius of the needle. Then
2.4. Analytic Solutions of Part (1)
Let ; (6a) is expressed asLet ; then we getEquation (17b) is a 1order Bessel equation, and the solutions arewhere J_{1}, Y_{1} are 1order Bessel function and Neumann function, respectively.
Inputting , we get . Because and , based on the boundary condition , we get ; thereforeBecause , then .
Put into (6b), Let ; then we getEquation (21b) is a 0order 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, ; thereforePutting (19) and (23) into (24), we getBased on the recurrence relation of Bessel’s functionwe 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 onedimensional question [15]; therefore (6a) and (6b) is changed toBased 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 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)
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_{0}). It shows that r_{0} has a minimal effect on U.
(a)
(b)
(a)
(b)
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)
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)
(a)
(b)
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.

(a)
(b)
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)
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_{0} 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_{0}) 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_{0} (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 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 liftingthrusting 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 quasilinear 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. View at: Google Scholar
 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,” EvidenceBased Complementary and Alternative Medicine, vol. 2013, Article ID 160651, 4 pages, 2013. View at: Google Scholar
 Y. L. Goh, C. E. Ho, and B. Zhao, “Acupuncture and depth: future direction for acupuncture research,” Evidencebased Complementary and Alternative Medicine, vol. 2014, Article ID 871217, 5 pages, 2014. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 H. M. Langevin, “The science of stretch,” The Scientist, vol. 27, no. 5, p. 32, 2013. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 A. Mvogo, G. H. BenBolie, and T. C. Kofané, “Energy transport in the three coupled α−polypeptide chains of collagen molecule with longrange interactions effect,” Chaos, Article ID 063115, pp. 1–9, 2015. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 S. Nicolle, J. Decorps, B. Fromy, and J. Palierne, “New regime in the mechanical behavior of skin: strainsoftening occurring before strainhardening,” Journal of the Mechanical Behavior of Biomedical Materials, vol. 69, pp. 98–106, 2017. View at: Publisher Site  Google Scholar
 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). View at: Google Scholar
 Y. C. Fung, Mechanical Properties of Living Tissues, Springer, Berlin, Germany, 1993.
 X.M. Wang, “Liftingthrusting and rotating manipulations :a comparison on energy input,” Chinese Acupuncture & Moxibustion, vol. 31, no. 1, pp. 71–74, 2011 (Chinese). View at: Google Scholar
 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. View at: Google Scholar
 G. H. Ding, X. Y. Shen, J. H. Dai et al., “The difference of acupuncture manipulation frequencies of liftingthrusting and twirlingrotating at deqi and nondeqi,” Chinese Acupuncture & Moxibustion, vol. 22, no. 10, pp. 679–681, 2002 (Chinese). View at: Google Scholar
 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). View at: Google Scholar
 Y. L. Ren, T. P. Guo, H. B. Du et al., “A survey of the practice and perspectives of chinese acupuncturists on deqi,” EvidenceBased Complementary and Alternative Medicine, vol. 2015, Article ID 684708, 8 pages, 2015. View at: Google Scholar
 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. View at: Google Scholar
 L. FallandCheung, 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
Copyright
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.