arguments
fname=strip.msh
degree=1
Δuload=0.004999999999999999mm
nsteps=11
E=70GPa
nu=0.2
σyield=243MPa
hardening=2.25GPa
referencedata=zienkiewicz.csv
testing=False
gmsh
loaded 2d gmsh topology consisting of #328 elements
- volume groups: interior #328
- boundary groups: hsymmetry #6, vsymmetry #14, load #5
solve
solving 356 dof system to machine precision using arnoldi solver
constructing direct preconditioner
residual decreased by 14.9 orders using 0 krylov vectors
residual decreased by 0.2 orders using 1 krylov vectors
solver returned with residual 8e-13
step 1/11
NutilsDeprecationWarning: function.sqrt is deprecated; use numpy.sqrt instead
In /usr/local/lib/python3.11/dist-packages/nutils/function.py:1226
NutilsDeprecationWarning: function.inverse is deprecated; use numpy.linalg.inv instead
In /usr/local/lib/python3.11/dist-packages/nutils/function.py:1226
NutilsDeprecationWarning: function.greater is deprecated; use numpy.greater or the > operator instead
In /usr/local/lib/python3.11/dist-packages/nutils/function.py:1226
newton
iter 0
loop _sample_X 62%
rmap 1/20
rmap residual = 0.36924679291437684
rmap 2/20
rmap residual = 0.0577093254594627
rmap 3/20
rmap residual = 0.0009115452897854283
rmap 4/20
rmap residual = 4.053733565830759e-06
rmap 5/20
rmap residual = 5.803540863421531e-09
rmap 1/20
rmap residual = 0.36924679291437684
rmap 2/20
rmap residual = 0.0577093254594627
rmap 3/20
rmap residual = 0.0009115452897854283
rmap 4/20
rmap residual = 4.053733565830759e-06
rmap 5/20
rmap residual = 5.803540863421531e-09
rmap 1/20
rmap residual = 0.36924679291437684
rmap 2/20
rmap residual = 0.0577093254594627
rmap 3/20
rmap residual = 0.0009115452897854283
rmap 4/20
rmap residual = 4.053733565830759e-06
rmap 5/20
rmap residual = 5.803540863421531e-09
loop _sample_X 73%
rmap 1/20
rmap residual = 0.03494458996042111
rmap 2/20
rmap residual = 0.000524271305684226
rmap 3/20
rmap residual = 1.7752035121016353e-07
rmap 4/20
rmap residual = 1.3515167103687189e-11
rmap 1/20
rmap residual = 0.03494458996042111
rmap 2/20
rmap residual = 0.000524271305684226
rmap 3/20
rmap residual = 1.7752035121016353e-07
rmap 4/20
rmap residual = 1.3515167103687189e-11
rmap 1/20
rmap residual = 0.03494458996042111
rmap 2/20
rmap residual = 0.000524271305684226
rmap 3/20
rmap residual = 1.7752035121016353e-07
rmap 4/20
rmap residual = 1.3515167103687189e-11
loop _sample_X 81%
rmap 1/20
rmap residual = 0.17059004666396904
rmap 2/20
rmap residual = 0.012217587660233401
rmap 3/20
rmap residual = 4.5525090080967413e-05
rmap 4/20
rmap residual = 9.611667300952833e-08
rmap 5/20
rmap residual = 8.513752311591405e-11
rmap 1/20
rmap residual = 0.17059004666396904
rmap 2/20
rmap residual = 0.012217587660233401
rmap 3/20
rmap residual = 4.5525090080967413e-05
rmap 4/20
rmap residual = 9.611667300952833e-08
rmap 5/20
rmap residual = 8.513752311591405e-11
rmap 1/20
rmap residual = 0.17059004666396904
rmap 2/20
rmap residual = 0.012217587660233401
rmap 3/20
rmap residual = 4.5525090080967413e-05
rmap 4/20
rmap residual = 9.611667300952833e-08
rmap 5/20
rmap residual = 8.513752311591405e-11
iter 1 (0%)
solve
solving 356 dof system to tolerance 1e+00 using arnoldi solver
constructing direct preconditioner
residual decreased by 14.5 orders using 0 krylov vectors
solver returned with residual 3e-12
loop _sample_X 62%
rmap 1/20
rmap residual = 1.4129548166806147
rmap 2/20
rmap residual = 0.6715264478217722
rmap 3/20
rmap residual = 0.07058498406504249
rmap 4/20
rmap residual = 0.002620576808288801
rmap 5/20
rmap residual = 1.6244913291964725e-05
rmap 6/20
rmap residual = 3.965055771777252e-08
rmap 7/20
rmap residual = 1.487986432796612e-10
rmap 1/20
rmap residual = 1.4129548166806147
rmap 2/20
rmap residual = 0.6715264478217722
rmap 3/20
rmap residual = 0.07058498406504249
rmap 4/20
rmap residual = 0.002620576808288801
rmap 5/20
rmap residual = 1.6244913291964725e-05
rmap 6/20
rmap residual = 3.965055771777252e-08
rmap 7/20
rmap residual = 1.487986432796612e-10
rmap 1/20
rmap residual = 1.4129548166806147
rmap 2/20
rmap residual = 0.6715264478217722
rmap 3/20
rmap residual = 0.07058498406504249
rmap 4/20
rmap residual = 0.002620576808288801
rmap 5/20
rmap residual = 1.6244913291964725e-05
rmap 6/20
rmap residual = 3.965055771777252e-08
rmap 7/20
rmap residual = 1.487986432796612e-10
loop _sample_X 73%
rmap 1/20
rmap residual = 0.6994190001982555
rmap 2/20
rmap residual = 0.1096905118114865
rmap 3/20
rmap residual = 0.002058829116860351
rmap 4/20
rmap residual = 8.675800768209244e-06
rmap 5/20
rmap residual = 4.517051176419382e-09
rmap 1/20
rmap residual = 0.6994190001982555
rmap 2/20
rmap residual = 0.1096905118114865
rmap 3/20
rmap residual = 0.002058829116860351
rmap 4/20
rmap residual = 8.675800768209244e-06
rmap 5/20
rmap residual = 4.517051176419382e-09
rmap 1/20
rmap residual = 0.6994190001982555
rmap 2/20
rmap residual = 0.1096905118114865
rmap 3/20
rmap residual = 0.002058829116860351
rmap 4/20
rmap residual = 8.675800768209244e-06
rmap 5/20
rmap residual = 4.517051176419382e-09
loop _sample_X 78%
rmap 1/20
rmap residual = 0.45756031087925575
rmap 2/20
rmap residual = 0.04177240681320414
rmap 3/20
rmap residual = 0.0002985336859977726
rmap 4/20
rmap residual = 7.43306777751856e-07
rmap 5/20
rmap residual = 5.170126667152799e-10
rmap 1/20
rmap residual = 0.45756031087925575
rmap 2/20
rmap residual = 0.04177240681320414
rmap 3/20
rmap residual = 0.0002985336859977726
rmap 4/20
rmap residual = 7.43306777751856e-07
rmap 5/20
rmap residual = 5.170126667152799e-10
rmap 1/20
rmap residual = 0.45756031087925575
rmap 2/20
rmap residual = 0.04177240681320414
rmap 3/20
rmap residual = 0.0002985336859977726
rmap 4/20
rmap residual = 7.43306777751856e-07
rmap 5/20
rmap residual = 5.170126667152799e-10
loop _sample_X 81%
rmap 1/20
rmap residual = 0.967094219061837
rmap 2/20
rmap residual = 0.2794400608520845
rmap 3/20
rmap residual = 0.013934450010759493
rmap 4/20
rmap residual = 0.00014438723010629637
rmap 5/20
rmap residual = 2.974739991923109e-07
rmap 6/20
rmap residual = 7.738348426686146e-10
rmap 1/20
rmap residual = 0.967094219061837
rmap 2/20
rmap residual = 0.2794400608520845
rmap 3/20
rmap residual = 0.013934450010759493
rmap 4/20
rmap residual = 0.00014438723010629637
rmap 5/20
rmap residual = 2.974739991923109e-07
rmap 6/20
rmap residual = 7.738348426686146e-10
rmap 1/20
rmap residual = 0.967094219061837
rmap 2/20
rmap residual = 0.2794400608520845
rmap 3/20
rmap residual = 0.013934450010759493
rmap 4/20
rmap residual = 0.00014438723010629637
rmap 5/20
rmap residual = 2.974739991923109e-07
rmap 6/20
rmap residual = 7.738348426686146e-10
estimated residual minimum at 46% of update vector
estimated residual minimum at 89% of update vector
update accepted at relaxation 0.46479
iter 2 (4%)
solve
solving 356 dof system to tolerance 5e-01 using arnoldi solver
constructing direct preconditioner
residual decreased by 15.2 orders using 0 krylov vectors
solver returned with residual 3e-13
estimated residual minimum at 241% of update vector
update accepted at relaxation 0.41552
iter 3 (6%)
solve
solving 356 dof system to tolerance 3e-01 using arnoldi solver
constructing direct preconditioner
residual decreased by 15.2 orders using 0 krylov vectors
solver returned with residual 2e-13
estimated residual minimum at 120% of update vector
update accepted at relaxation 0.83104
iter 4 (15%)
solve
solving 356 dof system to tolerance 4e-02 using arnoldi solver
constructing direct preconditioner
residual decreased by 15.2 orders using 0 krylov vectors
solver returned with residual 3e-14
estimated residual minimum at 100% of update vector
update accepted at relaxation 1
converged in 4 iterations to residual 3.1e-14
NutilsDeprecationWarning: function.sqrt is deprecated; use numpy.sqrt instead
In /usr/local/lib/python3.11/dist-packages/nutils/function.py:1226
step 2/11
NutilsDeprecationWarning: function.inverse is deprecated; use numpy.linalg.inv instead
In /usr/local/lib/python3.11/dist-packages/nutils/function.py:1226
NutilsDeprecationWarning: function.greater is deprecated; use numpy.greater or the > operator instead
In /usr/local/lib/python3.11/dist-packages/nutils/function.py:1226
newton
iter 0
loop _sample_X 62%
rmap 1/20
rmap residual = 0.4011851505292237
rmap 2/20
rmap residual = 0.06823230456152028
rmap 3/20
rmap residual = 0.001250840168702049
rmap 4/20
rmap residual = 6.0121539269803345e-06
rmap 5/20
rmap residual = 8.437954305605043e-09
rmap 1/20
rmap residual = 0.4011851505292237
rmap 2/20
rmap residual = 0.06823230456152028
rmap 3/20
rmap residual = 0.001250840168702049
rmap 4/20
rmap residual = 6.0121539269803345e-06
rmap 5/20
rmap residual = 8.437954305605043e-09
rmap 1/20
rmap residual = 0.4011851505292237
rmap 2/20
rmap residual = 0.06823230456152028
rmap 3/20
rmap residual = 0.001250840168702049
rmap 4/20
rmap residual = 6.0121539269803345e-06
rmap 5/20
rmap residual = 8.437954305605043e-09
loop _sample_X 73%
rmap 1/20
rmap residual = 0.0702199036908595
rmap 2/20
rmap residual = 0.0021164225802312122
rmap 3/20
rmap residual = 1.6711226133461008e-06
rmap 4/20
rmap residual = 1.0531053659972039e-09
rmap 1/20
rmap residual = 0.0702199036908595
rmap 2/20
rmap residual = 0.0021164225802312122
rmap 3/20
rmap residual = 1.6711226133461008e-06
rmap 4/20
rmap residual = 1.0531053659972039e-09
rmap 1/20
rmap residual = 0.0702199036908595
rmap 2/20
rmap residual = 0.0021164225802312122
rmap 3/20
rmap residual = 1.6711226133461008e-06
rmap 4/20
rmap residual = 1.0531053659972039e-09
loop _sample_X 81%
rmap 1/20
rmap residual = 0.20395102861814213
rmap 2/20
rmap residual = 0.017522557406795414
rmap 3/20
rmap residual = 9.217131836468041e-05
rmap 4/20
rmap residual = 2.3420716619012625e-07
rmap 5/20
rmap residual = 2.4306491665668163e-10
rmap 1/20
rmap residual = 0.20395102861814213
rmap 2/20
rmap residual = 0.017522557406795414
rmap 3/20
rmap residual = 9.217131836468041e-05
rmap 4/20
rmap residual = 2.3420716619012625e-07
rmap 5/20
rmap residual = 2.4306491665668163e-10
rmap 1/20
rmap residual = 0.20395102861814213
rmap 2/20
rmap residual = 0.017522557406795414
rmap 3/20
rmap residual = 9.217131836468041e-05
rmap 4/20
rmap residual = 2.3420716619012625e-07
rmap 5/20
rmap residual = 2.4306491665668163e-10
iter 1 (0%)
solve
solving 356 dof system to tolerance 1e+00 using arnoldi solver
constructing direct preconditioner
residual decreased by 14.6 orders using 0 krylov vectors
solver returned with residual 3e-12
loop _sample_X 62%
rmap 1/20
rmap residual = 1.3081498724072307
rmap 2/20
rmap residual = 0.5802238612836635
rmap 3/20
rmap residual = 0.05490555912927138
rmap 4/20
rmap residual = 0.0016473172692975833
rmap 5/20
rmap residual = 9.275756522223444e-06
rmap 6/20
rmap residual = 2.4960092494488014e-08
rmap 7/20
rmap residual = 9.020091486766619e-11
rmap 1/20
rmap residual = 1.3081498724072307
rmap 2/20
rmap residual = 0.5802238612836635
rmap 3/20
rmap residual = 0.05490555912927138
rmap 4/20
rmap residual = 0.0016473172692975833
rmap 5/20
rmap residual = 9.275756522223444e-06
rmap 6/20
rmap residual = 2.4960092494488014e-08
rmap 7/20
rmap residual = 9.020091486766619e-11
rmap 1/20
rmap residual = 1.3081498724072307
rmap 2/20
rmap residual = 0.5802238612836635
rmap 3/20
rmap residual = 0.05490555912927138
rmap 4/20
rmap residual = 0.0016473172692975833
rmap 5/20
rmap residual = 9.275756522223444e-06
rmap 6/20
rmap residual = 2.4960092494488014e-08
rmap 7/20
rmap residual = 9.020091486766619e-11
loop _sample_X 73%
rmap 1/20
rmap residual = 0.6085048889215479
rmap 2/20
rmap residual = 0.08402669838224323
rmap 3/20
rmap residual = 0.001248682903179151
rmap 4/20
rmap residual = 4.6635796110793005e-06
rmap 5/20
rmap residual = 3.5134854714843127e-09
rmap 1/20
rmap residual = 0.6085048889215479
rmap 2/20
rmap residual = 0.08402669838224323
rmap 3/20
rmap residual = 0.001248682903179151
rmap 4/20
rmap residual = 4.6635796110793005e-06
rmap 5/20
rmap residual = 3.5134854714843127e-09
rmap 1/20
rmap residual = 0.6085048889215479
rmap 2/20
rmap residual = 0.08402669838224323
rmap 3/20
rmap residual = 0.001248682903179151
rmap 4/20
rmap residual = 4.6635796110793005e-06
rmap 5/20
rmap residual = 3.5134854714843127e-09
loop _sample_X 78%
rmap 1/20
rmap residual = 0.39555150036823317
rmap 2/20
rmap residual = 0.029957977390283105
rmap 3/20
rmap residual = 0.00015136804045934633
rmap 4/20
rmap residual = 3.1654015281060113e-07
rmap 5/20
rmap residual = 1.9162197358516036e-10
rmap 1/20
rmap residual = 0.39555150036823317
rmap 2/20
rmap residual = 0.029957977390283105
rmap 3/20
rmap residual = 0.00015136804045934633
rmap 4/20
rmap residual = 3.1654015281060113e-07
rmap 5/20
rmap residual = 1.9162197358516036e-10
rmap 1/20
rmap residual = 0.39555150036823317
rmap 2/20
rmap residual = 0.029957977390283105
rmap 3/20
rmap residual = 0.00015136804045934633
rmap 4/20
rmap residual = 3.1654015281060113e-07
rmap 5/20
rmap residual = 1.9162197358516036e-10
loop _sample_X 81%
rmap 1/20
rmap residual = 0.8709816974522826
rmap 2/20
rmap residual = 0.22962692047591995
rmap 3/20
rmap residual = 0.009797148526775342
rmap 4/20
rmap residual = 8.399142916867144e-05
rmap 5/20
rmap residual = 9.878729317815382e-08
rmap 6/20
rmap residual = 2.447624983762004e-10
rmap 1/20
rmap residual = 0.8709816974522826
rmap 2/20
rmap residual = 0.22962692047591995
rmap 3/20
rmap residual = 0.009797148526775342
rmap 4/20
rmap residual = 8.399142916867144e-05
rmap 5/20
rmap residual = 9.878729317815382e-08
rmap 6/20
rmap residual = 2.447624983762004e-10
rmap 1/20
rmap residual = 0.8709816974522826
rmap 2/20
rmap residual = 0.22962692047591995
rmap 3/20
rmap residual = 0.009797148526775342
rmap 4/20
rmap residual = 8.399142916867144e-05
rmap 5/20
rmap residual = 9.878729317815382e-08
rmap 6/20
rmap residual = 2.447624983762004e-10
estimated residual minimum at 44% of update vector
estimated residual minimum at 93% of update vector
update accepted at relaxation 0.4442
iter 2 (4%)
solve
solving 356 dof system to tolerance 4e-01 using arnoldi solver
constructing direct preconditioner
residual decreased by 15.1 orders using 0 krylov vectors
solver returned with residual 3e-13
estimated residual minimum at 241% of update vector
update accepted at relaxation 0.41428
iter 3 (7%)
solve
solving 356 dof system to tolerance 2e-01 using arnoldi solver
constructing direct preconditioner
residual decreased by 15.1 orders using 0 krylov vectors
solver returned with residual 2e-13
estimated residual minimum at 121% of update vector
update accepted at relaxation 0.82855
iter 4 (15%)
solve
solving 356 dof system to tolerance 4e-02 using arnoldi solver
constructing direct preconditioner
residual decreased by 15.1 orders using 0 krylov vectors
solver returned with residual 3e-14
estimated residual minimum at 100% of update vector
update accepted at relaxation 1.0
converged in 4 iterations to residual 4.4e-13