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

Authors: JunHyung Jo; Yun Lee
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

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.


---

## Body

## 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 [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.


## 2. Hat Wavelets Series

### 2.1.L2 Function and the Space HmΩ

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 Figures1 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<∞ [9].Figure 1
A scaling function(φ); wavelet functions (Φ,V1=V0+W0).Figure 2
Derivatives of a scaling function and wavelet functions.The definition of Sobolev classes may be summarized as follows. A functionf 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Ω=f∣f,∂f∂x,∂f∂y,…,∂mf∂zm∈L2Ω.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=j0∞∑k=-∞∞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; j∈Z is the scaling parameter; k∈Z 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=0∞∑k=-∞∞dj,kΦj,kxThe subspace ofL2R spanned by the scaling functions φkx=φx-k is obtained as V0=Spankφk(x)¯. The different subspaces are related as shown in(5)V0⊂V1⊂V2⊂V3⊂⋯⊂L2Rwithu∈Vj.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=0∞∑k=-∞∞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 forj. 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 (inV0) are(9)φ0x=x,x∈0,12-x,x∈1,20,otherwise.Wavelet functions (inW0) are(10)Φ0x=-x,x∈0,0.53x-2,x∈0.5,1-3x+4,x∈1,1.5x-2,x∈1.5,20,otherwise,(11)Φ1x=x,x∈0,0.5-2x+32,x∈0.5,1x-32,x∈1,1.50,otherwise,(12)Φ2x=x,x∈0,0.51-x,x∈0.5,10,otherwise.


## 2.1.L2 Function and the Space HmΩ

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 Figures1 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<∞ [9].Figure 1
A scaling function(φ); wavelet functions (Φ,V1=V0+W0).Figure 2
Derivatives of a scaling function and wavelet functions.The definition of Sobolev classes may be summarized as follows. A functionf 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Ω=f∣f,∂f∂x,∂f∂y,…,∂mf∂zm∈L2Ω.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=j0∞∑k=-∞∞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; j∈Z is the scaling parameter; k∈Z 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=0∞∑k=-∞∞dj,kΦj,kxThe subspace ofL2R spanned by the scaling functions φkx=φx-k is obtained as V0=Spankφk(x)¯. The different subspaces are related as shown in(5)V0⊂V1⊂V2⊂V3⊂⋯⊂L2Rwithu∈Vj.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=0∞∑k=-∞∞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 forj. 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 (inV0) are(9)φ0x=x,x∈0,12-x,x∈1,20,otherwise.Wavelet functions (inW0) are(10)Φ0x=-x,x∈0,0.53x-2,x∈0.5,1-3x+4,x∈1,1.5x-2,x∈1.5,20,otherwise,(11)Φ1x=x,x∈0,0.5-2x+32,x∈0.5,1x-32,x∈1,1.50,otherwise,(12)Φ2x=x,x∈0,0.51-x,x∈0.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 functionsw1 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)findw∈H1such  thatBw1,v1=lv1,v1∈H1finds∈H1such  thatBs1,v2=lv2,v2∈H1.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Ψ+dwdx≡V,(18)EIdΨdx≡M.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=0∞∑k=-∞∞dj,kΦj,kx,w2x=∑k=-∞∞ckφkx+∑j=0∞∑k=-∞∞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=wi1∫0Ωχ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χi1l⋯Q10+⋯+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χ2dx⋯dχndxχn+1χn+2⋯χ2n(32)Λ2Ω≔0102⋯0nχ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 [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 Ωs∈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.


## 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 Figure3, 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.Figure 3
The fixed end supported beam applied by a Dirac delta function.Figure 4
7 functions and 8 sections (φx: scaling functions, Φx: wavelet functions).Figure 5
15 functions and 16 sections (φx: scaling functions, Φx: wavelet functions).As shown in Figure6, 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. 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.Figure 6
Deflection with the exact and reduced integration (h=0.2,EI = 6.667E-04, and GAKs = 1.667E-01).Figure 7
Deflection with the exact and reduced integration (h=0.4,EI = 0.533E-02, and GAKs = 3.333E-01).Figure 8
Deflection with the exact and reduced integration (h=0.6,EI = 1.800E-02, and GAKs = 5.000E-01).Figure 9
Deflection with the exact and reduced integration (h=0.8,EI = 4.267E-02, and GAKs = 6.667E-01).Figure 10
Deflection with the exact and reduced integration (h=1.0,EI = 8.333E-02, and GAKs = 8.333E-01).As shown in Figure11, 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. Figures 11–15 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.Figure 11
Rotations with the exact and reduced integration (h=0.2,EI = 6.667E-04, and GAKs = 1.667E-01).Figure 12
Rotations with the exact and reduced integration (h=0.4,EI = 0.533E-02, and GAKs = 3.333E-01).Figure 13
Rotations with the exact and reduced integration (h=0.6,EI = 1.800E-02, and GAKs = 5.000E-01).Figure 14
Rotations with the exact and reduced integration (h=0.8,EI = 4.267E-02, and GAKs = 6.667E-01).Figure 15
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.Figure 16
Condition numbers of[k](h=0.2).Figure 17
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)x∈0,1:EI1=6.667E-04,GAKs=1.667E-01(43)x∈1,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 Figure18, 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.Figure 18
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. Figure19 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.Figure 19
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 domainx∈[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).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).


## 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 Figure3, 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.Figure 3
The fixed end supported beam applied by a Dirac delta function.Figure 4
7 functions and 8 sections (φx: scaling functions, Φx: wavelet functions).Figure 5
15 functions and 16 sections (φx: scaling functions, Φx: wavelet functions).As shown in Figure6, 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. 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.Figure 6
Deflection with the exact and reduced integration (h=0.2,EI = 6.667E-04, and GAKs = 1.667E-01).Figure 7
Deflection with the exact and reduced integration (h=0.4,EI = 0.533E-02, and GAKs = 3.333E-01).Figure 8
Deflection with the exact and reduced integration (h=0.6,EI = 1.800E-02, and GAKs = 5.000E-01).Figure 9
Deflection with the exact and reduced integration (h=0.8,EI = 4.267E-02, and GAKs = 6.667E-01).Figure 10
Deflection with the exact and reduced integration (h=1.0,EI = 8.333E-02, and GAKs = 8.333E-01).As shown in Figure11, 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. Figures 11–15 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.Figure 11
Rotations with the exact and reduced integration (h=0.2,EI = 6.667E-04, and GAKs = 1.667E-01).Figure 12
Rotations with the exact and reduced integration (h=0.4,EI = 0.533E-02, and GAKs = 3.333E-01).Figure 13
Rotations with the exact and reduced integration (h=0.6,EI = 1.800E-02, and GAKs = 5.000E-01).Figure 14
Rotations with the exact and reduced integration (h=0.8,EI = 4.267E-02, and GAKs = 6.667E-01).Figure 15
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.Figure 16
Condition numbers of[k](h=0.2).Figure 17
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)x∈0,1:EI1=6.667E-04,GAKs=1.667E-01(43)x∈1,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 Figure18, 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.Figure 18
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. Figure19 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.Figure 19
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 domainx∈[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).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.



---

*Source: 4868050-2017-08-27.xml*