Darcy-Stokes问题的一个非协调有限元分析

本文针对Darcy-Stokes问题构造一个新的非协调的矩形单元,证明该单元的收敛性,并得到速度能量模的二阶收敛结果,数值算例验证了理论分析结果.     

多变量有限元:相容性与模式优化

本文给出基于非协调试解函数的多变量有限元的非线性相容分析,导出能量相容条件和单元优化条件.提出了实施这些条件的杂交元的优化方法,建立了单元优化格式和多变量参数匹配原理.研究结果已应用于弹性力学等问题.     

柱塞泵油膜压力分析的高效有限单元法

利用有限单元法,数值分析轴向柱塞泵缝隙高度可变时阀门板上油膜的压力分布.缝隙中油压变化的确定,通常是应用Pasynkow 润滑理论中的Reynolds方程.使用基于有限单元法的自编程序,数值地求解该方程.为了提高结果的精度,基于解的残差估计,应用了加密的自适应网格.给出了依赖于泵的几何条件和工作参数的计算结果.     

复合材料层合板的弯曲与振动

本文应用应力杂交有限元方法分析了复合材料层合板的弯曲与振动.在本文中,首先根据修正的余能变分原理,构造了一个适合于复合材料层合板特点的矩形应力杂交板弯曲单元.在单元内,分层假设应力参数,在单元的边界上,根据YNS理论的假设确定边界位移场.这样使得构造出来的单元不仅能够考虑横向剪切变形的影响和局部扭曲效应,而且具有较少的自由度数.其次,用此单元求解了层合板的弯曲与振动问题,并将计算结果与精确解进行了比较,比较表明二者非常接近.这说明了在计算方面本文单元具有较高的精确度.     

改进的非协调广义混合单元及性能分析

非协调广义混合单元最突出的特点是避免了传统混合单元中系数矩阵主对角线上存在零元素的问题，因此位移和应力结果的收敛是稳定的.以最小势能原理和H R变分原理为基础，联合增强假设应变理论建立了新的8结点非协调广义混合单元.一方面，该单元保持了已有非协调广义混合单元的全部优点；另一方面，该单元简化了积分计算.数值实例表明，改进的非协调广义混合单元的数值结果精度高，计算速度快并且对单元的几何扭曲敏感度低.     

基于扩展多面体组合单元的非规则颗粒材料离散元方法北大核心CSCD

非规则颗粒材料广泛地存在于自然界和工业生产中,其复杂的几何形态对力学性质有显著的影响.为构建更接近真实颗粒形态的理论模型,以扩展多面体为基本单元,发展了扩展多面体组合单元.为验证扩展多面体组合单元的可靠性,分别对凸形三棱柱单元、凹形正倒锥体单元在平底漏斗中的卸料过程进行了离散元模拟,并与试验结果进行比较分析,得到其具有较好的一致性.在此基础上,对不同形态的组合单元进行堆积和卸料离散元模拟,研究了颗粒形状对堆积分数、卸料流量和休止角的影响.结果表明,颗粒形状越复杂,颗粒之间的互锁效应越显著,颗粒系统更加稳定.扩展多面体组合单元的有效应用,为离散元数值模拟描述任意形态颗粒材料提供了一种新的构建方法.     

带无序保护单元的系统可靠性分析

在寿命分布为离散型和连续型的情形下,对带无序保护单元的系统进行可靠性分析,给出了此类系统部分可靠性参数的计算方法和计算结果.同时还给出了一个计算实例,用以说明上述分析结果的具体应用.     

各向异性网格下Wilson元的超收敛性分析

在各向异性网格下研究了二阶椭圆边值问题的Wilson有限元方法,利用单元构造的特殊性和一些新的技巧得到相应的超逼近和超收敛结果.数值算例的结果与理论分析是相吻合的.     

基于带有虚拟单元聚类分析的广义DEA方法

广义DEA方法是相对于样本单元对决策单元进行评价,如何选择样本单元对于评价的结果有一定的影响.有时可以选择决策单元集的子集作为样本单元集,这样更容易被决策者所接受.通过构造符合一定条件的虚拟决策单元与所有决策单元一起采取最短距离法进行聚类分析,用包含虚拟决策单元的类作为样本单元集,可以使评价结果与传统DEA方法相比能够进一步体现出相对效率之间的差异.     

Stokes问题的非协调有限元分析

