format: text

text:
# Scattering of the acoustic waves on finite cylindrical covers

The interaction problems of wave fields with resilient bounded bodies arise in the different fields of mechanics, hydro acoustics, geophysics and seismology.

The analytical solutions of the problems of hydro acoustic waves scattering by resilient bodies are directly connected with the methods development of modeling wave fields interaction processes with underwater objects. Increasing the complexity of underwater objects geometry we get the increase of mathematical complexity of the pressure field determination problem solution. In this connection the problems are solved with simplified mathematical statement, which analysis helps to indicate the significant features of geometry and structure of modeling object. Up till now number of analytical and experimental researches were implemented according to the statement of the problem. These researches were directed to the conformities determination of influence of any construction or technologic property on structure parameters of reflected signals.

The problems of scattering of acoustic pressure waves by resilient bodies with complex geometry have no analytical solution.

Applying the direct numerical methods (e.g. finite difference method) to solve partial differential equations of these problems is not expedient even using powerful computers due to multi-dimensional type of the problem. For this reason high-frequency and low-frequency methods are often used to solve such problems. But when wavelength is comparable by order with the characteristic size of disperse body the frequency lies in moderate range. Exactly this range as shown in  is highly informational.

To solve the scattering problems in the moderate frequency range the limited integral equations method is applied.

Applying Grin theorem let us reduce the solution of stationary problem of resilient body dispersion to the solution of combined singular equations for the system of potential displacements. The solution of the above determines the surface potentials, which are used to define the outside field. Thus, the number of dimension variables is reduced and it becomes possible to use numerical methods.

In works  the problem of sound diffraction on the finite cylinder with mixed boundary conditions is solved by Grin’s function method. In  the axis-symmetrical dispersed pressure of resilient disk is computed using finite element method and HelmHoltz integral equation.

As is well known, there is a tight relation between the solution of scattering and emission problems. It could be defined by reciprocity theorem . This relation could be used to transfer emission problem results to scattering problems and vise versa .

While solving complex problems, for example three-dimensional problems of acoustic waves scattering on the resilient bodies or two-dimensional problems of scattering on the piece-wise-smooth resilient bodies and problems of scattering of acoustic waves on the thin-wall covers, troubles with definite satisfaction of the appropriate boundary conditions arise.

The considered work is devoted to the research of scattered and emissed acoustic waves by finite cylindrical covers, resilient cylinders and round plates.

The presence of ribs fracture line in bottoms interface and cylindrical covers on the face of resilient cylinder and on edge of round plates inputs the additional troubles of evaluating of such problems.

During last years for the solution of space problems of scattering of acoustic waves numerical and integral equations methods are effectively applied. Despite the achieved researches of applying these methods to three-dimensional problems, questions of efficiency increase of existing methods and their PC-realization are still actual. In particular, one of the hydro-acoustic problems is the evaluation of scattered and emitted pressure field by three-dimensional bodies of cylindrical form. But during the scattered pressure field exploration near bounded disperser with edges and also internal features of dissipater troubles arise, which are hardly solved without analytical research. That is why the important problem of analytic-numerical methods development and development of software for acoustic space interaction with finite bodies of cylindrical form are actual.

During the solution of problems of emission and scattering of acoustic waves by bounded bodies of cylindrical form we use the strategy based on applying of bounded integral equations method combined with series method and later using the improvement of convergence method, which consider particularities of desired functions. The proposed strategy is authorized by the fact that the series method is effective to solution of equations of resiliency theory, thin covers and plates, and the bounded integral equations method is effective to research of infinite acoustic environment vibrations. Using this methodology let us describe the behavior of bounded resilient bodies near salient line, explore the pressure field in near and far zone, obtain the solution of the problem with defined precision in wide frequency range.

In the first chapter basic linear equations of motion and initial relations are considered, statement and solution of problem of scattering of flat acoustic wave by finite cylindrical cover bounded at the ends by resilient bottoms is done.

The improvement of convergence of obtained infinite systems of linear algebraic equations was held, using asymptotic properties of Fourier decomposition coefficients of the desired functions while solving the equations of cover theory.

The solutions of the problem of scattering of flat pressure wave by open-end free cylindrical cover, and also by tension cylindrical cover with defined inner tensions. In case of thin resilient cylindrical cover and round plates its motion is modeled by the linear Thimoshenko model.

Today the problem of learning of dynamical processes of acoustic fields interaction with thin-walled structures in liquids is of great interest. The theoretical exploration of these processes is strongly associated with building of effective analytic-numerical methods of mathematical physics equation solving, which are realized on today’s PC. For example, while analyzing wave diffraction on solid bodies and covers of revolution with fixed ends in  finite-element and spline-functions methods based on function approximation theory were applied.

We assume that scattering objects could easily move in boundless liquid environment or to be deformless.

In this work we propose the analytical method of dynamical problems solution of fixed oscillation of infinite acoustic environment interacted with resilient finite cylindrical bodies.

In the nowadays technics problems of oscillation of finite cylindrical covers, filled with liquid, their sound waves emission at vibration in infinite acoustic environment are actual.

While interaction of wave fields with resilient bounded bodies problems arise in different fields of technics, medicine, hydro-acoustics and geo-physics.

The analytical solutions of problems of scattering of acoustic waves by resilient bodies are directly related with developing methods of modelling of interaction processes of wave fields with underwater objects.

Let in the infinite ideal compressible rigid media the cylindrical cover of finite length bounded at the ends by bottoms is placed. Cover is filled with ideal compressible fluid inside with density $`\rho _0`$ and sound velocity $`c_0`$.

The considered cover is referred to cylindrical system of coordinates $`r^{},z^{},\mathrm{\Theta }`$ where axis $`oz^{}`$ is aligned with axis of cylindrical cover and plane $`z^{}=0`$ is equidistant to ends (see Fig. 1). Flat harmonic in time pressure wave falls to the cover (see Fig. 1):

$$P^\alpha (r,z,\mathrm{\Theta },\tau )=P_0\mathrm{exp}\left(i\omega \left(r\mathrm{cos}\mathrm{\Theta }\mathrm{sin}\phi ^{}+z\mathrm{cos}\phi ^{}\tau \right)\right)$$
(1)
In formula (1) and below following symbolism is used: $`r=\frac{r^{}}{R_0};z=\frac{z^{}}{R_0}`$ are nondimentional coordinates, $`\tau =\frac{ct}{R_0};\omega =\frac{\mathrm{\Omega }R_0}{c}`$ are nondimentional time and angular frequency, $`c`$ is sound velocity in the dissipating body environment, $`R_0`$ is radius of middle surface of cylindrical cover, $`t,\mathrm{\Omega }`$ are time and angular frequency, $`\phi ^{}`$ is angle between wave direction and axis of the cylindrical cover, $`P_0`$ is constant with the dimension of pressure. Time factor $`\mathrm{exp}\left(i\omega \tau \right)`$ is omitted below.

The problem of definition of dissipated pressure field outside the cover $`P^e(r,z,\mathrm{\Theta })`$ and in filler $`P^0(r,z,\mathrm{\Theta })`$ and components of cover movement vector $`\stackrel{}{U}`$ and plates $`\stackrel{}{W}`$ is reduced to solution of system of differential equations.

Dissipated acoustic pressure $`P^e(r,z,\mathrm{\Theta })`$ in unbounded acoustic medium is described by wave equation

$$\left(\mathrm{\Delta }+\varpi ^2\right)P^e(r,z,\theta )=0$$
(2)
where $`\mathrm{\Delta }`$ is three-dimensional Laplace operator.

Dissipated pressure in fluid which fills the cover $`P^0(r,z,\mathrm{\Theta })`$ satisfies the wave equation

$$\left(\mathrm{\Delta }+\frac{\varpi ^2}{\beta _0^2}\right)P^o(r,z,\theta )=0$$
(3)
where $`\beta _0^2=\frac{c^2}{c_0^2}.`$

The dynamics of thin resilient cylindrical cover is described by linear equation of Timoshenko cover theory which takes into account rotary inertia and deformation of transversed shift

$$L_{ij}u_j=g_1\delta _{i3};\left(L_{ij}=L_{ji};i,j=1,2,3,4,5\right)$$
(4)
where $`\delta _{ij}`$ is the Cronicler’s symbol,

$`L_{11}`$ $`=`$ $`\left(1+a^2\right)\left({\displaystyle \frac{^2}{\theta ^2}}\varkappa _1^2\right)+\sigma _1{\displaystyle \frac{^2}{z^2}}+{\displaystyle \frac{\varpi ^2}{\beta ^2}};`$
$`L_{12}`$ $`=`$ $`a^2\left({\displaystyle \frac{^2}{\theta ^2}}+\varkappa _1^2+\sigma _1{\displaystyle \frac{^2}{z^2}}+{\displaystyle \frac{\varpi ^2}{\beta ^2}}\right)+\varkappa _1^2`$
$`L_{22}`$ $`=`$ $`a^2\left({\displaystyle \frac{^2}{\theta ^2}}\varkappa _1^2+\sigma _1{\displaystyle \frac{^2}{z^2}}+{\displaystyle \frac{\varpi ^2}{\beta ^2}}\right)\varkappa _1^2;`$
$`L_{13}`$ $`=`$ $`\left(1+a^2\right)\left(1+\varkappa _1^2\right){\displaystyle \frac{}{\theta }};L_{14}=\sigma _2{\displaystyle \frac{^2}{z\theta }};L_{15}=L_{24}=0;`$
$`L_{23}`$ $`=`$ $`\left[\left(1+a^2\right)\varkappa _1^2+a^2\right]{\displaystyle \frac{}{\theta }};L_{25}=a^2\sigma _2{\displaystyle \frac{^2}{z\theta }};`$
$`L_{33}`$ $`=`$ $`\left(1+a^2\right)\left(\varkappa _1^2{\displaystyle \frac{^2}{\theta ^2}}+1\right)\varkappa _1^2{\displaystyle \frac{^2}{z^2}}{\displaystyle \frac{\varpi ^2}{\beta ^2}};L_{34}=\nu _{01}{\displaystyle \frac{}{z}};`$
$`L_{35}`$ $`=`$ $`\varkappa _1^2{\displaystyle \frac{}{z}};L_{44}=\left(1+a^2\right)\sigma _1{\displaystyle \frac{^2}{\theta ^2}}+{\displaystyle \frac{^2}{z^2}}+{\displaystyle \frac{\varpi ^2}{\beta ^2}};`$
$`L_{45}`$ $`=`$ $`a^2\left(\sigma _1{\displaystyle \frac{^2}{\theta ^2}}+{\displaystyle \frac{^2}{z^2}}+{\displaystyle \frac{\varpi ^2}{\beta ^2}}\right);\sigma _1={\displaystyle \frac{1\nu _{01}}{2}};\sigma _2={\displaystyle \frac{1+\nu _{01}}{2}};`$
$`L_{55}`$ $`=`$ $`a^2\left(\sigma _1{\displaystyle \frac{^2}{\theta ^2}}+{\displaystyle \frac{^2}{z^2}}+{\displaystyle \frac{\varpi ^2}{\beta ^2}}\right);\beta ^2={\displaystyle \frac{c_{10}^2}{c^2}};\varkappa _1={\displaystyle \frac{c_{20}}{c_{10}}};`$
$`g_1`$ $`=`$ $`{\displaystyle \frac{g_r}{\xi \rho c^2}};\xi ={\displaystyle \frac{\beta ^2h_1\rho _1}{R_0\rho }};`$
$$a^2=\frac{h_1^2}{R_0^2};c_{10}^2=\frac{E_1}{\rho \left(1\nu _{01}^2\right)};c_{20}^2=\frac{E_1k_{T1}}{2\rho _1\left(1+\nu _{01}\right)}.$$

$`E_1,\nu _{01},\rho _1,h_1`$ is Young’s modulus, Poisson’s ratio, density and thickness of cover material, $`\rho `$ is density of outside liquid medium, $`U_2,U_4`$ are turning angle of normal of cover’s middle surface in planes $`r\mathrm{\Theta }`$ and $`rz`$ appropriately; $`g_r`$ is radial component of outside force falls on middle surface area unit.

Differential equations of Timoshenko-Mindlin which describe the transverse motion of round resilient plates with shifted and rotary inertia taken into account and characterizing few first asymmetrical thickness mod oscillations have the following form

$$\left(\mathrm{\Delta }_0^2+\alpha _1\mathrm{\Delta }_0\alpha _2\right)W_z(r,\theta )=\left(\alpha _2^{}\alpha _1^{}\mathrm{\Delta }_0\right)g_z$$
(5)
$$\left(\mathrm{\Delta }_0+k_3^2\right)\mathrm{\Phi }(r,\theta )=0$$
(6)
where

$`\mathrm{\Delta }_0`$ $`=`$ $`{\displaystyle \frac{^2}{r^2}}+{\displaystyle \frac{1}{r}}{\displaystyle \frac{}{r}}+{\displaystyle \frac{1}{r^2}}{\displaystyle \frac{^2}{\theta ^2}};\alpha _1={\displaystyle \frac{mc^2\omega ^2\left(S+I\right)}{D_2}};`$
$`\alpha _2`$ $`=`$ $`{\displaystyle \frac{mc^2\omega ^2}{D_2}}\left(1{\displaystyle \frac{mc^2\omega ^2SI}{D_2}}\right);\alpha _1^{}={\displaystyle \frac{S}{D_2}};`$
$`m`$ $`=`$ $`{\displaystyle \frac{h_2\rho _2}{R_0}};\alpha _2^{}={\displaystyle \frac{1}{D_2}}\left(1{\displaystyle \frac{mc^2\omega ^2SI}{D_2}}\right);S={\displaystyle \frac{E_2h_2^2}{G_2R_0^2k_{T2}\left(1\nu _{02}^2\right)}};`$
$`I`$ $`=`$ $`{\displaystyle \frac{h_2^2}{12R_0^2}};D_2={\displaystyle \frac{E_2h_2^3}{12R_0^2\left(1\nu _{02}^2\right)}};`$
$$k_3^2=\frac{\varpi ^2}{\beta _3^2}\frac{2}{S\left(1\nu _{02}\right)};\beta _3^2=\frac{c^2}{c_{22}^2};c_{22}^2=\frac{E_2k_{T2}}{2\rho _2\left(1+\nu _{02}\right)};$$

$`g_z`$ \- normal component of outside force, which falls falls on middle surface area unit; $`E_2,\nu _{02},G_2,\rho _2,h_2`$ is Young’s modulus, Poisson’s ratio, shift ratio, density and plate thickness appropriately, $`K_{T2}`$ is numerical shift coefficient; $`\mathrm{\Phi }`$ is auxiliary function; $`C_{22}`$ is velocity of transversed waves in plate.

Let us define the radial and tangent components of motion vector of plates $`W_r,W_\mathrm{\Theta }`$ which characterize the flat stress state as follows

$$W_z=\frac{\phi }{r}+\frac{1}{r}\frac{\mathrm{\Psi }}{r};W_\theta =\frac{1}{r}\frac{\phi }{\theta }\frac{\mathrm{\Psi }}{r};$$
(7)
Here the scalar potential $`\phi `$ and non-zero component of vector potential $`\psi `$ satisfy the wave equations

$$\left(\mathrm{\Delta }_0+\frac{\omega ^2}{\beta _1^2}\right)\phi (r,\theta )=0;\left(\mathrm{\Delta }_0+\frac{\omega ^2}{\beta _2^2}\right)\mathrm{\Psi }(r,\theta )=0$$
(8)
where

$$\beta _1=\frac{c_{11}}{c};\beta _2=\frac{c_{21}}{c};c_{11}^2=\frac{E_2\left(1\nu _{02}\right)}{\rho \left(1+\nu _{02}\right)\left(12\nu _{02}\right)};c_{21}^2=\frac{E_2}{2\rho \left(1\nu _{02}\right)}.$$

The solution of differential equations (2)-(8) must satisfy the following boundary conditions

a) continuity translation on middle cover surface and plates

$$\left[\frac{}{r}\left(P^a+P^e\right)\right]_{r=1}=\rho \omega ^2c^2u_3;g_r=\left[P^aP^e+P^o\right]_{r=1}$$
(9)
$$\left[\frac{}{r}P^o\right]_{r=1}=\xi _{01}\beta _0^1\rho \omega ^2c^2u_3;\xi _{01}=\frac{\beta _0^2h_1\rho _1}{R_0\rho };$$

