ALE FINITE ELEMENT ANALYSIS OF THE OPENING AND CLOSING PROCESS OF THE ARTIFICIAL MECHANICAL VALVE

ＡＬＥＦＩＮＩＴＥＥＬＥＭＥＮＴＡＮＡＬＹＳＩＳＯＦＴＨＥＯＰＥＮＩＮＧＡＮＤＣＬＯＳＩＮＧＰＲＯＣＥＳＳＯＦＴＨＥＡＲＴＩＦＩＣＩＡＬＭＥＣＨＡＮＩＣＡＬＶＡＬＶＥ（张建海）（陈大鹏）（邹盛铨）ＡＬＥＦＩＮＩＴＥＥＬＥＭＥＮＴＡＮＡＬＹＳＩＳＯＦＴ...     

A STUDY ON SHIP VIBRATION USING FINITE ELEMENT METHOD

when the beam theory was used to calculate ship hull vibration,greater discrepancies were found between theoretical cal-culations and actual measurements especially at higher modes.Thusthe beam model cannot be considered as a practical one for higher-mode calculations.This paper presents the application of two-dimensional finite element model for the calculation of ship ver-tical vibration.Using the multi-element structural dynamic a-nalysis program DDJ(DL)developed by ourselves,the hull vibra-tion analysis of two ships(vessel A and vessel B)was carriedout on the Model-709 Computer made in the People’s Republic ofChina.The results of the calculation,when compared with actualmeasurements.show that the two-dimensional model is much moreefficient than the traditional beam model.The agreement be-tween the calculations and measurements has been improved greatly,and this discrepancy at the4th·and5th-modes has decreased towithin5%as compared to that of more than20%in the traditionalmodel.Furthermore.the mode     

THE CONSTRUCTION OF LARGE ELEMENT IN FINITE ELEMENT METHOD

In the usual finite element method.the order of the inter-polation in an element is kept unchanged.and the accuracyis raised by subdividing the grid denser and denser.Alter-natively.in the large element method.the grid is kept un-changed.and the terms of approximate series in the elementare increased to raise the accuracy.In this paper,a method for constructing large elementsis presented.when using this method,two sets of variables.one set defined inside the element.and the other defined onthe boundary of the element.are adopted.Then,these twosets of variables are combined by the hybrid-penalty functionmethod.This method can be applied to any elliptic equationsin a domain With arbitrary shape and arbitrary complex boun-dary condition.It is proved with strict mathematical methodin this paper.that in general cases,the accuracy of this me-thod is much higher than that of the usual element and thelarge element method presented in[7].Therefore.the deqreesof freedom needed in this method are much fewer than t     

THE FINITE ELEMENT METHOD OF SINGULAR PERTURBATION PROBLEM

In this paper we construct new finite element subspace using polynomials of different degrees and the new finite element scheme is established. The convergence of the scheme and the stability of the reduced difference equation are proved.     

THE RATIONALISM THEORY AND ITS FINITE ELEMENT ANALYSIS METHOD OF SHELL STRUCTURES

In this paper,a kind of rationalism theory of shell is established which is of differentmechanic characters in tension and in compression,and the finite element numericalanalysis method is also described.     

FUZZY ARITHMETIC AND SOLVING OF THE STATIC GOVERNING EQUATIONS OF FUZZY FINITE ELEMENT METHOD

The key component of finite element analysis of structures with fuzzy parameters, which is associated with handling of some fuzzy information and arithmetic relation of fuzzy variables , was the solving of the governing equations of fuzzy finite element method. Based on a given interval representation of fuzzy numbers, some arithmetic rules of fuzzy numbers and fuzzy variables were developed in terms of the properties of interval arithmetic. According to the rules and by the theory of interval finite element method, procedures for solving the static governing equations of fuzzy finite element method of structures were presented. By the proposed procedure, the possibility distributions of responses of fuzzy structures can be generated in terms of the membership functions of the input fuzzy numbers. It is shown by a numerical example that the computational burden of the presented procedures is low and easy to implement. The effectiveness and usefulness of the presented procedures are also illustrated.     

ANALYSIS OF COMPOSITE LAMINATE BEAMS USING COUPLING CROSS-SECTION FINITE ELEMENT METHOD

Beams and plates manufactured from laminates of composite materials have distinct advantages in a significant number of applications. However, the anisotropy arising from these materials adds a significant degree of complexity, and thus time, to the stress and deformation analyses of such components, even using numerical approaches such as finite elements. The analysis of composite laminate beams subjected to uniform extension, bending, and/or twisting loads was performed by a novel implementation of the usual finite element method. Due to the symmetric features of the deformations, only a thin slice of the beam to be analysed needs to be modelled. Conventional three-dimensional solid finite elements were used for the structural discretization. The accurate deformation relationships were formulated and implemented through the coupling of nodal translational degrees of freedom in the numerical analysis. A sample solution for a rectangular composite laminate beam is presented to show the validity and accuracy of the proposed method.     

