The instabilities of soil specimens in laboratory or soil made geotechnical structures in field are always numerically simulated by the classical continuum mechanics-based constitutive models with finite element method. However, finite element mesh dependency problems are inevitably encountered when the strain localized failure occurs especially in the post-bifurcation regime. In this paper, an attempt is made to summarize several main numerical regularization techniques used in alleviating the mesh dependency problems, i.e., viscosity theory, nonlocal theory, high-order gradient and micropolar theory. Their fundamentals as well as the advantages and limitations are presented, based on which the combinations of two or more regularization techniques are also suggested. For all the regularization techniques, at least one implicit or explicit parameter with length scale is necessary to preserve the ellipticity of the partial differential governing equations. It is worth noting that, however, the physical meanings and their relations between the length parameters in different regularization techniques are still an open question, and need to be further studied. Therefore, the micropolar theory or its combinations with other numerical methods are promising in the future.

A large number of experimental studies on sand–clay mixtures are well documented in the literature; however, the preparation protocol is rarely clearly detailed or varies a lot according to the authors. Variations in the preparation technique obviously increases the challenge of comparing different test results. As a consequence, sample preparation for sand–clay mixtures should be kept as constant as possible to ensure homogeneity and uniformity of samples and limit result variability. This paper develops a detailed procedure on how to prepare sand–clay samples for interface direct shear tests. Sand–clay mixtures are prepared with Fontainebleau sand, kaolinite clay and distilled water by the S1 (sand–water–clay) protocol. The uniformity of the reconstituted specimens is assessed by measuring the water content and density on three slices from the top to the bottom across the specimens. The repeatability of the samples is checked with oedometer and interface direct shear tests. This sample preparation procedure can be used for preparing sand–clay mixture for interface direct shear tests to investigate the influence of clay content or other effects (e.g., temperature) on the mechanical behavior of soil–structure interface. It has demonstrated great performance in preparing samples with good homogeneity and shape, compared to other traditional reconstitution techniques. With the sample preparation procedure, we can obtain repeatable test results as well.

The mechanical behaviors of Martian regolith-structure interfaces are of great significance for the design of rover, development of excavation tools, and construction of infrastructure in Mars exploration. This paper presents an experimental investigation on the properties of a Martian regolith simulant (JEZ-1) through one-dimensional oedometer test, direct shear test, and interface direct shear tests between JEZ-1 and steel plates with different roughness. Oedometer result reveals that the compression and swelling indexes of the JEZ-1 are quite low, thus it is a less compressible and lower swelling soil. The cohesion and adhesion of JEZ-1 are lower than 5 kPa. The values of the internal friction angle range from 39.7° to 40.6°, and the interface friction angles are 16.7° to 36.2° for the smooth and rough interface. Furthermore, the direct shear and interface direct shear results indicate that the interface friction angles are lower than the internal friction angles of JEZ-1 and increase close to the internal friction angles with increasing interface roughness.

The expansion induced by sulfate attack on cement-treated aggregates (SACA) is a well-known problem that can be solved. It causes obvious heaves in road bases and railway subgrades. In this paper, the effects of the sodium sulfate content, cement content, degree of compaction, sulfate types, attack types, aluminum ion supply, concentration of curing sulfate solution, and temperature on the expansion behavior induced by SACA were investigated over 60 days in the laboratory. Based on the Sobol sensitivity analysis method, the influence of the sensitivity of each factor on the expansion was quantitatively analyzed, and the dominant factor of expansion was proposed. Results show that sulfate content is the domain factor of expansion that is induced by SACA, and it presents a logarithmic function relationship with strain. The 0.5% sodium sulfate content is the minimum sulfate content which causes the expansion that is induced by SACA. When the sulfate content is less than 1%, the expansion induced by SACA is minor. When the sulfate content is between 1% and 3%, the expansion behavior is expressed in four stages as follows: rapid strain increase, followed by a short stagnation period, then a significant strain increase and, finally, constant strain. When the sulfate content is greater than 5%, there are two stages comprising the expansion behavior as follows: the rapid strain increases and constant strain occurs. Greater sulfate content, greater degree of compaction, and lower temperature have positive effects on the expansion induced by SACA. The cement content does not have a consistent effect on expansion behavior. Compared with a sodium sulfate attack, both the reaction rate and expansion of cement-treated aggregates that are attacked by gypsum are smaller, and the attack period is also longer. When the sulfate content is greater than 1%, the addition of kaolin promotes the progression of the expansion induced by SACA. A small amount of water is sufficient for the demand for the sulfate attack. When the sulfate content is at a certain level, the expansion induced by SACA that is under external attack is much smaller than the expansion that is under internal attack. This study is expected to serve as a reference for future research on the mechanics of SACA, and it attempts to provide theoretical support for amending expansions that are induced by SACA.

In traditional finite element failure analyses of geotechnical structures, the micro grain rotations cannot be modelled and numerical solutions are mesh dependent. In this study, a user element including rotational degree of freedom has been developed based on micropolar theory (Cosserat theory), then an enhanced non-associated sand model is calibrated with laboratory data and used to model the plane strain tests. The simulated results demonstrate the polarized model is able to model reasonably the sand behavior as well as the grain rotations in the localized region. What’s more, with this enhanced model, the mesh independent numerical solutions in terms of mechanical responses, shear bands thickness and orientations have been obtained. Article highlights: In failure analysis of geostructures, significant rotations of soil grains have been observed to occur in the strain localized regions, but the current commercial Finite Element tools cannot model the micro rotations. Therefore, a user defined element must be developed to include the rotational degree of freedom. The micropolar approach is proven to be effective to model the grain rations in present paper. More suitable than other classical soil or sand constitutive models, the selected non-associated sand model in present paper is capable of describing well the contraction and shear dilatancy behaviors of sand. Then the model has been enhanced by means of micropolar technique, in this way, the reasonable strain localization phenomena in laboratory tests could be predicted well. There are always the mesh dependent problems for traditional simulations of the strain localization phenomena in finite element analysis. It can be found in present paper that the mesh independent numerical solutions are obtained by means of micropolar technique. Furthermore, the micropolar approach can obviously improve convergence difficulties in finite element analyses.