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.

SPACE-TIME FINITE ELEMENT METHOD FOR SCHRDINGER EQUATION AND ITS CONSERVATION

Energy conservation of nonlinear Schrodinger ordinary differential equation was proved through using continuous finite element methods of ordinary differential equation; Energy integration conservation was proved through using space-time continuous fully discrete finite element methods and the electron nearly conservation with higher order error was obtained through using time discontinuous only space continuous finite element methods of nonlinear Schrodinger partial equation. The numerical results are in accordance with the theory.

Structural Integrity Analysis of Containers Lost at Sea Using Finite Element Method

Unlike traditional transportation,container transportation is a relatively new logistics transportation mode.Shipping containers lost at sea have raised safety concerns.In this study,finite element analysis of containers subjected to hydrostatic pressure,using commercial software ANSYS APDL was performed.A computer model that can reasonably predict the state of an ISO cargo shipping container was developed.The von Mises stress distribution of the container was determined and the yield strength was adopted as the failure criterion.Numerical investigations showed that the conventional ship container cannot withstand hydrostatic pressure in deep water conditions.A strengthened container option was considered for the container to retain its structural integrity in water conditions.

Design and Development of Small Modular Courses for“Education-Training Integration”in Vocational Colleges:A Case Study

Vocational education plays a vital role in the development of skilled technical professionals and the advancement of the economy.However,the emphasis on campus education often neglects the importance of external training,hindering the overall development of vocational education.This study aims to address this issue by exploring the design and development of small modular courses that integrate training and education in vocational colleges,focusing on the mechanics course as a case study.The research methods employed in this study include an in-depth analysis of enterprise training needs,the development of digital teaching resources utilizing the finite element method(FEM),and the implementation of a small modular course integrating education and training.The data analysis reveals positive outcomes in terms of learners’comprehension and engagement with complex mechanics formulas through the use of stress nephograms and other digital resources.This research provides a new perspective on curriculum design and offers insights into the integration of training and education in vocational colleges.The findings underscore the significance of incorporating innovative teaching methodologies and digital resources in enhancing the quality and relevance of vocational education,ultimately contributing to the cultivation of skilled professionals and the growth of the vocational education sector.

Evaluation of mixed-mode stress intensity factors by extended finite element method

Extended finite element method （XFEM） implementation of the interaction integral methodology for evaluating the stress intensity factors （SIF） of the mixed-mode crack problem is presented. A discontinuous function and the near-tip asymptotic function are added to the classic finite element approximation to model the crack behavior. Two-state integral by the superposition of actual and auxiliary fields is derived to calculate the SIFs. Applications of the proposed technique to the inclined centre crack plate with inclined angle from 0° to 90° and slant edge crack plate with slant angle 45°, 67.5° and 90° are presented, and comparisons are made with closed form solutions. The results show that the proposed method is convenient, accurate and computationallv efficient.

CONSERVATION LAWS IN FINITE MICROCRACKING BRITTLE SOLIDS

This paper addresses the conservation laws in finite brittle solids with microcracks. The discussion is limited to the 2-D cases. First, after considering the combination of the Pseudo-Traction Method and the indirect Boundary Element Method, a versatile method for solving multi-crack interacting problems in finite plane solids is proposed, by which the fracture parameters （SIF and path-independent integrals） can be calculated with a desirable accuracy. Second, with the aid of the method proposed, the roles the conservation laws play in the fracture analysis for finite microcracking solids are studied. It is concluded that the conservation laws do play important roles in not only the fracture analysis but also the analysis of damage and stability for the finite microcracking system. Finally, the physical interpretation of the M-integral is discussed further. An explicit relation between the M-integral and the crack face area, i.e., M = GS, has been discovered using the analytical method, which can shed some light on the Damage Mechanics issues from a different perspective.

Continuous finite element methods for Hamiltonian systems