SOLUTION OF THE CONNECTION PROBLEMS BETWEEN A FINITE HOLED PLATE AND A STIFFENER BY USING THE PARTITIONING CONCEPT OF THE GENERALIZED VARIATIONAL METHOD

In this paper a new finite element method is presented.in which complex functions arechosen to be the finite element model and the partitioning concept of the generalizedvariational method is utilized.The stress concentration factors for a finite holed platewelded by a stiffener are calculated and the analytical solutions in series form are obtained.From some computer trials it is demonstrated that the problem of displacementcompatibility and continuity of tractions between the holed plate and the stiffener issuccessfully analysed by using this method.Since only three elements need to beformulated.relatively less storage is required than the usual finite element methods.Furthermore,the accurary of solutions is improved and the computer time requirements areconsiderably reduced.Numerical results of stress concentration factors and stresses alongthe welded-line which may be referential to engineers are shown in tables.     

ON THE COMPACTNESS OF QUASI-CONFORMING ELEMENT SPACES AND THE CONVERGENCE OF QUASI-CONFORMING ELEMENT METHOD

In this paper, the compactness of quasi-conforming element spaces and theconvergence of quasi-conforming element method are discussed The well-known Rellichcompactness theorem is generalized to the sequences of quasi-conforming element spaceswith certain properties. and the generalized Poincare inequality.The generalized Friedrichsinequality and the generalized inequality of Poincare-Friedrichs are proved true for them.The error estimates are also given. It is shown that the quasi-conforming element method isconvergent if the quasi-conforming element spaces have the approximability and the strongcontinuity, and satisfy the rank condition of element and pass the test IPT As practicalexamples, 6-paramenter. 9-paramenter, 12-paramenter, 15-parameter, 18-parameter and21-paramenter quasi-conforming elements are shown to be convergent, and theirL~(2‘2)(Ω)-errors are O(h_τ)、 O(h_τ)、O(h_τ~2)、O(h_τ~2), O(h_τ~3), and O(h_τ~4) respectively.     

A 3D finite element ALE method using an approximate Riemann solution

Arbitrary Lagrangian–Eulerian finite volume methods that solve a multidimensional Riemann‐like problem at the cell center in a staggered grid hydrodynamic (SGH) arrangement have been proposed. This research proposes a new 3D finite element arbitrary Lagrangian–Eulerian SGH method that incorporates a multidimensional Riemann‐like problem. Two different Riemann jump relations are investigated. A new limiting method that greatly improves the accuracy of the SGH method on isentropic flows is investigated. A remap method that improves upon a well‐known mesh relaxation and remapping technique in order to ensure total energy conservation during the remap is also presented. Numerical details and test problem results are presented. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical simulation of impact forces by a falling sphere in air using ALE finite element method

ABSTRACT

A numerical method is developed to simulate the process that a falling rigid sphere hits rigid ground and bounces back in air. The problem is treated as fluid-structure interaction problem based on the ALE finite element flow analysis. In order to introduce the numerical process of impact into the present staggered fluid-structure time marching algorithm, the impact force is applied to the equation of motion of the sphere. The magnitude of the impact force is determined by iteration so that the velocity of the sphere after impact converges to zero. Application of the impact force at a single time instant causes unphysical pressure oscillation. This has been suppressed by applying the impact force smoothly over multiple short time steps. In the present method impulse is evaluated instead of impact force. Computations with different density ratio of the sphere to air showed effect of the air on the sphere motion.     

Dynamic analysis of baffled fuel‐storage tanks using the ALE finite element method

This paper is concerned with the parametric investigation on the structural dynamic response of moving fuel‐storage tanks with baffles. Since the structural dynamic behaviour is strongly coupled with interior liquid motion, the design of a fuel‐storage tank securing the structural stability becomes the appropriate suppression of the flow motion, which is in turn related to the baffle design. In order to numerically investigate the parametric dynamic characteristics of moving tanks, we employ the arbitrary Lagrangian–Eulerian (ALE) finite element method that is widely being used to deal with the problems with free surface, moving boundary, large deformation and interface contact. Following the theoretical and numerical formulations of fluid‐structure interaction problems, we present parametric numerical results of a cylindrical fuel‐storage tank moving with uniform vertical acceleration, with respect to the baffle number and location, and the inner‐hole diameter. Copyright © 2003 John Wiley & Sons, Ltd.     

Effects of the particle Stokes number on wind turbine airfoil erosion

Under natural conditions, wind turbines are inevitably eroded by the action of sand-wind flow. To further investigate the effects of dust drift on the erosion of the wind turbine blades in sand-wind environments, the effects of the wind velocity, particle diameter, and particle density on the erosion of wind turbine airfoils are studied, and the effects of the particle Stokes number on the airfoil erosion are discussed. The results show that, when the angle of attack(AOA) is 6.1°, there will be no erosion on the airfoil surface if the particle Stokes number is lower than 0.013 5, whereas erosion will occur if the particle Stokes number is higher than 0.015 1. Therefore, there exists a critical range for the particle Stokes number. When the particle Stokes number is higher than the maximum value in the critical range, airfoil erosion will occur. The result is further confirmed by changing the particle diameter, particle density, and inflow speed. It is shown that the erosion area on the airfoil and the maximum erosion rate are almost equal under the same particle Stokes number and AOA. The extent of airfoil erosion increases when the particle Stokes number increases, and the critical particle Stokes number increases when the AOA increases. Moreover, the geometric shape of the airfoil pressure surface greatly affects the airfoil erosion, especially at the curvature near the leading edge.     

