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.

Highly Accurate Golden Section Search Algorithms and Fictitious Time Integration Method for Solving Nonlinear Eigenvalue Problems

This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.

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.

Application of Choquet Integral-Importance-Performance Analysis and TOPSIS Methods in Approaching the Preference Factors of Calligraphy and Seal Engraving Imagery

Classical Chinese characters,presented through calligraphy,seal engraving,or painting,can exhibit different aesthetics and essences of Chinese characters,making them the most important asset of the Chinese people.Calligraphy and seal engraving,as two closely related systems in traditional Chinese art,have developed through the ages.Due to changes in lifestyle and advancements in modern technology,their original functions of daily writing and verification have gradually diminished.Instead,they have increasingly played a significant role in commercial art.This study utilizes the Evaluation Grid Method(EGM)and the Analytic Hierarchy Process(AHP)to research the key preference factors in the application of calligraphy and seal engraving imagery.Different from the traditional 5-point equal interval semantic questionnaire,this study employs a non-equal interval semantic questionnaire with a golden ratio scale,distinguishing the importance ratio of adjacent semantic meanings and highlighting the weighted emphasis on visual aesthetics.Additionally,the study uses Importance-Performance Analysis(IPA)and Technique for Order Preference by Similarity to Ideal Solution(TOPSIS)to obtain the key preference sequence of calligraphy and seal engraving culture.Plus,the Choquet integral comprehensive evaluation is used as a reference for IPA comparison.It is hoped that this study can provide cultural imagery references and research methods,injecting further creativity into industrial design.

Teaching Method for Comprehensive English Integration of Grammar-translation Method and Communicative Approach

Comprehensive English is a very basic and important course for English majors,according to the features of the text-book A New English Course,teachers should adopt the integration of grammar translation method and communicative approach to improve students' linguistic competence and communicative competence.

Study on Characteristics of 3-D Translating-Pulsating Source Green Function of Deep-Water Havelock Form and Its Fast Integration Method

The singularities, oscillatory performances and the contributing factors to the 3-＇D translating-pulsating source Green function of deep-water Havelock form which consists of a local disturbance part and a far-field wave-like part, are analyzed systematically. Relative numerical integral methods about the two parts are presented in this paper. An improved method based on LOBATTO rule is used to eliminate singularities caused respectively by infinite discontinuity and jump discontinuous node from the local disturbance part function, which makes the improvement of calculation efficiency and accuracy possible. And variable substitution is applied to remove the singularity existing at the end of the integral interval of the far-field wave-like part function. Two auxiliary techniques such as valid interval calculation and local refinement of integral steps technique in narrow zones near false singularities are applied so as to avoid unnecessary integration of invalid interval and improve integral accordance. Numerical test results have proved the efficiency and accuracy in these integral methods that thus can be applied to calculate hydrodynamic performance of floating structures moving in waves.

PARAMETRIC EQUATIONS OF NONHOLONOMIC NONCONSERVATIVE SYSTEMS IN THE EVENT SPACE AND THE METHOD OF THEIR INTEGRATION

In this paper,the parametric equations with multipliers of nonholonomic nonconservative sys- tems in the event space are established,their properties are studied,and their explicit formulation is obtained. And then the field method for integrating these equations is given.Finally,an example illustrating the appli- cation of the integration method is given.

Study of probability integration method parameter inversion by the genetic algorithm

In order to obtain accurate probability integration method(PIM) parameters for surface movement of multi-panel mining, a genetic algorithm(GA) was used to optimize the parameters. As the measured surface movement is affected by more than one mining panel, traditional PIM parameter inversion model is difficult to ensure the reliability of the results due to the complexity of rock movement. With crossover,mutation and selection operators, GA can perform a global optimization search and has high computation efficiency. Compared with the pattern search algorithm, the fitness function can avoid falling into local minima traps. GA reduces the risk of local minima traps which improves the accuracy and reliability with the mutation mechanism. Application at Xuehu colliery shows that GA can be used to inverse the PIM parameters for multi-panel surface movement observation, and reliable results can be obtained. The research provides a new way for back-analysis of PIM parameters for mining subsidence under complex conditions.

Model for the Whole Roller Leveling Process of Plates with Random Curvature Distribution Based on the Curvature Integration Method

