Taylor bar example#

This example is inspired by the “Taylor Bar” example on the LS-DYNA Knowledge Base site. It shows how to use PyDyna to create a keyword file for LS-DYNA and solve it within a Pythonic environment.

Perform required imports#

Import required packages, including those for the keywords, deck, and solver.

import os
import pathlib
import shutil
import tempfile

import ansys.dpf.core as dpf
import matplotlib.pyplot as plt

from ansys.dyna.core import Deck, keywords as kwd
from ansys.dyna.core.run import run_dyna
from ansys.dyna.core.utils.download_utilities import EXAMPLES_PATH, DownloadManager

workdir = tempfile.TemporaryDirectory()

mesh_file_name = "taylor_bar_mesh.k"
mesh_file = DownloadManager().download_file(
    mesh_file_name, "ls-dyna", "Taylor_Bar", destination=os.path.join(EXAMPLES_PATH, "Taylor_Bar")
)

Create a deck and keywords#

Create a deck, which is the container for all the keywords. Then, create and append individual keywords to the deck.

def create_input_deck(initial_velocity):
    deck = Deck()
    deck.title = "Taylor-Bar Velocity - %s - Unit: t-mm-s" % initial_velocity

    # Define material
    mat_1 = kwd.Mat003(mid=1)
    mat_1.ro = 7.85000e-9
    mat_1.e = 150000.0
    mat_1.pr = 0.34
    mat_1.sigy = 390.0
    mat_1.etan = 90.0

    # Define section
    sec_1 = kwd.SectionSolid(secid=1)
    sec_1.elform = 1

    # Define part
    part_1 = kwd.Part(pid=1, mid=mat_1.mid, secid=sec_1.secid)

    # Define coordinate system
    cs_1 = kwd.DefineCoordinateSystem(cid=1)
    cs_1.xl = 1.0
    cs_1.yp = 1.0

    # Define initial velocity
    init_vel = kwd.InitialVelocityGeneration()
    init_vel.id = part_1.parts["pid"][0]
    init_vel.styp = 2
    init_vel.vy = initial_velocity
    init_vel.icid = cs_1.cid

    # Define box for node set
    box_1 = kwd.DefineBox(boxid=1, xmn=-500, xmx=500, ymn=39.0, ymx=40.1, zmn=-500, zmx=500)

    # Create node set
    set_node_1 = kwd.SetNodeGeneral()
    set_node_1.sid = 1
    set_node_1.option = "BOX"
    set_node_1.e1 = box_1.boxid

    # Define rigid wall
    rw = kwd.RigidwallPlanar(id=1)
    rw.nsid = set_node_1.sid
    rw.yt = box_1.ymx
    rw.yh = box_1.ymn

    # Define control termination
    control_term = kwd.ControlTermination(endtim=8.00000e-5, dtmin=0.001)

    # Define database cards
    deck_dt_out = 8.00000e-8
    deck_glstat = kwd.DatabaseGlstat(dt=deck_dt_out, binary=3)
    deck_matsum = kwd.DatabaseMatsum(dt=deck_dt_out, binary=3)
    deck_nodout = kwd.DatabaseNodout(dt=deck_dt_out, binary=3)
    deck_elout = kwd.DatabaseElout(dt=deck_dt_out, binary=3)
    deck_rwforc = kwd.DatabaseRwforc(dt=deck_dt_out, binary=3)
    deck_d3plot = kwd.DatabaseBinaryD3Plot(dt=4.00000e-6)

    # Define deck history node
    deck_hist_node_1 = kwd.DatabaseHistoryNodeSet()
    deck_hist_node_1.id1 = set_node_1.sid

    # Append all cards to input deck
    deck.extend(
        [
            deck_glstat,
            deck_matsum,
            deck_nodout,
            deck_elout,
            deck_rwforc,
            deck_d3plot,
            set_node_1,
            control_term,
            rw,
            box_1,
            init_vel,
            cs_1,
            part_1,
            mat_1,
            sec_1,
            deck_hist_node_1,
        ]
    )

    return deck


def write_input_deck(**kwargs):
    initial_velocity = kwargs.get("initial_velocity")
    wd = kwargs.get("wd")
    if not all((initial_velocity, wd)):
        raise Exception("Missing input!")
    deck = create_input_deck(initial_velocity)
    # Import mesh
    deck.append(kwd.Include(filename=mesh_file_name))

    # Write LS-DYNA input deck
    os.makedirs(wd, exist_ok=True)
    deck.export_file(os.path.join(wd, "input.k"))
    shutil.copyfile(mesh_file, os.path.join(wd, mesh_file_name))

Define the Dyna solver function#

def run(directory):
    run_dyna("input.k", working_directory=directory, stream=False)
    assert os.path.isfile(os.path.join(directory, "d3plot")), "No result file found"

Define the DPF output function#

def get_global_ke(directory):
    ds = dpf.DataSources()
    result_file = os.path.join(directory, "d3plot")
    assert os.path.isfile(result_file)
    ds.set_result_file_path(result_file, "d3plot")
    model = dpf.Model(ds)

    gke_op = dpf.operators.result.global_kinetic_energy()
    gke_op.inputs.data_sources.connect(ds)
    gke = gke_op.eval()
    field = gke.get_field(0)
    ke_data = field.data

    time_data = model.metadata.time_freq_support.time_frequencies.data_as_list
    return time_data, ke_data

View the model#

etc etc

deck_for_graphic = create_input_deck(300e3)
deck_for_graphic.append(kwd.Include(filename=mesh_file))
deck_for_graphic.plot()
plot taylor bar
/home/runner/work/pydyna/pydyna/.venv/lib/python3.13/site-packages/ansys/dyna/core/lib/deck_plotter.py:567: PyVistaFutureWarning: The default value of `algorithm` for the filter
`UnstructuredGrid.extract_surface` will change in the future. It currently defaults to
`'dataset_surface'`, but will change to `None`. Explicitly set the `algorithm` keyword to
silence this warning.
  plot_data = plot_data.extract_surface()

Run a parametric solve#

etc etc

# Define base working directory

color = ["b", "r", "g", "y"]
# Specify different velocities in mm/s
initial_velocities = [275.0e3, 300.0e3, 325.0e3, 350.0e3]

for index, initial_velocity in enumerate(initial_velocities):
    # Create a folder for each parameter
    wd = os.path.join(workdir.name, "tb_vel_%s" % initial_velocity)
    pathlib.Path(wd).mkdir(exist_ok=True)
    # Create LS-Dyna input deck
    write_input_deck(initial_velocity=initial_velocity, wd=wd)
    # Run Solver
    run(wd)
    # Run PyDPF Post
    time_data, ke_data = get_global_ke(wd)
    # Add series to the plot
    plt.plot(time_data, ke_data, color[index], label="KE at vel. %s mm/s" % initial_velocity)