$$\left[\frac{}{z}\left(P^a+P^e\right)\right]_{z=l}=\rho \omega ^2c^2W_z;g_z=\left[P^aP^e+P^o\right]_{z=_{}^+l};$$

$$\left[\frac{}{z}P^o\right]_{z=l}=\xi _{01}\beta _0^1\rho \omega ^2c^2W_z;\xi _{02}=\frac{\beta _0^2h_2\rho _2}{R_0\rho };$$

b) continuity translation in junction place of cover and plates

$$\left[u_3\right]_{z=l}=\left[W_r\right]_{r=1};\left[u_5\right]_{z=l}=\left[W_z\right]_{r=1};\left[u_1\right]_{z=l}=\left[W_\theta \right]_{r=1};$$
(10)
c) continuity of turning angle $`\left(r=1,z=\pm l\right)`$

$$\left[u_2\right]_{z=l}=\left[\frac{1}{R_0}\frac{W_r}{r}\frac{W_\theta }{R_0}\right]_{r=1};\left[u_4\right]_{z=l}=\left[\frac{G^{}}{r}+\frac{1}{r}\frac{\mathrm{\Phi }}{r}\right]_{r=1}$$
(11)
d) continuity of forces $`\left(r=1,z=l\right)`$

$$B_1\left[\frac{u_5}{z}+\frac{\nu _{01}}{R_0}\left(\frac{u_1}{\theta }+u_3\right)\right]=\mathrm{\Lambda }_1^\stackrel{´}{}\left[\frac{}{r}\left(G^{}W_z\right)+\frac{1}{r}\frac{\mathrm{\Phi }}{r}\right];$$

$$\mathrm{\Lambda }_1^\stackrel{´}{}\left(u_4+\frac{u_3}{z}\right)=\frac{E_2h_2}{\left(1\nu _{02}^2\right)R_0}\left[\frac{W_r}{r}+\frac{\nu _{02}}{r}\left(W_z+\frac{W_\theta }{\theta }\right)\right];$$
(12)
$$B_1\frac{1\nu _{01}}{2}\left(\frac{1}{R_0}\frac{u_5}{\theta }+\frac{u_1}{z}\right)=\frac{E_2h_2}{\left(1\nu _{02}\right)R_0}\left(\frac{1}{r}\frac{W_z}{r}+\frac{W_\theta }{r}\frac{W_\theta }{r}\right);$$

e) continuity of moments $`\left(r=1,z=l\right)`$

$$D_1\left(\frac{u_4}{z}+\frac{\nu _{01}}{R_0}\frac{u_2}{\theta }\right)=D_2\left[\frac{^2G^{}}{r^2}+\nu _{02}\left(\frac{1}{r}\frac{G^{}}{r}+\frac{1}{r^2}\frac{^2G^{}}{r^2}\right)\frac{1\nu _{02}}{2}\frac{\mathrm{\Phi }}{r}\right];$$
(13)
$$D_1\frac{1\nu _{01}}{2}\left(\frac{u_2}{z}+\frac{1}{R_0}\frac{u_4}{\theta }+\frac{u_1}{z}\right)=$$

$$D_2\left(1\nu _{02}\right)\left\{\frac{^2G^{}}{r\theta }+\frac{1}{2}\left[\mathrm{\Delta }_0\mathrm{\Phi }2\left(\frac{1}{r}\frac{\mathrm{\Phi }}{r}+\frac{1}{r^2}\frac{^2\mathrm{\Phi }}{\theta ^2}\right)\right]\right\}$$

where

$$B_1=\frac{E_1h_1}{\left(1\nu _{02}^2\right)R_0};G^{}=W_z+S\mathrm{\Delta }_0W_z+\frac{S^2}{D_2}g_z;$$

$$\mathrm{\Lambda }_\sigma ^\stackrel{´}{}=\frac{k_{t\sigma }h_\sigma }{R_0};D_\sigma =\frac{E_\sigma h_\sigma ^3}{12\left(1\nu _{0\sigma }^2\right)R_0^3};\sigma =1,2.$$

Taking into account the symmetry relative to plane $`\mathrm{\Theta }=0`$ we find the solution of equations (2)-(8) in the Fourier series form:

$$u_1(z,\theta )=\underset{n=0}{\overset{\mathrm{}}{}}\underset{\nu =0}{\overset{\mathrm{}}{}}\left(u_{1n\nu }^+\mathrm{cos}\beta _\nu z+u_{1n\nu }^{}\mathrm{sin}\beta _\nu z\right)\mathrm{sin}n\theta ;$$

$$u_2(z,\theta )=\underset{n=0}{\overset{\mathrm{}}{}}\underset{\nu =0}{\overset{\mathrm{}}{}}\left(u_{2n\nu }^+\mathrm{cos}\beta _\nu z+u_{2n\nu }^{}\mathrm{sin}\beta _\nu z\right)\mathrm{sin}n\theta ;$$

$$u_3(z,\theta )=\underset{n=0}{\overset{\mathrm{}}{}}\underset{\nu =0}{\overset{\mathrm{}}{}}\left(u_{3n\nu }^+\mathrm{cos}\beta _\nu z+u_{3n\nu }^{}\mathrm{sin}\beta _\nu z\right)\mathrm{cos}n\theta ;$$

$$u_4(z,\theta )=\underset{n=0}{\overset{\mathrm{}}{}}\underset{\nu =0}{\overset{\mathrm{}}{}}\left(u_{4n\nu }^+\mathrm{sin}\beta _\nu z+u_{4n\nu }^{}\mathrm{cos}\beta _\nu z\right)\mathrm{cos}n\theta ;$$

$$u_5(z,\theta )=\underset{n=0}{\overset{\mathrm{}}{}}\underset{\nu =0}{\overset{\mathrm{}}{}}\left(u_{5n\nu }^+\mathrm{sin}\beta _\nu z+u_{5n\nu }^{}\mathrm{cos}\beta _\nu z\right)\mathrm{cos}n\theta $$
(14)
$$W_z(r,\theta )=\underset{n=0}{\overset{\mathrm{}}{}}W_{z,n}\left(r\right)\mathrm{cos}n\theta ;W_r(r,\theta )=\underset{n=0}{\overset{\mathrm{}}{}}W_{r,n}\left(r\right)\mathrm{cos}n\theta ;$$

$$\phi (r,\theta )=\underset{n=0}{\overset{\mathrm{}}{}}\phi _n\left(r\right)\mathrm{cos}n\theta ;\mathrm{\Psi }(r,\theta )=\underset{n=0}{\overset{\mathrm{}}{}}\mathrm{\Psi }_n\left(r\right)\mathrm{sin}n\theta $$

$$W_{\theta n}(r,\theta )=\underset{n=0}{\overset{\mathrm{}}{}}W_{\theta ,n}\left(r\right)\mathrm{sin}n\theta ;\mathrm{\Phi }(r,\theta )=\underset{n=0}{\overset{\mathrm{}}{}}\mathrm{\Phi }_n\left(r\right)\mathrm{sin}n\theta $$

$$P^e(r,z,\theta )=\underset{n=0}{\overset{\mathrm{}}{}}P_n^e(r,z)\mathrm{cos}n\theta ;P^0(r,z,\theta )=\underset{n=0}{\overset{\mathrm{}}{}}P_n^0(r,z)\mathrm{cos}n\theta $$

where

$$\beta _\nu =\frac{\nu \pi }{l};l=\frac{L}{R_0};$$

In relations (14) ”+” is the symmetric and ”-” is the non-symmetric components of general problem solution.

Let us exemplify the solution of wave equation (2) which satisfies the Zommerfield’s condition in the form of Helmholtz-Huygens integral

$$P^e(r,z,\theta )=_{\sigma _0}\left(P^e_{\sigma _0}\frac{G}{n_0}\frac{P^e}{n}_{\sigma _0}G\right)𝑑\sigma _0$$
(15)
where $`d\sigma _0`$ is element of dissipating surface area $`\sigma _0;`$

$`G=\left(4\pi R^{}\right)^1\mathrm{exp}\left(i\omega R^{}\right)`$ is the fundamental solution of Helmholtz equation; $`R^2=r^2+r_0^22rr_0\mathrm{cos}\left(\mathrm{\Theta }\mathrm{\Theta }_0\right)+\left(zz_0\right)^2`$ is the distance between point of observation $`(r,z,\mathrm{\Theta })`$ and point $`(r_0,z_0,\mathrm{\Theta }_0)`$ on the surface $`\sigma _0;`$ $`\frac{}{n_0}`$ is the derivative with respect to outside normal of surface $`\sigma _0`$ in point $`(r,z,\mathrm{\Theta });`$ $`P^e(r_0,z_0,\mathrm{\Theta }_0)`$ is the value of dissipated pressure on surface $`\sigma _0.`$

Let us arrange the defined representation of function $`G`$ in cylindrical coordinates  with divided coordinate $`z`$

$$G=\frac{1}{4\pi }\underset{m=0}{\overset{\mathrm{}}{}}ϵ_m\mathrm{cos}\left[m\left(\theta \theta _0\right)\right]_0^{\mathrm{}}\lambda J_m\left(\lambda r_0\right)J_m\left(\lambda r\right)\frac{e^{\varkappa zz_0}}{\varkappa }𝑑\lambda $$
(16)
where æ=$`\{\begin{array}{c}\sqrt{\lambda ^2\omega ^2},\text{ if }\lambda >\omega \\ i\sqrt{\omega ^2\lambda ^2},\text{ if }\omega >\lambda \end{array},`$ $`J_m\left(\right)`$ is Bessel’s function of the first kind; $`\epsilon _0=1;\epsilon _m=2;m1.`$

Let us assume on the surface of closed cylindrical cover the distributed pressure $`P^e(r_0,z_0,\mathrm{\Theta }_0)`$ and its formal derivatives are presented in the following form:

a) by Fourier’s series on the side surface of cylinder

$$P^e(z_0,\theta _0)=\underset{n=0}{\overset{\mathrm{}}{}}\underset{\nu =0}{\overset{\mathrm{}}{}}\left(f_{n\nu }^{+e}\mathrm{cos}\beta _\nu z_0+f_{n\nu }^e\mathrm{sin}\beta _\nu z_0\right)\mathrm{cos}n\theta $$
(17)
$$\frac{}{r}P^e(z_0,\theta _0)=\underset{n=0}{\overset{\mathrm{}}{}}\underset{\nu =0}{\overset{\mathrm{}}{}}\left(f_{n\nu }^{+e}\mathrm{cos}\beta _\nu z_0+f_{n\nu }^e\mathrm{sin}\beta _\nu z_0\right)\mathrm{cos}n\theta $$

b) by Fourier-Bessel’s series on plates ends $`z_0=\pm l`$

$$P^e(z_0,\theta _0)=\underset{n=0}{\overset{\mathrm{}}{}}\underset{\nu =0}{\overset{\mathrm{}}{}}g_{nj}^{e\pm }J_n\left(\gamma _{jn}r_0\right)\mathrm{cos}n\theta _0$$

$$\frac{}{z}P^e(z_0,\theta _0)=\underset{n=0}{\overset{\mathrm{}}{}}\underset{\nu =0}{\overset{\mathrm{}}{}}g_{nj}^{e\pm }J_n\left(\gamma _{jn}r_0\right)\mathrm{cos}n\theta _0$$
(18)
where $`\gamma _{jn}`$ is the zeroes of Bessel’s function derivative $`\left(J_n^{}\left(\theta _{jn}\right)=0\right)`$

$$g_{njz_{0=l}}^{e+}=g_{njz_{0=l}}^{e+};g_{njz_{0=l}}^e=g_{njz_{0=l}}^e;g_{njz_{0=l}}^{e+}=g_{njz_{0=l}}^{e+};g_{njz_{0=l}}^e=g_{njz_{0=l}}^e$$

Let us write the distributed pressure $`P^e`$ as the sum of symmetrical and antisymmetric components relative to $`z=0`$

$$P^e(r,z,\theta )=P^{e+}(r,z,\theta )+P^e(r,z,\theta )$$

Substituting the expansions (16), (17), (18) into (15) and integrating using (14) we get the following expressions for the distributed pressure to region

a)symmetrical component

$$P_n^{e+}(r,z)=\frac{\epsilon _n}{2}\underset{\nu =0}{\overset{\mathrm{}}{}}f_{n\nu }^{e+}_0^{\mathrm{}}\frac{\lambda ^2j_n^{}\left(\lambda \right)J_n\left(\lambda r\right)}{\varkappa ^2+\beta _\nu ^2}\left[\mathrm{cos}\beta _\nu z\left(1\right)^\nu ch\left(\varkappa z\right)e^{\varkappa l}\right]𝑑\lambda $$

$$\frac{\epsilon _n}{2}\underset{\nu =0}{\overset{\mathrm{}}{}}f_{n\nu }^{e+}_0^{\mathrm{}}\frac{\lambda j_n\left(\lambda \right)J_n\left(\lambda r\right)}{\varkappa ^2+\beta _\nu ^2}\left[\mathrm{cos}\beta _\nu z\left(1\right)^\nu ch\left(\varkappa z\right)e^{\varkappa l}\right]𝑑\lambda $$
(19)
$$\frac{\epsilon _n}{2}\underset{\nu =0}{\overset{\mathrm{}}{}}g_{nj}^{e+}J_n\left(\gamma _{jn}\right)_0^{\mathrm{}}\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda r\right)}{\gamma _{jn}^2\lambda ^2}e^{\varkappa l}sh(\varkappa z)d\lambda $$

$$\frac{\epsilon _n}{2}\underset{\nu =0}{\overset{\mathrm{}}{}}g_{nj}^{e+}J_n\left(\gamma _{jn}\right)_0^{\mathrm{}}\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda r\right)}{\varkappa \left(\gamma _{jn}^2\lambda ^2\right)}e^{\varkappa l}sh(\varkappa z)d\lambda .$$

b)antisymmetric component

$$P_n^e(r,z)=\frac{\epsilon _n}{2}\underset{\nu =0}{\overset{\mathrm{}}{}}f_{n\nu }^e_0^{\mathrm{}}\frac{\lambda ^2j_n^{}\left(\lambda \right)J_n\left(\lambda r\right)}{\varkappa \left(\varkappa ^2+\beta _\nu ^2\right)}\left[\varkappa \mathrm{sin}\beta _\nu z\left(1\right)^\nu \beta _\nu sh\left(\varkappa z\right)e^{\varkappa l}\right]𝑑\lambda $$

$$\frac{\epsilon _n}{2}\underset{\nu =0}{\overset{\mathrm{}}{}}f_{n\nu }^e_0^{\mathrm{}}\frac{\lambda J_n\left(\lambda \right)J_n\left(\lambda r\right)}{\varkappa \left(\varkappa ^2+\beta _\nu ^2\right)}\left[\varkappa \mathrm{sin}\beta _\nu z\left(1\right)^\nu \beta _\nu sh\left(\varkappa z\right)e^{\varkappa l}\right]𝑑\lambda \pm $$

$$\pm \frac{\epsilon _n}{2}\underset{\nu =0}{\overset{\mathrm{}}{}}g_{nj}^eJ_n\left(\gamma _{jn}\right)_0^{\mathrm{}}\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda r\right)}{\gamma _{jn}^2\lambda ^2}e^{\varkappa l}ch\left(\varkappa z\right)𝑑\lambda \pm $$
(20)
$$\pm \frac{\epsilon _n}{2}\underset{\nu =0}{\overset{\mathrm{}}{}}g_{nj}^eJ_n\left(\gamma _{jn}\right)_0^{\mathrm{}}\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda r\right)}{\varkappa \left(\gamma _{jn}^2\lambda ^2\right)}e^{\varkappa l}sh(\varkappa z)d\lambda .$$

where $`\epsilon _0=1;\epsilon _n=2;n1.`$

For the region $`\left(r0,\left|z\right|l,0\mathrm{\Theta }2\pi \right)`$ expression for $`P_n^{e+}(r,z)`$ and $`P_n^e(r,z)`$ will look

