Propagations of Rayleigh and Love waves in ZnO films/glass substrates analyzed by three-dimensional finite element method

Propagation characteristics of surface acoustic waves（SAWs） in ZnO films/glass substrates are theoretically investigated by the three-dimensional（3D） finite element method. At first, for（11ˉ20） ZnO films/glass substrates, the simulation results confirm that the Rayleigh waves along the [0001] direction and Love waves along the [1ˉ100] direction are successfully excited in the multilayered structures. Next, the crystal orientations of the ZnO films are rotated, and the influences of ZnO films with different crystal orientations on SAW characterizations, including the phase velocity, electromechanical coupling coefficient, and temperature coefficient of frequency, are investigated. The results show that at appropriate h/λ, Rayleigh wave has a maximum k^2 of 2.4% in（90°, 56.5°, 0°） ZnO film/glass substrate structure; Love wave has a maximum k^2 of 3.81% in（56°, 90°, 0°） ZnO film/glass substrate structure. Meantime, for Rayleigh wave and Love wave devices, zero temperature coefficient of frequency（TCF） can be achieved at appropriate ratio of film thickness to SAW wavelength. These results show that SAW devices with higher k^2 or lower TCF can be fabricated by flexibly selecting the crystal orientations of ZnO films on glass substrates.

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.

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 Deep Learning Approach to Shape Optimization Problems for Flexoelectric Materials Using the Isogeometric Finite Element Method

A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization.

A New Isogeometric Finite Element Method for Analyzing Structures

High-performance finite element research has always been a major focus of finite element method studies.This article introduces isogeometric analysis into the finite element method and proposes a new isogeometric finite element method.Firstly,the physical field is approximated by uniform B-spline interpolation,while geometry is represented by non-uniform rational B-spline interpolation.By introducing a transformation matrix,elements of types C^(0)and C^(1)are constructed in the isogeometric finite element method.Subsequently,the corresponding calculation formats for one-dimensional bars,beams,and two-dimensional linear elasticity in the isogeometric finite element method are derived through variational principles and parameter mapping.The proposed method combines element construction techniques of the finite element method with geometric construction techniques of isogeometric analysis,eliminating the need for mesh generation and maintaining flexibility in element construction.Two elements with interpolation characteristics are constructed in the method so that boundary conditions and connections between elements can be processed like the finite element method.Finally,the test results of several examples show that:(1)Under the same degree and element node numbers,the constructed elements are almost consistent with the results obtained by traditional finite element method;(2)For bar problems with large local field variations and beam problems with variable cross-sections,high-degree and multi-nodes elements constructed can achieve high computational accuracy with fewer degrees of freedom than finite element method;(3)The computational efficiency of isogeometric finite element method is higher than finite element method under similar degrees of freedom,while as degrees of freedom increase,the computational efficiency between the two is similar.

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.

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.

A Computational Framework for Parachute Inflation Based on Immersed Boundary/Finite Element Approach

A computational framework for parachute inflation is developed based on the immersed boundary/finite element approach within the open-source IBAMR library.The fluid motion is solved by Peskin's diffuse-interface immersed boundary(IB)method,which is attractive for simulating moving-boundary flows with large deformations.The adaptive mesh refinement technique is employed to reduce the computational cost while retain the desired resolution.The dynamic response of the parachute is solved with the finite element approach.The canopy and cables of the parachute system are modeled with the hyperelastic material.A tether force is introduced to impose rigidity constraints for the parachute system.The accuracy and reliability of the present framework is validated by simulating inflation of a constrained square plate.Application of the present framework on several canonical cases further demonstrates its versatility for simulation of parachute inflation.

Modularized and Parametric Modeling Technology for Finite Element Simulations of Underground Engineering under Complicated Geological Conditions

The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.

Probabilistic Analysis of Slope Using Finite Element Approach and Limit Equilibrium Approach around Amalpata Landslide of West Central, Nepal

