一种高精度的裂纹奇异单元

基于广义伽辽金法的多变量有限元算法,增加了连续体力学有限元模型建立的灵活性.本文利用它,通过数值试验的对比建立了一种高精度的含奇异性的裂纹单元,并对多变量奇异元的构成进行了探讨.     

求解Stokes问题的一类新混合元格式

利用Specht九参数元的分析技巧 ,通过对Adini元增加一个自由度 ,本文构造一类求解Stokes问题的新二阶混合有限元格式 .该格式简单实用且自由度较已有的许多二阶格式少 .     

A stabilized mixed formulation for finite strain deformation for low-order tetrahedral solid elements

A stabilized mixed finite element formulation for four-noded tetrahedral elements is introduced for robustly solving small and large deformation problems. The uniqueness of the formulation lies within the fact that it is general in that it can be applied to any type of material model without requiring specialized geometric or material parameters. To overcome the problem of volumetric locking, a mixed element formulation that utilizes linear displacement and pressure fields was implemented. The stabilization is provided by enhancing the rate of deformation tensor with a term derived using a bubble function approach. The element was implemented through a user-programmable element of the commercial finite element software ANSYS. Using the ANSYS platform, the performance of the element was evaluated by comparing the predicted results with those obtained using mixed quadratic tetrahedral elements and hexahedral elements with a B-bar formulation. Based on the quality of the results, the new element formulation shows significant potential for use in simulating complex engineering processes.     

Rail-bridge coupling element of unequal lengths for analysing train-track-bridge interaction systems

This paper presents a rail-bridge coupling element of unequal lengths, in which the length of a bridge element is longer than that of a rail element, to investigate the dynamic problem of train-track-bridge interaction systems. The equation of motion in matrix form is given for a train-track-bridge interaction system with the proposed element. The first two numerical examples with two types of bridge models are chosen to illustrate the application of the proposed element. The results show that, for the same length of rail element, (1) the dynamic responses of train, track and bridge obtained by the proposed element are almost identical to those obtained by the rail-bridge coupling element of equal length, and (2) compared with the rail-bridge coupling element of equal length, the proposed element can help to save computer time. Furthermore, the influence of the length of rail element on the dynamic responses of rail is significant. However, the influence of the length of rail element on the dynamic responses of bridge is insignificant. Therefore, the proposed element with a shorter rail element and a longer bridge element may be adopted to study the dynamic responses of a train-track-bridge interaction system. The last numerical example is to investigate the effects of two types of track models on the dynamic responses of vehicle, rail and bridge. The results show that: (1) there are differences of the dynamic responses of vehicle, rail and bridge based on the single-layer and double-layer track models, (2) the maximum differences increase with the increase of the mass of sleeper, (3) the double-layer track model is more accurate.     

定常不可压Navier—Stokes方程迎风格式

本文利用齐次定解条件对定常不可压Navier—Stokes方程的非线性项进行处理,给出了相应的一种迎风Galerkin有限元算法;针对这种迎风Galerkin有限元算法,在迎风参数满足一定条件下,利用其三项式具有的一些很好性质,更简单地证明了该问题解的存在唯一性。     

Explicit error estimate for the nonconforming wilson's element

In this article, we study the explicit expressions of the constants in the error estimate of the nonconforming finite element method. We explicitly obtain the approximation error estimate and the consistency error estimate for the Wilson's element without the regular assumption, respectively, which implies the final finite element error estimate. Such explicit a priori error estimates can be used as computable error bounds.     

三维Poisson方程特征值的四种有限元解及比较

应用三维EQ1rot元、三维Crouzeix-Raviart元、八节点等参数元、四面体线性元计算三维Poisson方程的近似特征值.计算结果表明:三维EQ1rot元和三维Crouzeix-Raviart元特征值下逼近准确特征值,八节点等参数元、四面体线性元特征值上逼近准确特征值,三维EQr1ot元和三维Crouzeix-Raviart元外推特征值下逼近准确特征值.计算结果还表明三维Crouzeix-Raviart元是一种计算效率较高的非协调元.     

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.     

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

本文在阶梯折算法和精确解析法的基础上,提出构造有限元的新方法——精确元法.该方法不用变分原理,可适用于任意变系数正定和非正定偏微分方程.利用该方法,得到Reissuer板弯曲的一个非协调单元,它具有十五个自由度.由于节点位移参数仅含有挠度和转角,因此处理任意边界条件非常容易.文中给出证明,位移和内力均收敛于精确解.由精确元法所得到的单元不仅能用于厚板,也可用于薄板.文末给出四个算例.算例表明,利用本文的方法,可获得满意的结果,并有较高的数值精度.     

基于剪切效应纤维梁单元的结构非线性有限元数值模拟

