A stable implicit nodal integration-based particle finite element method(N-PFEM)for modelling saturated soil dynamics

In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.

Crystal plasticity finite element simulations on extruded Mg-10Gd rod with texture gradient

The mechanical properties of an extruded Mg-10Gd sample, specifically designed for vascular stents, are crucial for predicting its behavior under service conditions. Achieving homogeneous stresses in the hoop direction, essential for characterizing vascular stents, poses challenges in experimental testing based on standard specimens featuring a reduced cross section. This study utilizes an elasto-visco-plastic self-consistent polycrystal model(ΔEVPSC) with the predominant twinning reorientation(PTR) scheme as a numerical tool, offering an alternative to mechanical testing. For verification, various mechanical experiments, such as uniaxial tension, compression, notched-bar tension, three-point bending, and C-ring compression tests, were conducted. The resulting force vs. displacement curves and textures were then compared with those based on the ΔEVPSC model. The computational model's significance is highlighted by simulation results demonstrating that the differential hardening along with a weak strength differential effect observed in the Mg-10Gd sample is a result of the interplay between micromechanical deformation mechanisms and deformation-induced texture evolution. Furthermore, the study highlights that incorporating the axisymmetric texture from the as-received material incorporating the measured texture gradient significantly improves predictive accuracy on the strength in the hoop direction. Ultimately, the findings suggest that the ΔEVPSC model can effectively predict the mechanical behavior resulting from loading scenarios that are impossible to realize experimentally, emphasizing its valuable contribution as a digital twin.

THE SUPERCLOSENESS OF THE FINITE ELEMENT METHOD FOR A SINGULARLY PERTURBED CONVECTION-DIFFUSION PROBLEM ON A BAKHVALOV-TYPE MESH IN 2D

For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of the mesh in the layer adjacent to the transition point,resulting in a suboptimal estimate for convergence.Existing analysis techniques cannot handle these difficulties well.To fill this gap,here a novel interpolation is designed delicately for the smooth part of the solution,bringing about the optimal supercloseness result of almost order 2 under an energy norm for the finite element method.Our theoretical result is uniform in the singular perturbation parameterεand is supported by the numerical experiments.

A Wave Superposition-Finite Element Method for Calculating the Radiated Noise Generated by Volumetric Targets in Shallow Water

A combined method of wave superposition and finite element is proposed to solve the radiation noise of targets in shallow sea.Taking the sound propagation of spherical sound source in shallow sea as an example,the radiation sound field of the spherical sound source is equivalent to the linear superposition of the radiation sound field of several internal point sound sources,and then the radiated noise induced by spherical sound source can be predicted quickly.The accuracy and efficiency of the method are verified by comparing with the numerical results of finite element method,and the rapid prediction of underwater radiated noise of cylindrical shell is carried out based on the method.The results show that compared with the finite element method,the relative error of the calculation results under different simulation conditions does not exceed 0.1%,and the calculation time is about 1/10 of the finite element method,so this method can be used to solve the radiated noise of shallow underwater targets.

Multi-scale Modeling and Finite Element Analyses of Thermal Conductivity of 3D C/SiC Composites Fabricating by Flexible-Oriented Woven Process

Thermal conductivity is one of the most significant criterion of three-dimensional carbon fiber-reinforced SiC matrix composites(3D C/SiC).Represent volume element(RVE)models of microscale,void/matrix and mesoscale proposed in this work are used to simulate the thermal conductivity behaviors of the 3D C/SiC composites.An entirely new process is introduced to weave the preform with three-dimensional orthogonal architecture.The 3D steady-state analysis step is created for assessing the thermal conductivity behaviors of the composites by applying periodic temperature boundary conditions.Three RVE models of cuboid,hexagonal and fiber random distribution are respectively developed to comparatively study the influence of fiber package pattern on the thermal conductivities at the microscale.Besides,the effect of void morphology on the thermal conductivity of the matrix is analyzed by the void/matrix models.The prediction results at the mesoscale correspond closely to the experimental values.The effect of the porosities and fiber volume fractions on the thermal conductivities is also taken into consideration.The multi-scale models mentioned in this paper can be used to predict the thermal conductivity behaviors of other composites with complex structures.

Extended finite element-based cohesive zone method for modeling simultaneous hydraulic fracture height growth in layered reservoirs

In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.

