Added-Mass Based Efficient Fluid–Structure Interaction Model for Dynamics of Axially Moving Panels with Thermal Expansion

This is a self-archived version of an original article. This version may differ from the original in pagination and typographic details. Author(s): Title: Year: Version: Copyright: Rights: Rights url: Please cite the original version: CC BY 4.0 https://creativecommons.org/licenses/by/4.0/ Added-Mass Based Efficient Fluid–Structure Interaction Model for Dynamics of Axially Moving Panels with Thermal Expansion © 2020 by the authors. Licensee MDPI, Basel, Switzerland. Published version Banichuk, Nikolay; Ivanova, Svetlana; Makeev, Evgeny; Jeronen, Juha; Tuovinen, Tero Banichuk, N., Ivanova, S., Makeev, E., Jeronen, J., & Tuovinen, T. (2020). Added-Mass Based Efficient Fluid–Structure Interaction Model for Dynamics of Axially Moving Panels with Thermal Expansion. Mathematical and Computational Applications, 25(1), Article 9. https://doi.org/10.3390/mca25010009 2020 Mathematical and Computational Applications Article Added-Mass Based Efficient Fluid–Structure Interaction Model for Dynamics of Axially Moving Panels with Thermal Expansion Nikolay Banichuk 1, Svetlana Ivanova 1, Evgeny Makeev 1, Juha Jeronen 2 and Tero Tuovinen 2,* 1Ishlinsky Institute for Problems in Mechanics of the Russian Academy of Sciences, Prospekt Vernadskogo 101-1, Moscow, 119526, Russia; [email protected] (N.B.); [email protected]u (S.I.); [email protected] (E.M.) 2Faculty of Information Technology, University of Jyväskylä, Mattilanniemi 2, 40014 Jyväskylä, Finland; [email protected] *Correspondence: ter[email protected]; Tel.: +358504413685 Received: 13 December 2019; Accepted: 20 January 2020; Published: 22 January 2020 Abstract: The paper considers the analysis of a traveling panel, submerged in axially flowing fluid. In order to accurately model the dynamics and stability of a lightweight moving material, the interaction between the material and the surrounding air must be taken into account. The lightweight material leads to the inertial contribution of the surrounding air to the acceleration of the panel becoming significant. This formulation is novel and the case complements our previous studies on the field. The approach described in this paper allows for an efficient semi-analytical solution, where the reaction pressure of the fluid flow is analytically represented by an added-mass model in terms of the panel displacement. Then, the panel displacement, accounting also for the fluid–structure interaction, is analyzed with the help of the weak form of the governing partial differential equation, using a Galerkin method. In the first part of this paper, we represent the traveling panel by a single partial differential equation in weak form, using an added-mass approximation of the exact fluid reaction. In the second part, we apply a Galerkin method for dynamic stability analysis of the panel, and present an analytical investigation of static stability loss (divergence, buckling) based on the added-mass model. Keywords: mathematical modeling; axially moving panel; fluid–elastic interaction; aerothermoelastic stability; divergence; added-mass approximation 1. Introduction In this study, we develop a mathematical model representing the behavior of an open draw in the paper making process, taking into account the effects of fluid–structure interaction and of thermal expansion of the material. Open draws are the spans, where the web travels without support and interacts with gas (fluid). Within these open draws, the paper web is subject to temperature forces. Stability of these parts is therefore an interesting research question for the applications. As for applied aspects of this study, axially moving materials have many applications in industry. Examples of systems with axially moving materials are band saws, power transmission belts, paper making processes, printing presses, manufacturing of plastic films and sheets. The model can also be applied to other open draws involving lightweight materials, such as the vibrations of a tape in a tape drive. In our paper, we however concentrate on paper making processes. Examples including transportation from the press nip to the drying section. We have simplified the problem, while capturing relevant aspects of the physical phenomena. A key point is that the behavior of the moving material is strongly dependent on problem parameters such Math. Comput. Appl. 2020,25, 9; doi:10.3390/mca25010009 www.mdpi.com/journal/mca Math. Comput. Appl. 2020,25,9 2 of 16 as panel and fluid velocities, densities of the fluid and the panel, bending rigidity, temperature, in-plane tension, and also external forces and geometric factors. For lightweight materials, the interaction between the material and the surrounding air must be taken into account because it significantly changes the dynamic response, as was first observed by [1]. In many previous investigations, axially moving materials have been modeled as traveling flexible strings, beams, panels and plates. Some classical articles have been made by Archibald and Emslie [ 2 ], Mote [ 3 ], Wickert and Mote [ 4 ], Wickert [ 5 ], Parker [ 6 ] and Wang et al. [ 7 ]. The majority of the performed studies have been devoted to isotropic moving materials. As for orthotropic moving materials, many important results, concerning modeling and the analysis of behavior, can be found in the monographs by Marynowski [ 8 ], Banichuk et al. [ 9 ], Tuovinen [ 10 ], and again Banichuk et al. [ 11 ]. In the article by Liu et al. [ 12 ], free vibrations of nonuniform rectangular plates varying in one or two directions are investigated. The plate was orthotropic and resting upon a Winkler-type spring foundation. In the article by Ghayesh and Amabili [ 13 ], the investigation of the transverse vibration response, stability, and bifurcations of an axially moving viscoelastic beam with time-dependent axial speed was considered. Research by Hozhabrossadati and Sani [ 14 ] deals with the determination of the static displacement function of a Euler–Bernoulli beam with two guided supports. An extensive literature review about mechanics of axially moving continua can also be found in Marynowski & Kapitaniak [ 15 ]. Note also that the papers by Banichuk et al. [ 16 ] and Banichuk and Ivanova are thematically adjointed [17]. Fundamental problems of fluid–structure interaction have been studied (Banichuk and Mironov [ 18 , 19 ], Ashley and Landahl [ 20 ], Anderson [ 21 ] and analyzed Bisplinghoff and Ashley [ 22 ], Ashley and McIntosh [23], Kornecki et al. [24], Vaughan and Raman [25], and Banichuk et al. [9]). In a book by Banichuk et al. [ 9 ], following an approach similar to Ashley and Landahl [ 20 ], in the case of potential flow, the aerodynamic reaction was found as a functional fluid–structure interaction problem involving the solution of integrodifferential equation with a singular kernel. The problem can be further simplified by approximating the singular kernel. For a review on studies concerning circular cylindrical shells including fluid–structure interaction, see Amabili and Paidoussis [ 26 ]. For shells containing flowing fluid, see Amabili et al. [ 27 ] and [ 28 ], Amabili and Garziera [ 29 ] and [ 30 ]. Periodically supported plates subjected to axial flow have been considered in Tubaldi et al. [31]. In Liu and Chang [ 32 ], a fluid–structure interaction problem dealing with free vibration analysis of an axisymmetrical, nonuniform circular plate with quadratic change of thickness in contact with fluid was investigated. Moreover, taking into account thermoelastic behaviour, the forced vibration analysis of nonhomogeneous thermoelastic, isotropic, thin annular disk under periodic and exponential types of axisymmetric dynamic pressures applied on its inner boundary has been performed and an analytical benchmark solution has been obtained, e.g., Mishra et al. [ 33 ]. In the article by Mao et al. [ 34 ], the authors have studied the thermoelastic instability (TEI) of a functionally graded material (FGM) half-plane sliding against a homogeneous half-plane at the in-plane direction. The interaction of the frictional heat and thermal contact resistance is taken into account in the TEI analysis. In this study, we formulate the dynamic aerothermomechanical problem for an axially moving panel, interacting with potential flow and subject to thermal expansion. This formulation is novel and the case complements our previous studies on the field. As for simplification of the model, we approximate the rigorous solution for the fluid reaction via an added-mass approach. As a result, the integrodifferential governing equation reduces to a partial differential equation, where we use added-mass approximations of the original inertial, Coriolis and centrifugal terms. In a steady state, the equation is further reduced into an ordinary differential equation. The resulting differential equations can be effectively solved for various initial and boundary conditions. As a result, we find solutions for some stationary and time-dependent problems. In the stationary case, we find the analytical solution of the stability problem and present expressions for the critical divergence velocity and for critical temperature of the loss of stability. We also study the dynamic stability of the system. Math. Comput. Appl. 2020,25,9 3 of 16 The proposed formulation of problem and method of analysis are simplified to give clear physical motivation and realize effective solutions. First, simplification consists of representation of moving elastic plate (two-dimensional model) as moving in the axial direction elastic panel, describing by a one-dimensional deflecting function w=w(x , t) . The travelling panel is interacting with a two-dimensional model of fluid moving in axial and normal (transverse) directions. These assumptions give us the possibilities to apply complex variable techniques for finding analytical solutions of a two-dimensional hydrodynamical (plane) problem with arbitrary distribution of elastic deflections w(x , t) . Thus, the original coupled problem of fluid–structure interaction admits the reduction to a one-dimensional integro-differential equation. This equation can be used as a test for some numerical solutions with 3D hydrodynamical and 2D elastodynamical parts of the original problem. The analytical approach, proposed in this paper, has some constraints. The model does not cover all possible aspects of the physical situation. The one-dimensional geometry of the panel limits the applicability to long and narrow draws. For short and wide draws, a plate model should be used, due to the localization of deflections near the free edges (Banichuk et al. [ 9 ]). In addition, paper as a material is slightly viscoelastic. It is known that the presence of viscosity, no matter how small, will qualitatively change the behavior of the stability exponents. However, the changes appear in the post-critical range, and the first critical velocity remains almost the same, so, for predicting the initial loss of stability, the elastic material model is sufficient. Moreover, the moving fluid is supposed to be ideal and the property of viscosity is excluded. If in reality liquid or gas are not ideal, but the viscosity is small, then the obtained ideal solution can be considered as a zero-order approximation in the small parameter method. Perturbation with respect to small viscosity parameter µ will be realized by the first order approximation. Finally, we have not accounted for the temperature dependence of the Young’s modulus. This is expected not to have a significant effect on the results. The bending rigidity parameter is small, so it mainly has a singular perturbation effect on the partial differential equation analyzed. For small but finite bending rigidity, the behavior of the stability exponents describing free vibrations is qualitatively different from the case with no bending rigidity (Wang et al. [ 7 ] and Jeronen [ 35 ]), but the first critical velocity remains almost the same. Thus, even if the bending rigidity slightly varies due to the temperature dependence of the Young’s modulus, it will not have a major effect on the critical velocity predicted by the model. 2. Basic Relations and Aerothermoelastic Model Consider a traveling thermoelastic panel subjected to a potential flow, where the free stream flows toward the right at velocity n• (with respect to the laboratory frame), see Figure 1. The governing equation describing small transverse vibrations of the panel is m∂2w ∂t2+2mV0 ∂2w ∂x∂t+mV2 0 ∂2w ∂x2T∂2w ∂x2+D∂4w ∂x4=qf+g, (1) where w is the transverse displacement, qf⌘qf(w) is the fluid reaction pressure, and g⌘g(x , t) represents external forces inside the domain. The mass per unit area is m . The panel travels axially at the velocity V0 and is subjected to in-plane axial mechanical tension Tm and thermal compression Tq . The unit of Tmand Tqis force per unit length. The bending rigidity of the panel is D. We have T=TmTq,Tm=T0, Tq=Eh 1n#q,D=Eh3 12(1n2). (2) Here, h is the thickness of the panel, E is the Young’s modulus of the panel material, n is its Poisson ratio, T0is a given constant tension, and #qis the generalized thermal strain defined by #q=1 hZh/2 h/2 aqq(z)dz. (3) Math. Comput. Appl. 2020,25,9 4 of 16 Figure 1. One-dimensional model of the moving panel. The coefficient of linear thermal expansion aq is constant, and the function q takes the form q(z)=qabs(z)q0 , where qabs and the reference temperature q0 are given in Kelvin. If the panel is at a uniform temperature of q0 , then q(z)= 0, and consequently #q= 0. In the case where aq and q are constant, we have #q=a0q. The reaction pressure qf is unknown, to be solved from the flow model in terms of w . The reaction pressure describes how the surrounding fluid pushes back on the panel, in reaction to the panel vibration. The external forces g are considered known, and are allowed to vary dynamically. For finding the fluid reaction pressure, we apply techniques from the aerodynamics of thin aerofoils, constructing a Green’s function type solution, as was made in [ 18 , 19 ] and presented in the books [ 9 , 20 – 22 , 36 ]. In non-dimensional coordinates x0=x/`,t0=t/twe have (primes omitted) qf(x,t)=rf `✓` t ∂ ∂t+v•∂ ∂x◆Z1 1N(x,x)✓` t ∂ ∂t+v•∂ ∂x◆w(x,x)dx, (4) N(x,x)⌘1 pln  1+L(x,x) 1L(x,x) ,L(x,x)=s(1x)(1+x) (1x)(1+x), (5) (x,t)2W,W={(x,t):1<x<1, 0 <t<•}, (6) where v•= 0 means no free-stream flow in laboratory coordinates and v•=V0 means that the whole air mass moves with the panel. Briefly stated, this Green’s function solution describes two-dimensional potential flow in the plane excluding a linear cut. To make an analytical solution possible, since we consider small vibrations around the trivial equilibrium position only, we geometrically approximate the panel as a straight line segment. The region z= 0,  1 <x< 1 is cut out from the plane, approximating the space occupied by the panel. The vibrations of the panel are taken into account in the no-penetration boundary condition on the velocity of the flow. This leads to Equations (4)–(6). A detailed treatment can be found in the book [9]. In (4), rf is the density ( kg/ m 3 ) of the surrounding medium. The spatial scaling factor ` is the half-length of the span, and t is an arbitrary scaling factor for time coordinate. For a physically meaningful time scaling, one can choose, e.g., t=`/C , where C=pT/m , the critical velocity of a traveling ideal string or membrane. The unit of tis [t]=s. Added-mass models are sometimes used for taking into account the inertial effects of fluid–structure interaction, expressed by an inertial operator in a simplified approximate setting. In this study, we consider an added-mass approximation for the above model of fluid–structure interaction of an axially moving panel. We will use the constructed approximation of the total inertial term md2w/dt2 for the solution of the problems of vibration and stability. As was shown in the book [9], we may use the following approximation: Z1 1N(x,x)f(x)dx⇡µZ1 1d(x,x)f(x)dx=µf(x),x2[1, 1], (7) Math. Comput. Appl. 2020,25,9 5 of 16 where d(x,x)is the Dirac delta distribution, and µis the constant µ⌘mean x2(1,1)Z1 1N(x,x)dx=1 2Z1 1Z1 1N(x,x)dxdx⇡p 4. (8) The symbol mean(·) denotes the average of the indicated quantity over the indicated interval. The idea behind the approximation (7) is that N(x , x) is a singular but integrable ( L1 ) kernel, which decays quickly as |xx| increases; it can be shown that N(x , x)<ln( 4 /|xx|) . Hence, for any sufficiently well-behaved f , most of the contribution to the integral on the left-hand side comes from the singularity. In physical terms, we approximate the fluid reaction pressure as a pointwise local effect, whereas in reality a disturbance in the flow propagates along the fluid, causing the pressure field to change globally (limited by the speed of sound in the medium). Especially from a mathematical modeling perspective, recall that a potential flow will instantly reconfigure itself to satisfy the Laplace equation. Any disturbance anywhere in the fluid domain will instantly cause global (even if small) effects. These effects are ignored in the present approach; this makes the solution process significantly simpler at the cost of some accuracy. Another approach is to treat the singular integrals numerically, which has been performed earlier in [9] and [35]. Let us now apply (7) and (8) to Equation (1). For now, consider only the inertial terms and qf . In other words, in Equation (1), let T= 0, D= 0 and g⌘ 0, as we may easily add these terms back when finished. We approximate qf by inserting Equation (7) and µ=p/ 4 from Equation (8) into Equation (4). We have m t2 ∂2w ∂t2+2mV0 `t ∂2w ∂x∂t+mV2 0 `2 ∂2w ∂x2=rf p 4✓` t2 ∂2w ∂t2+2v• t ∂2w ∂x∂t+v2 • ` ∂2w ∂x2◆. (9) Collecting terms that have the same derivative of w, we obtain 1 t2(m+ma)∂2w ∂t2+2V0 `t(m+marv)∂2w ∂x∂t+V2 0 `2⇣m+mar2 v⌘∂2w ∂x2=0, (10) where we have defined ma=p`rf 4,rv=v• V0 . (11) The expressions (10) and (11) reduce to a classical one- or three-term single-parameter added-mass model by choosing rv= 0 (i.e., v•= 0, no free-stream flow in laboratory coordinates) or rv= 1 (i.e., v•=V0, the whole air mass moves with the panel). 3. Non-Dimensional Initial Boundary Value Problem Using (10) and (11), the governing equation for dynamical behavior can be written in non-dimensional form as 1 t2(m+ma)∂2w ∂t2+2V0 `t(m+marv)∂2w ∂x∂t+V2 0 `2⇣m+mar2 v⌘∂2w ∂x2 T1 `2 ∂2w ∂x2+D1 `4 ∂4w ∂x4=1 `g(x,t). (12) Define C=pT/m, (13) which is the critical velocity of a traveling string or membrane with no bending rigidity. Let us define also the non-dimensional quantities g=` mrf,b=D mC2`2=D `2T, Math. Comput. Appl. 2020,25,9 6 of 16 c=V0 C,k=v• C,a=` tC,rm=ma m. (14) We will also use the non-dimensional external load g0x0,t0=` mC2g`x0,tt0=` Tg(`x0,tt0), (15) where g is the original dimensional load function. Then, we multiply Equation (12) by the factor `2/mC2, and use the definitions (13)–(15). We obtain a2(1+rm)∂2w ∂t2+2ac(1+rvrm)∂2w ∂x∂t+⇣c21+c2r2 vrm⌘∂2w ∂x2+b∂4w ∂x4=g(x,t), (16) where all the variables and functions are in non-dimensional form, and the domain is W={1<x<1, t>0}. We will use Equation (16) to determine the dynamical behavior and static loss of stability. As for the boundary conditions, we will use the simply supported conditions: (w)x=1=(w)x=1=0, (17) ✓∂2w ∂x2◆x=1 =✓∂2w ∂x2◆x=1 =0. (18) A direct simulation of dynamical behavior requires also two initial conditions (w)t=0=g1(x), (19) ✓∂w ∂t◆t=0 =g2(x), (20) where g1(x) and g2(x) are known. For determining free vibration modes, no initial conditions are needed. 4. Static Problem of Aerothermoelastic Stability In the context of static analysis, we set ∂2w ∂t2(x,t)=0, ∂2w ∂x∂t=0, g(x,t)=0 (21) in Equation (16). For the static stability problem, we use the Euler approach of determining a nontrivial steady-state solution, and the associated critical values of the problem parameters of interest. We concentrate on finding the critical panel velocity and critical temperature, which play a role analogous to Euler’s critical compression force for the axially compressed beam. Using (21) in (16), the steady-state equation is hc2⇣1+rmr2 v⌘1i∂2w ∂x2+b∂4w ∂x4=0, (22) with boundary conditions (19) and (20). Equation (22) describes the buckling of a traveling panel submerged in potential flow. Math. Comput. Appl. 2020,25,9 7 of 16 It is convenient to transform the original problem (22) with boundary conditions (17) and (18) into a sequence of two problems. Let us introduce a new variable y: y(x)=∂2w ∂x2,1x1. (23) This transforms Equation (22) into ∂2y ∂x2+ly =0, 1x1, (24) y(1)=0, y(1)=0, (25) where the introduced eigenvalue l(the load parameter of the stability problem) is l=1 bhc2⇣1+rmr2 v⌘1i. (26) The solution of (24) is y(x)=A1sin ✓pl✓x+1 2◆◆+A2✓cos ✓pl✓x+1 2◆◆◆, where A1 and A2 are arbitrary constants. With the help of boundary conditions (25), the solution is reduced to y(x)=A1sin ✓pl✓x+1 2◆◆, (27) l=✓jp 2◆2 ,j=1, 2, . . . (28) The minimal value of l corresponds to j= 1, and consequently lmin =p2/ 4 . Thus, we have the following basic relation for the critical parameters of static loss of stability: p2 4=`2 D⇣mV2 0+mav2 •T⌘=`2 D✓mV2 0+mav2 •T0+Eh 1n#0◆. (29) As follows from (29) and the expressions for #q and ma , the critical velocity of static loss of stability (divergence, buckling) and the critical temperature are expressed as ⇣V2 0⌘div =1 m✓p2D 4`2+T0Eh 1na0qp` 4rfv2 •◆, (30) qdiv =1n Ehaq✓p2D 4`2+T0mV2 0p` 4rfv2 •◆. (31) The influence of the added mass of the fluid on the critical velocity is represented by ✓V0 v•◆div =rKma m,K=1 v•✓p2D 4m`2+T m◆, (32) and is shown in Figure 2by the solid lines corresponding to K=1, 2, 3, 4. Math. Comput. Appl. 2020,25,9 8 of 16 Figure 2. Critical ratio V0/n•as a function of the quotient ma/m. 5. Dynamic Problem of Free Vibrations For numerical analysis, a discrete approximation will be used for the partial differential equation (16) with boundary conditions (17) and (18) and initial conditions (19) and (20). In this context, Galerkin methods are especially convenient and well-known techniques to solve the problem. We will space-discretize using the finite element method with C2continuous Hermite elements. To characterize the dynamics of the system, let us perform a numerical dynamic stability analysis. Following [ 37 ], we set the load as g⌘ 0, and formulate an eigenvalue problem, using the time-harmonic trial function w(x,t)=exp(st)W(x), (33) where s is the stability exponent (a complex number) and W(x) is the vibration mode. We will solve the problem for eigenvalue–eigenfunction pairs (s , W) . Loss of stability occurs at such values of the axial drive velocity V0 , where at least one eigenvalue s transitions to the positive half-plane (i.e., where the real part Re s becomes positive). The critical velocity is the smallest positive V0 such that stability is lost. In the original continuum problem, there will be a countably infinite number of solutions, corresponding to the spectrum of vibration modes. When the problem is discretized, if there are N degrees of freedom in the discretization, it will have 2 N solutions because, upon inserting (33) into (16), the resulting ordinary differential equation will be a quadratic polynomial in s . This results in a quadratic eigenvalue problem, which can be converted into a twice larger generalized linear eigenvalue problem by the companion form technique, see [ 38 ]. However, not all of the numerical solutions will be solutions of the continuum problem. Attempting to represent an infinite spectrum using a finite basis leads to an aliasing phenomenon. Best accuracy will be obtained for the first few of the lowest modes (slowest vibration, smallest imaginary part |Im s| ). Below, we will investigate the four lowest modes. Because the material is fully elastic, and hence the model ignores dissipation, then if s is an eigenvalue, s is also. This is in addition to the general symmetry (also valid for viscoelastic models) that, if s is an eigenvalue, then also conj(s) is. This latter property only reflects the physical symmetry with respect to the xaxis (z=0) in the problem setup. Math. Comput. Appl. 2020,25,9 15 of 16 17. Banichuk, N.; Ivanova, S.Y. Mathematical modelling of the axially moving panels subjected to thermomechanical actions. Mech. Based Des. Struct. Mach. 2018,46, 101–109. 18. Banichuk, N.; Mironov, A. Problems in optimization of vibrating plates submerged in an ideal fluid. PMM 1976,4, 520–527. 19. Banichuk, N.; Mironov, A. A scheme of a jet flow for investigating the equilibrium shapes of elastic plates submerged in a stream of a fluid and some optimization problems. PMM 1979,43, 83–90. 20. Ashley, H.; Landahl, M. Aerodynamics of Wings and Bodies; Dover: New York, NY, USA, 1985. 21. Anderson, Jr., J.D. Fundamentals of Aerodynamics; McGraw-Hill: New York, NY, USA, 1985. 22. Bisplinghoff, R.L.; Ashley, H. Principles of Aeroelasticity; Dover Publications, Inc.: New York, NY, USA, 1962. 23. Ashley, H.; McIntosh, S.C. Applications of aeroelastic constraints on structural optimization. In Proceedings of the 12th International Congress of Theoretical and Applied Mechanics, Stanford, CA, USA, 26–31 August 1968; Springer-Verlag: Berlin, Germany, 1969; pp. 100–113. 24. Kornecki, A.; Dowell, E.H.; O’Brien, J. On the aeroelastic instability of two-dimensional panels in uniform incompressible flow. J. Sound Vib. 1976,47, 163–178. 25. Vaughan, M.; Raman, A. Aeroelastic Stability of Axially Moving Webs Coupled to Incompressible Flows. ASME J. Appl. Mech. 2010,77, 021001. doi:10.1115/1.2910902. 26. Amabili, M.; Païdoussis, M.P. Review of studies on geometrically nonlinear vibrations and dynamics of circular cylindrical shells and panels, with and without fluid–structure interaction. Appl. Mech. Rev. 2003 , 56, 349–381. doi:10.1115/1.1565084. 27. Amabili, M.; Pellicano, F.; Païdoussis, M.P. Non-linear dynamics and stability of circular cylindrical shells containing flowing fluid. Part I: Stability. J. Sound Vib. 1999,225, 655–699. doi:10.1006/jsvi.1999.2255. 28. Amabili, M.; Karagiozis, K.; Païdoussis, M.P. Effect of geometric imperfections on non-linear stability of circular cylindrical shells conveying fluid. Int. J. Non-Linear Mech. 2009 ,44, 276–289. doi:10.1016/j.ijnonlinmec.2008.11.006. 29. Amabili, M.; Garziera, R. Vibrations of circular cylindrical shells with nonuniform constraints, elastic bed and added mass. Part II: Shells containing or immersed in axial flow. J. Fluids Struct. 2002 ,16, 31–51. doi:10.1006/jfls.2001.0402. 30. Amabili, M.; Garziera, R. Vibrations of circular cylindrical shells with nonuniform constraints, elastic bed and added mass. Part III: Steady viscous effects on shells conveying fluid. J. Fluids Struct. 2002 ,16, 795–809. doi:10.1006/jfls.2002.0446. 31. Tubaldi, E.; Alijani, F.; Amabili, M. Non-linear vibrations and stability of a periodically supported rectangular plate in axial flow. Int. J. Non-Linear Mech. 2014 ,66, 54–65. doi:10.1016/j.ijnonlinmec.2013.12.004. 32. Liu, M.F.; Chang, T.P. Axisymmetric Vibration of a Varying-Thickness Circular Plate in Contact with Fluid. Mech. Based Des. Struct. Mach. 2004,32, 39–56. 33. Mishra, K.C.; Sharma, J.N.; Sharma, P.K. Analysis of vibrations in a nonhomogeneous thermoelastic thin annular disk under dynamic pressure. Mech. Based Des. Struct. Mach. 2017,45, 207–218. 34. Mao, J.J.; Ke, L.L.; Yang, J.; Kitipornchai, S.; Wang, Y.S. Thermoelastic instability of functionally graded materials with interaction of frictional heat and contact resistance. Mech. Based Des. Struct. Mach. 2018 , 46, 139–156. 35. Jeronen, J. On the mechanical stability and out-of-plane dynamics of a travelling panel submerged in axially flowing ideal fluid: A study into paper production in mathematical terms. Ph.D. Thesis, Department of Mathematical Information Technology, University of Jyväskylä, 2011. 36. Banichuk, V.V. Problems and Methods of Optimal Structural Design; Plenum Press: New York, NY, USA, 1983. 313 Pages. 37. Bolotin, V.V. Nonconservative Problems of the Theory of Elastic Stability; Pergamon Press: New York, NY, USA, 1963. 38. Tisseur, F.; Meerbergen, K. The quadratic eigenvalue problem. SIAM Rev. 2001,43, 235–286. 39. Adams, R.A. Sobolev Spaces; Academic Press: Cambridge, MA, USA, 1975. Math. Comput. Appl. 2020,25,9 16 of 16 40. Kouko, J. Effects of Heating, Drying and Straining on the Relaxation and Tensile Properties of Wet Paper. Ph.D. Thesis, Department of Physics, University of Jyväskylä, 2014. 41. Bender, C.M.; Orszag, S.A. Advanced Mathematical Methods for Scientists and Engineers I: Asymptotic Methods and Perturbation Theory; Springer-Verlag: New York, NY, USA, 1978. 1999 Reprint: ISBN 978-0-387-98931-0. c 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).