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.

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.

Multiscale Finite Element Method for Coupling Analysis of Heterogeneous Magneto-Electro-Elastic Structures in Thermal Environment

Magneto-electro-elastic (MEE) materials, a new type of composite intelligent materials, exhibit excellent multifield coupling effects. Due to the heterogeneity of the materials, it is challenging to use the traditional finite element method (FEM) for mechanical analysis. Additionally, the MEE materials are often in a complex service environment, especially under the influence of the thermal field with thermoelectric and thermomagnetic effects, which affect its mechanical properties. Therefore, this paper proposes the efficient multiscale computational method for the multifield coupling problem of heterogeneous MEE structures under the thermal environment. The method constructs a multi-physics field with numerical base functions (the displacement, electric potential, and magnetic potential multiscale base functions). It equates a single cell of heterogeneous MEE materials to a macroscopic unit and supplements the macroscopic model with a microscopic model. This allows the problem to be solved directly on a macroscopic scale. Finally, the numerical simulation results demonstrate that compared with the traditional FEM, the multiscale finite element method (MsFEM) can achieve the purpose of ensuring accuracy and reducing the degree of freedom, and significantly improving the calculation efficiency.

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.

High Accuracy Finite Volume Method for Solving Nonlinear Aeroacoustics Problems on Unstructured Meshes

The paper presents a finite volume numerical method universally applicable for solving both linear and nonlinear aeroacoustics problems on arbitrary unstructured meshes. It is based on the vertexcentered multi-parameter scheme offering up to the 6th accuracy order achieved on the Cartesian meshes. An adaptive dissipation is added for the numerical treatment of possible discontinuities. The scheme properties are studied on a series of test cases, its efficiency is demonstrated at simulating the noise suppression in resonance-type liners.

A 3D Nonhydrostatic Compressible Atmospheric Dynamic Core by Multi-moment Constrained Finite Volume Method

A 3D compressible nonhydrostatic dynamic core based on a three-point multi-moment constrained finite-volume (MCV) method is developed by extending the previous 2D nonhydrostatic atmospheric dynamics to 3D on a terrainfollowing grid. The MCV algorithm defines two types of moments: the point-wise value (PV) and the volume-integrated average (VIA). The unknowns (PV values) are defined at the solution points within each cell and are updated through the time evolution formulations derived from the governing equations. Rigorous numerical conservation is ensured by a constraint on the VIA moment through the flux form formulation. The 3D atmospheric dynamic core reported in this paper is based on a three-point MCV method and has some advantages in comparison with other existing methods, such as uniform third-order accuracy, a compact stencil, and algorithmic simplicity. To check the performance of the 3D nonhydrostatic dynamic core, various benchmark test cases are performed. All the numerical results show that the present dynamic core is very competitive when compared to other existing advanced models, and thus lays the foundation for further developing global atmospheric models in the near future.

Characteristic-Based Time Domain Method for Cylindrically Conformal Microstrip Patch Antennas

The characteristics of a cylindrical conformal microstrip patch antenna are analyzed by using the characteristic-based time domain （CBTD） method. A governing equation in the cylindrical coordinate system is formulated directly to facilitate the analysis of cylindrically conformal microstrip patch antennas. The algorithm has second-order accuracy both in time and space domain and has the potential to eliminate the spurious wave reflection from the numerical boundaries of the computational domain, Numerical results demonstrate the important merits and accuracy of the proposed technique in computational electromagnetics,

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

裂纹缺陷限制了压电材料更广泛的应用和相关器件性能的进一步提升.为提高求解压电材料断裂参数的精度和效率,将复势函数法和Cell-based光滑有限元法引入到含裂纹平面压电材料问题中,提出了含裂纹压电材料的Cell-based光滑有限元法.以含中心裂纹压电平板的问题为例,对不同材料、裂纹长度、网格和光滑子元数下的正则强度因子进行了讨论,并与FEM作了对比,数值算例结果表明,Cell-based光滑有限元法改善了有限元法刚度"偏硬"的缺点,具有高精度、高效率的优点.