A model based on the curvature integration method has been applied in an online plate leveling system. However, there are some shortcomings in the current leveling models. On the one hand, the models cannot deal with the leveling process of plates with a random curvature distribution. On the other hand, the current models are suitable only for stable leveling processes and ignore the biting in and tailing out stages. This study presents a new plate-leveling model based on the curvature integration method, which can describe the leveling process of plates with random curvature distribution. Further, the model is solved in two cases in order to take the biting in and tailing out stages into consideration. The proposed model is evaluated by comparing with a plate leveling experiment. Finally, the leveling process of a plate with a wave bent is studied using the proposed model. It is found that the contact angles vary greatly during the biting in and tailing out stages. However, they are relatively steady during the 5 roller leveling stage. In addition, the contact angle of roller No. 2 is the smallest, which is close to 0. Roller leveling can effectively eliminate bending in the plate, but there are regions in the head and tail of the plate, where roller leveling is not effective. The non-leveling region length is about 2 times that of the roller space. This study proposes a quasi-static plate-leveling model, which makes it possible to analyze the dynamic straightening process using a curvature integration method. It also makes it possible to analyze the straightening process of a plate with random curvature distribution.

Modified precise time step integration method of structural dynamic analysis

The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method （e.g., an algorithm with an arbitrary order of accuracy） is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.

Poisson theory and integration method of Birkhoffian systems in the event space

This paper focuses on studying the Poisson theory and the integration method of a Birkhoffian system in the event space. The Birkhoff＇s equations in the event space are given. The Poisson theory of the Birkhoffian system in the event space is established. The definition of the Jacobi last multiplier of the system is given, and the relation between the Jacobi last multiplier and the first integrals of the system is discussed. The researches show that for a Birkhoffian system in the event space, whose configuration is determined by （2n ＋ 1） Birkhoff＇s variables, the solution of the system can be found by the Jacobi last multiplier if 2n first integrals are known. An example is given to illustrate the application of the results.

Fast precise integration method for hyperbolic heat conduction problems

A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to the fact that the matrix exponential is sparse. The presented method employs the sparsity of the matrix exponential to improve the original precise integration method. The merits are that the proposed method is suitable for large hyperbolic heat equations and inherits the accuracy of the original version and the good computational efficiency, which are verified by two numerical examples.

Response Sensitivity Analysis of the Dynamic Milling Process Based on the Numerical Integration Method

As one of the bases of gradient-based optimization algorithms, sensitivity analysis is usually required to calculate the derivatives of the system response with respect to the machining parameters. The most widely used approaches for sensitivity analysis are based on time-consuming numerical methods, such as finite difference methods. This paper presents a semi-analytical method for calculation of the sensitivity of the stability boundary in milling. After transforming the delay-differential equation with time-periodic coefficients governing the dynamic milling process into the integral form, the Floquet transition matrix is constructed by using the numerical integration method. Then, the analytical expressions of derivatives of the Floquet transition matrix with respect to the machining parameters are obtained. Thereafter, the classical analytical expression of the sensitivity of matrix eigenvalues is employed to calculate the sensitivity of the stability lobe diagram. The two-degree-of-freedom milling example illustrates the accuracy and efficiency of the proposed method. Compared with the existing methods, the unique merit of the proposed method is that it can be used for analytically computing the sensitivity of the stability boundary in milling, without employing any finite difference methods. Therefore, the high accuracy and high efficiency are both achieved. The proposed method can serve as an effective tool for machining parameter optimization and uncertainty analysis in high-speed milling.

Evaluation of integration methods for hybrid simulation of complex structural systems through collapse

