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.

TWO-LEVEL MULTISCALE FINITE ELEMENT METHODS FOR THE STEADY NAVIER-STOKES PROBLEM

In this article, on the basis of two-level discretizations and multiscale finite element method, two kinds of finite element algorithms for steady Navier-Stokes problem are presented and discussed. The main technique is first to use a standard finite element discretization on a coarse mesh to approximate low frequencies, then to apply the simple and Newton scheme to linearize discretizations on a fine grid. At this process, multiscale finite element method as a stabilized method deals with the lowest equal-order finite element pairs not satisfying the inf-sup condition. Under the uniqueness condition, error analyses for both algorithms are given. Numerical results are reported to demonstrate the effectiveness of the simple and Newton scheme.

THE TWO-LEVEL STABILIZED FINITE ELEMENT METHOD BASED ON MULTISCALE ENRICHMENT FOR THE STOKES EIGENVALUE PROBLEM

In this paper,we first propose a new stabilized finite element method for the Stokes eigenvalue problem.This new method is based on multiscale enrichment,and is derived from the Stokes eigenvalue problem itself.The convergence of this new stabilized method is proved and the optimal priori error estimates for the eigenfunctions and eigenvalues are also obtained.Moreover,we combine this new stabilized finite element method with the two-level method to give a new two-level stabilized finite element method for the Stokes eigenvalue problem.Furthermore,we have proved a priori error estimates for this new two-level stabilized method.Finally,numerical examples confirm our theoretical analysis and validate the high effectiveness of the new methods.

The Numerical Accuracy Analysis of Asymptotic Homogenization Method and Multiscale Finite Element Method for Periodic Composite Materials

In this paper,we discuss the numerical accuracy of asymptotic homogenization method(AHM)and multiscale finite element method(MsFEM)for periodic composite materials.Through numerical calculation of the model problems for four kinds of typical periodic composite materials,the main factors to determine the accuracy of first-order AHM and second-order AHM are found,and the physical interpretation of these factors is given.Furthermore,the way to recover multiscale solutions of first-order AHM and MsFEM is theoretically analyzed,and it is found that first-order AHM and MsFEM provide similar multiscale solutions under some assumptions.Finally,numerical experiments verify that MsFEM is essentially a first-order multiscale method for periodic composite materials.

Reduced-order finite element method based on POD for fractional Tricomi-type equation

The reduced-order finite element method （FEM） based on a proper orthogo- nal decomposition （POD） theory is applied to the time fractional Tricomi-type equation. The present method is an improvement on the general FEM. It can significantly save mem- ory space and effectively relieve the computing load due to its reconstruction of POD basis functions. Furthermore, the reduced-order finite element （FE） scheme is shown to be un- conditionally stable, and error estimation is derived in detail. Two numerical examples are presented to show the feasibility and effectiveness of the method for time fractional differential equations （FDEs）.

HIGH-ORDER RUNGE-KUTTA DISCONTINUOUS GALERKIN FINITE ELEMENT METHOD FOR 2-D RESONATOR PROBLEM

The Runge-Kutta discontinuous Galerkin finite element method （RK-DGFEM） is introduced to solve the classical resonator problem in the time domain. DGFEM uses unstructured grid discretization in the space domain and it is explicit in the time domain. Consequently it is a best mixture of FEM and finite volume method （FVM）. RK-DGFEM can obtain local high-order accuracy by using high-order polynomial basis. Numerical experiments of transverse magnetic （TM） wave propagation in a 2-D resonator are performed. A high-order Lagrange polynomial basis is adopted. Numerical results agree well with analytical solution. And different order Lagrange interpolation polynomial basis impacts on simulation result accuracy are discussed. Computational results indicate that the accuracy is evidently improved when the order of interpolation basis is increased. Finally, L^2 errors of different order polynomial basis in RK-DGFEM are presented. Computational results show that L^2 error declines exponentially as the order of basis increases.

New nonconforming finite element method for solving transient Naiver-Stokes equations

For transient Naiver-Stokes problems, a stabilized nonconforming finite element method is presented, focusing on two pairs inf-sup unstable finite element spaces, i.e., pNC/pNC triangular and pNQ/pNQ quadrilateral finite element spaces. The semi- and full-discrete schemes of the stabilized method are studied based on the pressure projection and a variational multi-scale method. It has some attractive features： avoiding higher-order derivatives and edge-based data structures, adding a discrete velocity term only on the fine scale, being effective for high Reynolds number fluid flows, and avoiding increased computation cost. For the full-discrete scheme, it has second-order estimations of time and is unconditionally stable. The presented numerical results agree well with the theoretical results.

Meshfree Finite Volume Element Method for Constrained Optimal Control Problem Governed by Random Convection Diffusion Equations

