三角形REISSNER-MINDLIN板元

本文提出构造无自锁现象的Reissuer－Mindlin板元的一个一般性方法.此方法将剪切应变用它的适当的插值多项式代替,当板厚趋于零时这对应于插值点的Kirchhoff条件,因而单元无自锁现象.根据这种方法我们构造两个三角形元──一个3节点元和一个6节点元,并给出数值结果.     

基于台劳展式的矩形Reissner-Mindlin板元

1.引言 Rdissner-Mindlin板模型放弃了经典板模型的Kirchhoff假说,考虑了剪切变形,能应用于更广泛的板问题。Reissner-Mindlin板模型的挠度与转角是相互独立的,单元只需具有c~0连续性,这一点优于需要具有c~1连续性的Kirchhoff板元,但一个严重困难是普通c~0元,尤其是低阶c~0元,当板厚趋于零时不收敛,这就是所谓的自锁现象(locking)。近年来,研究避免自锁现象的Reissner-Mindlin模型板元吸引了不少的注     

具有几何对称性的12参数矩形板元

1 引言 三角形板元中,形式最简单的是九参数元,节点参数是单元三个顶点上的函数值和两个一阶偏导数值,非协调九参三角形板元的研究取得了丰硕成果,根据不同方法已构造出多种收敛性能好的单元.相比之下,矩形板元的研究较少见报道.矩形板元中形式最简单的是12参元,节点参数是单元4个顶点上的函数值和两个一阶偏导数值,这类似于九参三角形板元.常见的12参矩形板元是ACM元,其形函数空间是完整3次多项式空间加上两个4次多项式的基函数,ACM元是C°元,但位移形函数的外法向导数平均值在单元间不连续,这类似于Zienkiewicz九参三角形板元,但由于矩形单元的特殊形状,ACM元是收敛的.龙驭球教授等在[1]中提出一种12参矩形广义协调元,其位移形函数的外法向导数平均值在     

Reissner-Mindlin板问题带约束非协调旋转Q_1有限元方法

本文利用带约束非协调旋转Q_1元逼近Reissner-Mindlin板问题中旋度的两个分量.并分别选择Wilson元、双线性元和带约束非协调旋转Q_1元逼近挠度,相应地选取不连续的矢量值分片线性函数空间、最低阶旋转Raviart-Thomas元空间和矢量值分片常数函数空间为离散的剪应力空间,在矩形网格上构造了三个板元.通过证明一个离散的Korn不等式,并借助MITC4元的解构造了旋度、挠度和剪应力一个具有某种特殊且关键的可交换性的插值.再利用Helmholtz分解分析相容性误差.我们证明了这三个矩形元在能量范数意义下与板厚无关的一致最优收敛性.数值算例验证了我们的理论结果.     

高维空间中一类Ginzburg-Landau型泛函的极小元的渐近分析

在空间H1,pg(Ω,Rn)中讨论如下一类变系数Ginzburg-Landau型泛函Eε(Ω)=∫Ωa(x)p|Δu|p+14εpb(x)(|u|2-β2(x))2dx的极小元列的渐近性质.这里2≤p0,m≤a(x),b(x),β(x)≤M,且a(x),b(x),β(x)是光滑函数.研究了当ε→0时极小元的渐近性态,证明了极小元列在H1,pg(Ω,Rn)中强收敛于某个元素,且得到了该元素所满足的微分方程边值问题.     

量子群主Tilting模的张量积及其滤过

A=z[υ]Ω,Ω是Z[υ]的由υ-1和奇素数p生成的理想.U是A上的量子代数.令φp是p次分圆多项式,B=A/(φp),Γ是商代数B关于理想(ξ-1)的完备化,式中ξ是p次本原根.对λ∈X+,Mr(λ)表首权为λ的不可分解Uг-Tilting模(称为主Ur模).本文给出了量子群主Ur模的张量积定理.对p≥2(h-1),在p2室中描述了量子群主Ur模好滤过滤过商之首权的分布状态及其滤过重数.作为例子,对秩1型和A2型的量子群情形给出了p2室中一般位置室主Ur模好滤过的分解模式.     

关于非线性椭圆边值问题解的存在性的注

