一类Kirchhoff 板弯问题自适应高阶混合元方法及理论分析

本文针对Kirchhoff 板弯问题提出了一个基于高阶Hellan-Herrmann-Johnson (简记为H-H-J)方法的自适应有限元算法, 分析了它的收敛性和计算复杂度. 证明了算法在执行过程中, 相应的拟能量误差会以几何级数单调衰减, 从而得到收敛性. 利用此单调下降性质, 进一步给出了算法的计算复杂度. 推导过程中的一个关键步骤是建立基于平衡方程的单元误差表示(error indicator) 与平衡方程右端载荷震荡项(data oscillation) 的局部等价关系.     

A compact C0 discontinuous Galerkin method for Kirchhoff plates

A compact C⁰ discontinuous Galerkin (CCDG) method is developed for solving the Kirchhoff plate bending problems. Based on the CDG (LCDG) method for Kirchhoff plate bending problems, the CCDG method is obtained by canceling the term of global lifting operator and enhancing the term of local lifting operator. The resulted CCDG method possesses the compact stencil, that is only the degrees of freedom belonging to neighboring elements are connected. The advantages of CCDG method are: (1) CCDG method just requires C⁰ finite element spaces; (2) the stiffness matrix is sparser than CDG (LCDG) method; and (3) it does not contain any parameter which can not be quantified a priori compared to C⁰ interior penalty (IP) method. The optimal order error estimates in certain broken energy norm and H¹‐norm for the CCDG method are derived under minimal regularity assumptions on the exact solution with the help of some local lower bound estimates of a posteriori error analysis. Some numerical results are included to verify the theoretical convergence orders. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1265–1287, 2015     

A posteriori error estimates for the Morley plate bending element

A local a posteriori error indicator for the well known Morley element for the Kirchhoff plate bending problem is presented. The error indicator is proven to be both reliable and efficient. The technique applied is general and it is shown to have also other applications.     

混合杂交罚函数有限元方法及其应用

对一般非协调有限元,目前采用最多的两种方法是罚函数法和混合、杂交法.前一种方法总能保证收敛,但精度差,条件数和稀疏性不好;后一种方法则要满足“秩条件”才能保证收敛,故单元的构造受到很大的限制.本文提出把这两种方法结合一起的有限元方法——混合杂交罚函数法.从理论上严格证明了(在非常一般的条件下)这种新方法总是收敛的,并且其精度、条件数以及稀疏性等皆与协调元相同,也就是说都是最优的. 最后应用这一方法具体构造了一个新的九自由度任意三角形弯板单元(每个顶点给三个自由度——一个位移和两个转角),其单元刚度矩阵计算公式与旧的九自由度三角形弯板单元的计算公式相差不多.但它对任意几何形状的平板都收敛于真解,如果真解u∈H3的话,它的三个弯矩具有一阶精度,位移及两个转角均具有二阶精度.     

Optimal Control of a Variational Inequality with Application to Equilibrium Problem of an Elastic Nonhomogeneous and Anisotropic Plate Resting on Unilateral Elastic Foundation

Several optimal control problems with the same state problem—a variational inequality with a monotone operator—are considered. The inequality represents bending of an elastic, nonhomogeneous, anisotropic Kirchhoff plate resting on some unilateral elasto-rigid foundation and point supports. Both the thickness of the plate and the coefficient of the unilateral elastic foundation play the role of design variables. Cost functionals include the work of external forces (compliance), total reaction forces of the foundation or the weight of the plate. The solvability of all the problems is proved. Moreover, approximate methods for the optimal control and weight minimization problems are proposed, making use of finite elements. The design variables are approximated by piecewise affine functions. The solvability of the approximate problems is proved and some convergence analysis is presented.     

Vibration analysis of functionally graded plate with a moving mass

A hybrid method is proposed to predict the dynamic behavior of functionally graded (FG) plate subjected to a moving mass. The governing equations of motion of FG plate are derived using the Kirchhoff plate theory and Lagrange equation. Improved Rayleigh–Ritz solution is used to treat the spatial partial derivatives. Penalty method is employed to deal with the constraints, and the energy terms due to boundary conditions are included in Lagrange, hence it is not necessary to particularly consider the constraints in the modeling process. And the combination of simple polynomials and trigonometric functions is selected as the admissible functions. The advantage of this improvement in Rayleigh–Ritz method is that it is not needed to find satisfied admissible functions for different boundary conditions while the convergence of the solution is improved. Meanwhile, the method can be used to handle the versatile boundary conditions. Differential quadrature method (DQM) as a step-by-step time integration scheme is employed for discretization of temporal derivatives. The validated results show that the presented method is very reliable and efficient, and its convergence and accuracy are also better compared to finite element method for solving the dynamic problems of FG plate with moving loads (force and mass). Moreover, the influences of material properties and boundary conditions on maximum dynamic deflections are investigated, as well as moving speeds and inertial effects of loads (mass and force). Although only four edge boundary conditions are addressed in the present work, the proposed procedure is applicable for any arbitrary edge boundary conditions.     