A finite element model for the simulation of lost foam casting

In this paper, we present a numerical model to simulate the lost foam casting process. We introduce this particular casting first in order to capture the different physical processes in play during a casting. We briefly comment on the possible physical and numerical models used to envisage the numerical simulation. Next we present a model which aims to solve ‘part of’ the complexities of the casting, together with a simple energy budget that enables us to obtain an equation for the velocity of the metal front advance. Once the physical model is established we develop a finite element method to solve the governing equations. The numerical and physical methodologies are then validated through the solution of a two‐ and a three‐dimensional example. Finally, we discuss briefly some possible improvements of the numerical model in order to capture more physical phenomena. Copyright © 2004 John Wiley & Sons, Ltd.     

Effective electroelastic constants for three-phase confocal elliptical cylinder model in piezoelectric quasicrystal composites

A three-phase confocal elliptical cylinder model is proposed to analyze micromechanics of one-dimensional hexagonal piezoelectric quasicrystal (PQC) composites. Exact solutions of the phonon, phason, and electric fields are obtained by using the conformal mapping combined with the Laurent expansion technique when the model is subject to far-field anti-plane mechanical and in-plane electric loadings. The effective electroelastic constants of several different composites made up of PQC, quasicrystal (QC), and piezoelectric (PE) materials are predicted by the generalized self-consistent method. Numerical examples are conducted to show the effects of the volume fraction and the cross-sectional shape of inclusion (or fiber) on the effective electroelastic constants of these composites. Compared with other micromechanical methods, the generalized selfconsistent and Mori-Tanaka methods can predict the effective electroelastic constants of the composites consistently.     

A finite element formulation satisfying the discrete geometric conservation law based on averaged Jacobians

In this article, a new methodology for developing discrete geometric conservation law (DGCL) compliant formulations is presented. It is carried out in the context of the finite element method for general advective–diffusive systems on moving domains using an ALE scheme. There is an extensive literature about the impact of DGCL compliance on the stability and precision of time integration methods. In those articles, it has been proved that satisfying the DGCL is a necessary and sufficient condition for any ALE scheme to maintain on moving grids the nonlinear stability properties of its fixed‐grid counterpart. However, only a few works proposed a methodology for obtaining a compliant scheme. In this work, a DGCL compliant scheme based on an averaged ALE Jacobians formulation is obtained. This new formulation is applied to the θ family of time integration methods. In addition, an extension to the three‐point backward difference formula is given. With the aim to validate the averaged ALE Jacobians formulation, a set of numerical tests are performed. These tests include 2D and 3D diffusion problems with different mesh movements and the 2D compressible Navier–Stokes equations. Copyright © 2011 John Wiley & Sons, Ltd.     

Reliability of elasto-plastic structure using finite element method

A solution of probabilistic FEM for elastic-plastic materials is presented based on the incremental theory of plasticity and a modified initial stress method. The formulations are deduced through a direct differentiation scheme. Partial differentiation of displacement, stress and the performance function can be iteratively performed with the computation of the mean values of displacement and stress. The presented method enjoys the efficiency of both the perturbation method and the finite difference method, but avoids the approximation during the partial differentiation calculation. In order to improve the efficiency, the adjoint vector method is introduced to calculate the differentiation of stress and displacement with respect to random variables. In addition, a time-saving computational method for reliability index of elastic-plastic materials is suggested based upon the advanced First Order Second Moment (FOSM) and by the usage of Taylor expansion for displacement. The suggested method is also applicable to 3-D cases. The project supported by the Research Grant Council of Hong Kong (HKUST 722196E, 6039197E), the National Natural Science Foundation of China(59809003) and the Foundation of University Key Teacher by the Chinese Ministry of Education     

Fractional four-step finite element method for analysis of thermally coupled fluid-solid interaction problems

An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal stress in the solid. The fractional four-step finite element method and the streamline upwind Petrov-Galerkin (SUPG) method are used to analyze the viscous thermal flow in the fluid. Analyses of the heat transfer and the thermal stress in the solid are performed by the Galerkin method. The second-order semiimplicit Crank-Nicolson scheme is used for the time integration. The resulting nonlinear equations are linearized to improve the computational efficiency. The integrated analysis method uses a three-node triangular element with equal-order interpolation functions for the fluid velocity components, the pressure, the temperature, and the solid displacements to simplify the overall finite element formulation. The main advantage of the present method is to consistently couple the heat transfer along the fluid-solid interface. Results of several tested problems show effiectiveness of the present finite element method, which provides insight into the integrated fluid-thermal-structural interaction phenomena.