1 st-Order Shear Deformable Beam Formulation Based on Meshless Wavelet Galerkin Method

This paper examined and discussed aMeshlessWavelet GalerkinMethod (MWGM) formulation for a first-order shear deformable beam, the properties of the MWGM, the differences between the MWGM and EFG, and programming methods for the MWGM. The first-order shear deformable beam (FSDB) consists of a pair of second-order elliptic differential equations. The weak forms of two differential equations are deduced using Hat wavelet series.The exact integration and reduced integration were used to analyze the problems. Some indeterminate beam problems are considered. Condition numbers of the stiffness matrix were analyzed with exact integration and reduced integration for two cases of these problems. Consequently, the results were converged on the analytic solutions. The shear-locking phenomenon also occurred in the MWGM as it occurs in the conventional FEM.The stiffness matrix calculated from the reduced integration causes a similar numerical error to the stiffnessmatrix calculated from the exact integration in the MWGM.The MWGM showed desirable results in the examples.


Introduction
The aim of this study was to formulate a Meshless Wavelet Galerkin Method (MWGM) to solve a pair of 2nd-order elliptic differential equations, in other words, the Timoshenko beam differential equation using Hat wavelets.Usually the formulation of a MWGM is similar to the conventional formulation of FEM.However, some differences exist and this paper is focused on these differences.Additionally, the property of the 1st-order shear deformable theory, which is the locking phenomenon, is discussed.
Meanwhile, the wavelet theory has evolved over the last few decades.Similar to the Fourier series expansion, wavelets are used to express arbitrary signals.Even though wavelets so far have been used to process signals, they can be used to solve differential equations.Moreover, wavelet numerical methods have been investigated by many researchers in numerical analysis as well as structural analysis fields.For instance, the B-spline wavelet finite element method was developed to analyze vibrations in structures [1].The wavelet element method was used to combine biorthogonal wavelet systems using the spectral element method [2].A study that used Daubechies wavelets as the shape functions was performed to solve the Euler beam bending problem [3].Additionally, a study using B-spline wavelet as the shape functions was conducted [4].Basu et al. did a study comparing the FEM, BEM, Meshless Method, and Wavelet Methods [5,6].Meanwhile, thin plate problems were analyzed with Daubechies wavelets [7].Recently Daubechies wavelets and Deslauriers-Dubuc interpolating functions were used to analyze the Euler beam on elastic foundation and thin plate problem [8].Most of the studies in the literature have been about spline wavelet and Daubechies wavelet.
This study used the Hat wavelet series for expandability and for the ability of inosculation to other numerical methods.The possibility of the inosculation of numerical methods has become important recently because existing methods have already come into wide use such as the standard h-FEM, BEM, and Meshless Method.It is better if expansion and inosculation as well as realization are easy.

Hat Wavelets Series
2.1. 2 Function and the Space   (Ω).Usually, the first derivative of the shape functions of FEM must be square integrable function in the Lévesque sense.In most FEMs, approximations of variational boundary-value problems are concerned with the square integral of the first derivative of functions.As shown in Figures 1 and 2, the wavelet series (scaling functions and wavelet functions) belong to  0 (Ω) and  2 (Ω).The first derivatives of the wavelet series are discontinuous but well defined members of  2 (Ω), in other words, those with the property for which ∫  2  < ∞ [9].
The definition of Sobolev classes may be summarized as follows.A function  belongs to the Sobolev classes   (Ω) if  and all of its partial derivatives of order less than or equal to  ( being to a nonnegative integer) are members of  2 (Ω).In (1),  usually is set to 1 because the first derivatives of functions are considered in the variational boundary-value problem.
With the above considerations, the following properties should be satisfied to use these functions as basic functions for MWGM.
and Φ should be sufficiently differentiable and have a square integrable property. and Φ should satisfy the boundary conditions of the problem. and Φ should have a completeness property.
To summarize,  and Φ should exist in the  1 - and have the same continuity requirement.