复合材料层合板的摄动随机Cell-Based光滑有限元

随机性是实际工程结构的固有特性,如何更真实地描述含随机参数结构的随机响应及统计特性,对工程结构的可靠性设计具有非常重要的意义。本文基于Cell-Based光滑有限元,采用四边形单元,推导了基于一阶剪切变形理论的复合材料层合板的光滑有限元公式,降低了网格划分要求,适应不规则网格,并采用离散剪切间隙有效地消除了剪切自锁;结合摄动法和随机场理论,导出了复合材料层合板的摄动随机光滑有限元平衡方程,并给出了结构随机响应数字特征的计算公式,求解了材料属性含随机性的复合材料层合板的随机响应问题,数值算例结果表明了本方法的有效性和准确性。

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.

Investigation of Coulomb force effects on ethylene glycol based nanofluid laminar flow in a porous enclosure

Forced convection heat transfer of ethylene glycol based nanofluid with FeOinside a porous medium is studied using the electric field. The control volume based finite element method(CVFEM) is selected for numerical simulation. The impact of the radiation parameter(R), the supplied voltage(?φ), the volume fraction of nanofluid(?), the Darcy number(Da), and the Reynolds number(Re) on nanofluid treatment is demonstrated. Results prove that thermal radiation increases the temperature gradient near the positive electrode. Distortion of isotherms increases with the enhance of the Darcy number and the Coulomb force.

Numerical Simulation for the External Combustion of Base-Bleed Projectile Using Gridless Method

The gridless method coupled with finite rate chemistry model is employed to simulate the external combustion flow fields of M864 base bleed projectile. The fluid dynamics process is described by Euler Equation in 2-D axisymmetric coordinate. The numerical method is based on least-square gridless method,and the inviscid flux is calculated by multi-component HLLC( Harten-Lax-van Leer-Contact) scheme,and a H2-CO reaction mechanism involving 9 species and 11 reactions is used. The computations are performed for the full projectile configuration of Ma = 1. 5,2,and 3. The hot air injection cases and inert cases are simulated for comparison. The numerical results show that due to the combustion in the weak region,the recirculation zone enlarges and moves downstream,the base pressure increases and the total drag force coefficient decreases. At Ma = 3. 0,the rear stagnation point shifts downstream approximate 0. 26 caliber,and the base pressure increases about 53. 4%,and the total drag force coefficient decreases to 0. 182 which agrees well with the trajectory model prediction. Due to neglecting the effects of viscosity and turbulence,there exists a certain difference at Ma = 1. 5,2. 0.

力电湿多物理场耦合Cell-Based光滑有限元法研究

为了解决压电结构在实际工程应用中受到空气湿度影响的问题,基于压电材料的本构方程、几何方程和平衡方程提出了力电湿多物理场耦合Cell-based光滑有限元模型.建立了一种新型压电俘能器,并对其简化模型的静力学性能进行求解,对比了不同湿度变化对结构性能的影响.与有限元法的计算结果对比,Cell-based光滑有限元法采用较少的网格就能够达到与有限元法采用较多的网格相同的精度,验证了本方法的正确性和有效性.研究发现,较小的湿度变化给结构的广义位移带来了较大的影响.压电结构在潮湿环境下的性能研究为复合材料产品化过程提供了理论依据.

Test Selection on Extended Finite State Machines with Provable Guarantees

