Implementing the Node Based Smoothed Finite Element Method as User Element in Abaqus for Linear and Nonlinear Elasticity

In this paper,the node based smoothed-strain Abaqus user element(UEL)in the framework of finite element method is introduced.The basic idea behind of the node based smoothed finite element(NSFEM)is that finite element cells are divided into subcells and subcells construct the smoothing domain associated with each node of a finite element cell[Liu,Dai and Nguyen-Thoi(2007)].Therefore,the numerical integration is globally performed over smoothing domains.It is demonstrated that the proposed UEL retains all the advantages of the NSFEM,i.e.,upper bound solution,overly soft stiffness and free from locking in compressible and nearly-incompressible media.In this work,the constant strain triangular(CST)elements are used to construct node based smoothing domains,since any complex two dimensional domains can be discretized using CST elements.This additional challenge is successfully addressed in this paper.The efficacy and robustness of the proposed work is obtained by several benchmark problems in both linear and nonlinear elasticity.The developed UEL and the associated files can be downloaded from https://github.com/nsundar/NSFEM.

Nodeless variable finite element method for heat transfer analysis by means of flux-based formulation and mesh adaptation

Based on flux-based formulation, a nodeless variable element method is developed to analyze two-dimensional steady-state and transient heat transfer problems. The nodeless variable element employs quadratic interpolation functions to provide higher solution accuracy without necessity to actually generate additional nodes. The flux-based formulation is applied to reduce the complexity in deriving the finite element equations as compared to the conventional finite element method, The solution accuracy is further improved by implementing an adaptive meshing technique to generaie finite element mesh that can adapt and move along corresponding to the solution behavior. The technique generates small elements in the regions of steep solution gradients to provide accurate solution, and meanwhile it generates larger elements in the other regions where the solution gradients are slight to reduce the computational time and the computer memory. The effectiveness of the combined procedure is demonstrated by heat transfer problems that have exact solutions. These problems tire： （a） a steady-state heat conduction analysis in a square plate subjected to a highly localized surface heating, and （b） a transient heat conduction analysis in a long plate subjected to moving heat source.

SPLITTING MODULUS FINITE ELEMENT METHOD FOR ORTHOGONAL ANISOTROPIC PLATE BENGING

Splitting modulus variational principle in linear theory of solid mechanics was introduced, the principle for thin plate was derived, and splitting modulus finite element method of thin plate was established too. The distinctive feature of the splitting model is that its functional contains one or more arbitrary additional parameters, called splitting factors, so stiffness of the model can be adjusted by properly selecting the splitting factors. Examples show that splitting modulus method has high precision and the ability to conquer some ill-conditioned problems in usual finite elements. The cause why the new method could transform the ill-conditioned problems into well-conditioned problem, is analyzed finally.

MIXED FINITE ELEMENT METHODS FOR THE SHALLOW WATER EQUATIONS INCLUDING CURRENT AND SILT SEDIMENTATION (Ⅱ)——THE DISCRETE-TIME CASE ALONG CHARACTERISTICS

The mixed finite element(MFE) methods for a shallow water equation system consisting of water dynamics equations,silt transport equation,and the equation of bottom topography change were derived.A fully discrete MFE scheme for the discrete_time along characteristics is presented and error estimates are established.The existence and convergence of MFE solution of the discrete current velocity,elevation of the bottom topography,thickness of fluid column,and mass rate of sediment is demonstrated.

Least Square Finite Element Method for Viscous Splitting of Unsteady Incompressible Navier–Stokes Equations

In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split into diffusive and convective parts in each time step. The diffusive part is discretized by the backward difference method in time and discretized by the standard Galerkin method in space. The convective part is a first-order nonlinear equation.After the linearization of the nonlinear part by Newton’s method, the convective part is also discretized by the backward difference method in time and discretized by least square scheme in space. C0-type element can be used for interpolation of the velocity and pressure in the present model. Driven cavity flow and flow past a circular cylinder are conducted to validate the present model. Numerical results agree with previous numerical results, and the model has high accuracy and can be used to simulate problems with complex geometry.

A Characteristics-Mix Stabilized Finite Element Method for Variable Density Incompressible Navier-Stokes Equations

This paper describes a characteristics-mix finite element method for the computation of incompressible Navi-er-Stokes equations with variable density. We have introduced a mixed scheme which combines a characteristics finite element scheme for treating the mass conservation equation and a finite element method to deal with the momentum equation and the divergence free constraint. The proposed method has a lot of attractive computational properties: parameter-free, very flexible, and averting the difficulties caused by the original equations. The stability of the method is proved. Finally, several numerical experiments are given to show that this method is efficient for variable density incompressible flows problem.

Finite Element Method Computations of the Acoustics of the Human Head Based on the Projection Based Interpolation Data

