^{1}

^{2}

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 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.

As is well known (e.g., [

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.

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., [

Remark 1.1. The expressions for the elastic components of the entries of the Cauchy stress matrix in a solid material (e.g., [

E = 3 K ( 1 − 2 ν ) = 2 G ( 1 + ν ) , ν = [ 3 − 2 ( G / K ) ] / [ 6 + 2 ( G / K ) ] , (1.1)

and determine the volume and shear viscosities η and μ ,

η = K θ , μ = G θ (1.2)

where θ ≥ 0 is the stress-relaxation time (SRT).

As follows from the experimental data in

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., [

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., [

Soft tissue | Bulk modulus, K [7, | Shear modulus, G [7, | Ratio of the shear modulus to the bulk one, G/K -- |
---|---|---|---|

1 | 2 | 3 | |

Stomach | 480 | 8 - 45 | ( 0.017 - 0.094 ) × 10 − 3 |

Liver | 280 | 37 - 340 | ( 0.132 - 1.214 ) × 10 − 3 |

Heart | 490 | 60 - 48 | ( 0.122 - 0.302 ) × 10 − 3 |

Lung | 150 | 10 - 54 | ( 0.067 - 0.36 ) × 10 − 3 |

elastic moduli (1.1) but also viscosities (1.2). In particular, a version of this model is used as the linearized core ( [

Remark 1.2. According to theoretical physics (e.g., [

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 [

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.

In order to fill the gap indicated in the last sentence of the previous subsection, paper [

There are more complex approximations for the NSRF, for instance, the Prony series (e.g., [

The following four facts should also be noted.

Firstly, in a number of cases (e.g., [^{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., [

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 [

Fourthly, the computer simulations in the aforementioned papers [

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 one- dimensional 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., [

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.

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).

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 the first bullet of Section 10.

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

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 [

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.

The model described below assumes that SLTs are isotropic, isothermal, spatially homogeneous at equilibrium, and, according to the SLT data in

∂ 2 Δ P ∂ t 2 = s 2 ∇ 2 [ Δ P + ∫ 0 ∞ Ψ ( ψ ) ∂ Δ P ( t − ψ ) ∂ t ∂ ψ ] , t ∈ ℝ , (2.1)

where s > 0 is the speed of the bulk acoustic waves, ∇ 2 is the Laplace

differential expression with respect to the three spatial coordinates, ψ ≥ 0 is the time separation, and Ψ ( ψ ) is the NSRF. The term “normalized’’ denoted with the letter “N’’ in the abbreviation “NSRF’’ (introduced in Subsection 1.2) means that Ψ ( 0 ) = 1 .

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 [

s = K / ρ (2.2)

where ρ > 0 is the volumetric mass density of an SLT. The integral in (2.1) is known as the Boltzmann superposition integral. Due to the presence of this integral in Equation (2.1), the equation is suitable for viscoelastic rather than elastic-only materials.

The basic parameter of NSRF Ψ ( ψ ) is SRT θ (already mentioned in Remark 1.1), which is determined with the well-known relation (e.g., [

θ = ∫ 0 ∞ Ψ ( ψ ) d ψ (2.3)

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, η = 3.09 × 10 − 3 and K = 2.2 × 10 9 (e.g., [

A number of works on non-invasive measurements in SLTs available in the literature indicate that in various SLTs, μ = 0.1 - 20 . As follows from Column 2 of ^{3}. Application of the values of μ 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 0.3 × 10 − 6 and 2.5 × 10 − 3 .

Equation (2.1) was derived under the condition that (e.g., [

| Δ p | / K ≪ 1 (2.4)

(see Remark 1.1 for Δ p ). Condition (2.4) 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 [

1 K ∂ Δ p ( x , t ) ∂ t = − tr { ∂ ε ( x , t ) ∂ t [ I + ε ( x , t ) ] − 1 } (2.5)

where ε ( x , t ) is the strain matrix, I is the 3 × 3 identity matrix, and tr ( ⋅ ) is the trace of a matrix. Equality (2.5) enables one to formulate condition (2.4) in terms of the strain matrix.

If one prefers the frequency-domain paradigm rather than the time-domain one, one can use a consistently derived frequency-domain counterpart of relation (2.5) in order to check condition (2.4). ,

In the simplest, exponential approximation for NSRF Ψ ( ψ ) , i.e., under the assumption that Ψ ( ψ ) = e x p ( − ψ / θ ) at all ψ ≥ 0 , Equation (2.1) is specified to (see (2.12) in [

Paper [

θ ∂ 3 Δ P / ∂ t 3 + ∂ 2 Δ P / ∂ t 2 = s 2 ∇ 2 ( Δ P + 2 θ ∂ Δ P / ∂ t ) (2.7)

regarded as an ordinary differential equation (ODE) of the first order for variable ∂ 2 Δ P / ∂ t 2 − s 2 ∇ 2 Δ P . Equation (2.7) is formally a PDE of the Zener type ( [

Each of the PIDEs (2.1) and (2.6), and, thus, PDE (2.7) is valid only if condition (2.4) holds. In order to verify it, Δ p should be available. Work [

θ ∂ Δ P / ∂ t + Δ P = Δ p + 2 θ ∂ Δ p / ∂ t (2.8)

As is noted on p. 512 in [

Remark 2.3. In view of (2.3), the asymptotic representation of PIDE (2.1) in the limit case as θ ↓ 0 is PDE (2.1) in [

∂ 2 Δ P / ∂ t 2 = s 2 ∇ 2 ( Δ P + θ ∂ Δ P / ∂ t ) (2.9)

that, at θ = 0 , is reduced to

∂ 2 Δ P / ∂ t 2 = s 2 ∇ 2 Δ P (2.10)

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.

Equation (2.9) was derived in [

Equation (2.10) presents the zero-SRT case of (2.1) and corresponds to a purely elastic (or inviscid) material.

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) or (2.7) is that any of them eliminates the zero-SRT asymptote limitation of the Stokes-type equation indicated in Remark 2.3.

Remark 2.4. The input data for PDE (2.7) are the following three parameters: ρ , K , and θ . Number s used in (2.7) is determined in terms of ρ and K by means of (2.2).

The output data of LNS-PDE (2.7) include (but need not be limited to) the following two space-time-dependent variables of an SLT, acoustic ANS Δ P and acoustic pressure Δ p (see (2.8)). ,

Expression (2.8) not only describes pressure Δ p in terms of the Δ P but also enables one to answer the question on if the results of the above linear models are acceptable. For example, any solution of PDE (2.7), which, by means of (2.8), provides the values of Δ p such that (2.4) holds, presents an adequate result. This, so to say, self-testing allows to avoid a use of much more complex, nonlinear models (e.g., [

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, x = 0 and x = h where h > 0 is the thickness of the layer. Consequently, PDE (2.7) is reduced to

θ ∂ 3 Δ P / ∂ t 3 + ∂ 2 Δ P / ∂ t 2 = s 2 ∂ 2 ( Δ P + 2 θ ∂ Δ P / ∂ t ) / ∂ x 2 , 0 ≤ x ≤ h . (3.1)

One of the basic settings for this equation presumes that the external acceleration acting on the SLT layer is present at plane x = 0 , whereas the other plane, the one at x = h , contacts air. The corresponding boundary conditions are

a ( t ) = − ρ − 1 ∂ Δ P ( 0 , t ) / ∂ x (3.2)

where a ( t ) is the acceleration and

Δ P ( h , t ) = 0 . (3.3)

Note that term − ∂ Δ P ( x , t ) / ∂ x 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 a ( ⋅ ) .

The output data of the mentioned problem are the versions of the two space- time-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, 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.

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., [

One introduces the change of variable

Δ P ( x , t ) = ρ [ w ( x , t ) − a ( t ) ( x − h ) ] (4.1)

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 w ( x , t ) ,

θ ∂ 3 w / ∂ t 3 + ∂ 2 w / ∂ t 2 = s 2 ∂ 2 ( w + 2 θ ∂ w / ∂ t ) / ∂ x 2 − a ¯ ( t ) ( h − x ) , 0 ≤ x ≤ h (4.2)

∂ w ( 0 , t ) / ∂ x = 0 (4.3)

w ( h , t ) = 0 (4.4)

where

a ¯ ( t ) = θ d 3 a ( t ) / d t 3 + d 2 a ( t ) / d t 2 (4.5)

The last term on the right-hand side of (4.2) presents the source term. It results from change of variable (4.1). In connection with the source terms, we accept the following.

Assumption 4.1. Function d 3 a ( t ) / d t 3 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.

The present section considers solution of boundary-value problem (4.2)-(4.4) by means of the Fourier method (e.g., [

・ 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., [

w ( x , t ) = ∑ i = 1 ∞ f i ( t ) X i ( x ) , 0 ≤ x ≤ h (5.1)

where X i ( x ) are the orthonormed eigenfunctions of the Laplace operator in PDE (4.2). This operator is differential expression ∂ 2 / ∂ x 2 endowed with boundary conditions (4.3) and (4.4). As one can easily check, the eigenvalues and eigenfunctions of the operator are numbers − κ i 2 , i = 1,2, ⋯ , where

κ i = π h 2 i − 1 2 , i = 1 , 2 , ⋯ (5.2)

and

X i ( x ) = 2 / h cos ( κ i x ) , 0 ≤ x ≤ h , i = 1 , 2 , ⋯ . (5.3)

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., [

h − x = ∑ i = 1 ∞ [ ∫ 0 h ( h − x * ) X i ( x * ) d x * ] X i ( x ) = 2 / h ∑ i = 1 ∞ κ i − 2 X i ( x ) , 0 ≤ x ≤ h . (5.4)

Application of (5.1) and (5.4) to PDE (4.2) results in

θ d 3 f i ( t ) / d t 3 + d 2 f i ( t ) / d t 2 + 2 θ s 2 κ i 2 d f i ( t ) / d t + s 2 κ i 2 f i ( t ) = − 2 / h a ¯ ( t ) κ i − 2 , i = 1 , 2 , ⋯ (5.5)

Relations (5.5) present the countably infinite set of linear non-autonomous ODEs of the third order for time-dependent coefficients f i ( t ) in the Fourier expansion (5.1). The forms of these ODEs show that they can be mutually independent.

The third-order ODEs (5.5) can be represented in the form of the first-order ODE systems

d F i ( t ) / d t = Γ i F i ( t ) − 2 / h a ¯ ( t ) ϕ i , i = 1 , 2 , ⋯ , (5.6)

where vectors F i ( t ) and ϕ i , and matrixes Γ i are described with the following expressions

F i ( t ) = ( f i ( t ) d f i ( t ) / d t d 2 f i ( t ) / d t 2 ) , i = 1 , 2 , ⋯ , (5.7)

ϕ i = ( 0 0 θ − 1 κ i − 2 ) , i = 1 , 2 , ⋯ , (5.8)

Γ i = ( 0 1 0 0 0 1 − θ − 1 − 2 s 2 κ i 2 − θ − 1 s 2 κ i 2 ) , i = 1 , 2 , ⋯ . (5.9)

As is well known in theory of systems of linear ODEs, if systems (5.6) are asymptotically stable, then:

・ the steady-state solutions of these systems (under Assumption 4.1) are described as

F i ( t ) = − 2 / h ∫ − ∞ t a ¯ ( t * ) Φ i ( t − t * ) ϕ i d t * , i = 1 , 2 , ⋯ (5.10)

where

Φ i ( t − t * ) = exp [ ( t − t * ) Γ i ] , i = 1 , 2 , ⋯ (5.11)

are the Cauchy matrixes of the systems

・ and all entries of these matrixes exponentially tend to zero as t − t * → ∞ .

Note that column-vectors Φ i .1 ( t − t * ) , Φ i .2 ( t − t * ) , and Φ i .3 ( t − t * ) of matrixes (5.11) are the solutions of homogeneous version

d F i ( t ) / d t = Γ i F i ( t ) , i = 1 , 2 , ⋯ (5.12)

of the ODE system (5.6), which are such that (see (5.11)) Φ i ( 0 ) = I , i = 1 , 2 , ⋯ , where I is the identity 3 × 3 -matrix, or, equivalently,

Φ i . l ( 0 ) = e l , l = 1 , 2 , 3 , i = 1 , 2 , ⋯ (5.13)

where e l is the l th column of matrix I . Due to (5.8), relations (5.10) are reduced to expressions

F i ( t ) = − 2 / h θ κ i 2 ∫ − ∞ t a ¯ ( t * ) Φ i .3 ( t − t * ) d t * , i = 1 , 2 , ⋯ , (5.14)

which do not include matrixes (5.11) and emphasize a special role of their third columns Φ i .3 ( t − t * ) . In view of the structure of vectors (5.7), it is sufficient to consider the first entries of vectors (5.14) only. The versions of (5.14) for these entries are

f i ( t ) = − 2 / h θ κ i 2 ∫ − ∞ t Φ i .13 ( t − t * ) a ¯ ( t * ) d t * , i = 1 , 2 , ⋯ (5.15)

where Φ i .13 ( t − t * ) is the first entry of vector Φ i .3 ( t − t * ) . Thus, the analytical expressions for the time-dependent coefficients f i ( t ) in expansion (5.1) are (5.15). Note that the physical dimension of Φ i .3 ( t − t * ) is the squared time.

Assumption 5.1. Each of the functions d 2 a ( t ) / d t 2 and d 3 a ( t ) / d t 3 , 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

Δ P ( x , t ) = − ρ { a ( t ) ( x − h ) + 2 h 1 θ ∑ i = 1 ∞ [ 1 κ i 2 ∫ − ∞ t Φ i .13 ( t − t * ) a ¯ ( t * ) d t * ] cos ( κ i x ) } , 0 ≤ x ≤ h (5.16)

Remark 5.1. The exact analytical expression for the steady-state solution of boundary-value problem (3.1)-(3.3) is described with expression (5.16) (see also (4.5) and (5.2)) in terms of entries Φ i .13 ( t − t * ) of the Cauchy matrixes (5.11) of ODE systems (5.6) under Assumptions 4.1 and 5.1, as well as the assumption that all of the ODE systems (5.6) are asymptotically stable. ,

The next section establishes asymptotic stability of systems (5.6) and derives the exact analytical expressions for the entries mentioned in Remark 5.1.

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

θ λ 3 + λ 2 + 2 θ ( s κ i ) 2 λ + ( s κ i ) 2 = 0 , i = 1 , 2 , ⋯ , (6.1)

which are the characteristic equations for homogeneous versions

θ d 3 f i ( t ) / d t 3 + d 2 f i ( t ) / d t 2 + 2 θ s 2 κ i 2 d f i ( t ) / d t + s 2 κ i 2 f i ( t ) = 0 , i = 1 , 2 , ⋯ , (6.2)

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 [

(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 .

(ii) 1 / θ i + 2 / Θ i = 1 / θ .

(iii) lim s κ i → 0 θ i = θ , lim s κ i → ∞ θ i = 2 θ , and θ i strictly monotonically increases as s κ i increases from zero to infinity.

(iv) Asymptotic representation for 1 / Θ i in the limit case as s κ i → 0 is 1 / Θ i = θ ( s κ i ) 2 / 2 .

(v) ω i = Θ i ( s κ i ) 2 2 θ − ( 1 θ − 1 Θ i ) 2 .

(vi) Asymptotic representations for ω i in the limit cases as s κ i → 0 and s κ i → ∞ are ω i = s κ i and ω i = 2 s κ i , respectively.

(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 s κ i , 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), λ = − 1 / θ 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., [

1 θ i = 1 − M ( θ s κ i ) θ , i = 1 , 2 , ⋯ (6.3)

where

M ( χ ) = M + ( χ ) + M − ( χ ) (6.4)

and

M ± ( χ ) = 1 3 [ 1 − 1 + 9 χ 2 ± 3 | χ | 3 ( 32 χ 4 − 13 χ 2 + 4 ) 2 3 ] . (6.5)

Continuous variable χ in (6.4) and (6.5) generalizes discrete variable θ s κ i > 0 in (6.3). Numerical values of function (6.4) are exemplified in

In view of (6.3), the equality in statement (ii) is equivalent to

1 Θ i = M ( θ s κ i ) 2 θ , i = 1 , 2 , ⋯ , (6.6)

which, in turn, transforms the expression in statement (v) into

ω i = 1 θ ( θ s κ i ) 2 M ( θ s κ i ) − [ 1 − M ( θ s κ i ) 2 ] 2 = 1 θ ( θ s κ i ) 2 M ( θ s κ i ) − 1 + M ( θ s κ i ) − [ M ( θ s κ i ) ] 2 4 , i = 1 , 2 , ⋯ . (6.7)

Remark 6.1. As follows from statements (iii) and (iv), as well as expressions

χ | M ( χ ) | χ | M ( χ ) | χ | M ( χ ) | χ | M ( χ ) | χ | M ( χ ) | χ | M ( χ ) |
---|---|---|---|---|---|---|---|---|---|---|---|

0.0 | 0 | 2.0 | 0.48389 | 4.0 | 0.49606 | 6.0 | 0.49826 | 8.0 | 0.49902 | 10.0 | 0.49937 |

0.2 | 0.0399 | 2.2 | 0.48675 | 4.2 | 0.49643 | 6.2 | 0.49837 | 8.2 | 0.49907 | 10.2 | 0.4994 |

0.4 | 0.15465 | 2.4 | 0.48891 | 4.4 | 0.49675 | 6.4 | 0.49847 | 8.4 | 0.49911 | 10.4 | 0.49942 |

0.6 | 0.29622 | 2.6 | 0.49058 | 4.6 | 0.49703 | 6.6 | 0.49856 | 8.6 | 0.49915 | 10.6 | 0.49944 |

0.8 | 0.38634 | 2.8 | 0.49190 | 4.8 | 0.49727 | 6.8 | 0.49865 | 8.8 | 0.49919 | 10.8 | 0.49946 |

1.0 | 0.43016 | 3.0 | 0.49296 | 5.0 | 0.49749 | 7.0 | 0.49872 | 9.0 | 0.49923 | 11.0 | 0.49948 |

1.2 | 0.45293 | 3.2 | 0.49382 | 5.2 | 0.49768 | 7.2 | 0.49879 | 9.2 | 0.49926 | 11.2 | 0.49950 |

1.4 | 0.46611 | 3.4 | 0.49454 | 5.4 | 0.49785 | 7.4 | 0.49886 | 9.4 | 0.49929 | 11.4 | 0.49952 |

1.6 | 0.4744 | 3.6 | 0.49513 | 5.6 | 0.49800 | 7.6 | 0.49892 | 9.6 | 0.49932 | 11.6 | 0.49954 |

1.8 | 0.47997 | 3.8 | 0.49563 | 5.8 | 0.49814 | 7.8 | 0.49897 | 9.8 | 0.49935 | 11.8 | 0.49955 |

(6.3) and (6.6), the properties below hold.

・ l i m χ → 0 M ( χ ) = 0

・ l i m χ → ∞ M ( χ ) = 1 / 2

・ Function M ( χ ) strictly monotonically increases as χ increases from zero to infinity.

・ The asymptotic representation for M ( χ ) in the limit case as χ → 0 is M ( χ ) = χ 2 . ,

Remark 6.2. As follows from the fist two bullets in Remark 6.1 and expression (6.7), relations

lim θ ↓ 0 ω i = s κ i , lim θ → ∞ ω i = 2 s κ i , i = 1 , 2 , ⋯

hold. ,

Remark 6.3. The above consideration enables one to express three sequences { θ i } i ≥ 1 , { Θ i } i ≥ 1 , and { ω i } i ≥ 1 of the parameters of the real root and complex conjugate roots in statement (i) in terms of parameter θ , physically dimensionless function M ( ⋅ ) , and physically dimensionless sequence { θ s κ i } i ≥ 1 where s and κ i are determined with (2.2) and (5.2) by means of (6.3), (6.6), and (6.7). ,

According to theory of linear ODEs with time-independent coefficients, the roots of Equation (6.1), which are indicated in Remark 6.2, enable one to specify the general solutions of the third-order scalar ODEs (6.2) at the initial time point, s , which was already used in expressions (5.11). These solutions are

g ^ i ( t − t * ) = E i exp [ − ( t − t * ) / θ i ] + { S i sin [ ω i ( t − t * ) ] + C i cos [ ω i ( t − t * ) ] } exp [ − ( t − t * ) / Θ i ] , i = 1 , 2 , ⋯ , (6.8)

where E i , S i , and C i are the arbitrary parameters independent of time. This, in particular, means that

θ d 3 g ^ i ( t − t * ) / d t 3 + d 2 g ^ i ( t − t * ) / d t 2 + 2 θ s 2 κ i 2 d g ^ i ( t − t * ) / d t + s 2 κ i 2 g ^ i ( t − t * ) = 0 , i = 1 , 2 , ⋯ . (6.9)

The forms of the corresponding solutions of the first-order ODE systems (5.12) are

G ^ i ( t − t * ) = ( g ^ i ( t − t * ) d g ^ i ( t − t * ) / d t d 2 g ^ i ( t − t * ) / d t 2 ) , i = 1 , 2 , ⋯ . (6.10)

These vector-functions describe all of the solutions of the systems specifying them with individual values of parameters A i , B i , and C i . In particular, vectors (6.10) describe the three columns of each of the Cauchy matrixes (5.11) by means of conditions (5.13) with expressions Φ i . l ( t − t * ) = G ^ i ( t − t * ) | G ^ i ( 0 ) = e l , l = 1 , 2 , 3 , i = 1 , 2 , ⋯ . Consequently, the third columns, which include entries Φ i .13 ( t − t * ) used in (5.16), are determined as Φ i .3 ( t − t * ) = G ^ i ( t − t * ) | G ^ i ( 0 ) = e 3 , i = 1 , 2 , ⋯ , or, in more detail,

Φ i .3 ( t − t * ) = ( g ^ i ( t − t * ) d g ^ i ( t − t * ) / d t d 2 g ^ i ( t − t * ) / d 2 t ) | g ^ i ( 0 ) = 0 , d g ^ i ( 0 ) / d t = 0 , d 2 g ^ i ( 0 ) / d 2 t = 1 , i = 1 , 2 , ⋯ . (6.11)

In order to carry out the calculations presumed by (6.11), one needs to evaluate the first and second time derivatives of function g ^ i ( t − s ) . As follows from (6.8),

d g ^ i ( t − t * ) d t = − E i θ i exp ( − t − t * θ i ) + { ( − S i − Θ i ω i C i ) sin [ ω i ( t − t * ) ] + ( Θ i ω i S i − C i ) cos [ ω i ( t − t * ) ] } 1 Θ i exp ( − t − t * Θ i ) , i = 1 , 2 , ⋯ , (6.12)

d 2 g ^ i ( t − t * ) d t 2 = E i θ i 2 e x p ( − t − t * θ i ) + { { [ 1 − ( Θ i ω i ) 2 ] S i + 2 Θ i ω i C i } sin [ ω i ( t − t * ) ] + { − 2 Θ i ω i S i + [ 1 − ( Θ i ω i ) 2 ] C i } cos [ ω i ( t − t * ) ] } × 1 Θ i 2 e x p ( − t − t * Θ i ) , i = 1 , 2 , ⋯ . (6.13)

Entries Φ i .13 ( 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

g ^ i ( 0 ) = E i + C i = 0 , i = 1 , 2 , ⋯ ,

d g ^ i ( 0 ) d t = − E i θ i + Θ i ω i S i − C i Θ i = 0 , i = 1 , 2 , ⋯ ,

d 2 g ^ i ( 0 ) d t 2 = E i θ i 2 + − 2 Θ i ω i S i + [ 1 − ( Θ i ω i ) 2 ] C i Θ i 2 = 1 , i = 1 , 2 , ⋯ .

These systems are the ones for determination of parameters E i , S i , and C i . One can readily check that the corresponding solutions are

E i = − C i = 1 ( 1 / θ i − 1 / Θ i ) 2 + ω i 2 , i = 1 , 2 , ⋯ ,

S i = 1 ω i ( 1 / θ i − 1 / Θ i ) A i , i = 1 , 2 , ⋯ ,

or, in view of (6.3) and (6.6),

E i = − C i = θ 2 / D i , i = 1 , 2 , ⋯ , (6.14)

S i = γ i θ ω i E i , i = 1 , 2 , ⋯ , (6.15)

where

D i = γ i 2 + θ 2 ω i 2 , i = 1 , 2 , ⋯ , (6.16)

γ i = 1 − ( 3 / 2 ) ( θ s κ i ) , i = 1 , 2 , ⋯ . (6.17)

Note that parameters (6.16) and (6.17) are physically dimensionless. The remark below is similar to Remark 6.2.

Remark 6.4. As follows from the fist two bullets in Remark 6.1, Remark 6.2, and express- ions (6.15)-(6.17), relations

lim θ ↓ 0 γ i = 1 , lim θ → ∞ γ i = 1 / 4 , i = 1 , 2 , ⋯ ,

lim θ ↓ 0 ( θ 2 / D i ) = 0 , lim θ → ∞ ( θ 2 / D i ) = ( 1 / 2 ) / ( s κ i ) 2 , i = 1 , 2 , ⋯ ,

hold. ,

Remark 6.5. Three sequences { E i } i ≥ 1 , { S i } i ≥ 1 , and { C i } i ≥ 1 are expressed in terms of SRT θ , physically dimensionless function M ( ⋅ ) (see (6.4)), and physically dimensionless sequence { θ s κ i } i ≥ 1 by means of (6.7) and (6.14)-(6.17). ,

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 D i in the limit case as i → ∞ is D i = 2 ( θ s κ i ) 2 . ,

As is well known (e.g., [

sin [ ω i ( t − t ∗ ) ] = sin ( ω i t ) cos ( ω i t ∗ ) − cos ( ω i t ) sin ( ω i t ∗ ) , i = 1 , 2 , ⋯ ,

cos [ ω i ( t − t ∗ ) ] = cos ( ω i t ) cos ( ω i t ∗ ) + sin ( ω i t ) sin ( ω i t ∗ ) , i = 1 , 2 , ⋯ .

Substituting these expressions into (6.8) and taking into account the first entries of vectors (6.11), one obtains

Φ i .13 ( t − t * ) = E i exp [ − ( t − t * ) / θ i ] + { [ C i sin ( ω i t ∗ ) + S i cos ( ω i t ∗ ) ] sin ( ω i t ) + [ − S i sin ( ω i t ∗ ) + C i cos ( ω i t ∗ ) ] cos ( ω i t ) } × e x p [ − ( t − t ∗ ) / Θ i ] , i = 1 , 2 , ⋯ ,

which, after application of (6.14)-(6.17), become

Φ i .13 ( t − t * ) = ( θ / D i ) × { θ e x p ( − t − t * Θ i ) − { [ θ sin ( ω i t ∗ ) − ( γ i / ω i ) cos ( ω i t ∗ ) ] sin ( ω i t ) + [ ( γ i / ω i ) sin ( ω i t ∗ ) + θ cos ( ω i t ∗ ) ] cos ( ω i t ) } × e x p ( − t − t * Θ i ) } , i = 1 , 2 , ⋯ .

These expressions enable one to calculate the integrals in (5.16) as

∫ − ∞ t Φ i .13 ( t − t * ) a ¯ ( s ) d s = θ D i [ E i ∨ ( t , a ¯ ( ⋅ ) ) + S i ∨ ( t , a ¯ ( ⋅ ) ) sin ( ω i t ) + C i ∨ ( t , a ¯ ( ⋅ ) ) cos ( ω i t ) ] , i = 1 , 2 , ⋯ , (6.18)

where

E i ∨ ( t , a ¯ ( ⋅ ) ) = θ exp ( − t / θ i ) ∫ − ∞ t exp ( t ∗ / θ i ) a ¯ ( t * ) d t * , i = 1 , 2 , ⋯ , (6.19)

S i ∨ ( t , a ¯ ( ⋅ ) ) = − exp ( − t / Θ i ) ∫ − ∞ t exp ( t ∗ / Θ i ) [ θ sin ( ω i t ∗ ) − ( γ i / ω i ) cos ( ω i t ∗ ) ] a ¯ ( t * ) d t * , i = 1 , 2 , ⋯ , (6.20)

C i ∨ ( t , a ¯ ( ⋅ ) ) = − exp ( − t / Θ i ) ∫ − ∞ t exp ( t ∗ / Θ i ) [ ( γ i / ω i ) sin ( ω i t ∗ ) + θ cos ( ω i t ∗ ) ] a ¯ ( t * ) d t * , i = 1 , 2 , ⋯ . (6.21)

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

Δ P ( x , t ) = − ρ × { a ( t ) ( x − h ) + 2 h ∑ i = 1 ∞ E i ∨ ( t , a ¯ ( ⋅ ) ) + S i ∨ ( t , a ¯ ( ⋅ ) ) sin ( ω i t ) + C i ∨ ( t , a ¯ ( ⋅ ) ) cos ( ω i t ) κ i 2 D i cos ( κ i x ) } , 0 ≤ x ≤ h . (6.22)

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.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).

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

[ S i ∨ ( t , a ¯ ( ⋅ ) ) ] 2 + [ C i ∨ ( t , a ¯ ( ⋅ ) ) ] 2 ,

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.

Boundary-value problem (3.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.

The exact analytical expression for the steady-state solution of the problem is, according to Remark 5.1, expression (6.22) (see also (4.5), (5.2), (6.16), and (6.17)) under Assumptions 4.1 and 5.1, and with coefficients (6.19)-(6.21). Parameters s , κ i , θ i , Θ i , and ω i are determined with (2.2), (5.2), (6.3), (6.6), and (6.7), respectively, in terms of the input parameters ρ , K , h , and θ .

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., [

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.

As follows from Remark 3.1, in addition to solution (6.22), there is one more space-time-dependent output variable: the acoustic pressure. The corresponding expressions are derived in this section.

Expression (6.22) presents the steady-state acoustic ANS. The related acoustic pressure Δ p ( t ) can be obtained as the steady-state solution of ODE (2.8) regarded as the equation for Δ p ( t ) .

The steady-state solution of ODE (2.8) is

Δ p ( x , t ) = 1 2 θ ∫ − ∞ t exp ( − t − t * 2 θ ) Δ P ( x , t * ) d t * + 1 2 ∫ − ∞ t exp ( − t − t * 2 θ ) ∂ Δ P ( x , t * ) ∂ t * d t * , 0 ≤ x ≤ h . (8.1)

We accept the following assumption.

Assumption 8.1. Function Δ P ( x , t ) , as a function of t at every fixed x ∈ [ 0, h ] , is uni- formly bounded in the entire t-axis. ,

Integrating the second integral by parts and taking into account Assumption 8.1, one obtains that

∫ − ∞ t exp ( − t − t * 2 θ ) ∂ Δ P ( x , t * ) ∂ t * d t * = Δ P ( x , t ) − 1 2 θ ∫ − ∞ t exp ( − t − t * 2 θ ) Δ P ( x , t * ) d t * , 0 ≤ x ≤ h . (8.2)

Substitution of (8.2) into (8.1) leads to

Δ p ( x , t ) = 1 2 [ Δ P ( x , t ) + 1 2 θ ∫ − ∞ t exp ( − t − t * 2 θ ) Δ P ( x , t * ) d t * ] , 0 ≤ x ≤ h . (8.3)

In the limit case as θ ↓ 0 , the integral in (8.3) tends to Δ P ( x , t ) .

Term Δ P ( x , t ) in (8.3) is determined with (6.22) (see also (6.19)-(6.21)). Regarding the convergence of the series in (6.22), one can note the following.

Remark 8.1. One can readily check that each of the coefficients (6.19)-(6.21) of the series in (6.22) is uniformly bounded in x ∈ [ 0, h ] and t ∈ ( − ∞ , ∞ ) . Moreover, the i th term, i = 1 , 2 , ⋯ , of the series in the limit case as i → ∞ , is proportional to κ i − 4 or, equivalently (see (5.2)), to ( 2 i − 1 ) − 4 . 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 ∑ i = 1 ∞ ( 2 i − 1 ) − 4 converges (e.g., [

Remark 8.1 shows that representation (6.22) is properly defined.

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 θ : h = 10 − 2 , ρ = 1.04 × 10 3 , K = 2.2 × 10 9 , and θ = θ * where

θ * = 0.95 × 10 − 4 (9.1)

The value of θ is determined as described below.

As follows from the discussion on pp. 1961 and 1966 in [

The present work considers the three examples. They are based on the above values of h , ρ , and K at different values of θ indicated in the captions for Figures 1-12. The examples also differ with different external-acceleration functions of the following two forms

a ( t ) = { 0 , t < 0 or t > T p 10 4 { 1 − cos [ ( 2 π / T p ) t ] } / 2 , 0 ≤ t ≤ T p (9.2)

a ( t ) = { 0 , t < 0 10 4 { 1 − cos [ ( 2 π / T p ) t ] } / 2 , 0 ≤ t (9.3)

where T p > 0 is the parameter of function a ( t ) . More specifically, function a ( t ) is: (9.2) at T p = 10 − 5 in Example 1, (9.3) at T p = 10 − 5 in Example 2, and (9.2) at T p = 10 − 6 in Example 3. In all of the examples, the pressures in Figures 2-10 and

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 [

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 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 bottom.

Analysis of Figures 1-12 enables one to indicate a few features.

1) Comparison of

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

The influence of the SRT upon the pressure is more pronounced in

Figures 5-9 corresponding to the SRT values θ = 0.010 × θ * , θ = 0.030 × θ * , θ = 0.100 × θ * , θ = 0.300 × θ * , and θ = 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.,

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 T p = 10 − 5 , i.e., with the angular frequency

ω = 2 π / T p = 2 π × 10 5 ≈ 628319 (9.4)

As is well known (e.g., [

Note that, in the case of Example 3 in

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

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.7 emphasizes the new components of the solution, which are due to the exponential-NSRF- originated 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

SRT θ | Pair of the angular natural frequencies | Is the smallest of the two frequencies the principal one? |
---|---|---|

1 0.000 × θ * | 228,463, 685,389 | yes |

2 0.001 × θ * | 228,499, 686,478 | yes |

3 0.003 × θ * | 228,826, 695,624 | yes |

4 0.010 × θ * | 232,716, 816,686 | yes |

5 0.030 × θ * | 272,229, 949,518 | yes |

6 0.100 × θ * | 317,753, 967,498 | yes |

7 0.300 × θ * | 322,499, 969,085 | yes |

8 1.000 × θ * | 323,041, 969,266 | yes |

amplitudes). These effects are especially clearly seen in Figures 5-9. The presented behaviors and values of Δ p ( x , t ) unavailable with common wave PDE (2.10) or the Stokes-type PDE (2.9).

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.

5) The above results are obtained in the exponential approximation for NSRF Ψ ( ψ ) (see the text above (2.6)). The exponential dependence qualitatively correctly represents the behaviors of the NSRF of more complex shapes. Consequently, the qualitative aspects of the obtained results remain valid in the general case where the NSRF need not be exponential. Also, since NSRF Ψ ( ψ ) enters the general PIDE (2.1) under the sign of integral, the errors of one or another approximation for the NSRF affect the quantitative behaviors of the acoustic ANS Δ P ( x , t ) and, thus (see (8.3)), the acoustic pressure Δ p ( x , t ) not very significantly. We also note the following.

Remark 9.1. In view of the value K = 2.2 × 10 9 (see above) and relation (2.4), the pressure values in Figures 1-12 meet requirement (2.4). This shows that, in all of the mentioned cases, the linearity of the present model is adequate. ,

Remark 9.2. The present work assumes that the soft tissue under consideration is linearly viscoelastic. The work emphasizes condition (2.4) for the linear-model adequacy and check it with the obtained results. It appears that the linear-model results are adequate.

In contrast to that, work [

Remark 9.3. Work [

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., [

One can also note that paper [^{5} and 10^{6} used in the above Examples 1 - 3). The conclusion of [^{−3}) at all frequencies in the range ( 10 3 , ∞ ) (see [

The results of the present works are summarized in the next section.

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 (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 [

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 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.

The authors are grateful to Vinnova, Sweden’s Innovation Agency, for a partial support of the present work from the project 2015-00352 “Zero Vibration Injuries, Phase 2’’. The authors thank Hans Lindell, Swerea IVF AB, Mölndal, Sweden, for valuable information on hand-held impact machines design and performance, numerous discussions, and reading the manuscript. The authors also thank Dr. Henrik Enquist, Medical Services, Laboratory Medicine, Region Skåne, Lund, Sweden, for drawing their attention to reference [

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

The VWF injuries or, in a more professional language, lesions in human fingers are well known. They are described in work [

Paper [

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

this respect, the mentioned populations are similar to other long-lasting lesions such as scars or other infiltrates in chronic inflammation, and lesions caused by homeorhetic dysfunctions. The infiltrates and the notion of a homeorhetic dysfunction are discussed, for example, in [

The skin of the VWF patients also contains lesions. The experimental studies of [

As is well known, fingers, in contrast to palms, do not contain muscles (other than arrector pili). The soft tissues of a finger comprise the epidermis, dermis, and subcutaneous tissue. The soft tissues of a palm include the muscle tissue as well. The lesions in fingers of the VWF patients are also present in palms, first of all, in the palm regions neighboring the fingers (e.g., [

Thus, the vibration injuries to be prevented are VSMCs with noticeably increased volumes or, briefly, enlarged VSMCs. Enlargement of a cell is literally a geometric effect. It has a number of mechanical implications. However, the fact that, as is noted above, enlargement of the cells was emphasized by histopathologists (e.g., [

In order to provide the aforementioned prevention, it is necessary to develop a quantitative acoustic model that would include the influence of the acceleration of the vibrating hand-held tool upon the volumes of VSMCs. The questions to be answered on the way to this model are the following.

1) What characteristics of an SLT directly influence upon enlargements of the cells in the SLT?

2) How specifically is the mentioned influence implemented?

3) How can one reduce the resulting enlargements of the cells?

The answers to all of these questions are suggested in the present section. It deals with a cell, which is located at fixed spatial point x in an SLT, and acoustic pressure Δ p = Δ p ( x , t ) in the SLT (see (8.3)). In what follows, we, for brevity, do not explicitly indicate the dependences on x but acknowledge their presence by denoting the time differential as the partial differential ∂ t .

Let V c = V c ( t ) and K c be the volume and bulk modulus of the cell, respectively. We also consider the pressures p c = p c ( x , t ) and p i f = p i f ( x , t ) in the cell and the interstitial fluid (IF) surrounding the cell in the SLT. The equilibrium versions of these pressures are denoted with p ¯ c = p ¯ c ( x ) and p ¯ i f = p ¯ i f ( x ) . Then the corresponding acoustic pressures are expressed as

Δ p c = p c − p ¯ c , Δ p i f = p i f − p ¯ i f (A.1)

where the sign “overline’’ denotes the equilibrium version of the corresponding quantity. This convention is used down to the end of the present section.

As is assumed above (see (2.4)), the SLT is a linear solid. Similarly, we regard the cell as a solid body and assume that it is also a linear solid, i.e.,

| Δ p c | / K c ≪ 1 (A.2)

Note that the property (A.2) also means (e.g., [

| ε c | ≪ 1 (A.3)

where

ε c = ( V c − V ¯ c ) / V ¯ c (A.4)

is the relative change of the cell volume.

Taking into account the linearity and the well-known continuum-mechanics result ( [

∂ V c ∂ t + 1 K c ∂ Δ p c ∂ t V c = − Φ c (A.5)

where Φ c is the cell-volume rate due to the phenomena that are not represented with the terms on the left-hand side. This rate is expressed with the well-known Starling equation (e.g., [

Φ c = H c S c [ p c − p i f − σ c k B T ( n c * − n i f * ) ] (A.6)

where H c is the hydraulic conductance of the semi-permeable membrane of the cell, S c is the area of the cell-membrane surface, σ c is the reflection coefficient of the cell membrane (e.g., [

Note that the sign “minus’’ in front of the term (A.6) in equation (A.5) is due to the well-known facts (e.g., [

If the term in the brackets in (A.6) is positive, there is a flow of fluid from the ICF into the IF. This flow decreases the cell volume.

If the term in the brackets in (A.6) is negative, there is a flow of fluid from the IF into the ICF. This flow increases the cell volume.

At equilibrium, the term in the brackets in (A.6) is zero, i.e., p ¯ c − p ¯ i f − σ c k B T ( n ¯ c * − n ¯ i f * ) = 0 . Rewriting (A.6) in terms of (A.1) and allowing for the latter equality, one obtains

Φ c = H c S c { Δ p c − Δ p i f − σ c k B T [ ( n c * − n ¯ c * ) − ( n i f * − n ¯ i f * ) ] } . (A.7)

We assume that concentration n i f * remains unchanged, i.e.,

n i f * = n ¯ i f * (A.8)

and the number of the solute molecules in the ICF, n ¯ c * V ¯ c , remains unchanged as well that results in

n c * = n ¯ c * ( V ¯ c / V c ) . (A.9)

Combination of (A.7)-(A.9) enables one to present Equation (A.5) as follows

∂ V c ∂ t + 1 K c ∂ Δ p c ∂ t V c = − H c S c { Δ p c − Δ p i f − σ c k B T n ¯ c * [ ( V ¯ c / V c ) − 1 ] } . (A.10)

This equation includes three acoustic pressures: Δ p , which is modeled in the main text of the present work, as well as Δ p c and Δ p i f , which are not modeled in the work. In order to simplify Equation (A.10) to the form where the available pressure is the only one, which is used, we assume that Δ p c = Δ p and Δ p c = Δ p i f . These relations transform (A.10) into

∂ V c ∂ t + 1 K c ∂ Δ p c ∂ t V c = H c σ c k B T n ¯ c * S c [ ( V ¯ c / V c ) − 1 ] (A.11)

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 ( [

S c = S ¯ c ( 1 + ε c 2 / 3 ) (A.12)

Substituting (A.4) into (A.11) and (A.12), one obtains

V ¯ c ∂ ε c ∂ t + 1 K c ∂ Δ p ∂ t V ¯ c ( 1 + ε c ) = − H c σ c k B T n ¯ c * S ¯ c ( 1 + ε c 2 / 3 ) ε c 1 + ε c

or, equivalently,

∂ ε c ∂ t = − 1 τ c 1 + ε c 2 / 3 1 + ε c ε c − 1 K c ∂ Δ p ∂ t ε c − 1 K c ∂ Δ p ∂ t (A.13)

where

τ c − 1 = H c σ c k B T n ¯ c * S ¯ c / V ¯ c (A.14)

is the relaxation time induced by the cell osmosis. We consider ODE (A.13) in a semi-infinite time interval t > t o under the initial condition that the cell volume is at equilibrium at point t = t o , i.e. (see (A.4)),

ε c | t = t o = 0. (A.15)

Example A.1. According to experimental study in [

τ c − 1 ≈ H c σ c k B T n ¯ c * × 1.44 × 10 6 . (A.16)

Next, at normal human temperature of about 37˚C, T = 310 . This results in k B T ≈ 4.28 × 10 − 21 that transforms (A.16) into

τ c − 1 ≈ H c σ c n ¯ c * × 6.16 × 10 − 15 . (A.17)

Parameter H c , the hydraulic conductance of the VSMC membrane, can be taken from [

τ c − 1 ≈ σ c n ¯ c * × 7.39 × 10 − 25 (A.18)

We assume that the well-known typical value of σ c n ¯ c for an intracellular fluid (e.g., [

Remark A.1. Multiplying Equation (A.13) considered in interval t > t o by to τ c and passing in the resulting equation to the limit as τ c ↓ 0 , one obtains that ε c = 0 at all t > t o . This result and (A.15) lead to property

lim τ c ↓ 0 ε c = 0 , t ≥ t o (A.19)

It indicates that the cell osmosis plays a mechanoprotective role in the cell-volume maintenance. ,

We note note a few features of the problem under consideration.

Vibrating hand-held tools (e.g., impact wrenches outlined in Section 1) are usually used in such a way that the periods when the tool is active (continuously functions) alternate with the periods when the tool is inactive (does not function). Any of the activity periods is bounded and the total number of them is finite. In each of the inactivity periods, acceleration a ( t ) in (6.22) is a function with bounded support, and thus Δ P , and, in view of (8.3), Δ p exponentially tend to zero (cp., (6.19)-(6.22)). In view of these features and the fact that Δ p , as function of time, is uniformly bounded (e.g., see (A.2)), property

lim t → ∞ ∫ t o t | Δ p ( t * ) | d t * exists and is finite, (A.20)

holds. It in particular implies that

lim t → ∞ Δ p ( t ) = 0 (A.21)

Vibrations, which are produced by the handles of active hand-held tools on the fingers of the users, present series of more or less regular pulses of the acoustic pressure in the finger SLT Δ p ( t ) , which are compressive, i.e.,

values of Δ p ( t ) at all t ≥ t o are predominantly positive; this feature is illustrated with any of Figures 3-10. (A.22)

As is noted in the second paragraph of the present section, an enlarged cell is a significant lesion if its volume is about 1.5 greater than the volume of a normal cell. This figure corresponds to ε c ≈ 0.5 . However, such lesions become measurable only after a sufficiently long time of use of the vibrating tools. More specifically, the cell enlargement is an additive effect of a large number of the activity periods, i.e., a significant relative enlargement is achieved as the result of a long series of enlargements for small amounts. This means that it is sufficient to consider ODE (A.13) at small values of | ε c | prescribed by (A.3) in any period of the tool activity. This leads to reduction (A.13) to equation

∂ ε c ∂ t = − ( 1 τ c + 1 K c ∂ Δ p ∂ t ) ε c − 1 K c ∂ Δ p ∂ t (A.23)

which is a linear ODE.

The solution on the initial-value problem (A.23), (A.15) are studied in the following theorem.

Theorem A.1. Let R ( t ) be the range of function Δ p ( t ) in interval [ t o , t ] , i.e.,

R ( t ) = M ( t ) − m ( t ) , t ≥ t o , (A.24)

where

m ( t ) = min t * ∈ [ t o , t ] Δ p ( t * ) , M ( t ) = max t * ∈ [ t o , t ] Δ p ( t * ) , t ≥ t o . (A.25)

It is assumed that

R ( t ) / K c ≪ 1 , t ≥ t o . (A.26)

We also accept notation

ε c . ∞ = lim t → ∞ ε c (A.27)

Then the following assertions are valid.

(i) The solution of initial-value problem (A.23), (A.15) is

ε c = 1 τ c ∫ t o t exp { − [ t − t * τ c + Δ p ( t ) − Δ p ( t ∗ ) K c ] } ∂ t * + exp { − [ t − t o τ c + Δ p ( t ) − Δ p ( t o ) K c ] } − 1 , t > t o (A.28)

(ii) Inequalities

exp [ m ( t ) − Δ p ( t ) K c ] − 1 ≤ ε c ≤ exp [ M ( t ) − Δ p ( t ) K c ] − 1 , t ≥ t o (A.29)

exp [ − R ( t ) / K c ] − 1 ≤ ε c ≤ exp [ − R ( t ) / K c ] − 1 , t ≥ t o (A.30)

hold.

(iii) Linear ODE (A.23) is applicable.

(iv) Relation

ε c . ∞ = lim t → ∞ { 1 τ c ∫ t o t exp [ − t − t * τ c + Δ p ( t ∗ ) K c ] ∂ t ∗ } − 1 (A.31)

is valid.

Proof. According to theory of non-autonomous linear ODEs, the solution of initial-value problem (A.23), (A.15) is

ε c = − exp { − [ t − t o τ c + 1 K c ∫ t o t Δ p ( t ∗ ) K c ∂ t * ] } × ∫ t o t exp [ t * − t o τ c + 1 K c ∫ t o t ∗ ∂ Δ p ( t * ) ∂ t * ∂ t * ] ∂ Δ p ( t ∗ ) ∂ t * ∂ t * K c , t > t o ,

or, equivalently,

ε c = − exp { − [ t − t o τ c + Δ p ( t ) − Δ p ( t o ) K c ] } × ∫ t o t exp ( t * − t o τ c ) ∂ ∂ t * exp [ 1 K c ∫ t o t ∗ ∂ Δ p ( t * ) ∂ t * ∂ t * ] ∂ t * , t > t o (A.32)

Calculating the integral, which is the second multiplier in (A.32), by parts, one obtains

∫ t o t exp ( t * − t o τ c ) ∂ ∂ t * exp [ 1 K c ∫ t o t ∗ ∂ Δ p ( t * ) ∂ t * ∂ t * ] ∂ t * = { exp ( t * − t o τ c ) exp [ 1 K c ∫ t o t ∗ ∂ Δ p ( t * ) ∂ t * ∂ t * ] } | t o t − ∫ t o t ∂ ∂ t * [ exp ( t * − t o τ c ) ] exp [ 1 K c ∫ t o t ∗ ∂ Δ p ( t * ) ∂ t * ∂ t * ] ∂ t * , t > t o

that is

∫ t o t exp ( t * − t o τ c ) ∂ ∂ t * exp [ 1 K c ∫ t o t ∗ ∂ Δ p ( t * ) ∂ t * ∂ t * ] ∂ t * = exp [ t − t o τ c + Δ p ( t ) − Δ p ( t o ) K c ] − 1 − 1 τ c ∫ t o t exp [ t * − t o τ c + Δ p ( t * ) − Δ p ( t o ) K c ] ∂ t * , t > t o (A.33)

Substitution of (A.33) into (A.32) results in (A.28). This proves assertion (i).

Inequalities (A.29) follow from the from of the right-hand side of (A.28) and notations (A.25), (A.24). Inequalities (A.30) follow from (A.29) and also notations (A.25), (A.24). This proves assertion (ii).

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 (A.26), notations (A.24), (A.25), and limit relation account (A.21), one obtains (A.31). This completes both the proofs of assertion (iv) and the theorem. ,

In view of (A.27), the infinite-time value ε c . ∞ (see (A.31)) presents the infinite-time residual value of the cell volume. In view of (A.22), one should expect that the limit in (A.31) is greater than unit and, thus,

ε c . ∞ > 0 (A.34)

i.e., the above residual value corresponds to the relative enlargement of a cell.

Remark A.2. The main result of the present section is the explicit expression (A.31) (see also (A.34)) for the infinite-time value ε c . ∞ of the cell-volume relative enlargement in terms of:

・ τ c , the relaxation time induced by the cell osmosis and described with (A.14);

・ K c , the bulk modulus of the cell;

・ Δ p ( t ) , the acoustic pressure (see (8.3)) in the SLT, which includes the cell under consideration.

Parameter K c and the parameters underlying Δ p ( t ) are mechanical. Concerning τ c , one can note that definition (A.14) of this parameter also includes mechanical parameters, namely V ¯ c , S ¯ c , k B , and T . However, the three other parameters in (A.14), namely

・ H c , the hydraulic conductance of the semi-permeable membrane of the cell;

・ σ c , the reflection coefficient of the cell membrane;

・ n ¯ c * , the equilibrium value the total concentration of the solutes in the ICF of the cell;

are biological. Thus, expression (A.31) provides an explicit picture of the combined, mechanobiological effect on the infinite-time value ε c . ∞ of the cell-volume relative enlargement. This picture indicates that ε c . ∞ is directly proportional to both the ratio Δ p ( t ) / K c (see (A.26) and τ c (cp., Remark A.1).

The developed mechnobiological picture for the cell enlargement can be summarized in terms of the answers to questions (1)-(3) formulated in the text above (A.1) as follows.

According to (A.31), an acoustic variable in an SLT that directly affects the cell enlargement is acoustic pressure Δ p in the SLT. This is the answer to quation (1).

Expression (A.31) also provides the answer to question (2). Indeed, the role of Δ p is specified with (A.31).

The corresponding answer to question (3) is the following. In order to reduce the cell 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. In particular, in the idealized case as this relaxation time tends to zero, the cell volume remains unchanged (see Remark A.1). ,

Design of future vibtrating hand-held tools often focuses on improvement of the vibration characteristics, for example, on reduction of Δ p ( t ) . However, these activities cannot be regarded successful until the effect of the improvement is estimated from either biopsy tests of the fingers of the persons who used the advanced tools (which is highly resource and time consuming) or, at least, from quantitative simulation in accordance with the summary formulated in Remark A.2 (which is much more affordable). ,

ANS average normal stress

HAVS hand-arm vibration syndrome

KV Kelvin-Voigt

LNS linear nonstationary

NSRF normalized stress-relaxation function

ODE ordinary differential equation

PDE partial differential equation

PIDE partial integro-differential equation

SLT soft living tissue

SRT stress-relaxation time

VWF vibration-induced white finger