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

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

非协调曲边四边形十二自由度平板弯曲单元*

本文利用精确元法[1],给出一个十二自由度曲边四边形板弯曲单元.该方法不需要变分原理,适用于任意正定和非正定偏微分方程.利用这个方法,单元之间的协调条件很容易满足,仅须位移和内力在单元节点上连续,即可保证所得到的解收敛于精确解.利用本文方法所获得的解,无论是位移还是内力可同时有二阶收敛精度.文末给出数值算例.表明了本文所得到的单元有非常好的精度.     

弹性力学平面问题的精确元法

本文在阶梯折算法[1]和精确解析法[2]的基础上,提出构造有限元的新方法——精确元法.该方法不用一般变分原理,可适用于任意变系数正定和非正定偏微分方程.利用该方法得到弹性力学平面问题的一个非协调任意四边形单元.它具有八个自由度.由于没有采用雅可比变换,该单元可以蜕化为三角形单元,在工程中使用起来较为方便.文中给出收敛性证明.文末给出算例,位移和应力均给出较好的结果,在单元的节点上有较好的数值精度.     

非协调曲边四过形十二自由度平板变曲单元

本利用精确元法^[1]，给出一个十二自由曲边四边形板弯曲单元，该方法不需要变分原理，适用于任意正定和非正定偏微分方程。利用这个方法，单元之间的协调条件很容易满足，仅须位移和内力的单元节点上连续，即可保证所得的解收敛于精确解。利用本方法所获得的解，无论是位移还是内可力同时有二阶收敛精度。末给出数值算例，表明了本所得到的单元有非常好的精度。     

精确有限元法

本文提出构造有限单元的新方法——精确有限元法.它可以求解在任意边界条件下任意变系数正定或非正定偏微分方程。文中给出它的收敛性证明和计算偏微分方程的一般格式。用精确元法所得到的单元是一个非协调元,单元之间的相容条件容易处理.与相同自由度普通有限元相比,由精确元法所得到的解的高阶导数具有较高的收敛精度.文末给出数值算例,所得到的结果均收敛于精确解,并有较好的数值精度.     

非均匀弹性地基圆薄板大挠度问题的一般解*

本文在文[1]的基础上提出了一个新的方法可用于求解任意变系数非线性常微分方程组.文中导出了任意轴对称载荷和不同边界条件下的非均匀弹性地基圆薄板大变形的一般解,并给出了收敛于精确解的证明.问题最后可归结为求解一个仅含有三个未知量的非线性代数方程组.该方法和其它方法比较,具有收敛范围大,计算简便迅速等特点.文末给出算例表明内力和位移均可得到满意的结果,验证了本文理论的正确性.     

非均匀圆柱壳非线性轴对称变形一般解

本文利用阶梯折算法[1],得到了非均匀圆柱壳非线性轴对称变形的一般解.文中导出了在任意轴对称载荷下求解非均匀圆柱壳非线性弯曲的位移和内力的一般公式,并给出一致收敛于精确解的证明.问题最后归结为求解二元一次代数方程组,文末给出算例.算例表明,无论内力和位移都可得到满意的结果,并收敛于精确解.     

非均匀变截面梁动力响应的一般解

本文利用精确解析法[1]给出非均匀变截面梁在任意谐振荷载和边界条件下的动力响应的一般解.问题最后归结为求解一个非正定微分方程.对于这一问题用一般变分法求解失败,利用本文的方法,这个问题的一般解可以写成解析的形式.因此对优化特别方便.本文给出收敛性证明.文末给出的算例表明本文的方法可获得满意的结果.     

非均匀双向加肋圆柱壳非线性轴对称变形的精确解析法

本文首次利用精确解析法分析了环向和纵向加肋非均匀圆柱壳在任意载荷和边界条件下非线性轴对称变形问题.导出了一致收敛于精确解的位移和内力解析表达式,文中给出收敛性问题.问题最后归结为求解二元一次代数方程组,计算既简便又迅速.文末给出四个数值算例表明,本文提出的方法,可以得到满意的结果.     

任意变系数微分方程的精确解析法

工程中的许多问题归结为求解任意变系数微分方程的解.本文首次提出精确解析法,用以求解任意变系数微分方程在任意边界条件下的解.文中还给出精确解析法的一般计算格式,得到了一致收敛于精确解及其任意阶导数的解析表达式,并给出收敛性证明.文末给出四个算例,均得到较好的结果,证明了本文理论的正确性.     

弱有限元方法在线弹性问题中的应用

