key: cord-0906766-yl8zq78f
authors: Panchore, Vijay
title: Selection of Radial Basis Functions for the Accuracy of Meshfree Galerkin Method in Rotating Euler–Bernoulli Beam Problem
date: 2022-04-30
journal: Int J Appl Comput Math
DOI: 10.1007/s40819-022-01327-z
sha: 34f7003cadf0c4701587970399078cca11449ce8
doc_id: 906766
cord_uid: yl8zq78f

In this work, the radial basis function approximations are used to improve the accuracy of meshfree Galerkin method. The method is applied to the free vibration problems of non-rotating and rotating Euler–Bernoulli beams. The stiffness and mass matrices are derived by using conventional methods. In this meshfree method, only six nodes are considered within a single sub-domain. The parameters are varied for different approximations; the results are obtained with different approximations and found accurate. Two new basis function have been developed which are relatively accurate than conventional basis function: the first new basis function is obtained by multiplication of linear function to radial basis function and second new basis function is obtained by multiplying cubuic radial basis function to Gaussian radial basis function. The first few modes show same result that is available in literature using finite element method and higher modes are found very accurate as well. The result are found to be more accurate for first three modes of non-rotating and rotating Euler–Bernoulli beams where the cantilever beam boundary conditions are used; the first three modes do not change with the change in the parameter c of radial basis function.

Radial basis functions are frequently used in numerical methods; the meshfree method is used for structural problems where it is superior to conventional finite element method. Generally, moving least squares functions are used in formulation of meshfree method but it does not satisfy the Kronecker delta property like radial basis functions; it is easier to apply the essential boundary conditions with radial basis functions. The selection of parameters for radial basis functions plays an important role for the accuracy of method. Also, it depends on the problem which we solve. Number of solutions can be obtained with changing the parameters of radial basis function: in meshfree methods the length of the subdomain may change and accurate results can be obtained with only few nodes. The meshfree Galerkin method can be used to conclude the results using radial basis function.

The traditional Galerkin method and Galerkin finite element method are generally used. The governing partial differential equation of rotating Euler-Bernoulli beam has been solved with traditional Galerkin method, Galerkin finite element method and Galerkin meshfree method. In Galerkin meshfree method, the subdomain for trial and test functions remains same; we have to be careful with selection of basis functions. The Petrov-Galerkin method which is a truly meshless method has explored in literature as well with radial basis functions. It is very essential to select the parameters of radial basis functions based on the problem. In literature, it is also found that summation of two radial basis functions increases the subdomain length: increase in number of nodes within the subdomain. The stiffness matrix of a rotating Euler-Bernoulli beam has been converted in symmetric form in recent literature.

The radial basis functions are generally used in neural networks and other computations; the radial basis functions are used to solve the differential equations [1] . The collective approach of radial basis functions and finite element method has been used to solve the electromagnetic problems [2] . The radial basis functions which are self-organizing are used for the development of neural network [3] . These entire radial basis functions are assumed in its original form while the improvements can be done based on the nature of the problem. In numerical problems, radial basis functions are assumed as basis functions and results obtained are accurate; the higher order derivatives of a radial basis function can be computed easily. The radial basis function has been frequently used with meshfree methods as well where the shape functions are defined on the nodes.

The meshless methods have relative advantage in few areas to finite element method: In problems of structures having cracks, the nodes can be placed at the location of crack and the analysis can be done. The meshfree methods have been used with different basis function: (1) moving least squares basis function, (2) basis functions derived from Kriging polynomial and (3) radial basis function. The advantage of radial basis function is that it does satisfy the Kronecker delta property and the essential boundary conditions can be applied easily. The formulation involves mostly Galerkin method or Petrov-Galerkin method. To check the initial results of the method Galerkin method can be used where the subdomain of trial and test function is same.

The basic accuracy has been discussed with the combination of meshless method and radial basis function [4] . The Dirichlet problems have been solved using this combination [5] . The Navier-Stokes equation has also been solved where results show good agreement with literature [6] . Partial differential equations have been solved with this combination [7] . The collocation method has been used to solve the partial differential equation with similar combination of meshless method and radial basis function [8] . The radial basis function interpolation along with the barycentric rational interpolation has been used to 2D viscoelastic wave equation [9] . The global and local meshless method has been used for solving problems of incompressible fluid flow with heat transfer where radial basis functions have been used [10] . Kernel Techniques have also been used for meshless methods [11] .

The ordinary differential equations followed by high order L-stable method have been discussed in literature [12] . The partial differential equations have been solved using the integrated multi quadric radial basis function approximation methods [13] . The radial basis functions with finite difference method have been discussed in literature [14] . The unsteady Burgers' equation has been solved using the meshless methods where radial basis functions are used as basis functions [15] . On arbitrary surfaces, the partial differential equations are solved by orthogonal gradients method where basis functions are radial basis functions [16] . A unique stable basis for the interpolation using radial basis function has been discussed as well [17] . The advection-diffusion-reaction equations have been solved using the local radial basis function method [18] . Non local boundary value problems with Neumann's boundary conditions are solved using the radial basis function method [19] . In literature, the orthogonal polynomials have also been used as basis functions in meshless method approximation.

