Abstract

A stochastic finite element method based on B-spline wavelet on the interval (BSWI-SFEM) is presented for static analysis of 1D and 2D structures in this paper. Instead of conventional polynomial interpolation, the scaling functions of BSWI are employed to construct the displacement field. By means of virtual work principle and BSWI, the wavelet finite elements of beam, plate, and plane rigid frame are obtained. Combining the Monte Carlo method and the constructed BSWI elements together, the BSWI-SFEM is formulated. The constructed BSWI-SFEM can deal with the problems of structural response uncertainty caused by the variability of the material properties, static load amplitudes, and so on. Taking the widely used Timoshenko beam, the Mindlin plate, and the plane rigid frame as examples, numerical results have demonstrated that the proposed method can give a higher accuracy and a better constringency than the conventional stochastic finite element methods.

1. Introduction

The finite element method (FEM) is well accepted now in numerous industrial areas [14]. However, the validity of the results obtained by FEM can be drastically limited by the uncertainties of the parameters introduced in the model. Many factors may cause the uncertainties, such as a consequence of ignorance and spatial variability [5, 6].

On the other hand, the reliability design and analysis for engineering structures also require a kind of FE model which could deal with the variability and uncertainty. For example, a slight change in design variables may result in a poignant change of structural properties. As a powerful tool in computational stochastic mechanics, the stochastic finite element method (SFEM) appeared.

SFEM is an extension of the classical FEM to stochastic framework. The most important factor of SFEM is the uncertainty modeling. Over the past several decades, many methods were combined with FEM to approximate the uncertainties mentioned above. Chakraborty and Dey investigated the stochastic finite element simulation based on the Neumann expansion of uncertain structures subjected to earthquake [7]. A method named first-order, second-moment (FOSM) method was proposed by Cornell [8]; in fact it is a special case of perturbation method [9, 10]. Valdebenito et al. studied the application of first-order expansion considering intervening variables for estimating second-order statistics for stochastic finite element analysis [11]. The Karhunen-Loeve [12] expansion, which can be seen as a kind of the orthogonal series expansion, is also an important representation of SFEM fields. Sakamoto and Ghanem [13] proposed a method based on polynomial chaos expansion for presenting the variability. Proppe implemented multiresolution analysis for stochastic finite element problems with Karhunen-Loeve expansion [14]. Based on Karhunen-Loeve series expansion and hierarchical approach, Yang et al. applied the hierarchical stochastic finite element method for structural analysis [15]. Another method is Monte Carlo simulation (MCS) [16]; the name “Monte Carlo” was coined in the 1940s by scientists working in Los Alamos to designate a class of numerical methods based on the use of random numbers. Compared with other technologies, MCS is more direct and widely used. Taking a certain amount of samples as the priori testing data, MCS can give an accurate approximation for the stochastic field [17]. Besides, Diazdelao and Adhikari explored a method to reduce the computation cost of stochastic finite element codes based on the Gaussian process emulation [18]. Han and Bang performed the probabilistic optimal safety assessment of steel cable-stayed bridges using SFEM and expected life cycle cost concept [19]. Lang et al. investigated the prediction of heat transfer in composite materials with uncertain geometry based on the combination of the extended finite element method and the spectral stochastic finite element method [20]. Augustin and Rentrop constructed the stochastic Galerkin techniques for the random ordinary differential equations and demonstrated the advantages through numerical examples [21]. Shang and Yun proposed a stochastic finite element within a general purpose finite element analysis program ABAQUS to simulate the probabilistic structural response of stochastic materials [22].

Another important factor of SFEM is the properties of FEM element which influence the presentation directly. Monte Carlo simulation is widely used to solve complex mechanics and physics problems, particularly those involving multiple independent stochastic variables. The essential characteristic of Monte Carlo method is the use of random numbers and random variables. A random variable is a quantity that results from a repeatable process and whose actual values cannot be predicted with certainty [23]. The occurrence frequency of a sample value would be deemed to be the approximation to its probability. With the famous “law of large number,” the correctness of Monte Carlo simulation can be understood easily. Usually, a sample size of the order 30–100 could already provide estimates with acceptable accuracy [24]. On the other hand, the law makes the exactness of Monte Carlo simulation heavily depend on the sampling counter. With the growing number of sampling counters, the calculation efficiency seems to be particularly important. That is to say, the realization of MCS should be based not only on an effective computer system but also on an efficient algorithm. Therefore, a kind of new stochastic finite element method based on B-spline wavelet on the interval and MCS (BSWI-SFEM) is constructed to deal with the current problems in SFEM analysis.

