Nonlinear quasi-conforming finite element method

The nonlinear quasi-conforming FEM is presented based on the basic concept of the quasiconforming finite element. First, the incremental principle of stationary potential energy is discussed. Then, the formulation process of the nonlinear quasi-conforming FEM is given. Lastly, two computational examples of shells are given.     

The displacement function of quasi-conforming element and its node error

Based on the strain formulation of the quasi-conforming finite element, displacement functions are constructed which have definite physical meaning, and a conclusion can be obtained that the coefficients of the constant and the linear strain are uniquely determined, and the quasi-conforming finite element method is convergent to constant strain. There are different methods for constructing the rigid displacement items, and different methods correspond to different order node errors, and this is different from ordinary displacement method finite element. Contributed by TANG Li-min Biography: HE Dong-sheng (1964-)     

Boundary-type finite element method for wave propagation analysis

The boundary-type finite element method has been investigated and applied to the Helmholz and mild-slope equations. Four types of interpolation function are examined based on trigonometric function series. Three-node triangular, four-node quadrilateral, six-node triangular and eight-node quadrilateral elements are tested; these are all non-conforming elements. Three types of numerical example show that the three-node triangular and four-node quadrilateral elements are useful for practical analysis.     

随机有限元分析的二级区域分解并行求解算法

采用蒙特卡罗模拟(MCS)和加权积分法对二维问题进行随机有限元分析。尽管MCS方法对任何有确定解的问题都具有求解精度高的优点,但由于求解所需的计算量巨大使其应用受到限制。利用并行求解技术可有效地处理这种密集型计算问题。基于有限元分裂对接法(FETI)的并行特性并利用预处理共轭梯度法(PCG)的求解高效性,结合整体子区域实现(GSI-PCG)和FETI法,提出二级求解算法,并在工作站机群上实现了数值算例。算例计算结果表明本文GSI(PCG)-FETI算法具有较高的并行加速比和并行效率,具有良好的性能,可有效地进行二维问题的随机有限元分析。     

Exact geometry based quasi-conforming analysis for Euler-Bernoulli beam

The problem of quick analysis using exact geometry data was proposed by Hughes et al. and the isogeometric analysis framework was introduced as a solution. In this letter, the exact geometry concept is combined into the quasi-conforming framework and a novel method, i.e., the exact geometry based quasi-conforming analysis is proposed. In present method the geometry is exactly described by non-uniform rational B-spline bases, while the solution space by traditional polynomial bases. Present method combines the merits of both isogeometric analysis and quasi-conforming finite element method. In this letter Euler-Bernoulli beam problem is solved as an example and the results show that the present method is effective and promising.     

A moving finite element method for the population balance equation

A moving finite element algorithm has been compared against the upwind-differencing and Smolarkiewicz methods for the population balance equation of multicomponent particle growth processes. Analytical solutions and an error function have been used to test the numerical methods. The moving finite elements technique is much more accurate than other methods for a wide range of parameters. Since this method uses moving grids, it is able to model very narrow particle size distributions. It is also shown that the method can be extended to solve condensational growth problems which include particle curvature and non-continuum mass transfer effects.     

A boundary-type finite element model for water surface wave problems

A new combinative method of boundary-type finite elements and boundary solutions is presented to study wave diffraction-refraction and harbour oscillation problems. The numerical model is based on the mild-slope equation. The key feature of this method is that the discretized matrix equation can be formulated only by the calculation of a line integral, since the interpolation equation which satisfies the governing equation in each element is used. The numerical solutions are compared with existing analytical, experimental, observed and other numerical results. The present method is shown to be an effective and accurate method for water surface wave problems.     

Complex wavenumber Fourier analysis of the B-spline based finite element method

We present the results of one-dimensional complex wavenumber Fourier analysis of the B-spline variant of Finite Element Method (FEM). Generally, numerical results of elastic wave propagation in solids obtained by FEM are polluted by dispersion and attenuation. It was shown for the higher-order B-spline based FEM, that the optical modes did not occur in the case of infinite domains, unlike the higher-order Lagrangian and Hermitian finite elements, and also the dispersion errors are smaller. The paper’s main focus is on the wave propagation through B-spline multi-patch/segment discretization with the C⁰$C^0$ connection of B-spline segments and, chiefly, to the determining of dispersion and attenuation dependences. The numerical approach employed leads to substantial minimization of dispersion errors. Furthermore, the errors decrease in line with the increasing order of the B-spline elements/segments, with the local refinement, and also by the particular choice of the positions of control points through the optimizing procedure.     

薄膜横向振动的有限元分析

构造了一种3节点三角形膜单元,以适用于平面薄膜横向振动的有限元分析.在给出单元形函数的基础上,根据最小势能原理建立了薄膜自由振动方程,并推导了单元刚度矩阵和单元质量矩阵.研究结果表明,单元刚度矩阵和单元质量矩阵形式简单,且自由度少;通过两个典型算例,证明3节点三角形膜单元的计算结果非常接近理论解,同时可以达到很高的精度...     

Three dimensional finite element analysis for Rayleight-Benard convection in rectangular enclosures

