薄板弯曲单元的精化矩阵

本文提出一个改进薄板弯曲单元计算精度的一般性方法：用精化矩阵对原有的单元刚度阵进行修正。实践表明，这种方法可以有效地改变单元的刚柔特性，从而达到提高单元精度的目的。文中证明了，单刚中加进精化矩阵不会破坏原来单元的收敛性质。     

薄板几何非线性中的精化元方法及膜闭锁问题

基于放松单元间协调条件的大变形分原理和全局拉拉格朗日方法，推导出几何非线性精化三角形薄板单元，对几何刚度矩阵，通过引入特殊的单元位移函数，有效地消除了薄板弯曲问题中伴生的膜闭锁现象，数值结果表明该单元在几何非线性分析中既能消除膜闭锁又具有较高精度。     

薄板几何非线性中的精化元方法及膜闭锁问题

基于放松单元间协调条件的大变形变分原理和全局拉格朗日方法，推导了几何非线性精化三角形薄板单元。对几何刚度矩阵，通过引入特殊的单元位移函数，有效地消除了薄板弯曲问题中伴生的膜闭锁现象。数值结果表明该单元在几何非线性分析中既能消除膜闭锁又具有较高精度。     

8—21节点块体精化不协调元

本文基于文（２）提出的精化直接刚度法的变分原理，构造了８－２１节点元的精化不协调位移模式，并装入ＳＡＰ５程序中，该单元能通过分片试验数据结果表明：采用精化直接刚度法，可以有效地改装单元特性，使新单元在具原单元的简单、高效、可靠性好的同时，在精度，适应性等方面有大幅度的提高。     

精化杂交压电固体单元的研究

基于力、电耦合问题的三类交量广义交分原理，提出了广义杂交压电单元列式。为了进一步改进单元的性能和保证单元能够通过分片检验，通过引入非协调模式、放松电学方程约束条件和单元间的弱连续性条件，建立了新的、修正的广义交分原理，在此基础上成功地引入了应力、应交的正交化插值模式，从而建立了精化杂交压电单元法，它继承了常规精化杂交单元的全部优点。文中所推导的八节点精化杂交压电固体单元列式完全避免了矩阵求逆运算，较广义杂交压电单元和杂交应力压电单元均显著提高了计算效率。数值算例表明，与同类型其他单元相比，该单元明显具有更好的对歪斜网格的适应性。     

一种提高薄板稳定分析精度的方法

在薄板稳定分析中，九参数三角形薄板单元因其形状简单，使用方便，在实际工程中得到了广泛应用。本文基于参数调正的几何刚度矩阵，对九参数三角形薄板单元的一致刚度矩阵进行了修正，为提高计算精度提供了一种有效方法。     

大型结构分析程序系统单元库更新

阐述了大型结构分析程序系统单元库更新的概念及其意义，报告了装入ＡＰ－１程序系统单元库的精化四节点及八节点平面等参元的情况，给出了更新前、后单元通过考题库测试的结果及与其它程序系统结果的比较。数值结果表明通过更新单元来提高程序系统的分析精度与计算效率从而增强其竞争力是可行的。     

复合材料层合板精化高阶理论及其精化三角形板单元

提出一种新的精化高阶理论，该理论满足层间位移、应力连续条件，由此建立了三角形精化板单元。该单元满足单元间Ｃ１类弱连续条件，其收敛性得到保证，且具有简单、高效率的优点。     

大型结构分析程序系统单元库更新

阐述了大型结构分析系统单元库更新的概念及含意义，报告了装入ＡＰ－１程序系统单元库的精化四节点及八节点平面等参元的情况，给出了更新前、后单元通过考题库测试的结果及与其它程序系统结果的比较。数值结果表明通过更新单元来提高程序系统的分析精度与计算效率从而增强其竞争力是可行的。     

埃尔米特梁单元的分块对角与高阶质量矩阵