It is necessary to recommend wavelet analysis and BSWI briefly before we introduce BSWI-FEM. Wavelet analysis is a new method developed in recent years [25]. The wavelet method can be viewed as a method in which the approximating function is defined by a multiresolution technique based on scaling or wavelet functions, similar to those used in signal and image processing. With its desirable characteristics, such as multiresolution properties and various basis functions for structural analysis, wavelet method is well argued by many researchers not only in numerical analysis domains [2628] but also in structural analysis fields [2931]. As a kind of wavelet, BSWI basis has the good characteristics of compact support, smoothness, and symmetry besides the multiresolution analysis. Moreover, it has the explicit expression, which will not lead to any trouble for differentiation and integration. So it is the best one among all existing wavelets in approximation of numerical calculation [32]. Based on the good characteristics of BSWI, several investigators have tried to combine it with FEM to solve practical problems [3337]. However, the above-mentioned methods, which are all constructed in a determined wavelet space, can hardly make a satisfying approximation for the influence caused by structural uncertainties.

In this paper, a new BSWI-SFEM is presented. Combining BSWI element with Monte Carlo method, a BSWI-SFEM formulation is constructed to deal with the static and vibration analysis of 1D and 2D structures. Taking the Timoshenko beam, the Mindlin plate, and the plane rigid frame as examples, numerical results indicate that the proposed method has a higher accuracy and a better constringency than the conventional SFEM.

2. The Construction of BSWI Element

2.1. One-Dimensional BSWI on the Interval

Goswami et al. constructed the scaling functions of BSWI and gave the th order wavelet formulations at scale in 1995 [39]. To any scale , the th order BSWI scaling functions must satisfy the following condition [39] in order to have at least one inner wavelet on the interval :

Then scale th order BSWI (BSWI) scaling functions and the corresponding wavelets can be evaluated by the following formulas:

Let be the scale which satisfies (1); then let in (2), for any ; according to the 0 scale th order scaling functions and wavelets given in [39], we can get the scale th order scaling functions and wavelets. There are boundary scaling functions and wavelets, inner scaling functions, and inner wavelets. The scaling functions and the wavelets on the interval can be written in vector form as follows: The corresponding eleven scaling functions and eleven wavelets of BSWI43 are given in Figure 1.

2.2. Two-Dimensional Tensor Product BSWI on the Interval

By means of tensor product, 2D BSWI can be easily extended from 1D BSWI. A separable wavelet semiorthonormal basis at scale of is constructed with tensor products of the multiresolution approximation spaces and , so the tensor product subspace can be expressed as , and the scaling functions can be expressed as where , indicating one row vector combined by the th order scaling functions at scale . , indicating another row vector combined by the th order scaling functions at scale . The wavelets functions are where and are wavelets of BSWI in (4). Figure 2(a) shows all the scaling functions of 2D tensor products BSWI, and the wavelets of tensor products BSWI are shown in Figures 2(b), 2(c), and 2(d).

2.3. The BSWI Timoshenko Beam Element

Considering the influence of shear deformation, the generalized function of potential energy of the Timoshenko beam element can be given as [32] where is the bending rigidity, is cross section rotation, is shear modulus, is the area of cross section, is the lump force, is the lump bending moment, and is the shear deformation coefficient. According to the Timoshenko beam theory, and can be independently interpolated by BSWI scaling functions as

Translating (7) into interval and taking (8) into (7), let , and the element solving equation of bending problem can be obtained as follows: where

The element solving equation of vibration problem is where is the consistent mass matrix as follows: is the stiffness matrix as follows: is density and is the vibration engine value.

2.4. The BSWI Mindlin Plate Element

According to the Mindlin plate theory, the displacements and rotations are interpolated independently, then the generalized potential energy function for the Mindlin plate bending problem is [32] where is the concentrated load. Consider where and are rotations in directions and , is the thickness of the Mindlin plate, and is the Poisson ratio. Other symbols are the same as in (7). Taking BSWI scaling functions as interpolating function is as follows: Submitting (16) into (14), according to the generalized variational principle, the BSWI element formulation of the solving equation can be obtained as follows: where we denote the following:    is similar to if the lengths of element in -axis , , , and are replaced by , , , and , respectively. Then the terms of (17) can be written as

The element solving equation of vibration problem is where is the following consistent mass matrix: is the following stiffness matrix: is density and is the vibration engine value.