By applying the continuous finite element methods of ordinary differential equations, the linear element methods are proved having second-order pseudo-symplectic scheme and the quadratic element methods are proved having third-order pseudo- symplectic scheme respectively for general Hamiltonian systems, and they both keep energy conservative. The finite element methods are proved to be symplectic as well as energy conservative for linear Hamiltonian systems. The numerical results are in agree-ment with theory.

Finite Element Methods for Coupled Stokes and Darcy Problems

We derived and analyzed a new numerical scheme for the coupled Stokes and Darcy problems by using H（div） conforming elements in the entire domain. The approach employs the mixed finite element method for the Darcy equations and a stabilized H（div） finite element method for the Stokes equations. Optimal error estimates for the fluid velocity and pressure are derived. The finite element solutions from the new scheme not only feature a full satisfaction of the continuity equation, which is highly demanded in scientific computing, but also satisfy the mass conservation.

ELEMENT－BY－ELEMENT MATRIX DECOMPOSITION ANDSTEP－BY－STEP INTEGRATION METHOD FOR TRANSIENTDYNAMIC PROBLEMS

In this paper a general matrix decomposition scheme as well as an element-by-clement relaxation algorithm combined with step-by -step integration method is presented for transient dynamic problems thus the finite element method can be fromforming global stiffness matrix global mass matrix as well as solyin large scale sparse equations Theory analysis and numerical results show that the presented matrix decomposition scheme is the optimal one The presented algoithm has else physicalmeaning and can be busily applied to finite element codes

A Posteriori Error Estimates for Finite Element Methods for Systems of Nonlinear,Dispersive Equations

The present study regards the numerical approximation of solutions of systems of Korteweg-de Vries type,coupled through their nonlinear terms.In our previous work[9],we constructed conservative and dissipative finite element methods for these systems and presented a priori error estimates for the semidiscrete schemes.In this sequel,we present a posteriori error estimates for the semidiscrete and fully discrete approximations introduced in[9].The key tool employed to effect our analysis is the dispersive reconstruction devel-oped by Karakashian and Makridakis[20]for related discontinuous Galerkin methods.We conclude by providing a set of numerical experiments designed to validate the a posteriori theory and explore the effectivity of the resulting error indicators.

Research of the Higher-order Finite Element Arithmetic for Radiation Exchange

With the development of aeronautic and astronautic techniques, radiation becomes much more significant while the structure is exposed to the higher and higher temperature. Most of the current finite element software packages treat it using the net-radiation method or absorbed radiation method based on the assumption of isothermal surface with uniform radiation heat flux, which brings the conflict between the precision and the quantity of grids. Using integral method to compute the variable radiation heat flux in higher-order finite element, the precision can be improved greatly while using the same quantity of grids, because it is more consistent with the distribution of real temperature. In this paper, the integral is only processed on the same integral points as those used for solving the finite element equations, so it may be of high efficiency. In an academic testing model, the result is contrast to which get in ANSYS, proving the high precision of the method. Then an actual sandwich panel used in the thermal protection system is analyzed with the method, and the error is comparatively low to the analytical answer while the computation being of high efficiency.

CHARACTERISTIC GALERKIN METHOD FOR CONVECTION-DIFFUSION EQUATIONS AND IMPLICIT ALGORITHM USING PRECISE INTEGRATION

This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability.

The Numerical Integration of Discrete Functions on a Triangular Element

With the application of Hammer integral formulas of a continuous function on a triangular element, the numerical integral formulas of some discrete functions on the element are derived by means of decomposition and recombination of base functions. Hammer integral formulas are the special examples of those of the paper.

On the Total Dynamic Response of Soil-Structure Interaction System in Time Domain Using Elastodynamic Infinite Elements with Scaling Modified Bessel Shape Functions

