Skip to content

Upcoming

The release notes for the upcoming release are:

New features:

  • Truss elements: We are happy to announce that with this release numgeo now extends its finite element library by truss elements! Have a look at the Reference Manual, the Theory Manual or the Benchmark Manual for more information.

  • Beam elements: We are happy to announce that with this release numgeo now extends its finite element library by beam elements! Have a look at the Reference Manual, the Theory Manual or the Benchmark Manual for more information.

  • Flexible Boundary Conditions: You can now define Dirichlet/essential boundary conditions which change their status (value and active/inactive) with space and time within a calculation step. For more information see *Flexible-Boundary and the corresponding section in the Theory Manual.

  • Model change: You can now change a material assignment mid analysis! To find more about this option, have a look at the corresponding Sections in the Reference Manual or the Theory Manual.

  • Embedded region: You can now embed elements into a host mesh using the new *Embedded region feature! For details, see the corresponding Sections in the Reference Manual or the Theory Manual.

  • GiD output: With this release, field output can now also be requested for the software GiD, supplementing the previously available ParaView (VTK) output. This addition marks an important step toward the full integration of numgeo into the GiD software, now covering the entire workflow, including pre-processing, calculation, and post-processing. This enhancement significantly improves usability by allowing seamless data transfer between numgeo and GiD, facilitating a more streamlined and comprehensive analysis process.

  • Improved Hypoplasticity + ISA: the new release features an improved implementation of the Hypo-ISA model in numgeo. This update features an Euler-Richardson integration scheme, optional projection of inadmissible stress states on the Matsuoka-Nakai yield surface and more. For more details check out the Theory Manual and the Reference Manual

General modifications:

  • *Initial conditions, type = stress not mandatory anymore. Definition of an initial stress field is not mandatory anymore. It is the users responsibility to prescribe (or not) an initial (effective) stress field.

  • User boundary can now consider type (velocity/acceleration) similarly to regular boundary conditions, see UBoundary. In addition, the arguments of the user boundary suboutine have been adjusted, such that the full name of the degree of freedom is passed (with length 2, e.g. U1), allowing to define multiple degrees of freedom in one subroutine by defining a case.

  • New feature for SANISAND-MSf: Users can now choose between two formulations for \(M_b\)

  • If the order of slave/master set in the definition of Equations is incorrect (i.e. the master set has more nodes than the slave set) the sets are now automatically reordered instead of issuing an error message

  • The *Solid Section keyword to assign a material model to an element set is now renamed to *Assign Material (backward compatibility is ensured)

  • Performance

    • Considerable speed-up of some basic sorting algorithms
    • Considerable speed-up of contact routines
    • Speed-up of element assembly, especially in 3D simulations

  • Simplified the Modified Cam Clay material interface (removed one algorithmic parameter)

  • Two new element types for partially saturated soils:

    • `u8p8-stab-ax: 8 node, axisymmetric rectangular element
    • `u8p8-stab-ax-red: 8 node, axisymmetric rectangular element with reduced (4 point) integration

  • VTK/Paraview output: nodal averaged strains are now displayed as symmetric tensors (like stress)

Resolved bugs:

  • Amplitudes are now evaluated once per increment. This resolves an error in cases where the same lab-amplitude was referenced multiple times in one step. This issue only occurred for lab-amplitudes with an early termination criterion, other amplitudes were not affected.

  • An error is issued if *Surface, type=node is defined within the assembly, which is not allowed

  • Fixed potential of NaN in some cases if no gravity is defined by the user

  • Error message if initial conditions of same type are defined multiple times

  • Fixed extrapolation of integration point values to nodes for 10-node tetrahedral elements

Manuals: