Non-deterministic `AssertionError: !(any(isnan, component_data))` when running the `von_misses_truss` example with ONSAS v0.4.6 on Julia 1.11

[DilumAluthge]

I am able to reproduce this with ONSAS v0.4.6 on Julia 1.11.6.

MWE

Open Julia 1.11.6, and run the following code:

import Pkg

Pkg.add(name = "ONSAS", version = v"0.4.6")

Close Julia, and then open a new Julia session, and run the following code:

import ONSAS

julia_cmd = Base.julia_cmd()
script = pkgdir(ONSAS, "examples", "von_misses_truss", "von_misses_truss.jl")
cmd = `"$julia_cmd" "$script"`

for i in 1:100
    @info "Starting iteration $i"
    run(cmd)
    @info "Finished iteration $i"
end

In my experience, I can make it up to e.g. 5-10 iterations, but eventually I always get the failure, which looks like this:

Test Summary:                                                    | Pass  Total  Time
Analytic and numeric load factors case: RotatedEngineeringStrain |    1      1  0.1s
ERROR: LoadError: AssertionError: !(any(isnan, component_data))
Stacktrace:
  [1] write_cell_data(vtk::VTKMeshFile{WriteVTK.DatasetFile}, celldata::Vector{Matrix{Float64}}, name::String; component_names::Vector{String}, kwargs::@Kwargs{})
    @ ONSAS.VTK ~/.julia/packages/ONSAS/49upX/src/Interfaces/VTK.jl:137
  [2] write_cell_data
    @ ~/.julia/packages/ONSAS/49upX/src/Interfaces/VTK.jl:128 [inlined]
  [3] (::ONSAS.VTK.var"#24#25"{Dict{Symbol, Vector{Matrix{Float64}}}, Dict{Symbol, Vector{Float64}}})(vtx::VTKMeshFile{WriteVTK.DatasetFile})
    @ ONSAS.VTK ~/.julia/packages/ONSAS/49upX/src/Interfaces/VTK.jl:232
  [4] VTKMeshFile(::ONSAS.VTK.var"#24#25"{Dict{Symbol, Vector{Matrix{Float64}}}, Dict{Symbol, Vector{Float64}}}, ::String, ::Vararg{Any}; kwargs::@Kwargs{append::Bool})
    @ ONSAS.VTK ~/.julia/packages/ONSAS/49upX/src/Interfaces/VTK.jl:40
  [5] VTKMeshFile
    @ ~/.julia/packages/ONSAS/49upX/src/Interfaces/VTK.jl:37 [inlined]
  [6] #write_vtk#23
    @ ~/.julia/packages/ONSAS/49upX/src/Interfaces/VTK.jl:225 [inlined]
  [7] write_vtk
    @ ~/.julia/packages/ONSAS/49upX/src/Interfaces/VTK.jl:215 [inlined]
  [8] write_vtk(sol::Solution{StaticState, NonLinearStaticAnalysis{Structure{3, Mesh{3, Node{3, Float64}, Truss{3, GreenStrain, Float64, Node{3, Float64}, StaticArraysCore.SVector{2, Node{3, Float64}}}, AbstractFace, Nothing}, AbstractMaterial, Truss{3, GreenStrain, Float64, Node{3, Float64}, StaticArraysCore.SVector{2, Node{3, Float64}}}, AbstractBoundaryCondition, AbstractBoundaryCondition}, Float64, Vector{Float64}}, NewtonRaphson}, base_filename::String; fields::Vector{Symbol}, append::Bool)
    @ ONSAS.VTK ~/.julia/packages/ONSAS/49upX/src/Interfaces/VTK.jl:260
  [9] write_vtk(sol::Solution{StaticState, NonLinearStaticAnalysis{Structure{3, Mesh{3, Node{3, Float64}, Truss{3, GreenStrain, Float64, Node{3, Float64}, StaticArraysCore.SVector{2, Node{3, Float64}}}, AbstractFace, Nothing}, AbstractMaterial, Truss{3, GreenStrain, Float64, Node{3, Float64}, StaticArraysCore.SVector{2, Node{3, Float64}}}, AbstractBoundaryCondition, AbstractBoundaryCondition}, Float64, Vector{Float64}}, NewtonRaphson}, base_filename::String)
    @ ONSAS.VTK ~/.julia/packages/ONSAS/49upX/src/Interfaces/VTK.jl:251
 [10] run()
    @ Main ~/.julia/packages/ONSAS/49upX/examples/von_misses_truss/von_misses_truss.jl:135
 [11] top-level scope
    @ ~/.julia/packages/ONSAS/49upX/examples/von_misses_truss/von_misses_truss.jl:140
in expression starting at /Users/dilum/.julia/packages/ONSAS/49upX/examples/von_misses_truss/von_misses_truss.jl:140
ERROR: failed process: Process(`/path/to/julia/1.11.6/bin/julia -C native -J/path/to/julia/1.11.6/lib/julia/sys.dylib -g1 /Users/dilum/.julia/packages/ONSAS/49upX/examples/von_misses_truss/von_misses_truss.jl`, ProcessExited(1)) [1]

Stacktrace:
 [1] pipeline_error
   @ ./process.jl:598 [inlined]
 [2] run(::Cmd; wait::Bool)
   @ Base ./process.jl:513
 [3] run(::Cmd)
   @ Base ./process.jl:510
 [4] top-level scope
   @ ./REPL[5]:3

[mvanzulli]

So it seems we have some numerical instability in that example. One option is to remove the assertion so if a Nan is passed along then the solver errors but not the write vtk. Will remove the check from there for now.