Propagation of Acoustic Waves Caused by the Accelerations of Vibrating Hand-Held Tools in Viscoelastic Soft Tissues of Human Hands and a Mechanobiological Picture for the Related Injuries

As is well known, hand-arm vibration syndrome (HAVS), or vibration-induced white finger (VWF), which is a secondary form of Raynaud’s syndrome, is an industrial injury triggered by regular use of vibrating hand-held tools. According to the related biopsy tests, the main vibration-caused lesion is an increase in the thickness of the artery walls of the small arteries and arterioles resulted from enlarged vascular smooth muscle cells (VSMCs) in the wall layer known as tunica media. The present work develops a mechanobiological picture for the cell enlargement. The work deals with acoustic variables in solid materials, i.e., the non-equilibrium components of mechanical variables in the materials in the case where these components are weakly non-equilibrium. The work derives an explicit expression for the infinite-time cell-volume relative enlargement. This enlargement is directly affected by the acoustic pressure in the soft living tissue (SLT). In order to reduce the enlargement, one can reduce either the ratio of the acoustic pressure in the SLT to the cell bulk modulus or the relaxation time induced by the cell osmosis, or both the characteristics. Also, a mechanoprotective role of the above relaxation time in the cell-volume maintenance is noted. The above mechanobiological picture focuses attention on the pressure in an SLT and, thus, modeling of propagation of acoustic waves caused by the acceleration of a vibrating hand-held tool. The present work analyzes the propagation along the thickness of an infinite planar layer of an SLT. The work considers acoustic modeling. As a general viscoelastic acoustic model, the work suggests linear non-stationary partial How to cite this paper: Mamontov, E. and Berbyuk, V. (2017) Propagation of Acoustic Waves Caused by the Accelerations of Vibrating Hand-Held Tools in Viscoelastic Soft Tissues of Human Hands and a Mechanobiological Picture for the Related Injuries. Journal of Applied Mathematics and Physics, 5, 1997-2043. https://doi.org/10.4236/jamp.2017.510169 Received: July 7, 2017 Accepted: October 23, 2017 Published: October 26, 2017 Copyright © 2017 by authors and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/ Open Access E. Mamontov, V. Berbyuk DOI: 10.4236/jamp.2017.510169 1998 Journal of Applied Mathematics and Physics integro-differential equation (PIDE) for the weakly non-equilibrium component of the average normal stress (ANS) or, briefly, the acoustic ANS. The PIDE is, in the exponential approximation for the normalized stress-relaxation function (NSRF) reduced to the third-order linear non-stationary partial differential equation (PDE), which is of the Zener type. The unique advantage of the PIDE is that it presents a compact model for the acoustic ANS in an SLT, which explicitly includes the NSRF, thereby enabling a consistent description of the lossy-propagation effects inherent in SLTs. The one-spatial-coordinate version of this PDE in the planar SLT layer with the corresponding boundary conditions is considered. The relevance of these settings is motivated by a conclusion of other authors, which is based on the results of the frequency-domain simulation in three spatial coordinates. The boundary-value problem at arbitrary value of the stress-relaxation time (SRT) and arbitrary but sufficiently regular shape of the external acceleration is analytically solved by means of the Fourier method. The obtained solution is the steady-state acoustic ANS and allows calculation of the corresponding steady-state acoustic pressure as well. The derived analytical representations are computationally implemented. Propagation of the pressure waves in the SLT layer at zero and different nonzero values of the SRT, and the single-pulse external acceleration is presented. They complement the zero-SRT and zero-SRT-asymptote results with the results for various values of the SRT. The obtained pressure values are, at all of the space-time points under consideration, meeting the condition for the adequateness of the linear model. In the case where the SRT is zero, the results well agree with the ones obtained by using the simulation software package LS-DYNA. The dependence of the damping of acoustic variables in an SLT on the SRT in the present third-order case significantly generalizes the one in the second-order linear systems. The related resonance effect in the waves of the acoustic pressure propagating in an SLT is also discussed. The effects of the NSRF-originated memory function provided by the present third-order PDE model are necessary for proper simulation of the pressure, which is of special importance in the aforementioned mechanoboiological picture. The results obtained in the work present a viscoelastic acoustic framework for SLTs. These results open a way to quantitatively specific evaluation of technological strategies for reduction of the vibration-caused injuries or, loosely speaking, achieving “zero’’ injury.


Introduction
As is well known (e.g., [1]), vibration and compressive loads are mechanical stimuli that have a powerful influence on living tissues.Recent studies in the E. Mamontov 7 -9 in [1] and the papers in [2]).One of the most well-known negative implications is hand-arm vibration syndrome (HAVS), or vibration-induced white finger (VWF).It is a secondary form of Raynaud's syndrome, an industrial injury triggered by regular use of vibrating hand-held tools.
The vibrating tools generate the accelerations acting on, and perpendicular to the surface of human fingers and neighboring regions of palms.These accelerations create acoustic waves propagating in the human soft tissues and causing injuries.An acoustic variable in a solid material is the non-equilibrium component of a mechanical variable in the material in the case where this component is weakly non-equilibrium.
The Appendix develops a mechanobiological picture for the main vibration injuries.It indicates that there are two characteristics that directly affect the cell volume.One of them is the mechanical variable, which is the acoustic pressure in a SLT, whereas the other is the biophysical variable, which is the relaxation time induced by the cell osmosis.The main text of the present work focuses on modeling of the mentioned pressure.The purpose of the work is formulated in Subsection 1.2 on the basis of the discussion of the state of the art in the area in Subsection 1.1.Subsection 1.3 specifies the topics of the other sections and outlines the approach of the work.

State of the Art in Propagation of Acoustic Waves in Soft Solids
Propagation of acoustic waves in solids in the cases where the solids are elastic, i.e., inviscid, or, equivalently, conservative systems, is the well-studied topic.
There is a vast literature on it (e.g., [3] and the references therein) and a lot of corresponding software packages, commercial or research.Moreover, the topic is described and analyzed in many text-books.To mention a few, we note [4] [5] and [6].These materials enable one to solve virtually any elastic-solid acoustic problem.This is equally true with respect to soft solids as is explained in the remark below.Remark 1.1.The expressions for the elastic components of the entries of the Cauchy stress matrix in a solid material (e.g., [5], (1.3), (1.2)) in terms of the strains are known as the Hooke law (e.g., [5], (1.43)).These expressions are linear in the strains with the coefficients, which present the bulk and shear moduli of the material, K and G .The latter one is also known as the modulus of rigidity.These moduli are coupled with the Young modulus E and the Poisson coefficient ν as follows and determine the volume and shear viscosities η and µ , As follows from the experimental data in Table 1, shear moduli of SLTs are usually a few orders less than the corresponding bulk moduli.Table 1 specifies this feature in the case of SLTs in Column 3. Note that elastic properties and parameters of SLTs are discussed in considerable detail in the review work [9] (see also the references therein).
The mentioned feature enables one to neglect the shear-modulus-related terms in the aforementioned expressions for the elastic components of the entries of the Cauchy stress matrix, thereby reducing the stress matrix to the scalar matrix.The scalar in this matrix is P −∆ where P ∆ is the acoustic average normal stress (ANS).
Term P ∆ presents the sum of its elastic and viscous components.The elastic one is known as the acoustic pressure, p ∆ .Consequently, P ∆ coincides with p ∆ if and only if the material is inviscid, i.e., (see (1.2)), at 0 θ = .Thus, P p ∆ − ∆ is the viscous components of P ∆ .


The above discussion shows that acoustic problems in soft elastic solids can adequately be solved in terms of the well-established elastic-solid models (e.g., [3] and [9], and the references therein).This means that new research and development in acoustics of soft solids are associated with other, less studied problems.One of them is soft viscoelastic solids.An important example of viscoelastic solids is living tissues (e.g., p. R3 of [9], pp.194 and 196 of [3]), both SLTs and bones (e.g., [10]).
In the case of soft viscoelastic solids, the reduction to the shearless version (see Remark 1.1) is also applicable.However, key difficulties here are associated with allowing for non-nonzero values of bulk (or volume) viscosity η or, equivalently (see (1.2)), SRT θ .
The spatially heterogeneous viscoelastic model for solids, which is known since long ago, is the corresponding generalization (e.g., [5], (6.15), [11], Chapter 3 and Epilogue, [12], §34) of the spatially homogeneous Kelvin-Voigt (KV) model (e.g., [5], (6.8) or (6.9)).This generalization presents the system of linear nonstationary (LNS) partial differential equations (PDEs) for the entries of the three-dimensional displacement vector in a solid, and includes not only Table 1.Experimental data on the bulk and shear moduli of soft tissues [7].For comparison, one can obtain value for natural rubber from the Poisson coefficient reported in Table 1 in [8] and the first equality in (1.1).elastic moduli (1.1) but also viscosities (1.2).In particular, a version of this model is used as the linearized core ( [3], (5) or ( 6)) of the non-linear model in one spatial coordinate considered in [3].(Note that the equation of the form ( [3], (5)) was published by G.G. Stokes [13], and further analyzed in ( [11], Chapter 3 and Epilogue) and [14] long before work [3] was submitted.)Remark 1.2.According to theoretical physics (e.g., [15], Section 6 of Chapter II, [16], (8.6); see also the discussion in [17], Appendix), the acoustic stress in any viscoelastic solid is the sum of its elastic and viscous components where the viscous component is presented with the Boltzmann superposition integral.The integrand in this integral includes the normalized stress-relaxation function (NSRF).It is a function of the non-negative time separation.The integral of the NSRF in the time separation over the interval from zero to infinity is known as the SRT, θ . (The above term "normalized'' means that the NSRF at zero time separation is unit.) As one can show, the Boltzmann superposition integral, in the asymptotic limit case as 0 θ ↓ , is the product of θ and the time-derivative of the elastic component of the acoustic stress.This product presents the zero-SRT asymptote.Also, in the mentioned limit case, the expression for the total acoustic stress in a solid presents the KV model. One should note that the models based on the KV representation are, as models for the stress rather than the strain, applicable in the asymptotic limit indicated in Remark 1.2 only.The aforementioned work [3]  Along with this, the SRT values for the SLTs of the VWF-patient fingers and palms need not be that low.(These SLTs include the dermis, epidermis, subcutaneous tissue, and muscle tissue mentioned in Appendix.)For this very reason, the zero-SRT asymptote or the KV treatments are generally not suitable for the VWF problems.
However, to our knowledge, there is no any consistent acoustic model, which takes into account the NSRF at arbitrary values of the SRT rather than in the particular case of the zero-SRT asymptote only.

Beyond the State of the Art. Purpose of the Work
In order to fill the gap indicated in the last sentence of the previous subsection, paper [17] derives the general LNS partial integro-differential equation (PIDE) in the three spatial coordinates for the acoustic ANS (see [17], (2.10)).This equation takes into account the NSRF strictly in the way summarized in Remark 1.2.Moreover, in the NSRF exponential approximation determined by a single parameter, SRT, the derived equation is reduced to the third-order LNS PDE (see [17], (2.11)), which appears to be of the Zener type.
There are more complex approximations for the NSRF, for instance, the Prony series (e.g., [18]) and a number of other empirical and semi-empirical models (e.g., [19]).However, parameters of these models are determined from the NSRF measured values, which are difficult (if possible at all) to obtain by means of the patient-specific in vivo or ex vivo tests.In contrast to that, the single parameter of the exponential NSRF model, SRT θ , can be obtained from the measured values of one of the viscosities and the corresponding elastic modulus (see (1.2)).Another advantage of the exponential approximation is that it not only includes the basic parameter, SRT, but also allows to keep the number of the input parameters for the modeling at the minimum level.Also note that a consistent theoretical-physics model for the NSRF and the list of its parameters remains unknown.
The following four facts should also be noted.
Firstly, in a number of cases (e.g., [20], Figure 3), the accelerations generated by vibrating hand-held tools, can be outlined the following way.The acceleration is in the form of, loosely speaking, large periodic pulses (e.g., with the periods of 0.03 -0.02 s [20]).The pulse-peak values of these accelerations may reach ±10 5 m/s 2 .The corresponding displacements are, however, extremely small.All of the pulses are very similar to each other, roughly, the same.The shape of a pulse tends to zero as the time tends to infinity.The settling-to-zero time is usually noticeably less than the period.The corresponding acoustic signals in SLTs also well settle to zero within the period between the pulses.Consequently, it is usually sufficient to analyze the acoustic response to a single pulse only.
Secondly, works related to the areas described in Subsection 1.1 (e.g., [3] and [9] and the references therein) include numerous illustrations of computer-simulation results on propagation of acoustic waves in SLTs.They facilitate understanding of the acoustic phenomena.
Thirdly, a finger and a palm of the hand, which holds the handle of a vibrating tool, presents a complex geometrical system.Within a bounded area of the finger/palm-handle contact, the contact surface can be regarded as planar and the tool acceleration acting on the surface can be regarded as acting perpendicularly to the planar surface.These settings are adopted, for example, in paper [21] (where the planar contact region of the finger is termed the finger pad).The acoustic computer simulation of a finger model in three spatial dimensions presented in [21], which was preceded by the simulation in two spatial dimensions [22], comes to the following conclusion.The region of the soft tissue, which is the most prone to formation of vibration-induced injuries, is the tissue domain neighboring the contact, and this domain is the location of the greatest variations of acoustic variables (e.g., strains) during vibration exposure.
Consequently, analysis of the dominating acoustic phenomena in the finger or palm soft tissue can focus on a planar layer of the tissue between the contact plane and the finger/palm bones in the direction of the spatial coordinate perpendicular to the plane, i.e., along the thickness of the layer.One can add that, for simplicity, the planar layer can be regarded as infinite.
Fourthly, the computer simulations in the aforementioned papers [21] and [22] are based on nonlinear frequency-domain methods, which provide nonlinear response functions, and only deal with harmonic external time-varying signals.
In contrast to that, the approach developed in the present work is intended to deal not only with smooth signals (such as linear combinations of harmonics) but also the measured, highly irregular external time-varying signals.For this reason, it considers the time-domain models and methods only.
Taking into account the above four facts and the fact that the aforementioned third-order PDE open a way to consistent computer simulation of SLTs, the purpose of the present work is three-fold: 1) derivation of the answers to the three questions formulated in Appendix; 2) development of the Fourier-method steady-state solutions of the boundary-value problem for the third-order LNS PDE of the Zener type in an infinite planar layer of an SLT under the action of the acceleration produced by the vibrating hand-held tools; 3) simulation of propagation of the pressure waves in the SLT layer, which are caused by the single-pulse accelerations, at zero and different nonzero values of the SRT, and different settings for the acceleration pulse; the simulation is supposed to be implemented in two different ways: by means of the onedimensional version of the model in Point (1) and, at zero relaxation time, by means of LS-DYNA, the well-known purely numerical, finite-element simulator (e.g., [23]).
Achieving the above purpose will complement the aforementioned, zero-SRT and zero-SRT-asymptote results with the results for arbitrary values of the SRT.Task (1) is solved in Appendix.Tasks (2) and (3) are solved in the way described in the next subsection.

