多尺度六面体单元网格连接方法研究

提出多结点六面体单元的结点形函数构建方法,解决了数字化分析时两个多尺度有限元网格在接触界面由于单元结点不匹配而导致的结点属性不能连续传递问题.首先将全局坐标下的六面体单元及其表面上多结点通过等参逆变换转成局部坐标下的规则六面体单元及其表面上多结点;在规则单元中,以每个结点为基点,分别沿三个正交的局部坐标方向在单元内寻找...     

The improvement of numerical modeling in the solution of incompressible viscous flow problems using finite element method based on spherical Hankel shape functions

In this paper, the finite element method with new spherical Hankel shape functions is developed for simulating 2‐dimensional incompressible viscous fluid problems. In order to approximate the hydrodynamic variables, the finite element method based on new shape functions is reformulated. The governing equations are the Navier‐Stokes equations solved by the finite element method with the classic Lagrange and spherical Hankel shape functions. The new shape functions are derived using the first and second kinds of Bessel functions. In addition, these functions have properties such as piecewise continuity. For the enrichment of Hankel radial basis functions, polynomial terms are added to the functional expansion that only employs spherical Hankel radial basis functions in the approximation. In addition, the participation of spherical Bessel function fields has enhanced the robustness and efficiency of the interpolation. To demonstrate the efficiency and accuracy of these shape functions, 4 benchmark tests in fluid mechanics are considered. Then, the present model results are compared with the classic finite element results and available analytical and numerical solutions. The results show that the proposed method, even with less number of elements, is more accurate than the classic finite element method.     

基于平面偶应力-Reissner/Mindlin板比拟的偶应力有限元

偶应力理论的有限元列式面临本质性的C1连续性困难. 平面偶应力理论和Reissner/Mindlin板弯曲理论之间的比拟关系表明这两个理论系统的有 限元的同一性，而R/M板有限元并不存在C1连续性困难. 因此，研究将R/M板单元转化为具有一般位移自由度的平面偶应力单元的一般方法. 根据这一方法，将典型的8节点Serendipity型R/M板单元Q8S转化为一个4节点12 自由度的四边形平面偶应力单元，数值结果表明该单元具有良好的精度和收敛性     

基于等几何分析的比例边界有限元方法

提出了一种具有比例边界有限元的半解析特性和等几何分析的几何特性的新方法。该新方法是在比例边界有限元框架中用NURBS曲线或曲面精确描述域边界几何形状,同时域边界位移场采用描述几何形状的NURBS形函数等参构造。这种新方法具有比例边界有限元固有的径向解析特性和NURBS的高阶连续性的优点。数值算例显示,与传统的比例边界有限元相比,基于等几何分析的比例边界有限元方法提高了域边界单元和域内应力场的连续性,减少了计算自由度。应用此方法可以用较少的计算自由度获得更高连续阶和更高精度的位移、应力和应变场。     

有限元方法中的光滑积分伪弱形式

提出了一种光滑积分伪弱形式,将光滑积分拓展至被积函数非偏导项求解。结合光滑应变技术和伪弱形式,可实现有限元系统方程统一光滑积分求解,即对刚度矩阵和质量矩阵中的应变矩阵和形函数矩阵均可进行光滑积分处理,并转化为光滑子域的边界积分。光滑积分伪弱形式与光滑应变技术比较,增加了形函数矩阵不定积分处理过程,且没有降低有限元求解对形函数连续性的要求。不过,伪弱形式改变了单元积分的求解形式,连续质量矩阵求解也无需坐标映射和雅可比矩阵计算。以轴对称二维问题为研究对象,结果表明极度不规则三角形和四边形单元光滑积分伪弱形式在静态和动态有限元方程求解中也具有很好的精度。     

基于余弦样条函数的有限点阵法

应用节点影响域的概念,提出了基于余弦样条函数的有限点阵方法.利用余弦样条函数的性质,通过张量积的形式构造余弦样条函数正则解空间,用于逼近场函数,余弦样条基的线性组合使得边界条件处理如同有限元法一样方便.余弦样条与边界型方法结合,可用于求解不规则域问题.数值实例的计算结果表明,文中方法避免了高阶多项式构造形函数所带来的数值振荡,解的连续性不受限制,为改进计算精度而加密点阵所导致的计算量增加量较少,计算收敛快.     

节点梯度光滑有限元配点法