$$P_n^{e+}(r,z)=\frac{\epsilon _n}{2}\left\{\begin{array}{c}_{\nu =0}^{\mathrm{}}f_{n\nu }^{e+}_0^{\mathrm{}}\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda r\right)}{\varkappa ^2+\beta _\nu ^2}\left(1\right)^\nu sh\left(\varkappa l\right)e^{\varkappa z}𝑑\lambda \\ _{\nu =0}^{\mathrm{}}f_{n\nu }^{e+}_0^{\mathrm{}}\frac{\lambda J_n\left(\lambda \right)J_n\left(\lambda r\right)}{\varkappa ^2+\beta _\nu ^2}\left(1\right)^\nu sh\left(\varkappa l\right)e^{\varkappa z}𝑑\lambda +\\ +_{j=0}^{\mathrm{}}g_{nj}^{e+}J_n\left(\gamma _{jn}\right)_0^{\mathrm{}}\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda r\right)}{\gamma _{jn}{}_{}{}^{2}\lambda ^2}sh\left(\varkappa l\right)e^{\varkappa z}𝑑\lambda \\ _{j=0}^{\mathrm{}}g_{nj}^{e+}J_n\left(\gamma _{jn}\right)_0^{\mathrm{}}\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda r\right)}{\varkappa \left(\gamma _{jn}{}_{}{}^{2}\lambda ^2\right)}ch\left(\varkappa l\right)e^{\varkappa z}𝑑\lambda \end{array}\right\}$$
(21)
$$P_n^e(r,z)=\pm \frac{\epsilon _n}{2}\left\{\begin{array}{c}_{\nu =0}^{\mathrm{}}f_{n\nu }^e_0^{\mathrm{}}\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda r\right)}{\varkappa ^2+\beta _\nu ^2}\left(1\right)^\nu \beta _\nu sh\left(\varkappa l\right)e^{\varkappa z}𝑑\lambda \\ _{\nu =0}^{\mathrm{}}f_{n\nu }^e_0^{\mathrm{}}\frac{\lambda J_n\left(\lambda \right)J_n\left(\lambda r\right)}{\varkappa \left(\varkappa ^2+\beta _\nu ^2\right)}\left(1\right)^\nu \beta _\nu sh\left(\varkappa l\right)e^{\varkappa z}𝑑\lambda +\\ +_{j=0}^{\mathrm{}}g_{nj}^eJ_n\left(\gamma _{jn}\right)_0^{\mathrm{}}\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda r\right)}{\gamma _{jn}{}_{}{}^{2}\lambda ^2}ch\left(\varkappa l\right)e^{\varkappa z}𝑑\lambda \\ _{j=0}^{\mathrm{}}g_{nj}^eJ_n\left(\gamma _{jn}\right)_0^{\mathrm{}}\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda r\right)}{\varkappa \left(\gamma _{jn}{}_{}{}^{2}\lambda ^2\right)}sh\left(\varkappa l\right)e^{\varkappa z}𝑑\lambda \end{array}\right\}$$
(22)
\+ when $`z>l`$

\- when $`z<l.`$

Functions $`P_n^{e\pm }(r,z)`$ and $`P_n^{e\pm }(r_0,z_0)`$ and introduced in series (17)-(18) unknown coefficients $`f_{n\nu }^{e\pm },g_{nj}^{e\pm }`$ and functions $`P_n^{e\pm }(r_0,z_0)`$ are tied by relations

$$f_{n\nu }^{e+}=\frac{\epsilon _\nu }{l}_j^lP_n^{e+}(r,z)_{r=1}\mathrm{cos}\left(\beta _\nu z\right)dz;\epsilon _0=\frac{1}{2};\epsilon _\nu =1,\nu 1;$$

$$f_{n\nu }^e=\frac{1}{l}_j^lP_n^e(r,z)_{r=1}\mathrm{sin}\left(\beta _\nu z\right)dz;$$
(23)
$$g_{nj}^{e\pm }=\epsilon _{nj}_0^1P_n^{e\pm }(r,z)_{z=l}J_n\left(\gamma _{jn}r\right)rdr;ϵ_{nj}=\frac{2\gamma _{jn}^2}{\left(\gamma _{jn}^2n^2\right)J_n^2\left(\gamma _{jn}\right)}.$$

Substituting the expression (19),(20) into (23) we get algebraic equations

$$\underset{\nu =0}{\overset{\mathrm{}}{}}\left(f_{n\nu }^{e\pm }F_{n\nu \mu }^{1\pm }+f_{n\nu }^{e\pm }F_{n\nu \mu }^{1\pm }\right)+\underset{j=0}{\overset{\mathrm{}}{}}\left(g_{nj}^{e\pm }G_{nj\mu }^{1\pm }+g_{nj}^{e\pm }G_{nj\mu }^{1\pm }\right)=0;$$

$$\underset{\nu =0}{\overset{\mathrm{}}{}}\left(f_{n\nu }^{e\pm }F_{n\nu \mu }^{2\pm }+f_{n\nu }^{e\pm }F_{n\nu \mu }^{2\pm }\right)+\underset{j=0}{\overset{\mathrm{}}{}}\left(g_{nj}^{e\pm }G_{nj\mu }^{2\pm }+g_{nj}^{e\pm }G_{nj\mu }^{2\pm }\right)=0.$$
(24)
where

$$F_{N\nu \mu }^{1+}=_0^{\mathrm{}}\frac{\lambda J_n^2\left(\lambda \right)}{\left(\varkappa ^2+\beta _\nu ^2\right)}\left[\epsilon _\mu l\delta _{\nu \mu }\frac{\left(1\right)^{\nu +\mu }\varkappa \left(1e^{2\varkappa l}\right)}{\left(\varkappa ^2+\beta _\nu ^2\right)}\right]𝑑\lambda ;$$

$$F_{N\nu \mu }^{1+}=_0^{\mathrm{}}\frac{J_n^{}\left(\lambda \right)J_n\left(\lambda \right)}{\left(\varkappa ^2+\beta _\nu ^2\right)}\left[\epsilon _\mu l\delta _{\nu \mu }\left(\omega ^2\beta _\nu ^2\right)\frac{\left(1\right)^{\nu +\mu }\left(y_1y_2\right)}{\varkappa \left(\varkappa ^2+\beta _\nu ^2\right)}\right]𝑑\lambda $$

$$\epsilon _\mu l\delta _{\nu \mu }\left(\frac{n}{2}+1\right);y_1=\varkappa ^2\lambda ^2\left(1e^{2\varkappa l}\right);y_2=\left(\varkappa ^2+\beta _\nu ^2\right)^2;$$

$$G_{nj\mu }^{1+}=J_n\left(\gamma _{jn}\right)_0^{\mathrm{}}\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda \right)\left(1\right)^\mu \left(1e^{2\varkappa l}\right)}{\left(\gamma _{jn}^2\lambda ^2\right)\left(\varkappa ^2+\beta _\mu ^2\right)}𝑑\lambda $$

$$G_{nj\mu }^{1+}=J_n\left(\gamma _{jn}\right)_0^{\mathrm{}}\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda \right)\left(1\right)^\mu \varkappa \left(1e^{2\varkappa l}\right)}{\left(\gamma _{jn}^2\lambda ^2\right)\left(\varkappa ^2+\beta _\mu ^2\right)}𝑑\lambda ;$$

$$G_{nj\mu }^{2+}=J_n\left(\gamma _{jn}\right)J_n\left(\gamma _{\mu n}\right)_0^{\mathrm{}}\frac{\lambda ^3J_n^2\left(\lambda \right)\left(1e^{2\varkappa l}\right)}{\left(\gamma _{jn}^2\lambda ^2\right)\left(\gamma _{\mu n}^2\lambda ^2\right)}𝑑\lambda \frac{\delta _{j\mu }}{\epsilon _{n\mu }};$$

$$G_{nj\mu }^{2+}=J_n\left(\gamma _{jn}\right)J_n\left(\gamma _{\mu n}\right)_0^{\mathrm{}}\frac{\lambda ^3J_n^2\left(\lambda \right)\left(1+e^{2\varkappa l}\right)}{\varkappa \left(\gamma _{jn}^2\lambda ^2\right)\left(\gamma _{\mu n}^2\lambda ^2\right)}𝑑\lambda ;$$

$$F_{n\nu \mu }^{2+}=J_n\left(\gamma _{\mu n}\right)_0^{\mathrm{}}\frac{\lambda ^3J_n^2\left(\lambda \right)\left(1\right)^\nu \left(1e^{2\varkappa l}\right)}{\left(\varkappa ^2+\beta _\mu ^2\right)\left(\gamma _{\mu n}^2\lambda ^2\right)}𝑑\lambda ;$$

$$F_{n\nu \mu }^{2+}=J_n\left(\gamma _{\mu n}\right)_0^{\mathrm{}}\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda \right)\left(1\right)^\nu \left(1e^{2\varkappa l}\right)}{\left(\varkappa ^2+\beta _\mu ^2\right)\left(\gamma _{\mu n}^2\lambda ^2\right)}𝑑\lambda ;$$

$$F_{n\nu \mu }^1=_0^{\mathrm{}}\frac{J_n^{}\left(\lambda \right)J_n\left(\lambda \right)}{\varkappa \left(\varkappa ^2+\beta _\nu ^2\right)}\left[\varkappa l\delta _{\nu \mu }\left(\varpi ^2\beta _\nu ^2\right)\frac{\left(1\right)^{\nu +\mu }\beta _\nu \beta _\mu \lambda ^2\left(1e^{2\varkappa l}\right)}{\left(\varkappa ^2+\beta _\mu ^2\right)}\right]𝑑\lambda $$

$$\delta _{\nu \mu }l\left(\frac{n}{2}+1\right);$$

$$F_{n\nu \mu }^1=_0^{\mathrm{}}\frac{\lambda J_n^2\left(\lambda \right)}{\varkappa \left(\varkappa ^2+\beta _\nu ^2\right)}\left[\varkappa l\delta _{\nu \mu }+\frac{\left(1\right)^{\nu +\mu }\beta _\nu \left(1e^{2\varkappa l}\right)}{\left(\varkappa ^2+\beta _\mu ^2\right)}\right]𝑑\lambda ;$$

$$G_{nj\mu }^1=J_n\left(\gamma _{jn}\right)_0^{\mathrm{}}\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda \right)\left(1\right)^\mu \beta _\mu \left(1e^{2\varkappa l}\right)}{\left(\gamma _{jn}^2\lambda ^2\right)\left(\varkappa ^2+\beta _\mu ^2\right)}𝑑\lambda ;$$

$$G_{nj\mu }^1=J_n\left(\gamma _{jn}\right)_0^{\mathrm{}}\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda \right)\left(1\right)^\mu \beta _\mu \left(1e^{2\varkappa l}\right)}{\varkappa \left(\gamma _{jn}^2\lambda ^2\right)\left(\varkappa ^2+\beta _\mu ^2\right)}𝑑\lambda ;$$

$$F_{n\nu \mu }^2=J_n\left(\gamma _{\mu n}\right)_0^{\mathrm{}}\frac{\lambda ^3J_n^{}\left(\lambda \right)\left(1\right)^\nu \beta _\nu \left(1e^{2\varkappa l}\right)}{\varkappa \left(\varkappa ^2+\beta _\mu ^2\right)\left(\gamma _{\mu n}^2\lambda ^2\right)}𝑑\lambda ;$$

$$F_{nj\mu }^2=J_n\left(\gamma _{\mu n}\right)_0^{\mathrm{}}\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda \right)\left(1\right)^\mu \beta _\nu \left(1e^{2\varkappa l}\right)}{\varkappa \left(\gamma _{jn}^2\lambda ^2\right)\left(\varkappa ^2+\beta _\mu ^2\right)}𝑑\lambda ;$$

$$G_{nj\mu }^2=J_n\left(\gamma _{jn}\right)J_n\left(\gamma _{\mu n}\right)_0^{\mathrm{}}\frac{\lambda ^3J_n^2\left(\lambda \right)\left(1+e^{2\varkappa l}\right)}{\left(\gamma _{jn}^2\lambda ^2\right)\left(\gamma _{\mu n}^2\lambda ^2\right)}𝑑\lambda \frac{\delta _{j\mu }}{\epsilon _{n\mu }};$$

$$G_{nj\mu }^2=J_n\left(\gamma _{jn}\right)J_n\left(\gamma _{\mu n}\right)_0^{\mathrm{}}\frac{\lambda ^3J_n^2\left(\lambda \right)\left(1e^{2\varkappa l}\right)}{\varkappa \left(\gamma _{jn}^2\lambda ^2\right)\left(\gamma _{\mu n}^2\lambda ^2\right)}𝑑\lambda .$$

$`\delta _{\nu \mu }`$ is the Kronecker’s symbol.

Let us consider a question of derivatives representation from series (14) also by Fourier and Fourier-Bessel series, as during the solution of equations (3)-(8) we need to differentiate it.

Let us assume that even function $`f\left(z\right)`$ and odd function $`\phi \left(z\right)`$ and its derivatives could be expanded into Fourier series in $`[l,l]:`$

$$f^{\left(2k\right)}\left(z\right)=\underset{\nu =0}{\overset{\mathrm{}}{}}a_\nu ^{2k}\mathrm{cos}\beta _\nu z;\phi ^{\left(2k\right)}\left(z\right)=\underset{\nu =1}{\overset{\mathrm{}}{}}b_\nu ^{2k}\mathrm{sin}\beta _\nu z;$$

$$f^{\left(2k+1\right)}\left(z\right)=\underset{\nu =0}{\overset{\mathrm{}}{}}a_\nu ^{2k+1}\mathrm{sin}\beta _\nu z;\phi ^{\left(2k1\right)}\left(z\right)=\underset{\nu =1}{\overset{\mathrm{}}{}}b_\nu ^{2k1}\mathrm{cos}\beta _\nu z.$$
(25)
Coefficients $`a_\nu ^{2K},a_\nu ^{2K1},b_\nu ^{2K},b_\nu ^{2K1}`$ are defined by formulas

$$a_\nu ^{2k}=\frac{\epsilon _\nu }{2}_l^lf^{\left(2k\right)}\left(z\right)\mathrm{cos}\beta _\nu zdz;b_\nu ^{2k}=\frac{\epsilon _\nu }{2}_l^l\phi ^{\left(2k\right)}\left(z\right)\mathrm{sin}\beta _\nu zdz$$

$$a_\nu ^{2k+1}=\frac{\epsilon _\nu }{2}_l^lf^{\left(2k+1\right)}\left(z\right)\mathrm{sin}\beta _\nu zdz;b_\nu ^{2k1}=\frac{\epsilon _\nu }{2}_l^l\phi ^{\left(2k1\right)}\left(z\right)\mathrm{cos}\beta _\nu zdz$$
(26)
Integrating (26) by parts and taking into account the evenness and oddness of derivatives from $`f\left(z\right)`$ and $`\phi \left(z\right)`$ we find that

$$a_\nu ^{2k}=\frac{2\epsilon _\nu }{l}\left[\left(1\right)^{i+\nu +1}\beta _\nu ^{2i2}f^{\left(2k2i+1\right)}\left(l\right)+\left(1\right)^k\beta _\nu ^{2k}a_\nu ^0\right];$$

$$b_\nu ^{2k1}=\frac{2}{l}\left[\left(1\right)^{i+\nu +1}\beta _\nu ^{2i2}\phi ^{\left(2k2i\right)}\left(l\right)\left(1\right)^k\beta _\nu ^{2k1}b_\nu ^0\right];$$
(27)
$$a_\nu ^{2k+1}=\beta _\nu a_\nu ^{2k};b_\nu ^{2k}=\beta _\nu b_\nu ^{2k1};$$

During the solution of equations (3) we present the function $`P_n^0(r,z)`$ as the series

$$\underset{\nu =0}{\overset{\mathrm{}}{}}\left[P_{n\nu }^{o+}\left(r\right)\mathrm{cos}\beta _\nu z+P_{n\nu }^o\left(r\right)\mathrm{sin}\beta _\nu z\right]=P_n^o(r,z)$$
(28)
The second derivative from $`P_n^0(r,z)`$ taking into account (25), (27) we write as

$$\frac{^2}{z^2}P_n^o(r,z)=\underset{\nu =0}{\overset{\mathrm{}}{}}\left(\frac{}{z}P_n^{o+}(r,z)_{z=l}\frac{2\epsilon _\nu \left(1\right)^\nu }{l}\beta _\nu ^2P_{n\nu }^{o+}\left(r\right)\right)\mathrm{cos}\beta _\nu z+$$