Finite Element Analysis for Magneto-Convection Heat Transfer Performance in Vertical Wavy Surface Enclosure:Fin Size Impact

The goal of this paper is to represent a numerical study of magnetohydrodynamic mixed convection heat transfer in a lid-driven vertical wavy enclosure with a fin attached to the bottomwall.We use a finite elementmethod based on Galerkin weighted residual(GWR)techniques to set up the appropriate governing equations for the present flow model.We have conducted a parametric investigation to examine the impact of Hartmann and Richardson numbers on the flow pattern and heat transmission features inside a wavy cavity.We graphically represent the numerical results,such as isotherms,streamlines,velocity profiles,local and mean Nusselt numbers,and average surface temperature.Comparisons between the results of this work and previously published work in a literature review have been produced to examine the reliability and consistency of the data.The different sizes of the fin surface significantly impact flow creation and temperature fields.Additionally,the long fin size is necessary to enhance the heat transfer rate on the right surface at large Richardson numbers and low Hartmann numbers.Fin surfaces can significantly increase the mixing of fluid inside the enclosure,which can mean reductions in reaction times and operating costs,along with increases in heat transfer and efficiency.

Collapse Behavior of Pipe-Framed Greenhouses with and without Reinforcement under Snow Loading:A 3-D Finite Element Analysis

The present paper first investigates the collapse behavior of a conventional pipe-framed greenhouse under snow loading based on a 3-D finite element analysis,in which both geometrical and material non-linearities are considered.Three snow load distribution patterns related to the wind-driven snow particle movement are used in the analysis.It is found that snow load distribution affects the deformation and collapse behavior of the pipe-framed greenhouse significantly.The results obtained in this study are consistent with the actual damage observed.Next,discussion is made of the effects of reinforcements by adding members to the basic frame on the strength of the whole structure,in which seven kinds of reinforcement methods are examined.A buckling analysis is also carried out.The results indicate that the most effective reinforcement method depends on the snow load distribution pattern.

A Full Predictor-Corrector Finite Element Method for the One-Dimensional Heat Equation with Time-Dependent Singularities

The energy norm convergence rate of the finite element solution of the heat equation is reduced by the time-regularity of the exact solution. This paper presents an adaptive finite element treatment of time-dependent singularities on the one-dimensional heat equation. The method is based on a Fourier decomposition of the solution and an extraction formula of the coefficients of the singularities coupled with a predictor-corrector algorithm. The method recovers the optimal convergence rate of the finite element method on a quasi-uniform mesh refinement. Numerical results are carried out to show the efficiency of the method.

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.

Hermite Finite Element Method for Vibration Problem of Euler-Bernoulli Beam on Viscoelastic Pasternak Foundation

Viscoelastic foundation plays a very important role in civil engineering. It can effectively disperse the structural load into the foundation soil and avoid the damage caused by the concentrated load. The model of Euler-Bernoulli beam on viscoelastic Pasternak foundation can be used to analyze the deformation and response of buildings under complex geological conditions. In this paper, we use Hermite finite element method to get the numerical approximation scheme for the vibration equation of viscoelastic Pasternak foundation beam. Convergence and error estimation are rigourously established. We prove that the fully discrete scheme has convergence order O(τ2+h4), where τis time step size and his space step size. Finally, we give four numerical examples to verify the validity of theoretical analysis.

Gradient Recovery Based Two-Grid Finite Element Method for Parabolic Integro-Differential Optimal Control Problems

In this paper, the optimal control problem of parabolic integro-differential equations is solved by gradient recovery based two-grid finite element method. Piecewise linear functions are used to approximate state and co-state variables, and piecewise constant function is used to approximate control variables. Generally, the optimal conditions for the problem are solved iteratively until the control variable reaches error tolerance. In order to calculate all the variables individually and parallelly, we introduce a gradient recovery based two-grid method. First, we solve the small scaled optimal control problem on coarse grids. Next, we use the gradient recovery technique to recover the gradients of state and co-state variables. Finally, using the recovered variables, we solve the large scaled optimal control problem for all variables independently. Moreover, we estimate priori error for the proposed scheme, and use an example to validate the theoretical results.

Finite element analysis of stress at implant-bone interface of dental implants with different structures

