Multilevel preconditioners for a quadrature Galerkin solution of a biharmonic problem

Efficient multilevel preconditioners are developed and analyzed for the quadrature finite element Galerkin approximation of the biharmonic Dirichlet problem. The quadrature scheme is formulated using the Bogner–Fox–Schmit rectangular element and the product two‐point Gaussian quadrature. The proposed additive and multiplicative preconditioners are uniformly spectrally equivalent to the operator of the quadrature scheme. The preconditioners are implemented by optimal algorithms, and they are used to accelerate convergence of the preconditioned conjugate gradient method. Numerical results are presented demonstrating efficiency of the preconditioners. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

Heinz‐type inequality and bi‐Lipschitz continuity for quasiconformal mappings satisfying inhomogeneous biharmonic equations

Let be the class of all sense‐preserving homeomorphic self‐mappings of . The aim of this paper is twofold. First, we obtain Heinz‐type inequality for (K,K′)‐quasiconformal mappings satisfying inhomogeneous biharmonic equation Δ(Δω) = g in unit disk with associated boundary value conditions and . Second, we establish biLipschitz continuity for (K,K′)‐quasiconformal mappings satisfying aforementioned inhomogeneous biharmonic equation when and are small enough.     

Interpolating minimal energy C1‐Surfaces on Powell–Sabin Triangulations: Application to the resolution of elliptic problems

In this article, we present a method to obtain a C¹‐surface, defined on a bounded polygonal domain Ω, which interpolates a specific dataset and minimizes a certain “energy functional.” The minimization space chosen is the one associated to the Powell–Sabin finite element, whose elements are C¹‐quadratic splines. We develop a general theoretical framework for that, and we consider two main applications of the theory. For both of them, we give convergence results, and we present some numerical and graphical examples. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 798–821, 2015     

梁受正弦分布压力作用的解析解

对梁受正弦分布压力作用的情况，本文以矩形截面简支梁为例，通过构造满足所有边界条件和双调和方程的应力函数给出了相应的解析解，这为求解梁受任意分布压力作用的解的问题打下了基础。     

A quadrature finite element Galerkin scheme for a biharmonic problem on a rectangular polygon

A quadrature Galerkin scheme with the Bogner–Fox–Schmit element for a biharmonic problem on a rectangular polygon is analyzed for existence, uniqueness, and convergence of the discrete solution. It is known that a product Gaussian quadrature with at least three‐points is required to guarantee optimal order convergence in Sobolev norms. In this article, optimal order error estimates are proved for a scheme based on the product two‐point Gaussian quadrature by establishing a relation with an underdetermined orthogonal spline collocation scheme. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

A C0 finite element method for the biharmonic problem without extrinsic penalization

A symmetric C ⁰ finite element method for the biharmonic problem is constructed and analyzed. In our approach, we introduce one‐sided discrete second‐order derivatives and Hessian matrices to formulate our scheme. We show that the method is stable and converge with optimal order in a variety of norms. A distinctive feature of the method is that the results hold without extrinsic penalization of the gradient across interelement boundaries. Numerical experiments are given that support the theoretical results, and the extension to Kirchhoff plates is also discussed. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1254–1278, 2014     

ELASTIC CALCULATION OF STRESSES IN RINGS USING AIRY STRESS FUNCTION

Stress calculation formulae for a ring have been obtained by using Airy stress function of the plane strain field with the decomposition of the solutions for normal stresses of Airy biharmonic equation into two parts when it is loaded under two opposite inside forces along a diameter. One part should fulfill a constraint condition about normal stress distribution along the circumference at an energy valley to do the minimum work. Other part is a stress residue constant. In order to verify these formulae and the computed results, the computed contour lines of equi-maximal shear stresses were plotted and quite compared with that of photo-elasticity test results. This constraint condition about normal stress distribution along circumference is confirmed by using Greens’ theorem. An additional compression exists along the circumference of the loaded ring, explaining the divorcement and displacement of singularity points at inner and outer boundaries.     

带有临界Sobolev指数的奇异双调和问题

LetΩR~N be a smooth bounded domain such that 0∈Ω,N≥5,2~*:=(2N)/(N-4) is the critical Sobolev exponent,and f(x) is a given function.By using the variational methods, the paper proves the existence of solutions for the singular critical in the homogeneous problemΔ~u-μu/(|x|~4)=|u|~(2~*-2)u f(x) with Dirichlet boundary condition on Ωunder some assumptions on f(x) andμ.     

Solution of the planar Newtonian stick–slip problem with the singular function boundary integral method

A singular function boundary integral method (SFBIM) is proposed for solving biharmonic problems with boundary singularities. The method is applied to the Newtonian stick–slip flow problem. The streamfunction is approximated by the leading terms of the local asymptotic solution expansion which are also used to weight the governing biharmonic equation in the Galerkin sense. By means of the divergence theorem the discretized equations are reduced to boundary integrals. The Dirichlet boundary conditions are weakly enforced by means of Lagrange multipliers, the values of which are calculated together with the singular coefficients. The method converges very fast with the number of singular functions and the number of Lagrange multipliers, and accurate estimates of the leading singular coefficients are obtained. Comparisons with the analytical solution and results obtained with other numerical methods are also made. Copyright © 2005 John Wiley & Sons, Ltd.     

双调和算子特征值问题的混合三角谱元方法

本文针对双调和算子特征值问题设计了基于混合变分形式的三角谱元逼近格式,其基函数采用指标为(-1,-1,-1)的广义Koornwinder多项式.在H~1-及H_0~1-正交谱元投影的逼近理论基础上,我们建立了双调和算子特征值与特征函数的收敛性估计;它关于网格尺寸h是最优的,关于多项式次数M是次优的.然而,在H_0~2-正交谱元投影的最优估计假设前提下,关于M的次优收敛阶估计则提升为最优.此外,Koornwinder分片多项式逼近的结果还表明,在带权Besov空间范数的度量下,对于存在着区域角点奇性的双调和算子特征值问题,谱元方法的收敛阶能达到h-型有限元方法的2倍.最后,本文的数值实验结果展示了谱元逼近格式的高效性,同时也验证了相关理论的正确性.