An analytical solution has been obtained with the combination of orthogonal polynomials and Walsh function for fractional transport equation [20] . The general fractional model of COVID-19 has been comparatively studied with the effects of isolation and quarantine where general fractional model provides better results when compared to conventional methods [21] . A novel fractional chaotic system has been studied where quadratic and cubic nonlinearities are considered and a nonstandard finite difference scheme is used [22] . The motion of a beam on nanowire has been studied where Lagrangian and fractional Euler-Lagrangian equations [23] . A similar formulation has been used for capacitor microphone case study [24] . A fractional hybrid integro-differential equation with mixed hybrid integral boundary value conditioned was investigated which leads to less material consumption [25] . A new fractional hybrid value problem has been explored with examples [26] . The solution for fractional differential inclusions has been studied which is useful in solving real world problems [27] . The extensions of fractional thermostat model were provided using the solved examples [28] .

The Caputo-Fabrizio fractional modeling has been used to study the hearing loss due to Mumps virus [29] . It has also been used to analyze transient response of the parallel RCL circuit [30] . The generalized Caputo fractional derivatives were used to investigate the p-Laplacian nonperiodic nonlinear boundary value problem [31] . In the recent literature, the extended fractional Caputo-Fabrizio derivation has been used for fractional integrodifferential inclusions [32] . This derivative has also been applied to higher order fractional integro-differential equations and to analyze the model of HIV-1 infection of CD4 + T-cell [33, 34] . The new class of hybrid type fractional differential equations has been solved where the examples are provided as well [26] .

In this paper, we solve the free vibration problem of rotating Euler-Bernoulli beam using the meshless method where radial basis functions are used. Generally, the problem is solved using the finite element method [35] [36] [37] [38] . The non-rotating beam solutions are also obtained by using the meshless method where radial basis functions are used [39] . The parameter used in radial basis function approximation is never explored. In this study, we change the parameter of radial basis function to get accurate solutions of rotating beams. Also, based on the concept of centrifugal force which varies along the length of the beam, we develop a new basis function which is found to be more accurate for the first three modes. The results are discussed in detail where only one sub-domain is considered for trial and test functions.

The accurate values of this problem have been obtained in recent literature [40] . The results include the non-rotating and rotating Euler-Bernoulli beam and non-uniform rotating and non-uniform non-rotating beam as well. The solutions can be obtained with numbers of techniques where the radial basis functions are used for weak form of Petrov-Galerkin method. Here we choose this method and notice the change in parameters which influence the results.

The partial differential equation for a rotating Euler-Bernoulli beam is given by

For a fixed free beam the boundary conditions and initial conditions are given by

Here, the flexural stiffness is E I (x), mass per unit length is m(x). G(x) is the centrifugal force which is given by

where, is the angular velocity and R is the radius of the rotating beam. The centrifugal force for a rotating beam is maximum at the root while minimum at the free end. Figure 1 shows a rotating Euler-Bernoulli beam.

The problem will be an ordinary differential equation after substituting w(x, t) e iωt w(x) to get

The flexural stiffness and mass per unit length will be constant for a uniform beam and the equation can be written as

where, ω is the natural frequency. The weak form of Eq. (4) is given by

where, v(x) is the test function.

Writing

The terms outside the integral will be zero as we satisfy the natural boundary conditions. Then we get

The mass and stiffness matrices are given by

The test and trial functions are similar in Galerkin method and are having the same subdomain length (

s ). Figure 2 shows the distribution of the nodes in a rotating Euler-Bernoulli beam. In the current study only six nodes are considered for the subdomain of trial and test function.

The transverse displacement of an Euler-Bernoulli beam is given by

where, a 1 , b 1 , a 2 , b 2 ,…, a N , b N are arbitrary constants. The radial basis function is given by

Generally, the value of c is considered as one. In the present study we explore the values for c 1, 1 2 , 1 5 , 1 10 , 1 50 , 1 100

The derivative of radial basis function is given by

Values of c and s j are user defined. The slope is given by

We can then rewrite the transverse displacement as

where,

and

The slope can be given as

where,

Eventually the assembled equation can be written as where,

and

Here, w 1 , θ 1 , w 2 , θ 2 ,…, w N , θ N are the nodal degrees of freedom. Rewriting Eq. 21, we get

The transverse displacement can be given as

where, [H (x)] is the shape function vector.

i (x) are the shape functions associated with node i. We can then write transverse displacement as

In meshfree Galerkin method, we assume the test function similar to the trial function and it is given by The mass matrix and stiffness matrix for each subdomain is given by

In case of multiple subdomain the assembly is required. Figure 3 shows the multiple subdomains.