Basic Approximation Functions Based on Wavelets and
Their Construction Method.The basic scaling functions and wavelet functions should be combined linearly to construct the basis functions for the wavelet series analysis as in (2).In a one-dimensional problem domain, the unknown function  can be approximated by a series expansion.Generally, the series expansion can be defined in terms of wavelets as in (2) as follows: where () ∈  2 () is the scaling functions; Φ() is the wavelet functions; () and (, ) are the coefficients;  ∈  is the scaling parameter;  ∈  is the translation index (controls the spatial location);  0 is the starting scale parameter; and  is the space of all integers.In the expression for   0 , (), the factor 2  0 /2 helps maintain a constant norm independent of the scale  0 .The coarsest scale for the space spanned by   () depends on the parameter  0 .In the proposed research, it is assumed that  0 = 0; that is, the coarsest space is assumed for the initial resolution.This will lead to   0 , () = ( − ), and (2) becomes The subspace of  2 () spanned by the scaling functions   () = ( − ) is obtained as  0 = Span  {  ()}.The different subspaces are related as shown in The next equation is used for approximating signals or functions over the span (−∞, ∞); thus, the translation factor  is infinity covering all the integers from −∞ and ∞ shown in (6).The integer  will be finite to be used in the weak form of the MWGM because the infinity of  cannot be used in the numerical analysis.
Wavelet basis functions consist of scaling functions and wavelet functions as mentioned above.The generalized expression of extension is shown in (7) and is used as a refinement scheme.Equation (7) corresponds to the refinement of the mesh for the h-version finite element method.
The following expressions represent examples of the relationship for each level,  = 0, 1, 2, 3: The extension of the spaces shows that the spaces are hierarchical, and the refinement can be done by using a suitable value for .If  is taken to be ∞, the results of the analysis are exact solutions.

The Property of the Hat Wavelet
Series.There are many wavelet series including the Daubechies, Trigonometric, Hermite Cubic, Haar, and Hat.Because of its suitable properties, this study used the Hat wavelet series as basis functions for the numerical analysis of the weak form of the differential equation, in other words, test functions and trial functions.
Hat wavelet functions have a continuity property but not an orthonormality.After differentiation, the Hat wavelet functions are discontinuous and have a jump property.However the wavelet functions can be integrable after differentiation.The first derivative of the Hat wavelet functions is square integrable.This property is suitable for basis functions for the weak form of the 2nd-order elliptic differential equation.Equations ( 9)∼(12) represent the Hat wavelet functions.They are continuous but do not form an orthonormal basis.Figure 1 shows the relatively simple Hat wavelet functions, which are constructed with (9) and (11) and are used as the basis functions of the MWGM.Usually, ( 9) is called basic scaling functions and ( 11) is called wavelet functions.
Basic scaling functions (in  0 ) are Wavelet functions (in  0 ) are

