Modified mixed Ritz-DQ formulation for free vibration of thick rectangular and skew plates with general boundary conditions

Recently, the present authors proposed a simple mixed Ritz-differential quadrature (DQ) methodology for free and forced vibration, and buckling analysis of rectangular plates. In this technique, the Ritz method is first used to discretize the spatial partial derivatives with respect to a coordinate direction of the plate. The DQ method is then employed to analogize the resulting system of ordinary or partial differential equations. The mixed method was shown to work well for vibration and buckling problems of rectangular plates with simple boundary conditions. But, due to the use of conventional Ritz method in one coordinate direction of the plate, the geometric boundary conditions of the problem can only be satisfied in that direction. Therefore, the conventional mixed Ritz-DQ methodology may encounter difficulties when dealing with rectangular plates involving adjacent free edges and skew plates. To overcome this difficulty, this paper presents a modified mixed Ritz-DQ formulation in which all the natural boundary conditions are exactly implemented. The versatility, accuracy and efficiency of the proposed method for free vibration analysis of thick rectangular and skew plates are tested against other solution procedures. It is revealed that the proposed method can produce highly accurate solutions for the natural frequencies of thick rectangular plates involving adjacent free edges and skew plates using a small number Ritz terms and DQ sampling points.     

A quadrilateral plate bending element based on deformation modes

In this paper, a quadrilateral element is proposed for the analysis of thin plate bending. This element is non-conforming and consists of four-nodes and twelve degrees of freedom. A third-order field for the element displacement is written in terms of the deformation modes. Moreover, the rotational fields are obtained by utilizing the first-order Jacobean matrix. All interpolation functions are explicitly found by the presented formulation. The stiffness matrix of the element is then computed by using these functions. Finally, the accuracy of the suggested element is evaluated by solving some thin plate bending structures. Numerical findings reveal the new quadrilateral element MKQ12 is robust and accurate for analysis of thin plates.     

弹性厚板的分区广义变分原理

本文提出弹性厚板分区广义变分原理,其要点如下:1.各分区可任意定为势能区或余能区.分区势能、分区余能、分区混合变分原理是它的三种特殊形式.2.每个分区中独立变分变量的个数可任意规定.每个分区可定为单类变量区、二类变量区或三类变量区.3.每个交界线上的位移和力的连接条件可以放宽.这个原理为非协调元的厚板有限元法提供理论基础.各种厚板有限元模型可看作这个原理的特殊应用.特别是弹性厚板分区混合变分原理的提出为分区混合有限元法应用于厚板问题打下了基础.     

三维Poisson方程特征值的四种有限元解及比较

应用三维EQ1rot元、三维Crouzeix-Raviart元、八节点等参数元、四面体线性元计算三维Poisson方程的近似特征值.计算结果表明:三维EQ1rot元和三维Crouzeix-Raviart元特征值下逼近准确特征值,八节点等参数元、四面体线性元特征值上逼近准确特征值,三维EQr1ot元和三维Crouzeix-Raviart元外推特征值下逼近准确特征值.计算结果还表明三维Crouzeix-Raviart元是一种计算效率较高的非协调元.     

A simple nonconforming quadrilateral finite element

We introduce and analyze a simple nonconforming quadrilateral finite element and then we derive optimal a priori error estimates for arbitrary regular quadrilaterals. The idea of extension to some non-conforming elements for three-dimensional hexahedrons is also presented.     

位移障碍下四阶变分不等式问题的非协调有限元一般误差估计式

1.引 言 关于二阶变分不等式问题的非协调有限元逼近已有大量研究[1-5].但是,对于四阶变分不等式的研究相对而言较少[6-7].[8,9,10]给出了位移障碍问题的非协调有限元,包括C0元(如Zienkiewicz元及Adini元)和非C0元(如Morley元及De Veubeke元)逼近的理论分析及最优误差估计.经过仔细分析发现,其成功的关键技巧是充分利用上述单元的一个     

三维区域上Q_1~(rot)元的渐近展开及外推

本文研究了正方体区域上Q₁^rot非协调元渐近展开式.利用林群、吕涛等提出的有限元误差渐近展开法,获得了正方体区域上Q₁^rot非协调元特征值的误差渐近展开式.理论分析和数值实验结果表明三维Q₁^rot非协调元特征值外推公式是有效的,可以把特征值的精度从二阶提高到四阶.     

