On the Application of Modified Finite Sine Transform to Structural Mechanics

Themethod of modified finite sine transform (MFST) was introduced to solve fourth-order boundary value problems in structural mechanics. The analytical features and computational aspects of MFST were presented, including integral transforms of the derivatives of functions and the term-by-term differentiation of inverse MFST. For the simplicity of demonstrating the solution process, the cantilever beamwas used in the discussion.Numerical results for transverse vibration based on the presentmethodwere compared and validated with existing results for natural frequencies. Although the present method is not superior to conventional finite sine transform method in one-dimensional structures, it could be significantly efficient for two-dimensional structures with asymmetric boundary conditions, especially including cantilever plates and other plates with only one free edge or two adjacent free edges.


Introduction
The modified Fourier series has recently been introduced in [1] as an effective supplement of classical Fourier series for the approximation of nonperiodic functions, whereby in [−1, 1] the sin  functions are adjusted by sin( − 1/2),  ⩾ 1.The main immediate advantage of the modified Fourier series is the faster convergence rate [2], which can also be further accelerated [3].On the other hand, dating back to the late 1940s, the modified finite Fourier sine transform (MFST) was defined by Roettinger [4] to solve the second-order differential equations with mixed/asymmetric boundary conditions (e.g., (0) = 0 and   () = 0 in [0, ]), where the integral kernel is of sin( − 1/2),  = 1, 2, . . . in [0, ].Such transformation complements the finite Fourier sine transform (FST) which involves the Dirichlet boundary conditions, and the finite Fourier cosine transform (FCT) which is compatible with Neumann boundary conditions of second-order differential equations [5].
In the field of structural mechanics, Fourier sine/cosine series equipped with Stokes's transformation have been extensively applied to higher-order boundary value problems for beams [6][7][8][9][10], plates [11], and shells [12,13].Similarly, FST and FCT have recently been adopted to obtain the exact series solutions of bending of plates with complex boundary conditions [14][15][16][17][18].However, although MFST was introduced at the same time as FST and FCT [4,19], very little attention has been paid to the application of MFST to higherorder boundary value problems in structural mechanics.Furthermore, it can be found that the applications of FST and FCT to bending and vibration of plates are mostly reported for symmetry boundary conditions.Plates with asymmetric boundary conditions such as cantilever plates and other plates with only one free edge or two adjacent free edges are rarely addressed by using FST and FCT.One possible reason is that these two methods are not applicable.Instead, MFST can provide efficient solutions.
In the present study, MFST is introduced to solve the fourth-order boundary value problems with asymmetric boundary conditions.At first, MFST was defined and MFST of the function's derivatives were derived.Then, the Stokes's transformation was applied to obtain the term-by-term differentiations of inverse MFST (i.e., modified sine series).A simple example of free vibration of cantilever beam was analyzed by using MFST to present the solving procedure.At last, more detailed discussions involving MFST and conventional Fourier series as well as FST and FCT are presented.It can be concluded that MFST presents genuine advantages over classical FST and FCT in specific situations of asymmetric boundary conditions.

Modified Finite Sine Transform
MFST defined on the interval (0, ) can be given by [4] The inverse transformation is expressed as Similarly, the pair of modified finite cosine transform (MFCT) is defined by It can be found that MFCT and MFST are essentially the same [19] because Using integration by parts, MFST of derivatives of functions can be obtained, which is presented in Appendix A.

Term-By-Term Differentiations of Inverse MFST
The inverse MFST will be a modified Fourier sine series, which can be expressed by where   = ( − 1/2)/ and the Fourier coefficient   is given by Note that the end value (0) will be forced to be zero but not for () if they are included in the modified sine series in (5).Therefore, (0) on the interval [0, ] will be defined separately.
In this manner, the function of solution () is defined in two separate regions, one is  = 0 and the other is 0 <  ⩽ .Thus, Stokes's transformation [20, pp. 375] should be employed to differentiate these series properly.Care must be taken when finding the correct modified cosine series corresponding to   (), which should not include the constant term as that of Fourier cosine series.Thus, the first-order derivative   () can be formulated by in which The expressions of first four derivatives of () obtained by using Stokes' transformation can be found in Appendix B.

