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

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.

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

ASYMPTOTICAL STABILITY OF NEUTRAL REACTION-DIFFUSION EQUATIONS WITH PCAS AND THEIR FINITE ELEMENT METHODS

This paper focuses on the analytical and numerical asymptotical stability of neutral reaction-diffusion equations with piecewise continuous arguments.First,for the analytical solutions of the equations,we derive their expressions and asymptotical stability criteria.Second,for the semi-discrete and one-parameter fully-discrete finite element methods solving the above equations,we work out the sufficient conditions for assuring that the finite element solutions are asymptotically stable.Finally,with a typical example with numerical experiments,we illustrate the applicability of the obtained theoretical results.

A symplectic finite element method based on Galerkin discretization for solving linear systems

We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is the dissipative algorithm and cannot maintain long-term energy conservation.Thus,a symplectic finite element method with energy conservation is constructed in this paper.A linear elastic system can be discretized into multiple elements,and a Hamiltonian system of each element can be constructed.The single element is discretized by the Galerkin method,and then the Hamiltonian system is constructed into the Birkhoffian system.Finally,all the elements are combined to obtain the vibration equation of the continuous system and solved by the symplectic difference scheme.Through the numerical experiments of the vibration response of the Bernoulli-Euler beam and composite plate,it is found that the vibration response solution and energy obtained with the algorithm are superior to those of the Runge-Kutta algorithm.The results show that the symplectic finite element method can keep energy conservation for a long time and has higher stability in solving the dynamic responses of linear elastic systems.

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 Simulation of High Frequency Electromagnetics Wave Propagation on Vivaldi Antenna Using Finite Element Method

The simulation of the electromagnetic wave propagation plays an important role in predicting the performance of wireless transmission and communication systems. This research paper performs a numerical simulation using the finite element method (FEM) to study electromagnetic propagation through both conductive and dielectric media. The simulations are made using the COMSOL Multiphysics software which notably implements the finite element method. The microwave is produced by a Vivaldi antenna at the respective frequencies of 2.6 and 5 GHz and the propagation equation is formulated from Maxwell’s equations. The results obtained show that in the air, strong electric fields are observed in the slot and the micro-strip line for the two frequencies, they are even greater when the wave propagates in the glass and very weak for the copper. The 3D evolutions of the wave in air and glass present comparable values at equal frequencies, the curves being more regular in air (dielectric). The radiation patterns produced for air and glass are directional, with a large main lobe, which is narrower at 5 GHz. For copper, the wave propagation is quite uniform in space, and the radiation patterns show two main lobes with a much larger size at 2.6 GHz than at 5 GHz. The propagation medium would therefore influence the range of values of the gain of the antenna.

Modal and Thermal Analysis of a Modified Connecting Rod of an Internal Combustion Engine Using Finite Element Method

The connecting rod is one of the most important moving components in an internal combustion engine. The present work determined the possibility of using aluminium alloy 7075 material to design and manufacture a connecting rod for weight optimisation without losing the strength of the connecting rod. It considered modal and thermal analyses to investigate the suitability of the material for connecting rod design. The parameters that were considered under the modal analysis were: total deformation, and natural frequency, while the thermal analysis looked at the temperature distribution, total heat flux and directional heat flux of the four connecting rods made with titanium alloy, grey cast iron, structural steel and aluminium 7075 alloy respectively. The connecting rod was modelled using Autodesk inventor2017 software using the calculated parameters. The steady-state thermal analysis was used to determine the induced heat flux and directional heat flux. The study found that Aluminium 7075 alloy deformed more than the remaining three other materials but has superior qualities in terms of vibrational natural frequency, total heat flux and lightweight compared to structural steel, grey cast iron and titanium alloy.

A Posteriori Error Estimate of Two Grid Mixed Finite Element Methods for Semilinear Elliptic Equations

In this paper, we present the a posteriori error estimate of two-grid mixed finite element methods by averaging techniques for semilinear elliptic equations. We first propose the two-grid algorithms to linearize the mixed method equations. Then, the averaging technique is used to construct the a posteriori error estimates of the two-grid mixed finite element method and theoretical analysis are given for the error estimators. Finally, we give some numerical examples to verify the reliability and efficiency of the a posteriori error estimator.

PARTITION OF UNITY FINITE ELEMENT METHOD FOR SHORT WAVE PROPAGATION IN SOLIDS

A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.

Galerkin-based quasi-smooth manifold element(QSME)method for anisotropic heat conduction problems in composites with complex geometry

The accurate and efficient analysis of anisotropic heat conduction problems in complex composites is crucial for structural design and performance evaluation. Traditional numerical methods, such as the finite element method(FEM), often face a trade-off between calculation accuracy and efficiency. In this paper, we propose a quasi-smooth manifold element(QSME) method to address this challenge, and provide the accurate and efficient analysis of two-dimensional(2D) anisotropic heat conduction problems in composites with complex geometry. The QSME approach achieves high calculation precision by a high-order local approximation that ensures the first-order derivative continuity.The results demonstrate that the QSME method is robust and stable, offering both high accuracy and efficiency in the heat conduction analysis. With the same degrees of freedom(DOFs), the QSME method can achieve at least an order of magnitude higher calculation accuracy than the traditional FEM. Additionally, under the same level of calculation error, the QSME method requires 10 times fewer DOFs than the traditional FEM. The versatility of the proposed QSME method extends beyond anisotropic heat conduction problems in complex composites. The proposed QSME method can also be applied to other problems, including fluid flows, mechanical analyses, and other multi-field coupled problems, providing accurate and efficient numerical simulations.

Finite Element Simulation of Radial Tire Building and Shaping Processes Using an Elasto-Viscoplastic Model