plt.xlabel("Time (s)")
plt.ylabel("Energy (mJ)")
 License option : check ansys licenses only


 ***************************************************************
 *            ANSYS LEGAL NOTICES                              *
 ***************************************************************
 *                                                             *
 * Copyright 1971-2025 ANSYS, Inc.  All rights reserved.       *
 * Unauthorized use, distribution or duplication is            *
 * prohibited.                                                 *
 *                                                             *
 * Ansys is a registered trademark of ANSYS, Inc. or its       *
 * subsidiaries in the United States or other countries.       *
 * See the ANSYS, Inc. online documentation or the ANSYS, Inc. *
 * documentation CD or online help for the complete Legal      *
 * Notice.                                                     *
 *                                                             *
 ***************************************************************
 *                                                             *
 * THIS ANSYS SOFTWARE PRODUCT AND PROGRAM DOCUMENTATION       *
 * INCLUDE TRADE SECRETS AND CONFIDENTIAL AND PROPRIETARY      *
 * PRODUCTS OF ANSYS, INC., ITS SUBSIDIARIES, OR LICENSORS.    *
 * The software products and documentation are furnished by    *
 * ANSYS, Inc. or its subsidiaries under a software license    *
 * agreement that contains provisions concerning               *
 * non-disclosure, copying, length and nature of use,          *
 * compliance with exporting laws, warranties, disclaimers,    *
 * limitations of liability, and remedies, and other           *
 * provisions.  The software products and documentation may be *
 * used, disclosed, transferred, or copied only in accordance  *
 * with the terms and conditions of that software license      *
 * agreement.                                                  *
 *                                                             *
 * ANSYS, Inc. is a UL registered                              *
 * ISO 9001:2015 company.                                      *
 *                                                             *
 ***************************************************************
 *                                                             *
 * This product is subject to U.S. laws governing export and   *
 * re-export.                                                  *
 *                                                             *
 * For U.S. Government users, except as specifically granted   *
 * by the ANSYS, Inc. software license agreement, the use,     *
 * duplication, or disclosure by the United States Government  *
 * is subject to restrictions stated in the ANSYS, Inc.        *
 * software license agreement and FAR 12.212 (for non-DOD      *
 * licenses).                                                  *
 *                                                             *
 ***************************************************************

      Date: 03/05/2026      Time: 00:28:30
      ___________________________________________________
     |                                                   |
     |  LS-DYNA, A Program for Nonlinear Dynamic         |
     |  Analysis of Structures in Three Dimensions       |
     |  Date    : 10/01/2025    Time: 13:03:42           |
     |  Version : smp d R16                              |
     |  Revision: R16.1.1-20-g0c90cad538                 |
     |  AnLicVer: 2025 R2 (20250506+dl-73-g4dde854)      |
     |                                                   |
     |  Features enabled in this version:                |
     |    Shared Memory Parallel                         |
     |    CESE CHEMISTRY EM ICFD STOCHASTIC_PARTICLES    |
     |    FFTW (multi-dimensional FFTW Library)          |
     |    ARPACK (nonsymmetric eigensolver library)      |
     |    ANSYSLIC enabled                               |
     |                                                   |
     |  Platform   : Xeon64 System                       |
     |  OS Level   : Linux 3.10.0 uum                    |
     |  Compiler   : Intel Fortran Compiler 19.0 SSE2    |
     |  Hostname   : c5f74e11fadc                        |
     |  Precision  : Double precision (I8R8)             |
     |                                                   |
     |  Unauthorized use infringes Ansys Inc. copyrights |
     |___________________________________________________|

 > i=input.k memory=20m
