The use of the node-based smoothed finite element method to estimate static and seismic bearing capacities of shallow strip footings

The node-based smoothed finite element method(NS-FEM)is shortly presented for calculations of the static and seismic bearing capacities of shallow strip footings.A series of computations has been performed to assess variations in seismic bearing capacity factors with both horizontal and vertical seismic accelerations.Numerical results obtained agree very well with those using the slip-line method,revealing that the magnitude of the seismic bearing capacity is highly dependent upon the combinations of various directions of both components of the seismic acceleration.An upward vertical seismic acceleration reduces the seismic bearing capacity compared to the downward vertical seismic acceleration in calculations.In addition,particular emphasis is placed on a separate estimation of the effects of soil and superstructure inertia on each seismic bearing capacity component.While the effect of inertia forces arising in the soil on the seismic bearing capacity is non-trivial,and the superstructure inertia is the major contributor to reductions in the seismic bearing capacity.Both tables and charts are given for practical application to the seismic design of the foundations.

Implementing the Node Based Smoothed Finite Element Method as User Element in Abaqus for Linear and Nonlinear Elasticity

In this paper,the node based smoothed-strain Abaqus user element(UEL)in the framework of finite element method is introduced.The basic idea behind of the node based smoothed finite element(NSFEM)is that finite element cells are divided into subcells and subcells construct the smoothing domain associated with each node of a finite element cell[Liu,Dai and Nguyen-Thoi(2007)].Therefore,the numerical integration is globally performed over smoothing domains.It is demonstrated that the proposed UEL retains all the advantages of the NSFEM,i.e.,upper bound solution,overly soft stiffness and free from locking in compressible and nearly-incompressible media.In this work,the constant strain triangular(CST)elements are used to construct node based smoothing domains,since any complex two dimensional domains can be discretized using CST elements.This additional challenge is successfully addressed in this paper.The efficacy and robustness of the proposed work is obtained by several benchmark problems in both linear and nonlinear elasticity.The developed UEL and the associated files can be downloaded from https://github.com/nsundar/NSFEM.

Bubble-Enriched Smoothed Finite Element Methods for Nearly-Incompressible Solids

This work presents a locking-free smoothed finite element method(S-FEM)for the simulation of soft matter modelled by the equations of quasi-incompressible hyperelasticity.The proposed method overcomes well-known issues of standard finite element methods(FEM)in the incompressible limit:the over-estimation of stiffness and sensitivity to severely distorted meshes.The concepts of cell-based,edge-based and node-based S-FEMs are extended in this paper to three-dimensions.Additionally,a cubic bubble function is utilized to improve accuracy and stability.For the bubble function,an additional displacement degree of freedom is added at the centroid of the element.Several numerical studies are performed demonstrating the stability and validity of the proposed approach.The obtained results are compared with standard FEM and with analytical solutions to show the effectiveness of the method.

A Cell-Based Linear Smoothed Finite Element Method for Polygonal Topology Optimization

The aim of this work is to employ a modified cell-based smoothed finite element method(S-FEM)for topology optimization with the domain discretized with arbitrary polygons.In the present work,the linear polynomial basis function is used as the weight function instead of the constant weight function used in the standard S-FEM.This improves the accuracy and yields an optimal convergence rate.The gradients are smoothed over each smoothing domain,then used to compute the stiffness matrix.Within the proposed scheme,an optimum topology procedure is conducted over the smoothing domains.Structural materials are distributed over each smoothing domain and the filtering scheme relies on the smoothing domain.Numerical tests are carried out to pursue the performance of the proposed optimization by comparing convergence,efficiency and accuracy.

A novel twice-interpolation finite element method for solid mechanics problems

Formulation and numerical evaluation of a novel twice-interpolation finite element method （TFEM） is presented for solid mechanics problems. In this method, the trial function for Galerkin weak form is constructed through two stages of consecutive interpolation. The primary interpolation follows exactly the same procedure of standard FEM and is further reproduced according to both nodal values and averaged nodal gradients obtained from primary interpolation. The trial functions thus constructed have continuous nodal gradients and contain higher order polynomial without increasing total freedoms. Several benchmark examples and a real dam problem are used to examine the TFEM in terms of accuracy and convergence. Compared with standard FEM, TFEM can achieve significantly better accuracy and higher convergence rate, and the continuous nodal stress can be obtained without any smoothing operation. It is also found that TFEM is insensitive to the quality of the elemental mesh. In addition, the present TFEM can treat the incompressible material without any modification.

Optimization of buckling load for laminated composite plates using adaptive Kriging-improved PSO:A novel hybrid intelligent method

An effective hybrid optimization method is proposed by integrating an adaptive Kriging(A-Kriging)into an improved partial swarm optimization algorithm(IPSO)to give a so-called A-Kriging-IPSO for maximizing the buckling load of laminated composite plates(LCPs)under uniaxial and biaxial compressions.In this method,a novel iterative adaptive Kriging model,which is structured using two training sample sets as active and adaptive points,is utilized to directly predict the buckling load of the LCPs and to improve the efficiency of the optimization process.The active points are selected from the initial data set while the adaptive points are generated using the radial random-based convex samples.The cell-based smoothed discrete shear gap method(CS-DSG3)is employed to analyze the buckling behavior of the LCPs to provide the response of adaptive and input data sets.The buckling load of the LCPs is maximized by utilizing the IPSO algorithm.To demonstrate the efficiency and accuracy of the proposed methodology,the LCPs with different layers(2,3,4,and 10 layers),boundary conditions,aspect ratios and load patterns(biaxial and uniaxial loads)are investigated.The results obtained by proposed method are in good agreement with the literature results,but with less computational burden.By applying adaptive radial Kriging model,the accurate optimal resultsebased predictions of the buckling load are obtained for the studied LCPs.

A Smoothing Newton Method for the Box Constrained Variational Inequality Problems

The box constrained variational inequality problem can be reformulated as a nonsmooth equation by using median operator.In this paper,we present a smoothing Newton method for solving the box constrained variational inequality problem based on a new smoothing approximation function.The proposed algorithm is proved to be well defined and convergent globally under weaker conditions.

A Smoothing SAA Method for a Stochastic Linear Complementarity Problem

Utilizing the well-known aggregation technique, we propose a smoothing sample average approximation （SAA） method for a stochastic linear complementarity problem, where the underlying functions are represented by expectations of stochastic functions. The method is proved to be convergent and the preliminary numerical results are reported.

Numerical analysis of submarine landslides using a smoothed particle hydrodynamics depth integral model

Submarine landslides can cause severe damage to marine engineering structures. Their sliding velocity and runout distance are two major parameters for quantifying and analyzing the risk of submarine landslides.Currently, commercial calculation programs such as BING have limitations in simulating underwater soil movements. All of these processes can be consistently simulated through a smoothed particle hydrodynamics（SPH） depth integrated model. The basis of the model is a control equation that was developed to take into account the effects of soil consolidation and erosion. In this work, the frictional rheological mode has been used to perform a simulation study of submarine landslides. Time-history curves of the sliding body＇s velocity, height,and length under various conditions of water depth, slope gradient, contact friction coefficient, and erosion rate are compared; the maximum sliding distance and velocity are calculated; and patterns of variation are discussed.The findings of this study can provide a reference for disaster warnings and pipeline route selection.

GLOBAL LINEAR AND QUADRATIC ONE-STEP SMOOTHING NEWTON METHOD FOR VERTICAL LINEAR COMPLEMENTARITY PROBLEMS

A one_step smoothing Newton method is proposed for solving the vertical linear complementarity problem based on the so_called aggregation function. The proposed algorithm has the following good features: (ⅰ) It solves only one linear system of equations and does only one line search at each iteration; (ⅱ) It is well_defined for the vertical linear complementarity problem with vertical block P 0 matrix and any accumulation point of iteration sequence is its solution.Moreover, the iteration sequence is bounded for the vertical linear complementarity problem with vertical block P 0+R 0 matrix; (ⅲ) It has both global linear and local quadratic convergence without strict complementarity. Many existing smoothing Newton methods do not have the property (ⅲ).

The Smoothing Newton Method for Solving the Extended Linear Complementarity Problem

The extended linear complementarity problem(denoted by ELCP) can be reformulated as the solution of a nonsmooth system of equations. By the symmetrically perturbed CHKS smoothing function, the ELCP is approximated by a family of parameterized smooth equations. A one-step smoothing Newton method is designed for solving the ELCP. The proposed algorithm is proved to be globally convergent under suitable assumptions.

A Class of Smoothing Modulus-Based Iterative Method for Solving Implicit Complementarity Problems

In this paper, a class of smoothing modulus-based iterative method was presented for solving implicit complementarity problems. The main idea was to transform the implicit complementarity problem into an equivalent implicit fixed-point equation, then introduces a smoothing function to obtain its approximation solutions. The convergence analysis of the algorithm was given, and the efficiency of the algorithms was verified by numerical experiments.

2020—2025年广东省医疗机构床位需求预测

目的预测2020—2025年广东省医疗机构的床位需求总量。方法基于卫生服务需求法与Holt双参数指数平滑模型,结合年龄别人口数据预测床位需求。结果2025年,广东省住院人数为2425.11万人,床位需求数为70.04万张,每千常住人口床位需求数为5.63张、每千常住人口拥有床位数5.55张,供需比例为98.58%。预测模型的平均百分误差为1.63%(标准差=±1.90%,均方根误差=20.63)。结论结合人口的年龄结构进行预测结果更稳定、误差更小。2020年,广东省的床位配置量基本能满足床位需求,供需较为平衡。但2024年床位需求将超过床位配置总量。未来,广东省应加大床位资源的投入力度,提高基层卫生机构的床位利用率,全面落实分级诊疗制度。

SMOOTHING NEWTON ALGORITHM FOR THE CIRCULAR CONE PROGRAMMING WITH A NONMONOTONE LINE SEARCH

In this paper, we present a nonmonotone smoothing Newton algorithm for solving the circular cone programming(CCP) problem in which a linear function is minimized or maximized over the intersection of an affine space with the circular cone. Based on the relationship between the circular cone and the second-order cone(SOC), we reformulate the CCP problem as the second-order cone problem(SOCP). By extending the nonmonotone line search for unconstrained optimization to the CCP, a nonmonotone smoothing Newton method is proposed for solving the CCP. Under suitable assumptions, the proposed algorithm is shown to be globally and locally quadratically convergent. Some preliminary numerical results indicate the effectiveness of the proposed algorithm for solving the CCP.

AN IMPROVED MODAL ANALYSIS FOR THREE-DIMENSIONAL PROBLEMS USING FACE-BASED SMOOTHED FINITE ELEMENT METHOD

In this work, we further extended the face-based smoothed finite element method （FS-FEM） for modal analysis of three-dimensional solids using four-node tetrahedron elements. The FS-FEM is formulated based on the smoothed Calerkin weak form which employs smoothed strains obtained using the gradient smoothing operation on face-based smoothing domains. This strain smoothing operation can provide softening effect to the system stiffness and make the FSFEM provide more accurate eigenfrequency prediction than the FEM does. Numerical studies have verified this attractive property of FS-FEM as well as its ability and effectiveness on providing reliable eigenfrequency and eigenmode prediction in practical engineering application.

Multi-resolution technique integrated with smoothed particle element method (SPEM) for modeling fluid-structure interaction problems with free surfaces

Free-surface flows, especially those associated with fluid-structure interactions(FSIs), pose challenging problems in numerical simulations. The authors of this work recently developed a smoothed particle element method(SPEM) to simulate FSIs. In this method, both the fluid and solid regions are initially modeled using a smoothed finite element method(S-FEM) in a Lagrangian frame, whereas the fluid regions undergoing large deformations are adaptively converted into particles and modeled with an improved smoothed particle hydrodynamics(SPH) method. This approach greatly improves computational accuracy and efficiency because of the advantages of the S-FEM in efficiently treating solid/fluid regions showing small deformations and the SPH method in effectively modeling moving interfaces. In this work, we further enhance the efficiency of the SPEM while effectively capturing local fluid information by introducing a multi-resolution technique to the SPEM and developing an effective approach to treat multi-resolution element-particle interfaces. Various numerical examples demonstrate that the multiresolution SPEM can significantly reduce the computational cost relative to the original version with a constant resolution.Moreover, the novel approach is effective in modeling various incompressible flow problems involving FSIs.

Selective Smoothed Finite Element Method

The paper examines three selective schemes for the smoothed finite element method （SFEM） which was formulated by incorporating a cell-wise strain smoothing operation into the standard compatible finite element method （FEM）. These selective SFEM schemes were formulated based on three selective integration FEM schemes with similar properties found between the number of smoothing cells in the SFEM and the number of Gaussian integration points in the FEM. Both scheme 1 and scheme 2 are free of nearly incompressible locking, but scheme 2 is more general and gives better results than scheme 1. In addition, scheme 2 can be applied to anisotropic and nonlinear situations, while scheme 1 can only be applied to isotropic and linear situations. Scheme 3 is free of shear locking. This scheme can be applied to plate and shell problems. Results of the numerical study show that the selective SFEM schemes give more accurate results than the FEM schemes.

A stable CS-FEM for the static and seismic stability of a single square tunnel in the soil where the shear strength increases linearly with depth

A numerical procedure using a stable cell-based smoothed finite element method(CS-FEM)is presented for estimation of stability of a square tunnel in the soil where the shear strength increases linearly with depth.The kinematically admissible displacement fields are approximated by uniform quadrilateral elements in conjunction with the strain smoothing technique,eliminating volumetric locking issues and the singularity associated with the MohreCoulomb model.First,a rich set of simulations was performed to compute the static stability of a square tunnel with different geometries and soil conditions.The presented results are in excellent agreement with the upper and lower bound solutions using the standard finite element method(FEM).The stability charts and tables are given for practical use in the tunnel design,along with a newly proposed formulation for predicting the undrained stability of a single square tunnel.Second,the seismic stability number was computed using the present numerical approach.Numerical results reveal that the seismic stability number reduces with an increasing value of the horizontal seismic acceleration(a_(h)),for both cases of the weightless soil and the soil with unit weight.Third,the link between the static and seismic stability numbers is described using corrective factors that represent reductions in the tunnel stability due to seismic loadings.It is shown from the numerical results that the corrective factor becomes larger as the unit weight of soil mass increases;however,the degree of the reduction in seismic stability number tends to reduce for the case of the homogeneous soil.Furthermore,this advanced numerical procedure is straightforward to extend to three-dimensional(3D)limit analysis and is readily applicable for the calculation of the stability of tunnels in highly anisotropic and heterogeneous soils which are often encountered in practice.

Predicting LTE Throughput Using Traffic Time Series

Throughput prediction is essential for congestion control and LTE network management. In this paper, the autoregressive integrated moving average （ARIMA） model and exponential smoothing model are used to predict the throughput in a single cell and whole region in an LTE network. The experimental results show that these two models perform differently in both scenarios. The ARIMA model is better than the exponential smoothing model for predicting throughput on weekdays in a whole region. The exponential smoothing model is better than the ARIMA model for predicting throughput on weekends in a whole region. The exponential smoothing model is better than the ARIMA model for predicting throughput in a single cell. In these two LTE network scenarios, throughput prediction based on traffic time series leads to more efficient resource management and better QoS.

Numerical modeling of sediment transport based on unsteady and steady flows by incompressible smoothed particle hydrodynamics method

The purpose of the present paper is to introduce a simple two-part multi-phase model for the sediment transport problems based on the incompressible smoothed particle hydrodynamics（ISPH） method. The proposed model simulates the movement of sediment particles in two parts. The sediment particles are classified into three categories, including the motionless particles, moving particles behave like a rigid body, and moving particles with a pseudo fluid behavior. The criterion for the classification of sediment particles is the Bingham rheological model. Verification of the present model is performed by simulation of the dam break waves on movable beds with different conditions and the bed scouring under steady flow condition. Comparison of the present model results, the experimental data and available numerical results show that it has good ability to simulate flow pattern and sediment transport.