The effect of structure,elastic modulus and thickness of lower modulus layer in porous titanium implants on the stress distribution at the implant-bone interface was investigated.Three-dimensional finite element models of different titanium implants were constructed.The structures of the implants included the whole lower modulus style （No.1）,bio-mimetic style （No.2）,the whole lower modulus style in cancellous bone （No.3） and the whole dense style No.4.The stress distributions at bone-implant interface under static loading were analyzed using Ansys Workbench 10.0 software.The results indicated that the distribution of interface stress is strongly depended on the structure of the implants.The maximum stresses in cancellous bone and root region of implant No.2 are lower than those in the other three implants.A decrease in the modulus of the low modulus layer facilitates the interface stress transferring.Increasing the thickness of the low modulus layer can reduce the stress and induce a more uniform stress distribution at the interface.Among the four implants,biomimetic style implant No.2 is superior in transferring implant-bone interface stress to surrounding bones.

NUMERICAL SIMULATION OF UNSTEADY-STATE UNDEREXPANDED JET USING DISCONTINUOUS GALERKIN FINITE ELEMENT METHOD

A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.

EXPERIMENTAL RESEARCH AND NONLINEAR FINITE ELEMENT ANALYSIS ON NEW TYPE JOINT BETWEEN COLUMN AND STEEL BEAM OF CONCRETE-FILLED RECTANGULAR STEEL TUBULAR

Experimental results of new type joints between the column and the. steel beam of concrete-filled rectangular steel tubular （CFRT） under reversed cyclic loads are presented. The earthquake resistant capacity of the joint is influenced by infilled concrete, stiffener length and relative dimensions of column and beam. It is found that the hysteresis curves obtained in the experiment are full and the joints have a good energy dissipation capacity. The nonlinear finite element models are also used to analyze the hysteresis behavior of the joints under reversed cyclic loads using ANSYS 8.0. The influences of the stiffener length and the infilled concrete are analyzed. Analytical results show that the stiffener length and the infilled concrete are critical for the joints. Furthermore, the skeleton curves of the finite element models are in good agreement with those of experiments.

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.

Phase-field modeling of dendritic growth under forced flow based on adaptive finite element method

A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.

Dynamic finite element model updating using meta-model and genetic algorithm

Current dynamic finite element model updating methods are not efficient or restricted to the problem of local optima. To circumvent these, a novel updating method which integrates the meta-model and the genetic algorithm is proposed. Experimental design technique is used to determine the best sampling points for the estimation of polynomial coefficients given the order and the number of independent variables. Finite element analyses are performed to generate the sampling data. Regression analysis is then used to estimate the response surface model to approximate the functional relationship between response features and design parameters on the entire design space. In the fitness evaluation of the genetic algorithm, the response surface model is used to substitute the finite element model to output features with given design parameters for the computation of fitness for the individual. Finally, the global optima that corresponds to the updated design parameter is acquired after several generations of evolution. In the application example, finite element analysis and modal testing are performed on a real chassis model. The finite element model is updated using the proposed method. After updating, root-mean-square error of modal frequencies is smaller than 2%. Furthermore, prediction ability of the updated model is validated using the testing results of the modified structure. The root-mean-square error of the prediction errors is smaller than 2%.

Analysis of the Elastic-Plastic Dynamic Finite Element on a Jointed Rock Mass

Aim To study the elastic plastic dynamical constitutive relations about a jointed rock mass under explosion load and its computer simulation. Methods\ Stress history is taken into account and stresses will follow changes in time during a period of explosion load. According to the principle of static force balance, the corresponding nodal concentrated force is calculated and the nodal displacement is counted. The elastic plastic dynamic finite element equations are thus obtained. Results\ A finite element method is given for a jointed rock mass under explosion load. Conclusion\ The problem of large plastic deformation for jointed rock mass on blasting was efficiently resolved through dynamic finite element analysis and the range of damages by blasting simulated, and this pushes forward the problem to engineering practice.

FINITE ELEMENT-ARTIFICIAL TRANSMITTING BOUNDARY METHOD FOR WAVE SCATTERING FROM IRREGULAR CYLINDER

The finite element artificial transmitting boundary method is employed here to treat the near field scattering of a cylindrical wave from an irregular cylinder. A comparison is made between this method and the analytical one. And then examples are given to demonstrate the solution of several problems of the irregular object scattering. The method can not only produce clear physical pictures, but can efficiently handle many complicated scattering problems.