Messages file /ansys_inc/v252/licensingclient/language/en-us/ansysli_msgs.xml does not exist.
[license/info] Successfully checked out 1 of "dyna_solver_core".
[license/info] --> Checkout ID: c5f74e11fadc-root-13-000003 (days left: 320)
[license/info] --> Customer ID: 0
[license/info] Successfully started "LSDYNA (Core-based License)".

 Executing with ANSYS license

 Command line options: i=input.k
                       memory=20m

 Input file: input.k

 The native file format       : 64-bit small endian
 Memory size from command line:    20000000

 on UNIX computers note the following change:

 ctrl-c interrupts ls-dyna and prompts for a sense  switch.
 type the desired sense switch: sw1., sw2., etc. to continue
 the execution.  ls-dyna will respond as explained in the users manual

    type                      response
   -----   ------------------------------------------------------------
   quit    ls-dyna terminates.
   stop    ls-dyna terminates.
   sw1.    a restart file is written and ls-dyna terminates.
   sw2.    ls-dyna responds with time and cycle numbers.
   sw3.    a restart file is written and ls-dyna continues calculations.
   sw4.    a plot state is written and ls-dyna continues calculations.
   sw5.    ls-dyna enters interactive graphics phase.
   swa.    ls-dyna flushes all output i/o buffers.
   swb.    a dynain is written and ls-dyna continues calculations.
   swc.    a restart and dynain are written and ls-dyna continues calculations.
   swd.    a restart and dynain are written and ls-dyna terminates.
   swe.    stop dynamic relaxation just as though convergence
   endtime=time change the termination time
   lpri    toggle implicit lin. alg. solver output on/off.
   nlpr    toggle implicit nonlinear solver output on/off.
   iter    toggle implicit output to d3iter database on/off.
   prof    output timing data to messag and continue.
   conv    force implicit nonlinear convergence for current time step.
   ttrm    terminate implicit time step, reduce time step, retry time step.
   rtrm    terminate implicit at end of current time step.


 ********  notice  ********  notice  ********  notice  ********
 *                                                            *
 * This is the LS-DYNA Finite Element code.                   *
 *                                                            *
 * Neither LST nor the authors assume any responsibility for  *
 * the validity, accuracy, or applicability of any results    *
 * obtained from this system. Users must verify their own     *
 * results.                                                   *
 *                                                            *
 * LST endeavors to make the LS-DYNA code as complete,        *
 * accurate and easy to use as possible.                      *
 * Suggestions and comments are welcomed.  Please report any  *
 * errors encountered in either the documentation or results  *
 * immediately to LST through your site focus.                *
 *                                                            *
 * Copyright (C) 1990-2021                                    *
 * by Livermore Software Technology, LLC                      *
 * All rights reserved                                        *
 *                                                            *
 ********  notice  ********  notice  ********  notice  ********

 Beginning of keyword reader                                   03/05/26 00:28:36

                                                               03/05/26 00:28:36
 Open include file: taylor_bar_mesh.k

 Memory required to process keyword     :     195795
 Additional dynamic memory required     :    2200515


 input of data is completed

 initial kinetic energy = 0.11118392E+06

 Memory required to begin solution      :     196K
 Additional dynamically allocated memory:    2257K
                                   Total:    2453K

 initialization completed
       1 t 0.0000E+00 dt 6.72E-08 flush i/o buffers            03/05/26 00:28:36
       1 t 0.0000E+00 dt 6.72E-08 write d3plot file            03/05/26 00:28:36
 cpu time per zone cycle............       864 nanoseconds
 average cpu time per zone cycle....       911 nanoseconds
 average clock time per zone cycle..      1129 nanoseconds

 estimated total cpu time          =         0 sec (       0 hrs  0 mins)
 estimated cpu time to complete    =         0 sec (       0 hrs  0 mins)
 estimated total clock time        =         6 sec (       0 hrs  0 mins)
 estimated clock time to complete  =         1 sec (       0 hrs  0 mins)
 termination time                  = 8.000E-05
      60 t 3.9547E-06 dt 6.71E-08 write d3plot file            03/05/26 00:28:36
     136 t 7.9703E-06 dt 3.83E-08 write d3plot file            03/05/26 00:28:36
     272 t 1.1990E-05 dt 2.37E-08 write d3plot file            03/05/26 00:28:36
     493 t 1.5998E-05 dt 1.53E-08 write d3plot file            03/05/26 00:28:36
     812 t 1.9999E-05 dt 1.10E-08 write d3plot file            03/05/26 00:28:36
    1234 t 2.3997E-05 dt 8.71E-09 write d3plot file            03/05/26 00:28:36
    1734 t 2.7997E-05 dt 7.45E-09 write d3plot file            03/05/26 00:28:36
    2302 t 3.1995E-05 dt 6.80E-09 write d3plot file            03/05/26 00:28:36
    2905 t 3.5998E-05 dt 6.54E-09 write d3plot file            03/05/26 00:28:36
    3526 t 3.9999E-05 dt 6.52E-09 write d3plot file            03/05/26 00:28:37
    4157 t 4.4000E-05 dt 6.15E-09 write d3plot file            03/05/26 00:28:37
    4807 t 4.7995E-05 dt 6.15E-09 write d3plot file            03/05/26 00:28:37
    5000 t 4.9198E-05 dt 6.53E-09 flush i/o buffers            03/05/26 00:28:37
    5450 t 5.1999E-05 dt 6.15E-09 write d3plot file            03/05/26 00:28:37
    6095 t 5.5999E-05 dt 6.15E-09 write d3plot file            03/05/26 00:28:37
    6722 t 5.9997E-05 dt 6.15E-09 write d3plot file            03/05/26 00:28:38
    7340 t 6.3996E-05 dt 6.53E-09 write d3plot file            03/05/26 00:28:38
    7963 t 6.7995E-05 dt 6.53E-09 write d3plot file            03/05/26 00:28:38
    8579 t 7.1994E-05 dt 6.16E-09 write d3plot file            03/05/26 00:28:38
    9223 t 7.5998E-05 dt 6.57E-09 write d3plot file            03/05/26 00:28:38
    9835 t 7.9994E-05 dt 6.57E-09 write d3plot file            03/05/26 00:28:38

 *** termination time reached ***
    9835 t 8.0001E-05 dt 6.57E-09 write d3dump01 file          03/05/26 00:28:38
    9835 t 8.0001E-05 dt 6.57E-09 write d3plot file            03/05/26 00:28:38

 N o r m a l    t e r m i n a t i o n                          03/05/26 00:28:38

 Memory required to complete solution   :     196K
 Additional dynamically allocated memory:    2689K
                                   Total:    2884K

 T i m i n g   i n f o r m a t i o n
                        CPU(seconds)   %CPU  Clock(seconds) %Clock
  ----------------------------------------------------------------
  Keyword Processing ... 1.5214E-02    0.52     1.5224E-02    0.18
    KW Reading ......... 5.8835E-03    0.20     5.8940E-03    0.07
    KW Writing ......... 2.7163E-03    0.09     2.7170E-03    0.03
  Initialization ....... 1.1635E-01    4.01     5.7518E+00   67.32
    Init Proc Phase 1 .. 5.4813E-02    1.89     5.5690E-02    0.65
    Init Proc Phase 2 .. 1.1077E-03    0.04     1.6590E-03    0.02
  Init solver .......... 5.3880E-05    0.00     5.4000E-05    0.00
  Element processing ... 1.8239E+00   62.88     1.8239E+00   21.35
    Solids ............. 1.7313E+00   59.69     1.7290E+00   20.24
    Shells ............. 2.6205E-02    0.90     2.8477E-02    0.33
    ISO Shells ......... 9.4960E-03    0.33     9.5570E-03    0.11
    E Other ............ 4.2378E-02    1.46     4.2173E-02    0.49
  Binary databases ..... 1.6844E-02    0.58     1.7114E-02    0.20
  ASCII database ....... 4.0631E-01   14.01     4.0611E-01    4.75
  Contact algorithm .... 2.3115E-02    0.80     2.3339E-02    0.27
  Rigid Bodies ......... 4.4223E-02    1.52     4.4425E-02    0.52
    Mass Scaling ....... 1.0974E-02    0.38     1.1048E-02    0.13
    Update nodes ....... 1.0142E-02    0.35     1.0197E-02    0.12
  Time step size ....... 1.2269E-02    0.42     1.2259E-02    0.14
  Rigid wall ........... 3.8169E-02    1.32     3.7436E-02    0.44
  Group force file ..... 1.0388E-02    0.36     1.0567E-02    0.12
  Others ............... 2.6498E-02    0.91     2.6186E-02    0.31
  Misc. 1 .............. 1.9253E-01    6.64     2.0005E-01    2.34
    Scale Masses ....... 1.0028E-02    0.35     1.0023E-02    0.12
    Force Constraints .. 9.9771E-03    0.34     9.9240E-03    0.12
    Force to Accel ..... 3.8689E-02    1.33     3.8464E-02    0.45
  Misc. 2 .............. 9.0587E-02    3.12     9.0466E-02    1.06
    Geometry Update .... 5.5900E-02    1.93     5.5238E-02    0.65
  Misc. 3 .............. 2.0084E-02    0.69     2.0200E-02    0.24
  Misc. 4 .............. 6.3986E-02    2.21     6.4230E-02    0.75
    Timestep Init ...... 2.6463E-02    0.91     2.6300E-02    0.31
    Apply Loads ........ 1.6396E-02    0.57     1.6403E-02    0.19
  ----------------------------------------------------------------
  T o t a l s            2.9005E+00  100.00     8.5434E+00  100.00

 Problem time       =    8.0001E-05
 Problem cycle      =      9835
 Total CPU time     =         3 seconds (   0 hours  0 minutes  3 seconds)
 CPU time per zone cycle  =        366.110 nanoseconds
 Clock time per zone cycle=        367.095 nanoseconds

 Number of CPU's    1
 NLQ used/max     192/  192
 Start time   03/05/2026 00:28:36
 End time     03/05/2026 00:28:38
 Elapsed time       2 seconds for    9835 cycles using  1 SMP thread
             (      0 hour   0 minute   2 seconds)

 N o r m a l    t e r m i n a t i o n                          03/05/26 00:28:38
 License option : check ansys licenses only


 ***************************************************************
 *            ANSYS LEGAL NOTICES                              *
 ***************************************************************
 *                                                             *
 * Copyright 1971-2025 ANSYS, Inc.  All rights reserved.       *
 * Unauthorized use, distribution or duplication is            *
 * prohibited.                                                 *
 *                                                             *
 * Ansys is a registered trademark of ANSYS, Inc. or its       *
 * subsidiaries in the United States or other countries.       *
 * See the ANSYS, Inc. online documentation or the ANSYS, Inc. *
 * documentation CD or online help for the complete Legal      *
 * Notice.                                                     *
 *                                                             *
 ***************************************************************
 *                                                             *
 * THIS ANSYS SOFTWARE PRODUCT AND PROGRAM DOCUMENTATION       *
 * INCLUDE TRADE SECRETS AND CONFIDENTIAL AND PROPRIETARY      *
 * PRODUCTS OF ANSYS, INC., ITS SUBSIDIARIES, OR LICENSORS.    *
 * The software products and documentation are furnished by    *
 * ANSYS, Inc. or its subsidiaries under a software license    *
 * agreement that contains provisions concerning               *
 * non-disclosure, copying, length and nature of use,          *
 * compliance with exporting laws, warranties, disclaimers,    *
 * limitations of liability, and remedies, and other           *
 * provisions.  The software products and documentation may be *
 * used, disclosed, transferred, or copied only in accordance  *
 * with the terms and conditions of that software license      *
 * agreement.                                                  *
 *                                                             *
 * ANSYS, Inc. is a UL registered                              *
 * ISO 9001:2015 company.                                      *
 *                                                             *
 ***************************************************************
 *                                                             *
 * This product is subject to U.S. laws governing export and   *
 * re-export.                                                  *
 *                                                             *
 * For U.S. Government users, except as specifically granted   *
 * by the ANSYS, Inc. software license agreement, the use,     *
 * duplication, or disclosure by the United States Government  *
 * is subject to restrictions stated in the ANSYS, Inc.        *
 * software license agreement and FAR 12.212 (for non-DOD      *
 * licenses).                                                  *
 *                                                             *
 ***************************************************************

      Date: 03/05/2026      Time: 00:29:00
      ___________________________________________________
     |                                                   |
     |  LS-DYNA, A Program for Nonlinear Dynamic         |
     |  Analysis of Structures in Three Dimensions       |
     |  Date    : 10/01/2025    Time: 13:03:42           |
     |  Version : smp d R16                              |
     |  Revision: R16.1.1-20-g0c90cad538                 |
     |  AnLicVer: 2025 R2 (20250506+dl-73-g4dde854)      |
     |                                                   |
     |  Features enabled in this version:                |
     |    Shared Memory Parallel                         |
     |    CESE CHEMISTRY EM ICFD STOCHASTIC_PARTICLES    |
     |    FFTW (multi-dimensional FFTW Library)          |
     |    ARPACK (nonsymmetric eigensolver library)      |
     |    ANSYSLIC enabled                               |
     |                                                   |
     |  Platform   : Xeon64 System                       |
     |  OS Level   : Linux 3.10.0 uum                    |
     |  Compiler   : Intel Fortran Compiler 19.0 SSE2    |
     |  Hostname   : ee91f4dc9087                        |
     |  Precision  : Double precision (I8R8)             |
     |                                                   |
     |  Unauthorized use infringes Ansys Inc. copyrights |
     |___________________________________________________|

 > i=input.k memory=20m