The stability study of the ongoing and recurring Amalpata landslide in Baglung in Nepal’s Gandaki Province is presented in this research. The impacted slope is around 200 meters high, with two terraces that have different slope inclinations. The lower bench, located above the basement, consistently fails and sets others up for failure. The fluctuating water level of the slope, which travels down the slope masses, exacerbates the slide problem. The majority of these rocks are Amalpata landslide area experiences several structural disruptions. The area’s stability must be evaluated in order to prevent and control more harm from occurring to the nearby agricultural land and people living along the slope. The slopes’ failures increase the damages of house existing in nearby area and the erosion of the slope. Two modeling techniques the finite element approach and the limit equilibrium method were used to simulate the slope. The findings show that, in every case, the terrace above the basement is where the majority of the stress is concentrated, with a safety factor of near unity. Using probabilistic slope stability analysis, the failure probability was predicted to be between 98.90% and 100%.

Three-dimensional forward modeling for magnetotelluric sounding by finite element method

A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.

Alternating Direction Finite Volume Element Methods for Three-Dimensional Parabolic Equations

This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite difference scheme with second-order splitting error, the other two schemes have third-order splitting error, and the last one is an extended LOD scheme. The L2 norm and H1 semi-norm error estimates are obtained for the first scheme and second one, respectively. Finally, two numerical examples are provided to illustrate the efficiency and accuracy of the methods.

Simulation of three-dimensional tension-induced cracks based on cracking potential function-incorporated extended finite element method

In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination of the cracking direction constitutes a great challenge.In most cases,the local stress state provides the fundamental criterion to judge the presence of cracks and the direction of crack propagation.However,in the case of three-dimensional analysis,the coordination relationship between grid elements due to occurrence of cracks becomes a difficult problem for this method.In this paper,based on the extended finite element method,the stress-related function field is introduced into the calculation domain,and then the boundary value problem of the function is solved.Subsequently,the envelope surface of all propagation directions can be obtained at one time.At last,the possible surface can be selected as the direction of crack development.Based on the aforementioned procedure,such method greatly reduces the programming complexity of tracking the crack propagation.As a suitable method for simulating tension-induced failure,it can simulate multiple cracks simultaneously.

STUDY ON THREE-DIMENSIONAL FINITE BODIES CONTAINING CRACKS USING THE FINITE ELEMENT METHOD OF LINES

The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good accuracy can be obtained with relatively coarse girds. In particular, application to the tension specimen shows very good agreement with the evaluation of stress intensity factors, which is better than the results of other methods. This implies a considerable potential for using this method in the 3D analysis of finite geometry solids and suggests a possible extension of this technique to nonlinear material behavior.

Three-Dimensional Thermo-Elastic-Plastic Finite Element Method Modeling for Predicting Weld-Induced Residual Stresses and Distortions in Steel Stiffened-Plate Structures

The objective of the present paper is to develop nonlinear finite element method models for predicting the weld-induced initial deflection and residual stress of plating in steel stiffened-plate structures. For this purpose, three-dimensional thermo-elastic-plastic finite element method computations are performed with varying plate thickness and weld bead length (leg length) in welded plate panels, the latter being associated with weld heat input. The finite element models are verified by a comparison with experimental database which was obtained by the authors in separate studies with full scale measurements. It is concluded that the nonlinear finite element method models developed in the present paper are very accurate in terms of predicting the weld-induced initial imperfections of steel stiffened plate structures. Details of the numerical computations together with test database are documented.

Three-dimensional analysis of elastic stress distribution of indented ceramic surface by finite element method

The three-dimensional stress distributions in the area surrounding indentation pattern for three different materials, Al2O3, Si3N4 and SiC were analyzed by finite element method(FEM). Those theoretical results were also compared with the experimental ones by Rockwell hardness test. The effect of loading stress on the plastic deformation in specimens, surface was investigated on the assumption of shear strain energy theory by Huber-Mises when the materials were indented. The distributions of nomal stress, shear stress, and Mises stress were analysed with variations of loading conditions. It is clear that the analytical results for the stress distributions, the crack length and its density of probability are in good agreement with the experimental results.

CHEBYSHEV PSEUDOSPECTRAL-HYBRID FINITE ELEMENT METHOD FOR THREE-DIMENSIONAL VORTICITY EQUATION

In this paper,Chebyshev pseudospectral-finite element schemes are proposed for solving three dimensional vorticity equation.Some approximation results in nonisotropic Sobolev spaces are given.The generalized stability and the convergence are proved strictly.The numerical results show the advantages of this method.The technique in this paper is also applicable to other three-dimensional nonlinear problems in fluid dynamics.