In this paper,we investigate a stochastic meshfree finite volume element method for an optimal control problem governed by the convection diffusion equations with random coefficients.There are two contributions of this paper.Firstly,we establish a scheme to approximate the optimality system by using the finite volume element method in the physical space and the meshfree method in the probability space,which is competitive for high-dimensional random inputs.Secondly,the a priori error estimates are derived for the state,the co-state and the control variables.Some numerical tests are carried out to confirm the theoretical results and demonstrate the efficiency of the proposed method.

A Multiscale Multilevel Monte Carlo Method for Multiscale Elliptic PDEs with Random Coefficients

We propose a multiscale multilevel Monte Carlo(MsMLMC)method to solve multiscale elliptic PDEs with random coefficients in the multi-query setting.Our method consists of offline and online stages.In the offline stage,we construct a small number of reduced basis functions within each coarse grid block,which can then be used to approximate the multiscale finite element basis functions.In the online stage,we can obtain the multiscale finite element basis very efficiently on a coarse grid by using the pre-computed multiscale basis.The MsMLMC method can be applied to multiscale RPDE starting with a relatively coarse grid,without requiring the coarsest grid to resolve the smallestscale of the solution.We have performed complexity analysis and shown that the MsMLMC offers considerable savings in solving multiscale elliptic PDEs with random coefficients.Moreover,we provide convergence analysis of the proposed method.Numerical results are presented to demonstrate the accuracy and efficiency of the proposed method for several multiscale stochastic problems without scale separation.

时空降阶单元及其自适应分析初探

基于初值问题和边值问题中降阶单元的成功实现,进一步提出时空降阶单元,以弹性弦振动方程(波动方程)为例,构造了自适应有限元分析算法。时空降阶单元继承了初值和边值问题中降阶单元的所有优点,可以给出按最大模度量满足误差限的降阶单元的解答。该文对这一研究进展做一简要介绍,并给出包括若干强迫振动以及移动荷载反应的典型算例以展示本法的可行性、有效性和可靠性。

A New Formulation of the Scaled Boundary Finite Element Method for Heterogeneous Media:Application to Heat Transfer Problems

The solution to heat transfer problems in two-dimensional heterogeneous media is attended based on the scaled boundary finite element method(SBFEM)coupled with equilibrated basis functions(EqBFs).The SBFEM reduces the model order by scaling the boundary solution onto the inner element.To this end,tri-lateral elements are emanated from a scaling center,followed by the development of a semi-analytical solution along the radial direction and a finite element solution along the circumferential/boundary direction.The discretization is thus limited to the boundaries of the model,and the semi-analytical radial solution is found through the solution of an eigenvalue problem,which restricts the methods’applicability to heterogeneous media.In this research,we first extracted the SBFEM formulation considering the heterogeneity of the media.Then,we replaced the semi-analytical radial solution with the EqBFs and removed the eigenvalue solution step from the SBFEM.The varying coefficients of the partial differential equation(PDE)resulting from the heterogeneity of the media are replaced by a finite series in the radial and circumferential directions of the element.A weighted residual approach is applied to the radial equation.The equilibrated radial solution series is used in the new formulation of the SBFEM.

Calculation method of ship collision force on bridge using artificial neural network

Ship collision on bridge is a dynamic process featured by high nonlinearity and instantaneity. Calculating ship-bridge collision force typically involves either the use of design-specification-stipulated equivalent static load, or the use of finite element method (FEM) which is more time-consuming and requires supercomputing resources. In this paper, we proposed an alternative approach that combines FEM with artificial neural network (ANN). The radial basis function neural network (RBFNN) employed for calculating the impact force in consideration of ship-bridge collision mechanics. With ship velocity and mass as the input vectors and ship collision force as the output vector, the neural networks for different network parameters are trained by the learning samples obtained from finite element simulation results. The error analyses of the learning and testing samples show that the proposed RBFNN is accurate enough to calculate ship-bridge collision force. The input-output relationship obtained by the RBFNN is essentially consistent with the typical empirical formulae. Finally, a special toolbox is developed for calculation efficiency in application using MATLAB software.

An Inf-Sup Stabilized Finite Element Method by Multiscale Functions for the Stokes Equations

In the paper,an inf-sup stabilized finite element method by multiscale functions for the Stokes equations is discussed.The key idea is to use a PetrovGalerkin approach based on the enrichment of the standard polynomial space for the velocity component with multiscale functions.The inf-sup condition for P_(1)-P_(0)triangular element(or Q_(1)-P_(0)quadrilateral element)is established.The optimal error estimates of the stabilized finite element method for the Stokes equations are obtained.