$$+\underset{\nu =0}{\overset{\mathrm{}}{}}\left(\frac{2\left(1\right)^\nu \beta _\nu }{l}P_n^{o+}(r,z)_{z=l}\beta _\nu ^2P_{n\nu }^o\left(r\right)\right)\mathrm{sin}\beta _\nu z$$
(29)
Let the dissipated pressure and its derivatives on the interior side of the plates are defined as follows

$$P_n^{o\pm }(r,l)=\underset{j=0}{\overset{\mathrm{}}{}}g_{nj}^{o\pm }J_n\left(\gamma _{jn}r\right);\frac{}{z}P_n^{o\pm }(r,l)=\underset{j=0}{\overset{\mathrm{}}{}}g_{nj}^{o\pm }J_n\left(\gamma _{jn}r\right)$$
(30)
Taking into account (18)-(30) wave equation (3) for the symmetrical and antisymmetric components obtains the following form

$$\left[\frac{^2}{r^2}+\frac{1}{r}\frac{}{r}\left(\frac{n^2}{r^2}\varkappa _2^2\right)\right]P_{n\nu }^{o+}\left(r\right)=\underset{j=0}{\overset{\mathrm{}}{}}\frac{2\epsilon _\nu \left(1\right)^\nu }{l}g_{nj}^{o+}J_n\left(\gamma _{jn}r\right);$$
(31)
$$\left[\frac{^2}{r^2}+\frac{1}{r}\frac{}{r}\left(\frac{n^2}{r^2}\varkappa _2^2\right)\right]P_{n\nu }^o\left(r\right)=\underset{j=0}{\overset{\mathrm{}}{}}\frac{2\beta _\nu \left(1\right)^\nu }{l}g_{nj}^oJ_n\left(\gamma _{jn}r\right);$$

where æ$`{}_{2}{}^{2}=\frac{\omega ^2}{\beta _0^2\beta _\nu ^2};`$

General solutions of (31) are

$$P_{n\nu }^{o+}\left(r\right)=f_{n\nu }^{o+}J_n\left(\varkappa _2r\right)\underset{j=0}{\overset{\mathrm{}}{}}C_{n\nu j}^+g_{nj}^{o+}J_n\left(\gamma _{jn}r\right);$$
(32)
$$P_{n\nu }^o\left(r\right)=f_{n\nu }^oJ_n\left(\varkappa _2r\right)+\underset{j=0}{\overset{\mathrm{}}{}}C_{n\nu j}^{}g_{nj}^oJ_n\left(\gamma _{jn}r\right);$$

where

$$C_{n\nu j}^+=\frac{2\epsilon _\nu \left(1\right)^\nu }{l\left(\varkappa _2^2\gamma _{jn}^2\right)};C_{n\nu j}^{}=\frac{2\beta _\nu \left(1\right)^\nu }{l\left(\varkappa _2^2\gamma _{jn}^2\right)}$$

From (32) we obtain that pressure $`P_n^{0\pm }\left(r\right)`$ and $`\frac{}{r}P_n^{0\pm }\left(r\right)`$ on the interior surface of cylindrical cover $`\left(z=1\right)`$ are

$`{\displaystyle \frac{}{r}}P_n^{o\pm }\left(1\right)`$ $`=`$ $`f_{n\nu }^{o\pm }\left(\varkappa _2\right)J_n^{}\left(\varkappa _2\right).`$ (33)
$`P_{n\nu }^{o+}\left(1\right)`$ $`=`$ $`f_{n\nu }^{o+}J_n\left(\varkappa _2\right){\displaystyle \underset{j=0}{\overset{\mathrm{}}{}}}C_{n\nu j}^+g_{nj}^{o+}J_n\left(\gamma _{jn}\right)`$
$$P_{n\nu }^o\left(1\right)=f_{n\nu }^oJ_n\left(\varkappa _2\right)+\underset{j=0}{\overset{\mathrm{}}{}}C_{n\nu j}^{}g_{nj}^oJ_n\left(\gamma _{jn}\right)$$

Let us exemplify $`P_{n\nu }^{0\pm }`$ from (32) by Fourier-Bessel’s series and equate it to coefficients (30). We get

$$g_{nj}^o=\underset{\nu =0}{\overset{\mathrm{}}{}}f_{n\nu }^{o+}\frac{\left(1\right)^\nu \varkappa _2J_n^{}\left(\varkappa _2\right)J_n\left(\gamma _{jn}\right)\epsilon _{n\mu }}{\gamma _{jn}^2\varkappa _2^2}g_{nj}^{o+}\underset{\nu =0}{\overset{\mathrm{}}{}}\left(1\right)^\nu C_{n\nu \mu }^+;$$
(34)
$$g_{nj}^o=\underset{\nu =0}{\overset{\mathrm{}}{}}f_{n\nu }^{o+}\frac{\left(1\right)^\nu \varkappa _2J_n^{}\left(\varkappa _2\right)J_n\left(\gamma _{jn}\right)\epsilon _{n\mu }}{\gamma _{jn}^2\varkappa _2^2}g_{nj}^{o+}\underset{\nu =0}{\overset{\mathrm{}}{}}\left(1\right)^\nu C_{n\nu \mu }^{}.$$
(35)
Let us represent the falling pressure $`P^a(r,z,\mathrm{\Theta })`$ and ints normal derivatives on the surface of a cover and plates a the series analogical to (14), (18), (19):

$$P^a(r,z,\theta )=\underset{n=0}{\overset{\mathrm{}}{}}P_n^a(r,z)\mathrm{cos}n\theta ;$$

$$P_n^a(1,z)=\underset{\nu =0}{\overset{\mathrm{}}{}}\left(f_{n\nu }^{a+}\mathrm{cos}\beta _\nu z+f_{n\nu }^a\mathrm{sin}\beta _\nu z\right);P_n^{a\pm z}(r,l)=\underset{j=0}{\overset{\mathrm{}}{}}g_{nj}^{a\pm }J_n\left(\gamma _{jn}r\right);$$

$$\frac{}{r}P_n^a(r,z)_{r=1}=\underset{\nu =0}{\overset{\mathrm{}}{}}\left(f_{n\nu }^{a+}\mathrm{cos}\beta _\nu z+f_{n\nu }^a\mathrm{sin}\beta _\nu z\right);$$
(36)
$$\frac{}{z}P_n^{a\pm }(r,l)_{z=\pm j}=\underset{j=0}{\overset{\mathrm{}}{}}\left(g_{nj}^{a+}\pm g_{nj}^{a+}\right)J_n\left(\gamma _{jn}r\right);$$

Let us present the flat wave (1) as the expansion

$$P^a(r,z,\theta )=P_0\underset{j=0}{\overset{\mathrm{}}{}}\epsilon _ni^nJ_n\left(\omega \mathrm{sin}\phi ^{}r\right)\mathrm{exp}\left(i\omega z\mathrm{cos}\phi ^{}\right)\mathrm{cos}n\theta ;$$

and integrating using formulas that define the correspondent series components (36) we obtain

$$f_{n\nu }^{a+}=\frac{\left(1\right)^\nu 2\epsilon _ni^nJ_n\left(\omega \mathrm{sin}\phi ^{}\right)\omega \mathrm{cos}\phi ^{}}{\beta _\nu ^2\omega ^2\mathrm{cos}^2\phi ^{}}\mathrm{sin}\left(\omega l\mathrm{cos}\phi ^{}\right);$$

$$f_{n\nu }^a=\frac{2\left(1\right)^\nu \epsilon _ni^nJ_n\left(\omega \mathrm{sin}\phi ^{}\right)\beta _\nu }{\beta _\nu ^2\omega ^2\mathrm{cos}^2\phi ^{}}\mathrm{sin}\left(\omega l\mathrm{cos}\phi ^{}\right);$$

$$g_{n\nu }^a=\frac{\epsilon _{nj}\epsilon _ni^n\omega \mathrm{sin}\phi ^{}J_n^{}\left(\omega \mathrm{sin}\phi ^{}\right)J_n\left(\gamma _{jn}\right)}{\gamma _{jn}^2\omega ^2\mathrm{sin}^2\phi ^{}}\mathrm{cos}\left(\omega l\mathrm{cos}\phi ^{}\right);$$

$$g_{n\nu }^{a+}=\frac{\epsilon _ni^n\epsilon _{nj}\left(\omega \mathrm{sin}\phi ^{}\right)J_n\left(\gamma _{jn}\right)}{\gamma _{jn}^2\omega ^2\mathrm{sin}^2\phi ^{}}\mathrm{sin}\left(\omega l\mathrm{cos}\phi ^{}\right);$$

$$f_{n\nu }^{a+}=\epsilon _ni^n\omega \mathrm{sin}\phi ^{}J_n^{}\left(\omega \mathrm{sin}\phi ^{}\right)\frac{2\mathrm{cos}\phi ^{}\left(1\right)^\nu \omega }{\beta _\nu ^2\omega ^2\mathrm{cos}^2\phi ^{}}\mathrm{sin}\left(\omega l\mathrm{cos}\phi ^{}\right);$$

$$f_{n\nu }^a=\epsilon _ni^n\omega \mathrm{sin}\phi ^{}J_n^{}\left(\omega \mathrm{sin}\phi ^{}\right)\frac{2\beta _\nu }{\beta _\nu ^2\omega ^2\mathrm{cos}^2\phi ^{}}\mathrm{sin}\left(\omega l\mathrm{cos}\phi ^{}\right);$$

$$g_{n\nu }^{a+}=i\omega \mathrm{cos}\phi ^{}g_{nj}^{a+};g_{n\nu }^a=i\omega \mathrm{cos}\phi ^{}g_{nj}^a.$$

To find the solution of inhomogeneous equation (5) we define the auxiliary function

$$F^\pm =\mathrm{\Delta }_0W^\pm +\frac{S}{D_2}g_z^\pm ;$$
(37)
Taking into account (37) to define $`F^\pm (r,\mathrm{\Theta })`$ from (5) we get the following differential equation:

$$\left(\mathrm{\Delta }_0^2\alpha _1\mathrm{\Delta }_0\alpha _2\right)F^\pm (r,\theta )=\left(\beta _1^{}\beta _2^{}\right)g_z^\pm (r,z);$$
(38)
where

$$\beta _1^{}=\alpha _2^{}+\alpha _1^{}\frac{S}{D_2};\beta _2^{}=\alpha _2\frac{S}{D_2};$$

$$g_z^\pm (r,z)=\underset{n=0}{\overset{\mathrm{}}{}}\underset{j=0}{\overset{\mathrm{}}{}}\left(g_{nj}^{a\pm }+g_{nj}^{e\pm }g_{nj}^{o\pm }\right)J_n\left(\gamma _{jn}r\right)\mathrm{cos}n\theta ;$$

The solution of equation (38) we write as follows:

$$F^\pm (r,\theta )=\underset{n=0}{\overset{\mathrm{}}{}}F_n^\pm \left(r\right)\mathrm{cos}n\theta $$

$$F_n^\pm \left(r\right)=A_{1n}^\pm J_n\left(\gamma _1r\right)+A_{2n}^\pm J_n\left(\gamma _2r\right)+\underset{j=0}{\overset{\mathrm{}}{}}A_{3nj}g_{znj}J_n\left(\gamma _{jn}r\right)$$
(39)
$$\gamma _{1,2}=\sqrt{\frac{\alpha _1\pm \sqrt{\alpha _1^2+4\alpha _2}}{2}};A_{3nj}=\frac{\beta _1^{}\gamma _{jn}^2\beta _2^{}}{\gamma _{jn}^4\alpha _1\gamma _{jn}^2\alpha _2};$$

$$g_{znj}^\pm =g_{nj}^{a\pm }+g_{nj}^{e\pm }+g_{nj}^{o.\pm }.$$

$`A_{1n}^\pm ,A_{2n}^\pm `$ are unknown coefficients.

Using (39) and (37) to the normal shifts of the plates we obtain

$$W_n^\pm \left(r\right)=\frac{A_{1n}^\pm J_n\left(\gamma _1r\right)}{\gamma _1^2}+\frac{A_{2n}^\pm J_n\left(\gamma _2r\right)}{\gamma _2^2}+\underset{j=0}{\overset{\mathrm{}}{}}\frac{\left(A_{3nj}\frac{S}{D_2}\right)}{\gamma _{jn}^2}g_{znj}^\pm J_n\left(\gamma _{jn}r\right).$$
(40)
The solution of equations (6)-(8) will be rewritten as follows:

$$\mathrm{\Phi }_n^\pm \left(r\right)=\mathrm{\Phi }_{1n}^\pm J_n\left(k_1r\right);$$
(41)
$$\phi _n^\pm \left(r\right)=\phi _{1n}^\pm J_n\left(k_1r\right)$$
(42)
$$\mathrm{\Psi }_n^\pm \left(r\right)=\mathrm{\Psi }_{1n}^\pm J_n\left(k_2r\right)$$
(43)
Substituting expansions (14) into equation (4) taking into account (25), (27) we get the algebraic equations relatively to components of shift vector components:

$$a_{ij}^\pm u_{nj\nu }=b_{in}^\pm +\delta _{i3}g_{n\nu }^\pm ;a_{ij}^\pm =a_{ji}^\pm ;i,j=1,2,3,4,5.$$
(44)
where

$$a_{11}^\pm =\left(1+a^2\right)\left(n^2\sigma _1\beta _\nu ^2+\varkappa _1^2\right)\sigma _1\beta _\nu ^2\frac{\omega ^2}{\beta _\nu ^2};$$

$$a_{12}^\pm =a^2\left(n^2\sigma _1\beta _\nu ^2+\varkappa _1^2\right)\frac{\omega ^2}{\beta _\nu ^2};a_{13}^\pm =\left(1+a^2\right)\left(1+\varkappa _1^2\right)\left(\pm n\right);$$

$$a_{14}^\pm =\beta _\nu \sigma _2\left({}_{+}{}^{}n\right);a_{15}^\pm =0;a_{22}^\pm =a^2\left(n^2\sigma _1\beta _\nu ^2\varkappa _1^2+\frac{\omega ^2}{\beta _\nu ^2}\right)\varkappa _1^2;$$

$$a_{23}^\pm =\left[\left(1+a^2\right)\varkappa _1^2+a^2\right]\left(\pm n\right);a_{24}^\pm =0;a_{25}^\pm =a^2\sigma _2\beta _\nu ^2\left(\pm n\right);$$

$$a_{34}^\pm =\nu _{01}\beta _\nu ;a_{33}^\pm =\left(1+a^2\right)\left(1+\varkappa _1^2n^2\right)+\varkappa _1^2\beta _\nu ^2\frac{\omega ^2}{\beta _\nu ^2};a_{35}^\pm =_+^{}\beta _\nu \varkappa _1^2$$

$$a_{44}^\pm =\left(1+a^2\right)\sigma _1n^2\beta _\nu ^2+\frac{\omega ^2}{\beta _\nu ^2};a_{45}^\pm =a^2\left(n^2\sigma _1\beta _\nu ^2+\frac{\omega ^2}{\beta _\nu ^2}\right);$$

$$a_{55}^\pm =a^2\left(n^2\sigma _1\beta _\nu ^2+\frac{\omega ^2}{\beta _\nu ^2}\right);$$

$$b_{1n}^+=C_{11}^+u_{1n}^++C_{12}^+u_{2n}^++C_{14}^+u_{4n}^+;b_{2n}^+=C_{21}^+u_{1n}^++C_{22}^+u_{2n}^++C_{25}^+u_{5n}^+;$$

$$b_{3n}^+=C_{33}^+u_{3n}^++C_{34}^+u_{4n}^++C_{35}^+u_{5n}^++g_{rn\nu }^+;b_{4n}^+=C_{44}^+u_{4n}^++C_{45}^+u_{5n}^+;$$

$$b_{5n}^+=C_{54}^+u_{4n}^++C_{55}^+u_{5n}^+;b_{1n}^{}=C_{11}^{}u_{1n}^{}+C_{12}^{}u_{2n}^{}+C_{14}^{}u_{4n}^{};$$

$$b_{2n}^{}=C_{21}^{}u_{1n}^{}+C_{22}^{}u_{2n}^{}+C_{25}^{}u_{5n}^{};b_{3n}^{}=C_{33}^{}u_{3n}^{}+g_{rn\nu }^{};$$

