Self-adaptive one-dimensional nonlinear finite element method based on element energy projection method

The element energy projection （EEP） method for computation of super- convergent resulting in a one-dimensional finite element method （FEM） is successfully used to self-adaptive FEM analysis of various linear problems, based on which this paper presents a substantial extension of the whole set of technology to nonlinear problems. The main idea behind the technology transfer from linear analysis to nonlinear analysis is to use Newton＇s method to linearize nonlinear problems into a series of linear problems so that the EEP formulation and the corresponding adaptive strategy can be directly used without the need for specific super-convergence formulation for nonlinear FEM. As a re- sult, a unified and general self-adaptive algorithm for nonlinear FEM analysis is formed. The proposed algorithm is found to be able to produce satisfactory finite element results with accuracy satisfying the user-preset error tolerances by maximum norm anywhere on the mesh. Taking the nonlinear ordinary differential equation （ODE） of second-order as the model problem, this paper describes the related fundamental idea, the imple- mentation strategy, and the computational algorithm. Representative numerical exam- ples are given to show the efficiency, stability, versatility, and reliability of the proposed approach.

Self-adaptive strategy for one-dimensional finite element method based on EEP method with optimal super-convergence order

Based on the newly-developed element energy projection （EEP） method with optimal super-convergence order for computation of super-convergent results, an improved self-adaptive strategy for one-dimensional finite element method （FEM） is proposed. In the strategy, a posteriori errors are estimated by comparing FEM solutions to EEP super-convergent solutions with optimal order of super-convergence, meshes are refined by using the error-averaging method. Quasi-FEM solutions are used to replace the true FEM solutions in the adaptive process. This strategy has been found to be simple, clear, efficient and reliable. For most problems, only one adaptive step is needed to produce the required FEM solutions which pointwise satisfy the user specified error tolerances in the max-norm. Taking the elliptical ordinary differential equation of the second order as the model problem, this paper describes the fundamental idea, implementation strategy and computational algorithm and representative numerical examples are given to show the effectiveness and reliability of the proposed approach.

SELF-ADAPTIVE STRATEGY FOR ONE-DIMENSIONAL FINITE ELEMENT METHOD BASED ON ELEMENT ENERGY PROJECTION METHOD

Based on the newly-developed element energy projection （EEP） method for computation of super-convergent results in one-dimensional finite element method （FEM）, the task of self-adaptive FEM analysis was converted into the task of adaptive piecewise polynomial interpolation. As a result, a satisfactory FEM mesh can be obtained, and further FEM analysis on this mesh would immediately produce an FEM solution which usually satisfies the user specified error tolerance. Even though the error tolerance was not completely satisfied, one or two steps of further local refinements would be sufficient. This strategy was found to be very simple, rapid, cheap and efficient. Taking the elliptical ordinary differential equation of second order as the model problem, the fundamental idea, implementation strategy and detailed algorithm are described. Representative numerical examples are given to show the effectiveness and reliability of the proposed approach.

SOLVING 2-D NONLINEAR COUPLED HEAT AND MOISTURE TRANSFER PROBLEMS VIA A SELF-ADAPTIVE PRECISE ALGORITHM IN TIME DOMAIN

A self-adaptive precise algorithm in the time domain was employed to solve 2-D nonlinear coupled heat and moisture transfer problems. By expanding variables at a discretized time interval, the variations of variables can be described more precisely,and a nonlinear coupled initial and boundary value problem was converted into a series of recurrent linear boundary value problems which are solved by FE technique. In the computation, no additional assumption and the nonlinear iteration are required, and a criterion for self-adaptive computation is proposed to maintain sufficient computing accuracy for the change sizes of time steps. In the numerical comparison, the variations of material properties with temperature, moisture content, and both temperature and moisture content are taken into account, respectively. Satisfactory results have been obtained, indicating that the proposed approach is capable of dealing with complex nonlinear problems.

Low-velocity Impact Damage Analysis of Composite Laminates Using Self-adapting Delamination Element Method

On the basis ofa 2D 4-node Mindlin shell element method, a novel self-adapting delamination finite element method is presented, which is developed to model the delamination damage of composite laminates. In the method, the sublaminate elements are generated automatically when the delamination damage occurs or extends. Thus, the complex process and state of delamination damage can be simulated practically with high efficiency for both analysis and modeling. Based on the self-adapting delamination method, linear dynamic finite element damage analysis is performed to simulate the low-velocity impact damage process of three types of mixed woven composite laminates. Taking the frictional force among sublaminations during delaminating and the transverse normal stress into account, the analytical results are consistent with those of the experimental data.

REFINED HYBRID MINDLIN PLATE ELEMENT FOR GEOMETRICALLY NON LINEAR ANALYSIS