利用非线性增生映射值域的扰动理论,本文研究了与P拉普拉斯算子△p相关的非线性椭圆边值问题@在Ls(Ω)空间中解的存在性,其中2>sp>2nn+1且n1.@-Δpu+|u(x)|p-2u(x)+g(x,u(x))=fa.e.x∈Ω-〈υ,|u|p-2u〉=0a.e.x∈Γ其中f∈Ls(Ω)给定,ΩRn,n1,Δpu=div(|u|p-2u)为P拉普拉斯算子,υ为Γ的外法向导数,g∶Ω×R→R满足Caratheodory条件.本文所讨论的方程及所用的方法是对以往一些工作的补充和延续.     

数值求积与有限元外推

其中 Ω是 R~2中的有界区域.设(?)={T}是对Ω的单元剖分网格,V_h 是相应的有限元空间.问题(1.1)的标准有限元方法,是寻求 Ritz 投影 u~h∈V~h 满足     

有限元方法的渐近展开及其外推

首先将 Ω 剖分成大三角形域 Ω_k,Ω_k 走的顶点亦为Ω的角点.对诸Ω_k 进行一致剖分(参看下页图),设Ω_k~h={τ_(kl)},Ω~h=(?)Ω_k~h={τ_(kl)}k,l;h_k 表示Ω_k~h 中单元的直径,h=(?)(h_k),S_0~h(Ω~h)表示线性有限元空间;u~h 和 u~l 分别表示问题 (P) 的解 u 在 S_0~k(Ω~h)上的 Ritz 投影及其线性插值,G_(z_(?))~h(z) 表示问题 (P) 的 Green 函数 G_(z_0)(z) 在 S_0~h(Ω~h)上的 Ritz 投影.由[3]知     

简单自反MTS(υ,λ)的存在谱

