定常Satokes方程的窄边四边形有限元法

本文提出了一个试探函数不满足divν＝0的定常Stokes方程的窄边四边形有限元法，利用窄边四边形等参有限元插值定理和Falk的思想，得到了H^1(Ω）模误差的最优收敛估计阶O（h）。     

Nonlinear finite element reliability analysis of Concrete-Faced Rockfill (CFR) dams under static effects

The response of concrete slab on Concrete-Faced Rockfill (CFR) dams is very important. This study investigates the reliability of the concrete slab on a CFR dam by the improved Rackwitz–Fiessler method under static loads. For this purpose, ANSYS finite element analysis software and FERUM reliability analysis program are combined with direct coupled method and response surface method. Reliability index and probability of failure of the concrete are computed in the all critical points of the concrete slab by dam height. This study is also expanded for the reliability of CFR dams including different concrete slab thickness. In addition to the linear behavior, geometrically and materially non-linear responses of the dam are considered in the finite element analysis which is performed with reliability analysis. The Drucker–Prager method and the multi linear kinematic hardening method are, respectively, used for concrete slab and for rockfill and foundation rock. Finite element model used in the analyses includes dam–reservoir–foundation interaction. Reservoir water is modeled by the Lagrangian approach. Welded and friction contact based on the Coulomb’s friction law are considered in the joints of the dam. One-dimensional two noded contact elements are used to define friction. The self-weight of the dam and the hydrostatic pressure of the reservoir water are considered in the numerical solutions. According to this study, hydrostatic pressure, nonlinear response of the rockfill and the decrease in the concrete slab thickness reduce the reliability of the concrete slab of the CFR dam. Besides, the CFR dam models including friction are safer than the models including welded contact in the joints.     

Geometrically nonlinear analysis with a 4-node membrane element formulated by the quadrilateral area coordinate method

Recently, a 4-node quadrilateral membrane element AGQ6-I, has been successfully developed for analysis of linear plane problems. Since this model is formulated by the quadrilateral area coordinate method (QACM), a new natural coordinate system for developing quadrilateral finite element models, it is much less sensitive to mesh distortion than other 4-node isoparametric elements and free of various locking problems that arise from irregular mesh geometries. In order to extend these advantages of QACM to nonlinear applications, the total Lagrangian (TL) formulations of element AGQ6-I was established in this paper, which is also the first time that a plane QACM element being applied in the implicit geometrically nonlinear analysis. Numerical examples of geometrically nonlinear analysis show that the presented formulations can prevent loss of accuracy in severely distorted meshes, and therefore, are superior to those of other 4-node isoparametric elements. The efficiency of QACM for developing simple, effective and reliable serendipity plane membrane elements in geometrically nonlinear analysis is demonstrated clearly.     

Two different approaches for the treatment of boundary singularities

The finite element treatment of boundary singularities in elliptic problems has demanded that special techniques be developed. Many of these use some form of singular element in the neighborhood of a singular point. For a homogenous second-order problem, defined in a domain with a polygonal boundary, we studied three cases with different singularity orders. Some results about the accuracy of the solution are presented. Numerical results have been obtained using Akin singular triangular elements with three and six nodes and quadrilateral elements having four and eight nodes. The behavior of the gradient in each of these elements is also analyzed. However, these elements are not completely satisfactory, so that an alternative technique using curved isoparametric elements is given here. The results obtained with the two methods are compared. Conclusions about numerical accuracy of each method, the order of integration and the simplicity of application are made.     

Shape optimal design of arch dams including dam-water–foundation rock interaction using a grading strategy and approximation concepts

Optimal design of arch dams including dam-water–foundation rock interaction is achieved using the soft computing techniques. For this, linear dynamic behavior of arch dam-water–foundation rock system subjected to earthquake ground motion is simulated using the finite element method at first and then, to reduce the computational cost of optimization process, a wavelet back propagation neural network (WBPNN) is designed to predict the arch dam response instead of directly evaluating it by a time-consuming finite-element analysis (FEA). In order to enhance the performance generality of the neural network, a dam grading technique (DGT) is also introduced. To assess the computational efficiency of the proposed methodology for arch dam optimization, an actual arch dam is considered. The optimization is implemented via the simultaneous perturbation stochastic approximation (SPSA) algorithm for the various conditions of the interaction problem. Numerical results show the merits of the suggested techniques for arch dam optimization. It is also found that considering the dam-water–foundation rock interaction has an important role for safely designing an arch dam.     

<Emphasis Type="Italic">L</Emphasis><Superscript>2</Superscript> error estimate of the finite volume element methods on quadrilateral meshes

In this paper the optimal L ² error estimates of the finite volume element methods (FVEM) for Poisson equation are discussed on quadrilateral meshes. The trial function space is taken as isoparametric bilinear finite element space on quadrilateral partition, and the test function space is defined as piecewise constant space on dual partition. Under the assumption that all elements on quadrilateral meshes are O(h ²) quasi-parallel quadrilateral elements, we prove convergence rate to be O(h ²) in L ² norm.     

线性蠕变平面问题中的摄动有限元方法

本文提出一种利用摄动有限元求解线性蠕变问题的方法(PFMC).它可用于平面蠕变问题,诸如钢筋混凝土梁,预应力混凝土梁,钢筋混凝土圆筒或位于弹性或粘弹性介质中的钢筋混凝土隧洞及地下建筑等.本方法不采用一般增量法中关于在一个时段内各物理量保持不变的假设,提高了分析精度,加大了计算步长,减少了机器存储,提高了运算效率.本文构造了包含钢筋在内的四节点四边形等参单元的有限元摄动格式.并给出五个算例,与解析解相比,有令人满意的精度.     

The interpolation theorem for narrow quadrilateral isoparametric finite elements

Summary. The interpolation theorem for convex quadrilateral isoparametric finite elements is proved in the case when the condition is not satisfied, where is the diameter of the element and is the radius of an inscribed circle in . The interpolation error is in the -norm and in the -norm provided that the interpolated function belongs to . In the case when the long sides of the quadrilateral are parallel the constants appearing in the estimates are evaluated. Received September 1993 / Revised version received March 6, 1995     

Damage analysis of arch concrete dams subjected to underwater explosion

The nonlinear behavior of arch concrete dams subjected to underwater explosions was investigated. Pressure waves were spherically modeled using a 3-D finite element method. The nonlinear behavior of concrete was modeled using the concrete damage plasticity model. In addition, absolute wave formulation was used to develop a fully-coupled model incorporating the propagation of the shock wave and water–dam interaction. Analysis of an arch concrete dam subject to underwater explosion was performed upon model verification. The dynamic response of the dam subject to the underwater explosion was studied for different sizes of explosions at different depths. The results showed that the closer the point to the explosive source, the sooner the impression and damage was caused by the blast. It was clear that the middle part of the dam facing the explosion was the first location to incur damage and then cracks extended to the downstream face and abutments. Finally, cracks extended to the bottom of the dam. The results of the dynamic analysis and dissipated energy by the loss of elements’ resistance indicated that the time slice of 150 ms was appropriate to analyze arch concrete dam structure subject to an underwater explosion. In addition, for a given amount of explosives, an increase in the depth of explosion corresponded to an increase in the maximum displacement of the upstream face of the dam. The results also confirmed that the damage distribution to the arch concrete dam subject to an underwater explosion depended upon the amount and depth of explosive materials.     

非矩形区域上半导体问题的交替方向有限元方法

Transient behavior of semiconductor with heat-conduction on nonrectangular is studied using isoparametric elements and an approximation to the Jacobian of the isoparametric map.Concentration and heat-conduction equations are solved by alternating-direction methods and electric potential equation is approximated by finite element method.Optimal order error estimates in L^2 are demonstrated using the theory and technique of a prior estimate of differential equation.     

Transient dynamics of composite sandwich plates using 4-, 8-, 9-noded isoparametric quadrilateral elements

A finite element method based on Mindlin's theory is employed in the prediction of the dynamic transient response of multilayered composite sandwich plates. Numerical convergence and stability of 4-noded linear, 8-noded serendipity, and 9-noded Lagrangian elements are established using an explicit time integration technique. A special mass matrix diagonalization scheme is adopted which conserves the total mass of the element and includes the effects due to rotary inertia terms. The parametric effects of the time step, finite element mesh, lamination scheme, and orthotropy on the transient response are investigated. Numerical results for deflections and stresses are presented under various boundary conditions and loadings. The results presented herein should be of interest to composite-structure designers, experimentalists, and numerical analysists in verifying their results.     

Thermal analysis of hot strip rolling using finite element and upper bound methods

Thermal analysis of hot rolling process has been studied in this work. A finite element method has been coupled with an upper bound solution assuming, triangular velocity field, to predict temperature field during hot strip rolling operation. To do so, an Upwind Petrov–Galerkin scheme together with isoparametric quadrilateral elements has been employed to solve the steady-state heat transfer equation. A comparison has been made between the published and the model predictions and a good agreement was observed showing the accuracy of the proposed model.     

Analysis of reinforced concrete structures subjected to dynamic loads with a viscoplastic Drucker–Prager model

The behavior of reinforced concrete structures subjected to dynamic loads is analyzed. The concrete material is modelled by an elasto-viscoplastic law, whose inviscid counterpart is the Drucker–Prager model. A viscous regularization is introduced in order to avoid the mesh dependency effects that usually appear when strain softening occurs. The model is implemented in a general finite element computer code for fast transient analysis of fluid-structure systems, based on an explicit central difference scheme. The model is activated to both continuum elements and layered shell elements. So, realistic numerical analyses of complex 3-D engineering problems are simple and efficient. Three examples, two of which are modelled with layered shell elements, are presented below.     

Implementing and using high-order finite element methods

We present theoretical analyses of and detailed timings for two programs which use high-order finite element methods to solve the Navier- Strokes equations in two and three dimensions. The analyses show that algorithms popular in low-order finite element implementations are not always appropriate for high-order methods. The timings show that with the proper algorithms high-order finite element methods are viable for solving the Navier-Stokes equations. We show that it is more efficient, both in time and storage, not to precompute element matrices as the degree of approximating functions increases. We also study the cost of assembling the stiffness matrix versus directly evaluating bilinear forms in two and three dimensions. We show that it is more efficient not to assemble the full stiffness matrix for high-order methods in some cases. We consider the computational issues with regard to both Euclidean and isoparametric elements. We show that isoparametric elements may be used with higher-order elements without increasing the order of computational complexity.     

广义协调等参元

本文根据广义协调条件,对平面应力四边形单元提出一个广义协调等参元GC-Q6.单元GC-Q6是对Wilson非协调等参元Q6的一个改进方案:单元Q6只对平行四边形网格情况能通过分片检验,对一般四边形情况却不能通过,而本文的单元GC-Q6则对一般四边形情况也能通过.当单元为平行四边形时,GC-Q6单元即退化为QS单元.算例表明,广义协调等参元的应力精度高于文献中已有的单元,对不规则网格均能保持良好的性态.     

A flexible approach to the simulation of hybrid multibody systems

The present works deals with the incorporation of both flexible beam and shell structures into the realm of flexible multibody dynamics. Geometrically exact beam formulations based on classical Simo-Reissner kinematics are suitable for modelling beam-type flexible components in the context of finite-deformation multibody dynamics. So geometrically exact shell formulations are based on Reissner-Mindlin kinematics. In [2], a flexible framework for dealing with flexible structural elements in a multibody context is described. A specific isoparametric finite element discretization of a shell formulation leads to semi-discrete equations of motions assuming the form of differential-algebraic equations (DAEs). A compatible isoparametric formulation of beams has already been developed in [1]. The uniform DAE framework makes possible the incorporation of alternative finite element formulations. In addition to that, various time-stepping schemes such as energy-momentum methods or variational integrators can be applied. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Nonconforming finite element methods on quadrilateral meshes

It is well known that it is comparatively difcult to design nonconforming fnite elements on quadrilateral meshes by using Gauss-Legendre points on each edge of triangulations.One reason lies in that these degrees of freedom associated with these Gauss-Legendre points are not all linearly independent for usual expected polynomial spaces,which explains why only several lower order nonconforming quadrilateral fnite elements can be found in literature.The present paper proposes two families of nonconforming fnite elements of any odd order and one family of nonconforming fnite elements of any even order on quadrilateral meshes.Degrees of freedom are given for these elements,which are proved to be well-defned for their corresponding shape function spaces in a unifying way.These elements generalize three lower order nonconforming fnite elements on quadrilaterals to any order.In addition,these nonconforming fnite element spaces are shown to be full spaces which is somehow not discussed for nonconforming fnite elements in literature before.     

Seismic collapse capacity spectra for simple structures vulnerable to dynamic instabilities

In this study collapse capacity spectra based on various definitions of the seismic intensity are set in contrast and evaluated. The presented collapse capacities for highly inelastic non-deteriorating single-degree-of-freedom (SDOF) systems, which are vulnerable to the destabilizing effect of gravity loads, are derived for a near-fault set of ground motions with distinct pulse characteristics. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

惩罚函数法在三维等参杂交元的应用

提出在等参杂交元中用惩罚函数法引入平衡约束条件,具体讨论了惩罚函数法在三维等参杂交元中的运用,并提出采用分项罚数的方法,建立最佳的罚平衡杂交元模型．罚平衡法可以在不增加自由度的前提下,有效地扼制寄生应力．数值实验表明,新建立的单元,可以有效地抑制单元畸变对计算精度的影响,从而大幅度提高畸变网格下的计算精度,方法带有普遍性．