Note
Go to the end to download the full example code.
Example Thermal Transient Analysis#
This example shows how you can use PyMAPDL to input a time dependent temperature table to vary the temperature at a beam. This uses convection loads with independently varying convection coefficient and bulk temperature.
Example adapted from: https://www.simutechgroup.com/tips-and-tricks/fea-articles/97-fea-tips-tricks-thermal-transient
Thanks SimuTech!
import matplotlib.pyplot as plt
import numpy as np
from ansys.mapdl.core import launch_mapdl
mapdl = launch_mapdl(loglevel="ERROR")
mapdl.clear()
mapdl.prep7()
# Material properties-- 1020 steel in imperial
mapdl.units("BIN") # U.S. Customary system using inches (in, lbf*s2/in, s, °F).
mapdl.mp("EX", 1, 30023280.0)
mapdl.mp("NUXY", 1, 0.290000000)
mapdl.mp("ALPX", 1, 8.388888889e-06)
mapdl.mp("DENS", 1, 7.346344000e-04)
mapdl.mp("KXX", 1, 6.252196000e-04)
mapdl.mp("C", 1, 38.6334760)
# use a thermal element type
mapdl.et(1, "SOLID70")
1
Geometry and Mesh#
Create a block 5x1x1 inches in size and mesh it
mapdl.block(0, 5, 0, 1, 0, 1)
mapdl.lesize("ALL", 0.2, layer1=1)
mapdl.mshape(0, "3D")
mapdl.mshkey(1)
mapdl.vmesh(1)
mapdl.eplot()
Setup the Solution#
Solve a transient analysis while ramping the load up and down.
Note the solution time commands in the above code fragment. The final TIME is set to 1000 seconds. Time substep size is permitted to range from a minimum of 2 seconds to a maximum of 50 seconds in the DELTIM command. A first substep of 10 seconds is applied. Automatic time substep sizing will vary substeps between the extremes.
A Table Array is used for the time-dependent Convection Coefficient values. Times go in the Zeroth column, while associated Convection Coefficients go in the First column.
mapdl.run("/SOLU")
mapdl.antype(4) # transient analysis
mapdl.trnopt("FULL") # full transient analysis
mapdl.kbc(0) # ramp loads up and down
# Time stepping
end_time = 1500
mapdl.time(end_time) # end time for load step
mapdl.autots("ON") # use automatic time stepping
# setup where the subset time is 10 seconds, time
mapdl.deltim(10, 2, 25) # substep size (seconds)
# -- minimum value shorter than smallest
# time change in the table arrays below
# Create a table of convection times and coefficients and transfer it to MAPDL
my_conv = np.array(
[
[0, 0.001], # start time
[120, 0.001], # end of first "flat" zone
[130, 0.005], # ramps up in 10 seconds
[700, 0.005], # end of second "flat zone
[710, 0.002], # ramps down in 10 seconds
[end_time, 0.002],
]
) # end of third "flat" zone
mapdl.load_table("my_conv", my_conv, "TIME")
# Create a table of bulk temperatures for a given time and transfer to MAPDL
my_bulk = np.array(
[
[0, 100], # start time
[120, 100], # end of first "flat" zone
[500, 300], # ramps up in 380 seconds
[700, 300], # hold temperature for 200 seconds
[900, 75], # temperature ramps down for 200 seconds
[end_time, 75],
]
) # end of second "flat" zone
mapdl.load_table("my_bulk", my_bulk, "TIME")
The Transient Thermal Solve#
This model is to be solved in one time step. For this reason, a
TSRES
command is used for force the solver to include a
SOLVE
at every time point in the two Table Arrays above. This
ensures that the time-dependent curves are followed by the transient
analysis. Intermediate solutions between the TSRES
time points
will be included according to the DELTIM
command and the
automatic time stepping decisions of the ANSYS solver.
In this example, the times for the TSRES
array illustrated above
have been determined manually. A set of APDL commands could be used
to automate this process for chosen Table Array entries, in more
complex modeling situations, including checks that no time intervals
are too short.
Results at substeps will be wanted if the intermediate solutions of
the time-transient analysis are to be available for post-processing
review. The OUTRES
command is used to control how much is written to
the results file. In this example the OUTRES command will be used to
simply write out all results for all substeps. In work with large
models and may substeps, too much data will be written if such a
strategy is employed for OUTRES
, and other options will need to be
considered. Note that one option for the OUTRES
command is to
control times at which results are written with a Table Array, much
as is used in the TSRES
command, but typically for a larger number
of time points, although including those of the TSRES array.
The initial condition starting temperature is controlled for this
example with the TUNIF
command. Note that thermal transient
analyses can also have a starting temperature profile formed by a
static thermal SOLVE
. If a user neglects to set an initial
temperature in ANSYS Mechanical APDL, a value of zero will be used,
which is often not what is desired.
The thermal convective loads are applied with an SF family command—in this example a convective load is applied to the end face of the solid model by the SFA command, using the Table Array entries for convection and bulk temperature that were developed above. The Table Array names are surrounded with percent signs (%). A SOLVE is then performed.
# Force transient solve to include the times within the conv and bulk arrays
# my_tres = np.unique(np.vstack((my_bulk[:, 0], my_conv[:, 0])))[0] # same as
mapdl.parameters["my_tsres"] = [120, 130, 500, 700, 710, 900, end_time]
mapdl.tsres("%my_tsres%")
mapdl.outres("ERASE")
mapdl.outres("ALL", "ALL")
mapdl.eqslv("SPARSE") # use sparse solver
mapdl.tunif(75) # force uniform starting temperature (otherwise zero)
# apply the convective load (convection coefficient plus bulk temperature)
# use "%" around table array names
mapdl.sfa(6, 1, "CONV", "%my_conv%", " %my_bulk%")
# solve
mapdl.solve()
***** MAPDL SOLVE COMMAND *****
TRANSFER SOLID MODEL BOUNDARY CONDITIONS TO FINITE ELEMENT MODEL
SURFACE LOADS TRANSFERRED FROM AREAS = 25
*** NOTE *** CP = 0.000 TIME= 00:00:00
There is no title defined for this analysis.
*****MAPDL VERIFICATION RUN ONLY*****
DO NOT USE RESULTS FOR PRODUCTION
S O L U T I O N O P T I O N S
PROBLEM DIMENSIONALITY. . . . . . . . . . . . .3-D
DEGREES OF FREEDOM. . . . . . TEMP
ANALYSIS TYPE . . . . . . . . . . . . . . . . .TRANSIENT
SOLUTION METHOD. . . . . . . . . . . . . . .FULL
EQUATION SOLVER OPTION. . . . . . . . . . . . .SPARSE
GLOBALLY ASSEMBLED MATRIX . . . . . . . . . . .SYMMETRIC
*** NOTE *** CP = 0.000 TIME= 00:00:00
The conditions for direct assembly have been met. No .emat or .erot
files will be produced.
D I S T R I B U T E D D O M A I N D E C O M P O S E R
...Number of elements: 625
...Number of nodes: 936
...Decompose to 0 CPU domains
...Element load balance ratio = 0.000
L O A D S T E P O P T I O N S
LOAD STEP NUMBER. . . . . . . . . . . . . . . . 1
TIME AT END OF THE LOAD STEP. . . . . . . . . . 1500.0
AUTOMATIC TIME STEPPING . . . . . . . . . . . . ON
STARTING TIME STEP SIZE. . . . . . . . . . . 10.000
MINIMUM TIME STEP SIZE . . . . . . . . . . . 2.0000
MAXIMUM TIME STEP SIZE . . . . . . . . . . . 25.000
TIME STEP RESET ARRAY . . . . . . . . . . . . . MY_TSRES
STEP CHANGE BOUNDARY CONDITIONS . . . . . . . . NO
TRANSIENT (INERTIA) EFFECTS
THERMAL DOFS . . . . . . . . . . . . . . . . ON
TRANSIENT INTEGRATION PARAMETERS
THETA. . . . . . . . . . . . . . . . . . . . 1.0000
OSCILLATION LIMIT CRITERION. . . . . . . . . 0.50000
TOLERANCE. . . . . . . . . . . . . . . . . . 0.0000
PRINT OUTPUT CONTROLS . . . . . . . . . . . . .NO PRINTOUT
DATABASE OUTPUT CONTROLS
ITEM FREQUENCY COMPONENT
ALL ALL
**** CENTER OF MASS, MASS, AND MASS MOMENTS OF INERTIA ****
CALCULATIONS ASSUME ELEMENT MASS AT ELEMENT CENTROID
TOTAL MASS = 0.36732E-02
MOM. OF INERTIA MOM. OF INERTIA
CENTER OF MASS ABOUT ORIGIN ABOUT CENTER OF MASS
XC = 2.5000 IXX = 0.2424E-02 IXX = 0.5877E-03
YC = 0.50000 IYY = 0.3181E-01 IYY = 0.7934E-02
ZC = 0.50000 IZZ = 0.3181E-01 IZZ = 0.7934E-02
IXY = -0.4591E-02 IXY = 0.8674E-17
IYZ = -0.9183E-03 IYZ = 0.2819E-17
IZX = -0.4591E-02 IZX = 0.1475E-16
*** MASS SUMMARY BY ELEMENT TYPE ***
TYPE MASS
1 0.367317E-02
Range of element maximum matrix coefficients in global coordinates
Maximum = 5.168130667E-05 at element 0.
Minimum = 4.168130667E-05 at element 0.
*** ELEMENT MATRIX FORMULATION TIMES
TYPE NUMBER ENAME TOTAL CP AVE CP
1 625 SOLID70 0.000 0.000000
Time at end of element matrix formulation CP = 0.
ALL CURRENT MAPDL DATA WRITTEN TO FILE NAME=
FOR POSSIBLE RESUME FROM THIS POINT
DISTRIBUTED SPARSE MATRIX DIRECT SOLVER.
Number of equations = 936, Maximum wavefront = 0
Memory available (MB) = 0.0 , Memory required (MB) = 0.0
Distributed sparse solver maximum pivot= 0 at node 0 .
Distributed sparse solver minimum pivot= 0 at node 0 .
Distributed sparse solver minimum pivot in absolute value= 0 at node 0
.
*** ELEMENT RESULT CALCULATION TIMES
TYPE NUMBER ENAME TOTAL CP AVE CP
1 625 SOLID70 0.000 0.000000
*** NODAL LOAD CALCULATION TIMES
TYPE NUMBER ENAME TOTAL CP AVE CP
1 625 SOLID70 0.000 0.000000
*** LOAD STEP 1 SUBSTEP 1 COMPLETED. CUM ITER = 1
*** TIME = 10.0000 TIME INC = 10.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2324 OSCILLATION LIMIT = 2.324
*** AUTO TIME STEP: NEXT TIME INC = 3.3333 DECREASED (FACTOR = 0.3333)
*** LOAD STEP 1 SUBSTEP 2 COMPLETED. CUM ITER = 2
*** TIME = 13.3333 TIME INC = 3.33333 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.8326E-01 OSCILLATION LIMIT = 0.2775
*** AUTO STEP TIME: NEXT TIME INC = 3.3333 UNCHANGED
*** LOAD STEP 1 SUBSTEP 3 COMPLETED. CUM ITER = 3
*** TIME = 16.6667 TIME INC = 3.33333 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.5387E-01 OSCILLATION LIMIT = 0.1796
*** AUTO TIME STEP: NEXT TIME INC = 9.2821 INCREASED (FACTOR = 2.7846)
*** LOAD STEP 1 SUBSTEP 4 COMPLETED. CUM ITER = 4
*** TIME = 25.9488 TIME INC = 9.28214 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.3279E-01 OSCILLATION LIMIT = 0.3044
*** AUTO TIME STEP: NEXT TIME INC = 15.248 INCREASED (FACTOR = 1.6427)
*** LOAD STEP 1 SUBSTEP 5 COMPLETED. CUM ITER = 5
*** TIME = 41.1965 TIME INC = 15.2477 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2082E-01 OSCILLATION LIMIT = 0.3174
*** AUTO TIME STEP: NEXT TIME INC = 24.019 INCREASED (FACTOR = 1.5753)
*** LOAD STEP 1 SUBSTEP 6 COMPLETED. CUM ITER = 6
*** TIME = 65.2156 TIME INC = 24.0191 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1335E-01 OSCILLATION LIMIT = 0.3207
*** AUTO TIME STEP: NEXT TIME INC = 25.000 INCREASED (FACTOR = 1.0408)
*** LOAD STEP 1 SUBSTEP 7 COMPLETED. CUM ITER = 7
*** TIME = 90.2156 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.9188E-02 OSCILLATION LIMIT = 0.2297
*** AUTO TIME STEP: NEXT TIME INC = 14.892 DECREASED (FACTOR = 0.5957)
*** LOAD STEP 1 SUBSTEP 8 COMPLETED. CUM ITER = 8
*** TIME = 105.108 TIME INC = 14.8922 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.7332E-02 OSCILLATION LIMIT = 0.1092
*** AUTO STEP TIME: NEXT TIME INC = 14.892 UNCHANGED
*** LOAD STEP 1 SUBSTEP 9 COMPLETED. CUM ITER = 9
*** TIME = 120.000 TIME INC = 14.8922 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.5920E-02 OSCILLATION LIMIT = 0.8817E-01
*** AUTO TIME STEP: NEXT TIME INC = 10.000 DECREASED (FACTOR = 0.6715)
*** LOAD STEP 1 SUBSTEP 10 COMPLETED. CUM ITER = 10
*** TIME = 130.000 TIME INC = 10.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.6812 OSCILLATION LIMIT = 6.812
*** AUTO STEP TIME: NEXT TIME INC = 10.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 11 COMPLETED. CUM ITER = 11
*** TIME = 140.000 TIME INC = 10.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.3243 OSCILLATION LIMIT = 3.243
*** AUTO TIME STEP: NEXT TIME INC = 3.3333 DECREASED (FACTOR = 0.3333)
*** LOAD STEP 1 SUBSTEP 12 COMPLETED. CUM ITER = 12
*** TIME = 143.333 TIME INC = 3.33333 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2992 OSCILLATION LIMIT = 0.9973
*** AUTO STEP TIME: NEXT TIME INC = 3.3333 UNCHANGED
*** LOAD STEP 1 SUBSTEP 13 COMPLETED. CUM ITER = 13
*** TIME = 146.667 TIME INC = 3.33333 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2844 OSCILLATION LIMIT = 0.9480
*** AUTO STEP TIME: NEXT TIME INC = 3.3333 UNCHANGED
*** LOAD STEP 1 SUBSTEP 14 COMPLETED. CUM ITER = 14
*** TIME = 150.000 TIME INC = 3.33333 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2723 OSCILLATION LIMIT = 0.9075
*** AUTO TIME STEP: NEXT TIME INC = 2.0000 DECREASED (FACTOR = 0.6000)
*** LOAD STEP 1 SUBSTEP 15 COMPLETED. CUM ITER = 15
*** TIME = 152.000 TIME INC = 2.00000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2656 OSCILLATION LIMIT = 0.5313
*** AUTO STEP TIME: NEXT TIME INC = 2.0000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 16 COMPLETED. CUM ITER = 16
*** TIME = 154.000 TIME INC = 2.00000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2595 OSCILLATION LIMIT = 0.5189
*** AUTO STEP TIME: NEXT TIME INC = 2.0000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 17 COMPLETED. CUM ITER = 17
*** TIME = 156.000 TIME INC = 2.00000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2537 OSCILLATION LIMIT = 0.5074
*** AUTO STEP TIME: NEXT TIME INC = 2.0000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 18 COMPLETED. CUM ITER = 18
*** TIME = 158.000 TIME INC = 2.00000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2483 OSCILLATION LIMIT = 0.4967
*** AUTO TIME STEP: NEXT TIME INC = 2.0134 INCREASED (FACTOR = 1.0067)
*** LOAD STEP 1 SUBSTEP 19 COMPLETED. CUM ITER = 19
*** TIME = 160.013 TIME INC = 2.01340 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2432 OSCILLATION LIMIT = 0.4897
*** AUTO TIME STEP: NEXT TIME INC = 2.0556 INCREASED (FACTOR = 1.0209)
*** LOAD STEP 1 SUBSTEP 20 COMPLETED. CUM ITER = 20
*** TIME = 162.069 TIME INC = 2.05556 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2384 OSCILLATION LIMIT = 0.4899
*** AUTO TIME STEP: NEXT TIME INC = 2.0977 INCREASED (FACTOR = 1.0205)
*** LOAD STEP 1 SUBSTEP 21 COMPLETED. CUM ITER = 21
*** TIME = 164.167 TIME INC = 2.09773 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2337 OSCILLATION LIMIT = 0.4901
*** AUTO TIME STEP: NEXT TIME INC = 2.1399 INCREASED (FACTOR = 1.0201)
*** LOAD STEP 1 SUBSTEP 22 COMPLETED. CUM ITER = 22
*** TIME = 166.307 TIME INC = 2.13991 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2291 OSCILLATION LIMIT = 0.4903
*** AUTO TIME STEP: NEXT TIME INC = 2.1821 INCREASED (FACTOR = 1.0197)
*** LOAD STEP 1 SUBSTEP 23 COMPLETED. CUM ITER = 23
*** TIME = 168.489 TIME INC = 2.18210 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2248 OSCILLATION LIMIT = 0.4905
*** AUTO TIME STEP: NEXT TIME INC = 2.2243 INCREASED (FACTOR = 1.0193)
*** LOAD STEP 1 SUBSTEP 24 COMPLETED. CUM ITER = 24
*** TIME = 170.713 TIME INC = 2.22429 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2206 OSCILLATION LIMIT = 0.4907
*** AUTO TIME STEP: NEXT TIME INC = 2.2665 INCREASED (FACTOR = 1.0190)
*** LOAD STEP 1 SUBSTEP 25 COMPLETED. CUM ITER = 25
*** TIME = 172.979 TIME INC = 2.26648 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2166 OSCILLATION LIMIT = 0.4909
*** AUTO TIME STEP: NEXT TIME INC = 2.3087 INCREASED (FACTOR = 1.0186)
*** LOAD STEP 1 SUBSTEP 26 COMPLETED. CUM ITER = 26
*** TIME = 175.288 TIME INC = 2.30867 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2127 OSCILLATION LIMIT = 0.4910
*** AUTO TIME STEP: NEXT TIME INC = 2.3509 INCREASED (FACTOR = 1.0183)
*** LOAD STEP 1 SUBSTEP 27 COMPLETED. CUM ITER = 27
*** TIME = 177.639 TIME INC = 2.35087 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2089 OSCILLATION LIMIT = 0.4912
*** AUTO TIME STEP: NEXT TIME INC = 2.3931 INCREASED (FACTOR = 1.0180)
*** LOAD STEP 1 SUBSTEP 28 COMPLETED. CUM ITER = 28
*** TIME = 180.032 TIME INC = 2.39307 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2053 OSCILLATION LIMIT = 0.4913
*** AUTO TIME STEP: NEXT TIME INC = 2.4353 INCREASED (FACTOR = 1.0176)
*** LOAD STEP 1 SUBSTEP 29 COMPLETED. CUM ITER = 29
*** TIME = 182.467 TIME INC = 2.43529 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2018 OSCILLATION LIMIT = 0.4915
*** AUTO TIME STEP: NEXT TIME INC = 2.4775 INCREASED (FACTOR = 1.0173)
*** LOAD STEP 1 SUBSTEP 30 COMPLETED. CUM ITER = 30
*** TIME = 184.945 TIME INC = 2.47752 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1984 OSCILLATION LIMIT = 0.4916
*** AUTO TIME STEP: NEXT TIME INC = 2.5198 INCREASED (FACTOR = 1.0171)
*** LOAD STEP 1 SUBSTEP 31 COMPLETED. CUM ITER = 31
*** TIME = 187.465 TIME INC = 2.51977 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1952 OSCILLATION LIMIT = 0.4917
*** AUTO TIME STEP: NEXT TIME INC = 2.5621 INCREASED (FACTOR = 1.0168)
*** LOAD STEP 1 SUBSTEP 32 COMPLETED. CUM ITER = 32
*** TIME = 190.027 TIME INC = 2.56205 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1920 OSCILLATION LIMIT = 0.4919
*** AUTO TIME STEP: NEXT TIME INC = 2.6044 INCREASED (FACTOR = 1.0165)
*** LOAD STEP 1 SUBSTEP 33 COMPLETED. CUM ITER = 33
*** TIME = 192.631 TIME INC = 2.60438 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1889 OSCILLATION LIMIT = 0.4920
*** AUTO TIME STEP: NEXT TIME INC = 2.6467 INCREASED (FACTOR = 1.0163)
*** LOAD STEP 1 SUBSTEP 34 COMPLETED. CUM ITER = 34
*** TIME = 195.278 TIME INC = 2.64674 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1859 OSCILLATION LIMIT = 0.4921
*** AUTO TIME STEP: NEXT TIME INC = 2.6892 INCREASED (FACTOR = 1.0160)
*** LOAD STEP 1 SUBSTEP 35 COMPLETED. CUM ITER = 35
*** TIME = 197.967 TIME INC = 2.68917 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1830 OSCILLATION LIMIT = 0.4922
*** AUTO TIME STEP: NEXT TIME INC = 2.7317 INCREASED (FACTOR = 1.0158)
*** LOAD STEP 1 SUBSTEP 36 COMPLETED. CUM ITER = 36
*** TIME = 200.699 TIME INC = 2.73166 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1802 OSCILLATION LIMIT = 0.4923
*** AUTO TIME STEP: NEXT TIME INC = 2.7742 INCREASED (FACTOR = 1.0156)
*** LOAD STEP 1 SUBSTEP 37 COMPLETED. CUM ITER = 37
*** TIME = 203.473 TIME INC = 2.77422 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1775 OSCILLATION LIMIT = 0.4924
*** AUTO TIME STEP: NEXT TIME INC = 2.8169 INCREASED (FACTOR = 1.0154)
*** LOAD STEP 1 SUBSTEP 38 COMPLETED. CUM ITER = 38
*** TIME = 206.290 TIME INC = 2.81688 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1748 OSCILLATION LIMIT = 0.4925
*** AUTO TIME STEP: NEXT TIME INC = 2.8596 INCREASED (FACTOR = 1.0152)
*** LOAD STEP 1 SUBSTEP 39 COMPLETED. CUM ITER = 39
*** TIME = 209.149 TIME INC = 2.85963 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1723 OSCILLATION LIMIT = 0.4926
*** AUTO TIME STEP: NEXT TIME INC = 2.9025 INCREASED (FACTOR = 1.0150)
*** LOAD STEP 1 SUBSTEP 40 COMPLETED. CUM ITER = 40
*** TIME = 212.052 TIME INC = 2.90250 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1698 OSCILLATION LIMIT = 0.4927
*** AUTO TIME STEP: NEXT TIME INC = 2.9455 INCREASED (FACTOR = 1.0148)
*** LOAD STEP 1 SUBSTEP 41 COMPLETED. CUM ITER = 41
*** TIME = 214.997 TIME INC = 2.94549 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1673 OSCILLATION LIMIT = 0.4928
*** AUTO TIME STEP: NEXT TIME INC = 2.9886 INCREASED (FACTOR = 1.0146)
*** LOAD STEP 1 SUBSTEP 42 COMPLETED. CUM ITER = 42
*** TIME = 217.986 TIME INC = 2.98863 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1649 OSCILLATION LIMIT = 0.4929
*** AUTO TIME STEP: NEXT TIME INC = 3.0319 INCREASED (FACTOR = 1.0145)
*** LOAD STEP 1 SUBSTEP 43 COMPLETED. CUM ITER = 43
*** TIME = 221.018 TIME INC = 3.03192 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1626 OSCILLATION LIMIT = 0.4929
*** AUTO TIME STEP: NEXT TIME INC = 3.0754 INCREASED (FACTOR = 1.0143)
*** LOAD STEP 1 SUBSTEP 44 COMPLETED. CUM ITER = 44
*** TIME = 224.093 TIME INC = 3.07537 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1603 OSCILLATION LIMIT = 0.4930
*** AUTO TIME STEP: NEXT TIME INC = 3.1190 INCREASED (FACTOR = 1.0142)
*** LOAD STEP 1 SUBSTEP 45 COMPLETED. CUM ITER = 45
*** TIME = 227.212 TIME INC = 3.11902 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1581 OSCILLATION LIMIT = 0.4931
*** AUTO TIME STEP: NEXT TIME INC = 3.1629 INCREASED (FACTOR = 1.0141)
*** LOAD STEP 1 SUBSTEP 46 COMPLETED. CUM ITER = 46
*** TIME = 230.375 TIME INC = 3.16286 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1559 OSCILLATION LIMIT = 0.4931
*** AUTO TIME STEP: NEXT TIME INC = 3.2069 INCREASED (FACTOR = 1.0139)
*** LOAD STEP 1 SUBSTEP 47 COMPLETED. CUM ITER = 47
*** TIME = 233.582 TIME INC = 3.20692 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1538 OSCILLATION LIMIT = 0.4932
*** AUTO TIME STEP: NEXT TIME INC = 3.2512 INCREASED (FACTOR = 1.0138)
*** LOAD STEP 1 SUBSTEP 48 COMPLETED. CUM ITER = 48
*** TIME = 236.833 TIME INC = 3.25122 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1517 OSCILLATION LIMIT = 0.4932
*** AUTO TIME STEP: NEXT TIME INC = 3.2958 INCREASED (FACTOR = 1.0137)
*** LOAD STEP 1 SUBSTEP 49 COMPLETED. CUM ITER = 49
*** TIME = 240.129 TIME INC = 3.29577 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1497 OSCILLATION LIMIT = 0.4933
*** AUTO TIME STEP: NEXT TIME INC = 3.3406 INCREASED (FACTOR = 1.0136)
*** LOAD STEP 1 SUBSTEP 50 COMPLETED. CUM ITER = 50
*** TIME = 243.470 TIME INC = 3.34059 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1477 OSCILLATION LIMIT = 0.4933
*** AUTO TIME STEP: NEXT TIME INC = 3.3857 INCREASED (FACTOR = 1.0135)
*** LOAD STEP 1 SUBSTEP 51 COMPLETED. CUM ITER = 51
*** TIME = 246.855 TIME INC = 3.38569 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1457 OSCILLATION LIMIT = 0.4934
*** AUTO TIME STEP: NEXT TIME INC = 3.4311 INCREASED (FACTOR = 1.0134)
*** LOAD STEP 1 SUBSTEP 52 COMPLETED. CUM ITER = 52
*** TIME = 250.286 TIME INC = 3.43110 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1438 OSCILLATION LIMIT = 0.4934
*** AUTO TIME STEP: NEXT TIME INC = 3.4768 INCREASED (FACTOR = 1.0133)
*** LOAD STEP 1 SUBSTEP 53 COMPLETED. CUM ITER = 53
*** TIME = 253.763 TIME INC = 3.47684 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1419 OSCILLATION LIMIT = 0.4935
*** AUTO TIME STEP: NEXT TIME INC = 3.5229 INCREASED (FACTOR = 1.0133)
*** LOAD STEP 1 SUBSTEP 54 COMPLETED. CUM ITER = 54
*** TIME = 257.286 TIME INC = 3.52291 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1401 OSCILLATION LIMIT = 0.4935
*** AUTO TIME STEP: NEXT TIME INC = 3.5694 INCREASED (FACTOR = 1.0132)
*** LOAD STEP 1 SUBSTEP 55 COMPLETED. CUM ITER = 55
*** TIME = 260.856 TIME INC = 3.56935 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1383 OSCILLATION LIMIT = 0.4935
*** AUTO TIME STEP: NEXT TIME INC = 3.6162 INCREASED (FACTOR = 1.0131)
*** LOAD STEP 1 SUBSTEP 56 COMPLETED. CUM ITER = 56
*** TIME = 264.472 TIME INC = 3.61617 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1365 OSCILLATION LIMIT = 0.4936
*** AUTO TIME STEP: NEXT TIME INC = 3.6634 INCREASED (FACTOR = 1.0131)
*** LOAD STEP 1 SUBSTEP 57 COMPLETED. CUM ITER = 57
*** TIME = 268.135 TIME INC = 3.66338 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1347 OSCILLATION LIMIT = 0.4936
*** AUTO TIME STEP: NEXT TIME INC = 3.7110 INCREASED (FACTOR = 1.0130)
*** LOAD STEP 1 SUBSTEP 58 COMPLETED. CUM ITER = 58
*** TIME = 271.846 TIME INC = 3.71101 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1330 OSCILLATION LIMIT = 0.4936
*** AUTO TIME STEP: NEXT TIME INC = 3.7591 INCREASED (FACTOR = 1.0130)
*** LOAD STEP 1 SUBSTEP 59 COMPLETED. CUM ITER = 59
*** TIME = 275.605 TIME INC = 3.75907 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1313 OSCILLATION LIMIT = 0.4936
*** AUTO TIME STEP: NEXT TIME INC = 3.8076 INCREASED (FACTOR = 1.0129)
*** LOAD STEP 1 SUBSTEP 60 COMPLETED. CUM ITER = 60
*** TIME = 279.413 TIME INC = 3.80759 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1296 OSCILLATION LIMIT = 0.4936
*** AUTO TIME STEP: NEXT TIME INC = 3.8566 INCREASED (FACTOR = 1.0129)
*** LOAD STEP 1 SUBSTEP 61 COMPLETED. CUM ITER = 61
*** TIME = 283.269 TIME INC = 3.85658 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1280 OSCILLATION LIMIT = 0.4937
*** AUTO TIME STEP: NEXT TIME INC = 3.9061 INCREASED (FACTOR = 1.0128)
*** LOAD STEP 1 SUBSTEP 62 COMPLETED. CUM ITER = 62
*** TIME = 287.175 TIME INC = 3.90605 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1264 OSCILLATION LIMIT = 0.4937
*** AUTO TIME STEP: NEXT TIME INC = 3.9560 INCREASED (FACTOR = 1.0128)
*** LOAD STEP 1 SUBSTEP 63 COMPLETED. CUM ITER = 63
*** TIME = 291.131 TIME INC = 3.95603 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1248 OSCILLATION LIMIT = 0.4937
*** AUTO TIME STEP: NEXT TIME INC = 4.0065 INCREASED (FACTOR = 1.0128)
*** LOAD STEP 1 SUBSTEP 64 COMPLETED. CUM ITER = 64
*** TIME = 295.138 TIME INC = 4.00654 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1232 OSCILLATION LIMIT = 0.4937
*** AUTO TIME STEP: NEXT TIME INC = 4.0576 INCREASED (FACTOR = 1.0127)
*** LOAD STEP 1 SUBSTEP 65 COMPLETED. CUM ITER = 65
*** TIME = 299.196 TIME INC = 4.05758 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1217 OSCILLATION LIMIT = 0.4937
*** AUTO TIME STEP: NEXT TIME INC = 4.1092 INCREASED (FACTOR = 1.0127)
*** LOAD STEP 1 SUBSTEP 66 COMPLETED. CUM ITER = 66
*** TIME = 303.305 TIME INC = 4.10918 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1202 OSCILLATION LIMIT = 0.4937
*** AUTO TIME STEP: NEXT TIME INC = 4.1613 INCREASED (FACTOR = 1.0127)
*** LOAD STEP 1 SUBSTEP 67 COMPLETED. CUM ITER = 67
*** TIME = 307.466 TIME INC = 4.16135 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1186 OSCILLATION LIMIT = 0.4937
*** AUTO TIME STEP: NEXT TIME INC = 4.2141 INCREASED (FACTOR = 1.0127)
*** LOAD STEP 1 SUBSTEP 68 COMPLETED. CUM ITER = 68
*** TIME = 311.680 TIME INC = 4.21410 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1172 OSCILLATION LIMIT = 0.4937
*** AUTO TIME STEP: NEXT TIME INC = 4.2675 INCREASED (FACTOR = 1.0127)
*** LOAD STEP 1 SUBSTEP 69 COMPLETED. CUM ITER = 69
*** TIME = 315.948 TIME INC = 4.26746 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1157 OSCILLATION LIMIT = 0.4938
*** AUTO TIME STEP: NEXT TIME INC = 4.3214 INCREASED (FACTOR = 1.0126)
*** LOAD STEP 1 SUBSTEP 70 COMPLETED. CUM ITER = 70
*** TIME = 320.269 TIME INC = 4.32143 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1143 OSCILLATION LIMIT = 0.4938
*** AUTO TIME STEP: NEXT TIME INC = 4.3760 INCREASED (FACTOR = 1.0126)
*** LOAD STEP 1 SUBSTEP 71 COMPLETED. CUM ITER = 71
*** TIME = 324.645 TIME INC = 4.37603 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1128 OSCILLATION LIMIT = 0.4938
*** AUTO TIME STEP: NEXT TIME INC = 4.4313 INCREASED (FACTOR = 1.0126)
*** LOAD STEP 1 SUBSTEP 72 COMPLETED. CUM ITER = 72
*** TIME = 329.076 TIME INC = 4.43128 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1114 OSCILLATION LIMIT = 0.4938
*** AUTO TIME STEP: NEXT TIME INC = 4.4872 INCREASED (FACTOR = 1.0126)
*** LOAD STEP 1 SUBSTEP 73 COMPLETED. CUM ITER = 73
*** TIME = 333.564 TIME INC = 4.48717 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1100 OSCILLATION LIMIT = 0.4938
*** AUTO TIME STEP: NEXT TIME INC = 4.5437 INCREASED (FACTOR = 1.0126)
*** LOAD STEP 1 SUBSTEP 74 COMPLETED. CUM ITER = 74
*** TIME = 338.107 TIME INC = 4.54374 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1087 OSCILLATION LIMIT = 0.4938
*** AUTO TIME STEP: NEXT TIME INC = 4.6010 INCREASED (FACTOR = 1.0126)
*** LOAD STEP 1 SUBSTEP 75 COMPLETED. CUM ITER = 75
*** TIME = 342.708 TIME INC = 4.60097 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1073 OSCILLATION LIMIT = 0.4938
*** AUTO TIME STEP: NEXT TIME INC = 4.6589 INCREASED (FACTOR = 1.0126)
*** LOAD STEP 1 SUBSTEP 76 COMPLETED. CUM ITER = 76
*** TIME = 347.367 TIME INC = 4.65890 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1060 OSCILLATION LIMIT = 0.4938
*** AUTO TIME STEP: NEXT TIME INC = 4.7175 INCREASED (FACTOR = 1.0126)
*** LOAD STEP 1 SUBSTEP 77 COMPLETED. CUM ITER = 77
*** TIME = 352.085 TIME INC = 4.71751 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1047 OSCILLATION LIMIT = 0.4938
*** AUTO TIME STEP: NEXT TIME INC = 4.7768 INCREASED (FACTOR = 1.0126)
*** LOAD STEP 1 SUBSTEP 78 COMPLETED. CUM ITER = 78
*** TIME = 356.862 TIME INC = 4.77684 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1034 OSCILLATION LIMIT = 0.4938
*** AUTO TIME STEP: NEXT TIME INC = 4.8369 INCREASED (FACTOR = 1.0126)
*** LOAD STEP 1 SUBSTEP 79 COMPLETED. CUM ITER = 79
*** TIME = 361.698 TIME INC = 4.83687 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1021 OSCILLATION LIMIT = 0.4938
*** AUTO TIME STEP: NEXT TIME INC = 4.8976 INCREASED (FACTOR = 1.0126)
*** LOAD STEP 1 SUBSTEP 80 COMPLETED. CUM ITER = 80
*** TIME = 366.596 TIME INC = 4.89761 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1008 OSCILLATION LIMIT = 0.4938
*** AUTO TIME STEP: NEXT TIME INC = 4.9591 INCREASED (FACTOR = 1.0126)
*** LOAD STEP 1 SUBSTEP 81 COMPLETED. CUM ITER = 81
*** TIME = 371.555 TIME INC = 4.95908 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.9958E-01 OSCILLATION LIMIT = 0.4938
*** AUTO TIME STEP: NEXT TIME INC = 5.0213 INCREASED (FACTOR = 1.0125)
*** LOAD STEP 1 SUBSTEP 82 COMPLETED. CUM ITER = 82
*** TIME = 376.576 TIME INC = 5.02128 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.9834E-01 OSCILLATION LIMIT = 0.4938
*** AUTO TIME STEP: NEXT TIME INC = 5.0842 INCREASED (FACTOR = 1.0125)
*** LOAD STEP 1 SUBSTEP 83 COMPLETED. CUM ITER = 83
*** TIME = 381.661 TIME INC = 5.08421 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.9713E-01 OSCILLATION LIMIT = 0.4938
*** AUTO TIME STEP: NEXT TIME INC = 5.1479 INCREASED (FACTOR = 1.0125)
*** LOAD STEP 1 SUBSTEP 84 COMPLETED. CUM ITER = 84
*** TIME = 386.808 TIME INC = 5.14787 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.9593E-01 OSCILLATION LIMIT = 0.4938
*** AUTO TIME STEP: NEXT TIME INC = 5.2123 INCREASED (FACTOR = 1.0125)
*** LOAD STEP 1 SUBSTEP 85 COMPLETED. CUM ITER = 85
*** TIME = 392.021 TIME INC = 5.21227 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.9474E-01 OSCILLATION LIMIT = 0.4938
*** AUTO TIME STEP: NEXT TIME INC = 5.2774 INCREASED (FACTOR = 1.0125)
*** LOAD STEP 1 SUBSTEP 86 COMPLETED. CUM ITER = 86
*** TIME = 397.298 TIME INC = 5.27740 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.9358E-01 OSCILLATION LIMIT = 0.4938
*** AUTO TIME STEP: NEXT TIME INC = 5.3433 INCREASED (FACTOR = 1.0125)
*** LOAD STEP 1 SUBSTEP 87 COMPLETED. CUM ITER = 87
*** TIME = 402.641 TIME INC = 5.34327 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.9242E-01 OSCILLATION LIMIT = 0.4938
*** AUTO TIME STEP: NEXT TIME INC = 5.4099 INCREASED (FACTOR = 1.0125)
*** LOAD STEP 1 SUBSTEP 88 COMPLETED. CUM ITER = 88
*** TIME = 408.051 TIME INC = 5.40988 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.9129E-01 OSCILLATION LIMIT = 0.4939
*** AUTO TIME STEP: NEXT TIME INC = 5.4772 INCREASED (FACTOR = 1.0124)
*** LOAD STEP 1 SUBSTEP 89 COMPLETED. CUM ITER = 89
*** TIME = 413.528 TIME INC = 5.47722 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.9017E-01 OSCILLATION LIMIT = 0.4939
*** AUTO TIME STEP: NEXT TIME INC = 5.5453 INCREASED (FACTOR = 1.0124)
*** LOAD STEP 1 SUBSTEP 90 COMPLETED. CUM ITER = 90
*** TIME = 419.074 TIME INC = 5.54529 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.8906E-01 OSCILLATION LIMIT = 0.4939
*** AUTO TIME STEP: NEXT TIME INC = 5.6141 INCREASED (FACTOR = 1.0124)
*** LOAD STEP 1 SUBSTEP 91 COMPLETED. CUM ITER = 91
*** TIME = 424.688 TIME INC = 5.61408 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.8797E-01 OSCILLATION LIMIT = 0.4939
*** AUTO TIME STEP: NEXT TIME INC = 5.6836 INCREASED (FACTOR = 1.0124)
*** LOAD STEP 1 SUBSTEP 92 COMPLETED. CUM ITER = 92
*** TIME = 430.371 TIME INC = 5.68359 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.8690E-01 OSCILLATION LIMIT = 0.4939
*** AUTO TIME STEP: NEXT TIME INC = 5.7538 INCREASED (FACTOR = 1.0124)
*** LOAD STEP 1 SUBSTEP 93 COMPLETED. CUM ITER = 93
*** TIME = 436.125 TIME INC = 5.75382 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.8584E-01 OSCILLATION LIMIT = 0.4939
*** AUTO TIME STEP: NEXT TIME INC = 5.8248 INCREASED (FACTOR = 1.0123)
*** LOAD STEP 1 SUBSTEP 94 COMPLETED. CUM ITER = 94
*** TIME = 441.950 TIME INC = 5.82476 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.8480E-01 OSCILLATION LIMIT = 0.4939
*** AUTO TIME STEP: NEXT TIME INC = 5.8964 INCREASED (FACTOR = 1.0123)
*** LOAD STEP 1 SUBSTEP 95 COMPLETED. CUM ITER = 95
*** TIME = 447.846 TIME INC = 5.89639 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.8377E-01 OSCILLATION LIMIT = 0.4939
*** AUTO TIME STEP: NEXT TIME INC = 5.9687 INCREASED (FACTOR = 1.0123)
*** LOAD STEP 1 SUBSTEP 96 COMPLETED. CUM ITER = 96
*** TIME = 453.815 TIME INC = 5.96871 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.8276E-01 OSCILLATION LIMIT = 0.4940
*** AUTO TIME STEP: NEXT TIME INC = 6.0417 INCREASED (FACTOR = 1.0122)
*** LOAD STEP 1 SUBSTEP 97 COMPLETED. CUM ITER = 97
*** TIME = 459.857 TIME INC = 6.04171 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.8176E-01 OSCILLATION LIMIT = 0.4940
*** AUTO TIME STEP: NEXT TIME INC = 6.1154 INCREASED (FACTOR = 1.0122)
*** LOAD STEP 1 SUBSTEP 98 COMPLETED. CUM ITER = 98
*** TIME = 465.972 TIME INC = 6.11538 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.8078E-01 OSCILLATION LIMIT = 0.4940
*** AUTO TIME STEP: NEXT TIME INC = 6.1897 INCREASED (FACTOR = 1.0122)
*** LOAD STEP 1 SUBSTEP 99 COMPLETED. CUM ITER = 99
*** TIME = 472.162 TIME INC = 6.18970 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.7981E-01 OSCILLATION LIMIT = 0.4940
*** AUTO TIME STEP: NEXT TIME INC = 6.2647 INCREASED (FACTOR = 1.0121)
*** LOAD STEP 1 SUBSTEP 100 COMPLETED. CUM ITER = 100
*** TIME = 478.427 TIME INC = 6.26466 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.7886E-01 OSCILLATION LIMIT = 0.4940
*** AUTO TIME STEP: NEXT TIME INC = 6.3402 INCREASED (FACTOR = 1.0121)
*** LOAD STEP 1 SUBSTEP 101 COMPLETED. CUM ITER = 101
*** TIME = 484.767 TIME INC = 6.34024 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.7792E-01 OSCILLATION LIMIT = 0.4941
*** AUTO TIME STEP: NEXT TIME INC = 6.4164 INCREASED (FACTOR = 1.0120)
*** LOAD STEP 1 SUBSTEP 102 COMPLETED. CUM ITER = 102
*** TIME = 491.183 TIME INC = 6.41644 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.7700E-01 OSCILLATION LIMIT = 0.4941
*** AUTO TIME STEP: NEXT TIME INC = 6.4932 INCREASED (FACTOR = 1.0120)
*** LOAD STEP 1 SUBSTEP 103 COMPLETED. CUM ITER = 103
*** TIME = 497.676 TIME INC = 6.49323 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.7610E-01 OSCILLATION LIMIT = 0.4941
*** AUTO TIME STEP: NEXT TIME INC = 2.3235 DECREASED (FACTOR = 0.3578)
*** LOAD STEP 1 SUBSTEP 104 COMPLETED. CUM ITER = 104
*** TIME = 500.000 TIME INC = 2.32353 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.7578E-01 OSCILLATION LIMIT = 0.1761
*** AUTO TIME STEP: NEXT TIME INC = 10.000 INCREASED (FACTOR = 4.3038)
*** LOAD STEP 1 SUBSTEP 105 COMPLETED. CUM ITER = 105
*** TIME = 510.000 TIME INC = 10.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.9532E-02 OSCILLATION LIMIT = 0.9532E-01
*** AUTO TIME STEP: NEXT TIME INC = 25.000 INCREASED (FACTOR = 2.5000)
*** LOAD STEP 1 SUBSTEP 106 COMPLETED. CUM ITER = 106
*** TIME = 535.000 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.4495E-02 OSCILLATION LIMIT = 0.1124
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 107 COMPLETED. CUM ITER = 107
*** TIME = 560.000 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.3177E-02 OSCILLATION LIMIT = 0.7942E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 108 COMPLETED. CUM ITER = 108
*** TIME = 585.000 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2619E-02 OSCILLATION LIMIT = 0.6549E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 109 COMPLETED. CUM ITER = 109
*** TIME = 610.000 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2350E-02 OSCILLATION LIMIT = 0.5874E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 110 COMPLETED. CUM ITER = 110
*** TIME = 635.000 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2213E-02 OSCILLATION LIMIT = 0.5533E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 111 COMPLETED. CUM ITER = 111
*** TIME = 660.000 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2143E-02 OSCILLATION LIMIT = 0.5358E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 112 COMPLETED. CUM ITER = 112
*** TIME = 685.000 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2107E-02 OSCILLATION LIMIT = 0.5268E-01
*** AUTO TIME STEP: NEXT TIME INC = 15.000 DECREASED (FACTOR = 0.6000)
*** LOAD STEP 1 SUBSTEP 113 COMPLETED. CUM ITER = 113
*** TIME = 700.000 TIME INC = 15.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2094E-02 OSCILLATION LIMIT = 0.3141E-01
*** AUTO TIME STEP: NEXT TIME INC = 10.000 DECREASED (FACTOR = 0.6667)
*** LOAD STEP 1 SUBSTEP 114 COMPLETED. CUM ITER = 114
*** TIME = 710.000 TIME INC = 10.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2760 OSCILLATION LIMIT = 2.760
*** AUTO STEP TIME: NEXT TIME INC = 10.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 115 COMPLETED. CUM ITER = 115
*** TIME = 720.000 TIME INC = 10.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1647 OSCILLATION LIMIT = 1.647
*** AUTO TIME STEP: NEXT TIME INC = 3.3333 DECREASED (FACTOR = 0.3333)
*** LOAD STEP 1 SUBSTEP 116 COMPLETED. CUM ITER = 116
*** TIME = 723.333 TIME INC = 3.33333 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1550 OSCILLATION LIMIT = 0.5168
*** AUTO STEP TIME: NEXT TIME INC = 3.3333 UNCHANGED
*** LOAD STEP 1 SUBSTEP 117 COMPLETED. CUM ITER = 117
*** TIME = 726.667 TIME INC = 3.33333 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1495 OSCILLATION LIMIT = 0.4983
*** AUTO TIME STEP: NEXT TIME INC = 3.3446 INCREASED (FACTOR = 1.0034)
*** LOAD STEP 1 SUBSTEP 118 COMPLETED. CUM ITER = 118
*** TIME = 730.011 TIME INC = 3.34458 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1451 OSCILLATION LIMIT = 0.4854
*** AUTO TIME STEP: NEXT TIME INC = 3.4455 INCREASED (FACTOR = 1.0302)
*** LOAD STEP 1 SUBSTEP 119 COMPLETED. CUM ITER = 119
*** TIME = 733.457 TIME INC = 3.44549 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1412 OSCILLATION LIMIT = 0.4864
*** AUTO TIME STEP: NEXT TIME INC = 3.5415 INCREASED (FACTOR = 1.0279)
*** LOAD STEP 1 SUBSTEP 120 COMPLETED. CUM ITER = 120
*** TIME = 736.998 TIME INC = 3.54150 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1375 OSCILLATION LIMIT = 0.4871
*** AUTO TIME STEP: NEXT TIME INC = 3.6354 INCREASED (FACTOR = 1.0265)
*** LOAD STEP 1 SUBSTEP 121 COMPLETED. CUM ITER = 121
*** TIME = 740.634 TIME INC = 3.63542 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1341 OSCILLATION LIMIT = 0.4876
*** AUTO TIME STEP: NEXT TIME INC = 3.7282 INCREASED (FACTOR = 1.0255)
*** LOAD STEP 1 SUBSTEP 122 COMPLETED. CUM ITER = 122
*** TIME = 744.362 TIME INC = 3.72823 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1309 OSCILLATION LIMIT = 0.4879
*** AUTO TIME STEP: NEXT TIME INC = 3.8203 INCREASED (FACTOR = 1.0247)
*** LOAD STEP 1 SUBSTEP 123 COMPLETED. CUM ITER = 123
*** TIME = 748.182 TIME INC = 3.82033 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1278 OSCILLATION LIMIT = 0.4883
*** AUTO TIME STEP: NEXT TIME INC = 3.9119 INCREASED (FACTOR = 1.0240)
*** LOAD STEP 1 SUBSTEP 124 COMPLETED. CUM ITER = 124
*** TIME = 752.094 TIME INC = 3.91194 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1249 OSCILLATION LIMIT = 0.4886
*** AUTO TIME STEP: NEXT TIME INC = 4.0032 INCREASED (FACTOR = 1.0233)
*** LOAD STEP 1 SUBSTEP 125 COMPLETED. CUM ITER = 125
*** TIME = 756.097 TIME INC = 4.00319 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1221 OSCILLATION LIMIT = 0.4889
*** AUTO TIME STEP: NEXT TIME INC = 4.0942 INCREASED (FACTOR = 1.0227)
*** LOAD STEP 1 SUBSTEP 126 COMPLETED. CUM ITER = 126
*** TIME = 760.192 TIME INC = 4.09422 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1195 OSCILLATION LIMIT = 0.4891
*** AUTO TIME STEP: NEXT TIME INC = 4.1851 INCREASED (FACTOR = 1.0222)
*** LOAD STEP 1 SUBSTEP 127 COMPLETED. CUM ITER = 127
*** TIME = 764.377 TIME INC = 4.18515 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1169 OSCILLATION LIMIT = 0.4894
*** AUTO TIME STEP: NEXT TIME INC = 4.2761 INCREASED (FACTOR = 1.0217)
*** LOAD STEP 1 SUBSTEP 128 COMPLETED. CUM ITER = 128
*** TIME = 768.653 TIME INC = 4.27609 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1145 OSCILLATION LIMIT = 0.4896
*** AUTO TIME STEP: NEXT TIME INC = 4.3672 INCREASED (FACTOR = 1.0213)
*** LOAD STEP 1 SUBSTEP 129 COMPLETED. CUM ITER = 129
*** TIME = 773.020 TIME INC = 4.36718 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1121 OSCILLATION LIMIT = 0.4898
*** AUTO TIME STEP: NEXT TIME INC = 4.4586 INCREASED (FACTOR = 1.0209)
*** LOAD STEP 1 SUBSTEP 130 COMPLETED. CUM ITER = 130
*** TIME = 777.479 TIME INC = 4.45858 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1099 OSCILLATION LIMIT = 0.4899
*** AUTO TIME STEP: NEXT TIME INC = 4.5504 INCREASED (FACTOR = 1.0206)
*** LOAD STEP 1 SUBSTEP 131 COMPLETED. CUM ITER = 131
*** TIME = 782.029 TIME INC = 4.55044 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1077 OSCILLATION LIMIT = 0.4900
*** AUTO TIME STEP: NEXT TIME INC = 4.6429 INCREASED (FACTOR = 1.0203)
*** LOAD STEP 1 SUBSTEP 132 COMPLETED. CUM ITER = 132
*** TIME = 786.672 TIME INC = 4.64292 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1056 OSCILLATION LIMIT = 0.4902
*** AUTO TIME STEP: NEXT TIME INC = 4.7362 INCREASED (FACTOR = 1.0201)
*** LOAD STEP 1 SUBSTEP 133 COMPLETED. CUM ITER = 133
*** TIME = 791.408 TIME INC = 4.73621 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1035 OSCILLATION LIMIT = 0.4902
*** AUTO TIME STEP: NEXT TIME INC = 4.8305 INCREASED (FACTOR = 1.0199)
*** LOAD STEP 1 SUBSTEP 134 COMPLETED. CUM ITER = 134
*** TIME = 796.239 TIME INC = 4.83050 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1015 OSCILLATION LIMIT = 0.4903
*** AUTO TIME STEP: NEXT TIME INC = 4.9260 INCREASED (FACTOR = 1.0198)
*** LOAD STEP 1 SUBSTEP 135 COMPLETED. CUM ITER = 135
*** TIME = 801.165 TIME INC = 4.92601 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.9954E-01 OSCILLATION LIMIT = 0.4904
*** AUTO TIME STEP: NEXT TIME INC = 5.0230 INCREASED (FACTOR = 1.0197)
*** LOAD STEP 1 SUBSTEP 136 COMPLETED. CUM ITER = 136
*** TIME = 806.188 TIME INC = 5.02296 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.9763E-01 OSCILLATION LIMIT = 0.4904
*** AUTO TIME STEP: NEXT TIME INC = 5.1216 INCREASED (FACTOR = 1.0196)
*** LOAD STEP 1 SUBSTEP 137 COMPLETED. CUM ITER = 137
*** TIME = 811.309 TIME INC = 5.12156 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.9575E-01 OSCILLATION LIMIT = 0.4904
*** AUTO TIME STEP: NEXT TIME INC = 5.2221 INCREASED (FACTOR = 1.0196)
*** LOAD STEP 1 SUBSTEP 138 COMPLETED. CUM ITER = 138
*** TIME = 816.531 TIME INC = 5.22208 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.9390E-01 OSCILLATION LIMIT = 0.4904
*** AUTO TIME STEP: NEXT TIME INC = 5.3248 INCREASED (FACTOR = 1.0197)
*** LOAD STEP 1 SUBSTEP 139 COMPLETED. CUM ITER = 139
*** TIME = 821.856 TIME INC = 5.32476 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.9208E-01 OSCILLATION LIMIT = 0.4903
*** AUTO TIME STEP: NEXT TIME INC = 5.4299 INCREASED (FACTOR = 1.0197)
*** LOAD STEP 1 SUBSTEP 140 COMPLETED. CUM ITER = 140
*** TIME = 827.286 TIME INC = 5.42987 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.9029E-01 OSCILLATION LIMIT = 0.4903
*** AUTO TIME STEP: NEXT TIME INC = 5.5377 INCREASED (FACTOR = 1.0199)
*** LOAD STEP 1 SUBSTEP 141 COMPLETED. CUM ITER = 141
*** TIME = 832.824 TIME INC = 5.53768 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.8852E-01 OSCILLATION LIMIT = 0.4902
*** AUTO TIME STEP: NEXT TIME INC = 5.6485 INCREASED (FACTOR = 1.0200)
*** LOAD STEP 1 SUBSTEP 142 COMPLETED. CUM ITER = 142
*** TIME = 838.472 TIME INC = 5.64849 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.8677E-01 OSCILLATION LIMIT = 0.4901
*** AUTO TIME STEP: NEXT TIME INC = 5.7626 INCREASED (FACTOR = 1.0202)
*** LOAD STEP 1 SUBSTEP 143 COMPLETED. CUM ITER = 143
*** TIME = 844.235 TIME INC = 5.76257 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.8503E-01 OSCILLATION LIMIT = 0.4900
*** AUTO TIME STEP: NEXT TIME INC = 5.8802 INCREASED (FACTOR = 1.0204)
*** LOAD STEP 1 SUBSTEP 144 COMPLETED. CUM ITER = 144
*** TIME = 850.115 TIME INC = 5.88023 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.8331E-01 OSCILLATION LIMIT = 0.4899
*** AUTO TIME STEP: NEXT TIME INC = 6.0018 INCREASED (FACTOR = 1.0207)
*** LOAD STEP 1 SUBSTEP 145 COMPLETED. CUM ITER = 145
*** TIME = 856.117 TIME INC = 6.00179 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.8160E-01 OSCILLATION LIMIT = 0.4897
*** AUTO TIME STEP: NEXT TIME INC = 6.1275 INCREASED (FACTOR = 1.0210)
*** LOAD STEP 1 SUBSTEP 146 COMPLETED. CUM ITER = 146
*** TIME = 862.244 TIME INC = 6.12754 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.7990E-01 OSCILLATION LIMIT = 0.4896
*** AUTO TIME STEP: NEXT TIME INC = 6.2578 INCREASED (FACTOR = 1.0213)
*** LOAD STEP 1 SUBSTEP 147 COMPLETED. CUM ITER = 147
*** TIME = 868.502 TIME INC = 6.25782 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.7821E-01 OSCILLATION LIMIT = 0.4894
*** AUTO TIME STEP: NEXT TIME INC = 6.3929 INCREASED (FACTOR = 1.0216)
*** LOAD STEP 1 SUBSTEP 148 COMPLETED. CUM ITER = 148
*** TIME = 874.895 TIME INC = 6.39295 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.7653E-01 OSCILLATION LIMIT = 0.4893
*** AUTO TIME STEP: NEXT TIME INC = 6.5333 INCREASED (FACTOR = 1.0219)
*** LOAD STEP 1 SUBSTEP 149 COMPLETED. CUM ITER = 149
*** TIME = 881.428 TIME INC = 6.53326 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.7486E-01 OSCILLATION LIMIT = 0.4891
*** AUTO TIME STEP: NEXT TIME INC = 6.6791 INCREASED (FACTOR = 1.0223)
*** LOAD STEP 1 SUBSTEP 150 COMPLETED. CUM ITER = 150
*** TIME = 888.107 TIME INC = 6.67908 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.7320E-01 OSCILLATION LIMIT = 0.4889
*** AUTO TIME STEP: NEXT TIME INC = 6.8308 INCREASED (FACTOR = 1.0227)
*** LOAD STEP 1 SUBSTEP 151 COMPLETED. CUM ITER = 151
*** TIME = 894.938 TIME INC = 6.83075 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.7154E-01 OSCILLATION LIMIT = 0.4887
*** AUTO TIME STEP: NEXT TIME INC = 5.0620 DECREASED (FACTOR = 0.7411)
*** LOAD STEP 1 SUBSTEP 152 COMPLETED. CUM ITER = 152
*** TIME = 900.000 TIME INC = 5.06196 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.7035E-01 OSCILLATION LIMIT = 0.3561
*** AUTO TIME STEP: NEXT TIME INC = 10.000 INCREASED (FACTOR = 1.9755)
*** LOAD STEP 1 SUBSTEP 153 COMPLETED. CUM ITER = 153
*** TIME = 910.000 TIME INC = 10.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2093E-01 OSCILLATION LIMIT = 0.2093
*** AUTO TIME STEP: NEXT TIME INC = 23.886 INCREASED (FACTOR = 2.3886)
*** LOAD STEP 1 SUBSTEP 154 COMPLETED. CUM ITER = 154
*** TIME = 933.886 TIME INC = 23.8864 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1164E-01 OSCILLATION LIMIT = 0.2781
*** AUTO TIME STEP: NEXT TIME INC = 25.000 INCREASED (FACTOR = 1.0466)
*** LOAD STEP 1 SUBSTEP 155 COMPLETED. CUM ITER = 155
*** TIME = 958.886 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.7726E-02 OSCILLATION LIMIT = 0.1931
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 156 COMPLETED. CUM ITER = 156
*** TIME = 983.886 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.5386E-02 OSCILLATION LIMIT = 0.1346
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 157 COMPLETED. CUM ITER = 157
*** TIME = 1008.89 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.3926E-02 OSCILLATION LIMIT = 0.9816E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 158 COMPLETED. CUM ITER = 158
*** TIME = 1033.89 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.3049E-02 OSCILLATION LIMIT = 0.7624E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 159 COMPLETED. CUM ITER = 159
*** TIME = 1058.89 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2544E-02 OSCILLATION LIMIT = 0.6361E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 160 COMPLETED. CUM ITER = 160
*** TIME = 1083.89 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2262E-02 OSCILLATION LIMIT = 0.5655E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 161 COMPLETED. CUM ITER = 161
*** TIME = 1108.89 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2107E-02 OSCILLATION LIMIT = 0.5268E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 162 COMPLETED. CUM ITER = 162
*** TIME = 1133.89 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.2023E-02 OSCILLATION LIMIT = 0.5059E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 163 COMPLETED. CUM ITER = 163
*** TIME = 1158.89 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1978E-02 OSCILLATION LIMIT = 0.4946E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 164 COMPLETED. CUM ITER = 164
*** TIME = 1183.89 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1954E-02 OSCILLATION LIMIT = 0.4885E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 165 COMPLETED. CUM ITER = 165
*** TIME = 1208.89 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1941E-02 OSCILLATION LIMIT = 0.4853E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 166 COMPLETED. CUM ITER = 166
*** TIME = 1233.89 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1934E-02 OSCILLATION LIMIT = 0.4835E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 167 COMPLETED. CUM ITER = 167
*** TIME = 1258.89 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1930E-02 OSCILLATION LIMIT = 0.4826E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 168 COMPLETED. CUM ITER = 168
*** TIME = 1283.89 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1928E-02 OSCILLATION LIMIT = 0.4821E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 169 COMPLETED. CUM ITER = 169
*** TIME = 1308.89 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1927E-02 OSCILLATION LIMIT = 0.4818E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 170 COMPLETED. CUM ITER = 170
*** TIME = 1333.89 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1927E-02 OSCILLATION LIMIT = 0.4817E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 171 COMPLETED. CUM ITER = 171
*** TIME = 1358.89 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1926E-02 OSCILLATION LIMIT = 0.4816E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 172 COMPLETED. CUM ITER = 172
*** TIME = 1383.89 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1926E-02 OSCILLATION LIMIT = 0.4816E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 173 COMPLETED. CUM ITER = 173
*** TIME = 1408.89 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1926E-02 OSCILLATION LIMIT = 0.4815E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 174 COMPLETED. CUM ITER = 174
*** TIME = 1433.89 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1926E-02 OSCILLATION LIMIT = 0.4815E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 175 COMPLETED. CUM ITER = 175
*** TIME = 1458.89 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1926E-02 OSCILLATION LIMIT = 0.4815E-01
*** AUTO STEP TIME: NEXT TIME INC = 25.000 UNCHANGED
*** LOAD STEP 1 SUBSTEP 176 COMPLETED. CUM ITER = 176
*** TIME = 1483.89 TIME INC = 25.0000 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1926E-02 OSCILLATION LIMIT = 0.4815E-01
*** AUTO TIME STEP: NEXT TIME INC = 16.114 DECREASED (FACTOR = 0.6445)
*** LOAD STEP 1 SUBSTEP 177 COMPLETED. CUM ITER = 177
*** TIME = 1500.00 TIME INC = 16.1136 NEW TRIANG MATRIX
*** RESPONSE EIGENVALUE = 0.1926E-02 OSCILLATION LIMIT = 0.3104E-01
Post-Processing#
Open up the result file using ansys.mapdl.reader
result = mapdl.thermal_result
mapdl.post1()
FINISH SOLUTION PROCESSING
***** ROUTINE COMPLETED ***** CP = 0.000
*****MAPDL VERIFICATION RUN ONLY*****
DO NOT USE RESULTS FOR PRODUCTION
***** MAPDL RESULTS INTERPRETATION (POST1) *****
USE LAST SUBSTEP ON RESULT FILE FOR LOAD CASE 0
SET COMMAND GOT LOAD STEP= 1 SUBSTEP= 177 CUMULATIVE ITERATION= 177
TIME/FREQUENCY= 1500.0
TITLE=
Visualize a Slice#
Visualize a slice through the dataset using pyvista
for more details visit PyVista documentation.
# get the temperature of the 30th result set
mapdl.set(1, 30)
temp = mapdl.post_processing.nodal_temperature()
# Load this result into the underlying VTK grid
grid = mapdl.mesh._grid
grid["temperature"] = temp
# generate a single horizontal slice slice along the XY plane
single_slice = grid.slice(normal=[0, 0, 1], origin=[0, 0, 0.5])
single_slice.plot(scalars="temperature")
Visualize Several Slices#
This shows how you can visualize a series of slices through a dataset
# get the temperature of a different result set
mapdl.set(1, 120)
temp = mapdl.post_processing.nodal_temperature()
# Load this result into the underlying VTK grid
grid = mapdl.mesh._grid
grid["temperature"] = temp
# generate a single horizontal slice slice along the XY plane
slices = grid.slice_along_axis(7, "y")
slices.plot(scalars="temperature", lighting=False, show_edges=True)
Temperature at a Single Point#
Extract the temperature at a single node and plot it with respect to
the input temperatures using ansys.mapdl
Here, we use the get_value
command which is very similar to the
*GET
command, except it immediately returns the value as a
python accessible variable, rather than storing it to a MAPDL value.
# for example, the temperature of Node 12 is can be retrieved simply with:
mapdl.get_value("node", 12, "TEMP")
# note that this is similar to # *GET, Par, NODE, N, Item1, IT1NUM, Item2, IT2NUM
# See the MAPDL reference for all the items you can obtain using *GET
251.05408986663753
Here, we extract the temperature of the node across for each solution
[76.27707675716567, 76.8179922435201, 77.4215154009758, 79.04838835561274, 81.1910830004936, 83.56709007016187, 85.35125295387681, 86.23266429275375, 86.97754017604889, 88.31275589884471]
Alternatively, you can simply grab the data for the node from the entire response. This is less efficient as the entire data set is sent back for each result.
# get the index of node 12 in MAPDL
idx = np.nonzero(mapdl.mesh.nnum == 12)[0][0]
# get the temperature at that index for each result
node_temp_from_post = []
for i in range(1, 1 + nsets):
mapdl.set(1, i)
node_temp_from_post.append(mapdl.post_processing.nodal_temperature()[idx])
# Again, the first 10 temperatures
node_temp_from_post[:10]
[np.float64(76.27707675716567), np.float64(76.8179922435201), np.float64(77.4215154009758), np.float64(79.04838835561274), np.float64(81.1910830004936), np.float64(83.56709007016187), np.float64(85.35125295387681), np.float64(86.23266429275375), np.float64(86.97754017604889), np.float64(88.31275589884471)]
Plot the temperature as a function of time
time_values = mapdl.post_processing.time_values
plt.plot(time_values, node_temp, label="Node 12")
plt.plot(my_bulk[:, 0], my_bulk[:, 1], ":", label="Input")
plt.legend()
plt.xlabel("Time (seconds)")
plt.ylabel("Temperature ($^\circ$F)")
plt.show()
Stop mapdl#
mapdl.exit()
Total running time of the script: (0 minutes 10.029 seconds)