2.5. BSWI Plane Rigid Frame Element

The BSWI plane rigid frame element can be seen as the simple superposition of the axial rod element and the Euler beam element. The displacement of element nodes, solving domain, and corresponding coordinates are shown in Figure 3.

It can be seen from Figure 3 that there are DOFs of one element. The physical DOFs of BSWI plane rigid element in local coordinate are while the corresponding DOFs in global coordinate () are

In order to construct the BSWI plane rigid frame element, the relationship between global and local coordinates should be obtained. Consider where Therefore,

Therefore, the BSWI plane rigid frame element for bending problem analysis in global coordinate is

The BSWI plane rigid frame element for vibration problem analysis in global coordinate is

Then, by introducing the BSWI axial rod element and the BSWI Euler beam element [40], the BSWI plane rigid frame can be obtained.

3. Stochastic Finite Element Based on B-Spline Wavelet on the Interval

The application of the MCS-FEM to large-scale structures was impossible due to its excessive computation cost a few years ago due to the huge calculation amount. Let us consider the event that the random variable can take the discrete values with probabilities , respectively. The corresponding probability distribution functions can be expressed as where is the Dirac function, and the corresponding cumulative distribution function can be expressed as

Since a discrete distribution is used to approximate the continuous distribution in the real events, the sampling counter should be big enough.

While BSWI-FEM has the property of high accuracy and efficiency [26], through inosculating this technology with MCS, the limitation of applying MCS-FEM to large-scale structures can be alleviated further.

Therefore, the stochastic wavelet finite element based on B-spline wavelet on the interval is proposed by the aid of combining Monte Carlo method with BSWI-FEM. The first step is to select the initialization parameters, such as element type and element number. Then the sampling counter should be decided, and the probability model of variables also needs to be selected. With these inputs, the calculation of BSWI-SFEM can be started as shown in Figure 4. Then the stiffness matrix and the force vector are calculated by BSWI-SFEM. Finally, through BSWI-SFEM, the results can be obtained. This process will be repeated again and again until converged. Considering the robustness of the procedure, the sampling counter is set as a variable in course; this can guarantee the convergence.

4. Numerical Examples

To stochastic analysis problem, it is hard to define an exact solution. But in one sense, the solution obtained by a large number of samples could be a good approximation of the real exact solution. In this paper, we use this kind of sample exact solution as an index to evaluate the property of the constructed BSWI-FEM. Furthermore, we will give a traditional result of finite element analysis before the BSWI-FEM analysis to validate the accuracy of the proposed method.

All the programs are implemented in the MATLAB environment. In all stochastic simulations, a computer with a 2.5 GHz Intel CPU and 2 gigabytes of memory is used, and the operating system is 32-bit Windows XP. We assume that all the parameters obey the normal distribution: the mean value of length is 1000 mm (1000 mm × 1000 mm when it is plate), the mean value of the cross section for beam is 20 mm × 12 mm, the mean value of plate thickness is 10 mm, the mean value amplitude of load is 100 N, the mean value of density is 7917 kg/m3, and the mean value of Young’s modulus is 206 MPa. The standard deviations of these parameters are 5% of their mean values.

4.1. Bending of the Timoshenko Beam

A clamp supported Timoshenko beam under a uniform load (as shown in Figure 5) is considered in this example. The parameters are mentioned above. It is demonstrated in Table 1 that the proposed BSWI element has a good agreement with the exact solution [38].

4.2. Stochastic Analysis of the Timoshenko Beam Bending

A stochastic analysis for a clamp supported Timoshenko beam under a uniform load as shown in Figure 5 is considered in this numerical example. The results of BSWI-SFEM are shown in Table 2 compared with conventional FEM. It is evident that 1 BSWI-SFEM element can achieve a similar accuracy as 320 conventional Timoshenko elements. However, the proposed method only needs 0.078 s to establish the computation, while the conventional method needs 154.875 s as shown in Figure 6.

Figure 7 presents the mean value of deflection at the center of beam which varies with the counter of sample. It can be seen that 1 BSWI-SFEM element obtains a good approximation of the center deflection compared with conventional method. Although the conventional method could approximate the sample exact solution with mesh refinement, the time consumption increases sharply. This kind of increase would result in unacceptable time consumption for the conventional method when it is used in large-scale structures. However, owing to the efficiency and accuracy of BSWI-SFEM, this problem can be greatly eased.

4.3. Stochastic Analysis of the Timoshenko Beam Vibration