Topics of the Subsequent Sections and an Outline of the Present Approach
The work presents and discusses the aforementioned third-order PDE for the acoustic ANS in an SLT (Section 2).The corresponding boundary-value problem is formulated in the case of an infinite planar layer of an SLT (Sections 3 and 4).
This formulation includes the external acceleration of a vibrating hand-held tool.
The exact analytical expression for the acoustic ANS in the general case of the arbitrary but sufficiently regular shape of the acceleration is derived in the form of the Fourier-method series (Sections 5 and 6).In general, exact analytical expressions for solutions of models provide the most transparent views of the structure of the solutions.These forms include the explicit dependences of the solutions on the parameters and other input variables of the models.This is an advantage of analytical solutions.
Another advantage is that analytical solutions can provide the modeling terms, which are inaccessible in purely numerical approaches.For instance, the latter do not allow to determine the natural angular frequencies, which, in particular, can play a noticeable role in interpretation of the features of the related computational results (see (6.7) in connection with the text in Table 3).
The present work only considers the steady-state solutions of the aforementioned boundary-value problem.Loosely speaking, these are the ones independent of the initial values (e.g., see [24] for further details).Due to this feature, the steady-state solutions present the core behaviors of the systems, i.e., the behaviors, which are not affected by specific initial values or initial time points.
Also note that all of the physical-quantity values specified below, with the exception of the ones in Figures 1-12 are presented in the SI.

Acoustic Equation for the Scalar Stress in Soft Living Tissues
The model described below assumes that SLTs are isotropic, isothermal, spatially homogeneous at equilibrium, and, according to the SLT data in where 0 s > is the speed of the bulk acoustic waves, 2 ∇ is the Laplace            The first and second terms in the brackets in (2.1) correspond to the elastic and viscous stresses in the way discussed in Section 2 and Appendix of work [17].
Parameters s and K are coupled by means of relation (e.g., [17] ( ) Parameter θ is included in (1.2).Any of the relations (1.2) can be used for evaluation of θ .Two of many examples are indicated in the remarks below.Remark 2.1.In the case of liquid water, µ and G to the second equality in (1.2) results in the range of the values of SRT θ , which appear to be typical for SLTs, namely between Equation (2.1) was derived under the condition that (e.g., [27], (A.1.13)) ) assures that the SLT under consideration can be regarded as a linear material.Remark 2.2.As follows from equations (A.1.14)and (B.8) in [27], that where ( ) Paper [17] shows that PIDE (2.6) expresses the steady-state solution of LNS PDE (see [17], (2.11)) ( ) As is noted on p. 512 in [17], ODE (2.8) is the ANS-member of the family of spatially homogeneous versions of the Zener acoustic models (e.g., [5], (6.25)), which corresponds to the exponential approximation for the NSRF.Notably, this family is also used for experimental determination of SRTs in SLTs.For example, work [28] proposes a framework of quantitative assessment of SLTs that can be used in clinical diagnostics.The work applies the Zener models "(e.g., see (1) in [28]) and determines the tissue SRT in terms of the Young modulus and the corresponding viscosity (see [28], the right column on p. 154).Since the Young modulus E is coupled with both bulk and shear moduli K and G with the well-known linear relations (see the first relation in (1.1)), the SRT determined in the above way can be used in the capacity of SRT θ regarded in the present work (in particular, in expressions (1.2) for the corresponding viscosities).
Remark 2.3.In view of (2.3), the asymptotic representation of PIDE (2.1) in the limit case as ( ) Both (2.9) and (2.10) do not include the integral indicated in (2.1).This integral where the kernel function is the NSRF is often termed the memory function.It provides the integral influence of the values of a solution of (2.1) at all t t ψ − < upon the value at t.In this sense, the value at t "remembers'' the values at all t t ψ − < .This effect significantly complicates the behaviors of the solutions.The NSRF-originated memory function corresponds to the friction distributed in time.
In the case where the solution variable is the velocity potential of a fluid, the version of (2.9) was derived by G.G. Stokes [13] (see the discussion in [27], p.

Comparison of particular version (2.6) of Equation (2.1) with Equation (2.7)
shows that, in the latter equation, the terms, which include SRT θ , originate from the exponential-NSRF-originated memory function in (2.6). The advantage of PIDE (2.1) and its particular cases (2.6)  ) holds, presents an adequate result.This, so to say, self-testing allows to avoid a use of much more complex, nonlinear models (e.g., [3]) in the cases where linear models suffice.