一个Mendelsohn三元系MTS(υ,λ)=(X,B)被称作是自反的,如果它与它的逆(X,B-1)是同构的,其中B-1=｛〈z,y,x〉;〈x,y,z〉∈B.在[2]中已给出了简单自反MTS((υ,1)的存在谱,即υ≡0,1(mod3),υ3且υ≠6.本文讨论一般λ的情况,并得到简单自反MTS(υ,λ)的存在谱是λυ(υ-1)≡0(mod3);λυ-2,υ3且(υ,λ)≠(6,1);(6,3).     

On the connection between some linear triangular Reissner-Mindlin plate bending elements

Summary. We consider three triangular plate bending elements for the Reissner-Mindlin model. The elements are the MIN3 element of Tessler and Hughes [19], the stabilized MITC3 element of Brezzi, Fortin and Stenberg [5] and the T3BL element of Xu, Auricchio and Taylor [2, 17, 20]. We show that the bilinear forms of the stabilized MITC3 and MIN3 elements are equivalent and that their implementation may be simplified by using numerical integration of reduced order. The T3BL element is shown to be essentially the same as the MIN3 and stabilized MITC3 elements with reduced integration. We finally introduce a general stabilized finite element formulation which covers all three methods. For this class of methods we prove the stability and optimal convergence properties. Received November 4, 1996 / Revised version received May 29, 1997 / Published online January 27, 2000     

Virtual Elements for the Reissner‐Mindlin plate problem

In this work, we present a virtual element method for the approximation of the plate bending problem in the Reissner‐Mindlin formulation. The proposed method follows the MITC approach of the FEM context. We construct a family of VEM spaces with arbitrary degree of accuracy that satisfies the conditions of the MITC philosophy. We perform some numerical tests which allow us to assess the convergence and the robustness of the method.     

Optimal $\mathcal{L}^2$ error bounds for MITC3 type element

Summary. In this paper, we derive the optimal error bounds for the stabilized MITC3 element [3], the MIN3 type element [7] and the T3BL element [8]. In this way we have solved the problem proposed recently in [5] in a positive manner. Moreover, we estimate the difference between stabilized MITC3 and MIN3 and show it is of order uniform in the plate thickness. Received May 31, 2000 / Revised version received April 2, 2001 / Published online September 19, 2001     

Upper and lower error bounds for plate-bending finite elements

Summary. We compare the robustness of three different low-order mixed methods that have been proposed for plate-bending problems: the so-called MITC, Arnold-Falk and Arnold-Brezzi elements. We show that for free plates, the asymptotic rate of convergence in the presence of quasiuniform meshes approaches the optimal O(h) for MITC elements as the thickness approaches 0, but only approaches for the latter two. We accomplish this by establishing lower bounds for the error in the rotation. The deterioration occurs due to a consistency error associated with the boundary layer – we show how a modification of the elements at the boundary can fix the problem. Finally, we show that the Arnold-Brezzi element requires extra regularity for the convergence of the limiting (discrete Kirchhoff) case, and show that it fails to converge in the presence of point loads. Received June 9, 1998 / Published online December 6, 1999     

A p-version MITC finite element method for Reissner–Mindlin plates with curved boundaries

We consider the approximation of Reissner–Mindlin plates with curved boundaries, using a p-version MITC finite element method. We describe in detail the formulation and implementation of the method, and emphasize the need for a Piola-type map in order to handle the curved geometry of the elements. The results of our numerical computations demonstrate the robustness of the method and suggest that it gives near exponential convergence when the error is measured in the energy norm. For the robust computation of quantities of engineering interest, such as the shear force, the proposed method yields very satisfactory results without the need for any additional post-processing. Comparisons are made with the standard finite element formulation, with and without post-processing.     

A posteriori error analysis for conforming MITC elements for Reissner-Mindlin plates

This paper establishes a unified a posteriori error estimator for a large class of conforming finite element methods for the Reissner-Mindlin plate problem. The analysis is based on some assumption (H) on the consistency of the reduction integration to avoid shear locking. The reliable and efficient a posteriori error estimator is robust in the sense that the reliability and efficiency constants are independent of the plate thickness . The presented analysis applies to all conforming MITC elements and all conforming finite element methods without reduced integration from the literature.

     

The hp-MITC finite element method for the Reissner–Mindlin plate problem

The popular MITC finite elements used for the approximation of the Reissner–Mindlin plate are extended to the case where elements of non-uniform degree p distribution are used on locally refined meshes. Such an extension is of particular interest to the hp-version and hp-adaptive finite element methods. A priori error bounds are provided showing that the method is locking-free. The analysis is based on new approximation theoretic results for non-uniform Brezzi–Douglas–Fortin–Marini spaces, and extends the results obtained in the case of uniform order approximation on globally quasi-uniform meshes presented by Stenberg and Suri (SIAM J. Numer. Anal. 34 (1997) 544). Numerical examples illustrating the theoretical results and comparing the performance with alternative standard Galerkin approaches are presented for two new benchmark problems with known analytic solution, including the case where the shear stress exhibits a boundary layer. The new method is observed to be locking-free and able to provide exponential rates of convergence even in the presence of boundary layers.     

Spectrum and states of the bcs hamiltonian in a finite domain. I. Spectrum

The BCS Hamiltonian in a finite cube with periodic boundary condition is considered. The special subspace of pairs of particles with opposite momenta and spin is introduced. It is proved that, in this subspace, the spectrum of the BCS Hamiltonian is defined exactly for one pair, and for n pairs the spectrum is defined through the eigenvalues of one pair and a term that tends to zero as the volume of the cube tends to infinity. On the subspace of pairs, the BCS Hamiltonian can be represented as a sum of two operators. One of them describes the spectra of noninteracting pairs and the other one describes the interaction between pairs that tends to zero as the volume of the cube tends to infinity. It is proved that, on the subspace of pairs, as the volume of the cube tends to infinity, the BCS Hamiltonian tends to the approximating Hamiltonian, which is a quadratic form with respect to the operators of creation and annihilation.     

Analysis of a 3‐node mixed‐shear‐projected triangular element method for Reissner–Mindlin plates

This article analyses an existing 3‐node hybrid triangular element, called MiSP3, for Reissner–Mindlin plates which behaves robustly in numerical benchmark tests (Ayad, Dhatt, and Batoz, Int J Numer Method Eng 42 (1998), 1149–1179). Based on Hellinger‐Reissner variational principle and the mixed shear interpolation/projection technique of MITC family, the MiSP3 element uses continuous piecewise linear polynomials for the approximations of displacements and a piecewise‐independent equilibrium mode for the approximations of bending moments/shear stresses. Due to local elimination of the parameters of moments/stresses, the element is almost of the same computational cost as the conforming linear triangular displacement element. We derive uniform stability and convergence results with respect to the plate thickness. The main tools of our analysis are the self‐equilibrium relation of the moments/stresses approximations, the properties of the mixed shear interpolation and the discrete Helmholtz decomposition of the shear stress approximation. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 241–258, 2017