$$b_{4n}^{}=C_{41}^{}u_{1n}^{}+C_{43}^{}u_{3n}^{}+C_{44}^{}u_{4n}^{}+C_{45}^{}u_{5n}^{};g_{rn\nu }^{}=f_{n\nu }^{e+}+f_{n\nu }^{a+}f_{n\nu }^{o+};$$

$$b_{5n}^{}=C_{52}^{}u_{2n}^{}+C_{53}^{}u_{3n}^{}+C_{54}^{}u_{4n}^{}+C_{55}^{}u_{5n}^{};C_{11}^+=\sigma _1\frac{2\epsilon _\nu \left(1\right)^\nu }{l};$$

$$C_{12}^+=a^2\sigma _1\frac{2\epsilon _\nu \left(1\right)^\nu }{l};C_{14}^+=\sigma _2n\frac{2\epsilon _\nu \left(1\right)^\nu }{l};C_{21}^+=C_{12}^+;$$

$$C_{22}^+=a^2\sigma _2\frac{2\epsilon _\nu \left(1\right)^\nu }{l};C_{25}^+=a^2\sigma _2n\frac{2\epsilon _\nu \left(1\right)^\nu }{l};C_{33}^+=\varkappa _1^2\frac{2\epsilon _\nu \left(1\right)^\nu }{l};$$

$$C_{34}^+=\nu _{01}\frac{2\epsilon _\nu \left(1\right)^\nu }{l};C_{35}^+=C_{33}^+;C_{44}^+=\beta _\nu \frac{2\epsilon _\nu \left(1\right)^\nu }{l};C_{45}^+=a^2\beta _\nu \frac{2\epsilon _\nu \left(1\right)^\nu }{l};$$

$$C_{54}^+=C_{55}^++C_{45}^+;C_{11}^{}=\sigma _1\frac{2\left(1\right)^\nu \beta _\nu }{l};C_{12}^{}=a^2\frac{2\left(1\right)^\nu \beta _\nu }{l};C_{21}^{}=C_{12}^{};$$

$$C_{33}^{}=\varkappa _1^2\frac{2\left(1\right)^\nu \beta _\nu }{l};C_{33}^{}=\sigma _2h\frac{2\left(1\right)^\nu \beta _\nu }{l};C_{53}^{}=C_{33}^{};C_{54}^{}=C_{55}^{}=C_{45}^{};$$

$$C_{52}^{}=a^2C_{41}^{};C_{45}^\pm =\pm C_{12}^\pm ;$$

From (44) we define $`U_{jn\nu }^\pm `$ as:

$$u_{nj\nu }^+=L_{1n\nu }^{j+}u_{1n}^++L_{2n\nu }^{j+}u_{2n}^++L_{3n\nu }^{j+}u_{3n}^++L_{4n\nu }^{j+}u_{4n}^++L_{5n\nu }^{j+}u_{5n}^++L_{6n\nu }^{j+}g_{nr\nu }^+;$$

$$u_{nj\nu }^{}=L_{1n\nu }^ju_{1n}^{}+L_{2n\nu }^ju_{2n}^{}+L_{3n\nu }^ju_{3n}^{}+L_{4n\nu }^ju_{4n}^{}+L_{5n\nu }^ju_{5n}^{}+L_{6n\nu }^jg_{nr\nu }^{};j=\overline{1,5};$$
(45)
where

$$L_{1n\nu }^{j+}=\frac{1}{\mathrm{\Delta }_5^+}\left(C_{11}^+A_{1j}^++C_{21}^+A_{2j}^+\right);L_{2n\nu }^{j+}=\frac{1}{\mathrm{\Delta }_5^+}\left(C_{12}^+A_{1j}^++C_{22}^+A_{2j}^+\right);$$

$$L_{3n\nu }^{j+}=\frac{1}{\mathrm{\Delta }_5^+}\left(C_{33}^+A_{3j}^+\right);L_{6n\nu }^{j+}=\frac{A_{3j}^+}{\mathrm{\Delta }_5^+};$$

$$L_{4n\nu }^{j+}=\frac{1}{\mathrm{\Delta }_5^+}\left(C_{44}^+A_{4j}^++C_{54}^+A_{5j}^++C_{14}^+A_{1j}^++C_{34}^+A_{3j}^+\right);$$

$$L_{5n\nu }^{j+}=\frac{1}{\mathrm{\Delta }_5^+}\left(C_{45}^+A_{4j}^++C_{55}^+A_{5j}^++C_{25}^+A_{2j}^++C_{35}^+A_{3j}^+\right);$$

$$L_{1n\nu }^j=\frac{1}{\mathrm{\Delta }_5^{}}\left(C_{11}^{}A_{1j}^{}+C_{21}^{}A_{2j}^{}+C_{41}^{}A_{2j}^{}\right);$$

$$L_{2n\nu }^j=\frac{1}{\mathrm{\Delta }_5^{}}\left(C_{12}^{}A_{1j}^{}+C_{22}^{}A_{2j}^{}+C_{52}^{}A_{5j}^{}\right);$$

$$L_{3n\nu }^j=\frac{1}{\mathrm{\Delta }_5^{}}\left(C_{33}^{}A_{3j}^{}+C_{43}^{}A_{4j}^{}+C_{53}^{}A_{5j}^{}\right);$$

$$L_{4n\nu }^j=\frac{1}{\mathrm{\Delta }_5^{}}\left(C_{44}^{}A_{4j}^{}+C_{54}^{}A_{5j}^{}\right);L_{5n\nu }^j=\frac{1}{\mathrm{\Delta }_5^{}}\left(C_{45}^{}A_{4j}^{}+C_{55}^{}A_{5j}^{}\right);$$

$$L_{6n\nu }^j=\frac{A_{3j}^{}}{\mathrm{\Delta }_5^{}};\mathrm{\Delta }_5^\pm =det\left|a_{ij}^\pm \right|;i,j=(1,2,3,4,5);$$

$`A_{i,j}^\pm `$ are algebraic adjuncts of elements $`a_{ij}^\pm .`$

To improve the convergence of the obtained solution (45) we use the asymptotic properties of expansion coefficients of Fourier series of the desired functions. Basing on (25), (27) we could define the following relations to the even functions $`f\left(z\right)`$ and odd functions $`\phi \left(z\right)`$ and its derivatives represented as Fourier series on $`[l,l],`$

a) for even function $`f\left(z\right)`$

$$f\left(z\right)=\underset{\nu =0}{\overset{N}{}}f_\nu \mathrm{cos}\beta _\nu z+\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{2\epsilon _\nu \left(1\right)^\nu }{l}\left[\frac{f^{}\left(l\right)}{\beta _\nu ^2}\frac{f^{}\left(l\right)}{\beta _\nu ^4}+\frac{f^\text{v}\left(l\right)}{\beta _\nu ^6}\right]\mathrm{cos}\beta _\nu z;$$
(46)
$$\frac{^2f\left(z\right)}{z^2}=\underset{\nu =0}{\overset{N}{}}\left[\beta _\nu f_\nu +\frac{2\epsilon _\nu \left(1\right)^\nu }{l}f^{}\left(l\right)\right]\mathrm{cos}\beta _\nu z+\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{2\epsilon _\nu \left(1\right)^\nu }{l}\left[\frac{f^{}\left(l\right)}{\beta _\nu ^4}\frac{f^\text{v}\left(l\right)}{\beta _\nu ^6}\right]\mathrm{cos}\beta _\nu z$$

$$\frac{^4f\left(z\right)}{z^4}=\underset{\nu =0}{\overset{N}{}}\left[\beta _\nu ^4f_\nu +\frac{2\epsilon _\nu \left(1\right)^\nu }{l}\left(\beta _\nu ^2f^{}\left(l\right)f^{}\left(l\right)\right)\right]\mathrm{cos}\beta _\nu z+$$

$$+\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{2\epsilon _\nu \left(1\right)^\nu }{l}\left[\frac{f^\text{v}\left(l\right)}{\beta _\nu ^2}\frac{f^\text{v}\left(l\right)}{\beta _\nu ^6}\right]\mathrm{cos}\beta _\nu z;f^{}\left(l\right)=\frac{^{}f\left(z\right)}{z}_{z=l};$$

$$f^{}\left(l\right)=\frac{^3f\left(z\right)}{z^3}_{z=l};f^\text{v}\left(l\right)=\frac{^5f\left(z\right)}{z^5}_{z=l};f^\text{v}\left(l\right)=\frac{^7f\left(z\right)}{z^7}_{z=l};$$

etc.

b) for odd function

$$\phi \left(z\right)=\underset{\nu =0}{\overset{N}{}}\phi _\nu \mathrm{sin}\beta _\nu z;$$

$$\frac{\phi \left(z\right)}{z}=\underset{\nu =1}{\overset{N}{}}\left[\beta _\nu \phi _\nu +\frac{2\left(1\right)^\nu }{l}\phi \left(l\right)\right]\mathrm{cos}\beta _\nu z+\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{2\left(1\right)^\nu }{l}\left[\frac{\phi ^{}\left(l\right)}{\beta _\nu }\frac{\phi ^\text{v}\left(l\right)}{\beta _\nu ^4}\right]\mathrm{cos}\beta _\nu z$$
(47)
$$\frac{^3\phi \left(z\right)}{z^3}=\underset{\nu =1}{\overset{N}{}}\left[\beta _\nu ^3\phi _\nu \frac{2\left(1\right)^\nu }{l}\left(\beta _\nu ^2\phi \left(l\right)\phi ^{}\left(l\right)\right)\right]\mathrm{cos}\beta _\nu z+$$

$$+\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{2\left(1\right)^\nu }{l}\left[\frac{\phi ^\text{v}\left(l\right)}{\beta _\nu ^2}\frac{\phi ^\text{v}\left(l\right)}{\beta _\nu ^4}\right]\mathrm{cos}\beta _\nu z;\frac{^2\phi \left(z\right)}{z^2}_{z=l}=\phi ^{}\left(l\right);$$

$$\frac{^5\phi \left(z\right)}{z^5}=\underset{\nu =1}{\overset{N}{}}\left[\beta _\nu ^5\phi _\nu +\frac{2\left(1\right)^\nu }{l}\left(\beta _\nu ^4\phi \left(l\right)\beta _\nu \phi ^{}\left(l\right)\right)\right]\mathrm{cos}\beta _\nu z+$$

$$+\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{2\left(1\right)^\nu }{l}\left[\frac{\phi ^\text{v}\left(l\right)}{\beta _\nu ^2}\frac{\phi ^\text{v}\left(l\right)}{\beta _\nu ^4}\right]\mathrm{cos}\beta _\nu z;\frac{^4\phi \left(z\right)}{z^4}_{z=l}=\phi ^\text{v}\left(l\right);\frac{^6\phi \left(z\right)}{z^6}_{z=l}=\phi ^\text{v}\left(l\right).$$

Let us introduce the symbolism

$$X_1^+=u_{1n}^+;X_2^+=u_{2n}^+;X_3^+=u_{3n}^+;X_4^+=u_{4n}^+;X_5^+=u_{5n}^+;$$

$$X_1^{}=u_{1n}^{};X_2^{}=u_{2n}^{};X_3^{}=u_{3n}^{};X_4^{}=u_{4n}^{};X_5^{}=u_{5n}^{};$$
(48)
$$X_6^\pm =A_{1n}^\pm ;X_7^\pm =A_{2n}^\pm ;X_8^\pm =\mathrm{\Phi }_{1n}^\pm ;X_9^\pm =\phi _{1n}^\pm ;X_{10}^\pm =\mathrm{\Psi }_{1n}^\pm .$$

Substituting obtained solutions of differential equations (2)-(8) into the boundary conditions (9)-(13) we get the system of equations relative to unknown constants (48) with asymptotic (46), (47) for $`U_{jn\nu }^\pm `$ taken into account:

$$\left(a^\pm \right)_{ij}\left\{X^\pm \right\}_i=\left\{F^\pm \right\}_i;$$
(49)
where

$$a_{11}^+=\underset{\nu =0}{\overset{N}{}}L_{1n\nu }^{3+}\left(1\right)^\nu ;a_{12}^+=\underset{\nu =0}{\overset{N}{}}L_{2n\nu }^{3+}\left(1\right)^\nu ;a_{13}^+=\underset{\nu =0}{\overset{N}{}}L_{3n\nu }^{3+}\left(1\right)^\nu +B_1^{};a_{14}^+=\underset{\nu =0}{\overset{N}{}}L_{4n\nu }^{3+}\left(1\right);$$

$$a_{15}^+=\underset{\nu =0}{\overset{N}{}}L_{5n\nu }^{3+}\left(1\right)^\nu ;a_{25}^+=1;a_{19}^\pm =k_1J_n^{}\left(k_1\right);a_{110}^\pm =nJ_n\left(k_2\right);B_1^{}=\frac{2l}{\pi }\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{1}{\nu ^2};$$

$$a_{31}^+=\underset{\nu =0}{\overset{N}{}}L_{1n\nu }^{1+}\left(1\right)^\nu +B_1^{};a_{32}^+=\underset{\nu =0}{\overset{N}{}}L_{2n\nu }^{1+}\left(1\right)^\nu ;a_{33}^+=\underset{\nu =0}{\overset{N}{}}L_{3n\nu }^{1+}\left(1\right)^\nu ;a_{34}^+=\underset{\nu =0}{\overset{N}{}}L_{4n\nu }^{1+}\left(1\right)^\nu ;$$

$$a_{26}^\pm =\frac{J_n\left(\gamma _1\right)}{\gamma _1};a_{27}^\pm =\frac{J_n\left(\gamma _2\right)}{\gamma _2};a_{310}^\pm =k_2J_n^{}\left(k_2\right);a_{39}^\pm =nJ_n\left(k_1\right);$$

$$a_{35}^+=\underset{\nu =0}{\overset{N}{}}L_{5n\nu }^{1+}\left(1\right)^\nu ;a_{41}^+=\underset{\nu =0}{\overset{N}{}}L_{1n\nu }^{2+}\left(1\right)^\nu ;a_{42}^+=\underset{\nu =0}{\overset{N}{}}L_{2n\nu }^{2+}\left(1\right)^\nu +B_1^{};a_{43}^+=\underset{\nu =0}{\overset{N}{}}L_{3n\nu }^{2+}\left(1\right)^\nu ;$$

$$a_{44}^+=\underset{\nu =0}{\overset{N}{}}L_{4n\nu }^{2+}\left(1\right)^\nu ;a_{45}^+=\underset{\nu =0}{\overset{N}{}}L_{3n\nu }^{2+}\left(1\right)^\nu ;a_{54}^+=1;a_{49}^\pm =\frac{n}{R_0}\left(k_1J_n^{}\left(k_1\right)J_n\left(k_1\right)\right);$$

$$a_{410}^\pm =\frac{nJ_n\left(k_2\right)}{R_0}\frac{k_2J_n^{}\left(k_2\right)}{R_0};a_{56}^\pm =\left(\frac{1}{\gamma _1}+S\gamma _1\right)J_n^{}\left(\gamma _1\right);a_{57}^\pm =\left(\frac{1}{\gamma _2}+S\gamma _2\right)J_n^{}\left(\gamma _2\right);$$

$$a_{61}^+=B_1\left\{\underset{\nu =1}{\overset{N}{}}L_{1n\nu }^{1+}\left(1\right)^\nu \beta _\nu +\frac{\nu _{01}}{R_0}n\left[\underset{\nu =0}{\overset{N}{}}L_{1n\nu }^{1+}\left(1\right)^\nu +B_1^{}\right]+\frac{\nu _{01}}{R_0}\underset{\nu =0}{\overset{N}{}}L_{1n\nu }^{3+}\left(1\right)^\nu \right\}$$

$$a_{58}^\pm =k_3J_n^{}\left(k_3\right);a_{62}^+=B_1\left\{\underset{\nu =1}{\overset{N}{}}L_{2n\nu }^{1+}\left(1\right)^\nu \beta _\nu +\frac{\nu _{01}}{R_0}n\underset{\nu =0}{\overset{N}{}}L_{1n\nu }^{1+}\left(1\right)^\nu +\frac{\nu _{01}}{R_0}\underset{\nu =0}{\overset{N}{}}L_{1n\nu }^{3+}\left(1\right)^\nu \right\}$$