This paper is devoted to a new approach—the dynamic response of Soil-Structure System (SSS), the far field of which is discretized by decay or mapped elastodynamic infinite elements, based on scaling modified Bessel shape functions are to be calculated. These elements are appropriate for Soil-Structure Interaction problems, solved in time or frequency domain and can be treated as a new form of the recently proposed elastodynamic infinite elements with united shape functions (EIEUSF) infinite elements. Here the time domain form of the equations of motion is demonstrated and used in the numerical example. In the paper only the formulation of 2D horizontal type infinite elements (HIE) is used, but by similar techniques 2D vertical (VIE) and 2D corner (CIE) infinite elements can also be added. Continuity along the artificial boundary (the line between finite and infinite elements) is discussed as well and the application of the proposed elastodynamical infinite elements in the Finite element method is explained in brief. A numerical example shows the computational efficiency and accuracy of the proposed infinite elements, based on scaling Bessel shape functions.

Finite element analysis of stresses and interface crack in TBC systems

Thermal barrier coatings have been used on high temperature components. Due to high stresses leading to unpredictable failure, a transient thermal-structural finite element solution was employed to analyze the stress distribution and J-integral at the interface between the bond coating and thermally growing oxide(TGO) in the EB-PVD thermal barrier coatings subjected to thermal loadings. The effects of some environmental and material parameters were studied, such as thermal convection coefficient, ceramic elastic modulus and TGO thickness. The results show that the stresses and J-integral values are impacted by these parameters.

A Conservative SAV-RRK Finite Element Method for the Nonlinear Schrodinger Equation

Abstract.In this paper,we propose,analyze and numerically validate a conservative finite element method for the nonlinear Schrodinger equation.A scalar auxiliary variable(SAV)is introduced to reformulate the nonlinear Schrodinger equation into an equivalent system and to transform the energy into a quadratic form.We use the standard continuous finite element method for the spatial discretization,and the relaxation Runge-Kutta method for the time discretization.Both mass and energy conservation laws are shown for the semi-discrete finite element scheme,and also preserved for the full-discrete scheme with suitable relaxation coefficient in the relaxation Runge-Kutta method.Numerical examples are presented to demonstrate the accuracy of the proposed method,and the conservation of mass and energy in long time simulations.

Unconditional Superconvergence Analysis of Energy Conserving Finite Element Methods for the Nonlinear Coupled Klein-Gordon Equations

In this paper,we consider the energy conserving numerical scheme for coupled nonlinear Klein-Gordon equations.We propose energy conserving finite element method and get the unconditional superconvergence resultO(h^(2)+Dt^(2))by using the error splitting technique and postprocessing interpolation.Numerical experiments are carried out to support our theoretical results.

CONVERGENCE OF AN EXPLICIT UPWIND FINITE ELEMENT METHOD TO MULTI-DIMENSIONAL CONSERVATION LAWS

Investigates the convergence of an explicit upwind finite element method, the edge-averaged finite element method to multidimensional scalar conservation laws. Discussion on the finite element scheme and main convergence theorem; Proofs of the technical lemmas.

Analysis of multiple interfacial cracks in three-dimensional bimaterials using hypersingular integro-differential equation method

By using the concept of finite-part integral, a set of hypersingular integro-differential equations for multiple interracial cracks in a three-dimensional infinite bimaterial subjected to arbitrary loads is derived. In the numerical analysis, unknown displacement discontinuities are approximated with the products of the fundamental density functions and power series. The fundamental functions are chosen to express a two-dimensional interface crack rigorously. As illustrative examples, the stress intensity factors for two rectangular interface cracks are calculated for various spacing, crack shape and elastic constants. It is shown that the stress intensity factors decrease with the crack spacing.

A New Iterative Method for Multi-Moving Boundary Problems Based Boundary Integral Method

The present paper deals with very important practical problems of wide range of applications. The main target of the present paper is to track all moving boundaries that appear throughout the whole process when dealing with multi-moving boundary problems continuously with time up to the end of the process with high accuracy and minimum number of iterations. A new numerical iterative scheme based the boundary integral equation method is developed to track the moving boundaries as well as compute all unknowns in the problem. Three practical applications, one for vaporization and two for ablation were solved and their results were compared with finite element, heat balance integral and the source and sink results and a good agreement were obtained.