Messages file /ansys_inc/v252/licensingclient/language/en-us/ansysli_msgs.xml does not exist.
[license/info] Successfully checked out 1 of "dyna_solver_core".
[license/info] --> Checkout ID: ee91f4dc9087-root-13-000003 (days left: 320)
[license/info] --> Customer ID: 0
[license/info] Successfully started "LSDYNA (Core-based License)".

 Executing with ANSYS license

 Command line options: i=input.k
                       memory=20m

 Input file: input.k

 The native file format       : 64-bit small endian
 Memory size from command line:    20000000

 on UNIX computers note the following change:

 ctrl-c interrupts ls-dyna and prompts for a sense  switch.
 type the desired sense switch: sw1., sw2., etc. to continue
 the execution.  ls-dyna will respond as explained in the users manual

    type                      response
   -----   ------------------------------------------------------------
   quit    ls-dyna terminates.
   stop    ls-dyna terminates.
   sw1.    a restart file is written and ls-dyna terminates.
   sw2.    ls-dyna responds with time and cycle numbers.
   sw3.    a restart file is written and ls-dyna continues calculations.
   sw4.    a plot state is written and ls-dyna continues calculations.
   sw5.    ls-dyna enters interactive graphics phase.
   swa.    ls-dyna flushes all output i/o buffers.
   swb.    a dynain is written and ls-dyna continues calculations.
   swc.    a restart and dynain are written and ls-dyna continues calculations.
   swd.    a restart and dynain are written and ls-dyna terminates.
   swe.    stop dynamic relaxation just as though convergence
   endtime=time change the termination time
   lpri    toggle implicit lin. alg. solver output on/off.
   nlpr    toggle implicit nonlinear solver output on/off.
   iter    toggle implicit output to d3iter database on/off.
   prof    output timing data to messag and continue.
   conv    force implicit nonlinear convergence for current time step.
   ttrm    terminate implicit time step, reduce time step, retry time step.
   rtrm    terminate implicit at end of current time step.


 ********  notice  ********  notice  ********  notice  ********
 *                                                            *
 * This is the LS-DYNA Finite Element code.                   *
 *                                                            *
 * Neither LST nor the authors assume any responsibility for  *
 * the validity, accuracy, or applicability of any results    *
 * obtained from this system. Users must verify their own     *
 * results.                                                   *
 *                                                            *
 * LST endeavors to make the LS-DYNA code as complete,        *
 * accurate and easy to use as possible.                      *
 * Suggestions and comments are welcomed.  Please report any  *
 * errors encountered in either the documentation or results  *
 * immediately to LST through your site focus.                *
 *                                                            *
 * Copyright (C) 1990-2021                                    *
 * by Livermore Software Technology, LLC                      *
 * All rights reserved                                        *
 *                                                            *
 ********  notice  ********  notice  ********  notice  ********

 Beginning of keyword reader                                   03/05/26 00:29:06

                                                               03/05/26 00:29:06
 Open include file: taylor_bar_mesh.k

 Memory required to process keyword     :     195795
 Additional dynamic memory required     :    2200515


 input of data is completed

 initial kinetic energy = 0.13231805E+06

 Memory required to begin solution      :     196K
 Additional dynamically allocated memory:    2257K
                                   Total:    2453K

 initialization completed
       1 t 0.0000E+00 dt 6.72E-08 flush i/o buffers            03/05/26 00:29:06
       1 t 0.0000E+00 dt 6.72E-08 write d3plot file            03/05/26 00:29:06
 cpu time per zone cycle............       855 nanoseconds
 average cpu time per zone cycle....       916 nanoseconds
 average clock time per zone cycle..      1120 nanoseconds

 estimated total cpu time          =         0 sec (       0 hrs  0 mins)
 estimated cpu time to complete    =         0 sec (       0 hrs  0 mins)
 estimated total clock time        =         6 sec (       0 hrs  0 mins)
 estimated clock time to complete  =         1 sec (       0 hrs  0 mins)
 termination time                  = 8.000E-05
      60 t 3.9564E-06 dt 6.71E-08 write d3plot file            03/05/26 00:29:06
     145 t 7.9823E-06 dt 3.30E-08 write d3plot file            03/05/26 00:29:06
     308 t 1.1998E-05 dt 1.92E-08 write d3plot file            03/05/26 00:29:06
     587 t 1.6000E-05 dt 1.18E-08 write d3plot file            03/05/26 00:29:06
    1009 t 1.9996E-05 dt 7.72E-09 write d3plot file            03/05/26 00:29:06
    1582 t 2.3999E-05 dt 6.32E-09 write d3plot file            03/05/26 00:29:06
    2293 t 2.7997E-05 dt 5.27E-09 write d3plot file            03/05/26 00:29:06
    3109 t 3.1998E-05 dt 4.69E-09 write d3plot file            03/05/26 00:29:07
    3996 t 3.5997E-05 dt 4.39E-09 write d3plot file            03/05/26 00:29:07
    4921 t 3.9998E-05 dt 4.29E-09 write d3plot file            03/05/26 00:29:07
    5000 t 4.0337E-05 dt 4.29E-09 flush i/o buffers            03/05/26 00:29:07
    5855 t 4.3996E-05 dt 4.29E-09 write d3plot file            03/05/26 00:29:07
    6806 t 4.7998E-05 dt 4.29E-09 write d3plot file            03/05/26 00:29:07
    7782 t 5.1998E-05 dt 4.04E-09 write d3plot file            03/05/26 00:29:08
    8766 t 5.5998E-05 dt 4.04E-09 write d3plot file            03/05/26 00:29:08
    9747 t 5.9999E-05 dt 4.05E-09 write d3plot file            03/05/26 00:29:08
   10000 t 6.1023E-05 dt 4.05E-09 flush i/o buffers            03/05/26 00:29:08
   10699 t 6.3996E-05 dt 4.30E-09 write d3plot file            03/05/26 00:29:08
   11643 t 6.7999E-05 dt 4.30E-09 write d3plot file            03/05/26 00:29:09
   12579 t 7.1996E-05 dt 4.30E-09 write d3plot file            03/05/26 00:29:09
   13556 t 7.6000E-05 dt 4.06E-09 write d3plot file            03/05/26 00:29:09
   14494 t 7.9998E-05 dt 4.32E-09 write d3plot file            03/05/26 00:29:09

 *** termination time reached ***
   14494 t 8.0002E-05 dt 4.32E-09 write d3dump01 file          03/05/26 00:29:09
   14494 t 8.0002E-05 dt 4.32E-09 write d3plot file            03/05/26 00:29:09

 N o r m a l    t e r m i n a t i o n                          03/05/26 00:29:09

 Memory required to complete solution   :     196K
 Additional dynamically allocated memory:    2689K
                                   Total:    2884K

 T i m i n g   i n f o r m a t i o n
                        CPU(seconds)   %CPU  Clock(seconds) %Clock
  ----------------------------------------------------------------
  Keyword Processing ... 1.4786E-02    0.37     1.4786E-02    0.16
    KW Reading ......... 5.7160E-03    0.14     5.7170E-03    0.06
    KW Writing ......... 2.7155E-03    0.07     2.7150E-03    0.03
  Initialization ....... 1.1659E-01    2.94     5.5482E+00   59.00
    Init Proc Phase 1 .. 5.4893E-02    1.39     5.5989E-02    0.60
    Init Proc Phase 2 .. 1.1118E-03    0.03     1.6250E-03    0.02
  Init solver .......... 5.5123E-05    0.00     5.5000E-05    0.00
  Element processing ... 2.6498E+00   66.92     2.6499E+00   28.18
    Solids ............. 2.5158E+00   63.54     2.5125E+00   26.72
    Shells ............. 3.8102E-02    0.96     4.1403E-02    0.44
    ISO Shells ......... 1.3977E-02    0.35     1.3998E-02    0.15
    E Other ............ 6.0794E-02    1.54     6.0562E-02    0.64
  Binary databases ..... 2.1889E-02    0.55     2.1979E-02    0.23
  ASCII database ....... 4.2168E-01   10.65     4.2148E-01    4.48
  Contact algorithm .... 3.3512E-02    0.85     3.3724E-02    0.36
  Rigid Bodies ......... 6.4577E-02    1.63     6.4761E-02    0.69
    Mass Scaling ....... 1.6086E-02    0.41     1.6091E-02    0.17
    Update nodes ....... 1.4845E-02    0.37     1.4911E-02    0.16
  Time step size ....... 1.7708E-02    0.45     1.7756E-02    0.19
  Rigid wall ........... 5.4297E-02    1.37     5.3348E-02    0.57
  Group force file ..... 1.5167E-02    0.38     1.5400E-02    0.16
  Others ............... 3.7891E-02    0.96     3.7767E-02    0.40
  Misc. 1 .............. 2.6035E-01    6.58     2.7356E-01    2.91
    Scale Masses ....... 1.4772E-02    0.37     1.4783E-02    0.16
    Force Constraints .. 1.4420E-02    0.36     1.4458E-02    0.15
    Force to Accel ..... 5.6725E-02    1.43     5.6553E-02    0.60
  Misc. 2 .............. 1.3285E-01    3.36     1.3253E-01    1.41
    Geometry Update .... 8.2007E-02    2.07     8.0864E-02    0.86
  Misc. 3 .............. 2.6189E-02    0.66     2.6340E-02    0.28
  Misc. 4 .............. 9.2166E-02    2.33     9.2599E-02    0.98
    Timestep Init ...... 3.8370E-02    0.97     3.8213E-02    0.41
    Apply Loads ........ 2.3136E-02    0.58     2.2931E-02    0.24
  ----------------------------------------------------------------
  T o t a l s            3.9595E+00  100.00     9.4042E+00  100.00

 Problem time       =    8.0002E-05
 Problem cycle      =     14494
 Total CPU time     =         4 seconds (   0 hours  0 minutes  4 seconds)
 CPU time per zone cycle  =        343.458 nanoseconds
 Clock time per zone cycle=        344.633 nanoseconds

 Number of CPU's    1
 NLQ used/max     192/  192
 Start time   03/05/2026 00:29:06
 End time     03/05/2026 00:29:09
 Elapsed time       3 seconds for   14494 cycles using  1 SMP thread
             (      0 hour   0 minute   3 seconds)

 N o r m a l    t e r m i n a t i o n                          03/05/26 00:29:09
 License option : check ansys licenses only


 ***************************************************************
 *            ANSYS LEGAL NOTICES                              *
 ***************************************************************
 *                                                             *
 * Copyright 1971-2025 ANSYS, Inc.  All rights reserved.       *
 * Unauthorized use, distribution or duplication is            *
 * prohibited.                                                 *
 *                                                             *
 * Ansys is a registered trademark of ANSYS, Inc. or its       *
 * subsidiaries in the United States or other countries.       *
 * See the ANSYS, Inc. online documentation or the ANSYS, Inc. *
 * documentation CD or online help for the complete Legal      *
 * Notice.                                                     *
 *                                                             *
 ***************************************************************
 *                                                             *
 * THIS ANSYS SOFTWARE PRODUCT AND PROGRAM DOCUMENTATION       *
 * INCLUDE TRADE SECRETS AND CONFIDENTIAL AND PROPRIETARY      *
 * PRODUCTS OF ANSYS, INC., ITS SUBSIDIARIES, OR LICENSORS.    *
 * The software products and documentation are furnished by    *
 * ANSYS, Inc. or its subsidiaries under a software license    *
 * agreement that contains provisions concerning               *
 * non-disclosure, copying, length and nature of use,          *
 * compliance with exporting laws, warranties, disclaimers,    *
 * limitations of liability, and remedies, and other           *
 * provisions.  The software products and documentation may be *
 * used, disclosed, transferred, or copied only in accordance  *
 * with the terms and conditions of that software license      *
 * agreement.                                                  *
 *                                                             *
 * ANSYS, Inc. is a UL registered                              *
 * ISO 9001:2015 company.                                      *
 *                                                             *
 ***************************************************************
 *                                                             *
 * This product is subject to U.S. laws governing export and   *
 * re-export.                                                  *
 *                                                             *
 * For U.S. Government users, except as specifically granted   *
 * by the ANSYS, Inc. software license agreement, the use,     *
 * duplication, or disclosure by the United States Government  *
 * is subject to restrictions stated in the ANSYS, Inc.        *
 * software license agreement and FAR 12.212 (for non-DOD      *
 * licenses).                                                  *
 *                                                             *
 ***************************************************************

      Date: 03/05/2026      Time: 00:29:32
      ___________________________________________________
     |                                                   |
     |  LS-DYNA, A Program for Nonlinear Dynamic         |
     |  Analysis of Structures in Three Dimensions       |
     |  Date    : 10/01/2025    Time: 13:03:42           |
     |  Version : smp d R16                              |
     |  Revision: R16.1.1-20-g0c90cad538                 |
     |  AnLicVer: 2025 R2 (20250506+dl-73-g4dde854)      |
     |                                                   |
     |  Features enabled in this version:                |
     |    Shared Memory Parallel                         |
     |    CESE CHEMISTRY EM ICFD STOCHASTIC_PARTICLES    |
     |    FFTW (multi-dimensional FFTW Library)          |
     |    ARPACK (nonsymmetric eigensolver library)      |
     |    ANSYSLIC enabled                               |
     |                                                   |
     |  Platform   : Xeon64 System                       |
     |  OS Level   : Linux 3.10.0 uum                    |
     |  Compiler   : Intel Fortran Compiler 19.0 SSE2    |
     |  Hostname   : 6f39724c9d35                        |
     |  Precision  : Double precision (I8R8)             |
     |                                                   |
     |  Unauthorized use infringes Ansys Inc. copyrights |
     |___________________________________________________|

 > i=input.k memory=20m