$$a_{63}^+=B_1\left\{\underset{\nu =1}{\overset{N}{}}L_{3n\nu }^{5+}\left(1\right)^\nu \beta _\nu +\frac{\nu _{01}}{R_0}n\underset{\nu =0}{\overset{N}{}}L_{3n\nu }^{1+}\left(1\right)^\nu +\frac{\nu _{01}}{R_0}\left[\underset{\nu =0}{\overset{N}{}}L_{1n\nu }^{1+}\left(1\right)^\nu +B_1^{}\right]\right\}$$

$$a_{64}^+=B_1\left\{\underset{\nu =1}{\overset{N}{}}L_{4n\nu }^{5+}\left(1\right)^\nu \beta _\nu +\frac{\nu _{01}}{R_0}n\underset{\nu =0}{\overset{N}{}}L_{4n\nu }^{1+}\left(1\right)^\nu +\frac{\nu _{01}}{R_0}\underset{\nu =0}{\overset{N}{}}L_{4n\nu }^{3+}\left(1\right)^\nu \right\}$$

$$a_{65}^+=B_1\left\{\underset{\nu =1}{\overset{N}{}}\left[L_{5n\nu }^{5+}\left(1\right)^\nu \beta _\nu +\frac{2N}{l}\right]+\frac{\nu _{01}}{R_0}n\underset{\nu =0}{\overset{N}{}}L_{5n\nu }^{1+}\left(1\right)^\nu +\frac{\nu _{01}}{R_0}\underset{\nu =0}{\overset{N}{}}L_{5n\nu }^{3+}\left(1\right)^\nu \right\}$$

$$a_{66}^\pm =\mathrm{\Lambda }^{}S\gamma _1J_n^{}\left(\gamma _1\right);a_{67}^\pm =\mathrm{\Lambda }^{}S\gamma _2J_n^{}\left(\gamma _2\right);a_{68}^\pm =nJ_n\left(k_3\right);a_{74}^+=a_{75}^+=\mathrm{\Lambda }^{};$$

$$a_{79}^\pm =\frac{E_2h_2}{\left(1\nu _{02}^2\right)R_0}\left[k_1^2J_n^{}\left(k_1\right)+\nu _{02}k_1J_n^{}\left(k_1\right)\nu _{02}n^2J_n\left(k_1\right)\right];$$

$$a_{710}^\pm =\frac{E_2h_2}{\left(1\nu _{02}^2\right)R_0}\left[nk_2J_n^{}\left(k_2\right)+\nu _{02}nJ_n\left(k_2\right)\nu _{02}nJ_n^{}\left(k_2\right)\right];$$

$$a_{81}^+=B_1\frac{1\nu _{01}}{2};a_{85}^+=B_1\frac{1\nu _{01}}{2}\frac{n}{R_0};$$

$$a_{89}^\pm =\frac{E_2h_2}{\left(1\nu _{02}^2\right)R_0}\left[nk_1J_n^{}\left(k_1\right)nk_1J_n^{}\left(k_1\right)+nJ_n\left(k_1\right)\right];$$

$$a_{810}^\pm =\frac{E_2h_2}{\left(1\nu _{02}^2\right)R_0}\left[n^2J_n\left(k_2\right)k_2^2J_n^{}\left(k_2\right)+k_2J_n^{}\left(k_1\right)\right];$$

$$a_{92}^+=D_1\left\{\underset{\nu =1}{\overset{N}{}}L_{2n\nu }^{4+}\left(1\right)^\nu \beta _\nu +\frac{n\nu _{01}}{R_0}\left[\underset{\nu =0}{\overset{N}{}}L_{4n\nu }^{2+}\left(1\right)^\nu +B_1^{}\right]\right\};$$

$$a_{91}^+=D_1\left\{\underset{\nu =1}{\overset{N}{}}L_{1n\nu }^{4+}\left(1\right)^\nu \beta _\nu +\frac{n\nu _{01}}{R_0}\underset{\nu =0}{\overset{N}{}}L_{1n\nu }^{2+}\left(1\right)^\nu \right\};$$

$$a_{93}^+=D_1\left\{\underset{\nu =1}{\overset{N}{}}L_{3n\nu }^{4+}\left(1\right)^\nu \beta _\nu +\frac{n\nu _{01}}{R_0}\underset{\nu =0}{\overset{N}{}}L_{3n\nu }^{2+}\left(1\right)^\nu \right\};$$

$$a_{94}^+=D_1\left\{\underset{\nu =1}{\overset{N}{}}\left[L_{4n\nu }^{4+}\left(1\right)^\nu \beta _\nu +\frac{2N}{l}\right]+\frac{\nu _{01}}{R_0}n\underset{\nu =0}{\overset{N}{}}L_{4n\nu }^{2+}\left(1\right)^\nu \right\};$$

$$a_{95}^+=D_1\left\{\underset{\nu =1}{\overset{N}{}}L_{5n\nu }^{4+}\left(1\right)^\nu \beta _\nu +\frac{n\nu _{01}}{R_0}\underset{\nu =0}{\overset{N}{}}L_{5n\nu }^{2+}\left(1\right)^\nu \right\}$$

$$a_{96}^\pm =D_1\left[\left(\gamma _1^2S1\right)J_n^{}\left(\gamma _1\right)+\nu _{02}\left(S\gamma _1\frac{1}{\gamma _1}\right)J_n^{}\left(\gamma _1\right)+\nu _{02}n^2\left(\frac{1\gamma _1^2}{\gamma _1^2}\right)\right];$$

$$a_{97}^\pm =D_1\left[\left(\gamma _1^2S1\right)J_n^{}\left(\gamma _2\right)+\nu _{02}\left(S\gamma _2\frac{1}{\gamma _2}\right)J_n^{}\left(\gamma _2\right)+\nu _{02}n^2\left(\frac{1\gamma _2^2}{\gamma _2^2}\right)\right];$$

$$a_{98}^+=D_2\left(\frac{1\nu _{02}}{2}\right)k_3J_n\left(k_3\right);a_{101}^+=D_2\left(\frac{1\nu _{02}}{2R_0}\right);a_{102}^+=D_1\left(\frac{1\nu _{02}}{2}\right);$$

$$a_{104}^+=D_1\left(\frac{1\nu _{02}}{2}\right)\frac{n}{R_0};a_{106}^+=D_2n\left(1\nu _{02}\right)\left(S\gamma _1\frac{1}{\gamma _1}\right)J_n^{}\left(\gamma _1\right);$$

$$a_{107}^+=D_2n\left(1\nu _{02}\right)\left(S\gamma _2\frac{1}{\gamma _2}\right)J_n^{}\left(\gamma _2\right);a_{13}^{}=1;$$

$$a_{98}^+=D_2\left(\frac{1\nu _{02}}{2}\right)\left[2n^2J_n\left(k_3\right)+2k_3J_n^{}\left(k_3\right)k_3^2J_n\left(k_3\right)\right];$$

$$a_{21}^{}=\underset{\nu =1}{\overset{N}{}}L_{1n\nu }^5\left(1\right)^\nu ;a_{22}^{}=\underset{\nu =1}{\overset{N}{}}L_{2n\nu }^5\left(1\right)^\nu ;a_{23}^{}=\underset{\nu =1}{\overset{N}{}}L_{3n\nu }^5\left(1\right)^\nu ;a_{24}^{}=\underset{\nu =1}{\overset{N}{}}L_{4n\nu }^5\left(1\right)^\nu ;$$

$$a_{25}^{}=\underset{\nu =1}{\overset{N}{}}L_{5n\nu }^5\left(1\right)^\nu +B_1^{};a_{51}^{}=\underset{\nu =1}{\overset{N}{}}L_{1n\nu }^4\left(1\right)^\nu ;a_{52}^{}=\underset{\nu =1}{\overset{N}{}}L_{2n\nu }^4\left(1\right)^\nu ;a_{31}^{}=a_{42=}^{}1;$$

$$a_{51}^{}=\underset{\nu =1}{\overset{N}{}}L_{1n\nu }^4\left(1\right)^\nu ;a_{52}^{}=\underset{\nu =1}{\overset{N}{}}L_{2n\nu }^4\left(1\right)^\nu ;a_{53}^{}=\underset{\nu =1}{\overset{N}{}}L_{3n\nu }^4\left(1\right)^\nu ;a_{54}^{}=\underset{\nu =1}{\overset{N}{}}L_{4n\nu }^4\left(1\right)^\nu +B_1^{};$$

$$a_{55}^{}=\underset{\nu =1}{\overset{N}{}}L_{5n\nu }^4\left(1\right)^\nu ;a_{61}^{}=B_1\frac{\nu _{01}n}{R_0};a_{63}^{}=B_1\frac{\nu _{01}}{R_0};a_{65}^{}=B_1;$$

$$a_{71}^{}=\mathrm{\Lambda }_1\left\{\underset{\nu =1}{\overset{N}{}}L_{1n\nu }^3\left(1\right)^\nu \beta _\nu +\underset{\nu =0}{\overset{N}{}}L_{1n\nu }^4\left(1\right)^\nu \right\};$$

$$a_{72}^{}=\mathrm{\Lambda }_1\left\{\underset{\nu =1}{\overset{N}{}}L_{2n\nu }^3\left(1\right)^\nu \beta _\nu +\frac{n\nu _{01}}{R_0}\underset{\nu =0}{\overset{N}{}}L_{2n\nu }^4\left(1\right)^\nu \right\};$$

$$a_{73}^{}=\mathrm{\Lambda }_1\left\{\underset{\nu =1}{\overset{N}{}}L_{3n\nu }^3\left(1\right)^\nu \beta _\nu +\frac{2N}{l}+\underset{\nu =0}{\overset{N}{}}L_{3n\nu }^4\left(1\right)^\nu \right\}$$

$$a_{74}^{}=\mathrm{\Lambda }_1\left\{\underset{\nu =1}{\overset{N}{}}L_{4n\nu }^3\left(1\right)^\nu \beta _\nu +\left[\underset{\nu =0}{\overset{N}{}}L_{3n\nu }^4\left(1\right)^\nu +B_1^{}\right]\right\}$$

$$a_{75}^{}=\mathrm{\Lambda }_1\left\{\underset{\nu =1}{\overset{N}{}}L_{5n\nu }^3\left(1\right)^\nu \beta _\nu +\underset{\nu =0}{\overset{N}{}}L_{5n\nu }^4\left(1\right)^\nu \right\};$$

$$a_{76}^{}=\mathrm{\Lambda }_1\left\{\underset{\nu =1}{\overset{N}{}}\left[L_{5n\nu }^3\left(1\right)^\nu \beta _\nu +L_{5n\nu }^4\left(1\right)^\nu \right]g_{rn\nu }^{}\right\}$$

$$a_{81}^{}=B_1\frac{1\nu _{02}}{2}\left\{\frac{n}{R_0}\underset{\nu =1}{\overset{N}{}}L_{1n\nu }^5\left(1\right)^\nu +\left[\underset{\nu =0}{\overset{N}{}}L_{1n\nu }^1\left(1\right)^\nu \beta _\nu +\frac{2N}{l}\right]\right\};$$

$$a_{82}^{}=B_1\frac{1\nu _{02}}{2}\left\{\frac{n}{R_0}\underset{\nu =1}{\overset{N}{}}L_{2n\nu }^5\left(1\right)^\nu +\underset{\nu =0}{\overset{N}{}}L_{2n\nu }^1\left(1\right)^\nu \beta _\nu \right\};$$

$$a_{83}^{}=B_1\frac{1\nu _{02}}{2}\left\{\frac{n}{R_0}\underset{\nu =1}{\overset{N}{}}L_{3n\nu }^5\left(1\right)^\nu +\underset{\nu =0}{\overset{N}{}}L_{3n\nu }^1\left(1\right)^\nu \beta _\nu \right\};$$

$$a_{84}^{}=B_1\frac{1\nu _{02}}{2}\left\{\frac{n}{R_0}\underset{\nu =1}{\overset{N}{}}L_{4n\nu }^5\left(1\right)^\nu +\underset{\nu =1}{\overset{N}{}}L_{4n\nu }^1\left(1\right)^\nu \beta _\nu \right\};$$

$$a_{85}^{}=B_1\frac{1\nu _{02}}{2}\left\{\frac{n}{R_0}\left[\underset{\nu =1}{\overset{N}{}}L_{5n\nu }^5\left(1\right)^\nu +B_1^{}\right]+\underset{\nu =0}{\overset{N}{}}L_{5n\nu }^1\left(1\right)^\nu \beta _\nu \right\}$$

$$a_{101}^{}=D_1\frac{1\nu _{02}}{2}\left\{\frac{1}{R_0}\left[\underset{\nu =1}{\overset{N}{}}L_{1n\nu }^1\left(1\right)^\nu \beta _\nu +\frac{2N}{l}\right]+\underset{\nu =0}{\overset{N}{}}L_{1n\nu }^2\left(1\right)^\nu \beta _\nu \frac{n}{R_0}\underset{\nu =1}{\overset{N}{}}L_{4n\nu }^4\left(1\right)^\nu \right\}$$

$$a_{102}^{}=D_1\frac{1\nu _{02}}{2}\left\{\begin{array}{c}\left[_{\nu =1}^NL_{2n\nu }^2\left(1\right)^\nu \beta _\nu \frac{2N}{l}\right]\\ \frac{n}{R_0}_{\nu =0}^NL_{4n\nu }^4\left(1\right)^\nu \frac{1}{R_0}_{\nu =1}^NL_{4n\nu }^4\left(1\right)^\nu \beta _\nu ;\end{array}\right\}$$

$$a_{92}^{}=D_1\frac{\nu _{01}}{R_0}n;a_{94}^{}=D_1;$$

$$a_{103}^{}=D_1\frac{1\nu _{02}}{2}\left\{\underset{\nu =1}{\overset{N}{}}L_{3n\nu }^2\left(1\right)^\nu \beta _\nu +\frac{1}{R_0}\underset{\nu =0}{\overset{N}{}}L_{3n\nu }^1\left(1\right)^\nu \beta _\nu \frac{n}{R_0}\underset{\nu =1}{\overset{N}{}}L_{3n\nu }^4\left(1\right)^\nu \right\};$$

$$a_{104}^{}=D_1\frac{1\nu _{02}}{2}\{\frac{n}{R_0}[\underset{\nu =1}{\overset{N}{}}L_{4n\nu }^4(1)^\nu +B_1^{}]+\underset{\nu =1}{\overset{N}{}}L_{4n\nu }^2(1)^\nu \beta _\nu \frac{1}{R_0}\underset{\nu =1}{\overset{N}{}}L_{4n\nu }^1(1)^\nu \beta _\nu ;\}$$

$$a_{105}^{}=D_1\frac{1\nu _{02}}{2}\left\{\underset{\nu =1}{\overset{N}{}}L_{5n\nu }^2\left(1\right)^\nu \beta _\nu +\frac{1}{R_0}\underset{\nu =0}{\overset{N}{}}L_{5n\nu }^1\left(1\right)^\nu \beta _\nu \frac{n}{R_0}\underset{\nu =1}{\overset{N}{}}L_{5n\nu }^4\left(1\right)^\nu \right\};$$

$$F_i^\pm =\underset{\nu =0}{\overset{N}{}}\overline{C}_{in\nu }^\pm f_{rn\nu }^\pm +\underset{\nu =0}{\overset{N}{}}\overline{C}_{in\nu }^\pm g_{znj}^\pm ;i=1,2,3\mathrm{}\mathrm{}10$$

$$\overline{C}_{1n\nu }^+=L_{6n\nu }^{3+}\left(1\right)^\nu ;\overline{C}_{2nj}^\pm =A_{3nj}^{};\overline{C}_{3n\nu }^+=L_{6n\nu }^{1+}\left(1\right)^\nu ;\overline{C}_{4n\nu }^+=L_{6n\nu }^{2+}\left(1\right)^\nu ;$$

$$A_{3nj}^{}=\frac{A_{3nj}\frac{S}{D_2}}{\gamma _{jn}^2}J_n\left(\gamma _{jn}\right);\overline{C}_{6n\nu }^+=B_1\left[L_{6n\nu }^{5+}+\frac{\nu _{01}n}{R_0}L_{6n\nu }^{1+}+\frac{\nu _{01}}{R_0}L_{6n\nu }^{3+}\right];$$