配点法构造简单、计算高效,但需要用到数值离散形函数的高阶梯度,而传统有限元形函数的梯度在单元边界处通常仅具有C0连续性,因此无法直接用于配点法分析.本文通过引入有限元形函数的光滑梯度,提出了节点梯度光滑有限元配点法.首先基于广义梯度光滑方法,定义了有限元形函数在节点处的一阶光滑梯度值,然后以有限元形函数为核函数构造了有...     

平面广义四节点等参元GQ4及其性能探讨

广义节点有限元是将传统有限元方法中的节点广义化,在不增加节点个数的前提下,仅通过提高广义节点的插值函数的阶次,从而达到提高有限元解精度的目的.与现有的p型和hp型有限元不同,在这种新的有限元中,节点自由度全部定义在节点处,在理论与程序实现上与传统有限元方法具有很好的相容性,传统有限元方法是这种新方法的广义节点退化为0阶时的特殊情形.文中主要讨论了这一新方法的四节点等参元（记为GQ4）的形式.对GQ4进行的各种数值试验表明,所发展的广义四节点等参单元具有精度高且无剪切自锁与体积自锁等的特点.     

Non-local finite element analysis of damped beams

Non-local viscoelastic beam models are used to analyse the dynamics of beams with different boundary conditions using the finite element method. Unlike local damping models the internal force of the non-local model is obtained as weighted average of state variables over a spatial domain via convolution integrals with spatial kernel functions that depend on a distance measure. In the finite element analysis, the interpolating shape functions of the element displacement field are identical to those of standard two-node beam elements. However, for non-local damping, nodes remote from the element do have an effect on the energy expressions, and hence on the damping matrix. The expressions of these direct and cross damping matrices may be obtained explicitly for some common spatial kernel functions and Euler–Bernoulli beam theory. Alternatively numerical integration may be applied to obtain solutions. Examples are given where the eigenvalues are compared to the exact solution for a pinned–pinned beam to demonstrate the convergence of the finite element method. The results for beams with other boundary conditions are used to demonstrate the versatility of the finite element technique.     

MODE III 2-D FRACTURE ANALYSIS BY THE SCALED BOUNDARY FINITE ELEMENT METHOD

The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique that combines the advantages of the finite element method and the boundary element method with unique properties of its own. This method has proven very efficient and accurate for determining the stress intensity factors （SIFs） for mode I and mode II two-dimensional crack problems. One main reason is that the SBFEM has a unique capacity of analytically representing the stress singularities at the crack tip. In this paper the SBFEM is developed for mode III （out of plane deformation） two-dimensional fracture anMysis. In addition, cubic B-spline functions are employed in this paper for constructing the shape functions in the circumferential direction so that higher continuity between elements is obtained. Numerical examples are presented at the end to demonstrate the simplicity and accuracy of the present approach for mode Ⅲ two-dimensional fracture analysis.     

The Edge-Based Face Element Method for 3D-Stream Function and Flux Calculations in Porous Media Flow

We present a velocity-oriented discrete analog of the partial differential equations governing porous media flow: the edge-based face element method. Conventional finite element techniques calculate pressures in the nodes of the grid. However, such methods do not satisfy the requirement of flux continuity at the faces. In contrast, the edge-based method calculates vector potentials along the edges, leading to continuity of fluxes. The method is algebraically equivalent with the popular block-centered finite difference method and with the mixed-hybrid finite element method, but is algorithmically different and has the same robustness as the more conventional node-based velocity-oriented method. The numerical examples are computed analytically and may, therefore, be considered as an 'heuristic proof' of the theory and its practical applicability for reservoir engineering and geohydrology.     

Hermitian space-time finite elements for estuarine mass transport

Space-time finite element solutions of the convection–dispersion equation using higher-order nodal continuity and Hermitian polynomial shape functions are described. Five separate elements ranging from a complete linear element with C^0,0 nodal continuity to a complete first-order Hermitian element with C^1,1 nodal continuity are subjected to detailed analysis. Wave deformation analyses identify the source of leading or trailing edge oscillations, trailing edge oscillations being the major source of difficulty. These observations are confirmed by numerical experiments which further demonstrate the potential of higher-order nodal continuity. The performance of the complete first-order Hermitian element is quite satisfactory and measurably superior to the linear element.     

基于局部插值的结构动力模型降阶方法

提出了一种基于局部插值对大型结构有限元模型的特征值问题进行降阶的方法. 该方法通过局部插值将复杂结构的有限元模型中节点的位移用凝聚点的位移插值来表示, 从而得到用插值函数表示的简化基向量, 实现对结构广义特征值问题的降阶. 为了提高降阶模型的精度, 采用非协调元的插值函数作为局部插值函数来弱化凝聚后的结构刚度, 并且在有限元模型上进行逆迭代, 对得到的降阶后的广义特征值问题的特征值和特征向量进行改善. 为了提高模型降阶的效率, 采用规整网格包围整个结构生成均匀的凝聚点, 高效地确定了有限元模型中节点所依附的凝聚点. 最后, 对 3 个机床部件的模态分析验证了提出的简化方法的高效性和可行性.      