Messages file /ansys_inc/v252/licensingclient/language/en-us/ansysli_msgs.xml does not exist.
[license/info] Successfully checked out 1 of "dyna_solver_core".
[license/info] --> Checkout ID: 6f39724c9d35-root-13-000003 (days left: 320)
[license/info] --> Customer ID: 0
[license/info] Successfully started "LSDYNA (Core-based License)".

 Executing with ANSYS license

 Command line options: i=input.k
                       memory=20m

 Input file: input.k

 The native file format       : 64-bit small endian
 Memory size from command line:    20000000

 on UNIX computers note the following change:

 ctrl-c interrupts ls-dyna and prompts for a sense  switch.
 type the desired sense switch: sw1., sw2., etc. to continue
 the execution.  ls-dyna will respond as explained in the users manual

    type                      response
   -----   ------------------------------------------------------------
   quit    ls-dyna terminates.
   stop    ls-dyna terminates.
   sw1.    a restart file is written and ls-dyna terminates.
   sw2.    ls-dyna responds with time and cycle numbers.
   sw3.    a restart file is written and ls-dyna continues calculations.
   sw4.    a plot state is written and ls-dyna continues calculations.
   sw5.    ls-dyna enters interactive graphics phase.
   swa.    ls-dyna flushes all output i/o buffers.
   swb.    a dynain is written and ls-dyna continues calculations.
   swc.    a restart and dynain are written and ls-dyna continues calculations.
   swd.    a restart and dynain are written and ls-dyna terminates.
   swe.    stop dynamic relaxation just as though convergence
   endtime=time change the termination time
   lpri    toggle implicit lin. alg. solver output on/off.
   nlpr    toggle implicit nonlinear solver output on/off.
   iter    toggle implicit output to d3iter database on/off.
   prof    output timing data to messag and continue.
   conv    force implicit nonlinear convergence for current time step.
   ttrm    terminate implicit time step, reduce time step, retry time step.
   rtrm    terminate implicit at end of current time step.


 ********  notice  ********  notice  ********  notice  ********
 *                                                            *
 * This is the LS-DYNA Finite Element code.                   *
 *                                                            *
 * Neither LST nor the authors assume any responsibility for  *
 * the validity, accuracy, or applicability of any results    *
 * obtained from this system. Users must verify their own     *
 * results.                                                   *
 *                                                            *
 * LST endeavors to make the LS-DYNA code as complete,        *
 * accurate and easy to use as possible.                      *
 * Suggestions and comments are welcomed.  Please report any  *
 * errors encountered in either the documentation or results  *
 * immediately to LST through your site focus.                *
 *                                                            *
 * Copyright (C) 1990-2021                                    *
 * by Livermore Software Technology, LLC                      *
 * All rights reserved                                        *
 *                                                            *
 ********  notice  ********  notice  ********  notice  ********

 Beginning of keyword reader                                   03/05/26 00:29:38

                                                               03/05/26 00:29:38
 Open include file: taylor_bar_mesh.k

 Memory required to process keyword     :     195795
 Additional dynamic memory required     :    2200515


 input of data is completed

 initial kinetic energy = 0.15528993E+06

 Memory required to begin solution      :     196K
 Additional dynamically allocated memory:    2257K
                                   Total:    2453K

 initialization completed
       1 t 0.0000E+00 dt 6.72E-08 flush i/o buffers            03/05/26 00:29:38
       1 t 0.0000E+00 dt 6.72E-08 write d3plot file            03/05/26 00:29:38
 cpu time per zone cycle............       850 nanoseconds
 average cpu time per zone cycle....       927 nanoseconds
 average clock time per zone cycle..      1126 nanoseconds

 estimated total cpu time          =         0 sec (       0 hrs  0 mins)
 estimated cpu time to complete    =         0 sec (       0 hrs  0 mins)
 estimated total clock time        =         6 sec (       0 hrs  0 mins)
 estimated clock time to complete  =         1 sec (       0 hrs  0 mins)
 termination time                  = 8.000E-05
      60 t 3.9415E-06 dt 6.22E-08 write d3plot file            03/05/26 00:29:38
     155 t 7.9752E-06 dt 2.87E-08 write d3plot file            03/05/26 00:29:38
     349 t 1.1998E-05 dt 1.58E-08 write d3plot file            03/05/26 00:29:38
     697 t 1.5994E-05 dt 9.29E-09 write d3plot file            03/05/26 00:29:38
    1238 t 1.9994E-05 dt 6.27E-09 write d3plot file            03/05/26 00:29:38
    1998 t 2.3998E-05 dt 4.74E-09 write d3plot file            03/05/26 00:29:38
    2943 t 2.7998E-05 dt 3.88E-09 write d3plot file            03/05/26 00:29:38
    4058 t 3.1998E-05 dt 3.18E-09 write d3plot file            03/05/26 00:29:39
    5000 t 3.5037E-05 dt 2.96E-09 flush i/o buffers            03/05/26 00:29:39
    5313 t 3.5998E-05 dt 3.09E-09 write d3plot file            03/05/26 00:29:39
    6657 t 3.9998E-05 dt 2.95E-09 write d3plot file            03/05/26 00:29:39
    8032 t 4.3999E-05 dt 2.91E-09 write d3plot file            03/05/26 00:29:40
    9438 t 4.7999E-05 dt 2.91E-09 write d3plot file            03/05/26 00:29:40
   10000 t 4.9634E-05 dt 2.91E-09 flush i/o buffers            03/05/26 00:29:40
   10861 t 5.1999E-05 dt 2.74E-09 write d3plot file            03/05/26 00:29:40
   12312 t 5.5999E-05 dt 2.74E-09 write d3plot file            03/05/26 00:29:41
   13771 t 5.9998E-05 dt 2.74E-09 write d3plot file            03/05/26 00:29:41
   15000 t 6.3551E-05 dt 2.74E-09 flush i/o buffers            03/05/26 00:29:42
   15163 t 6.3998E-05 dt 2.74E-09 write d3plot file            03/05/26 00:29:42
   16550 t 6.7998E-05 dt 2.91E-09 write d3plot file            03/05/26 00:29:42
   17924 t 7.1999E-05 dt 2.91E-09 write d3plot file            03/05/26 00:29:42
   19360 t 7.6000E-05 dt 2.75E-09 write d3plot file            03/05/26 00:29:43
   20000 t 7.7767E-05 dt 2.76E-09 flush i/o buffers            03/05/26 00:29:43
   20766 t 7.9998E-05 dt 2.93E-09 write d3plot file            03/05/26 00:29:43

 *** termination time reached ***
   20766 t 8.0001E-05 dt 2.93E-09 write d3dump01 file          03/05/26 00:29:43
   20766 t 8.0001E-05 dt 2.93E-09 write d3plot file            03/05/26 00:29:43

 N o r m a l    t e r m i n a t i o n                          03/05/26 00:29:43

 Memory required to complete solution   :     196K
 Additional dynamically allocated memory:    2689K
                                   Total:    2884K

 T i m i n g   i n f o r m a t i o n
                        CPU(seconds)   %CPU  Clock(seconds) %Clock
  ----------------------------------------------------------------
  Keyword Processing ... 1.4888E-02    0.28     1.4889E-02    0.14
    KW Reading ......... 5.8295E-03    0.11     5.8300E-03    0.05
    KW Writing ......... 2.7425E-03    0.05     2.7430E-03    0.03
  Initialization ....... 1.1802E-01    2.18     5.6440E+00   51.51
    Init Proc Phase 1 .. 5.4940E-02    1.02     5.6108E-02    0.51
    Init Proc Phase 2 .. 1.2402E-03    0.02     1.8840E-03    0.02
  Init solver .......... 5.6085E-05    0.00     5.7000E-05    0.00
  Element processing ... 3.7782E+00   69.83     3.7782E+00   34.48
    Solids ............. 3.5886E+00   66.32     3.5843E+00   32.71
    Shells ............. 5.3571E-02    0.99     5.7701E-02    0.53
    ISO Shells ......... 1.9960E-02    0.37     2.0063E-02    0.18
    E Other ............ 8.6090E-02    1.59     8.5527E-02    0.78
  Binary databases ..... 2.9106E-02    0.54     2.9328E-02    0.27
  ASCII database ....... 4.4935E-01    8.30     4.4926E-01    4.10
  Contact algorithm .... 4.8243E-02    0.89     4.8510E-02    0.44
  Rigid Bodies ......... 9.2604E-02    1.71     9.2900E-02    0.85
    Mass Scaling ....... 2.3067E-02    0.43     2.3264E-02    0.21
    Update nodes ....... 2.1190E-02    0.39     2.1269E-02    0.19
  Time step size ....... 2.5224E-02    0.47     2.5312E-02    0.23
  Rigid wall ........... 7.6503E-02    1.41     7.5217E-02    0.69
  Group force file ..... 2.1570E-02    0.40     2.1827E-02    0.20
  Others ............... 5.3743E-02    0.99     5.3344E-02    0.49
  Misc. 1 .............. 3.5394E-01    6.54     3.7482E-01    3.42
    Scale Masses ....... 2.1019E-02    0.39     2.1058E-02    0.19
    Force Constraints .. 2.0751E-02    0.38     2.0690E-02    0.19
    Force to Accel ..... 8.0459E-02    1.49     8.0190E-02    0.73
  Misc. 2 .............. 1.8366E-01    3.39     1.8324E-01    1.67
    Geometry Update .... 1.1197E-01    2.07     1.1070E-01    1.01
  Misc. 3 .............. 3.4902E-02    0.65     3.4652E-02    0.32
  Misc. 4 .............. 1.3078E-01    2.42     1.3119E-01    1.20
    Timestep Init ...... 5.4314E-02    1.00     5.4054E-02    0.49
    Apply Loads ........ 3.2913E-02    0.61     3.2808E-02    0.30
  ----------------------------------------------------------------
  T o t a l s            5.4108E+00  100.00     1.0957E+01  100.00

 Problem time       =    8.0001E-05
 Problem cycle      =     20766
 Total CPU time     =         5 seconds (   0 hours  0 minutes  5 seconds)
 CPU time per zone cycle  =        330.506 nanoseconds
 Clock time per zone cycle=        331.755 nanoseconds

 Number of CPU's    1
 NLQ used/max     192/  192
 Start time   03/05/2026 00:29:38
 End time     03/05/2026 00:29:43
 Elapsed time       5 seconds for   20766 cycles using  1 SMP thread
             (      0 hour   0 minute   5 seconds)

 N o r m a l    t e r m i n a t i o n                          03/05/26 00:29:43
 License option : check ansys licenses only


 ***************************************************************
 *            ANSYS LEGAL NOTICES                              *
 ***************************************************************
 *                                                             *
 * Copyright 1971-2025 ANSYS, Inc.  All rights reserved.       *
 * Unauthorized use, distribution or duplication is            *
 * prohibited.                                                 *
 *                                                             *
 * Ansys is a registered trademark of ANSYS, Inc. or its       *
 * subsidiaries in the United States or other countries.       *
 * See the ANSYS, Inc. online documentation or the ANSYS, Inc. *
 * documentation CD or online help for the complete Legal      *
 * Notice.                                                     *
 *                                                             *
 ***************************************************************
 *                                                             *
 * THIS ANSYS SOFTWARE PRODUCT AND PROGRAM DOCUMENTATION       *
 * INCLUDE TRADE SECRETS AND CONFIDENTIAL AND PROPRIETARY      *
 * PRODUCTS OF ANSYS, INC., ITS SUBSIDIARIES, OR LICENSORS.    *
 * The software products and documentation are furnished by    *
 * ANSYS, Inc. or its subsidiaries under a software license    *
 * agreement that contains provisions concerning               *
 * non-disclosure, copying, length and nature of use,          *
 * compliance with exporting laws, warranties, disclaimers,    *
 * limitations of liability, and remedies, and other           *
 * provisions.  The software products and documentation may be *
 * used, disclosed, transferred, or copied only in accordance  *
 * with the terms and conditions of that software license      *
 * agreement.                                                  *
 *                                                             *
 * ANSYS, Inc. is a UL registered                              *
 * ISO 9001:2015 company.                                      *
 *                                                             *
 ***************************************************************
 *                                                             *
 * This product is subject to U.S. laws governing export and   *
 * re-export.                                                  *
 *                                                             *
 * For U.S. Government users, except as specifically granted   *
 * by the ANSYS, Inc. software license agreement, the use,     *
 * duplication, or disclosure by the United States Government  *
 * is subject to restrictions stated in the ANSYS, Inc.        *
 * software license agreement and FAR 12.212 (for non-DOD      *
 * licenses).                                                  *
 *                                                             *
 ***************************************************************

      Date: 03/05/2026      Time: 00:30:05
      ___________________________________________________
     |                                                   |
     |  LS-DYNA, A Program for Nonlinear Dynamic         |
     |  Analysis of Structures in Three Dimensions       |
     |  Date    : 10/01/2025    Time: 13:03:42           |
     |  Version : smp d R16                              |
     |  Revision: R16.1.1-20-g0c90cad538                 |
     |  AnLicVer: 2025 R2 (20250506+dl-73-g4dde854)      |
     |                                                   |
     |  Features enabled in this version:                |
     |    Shared Memory Parallel                         |
     |    CESE CHEMISTRY EM ICFD STOCHASTIC_PARTICLES    |
     |    FFTW (multi-dimensional FFTW Library)          |
     |    ARPACK (nonsymmetric eigensolver library)      |
     |    ANSYSLIC enabled                               |
     |                                                   |
     |  Platform   : Xeon64 System                       |
     |  OS Level   : Linux 3.10.0 uum                    |
     |  Compiler   : Intel Fortran Compiler 19.0 SSE2    |
     |  Hostname   : d78eb0e84528                        |
     |  Precision  : Double precision (I8R8)             |
     |                                                   |
     |  Unauthorized use infringes Ansys Inc. copyrights |
     |___________________________________________________|

 > i=input.k memory=20m