$$\overline{C}_{9n\nu }^+=D_1\left[L_{6n\nu }^{4+}+\frac{\nu _{01}n}{R_0}L_{6n\nu }^{2+}\right];\overline{C}_{2n\nu }^{}=L_{6n\nu }^5\left(1\right)^\nu ;\overline{C}_{5n\nu }^{}=L_{6n\nu }^4\left(1\right)^\nu $$

$$\overline{C}_{9n\nu }^+=D_2\left[\frac{A_{3nj}\gamma _{jn}^2SA_{3nj}\frac{S}{D_2}}{\gamma _{jn}^2}\left(\gamma _{jn}^2n^2+1\right)\right]J_n\left(\gamma _{jn}\right);$$

$$\overline{C}_{7n\nu }^{}=\mathrm{\Lambda }_1\left\{L_{6n\nu }^3\left(1\right)^\nu \beta _\nu +\underset{\nu =0}{\overset{N}{}}L_{6n\nu }^4\left(1\right)^\nu \right\};$$

$$\overline{C}_{8n\nu }^{}=B_1\frac{1\nu _{02}}{2}\left\{\frac{n}{R_0}L_{6n\nu }^5\left(1\right)^\nu +L_{6n\nu }^1\left(1\right)^\nu \beta _\nu \right\};$$

$$\overline{C}_{10n\nu }^{}=D_1\frac{1\nu _{02}}{2}\left\{\frac{n}{R_0}L_{6n\nu }^4\left(1\right)^\nu +L_{6n\nu }^2\left(1\right)^\nu \beta _\nu +\frac{1}{R_0}L_{6n\nu }^1\left(1\right)^\nu \beta _\nu \right\};$$

Solution (49) will be

$$X_i=\underset{\nu =0}{\overset{N}{}}Z_{n\nu j}^\pm f_{rn\nu }^\pm +\underset{j=0}{\overset{N}{}}Z_{n\nu j}^\pm g_{rnj}^\pm ;\left(i=1,2,3\mathrm{}.10\right)$$
(50)
where

$$Z_{n\nu j}^\pm =\underset{k=1}{\overset{10}{}}\frac{\overline{C}_{kn\nu }^\pm \mathrm{\Delta }_{9ki}^\pm }{\mathrm{\Delta }_{10}^\pm };Z_{n\nu j}^\pm =\underset{k=1}{\overset{10}{}}\frac{\overline{C}_{kn\nu }^\pm \mathrm{\Delta }_{9ki}^\pm }{\mathrm{\Delta }_{10}^\pm }.$$

For $`U_{3n\nu }`$ and $`U_{zn\mu }`$ we obtain

$$u_{3n\mu }^\pm =L_{1n\mu }^{3\pm }X_1^\pm +L_{2n\mu }^{3\pm }X_2^\pm +L_{3n\mu }^{3\pm }X_3^\pm +L_{4n\mu }^{3\pm }X_4^\pm +L_{5n\mu }^{3\pm }X_5^\pm +L_{6n\mu }^{3\pm }f_{rn\mu }^\pm ;$$
(51)
$$W_{zn\mu }^\pm =L_{7n\mu }^\pm X_6^\pm +L_{8n\mu }^\pm X_7^\pm +L_{9n\mu }^\pm g_{zn\mu }^\pm ;$$

To improve the convergence of the obtained solution we will single-out the characteristic features in unknown functions. It is known that  approaching to angle circles on surface $`r=1`$ and $`z=l`$ expression for the pressure $`P^{e\pm }`$ and its derivatives must have the following singularities

a) for the pressure

$$P_n^{e+}\left(z_0\right)=A_n\sqrt[3]{\left(l^2z_0^2\right)^2};P_n^e\left(z_0\right)=B_nz_0\sqrt[3]{\left(l^2z_0^2\right)};P_n^{e\pm }\left(r_0\right)=D_n^\pm \sqrt[3]{(1r_0)^2;}$$
(52)
b) for the derivatives from pressure

$$\frac{P_n^{e+}\left(z_0r\right)}{r}_{r=1}=\frac{A_n^{}}{\sqrt[3]{l^2z_0^2}};\frac{P_n^e\left(z_0r\right)}{r}_{r=1}=\frac{B_n^{}z_0}{\sqrt[3]{l^2z_0^2}};$$

$$\frac{P_n^e\left(z_0r\right)}{z}_{z=l}=\frac{D_n^\pm r^n}{\sqrt[3]{1r_0^2}}.$$
(53)
Expanding (52), (53) into Fourier’s, Fourier-Bessel’s series corresponding to (15) we obtain

$$A_n\sqrt{\left(l^2z_0^2\right)^2}=\underset{\nu =0}{\overset{\mathrm{}}{}}f_{n\nu }^{s+}\mathrm{cos}\beta _\nu z_0;B_nz_0\sqrt[3]{\left(l^2z_0^2\right)}=\underset{\nu =0}{\overset{\mathrm{}}{}}f_{n\nu }^s\mathrm{sin}\beta _\nu z_0;$$

$$D_n^\pm \sqrt[3]{\left(1r_0^2\right)^2}=\underset{\nu =0}{\overset{\mathrm{}}{}}g_{nj}^sJ_n\left(\gamma _{jn}r\right);\frac{A_n^{}}{\sqrt[3]{l^2z_0^2}}=\underset{\nu =0}{\overset{\mathrm{}}{}}f_{n\nu }^{s+}\mathrm{cos}\beta _\nu z_0;$$
(54)
$$\frac{B_n^{}z_0}{\sqrt[3]{l^2z_0^2}}=\underset{\nu =0}{\overset{\mathrm{}}{}}f_{n\nu }^s\mathrm{sin}\beta _\nu z_0;\frac{D_n^\pm r^n}{\sqrt[3]{1r_0^2}}=\underset{\nu =0}{\overset{\mathrm{}}{}}g_{nj}^{s\pm }J_n\left(\gamma _{jn}r\right)$$

where

$$f_{n\nu }^{s+}=\frac{A_nl^{1/6}\mathrm{\Gamma }\left(5/3\right)\mathrm{\Gamma }\left(1/2\right)\mathrm{\Gamma }_{7/6}\left(\beta _\nu l\right)2^{7/6}}{\beta _\nu ^{7/6}};g_{nj}^{s\pm }=\frac{D_n^\pm \epsilon _{nj}\mathrm{\Gamma }\left(2/3\right)\mathrm{\Gamma }_{n+7/6}\left(\gamma _{jn}\right)}{\sqrt[3]{2\gamma _{jn}^2}};$$

$$g_{nj}^{s\pm }=\frac{D_n^\pm \epsilon _{nj}2^{2/3}\mathrm{\Gamma }\left(5/3\right)\mathrm{\Gamma }_{n+5/3}\left(\gamma _{jn}\right)}{\gamma _{jn}^{5/3}};f_{n\nu }^{s+}=\frac{A_n\mathrm{\Gamma }\left(2/3\right)\mathrm{\Gamma }\left(1/2\right)J_{1/6}\left(\beta _\nu l\right)2^{1/6}}{l^{5/6}\beta _\nu ^{7/6}};$$

$$f_{n\nu }^s=\frac{A_n\mathrm{\Gamma }\left(5/3\right)\mathrm{\Gamma }\left(1/2\right)J_{13/6}\left(\beta _\nu l\right)2^{7/6}l^{7/6}}{l^{5/6}\beta _\nu ^{7/6}};f_{n\nu }^s=\frac{A_n\mathrm{\Gamma }\left(2/3\right)\mathrm{\Gamma }\left(1/2\right)J_{1/6}\left(\beta _\nu l\right)2^{1/6}\beta _\nu l}{\beta _\nu ^{7/6}}.$$

Unknown coefficients $`f_{n\nu }^{e\pm },`$ $`f_{n\nu }^{e\pm },`$ $`g_{nj}^{e\pm },`$ $`g_{nj}^{e\pm }`$ at large $`\nu ,j\left(\nu ,j>N\right)`$ behave itself like coefficients of known expansions (54), i.e. we could write

$$f_{n\nu }^{e\pm }=\frac{f_{nN}^{e\pm }}{f_{nN}^{s\pm }}f_{n\nu }^{s\pm };f_{n\nu }^{e\pm }=\frac{f_{nN}^{e\pm }}{f_{nN}^\pm }f_{n\nu }^{s\pm };g_{n\nu }^{e\pm }=\frac{g_{nN}^{e\pm }}{g_{nN}^{s\pm }}g_{n\nu }^{s\pm };g_{n\nu }^{e\pm }=\frac{g_{nN}^{e\pm }}{g_{nN}^\pm }g_{n\nu }^{s\pm }.$$
(55)
Thus, if in the unlimited systems of algebraic linear equations (24) unknown coefficients are replaced with formulas (55) than elements of $`N`$-th columns in matrixes will be:

$$F_{n\mu N}^{1+}=\frac{\left(1\right)^\mu }{f_{nN}^{s+}}\left\{\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda \right)}{\left(\varkappa ^2+\beta _\mu ^2\right)}\varkappa \left(1e^{2\varkappa l}\right)\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{\left(1\right)^\nu f_{n\nu }^{s+}}{\varkappa ^2+\beta _\nu ^2}\right\}𝑑\lambda ;$$

$$F_{n\mu N}^{1+}=\frac{\left(1\right)^\mu }{f_{nN}^{s+}}\left\{\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda \right)}{\left(\varkappa ^2+\beta _\mu ^2\right)}\varkappa \left(1e^{2\varkappa l}\right)\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{\left(1\right)^\nu f_{n\nu }^{s+}}{\varkappa ^2+\beta _\nu ^2}\right\}𝑑\lambda ;$$

$$G_{n\mu N}^{1+}=\frac{\left(1\right)^\mu }{g_{nN}^{s+}}\left\{\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda \right)}{\left(\varkappa ^2+\beta _\mu ^2\right)}\varkappa \left(1e^{2\varkappa l}\right)\underset{j=N+1}{\overset{\mathrm{}}{}}\frac{g_{nj}^{s+}J_n\left(\gamma _{jn}\right)}{\gamma _{jn}^2\lambda ^2}\right\}𝑑\lambda ;$$

$$G_{n\mu N}^{1+}=\frac{\left(1\right)^\mu }{g_{nN}^{s+}}\left\{\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda \right)}{\left(\varkappa ^2+\beta _\mu ^2\right)}\varkappa \left(1e^{2\varkappa l}\right)\underset{j=N+1}{\overset{\mathrm{}}{}}\frac{g_{nj}^{s+}J_n\left(\gamma _{jn}\right)}{\gamma _{jn}^2\lambda ^2}\right\}𝑑\lambda ;$$

$$F_{n\mu N}^{2+}=\frac{J_n\left(\gamma _{jn}\right)}{f_{nN}^{s+}}\left\{\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda \right)}{\gamma _{jn}^2\lambda ^2}\left(1e^{2\varkappa l}\right)\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{\left(1\right)^\nu f_{n\nu }^{s+}}{\varkappa ^2+\beta _\nu ^2}\right\}𝑑\lambda ;$$

$$F_{n\mu N}^{2+}=\frac{J_n\left(\gamma _{jn}\right)}{f_{nN}^{s+}}\left\{\frac{\lambda ^3J_n^2\left(\lambda \right)}{\gamma _{jn}^2\lambda ^2}\left(1e^{2\varkappa l}\right)\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{\left(1\right)^\nu f_{n\nu }^{s+}}{\varkappa ^2+\beta _\nu ^2}\right\}𝑑\lambda ;$$
(56)
$$G_{n\mu N}^{2+}=\frac{J_n\left(\gamma _{jn}\right)}{g_{nN}^{s+}}\left\{\frac{\lambda ^3J_n^2\left(\lambda \right)}{\varkappa \left(\gamma _{jn}^2\lambda ^2\right)}\left(1+e^{2\varkappa l}\right)\underset{j=N+1}{\overset{\mathrm{}}{}}\frac{g_{nj}^{s+}J_n\left(\gamma _{jn}\right)}{\gamma _{jn}^2\lambda ^2}\right\}𝑑\lambda ;$$

$$G_{n\mu N}^{2+}=\frac{J_n\left(\gamma _{jn}\right)}{g_{nN}^{s+}}\left\{\frac{\lambda ^3J_n^2\left(\lambda \right)}{\left(\gamma _{jn}^2\lambda ^2\right)}\left(1e^{2\varkappa l}\right)\underset{j=N+1}{\overset{\mathrm{}}{}}\frac{g_{nj}^{s+}J_n\left(\gamma _{jn}\right)}{\gamma _{jn}^2\lambda ^2}\right\}𝑑\lambda ;$$

$$F_{n\mu N}^1=\frac{\left(1\right)^\mu \beta _\mu }{f_{nN}^s}\left\{\frac{\lambda J_n^2\left(\lambda \right)}{\varkappa \left(\varkappa ^2+\beta _\mu ^2\right)}\left(1e^{2\varkappa l}\right)\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{\left(1\right)^\nu \beta _\nu f_{n\nu }^s}{\varkappa ^2+\beta _\nu ^2}\right\}𝑑\lambda ;$$

$$F_{n\mu N}^1=\frac{\left(1\right)^\mu \beta _\mu }{f_{nN}^s}\left\{\frac{\lambda J_n^{}\left(\lambda \right)J_n\left(\lambda \right)}{\varkappa \left(\varkappa ^2+\beta _\mu ^2\right)}\left(1e^{2\varkappa l}\right)\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{\left(1\right)^\nu \beta _\nu f_{n\nu }^s}{\varkappa ^2+\beta _\nu ^2}\right\}𝑑\lambda ;$$

$$G_{n\mu N}^1=\frac{\left(1\right)^\mu \beta _\mu }{g_{nN}^s}\left\{\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda \right)}{\varkappa \left(\varkappa ^2+\beta _\mu ^2\right)}\left(1e^{2\varkappa l}\right)\underset{j=N+1}{\overset{\mathrm{}}{}}\frac{g_{nj}^sJ_n\left(\gamma _{jn}\right)}{\gamma _{jn}^2\lambda ^2}\right\}𝑑\lambda ;$$

$$G_{n\mu N}^1=\frac{\left(1\right)^\mu \beta _\mu }{g_{nN}^s}\left\{\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda \right)}{\left(\varkappa ^2+\beta _\mu ^2\right)}\left(1e^{2\varkappa l}\right)\underset{j=N+1}{\overset{\mathrm{}}{}}\frac{g_{nj}^sJ_n\left(\gamma _{jn}\right)}{\gamma _{jn}^2\lambda ^2}\right\}𝑑\lambda ;$$

$$F_{n\mu N}^2=\frac{J_n\left(\gamma _{jn}\right)}{f_{nN}^s}\left\{\frac{\lambda ^2J_n^{}\left(\lambda \right)J_n\left(\lambda \right)}{\varkappa \left(\gamma _{jn}^2\lambda ^2\right)}\left(1e^{2\varkappa l}\right)\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{\left(1\right)^\nu \beta _\nu f_{n\nu }^s}{\varkappa ^2+\beta _\nu ^2}\right\}𝑑\lambda ;$$

$$F_{n\mu N}^2=\frac{J_n\left(\gamma _{jn}\right)}{f_{nN}^s}\left\{\frac{\lambda ^3J_n^2\left(\lambda \right)}{\varkappa \left(\gamma _{jn}^2\lambda ^2\right)}\left(1e^{2\varkappa l}\right)\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{\left(1\right)^\nu \beta _\nu f_{n\nu }^s}{\varkappa ^2+\beta _\nu ^2}\right\}𝑑\lambda ;$$

$$G_{n\mu N}^2=\frac{J_n\left(\gamma _{jn}\right)}{g_{nN}^s}\left\{\frac{\lambda ^3J_n^2\left(\lambda \right)}{\varkappa \left(\gamma _{jn}^2\lambda ^2\right)}\left(1e^{2\varkappa l}\right)\underset{j=N+1}{\overset{\mathrm{}}{}}\frac{g_{nj}^sJ_n\left(\gamma _{jn}\right)}{\gamma _{jn}^2\lambda ^2}\right\}𝑑\lambda ;$$