Eigen‐frequencies in thin elastic 3‐D domains and Reissner–Mindlin plate models

The eigen‐frequencies of elastic three‐dimensional thin plates are addressed and compared to the eigen‐frequencies of two‐dimensional Reissner–Mindlin plate models obtained by dimension reduction. The qualitative mathematical analysis is supported by quantitative numerical data obtained by the p‐version finite element method. The mathematical analysis establishes an asymptotic expansion for the eigen‐frequencies in power series of the thickness parameter. Such results are new for orthotropic materials and for the Reissner–Mindlin model. The 3‐D and R–M asymptotics have a common first term but differ in their second terms. Numerical experiments for clamped plates show that for isotropic materials and relatively thin plates the Reissner–Mindlin eigen‐frequencies provide a good approximation to the three‐dimensional eigen‐frequencies. However, for some anisotropic materials this is no longer the case, and relative errors of the order of 30 per cent are obtained even for relatively thin plates. Moreover, we showed that no shear correction factor is known to be optimal in the sense that it provides the best approximation of the R–M eigen‐frequencies to their 3‐D counterparts uniformly (for all relevant thicknesses range). Copyright © 2002 John Wiley & Sons, Ltd.     

非协调区域分解Lagrange乘子法的超收敛

本文讨论分析非协调区域分解Ｌａｇｒａｎｇｅ乘子法对二阶椭圆型方程Ｄｉｒｉｃｈｌｅｔ问题的有限元超收敛现象。文中通过利用积分恒等式，适宜地引进Ｌ２投影过渡以及高次插值后处理等技巧，经过一系列误差分析及估计，得到了高出半阶的超收敛结果，实现了非协调区域分解法与高精度算法的结合。     

Eigenvalue approximation from below using non-conforming finite elements

This is a survey article about using non-conforming finite elements in solving eigenvalue problems of elliptic operators,with emphasis on obtaining lower bounds. In addition,this article also contains some new materials for eigenvalue approximations of the Laplace operator,which include:1) the proof of the fact that the non-conforming Crouzeix-Raviart element approximates eigenvalues associated with smooth eigenfunctions from below;2) the proof of the fact that the non-conforming EQ rot1 element approximates eigenvalues from below on polygonal domains that can be decomposed into rectangular elements;3) the explanation of the phenomena that numerical eigenvalues λ 1,h and λ 3,h of the non-conforming Q rot1 element approximate the true eigenvalues from below for the L-shaped domain. Finally,we list several unsolved problems.     

分析各向异性板的各向同性化样条积分方程法

本文分别按Reissner理论和Kirchhoff理论导出各向导性板的各向同性化控制方程,并论证了它们间在正交各向异性简支矩形板中的相通性.在用样条积分方程法求解中采用的只是些简单的各向同性板基本解,在稀疏剖分下也能有良好的计算精度.对双参数弹性地基上的板也只需在板上虚载的取值上附加某些项而不致增加多大的工作量.     

Optimization of the static and dynamic characteristics of plates with isogrid stiffeners

In this paper, the orientation angles of stiffeners arranged in the form of isogrid configuration over a flat plate are selected to optimize the static and dynamic characteristics of these plates/stiffeners assemblies. The static characteristics are optimized by maximizing the critical buckling loads of the isogrid plate, while the dynamic characteristics are optimized by maximizing multiple natural frequencies of the stiffened plate.

A finite element model is developed to describe the statics and dynamics of Mindlin plates which are stiffened with arbitrarily oriented stiffeners. The model is used as a basis for optimizing separately or simultaneously the critical buckling loads and natural frequencies of the plates per unit volume of the plates/stiffeners assemblies.

Numerical examples are presented to demonstrate the utility of the developed model and optimization procedures. The presented approach can be invaluable in the design of plates with isogrid stiffeners for various vibration and noise control applications.     

Analysis of laminated composite plates using higher-order shear deformation plate theory and node-based smoothed discrete shear gap method

This paper presents a novel finite element formulation for static, free vibration and buckling analyses of laminated composite plates. The idea relies on a combination of node-based smoothing discrete shear gap method with the higher-order shear deformation plate theory (HSDT) to give a so-called NS-DSG3 element. The higher-order shear deformation plate theory (HSDT) is introduced in the present method to remove the shear correction factors and improve the accuracy of transverse shear stresses. The formulation uses only linear approximations and its implementation into finite element programs is quite simple and efficient. The numerical examples demonstrated that the present element is free of shear locking and shows high reliability and accuracy compared to other published solutions in the literature.     