埃尔米特梁单元常用的集中质量矩阵,是由挠度自由度对应的一致质量矩阵元素通过行求和或节点积分构造。然而,数值结果表明该集中质量矩阵在求解包含自由端的梁振动问题时,会出现频率精度掉阶现象。本文首先从保障质量矩阵最优收敛性的数值积分精度出发,分别针对三次和五次梁单元,发展了质量矩阵的梯度增强节点积分方案。利用梯度增强节点积分方案,可以得到具有分块对角形式的单元质量矩阵,而其组装的整体质量矩阵除边界节点外仍然呈现对角形式。对于两种单元,其分块对角质量矩阵分别具有4阶最优精度和6阶次优精度。再者,将标准一致质量矩阵和具有同阶精度的梯度增强节点积分质量矩阵进行优化组合,建立了具有超收敛特性的高阶质量矩阵。最后,通过数值算例系统验证了三次和五次单元的分块对角与高阶质量矩阵的频率计算精度。     

The construction of large element in finite element method

In the usual finite element method, the order of the interpolation in an element is kept unchanged, and the accuracy is raised by subdividing the grid denser and denser. Alternatively, in the large element method, the grid is kept unchanged, and the terms of approximate series in the element are increased to raise the accuracy.In this paper, a method for constructing large elements is presented. When using this method, two sets of variables, one set defined inside the element, and the other defined on the boundary of the element, are adopted. Then, these two sets of variables are combined by the hybrid-penalty function method. This method can be applied to any elliptic equations in a domain with arbitrary shape and arbitrary complex boundary condition. It is proved with strict mathematical method in this paper, that in general cases, the accuracy of this method is much higher than that of the usual element and the large element method presented in [7]. Therefore, the degrees of freedom needed in this method are much fewer than those in the two methods if the same accuracy is preserved.     

A simple method for simulating viscoelastic fluid flows in slit channel

A low-cost semi-analysis finite element technique, named the finite piece method (FPM) is presented in this article. It aims to solve three-dimensional (3D) viscoelastic slit flows. The viscoelastic stress of the fluid is modelled using an K-BKZ integral constitutive equation of the Wagner type. Picard iteration is used to solve non-linear equations. The FPM is tested on flow problems in both planar and contraction channels. The accuracy of the method is assessed by comparing flow distributions and pressure with results obtained by 3D finite element method (FEM). It shows that the solution accuracy is excellent and a substantial amount of computing time and memory requirement can be saved.     

Geometrically non-linear generalized hybrid element and refined element method of non-conforming modes

A generalized hybrid method of non-conforming modes based on a non-linear generalized variational principles with relaxed interelement continuity requirements, is developed, and the plane quadrilateral geometrically non-linear element is presented, furthermore, non-linear refined element method is devised by orthogonal approach. It is shown that the refined element can improve the computational accuracy for non-conforming modes. The project supported by the National Natural Science Foundation of China     

A 3D pyramid spline element

In this paper,a 13-node pyramid spline element is derived by using the tetrahedron volume coordinates and the B-net method,which achieves the second order completeness in Cartesian coordinates.Some appropriate examples were employed to evaluate the performance of the proposed element.The numerical results show that the spline element has much better performance compared with the isoparametric serendipity element Q20 and its degenerate pyramid element P13 especially when mesh is distorted,and it is comparable to the Lagrange element Q27.It has been demonstrated that the spline finite element method is an efficient tool for developing high accuracy elements.     

A conformal decomposition finite element method for modeling stationary fluid interface problems

A method is developed for modeling fluid transport in domains that do not conform to the finite element mesh. One or more level set functions are used to describe the fluid domain. A background, non‐conformal mesh is decomposed into elements that conform to the level set interfaces. Enrichment takes place by adding nodes that lie on the interfaces. Unlike other enriched finite element methods, the proposed technique requires no changes to the underlying element assembly, element interpolation, or element quadrature. The complexity is entirely contained within the element decomposition routines. It is argued that the accuracy of the method is no less than that for eXtended Finite Element Methods (XFEM) with Heaviside enrichment. The accuracy is demonstrated using multiple numerical tests. In all cases, optimal rates of convergence are obtained for both volume and surface quantities. Jacobi preconditioning is shown to remove the ill‐conditioning that may result from the nearly degenerate conformal elements. Copyright © 2009 John Wiley & Sons, Ltd.     