$$G_{n\mu N}^2=\frac{J_n\left(\gamma _{jn}\right)}{g_{nN}^s}\left\{\frac{\lambda ^3J_n^2\left(\lambda \right)}{\left(\gamma _{jn}^2\lambda ^2\right)}\left(1+e^{2\varkappa l}\right)\underset{j=N+1}{\overset{\mathrm{}}{}}\frac{g_{nj}^sJ_n\left(\gamma _{jn}\right)}{\gamma _{jn}^2\lambda ^2}\right\}𝑑\lambda .$$

where

$$\frac{2^{1/6}l^{1/6}\mathrm{\Gamma }\left(2/3\right)\mathrm{\Gamma }\left(1/2\right)I_{1/6}\left(\varkappa l\right)}{\varkappa ^{7/6}sh\left(\varkappa l\right)}\underset{\nu =0}{\overset{N1}{}}\frac{f_{n\nu }^{s+}\left(1\right)^\nu }{\varkappa ^2+\beta _\nu ^2}=\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{f_{n\nu }^{s+}\left(1\right)^\nu }{\varkappa ^2+\beta _\nu ^2};$$

$$\frac{2^{7/6}l^{7/6}\mathrm{\Gamma }\left(5/3\right)\mathrm{\Gamma }\left(1/2\right)I_{7/6}\left(\varkappa l\right)}{\varkappa ^{13/6}sh\left(\varkappa l\right)}\underset{\nu =0}{\overset{N1}{}}\frac{f_{n\nu }^{s+}\left(1\right)^\nu }{\varkappa ^2+\beta _\nu ^2}=\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{f_{n\nu }^{s+}\left(1\right)^\nu }{\varkappa ^2+\beta _\nu ^2};$$

$$\frac{2^{2/3}\lambda ^{1/3}\mathrm{\Gamma }\left(2/3\right)J_{n+5/3}\left(\lambda \right)}{\left(\lambda ^2n^2\right)J_n^2\left(\lambda \right)J_n^{}\left(\lambda \right)}\underset{\nu =0}{\overset{N1}{}}\frac{g_{n\nu }^{s\pm }J_n\left(\gamma _{jn}\right)}{\gamma _{jn}^2\lambda ^2}=\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{g_{n\nu }^{s\pm }J_n\left(\gamma _{jn}\right)}{\gamma _{jn}^2\lambda ^2};$$

$$\frac{2^{2/3}\lambda ^{2/3}\mathrm{\Gamma }\left(5/3\right)\mathrm{\Gamma }\left(1/2\right)J_{n+2/3}\left(\lambda \right)}{\left(\lambda ^2n^2\right)J_n^2\left(\lambda \right)J_n^{}\left(\lambda \right)}\underset{\nu =0}{\overset{N1}{}}\frac{g_{n\nu }^{s\pm }J_n\left(\gamma _{jn}\right)}{\gamma _{jn}^2\lambda ^2}=\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{g_{n\nu }^{s\pm }J_n\left(\gamma _{jn}\right)}{\gamma _{jn}^2\lambda ^2};$$

$$\frac{2^{1/3}\mathrm{\Gamma }\left(2/3\right)\mathrm{\Gamma }\left(1/2\right)\left[\varkappa lI_{1/6}^{}\left(\varkappa l\right)1/6I_{1/6}\left(\varkappa l\right)\right]}{\varkappa ^{1/6}sh\left(\varkappa l\right)}\underset{\nu =0}{\overset{N1}{}}\frac{f_{n\nu }^s\left(1\right)^\nu \beta _\nu }{\varkappa ^2+\beta _\nu ^2}=\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{f_{n\nu }^s\left(1\right)^\nu \beta _\nu }{\varkappa ^2+\beta _\nu ^2}$$

$$\frac{2^{7/6}l^{13/6}\mathrm{\Gamma }\left(5/3\right)\mathrm{\Gamma }\left(1/2\right)\left[\varkappa lI_{7/6}^{}\left(\varkappa l\right)7/6I_{7/6}\left(\varkappa l\right)\right]}{\varkappa ^{13/6}sh\left(\varkappa l\right)}\underset{\nu =0}{\overset{N1}{}}\frac{f_{n\nu }^{s+}\left(1\right)^\nu \beta _\nu }{\varkappa ^2+\beta _\nu ^2}=\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{f_{n\nu }^{s+}\left(1\right)^\nu \beta _\nu }{\varkappa ^2+\beta _\nu ^2}$$

Substituting the obtained values into condition of hydro-resilient contact (9) and taking into account (24), (34), (35), (56) we obtain such unlimited systems written in matrix form to evaluate the unknown coefficients

$$f_{n\nu }^{e\pm },f_{n\nu }^{e\pm },g_{nj}^{e\pm },g_{nj}^{e\pm },f_{n\nu }^{o\pm },f_{n\nu }^{o\pm },g_{nj}^{o\pm },g_{nj}^{o\pm };(\nu ,j=0,1,2\mathrm{}\mathrm{}N);$$

$$\underset{N\times N}{\left(F_{1n}^{1\pm }\right)_{\nu \mu }}\underset{1\times N}{\left\{f_n^{e\pm }\right\}_\nu }+\underset{N\times N}{\left(F_{2n}^{1\pm }\right)_{\nu \mu }}\underset{1\times N}{\left\{f_n^{e\pm }\right\}_\nu }+\underset{N\times N}{\left(G_{1n}^{1\pm }\right)_{\nu \mu }}\underset{1\times N}{\left\{g_n^{e\pm }\right\}_j}+\underset{N\times N}{\left(G_{2n}^{1\pm }\right)_{\nu \mu }}\underset{1\times N}{\left\{g_n^{e\pm }\right\}_j}=0;$$
(57)
$$\underset{N\times N}{\left(F_{1n}^{2\pm }\right)_{\nu \mu }}\underset{1\times N}{\left\{f_n^{e\pm }\right\}_\nu }+\underset{N\times N}{\left(F_{2n}^{2\pm }\right)_{\nu \mu }}\underset{1\times N}{\left\{f_n^{e\pm }\right\}_\nu }+\underset{N\times N}{\left(G_{1n}^{2\pm }\right)_{\nu \mu }}\underset{1\times N}{\left\{g_n^{e\pm }\right\}_j}+\underset{N\times N}{\left(G_{2n}^{2\pm }\right)_{\nu \mu }}\underset{1\times N}{\left\{g_n^{e\pm }\right\}_j}=0$$
(58)
$$\underset{N\times N}{\left(R_{1n}^\pm \right)_{\nu \mu }}\underset{1\times N}{\left\{f_{rn}^\pm \right\}_\nu }+\underset{}{\underset{N\times N}{\left(R_{2n}^\pm \right)_{\nu \mu }}\underset{1\times N}{\left\{g_{zn}^\pm \right\}_j}}=\underset{1\times N}{\left\{f_n^{e\pm }\right\}_\nu }+\underset{1\times N}{\left\{f_n^{a\pm }\right\}_\nu };$$
(59)
$$\underset{N\times N}{\left(R_{3n}^\pm \right)_{\nu \mu }}\underset{1\times N}{\left\{f_{rn}^\pm \right\}_\nu }+\underset{}{\underset{N\times N}{\left(R_{4n}^\pm \right)_{\nu \mu }}\underset{1\times N}{\left\{g_{zn}^\pm \right\}_j}}=\underset{1\times N}{\left\{g_n^{e\pm }\right\}_j}+\underset{1\times N}{\left\{g_n^{a\pm }\right\}_j};$$
(60)
$$\underset{N\times N}{\left(R_{1n}^\pm \right)_{\nu \mu }}\underset{1\times N}{\left\{f_{rn}^\pm \right\}_\nu }+\underset{}{\underset{N\times N}{\left(R_{2n}^\pm \right)_{\nu \mu }}\underset{1\times N}{\left\{g_{zn}^\pm \right\}_j}}=\underset{N\times N}{\left(E_{1n}^\pm \right)_{\nu \mu }}\underset{1\times N}{\left\{f_n^{o\pm }\right\}_\nu };$$
(61)
$$\underset{N\times N}{\left(R_{3n}^\pm \right)_{\nu \mu }}\underset{1\times N}{\left\{f_{rn}^\pm \right\}_\nu }+\underset{}{\underset{N\times N}{\left(R_{4n}^\pm \right)_{\nu \mu }}\underset{1\times N}{\left\{g_{zn}^\pm \right\}_j}}=\underset{N\times N}{\left(E_{2n}^\pm \right)_{\nu \mu }}\underset{1\times N}{\left\{f_n^{o\pm }\right\}_\nu }$$
(62)
$$\underset{N\times N}{\left(E_{3n}^\pm \right)_{\nu \mu }}\underset{1\times N}{\left\{f_n^{o\pm }\right\}_\nu }+\underset{N\times N}{\left(M_{1n}^\pm \right)_{j\mu }}\underset{1\times N}{\left\{g_n^{o\pm }\right\}_j}=\underset{1\times N}{\left\{f_n^{o\pm }\right\}_\nu };$$
(63)
$$\underset{N\times N}{\left(E_{4n}^\pm \right)_{\nu \mu }}\underset{1\times N}{\left\{f_n^{o\pm }\right\}_\nu }+\underset{N\times N}{\left(M_{2n}^\pm \right)_{j\mu }}\underset{1\times N}{\left\{g_n^{o\pm }\right\}_j}=\underset{1\times N}{\left\{g_n^{o\pm }\right\}_\nu }$$
(64)
where

$$\overline{R}_{in\nu \mu }=\rho \omega ^2c^2R_{in\nu \mu };i=1,2,3,4.$$

$$f_{rn\nu }^\pm =f_{n\nu }^{e\pm }+f_{n\nu }^{a\pm }f_{n\nu }^{o\pm };g_{znj}^\pm =g_{nj}^{e\pm }+g_{nj}^{a\pm }g_{nj}^{o\pm };$$

$$E_{4n\mu \nu }^{}=\left[\underset{\nu =1}{\overset{N}{}}C_{n\nu \mu }^{}J_n\left(\gamma _{jn}\right)\left(1\right)^\nu \beta _\nu +\frac{2N}{l}\right]\delta _{\nu \mu };$$

$$E_{3n\mu \nu }^\pm =J_n\left(\varkappa _2\right)\delta _{\nu \mu };E_{2n\mu \nu }^\pm =\frac{\beta _0}{\xi _{02}\epsilon _{n\mu }}\delta _{\nu \mu };E_{1n\mu \nu }^\pm =\frac{\varkappa _2J_n^{}\left(\varkappa _2\right)}{\xi _{01}\beta _0^1};$$

$$E_{4n\mu \nu }^+=\left\{\left[\underset{\nu =1}{\overset{N}{}}C_{n\nu \mu }^+\epsilon _{\mu n}\left(1\right)^\nu \right]+\frac{2}{\pi l}\underset{\nu =N+1}{\overset{\mathrm{}}{}}\frac{1}{\nu ^2}\right\}\delta _{\nu \mu }$$

$$M_{1n\mu \nu }^\pm =C_{n\nu \mu }^\pm J_n\left(\gamma _{jn}\right);M_{2n\mu \nu }^+=\frac{\left(1\right)^\nu \varkappa _2J_n^{}\left(\varkappa _2\right)J_n\left(\gamma _{jn}\right)\epsilon _{\mu n}}{\gamma _{\mu n}^2\varkappa _2^2};$$

$$M_{2n\mu \nu }^+=\frac{\left(1\right)^\nu \varkappa _2J_n^{}\left(\varkappa _2\right)J_n\left(\gamma _{jn}\right)\epsilon _{\mu n}\beta _\nu }{\gamma _{\mu n}^2\varkappa _2^2}.$$

Thus we obtain the closed finite system of $`8\left(2N+1\right)`$ linear algebraic equations for the symmetrical and anti-symmetrical components of the solution of given problem. From the system (57) -(64) we define the unknown coefficients

$$f_{n\nu }^{e\pm },f_{n\nu }^{e\pm },g_{nj}^{e\pm },g_{nj}^{e\pm },f_{n\nu }^{o\pm },f_{n\nu }^{o\pm },g_{nj}^{o\pm },g_{nj}^{o\pm };\nu ,j=0,1,2\mathrm{}.N$$

Let us consider some particular cases, which are directly obtained from (57):

1) If liquid filler is absent inside the cylindrical cover than we obtain the soluted system of first four matrix equations (57) -(60) with the following modifications

$$f_{n\nu }^{o\pm }=0;g_{nj}^{o\pm }=0;$$

2) For absolute solid (soft) cylinder the soluted system will be formed of matrix equations (57), (58) and such two equations

$$\underset{N\times 1}{\left\{f_n^{e\pm }\right\}_\nu =}\underset{N\times 1}{\left\{f_n^{a\pm }\right\}_\nu ;}\underset{N\times 1}{\left\{g_n^{e\pm }\right\}_j=}\underset{N\times 1}{\left\{g_n^{a\pm }\right\}_j;}\underset{N\times 1}{\left\{f_n^{e\pm }\right\}_\nu =}\underset{N\times 1}{\left\{f_n^{a\pm }\right\}_\nu ;}\underset{N\times 1}{\left\{g_n^{e\pm }\right\}_j=}\underset{N\times 1}{\left\{g_n^{a\pm }\right\}_j;}$$
(65)
3) If the cylindrical cover is replaced with liquid cylinder we obtain the soluted system containing (57), (58), (63), (64) matrix equations and such matrix equations:

$$\underset{N\times 1}{\left\{f_n^{e\pm }\right\}_\mu +}\underset{N\times 1}{\left\{f_n^{a\pm }\right\}_\mu =}\underset{}{\underset{N\times N}{\left(E_{1N}^\pm \right)_{\nu \mu }}\underset{N\times 1}{\left\{f_n^{o\pm }\right\}}_{\mu \nu };}\underset{N\times 1}{\left\{g_n^{e\pm }\right\}_\mu +}\underset{N\times 1}{\left\{g_n^{a\pm }\right\}_\mu =}\underset{}{\underset{N\times N}{\left(E_{2N}^\pm \right)_{j\mu }}\underset{N\times 1}{\left\{g_n^{o\pm }\right\}_\mu };}$$
(66)
$$\underset{N\times 1}{\left\{f_n^{e\pm }\right\}_\mu +}\underset{N\times 1}{\left\{f_n^{a\pm }\right\}_\mu }\underset{}{=}\underset{N\times 1}{\left\{f_n^{o\pm }\right\}_\mu ;}\underset{N\times 1}{\left\{g_n^{e\pm }\right\}_\mu +}\underset{N\times 1}{\left\{g_n^{a\pm }\right\}_\mu }\underset{}{=}\underset{N\times 1}{\left\{g_n^{o\pm }\right\}_\mu ;}$$

4) If point pressure source

$$P^o=P_0\delta \left(rr_0\right)\delta \left(\theta \right)/4\pi r;$$

is placed inside the cylindrical cover than the system of matrix equations (57) - (64) must be changed as follows

$$f_{n\nu }^{a+}=f_{n\nu }^a=g_{nj}^a=g_{nj}^{a+}=0;$$
(67)
$$g_{nj}^{a+}=\underset{\nu =0}{\overset{\mathrm{}}{}}\frac{P_0\epsilon _\nu \epsilon _{nj}\left(1\right)^\nu J_n\left(\gamma _{jn}\right)\mathrm{cos}\beta _\nu z_0}{\pi l};g_{nj}^a=\underset{\nu =0}{\overset{\mathrm{}}{}}\frac{P_0\beta _\nu \epsilon _{nj}\left(1\right)^\nu J_n\left(\gamma _{jn}\right)\mathrm{sin}\beta _\nu z_0}{\pi l};$$

$$f_{nj}^a=\underset{j=0}{\overset{\mathrm{}}{}}\frac{P_0\epsilon _{nj}J_n^2\left(\gamma _{jn}\right)\mathrm{sin}\beta _\nu z_0}{\pi l};f_{nj}^{a+}=\underset{\nu =0}{\overset{\mathrm{}}{}}\frac{P_0\epsilon _\nu \epsilon _{nj}J_n\left(\gamma _{jn}\right)\mathrm{cos}\beta _\nu z_0}{\pi l}.$$

Thus, in this work a new problem solution methodology of contact interaction of acoustic medium with resilient finite bodies of cylindrical form, based on application of boundary integral equations method in conjunction with series method with later use of series convergence improvement taking into account the particularities of the determined functions.