Figure 8 shows the vibration analysis results of the clamp supported Timoshenko beam. BSWI results are compared with traditional finite element method. It can be seen that 1 BSWI-SFEM element can obtain a good approximation of natural frequencies compared with traditional method. Although the conventional method could approximate the sample exact solution with mesh refinement, the time consumption increases sharply. In this numerical example, the proposed BSWI method only needs 6.005 s to establish the computation with 1 element, while the conventional method needs 240.758 s to obtain the results in a similar accuracy with 320 elements.

4.4. Bending of the Mindlin Plate

In this part, the Mindlin plate bending problem as presented in Figure 9 is used to test the accuracy of the BSWI Mindlin plate element. The corresponding parameters are chosen as the mean value of these stochastic parameters mentioned above. The numerical results given in Table 3 indicate that the proposed method has a good agreement with the exact solution [38] once more.

4.5. Stochastic Analysis of the Mindlin Plate Bending

Stochastic analysis of a square Mindlin plate with clamp and simply supported boundary conditions subjected to a uniform load as shown in Figure 9 is considered in this numerical example. The stochastic bending analysis of this plate is carried out by BSWI-SFEM incorporated with the Monte Carlo simulation. The center deflection of the two kinds of boundary conditions is presented in Tables 4 and 5. The results show that the BSWI-SFEM can give a better approximation for the sample exact solution compared with the conventional method. Furthermore, the time consumption as shown in Figures 10 and 11 indicates that BSWI-SFEM is more efficient than the conventional method, which allows the proposed approach to be used for the large-scale structures. The further comparisons are presented in Figures 12 and 13. Using the mean value of output variables, the reliability calculations can be thereafter performed.

4.6. Stochastic Analysis of the Mindlin Plate Vibration

The vibration analysis results of the clamp supported Mindlin plate are shown in Figure 14. The first three natural frequency results of BSWI element are compared with conventional element. It can be seen from the results that 1 BSWI element can obtain the results in a similar accuracy with 320 traditional elements. However, 1 BSWI method only needs 10.397 s to establish the computation, while the conventional method needs 4320.385 s to obtain the results in a similar accuracy with 320 elements.

4.7. Stochastic Analysis of Plane Rigid Frame Bending

As shown in Figure 15, the plane rigid frame combined by 15 components is considered in this example. We assume that all the parameters obey the normal distribution. The corresponding material parameters are elastic modulus  N/m2, width  m, height  m, and density  kg/m3, respectively. The standard deviations of these parameters are 5% of their mean values.

As shown in Figure 16, the bending analysis results of BSWI element on points A and B are compared with BEAM3 element. Although BEAM3 element can obtain the results with high accuracy by grid refinement, it will spend much more time on computation. However, the same solving domain of plane rigid frame is discreted by only 15 BSWI elements; it can achieve the results in a similar accuracy with 480 BEAM3 elements. So the efficiency and solving precision of BSWI-SFEM plane rigid element are verified.

4.8. Stochastic Analysis of Plane Rigid Frame Vibration

The stochastic vibration analysis of plane rigid frame is carried out by BSWI-SFEM incorporated with the Monte Carlo simulation in this example. The corresponding analysis results are shown in Figure 17. BSWI element results are compared with BEAM3 element. It can be seen from the comparison that, although BEAM3 element can acquire the results with high accuracy by grid refinement, much more time is spent. BSWI-SFEM is superior, with only 15 elements, and the results with high efficiency and accuracy can be achieved.

5. Conclusions

Based on B-spline wavelet on the interval and the Monte Carlo simulation, the stochastic finite element method is developed in the paper. The proposed method can deal with the response variability resulting from the spatial variability of the structural material properties. Taking the Timoshenko beam and the Mindlin plate as examples, the corresponding model are constructed and solved. The numerical examples illustrate that the proposed method can achieve a high analytical accuracy and efficient convergence in solving the stochastic problems. Although only three kinds of stochastic BSWI elements are constructed in the paper, further extensive work is easy to be established for solving the buckling and vibration problems of structures based on the same interpolation functions and approaches. The BSWI-SFEM is feasible to carry out the reliability analysis of large-scale engineering structures.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

Special thanks should be expressed to the editor and the referees. This work was supported by the National Natural Science Foundation of China (nos. 51335006 and 51225501), the Project funded by China Postdoctoral Science Foundation (no. 2014M552432), the Fundamental Research for the Central Universities, the National Science and Technology Major Project of China (no. 2012ZX04002071), and the Program for Changjiang Scholars and Innovative Research Team in University.