This study examines the performance of integration methods for hybrid simulation of large and complex structural systems in the context of structural collapse due to seismic excitations. The target application is not necessarily for real-time testing, but rather for models that involve large-scale physical sub-structures and highly nonlinear numerical models. Four case studies are presented and discussed. In the first case study, the accuracy of integration schemes including two widely used methods, namely, modified version of the implicit Newmark with fixed-number of iteration （iterative） and the operator-splitting （non-iterative） is examined through pure numerical simulations. The second case study presents the results of 10 hybrid simulations repeated with the two aforementioned integration methods considering various time steps and fixed-number of iterations for the iterative integration method. The physical sub-structure in these tests consists of a single-degree-of-freedom （SDOF） cantilever column with replaceable steel coupons that provides repeatable highly- nonlinear behavior including fracture-type strength and stiffness degradations. In case study three, the implicit Newmark with fixed-number of iterations is applied for hybrid simulations of a 1：2 scale steel moment frame that includes a relatively complex nonlinear numerical substructure. Lastly, a more complex numerical substructure is considered by constructing a nonlinear computational model of a moment frame coupled to a hybrid model ofa 1：2 scale steel gravity frame. The last two case studies are conducted on the same porotype structure and the selection of time steps and fixed number of iterations are closely examined in pre-test simulations. The generated unbalance forces is used as an index to track the equilibrium error and predict the accuracy and stability of the simulations.

A COMBINED PARAMETRIC QUADRATIC PROGRAMMING AND PRECISE INTEGRATION METHOD BASED DYNAMIC ANALYSIS OF ELASTIC-PLASTIC HARDENING/SOFTENING PROBLEMS

The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcome the result mesh-sensitivity problem in the dynamic strain softening or strain localization analysis. The equations for the dynamic elastic-plastic problems are derived in terms of the parametric variational principle, which is valid for associated, non-associated and strain softening plastic constitutive models in the finite element analysis. The precise integration method, which has been widely used for discretization in time domain of the linear problems, is introduced for the solution of dynamic nonlinear equations. The new algorithm proposed is based on the combination of the parametric quadratic programming method and the precise integration method and has all the advantages in both of the algorithms. Results of numerical examples demonstrate not only the validity, but also the advantages of the algorithm proposed for the numerical solution of nonlinear dynamic problems.

A potential integration method for Birkhoffian system

This paper is intended to apply the potential integration method to the differential equations of the Birkhoffian system. The method is that, for a given Birkhoffian system, its differential equations are first rewritten as 2n first-order differential equations. Secondly, the corresponding partial differential equations are obtained by potential integration method and the solution is expressed as a complete integral. Finally, the integral of the system is obtained.

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.

ADAPTIVE INTERVAL WAVELET PRECISE INTEGRATION METHOD FOR PARTIAL DIFFERENTIAL EQUATIONS

The quasi-Shannon interval wavelet is constructed based on the interpolation wavelet theory, and an adaptive precise integration method, which is based on extrapolation method is presented for nonlinear ordinary differential equations ( ODEs). And then, an adaptive interval wavelet precise integration method (AIWPIM) for nonlinear partial differential equations(PDEs) is proposed. The numerical results show that the computational precision of AIWPIM is higher than that of the method constructed by combining the wavelet and the 4th Runge-Kutta method, and the computational amounts of these two methods are almost equal. For convenience, the Burgers equation is taken as an example in introducing this method, which is also valid for more general cases.

Selection of the best initial orbital elements of satellite based on fuzzy integration evaluation method

The fuzzy integration evaluation method （FIEM） is studied in order to select the best orbital elements from the multi-group initial orbits determined by a satellite TT＆C （Tracking, Telemetry and Control） center with all kinds of data sources. By employing FIEM together with the experience of TT＆C experts, the index system to evaluate the selection of the best initial orbits is established after the data sources and orbit determination theories are studied. Besides, the concrete steps in employing the method are presented. Moreover, by taking the objects to be evaluated as evaluation experts, the problem of how to generate evaluation matrices is solved. Through practical application, the method to select the best initial orbital elements has been proved to be flexible and effective The originality of the method is to find a new evaluation criterion （comparing the actually tracked orbits） replacing the traditional one （comparing the nominal orbits） for selecting the best orbital elements.

Comparison of capability of time integration methods in capturing dynamic loading

Numerical properties of the time integration method proposed by the first author of this paper in 2007 are the same as those of the constant average acceleration method （AAM） for linear elastic systems, except that the capability to capture dynamic loading was not explored. It was found that there were different quadrature equations to predict the next step displacement increment. A modified quadrature equation of this method was derived so that the equation to determine the next step displacement was numerically equivalent to the equation used in the constant AAM. It was verified that the original form of this method, in general, had a better capability to capture dynamic loadings than the constant AAM. This excellent property, in addition to computational efficiency, will help to make this method competitive with general secondorder accurate integration methods.