ASYMPTOTIC SOLUTIONS FOR PREDICTED NATURAL FREQUENCIES OF TWO-DIMENSIONAL ELASTIC SOLID VIBRATION PROBLEMS IN FINITE ELEMENT ANALYSIS

In order to assess the discretization error of a finite element solution, asymptotic solutions for predicted natural frequencies of two-dimensional elastic solid vibration problems in the finite element analysis are presented in this paper. Since the asymptotic solution is more accurate than the original finite element solution, it can be viewed as an alternative solution against which the original finite element solution can be compared. Consequently, the discretization error of the finite element solution can be evaluated. Due to the existence of two kinds of two-dimensional problems in engineering practice, both the plane stress problem and the plane strain problem have been considered and the corresponding asymptotic formulae for predicted natural frequencies of two-dimensional solids by the finite element method have been derived from the fact that a discretized finite element system approaches a continuous one if the finite element size approaches zero. It has been demonstrated, from the related numerical results of three examples, that the present asymptotic solution, which can be obtained by simply using the corresponding formula without any further finite element calculation, is indeed more accurate than the original finite element solution so that it can be considered as a kind of corrected solution for the discretization error estimation of a finite element solution.     

The domain composition method applied to Poisson's equation in two dimensions

Domain composition, a recently described method for formulating continuum field problems, removes certain restrictions on the construction of finite element models such that it is possible to solve a finite element problem without using a global compatible mesh. The domain composition method couples or otherwise constrains meshes in local regions to obtain a solution equivalent to that produced by conventional finite element methods. In particular, the domain composition method enables finite element models to be formulated with overlapping elements. Several advantages come from this, including an ability to automatically generate a finite element model from a solid geometric model, an ability to use a variety of element types in a single finite element model and an ability to exactly match element boundaries to the local geometry. This paper shows in detail a finite element formulation of Poisson's equation using domain composition and presents certain key algorithms that incorporate the domain composition method into well-established finite element procedures.     

Stress recovery and error estimation for the scaled boundary finite‐element method

The scaled boundary finite‐element method is a novel semi‐analytical technique, combining the advantages of the finite element and the boundary element methods with unique properties of its own. This paper develops a stress recovery procedure based on a modal interpretation of the scaled boundary finite‐element method solution process, using the superconvergent patch recovery technique. The recovered stresses are superconvergent, and are used to calculate a recovery‐type error estimator. A key feature of the procedure is the compatibility of the error estimator with the standard recovery‐type finite element estimator, allowing the scaled boundary finite‐element method to be compared directly with the finite element method for the first time. A plane strain problem for which an exact solution is available is presented, both to establish the accuracy of the proposed procedures, and to demonstrate the effectiveness of the scaled boundary finite‐element method. The scaled boundary finite‐element estimator is shown to predict the true error more closely than the equivalent finite element error estimator. Unlike their finite element counterparts, the stress recovery and error estimation techniques work well with unbounded domains and stress singularities. Copyright © 2002 John Wiley & Sons, Ltd.     

An enhanced layerwise finite element using a robust solid-shell formulation approach

The application of layerwise theories to correctly model the displacement field of sandwich structures or laminates with high modulus ratios, usually employs plate or facet shell finite element formulations to compute the element stiffness and mass matrices for each layer. In this work, a different approach is proposed, using a high performance hexahedral finite element to represent the individual layer mass and stiffness. This 8-node hexahedral finite element is formulated based on the application of the enhanced assumed strain method (EAS) to resolve several locking pathologies coming from the high aspect ratios of the finite element and the usual incompressibility condition of the core materials. The solid-shell finite element formulation is introduced in the layerwise theory through the definition of a projection operator, which is based on the finite element variables transformation matrix. The new finite element is tested and the implemented numerical remedies are verified. The results for a soft core sandwich plate are hereby presented to demonstrate the proposed finite element applicability and robustness.     

含裂纹压电材料的Cell-Based光滑扩展有限元法

为精确而有效地求解机电耦合作用下含裂纹压电材料的断裂参数,首先,通过将复势函数法、扩展有限元法和光滑梯度技术引入到含裂纹压电材料的断裂机理问题中,提出了含裂纹压电材料的Cell-Based光滑扩展有限元法;然后,对含中心裂纹的压电材料强度因子进行了模拟,并将模拟结果与扩展有限元法和有限元法的计算结果进行了对比。数值算例结果表明:Cell-Based光滑扩展有限元法兼具扩展有限元法和光滑有限元法的特点,不仅单元网格与裂纹面相互独立,且裂尖处单元不需精密划分,与此同时,Cell-Based光滑扩展有限元法还具有形函数简单且不需求导、对网格质量要求低且求解精度高等优点。所得结论表明Cell-Based光滑扩展有限元法是压电材料断裂分析的有效数值方法。      