关于TRUNC板元收敛性的注记

1.引言 自从Bergan,Argyris等提出并发展了一种称之为TRUNC三角形元的非常规板元以来,在工程界得到了广泛的应用,也在数值分析方面引起了很大的兴趣。实际计算表明,它克服了Zienkiewicz元的三平行方向剖分的限制,而且简化了刚度矩阵的形成,因此,工程界很重视。不久前石钟慈对这种板元进行了细致的数学分析,给出了很     

Model for Contact of Crack Boundaries in a Bending Plate

Closed cracks in a bending plate are studied in a two-dimensional formulation based on the Kirchhoff theory and the theory of the two-dimensional stressed state. A generalization of the end contact model is discussed which makes it possible formulate a geometrically linear contact bending problem for a plate with a rectilinear cut including the interaction of symmetric and antisymmetric deformation modes. A procedure for reducing the problem to a system of singular integral equations is described. An example of the calculations is given.     

Solution of some plate bending problems using the boundary element method

Some fundamental aspects of the boundary element method of the Kirchhoff theory of thin plate flexure are given. The direct boundary integral equation method with higher conforming properties (using first-order Hermitian interpolation for plate displacement ω, and zero-order Hermitian interpolation for angle of rotation θ, the moment M andthe equivalent shear V) are used for several computational examples. They are: square plate with simply-supported or clamped edges, the same square plate with square central opening and the cantilevered triangular plates. The results of computation as compared with some known experimental and theoritical results showed that the numerical schemes seemed to be satisfactory for the practical applications.     

A geometrically nonlinear (3,2) reduced degree-of-freedom unified zigzag laminated beam element for large deformation analysis

A geometrically nonlinear (3,2) unified zigzag beam element is developed with a reduced number of degree-of-freedom for the large deformation analysis. The main merit of the beam element model is the Kirchhoff and Cauchy shear stress solution for large deformation and large strain analysis is more accurate. The geometrically nonlinearity is considered in the calculation of the zigzag coefficients. Thus, the results of shear Cauchy stress are matching well with solid element analysis in case of the beam with aspect ratio greater than 20 under large deformation. The zigzag coefficients are derived explicitly. The Green strain and the second Piola Kirchhoff stress are used. The second Piola Kirchhoff shear stress is continuous at the interface between adjacent layers priori. The bottom surface second Piola Kirchhoff shear stress condition is used to determine the zigzag coefficient and the top surface second Piola Kirchhoff shear stress condition is used to reduce one degree-of-freedom. The nonlinear finite element equations are derived. In the numerical tests, several benchmark problems with large deformation are solved to verify the accuracy. It is observed that the proposed beam has accurate solution for beam with aspect ratio greater than 20. The second Piola Kirchhoff and Cauchy shear stress accuracy is also good. A convergence study is also presented.     

A COMBINED HYBRID FINITE ELEMENT METHOD FOR PLATE BENDING PROBLEMS

In this paper,a combined hybrid method is applied to finite element discretization of plate bending problems.It is shown that the resultant schemes are stabilized,i.e., the convergence of the schemes is independent of inf-sup conditions and any other patch test.Based on this,two new series of plate elements are proposed.     

基于膜板比拟理论的一个新的四边形薄板单元

平面弹性与板弯曲的相似理论为构造薄板单元提供了一条有效的新途径。由于避开了c1连续性的困难,使得薄板单元的构造变得简单。更为深入地讨论了该相似性理论的应用,并构造了一个新的四节点十六自由度薄板单元。数值结果表明,该单元能通过分片试验,具有良好的收敛性和精度。     

TRAPEZOIDAL PLATE BENDING ELEMENT WITH DOUBLE SET PARAMETERS

Using double set parameter method, a 12-parameter trapezoidal plate bending element is presented. The first set of degrees of freedom, which make the element convergent, are the values at the four vertices and the middle points of the four sides together with the mean values of the outer normal derivatives along four sides. The second set of degree of freedom, which make the number of unknowns in the resulting discrete system small and computation convenient are values and the first derivatives at the four vertices of the element. The convergence of the element is proved.     

A reduced local C0 discontinuous galerkin method for Kirchhoff plates

