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 coupled Legendre-Laguerre polynomial method with analytical integration for the Rayleigh waves in a quasicrystal layered half-space with an imperfect interface

The Laguerre polynomial method has been successfully used to investigate the dynamic responses of a half-space.However,it fails to obtain the correct stress at the interfaces in a layered half-space,especially when there are significant differences in material properties.Therefore,a coupled Legendre-Laguerre polynomial method with analytical integration is proposed.The Rayleigh waves in a one-dimensional(1D)hexagonal quasicrystal(QC)layered half-space with an imperfect interface are investigated.The correctness is validated by comparison with available results.Its computation efficiency is analyzed.The dispersion curves of the phase velocity,displacement distributions,and stress distributions are illustrated.The effects of the phonon-phason coupling and imperfect interface coefficients on the wave characteristics are investigated.Some novel findings reveal that the proposed method is highly efficient for addressing the Rayleigh waves in a QC layered half-space.It can save over 99%of the computation time.This method can be expanded to investigate waves in various layered half-spaces,including earth-layered media and surface acoustic wave(SAW)devices.

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.

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.

Least square method based on Haar wavelet to solve multi-dimensional stochastic Ito-Volterra integral equations

This paper proposes a method combining blue the Haar wavelet and the least square to solve the multi-dimensional stochastic Ito-Volterra integral equation.This approach is to transform stochastic integral equations into a system of algebraic equations.Meanwhile,the error analysis is proven.Finally,the effectiveness of the approach is verified by two numerical examples.

Research and application of sidewall stability predic-tion method based on analytic hierarchy process and fuzzy integrative evaluation method

As a difficult problem, sidewall instability has been beset drilling workers all the time. Not only does it cause huge economic losses, but also it determines the success or failure of drilling engineering. Due to complex relationship between various factors which influence sidewall stability, it hasn’t been found a widely applied method to predicate sidewall stability so far. Therefore, in order to formulate corresponding measures to ensure successful drilling, searching for a kind of better method to forecast sidewall stability before drilling becomes an imperative and significant topic for drilling engineering. On the basis of traditional sidewall stability analytical method, we have put forward the Fuzzy Comprehensive Evaluation Method to forecast sidewall stability regulation using physico-chemical performance parameters of the clay mineral. This method has been improved by introducing the Analytic Hierarchy Process (AHP) and the Maximum Subjection Principle in the application process. After introducing Analytic Hierarchy Process to identify weight, and Maximum Subjection Principle to obtain evaluation results, it has reduced the influence of human factors and enhanced the accuracy of the fuzzy evaluation results. The application in Hailaer Area indicates that this method can predict sidewall stability of gas-oil well with high credibility and strong practicability.

Exponential Time Differencing Method for a Reaction-Diffusion System with Free Boundary

For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geometries may lead to difficulties in the accuracy when discretizing the high-order derivatives on grid points near the boundary.It is very challenging to design numerical methods that can efficiently and accurately handle both difficulties.Applying an implicit scheme may be able to remove the stability constraints on the time step,however,it usually requires solving a large global system of nonlinear equations for each time step,and the computational cost could be significant.Integration factor(IF)or exponential time differencing(ETD)methods are one of the popular methods for temporal partial differential equations(PDEs)among many other methods.In our paper,we couple ETD methods with an embedded boundary method to solve a system of reaction-diffusion equations with complex geometries.In particular,we rewrite all ETD schemes into a linear combination of specificФ-functions and apply one state-of-the-art algorithm to compute the matrix-vector multiplications,which offers significant computational advantages with adaptive Krylov subspaces.In addition,we extend this method by incorporating the level set method to solve the free boundary problem.The accuracy,stability,and efficiency of the developed method are demonstrated by numerical examples.

On power series statistical convergence and new uniform integrability of double sequences

In the present paper,we mostly focus on P_(p)^(2)-statistical convergence.We will look into the uniform integrability via the power series method and its characterizations for double sequences.Also,the notions of P_(p)^(2)-statistically Cauchy sequence,P_(p)^(2)-statistical boundedness and core for double sequences will be described in addition to these findings.

Development Level and Dynamic Analysis of the Integration of Three Rural Industries in the Context of Rural Revitalization

As a key measure to comprehensively promote the strategy of rural revitalization,especially in the field of industrial prosperity,the integration of rural three industries is of great strategic significance.In this study,entropy method and TOPSIS method were employed to calculate the comprehensive evaluation index of integration of three industries and the ideal development value respectively.The development status of integration of three rural industries was systematically evaluated,and the development trend of different regions was compared and analyzed.The results indicate that the development index of integration of three industries showed a steady growth trend,among which the development level of Jiangsu Province ranked first,followed by Hubei Province,while the development level of Shaanxi Province was relatively low.When analyzing the Level I indicators of integration of three rural industries,the contribution of industrial integration behavior exceeded the performance of industrial integration.Among the behavioral indicators of industrial integration,the weight of agricultural multi-functionality and service integration was relatively large,which plays a significant role in promoting the development of integration of three rural industries.The scores and growth rates of Jiangsu Province in increasing farmers income,increasing agricultural production and rural economic development were higher than those of other regions,while Shaanxi Province still had a certain gap in rural industrial integration indicators compared with Jiangsu Province and Hubei Province.In view of this,we came up with some strategic recommendations for further promoting the development of integration of three rural industries.