Comparative study of different nonconserving time integrators for wave propagation in hyperelastic waveguides

This paper addresses the formulation and numerical efficiency of various numerical models of different nonconserving time integrators for studying wave propagation in nonlinear hyperelastic waveguides. The study includes different nonlinear finite element formulations based on standard Galerkin finite element model, time domain spectral finite element model, Taylor–Galerkin finite element model, generalized Galerkin finite element model and frequency domain spectral finite element model. A comparative study on the computational efficiency of these different models is made using a hyperelastic rod model, and the optimal computational scheme is identified. The identified scheme is then used to study the propagation of transverse and longitudinal waves in a Timoshenko beam with Murnaghan material nonlinearity.     

New basis functions and computational procedures for p-version finite element analysis

New basis functions and solution procedures for p-version finite element analysis are described. They are used in a highly efficient p-version finite element solver for linear elastostatics and dynamics, which has been used in an industrial environment for over two years. Using two sample applications it is shown that, using the techniques proposed here, p-version finite element analysis can have a substantially lower computational cost, for given accuracy, than standard finite element methods. This makes the industrial applicability of p-version finite element analysis much wider than is commonly believed.     

Moving particle finite element method with global smoothness

We describe a new version of the moving particle finite element method (MPFEM) that provides solutions within a C⁰ finite element framework. The finite elements determine the weighting for the moving partition of unity. A concept of ‘General Shape Function’ is proposed which extends regular finite element shape functions to a larger domain. These are combined with Shepard functions to obtain a smooth approximation. The Moving Particle Finite Element Method combines desirable features of finite element and meshfree methods. The proposed approach, in fact, can be interpreted as a ‘moving partition of unity finite element method’ or ‘moving kernel finite element method’. This method possesses the robustness and efficiency of the C⁰ finite element method while providing at least C¹ continuity. Copyright © 2004 John Wiley & Sons, Ltd.     

A solid-shell layerwise finite element for non-linear geometric and material analysis

The application of layerwise theories to correctly model the displacement field of sandwich structures or laminates with high modulus ratios usually employs plate or facet-shell finite element formulations to compute the element stiffness and mass matrices for each layer. In this work an alternative approach is proposed, using a high performance hexahedral finite element to represent the individual layer mass and stiffness. This eight-node hexahedral finite element is formulated based on the application of the enhanced assumed strain method (EAS) to solve several locking pathologies coming from the high aspect ratio of the finite element and the usual incompressibility condition of the core materials. The solid-shell finite element formulation is introduced in the layerwise theory through the definition of a projection operator, based on the finite element variables transformation matrix. The non-linear geometric and material capabilities are introduced into the finite element formulation, allowing for the representation of large displacements, large deformation and material non-linear behaviors. The developed formulation is numerically tested and benchmarked, being validated by using published experimental results obtained from sandwich specimens.     

Fracture analysis of interfacial crack by global-local finite element

An eigenfunction expansion is used to formulate the global element on the crack tip. The global-local finite element method employs both conventional finite element and classical Rayleigh-Ritz kinematic approach. The hybrid Ritz method not only preserves the finite element modelling capability but adds the advantage of using prior information regarding the anticipate behaviour of the particular problem. Thus, it is able to achieve better accuracy with fewer elements in comparison with conventional finite element. Several examples relative to crack problems are presented to demonstrate the global-local finite element method.     

Finite element response sensitivity analysis using three‐field mixed formulation: general theory and application to frame structures

This paper presents a method to compute response sensitivities of finite element models of structures based on a three‐field mixed formulation. The methodology is based on the direct differentiation method (DDM), and produces the response sensitivities consistent with the numerical finite element response. The general formulation is specialized to frame finite elements and details related to a newly developed steel–concrete composite frame element are provided. DDM sensitivity results are validated through the forward finite difference method (FDM) using a finite element model of a realistic steel–concrete composite frame subjected to quasi‐static and dynamic loading. The finite element model of the structure considered is constructed using both monolithic frame elements and composite frame elements with deformable shear connection based on the three‐field mixed formulation. The addition of the analytical sensitivity computation algorithm presented in this paper extends the use of finite elements based on a three‐field mixed formulation to applications that require finite element response sensitivities. Such applications include structural reliability analysis, structural optimization, structural identification, and finite element model updating. Copyright © 2006 John Wiley & Sons, Ltd.     

CLP(ℛ) as a general finite element model definition language