Boundary Conditions for the Stress in the Case of an Infinite Planar Layer of a Soft Living Tissue
The present section considers the geometrically simplest case of the configuration of an SLT, namely where the SLT is an infinite planar layer.
Without a loss of generality, one can assume that the coordinate axis perpendicular to the layer is the x-axis and the two planar surfaces of the layer correspond to two different values of coordinate x, say, 0 x = and x h = where 0 h > is the thickness of the layer.Consequently, PDE (2.7) is reduced to ( ) One of the basic settings for this equation presumes that the external acceleration acting on the SLT layer is present at plane 0 x = , whereas the other plane, the one at x h = , contacts air.The corresponding boundary conditions are where ( ) a t is the acceleration and ( ) Note that term The output data of the mentioned problem are the versions of the two spacetime-dependent variables indicated in Remark 2.4 in one spatial coordinate, x. Remark 3.2.The above single-layer system can be generalized to a multi-layer system.In this case, one should include the boundary conditions at the inter-layer surfaces.Each of these conditions presents the equality of the acoustic ANS and the equality of the accelerations in the neighboring layers.
Multi-layer systems can substantially extend the scope of the present approach.
For example, in a three-layer system, one of the layers can be the muscle tissue of a palm, whereas the other two layers can be the skin and a layer of a protective material intended for a reduction of the amplitude or altering the frequency content of the acoustic waves penetrating into the soft tissues.Another example is a three-layer system where the layers represent the stratum corneum,