In this paper we present the Projection Based Interpolation (PBI) technique for construction of continuous approximation of MRI scan data of the human head. We utilize the result of the PBI algorithm to perform three dimensional (3D) Finite Element Method (FEM) simulations of the acoustics of the human head. The computational problem is a multi-physics problem modeled as acoustics coupled with linear elasticity. The computational grid contains tetrahedral finite elements with the number of equations and polynomial orders of approximation varying locally on finite element edges, faces, and interiors. We utilize our own out-of-core parallel direct solver for the solution of this multi-physics problem. The solver minimizes the memory usage by dumping out all local systems from all nodes of the entire elimination tree during the elimination phase.

THE ANALYTICAL SOLUTIONS BASED ON THE CONCEPT OF FINITE ELEMENT METHODS

On the basis of the concept of finite element methods, the rigorous analytical solutions of structural response in terms of the design variables are researched in this paper. The spatial trusses are taken as an example for the solution of the analytical expressions of the explicit displacements which are proved mathematically; then some conclusions are reached that are useful to structural sensitivity analysis and optimization. In the third part of the paper, a generalized geometric programming method is sugguested for the optimal model with the explicit displacement. Finally, the analytical solutions of the displacements of three trusses are given as examples.

TWO-GRID METHOD FOR CHARACTERISTICS FINITE-ELEMENT SOLUTION OF 2D NONLINEAR CONVECTION-DOMINATED DIFFUSION PROBLEM

For two-dimension nonlinear convection diffusion equation, a two-grid method of characteristics finite-element solution was constructed. In this method the nonlinear iterations is only to execute on the coarse grid and the fine-grid solution can be obtained in a single linear step. For the nonlinear convection-dominated diffusion equation, this method can not only stabilize the numerical oscillation but also accelerate the convergence and improve the computational efficiency. The error analysis demonstrates if the mesh sizes between coarse-grid and fine-grid satisfy the certain relationship, the two-grid solution and the characteristics finite-element solution have the same order of accuracy. The numerical is more efficient than that of characteristics example confirms that the two-grid method finite-element method.

Extended multiscale finite element method for mechanical analysis of heterogeneous materials

An extended multiscale finite element method （EMsFEM） is developed for solving the mechanical problems of heterogeneous materials in elasticity.The underlying idea of the method is to construct numerically the multiscale base functions to capture the small-scale features of the coarse elements in the multiscale finite element analysis.On the basis of our existing work for periodic truss materials, the construction methods of the base functions for continuum heterogeneous materials are systematically introduced. Numerical experiments show that the choice of boundary conditions for the construction of the base functions has a big influence on the accuracy of the multiscale solutions, thus,different kinds of boundary conditions are proposed. The efficiency and accuracy of the developed method are validated and the results with different boundary conditions are verified through extensive numerical examples with both periodic and random heterogeneous micro-structures.Also, a consistency test of the method is performed numerically. The results show that the EMsFEM can effectively obtain the macro response of the heterogeneous structures as well as the response in micro-scale,especially under the periodic boundary conditions.

Research on Analytical Method of Fatigue Characteristics of Soft Yoke Mooring System Based on Full-Scale Measurement

By focusing on the vulnerability of the structure of marine equipments, together with considering the randomness of meta-ocean load in statistics, a kind of analytical method of fatigue characteristics of marine structure based on full- scale and actual measurement of prototype is proposed. On the basis of short-term field measurement results of structural response, research is carried out on the fatigue analysis of hinge joints at the upper part of the Soft Yoke single point Mooring System （SYMS） by simultaneously monitoring the environmental load and considering the design criteria of offshore structure. Through analysis of finite element modeling, the time-histories of typical stress response are obtained, and then the assessment of fatigue damage at key hinge joints is conducted. The simulation results indicate that the proposed method can accurately analyze the fatigue damage of offshore engineering structure caused by the effect of wave load. The present analytical method of fatigue characteristics can be extended on other offshore engineering structures subjected to meta-ocean load.

Calculation and Analysis of Acoustic Characteristics of Straight-Through Perforated Pipe Muffler Based on Multilayer Sound Absorbing Material

Using the multi-physical field simulation software COMSOL,the acoustic characteristics of the multilayer sound absorbing material straight-through perforated pipe muffler are studied by the finite element method.The results show that the finite element calculation of the multilayer sound absorbing material straight-through the perforated pipe muffler agrees well with the experimental measurement results.The reliability of the finite element method for studying the acoustic performance of the straight-through perforated pipe muffler with multilayer sound absorbing materials is shown.Furthermore,the influence of some structural parameters of porous sound absorbing material and micro-perforated plate on the acoustic performance of the multilayer sound absorbing material straight-through perforated pipe muffler is analyzed.The muffler based on multilayer sound absorbing material is different from the traditional muffler.After applying the multilayer sound absorbing material to the straight-through perforated pipe muffler,the transmission loss value greatly increases,which provides new ideas and directions for future research on the muffler.