We propose and analyze a reduced local C⁰ discontinuous Galerkin (reduced LCDG) method with minimal penalization for Kirchhoff plate bending problems. The resulting linear system of the method can be solved efficiently by Gaussian elimination. Based on the key observation that the reduced LCDG method can be viewed as the localization of Hellan–Herrmann–Johnson method, we are inspired to use the techniques for dealing with the latter method to derive the well‐posedness and a priori error estimates of the reduced LCDG method. With the help of Zienkiewicz–Guzmán–Neilan element space, the a posteriori error analysis is also developed. Some numerical results are provided to demonstrate the theoretical estimates.© 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1902–1930, 2014     

A wavelet-based stochastic finite element method of thin plate bending

A wavelet-based stochastic finite element method is presented for the bending analysis of thin plates. The wavelet scaling functions of spline wavelets are selected to construct the displacement interpolation functions of a rectangular thin plate element and the displacement shape functions are expressed by the spline wavelets. A new wavelet-based finite element formulation of thin plate bending is developed by using the virtual work principle. A wavelet-based stochastic finite element method that combines the proposed wavelet-based finite element method with Monte Carlo method is further formulated. With the aid of the wavelet-based stochastic finite element method, the present paper can deal with the problem of thin plate response variability resulting from the spatial variability of the material properties when it is subjected to static loads of uncertain nature. Numerical examples of thin plate bending have demonstrated that the proposed wavelet-based stochastic finite element method can achieve a high numerical accuracy and converges fast.     

样条有限元

本文用三次B样条变分方法解规则区域上板梁组合弹性结构的平衡问题.推导出了适用于各种边界条件的统一计算格式,便于在计算机上实现.与通常有限元相比,具有计算量少、精确度高等显著特点.文中对自然边界条件作为约束条件的影响给予了考虑,并以板的弯曲问题为例说明影响极微.给出了几个数值的例子.     

非均匀薄板弯曲的精确元法

本文在阶梯折算法的基础上,提出构造有限元的新方法——精确元法.它不用一般变分原理,可适用于任意变系数正定和非正定偏微分方程.利用该方法,得到薄板弯曲一个非协调三角形单元,它具有6个自由度.文中给出证明,位移和内力均收敛于精确解,并有很好的精度.文末给出算例.算例表明利用本文的方法,内力和位移均可获得满意的结果.     

A POSTERIORI ERROR ESTIMATES FOR LOCAL CO DISCONTINUOUS GALERKIN METHODS FOR KIRCHHOFF PLATE BENDING PROBLEMS

We derive some residual-type a posteriori error estimates for the local CO discontinuous Galerkin （LCDG） approximations （[31]） of the Kirchhoff bending plate clamped on the boundary. The estimator is both reliable and efficient with respect to the moment-field approximation error in an energy norm. Some numerical experiments are reported to demonstrate theoretical results.     

Bilayer Plates: Model Reduction, Γ‐Convergent Finite Element Approximation,and Discrete Gradient Flow

The bending of bilayer plates is a mechanism that allows for large deformations via small externally induced lattice mismatches of the underlying materials. Its mathematical modeling, discussed herein, consists of a nonlinear fourth‐order problem with a pointwise isometry constraint. A discretization based on Kirchhoff quadrilaterals is devised and its Γ‐convergence is proved. An iterative method that decreases the energy is proposed, and its convergence to stationary configurations is investigated. Its performance, as well as reduced model capabilities, are explored via several insightful numerical experiments involving large (geometrically nonlinear) deformations.© 2015 Wiley Periodicals, Inc.     

On the energy flux vector for bending vibrations of a plate : PMM vol. 40, n 6, 1976, pp. 1131–1135

An expression for the energy flux vector of plate bending vibrations is obtained in invariant form. The derivation of expressions for the transverse force, bending and twisting moments in an arbitrary orthogonal coordinate system and the derivation of an orthogonality type condition for normal waves being propagated in a thin elastic strip with free edges are considered as applications.In a number of cases it turns out to be useful to consider the energy flux vector in analyzing the vibrations in systems with distributed parameters. The expressions for the Umov-Poynting vector in electrodynamics and for the energy flux vector in acoustics are well-known. An analogous vector for the bending Vibrations of a plate was mentioned only in [1 – 3], This vector is used in [1] to prove a uniqueness theorem for a two-component acoustic model consisting of an ideal compressible fluid and elastic plates in contact with it. However, the expression for the energy flux in [1] (it was later cited in [2, 3] with a reference to [1]) is erroneous. An exact expression (within the framework of the applicability of the Kirchhoff equation) is found below for the energy flux vector of the bending vibrations of a plate and some applications of the formulas obtained are mentioned.