Transformation of the Boundary-Value Problem for the Stress to the One with Homogeneous Boundary Conditions
As is well known, the Fourier method is applicable to boundary-value problems if the boundary conditions are homogeneous.Problem (3.1)-(3.3)can be transformed into this form in the following way (e.g., [4], Section 4 in Chapter IX).
One introduces the change of variable where the physical dimension of variable ( ) , w x t is the squared velocity.This change allows to substitute (4.1) into (3.1)-(3.3)and thereby obtain the PDE and homogeneous boundary conditions for ( ) The last term on the right-hand side of (4.

Exact Analytical Expression for the Steady-State Stress
The present section considers solution of boundary-value problem (4.2)-(4.4)by means of the Fourier method (e.g., [4], Chapters VIII and IX).It is applicable to the PDEs, which are based on the Laplace differential operators.The Fourier method provides expansions of solutions of these PDEs in the space-dependent operator eigenfunctions with time-dependent coefficients (e.g., [29], Chapter V) and formulates the ODE systems for these coefficients.The main advantages of the Fourier method are the following.
• The above expansions enable the exact analytical solutions.
• The exact analytical expressions provide an explicit insight in the structure of the solution in the terms, which are also meaningful physically.• The method allows a variety of approximations of different accuracy and complexity.
The consideration in this section down (5.5) is based on the well-known ideas (e.g., [4] Since the source term in PDE (4.2) includes linear function h x − , it is necessary to involve the corresponding expansion for this function.By using (5.2), (5.3), and the well-known results (e.g., [29], Sections 21.4 and 22.3 of Chapter V, [30], 440.11), one can show that the expansion is Application of (5.1) and (5.4) to PDE (4.2) results in The third-order ODEs (5.5) can be represented in the form of the first-order ODE systems where vectors ( ) i F t and i φ , and matrixes i Γ are described with the following expressions (5.7) (5.9) ( ) , where I is the identity 3 3 where l e is the l th column of matrix I .Due to (5.8), relations (5.10) are reduced to expressions   The next section establishes asymptotic stability of systems (5.6) and derives the exact analytical expressions for the entries mentioned in Remark 5.1.

Expression for the Stress in Terms of the Input Data Only
The present section establishes asymptotic stability of the ODE systems for the time-coefficients in the Fourier representation for the scalar stress and derives expressions for the entries of the Cauchy matrixes of these systems.On the basis of this, the section obtains the exact analytical representation for the steady-state stress in terms of the input data.
In order to investigate the asymptotic stability of ODE systems (5.6), one needs to consider the real parts of the eigenvalues of matrixes (5.9) in (5.6).Due to the connection of the first-order systems (5.6) to the third-order scalar Equation (5.5), these eigenvalues are the roots λ of equations ( ) ( ) which are the characteristic equations for homogeneous versions of ODEs (5.5).Note that i κ in (6.1) are determined with (5.2).In view of the results developed in Theorem 2.1 of [17] on characteristic Equation (2.13), the following statements are valid.
(i) Each of the characteristic Equation ( 6.1) has one real root, 1 i θ − , and a pair of the complex conjugate roots 1 i i ιω − Θ ± , where ι is the imaginary unit, i.e., 2  1 ι = − .(iv) Asymptotic representation for 1 i Θ in the limit case as 0 (vi) Asymptotic representations for i ω in the limit cases as (vii) Each of ODEs (5.5) is asymptotically stable.
Statement (vii) implies that each of the ODE systems (5.6) is asymptotically stable and, thus, the analysis presented in Section 5 below (5.9) is relevant.As follows from statements (ii) and (v), the parameters i Θ and i ω of the complex conjugate root in statement (i) are functions of not only SRT θ and angular frequencies i sκ , which directly follow from the input data (see Remark 3.1, (2.2), and (5.2)), but also i θ , which is not quantified above.
As is indicated in statement (i), is the only real root of cubic Equation (6.1).It can be determined with the help of the well-known method of S. del Ferro and N. F. Tartaglia for solving cubic equations.Nowadays, it is known as the Cardano solution (e.g., [31], 1.8-3).One can show that, according to it, the mentioned root is determined as follows ( ) where and ( ) ( ) Continuous variable χ in (6.4) and (6.5) generalizes discrete variable 3).Numerical values of function (6.4) are exemplified in Table 2.
In view of (6.3), the equality in statement (ii) is equivalent to which, in turn, transforms the expression in statement (v) into Remark 6.1.As follows from statements (iii) and (iv), as well as expressions where i E , i S , and i C are the arbitrary parameters independent of time.This, in particular, means that .13 i t s Φ − are specified with the condition indicated on the right-hand side of (6.11).On the strength of (6.8), (6.12), and (6.13), this condition is equivalent to three equalities, which are the second ones in the following three relations These systems are the ones for determination of parameters i E , i S , and i C .


We also not the following feature.Remark 6.6.As follows from (6.16), (6.2), the properties of function ( ) M ⋅ in Remark 6.1, and statement (vi), the asymptotic representation for i D in the limit case as i → ∞ is ( ) As is well known (e.g., [30] which, after application of (6.14)-(6.17),become These expressions enable one to calculate the integrals in (5.16) as In view of (4.5), the physical dimension of each of the coefficients (6.19)-(6.21) is acceleration.Relation (6.18) transforms (5.16) into The derived analytical expression (6.22) explicitly shows the wave nature of the scalar stress propagating in layers of SLTs.The scalar stress presents the superposition of the space-or time-dependent sine and cosine waves with coefficients (6.19)-( 6 Remark 6.7.The numerator in each of the ratios in (6.22) comprises three additive terms.Expressions (6.19)-(6.21),as well as (6.3), (6.6), and (6.17) show the following.
The sum of the second and third of the three terms presents the harmonic-type oscillations with frequency i ω , time-dependent amplitude , , and a time-dependent phase around zero.These oscillations are present, no matter if SRT θ is positive (damped oscillations) or zero (undamped oscillations described with (2.10)).
In contrast to that, the first of the three terms, namely (6.19), is present only if SRT θ is positive.If the latter is zero, the term is also zero.This means that the first term is present if and only if the exponential-NSRF-originated memory function noted in Remark 2.3 is not neglected as, for instance, in the zero-SRT and zero-SRT-asymptote cases also noted in Remark 2.3.Notably, expression (6.19) indicates that this term is, in contrast to the above harmonic, of the non-oscillatory form.As follows from its role in the aforementioned numerator, it creates a time-dependent shift of the above harmonic along the P ∆ -axis, thereby producing the behavior and values of P ∆ unavailable in the harmonic-type oscillations. The effect of the exponential-NSRF-originated memory-function terms discussed in Remark 6.7 is quantitatively illustrated with computer-simulation results in Section 9.
Section 7 summarizes the exact analytical solution derived in the present section.
As follows from (6.22), the derived solution presents the superposition of the space-and time-dependent sine and cosine waves with the coefficients, which depend on the time and the external-acceleration function.Expression (6.22) is similar to the ones, which are well known in analysis of wave equations in acoustics (e.g., [4], Section 4 of Chapter IX).Thus, solution (6.22) explicitly shows the wave nature of the signals propagating in SLTs.
The above analytical expression provides the explicit dependences of the solution not only on the space-time point ( ) , x t but also on each of the input parameters and the external-acceleration function.

Expressions for the Steady-State Pressure. Convergences of the Corresponding Series
As follows from Remark 3.1, in addition to solution (6.22), there is one more space-time-dependent output variable: the acoustic pressure.
Substitution of (8.2) into (8.1)leads to In the limit case as  ( ) , of the series in the limit case as i → ∞ , is proportional to . This follows from the denominator in (6.22), expression (6.16) and (6.17), the form of sine and cosine coefficients (6.20) and (6.21), statement (vi) in Section 6, and Remark 6.6.
Since series ( ) converges (e.g., [30], 48.004), the series in (6.The value of θ is determined as described below. As follows from the discussion on pp.1961 and 1966 in [7], the tissues indicated in Table 1 are generally anisotropic materials and the ranges of the values in this table can be due to this feature.Since the heart tissues are muscle tissues, the lowest value of G for heart in Table 1, i.e.,   (6.22).The countably infinite sum in (6.22) is replaced with the finite sum where the number of the terms is on the order of a few tens or hundreds, depending on specific example.
For comparison of the quantitative results of the present modeling, the simulation with the multi-physics LS-DYNA software package is also carried out.
According to [33], the corresponding settings are the following.
The prediction of wave propagation in the SLT material is modeled with a two-dimensional plane strain finite-element model.The central-difference method is adopted.The numerical simulation model consists of a rectangular plate of the material with a height of 10 −2 , discretized with two-dimensional plane strain continuum elements.In order to make the simulation model comparable to the above one-dimensional simulation, a lateral displacement boundary condition is introduced, which reduces the model to a one-dimensional problem.In addition, a non-reflecting boundary condition is used at the lateral boundaries of the simulation model in order to prevent stress wave reflections, which otherwise would be generated at these boundaries.At the upper boundary, no boundary condition is applied, so that approaching pressure waves are completely reflected.The material-layer response is time and history dependent and is usually described by viscoelastic constitutive models, based on, in the simplest case, exponential stress relaxation functions.In the present one-dimensional version, shear relaxation phenomena are not present and such physical phenomena as strain rate dependence and bulk viscosity are not taken into account.Therefore, the material is described with a linear elastic isotropic material model with the aforementioned values of ρ and K (see the text 1) Comparison of Figure 1 and Figure 2, as well as Figure 11 and Figure 12 shows that, in the cases of Examples 1 and 3 and at zero SRT, the computational results of the present work well agree with the purely numerical results obtained with the help of the LS-DYNA simulation package.
2) The influence of different values of the SRT scaled with respect to value (9.1) is exemplified with the present results in Figures 3-10.
Comparison of Figure 3 and Figure 4 indicate that the SRT value * 0.001 is that small that the pressure is very weakly affected by the SRT.
1.000 θ θ = × , respectively, show that these values indicate that the pressures well "feel'' them.In Figures 5-7, one can even see how the pressures settle to zero.However, this settling cannot be seen in Figures 3-9 show that the damping of the pressure waves in an SLT is inversely proportional to SRT θ at not very small θ but is directly proportional to θ if θ is sufficiently small (cp., Figure 2 and Figure 3), and is not present at all if 0 θ = (see Figure 2).This picture is physically sound and remains valid not only for the pressure but also for other acoustic variables because it reflects the properties of the eigenvalues of matrixes (5.9) in (5.6) (see (6.6) and Remark 6.1).The mentioned dependence of the damping on SRT θ means that this dependence is nonlinear.
3) Concerning the high maximums of the pressures in Figures 2-9, one can note the following.
In the case of Figures 1-9, the external acceleration is (9.2) at As is well known (e.g., [4], (66) in Chapter IX), in the case where 0 θ = or, equivalently, PDE (2.7) is of the second rather than third order (e.g., see Figure 1), the Fourier series for the solution includes the resonance terms proportional to ( ) where ω is the angular frequency of the external signal (cp., (9.4)) and , 1, 2, i i ω = , are the natural angular frequencies of the system (cp., (6.7)).If there is such = , then the resonance occurs and the solution is infinite.None of these equalities takes place in the case corresponding to Figure 1, according to the present computational results.However, there are two natural angular frequencies, which are rather close to ω (see Row 1 of Table 3).Consequently, the corresponding versions of ( ) Note that, in the case of Example 3 in Figure 11 or Figure 12, the externalacceleration angular frequency is an order higher than the one in the case of Example 1 in Figure 1 or Figure 2. Consequently, the frequency in Example 3 is an order away from the principal frequency and the one closest to it in Row 1 of Table 3.This feature contributes to the fact that the pressure amplitude in Example 3 (see Figure 11 or Figure 12) is significantly lower than the one in Example 2 (see Figure 1 or Figure 2) due to the reduced resonance effect.
In the cases of Figures 2-9, SRT θ is positive rather than zero and, thus, PDE (2.7) is of the third order.In each of these cases, there are also two natural angular frequencies, which are rather close to ω (see Rows 2-8 of Table 3).
However, the non-zero SRT values preclude the solution to become infinite, and the third order of the equation makes the resonance-effect dependences more complex.When the pressure weakly "feels'' the SRT value, as is the cases of the behaviors in Figure 3 and Figure 4, the increase in the pressure amplitude is fairly small.When the pressure "feels'' the SRT value well, as is the cases of the behaviors in Figures 5-9, the increase in the pressure amplitude is quite pronounced.
Moreover, the fact that PDE (2.7) at positive SRT ( 0 θ > ) is of the third rather than second order creates the time behaviors, which do not exist in the second-order, i.e., zero-SRT ( 0 θ = ) case.Indeed, Remark 6.One should also note that the compressions and rarefactions in an SLT, which are associated with an oscillatory pressure, can cause lesions in it.The above shifts in the pressure oscillations lead to the pressure values that cannot be described in terms of familiar wave equations.Thus, the effects of the NSRF-originated memory function provided by the present, third-order PDE model are of importance for proper evaluation of the risks of formation of lesions in SLTs.
4) The aforementioned figures illustrate how the pressure damps to zero if the external acceleration damps to zero (see (9.2)).If the acceleration does not damp to zero, for example, is periodic at positive values of the time (see (9.3)), the pressure does not damp to zero either, and the shape of its damping is more complex.In contrast to that, work [22] assumes that the viscoelastic soft tissues, which are under consideration in it, are non-linearly elastic.The work calculates all of the entries of the strain matrix in the frequency domain.This matrix is used in the frequency-domain version of relation (2.5).However, work [22] does not check, e.g., in the way specified in Remark 2.2, if the results of [22] violate condition (2.4) for the linear-model adequacy at least at a single space-time point.This leaves the need in the elastic model of work [22] to be nonlinear beyond the scope.Consequently, the purpose of the nonlinear model underlying [22] remains elusive.As is well known, nonlinearity of models significantly complicates not only both qualitative and quantitative analyses but also a meaningful interpretation of their results. Remark 9.3.Work [22] dealing with a nonlinear acoustic model in the frequency domain notes a few peaks of the displacement in the finger soft tissue and interprets them as the resonance ones.However, the work does not include analysis sufficient for this interpretation.
Indeed, the nature of nonlinear systems is much more complex than the one of linear systems.In general, a peak in the frequency behavior of a scalar variable of a non-linear system need not be a manifestation of resonance.From the point of view of physics, resonances in both linear and nonlinear systems are defined by whether or not the external-signal frequency coincides with, or close to, one of the (generalized) natural frequencies of the system (e.g., [34], Section 29, [35]).
However, work [22] does not show that the above peaks result from the coincidence (or proximity) of the external-signal frequency and one of the generalized natural frequencies.In other words, no natural frequency is derived or at least approximately estimated.This is in a striking contrast with what is discussed in the above Point 3 in connection with the natural frequencies listed in Table 3.Thus, the resonance reading in [22] needs a substantial extra devel-

opment. 
One can also note that paper [22] formulates other conclusions, which are in fact at the hypothesis level.For example, this paper, by means of computer simulation, studies the behavior of the finger-tissue strain in the frequency interval 3 16, 2 10   ×   (that does not contain values 10 5 and 10 6 used in the above Examples 1 -3).The conclusion of [22] is that the strain is concentrated in the skin tissue (at the soft-tissue depth less than 10 −3 ) at all frequencies in the range ( ) 3 10 , ∞ (see [22], p. 726).The reason why this interval includes not only interval ( 3 3 10 , 2 10  ×  , which is a part of the above range , × ∞ , which is not considered in [22] at all, is unknown.
The results of the present works are summarized in the next section.

Conclusion
As is well known, hand-arm vibration syndrome (HAVS), or vibration-induced white finger (VWF), which is a secondary form of Raynaud's syndrome, is an industrial injury triggered by regular use of vibrating hand-held tools.According well agree with the ones obtained by numerical simulation with the help of the LS-DYNA software package.
• The influence of values of SRT θ is discussed in terms of the obtained quantitative results and graphic illustrations.It is shown that dependence of the damping of acoustic variables in an SLT on θ in the present, third-order case significantly generalizes the one in the second-order linear systems.This picture is consistent and physically sound.Also, the related resonance effect in the waves of the acoustic pressure propagating in an SLT is discussed.It is explained in terms of the angular frequency of the external acceleration and specific values of the angular natural frequencies.The latter are computed in the present approach but are difficult to be determined within the frame of purely numerical simulation software.• There are the solution components, which are new with respect to the zero-SRT and zero-SRT-asymptote cases and are due to the effects of the NSRF-originated memory-function.The role of the new components is discussed in connection with the shifts of the pressure plots along the pressure axis.The memory-function effects provided by the present, third-order PDE model are of importance for proper evaluation of acoustic variables in SLTs.
• It is found that, in general, a combination of the resonance effects and the effects of the above damping-only components complicate the space-time behaviors of the acoustic ANS and acoustic pressure.
The results mentioned in the above bullets, with the exception of the first one, are obtained in the exponential approximation for the NSRF.More general cases where the NSRF shape need not be exponential can also be targeted.However, it is explained that the results based on the exponential approximation are correct qualitatively and can be acceptable quantitatively.
The outcomes of the main part of the work, in conjunction with the mechanobiological picture developed in Appendix, present a viscoelastic acoustic framework for SLTs.This framework opens a way to quantitatively specific evaluation of technological strategies for reduction of the vibration-caused injuries or, loosely speaking, for achieving "zero'' injury.This evaluation is suggested as the topic for the nearest future research.

Appendix. Main Vibration Injuries in an SLT, the
Corresponding Mechanobiological Picture, and the Role of a Time-Varying Acoustic Pressure The VWF injuries or, in a more professional language, lesions in human fingers are well known.They are described in work [36] in terms of the data obtained by means of biopsy tests.This work notes ( [36], p. 280) that the main lesion is an increase in the thickness of the artery walls of the small arteries and arterioles resulted from muscular hypertrophy.As is shown in Figure 13, the muscular content of the artery wall is the wall layer known as tunica media and formed by vascular smooth muscle cells (VSMCs).More specifically ( [36], pp.280-281), the thickening of the tunica media is identified by the presence of hypertrophied nuclei and an increased cytoplasmic bulk of the cells, i.e., enlargement of VSMCs.
Paper [37], which is also based on the biopsy-test data, studies the enlargement of VSMCs in the small arteries as well.It deals with the arteries in the inner layer of the skin, i.e., the layer of subcutaneous fat.As follows from the measurement data in ([37], Table 2), the enlarged VSMC presents a significant lesion when the volume of an enlarged cell 1.5 times greater than the volume of a normal cell.
The mentioned results of biopsy tests mean that populations of the enlarged VSMCs present long-lasting lesions.To the best of our knowledge, the literature does not provide facts on that these lesions, in the course of time, return to normal states (i.e., the enlarged-cell volumes return to their normal values).In The cell enlargement develops approximately at the same proportion along each of the three spatial axes (e.g., see the experimental data on the VSMC enlargement in ([37], Table 2)).This means that ( ) ( ) or, in terms of (A.4), ( ) Example A.1.According to experimental study in [37], VSMCs have the shapes that can approximately be described as elongated cylinders (e.g., see [37], Table 2).In the case of a normal VSMC, the cylinder cell volume is Parameter c H , the hydraulic conductance of the VSMC membrane, can be taken from [44].Indeed, according to Fig. 2(b) in [44] on the basal-cell case, We assume that the well-known typical value of c c n σ for an intracellular fluid (e.g., [42], the last column in Table 3 and Section Substitution of (A.33) into (A.32)results in (A.28).This proves assertion (i).
In order to prove assertion (iii), it is sufficient to note that application of assumption (A.26) to inequalities (A.30) assures relation (A.3), which, as is explained in the text above equation (A.23), makes this equation valid.
Passing to the limit as t tends to infinity in (A.28) and taking into assumption the stress-relaxation time (SRT).

The above results are summarized in Section 7 .
The analytical representation for the acoustic pressure in an SLT and the convergence of the series are considered in Section 8.The developed analytical representations are computationally implemented in computer simulation.The corresponding simulation results along the thickness of the muscle-tissue layer of a palm are exemplified in Section 9. Section 10 concludes the work.It can also serve as an executive summary.The role of Appendix is mentioned in Subsection 1.2 and specified at the beginning and in E. Mamontov, V. Berbyuk DOI: 10.4236/jamp.2017.5101692004 Journal of Applied Mathematics and Physics the first bullet of Section 10.

Figure 1 .
Figure 1.Propagation of the pressure waves in the 10-mm layer of the muscle tissue in the case of Example 1 at zero value of SRT θ.The pressure unit is MPa.The data were obtained with the help of the LS-DYNA software package [33].

Figure 2 .
Figure 2. Propagation of the pressure waves in the 10-mm layer of the muscle tissue in the case of Example 1 at zero value of SRT θ.The data were obtained in the present treatment.

Figure 3 .
Figure 3. Propagation of the pressure wave in the 10-mm layer of the muscle tissue in the case of Example 1 at value 0.001 × θ * = 0.95 × 10 −7 of SRT θ (see (9.1) for θ * ).The data were obtained in the present treatment.

Figure 4 .
Figure 4. Propagation of the pressure waves in the 10-mm layer of the muscle tissue in the case of Example 1 at value 0.003 × θ * = 0.285 × 10 −6 of SRT θ (see (9.1) for θ * ).The data were obtained in the present treatment.

Figure 5 .
Figure 5. Propagation of the pressure waves in the 10-mm layer of the muscle tissue in the case of Example 1 at value 0.010 × θ * = 0.95 × 10 −6 of SRT θ (see (9.1) for θ * ).The data were obtained in the present treatment.

Figure 6 .
Figure 6.Propagation of the pressure waves in the 10-mm layer of the muscle tissue in the case of Example 1 at value 0.030 × θ * = 0.285 × 10 −5 of SRT θ (see (9.1) for θ * ).The data were obtained in the present treatment.

Figure 7 .
Figure 7. Propagation of the pressure waves in the 10-mm layer of the muscle tissue in the case of Example 1 at value 0.100 × θ * = 0.95 × 10 −5 of SRT theta (see (9.1) for θ * ).The data were obtained in the present treatment.

Figure 8 .
Figure 8. Propagation of the pressure waves in the 10-mm layer of the muscle tissue in the case of Example 1 at value 0.300 × θ * = 0.285 × 10 −4 of SRT θ (see (9.1) for θ * ).The data were obtained in the present treatment.

Figure 9 .
Figure 9. Propagation of the pressure waves in the 10-mm layer of the muscle tissue in the case of Example 1 at value 1.000 × θ * = 0.95 × 10 −4 of SRT theta (see (9.1) for θ * ).The data were obtained in the present treatment.

Figure 10 .
Figure 10.Propagation of the pressure waves in the 10-mm layer of the muscle tissue in the case of Example 2 at value 0.030 × θ * = 0.285 × 10 −5 of SRT theta (see (9.1) for θ * ).The data were obtained in the present treatment.

Figure 11 .
Figure 11.Propagation of the pressure waves in the 10-mm layer of the muscle tissue in the case of Example 3 at zero value of SRT θ.The pressure unit is kPa.The data were obtained with the help of the LS-DYNA software package[33].

Figure 12 .
Figure 12.Propagation of the pressure waves in the 10-mm layer of the muscle tissue in the case of Example 3 at zero value of SRT θ.The data were obtained in the present treatment.
presents the volumetric density of the force acting at space-time point ( ) , x t .Remark 3.1.The input data for boundary-value problem (3.1)-(3.3)are the three parameters indicated in Remark 2.4, as well as parameter h and function ( ) epidermis, dermis, and subcutaneous tissue of a finger.The list of the related examples can easily be continued. Individual solutions of boundary-value problem (3.1)-(3.3)are specified with the initial conditions and can be obtained by the Fourier method.
2) presents the source term.It results from change of variable (4.1).In connection with the source terms, we accept the following.and, thus (see (4.5)), function ( ) a t are piecewise differentiable in the entire time axis. Assumption 4.1 is used in the consideration below.

. 5 )
Relations (5.5) present the countably infinite set of linear non-autonomousODEs of the third order for time-dependent coefficients ( ) i f t in the Fourier expansion(5.1).The forms of these ODEs show that they can be mutually independent.
, and, thus (see (4.5)), function ( ) a t is uniformly bounded in the entire time axis. Taking into account the exponential tendency to zero noted in the text on (5.12) as well as Assumptions 4.1 and 5.1, one can show that all of the coefficients (5.1) are uniformly bounded in the entire time axis.Combination of (4.1), (5.1), (5.3), and (5.15) results in ( ) zero to infinity.
.21), which depend on time and the external-acceleration function with coefficients.A special advantage of expressions (6.22) and (6.19)-(6.21) is that they deal with the general case of the involved variables and parameters rather than selected particular examples.The general case enables any number of particular examples, which can present individual settings and be quantified by computing, thereby providing the input data for the corresponding visualization for illustrative purposes.Other related output acoustic variables (such as the pressure) have similar properties.The expression for the pressure based on solution (6.22) is discussed in Section 8, which also explains the convergences of the series in (6.22).

 9 .