The Weak Form of the First-Order Shear Deformable Beam Using Hat Wavelets
The first-order shear deformable beam (FSDB) consists of a pair of second-order elliptic differential equations (13).Each differential equation consists of two independent variables  and Ψ, where  is the transverse function and Ψ is the rotation function.The rotation of the -axis is denoted by an independent function Ψ.The governing equilibrium equations of the FSDB, namely, the Timoshenko beam theory, are shown by where  is the modulus of elasticity,  is the moment of inertia,  is the shear modulus, and   is the shear correction coefficient.
The weak forms of two differential equations are deduced by multiplying each differential equation with the test functions  1 and  2 [10,11].Each of the test functions  1 and  2 corresponds to  and Ψ.The test functions of the weak forms belong to  0 (Ω) and exist in  2 (Ω); consequently, the abovementioned Hat wavelet functions can be used as the test functions.Moreover, if the Hat wavelet functions can be used for the trial functions of  and Ψ, the symmetry property of the stiffness matrix still remains in the MWGM as the standard FEM.The weak form over the domain is developed using the usual procedure.Equation ( 14) is multiplied by each of the weight functions − 1 and − 2 and integrated over the domain.
It is clear that ( 14) is as follows: If the integrands are integrable on the problem domain Ω, (16) are obtained by integrating the first of each integral once by parts.The integral domain is taken over Ω in (14).Generally, after calculating the element stiffness matrices, the assembling procedure of the element stiffness matrices should be performed.However, nonzero values for the scaling functions and wavelet functions exist on several elements continuously in the MWGM.Accordingly, it is very easy to calculate the assembled global stiffness matrix if the integral domain is taken over the whole domain where the problem is defined.And in doing so, the procedures for assembling the element equations, that is, the stiffness matrices and load vectors, can be circumvented.In other words, the reason for not using the element but instead the domain to integrate the equations is that a global stiffness matrix of the problem domain can be calculated without assembling the element stiffness matrices.In fact, only integration sections are used to integrate a square of the first derivatives of the shape functions over the domain in the MWGM instead of the above-mentioned "element."The properties of the sections are discussed next.
As in (16), second terms of the right hand sides of these equations except for the weight functions  1 and  2 are the shear force and the bending moment shown in ( 17) and (18), respectively.
Equation ( 19) is the approximation functions of  and Ψ.If the scaling functions and wavelet functions can be considered as the same kind of functions,  1 and  2 can be concisely written as (20).The forms of (20) are similar to the standard FEM.The numbers  1 and  2 are the number of polynomial terms.If  1 and  2 are not equal, consistent interpolation elements appear in the equations.In this case, a reduced integration does not need to be used.If  1 and  2 are equal, reduced integration elements with equal interpolation appear in the equations.These attributes are equal to the standard FEM.
(1)    (1)   () , Equation ( 21) follows typical procedures for developing linear equations of the weak form of the partial differential equation.Some differences exist in (21) because of the shape functions that are used in the MWGM.First, (21) has [0, ] as integral domains.These are due to the property of the scaling functions and wavelet functions.This property is that nonzero values for the scaling functions and wavelet functions exist over sections that are introduced to integrate the equations.The sections should be determined to be small enough to integrate the equation on the whole domain.If the sections are determined to be large, the integration cannot be calculated because the scaling functions and wavelet functions are folded in the sections.Second, the solutions of (21) are not the values for the node but for the coefficients (or parameters) of the approximation functions.This also is due to the property of the scaling functions and wavelet functions.Bases of wavelet functions have coefficients of functions whose locations are not determined as unknowns because they have degree of freedoms of functions without fixed positions.Usually, standard FEM has the discontinuous property of the solution derivative in the boundary of the elements.MWGM has the same property also.The derivatives of solutions in the section boundary have the discontinuous property, in other words,  0 continuity, because scaling functions and wavelet functions are only piecewise basis functions for the approximation.

The Comparison of the EFG Method and MWGM
At this point, the similarities and difference of the EFG and MWGM are discussed.As for the similarities, first, the EFG like the MWGM also uses basis functions that exist over the cell which is introduced to integrate the basis functions in the problem domain.Thus, integration and assembly of the stiffness matrix for these two methods can be performed very similarly.As previously stated, sections are introduced to integrate the problem domain in the MWGM.These sections have an effect on the procedures of the MWGM as cells affect the procedures of the EFG.Second, EFG cannot obtain the solutions of the nodes directly; thus, it has to perform additional calculations with the parameters of the basis functions.MWGM also follows the same procedures using the parameters of the basis functions.
Finally, sections have also a property similar to that of the standard FEM elements, which is a discontinuous attribute in the boundary of the elements.This property stems from the folded shape of the wavelet basis functions.If the folded wavelet basis functions exist in one section, sections as many as wavelet basis functions that are folded should be inserted into a preexisting section to integrate wavelet basis functions.This is equivalent to the mesh refinement of standard FEM and to extending subspaces, which consist of piecewise polynomials.
As for the differences, the EFG performs a Gauss elimination to calculate the basis functions value of the integration points.However, the MWGM does not perform the gauss elimination to obtain the basis functions value of the integration points.The MWGM only follows the method of the standard FEM, which uses simplified procedures for constructing stiffness matrices.
A similar case of sections that are used in this paper is found in [4].Elements are defined and divided into segments that include nodes in [4].In other words, the elements are composed of many segments assigned with nodes.Thereafter, the elements are mapped into the standard element domain Ω  ∈ [0, 1].As a result, many scaling functions and wavelet functions exist in the each elements and the problem domain is the collection of those elements.However, Hat wavelet functions were generated without explicit elements in the problem domain in this paper.The problem domain is composed of the collection of Hat wavelet functions.This is the difference of the way of using basis functions and calculation.We considered the indeterminate beam problem shown in Figure 3, which is a fixed end supported beam applied by a Dirac delta function [12,13].The exact integration and reduced integration were used to analyze the problems.Analytic solutions were also obtained for each case.The heights of the beam vary between 0.2 and 1.0.The length of the beam is 2(= 2.0).Figures 4 and 5 show the basis functions of the MWGM.The shear correction factor   was set to 5/6.As shown in Figure 6, the shear-locking phenomenon definitely occurs in the exact integration when using a small number of functions in the solutions for deflection.As the analyses progress, the shear-locking phenomenon disappears in Figure 6. Figure 6 shows that the reduced integration  softens the beam stiffness in the MWGM.As the height of the beam cross section increases progressively from Figures 6-10, the shear-locking phenomenon disappears in the analysis using fewer functions.This is because the deep beam is stiff enough to counteract the shear-locking phenomenon.

