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.

ANALYSIS OF THE THERMOPHYSICAL PARAMETERS OF MOIST WOOD PARTICLE MATERIAL IN A COUPLED HEAT AND MASS TRANSFER PROCESS OF FREEZING BY USING FINITE ELEMENT METHOD

The coupled heat and moisture transfer in a freezing process of wood particle material was mathematically modeled in the paper. The models were interactively solved by using the numerical method(the finite element method and the finite difference method). By matching the theoretical calculation to an experiment, the nonlinear problem was analyzed and the variable thermophysical parameters concerned was evaluated. The analysis procedure and the evaluation of the parameters were presented in detail. The result of the study showed that by using the method as described in the paper, it was possible to determine the variable (with respect to temperature, moisture content and freezing state) thermophysical parameters which were unknown or difficult to measure as long as the governing equations for a considered process were available. The method can significantly reduces the experiment efforts for determining thermophysical parameters which arc very complicated to measure. The determined variable of the effective heat conductivity of wood particle material was given in the paper. The error of the numerical calculation was also estimated by the comparison with a matched experiment.

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.

Circuit-field coupled finite element analysis method for an electromagnetic acoustic transducer under pulsed voltage excitation

This paper presents an analytical method for electromagnetic acoustic transducers （EMATs） under voltage excitation and considers the non-uniform distribution of the biased magnetic field. A complete model of EMATs including the non-uniform biased magnetic field, a pulsed eddy current field and the acoustic field is built up. The pulsed voltage excitation is transformed to the frequency domain by fast Fourier transformation （FFT）. In terms of the time harmonic field equations of the EMAT system, the impedances of the coils under different frequencies are calculated according to the circuit-field coupling method and Poynting＇s theorem. Then the currents under different frequencies are calculated according to Ohm＇s law and the pulsed current excitation is obtained by inverse fast Fourier transformation （IFFT）. Lastly, the sequentially coupled finite element method （FEM） is used to calculate the Lorentz force in the EMATs under the current excitation. An actual EMAT with a two-layer two-bundle printed circuit board （PCB） coil, a rectangular permanent magnet and an aluminium specimen is analysed. The coil impedances and the pulsed current are calculated and compared with the experimental results. Their agreement verified the validity of the proposed method. Furthermore, the influences of lift-off distances and the non-uniform static magnetic field on the Lorentz force under pulsed voltage excitation are studied.

Numerical analysis of coupled finite element with element-free Galerkin in sheet flexible-die forming

A numerical method for coupled deformation between sheet metal and flexible-die was proposed. Based on the updated Lagrangian （UL） formulation, the elastoplastic deformation of sheet metal was analyzed with finite element method （FEM） and the bulk deformation of flexible-die was analyzed with element-free Galerkin method （EFGM）. The frictional contact between sheet metal and flexible-die was treated by the penalty function method. The sheet elastic flexible-die bulging process was analyzed with the FEM-EFGM program for coupled deformation between sheet metal and bulk flexible-die, called CDSB-FEM-EFGM for short. Compared with finite element code DEFORM-2D and experiment results, the CDSB-FEM-EFGM program is feasible. This method provides a suitable numerical method to analyze sheet flexible-die forming.

RANDOM MICROSTRUCTURE FINITE ELEMENT METHOD FOR EFFECTIVE NONLINEAR PROPERTIES OF COMPOSITE MATERIALS

Some theoretical methods have been reported to deal with nonlinear problems of composite materials but the accuracy is not so good. In the meantime, a lot of linear problems are difficult to be managed by the theoretical methods. The present study aims to use the developed method, the random microstructure finite element method, to deal with these nonlinear problems. In this paper, the random microstructure finite element method is used to deal with all three kinds of nonlinear property problems of composite materials. The analyzed results suggest the influences of the nonlinear phenomena on the effective properties of composite materials are significant and the random microstructure finite element method is an effective tool to investigate the nonlinear problems.

A New Unified Stabilized Mixed Finite Element Method of the Stokes-Darcy Coupled Problem: Isotropic Discretization

In this paper, we develop an a-priori error analysis of a new unified mixed finite element method for the coupling of fluid flow with porous media flow in RN, N ∈ {2,3}, on isotropic meshes. Flows are governed by the Stokes and Darcy equations, respectively, and the corresponding transmission conditions are given by mass conservation, balance of normal forces, and the Beavers-Joseph-Saffman law. The approach utilizes a modification of the Darcy problem which allows us to apply a variant nonconforming Crouzeix-Raviart finite element to the whole coupled Stokes-Darcy problem. The well-posedness of the finite element scheme and its convergence analysis are derived. Finally, the numerical experiments are presented, which confirm the excellent stability and accuracy of our method.