A Direct Method for the Generation of Shape Functions of Any Desired Continuity for Rectangular Finite Elements

Abstract

A direct method is presented for the generation of shape functions for rectangular finite elements, with any desired continuity of the shape functions across the interelement boundaries, i.e., the shape functions can be of any prescribed class C^N, N = 0, 1, 2, 3 …. The method is illustrated with seven examples, four of which are previously published elements. The other three examples represent new elements, two of which are conforming C¹ elements, and the third one is of class C². One of the C1 elements and the C² element reproduce exactly all the terms up to the sextic polynomial.     

US-FE-LSPIM四边形单元及其在几何非线性问题中的应用

为了提高在网格畸变时的数值计算精度,基于非对称有限单元的概念,提出US-FE-LSPIM四边形单元。该单元是利用传统的四节点等参元形函数集和FE—LSPIM四边形单元形函数集分别作为检验函数和试函数而构成。前者用于满足单元间和单元内的位移连续性要求,后者用于满足位移完备性要求。该单元结合了有限单元法和无网格法的优点,能...     

Numerical analysis of Biot's consolidation process by radial point interpolation method

An algorithm is proposed to solve Biot's consolidation problem using meshless method called a radial point interpolation method (radial PIM). The radial PIM is advantageous over the meshless methods based on moving least-square (MLS) method in implementation of essential boundary condition and over the original PIM with polynomial basis in avoiding singularity when shape functions are constructed. Two variables in Biot's consolidation theory, displacement and excess pore water pressure, are spatially approximated by the same shape functions through the radial PIM technique. Fully implicit integration scheme is proposed in time domain to avoid spurious ripple effect. Some examples with structured and unstructured nodes are studied and compared with closed-form solution or finite element method solutions.     

非匹配结点的有限单元等参插值方法及其应用研究

针对两零件的异构网格单元结点在接触界面不能相互匹配导致结点属性不能连续过渡和传递的问题,提出非匹配结点的有限单元等参插值方法,通过构建非匹配结点的形函数和修正原结点的形函数,将结点属性值的影响范围限制在可控的局部区域,从而实现两异构网格结点属性在接触界面的连续过渡和传递。通过两个异构的四边形单元网格的结点属性在接触界面的过渡实例和啮合齿轮的接触分析应用,验证了该方法的正确性和有效性。     

多边形有限单元形函数的比较研究

多边形有限单元形函数有wachspress插值、Laplace插值和平均值插值三种类型.本文对三种多边形有限单元形函数的性质作了比较研究,给出了三种形函数各自的优点和局限性.Waclaspress和Laplace形函数是有理函数形式,而平均值形函数是无理函数形式.三种形函数均满足单位分解性、线性完备性,且在单元边界上呈线性.在三角形单元上,它们都等价于三角形面积坐标插值.在矩形单元上,Wachspress和Laplace形函数等价于双线性多项式插值形函数.Wachspress和平均值形函数适用于任意凸多边形单元,Laplace形函数更适用于圆内接多边形单元.Wachspress形函数不能推广到含有边节点的单元,平均值形函数可以直接推广到含有边节点的单元.数值试验,验证了本文理论分析的结论.     

多边形有限元研究进展

有限元法是数值求解偏微分方程边值问题的重要方法,采用不规则多边形单元网格, 可以方便有效地模拟材料的力学性能, 又使得区域网格剖分变得灵活方便. 特别是对于复杂的几何形状, 多边形单元网格具有更大的优势. 本文对国内外有关多边形有限元法的最新进展作了初步的总结和评述, 主要以基于位移法的多边形有限元为主.论述了多边形有限元的发展历史, 给出了多边形单元上的Wachspress插值、Laplace插值和重心坐标的一些最新研究成果. 与经典有限元法形函数为多项式形式不同, 多边形单元的形函数为有理函数或者无理函数形式. 多边形单元插值形函数满足线性完备性, 可以再现线性位移场, 像经典有限元法一样直接施加本质边界条件; 插值函数在多边形的边界上是线性的,确保不同单元间的自动协调. 不同单元的插值形函数表达公式形式统一, 方便混合单元网格计算的程序编写. 提出了多边形有限元法今后需要研究的问题.