A precise singular finite element for crack analysis

The multi-variable finite element algorithm based on the generalized Galerkin’smethod is more flexible to establish a finite element model in the continuum mechanies.Byusing this algorithm and numerical tests a new singular finite element for elasto-plasticfracture analysis has been formulated.The results of numerical tests show that the newelement possesses high accuracy and good performance.Some rules for formulating amulti-variable singular finite element are also discussed in this paper.     

A 3D multi-field element for simulating the electromechanical coupling behavior of dielectric elastomers

We propose a multi-field implicit finite element method for analyzing the electromechanical behavior of dielectric elastomers. This method is based on a four-field variational principle, which includes displacement and electric potential for the electromechanical coupling analysis, and additional independent fields to address the incompressible constraint of the hyperelastic material. Linearization of the variational form and finite element discretization are adopted for the numerical implementation. A general FEM program framework is developed using C ++ based on the open-source finite element library deal.II to implement this proposed algorithm. Numerical examples demonstrate the accuracy, convergence properties, mesh-independence properties, and scalability of this method. We also use the method for eigenvalue analysis of a dielectric elastomer actuator subject to electromechanical loadings. Our finite element implementation is available as an online supplementary material.     

RTM充模过程数值模拟的隐式有限元算法

建立了基于欧拉方法描述树脂传递模塑(RTM)工艺充模过程的基本数学方程,并采用有限元隐式时间积分方法对基本方程进行了数值求解．编制了基于隐式有限元算法及传统有限元控制体算法的程序,通过具体算例比较了这两种算法的优缺点．与传统的有限元控制体法相比,该文提出的隐式有限元算法能节省计算时间,特别适合于单元、节点数目多的情况．隐式有限元算法是一种纯有限元方法,不需要使用控制体积技术,采用该算法计算出的流动前沿与时间步长无关。     

岩质圆形隧洞围岩应力场弹塑性新解

针对动态接触问题的有限元并行计算，提出了一种新的接触算法. 新算法引入局部拉氏 乘子技术来计算接触力. 由于同时考虑了无穿透的接触约束条件和相邻接触对的相互影响， 较之广泛使用的罚参数法，新算法使接触约束条件和系统平衡方程得到更充分的满足. 虽然 为提高接触计算精度而在局部采用了迭代技术，但算法仍然具有较高的效率，且与显式时间 积分方案完全相容. 此外，通过构造专门的区域分解方案，实现了将现有为串行程序开发的 搜索算法平滑移植到并行环境的目标. 数值算例表明，所提出的接触算法具有很好的并行性， 在保证了接触问题并行计算精度的同时，取得了满意的并行效率.     

框架结构屈曲的精确有限元求解

基于屈曲微分控制方程的一般解,构造了Euler梁在轴力作用下的精确形函数,建立了用于框架结构屈曲分析的精确有限单元,得到了单元刚度矩阵和几何刚度矩阵的显式表达,并提出了基于常规特征值计算的迭代算法以确定屈曲载荷及相应失稳模态的精确解. 研究表明, 对于线性稳定性分析而言,常规框架有限单元可视为精确有限单元的一种近似. 若采用精确单元,无需进行网格细分就可以获得精确的屈曲载荷和失稳模态. 数值算例证明了新单元和算法的效率和精度.      

An iterative parallel algorithm of finite element method

In this paper, a parallel algorithm with iterative form for solving finite element equation is presented. Based on the iterative solution of linear algebra equations, the parallel computational steps are introduced in this method. Also by using the weighted residual method and choosing the appropriate weighting functions, the finite element basic form of parallel algorithm is deduced. The program of this algorithm has been realized on the ELXSI-6400 parallel computer of Xi'an Jiaotong University. The computational results show the operational speed will be raised and the CPU time will be cut down effectively. So this method is one kind of effective parallel algorithm for solving the finite element equations of large-scale structures.     

A Eulerian approach to the finite element modelling of neo-Hookean rubber material

A Eulerian approach is applied to the finite element modelling of neo-Hookean rubber material. Two major problems are encountered. The first problem is the construction of an algorithm to calculate stresses in the rubber material from velocities instead of displacements. This problem is solved with an algorithm based on the definition of the velocity gradient. The second problem is the convection of stresses through the finite element mesh. This problem is solved by adapting the so-called Taylor-Galerkin technique. Solutions for both problems are implemented in a finite element program and their validity is shown by test problems. Results of these implementations are compared with results obtained by a standard Lagrangian approach finite element package and good agreement has been found.     

无穷域势流问题的有限元／差分线法混合求解

