MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi 10.1155/2017/4868050 4868050 Research Article 1st-Order Shear Deformable Beam Formulation Based on Meshless Wavelet Galerkin Method http://orcid.org/0000-0001-7356-2475 Jo JunHyung 1 http://orcid.org/0000-0002-9855-679X Lee Yun 2 Markovic Damijan 1 Korea Electric Power Corporation Daegu Kyungbuk Construction Branch Gong Pyeung Road 96 Junggu Daegu Republic of Korea kepco.co.kr 2 Department of Civil Engineering Daejeon University 96-3 Yongun-dong Dong-gu Daejeon 300-716 Republic of Korea dju.ac.kr 2017 2782017 2017 23 05 2017 02 07 2017 2782017 2017 Copyright © 2017 JunHyung Jo and Yun Lee. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This paper examined and discussed a Meshless Wavelet Galerkin Method (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 stiffness matrix calculated from the exact integration in the MWGM. The MWGM showed desirable results in the examples.

Ministry of Education, Science and Technology NRF-2016R1A2B4016632
1. 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 . The wavelet element method was used to combine biorthogonal wavelet systems using the spectral element method . A study that used Daubechies wavelets as the shape functions was performed to solve the Euler beam bending problem . Additionally, a study using B-spline wavelet as the shape functions was conducted . 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 . Recently Daubechies wavelets and Deslauriers-Dubuc interpolating functions were used to analyze the Euler beam on elastic foundation and thin plate problem . 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.

2. Hat Wavelets Series 2.1. <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M1"><mml:mrow><mml:msup><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> Function and the Space <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M2"><mml:msup><mml:mrow><mml:mi>H</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi></mml:mrow></mml:msup><mml:mfenced separators="|"><mml:mrow><mml:mi>Ω</mml:mi></mml:mrow></mml:mfenced></mml:math></inline-formula>

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 C0(Ω) and L2(Ω). The first derivatives of the wavelet series are discontinuous but well defined members of L2(Ω), in other words, those with the property for which f2dx< .

A scaling function (φ); wavelet functions (Φ,V1=V0+W0).

Derivatives of a scaling function and wavelet functions.

The definition of Sobolev classes may be summarized as follows. A function f belongs to the Sobolev classes Hm(Ω) if f and all of its partial derivatives of order less than or equal to m (m being to a nonnegative integer) are members of L2(Ω). In (1), m usually is set to 1 because the first derivatives of functions are considered in the variational boundary-value problem.(1)HmΩ=ff,fx,fy,,mfzmL2Ω.

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 H1-space and have the same continuity requirement.

2.2. 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 u can be approximated by a series expansion. Generally, the series expansion can be defined in terms of wavelets as in (2) as follows:(2)ux=k=-cj0kφj0,kx+j=j0k=-dj,kΦj,kx(3)Withφj0,kx=2j0/2φ2j0/2x-k,Φj,kx=2j0/2Φ2jx-k,where φ(x)L2R is the scaling functions; Φ(x) is the wavelet functions; c(k) and d(j,k) are the coefficients; jZ is the scaling parameter; kZ is the translation index (controls the spatial location); j0 is the starting scale parameter; and Z is the space of all integers.

In the expression for φj0,kx, the factor 2j0/2 helps maintain a constant norm independent of the scale j0. The coarsest scale for the space spanned by φkx depends on the parameter j0. In the proposed research, it is assumed that j0=0; that is, the coarsest space is assumed for the initial resolution. This will lead to φj0,kx=φx-k, and (2) becomes (4)ux=k=-ckφkx-k+j=0k=-dj,kΦj,kx

The subspace of L2R spanned by the scaling functions φkx=φx-k is obtained as V0=Spankφk(x)¯. The different subspaces are related as shown in(5)V0V1V2V3L2Rwith  uVj.

The next equation is used for approximating signals or functions over the span (-,); thus, the translation factor k is infinity covering all the integers from - and shown in (6). The integer k will be finite to be used in the weak form of the MWGM because the infinity of k cannot be used in the numerical analysis.(6)ux=k=-ckφkx+j=0k=-dj,kΦj,kx.

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.(7)Vj+1=V1+V2+V3++Wj,j=0,.