Generalized Multiscale Finite Element Methods.Nonlinear Elliptic Equations

In this paper we use the Generalized Multiscale Finite Element Method(GMsFEM)framework,introduced in[26],in order to solve nonlinear elliptic equations with high-contrast coefficients.The proposed solution method involves linearizing the equation so that coarse-grid quantities of previous solution iterates can be regarded as auxiliary parameters within the problem formulation.With this convention,we systematically construct respective coarse solution spaces that lend themselves to either continuous Galerkin(CG)or discontinuous Galerkin(DG)global formulations.Here,we use Symmetric Interior Penalty Discontinuous Galerkin approach.Both methods yield a predictable error decline that depends on the respective coarse space dimension,and we illustrate the effectiveness of the CG and DG formulations by offering a variety of numerical examples.

Gaussian Radial Basis Function interpolation in vertical deformation analysis

In many deformation analyses,the partial derivatives at the interpolated scattered data points are required.In this paper,the Gaussian Radial Basis Functions(GRBF)is proposed for the interpolation and differentiation of the scattered data in the vertical deformation analysis.For the optimal selection of the shape parameter,which is crucial in the GRBF interpolation,two methods are used:the Power Gaussian Radial Basis Functions(PGRBF)and Leave One Out Cross Validation(LOOCV)(LGRBF).We compared the PGRBF and LGRBF to the traditional interpolation methods such as the Finite Element Method(FEM),polynomials,Moving Least Squares(MLS),and the usual GRBF in both the simulated and actual Interferometric Synthetic Aperture Radar(InSAR)data.The estimated results showed that the surface interpolation accuracy was greatly improved by LGRBF and PGRBF methods in comparison withFEM,polynomial,and MLS methods.Finally,LGRBF and PGRBF interpolation methods are used to compute invariant vertical deformation parameters,i.e.,changes in Gaussian and mean Curvatures in the Groningen area in the North of Netherlands.

An Improved Model Reduction Method on AIMs for N-S Equations Using Multilevel Finite Element Method and Hierarchical Basis

A numerical method is proposed to approach the Approximate Inertial Man-ifolds(AIMs)in unsteady incompressible Navier-Stokes equations,using multilevel fi-nite element method with hierarchical basis functions.Following AIMS,the unknown variables,velocity and pressure in the governing equations,are divided into two com-ponents,namely low modes and high modes.Then,the couplings between low modes and high modes,which are not accounted by standard Galerkin method,are consid-ered by AIMs,to improve the accuracy of the numerical results.Further,the multilevel finite element method with hierarchical basis functions is introduced to approach low modes and high modes in an efficient way.As an example,the flow around airfoil NACA0012 at different angles of attack has been simulated by the method presented,and the comparisons show that there is a good agreement between the present method and experimental results.In particular,the proposed method takes less computing time than the traditional method.As a conclusion,the present method is efficient in numer-ical analysis of fluid dynamics,especially in computing time.

降阶单元的新进展:内置了最大模误差估计器的自适应有限元法初探

该文基于初值问题中降阶单元的成功实践,进一步对一般边值问题提出无需超收敛计算、无需结构化网格、无需结点位移修正的降阶单元;进而提出以降阶单元作为最终解,且内置了最大模误差估计器的自适应有限元算法。该文对这一研究进展做简要介绍,并给出一维和二维边值问题的初步算例,以展示本法的可行性、有效性和可靠性。

Sub-region Alopex Optimization Method with RSM for Design of Permanent Magnet Machines

Based on an Alopex optimization algorithm and a response surface model(RSM),a hybrid sub-region methodology is presented to solve the optimal design problems of permanent magnet(PM)machines.The Alopex optimization method is processed both in subspace and in global solution space.In order to decrease the computing time,a multi quadric radial basis function(MQRBF)is embedded in the optimization.The proposed method speeds up the convergence rate while keeps the accuracy of the solution.A numerical experiment is given to validate the efficiency and effectiveness of the method.

杯形谐波减速器齿廓修形方法及寿命预测分析

杯形谐波减速器柔轮与波发生器装配时产生的径向变形量理论与实际不符,柔轮与刚轮轮齿出现啮合干涉,基于柔轮装配的变形情况提出了柔轮齿廓修形方法。通过有限元法装配仿真确定了合适的修形量,修形后柔轮经过装配和运转有限元分析,最后基于Miner线性疲劳理论采用Fe-safe软件对柔轮寿命进行了预测。结果表明,修形后柔轮装配时最大应力从962.2 MPa降低至532.7 MPa,负载运转时应力为609.9 MPa,轮齿应力有较大改善且啮合干涉情况得到解决,柔轮寿命循环次数为4.28×10^(6)。