PREDICTION OF ACOUSTIC PERFORMANCE IN EXPANSION CHAMBER MUFFLERS WITH MEAN FLOW BY FINITE ELEMENT METHOD

The equation of wave propagation in a circular chamber with mean flow is obtained. Computational solution based on finite element method is employed to determine the transmission loss of expansive chamber. The effect of the mean flow and geometry (length of expansion chamber and expansion ratio)on acoustic attenuation performance is discussed, the predicted values of transmission loss of expansion chamber without and with mean flow are compared with those reported in the literature and they agree well. The accuracy of the prediction of transmission loss implies that finite element approximations are applicable to a lot of practical applications.

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.

The dimension split element-free Galerkin method for three-dimensional potential problems

This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-dimensional potential problem can be transformed into a series of two-dimensional problems. For these two-dimensional problems, the improved moving least-squares （IMLS） approximation is applied to construct the shape function, which uses an orthogonal function system with a weight function as the basis functions. The Galerkin weak form is applied to obtain a discretized system equation, and the penalty method is employed to impose the essential boundary condition. The finite difference method is selected in the splitting direction. For the purposes of demonstration, some selected numerical examples are solved using the DSEFG method. The convergence study and error analysis of the DSEFG method are presented. The numerical examples show that the DSEFG method has greater computational precision and computational efficiency than the IEFG method.

A Second Order Characteristic Mixed Finite Element Method for Convection Diffusion Reaction Equations

A combined approximate scheme is defined for convection-diffusion-reaction equations. This scheme is constructed by two methods. Standard mixed finite element method is used for diffusion term. A second order characteristic finite element method is presented to handle the material derivative term, that is, the time derivative term plus the convection term. The stability is proved and the L2-norm error estimates are derived for both the scalar unknown variable and its flux. The scheme is of second order accuracy in time increment, symmetric, and unconditionally stable.

Characteristics of permanent magnet linear synchronous motor fed by spwm inverter based on field-circuit coupled method

Presented field-circuit coupled adaptive time-stepping finite element method to study on permanent magnet linear synchronous motor (PMLSM) characteristics fed by SPWM voltage source inverter.In air-gap field where the direction or magnitude of the field is changing rapidly,the smallest elements are demanded due to high accuracy to use adaptive meshing technique.The co-simulation was used with the status space functions and time-step finite element functions,in which time-step of the status space functions was the smallest than finite element functions'.The magnitude relation of the normal elec- tromagnetic force and tangential electromagnetic force and the period were attained,and current curve was very abrupt at current zero area due to the bigger resistance and leak- age reactance,including main characteristics of motor voltage and velocity.The simulation results compare triumphantly with the experiments results.

Two-grid method for characteristic mixed finite-element solutions of nonlinear convection-diffusion equations

A two-grid method for solving nonlinear convection-dominated diffusion equations is presented. The method use discretizations based on a characteristic mixed finite-element method and give the linearization for nonlinear systems by two steps. The error analysis shows that the two-grid scheme combined with the characteristic mixed finite-element method can decrease numerical oscillation caused by dominated convections and solve nonlinear advection-dominated diffusion problems efficiently.

Combined adaptive meshing technique and characteristic-based split algorithm for viscous incompressible flow analysis

A combined characteristic-based split algorithm and all adaptive meshing technique for analyzing two-dimensional viscous incompressible flow are presented. Tile method uses the three-node triangular element with equal-order interpolation functions for all variables of tile velocity components and pressure. The main advantage of the combined nlethod is that it inlproves the sohltion accuracy by coupling an error estinla- tion procedure to an adaptive meshing technique that generates small elements in regions with a large change ill sohmtion gradients, mid at the same time, larger elements in the other regions. The performance of the combined procedure is evaluated by analyzing one test case of the flow past a cylinder, for their transient and steady-state flow behaviors.

Shape Optimal Design of Hydrodynamic Bearing Based on Dynamic Characteristics

ＳｈａｐｅＯｐｔｉｍａｌＤｅｓｉｇｎｏｆＨｙｄｒｏｄｙｎａｍｉｃＢｅａｒｉｎｇＢａｓｅｄｏｎＤｙｎａｍｉｃＣｈａｒａｃｔｅｒｉｓｔｉｃｓＣｈｅｎＹｕｎｆｅｉ（陈云飞）ＱｉａｎＲｕｉｍｉｎｇ（钱瑞明）（ＤｅｐａｒｔｍｅｎｔｏｆＭｅｃｈａｎｉｃａｌＥｎｇ...