The numerical implementation of a recently developed thermomechanical constitutive model for fine-grained soils based on hyperelasticity-hyperplasticity theory (Golchin et al. 2020), is presented. A new unconventional implicit stress return mapping algorithm, compatible with elasticity derived from Gibbs (complementary) energy potential, in strain invariant space, is designed and the consistent tangent operator for use in boundary value problems (such as in the finite element method) is derived. It is shown that the rate of convergence of the stress integration algorithm is quadratic. The numerical results are in good agreement with available data from thermomechanical element tests found in literature.

@en

In the Netherlands and other countries around the world there are extensive canals and rivers with banks supported by walls. Especially in the historic inner cities of the Netherlands, quay walls currently need to be replaced. During construction, temporary valuable underground space is created that can be used to accommodate energy installations for the benefit of the surrounding buildings. The quay wall itself can also be a thermal active structure relying on a combination of aquatic and geothermal energy. Two different designs for quay walls incorporating thermal energy infrastructure have been designed: (i) based around a sheet pile; and (ii) based around a concrete pile with a ‘L’ wall on top. An evaluation of the energy, mechanics, finances and CO2 saving is presented for a case study in this paper. It is seen that for the case studies elaborated, both walls can supply sufficient thermal energy for the houses along the canals improved by good use of surface water, reduce CO2 emissions by more than 1 tonne per meter length and can be financially positive. The influence on the mechanical performance of the energy quay wall (settlements, bearing capacity and bending) has been evaluated and found to be reassuringly small.@en

In geotechnical engineering, proper design of retaining structures is of great importance, since failure of these structures can lead to catastrophic consequences. Nowadays, the finite element method is seen as a reliable numerical technique to analyze soil behaviour and is widely used to assess the interaction be-tween soil and rigid structures. However, a disadvantage of this method is the difficulty of simulating contact between separate bodies. Because of this, the event of a slope failing and colliding with a rigid body cannot be analyzed, so that the additional forces acting against the rigid body caused by the motion of the ground are neglected. With the recent development of the material point method (MPM), this limitation has been over-come and problems involving large deformation and multiple bodies in contact can be analyzed. In this paper, the effect of a landslide colliding with a rigid wall has been studied, and multiple initial conditions have been considered in order to identify the critical case.@en