Skip to content

December 2024

The next upcoming release will include:

Graphical User Interface

We are proud to announce that numgeo now integrates into the Pre-Processor GiD! See GUI

Features:

  • Initial state
    • A new initial stress type type=full is now available allowing to prescribe all components of the stress tensor individually
  • Solver
    • It is now possible to omit the use of a linear equation solver (see *Step definition).
  • Early termination:
    • Define target stresses/strains for specific elements that lead to early termination of the step once they are exceeded. See Early termination
  • IMEX scheme:
    • A so-called mixed IMplicit-EXplicit (IMEX) scheme is available. It allows to mix the two time integrations spatially. See *Time integration
  • New boundary condition:
    • You can now specify a hydrostatic distribution with changing zero-level with time using the new boundary condition type=moving-hydrostatic, see *Boundary. A typical use case is to prescribe a changing water level to some parts of your model.
  • Improved contact behaviour of quadratic tetrahedra (10 nodes) and hexahedra (20 nodes) serendipity elements.
  • New elements:
  • New Bishop effective stress model:
    • A new approach to account for the influence of suction on the effective stress. The new approach considers the Bishop variable \(\chi\) as a simple power function of the effective degree of saturation, i.e. \(\chi=(S^e)^\kappa\), see Bishop effective stress. By appropriate choice of \(\kappa\) the modified Bishop stress approach from the software ZSoil can be recovered.
  • Line search:
    • In simulations involving multiple active physical fields, such as a consolidation analysis using coupled two-phase elements, the scaling factor \(\lambda\) is calculated separately for each active field.
  • Node set definition based on the geometry using circles/cylinders
  • Constitutive models
    • The Karlsruhe hydraulic model is now formulated in terms of the degree of saturation \(S^w\) and takes two residual degrees of saturation as input parameters: \(S_{r,d}\) for the main drying curve and \(S_{r,w}\) for the main wetting curve
    • Mohr-Coulomb implementation using the exact yield surface and an elastic tangent stiffness, see Mohr-Coulomb-2

Manuals:

General modifications:

  • "Static" steps now calculate the steady state solution in coupled simulations (e.g. using coupled u-p elements)

  • Behaviour of print output changed: The first 300 files will remain open throughout the simulation, while any additional files will be closed after each write operation and reopened when needed. Find more about it here

  • User files are compiled using ifx, since support of ifort stopped in recent OneAPI versions

Resolved bugs:

  • Bugfix in Karlsruhe Hydraulic model for update of the degree of saturation
  • In previous versions of numgeo, only one optional argument - maxiter, miniter, or inc - from the step definition was correctly recognized, and the arguments had to be provided in a specific order. This issue has now been resolved, allowing these arguments to be passed in any order.
  • Bugfix in calculation of friction in case of 3D