网架结构拟夹层板法的有限元验证

用拟夹层板法和有限元法对网架结构进行分析,对三类屋面网架（正放四角锥网架、两向正交正放网架和正放抽空四角锥网架）进行了均布荷载、局部荷载（半跨均布荷载）作用下的静力分析以及固有振动分析,对三类竖向承重网架墙体进行了稳定性分析。通过与有限元法分析结果的对比,表明了拟夹层板法作为一种简化的计算方法,其精度是比较高的,绝大多数的误差都小于5％,是可以直接用于工程结构设计的一种有效方法。此外,拟夹层板法还可作为一种在宏观上检验有限元建模正确与否的实用方法。     

The perturbation finite element method for solving the plane problem in consideration of linear creep

In this paper, a method (PFMC) for solving plane problem of linear creep is presented by using perturbation finite element. It can be used in plane problem in consideration of creep, such as reinforced concrete beam, presiressed concrete beam, reinforced concrete cylinder and reinforced concrete tunnel in elastic or visco-elastic medium, as well as underground building and so on.In the presented method, the assumption made in the general increment method that variables remain constant in a divided time interval is not taken. The accuracy is improved and the length of time step becomes larger. The computer storage can be reduced and the calculating efficiency can be increased.Perturbation finite element formulae for four-node quadrilateral isoparametric element including reinforcement are established and five numerical examples are given. As contrasted with the analytical solution, the accuracy is satisfactory.     

COARSE-MESH-ACCURACY IMPROVEMENT OF BILINEAR Q4-PLANE ELEMENT BY THE COMBINED HYBRID FINITE ELEMENT METHOD

IntroductionByadding‘nodeless’incompatiblebubblemodesandpreservinggeometriccharacteristicofthevariationalprincipleinmechanics,thecombinedhybridmethod[1～4 ]remarkablelyenhancedcoarse_mesh_accuracyofconventionalquadrilateralelementsoflowerorder.ThequadrilateralplaneelasticelementCH(0_1 )proposedinRef.[3 ]isasuccessfulexample.Followingthegeometricpointofviewinmechanicscombinedwithmathematicalanalysis,anovelexpressionofthecombinedhybridvariationalprincipleisintroducedtoclarifyitsintrinsicmecha…     

各向异性弹塑性中厚度板壳的有限元分析

本文在文献［２，３］的基础上，提出了一个解各向异性弹塑性中厚度板壳问题的有限元方法。考虑材料各向异性的特点，采用了Ｈｉｌｌ推广的Ｈｕｂｅｒ－Ｍｉｓｅｓ屈服准则；借用Ｏｗｅｎ的剪切修正系数，正确计及了叠层复合材料壳体的横向剪切效应；为了避免“自锁”现象，文中采用了９节点的Ｈｅｔｅｒｏｓｉｓ二次壳单元；特别是本文利用插值外推的思想，提出了一个带预测的弧长增量控制法，显著提高了确定变形路径的计算效率。几个数值算例表明本文给出的有限元方法对于各向异性中厚度板壳的弹塑性分析有较好的精度，尤其是对具有复杂变形路径的结构计算，收敛速度提高更快。     

Velocity–pressure integrated versus penalty finite element methods for high-Reynolds-number flows

Velocity–pressure integrated and consistent penalty finite element computations of high-Reynolds-number laminar flows are presented. In both methods the pressure has been interpolated using linear shape functions for a triangular element which is contained inside the biquadratic flow element. It has been shown previously that the pressure interpolation method, when used in conjunction with the velocity-pressure integrated method, yields accurate computational results for high-Reynolds-number flows. It is shown in this paper that use of the same pressure interpolation method in the consistent penalty finite element method yields computational results which are comparable to those of the velocity–pressure integrated method for both the velocity and the pressure fields. Accuracy of the two finite element methods has been demonstrated by comparing the computational results with available experimental data and/or fine grid finite difference computational results. Advantages and disadvantages of the two finite element methods are discussed on the basis of accuracy and convergence nature. Example problems considered include a lid-driven cavity flow of Reynolds number 10 000, a laminar backward-facing step flow and a laminar flow through a nest of cylinders.