本文考虑弱有限元（简称WG）方法在线弹性问题中的应用.WG方法是传统有限元方法的推广,用于偏微分方程的数值求解.和传统有限元一样,它的基本思想源于变分原理.WG方法的特点是使用在剖分单元内部和剖分单元边界上分别有定义的分片多项式函数（即弱函数）作为近似函数来逼近真解,并针对弱函数定义相应的弱微分算子代入数值格式进行计算.除此之外,WG方法允许在数值格式中引进稳定子以实现近似函数的弱连续性.WG方法具有允许使用任意多边形或多面体剖分,数值格式与逼近函数构造简单,易于满足相应的稳定性条件等优点.本文考虑WG方法在求解线弹性问题中的应用.围绕线弹性问题数值求解中常见的三个问题,即：数值格式的强制性,闭锁性,应力张量的对称性介绍WG方法在线弹性问题求解中的应用.     

A multilevel correction method for Stokes eigenvalue problems and its applications

In this paper, a new multilevel correction scheme is proposed to solve Stokes eigenvalue problems by the finite element method. This new scheme contains a series of correction steps, and the accuracy of eigenpair approximation can be improved after each step. In each correction step, we only need to solve a Stokes problem on the corresponding fine finite element space and a Stokes eigenvalue problem on the coarsest finite element space. This correction scheme can improve the efficiency of solving Stokes eigenvalue problems by the finite element method. As applications of this multilevel correction method, a multigrid method and an adaptive finite element technique are introduced for Stokes eigenvalue problems. Some numerical results are given to validate our schemes. Copyright © 2013 John Wiley & Sons, Ltd.     

ERROR ESTIMATES FOR SPARSE OPTIMAL CONTROL PROBLEMS BY PIECEWISE LINEAR FINITE ELEMENT APPROXIMATION

Optimization problems with L¹-control cost functional subject to an elliptic partial differential equation(PDE)are considered.However,different from the finite dimensiona l¹-regularization optimization,the resulting discretized L¹norm does not have a decoupled form when the standard piecewise linear finite element is employed to discretize the continuous problem.A common approach to overcome this difficulty is employing a nodal quadrature formula to approximately discretize the L¹-norm.In this paper,a new discretized scheme for the L¹-norm is presented.Compared to the new discretized scheme for L¹-norm with the nodal quadrature formula,the advantages of our new discretized scheme can be demonstrated in terms of the order of approximation.Moreover,finite element error estimates results for the primal problem with the new discretized scheme for the L¹-norm are provided,which confirms that this approximation scheme will not change the order of error estimates.To solve the new discretized problem,a symmetric Gauss-Seidel based majorized accelerated block coordinate descent(sGS-mABCD)method is introduced to solve it via its dual.The proposed sGS-mABCD algorithm is illustrated at two numerical examples.Numerical results not only confirm the finite element error estimates,but also show that our proposed algorithm is efficient.     

非均匀弹性地基圆板轴对称弯曲的一般解

本文在阶梯折算法的基础上,提出一个新的方法——精确解析法,得到了非均匀弹性地基圆板弯曲的一般解.文中导出了在任意轴对称载荷和边界条件下求解非均匀弹性地基圆板和中心带孔圆板弯曲的一般公式,并给出一致收敛于精确解的证明.文中得到的一般解可直接计算无弹性地基圆板的弯曲问题.问题最后归结为求解一个二元一次代数方程.文末给出算例,算例表明无论内力和位移均可得到满意的结果.     

Free vibration and buckling analysis of plates using B-spline wavelet on the interval Mindlin element

A finite element method (FEM) of B-spline wavelet on the interval (BSWI) is used in this paper to solve the free vibration and buckling problems of plates based on Reissner–Mindlin theory. By aid of the high accuracy of B-spline functions approximation for structural analysis, the proposed method could obtain a fast convergence and a satisfying numerical accuracy with fewer degrees of freedoms (DOF). The numerical examples demonstrate that the present BSWI method achieves the high accuracy compared to the exact solution and others existing approaches in the literatures. The BSWI finite element has potential to be used as a numerical method in analysis and design.     

A finite volume spectral element method for solving magnetohydrodynamic (MHD) equations

In this paper, the coupled equations in velocity and magnetic field for unsteady magnetohydrodynamic (MHD) flow through a pipe of rectangular section are solved using combined finite volume method and spectral element technique, improved by means of Hermit interpolation. The transverse applied magnetic field may have an arbitrary orientation relative to the section of the pipe. The velocity and induced magnetic field are studied for various values of Hartmann number, wall conductivity and orientation of the applied magnetic field. Comparisons with the exact solution and also some other numerical methods are made in the special cases where the exact solution exists. The numerical results for these sample problems compare very well to analytical results.