Insight into finite element shell discretizations by use of the “basic shell mathematical model”

The objective of this paper is to gain insight into finite element discretizations of shells using the basic shell mathematical model and, in particular, regarding the sources of “locking”. We briefly review the “basic shell mathematical model” and present a formulation of shell finite elements based on this model. These shell finite elements are equivalent to the widely-used continuum mechanics based shell finite elements. We consider a free hyperboloid shell problem, which is known to be difficult to solve accurately. Using a fine mesh of MITC9 elements based on the basic shell mathematical model, a detailed analysis is performed giving the distributions of all strain terms. A similar analysis using the MITC6 shell element shows why this element locks when the shell thickness is very small.     

基于边光滑有限元法的剪切变形板几何非线性分析

为改善在计算板的几何非线性问题时有限元法系统过硬的数值缺陷,提高计算精度,在考虑剪切变形的yon Karman假设下,基于全拉格朗日描述方法,将边光滑有限元法应用于板的几何非线性分析.计算公式基于1阶剪切变形理论,并采用离散剪切间隙有效地消除剪切自锁.在三角形单元的基础上进一步形成边界光滑域,在每个光滑域内对应变进行光...     

A subgrid stabilization finite element method for incompressible magnetohydrodynamics

This paper studies a numerical scheme for approximating solutions of incompressible magnetohydrodynamic (MHD) equations that uses eddy viscosity stabilization only on the small scales of the fluid flow. This stabilization scheme for MHD equations uses a Galerkin finite element spatial discretization with Scott-Vogelius mixed finite elements and semi-implicit backward Euler temporal discretization. We prove its unconditional stability and prove how the coarse mesh can be chosen so that optimal convergence can be achieved. We also provide numerical experiments to confirm the theory and demonstrate the effectiveness of the scheme on a test problem for MHD channel flow.     

An efficient 3D stochastic finite element method for failure probability analysis of high temperature components

The loadings on high temperature components are generally complex and the discreteness of the material strength is usually great. Therefore, the two-dimensional (2D) failure probability analysis model and the deterministic finite element method (DFEM) cannot be applied to evaluate the failure probability of asymmetrical three-dimensional (3D) components. To overcome the drawbacks of the 2D model and the DFEM, an efficient 3D stochastic finite element method (SFEM) is proposed in this paper. With this method, the failure probability of components subjected to complex loadings can be estimated by using the statistical analysis of the Von Mises stresses of element nodes. Meanwhile, ANSYS and MATLAB were employed to carry out 3D parametric modeling, solving and statistical analysis. The proposed method is efficient, as is verified for two cases, and it can also be easily applied in practical engineering.     

Recent developments in the Trefftz method for finite element and boundary element applications

In 1926 E. Trefftz published a paper about a variational formulation which utilizes boundary integrals. Almost half a century later researchers became interested again in the ideas of Trefftz when the potential advantage of the Trefftz-method for an efficient use in numerical application on a computer was recognized. The concept of Trefftz can be used both for finite element and boundary element applications. A crucial ingredient of the Trefftz- method is a set of linearly independent trial functions which a priori satisfy the governing differential equations under consideration. In this paper an overview of some recent developments to construct trial functions for the Trefftz-method in a systematic manner is given. Using different types of approximation functions (singular or non-singular) we can obtain very accurate finite element and boundary element algorithms.     

Mixed least-squares finite element model for the static analysis of laminated composite plates

This paper presents a mixed finite element model for the static analysis of laminated composite plates. The formulation is based on the least-squares variational principle, which is an alternative approach to the mixed weak form finite element models. The mixed least-squares finite element model considers the first-order shear deformation theory with generalized displacements and stress resultants as independent variables. Specifically, the mixed model is developed using equal-order C⁰ Lagrange interpolation functions of high p-levels along with full integration. This mixed least-squares-based discrete model yields a symmetric and positive-definite system of algebraic equations. The predictive capability of the proposed model is demonstrated by numerical examples of the static analysis of four laminated composite plates, with different boundary conditions and various side-to-thickness ratios. Particularly, the mixed least-squares model with high-order interpolation functions is shown to be insensitive to shear-locking.     

A geometrically nonlinear shell finite element for tire vibration analysis

An axisymmetric finite element is developed which includes such features as orthotropic material properties, doubly curved geometry, and both the first and second order nonlinear stiffness terms. This element can be used to predict the equilibrium state of an axisymmetric shell structure with geometrically nonlinear large displacements. Small amplitude vibration analysis can then be performed based on this equilibrium state. The nonlinear path is predicted by using the self-correcting incremental procedure and any point on the path can be checked by using the Newton-Raphson iterative scheme. The present formulation and solution procedure are evaluated by analyzing a series of examples with results compared with alternative known solutions. Examples include: free vibration of an isotropic cylindrical shell, a conical frustum, and an orthotropic cylindrical shell; buckling of a cylindrical shell; large deflection of a clamped disk, a spherical cap, and a steel belted radial tire. The final example is a free vibration analysis of the inflated tire and the natural frequencies obtained compared well with published experimental data.     

Optimal blank design based on finite element method for blades of large Francis turbines

Metal pressing process that is widely used in industries has advantages over casting process for producing large Francis turbine blades from thick plates. Prior to the pressing process, blank design is firstly performed to determine flat blanks. The traditional trial and error approach is not applicable to blade design for Francis turbines that are not standard due to hydraulic characteristics of power plant sites. The rapid development of computing technology makes it possible to obtain optimal flat blanks by numerical modelling and simulation. In this paper, inverse finite element approach is investigated for blank design and an elasto-plastic model has been built using the well-known commercial software ANSYS. Numerical simulations for blade unfolding models with thick shell elements, solid elements and shell elements have given results with negligible differences. Unfolding tests with simple geometries have been carried out and the numerical results agree well with the analytical solutions. A large and thick shape of a Francis turbine blade for a hydropower plant has been successfully unfolded by inverse FE model. Sensibility analysis shows that the middle surface of the flat blank is independent of blade thickness. For ensuring the machining of the blade after the pressing process, a new contour is obtained by extending the boundary of the flat blank provided by the numerical model. This research may provide a useful tool for optimal blank design of Francis turbine blades.