In case of rotating Euler-Bernoulli beam, the centrifugal force varies along the length. To capture this, the radial basis function is multiplied with function x. The resulting function is given by

The results show that the first, second, and third mode can be relatively accurate when we use above basis functions. The rest of the procedure is similar to radial basis function approximation.

The Gaussian radial basis function when multiplied to cubic radial basis function the approximation is give by Tables 3 and 4 . Figure 4 shows the respective mode shapes.

The results of new basis function Q j (x) are given in Tables 5 and 6 . Figure 5 show the respective mode shapes. Here in Table 5 , the results are not accurate for c 1 and c 0.5 which clearly show the importance of approximation. The similar observation can be seen in Table 6 .

The parameter (c) has been varied to get the approximate solutions; we observed that for better accuracy the values of c are in between 1 and 0.01. The results get accuracy for the decreasing value of c. The results are obtained for non-rotating and rotating Euler-Bernoulli beams. The two new basis functions which have been developed provide better results for the first three modes when compared to conventional method. The stiffness and mass matrices have been derived for meshfree method. The results are found for a fixed-free beam, which are accurate while considering only six nodes within the sub-domain of trial and test function. The new basis function is introduced to see the effect of centrifugal force in the solution. The results clearly show the importance of approximation in numerical methods: finite element method and meshfree methods. These methods can be applied to a numbers of problems such as Rayleigh beam theory and Timoshenko beam theory where the results may or may not vary with change in radial basis function paramter c. The results will also vary with numebrs on nodes within the sub-domain of trial function. One more important observation is for second new basis function the results are not accurate for c 1 and c 0.5. The radial basis function is having numbers of approximations: Gaussian functions, multi-quadric functions and thin plate spline functon. These approximations provides results which are not similar as we change the parameters.

Funding No funding is involved.

The mathematical codes will not be available while other datas can be shared by author.

Conflict of interest On behalf of all authors, the corresponding author states that there is no conflict of interest.

Multilayer perceptrons and radial basis function neural network methods for the solution of differential equations: a survey

Hybrid approach of radial basis function and finite element method for electromagnetic problems

Research on an online self-organizing radial basis function neural network

A note on the meshless method using radial basis functions

A meshless Galerkin method for Dirichlet problems using radial basis functions

Radial basis function meshless method for the steady incompressible Navier-Stokes equations

Meshless Galerkin methods using radial basis functions

Mehsless methods based on collocation with radial basis functions

Two meshless methods based on local radial basis function and barycentric rational interpolation for solving 2D viscoelastic wave equation

Global versus localized RBF meshless methods for solving incompressible fluid flow with heat transfer

Kernel techniques: from machine learning to meshless methods

Application of radial basis function with L-stable Padé time marching scheme for pricing exotic option

Integrated multiquadric radial basis function approximation methods

Stable calculation of Gaussian-based RBF-FD stencils

A meshless RBF method for computing a numerical solution of unsteady Burgers'-type equations

The orthogonal gradients method: a radial basis functions method for solving partial differential equations on arbitrary surfaces

A new stable basis for radial basis function interpolation

A local radial basis function method for advection-diffusion-reaction equations on complexly shaped domains

Radial basis functions method for solving of a non-local boundary value problem with Neumann's boundary conditions

Analytical method based on Walsh function combined with orthogonal polynomial for fractional transport equation

A new comparative study on the general fractional model of COVID-19 with isolation and quarantine effects

A nonstandard finite difference scheme for the modelling and nonidentical synchronization of a novel fractional chaotic system

Novel fractional-order Lagrangian to describe motion of beam on nanowire

A new and general fractional Lagrangian approach: a capacitor microphone case study

On a fractional hybrid integro-differential equation with mixed hybrid integral boundary value conditions by using three operators

On the new fractional hybrid boundary value problems with three-point integral hybrid conditions

On two fraction differential inclusions

A hybrid Caputo fractional modelling for thermostat with hybrid boundary value conditions

A theoretical study of the Caputo-Fabrizio fractional modelling for hearing loss due to Mumps virus with optimal control

Analyzing transient response of the parallel RCL circuit by using the Caputo-Fabrizio fractional derivative

Investigation of the p-Laplacian nonperiodic nonlinear boundary value problem via generalized Caputo fractional derivatives

On fractional integro-differential inclusions via the extended fractional Caputo-Fabrizio derivation

On higher order integro-differential equations including the Caputo-Fabrizio derivative

Analysis of the model of HIV-1 infection of CD4+ T-Cell with a new approach of fractional derivative

Free-vibration analysis of rotating beams by a variable-order finite element method

Rotor blade vibration by the Galerkin finite element method

Finite element approach to rotor blade modeling

Vibration of a rotating beam with tip mass

A radial basis function approach in the meshless local Petrov-Galerkin method for Euler-Bernoulli beam problems

Free vibration analysis of rotating blades with uniform tapers

Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations