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.

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.

Obstacles to the Integrated Care of Children Aged 6 to 59 Months Suffering from Acute and Severe Malnutrition at the Chu-Me and HATC of N’Djamena/Chad

Introduction: the problem of Severe Acute Malnutrition (SAM) and its consequences concern and challenge all social actors worldwide. The objective of the study was to identify the real obstacles that actors face in the integrated management of children aged 6 - 59 months suffering from SAM at the University Hospital of Mother and Child (CHUME) and the Chad-China Friendship Hospital (HATC) of N’Djamena. Methodology: this is a cross-sectional descriptive and interpretative study conducted from January to October 2023 at the CHUME and HATC of N’Djamena. It is based on questionnaires and direct observation of two hundred and thirty-six (236) parents or guardians of malnourished children admitted to these health facilities. All 236 agreed to participate freely in this study. Results: it appears from this study that the obstacles to the integrated management of malnourished children were: poor care-care relationship (26.7%) permanent product breaks (22.9%), the direct costs of services ranged from 30,000 to 60,000 FCFA (41.53%), drugs used for the treatment of patients came from the street (66.10%), malaria and anemia occupy an important place (27.5%), 57.2% respondents had monthly income below 60,000 CFAF. Conclusion: the present study carried out has identified the real bottlenecks to the integrated management of children suffering from SAM in the therapeutic nutritional units of the city of N’Djamena. .

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.

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.

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.

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.

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.

Precise integration method for a class of singular two-point boundary value problems

In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the method of variable coefficient dimensional expanding,the non-homogeneous ordinary differential equations（ODEs） are transformed into homogeneous ODEs.Then the interval is divided evenly,and the transfer matrix in every subinterval is worked out using the high order multiple perturbation method,and a set of algebraic equations is given in the form of matrix by the precise integration relation for each segment,which is worked out by the reduction method.Finally numerical examples are elaboratedd to validate the present method.

Self-Adaptive Filon’s Integration Method and Its Application to Computing Synthetic Seismograms

Based on the principle of self-adaptive Simpson integration method, and by incorporating the‘fifth-order’Filon’s integration algorithm [Bull. Seism. Soc. Am. 73(1983) 913], we have proposed a simple and efficient numerical integration method, i.e., the self-adaptive Filon’s integration method (SAFIM), for computing the synthetic seismograms at large epicentral distances. With numerical examples, we have demonstrated that the SAFIM is not only accurate but also very efficient. This new integration method is expected to be very useful in seismology, as well as in computing similar oscillatory integrals in other branches of physics.

Compact implicit integration factor methods for some complex-valued nonlinear equations

The compact implicit integration factor （cIIF） method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF method to some complex-valued nonlinear evolutionary equations such as the nonlinear SchrSdinger （NLS） equation and the complex Ginzburg-Landau （GL） equation. Detailed algorithm formulation and practical implementation of cIIF method are performed. The numerical results indicate that this method is very accurate and efficient.

The efficiency of direct integration methods in elastic contact-impact problems

A comparison of direct integration methods is madeand their efficiency is investigated for impact problems.New-mark,Wilson-θ,Central Difference and Houbolt Methodsare used as direct integration methods.Impact analysisincludes that of elastic and large deformation based uponupdated Lagrangian including buckling check.The resultsshow that the direct integration methods give differentresults in different contact-impact cases.

Short-term solar flare prediction using multi-model integration method

A multi-model integration method is proposed to develop a multi-source and heterogeneous model for short-term solar flare prediction. Different prediction models are constructed on the basis of extracted predictors from a pool of observation databases. The outputs of the base models are normal- ized first because these established models extract predictors from many data resources using different prediction methods. Then weighted integration of the base models is used to develop a multi-model integrated model （MIM）. The weight set that single models assign is optimized by a genetic algorithm. Seven base models and data from Solar and Heliospheric Observatory/Michelson Doppler Imager lon- gitudinal magnetograms are used to construct the MIM, and then its performance is evaluated by cross validation. Experimental results showed that the MIM outperforms any individual model in nearly every data group, and the richer the diversity of the base models, the better the performance of the MIM. Thus, integrating more diversified models, such as an expert system, a statistical model and a physical model, will greatly improve the performance of the MIM.

基于TFN-AHP熵权的铁钻工工业设计方案评价

针对铁钻工工业设计方案评价综合性强、存在一定主观性、评价影响因素多的问题,提出一种基于三角模糊数层次分析法和熵权法主、客观方法相结合的铁钻工方案评价模型。该模型在层次分析法框架下,根据专家意见构建铁钻工工业设计方案评价指标体系,划分4项1级指标和12项2级指标,以模糊决策矩阵为评价信息载体,引入三角模糊数优化传统层次分析法不能准确表征带有模糊不确定性特征的评价信息的问题,并使用客观方法熵权法组合赋权确定综合权重。将此评价模型运用于实际铁钻工设计项目方案评价,将评价结果与相关传统评价法进行对比,并将筛选出来的方案进行深入优化设计,证明了评价模型的有效性。研究结果能够有效辅助企业进行铁钻工工业设计方案决策。

Dynamic stability analysis of porous functionally graded beams under hygro-thermal loading using nonlocal strain gradient integral model

We present a study on the dynamic stability of porous functionally graded(PFG)beams under hygro-thermal loading.The variations of the properties of the beams across the beam thicknesses are described by the power-law model.Unlike most studies on this topic,we consider both the bending deformation of the beams and the hygro-thermal load as size-dependent,simultaneously,by adopting the equivalent differential forms of the well-posed nonlocal strain gradient integral theory(NSGIT)which are strictly equipped with a set of constitutive boundary conditions(CBCs),and through which both the stiffness-hardening and stiffness-softening effects of the structures can be observed with the length-scale parameters changed.All the variables presented in the differential problem formulation are discretized.The numerical solution of the dynamic instability region(DIR)of various bounded beams is then developed via the generalized differential quadrature method(GDQM).After verifying the present formulation and results,we examine the effects of different parameters such as the nonlocal/gradient length-scale parameters,the static force factor,the functionally graded(FG)parameter,and the porosity parameter on the DIR.Furthermore,the influence of considering the size-dependent hygro-thermal load is also presented.

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.