Example: Transverse Vibration of Cantilever Beams
In this section, MFST is applied to obtain the series solution of free vibration of a cantilever beam as shown in Figure 1.Such simple example is presented solely for illustrating the procedure and validating the method.The governing equation of free vibration of Euler beam can be expressed as [21]   4   4 + where (, ) is the displacement function,  is the flexural rigidity,  is the density, and  is the cross-sectional area of the beam.Using the method of separation of variables, the solution can be written as Substituting ( 10) into (9) yields where ,  =  4 / 4 and  4 =  2 /().The boundary conditions for the cantilever beam as shown in Figure 1  Taking MFST to both sides of (11) and substituting (A.3d) and boundary conditions result in where If ,  () and ,  (0) in ( 13) can be expressed by  *  (), (13) will become the characteristic equations for the frequency parameter , which can easily be solved as eigenvalue problems.First, ,  () can be expressed by  *  ().Based on (B.3), the second derivative of () can be obtained by Substituting the boundary conditions By truncating the infinite series to finite terms of , ,  () can be obtained by However, ,  (0) can not be directly expressed by  *  ().A slight change needs to be made on the boundary conditions.The cantilever beam in Figure 1 is modified to a rotationally restrained beam with a free edge as shown in Figure 2. Thus, the boundary condition at  = 0 will be The rotational stiffness  can be replaced by the end-fixity factor [15] as follows: For cantilever beam,  = ∞ and  = 1.By using (19), the cantilever condition can be approximated by a value of  close to 1, i.e.,  = 0.9999.Based on the Stokes' transformation for ,  () in (B.2), ,  (0) can be obtained by Substituting ( 19) and ( 20) into (18) yields At last, by substituting ,  () in ( 17) and ,  (0) in ( 21) into (13), a characteristic equation for the frequency parameter  can be obtained as Substituting   = ( − 1/2)/ into (22) results in where Equation ( 23) can be expressed in the following matrix form: where  = [ *  (1),  *  (2), . . .,  *  ()] and A is the corresponding coefficient matrix of size  ×  whose elements can be expressed as in which  = 1, 2, . . ., ,  = 1, 2, . . ., .
The natural frequencies can be obtained by solving the eigenvalue problem of (25) and the corresponding mode shape can be determined in terms of modified Fourier sine series by substituting  *  () of the eigenvector into (14b).The convergence of the first frequency parameter  is shown in Figure 3 for  = 0.9999.From the comparison with the exact value of [21], it can be found that the present result of  = 500 is close to the exact value and the convergence rate is not uniform after  = 1000 because overflow problems occurred.Thus, in present research,  = 500 is used to obtain the numerical results of  for validation.
As illustrated in Table 1, numerical results for the first six frequency parameters obtained by using MFST for rotationally restrained beams with  = 0.9999 are listed and compared with exact values of a cantilever beam in [21].Excellent agreements between the two methods can be observed in Table 1.

Discussion and Remarks
On the other hand, by using Fourier cosine series reported in [9] (i.e., inverse process of FCT), the frequency equation for cantilever beams can be expressed by where 4  (28a) The frequency parameter  can be obtained by solving the nonlinear equation of (27).By comparing with (25); one main disadvantage for solving (27) is that its procedure involves root-finding of nonlinear functions, which inevitably needs initial trial values and is often inconvenient and cumbersome.Furthermore, when extending such solution process to dynamic analysis of two-dimensional structures such as previous attempts made by [22,23] for completely free plates, the resulting nonlinear functions are highly complex and their root-finding for the frequency parameter is extremely difficult as pointed out by [17].
In addition, the solution process of MFST can be easily extended to two-dimensional structural mechanics problems such as plates or shells with one free edge, two adjacent free edges, and three free edges, because all these twodimensional structures involve so-called asymmetric boundary conditions for two opposite edges (i.e., one edge free and the opposite edge simply supported or clamped).For instance, the transverse vibration of cantilever plates can be solved by combining MFST and FCT for two-dimensional functions.For plates with two adjacent free edges, taking MFST on both directions of plates can be adopted.If plates have only one free edge, the joint transforms of MFST and FST would be used.Conversely, the conventional FST and FCT might not be efficiently applied to dynamic analysis of these plates with asymmetric boundary conditions.Therefore, although MFST has not shown significant advantages over conventional Fourier sine/cosine series, FST, or FCT in one-dimensional problems, it exhibits genuine advantages in two-dimensional structures with free edges.

Conclusions
This paper is devoted to analytical features, computational aspects, and applications of modified finite sine transform (MFST).Special care must be taken when finding derivatives of inverse MFST by using Stokes's transformation since the constant term associated with cosine series should not be included.Although only free vibration of cantilever beam was treated for demonstration, MFST can be extended for two-dimensional structures with asymmetric boundary conditions.Despite no significant advantages over conventional Fourier sine/cosine series, FST, or FCT in onedimensional problems, MFST exhibits genuine advantages in two-dimensional structures with asymmetric boundary conditions such as cantilever plates and other plates with only one free edge or two adjacent free edges.

Figure 3 :
Figure 3: Convergence of the fundamental frequency parameter  with  = 0.9999.

Table 1 :
Comparisons of frequency parameters, , between rotationally restrained beam with  0 = 0.9999 and cantilever beam.