Hi everyone,
In my element routine, I am solving a few nonlinear equations locally using Newton's method. In the case that the local Newton doesn't converge or if any other error occurs at the local level, I raise a flag by setting rmeas = 6 and forcing Feap to return from the element file.
I am also trying to output Stre values using Tplot. Feap outputs Stre values when the time command is used.
.pcd
batch
dt,,5.0
AUTO,MATE,0,0,60.0
loop,infi
time,,60.0
Loop iter 30
utang,,1
Next iter
next
endSo I started my analysis using dt,,5.0 and during the run time, when a Nan appeared at local level, a new attempt with a reduced time step (0.85*dt/rmeas = 0.70833) was started and that step converged well locally and globally. Everything was good till here. For the next step, to update the time step, Time command is used. Here, interestingly FEAP is using dt = 5, to compute the Stre values for Tplot. Shouldn't the dt here be 0.70833, the dt_old ?
After Time, in the Utang,,1, Feap did use the correct dt (dt_new = dt_old/rmeas, I had rmeas default to 0.9, when it is not 6) So the new dt was 0.70833/0.9
As a result of the wrong dt being used during the Time command, the Tplot output of stre values is incorrect. Can you help me out here!
I have attached the relevant output of the terminal window below. Thank you
*Command 7 * next iter v: 2.00 6.00 2.00
t= 0.03 0.00
*Command 6 * utan 0v: 1.00 0.00 0.00
t= 0.03 0.00
dt 5.000, rmeas 0.90
Nan! exiting material, elem 1 gauss 1 rmeas 6.00
exiting element (pos2) 1 gauss 1 rmeas 6.00
Residual norm = 0.0000000E+00 0.0000000E+00 t= 0.03 0.00
*D4TRI WARNING* Reduced diagonal is zero in 7 equations.
Step: 4 Iteration: 2
Condition check: D-max 0.0000E+00; D-min 0.0000E+00; Ratio 0.0000E+00
Maximum no. diagonal digits lost: 1
End Triangular Decomposition t= 0.03 0.00
*WARNING* Zero right-hand-side vector
Number of operations = 210 plus 0 Mega-ops
Time: CPU = 0.03 , System = 0.00
--> SOLVE AT 32.30 Mflops. Time= 0.00
Energy convergence test
Initial = 1.509213712200664E+02 Current = 0.000000000000000E+00
Relative = 0.000000000000000E+00 Tolerance = 1.000000000000000E-16
*Command 7 * next iter v: 30.0 6.00 3.00
t= 0.03 0.00
Computing solution at time 1.5000E+01: Total proportional load 1.5000E-01
Computing solution at time 1.5708E+01: Total proportional load 1.5708E-01
*Command 6 * utan 0v: 1.00 0.00 0.00
t= 0.03 0.00
dt 0.708, rmeas 0.90
Residual norm = 4.1266482E+02 1.0000000E+00 t= 0.03 0.00
Condition check: D-max 0.4104E+05; D-min 0.3221E+05; Ratio 0.1274E+01
Maximum no. diagonal digits lost: 1
End Triangular Decomposition t= 0.03 0.00
Number of operations = 210 plus 0 Mega-ops
Time: CPU = 0.03 , System = 0.00
--> SOLVE AT 42.01 Mflops. Time= 0.00
Energy convergence test
Initial = 2.143518501238372E+00 Current = 2.143518501238372E+00
Relative = 1.000000000000000E+00 Tolerance = 1.000000000000000E-16
*Command 7 * next iter v: 1.00 6.00 1.00
t= 0.03 0.00
*Command 6 * utan 0v: 1.00 0.00 0.00
t= 0.03 0.00
dt 0.708, rmeas 0.90
Residual norm = 4.0199588E+01 9.7414624E-02 t= 0.04 0.00
Condition check: D-max 0.2575E+05; D-min 0.1684E+05; Ratio 0.1529E+01
Maximum no. diagonal digits lost: 1
End Triangular Decomposition t= 0.04 0.00
Number of operations = 210 plus 0 Mega-ops
Time: CPU = 0.04 , System = 0.00
--> SOLVE AT 52.51 Mflops. Time= 0.00
Energy convergence test
Initial = 2.143518501238372E+00 Current = 2.162386395299861E-02
Relative = 1.008802300540251E-02 Tolerance = 1.000000000000000E-16
*Command 7 * next iter v: 2.00 6.00 2.00
t= 0.04 0.00
*Command 6 * utan 0v: 1.00 0.00 0.00
t= 0.04 0.00
dt 0.708, rmeas 0.90
Residual norm = 9.4240340E-03 2.2837018E-05 t= 0.04 0.00
Condition check: D-max 0.2602E+05; D-min 0.1680E+05; Ratio 0.1549E+01
Maximum no. diagonal digits lost: 1
End Triangular Decomposition t= 0.04 0.00
Number of operations = 210 plus 0 Mega-ops
Time: CPU = 0.04 , System = 0.00
--> SOLVE AT 42.01 Mflops. Time= 0.00
Energy convergence test
Initial = 2.143518501238372E+00 Current = 1.157288234726131E-09
Relative = 5.399012110497445E-10 Tolerance = 1.000000000000000E-16
*Command 7 * next iter v: 3.00 6.00 3.00
t= 0.04 0.00
*Command 6 * utan 0v: 1.00 0.00 0.00
t= 0.04 0.00
dt 0.708, rmeas 0.90
Residual norm = 1.7817941E-09 4.3177756E-12 t= 0.04 0.00
Condition check: D-max 0.2602E+05; D-min 0.1680E+05; Ratio 0.1549E+01
Maximum no. diagonal digits lost: 1
End Triangular Decomposition t= 0.04 0.00
Number of operations = 210 plus 0 Mega-ops
Time: CPU = 0.04 , System = 0.00
--> SOLVE AT 46.67 Mflops. Time= 0.00
Energy convergence test
Initial = 2.143518501238372E+00 Current = 4.240453396619205E-23
Relative = 1.978267691260594E-23 Tolerance = 1.000000000000000E-16
*Command 7 * next iter v: 30.0 6.00 4.00
t= 0.04 0.00
*Command 8 * next v: 1.00 4.00 1.00
t= 0.04 0.00
*Command 4 * time 0v: 60.0 0.00 0.00
t= 0.04 0.00
dt 5.000, rmeas 0.90
Computing solution at time 1.6495E+01: Total proportional load 1.6495E-01
*Command 5 * loop iter v: 30.0 8.00 0.00
t= 0.04 0.00
*Command 6 * utan 0v: 1.00 0.00 0.00
t= 0.04 0.00
dt 0.787, rmeas 0.90
Residual norm = 7.5256293E+02 1.0000000E+00 t= 0.04 0.00
*D4TRI WARNING* Sign of diagonal changed when reducing 2 equations.
Step: 5 Iteration: 0
Condition check: D-max 0.8777E+05; D-min 0.4084E+04; Ratio 0.2149E+02
Maximum no. diagonal digits lost: 1
End Triangular Decomposition t= 0.04 0.00
Number of operations = 210 plus 0 Mega-ops
Time: CPU = 0.04 , System = 0.00
--> SOLVE AT 38.19 Mflops. Time= 0.00
Energy convergence test
Initial = 7.369135282737404E+00 Current = 7.369135282737404E+00
Relative = 1.000000000000000E+00 Tolerance = 1.000000000000000E-16