提出了一种将有限元和差分线法相结合求解无穷域势流问题的算法。用两同心圆将求解域划分为存在重叠的有限和无限两个区域,在有限和无限域上分别用有限元和差分线法求解Laplace方程边值问题。用差分线法推导出的关系式修正有限元方程,求解该方程组从而得到原问题的解。本算法将求解无穷域问题转化为代数特征值问题和有限域内线性方程组的...     

一类针对高阶单元接触搜索算法：系统搜索算法

提出一类基于非线笥约束优化理论的高副搜索算法。通用的一些接触搜索算法一般针对低阶单元,它们具有较高的计算效率,但不适应于高阶单元。随着工整中对计算精度要求不断提高,高阶单元逐步被引入到接触问题计算程序中。因此,给出针对亢介单元而人有与现有的低阶算法相比扩计算效率的接触索算法具有现实意义,本文采用系统索方法和非线性约束优化理论,建立一种高阶单元接触点副搜索算法。计算表明,本算法具有较高的计算疚和良好     

Iterative methods for scattering problems in isotropic or anisotropic elastic waveguides

We consider the time-harmonic problem of the diffraction of an incident propagative mode by a localized defect, in an infinite elastic waveguide. We propose several iterative algorithms to compute an approximate solution of the problem, using a classical finite element discretization in a small area around the perturbation, and a modal expansion in the unbounded straight parts of the guide. Each algorithm can be related to a so-called domain decomposition method, with an overlap between the domains. Specific transmission conditions are used, so that at each step of the algorithm only the sparse finite element matrix has to be inverted, the modal expansion being obtained by a simple projection, using a bi-orthogonality relation. The benefit of using an overlap between the finite element domain and the modal domain is emphasized. An original choice of transmission conditions is proposed which enhances the effect of the overlap and allows us to handle arbitrary anisotropic materials. As a by-product, we derive transparent boundary conditions for an arbitrary anisotropic waveguide. The transparency of these new boundary conditions is checked for two- and three-dimensional anisotropic waveguides. Finally, in the isotropic case, numerical validation for two- and three-dimensional waveguides illustrates the efficiency of the new approach, compared to other existing methods, in terms of number of iterations and CPU time.     

On a finite element CFD algorithm for compressible,viscous and turbulent aerodynamic flows

This paper develops and analyses individual construction aspects of an efficient and accurate finite element algorithm for prediction of viscous and turbulent flow fields of impact in aerodynamics. The theoretical construction employs a Taylor weak statement (TWS) for coincident embedding of stability mechanisms within a classic Galerkin finite element formulation of semi-discrete approximation error orthogonalization. A wide variety of the stabilizing mechanisms of independently derived CFD algorithms are contained within the TWS theory. An implicit construction that meets the requirement of efficient convergence to steady state is developed. The theoretical asymptotic error estimates of the TWS finite element algorithm for supersonic and viscous boundary layer flows are verified. Application to a three-dimensional turbulent flow is cited.     

A DOMAIN DECOMPOSITION ALGORITHM WITH FINITE ELEMENT-BOUNDARY ELEMENT COUPLING

A domain decomposition algorithm coupling the finite element and the boundary element was presented. It essentially involves subdivision of the analyzed domain into sub-regions being independently modeled by two methods, i.e., the finite element method (FEM) and the boundary element method (BEM). The original problem was restored with continuity and equilibrium conditions being satisfied on the interface of the two sub-regions using an iterative algorithm. To speed up the convergence rate of the iterative algorithm, a dynamically changing relaxation parameter during iteration was introduced. An advantage of the proposed algorithm is that the locations of the nodes on the interface of the two sub-domains can be inconsistent. The validity of the algorithm is demonstrated by the consistence of the results of a numerical example obtained by the proposed method and those by the FEM, the BEM and a present finite element-boundary element (FE-BE) coupling method.     

On Taylor weak statement finite element methods for computational fluid dynamics

A Taylor series augmentation of a weak statement (a ‘Taylor weak statement’ or ‘Taylor-Galerkin’ method) is used to systematically reduce the dispersion error in a finite element approximation of the one-dimensional transient advection equation. A frequency analysis is applied to determine the phase velocity of semi-implicit linear, quadratic and cubic basis one-dimensional finite element methods and of several comparative finite difference/finite volume algorithms. The finite element methods analysed include both Galerkin and Taylor weak statements. The frequency analysis is used to obtain an improved linear basis Taylor weak statement finite element algorithm. Solutions are reported for verification problems in one and two dimensions and are compared with finite volume solutions. The improved finite element algorithms have sufficient phase accuracy to achieve highly accurate linear transient solutions with little or no artificial diffusion.     

Local and parallel finite element algorithms based on two-grid discretization for steady Navier-Stokes equations