Application of integrated geophysical techniques in geothermal exploration in Binhai County, Jiangsu Province

Integrated geophysical technology is a necessary and effective means for geothermal exploration.However,integration of geophysical technology for large‐scale surveys with those for geothermal reservoir localization is still in development.This study used the controlled source audio‐frequency magnetotelluric method technology for large‐scale exploration to obtain underground electrical structure information and micromotion detection technology to obtain underground wave velocity structure information.The combination of two detection technologies was used for local identification of geothermal reservoirs.Further,auxiliary correction and inversion constraint were implemented through the audio magnetotelluric sounding technology for maximum authenticity restoration of the near‐and transition‐field data.Through these technology improvements,a geothermal geological model was established for the Binhai County of Jiangsu Province in China and potential geothermal well locations were identified.On this basis,a geothermal well was drilled nearly 3000m deep,with a daily water volume of over 2000m3/day and a geothermal water temperature of 51°C at the well head.It is found that predictions using the above integrated geophysical exploration technology are in good agreement with the well geological formation data.This integrated geophysical technology can be effectively applied for geothermal exploration with high precision and reliability.

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.

A model for extracting large deformation mining subsidence using D-InSAR technique and probability integral method

Due to the difficulties in obtaining large deformation mining subsidence using differential Interferometric Synthetic Aperture Radar （D-InSAR） alone, a new algorithm was proposed to extract large deformation mining subsidence using D-InSAR technique and probability integral method. The details of the algorithm are as follows：the control points set, containing correct phase unwrapping points on the subsidence basin edge generated by D-InSAR and several observation points （near the maximum subsidence and inflection points）, was established at first; genetic algorithm （GA） was then used to optimize the parameters of probability integral method; at last, the surface subsidence was deduced according to the optimum parameters. The results of the experiment in Huaibei mining area, China, show that the presented method can generate the correct mining subsidence basin with a few surface observations, and the relative error of maximum subsidence point is about 8.3%, which is much better than that of conventional D-InSAR （relative error is 68.0%）.

A new method for integration of a Birkhoffian system

The idea of the gradient method for integrating the dynamical equations of a nonconservative system presented by Vujanovic is transplanted to a Birkhoffian system, which is a new method for the integration of Birkhoff＇s equations. First, the differential equations of motion of the Birkhoffian system are written out. Secondly, 2n Birkhoff＇s variables are divided into two parts, and assume that a part of the variables is the functions of the remaining part of the variables and time. Thereby, the basic quasi-linear partial differential equations are established. If a complete solution of the basic partial differential equations is sought out, the solution of the problem is given by a set of algebraic equations. Since one can choose n arbitrary Birkhoff＇s variables as the functions of n remains of variables and time in a specific problem, the method has flexibility. The major difficulty of this method lies in finding a complete solution of the basic partial differential equation. Once one finds the complete solution, the motion of the systems can be obtained without doing further integration. Finally, two examples are given to illustrate the application of the results.

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.

3D elastic wave equation forward modeling based on the precise integration method

The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data. However, because of its discretization, the FD method is only stable under certain conditions. The Arbitrary Difference Precise Integration （ADPI） method is based on the FD method and adopts an integration scheme in the time domain and an arbitrary difference scheme in the space domain. Therefore, the ADPI method is a semi-analytical method. In this paper, we deduce the formula for the ADPI method based on the 3D elastic equation and improve its stability. In forward modeling cases, the ADPI method was implemented in 2D and 3D elastic wave equation forward modeling. Results show that the travel time of the reflected seismic wave is accurate. Compared with the acoustic wave field, the elastic wave field contains more wave types, including PS- and PP- reflected waves, transmitted waves, and diffracted waves, which is important to interpretation of seismic data. The method can be easily applied to elastic wave equation numerical simulations for eoloical models.

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.

An integrated method of selecting environmental covariates for predictive soil depth mapping

Environmental covariates are the basis of predictive soil mapping.Their selection determines the performance of soil mapping to a great extent,especially in cases where the number of soil samples is limited but soil spatial heterogeneity is high.In this study,we proposed an integrated method to select environmental covariates for predictive soil depth mapping.First,candidate variables that may influence the development of soil depth were selected based on pedogenetic knowledge.Second,three conventional methods(Pearson correlation analysis(PsCA),generalized additive models(GAMs),and Random Forest(RF))were used to generate optimal combinations of environmental covariates.Finally,three optimal combinations were integrated to produce a final combination based on the importance and occurrence frequency of each environmental covariate.We tested this method for soil depth mapping in the upper reaches of the Heihe River Basin in Northwest China.A total of 129 soil sampling sites were collected using a representative sampling strategy,and RF and support vector machine(SVM)models were used to map soil depth.The results showed that compared to the set of environmental covariates selected by the three conventional selection methods,the set of environmental covariates selected by the proposed method achieved higher mapping accuracy.The combination from the proposed method obtained a root mean square error(RMSE)of 11.88 cm,which was 2.25–7.64 cm lower than the other methods,and an R^2 value of 0.76,which was 0.08–0.26 higher than the other methods.The results suggest that our method can be used as an alternative to the conventional methods for soil depth mapping and may also be effective for mapping other soil properties.

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.