Finite element method simulation for tensile process of sintered iron-base material

Different material properties leads to different metal fracture behaviors. Even if the powder material is composed of plastic metal, the fracture still does not show macroscopic plastic deformation characteristics if the material contains a large number of voids. Eight node isoparametric elastic plastic finite element method was used to simulate the tensile process of sintered powder material. By setting a number of voids in the analyzed metal cuboid, the initial density was taken into consideration. The material properties of the three dimensional solid for the tensile simulation were defined with reference to the known pure iron material parameters. The load displacement curves during elongation were obtained with a universal testing machine, and then the simulated curves were compared with the experimental results. The factors that cause the stress concentration and strength decrease were analyzed according to the simulated equivalent von Mises stress distribution.

Finite Element Analysis on Plane Stress Crack Growth in a Power-Law Hardening Material

Aim To determine numerically the field characteristics in the vied at the tip of a place crack growing steadily in a power-law hardening material. Meteods. Methods on the Euler mode and small-scale yield assumption, the numerical results were given by nonlinear finite element analysis. Results The numerical results of the shape of the active plastic sone, the angular distribution of stresseses and Clack tip opening displacement (CTOD) in the vicinity at the hp of the steadily groWing CraCk are determined. Conclusion The comparison between the numerical results given by the present wort and those given by analytic asymptotic analysis shows that the present work reached a very high accuracy.

A dynamic large-deformation particle finite element method for geotechnical applications based on Abaqus

In this paper,the application of Abaqus-based particle finite element method(PFEM)is extended from static to dynamic large deformation.The PFEM is based on periodic mesh regeneration with Delaunay triangulation to avoid mesh distortion.Additional mesh smoothing and boundary node smoothing techniques are incorporated to improve the mesh quality and solution accuracy.The field variables are mapped from the old to the new mesh using the closest point projection method to minimize the mapping error.The procedures of the proposed Abaqus-based dynamic PFEM(Abaqus-DPFEM)analysis and its implementation in Abaqus are detailed.The accuracy and robustness of the proposed approach are examined via four illustrative numerical examples.The numerical results show a satisfactory agreement with published results and further confirm the applicability of the Abaqus-DPFEM to solving dynamic large-deformation problems in geotechnical engineering.

Modeling and finite element analysis of transduction process of electromagnetic acoustic transducers for nonferromagnetic metal material testing

Facing the problems lack of considering the non-uniform distribution of the static bias magnetic field and computing the panicle displacements in the simulation model of electromagnetic acoustic transducer （EMAT）, a multi-field coupled model was established and the finite element method （FEM） was presented to calculate the entire transduction process. The multi-field coupled model included the static magnetic field, pulsed eddy current field and mechanical field. The FEM equations of the three fields were derived by Garlerkin FEM method. Thus, the entire transduction process of the EMAT was calculated through sequentially coupling the three fields. The transduction process of a Lamb wave EMAT was calculated according to the present model and method. The results show that, by the present method, it is valid to calculate the particle displacement under the given excitation signal and non-uniformly distributed static magnetic field. Calculation error will be brought about if the non-uniform distribution of the static bias magnetic field is neglected.

Finite Element Simulation of Flexible Roll Forming with Supplemented Material Data and the Experimental Verification

Flexible roll forming is a promising manufacturing method for the production of variable cross section products. Considering the large plastic strain in this forming process which is much larger than that of uniform deformation phase of uniaxial tensile test, the widely adopted method of simulating the forming processes with non-supplemented material data from uniaxial tensile test will certainly lead to large error. To reduce this error, the material data is supplemented based on three constitutive models. Then a finite element model of a six passes flexible roll forming process is established based on the supplemented material data and the original material data from the uniaxial tensile test. The flexible roll forming experiment of a B pillar reinforcing plate is carried out to verify the proposed method. Final cross section shapes of the experimental and the simulated results are compared. It is shown that the simulation calculated with supplemented material data based on Swift model agrees well with the experimental results, while the simulation based on original material data could not predict the actual deformation accurately. The results indicate that this material supplement method is reliable and indispensible, and the simulation model can well reflect the real metal forming process. Detailed analysis of the distribution and history of plastic strain at different positions are performed. A new material data supplement method is proposed to tackle the problem which is ignored in other roll forming simulations, and thus the forming process simulation accuracy can be greatly improved.