Local and parallel finite element algorithms based on two-grid discretization for Navier-Stokes equations in two dimension are presented. Its basis is a coarse finite element space on the global domain and a fine finite element space on the subdomain. The local algorithm consists of finding a solution for a given nonlinear problem in the coarse finite element space and a solution for a linear problem in the fine finite element space, then droping the coarse solution of the region near the boundary. By overlapping domain decomposition, the parallel algorithms are obtained. This paper analyzes the error of these algorithms and gets some error estimates which are better than those of the standard finite element method. The numerical experiments are given too. By analyzing and comparing these results, it is shown that these algorithms are correct and high efficient.     

算子自定义小波弹性板单元构造及自适应分析

提出一类用于分析弹性板问题的算子自定义小波弹性板单元构造方法。该方法的优点在于根据工程问题的求解需要灵活构造具有解耦特性的算子自定义小波基,使得系统多尺度刚阵具有沿对角线的强稀疏性,从而实现了该算法在每个尺度上独立、快速求解,系统方程的求解效率得到较大提高。建立多分辨Lagrange有限元空间和多尺度计算理论,提出基于稳定完备法的算子自定义小波弹性板单元构造方法及解耦条件。依据两尺度相对误差估计,提出自适应算子自定义小波有限元算法。数值算例证明,算子自定义小波弹性板单元具有求解精度与计算效率高等特点。     

复合材料弹塑性多尺度分析模型与算法

对材料非线性多尺度分析的计算模型与算法进行研究.在构建周期分布单胞分析算法的基础上,发展针对复合材料结构材料非线性多尺度分析的一般有限元方法.方法的特点是将所建立的单胞分析过程作为有限元分析的子程序嵌入到总程序系统当中,完成对应的高斯点应力计算,因而使所发展的方法具有实现方便的特点.给出数值计算结果,验证了方法与所发展的多尺度有限元分析程序的正确与有效性.     

FINITE ELEMENT RESOLUTION OF CONVECTION–DIFFUSION EQUATIONS WITH INTERIOR AND BOUNDARY LAYERS

The purpose of this paper is to present a new algorithm for the resolution of both interior and boundary layers present in the convection–diffusion equation in laminar regimes, based on the formulation of a family of polynomial– exponential elements. We have carried out an adaptation of the standard variational methods (finite element method and spectral element method), obtaining an algorithm which supplies non-oscillatory and accurate solutions. The algorithm consists of generating a coupled grid of polynomial standard elements and polynomial–exponential elements. The latter are able to represent the high gradients of the solution, while the standard elements represent the solution in the areas of smooth variation.     

基于CUDA的有限元矩阵并行装配算法研究

构建航天飞行器的结构有限元模型是准确模拟飞行仿真、完成飞行器在轨飞行阶段结构故障监测和诊断的基础。采用细长体飞行器简化梁模型，提出新的基于CUDA（Compute Unified Device Architecture）的有限元单元刚度矩阵生成和总刚度矩阵组装算法。依据梁单元矩阵的对称性，结合GPU硬件架构提出并行生成算法并进行改进。为有效减少装配时间，在装配过程中采用着色算法，提出了基于GPU（Graphics Processing Unit）共享内存的非零项组装策略，通过在不同计算平台下算例对比，验证了新算法的快速性。数值算例表明，本文算法的求解效率较高，针对一定计算规模内的模型可满足快速计算与诊断的实时性要求。     

基于改进位移模式的一维有限元超收敛算法

将多尺度方法的思想与超收敛计算的解析公式结合起来,提出了改进有限元位移模式的算法。利用超收敛计算的解析公式,将高阶有限元解的位移模式用常规有限元解的位移模式表示。用常规有限元解的位移模式与高阶有限元解的位移模式之和构造新的位移模式,采用积分形式推导了单元刚度矩阵。该算法在前处理和后处理两个阶段都使用超收敛计算公式,在常规试函数的基础上,增加了高阶试函数,使得单元内平衡方程的残差减少,从而达到提高精度的目标。对于线性单元,本文结点和单元的位移、导数都达到了h4阶的超收敛精度。     

铸件凝固温度场有限元分析中界面热阻的处理

提出一种处理铸件与铸型界面热阻问题的虚拟界面单元法，并给出了有限元计算公式。由于该公式不显含单元厚度(△l)，故该单元厚度△l可取任意值。当△l取为零时，使问题处理变得极为方便。针对某一具体金属型铝合金活塞的铸造凝固过程，按考虑和不考虑铸件与铸型间热阻影响两种方法作了有限元计算，通过与实测值相比较，本文提出的算法其计算精度远高于不考虑铸件与铸型间热阻影响的计算结果。另外．该方法使有限元建模方便、通用性强。