The comprehensive tire building and shaping processes are investigated through the finite element method(FEM)in this article.The mechanical properties of the uncured rubber from different tire components are investigated through cyclic loading-unloading experiments under different strain rates.Based on the experiments,an elastoviscoplastic constitutive model is adopted to describe themechanical behaviors of the uncured rubber.The distinct mechanical properties,including the stress level,hysteresis and residual strain,of the uncured rubber can all be well characterized.The whole tire building process(including component winding,rubber bladder inflation,component stitching and carcass band folding-back)and the shaping process are simulated using this constitutive model.The simulated green tire profile is in good agreement with the actual profile obtained through 3D scanning.The deformation and stress of the rubber components and the cord reinforcements during production can be obtained fromthe FE simulation,which is helpful for judging the rationality of the tire construction design.Finally,the influence of the parameter“drum width”is investigated,and the simulated result is found to be consistent with the experimental observations,which verifies the effectiveness of the simulation.The established simulation strategy provides some guiding significance for the improvement of tire design parameters and the elimination of tire production defects.

Study of Axisymmetric Infinite Guide Lined with Locally Reacting Material without Flow Using DtN Operators

The present work proposed a new method for the modeling by the finite element method of the acoustic propagation problems in infinite axisymmetric cylindrical guides lined with locally reacting absorbent materials without flow. The method deals with the development of an efficient transparent boundary condition based on DtN operators. The method developed in this study is successfully applied to a straight axisymmetric lined guide by imposing a mode on one of the artificial boundaries of the truncated guide. The results are in good agreement with analytical solutions. Applying the method for a non-uniform axisymmetric lined guide which is a complex case, proved its effectiveness and the results compared to those of PML layers are in very good agreement.

A multithreaded parallel upwind sweep algorithm for the S_(N) transport equations discretized with discontinuous finite elements

The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can overcome the deficiencies of conventionally structured meshes in complex geometry modeling. A multithreaded parallel upwind sweep algorithm for S_(N) transport was proposed to achieve a more accurate geometric description and improve the computational efficiency. The spatial variables were discretized using the standard discontinuous Galerkin finite-element method. The angular flux transmission between neighboring meshes was handled using an upwind scheme. In addition, a combination of a mesh transport sweep and angular iterations was realized using a multithreaded parallel technique. The algorithm was implemented in the 2D/3D S_(N) transport code ThorSNIPE, and numerical evaluations were conducted using three typical benchmark problems:IAEA, Kobayashi-3i, and VENUS-3. These numerical results indicate that the multithreaded parallel upwind sweep algorithm can achieve high computational efficiency. ThorSNIPE, with a multithreaded parallel upwind sweep algorithm, has good reliability, stability, and high efficiency, making it suitable for complex shielding calculations.

Finite Element Orthogonal Collocation Approach for Time Fractional Telegraph Equation with Mamadu-Njoseh Polynomials

Finite element method (FEM) is an efficient numerical tool for the solution of partial differential equations (PDEs). It is one of the most general methods when compared to other numerical techniques. PDEs posed in a variational form over a given space, say a Hilbert space, are better numerically handled with the FEM. The FEM algorithm is used in various applications which includes fluid flow, heat transfer, acoustics, structural mechanics and dynamics, electric and magnetic field, etc. Thus, in this paper, the Finite Element Orthogonal Collocation Approach (FEOCA) is established for the approximate solution of Time Fractional Telegraph Equation (TFTE) with Mamadu-Njoseh polynomials as grid points corresponding to new basis functions constructed in the finite element space. The FEOCA is an elegant mixture of the Finite Element Method (FEM) and the Orthogonal Collocation Method (OCM). Two numerical examples are experimented on to verify the accuracy and rate of convergence of the method as compared with the theoretical results, and other methods in literature.

Seismic Wavelet Analysis Based on Finite Element Numerical Simulation

The practice of exploration and production has proved that explosives are excited in different surrounding rocks and the seismic wavelets collected have different characteristics. In this paper, by establishing a numerical model of the explosion in the well, using finite element analysis technology for numerical simulation, the simulation calculated the stress structure in the near-source area of the earthquake excitation, and extracted the seismic wavelet. The results show that the simulation seismic wavelet characteristics of different thin interbedded sand and mudstone structures have changed significantly. Through excitation simulation, the amplitude and spectrum information of seismic wavelets can be compared and analyzed, and the excitation parameters can be optimized. .

Finite Element Analysis of the Influence of Artificial Cementation on the Strength Parameters and Bearing Capacity of Sandy Soil under a Strip Footing

Artificial cementation is a method commonly used to enhance and improve soil properties. This paper investigates the effect of using different amounts of cement on soil strength parameters and soil bearing capacity, using the finite element method. Experimental tests are conducted on soil samples with different amounts of Portland cement. A 2-D numerical model is created and validated using the numerical modelling software, COMSOL Multiphysics 5.6 software. The study finds that the cohesion, and the angle of the internal friction of the soil samples increase significantly as a result of adding 1%, 2%, and 4% of Portland cement. The results demonstrate that the stresses and strain under the strip footing proposed decrease by 3.24% and 7.42%. Moreover, the maximum displacement also decreases by 1.47% and 2.97%, as a result of adding cements of 2% and 4%. The bearing capacity values obtained are therefore excellent, especially when using the 2% and 4% cement. The increase identified is due to the increased values of the bearing capacity factors. It is concluded that from an economic viewpoint, using 2% cement is the best option.

Error estimates of H^1-Galerkin mixed finite element method for Schrdinger equation

An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same time, optimal error estimates are derived for fully discrete schemes. And it is showed that the H1-Galerkin mixed finite element approximations have the same rate of convergence as in the classical mixed finite element methods without requiring the LBB consistency condition.