Remark 8.1 shows that representation (6.22) is properly defined.Examples of the Computer Simulation Results along the Thickness of an Infinite Planar Layer of the Muscle Tissue of a Palm This section reports examples of computer simulation results obtained by means of computational implementation of the present modeling.The examples are based on the following values of the input parameters (see Remark 3.1) h , ρ , K , and θ :

Figures 1 - 12 .
Figures 1-12.The examples also differ with different external-acceleration functions of the following two forms

3 .
In all of the examples, the pressures in Figures 2-10 and Figure 12 are calculated according to (8.3) where term

above ( 9
.1)).The tissue layer is accelerated at the bottom boundary with a single sinusoidal pulse (9.1).The pressure is evaluated at five discrete spatial points throughout along the height of the tissue layer to capture the transient propagation and subsequent reflection of the pressure wave at the top and E. Mamontov, V. Berbyuk DOI: 10.4236/jamp.2017.5101692026 Journal of Applied Mathematics and Physics bottom.Analysis of Figures 1-12 enables one to indicate a few features.

Figure 8 and
Figure 8 and Figure 9 because it takes place at the time values beyond the interval shown in the figures due to rather high SRTs, the resonance-related increase in the amplitude of the solution.E. Mamontov, V. Berbyuk DOI: 10.4236/jamp.2017.5101692027 Journal of Applied Mathematics and Physics 7 emphasizes the new components of the solution, which are due to the exponential-NSRForiginated memory function in (2.6) and time-dependently shift the oscillations along the P ∆ -axis.These shifts, in view of (8.3), provide the time-dependent shifts along the p ∆ -axis as well.Since the single pulse in Example 1 (see (9.2)) and each pulse in the pulse sequence in Example 2 is positive, the shifts are in the positive direction of the p ∆ -axis (e.g., they, in particular, increase positive

Figure 10
Figure 10 corresponds to the SRT value of to the related biopsy tests (see the corresponding references in Appendix), the main vibration-caused lesion is an increase in the thickness of the artery walls of the small arteries and arterioles resulted from enlarged vascular smooth muscle cells (VSMCs) in the wall layer known as tunica media.The threefold purpose of the work formulated in Points (1)-(3) in Section 1.2 is achieved.More specifically, the outcomes of the work include the following key components.•The present work develops a mechanobiological picture for the cell enlargement.The work deals with acoustic variables in solid materials, i.e., the non-equilibrium components of mechanical variables in the materials in the case where these components are weakly non-equilibrium.The work derives an explicit expression for the infinite-time cell-volume relative enlargement.This enlargement is directly affected by the acoustic pressure in the SLT.In order to reduce the enlargement, one can reduce either the ratio of the acoustic pressure in the SLT to the cell bulk modulus or the relaxation time induced by the cell osmosis, or both the characteristics.Also, the mechanobiological picture shows a mechanoprotective role of the above relaxation time in the cell-volume maintenance.The mentioned results focus attention on the acoustic pressure and, thus, modeling of propagation of acoustic waves caused by the acceleration of a vibrating hand-held tool.This topic is addressed by the results listed below.• The present work analyzes the propagation along the thickness of an infinite planar layer of an SLT.The work considers the modeling in the time domain (rather than the frequency domain) only.As a general viscoelastic acoustic model, the work suggests linear non-stationary partial integro-differential equation (PIDE) (2.1) for the weakly non-equilibrium component of the average normal stress (ANS) or, briefly, the acoustic ANS that was derived in [17].This equation explicitly includes the normalized stress-relaxation function (NSRF).The unique advantage of Equation (2.1) is that it presents a compact model for the acoustic ANS in an SLT, which explicitly includes the NSRF, thereby enabling a consistent description of the lossy-propagation effects inherent in SLTs.A more simple equation, the third-order linear non-stationary partial differential equation (PDE) (2.7) of the Zener type, has the same advantage, however, only in the case where the NSRF is of the exponential type.• The one-spatial-coordinate version (3.1) of PDE (2.7) in the SLT layer with boundary conditions (3.2) and (3.3) is considered.The relevance of these settings is motivated by a conclusion of other authors on the results of the frequency-domain simulation in three spatial coordinates.The boundary-value problem (3.1)-(3.3)at arbitrary value of SRT 0 θ ≥ and arbitrary but sufficiently regular shape of the external acceleration is analytically solved by means of the Fourier method.The obtained solution is the steady-state acoustic ANS.It allows calculation of the corresponding steady-state acoustic pressure with (8.3).• The derived analytical representations are computationally implemented in simulation software.Propagation of the pressure waves in the SLT layer at zero and different nonzero values of the SRT and the single-pulse external acceleration is illustrated by the corresponding simulation results.They complement the zero-SRT and zero-SRT-asymptote results with the results for various values of the SRT.• The obtained pressure values are, at all of the space-time points under consideration, meeting the condition for the adequateness of the linearity of the applied model.• It is shown that, in the case where 0 θ = , the obtained quantitative results
time induced by the cell osmosis.We consider ODE (A.13) in a semi-infinite time interval o t t > under the initial condition that the cell volume is at equilibrium at point o t t = , i.e. (see (A.4)), normal human temperature of about 37˚C,

.
3.4.3) is applicable to the case of a VSMC.This value is 281This value can be used as a preliminary estimation. Remark A.1.Multiplying Equation (A.13) considered in interval