The following expressions represent examples of the relationship for each level, j=0,1,2,3:(8)Level  0:  V13=V11+W12j=0;Level  1:  V27=V01+W02+W14=V13+W14j=1;Level  2:  V315=V01+W02+W14+W28=V27+W28j=2;Level  3:  V431=V01+W02+W14+W28+W316=V315+W316j=3.

The extension of the spaces shows that the spaces are hierarchical, and the refinement can be done by using a suitable value for j. If j is taken to be , the results of the analysis are exact solutions.

2.3. 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 V0) are(9)φ0x=x,x0,12-x,x1,20,otherwise.

Wavelet functions (in W0) are(10)Φ0x=-x,x0,0.53x-2,x0.5,1-3x+4,x1,1.5x-2,x1.5,20,otherwise,(11)Φ1x=x,x0,0.5-2x+32,x0.5,1x-32,x1,1.50,otherwise,(12)Φ2x=x,x0,0.51-x,x0.5,10,otherwise.

3. 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 w and Ψ, where w is the transverse function and Ψ is the rotation function. The rotation of the y-axis is denoted by an independent function Ψ. The governing equilibrium equations of the FSDB, namely, the Timoshenko beam theory, are shown by (13)ddxGAKsΨ+dwdx+f=0,Ω0,L,ddxEIdΨdx+GAKsΨ+dwdx=0,Ω0,L,where E is the modulus of elasticity, I is the moment of inertia, G is the shear modulus, and Ks is the shear correction coefficient.

The weak forms of two differential equations are deduced by multiplying each differential equation with the test functions w1 and w2 [10, 11]. Each of the test functions w1 and w2 corresponds to w and Ψ. The test functions of the weak forms belong to C0(Ω) and exist in L2(Ω); consequently, the above-mentioned Hat wavelet functions can be used as the test functions. Moreover, if the Hat wavelet functions can be used for the trial functions of w 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 -w1 and -w2 and integrated over the domain.(14)0=0Ω-w1GAKsΨ+dwdx+fdx,0=0Ω-w2ddxEIdΨdx+GAKsΨ+dwdxdx.

It is clear that (14) is as follows:(15)find  wH1  such  that  Bw1,v1=lv1,  v1H1find  sH1  such  that  Bs1,v2=lv2,  v2H1.

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.(16)0=0Ωdw1dxGAKsΨ+dwdx+w1fdx-w1GAKsΨ+dwdxΩ,0=0Ωdw2dxEIdΨdx+w2GAKsΨ+dwdxdx-w2EIdΨdxΩ.

As in (16), second terms of the right hand sides of these equations except for the weight functions w1 and w2 are the shear force and the bending moment shown in (17) and (18), respectively.(17)GAKsΨ+dwdxV,(18)EIdΨdxM.

Equation (19) is the approximation functions of w and Ψ. If the scaling functions and wavelet functions can be considered as the same kind of functions, w1 and w2 can be concisely written as (20). The forms of (20) are similar to the standard FEM. The numbers n1 and n2 are the number of polynomial terms. If n1 and n2 are not equal, consistent interpolation elements appear in the equations. In this case, a reduced integration does not need to be used. If n1 and n2 are equal, reduced integration elements with equal interpolation appear in the equations. These attributes are equal to the standard FEM.(19)w1x=k=-ckφkx+j=0k=-dj,kΦj,kx,w2x=k=-ckφkx+j=0k=-dj,kΦj,kx,(20)w1=i=1n1wi1χi1x,w2=i=1n2wi2χi2x.

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,L] 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, C0 continuity, because scaling functions and wavelet functions are only piecewise basis functions for the approximation.(21)0ΩGAKsdχi1dxdχj1dxwi1wj1+GAKsdχi1dxχi2wi1wj2dx=0Ωχi1fwi1dx+w10Q10+w1LQ3L,0ΩEIdχi2dxdχj2dxwi2wj2+GAKsχi2χj2wi2wj2+GAKsχj2dχi1dxwi2wj1dx=w10Q20+w1LQ4L.