Messages file /ansys_inc/v252/licensingclient/language/en-us/ansysli_msgs.xml does not exist.
[license/info] Successfully checked out 1 of "dyna_solver_core".
[license/info] --> Checkout ID: d78eb0e84528-root-13-000003 (days left: 320)
[license/info] --> Customer ID: 0
[license/info] Successfully started "LSDYNA (Core-based License)".

 Executing with ANSYS license

 Command line options: i=input.k
                       memory=20m

 Input file: input.k

 The native file format       : 64-bit small endian
 Memory size from command line:    20000000

 on UNIX computers note the following change:

 ctrl-c interrupts ls-dyna and prompts for a sense  switch.
 type the desired sense switch: sw1., sw2., etc. to continue
 the execution.  ls-dyna will respond as explained in the users manual

    type                      response
   -----   ------------------------------------------------------------
   quit    ls-dyna terminates.
   stop    ls-dyna terminates.
   sw1.    a restart file is written and ls-dyna terminates.
   sw2.    ls-dyna responds with time and cycle numbers.
   sw3.    a restart file is written and ls-dyna continues calculations.
   sw4.    a plot state is written and ls-dyna continues calculations.
   sw5.    ls-dyna enters interactive graphics phase.
   swa.    ls-dyna flushes all output i/o buffers.
   swb.    a dynain is written and ls-dyna continues calculations.
   swc.    a restart and dynain are written and ls-dyna continues calculations.
   swd.    a restart and dynain are written and ls-dyna terminates.
   swe.    stop dynamic relaxation just as though convergence
   endtime=time change the termination time
   lpri    toggle implicit lin. alg. solver output on/off.
   nlpr    toggle implicit nonlinear solver output on/off.
   iter    toggle implicit output to d3iter database on/off.
   prof    output timing data to messag and continue.
   conv    force implicit nonlinear convergence for current time step.
   ttrm    terminate implicit time step, reduce time step, retry time step.
   rtrm    terminate implicit at end of current time step.


 ********  notice  ********  notice  ********  notice  ********
 *                                                            *
 * This is the LS-DYNA Finite Element code.                   *
 *                                                            *
 * Neither LST nor the authors assume any responsibility for  *
 * the validity, accuracy, or applicability of any results    *
 * obtained from this system. Users must verify their own     *
 * results.                                                   *
 *                                                            *
 * LST endeavors to make the LS-DYNA code as complete,        *
 * accurate and easy to use as possible.                      *
 * Suggestions and comments are welcomed.  Please report any  *
 * errors encountered in either the documentation or results  *
 * immediately to LST through your site focus.                *
 *                                                            *
 * Copyright (C) 1990-2021                                    *
 * by Livermore Software Technology, LLC                      *
 * All rights reserved                                        *
 *                                                            *
 ********  notice  ********  notice  ********  notice  ********

 Beginning of keyword reader                                   03/05/26 00:30:10

                                                               03/05/26 00:30:10
 Open include file: taylor_bar_mesh.k

 Memory required to process keyword     :     195795
 Additional dynamic memory required     :    2200515


 input of data is completed

 initial kinetic energy = 0.18009957E+06

 Memory required to begin solution      :     196K
 Additional dynamically allocated memory:    2257K
                                   Total:    2453K

 initialization completed
       1 t 0.0000E+00 dt 6.72E-08 flush i/o buffers            03/05/26 00:30:10
       1 t 0.0000E+00 dt 6.72E-08 write d3plot file            03/05/26 00:30:10
 cpu time per zone cycle............       854 nanoseconds
 average cpu time per zone cycle....       912 nanoseconds
 average clock time per zone cycle..      1106 nanoseconds

 estimated total cpu time          =         0 sec (       0 hrs  0 mins)
 estimated cpu time to complete    =         0 sec (       0 hrs  0 mins)
 estimated total clock time        =         6 sec (       0 hrs  0 mins)
 estimated clock time to complete  =         1 sec (       0 hrs  0 mins)
 termination time                  = 8.000E-05
      61 t 3.9737E-06 dt 5.73E-08 write d3plot file            03/05/26 00:30:10
     168 t 7.9894E-06 dt 2.48E-08 write d3plot file            03/05/26 00:30:10
     398 t 1.1989E-05 dt 1.29E-08 write d3plot file            03/05/26 00:30:10
     831 t 1.5993E-05 dt 6.92E-09 write d3plot file            03/05/26 00:30:10
    1528 t 1.9999E-05 dt 4.87E-09 write d3plot file            03/05/26 00:30:11
    2504 t 2.3998E-05 dt 3.41E-09 write d3plot file            03/05/26 00:30:11
    3765 t 2.7999E-05 dt 2.92E-09 write d3plot file            03/05/26 00:30:11
    5000 t 3.1350E-05 dt 2.55E-09 flush i/o buffers            03/05/26 00:30:11
    5259 t 3.2000E-05 dt 2.50E-09 write d3plot file            03/05/26 00:30:12
    7020 t 3.5999E-05 dt 2.12E-09 write d3plot file            03/05/26 00:30:12
    8909 t 3.9998E-05 dt 2.11E-09 write d3plot file            03/05/26 00:30:12
   10000 t 4.2229E-05 dt 2.06E-09 flush i/o buffers            03/05/26 00:30:13
   10865 t 4.3998E-05 dt 2.03E-09 write d3plot file            03/05/26 00:30:13
   12851 t 4.8000E-05 dt 2.01E-09 write d3plot file            03/05/26 00:30:13
   14898 t 5.1999E-05 dt 2.02E-09 write d3plot file            03/05/26 00:30:14
   15000 t 5.2205E-05 dt 2.02E-09 flush i/o buffers            03/05/26 00:30:14
   16946 t 5.5999E-05 dt 2.02E-09 write d3plot file            03/05/26 00:30:14
   19018 t 6.0000E-05 dt 2.02E-09 write d3plot file            03/05/26 00:30:15
   20000 t 6.1904E-05 dt 1.90E-09 flush i/o buffers            03/05/26 00:30:15
   21101 t 6.3999E-05 dt 1.90E-09 write d3plot file            03/05/26 00:30:15
   23112 t 6.7999E-05 dt 2.02E-09 write d3plot file            03/05/26 00:30:16
   25000 t 7.1737E-05 dt 1.90E-09 flush i/o buffers            03/05/26 00:30:16
   25130 t 7.1999E-05 dt 2.02E-09 write d3plot file            03/05/26 00:30:16
   27207 t 7.5999E-05 dt 2.02E-09 write d3plot file            03/05/26 00:30:17
   29235 t 8.0000E-05 dt 2.03E-09 write d3plot file            03/05/26 00:30:17

 *** termination time reached ***
   29235 t 8.0002E-05 dt 2.03E-09 write d3dump01 file          03/05/26 00:30:17
   29235 t 8.0002E-05 dt 2.03E-09 write d3plot file            03/05/26 00:30:17

 N o r m a l    t e r m i n a t i o n                          03/05/26 00:30:17

 Memory required to complete solution   :     196K
 Additional dynamically allocated memory:    2689K
                                   Total:    2884K

 T i m i n g   i n f o r m a t i o n
                        CPU(seconds)   %CPU  Clock(seconds) %Clock
  ----------------------------------------------------------------
  Keyword Processing ... 1.4835E-02    0.20     1.4846E-02    0.12
    KW Reading ......... 5.8983E-03    0.08     5.8990E-03    0.05
    KW Writing ......... 2.7242E-03    0.04     2.7240E-03    0.02
  Initialization ....... 1.1677E-01    1.57     5.5512E+00   43.02
    Init Proc Phase 1 .. 5.5072E-02    0.74     5.6267E-02    0.44
    Init Proc Phase 2 .. 1.0280E-03    0.01     1.6690E-03    0.01
  Element processing ... 5.3812E+00   72.27     5.3813E+00   41.71
    Solids ............. 5.1149E+00   68.70     5.1085E+00   39.59
    Shells ............. 7.5719E-02    1.02     8.2468E-02    0.64
    ISO Shells ......... 2.8215E-02    0.38     2.8305E-02    0.22
    E Other ............ 1.2058E-01    1.62     1.1984E-01    0.93
  Binary databases ..... 3.8845E-02    0.52     3.9037E-02    0.30
  ASCII database ....... 4.7116E-01    6.33     4.7123E-01    3.65
  Contact algorithm .... 6.7530E-02    0.91     6.8124E-02    0.53
  Rigid Bodies ......... 1.2953E-01    1.74     1.2980E-01    1.01
    Mass Scaling ....... 3.2062E-02    0.43     3.2182E-02    0.25
    Update nodes ....... 2.9779E-02    0.40     2.9931E-02    0.23
  Time step size ....... 3.5516E-02    0.48     3.5525E-02    0.28
  Rigid wall ........... 1.0609E-01    1.42     1.0401E-01    0.81
  Group force file ..... 3.0329E-02    0.41     3.0649E-02    0.24
  Others ............... 7.5878E-02    1.02     7.5412E-02    0.58
  Misc. 1 .............. 4.7778E-01    6.42     5.0146E-01    3.89
    Scale Masses ....... 2.9569E-02    0.40     2.9591E-02    0.23
    Force Constraints .. 2.9151E-02    0.39     2.9147E-02    0.23
    Force to Accel ..... 1.1395E-01    1.53     1.1376E-01    0.88
  Misc. 2 .............. 2.6893E-01    3.61     2.6869E-01    2.08
    Geometry Update .... 1.6873E-01    2.27     1.6675E-01    1.29
  Misc. 3 .............. 4.7689E-02    0.64     4.7357E-02    0.37
  Misc. 4 .............. 1.8338E-01    2.46     1.8359E-01    1.42
    Timestep Init ...... 7.6430E-02    1.03     7.6028E-02    0.59
    Apply Loads ........ 4.5229E-02    0.61     4.5058E-02    0.35
  ----------------------------------------------------------------
  T o t a l s            7.4455E+00  100.00     1.2902E+01  100.00

 Problem time       =    8.0002E-05
 Problem cycle      =     29235
 Total CPU time     =         7 seconds (   0 hours  0 minutes  7 seconds)
 CPU time per zone cycle  =        325.326 nanoseconds
 Clock time per zone cycle=        326.321 nanoseconds

 Number of CPU's    1
 NLQ used/max     192/  192
 Start time   03/05/2026 00:30:10
 End time     03/05/2026 00:30:17
 Elapsed time       7 seconds for   29235 cycles using  1 SMP thread
             (      0 hour   0 minute   7 seconds)

 N o r m a l    t e r m i n a t i o n                          03/05/26 00:30:17

Text(0, 0.5, 'Energy (mJ)')

Generate graphical output#

etc etc

plot taylor bar

Total running time of the script: (2 minutes 10.619 seconds)

Gallery generated by Sphinx-Gallery