There are today as many finite element input languages as there are finite element programs. These input languages were developed to serve one purpose: to provide the analyst with the means to convey the necessary analysis and design criteria to specific finite element programs. They were not developed as general finite element model-definition languages, and as a result their syntax and semantics are often unnatural and contrived. In this work we do not invent a new input language. Instead, we propose and evaluate the Contraint Logic Programming language CLP(?) as a general finite element model-definition and input language.     

Uncertainty analysis of elastostatic problems incorporating a new hybrid stochastic-spectral finite element method

This article aims to present a combination of stochastic finite element and spectral finite element methods as a new numerical tool for uncertainty quantification. One of the well-established numerical methods for reliability analysis of engineering systems is the stochastic finite element method. In this article, a commonly used version of the stochastic finite element method is combined with the spectral finite element method. Furthermore, the spectral finite element method is a numerical method employing special orthogonal polynomials (e.g., Lobatto) and quadrature schemes (e.g., Gauss-Lobatto-Legendre), leading to suitable accuracy, and much less domain discretization with excellent convergence as well. The proposed method of this article is a hybrid method utilizing efficiencies of both methods for analysis of stochastically linear elastostatic problems. Moreover, a spectral finite element method is proposed for numerical solution of a Fredholm integral equation followed by the present method, to provide further efficiencies to accelerate stochastic computations. Numerical examples indicate the efficiency and accuracy of the proposed method.     

Anomalous behavior of bilinear quadrilateral finite elements for modeling cohesive cracks with XFEM/GFEM

The performance of partition‐of‐unity based methods such as the generalized finite element method or the extended finite element method is studied for the simulation of cohesive cracking. The focus of investigation is on the performance of bilinear quadrilateral finite elements using these methods. In particular, the approximation of the displacement jump field, representing cohesive cracks, by extended finite element method/generalized finite element method and its effect on the overall behavior at element and structural level is investigated. A single element test is performed with two different integration schemes, namely the Newton‐Cotes/Lobatto and the Gauss integration schemes, for the cracked interface contribution. It was found that cohesive crack segments subjected to a nonuniform opening in unstructured meshes (or an inclined crack in a structured finite element mesh) result in an unrealistic crack opening. The reasons for such behavior and its effect on the response at element level are discussed. Furthermore, a mesh refinement study is performed to analyze the overall response of a cohesively cracked body in a finite element analysis. Copyright © 2013 John Wiley & Sons, Ltd.     

蜂窝纸板结构平压性能有限元分析

研究了蜂窝纸板结构的实体建模、单元属性、网格划分等有限元建模方法,并基于该有限元模型分析了蜂窝纸板结构的平压性能,有限元分析结果与试验结果吻合良好.     

Element by element a posteriori error estimation and improvement of stress solutions for two-dimensional elastic problems

A quantitative method of a posteriori error estimation based on finite element solutions is developed. The proposed method can estimate the error of solutions obtained by the finite element method for two-dimensional elastic problems with 4-node elements. Since the error is estimated element by element, the proposed method does not require a large memory or long computing time. Not only rectangular elements but also arbitrarily shaped 4-node elements can be used in this method for estimating the error in the finite element computation with high accuracy. The finite element solutions are improved by adding the estimated errors onto the original solutions. The proposed method can be utilized for any type of linear problem if an isoparametric finite element method is used.     

A nine node finite element for analysis of geometrically non-linear shells

A nine node finite element is presented for the analysis of thin shell structures undergoing large deflection. The finite element formulation is based on the concept of degenerate solid shell element and the Hellinger-Reissner principle with independent strain. Three versions of assumed independent strain are selected to suppress spurious kinematic modes. One version leads to a finite element model which is kinematically stable at element level while the other two give globally stable models. Numerical tests indicate that the finite element model which is stable at element level may reveal the locking effect in certain cases. However, the other two models are free of locking.     

Dynamic response in the time domain by coupled boundary and finite elements

This paper reports on the development of a finite element — boundary element coupling procedure for the analysis of arbitrary shaped elastic bodies subjected to dynamic loads. The coupling is accomplished through equilibrium and compatibility considerations along the boundary element — finite element interface.Several numerical studies are performed where one part of a uniform body is treated by finite elements, whereas the remaining region is descretized by boundary elements. The examples demonstrate the influence of different finite element approaches and the applicability and the accuracy of the proposed procedure.     

An isoparametric joint/interface element for finite element analysis

A generally applicable and simple joint/interface element for three- and two-dimensional finite element analysis is presented. The proposed element can model joints/interfaces between solid finite elements and shell finite elements. The derivation of the joint element stiffness is presented and algorithms for the treatment of nonlinear joint behaviour discussed. The performance of the element is tested on typical problems involving shell-to-shell and shell-to-solid interfaces.