Equation (22) is the equation for the matrix form. Finally, (23)–(25) are obtained, and (26) is the simplified matrix form.(22)0Ωwi1GAKsdχi1dxdχj1dxGAKsdχi1dxχi2dxwj1wj2=wi10Ωχi1fdx+wi1χi10Q10+wi1χi1LQ3L,0Ωwj2GAKsχj2dχi1dxwi2wj1EIdχi2dxdχj2dxwi2wj2+GAKsχi2χj2wi2wj2dxwj1wj2=wi2χi20Q20+wi2χi2LQ4L,(23)Kij11=0ΩGAKsdχi1dxdχj1dxdx,(24)Kij12=0ΩGAKsdχi1dxχj2dx,(25)Kij22=0ΩEIdχi2dxdχj2dx+GAKsχi2χj2dx,(26)K11K12K21K22ws=F1F2.

Integration should be performed with (27) over the whole domain to circumvent the assembly of the section stiffness matrix.(27)0l1Γ1xdx+l1lΓ1xdx++l1LΓ1xdxwj1wj2=0l1Γ2xdx+l1lΓ2xdx++l1LΓ2xdx+wi1χi10Q10+wi1χi1l1Q10+wi1χi1lQ10++wi1χi1LQ3L,(28)Γ1x=GAKsdχi1dxdχj1dxGAKsdχi1dxχi2,(29)Γ2x=χi1f.

Some sets of functions are introduced to code the above equations. Λ1Ω and Λ2Ω are sets of shape functions. If the shape functions and their derivatives are stored like in (31) and (32), the component matrices of the global matrix can be calculated by (33) and (34) without having to use complicated procedures. The symbol follows the dyadic product rule.(30)i,j=The  numer  of  wavelet  basis  functions(31)Λ1Ωdχ1dxdχ2dxdχndxχn+1χn+2χ2n(32)Λ2Ω01020nχn+1χn+2χ2n(33)K1Ω=0ΩΛ1ΩΛ1Ωdx(34)K2Ω=0ΩΛ2ΩΛ2Ωdx(35)KΩ=KΩ1+KΩ2.

4. 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 . Elements are defined and divided into segments that include nodes in . In other words, the elements are composed of many segments assigned with nodes. Thereafter, the elements are mapped into the standard element domain Ωs0,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.

5. Examples and Discussion 5.1. The Fixed End Supported Beam Applied by a Dirac Delta Function

Problem domain is (36)Ω0,2L,L=1.

Governing Equations are(37)ddxGAKsΨ+dwdx+f=δx-L(38)ddxEIdΨdx+GAKsΨ+dwdx=0.

Boundary Conditions are(39)w0=0,dw0dx=0(40)w2L=0,dw2Ldx=0.

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 2L(=2.0). Figures 4 and 5 show the basis functions of the MWGM. The shear correction factor ks was set to 5/6.

The fixed end supported beam applied by a Dirac delta function.

7 functions and 8 sections (φx: scaling functions, Φx: wavelet functions).

15 functions and 16 sections (φx: scaling functions, Φx: wavelet functions).

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 610, 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. 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.

Deflection with the exact and reduced integration (h=0.2, EI = 6.667E-04, and GAKs = 1.667E-01).

Deflection with the exact and reduced integration (h=0.4, EI = 0.533E-02, and GAKs = 3.333E-01).

Deflection with the exact and reduced integration (h=0.6, EI = 1.800E-02, and GAKs = 5.000E-01).

Deflection with the exact and reduced integration (h=0.8, EI = 4.267E-02, and GAKs = 6.667E-01).

Deflection with the exact and reduced integration (h=1.0, EI = 8.333E-02, and GAKs = 8.333E-01).

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 610. As the analyses progress, the shear-locking phenomenon disappears in Figure 11. Figures 1115 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.

Rotations with the exact and reduced integration (h=0.2, EI = 6.667E-04, and GAKs = 1.667E-01).

Rotations with the exact and reduced integration (h=0.4, EI = 0.533E-02, and GAKs = 3.333E-01).

Rotations with the exact and reduced integration (h=0.6, EI = 1.800E-02, and GAKs = 5.000E-01).

Rotations with the exact and reduced integration (h=0.8, EI = 4.267E-02, and GAKs = 6.667E-01).

Rotations with the exact and reduced integration (h=1.0, EI = 8.333E-02, and GAKs = 8.333E-01).