基于Euler-Bernoulli梁理论的经典纤维模型忽略了剪切变形给截面带来的影响,为了得到更加精确的梁单元模型,该文基于考虑剪切效应的纤维梁单元,根据Timoshenko梁理论,推导了该纤维梁单元的刚度矩阵,并结合弹塑性增量理论,同时考虑了几何非线性和材料非线性的双重影响,建立了压弯剪复杂应力状态下结构非线性有限元...     

A beam segment element for dynamic analysis of large aqueducts

Large aqueduct structure is a complex structure that is commonly used in hydraulic engineering, especially in large-scale water conveying projects. The analysis of dynamic response for an aqueduct structure is necessary if the aqueduct is built in an earthquake area. Traditional 3D finite element analysis is time consuming and the existing simplified response method cannot take into account all the effects, such as the bending-torsion coupling effect and the constrained torsion, of the deformations of the thin wall structure of the aqueduct body. For this special structure, a simple and yet accurate model for dynamic analysis is needed. In this paper, a beam segment element is developed and used for the calculation of dynamic response for aqueduct structures. With the frame of the aqueduct being modeled using beam element, the proposed model can calculate the dynamic response of the whole aqueduct structures. Results are compared with that of a general purpose finite element analysis software using 3D finite element model. Good agreement is achieved between the two models. However, the proposed model needs less elements and much less computing time.     

非自共轭和不定椭圆问题的有限体积元方法的一致收敛性

In this paper, we prove the existence, uniqueness and uniform convergence of the solution of finite volume element method based on the P1 conforming element for non-selfadjoint and indefinite elliptic problems under minimal elliptic regularity assumption.     

The Morley element for fourth order elliptic equations in any dimensions

In this paper, the well-known nonconforming Morley element for biharmonic equations in two spatial dimensions is extended to any higher dimensions in a canonical fashion. The general n-dimensional Morley element consists of all quadratic polynomials defined on each n-simplex with degrees of freedom given by the integral average of the normal derivative on each (n-1)-subsimplex and the integral average of the function value on each (n-2)-subsimplex. Explicit expressions of nodal basis functions are also obtained for this element on general n-simplicial grids. Convergence analysis is given for this element when it is applied as a nonconforming finite element discretization for the biharmonic equation. The work was supported by the National Natural Science Foundation of China (10571006). This work was supported in part by NSF DMS-0209497 and NSF DMS-0215392 and the Changjiang Professorship through Peking University.     

Coupled FEM-BEM for elastoplastic contact problems using Lagrange multipliers

The coupling of the elastoplastic finite element and elastic boundary element methods for two-dimensional frictionless contact stress analysis is presented. Interface traction matching (boundary element approach), which involves the force terms in the finite element analysis being transformed to tractions, is chosen for the coupling method. The analysis at the contact region is performed by the finite element method, and the Lagrange multiplier approach is used to apply the contact constraints. Since the analyses of elastoplastic problems are non-linear and involve iterative solution, the reduced size of the final system of equations introduced by combining the two methods is very advantageous, especially for contact problems where the nature of the problem also involves an iterative scheme.     

BCI-代数的换位元

本文在ＢＣＩ－代数中引入了换位元，证明了换位元的若干重要性质，讨论了换位元与极大元、极小元以及半群之间的关系．     

Decoupled Mixed Element Methods for Fourth Order Elliptic Optimal Control Problems with Control Constraints

In this paper, we study the finite element methods for distributed optimal control problems governed by the biharmonic operator. Motivated from reducing the regularity of solution space, we use the decoupled mixed element method which was used to approximate the solution of biharmonic equation to solve the fourth order optimal control problems. Two finite element schemes, i.e., Lagrange conforming element combined with full control discretization and the nonconforming Crouzeix-Raviart element combined with variational control discretization, are used to discretize the decoupled optimal control system. The corresponding a priori error estimates are derived under appropriate norms which are then verified by extensive numerical experiments.     

页岩气藏渗流模型的特征混合有限元数值模拟

将特征有限元方法和混合有限元方法进行耦合,对页岩气藏渗流模型进行了数值模拟,给出了详细的误差分析,得到了最优的L~2模误差估计,并用数值实验验证了方法的有效性.     

一类优化控制问题的有限元解的误差估计和超收敛

对积分微分方程的优化控制问题进行了介绍.讨论了积分微分方程的优化控制问题的混合有限元逼近,给出了优化控制问题的有限元逼近解的误差估计和超收敛性质.     

各向异性网格下特征值问题的有限元分析

主要目的是在各向异性网格下研究二阶椭圆特征值问题的两类非协调有限元—类Wilson矩形元和Carey三角形元—的收敛性分析.通过新的技巧和方法,得到了与传统有限元网格剖分下相同的特征对的最优误差估计.推广了已有的结果.     

A Superconvergence Analysis for Finite Element Solution by the Interpolant Postprocessing on Irregular Meshes for Smooth Problem

The post-processing procedure is given by an interpolant postprocessing of the finite element solution by appropriately-defined finite dimensional subspaces. The corresponding superconvergence is established on general quasi-regular finite element partitions.