τ
and passing in the resulting equation to the limit as0 c τ ↓ , one obtains that 0 c ε = at all o t t> .This result and (A.
, V. Berbyuk DOI: 10.4236/jamp.2017.5101691999 Journal of Applied Mathematics and Physics animal models demonstrate that certain types of mechanical loads regulate various tissues.However, "over exposure'' to mechanical stimuli is damaging for the tissues.It can cause symptoms, syndromes, and even diseases (e.g., see Refs.

Table 1 ,
shearless.In this case, one can use the LNS PIDE for acoustic ANS P [17]uation (2.7) is formally a PDE of the Zener type ([5], (6.40)) but, with respect to the origin and meaning, differs from the Zener PDEs noticeably (see the corresponding discussion in[17]).The fact that PIDE (2.6) is nothing but the aforementioned steady-state case of PDE(2.7)indicates that (2.7) is a more general description than (2.6).

Table 2 .
Examples of the numerical values of function )} 1)-(3.3) describes the acoustic ANS propagating in an SLT layer due to the external time-dependent acceleration.The input data for this problem are listed in Remark 3.1.

Table 3 .
The pairs of the angular natural frequencies of the solution of the model in the case of Example 1, which are closest to the external-acceleration angular frequency(9.4).The values of the SRT correspond to the ones for Figures1-12.These values present the fractions of the SRT value * θ (see (9.1)), which is discussed at the beginning of Section 9.The fractions are indicated in the left column.The data in the middle and right columns are obtained in the present approach.
SRT θ 15)lead to property Calculating the integral, which is the second multiplier in (A.32), by parts, one