广义协调等参元

本文根据广义协调条件,对平面应力四边形单元提出一个广义协调等参元GC-Q6.单元GC-Q6是对Wilson非协调等参元Q6的一个改进方案:单元Q6只对平行四边形网格情况能通过分片检验,对一般四边形情况却不能通过,而本文的单元GC-Q6则对一般四边形情况也能通过.当单元为平行四边形时,GC-Q6单元即退化为QS单元.算例表明,广义协调等参元的应力精度高于文献中已有的单元,对不规则网格均能保持良好的性态.     

A framework for obtaining guaranteed error bounds for finite element approximations

We give an overview of our recent progress in developing a framework for the derivation of fully computable guaranteed posteriori error bounds for finite element approximation including conforming, non-conforming, mixed and discontinuous finite element schemes. Whilst the details of the actual estimator are rather different for each particular scheme, there is nonetheless a common underlying structure at work in all cases. We aim to illustrate this structure by treating conforming, non-conforming and discontinuous finite element schemes in a single framework. In taking a rather general viewpoint, some of the finer details of the analysis that rely on the specific properties of each particular scheme are obscured but, in return, we hope to allow the reader to ‘see the wood despite the trees’.     

关于非协调位移元与杂交应力元的对应性

本文阐明了E.L.Wilson[3]等的非协调移位元与卞学鐄的杂交应力元之间所存在的对应性.     

A posteriori error estimators for the Stokes equations II non-conforming discretizations

Summary We present an a posteriori error estimator for the non-conforming Crouzeix-Raviart discretization of the Stokes equations which is based on the local evaluation of residuals with respect to the strong form of the differential equation. The error estimator yields global upper and local lower bounds for the error of the finite element solution. It can easily be generalized to the stationary, incompressible Navier-Stokes equations and to other non-conforming finite element methods. Numerical examples show the efficiency of the proposed error estimator.     

钢筋混凝土筏基的挠度分析

本假设筏基为一搁量在三参数地基上的钢筋混凝土薄板，利用伽辽金法求解了受均布荷裁的四边固定的钢筋混凝土板，所得结果颇具广泛性，且清度较高。     

Stability estimates for h-p spectral element methods for general elliptic problems on curvilinear domains

In this paper we show that the h-p spectral element method developed in [3,8,9] applies to elliptic problems in curvilinear polygons with mixed Neumann and Dirichlet boundary conditions provided that the Babuska-Brezzi inf-sup conditions are satisfied. We establish basic stability estimates for a non-conforming h-p spectral element method which allows for simultaneous mesh refinement and variable polynomial degree. The spectral element functions are non-conforming if the boundary conditions are Dirichlet. For problems with mixed boundary conditions they are continuous only at the vertices of the elements. We obtain a stability estimate when the spectral element functions vanish at the vertices of the elements, which is needed for parallelizing the numerical scheme. Finally, we indicate how the mesh refinement strategy and choice of polynomial degree depends on the regularity of the coefficients of the differential operator, smoothness of the sides of the polygon and the regularity of the data to obtain the maximum accuracy achievable. An erratum to this article is available at .     

非协调区域分解的Lagrangian乘子法

§1.引言 近年来随着并行计算机的迅速发展,求解椭圆型方程的区域分解法愈来愈引起人们的兴趣和重视.但是,目前能够见到的有限元区域分解法几乎都要求有限元空间在跨过子区域的边界时是协调的,必然限制有限元区域分解算法的优越性. [3]提出了一种非协凋区域分解法——非协调区域分解的杂交法.采用简化杂交法处理各子区域交界处的非协调性,这种方法在子区域的内部和边界采用两套不同的变量,允许内部变量在跨过各子区域的边界时不连续.但是这种方法有它的局限性,即要求边界变量在各子区域的顶点处必须保持连续性,这对推广到三维空间的情形带来很大的困难.本文提出一种非协调区域分解的Lagrangian乘子法,引进Lagrangian乘子来处理各子区域交界处的非协调性.这种方法也在子区域内部和边界采用两套不同的变量,它不仅允许内部变量在越过各子区域边界时的非协调性,并且还允许边界变量在各子区域的顶点处可以不连续,这就弥补了[3]的不足.同时,这种算法具有[3]的优点,即在不