Access the full text.
Sign up today, get DeepDyve free for 14 days.
Loading next page...
/lp/arxiv-cornell-university/constructing-separable-non-2-pi-periodic-solutions-to-the-navier-lam-e-QItTdITwg0
- ISSN
- 2070-0733
- eISSN
- ARCH-3337
- DOI
- 10.4208/aamm.OA-2019-0128
- Publisher site
- See Article on Publisher Site
Abstract
In a previous paper [Adv. Appl. Math. Mech. 10 (2018), pp. 1025-1056], we used the Buchwald representation to construct several families of separable cylindrical solutions to the Navier-Lam´e equation; these solutions had the property of being 2π-periodic in the circumferential coordinate. In this paper, we extend the analysis and obtain the complementary set of separable solutions whose circumferential parts are elementary 2π-aperiodic functions. Collectively, we construct eighteen distinct families of separable solutions; in each case, the circumferential part of the solution is one of three elementary 2π-aperiodic functions. These solutions are useful for solving a wide variety of dynamical problems that involve cylindrical geometries and for which 2π-periodicity in the angular coordinate is incompatible with the given boundary conditions. As illustrative examples, we show how the obtained solutions can be used to solve certain forced-vibration problems involving open cylindrical shells and open solid cylinders where (by virtue of the boundary conditions) 2π- periodicity in the angular coordinate is inappropriate. As an addendum to our prior work, we also include an illustrative example of a certain type of asymmetric problem that can be solved using the particular 2π-periodic subsolutions that ensue when there is no explicit dependence on the circumferential coordinate. arXiv:1906.11634v1 [math-ph] 25 Jun 2019 I. INTRODUCTION In linear elastodynamics [1], the displacement of a homogeneous, isotropic, and linearly- elastic solid is governed by the Navier-Lam´e (NL) equation. Over the last century, many specialized methods have been developed for obtaining solutions to the NL equation. One of the most popular and effective methods is the method of potentials [2–4], whereby the displacement vector is represented as a specific combination of one or more scalar and/or vector functions (called potentials). The scalar components of the NL equation are a set of coupled PDEs that cannot generally be solved in closed form. Representations in terms of displacement potentials may, once substituted into the NL equation, yield a system of PDEs that are uncoupled, less coupled, or at least simpler than the original component equations. Many different representations exist including the classical Helmholtz-Lam´e and Papkovich-Neuber representations. In applications, the analytical utility of a representation often depends on the choice of the working coordinate system. One representation that has proven to be analytically effective in anisotropic problems with cylindrical symmetry is the so-called Buchwald representation (see Refs. [5–7] and references therein). The Buch- wald representation is also applicable to problems involving isotropic media [8, 9], but its application to such problems is relatively rare [10–12]. The Buchwald representation involves three scalar potential functions (see Eq. (2) of Sec. II). It can be shown by various means (see, for example, Refs. [10, 13]) that the Buchwald representation reduces the original scalar components of the NL equation (which are a set of three coupled PDEs) to a set of two coupled PDEs involving two of the potentials and one separate decoupled PDE involving the remaining potential. By assuming separable product solutions for the three potentials and imposing certain conditions on their axial and temporal parts, it was shown in Ref. [13] that the coupled subsystem reduces to a homogeneous linear system that can be solved by elimination. The elimination procedure generates two independent fourth-order linear PDEs with constant coefficients, which can be subsequently solved using linear combinations of solutions to the two-dimensional Helmholtz equation in polar coordinates. A solution to the independent decoupled PDE can be obtained using separation of variables. In Ref. [13] (henceforth referred to as SC1), we specifically constructed particular solutions (to the PDEs governing the Buchwald potentials) possessing 2π-periodic angular parts. In the present paper, we consider the complementary class of 2 separable solutions whose angular parts are either aperiodic or periodic but not 2π-periodic. Although more involved, all three Buchwald potentials can again be completely determined (under the stipulated conditions), and following the approach of SC1 we can then construct eighteen distinct families of parametric solutions to the NL equation. The angular parts of these solutions (to the NL equation) are of course not 2π-periodic. One new complication that arises is that some solutions lacking 2π-periodicity in their angular parts are complex- valued, and in those cases, we discuss how to extract acceptable real-valued solutions. To our knowledge, the solutions obtained here constitute the first comprehensive set of 2π-aperiodic cylindrical solutions to the NL equation. Independent of the fact that they were derived using Buchwald potentials, these exact closed-form families of solutions are fundamentally important and useful in their own right. They can be directly applied to a fundamental set of linear-elastic boundary-value-problems in cylindrical coordinates where 2π-periodicity in the angular coordinate is incompatible with the given boundary conditions. In the second half of the paper, we provide three illustrative examples of such problems and detail how the obtained parametric solutions can be used to solve them. As an addendum to SC1, we also include an example that illustrates the use of a special set of 2π-periodic sub- solutions that results when there is no explicit functional dependence on the circumferential coordinate. II. EQUATIONS OF MOTION The NL equation can be written in vector form as [1, 4] ∂ u μ∇ u + (λ + μ)∇(∇ · u) + b = ρ , (1a) ∂t or (using the well-known identity for the vector Laplacian ) as ∂ u (λ + 2μ)∇(∇ · u) − μ∇ × (∇ × u) + b = ρ , (1b) ∂t where u is the displacement vector, λ and μ are the Lam´e constants from the classical linear theory of elasticity, b is the local body force, and ρ is the (constant) density of the material. In this paper, we consider the situation in which there are only surface forces acting on the 1 2 ∇ u = ∇(∇ · u) − ∇ × (∇ × u) 3 solid and no body forces (i.e., b = 0). In the circular cylindrical coordinate system, the displacement u ≡ u(r, θ, z, t) depends on the spatial coordinates (r, θ, z), which denote the radial, circumferential, and longitudinal coordinates, respectively, and on the time t. The (standard) Buchwald representation of the displacement field involves three scalar potentials Φ(r, θ, z, t), Ψ(r, θ, z, t), χ(r, θ, z, t), and is written as follows [5]: ∂Ψ ∂Φ u = ∇Φ + ∇ × (χˆz) + − ˆz, (2) ∂z ∂z where ˆz is the unit basis vector in the z-direction. It can be shown by various means [10, 13] that representation (2) satisfies Eq. (1) identically (assuming b = 0) provided the Buchwald potentials Φ, Ψ, and χ are solutions of the coupled PDE system: 2 2 2 ∂ Ψ ∂ Φ ∂ Φ (λ + 2μ)∇ Φ + (λ + μ) − (λ + μ) − ρ = 0, (3a) 2 2 2 ∂z ∂z ∂t 2 2 2 ∂ Φ ∂ Ψ ∂ Ψ 2 2 (λ + μ) ∇ Φ − + μ∇ Ψ + (λ + μ) − ρ = 0. (3b) 2 2 2 ∂z ∂z ∂t ∂ χ μ∇ χ − ρ = 0. (3c) ∂t In SC1 [13], we obtained, under certain conditions, particular separable solutions to subsystem (3a,3b) that are 2π-periodic in the circumferential coordinate θ. In the present paper, we seek to obtain the complementary set of separable solutions to subsystem (3a,3b) whose circumferential functions are 2π-aperiodic. By construction, the general method of solution will be guided by and closely follow the method developed in SC1. III. PARTICULAR SEPARABLE SOLUTIONS Assume that the Buchwald potentials Φ and Ψ have the separable form: Φ(r, θ, z, t) = φ (r, θ)φ (z)φ (t), (4a) ⊥ z t Ψ(r, θ, z, t) = ψ (r, θ)ψ (z)ψ (t), (4b) ⊥ z t which together with certain conditions on the functions {φ (z), φ (t), ψ (z), ψ (t)} (to be z t z t specified below) define the class of separable solutions to subsystem (3a,3b) that we seek to find in this paper. Using the notation 2 2 2 ∂ ∂ 1 ∂ 1 ∂ 2 2 ∇ ≡ ∇ − = + + (5) 2 2 2 2 ∂z ∂r r ∂r r ∂θ 4 to denote the two-dimensional Laplacian operator in polar coordinates, subsystem (3a,3b) can then be cast in the form: 2 2 2 ∂ Φ ∂ Ψ ∂ Φ (λ + 2μ)∇ Φ + μ + (λ + μ) − ρ = 0, (6a) 2 2 2 ∂z ∂z ∂t 2 2 ∂ Ψ ∂ Ψ 2 2 (λ + μ)∇ Φ + μ∇ Ψ + (λ + 2μ) − ρ = 0. (6b) ⊥ ⊥ 2 2 ∂z ∂t Substituting (4) into system (6) leads to the system 2 ′′ ′′ ′′ (λ + 2μ)φ φ ∇ + μφ φ − ρφ φ φ + [(λ + μ)ψ ψ ] ψ = 0, (7a) z t t z ⊥ t ⊥ ⊥ z t z 2 2 ′′ ′′ (λ + μ)φ φ ∇ φ + μψ ψ ∇ + (λ + 2μ)ψ ψ − ρψ ψ ψ = 0. (7b) z t ⊥ z t t z ⊥ ⊥ ⊥ z t A special class of solutions can be obtained by assuming common functional dependences in the axial and temporal parts of both potentials, that is, φ (z) = ψ (z), (8a) z z φ (t) = ψ (t), (8b) t t and by further assuming that the derivatives of the above functions are such that ′′ φ (z) = κφ (z), (8c) ′′ ψ (z) = κψ (z), (8d) ′′ φ (t) = τφ (t), (8e) ′′ ψ (t) = τψ (t), (8f) where κ ∈ R and τ ∈ R\{0} are arbitrary constants that we shall refer to as the axial and time constants, respectively. Ultimately, they will appear as parameters in the final solution. With these assumptions system (7) reduces to (λ + 2μ)∇ + μκ − ρτ φ + [(λ + μ)κ] ψ = 0, (9a) ⊥ ⊥ 2 2 (λ + μ)∇ φ + μ∇ + (λ + 2μ)κ − ρτ ψ = 0, (9b) ⊥ ⊥ ⊥ ⊥ which has the form of a homogeneous linear system: L φ + L ψ = 0, (10a) 1 ⊥ 2 ⊥ 5 L φ + L ψ = 0, (10b) 3 ⊥ 4 ⊥ where the linear operators {L , L , L , L } 1 2 3 4 L = (λ + 2μ)∇ + μκ − ρτ, (10c) L = (λ + μ)κ, (10d) L = (λ + μ)∇ , (10e) L = μ∇ + (λ + 2μ)κ − ρτ, (10f) are each first-order polynomial differential operators in ∇ with constant coefficients. Two distinct cases can be distinguished at this point depending on whether κ 6= 0 or κ = 0. The former is the general case while the latter is a special case for which the linear system (10) decouples. We shall consider these two cases separately. A. General Case: κ 6= 0 Under conditions (8) and assuming κ 6= 0, we can solve for φ (r, θ) and ψ (r, θ) using ⊥ ⊥ the method of elimination. It can be shown that φ (r, θ) and ψ (r, θ) independently satisfy ⊥ ⊥ the decoupled equations: (L L − L L )φ = 0, (11a) 1 4 2 3 ⊥ (L L − L L )ψ = 0, (11b) 1 4 2 3 ⊥ where (L L − L L ) is the operational determinant of system (10). By inspection, this 1 4 2 3 determinant is non-zero, and thus a unique solution to system (10) exists. We are free to solve Eq. (11a) to obtain a solution for φ and then use either of Eqs. (10a) or (10b) to solve for ψ , or alternatively, we can first solve Eq. (11b) for ψ and then use either of Eqs. (10a) ⊥ ⊥ or (10b) to solve for φ . For convenience, we choose the first option. After some straightforward algebra, it can be shown that Eq. (11a) assumes the form 4 2 a ∇ + a ∇ + a φ = 0, (12a) 2 1 0 ⊥ ⊥ ⊥ where a = μ(λ + 2μ), (12b) 6 a = 2μ(λ + 2μ)κ − (λ + 3μ)ρτ, (12c) a = (μκ − ρτ) [(λ + 2μ)κ − ρτ] . (12d) Since the operator in Eq. (12a) is a polynomial differential operator in ∇ with constant coefficients, we may look for separable solutions (denoted here by ϕ) of the same fundamental form as those of the two-dimensional Helmholtz equation in polar coordinates (∇ ϕ = Λϕ). This leads to the characteristic equation a Λ + a Λ + a = 0, (13) 2 1 0 whose roots, obtained from the quadratic formula −a ± a − 4a a 1 2 0 Λ = , (14a) 2a are ρτ ρτ Λ = − κ − , Λ = − κ − . (14b) − + (λ + 2μ) μ If ϕ (r, θ) and ϕ (r, θ) denote two independent separable solutions to Eq. (12) corre- − + sponding to the roots Λ and Λ , respectively, then φ (r, θ) = ϕ (r, θ)+ϕ (r, θ). We shall, − + ⊥ − + (1) (2) for notational convenience, let φ (r, θ) ≡ ϕ (r, θ) and φ (r, θ) ≡ ϕ (r, θ). In this more − + ⊥ ⊥ intuitive notation, 2 2 X X (s) (s) (s) φ (r, θ) = (r, θ) = (θ). (15) φ φ (r)φ ⊥ r θ s=1 s=1 (1) (2) The angular parts of ϕ (r, θ) and ϕ (r, θ) [i.e., φ (θ) and φ (θ) in Eq. (15), respectively] − + θ θ can be either periodic or non-periodic functions of θ. In SC1 [13], we considered exclusively the class of separable solutions to system (10) that possess 2π-periodic angular parts. In the present paper, we shall consider the complementary class of separable solutions to system (10) whose angular parts are either aperiodic or periodic but not 2π-periodic. Such solutions are useful for solving a wide variety of dynamical problems involving open cylindrical shells and panels as well as open solid cylinders, where 2π-periodicity in the angular coordinate is neither necessary nor (depending on the boundary conditions) appropriate [14, 15]. Utilizing the fundamental solutions to the 2D Helmholtz equation in polar coordinates (see Appendix A) with the Helmholtz constant taking the values Λ ≡ Λ and Λ ≡ Λ , we 1 − 2 + obtain the following particular solution to Eq. (12): (s) (s) φ (r, θ) = φ (r)φ (θ), (16a) s=1 7 where p p C exp − |η|θ + D exp |η|θ if η < 0 s s (s) φ (θ) = , (16b) C + D θ if η = 0 θ s s √ √ C cos ηθ + D sin ηθ if η > 0 s s p p √ √ A J |Λ |r + B Y |Λ |r if Λ < 0 and η < 0 s s s s s i |η| i |η| p p A cos |η| ln r + B sin |η| ln r if Λ = 0 and η < 0 s s s √ √ √ √ A I Λ r + B K Λ r if Λ > 0 and η < 0 s s s s s i |η| i |η| p p A J |Λ |r + B Y |Λ |r if Λ < 0 and η = 0 s 0 s s 0 s s (s) φ (r) = , (16c) A + B ln r if Λ = 0 and η = 0 s s s √ √ A I Λ r + B K Λ r if Λ > 0 and η = 0 s 0 s s 0 s s p p √ √ A J |Λ |r + B Y |Λ |r if Λ < 0 and η > 0 s η s s η s s √ √ η − η if Λ = 0 and η > 0 A r + B r s s s √ √ √ √ A I Λ Λ r + B K r if Λ > 0 and η > 0 s η s s η s s and {A , A , B , B , C , C , D , D } are arbitrary constants. Note that the above solution 1 2 1 2 1 2 1 2 involves three parameters: κ ∈ R\{0} and τ ∈ R\{0}, which occur in the (non-arbitrary) constants {Λ , Λ }, and η ∈ R. The parameter η in solution (16) is analogous to the 1 2 constant n ∈ N that appears in the solution to Eq. (12) when the circumferential functions n o (s) φ (θ) : s = 1, 2 are 2π-periodic [13]. Indeed, it can be readily verified that solution (16) reduces to the analogous 2π-periodic solution obtained in SC1 [13] when η ∈ {n : n ∈ N} with the caveat that D = 0 (s = 1, 2) when n = 0. Henceforth, we shall assume η ∈ 2 + R\ {N : N ∈ Z }. It is important to recognize that, when the parameter η < 0, the radial functions n o (s) φ (r) : s = 1, 2 are generally complex-valued functions. Real particular solutions can (s) however be extracted in these cases. When Λ > 0 and η < 0, φ (r) is a solution to (s) (s) d φ (r) dφ (r) r r 2 2 (s) r + r + |η| − Λ r φ (r) = 0, (17) dr dr and based on the discussion in the first paragraph of Appendix B, the real-valued functions ( ) n o p p √ √ Re I Λ r , K Λ r (18) s s i |η| i |η| 8 constitute a linearly independent set of solutions to Eq. (17). When Λ < 0 and η < 0, (s) φ (r) is a solution to (s) (s) d φ (r) dφ (r) r r 2 2 (s) r + r + |Λ |r + |η| φ (r) = 0, (19) dr dr and in this case (c.f., the discussion in the second paragraph of Appendix B), the functions n o p p π |η| ¯√ √ J |Λ |r ≡ sech Re J |Λ |r , (20a) s s |η| i |η| n o p p π |η| √ √ Y |Λ |r ≡ sech Re Y |Λ |r , (20b) s s |η| i |η| together constitute a real and linearly independent set of solutions to Eq. (19). To obtain the solution for ψ (r, θ), it is simplest to re-arrange Eq. (10a) and (using the solution for φ (r, θ)) directly solve for ψ (r, θ): ⊥ ⊥ [(λ + 2μ)∇ + (μκ − ρτ)] ψ (r, θ) = − φ (r, θ). (21) ⊥ ⊥ (λ + μ)κ In SC1 [13], it is shown that, given the solution for φ (r, θ) as expressed in Eq. (15), the solution for ψ (r, θ) can be expressed as: (s) (s) ψ (r, θ) = γ φ (r)φ (θ), (22a) ⊥ s r θ s=1 where the constant γ (s = 1, 2) is given by 1 if s = 1 γ = . (22b) 1 ρτ κ − if s = 2 κ μ n o (s) The solution for ψ (r, θ) is thus given by Eq. (22), where the functions φ (r) : s = 1, 2 ⊥ r n o (s) and φ (θ) : s = 1, 2 are now given by Eqs. (16c) and (16b), respectively. The only remaining task in solving system (6) is to determine the axial and temporal parts of the full potentials Φ and Ψ [see Eq. (4)]. The conditions listed in (8) provide the governing equations for determining these functions. Conditions (8a), (8c), and (8d) yield p p E cos |κ|z + F sin |κ|z if κ < 0 φ (z) = ψ (z) = , (23) z z √ √ E exp (− κz) + F exp ( κz) if κ > 0 9 and similarly, conditions (8b), (8e), and (8f) yield p p G cos |τ|t + H sin |τ|t if τ < 0 φ (t) = ψ (t) = , (24) t t √ √ G exp (− τt) + H exp ( τt) if τ > 0 where {E, F, G, H} are arbitrary constants. Particular separable solutions to subsystem (3c) can be obtained by assuming a product solution of the form χ(r, θ, z, t) = χ (r)χ (θ)χ (z)χ (t) (25) r θ z t and using separation of variables, as shown in Appendix C. The separation constants {υ , υ , υ , υ } defined in Appendix C are a set of independent parameters in their own t z θ r right, limited only by relation (C5). Thus, in theory, the full solution to system (3) can involve six independent parameters (e.g., {κ, τ, η, υ , υ , υ }). The number of parameters in t z θ the full solution can however be systematically reduced (yielding more compact particular solutions) by fixing one or more of the separation constants at prescribed value(s) that are linear combinations of the parameters {κ, τ, η}. Although not necessary, it is convenient (and often useful) to eliminate the separation constants {υ , υ , υ , υ } entirely by fixing t z θ r their values according to the following prescription: ρτ ρτ υ = , υ = κ, υ = η, υ = υ − υ = − κ ≡ Λ . (26) t z θ r t z 2 μ μ Upon doing so, one can immediately conclude (from inspection of the solutions given in Appendix C) that, apart from arbitrary constants, χ (z) and χ (t) are identical to φ (z) z t and φ (t), respectively, and that χ (r, θ) ≡ χ (r)χ (θ) is (apart from arbitrary constants) ⊥ r θ identical to the s = 2 term of Eq. (16). Explicitly: p p √ √ A J |Λ |r + B Y |Λ |r if Λ < 0 and η < 0 3 2 3 2 2 i |η| i |η| p p A cos |η| ln r + B sin |η| ln r if Λ = 0 and η < 0 3 3 2 √ √ √ √ A I Λ r + B K Λ r if Λ > 0 and η < 0 3 2 3 2 2 i |η| i |η| p p A J |Λ |r + B Y |Λ |r if Λ < 0 and η = 0 3 0 2 3 0 2 2 χ (r) = , (27a) A + B ln r if Λ = 0 and η = 0 3 3 2 √ √ A I Λ r + B K Λ r if Λ > 0 and η = 0 3 0 2 3 0 2 2 p p √ √ A J |Λ |r + B Y |Λ |r if Λ < 0 and η > 0 3 η 2 3 η 2 2 √ √ η − η A r + B r if Λ = 0 and η > 0 3 3 2 √ √ √ √ A I Λ r + B K Λ r if Λ > 0 and η > 0 3 η 2 3 η 2 2 10 p p C exp − |η|θ + D exp |η|θ if η < 0 3 3 χ (θ) = , (27b) C + D θ if η = 0 θ 3 3 √ √ C cos ηθ + D sin ηθ if η > 0 3 3 p p e e E cos |κ|z + F sin |κ|z if κ < 0 χ (z) = , (27c) z √ √ e e E exp (− κz) + F exp ( κz) if κ > 0 and p p e e G cos |τ|t + H sin |τ|t if τ < 0 χ (t) = , (27d) √ √ e e G exp (− τt) + H exp ( τt) if τ > 0 n o e e e e where A , B , C , D , E, F, G, H are arbitrary constants. The function χ (θ) given by 3 3 3 3 2 + Eq. (27b) is 2π-periodic when η ∈ {N : N ∈ Z } and in the special case η = D = 0. Since all three Buchwald potentials have now been completely determined (under con- ditions (8) and (26)), we can write down the corresponding particular solutions to Eq. (1) using the component form of Eq. (2): ∂Φ 1 ∂χ 1 ∂Φ ∂χ ∂Ψ ˆ ˆ u = u ˆr + u θ + u ˆz = + ˆr + − θ + ˆz. (28) r θ z ∂r r ∂θ r ∂θ ∂r ∂z Using Eq. (28), the cylindrical components of the displacement field are thus given by: (s) dφ (r) 1 dχ (θ) (s) u = φ (θ) φ (z)φ (t) + χ (r) χ (z)χ (t), (29a) r z t r z t dr r dθ s=1 2 (s) 1 dφ (θ) dχ (r) (s) r u = φ (r) φ (z)φ (t) − χ (θ)χ (z)χ (t), (29b) θ z t r θ z t r dθ dr s=1 dψ (z) (s) (s) u = γ φ (r)φ (θ) ψ (t), (29c) z s t r θ dz s=1 where: (s) (s) (i) φ (r) and φ (θ) in Eqs. (29a)-(29c) are given by Eqs. (16c) and (16b), respectively; (ii) φ (z) and φ (t) in Eqs. (29a)-(29b) are given by Eqs. (23) and (24), respectively; z t (iii) χ (r) and χ (θ) in Eqs. (29a)-(29b) are given by Eqs. (27a) and (27b), respectively; r θ (iv) χ (z) and χ (t) in Eqs. (29a)-(29b) are given by Eqs. (27c) and (27d), respectively; z t (v) ψ (z) and ψ (t) in Eq. (29c) are given by Eqs. (23) and (24), respectively; z t (vi) the constant γ in Eq. (29c) is given by Eq. (22b). 11 Note: Items (iii) and (iv) above are valid only under prescription (26). If and when prescription (26) is inapplicable (as in the example problem considered in Sec. V), then the functions {χ (r), χ (θ), χ (z), χ (t)} in Eqs. (29a)-(29b) are given by Eqs. (C9)-(C12) of r θ z t Appendix C. 1. Special 2π-Periodic Subsolution: No Dependence on the Circumferential Coordinate θ n o (s) When the circumferential functions φ (θ) : s = 1, 2 and χ (θ) are constant functions, θ θ the Buchwald potentials (and the resulting displacement field) are independent of the θ coordinate, and Eqs. (29a)-(29c) reduce to the following general form: (s) dφ (r) u = C φ (z)φ (t), (30a) r s z t dr s=1 dχ (r) u = −C χ (z)χ (t), (30b) θ 3 z t dr dψ (z) (s) u = γ C φ (r) ψ (t). (30c) z s s t dz s=1 In this case, the Buchwald potentials (and the resulting displacement field) are trivially 2π-periodic. It should however be emphasized that the absence of the θ coordinate in the solution does not mean that the displacement field is necessarily axisymmetric. Axisymmet- ric displacement fields must, in addition, have no circumferential component (i.e., u = 0). An example of an asymmetric problem having no explicit dependence on the θ coordinate is given in Sec. V. B. Special Case: κ = 0 When the parameter κ = 0, system (9) decouples and reduces to: (λ + 2μ)∇ − ρτ φ = 0, (31a) 2 2 (λ + μ)∇ φ + μ∇ − ρτ ψ = 0. (31b) ⊥ ⊥ ⊥ ⊥ Equation (31a) can be cast in the form of a 2D Helmholtz equation (in polar coordinates) (0) ρτ and immediately solved using separation of variables. Letting Λ ≡ , the following (λ+2μ) separable solution to Eq. (31a) is then obtained using the results of Appendix A: φ (r, θ) = R (r)Θ (θ), (32a) ⊥ 1 1 12 where p p C exp − |η|θ + D exp |η|θ if η < 0 1 1 Θ (θ) = , (32b) C + D θ if η = 0 1 1 √ √ C cos ηθ + D sin ηθ if η > 0 1 1 r r (0) (0) √ √ A J Λ r + B Y Λ r if τ < 0 and η < 0 1 1 1 1 i |η| i |η| q q (0) (0) √ √ A I Λ r + B K Λ r if τ > 0 and η < 0 1 1 1 1 i |η| i |η| r r (0) (0) A J Λ r + B Y Λ r if τ < 0 and η = 0 1 0 1 0 1 1 R (r) = q q , (32c) (0) (0) A I Λ r + B K Λ r if τ > 0 and η = 0 1 0 1 0 1 1 r r (0) (0) √ √ A J Λ r + B Y Λ r if τ < 0 and η > 0 1 η 1 η 1 1 q q (0) (0) A I Λ r + B K Λ r if τ > 0 and η > 0 1 η 1 η 1 1 where {A , B , C , D } are arbitrary constants. 1 1 1 1 2 2 To solve for ψ (r, θ), we substitute (λ + μ)∇ φ = − [μ∇ − ρτ] φ (which is obtained ⊥ ⊥ ⊥ ⊥ ⊥ from rearranging Eq. (31a)) into Eq. (31b). This yields the PDE: μ∇ − ρτ W = 0, where W ≡ (ψ − φ ) . (33) ⊥ ⊥ ⊥ ⊥ As before, Eq. (33) can be cast in the form of a 2D Helmholtz equation in polar coordinates. (0) ρτ Letting Λ ≡ , the following separable solution to Eq. (33) is then obtained using the results of Appendix A: W (r, θ) = R (r)Θ (θ), (34a) ⊥ 2 2 where p p C exp − |η|θ + D exp |η|θ if η < 0 2 2 Θ (θ) = , (34b) 2 C + D θ if η = 0 2 2 √ √ C cos ηθ + D sin ηθ if η > 0 2 2 13 r r (0) (0) √ √ A J Λ r + B Y Λ r if τ < 0 and η < 0 2 2 2 2 i |η| i |η| q q (0) (0) √ √ A I Λ r + B K Λ r if τ > 0 and η < 0 2 2 2 2 i |η| i |η| r r (0) (0) A J Λ r + B Y Λ r if τ < 0 and η = 0 2 0 2 0 2 2 R (r) = q , (34c) (0) (0) A I Λ r + B K Λ r if τ > 0 and η = 0 2 0 2 0 2 2 r r (0) (0) A J Λ r + B Y Λ r if τ < 0 and η > 0 2 η 2 η 2 2 q q (0) (0) √ √ A I Λ r + B K Λ r if τ > 0 and η > 0 2 η 2 η 2 2 where {A , B , C , D } are arbitrary constants. Since W ≡ (ψ − φ ), the solution for 2 2 2 2 ⊥ ⊥ ⊥ ψ (r, θ) is then the sum of the two solutions given by Eqs. (32) and (34): ψ (r, θ) = R (r)Θ (θ). (35) ⊥ s s s=1 The axial parts of the full potentials Φ and Ψ [c.f., Eqs. (4)], as determined from conditions (8a), (8c), and (8d) are: φ (z) = ψ (z) = E + Fz, (36) z z where {E, F} are arbitrary constants. The temporal parts of potentials Φ and Ψ are un- changed [c.f., Eq. (24)]. All that remains is to determine the Buchwald potential χ, which is governed by Eq. (3c), and by assumption, has the separable form (25). Under conditions (26), one can immediately conclude from inspection of the solutions given in Appendix C that, apart from arbitrary con- stants, the components χ (z) and χ (t) are identical to φ (z) and φ (t) as given by Eqs. (36) z t z t and (24), respectively, and that χ (r, θ) ≡ χ (r)χ (θ) is (apart from arbitrary constants) ⊥ r θ identical to the function W (r, θ) given by Eq. (34). Thus, analogous to Eqs. (27a)-(27c), the spatial components of χ are as follows: p p C exp − |η|θ + D exp |η|θ if η < 0 3 3 χ (θ) = , (37a) C + D θ if η = 0 3 3 √ √ C cos ηθ + D sin ηθ if η > 0 3 3 14 r r (0) (0) √ √ A J Λ r + B Y Λ r if τ < 0 and η < 0 3 3 2 2 i |η| i |η| q q (0) (0) √ √ A I Λ r + B K Λ r if τ > 0 and η < 0 3 3 2 2 i |η| i |η| r r (0) (0) A J Λ r + B Y Λ r if τ < 0 and η = 0 3 0 3 0 2 2 χ (r) = q , (37b) (0) (0) A I Λ r + B K Λ r if τ > 0 and η = 0 3 0 3 0 2 2 r r (0) (0) A J Λ r + B Y Λ r if τ < 0 and η > 0 3 η 3 η 2 2 q q (0) (0) √ √ A I Λ r + B K Λ r if τ > 0 and η > 0 3 η 3 η 2 2 e e χ (z) = E + Fz, (37c) n o e e where A , B , C , D , E, F are arbitrary constants. The temporal part of χ (i.e., χ (t)) is 3 3 3 3 again given by Eq. (27d). With all three Buchwald potentials now completely determined (under the stipulated conditions), we can write down the corresponding particular solutions to Eq. (1) using (28): dR (r) 1 dχ (θ) u = Θ (θ)φ (z)φ (t) + χ (r) χ (z)χ (t), (38a) r 1 z t r z t dr r dθ 1 dΘ (θ) dχ (r) u = R (r) φ (z)φ (t) − χ (θ)χ (z)χ (t), (38b) θ 1 z t θ z t r dθ dr dψ (z) u = R (r)Θ (θ) ψ (t), (38c) z s s t dz s=1 where: (i) R (r) and R (r) in Eqs. (38a)-(38c) are given by Eqs. (32c) and (34c), respectively; 1 2 (ii) Θ (θ) and Θ (θ) in Eqs. (38a)-(38c) are given by Eqs. (32b) and (34b), respectively; 1 2 (iii) φ (z) and φ (t) in Eqs. (38a)-(38b) are given by Eqs. (36) and (24), respectively; z t (iv) ψ (z) and ψ (t) in Eq. (38c) are given by Eqs. (36) and (24), respectively; z t (v) χ (r), χ (θ), χ (z), and χ (t) in Eqs. (38a)-(38b) are given by (37b), (37a), (37c), and r θ z t (27d), respectively. Note 1: Solution (38) is valid only when the parameter κ defined in Eqs. (8c)-(8d) is zero. Note 2: Item (v) above is valid only under prescription (26). If and when prescription (26) is inapplicable, then the functions {χ (r), χ (θ), χ (z), χ (t)} in Eqs. (38a)-(38c) are given r θ z t by Eqs. (C9)-(C12) of Appendix C. 15 IV. APPLICATIONS Our focus thus far has been on solving the PDE system (3) under assumptions (4) and (25) and conditions (8) and (26) and then using (28) to arrive at solutions (29) and (38). Independent of how these solutions were actually derived, we have at our disposal a set of particular solutions to the NL equation that can be directly applied to a fundamental set of linear-elastic boundary-value problems (BVPs) in cylindrical coordinates, namely, problems wherein the displacement, stress, and strain fields are: constant, linear, sinusoidal, expo- nential, or hyperbolic in the circumferential and longitudinal coordinates , and sinusoidal or exponential in the time coordinate. Such solutions are useful for solving dynamical problems involving cylindrical geometries where (rational) 2π-periodicity in the angular coordinate is incompatible with the given boundary conditions. In the remaining sections of the paper, we shall illustrate using simple (but mathematically non-trivial) examples how solutions (29) and (38) can be used to solve such BVPs. As an addendum to SC1, we first however present an example that illustrates the use of the special 2π-periodic solution (30). This example is conceptually useful since it also furnishes an explicit example of a situation where pre- scription (26) is inapplicable and one must use the more general solution for the Buchwald potential χ given in Appendix C. V. APPLICATION 1: FORCED VIBRATION OF A CLOSED SOLID CYLINDER Consider a closed solid elastic circular cylinder of finite length L and radius R subjected to the following harmonic stresses on its curved surface: kπ σ (R, θ, z, t) = A sin z sin(ωt), (39a) rr mπ σ (R, θ, z, t) = B sin z , (39b) rθ kπ σ (R, θ, z, t) = C cos z sin(ωt), (39c) rz where σ (r, θ, z, t) is a normal component of stress, σ (r, θ, z, t) and σ (r, θ, z, t) are shear rr rθ rz components of stress, {A, B, C} are prescribed constant stresses, {k, m} ∈ Z are prescribed dimensionless constants, and ω is the prescribed angular frequency of excitation. Note that Note that sinusoidal circumferential variations need not be 2π-periodic. 16 the circumferential stress given by Eq. (39b) is time independent. The cylinder is simply- supported at its flat ends, which are situated at z = 0 and z = L; the corresponding boundary conditions at the flat ends of the cylinder are thus: u (r, θ, 0, t) = u (r, θ, L, t) = 0, (40a) r r u (r, θ, 0, t) = u (r, θ, L, t) = 0, (40b) θ θ σ (r, θ, 0, t) = σ (r, θ, L, t) = 0, (40c) zz zz where σ (r, θ, z, t) is the normal component of stress along the axis of the cylinder. Condi- zz tions (39) are pure stress boundary conditions on the curved surface of the cylinder and must be satisfied for all θ ∈ [0, 2π], z ∈ (0, L), and arbitrary t. Conditions (40) are admissible mixed boundary conditions specifying orthogonal components of the stress and displacement at the flat ends of the cylinder; they must be satisfied for all r ∈ [0, R], θ ∈ [0, 2π], and arbitrary t. The condition that the displacement field is finite everywhere in the cylinder is implicit. For forced motion, at least one of {A, C} must be non-zero. (When A = C = 0, the cylin- der does not vibrate.) Since the stresses are independent of the circumferential coordinate θ in this problem, the displacement field of the cylinder must be so as well. Note however that when all of the stress constants {A, B, C} are non-zero, the motion is neither torsional (since both σ and σ are non-zero) nor axisymmetric (since σ is non-zero). rr rz rθ λ+2μ kπ Problem S: In the special case where the angular excitation frequency ω = , ρ L determine the displacement field at all points of the simply-supported cylinder subject to the non-uniform conditions (39) on its curved surface. The requirement of 2π-periodicity (since the cylinder is closed) and the fact that the displacement field is independent of θ (as deduced above) call upon the use of the special solution (30). Anticipating that the displacement component u will be time independent (due to the time independence of the stress component σ ) and that the displacement rθ components u and u must be time dependent implies that χ (t) in Eq. (30b) will not r z be of the same functional form as φ (t) and ψ (t) in Eqs. (30a) and (30c), respectively. t t Furthermore, unless the constants k and m in (39) are equal (which they need not be), χ (z) 17 must be of a different functional form than φ (z) and ψ (z). In such a situation, prescription z z (26) and the solution for the Buchwald potential χ that results from it [i.e., Eqs. (27a)-(27d)] are not valid and one must use the more general solution for χ given in Appendix C. The parameters associated with the solution for χ (i.e., the separation constants {υ , υ , υ , υ }) t z θ r must therefore be prescribed independently of the parameters associated with the solutions for Φ and Ψ (i.e., {κ, τ, η}). To obtain the general solutions for the u and u displacement components, fix the asso- r z ciated parameters κ, τ, and η as follows: kπ κ = − , τ = −ω , η = 0. (41) Then, choosing E = 0 and F = 1 in Eq. (23) and G = 0 and H = 1 in Eq. (24) yields: kπ φ (z) = ψ (z) = sin z , (42) z z φ (t) = ψ (t) = sin(ωt). (43) t t Choosing C = C = 1 in Eq. (16b) yields (noting that D = D = η = 0 is required for 1 2 1 2 2π-periodicity): (1) (2) φ (θ) = φ (θ) = 1. (44) θ θ The constants Λ and Λ , as determined from Eq. (14b), are: 1 2 2 2 2 2 kπ ρω kπ ρω Λ = − = 0, Λ = − < 0. (45) 1 2 L (λ + 2μ) L μ Since η = 0, Λ = 0, and Λ < 0, it then follows from Eq. (16c) that 1 2 (1) (2) φ (r) = A , φ (r) = A J |Λ |r , (46) 1 2 0 2 r r (1) where the ln r term generally present in the function φ (r) and the Y ( |Λ |r) term gen- 0 2 (2) erally present in the function φ (r) have been discarded (by prescribing B = B = 0 in 1 2 Eq. (16c)) in order to satisfy the finiteness condition at r = 0. To obtain the general solution for u , fix the parameters υ , υ , and υ associated with θ z t θ the solution for χ as follows: mπ υ = − , υ = 0, υ = 0. (47) z t θ The remaining parameter υ is then fixed by constraint (C5): mπ υ = υ − υ = > 0. (48) r t z 18 e e e e Then, choosing E = 0 and F = 1 in Eq. (C10) and G = 1 and H = 0 in Eq. (C11) yields: mπ χ (z) = sin z , (49) χ (t) = 1. (50) Choosing C = 1 in Eq. (C9) yields (noting that υ = D = 0 is required for 2π-periodicity): 3 θ 3 χ (θ) = 1. (51) Since υ = 0 and υ > 0, it then follows from Eq. (C12) that θ r mπ χ (r) = A I r , (52) 3 0 mπ where the K ( r) term has been discarded (by prescribing B = 0 in Eq. (C12)) in order 0 3 to again satisfy the finiteness condition. Inserting results (42)-(46) and (49)-(52) into (30) then yields the (general) displacement components: kπ u (r, z, t) = −A αJ (αr) sin z sin(ωt), (53a) r 2 1 mπ mπ u (r, z, t) = A I r sin z , (53b) θ 3 0 L L kπ kπ u (r, z, t) = [γ A + γ A J (αr)] cos z sin(ωt), (53c) z 1 1 2 2 0 L L where α ≡ |Λ |, and the constants γ (s = 1, 2) given by Eq. (22b) reduce (in this special 2 s case) to: 1 if s = 1 γ = . (53d) (λ+2μ) 1 − if s = 2 To complete the solution, we need only to determine the values of the constants {A , A , A } in Eqs. (53a)-(53c) that satisfy boundary conditions (39). Substituting 1 2 3 Eqs. (53) into the stress-displacement relations [Eqs. (D1) in Appendix D] yields (after some algebra) the required stress components as follows : ( ) kπ 1 kπ σ (r, z, t) = − A λ + 2A μα αJ (αr) − J (αr) sin z sin(ωt), (54a) rr 1 2 0 1 L r L Although not shown here, it can be verified that the stress component σ (r, z, t) automatically satisfies zz boundary condition (40c). mπ mπ 1 mπ mπ σ (r, z, t) = A μ I r − I r sin z , (54b) rθ 3 1 0 L L r L L kπ kπ σ (r, z, t) = A λ αJ (αr) cos z sin(ωt). (54c) rz 2 1 L L Application of boundary conditions (39) then yields the following 3 × 3 linear system: kπ ′ −λ −2μαJ (αR) 0 A A mπ mπ 1 mπ = , (55) 0 0 μ I R − I R A B 1 0 2 L L R L kπ 0 λ αJ (αR) 0 A C 1 3 where J (αR) ≡ αJ (αR) − J (αR) . A unique solution to system (55) exists only when 0 1 the following three conditions are satisfied: mπ mπ mπ λ 6= 0, J (αR) 6= 0, R I R 6= I R . (56) 1 1 0 L L L Otherwise the system is underdetermined and there exist an infinite number of solutions. Assuming conditions (56) are met, system (55) can be readily solved yielding the values: 1 2μJ (αR) A = − A + C , (57a) 2 kπ kπ λ J (αR) λ 1 A = , (57b) kπ λ αJ (αR) A = . (57c) mπ mπ 1 mπ μ I R − I R 1 0 L L R L The displacement field is therefore given by (53) with the constants {A , A , A } given by 1 2 3 (57). It is worthwhile to reiterate that when all three stress constants {A, B, C} are non- zero, the motion is neither torsional nor axisymmetric; the motion is axisymmetric only when B = 0. VI. APPLICATION 2: FORCED VIBRATION OF AN OPEN THICK CYLINDRI- CAL SHELL Consider an (elastic) open circular cylindrical shell whose domain is given by: R ≤ r ≤ R , 0 ≤ θ ≤ θ ≤ θ < 2π, and 0 ≤ z ≤ L. The geometrical parameters {L, R , R } are all 2 1 2 1 2 finite, but otherwise arbitrary. The flat ends of the shell, which are situated at z = 0 and z = L, are clamped: u (r, θ, 0, t) = u (r, θ, L, t) = 0, (58a) r r 20 u (r, θ, 0, t) = u (r, θ, L, t) = 0, (58b) θ θ u (r, θ, 0, t) = u (r, θ, L, t) = 0. (58c) z z Conditions (58a)-(58c) are pure displacement boundary conditions and must be satisfied for all r ∈ [R , R ], θ ∈ [θ , θ ], and arbitrary t. 1 2 1 2 In the following subsections, we shall specify several different sets of boundary condi- tions on the radial and circumferential faces of the shell and then obtain solutions to the corresponding boundary-value problems. A. Linear Circumferential Variations Suppose the (radial) curved surfaces of the shell are subjected to the following non- uniform stresses: kπ σ (R , θ, z, t) = (A + B θ) sin z sin(ωt), (59a) rr i i i kπ σ (R , θ, z, t) = C sin z sin(ωt), (59b) rθ i i kπ σ (R , θ, z, t) = (F + G θ) cos z sin(ωt), (59c) rz i i i where the index i = {1, 2}, {A , B , C , F , G : i = 1, 2} are prescribed constant stresses, i i i i i k ∈ Z is a prescribed constant, and ω is the prescribed angular excitation frequency. Conditions (59a)-(59c) are pure stress boundary conditions and must be satisfied for all θ ∈ [θ , θ ], z ∈ (0, L), and arbitrary t. Suppose furthermore that the following mixed 1 2 (non-uniform) conditions are prescribed on the circumferential faces of the shell: 1 kπ u (r, θ , z, t) = U R sin z sin(ωt), (60a) r i i r L u (r, θ , z, t) = 0, (60b) z i 1 kπ σ (r, θ , z, t) = S R sin z sin(ωt), (60c) θθ i i r L where again i = {1, 2}, {U : i = 1, 2} are prescribed constant displacements (U 6= U ), i 1 2 {S : i = 1, 2} are prescribed constant stresses (S 6= S ), and R ≡ (R + R )/2 is the i 1 2 1 2 mean radius of the shell. Conditions (60a)-(60c) are admissible mixed boundary conditions specifying orthogonal components of the stress and displacement on the circumferential faces of the shell and must be satisfied for all r ∈ [R , R ], z ∈ (0, L), and arbitrary t. 1 2 21 q μ kπ Problem A: In the special case where the excitation frequency ω = , determine ρ L the displacement field at all points of the clamped shell subject to conditions (59) on its radial surfaces and to conditions (60) on its circumferential faces. To obtain general solutions for {u , u , u } compatible with the given boundary condi- r θ z tions, fix the parameters κ, τ, and η (associated with the Buchwald potentials Φ and Ψ) as in (41) and employ prescription (26) for Buchwald potential χ. Then, choosing E = 0 and e e e F = 1 in Eq. (23), G = 0 and H = 1 in Eq. (24), E = 0 and F = 1 in Eq. (27c), and G = 0 and H = 1 in Eq. (27d) yields: kπ φ (z) = ψ (z) = χ (z) = sin z , (61) z z φ (t) = ψ (t) = χ (t) = sin(ωt). (62) t t Since η = 0 [c.f., (41)], it follows from Eqs. (16b) and (27b) that: (1) (2) φ (θ) = C + D θ, φ (θ) = C + D θ, χ (θ) = C + D θ. (63) 1 1 2 2 3 3 θ θ The constants Λ and Λ , as determined from Eq. (14b), are: 1 2 2 2 2 2 kπ ρω kπ ρω Λ = − > 0, Λ = − = 0. (64) 1 2 L (λ + 2μ) L μ Since η = 0, Λ > 0, and Λ = 0, it then follows from Eqs. (16c) and (27a) that: 1 2 p p (1) (2) φ (r) = A I Λ r + B K Λ r , φ (r) = A + B ln r, χ (r) = A + B ln r. 1 0 1 1 0 1 2 2 3 3 r r r (65) Furthermore, the constants γ (s = 1, 2) given by Eq. (22b) reduce (in this special case) to: 1 if s = 1 γ = . (66) 0 if s = 2 Inputting ingredients (61)-(66) into (29) and setting A = B = 0 in the result (since there 1 1 are no Bessel functions in the radial parts of (60)) yields the following general displacement field: B D kπ 2 3 u (r, θ, z, t) = (C + D θ) + (A + B ln r) sin z sin(ωt), (67a) r 2 2 3 3 r r L D B kπ 2 3 u (r, θ, z, t) = (A + B ln r) − (C + D θ) sin z sin(ωt), (67b) θ 2 2 3 3 r r L u (r, θ, z, t) = 0. (67c) 22 ln r Since there is no term in (60a), we choose B = 0 in (67a), which consequently eliminates the linear circumferential term in (67b). Furthermore, since U 6= U in (60a), the constant 1 2 B in (67a) must be non-zero. By defining a new set of arbitrary constants A ≡ A D , 2 2 2 2 ¯ ¯ C ≡ B C + A D , and D ≡ B D , (67) may then be compactly expressed as follows: 2 2 2 3 3 2 2 2 1 kπ ¯ ¯ u (r, θ, z, t) = C + D θ sin z sin(ωt), (68a) r 2 2 r L 1 kπ ¯ ¯ u (r, θ, z, t) = A + D ln r sin z sin(ωt), (68b) θ 2 2 r L u (r, θ, z, t) = 0. (68c) ¯ ¯ ¯ Note that each of the constants {A , C , D } has units of (length) . Application of boundary 2 2 2 condition (60a) using (68a) then yields: U θ − U θ U − U 1 2 2 1 2 1 ¯ ¯ C = R, D = R. (69) 2 2 θ − θ θ − θ 2 1 2 1 Condition (60b) is automatically satisfied by (68c). Substituting Eqs. (68) into the stress- displacement relations [Eqs. (D1) in Appendix D] yields the pertinent components of the stress field as follows: 2μ kπ ¯ ¯ σ (r, θ, z, t) = − C + D θ sin z sin(ωt), (70a) rr 2 2 r L 2μ kπ ¯ ¯ ¯ σ (r, θ, z, t) = − A − D + D ln r sin z sin(ωt), (70b) rθ 2 2 2 r L kπ μ kπ ¯ ¯ σ (r, θ, z, t) = C + D θ cos z sin(ωt), (70c) rz 2 2 L r L 2μ kπ ¯ ¯ σ (r, θ, z, t) = C + D θ sin z sin(ωt). (70d) θθ 2 2 r L Direct application of boundary conditions (59a)-(59c) and (60c) then yields a conditional ¯ ¯ solution to Problem A as follows. The displacement field (68) with constants C and D 2 2 ¯ ¯ given by (69) and A = D is a solution to Problem A provided the constant stresses 2 2 {A , B , C , F , G , S : i = 1, 2} are chosen as follows: i i i i i i 2μ 2μ ln R ¯ ¯ ¯ A = − C , B = − D , C = −2μ D , i 2 i 2 i 2 2 2 2 R R R i i i (71) μ kπ μ kπ 2μ ¯ ¯ ¯ ¯ F = C , G = D , S = C + D θ , i 2 i 2 i 2 2 i R L R L R i i ¯ ¯ where C and D are again given by (69). 2 2 23 B. Exponential Circumferential Variations Suppose instead that the radial surfaces of the shell are subjected to the following non- uniform stresses: kπ σ (R , θ, z, t) = A exp(−βθ) sin z sin(ωt), (72a) rr i i kπ σ (R , θ, z, t) = B exp(−βθ) sin z sin(ωt), (72b) rθ i i kπ σ (R , θ, z, t) = C exp(−βθ) cos z sin(ωt), (72c) rz i i where the index i = {1, 2}, {A , B , C : i = 1, 2} are prescribed constant stresses, k ∈ Z i i i and β > 0 are prescribed constants, and ω is the prescribed angular excitation frequency. Conditions (72a)-(72c) are pure stress boundary conditions and must be satisfied for all θ ∈ [θ , θ ], z ∈ (0, L), and arbitrary t. Suppose furthermore that the following (non- 1 2 uniform) conditions are prescribed on the circumferential faces of the shell: sin(β ln r) kπ u (r, θ , z, t) = D R sin z sin(ωt), (73a) r i i r L cos(β ln r) kπ u (r, θ , z, t) = D R sin z sin(ωt), (73b) θ i i r L u (r, θ , z, t) = 0, (73c) z i where again i = {1, 2}, {D : i = 1, 2} are prescribed constant displacements such that −β(θ −θ ) 2 1 D = D e and R ≡ (R + R )/2 is the mean radius of the shell. Conditions (73a)- 2 1 1 2 (73c) are pure displacement boundary conditions on the circumferential faces of the shell and must be satisfied for all r ∈ [R , R ], z ∈ (0, L), and arbitrary t. 1 2 μ kπ Problem B: In the special case where the excitation frequency ω = , determine ρ L the displacement field at all points of the clamped shell subject to conditions (72) on its radial surfaces and to conditions (73) on its circumferential faces. To obtain general solutions for {u , u , u } compatible with the given boundary condi- r θ z tions, fix the parameters κ and τ as in (41), fix η = −β , and employ prescription (26) for Buchwald potential χ. The axial and temporal parts of the Buchwald potentials are the same as those obtained for Problem A [c.f., Eqs. (61) and (62)]. Since η < 0, it follows from 24 Eqs. (16b) and (27b) that: (1) (2) −βθ +βθ −βθ +βθ −βθ +βθ φ (θ) = C e + D e , φ (θ) = C e + D e , χ (θ) = C e + D e . (74) 1 1 2 2 3 3 θ θ The constants Λ and Λ are the same as for Problem A [c.f., Eqs. (64)]. Since η < 0, 1 2 Λ > 0, and Λ = 0, it then follows from Eqs. (16c) and (27a) that: 1 2 p p (1) φ (r) = A I Λ r + B K Λ r , (75a) 1 iβ 1 1 iβ 1 (2) φ (r) = A cos(β ln r) + B sin(β ln r), (75b) 2 2 χ (r) = A cos(β ln r) + B sin(β ln r). (75c) 3 3 The constants γ (s = 1, 2) are the same as for Problem A [c.f., Eq. (66)]. Inputting the above ingredients into (29) and setting A = B = 0 in the result (since there are no Bessel 1 1 functions in the radial parts of (73)) and furthermore setting D = D = D = 0 (since 1 2 3 +βθ there are no e terms in (72)) yields the following general displacement field: u (r, θ, z, t) = {C [−A sin(β ln r) + B cos(β ln r)] − C [A cos(β ln r) + B sin(β ln r)]} r 2 2 2 3 3 3 kπ × exp(−βθ) sin z sin(ωt), (76a) u (r, θ, z, t) = − {C [A cos(β ln r) + B sin(β ln r)] + C [−A sin(β ln r) + B cos(β ln r)]} θ 2 2 2 3 3 3 kπ × exp(−βθ) sin z sin(ωt), (76b) u (r, θ, z, t) = 0. (76c) Since there is no cos(β ln r) term in (73a) and no sin(β ln r) term in (73b), we choose A = B = 0 in (76a) and (76b). For convenience, and without loss of generality, we may also set B = 0. By defining a new arbitrary constant C ≡ −A βC , (76) reduces to: 3 2 2 sin(β ln r) kπ u (r, θ, z, t) = C exp(−βθ) sin z sin(ωt), (77a) r L cos(β ln r) kπ u (r, θ, z, t) = C exp(−βθ) sin z sin(ωt), (77b) r L u (r, θ, z, t) = 0. (77c) Note that the constant C has units of (length) . Application of boundary condition (73a) using (77a) then yields: C = D exp(βθ )R = D exp(βθ )R. (78) 1 1 2 2 25 Condition (73b) is then automatically satisfied by (77b), and of course condition (73c) is trivially satisfied by (77c). Substituting Eqs. (77) into the stress-displacement relations [Eqs. (D1) in Appendix D] then yields the pertinent components of the stress field as follows: 2μC kπ σ (r, θ, z, t) = [β cos(β ln r) − sin(β ln r)] exp(−βθ) sin z sin(ωt), (79a) rr r L 2μC kπ σ (r, θ, z, t) = − [β sin(β ln r) + cos(β ln r)] exp(−βθ) sin z sin(ωt), (79b) rθ r L kπ sin(β ln r) kπ σ (r, θ, z, t) = μC exp(−βθ) cos z sin(ωt). (79c) rz L r L Direct application of boundary conditions (72a)-(72c) then yields a conditional solution to Problem B as follows. The displacement field (77) with constant C given by (78) is a solution to Problem B provided the constant stresses {A , B , C : i = 1, 2} are prescribed as follows: i i i 2μ A = [β cos(β ln R ) − sin(β ln R )] C, (80a) i i i 2μ B = − [β sin(β ln R ) + cos(β ln R )] C, (80b) i i i kπ sin(β ln R ) C = μ C, (80c) L R where C is again given by (78). VII. APPLICATION 3: FORCED VIBRATION OF AN OPEN SOLID CYLINDER Consider an (elastic) open circular solid cylinder whose domain is given by: 0 ≤ r ≤ R, 0 ≤ θ ≤ , and 0 ≤ z ≤ L. The flat ends of the cylinder, which are situated at z = 0 and z = L, are supported such that: σ (r, θ, 0, t) = σ (r, θ, L, t) = 0, (81a) rz rz σ (r, θ, 0, t) = σ (r, θ, L, t) = 0, (81b) θz θz u (r, θ, 0, t) = u (r, θ, L, t) = 0. (81c) z z Conditions (81) are admissible mixed boundary conditions specifying orthogonal components of the stress and displacement at the flat ends of the cylinder and must be satisfied for all 26 √ r ∈ [0, R], θ ∈ [0, π/ 101], and arbitrary t. The radial curved surface of the cylinder is subjected to the following non-uniform stresses: σ (R, θ, z, t) = A sin 101θ sin(ωt), (82a) rr σ (R, θ, z, t) = B cos 101θ sin(ωt), (82b) rθ σ (R, θ, z, t) = 0, (82c) rz where {A, B} are prescribed constant stresses and ω is the prescribed angular excitation frequency. Conditions (82) are pure stress boundary conditions and must be satisfied for all θ ∈ [0, π/ 101], z ∈ (0, L), and arbitrary t. The following uniform mixed conditions are furthermore prescribed on the circumferential faces of the cylinder: u (r, 0, z, t) = u (r, √ , z, t) = 0, (83a) r r σ (r, 0, z, t) = σ (r, √ , z, t) = 0, (83b) θθ θθ u (r, 0, z, t) = u (r, , z, t) = 0. (83c) z z Conditions (83) are admissible mixed boundary conditions specifying orthogonal components of the stress and displacement on the circumferential faces of the cylinder and must be satisfied for all r ∈ [0, R], z ∈ (0, L), and arbitrary t. Problem C: For arbitrary ω > 0, determine the displacement field at all points of the cylinder subject to: (i) conditions (81) at its flat ends; (ii) conditions (82) on its radial curved surface; and (iii) conditions (83) on its circumferential faces. The fact that conditions (82) are independent of z hints that we should employ the κ = 0 solution (38) with the axial parts of the Buchwald potentials equal to a constant. For n o e e convenience, we can choose {E = 1, F = 0} in Eq. (36) and E = 1, F = 0 in Eq. (37c), which yields: φ (z) = ψ (z) = χ (z) = 1. (84) z z The parameter τ can be set as in the previous two vibration problems (i.e., τ = −ω < 0) to generate the time dependence found in conditions (82). Choosing {G = 0, H = 1} in Eq. (24) n o e e and G = 0, H = 1 in Eq. (27d) then reproduces (62). Solution (38) subsequently reduces 27 to: dR (r) 1 dχ (θ) u = Θ (θ) + χ (r) sin(ωt), (85a) r 1 dr r dθ 1 dΘ (θ) dχ (r) u = R (r) − χ (θ) sin(ωt), (85b) θ 1 r dθ dr u = 0. (85c) By substituting (85) into (D1a) and comparing the result with (82a), it can be deduced that n o dχ (θ) d Θ (θ) the functions Θ (θ), , must each be of the form: sin 101θ . Setting the pa- 1 2 dθ dθ rameter η = 101 > 0, and then choosing {C = 0, D = 1} in Eq. (32b) and {C = 1, D = 0} 1 1 3 3 in Eq. (37a) gives: √ √ Θ (θ) = sin 101θ , χ (θ) = cos 101θ . (86) Then, choosing B = 0 in Eq. (32c) and B = 0 in Eq. (37b) (in order to satisfy the implicit 1 3 finiteness condition at r = 0) subsequently yields (since τ < 0 and η > 0): √ √ R (r) = A J (α r), χ (r) = A J (α r), (87) 1 1 1 3 2 101 r 101 q q 2 2 ρω ρω where α = and α = . Inputting (86) and (87) into (85) then yields the 1 2 λ+2μ μ displacement field: " # u = A J (α r) − A J (α r) sin 101θ sin(ωt), (88a) r 1 1 3 2 " # u = A J (α r) − A J (α r) cos 101θ sin(ωt), (88b) θ 1 1 3 2 u = 0, (88c) where the primes denote differentiation with respect to r. Note that displacement field (88) automatically satisfies conditions (81c), (83a), and (83c). Substituting (88) into the stress- displacement relations [Eqs. (D1) in Appendix D] then yields the pertinent components of 28 the stress field as follows: λ 101λ ′′ ′ √ √ σ (r, θ, z, t) = A (λ + 2μ)J (α r) + J (α r) − J (α r) rr 1 1 1 1 101 101 2 r r 2μ 101 1 + A J (α r) − J (α r) sin 101θ sin(ωt), (89a) 3 2 2 r r 1 101 ′′ ′ √ √ σ (r, θ, z, t) = μ A −J (α r) + J (α r) − J (α r) rθ 3 2 2 2 101 101 r r 2 101 1 + A J (α r) − J (α r) cos 101θ sin(ωt), (89b) 1 1 1 r r (λ + 2μ) 101 ′′ ′ √ √ σ (r, θ, z, t) = A λJ (α r) + J (α r) − J (α r) θθ 1 1 1 1 101 101 r r 2μ 101 1 + A J (α r) − J (α r) sin 101θ sin(ωt), (89c) 3 2 2 r r σ (r, θ, z, t) = 0, (89d) rz σ (r, θ, z, t) = 0, (89e) θz where the primes again denote differentiation with respect to r. By virtue of (89c) and (89e), conditions (83b) and (81b) are respectively satisfied. And by virtue of (89d), conditions (81a) and (82c) are naturally satisfied. Thus, the only condi- tions that remain to be satisfied are conditions (82a) and (82b). Application of these two remaining boundary conditions yields the following 2 × 2 linear system: a a A A 11 12 1 = , (90a) a a A B 21 22 3 where λ 101λ ′′ ′ √ √ a = (λ + 2μ)J (α R) + J (α R) − J (α R), (90b) 11 1 1 1 101 101 2 R R 2μ 101 1 a = J (α R) − J (α R) , (90c) 12 2 2 R R 2μ 101 1 a = J (α R) − J (α R) , (90d) 21 1 1 R R 1 101 ′′ ′ √ √ a = μ −J (α R) + J (α R) − J (α R) , (90e) 2 2 2 101 101 2 R R Although the matrix elements (90b)-(90e) can be written in several different forms free of derivatives, it is not relevant for our purposes here to do so. 29 and √ √ dJ (α r) d J (α r) 1/2 1/2 101 101 ′ ′′ √ √ J (α R) ≡ , J (α R) ≡ . (90f) 1/2 1/2 101 101 2 dr dr r=R r=R The solution to system (90) can then be readily obtained as follows: a A − a B a B − a A 22 12 11 21 A = , A = , (91) 1 3 a a − a a a a − a a 11 22 12 21 11 22 12 21 where {a , a , a , a } are as given by Eqs. (90b)-(90e). The solution to Problem C is 11 12 21 22 therefore given by the displacement field (88) with the constants {A , A } given by (91). 1 3 VIII. CONCLUSION In this paper, we have extended the work of SC1 to encompass a much broader class of separable solutions to the coupled equations of motion governing the Buchwald displacement potentials in cylindrical coordinates, namely, separable solutions whose circumferential parts are elementary 2π-aperiodic functions. Using the Buchwald representation (2), we then constructed eighteen distinct families of separable solutions to the NL equation; in each case, the circumferential part of the solution is one of three elementary 2π-aperiodic functions. To our knowledge, these solutions constitute the first comprehensive set of 2π-aperiodic cylindrical solutions to the NL equation. The solutions obtained in SC1 can reproduced by taking appropriate values of the angular parameter η and are thus special cases of the solutions obtained here. The obtained solutions can be directly applied to a fundamental set of linear-elastic boundary-value problems in cylindrical coordinates, in particular, problems wherein the displacement, stress, and strain fields are: constant, linear, sinusoidal, exponential, or hy- perbolic in the circumferential and longitudinal coordinates, and sinusoidal or exponential in the time coordinate. Such solutions are useful for solving a wide variety of dynamical prob- lems where 2π-periodicity in the angular coordinate is incompatible with the given boundary conditions. As illustrative examples, we showed how the obtained solutions can be used to solve various forced-vibration problems involving open thick cylindrical shells and open solid cylinders. As an addendum to SC1, we also included an example that illustrates the use of the special 2π-periodic solution (30), an important subsolution that is useful for problems where there is no explicit dependence on the circumferential coordinate. 30 The applications considered in this paper serve as a good starting point insofar as they are mathematically non-trivial but not overly complicated. More interesting and mathematically complex problems can be solved with the aid of the obtained solutions, but it is not within the scope of the present paper to do so. We hope to report on more advanced applications elsewhere. An important facet that we have not considered here is the rather large number of special solutions that derive from assuming one or more of the Buchwald potentials are zero, which as it turns out can yield more varieties of solutions and thus allow for more interesting applications. We hope to report on this topic in a future paper. Acknowledgments The authors acknowledge financial support from the Natural Sciences and Engineering Research Council (NSERC) of Canada and the Ontario Research Foundation (ORF). Appendix A: Fundamental Separable Solutions to the 2D Helmholtz Equation in Polar Coordinates In polar coordinates, the 2D Helmholtz equation is given by: 2 2 ∂ f 1 ∂f 1 ∂ f ∇ + Λ f = 0 ⇔ + + + Λf = 0, (A1) 2 2 2 ∂r r ∂r r ∂θ where f(r, θ) is some real-valued function of interest and Λ ∈ R is the Helmholtz constant. Fundamental separable solutions to (A1) can be obtained by assuming a solution of the form f(r, θ) = R(r)Θ(θ) and using separation of variables. Under the assumption that the angular function Θ(θ) is 2π-periodic, the solutions are well-known and readily available. The complementary solutions lacking 2π-periodicity do not however appear to be well-known nor are they (to our knowledge) readily found in the literature on PDEs. Since these complemen- tary (i.e., non-2π-periodic) solutions are fundamental building blocks in our mathematical construction, it is useful to establish them here for referential purposes. Substituting f(r, θ) = R(r)Θ(θ) into Eq. (A1) and then multiplying the result by r /R(r)Θ(θ) yields (after rearrangement of the radial and angular terms): ′′ ′ ′′ R (r) R (r) Θ (θ) 2 2 r + r + Λr = − . (A2) R(r) R(r) Θ(θ) 31 Since the left-hand side of Eq. (A2) depends only on r and the right-hand side depends only on θ, Eq. (A2) can be satisfied for all values of r and θ only if the terms on each side of the equation are equal to a constant. If we denote this constant by η, this leads to the equations: 2 ′′ ′ 2 r R (r) + rR (r) + Λr − η R(r) = 0, (A3) ′′ Θ (θ) + ηΘ(θ) = 0. (A4) The solution to Eq. (A4) is: p p C exp − |η|θ + D exp |η|θ if η < 0 Θ(θ) = , (A5) C + Dθ if η = 0 √ √ C cos ηθ + D sin ηθ if η > 0 where {C, D} are arbitrary constants. In the following, it will be assumed that η ∈ 2 + R\ {N : N ∈ Z }. The form of the solution to Eq. (A3) depends on the relative signs of Λ and η as follows. 1. Λ > 0 and η ≥ 0 In this case, Eq. (A3) is a standard (unmodified) Bessel equation with solution: √ √ √ √ R(r) = AJ Λr + BY Λr , (A6) η η where {A, B} are arbitrary constants. 2. Λ > 0 and η < 0 In this case, Eq. (A3) can be expressed in the form: 2 ′′ ′ 2 r R (r) + rR (r) + Λr + |η| R(r) = 0, (A7) which is again an unmodified Bessel equation whose solution is given by: √ √ √ √ R(r) = AJ Λr + BY Λr . (A8) i |η| i |η| Note that the Bessel functions in Eq. (A8) have purely imaginary orders, and are therefore complex-valued (see Ref. [16] for more information on such Bessel functions). 32 3. Λ < 0 and η ≥ 0 In this case, Eq. (A3) is of the form: 2 ′′ ′ 2 r R (r) + rR (r) − |Λ|r + η R(r) = 0, (A9) which is a modified Bessel equation with solution: p p √ √ R(r) = AI |Λ|r + BK |Λ|r . (A10) η η 4. Λ < 0 and η < 0 In this case, Eq. (A3) can be expressed in the form: 2 ′′ ′ 2 r R (r) + rR (r) − |Λ|r − |η| R(r) = 0, (A11) which is again a modified Bessel equation with solution given by: p p √ √ R(r) = AI |Λ|r + BK |Λ|r . (A12) i |η| i |η| Note that the modified Bessel functions in Eq. (A12) have purely imaginary orders; in this √ √ case, I (· · · ) is complex-valued, whereas K (· · · ) is real-valued. i |η| i |η| 5. Λ = 0 In this case, Eq. (A3) reduces to: 2 ′′ ′ r R (r) + rR (r) − ηR(r) = 0, (A13) which is a Cauchy-Euler (equidimensional) equation whose solution (depending on the value of η) is given by: √ √ η − η Ar + Br if η > 0 R(r) = . (A14) A + B ln r if η = 0 p p A cos |η| ln r + B sin |η| ln r if η < 0 33 Appendix B: A Brief on Bessel Functions with Purely Imaginary Orders For ν ∈ R and x ∈ (0, ∞), the functions {I (x) , K (x)} are linearly independent iν iν solutions of d y dy 2 2 2 x + x + ν − x y = 0. (B1) dx dx The functions I (x) and K (x) are complex- and real-valued, respectively. It can be shown iν iν that the Wronskian of I (x) ≡ Re{I (x)} and K (x) is non-zero [17], and thus these two ν iν iν functions together constitute a real and linearly independent set of solutions to Eq. (B1). For ν ∈ R and x ∈ (0, ∞), the functions {J (x) , Y (x)} are complex-valued linearly iν iν independent solutions of d y dy 2 2 2 x + x + x + ν y = 0. (B2) dx dx It can be shown that the Wronskian of the real functions πν J (x) ≡ sech Re{J (x)} (B3a) ν iν πν Y (x) ≡ sech Re{Y (x)} (B3b) ν iν is non-zero [17], and thus these functions together constitute a real and linearly independent set of solutions to Eq. (B2). Appendix C: Solution to Subsystem (3c) using Separation of Variables Substituting (25) into Eq. (3c) and then dividing the result by χ (r)χ (θ)χ (z)χ (t) yields r θ z t ′′ ′ ′′ ′′ ′′ χ (r) 1 χ (r) 1 χ (θ) χ (z) 1 χ (t) r r θ z t + + + = , (C1) χ (r) r χ (r) r χ (θ) χ (z) c χ (t) r r θ z T t where c ≡ μ/ρ. The first bracketed term depends only on r and θ, the second bracketed term only on z, and the right-hand side depends only on t. The above equation can be satisfied for all values of r, θ, z, and t only if each of the bracketed terms is a constant. This leads to the equations ′′ 2 χ (t) − υ c χ (t) = 0, (C2) t T t ′′ χ (z) − υ χ (z) = 0, (C3) z z and ′′ ′ ′′ χ (r) 1 χ (r) 1 χ (θ) r r θ + + = υ , (C4) χ (r) r χ (r) r χ (θ) r r θ 34 where the separation constants υ , υ , and υ obey (by virtue of Eq. (C1)) the relation t z r υ + υ = υ . (C5) r z t Equation (C4) can be written in the form ′′ ′ ′′ χ (r) χ (r) χ (θ) 2 r r 2 θ r + r − r υ = − . (C6) χ (r) χ (r) χ (θ) r r θ The left-hand side of Eq. (C6) depends only on r and the right-hand side only on θ. The two sides must therefore be equal to another constant, which we shall denote by υ . This second separation of variables yields the equations 2 ′′ ′ 2 r χ (r) + rχ (r) − r υ + υ χ (r) = 0, (C7) r θ r r r ′′ χ (θ) + υ χ (θ) = 0. (C8) θ θ Note that all real values of the separation constants υ , υ , υ , and υ are admissible but θ t z r that the latter three constants are not all independent. Upon independent specification of any two of {υ , υ , υ }, the third is automatically determined by relation (C5). Without loss t z r of generality, we may assign υ and υ to be independent parameters (that can take on any t z real values) whereupon υ = υ − υ . Solutions to Eqs. (C8), (C3), and (C2) are then given r t z by: p p C exp − |υ |θ + D exp |υ |θ if υ < 0 3 θ 3 θ θ χ (θ) = , (C9) C + D θ if υ = 0 3 3 θ √ √ C cos υ θ + D sin υ θ if υ > 0 3 θ 3 θ θ p p e e E cos |υ |z + F sin |υ |z if υ < 0 z z z e e χ (z) = , (C10) E + Fz if υ = 0 √ √ e e E exp − υ z + F exp υ z if υ > 0 z z z and p p e e G cos |υ |c t + H sin |υ |c t if υ < 0 t T t T t e e χ (t) = , (C11) G + Ht if υ = 0 √ √ e e G exp − υ c t + H exp υ c t if υ > 0 t T t T t n o e e e e respectively, where C , D , E, F, G, H are arbitrary constants. The function χ (θ) given 3 3 2 + by Eq. (C9) is 2π-periodic when υ ∈ {N : N ∈ Z } and in the special case υ = D = 0. θ θ 3 35 Using the results of Appendix A, solutions to Eq. (C7) can be immediately deduced; the solutions are as follows: p p √ √ A J |υ |r + B Y |υ |r if υ < 0 and υ < 0 3 r 3 r r θ i |υ | i |υ | θ θ p p A cos |υ | ln r + B sin |υ | ln r if υ = 0 and υ < 0 3 θ 3 θ r θ √ √ √ √ A I υ r + B K υ r if υ > 0 and υ < 0 3 r 3 r r θ i |υ | i |υ | θ θ p p A J |υ |r + B Y |υ |r if υ < 0 and υ = 0 3 0 r 3 0 r r θ χ (r) = , (C12) A + B ln r if υ = 0 and υ = 0 3 3 r θ √ √ A I υ r + B K υ r if υ > 0 and υ = 0 3 0 r 3 0 r r θ p p √ √ A J |υ |r + B Y |υ |r if υ < 0 and υ > 0 3 υ r 3 υ r r θ θ θ √ √ υ − υ θ θ A r + B r if υ = 0 and υ > 0 3 3 r θ √ √ √ √ A I υ r + B K υ r if υ > 0 and υ > 0 3 υ r 3 υ r r θ θ θ where {A , B } are arbitrary constants. 3 3 Appendix D: Stress-Displacement Relations For reference, we write here the linear elastic stress-displacement relations (in cylindrical coordinates), which provide the general mathematical connection between the components of the displacement and stress fields [18]: ∂u λ ∂u ∂u r θ z σ = (λ + 2μ) + + u + λ , (D1a) rr r ∂r r ∂θ ∂z ∂u (λ + 2μ) ∂u ∂u r θ z σ = λ + + u + λ , (D1b) θθ r ∂r r ∂θ ∂z ∂u λ ∂u ∂u r θ z σ = λ + + u + (λ + 2μ) , (D1c) zz r ∂r r ∂θ ∂z ∂u ∂u u r θ θ σ = μ + − = σ , (D1d) rθ θr r ∂θ ∂r r ∂u ∂u r z σ = μ + = σ , (D1e) rz zr ∂z ∂r ∂u 1 ∂u θ z σ = μ + = σ . (D1f) θz zθ ∂z r ∂θ [1] A. C. Eringen and E. S. Suhubi, Elastodynamics, Vol. 2, Linear Theory, Academic Press, New York, 1975. 36 [2] M. E. Gurtin, The Linear Theory of Elasticity, Mechanics of Solids, Vol. 2 (C. A. Truesdell, ed.), Springer, Berlin, 1972, pp. 1-295. [3] J. Miklowitz, The Theory of Elastic Waves and Wave Guides, Elsevier, Amsterdam, 1978. [4] V. B. Poruchikov, Methods of the Classical Theory of Elastodynamics, Springer, Berlin, 1993. [5] A. Rahman and F. Ahmad, Representation of the displacement in terms of scalar functions for use in transversely isotropic materials, J. Acoust. Soc. Am., 104 (1998), pp. 3675-3676. [6] F. Ahmad and A. Rahman, Acoustic scattering by transversely isotropic cylinders, Int. J. Eng. Sci., 38 (2000), pp. 325-335. [7] F. Ahmad, Guided waves in a transversely isotropic cylinder immersed in a fluid, J. Acoust. Soc. Am., 109 (2001), pp. 886-890. [8] F. A. Amirkulova, Acoustic and Elastic Multiple Scattering and Radiation from Cylindrical Structures, Ph.D. Thesis, Rutgers University, 2014. [9] J. Laurent, D. Royer, T. Hussain, F. Ahmad, and C. Prada, Laser induced zero-group velocity resonances in transversely isotropic cylinder, J. Acoust. Soc. Am., 137 (2015), pp. 3325-3334. [10] F. A. Amirkulova, Dispersion Relations for Elastic Waves in Plates and Rods, M.Sc. Thesis, Rutgers University, 2011. [11] P. Ponnusamy and R. Selvamani, Wave propagation in a homogeneous isotropic cylindrical panel, Adv. Theor. Appl. Math., 6 (2011), pp. 171-181. [12] R. Selvamani and P. Ponnusamy, Mathematical modelling of waves in a homogeneous isotropic rotating cylindrical panel, Rom. J. Math. Comp. Sci., 1 (2011), pp. 12-28. [13] J. Sakhr and B. A. Chronik, Solving the Navier-Lame equation in cylindrical coordinates using the Buchwald representation: some parametric solutions with applications, Adv. Appl. Math. Mech. 10 (2018), pp. 1025-1056. [14] D. Hui, Influence of geometric imperfections and in-plane constraints on nonlinear vibrations of simply supported cylindrical panels, J. Appl. Mech., 51 (1984), pp. 383-390. [15] C. W. Lim, K. M. Liew, and S. Kitipornchai, Vibration of open cylindrical shells: a three- dimensional elasticity approach, J. Acoust. Soc. Am., 104 (1998), pp. 1436-1443. [16] T. M. Dunster, Bessel Functions of Purely Imaginary Order with an Application to Second- Order Linear Differential Equations Having a Large Parameter, SIAM J. Math. Anal., 21 (1990), pp. 995-1018. [17] F. W. J. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark (Eds.), NIST Handbook of 37 Mathematical Functions, Cambridge University Press, Cambridge, 2010. [18] H. R. Hamidzadeh and R. N. Jazar, Vibrations of Thick Cylindrical Structures, Springer, New York, 2010.
Journal
Mathematical Physics – arXiv (Cornell University)
Published: Jun 25, 2019