Based upon a generalized variational principle, which relaxed the inter element continuity requirements, a novel refined hybrid Mindlin plate element is developed, its non linear element stiffness matrices are decomposed into a series of matrices with respect to the assumed strain modes. The formulation presented in this paper is different from any other non linear mixed/hybrid element formulation all successful experience of linear hybrid formulation is absorbed into the formulation(adding non conforming modes and realizing orthogonalization) Numerical results show that the present approach is more effective than any other non linear hybrid element formulation over the accuracy and computational efficiency. In addition, non conforming modes can also overcome the shear locking effect.

ANALYSIS OF NONLINEAR PIEZOELECTRIC CIRCULAR SHALLOW SPHERICAL SHELLS BY DIFFERENTIAL QUADRATURE ELEMENT METHOD

The static behavior of piezoelectric circular spherical shallow shells under both electrical and mechanical loads is studied by using the differential quadrature element method (DQEM). Geometrical nonlinearity effect is considered. Detailed formulations and procedures are given for the first time. Several examples are analyzed and accurate results are obtained by the DQEM. Based on the results in this paper, one may conclude that the DQEM is a useful tool for obtaining solutions of structural elements. It can be seen that the shell shape may be theore tically controlled and snap through may occur when the applied voltage reaches a critical value even without mechanical load for certain geometric configurations.

Linearity of Removal in the Proeess of Optical Element Polishing

In the light of some assumptions that are very close to the practical working conditions,a very complicated polishing process of optical element can be simplified as a linear and shift invariant system that is relatd only to the speed,pres- sure and time of processing.In polishing,the removed material can be represented and entreated by the convolution of the removal function of polishing head and the dwell function.The properties of removal function are presented.The assumptions and methods given by the author have been shown to be correct and applicable by experiments using a ring lap to polish the optical surfac.

A LOCKING-FREE ANISOTROPIC NONCONFORMING FINITE ELEMENT FOR PLANAR LINEAR ELASTICITY PROBLEM

The main aim of this article is to study the approximation of a locking-free anisotropic nonconforming finite element for the pure displacement boundary value problem of planar linear elasticity. The optimal error estimates are obtained by using some novel approaches and techniques. The method proposed in this article is robust in the sense that the convergence estimates in the energy and L^2-norms are independent-of the Lame parameter λ.

Optimal design of butterfly-shaped linear ultrasonic motor using finite element method and response surface methodology

A new method for optimizing a butterfly-shaped linear ultrasonic motor was proposed to maximize its mechanical output. The finite element analysis technology and response surface methodology were combined together to realize the optimal design of the butterfly-shaped linear ultrasonic motor. First, the operation principle of the motor was introduced. Second, the finite element parameterized model of the stator of the motor was built using ANSYS parametric design language and some structure parameters of the stator were selected as design variables. Third, the sample points were selected in design variable space using latin hypercube Design. Through modal analysis and harmonic response analysis of the stator based on these sample points, the target responses were obtained. These sample points and response values were combined together to build a response surface model. Finally, the simplex method was used to find the optimal solution. The experimental results showed that many aspects of the design requirements of the butterfly-shaped linear ultrasonic motor have been fulfilled. The prototype motor fabricated based on the optimal design result exhibited considerably high dynamic performance, such as no-load speed of 873 ram/s, maximal thrust of 27.5 N, maximal efficiency of 43%, and thrust-weight ratio of 45.8.

Finite element analysis of the effect of micro-pore defect on linear friction welding of medium carbon steel

Micro-pore is a very common material defect. In the present paper, the temperature fields of medium carbon steel joints with and without micro-pore defect during linear friction welding （LFW） were investigated by using finite element method. The effect of micro-pore defect on the axial shortening of joints during LFW was examined. The x- and y-direction displacements of micro-pore during the LFW process were also studied. In addition, the shape of micro-pore after LFW was observed. The heat conducted from the weld inteace to the specimen interior. The fluctuation range of the temperature curves for the joint with micro-pore is larger than that without micro-pore. Position of micro-pore changes with the change of the friction time. The circular shape of micro-pore becomes oval after welding.

Simple experimental scheme of teleporting a two-qubit state via linear optical elements

We propose a simple experimental scheme in which an unknown two-qubit state is faithfully and deterministically teleported from Alice to Bob. The scheme is constructed with four photons from parametric down conversion, linear optical elements, and conventional photon detectors, all of which are available in current technology. It is shown that the probability of successful teleportation ideally reaches 100% based on single-photon two-qubit-assisted Bell-state measurement, which can distinguish all four Bell-states simultaneously via conventional photon detectors. By generalizing the scheme, the teleportation of an unknown multi-qubit system can also be realized.

Multiple linear system techniques for 3D finite element method modeling of direct current resistivity