Examples and Discussion
Although the reduced integration and many functions are used simultaneously, the stiffness of the analysis domain is  not overly soft.All the analyses results do not exceed the analytic solution.
As shown in Figure 11, the shear-locking phenomenon definitely occurs in a solution of rotations with exact integration when using a small number of functions.As mentioned above, the solutions of rotations have a similar pattern shown in Figures 6-10.As the analyses progress, the shear-locking  phenomenon disappears in Figure 11.show the symmetry of the deformed domain.The reduced integration softens the beam stiffness in the MWGM with rotation solutions.It can be seen that everything mentioned above about deflection solutions also apply to rotation solutions.Additionally, condition numbers of the stiffness matrix were analyzed with exact integration and reduced integration  for two cases (ℎ = 0.2 and ℎ = 0.4) in Figures 16 and  17.The condition numbers increase linearly in both cases.Eventually, the ratios of the two values in each case approach 1.0.Thus, the stiffness matrix calculated from the reduced integration causes a similar numerical error as the stiffness matrix calculated from the exact integration in the MWGM.
Governing equations are as follows: Boundary conditions are as follows: Mathematical Problems in Engineering Next, we considered the indeterminate beam problem shown in Figure 18, which is a fixed end supported beam consisting of nonuniform sections.Additionally, the Dirac delta function is applied to this beam.The domain of the problem has a jump condition.The analytic solution cannot be obtained in a single domain but in multiple domains.Particularly, the domain  ∈ [1,2], for which the beam depth is large, has a high ratio of ℎ/2.It is difficult to analyze this type of model because of the domain jump condition.
As shown in the MWGM and analytic solution, maximum deflection occurs on the left side of the domain.Figure 19 clearly shows that the maximum deflection locations correspond with each solution.This problem also shows the shear-locking phenomenon, and, as the analyses progress, the shear-locking phenomenon disappears.
Both solutions with the exact integration and reduced integration agree well with each other in the previous problem.However, both solutions do not agree with each other in this problem.More specifically, the solutions with the exact integration correspond to the analytic solutions in domain  ∈ [1,2]; however, the solutions with the reduced integration correspond to analytic solutions in the domain  ∈ [0, 1].
It could be said that the reduced integration is effective for relatively thin beams, and the exact integration is effective for relatively thick beams.The reduced integration solution converges from a large value to a small value, in other words, an analytic solution.However, the exact integration converges from a small value to a large value, in other words, an analytic solution in the rotation solutions (Figure 20).

Conclusion
This paper examined and discussed the MWGM formulation for a first-order shear deformable beam, the properties of the MWGM, the differences between the MWGM and EFG, and programming methods for the MWGM.The MWGM has also the shear-locking phenomenon.Circumventing assembling the global stiffness matrix K enables the program to be written easily.Multiresolution and localization analysis, which is similar to adaptive analysis of the finite element method, can be performed with the MWGM.It is carefully predicted that all analysis methods could be used to solve one problem such as structure, fluid, soil, or multiphysics problems in the future.The MWGM could be used as part of the analyses.It is easy to perform multiresolution analysis, in other words, adaptive analysis using the MWGM, because integration and interpolation are separated, and linear functions are used as the bases.This study used the Hat wavelet series for expandability as well as for the ability of inosculation to other numerical methods.As mentioned in the introduction section of this paper, the MWGM is a good alternative and complementary method with traditional and recently developed methods.

Figure 3 :Figure 4 :Figure 5 :
Figure 3: The fixed end supported beam applied by a Dirac delta function.

1 Figure 18 :Figure 19 :
Figure 18: The fixed end supported beam with a nonuniform section applied by a Dirac delta function.