Additionally, condition numbers of the stiffness matrix were analyzed with exact integration and reduced integration for two cases (h=0.2 and h=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.

Condition numbers of [k](h=0.2).

Condition numbers of [k](h=0.4).

5.2. The Fixed End Supported Beam with a Nonuniform Section Applied by a Dirac Delta Function

Problem domain and properties of the cross sections are as follows: (41)Ω0,2L,L=1(42)x0,1:EI1=6.667E-04,GAKs=1.667E-01(43)x1,2:EI2=8EI=5.334E-03,GAKs=-3.5097E-01.

Governing equations are as follows: (44)ddxGAKsΨ+dwdx+f=δx-L(45)ddxEIdΨdx+GAKsΨ+dwdx=0.

Boundary conditions are as follows:(46)w0=0,dw0dx=0(47)w2L=0,dw2Ldx=0.

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 x[1,2], for which the beam depth is large, has a high ratio of h/2L. It is difficult to analyze this type of model because of the domain jump condition.

The fixed end supported beam with a nonuniform section applied by a Dirac delta function.

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.

Deflection with the exact and reduced integration (h=0.2,0.4, EI = 6.667E-04, GAKs = 1.667E-01, 8EI = 5.333E-03, and GAKs = -3.509E-01).

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 x[1,2]; however, the solutions with the reduced integration correspond to analytic solutions in the domain x[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).

Rotations with the exact and reduced integration (h=0.2, 0.4, EI=6.667E-04, GAKs = 1.667E-01, 8EI = 5.333E-03, and GAKs = -3.509E-01).

6. 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.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (NRF-2016R1A2B4016632).

Zhang X. Chen X. Yang Z. Li B. He Z. A stochastic wavelet finite element method for 1d and 2d structures analysis Shock and Vibration 2014 2014 2-s2.0-84901767048 10.1155/2014/104347 104347 Canuto C. Tabacco A. Urban K. The wavelet element method. I. Construction and analysis Applied and Computational Harmonic Analysis. Time-Frequency and Time-Scale Analysis, Wavelets, Numerical Algorithms, and Applications 1999 6 1 1 52 10.1006/acha.1997.0242 MR1664902 2-s2.0-0007372257 Ma J. Xue J. Yang S. He Z. A study of the construction and application of a Daubechies wavelet-based beam element Finite Elements in Analysis and Design 2003 39 10 965 975 2-s2.0-0037702901 10.1016/S0168-874X(02)00141-5 Xiang J. W. Chen X. F. He Z. J. Dong H. B. The construction of 1D wavelet finite elements for structural analysis Computational Mechanics 2007 40 2 325 339 10.1007/s00466-006-0102-5 MR2313467 Zbl1181.74143 2-s2.0-34247550722 Basu P. K. Jorge A. B. Badri S. Lin J. Higher-order modeling of continua by finite-element, boundary-element, meshless, and wavelet methods Computers & Mathematics with Applications 2003 46 1 15 33 10.1016/s0898-1221(03)90078-2 2-s2.0-0042825236 Zbl1288.76041 LIN J. H. Wavelets Modeling in Continuum Mechanics, Thesis Master of Science in civil engineering [M.s. thesis] 1999 Vanderbilt University Chen X. Yang S. Ma J. He Z. The construction of wavelet finite element and its application Finite Elements in Analysis and Design 2004 40 5-6 541 554 2-s2.0-0347600372 10.1016/S0168-874X(03)00077-5 Burgos R. B. Santos M. A. C. Silva R. R. e. Analysis of Beams and Thin Plates Using the Wavelet-Galerkin Method IACSIT International Journal of Engineering and Technology 7 2015 7 261 266 Becker E. B. Garey G. F. Oden J. T. Finte Elements An Introduction Volume 1 1981 Prentice-Hall Bathe K.-J. Finite element procedures Finite element procedures 1996 Reddy J. N. An introduction to the finite element method 1993 2nd McGRAW-HILL Augarde C. E. Deeks A. J. The use of Timoshenko's exact solution for a cantilever beam in adaptive analysis Finite Elements in Analysis and Design 2008 44 9-10 595 601 10.1016/j.finel.2008.01.010 2-s2.0-44249122351 Antes H. Fundamental solution and integral equations for Timoshenko beams Computers and Structures 2003 81 6 383 396 2-s2.0-0037370675 10.1016/S0045-7949(02)00452-2