The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and stored in two parts separately. One part is associated with the volume integral and the other is associated with the subsurface boundary integral. The equivalent multiple linear systems with closer right-hand sides than the original systems were constructed. A recycling Krylov subspace technique was employed to solve the multiple linear systems. The solution of the seed system was used as an initial guess for the subsequent systems. The results of two numerical experiments show that the improved algorithm reduces the iterations and CPU time by almost 50%, compared with the classical preconditioned conjugate gradient method.

A streamline diffusion nonconforming finite element method for the time-dependent linearized Navier-Stokes equations

A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretization. Crouzeix-Raviart nonconforming finite element approximation, namely, nonconforming （P1）2 - P0 element, is used for the velocity and pressure fields with the streamline diffusion technique to cope with usual instabilities caused by the convection and time terms. Stability and error estimates are derived with suitable norms.

High-efciency improved symmetric successive over-relaxation preconditioned conjugate gradient method for solving large-scale finite element linear equations

Fast solving large-scale linear equations in the finite element analysis is a classical subject in computational mechanics. It is a key technique in computer aided engineering （CAE） and computer aided manufacturing （CAM）. This paper presents a high-efficiency improved symmetric successive over-relaxation （ISSOR） preconditioned conjugate gradient （PCG） method, which maintains lelism consistent with the original form. Ideally, the by 50% as compared with the original algorithm. the convergence and inherent paralcomputation can It is suitable for be reduced nearly high-performance computing with its inherent basic high-efficiency operations. By comparing with the numerical results, it is shown that the proposed method has the best performance.

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.

Finite Element Analysis and Linear Regression of Maximum Temperature for Inner Wall of Chimney Foundation

The uniform design method was adopted and the twenty-four groups of different geometric and physical pa-rameters were chosen. The finite element model was built. Comparisons between the simulation results and the test re-sults prove that the simulation results are correct. The distribution of the temperature field of the chimney foundationwas analyzed. The multivariate linear regression of the hightest tomperature was performed on the inner wall of thechimney foundation by the numerical calculated results. The fitting property of the highest temperature with six influ-ence factors was obtained. A simple method for the calculation of the temperature field of the chimney foundation wasprovided.

C^1 natural element method for strain gradient linear elasticity and its application to microstructures

C^1 natural element method （C^1 NEM） is applied to strain gradient linear elasticity, and size effects on mi crostructures are analyzed. The shape functions in C^1 NEM are built upon the natural neighbor interpolation （NNI）, with interpolation realized to nodal function and nodal gradient values, so that the essential boundary conditions （EBCs） can be imposed directly in a Galerkin scheme for partial differential equations （PDEs）. In the present paper, C^1 NEM for strain gradient linear elasticity is constructed, and sev- eral typical examples which have analytical solutions are presented to illustrate the effectiveness of the constructed method. In its application to microstructures, the size effects of bending stiffness and stress concentration factor （SCF） are studied for microspeciem and microgripper, respectively. It is observed that the size effects become rather strong when the width of spring for microgripper, the radius of circular perforation and the long axis of elliptical perforation for microspeciem come close to the material characteristic length scales. For the U-shaped notch, the size effects decline obviously with increasing notch radius, and decline mildly with increasing length of notch.

Scheme for Concentration of Unknown Photonic GHZ Entangled States via Linear Optical Elements

We propose two schemes to concentrate unknown nonmaximally tripartite GHZ entangled states via linear optical elements. The finial maximally entangled states obtained from our schemes are shared by two or three parties. Our schemes only need polarizing beam splitters and single-photon detectors. In addition, the schemes can be demonstrated within current experimental technology.

Dynamic self-adaptive ANP algorithm and its application to electric field simulation of aluminum reduction cell

Region partition(RP) is the key technique to the finite element parallel computing(FEPC),and its performance has a decisive influence on the entire process of analysis and computation.The performance evaluation index of RP method for the three-dimensional finite element model(FEM) has been given.By taking the electric field of aluminum reduction cell(ARC) as the research object,the performance of two classical RP methods,which are Al-NASRA and NGUYEN partition(ANP) algorithm and the multi-level partition(MLP) method,has been analyzed and compared.The comparison results indicate a sound performance of ANP algorithm,but to large-scale models,the computing time of ANP algorithm increases notably.This is because the ANP algorithm determines only one node based on the minimum weight and just adds the elements connected to the node into the sub-region during each iteration.To obtain the satisfied speed and the precision,an improved dynamic self-adaptive ANP(DSA-ANP) algorithm has been proposed.With consideration of model scale,complexity and sub-RP stage,the improved algorithm adaptively determines the number of nodes and selects those nodes with small enough weight,and then dynamically adds these connected elements.The proposed algorithm has been applied to the finite element analysis(FEA) of the electric field simulation of ARC.Compared with the traditional ANP algorithm,the computational efficiency of the proposed algorithm has been shortened approximately from 260 s to 13 s.This proves the superiority of the improved algorithm on computing time performance.