THE COUPLING OF BOUNDARY ELEMENT AND FINITE ELEMENT METHODS FOR THE EXTERIOR NONSTATIONARY NAVIER-STOKES EQUATIONS

In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variational formulation and well posedness of the coupling method are obtained. The convergence and optimal estimates for the approximation solution are provided.

Calculation of the Coupling Coefficient of Twin-Core Fiber Based on the Supermode Theory with Finite Element Method

Currently, coupled mode theory (CMT) is widely used for calculating the coupling coefficient of twin-core fibers (TCFs) that are used in a broad range of important applications. This approach is highly accurate for scenarios with weak coupling between the cores but shows significant errors in the strong coupling scenarios, necessitating the use of a more accurate method for coupling coefficient calculations. Therefore, in this work, we calculate the coupling coefficients of TCFs using the supermode theory with finite element method (FEM) that has higher accuracy than CMT, particularly for the strong coupling TCF. To investigate the origin of the differences between the results obtained by these two methods, the modal field distributions of the supermodes of TCF are simulated and analyzed in detail.

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 threedimensional 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.

NEW ALGORITHM OF COUPLING ELEMENT-FREE GALERKIN WITH FINITE ELEMENT METHOD

Through the construction of a new ramp function, the element-flee Galerkin method and finite element coupling method were applied to the whole field, and was made fit for the structure of element nodes within the interface regions, both satisfying the essential boundary conditions and deploying meshless nodes and finite elements in a convenient and flexible way, which can meet the requirements of computation for complicated field. The comparison between the results of the present study and the corresponding analytical solutions shows this method is feasible and effective.

An explicit finite element method for dynamic analysis in three-medium coupling system and its application

In this paper, an explicit finite element method to analyze the dynamic responses of three-medium coupled systems with any terrain is developed on the basis of the numerical simulation of the continuous conditions on the bounda-ries among fluid saturated porous medium, elastic single-phase medium and ideal fluid medium. This method is a very effective one with the characteristic of high calculating speed and small memory needed because the formulae for this explicit finite element method have the characteristic of decoupling, and which does not need to solve sys-tem of linear equations. The method is applied to analyze the dynamic response of a reservoir with considering the dynamic interactions among water, dam, sediment and basement rock. The vertical displacement at the top point of the dam is calculated and some conclusions are given.

FINITE ELEMENT SOLUTION TO COUPLED THERMO－ELASTICCONTACT STRESS AND IMPACTRESPONSE OF MESHING GEARS

A mixed finite element solution of contact stresses in meshing gears is investigated with the consideration of coupled thermo-elastic deformation and impact behavior. A simulation procedure of finite element solution of meshing gears is developed. The versatility of the procedure for both numerical accuracy and computational efficiency is verified by numerical analysis of meshing gear teeth.

Existence and numerical approximation of a solution to frictional contact problem for electro-elastic materials

In this paper,a frictional contact problem between an electro-elastic body and an electrically conductive foundation is studied.The contact is modeled by normal compliance with finite penetration and a version of Coulomb’s law of dry friction in which the coefficient of friction depends on the slip.In addition,the effects of the electrical conductivity of the foundation are taken into account.This model leads to a coupled system of the quasi-variational inequality of the elliptic type for the displacement and the nonlinear variational equation for the electric potential.The existence of a weak solution is proved by using an abstract result for elliptic variational inequalities and a fixed point argument.Then,a finite element approximation of the problem is presented.Under some regularity conditions,an optimal order error estimate of the approximate solution is derived.Finally,a successive iteration technique is used to solve the problem numerically and a convergence result is established.

Friction Characteristics Between Marine Clay and Construction Materials

Structure-soil interface friction characteristics is of importance to investigate the interaction between engineering structures and soils,especially for offshore structures.The interface friction behavior between marine clay and structural materials with different roughness was studied in this paper by using 3D optical scanning tests,a modified direct shear device and numerical simulation.Relationships between the surface roughness of structures,water content and interface friction angle were presented by model tests.The increase of water contents decreased the interface friction angles.For interfaces with different roughness,the interface friction angles will be smaller than that of the soil when the water content exceeds a certain value.The roughness of the interface and the water content of the soil are mutually coupled to influence the coefficient of friction(COF).This paper proposed a Finite Element Method(FEM)to simulate the interface direct shear tests of structures with different roughness.The surface models with different roughness are established based on the structure data obtained by 3D scanning.The Coupled Eulerian-Lagrangian(CEL)approach was employed to analyse soils sheared by irregular surfaces.The interface behavior for interfaces with different roughness under cyclic shear stresses was analyzed by FEM.