The pattern of Rayleigh-Benard convection of air in a rectangular box heated-from-below is studied by numerically solving the three-dimensional time-dependent Navier-Stokes equations under the Boussinesq approximation. Slightly supercritical Rayleigh number was adopted to track the evolutions of flow structure as a function of enclosure's aspect ratio (A=L/H). The flow will asymptotically evolve to different patterns, among which, two possible types of flow pattern are found. One consists of the pair of straight vortex rolls and the other appears as closed vortex rings. The transition between the flow patterns indicates that there exists a flow bifurcation with the variation of container's aspect ratio. In addition, both steady and oscillatory flows have been observed, corresponding to the pair of straight vortex rolls and the vortex ring, respectively. The complexity of flow structure tends to increase with the increasing aspect ratio of the rectangular enclosure. The project supported by the National Natural Science Foundation of China (19889210), the National Distinguished Young Fund (10125210), the Hundred Talents Program of CAS, and the Training Program for the Trans-Century Outstanding Young of MOE     

A simple and efficient triangular finite element for plate bending

A simple and efficient triangular finite element is introduced for plate bending application. The element is a three-node triangular one with three basic degrees of freedom per node and two internal rotation degrees of freedom, using selective reduced integration. Numerical examples indicate that, despite its simplicity, the element is not only competitively accurate, but also useful as a thick/thin triangular plate bending element. It is also pointed out that this element using selective reduced integration is in fact a mixed element.     

A finite amplitude analysis of tunnel deformation by the finite element method with highly viscous fluid

A finite element method for highly viscous fluid is used to calculate the velocity and stress fields in the surrounding soft rock of a tunnel. In order to fit the calculated values with the measured displacement of tunnel wall, we inverted the boundary forces and the mechanical parameters of the surrounding rocks.     

A smeared / layered reinforced concrete model for finite element analysis

A new reinforced concrete model, in which the reinforcement steel is assumed as smeared / layered in concrete, is established and installed into a currently used finite element code for nonlinear analysis. It performs the nonlinear behaviors of both concrete and the reinforcement steel. The results of examples are in good agreement with the experimental data.     

Rigid finite element and limit analysis

According to the lower bound theorem of limit analysis the Rigid Finite Element Method (RFEM) is applied to structural limit analysis and the linear programmings for limit analysis are deduced in this paper. Moreover, the Thermo-Parameter Method (TPM) and Parametric Variational principles (PVP) are used to reduce the computational effort while maintaining the accuracy of solutions. A better solution is also obtained in this paper. The project supported by National Natural Science Foundation of China     

压电层合梁静动力学分析的高精度梁单元

给出了一个对复合材料压电层合梁进行数值分析的高精度压电层合梁单元。基于Shi三阶剪切变形板理论的位移场和Layer-wise理论的电势场,用力-电耦合的变分原理及Hamilton原理推导了压电层合梁单元列式。采用拟协调元方法推导了一个可显式给出单元刚度矩阵的两节点压电层合梁单元,并应用于压电层合梁的力-电耦合弯曲和自由振动分析。计算结果表明,该梁单元给出的梁挠度和固有频率与解析解吻合良好,并优于其它梁单元的计算结果,说明了本文所给压电层合梁单元的可靠性和准确性。研究结果可为力-电耦合作用下压电层合梁的力学分析提供一个简单、精确且高效的压电层合梁单元。     

A deforming finite element mesh for use in moving one-dimensional boundary wave problems

The finite element method is developed to solve the problem of wave run-up on a mild, plane slope. A novel approach to implementing a deforming mesh of one-dimensional, three-node, isoparametric elements is described and demonstrated. The discrete time interval (DTI), arbitrary Lagrangian–Eulerian (ALE) and space–time element (STE) methods are used to solve the unsteady one-dimensional shallow water wave equations. The boundary condition required is simply the seaward water surface elevation, and although the method has only been tested for monochromatic waves, it should be equally valid for any sea state which can be described as a water surface elevation as a function of time. All three solution methods are shown to given good results. Time histories of the terms of the governing equations are calculated and used to demonstrate how the ALE and STE methods account for mesh deformation. The model could be extended to two dimensions, which would have practical application to the run-up of obliquely incident waves. © 1997 John Wiley & Sons, Ltd.     

An explicit finite volume element method for solving characteristic level set equation on triangular grids

Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.     

Application of scaled boundary finite element method in static and dynamic fracture problems

The scaled boundary finite element method (SBFEM) is a recently developed numerical method combining advantages of both finite element methods (FEM) and boundary element methods (BEM) and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors (SIFs) directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.The project supported by the National Natural Science Foundation of China (50579081) and the Australian Research Council (DP0452681)The English text was polished by Keren Wang.     

Some applications of semi-analytical finite element method in static and dynamic analysis of structures

In this paper a kind of semi-analytical finite element method based on the transfer matrix analysis of shells of revolution is briefly formulated. Some recent investigations on its application are presented: (1) a reanalysis algorithm for improving the accuracy of free vibration analysis; (2) a kind of semi-analytical ring element for the stress analysis of a curved pipe of the slender torus shell type.