_{1}

A locking-free rectangular Mindlin plate element with a new multi-resolution analysis (MRA) is proposed and a new finite element method is hence presented. The MRA framework is formulated out of a mutually nesting displacement subspace sequence whose basis functions are constructed of scaling and shifting on the element domain of basic full node shape function. The basic full node shape function is constructed by extending the split node shape function of a traditional Mindlin plate element to other three quadrants around the coordinate zero point. As a result, a new rational MRA concept together with the resolution level (RL) is constituted for the element. The traditional 4-node rectangular Mindlin plate element and method is a mono-resolution one and also a special case of the proposed element and method. The meshing for the monoresolution plate element model is based on the empiricism while the RL adjusting for the multiresolution is laid on the rigorous mathematical basis. The analysis clarity of a plate structure is actually determined by the RL, not by the mesh. Thus, the accuracy of a plate structural analysis is replaced by the clarity, the irrational MRA by the rational and the mesh model by the RL that is the discretized model by the integrated.

Multi-resolution analysis (MRA) is a popular technique that has been applied in many domains such as the signal and image processing, the damage detection and health monitoring, the differential equation solution, etc. However, in the field of computational mechanics, the MRA has not been, in a real sense, fully utilized in the numerical solution of engineering problems either by the traditional finite element method (FEM) [

As is commonly known of the FEM, owing to the invariance of node number a single element contains the element that can be regarded as a monoresolution one from a MRA point of view and the FEM structural analysis is usually not associated with the MRA concept. The MRA seems to be rarely used when the FEM is employed to numerical analysis. However, it is, in fact, by means of meshing and re-meshing in which a cluster of monoresolution finite elements are assembled together artificially that the rough MRA is executed by the FEM. As we can see, in overall analysis process of a structure by the FEM, there is no mathematical foundation for the traditional finite element meshing and the finite elements are assembled together artificially. The traditional finite element model has to be re-meshed until sufficient accuracy is reached, which leads to the low computation efficiency or convergent rate. The deficiency of the FEM becomes much explicit in the accurate computation of structural problems with local steep gradient such as material nonlinear [

The great efforts have been made over the past thirty years to overcome the drawbacks of the FEM with many improved methods to come up, such as WFEM, MFM and NEM etc., which open up a transition from the monoresolution finite element method to the multiresolution finite element method featured with adjustable element node number. Although these MRA methods have illustrated their powerful capability and computational efficiency in dealing with some problems, they always have such major inherent deficiencies as the complexity of shape function construction, the absence of the Kronecker delta property of the shape function and the lack of a rigorous mathematical basis for the MRA, which make the treatment of element boundary condition complicated and the selection of element node layout empirical that substantially reduce computational efficiency. Hence, these MRA methods have never found a wide application in engineering practice just as the FEM. In fact, they can be viewed as the intermediate products in the transition of the FEM from the monoresolution to the multiresolution.

The deficiencies of all those MRA methods can be eliminated by the introduction of a new multiresolution finite element method in this paper. With respect to Mindlin plate element in the finite element stock, a new multiresolution locking-free rectangular Mindlin plate element is formulated by a new MRA, which is constituted by translated and scaled version as subspace basis functions of the basic node shape functions. The basic node shape function is then constructed from shifting to other three quadrants around a specific node of a basic element in one quadrant and joining the corresponding split node shape functions of four elements at the specific node. Hence, the node shape function construction is quite simple and clear. In addition, the proposed element method possesses a simple, clear and rigorous mathematical basis for MRA, which endows the proposed element with the resolution level (RL) that can be modulated to freely change the element node number and position in the element, adjusting structural analysis accuracy accordingly. As a result, the proposed element method can bring about substantial improvement of the computational efficiency in the structural analysis when compared with the corresponding FEM or other MRA methods.

A rectangular Mindlin plate element considered in this paper is shown as

The displacement of a classical rectangular Mindlin plate element shown in

where w,

The conventional shape functions

in which

Meanwhile, the shear angle

in which the parameters

of the plate,

Hence, substituting (6), (7) into (1), the following equations can be obtained as:

where

In addition, there exists the following relationship

As to the proposed element, the various shape functions, regarding to a node at the point coordinate of (0, 0) as shown in ^{2} or viewed as for a 1/4 split node, should be extended to the domain of

where the superscript I be denoted as b, s,

The Kronecker delta property holds for the basic full node shape functions

The basic full node shape functions

In order to carry out a MRA of a rectangular Mindlin plate structure, the mutual nesting displacement subspace sequence for a rectangular Mindlin plate element should be established. In this paper, a totally new technique is proposed to construct the MRA which is based on the concept that a subspace sequence (multi-resolution subspaces) can be formulated by subspace basis function vectors at different resolution levels whose elements- scaling function vector can be constructed by scaling and shifting on the domain

where

Here,

It is seen from Equation (14) that the nodes for the scaling process are equally spaced on the element domain of

Scaling of the basic full node shape functions on the domain of [−1, 1]^{2} (precisely on the domain of

^{2} will produce

the various node shape functions that are shown in the

Since the elements in the basis functions are linearly independent with the various scaling and the different shifting parameters, the subspaces in the subspace sequence can be established and are mutually nested, thus formulating a MRA framework, that is

where Z denoted as the positive integers,

Thus, it can be seen that the mutually nesting displacement subspace sequence

Based on the MRA, the deflection of a rectangular Mindlin plate element in the displacement subspace at the RL of

where

It can be seen that the proposed multi-resolution element is a meshfree one whose nodes are uniformly scattered at each coordinate respectively, node number and position fully determined by the RL. When the RL = 2 ´ 2, that is a traditional 4-node rectangular Mindlin plate element, Equation (16) will be reduced to Equation (1). Hence, the traditional 4-node rectangular Mindlin plate element can be regarded as a mono-resolution one and also a special case of the multi-resolution rectangular Mindlin plate element.

According to the classical assumption of a Mindlin plate theory, the generalized function of potential energy in a displacement subspace at the resolution level of

where

Substitute Equation (16) into Equation (17), the concise expression can be obtained after reassembling as follows:

where

According to the potential energy minimization principle, let

The element expressions of the stiffness matrix

where the superscript denoted as the row number of the matrix and the subscript as the aligned element node numbering (r, s). In terms of the properties of the full node shape functions, we have

in which

where

In order to carry out structural analysis, the element stiffness and mass matrices

where

where

The structural global stiffness, mass matrix

Example 1. A simply supported square plate with the geometric configuration of length L, the thickness h, the Poisson’s ratio

The displacement responses are found by the proposed element model (multiresolution element model) whose code flow is shown in

In

In light of tensor product formulation of the multi-dimensional MRA framework, the DOF of a multi-dimen- sional BSWI element will be so drastically increased from that of a one-dimensional element in an irrational way, resulting in substantial reduction of the computational efficiency. Secondly, there exists no such a parameter as the RL with a clear mathematical sense to adjust the element node number. As to the traditional monoresolution and the proposed multiresolution, the RL adjusting is more rationally and efficiently to be implemented

Four-side simply supported | Four-side clamped | ||||||||
---|---|---|---|---|---|---|---|---|---|

h/L | Analytic | Paper | BSWI23 | BSWI43 | Analytic | Paper | BSWI23 | BSWI43 | |

10^{−3} | 0.4062 | 0.4105 | 1.84E−4 | 0.4010 | 0.1265 | 0.1290 | 3.8E−5 | 0.1115 | |

10^{−2 } | 0.4062 | 0.4106 | 0.0173 | 0.4067 | 0.1265 | 0.1293 | 0.0037 | 0.1269 | |

10^{−1 } | 0.4273 | 0.4304 | 0.3510 | 0.4273 | 0.1499 | 0.1522 | 0.1152 | 0.1505 | |

0.15 | 0.4536 | 0.4566 | 0.4152 | 0.4536 | 0.1798 | 0.1802 | 0.1592 | 0.1788 | |

0.20 | 0.4906 | 0.4933 | 0.4678 | 0.4904 | 0.2167 | 0.2183 | 0.2047 | 0.2172 | |

0.30 | 0.5956 | 0.5968 | 0.5861 | 0.5957 | 0.3227 | 0.3235 | 0.3187 | 0.3246 | |

0.35 | 0.6641 | 0.6631 | 0.6579 | 0.6641 | 0.3951 | 0.3902 | 0.3898 | 0.3937 | |

than the meshing and the re-meshing to modulate element node number for the following two reasons. Firstly, the RLis based on the MRA framework which is constructed on a rigorous mathematical basis and the full node shape functions result in an integrated model, while the mesh, which resorts to the empiricism and the split node shape functions lead to a discretized model, has no MRA framework. Secondly, the stiffness matrix and the loading column vectors of the proposed element can be obtained automatically around the nodes while those of the traditional 4-node rectangular plate elements acquired by the artificially complex reassembling around the elements. It can be seen that the computational accuracy and efficiency of the proposed element model is higher than the other two.

Example 2. Two square plates are subjected to the uniform loading q with the geometric configuration of length L, thickness h, Poisson’s ratio ^{−3} and L/h = 0.3 respectively. Evaluate the deflections and moments at the center point of two plates under such boundary conditions as four sides simply supported or clamped.

The displacement responses and moments are evaluated by the proposed multiresolution element method, the traditional monoresolution element method. The central deflections and the bending moments for a thin and at hick plate with the boundary conditions as simply supported four edges (SS) and the clamped four edges (SC) under the different RLs and meshes are displayed respectively in

From the numerical example above, it is shown that based on the multi-resolution rectangular Mindlin plate element formulation, the multi-resolution finite element method is introduced, which incorporates such main steps as RL adjusting, element matrix formation, element matrix transformation from a local coordinate system to a global one and global structural matrix formation by splicing of the element matrices. Owing to the existence of the new MRA framework, the RL adjusting for the proposed method is more rationally and easily to be implemented than the meshing and re-meshing for the traditional 4-node rectangular Mindlin plate element method. Due to the basic full node shape function, the stiffness matrix and the loading column vectors of a proposed element can be automatically acquired through quadraturing around nodes in the element matrix formation step while those of the traditional 4-node rectangular Mindlin plate element obtained through complex artificially reassembling of the element matrix around the node-related elements in the re-meshing process for their 1/4 split nodes in a conventional element, which contributes a lot to computation efficiency improvement of the

RL | Mesh | Central deflection (ql^{4}/100C_{b}) | Bending moment (ql^{2}/10) | ||||||
---|---|---|---|---|---|---|---|---|---|

SS | SC | SS (central) | SC (side middle) | ||||||

Multi | Mono | Multi | Mono | Multi | Mono | Multi | Mono | ||

3 ´ 3 | 2 ´ 2 | 0.5063 | 0.5063 | 0.1480 | 0.1480 | 0.6602 | 0.6602 | −0.3551 | −0.3551 |

5 ´ 5 | 4 ´ 4 | 0.4328 | 0.4328 | 0.1403 | 0.1403 | 0.5217 | 0.5217 | −0.4761 | −0.4761 |

9 ´ 9 | 8 ´ 8 | 0.4129 | 0.4129 | 0.1304 | 0.1304 | 0.4892 | 0.4892 | −0.5028 | −0.5028 |

17 ´ 17 | 16 ´ 16 | 0.4079 | 0.4079 | 0.1275 | 0.1275 | 0.4814 | 0.4814 | −0.5104 | −0.5104 |

Analytical [ | 0.4062 | 0.1265 | 0.4789 | −0.5133 |

RL | Mesh | Central deflection (ql^{4}/100C_{b}) | Bending moment (ql^{2}/10) | ||||||
---|---|---|---|---|---|---|---|---|---|

SS | SC | SS (central) | SC (side middle) | ||||||

Multi | Mono | Multi | Mono | Multi | Mono | Multi | Mono | ||

3 ´ 3 | 2 ´ 2 | 0.7081 | 0.7081 | 0.3806 | 0.3806 | 0.8395 | 0.8395 | −0.4089 | −0.4089 |

5 ´ 5 | 4 ´ 4 | 0.6201 | 0.6201 | 0.3377 | 0.3377 | 0.5739 | 0.5739 | −0.4604 | −0.4604 |

9 ´ 9 | 8 ´ 8 | 0.5992 | 0.5992 | 0.3250 | 0.3250 | 0.5031 | 0.5031 | −0.4604 | −0.4604 |

17 ´ 17 | 16 ´ 16 | 0.5942 | 0.5942 | 0.3219 | 0.3219 | 0.4850 | 0.4850 | −0.4546 | −0.4546 |

Analytical [ | 0.5956 | 0.3227 | 0.4789 | −0.4260 |

proposed method. Moreover, since the multiresolution rectangular Mindlin plate element model of a structure usually contains much less elements than the traditional element model, thus requiring much less times of transformation matrix multiplying, the computation efficiency of the proposed method appears much higher than the traditional in the step of element matrix transformation. In addition, because of the simplicity and clarity of a full node shape function formulation with the Kronecker delta property and the solid mathematical basis for the new MRA framework, the proposed method is also superior to other corresponding MRA methods in terms of the computational efficiency, the application flexibility and extent. Hence, taking all those causes into account, the conclusion can be drawn that the multi-resolution Mindlin rectangular plate element method is more rationally, easily and efficiently to be executed, when compared with the traditional 4-node Mindlin rectangular plate element method or other corresponding MRA methods, and the proposed element would be the most accurate one formulated ever since.

In addition, multiresolution analysis (MRA) can be viewed as a technique by which amount of element details that are exposed can be modulated at a request. In the numerical analysis field, the node number a large-sized element contains could be adjusted respectively in various manners by different methods, such as the traditional FEM, the wavelet FEM (WFEM), the traditional mesh free method (MFM), the traditional natural element method (NEM), and the proposed multiresolution FEM (MFEM) etc. FEM applies the scheme of meshing and re-meshing, which is mainly relied on the empiricism, to adjust the element node number in a rough way, thus performing an irrational MRA; WFEM adopts the technique based on a mutually nested subspace sequence that is a weak mathematical basis and quite complex to be founded. MFM and NEM employ the strategy of prior artificial-selected element node layout which is also largely dependent on the empiricism. In a word, all those above or other methods are short of the parameter-resolution level (RL) with a clear mathematical sense that can be easily used to fully alter total element node number and locate element node because they do not have a simple, clear and solid mathematical basis. However, MFEM has such a simple, clear and rigorous mathematical basis that brings about the parameter RL to freely adjust total node number and locate nodes within the element. Hence, it can be said that WFEM, MFM, NEM etc are the intermediate products in the transition of the traditional FEM from the monoresolution to the multiresolution and MFEM consolidates all these irrational MRA approaches.

A new multiresolution finite element method that has both high power of resolution and strong flexibility of analysis accuracy is introduced into the field of numerical analysis. The method possesses such prominent features as follows:

・ A novel technique is proposed to construct the simple and clear basic full node shape function for a locking-free Mindlin rectangular plate element, which unveils the secret behind assembling artificially of node-related items in global matrix formation by the conventional FEM.

・ A mathematical basis for the MRA framework that is the mutually nesting displacement subspace sequence, is constituted out of the scaled and shifted version of the basic full node shape function, which brings about the rational MRA concept together with the RL.

・ The traditional 4-node Mindlin rectangular plate element and method is a monoresolution one and also a special case of the proposed. An element of superior analysis clarity surely contains more nodes when compared with that of the inferior.

・ The RL adjusting for the multiresolution Mindlinplate element model is laid on the rigorous mathematical basis while the meshing or remeshing for the monoresolution is based on the empiricism. The proposed element method can consolidate all corresponding irrational MRA approaches. Thus, the accuracy of a plate structural analysis is replaced by the clarity, the irrational MRA by the rational and the mesh by the RL that is the discretized model by the integrated.

・ A quite new concept is introduced into the FEM that the structural analysis clarity is actually determined by the RL―the density of node uniform distribution, not by the mesh.

・ With advent of the new finite element method [

The upcoming work will be focused on the treatment of interface between multiresolution elements of different RL. The interface may be extended to the bridging domain in which a transitional element can be used just as PS images of different RL. The transitional element could also be constructed by the technique of scaling and shifting of the basic full node shape function to virtual or real nodes.

Yiming Xia, (2016) Multiresolution Finite Element Method Based on a New Locking-Free Rectangular Mindlin Plate Element. World Journal of Mechanics,06,193-206. doi: 10.4236/wjm.2016.66016