Building high confidence regression test suites to validate new system versions is a challenging problem. A modelbased approach to build a regression test suite from a given test suite is described. The generated test suite includes every test that will traverse a change performed to produce the new version, and consists of only such tests to reduce the testing costs. Finite state machines extended with typed variables (EFSMs) are used to model systems and system changes are mapped to EFSM transition changes adding/deleting/replacing EFSM transitions and states. Tests are a sequence of input and expected output messages with concrete parameter values over the supported data types. An invariant is formulated to characterize tests whose runtime behavior can be accurately predicted by analyzing their descriptions along with the model. Incremental procedures to efficiently evaluate the invariant and to select tests for regression are developed. Overlaps among the test descriptions are exploited to extend the approach to simultaneously select multiple tests to reduce the test selection costs. Effectiveness of the approach is demonstrated by applying it to several protocols, Web services, and model programs extracted from a popular testing benchmark. Our experimental results show that the proposed approach is economical for regression test selection in all these examples. For all these examples, the proposed approach is able to identify all tests exercising changes more efficiently than brute-force symbolic evaluation.

Finite element model updating of a rod-type ultrasonicmotor based on response surface method

The conventional finite element model （FEM） of a rod-type ultrasonic motor is usually simplified by means of continuous composite structure. Because the actual contact characteristics between the parts of the ultrasonic motor is ignored, there is bigger error between the calculated values and experimental results. Aiming at solving problem, a new modeling method of a rod-type ultrasonic motor is presented to obtain a high-accuracy FEM. The bolt pretension and the normal contact stiffness and friction coefficient of the contact surface of ultrasonic motor are all considered in this method, and the significant parameters of working mode of the motor are selected by the response surface method, and the goal of calculating the structural response rapidly is realized by building the response surface model to replace the FEM. The result of finite element model updating shows that the average error of modal frequencies of updated model drops to 0.21% from 1.20%. The accuracy of FEM is obviously improved, which indicates that the FEM updating based on response surface method is of great application value on the design for a rod-type ultrasonic motor.

A novel simulation method for positive corona current pulses

A novel two-dimensional （2D） simulation method of positive corona current pulses is proposed. A control-volume- based finite element method （CV-FEM） is used to solve continuity equations, and the Galerkin finite element method （FEM） is used to solve Poisson＇s equation. In the proposed method, photoionization is considered by adopting an exact Helmholtz photoionization model. Furthermore, fully implicit discretization and variable time step are used to ensure the time-efficiency of the present method. Finally, the method is applied to a positive rod-plane corona problem. The numerical results are in agreement with the experimental results, and the validity of the proposed method is verified.

Structural characteristics of cement-stabilized soil bases with 3D finite element method

Cement-stabilized soil bases have been widely used in expressways due to its high strength,appropriate stiffness,good water resistance,and frost resistance.So far,the structural characteristics and mechanical behaviors of cement-stabilized soil bases were not investigated so much.In this paper,the 3D elastic-plastic finite element method(FEM)was used to analyze the mechanical behaviors and structural characteristics of cementstabilized soil bases from construction to operation.The pavement filling and the traffic loading processes were simulated,and a contact model was used to simulate the contact behavior between each layer of the pavement.Considering the construction process,the structural characteristics and mechanical behaviors of cementstabilized soil bases were studied under asphalt-concrete pavement conditions.Furthermore,the general rules of deformations and stresses in cement-stabilized soil bases under different conditions were discussed,and some suggestions were put forward for the design and construction of cement-stabilized soil bases.

基础激励下裂纹梁结构疲劳寿命预测模型

为了发展基于模型的裂纹梁结构疲劳寿命预测方法,提出了一种新的三维裂纹梁单元刚度矩阵。基于该刚度矩阵,建立基础激励下裂纹梁结构的三维有限元模型。推导三维裂纹梁单元的应力强度因子,建立三维裂纹梁单元的疲劳裂纹扩展模型。采用振动响应与裂纹扩展增量交互计算方法,计算循环载荷作用下裂纹梁结构的疲劳寿命。研究结果表明:裂纹梁结构寿命预测模型较好地体现了疲劳裂纹扩展寿命与基础激励参数、结构参数和裂纹参数之间的内在关系,反映了基础激励对裂纹梁结构疲劳寿命的影响。最后通过实验对模型进行了验证。