构造了两个非协调有限元,对Stokes问题满足离散B-B条件,单元对速度和压力具有二阶收敛性,数值实验验证了理论分析结果.     

A characterization of neutral elements by the exclusion of sublattices

The main result of this paper gives a characterization of neutral elements in lattices by the exclusion of seventeen types of sublattices containing it. The proof involves several generalizations of distributive, standard and neutral elements in lattices which are introduced using certain modular identities.     

An upwind-mixed method on changing meshes for two-phase miscible flow in porous media

Two-phase miscible flow in porous media is governed by a system of nonlinear partial differential equations. In this paper, the upwind-mixed method on dynamically changing meshes is presented for the problem in two dimensions. The pressure is approximated by a mixed finite element method and the concentration by a method which upwinds the convection and incorporates diffusion using an expanded mixed finite element method. The method developed is shown to obtain almost optimal rate error estimate. When the method is modified we can obtain the optimal rate error estimate that is well known for static meshes. The modification of the scheme is the construction of a linear approximation to the solution, which is used in projecting the solution from one mesh to another. Finally, numerical experiments are given.     

Three Dimension Quasi-Wilson Element for Flat Hexahedron Meshes

The well known Wilson‘s brick is only convergent for regular cuboid meshes. In this paper a quasi-Wilson element of three dimension is presented which is convergent for any hexahedron meshes. Meanwhile the element is anisotropic, that is it can be used to any flat hexahedron meshes for which the regular condition is unnecessary.     

MORTAR型旋转Q_1元的多重网格方法

本文讨论了mortar型旋转Q_1元的多重网格方法.证明了W循环的多重网格法是最优的,即收敛率与网格尺寸及层数无关.同时给出了一种可变的V循环多重网格算法,得到了一个条件数一致有界的预条件子.最后,数值试验验证了我们的理论结果.     

拟协调元的位移函数及节点误差

直接从拟协调元的应变关系式出发,构造具有明确物理意义的幂级数形式的位移函数,从而得出拟协调元的常应变和线性应变系数是唯一确定的,它只能收敛到常应变的结论;刚性位移项可采用多种构造方法,不同的方法得出的节点参数与单元的本身的节点参数存在不同阶次的误差,这与常规位移法有限元不同。     

An alternative formulation of the boundary element method

The paper suggests an alternative formulation of the Boundary Element Method, in which singular solutions generated by unit dislocations are required and moreover the stresses at the interior points of the body are directly computed from the boundary quantities, without passing through the displacements. Relationships between the singular solutions for unit dislocation and unit force are derived.     

2-Absorbing and Weakly 2-Absorbing Elements in Multiplicative Lattices

In this article we study weakly 2-absorbing elements and 2-absorbing elements in C-lattices. Next we characterize C-lattices in which every nonzero proper element is a (weakly) 2-absorbing element. We also establish a new characterization for principal element domains in terms of 2-absorbing elements.     

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

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

Anisotropic interpolation and quasi-Wilson element for narrow quadrilateral meshes

In this paper an anisotropic interpolation theorem is presentedthat can be easily used to check the anisotropy of an element.A kind of quasi-Wilson element is considered for second-orderproblems on narrow quadrilateral meshes for which the usualregularity condition _K/h_K c₀ > 0 is not satisfied, whereh_K is the diameter of the element K and _K is the radius of thelargest inscribed circle in K. Anisotropic error estimates ofthe interpolation error and the consistency error in the energynorm and the L²-norm are given. Furthermore, we give a Poincaréinequality on a trapezoid which improves a result of eniek.     

Guaranteed velocity error control for the pseudostress approximation of the Stokes equations

The pseudostress approximation of the Stokes equations rewrites the stationary Stokes equations with pure (but possibly inhomogeneous) Dirichlet boundary conditions as another (equivalent) mixed scheme based on a stress in H(div) and the velocity in L². Any standard mixed finite element function space can be utilized for this mixed formulation, e.g., the Raviart‐Thomas discretization which is related to the Crouzeix‐Raviart nonconforming finite element scheme in the lowest‐order case. The effective and guaranteed a posteriori error control for this nonconforming velocity‐oriented discretization can be generalized to the error control of some piecewise quadratic velocity approximation that is related to the discrete pseudostress. The analysis allows for local inf‐sup constants which can be chosen in a global partition to improve the estimation. Numerical examples provide strong evidence for an effective and guaranteed error control with very small overestimation factors even for domains with large anisotropy.© 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1411–1432, 2016