North American LS-DYNA User Forum 2023

More Details

Call for Papers: Join us at the 14th European LS-DYNA Conference in Germany!

More Details

Ansys LS-DYNA 2023 R1 Update Presentation – What’s New!

More Details

Community

Scripting Command Language LsPrePost – Python

Hi,

I would like to create some parametric models in LsPrePost using Python. I have seen some examples on https://ftp.lstc.com/anonymous/outgoing/lsprepost/SCLexamples/lsppscripting.pdf 

However, I cannot find any documentation for the modules (DataCenter, LsPrePost), help, and how to start creating my models.

If you could please help me with finding more resources to achieve my goal. I'd appreciate.

Thank you.

Avatar
Stefano Del Rosso October 1, 2023
Drilling SPH particles

Hi all.

I also simulate rock drilling using SPH particles.

I use the increased density to speed up the solution to the problem. Increasing the density for SPH is not recommended, but I did such tests, I also do tests without increasing the density, they have not yet completed. But I have questions.

The task is similar to LSTC

https://www.lstc-cmmg.org/drilling

Only I use the rock and my feed and drilling speeds. I use an average speed of 25 m/s feed speed and 27 rad/s drilling speed. While solving the problem, I encountered several problems. 

1. Particles scatter wildly. And when the particles fly apart, they pull the solid ring behind them and it is greatly distorted:

2.I reduced the density a little, the solution became better and I got good stress. I also replaced the EROSION contact with AUTOMATIC:

But the particles still fly away strongly and do not break off in large pieces. Why? how to deal with this? It happens like this: the drill pierces the particles, they begin to fly away more and more until the end of the calculation time.

I use concrete material JHC

Avatar
javat33489@gmail.com September 28, 2023
Drilling SPG particles

Hi all.

I'm simulating rock drilling using SPG particles.

I use the increased density to speed up the solution to the problem.

The task is similar to LSTC

https://www.lstc-cmmg.org/drilling

Only I use the rock and my feed and drilling speeds.

I use an average speed of 25 m/s feed speed and 27 rad/s drilling speed.

While solving the problem, I encountered several problems.

1.I get a negative volume error. And when the particles fly apart, they pull the solid ring behind them and it is greatly distorted:

At low feed speeds there are no such problems.

2.I reduced the density a little, the solution became better and I got good stress. I also replaced the EROSION contact with AUTOMATIC:

But now I get warnings:

Why? how to deal with this?

Avatar
September 28, 2023
AUTODYN Subroutine – Cycle limit

Dear all,

I am using an Autodyn subroutine for a project and the simulation stops at cycle 99,999, see the image below.

What I do not understand is that the criteria that could limit the number of cycles is much more important than 100,000:
- Controls > Wrapup Criteria: >100,000;
- Output > Interrupt: End cycle = 100,000,000;
- Output > Call EXEDIT: End cycle = 100,000,000;
In addition, the subroutine usrsub2.f90 calls the "WRAPUP" module but does not redefine it, so there is no arbitrary "NCYLIM" defined in the subroutine. Do you know where this limitation might come from?

Thank you in advance for your help,

Audrey

Avatar
Audrey PICCINI September 28, 2023
Make shear waves (S-waves) in cube

Hi!

I was trying to make shear waves in Ansys Explicit Dynamics. I am not sure what to do, and wants to listen to all suggestions. My initial plan was to make a box (5x5x10mm), and do something with (one/ some of) the nodes to make a shear wave propagating throughout the "box" by time.

This is my box so far:


Any suggestions?

Avatar
Magnus Norstenes September 27, 2023
problem about *part_ADPTIVE_FAILURE

i was trying to solve example of SPR which is in LS-Dyna example site.

https://www.dynaexamples.com/introduction/process_simulation/self-piercing-riveting

 

my status doesnt go further but in LS-run status is still running (cpu usage is amlost none at this moment)

message is end just like this 

'Cutting plane=      3053      3056   0.80254E-01'

 i cant solve whole simulation because of this problem

how can i solve this example????

 

Avatar
Yun Seung Min September 27, 2023
Stress data @ point

I want to see the stress variation over time at a specific point in LS dyna. Its a 3d model and ive attempted to use a stress probe with no luck so far. Any help would be appriecated! 

Avatar
Bill Tapp September 25, 2023
How to output temperature change of all the elements/nodes?

Hello,

I'm using SPH to simulate a high-velocity impact problem. 

The material type I used is MAT 15, whose history variable 5 is the temperature change. Now I would like to output all the temperature changes for all the SPH nodes in one part and it will take a very long time to read the result file and output the plot due to a large number of SPH nodes used. So, I would like to know whether it is possible to output all the element temperatures of the part when setting the keyword DATABASE. I've checked the manual but I cannot find which keyword to use. 

For instance, the nodal information can be output as nodout file and the energy information can be output as matsum.

Avatar
Rong Jin September 25, 2023
How to convert LS-DYNA reslut files(d3plot, binout) to ANSYS result file(rst)?

I used Lsdyna for vibration analysis and imported the response file into Virtual .Lab for noise analysis, but Virtual .Lab only supports rst files. So I need to convert the LSDYNA result file to rst file

Avatar
Cheng Xihao September 25, 2023
WB LS Dyna PreStress Bird Strike

Hi everyone,

I am trying to make a prestress bird strike analysis. In prestress (Static Structural) i have thermal, pressure and bolt preload loads. I am using SPH particle method to model bird, however, static structural doesn't support that. I couldn't find any way to solve it yet. Do you know how can constrain the SPH bird in static structural module? Or how can map the stresses which comes from WB Static structural to LS Dyna Explicit model ?

Regards.

Avatar
MEHMET AKKUŞ September 24, 2023
Model and Stress distribution from Ansys to LS Dyna

I am creating an additive manufactured Model with Iherent strain, LPBF function in Ansys and have created the model to study the residual stress. 

How can I import the model in LS Dyna for further research?

 

Thanks!

Avatar
Salem Obaid September 22, 2023
Blanks in cross sections using the section plane in LSPP

Hi,

When I observed the information of a cross-section using the section plane in LSPP, sometimes, some blank spots were observed. it seems that some elements at the selected cross-section are gone. I set the scale of displacement as zero.  What are the potential causes?  Thank you!

Kind regards,

Chayir

Avatar
chayir September 22, 2023
Point data

I have a thesis project using LS-DYNA 2023. I am attempting to use a single point of reference on a 3d solid plate and see the stress of it over time (cannot select surface). I have tried using named selections and coordinate system to highlight the point i want with no success. I have done similar things with fluent before but cannot seem to make this work. Sorry if my explaining is bad

 

Avatar
Bill Tapp September 22, 2023
The body does not rotate in LS-Dyna

Hi all.

A very strange situation. Previously I used a rigid body and set the translation and rotation using RIGID_BODY_MOTION. I made the body deformable and now set the rotation and translation using MOTION_SET. But rotation does not occur. The calculation gives an error, despite the fact that I have specified nodes for rotation

Avatar
javat33489@gmail.com September 21, 2023
Save Nodal Forces

I am trying to save the nodal forces of a DYNA simulation and can't get the NODFOR file to generate. My card deck is below. Does anyone know what is going on here? Thanks!

$ LS-DYNA Keyword file created by fem.dyna Python functions
$ Created on 05:53PM on July 24, 2023
*KEYWORD
*CONTROL_ENERGY
$ hgen, rwen, slnten, rylen
2,2,2,2
*CONTROL_MPP_IO_NOFULL
*CONTROL_TERMINATION
$ endtim, endcyc, dtmin, endeng, endmas
0.0045,0,0.0,0.0,0.0
*CONTROL_TIMESTEP
$ dtinit, tssfac, isdo, tslimt, dt2ms, lctm, erode, ms1st
1e-06,0.02,0,0.0,0.0,1,0,0
*DATABASE_NODOUT
$ dt, binary, lcur, ioopt, option1, option2
2e-05,0,0,1,0.0,0
*DATABASE_NODFOR
$ dt, binary, lcur, ioopt, option1, option2
1e-06,0,0,1,0.0,0
*DATABASE_EXTENT_BINARY
$ neiph, neips, maxint, strflg, sigflg, epsflg, rltflg, engflg
0,0,3,0,2,2,2,1
$ cmpflg, ieverp, beamip, dcomp, shge, stssz, n3thdt, ialemat
0,0,0,4,1,1,2,0
*DATABASE_FORMAT
$ iform, ibinary
0,1
*DATABASE_HISTORY_NODE_SET
$ id1, id2, id3, id4, id5, id6, id7, id8
1,0,0,0,0,0,0,0
*SET_NODE_GENERAL
1
ALL
*PART
$ Title
background
$ pid, secid, mid, eosid, hgid, grav, adpopt, tmid
1,1,1,0,0,0,0,0
*SECTION_SOLID
$ secid, elform, aet
1,2,0
*PART
$ Title
background pml
$ pid, secid, mid, eosid, hgid, grav, adpopt, tmid
2,2,2,0,0,0,0,0
*SECTION_SOLID
$ secid, elform, aet
2,2,0
*INCLUDE
../../nodes.dyn
*INCLUDE
../../elems.dyn
*INCLUDE
../NodalLoads.dyn
*INCLUDE
Materials.dyn
*INCLUDE
LoadCurves.dyn
*END

Avatar
Joey Richardson September 21, 2023
Shell Element with large distortions

Hello,

I am simulating the deformation of a metallic tube in the low velocity range.
Unfortunately, large distortions occur at the shell elements (see figure). 
These stretch to infinity and the voltage drops to 0.
These distortions occur during self-contact.
Attempts to counteract this by applying Hourglass were not really successful.

A short excerpt from the K-File:

*CONTROL_HOURGLASS
$#     ihq        qh  
         8      0.03
*CONTROL_SHELL
$#  wrpang     esort     irnxx    istupd    theory       bwc     miter      proj
      20.0         1        -1         0        16         2         1         0
$# rotascl    intgrd    lamsht    cstyp6    thshel      
       1.0         0         0         1         0
*CONTACT_AUTOMATIC_SINGLE_SURFACE_ID
$#     cid                                                                 title
         1Self_Contact
$#   surfa     surfb  surfatyp  surfbtyp   saboxid   sbboxid      sapr      sbpr
        13         0         3         0         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
       0.2      0.18       0.0       0.0       5.0         0       0.01.00000E20
$#    sfsa      sfsb      sast      sbst     sfsat     sfsbt       fsf       vsf
       2.0       1.0       0.0       0.0       1.0       1.0       1.0       1.0
$#    soft    sofscl    lcidab    maxpar     sbopt     depth     bsort    frcfrq
         1       0.2         0     1.025       3.0         2       100         1
$#  penmax    thkopt    shlthk     snlog      isym     i2d3d    sldthk    sldstf
       0.0         1         2         0         0         0       0.0       0.0
$#    igap    ignore    dprfac    dtstif    edgek     unused    flangl   cid_rcf
         1         0       0.0       0.0       0.0                 0.0         0
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_ID
$#     cid                                                                 title
         2dropmass
$#   surfa     surfb  surfatyp  surfbtyp   saboxid   sbboxid      sapr      sbpr
        13        10         3         3         0         0         2         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
       0.2      0.18       0.0       0.0                   0       0.01.00000E20
$#    sfsa      sfsb      sast      sbst     sfsat     sfsbt       fsf       vsf
       1.0       1.0                 0.0       1.0       1.0       1.0       1.0
$#    soft    sofscl    lcidab    maxpar     sbopt     depth     bsort    frcfrq
         2       0.1         0     1.025       1.0         2        10         1
$#  penmax    thkopt    shlthk     snlog      isym     i2d3d    sldthk    sldstf
       0.0         0         0         0         0         0       0.0       0.0
$#    igap    ignore    dprfac    dtstif    edgek     unused    flangl   cid_rcf
         1         0       0.0       0.0       0.0                 0.0         0
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_ID
$#     cid                                                                 title
         3bottom
$#   surfa     surfb  surfatyp  surfbtyp   saboxid   sbboxid      sapr      sbpr
        13        11         3         3         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
       0.2      0.18       0.0       0.0                   0       0.01.00000E20
$#    sfsa      sfsb      sast      sbst     sfsat     sfsbt       fsf       vsf
       1.0       1.0                 0.0       1.0       1.0       1.0       1.0
$#    soft    sofscl    lcidab    maxpar     sbopt     depth     bsort    frcfrq
         2       0.1         0     1.025       1.0         2        10         1
$#  penmax    thkopt    shlthk     snlog      isym     i2d3d    sldthk    sldstf
       0.0         0         0         0         0         0       0.0       0.0
$#    igap    ignore    dprfac    dtstif    edgek     unused    flangl   cid_rcf
         1         0       0.0       0.0       0.0                 0.0         0

 

Maybe someone has an idea why this happens?

Avatar
Noah September 21, 2023
Keep Track of Number of Eroded Elements

I'd like to plot Number of Eroded Elements vs Time in LSPP. Is this possible?

Avatar
Jacob Ellis September 20, 2023
About the d3plot files not fully generated

Hello,

There is an issue of the Ls-dyna SPH simulation. When I use the solver on my University HPC to do a SPH simulation, the d3plot cannot fully generated. The total set steps are about 102 d3plots for normal termination, while after generating about 80 d3plots, the solver got stuck and the task reported an error after running out of time. However, when I use the LS-Run on my Windows desktop, the job could end normally. The solver version on the HPC is R 10.1 and the local solver on my desktop is R 13.0. I’m not sure whether the solver version led to that issue.

The message file for the failed job:

 Date: 09/15/2023      Time: 15:40:54  

     ___________________________________________________
     |                                                 |
     |  Livermore  Software  Technology  Corporation   |
     |                                                 |
     |  7374 Las Positas Road                          |
     |  Livermore, CA 94551                            |
     |  Tel: (925) 449-2500  Fax: (925) 449-2507       |
     |  www.lstc.com                                   |
     |_________________________________________________|
     |                                                 |
     |  LS-DYNA, A Program for Nonlinear Dynamic       |
     |  Analysis of Structures in Three Dimensions     |
     |  Version : smp s R10.1.0   Date: 01/12/2018     |
     |  Revision: 123264          Time: 16:34:46       |
     |                                                 |
     |  Features enabled in this version:              |
     |    Shared Memory Parallel                       |
     |    ANSYS Database format                        |
     |    ANSYS License (ANSYS192)                     |
     |                                                 |
     |  Licensed to:                                   |
     |                                                 |
     |  Platform   : Xeon64 System                     |
     |  OS Level   : Linux 2.6.18 uum                  |
     |  Compiler   : Intel Fortran Compiler 16.0 SSE2  |
     |  Hostname   : cascade051                        |
     |  Precision  : Single precision (I4R4)           |
     |  SVN Version: 123300                            |
     |                                                 |
     |  Unauthorized use infringes LSTC copyrights     |
     |_________________________________________________|


 Input file: Run.k                                                                          

 The native file format       : 32-bit small endian
 Memory size from command line:  1000000000

 ********  notice  ********  notice  ********  notice  ********
 *                                                            *
 * This is the LS-DYNA Finite Element code.                   *
 *                                                            *
 * Neither LSTC 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.                                                   *
 *                                                            *
 * LSTC 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 LSTC through your site focus.               *
 *                                                            *
 * Copyright (C) 1990-2018                                    *
 * by Livermore Software Technology Corp. (LSTC)              *
 * All rights reserved                                        *
 *                                                            *
 ********  notice  ********  notice  ********  notice  ********
 
LSTC Environment Variables:
 
LSTC_LICENSE=ansys
 

 Beginning of keyword reader                                   09/15/23 15:40:54
 expanding   memory to  1000000000 d         600 beginning of keyword reader
 contracting memory to       21002 d         600  
 expanding   memory to       81130 d     2212719  
 expanding   memory to      111130 d     2212719  
 expanding   memory to      171130 d     2212719  
 expanding   memory to      191576 d     2212719  
 expanding   memory to      212264 d     2212719  
 expanding   memory to      233264 d     2212719  

 Open include file: SPH_Ball_15x15x15.k
 expanding   memory to      255440 d     2212719  

 Open include file: SPH_Plate_90x90x30.k
 expanding   memory to      277440 d     2212719  
 expanding   memory to      299040 d     2212719  
 expanding   memory to      320640 d     2212719  
 expanding   memory to      342240 d     2212719  
 expanding   memory to      363840 d     2212719  
 expanding   memory to      385440 d     2212719  
 expanding   memory to      407040 d     2212719  
 expanding   memory to      428640 d     2212719  
 expanding   memory to      450240 d     2212719  
 expanding   memory to      471840 d     2212719  
 expanding   memory to      493440 d     2212719  
 expanding   memory to      515040 d     2212719  
 expanding   memory to      536640 d     2212719  
 expanding   memory to      558240 d     2212719  
 expanding   memory to      579840 d     2212719  
 expanding   memory to      601440 d     2212719  
 expanding   memory to      623040 d     2212719  
 expanding   memory to      644640 d     2212719  
 expanding   memory to      666240 d     2212719  
 expanding   memory to      687840 d     2212719  
 expanding   memory to      709440 d     2212719  
 expanding   memory to      731040 d     2212719  
 expanding   memory to      752640 d     2212719  
 expanding   memory to      774240 d     2212719  
 expanding   memory to      795840 d     2212719  
 expanding   memory to      817440 d     2212719  
 expanding   memory to      839040 d     2212719  
 expanding   memory to      860640 d     2212719  
 expanding   memory to      882240 d     2212719  
 expanding   memory to      903840 d     2212719  
 expanding   memory to      925440 d     2212719  
 expanding   memory to      947040 d     2212719  
 expanding   memory to      968640 d     2212719  
 expanding   memory to      990240 d     2212719  
 expanding   memory to     1011840 d     2212719  
 expanding   memory to     1033440 d     2212719  
 expanding   memory to     1055040 d     2212719  
 expanding   memory to     1076640 d     2212719  
 expanding   memory to     1098240 d     2212719  
 expanding   memory to     1119840 d     2212719  
 expanding   memory to     1141440 d     2212719  
 expanding   memory to     1163040 d     2212719  
 expanding   memory to     1184640 d     2212719  
 expanding   memory to     1206240 d     2212719  
 expanding   memory to     1227840 d     2212719  
 expanding   memory to     1249440 d     2212719  
 expanding   memory to     1271040 d     2212719  
 expanding   memory to     1292640 d     2212719  
 expanding   memory to     1314240 d     2212719  
 expanding   memory to     1335840 d     2212719  
 expanding   memory to     1357440 d     2212719  
 expanding   memory to     1379040 d     2212719  
 expanding   memory to     1400640 d     2212719  
 expanding   memory to     1422240 d     2212719  
 expanding   memory to     1443840 d     2212719  
 expanding   memory to     1465440 d     2212719  
 expanding   memory to     1487040 d     2212719  
 expanding   memory to     1508640 d     2212719  
 expanding   memory to     1530240 d     2212719  
 expanding   memory to     1551840 d     2212719  
 expanding   memory to     1573440 d     2212719  
 expanding   memory to     1595040 d     2212719  
 expanding   memory to     1616640 d     2212719  
 expanding   memory to     1638240 d     2212719  
 expanding   memory to     1659840 d     2212719  
 expanding   memory to     1681440 d     2212719  
 expanding   memory to     1703040 d     2212719  

 Open include file: MatFile_Y5_PA.k
 expanding   memory to     1732040 d     2212719  
 expanding   memory to     1815114 d     2212719  
 expanding   memory to     1896712 d     2212719  
 expanding   memory to     1978310 d     2212719  
 expanding   memory to     2059908 d     2212719  
 expanding   memory to     2141506 d     2212719  
 expanding   memory to     2223104 d     2212719  
 expanding   memory to     2639248 d     2212719  
 expanding   memory to     2884040 d     2212719  
 expanding   memory to     3128832 d     2212719  
 expanding   memory to     5821534 d     2212719  
 expanding   memory to     6066326 d     2212719  
 expanding   memory to     6311118 d     2212719  
 expanding   memory to     6737848 d     2212725  
 expanding   memory to     7153992 d     2212725  

 Memory required to process keyword     :     7153992
 contracting memory to           1 d     2212725 end of keyword reader

 LS-DYNA will perform a structural only analysis

 Y5_PA_HV                                                                
                         ls-dyna smp.123264 s           date 01/12/2018
 expanding   memory to        2767 d         666 material, node and element data
 expanding   memory to        2845 d         668 mat add general damage 0
 expanding   memory to        2865 d         668 material added airbag porosity leakage
 expanding   memory to        2871 d         668 material properties
 expanding   memory to        2968 d      246487 element data for composite material
 expanding   memory to        2982 d     8323567 gradient of damage data
 expanding   memory to        3057 d     8323569 user control for contact interfaces
 contracting memory to        3039 d     8323569 thermal expansion 0
 expanding   memory to      247830 d     8323569 temp index array for node sorting
 contracting memory to        3039 d     8568364 for boundary cond. vectorization
 contracting memory to        3039 d     9057946 arbitrary element numbering input
 contracting memory to        3039 d     9057946 thick shell int. pt. strain
 expanding   memory to        3597 d    10281906 process element material data
 expanding   memory to        5598 d    10281904 load curve data
 contracting memory to        3604 d    10281904 user_define_loading_set
 expanding   memory to        3804 d    10281904 regularize time history data
 expanding   memory to        3808 d    10281904 *DEFINE_ELEMENT_DEATH
 contracting memory to        3808 d    10281904 rigidwall definitions
 expanding   memory to      248600 d    10281904 nodal constraint sets
 contracting memory to        3998 d    10281904 contact interfaces
 expanding   memory to      494580 d    10281904 implicit contact 1
 contracting memory to        4998 d    10282008 contact segment data 1
 expanding   memory to     1228953 d    10282008 contact segment data 3
 contracting memory to        3998 d    14493971 contact definition
 expanding   memory to      983386 d    14493971 rigid body
 contracting memory to        4638 d    14493995 contact friction
 expanding   memory to        4669 d    14493995 joints
 expanding   memory to    35744235 d    14494595 SPH elements
 expanding   memory to    36233817 d    14494595 discrete element
 contracting memory to    36233817 d    14739710 accelerometers
 expanding   memory to    36233917 d    14739710 rigid bodies for disp. termination
 contracting memory to    36233917 d    15474083 transformation to move nodes
 expanding   memory to    36233919 d    14739761 mat add cohesive 0
 contracting memory to    36233919 d    14739761 sensor_control for curve
 expanding   memory to    36233921 d    14739761 flexible body 2
 contracting memory to    36233921 d    14739761 implicit joints

 S t o r a g e   a l l o c a t i o n  

 Memory required to begin solution      :      37M
 Additional dynamically allocated memory:      15M
                                   Total:      51M
       1 t 0.0000E+00 dt 1.10E-12 flush i/o buffers            09/15/23 15:41:07
       1 t 0.0000E+00 dt 1.10E-12 write d3plot file            09/15/23 15:41:07
      96 t 9.4113E-08 dt 9.41E-09 write d3plot file            09/15/23 15:41:58
     103 t 1.8341E-07 dt 1.83E-08 write d3plot file            09/15/23 15:42:01
     108 t 2.9539E-07 dt 2.95E-08 write d3plot file            09/15/23 15:42:04
     111 t 3.8597E-07 dt 3.01E-08 write d3plot file            09/15/23 15:42:06
     114 t 4.8081E-07 dt 3.33E-08 write d3plot file            09/15/23 15:42:08
     117 t 5.8049E-07 dt 3.26E-08 write d3plot file            09/15/23 15:42:09
     120 t 6.7355E-07 dt 3.07E-08 write d3plot file            09/15/23 15:42:11
     123 t 7.7186E-07 dt 3.20E-08 write d3plot file            09/15/23 15:42:13
     126 t 8.6933E-07 dt 3.56E-08 write d3plot file            09/15/23 15:42:15
     129 t 9.7590E-07 dt 3.48E-08 write d3plot file            09/15/23 15:42:17
     132 t 1.0795E-06 dt 3.37E-08 write d3plot file            09/15/23 15:42:18
     135 t 1.1817E-06 dt 3.54E-08 write d3plot file            09/15/23 15:42:20
     138 t 1.2895E-06 dt 3.64E-08 write d3plot file            09/15/23 15:42:22
     141 t 1.3997E-06 dt 3.62E-08 write d3plot file            09/15/23 15:42:24
     143 t 1.4717E-06 dt 3.61E-08 write d3plot file            09/15/23 15:42:25
     146 t 1.5796E-06 dt 3.49E-08 write d3plot file            09/15/23 15:42:27
     149 t 1.6842E-06 dt 3.46E-08 write d3plot file            09/15/23 15:42:29
     152 t 1.7902E-06 dt 3.61E-08 write d3plot file            09/15/23 15:42:31
 
 *** Warning 40042 (SOL+42)
     Part ID = 5000001 plasticity algorithm did not converge for MAT15
     time = 1.8989E-06, ncycle = 154
     155 t 1.8989E-06 dt 3.67E-08 write d3plot file            09/15/23 15:42:33
 
 *** Warning 40042 (SOL+42)
     Part ID = 5000001 plasticity algorithm did not converge for MAT15
     time = 1.9356E-06, ncycle = 155
     157 t 1.9725E-06 dt 3.76E-08 write d3plot file            09/15/23 15:42:34
     160 t 2.0818E-06 dt 3.59E-08 write d3plot file            09/15/23 15:42:36
     163 t 2.1906E-06 dt 3.85E-08 write d3plot file            09/15/23 15:42:38
     165 t 2.2671E-06 dt 3.78E-08 write d3plot file            09/15/23 15:42:40
     168 t 2.3827E-06 dt 3.76E-08 write d3plot file            09/15/23 15:42:41
     171 t 2.4913E-06 dt 3.90E-08 write d3plot file            09/15/23 15:42:43
     173 t 2.5694E-06 dt 3.79E-08 write d3plot file            09/15/23 15:42:45
     176 t 2.6795E-06 dt 3.31E-08 write d3plot file            09/15/23 15:42:47
     179 t 2.7811E-06 dt 3.62E-08 write d3plot file            09/15/23 15:42:49
     182 t 2.8865E-06 dt 3.85E-08 write d3plot file            09/15/23 15:42:51
     184 t 2.9654E-06 dt 3.91E-08 write d3plot file            09/15/23 15:42:52
     187 t 3.0777E-06 dt 3.68E-08 write d3plot file            09/15/23 15:42:55
     190 t 3.1960E-06 dt 4.16E-08 write d3plot file            09/15/23 15:42:57
 
 *** Warning 40042 (SOL+42)
     Part ID = 5000001 plasticity algorithm did not converge for MAT15
     time = 3.2376E-06, ncycle = 190
     192 t 3.2788E-06 dt 3.88E-08 write d3plot file            09/15/23 15:42:58
     195 t 3.3979E-06 dt 3.85E-08 write d3plot file            09/15/23 15:43:00
     197 t 3.4761E-06 dt 4.04E-08 write d3plot file            09/15/23 15:43:02
     200 t 3.5958E-06 dt 4.16E-08 write d3plot file            09/15/23 15:43:04
     202 t 3.6748E-06 dt 3.24E-08 write d3plot file            09/15/23 15:43:05
     205 t 3.7821E-06 dt 4.01E-08 write d3plot file            09/15/23 15:43:08
     208 t 3.8980E-06 dt 3.79E-08 write d3plot file            09/15/23 15:43:10
     210 t 3.9736E-06 dt 3.74E-08 write d3plot file            09/15/23 15:43:11
     213 t 4.0845E-06 dt 4.02E-08 write d3plot file            09/15/23 15:43:13
     215 t 4.1662E-06 dt 4.12E-08 write d3plot file            09/15/23 15:43:15
     218 t 4.2846E-06 dt 4.04E-08 write d3plot file            09/15/23 15:43:17
     220 t 4.3617E-06 dt 3.91E-08 write d3plot file            09/15/23 15:43:18
 
 *** Warning 40042 (SOL+42)
     Part ID = 5000001 plasticity algorithm did not converge for MAT15
     time = 4.4007E-06, ncycle = 220
     223 t 4.4836E-06 dt 3.87E-08 write d3plot file            09/15/23 15:43:21
     225 t 4.5636E-06 dt 4.20E-08 write d3plot file            09/15/23 15:43:22
     228 t 4.6917E-06 dt 4.26E-08 write d3plot file            09/15/23 15:43:24
     230 t 4.7746E-06 dt 4.05E-08 write d3plot file            09/15/23 15:43:26
     233 t 4.8991E-06 dt 4.24E-08 write d3plot file            09/15/23 15:43:28
     235 t 4.9825E-06 dt 3.99E-08 write d3plot file            09/15/23 15:43:29
     237 t 5.0618E-06 dt 3.95E-08 write d3plot file            09/15/23 15:43:31
     240 t 5.1818E-06 dt 4.02E-08 write d3plot file            09/15/23 15:43:33
     242 t 5.2621E-06 dt 4.03E-08 write d3plot file            09/15/23 15:43:35
     245 t 5.3836E-06 dt 4.04E-08 write d3plot file            09/15/23 15:43:37
     247 t 5.4663E-06 dt 4.22E-08 write d3plot file            09/15/23 15:43:39
     250 t 5.5928E-06 dt 4.26E-08 write d3plot file            09/15/23 15:43:41
     252 t 5.6780E-06 dt 4.16E-08 write d3plot file            09/15/23 15:43:42
     254 t 5.7615E-06 dt 4.08E-08 write d3plot file            09/15/23 15:43:44
     257 t 5.8819E-06 dt 4.13E-08 write d3plot file            09/15/23 15:43:46
     259 t 5.9664E-06 dt 4.18E-08 write d3plot file            09/15/23 15:43:48
     262 t 6.0951E-06 dt 4.30E-08 write d3plot file            09/15/23 15:43:50
     264 t 6.1808E-06 dt 4.08E-08 write d3plot file            09/15/23 15:43:52
     266 t 6.2609E-06 dt 4.32E-08 write d3plot file            09/15/23 15:43:53
     269 t 6.3886E-06 dt 4.23E-08 write d3plot file            09/15/23 15:43:56
     271 t 6.4683E-06 dt 3.99E-08 write d3plot file            09/15/23 15:43:57
     274 t 6.5920E-06 dt 4.11E-08 write d3plot file            09/15/23 15:43:59
     276 t 6.6769E-06 dt 4.32E-08 write d3plot file            09/15/23 15:44:01
     278 t 6.7617E-06 dt 4.21E-08 write d3plot file            09/15/23 15:44:02
     281 t 6.8906E-06 dt 4.39E-08 write d3plot file            09/15/23 15:44:05
     283 t 6.9769E-06 dt 4.41E-08 write d3plot file            09/15/23 15:44:06
     285 t 7.0638E-06 dt 4.34E-08 write d3plot file            09/15/23 15:44:08
     288 t 7.1906E-06 dt 4.37E-08 write d3plot file            09/15/23 15:44:10
     290 t 7.2781E-06 dt 4.24E-08 write d3plot file            09/15/23 15:44:11
     292 t 7.3642E-06 dt 4.12E-08 write d3plot file            09/15/23 15:44:13
     295 t 7.4925E-06 dt 4.40E-08 write d3plot file            09/15/23 15:44:15
     297 t 7.5784E-06 dt 4.34E-08 write d3plot file            09/15/23 15:44:17
     299 t 7.6630E-06 dt 4.21E-08 write d3plot file            09/15/23 15:44:18
     302 t 7.7922E-06 dt 4.34E-08 write d3plot file            09/15/23 15:44:20
     304 t 7.8801E-06 dt 4.19E-08 write d3plot file            09/15/23 15:44:22
     306 t 7.9652E-06 dt 4.42E-08 write d3plot file            09/15/23 15:44:23
     309 t 8.0970E-06 dt 4.45E-08 write d3plot file            09/15/23 15:44:26
     311 t 8.1842E-06 dt 4.52E-08 write d3plot file            09/15/23 15:44:27
     313 t 8.2732E-06 dt 4.44E-08 write d3plot file            09/15/23 15:44:29
     315 t 8.3618E-06 dt 4.39E-08 write d3plot file            09/15/23 15:44:30
     318 t 8.4897E-06 dt 1.07E-08 write d3plot file            09/15/23 15:44:33
 
 *** Warning 40042 (SOL+42)
     Part ID = 5000001 plasticity algorithm did not converge for MAT15
     time = 8.5089E-06, ncycle = 1061
 
The message file for the successful job:

 

 ***************************************************************
 *            ANSYS 2023 R1          LEGAL NOTICES             *
 ***************************************************************
 *                                                             *
 * Copyright 1971-2023 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: 09/19/2023      Time: 15:42:25  
      ___________________________________________________ 
     |                                                   |
     |  LS-DYNA, A Program for Nonlinear Dynamic         |
     |  Analysis of Structures in Three Dimensions       |
     |  Date    : 11/15/2022    Time: 13:48:14           |
     |  Version : smp s R13                              |
     |  Revision: R13.1.1-13-g5480eca294                 |
     |                                                   |
     |  Features enabled in this version:                |
     |    Shared Memory Parallel                         |
     |    FFTW (multi-dimensional FFTW Library)          |
     |    Interactive Graphics                           |
     |    ANSYSLIC enabled                               |
     |    MADYMO indirect coupling (7.2 or higher)       |
     |                                                   |
     |  Licensed to:                                     |
     |  Issued by  : ANSYS                               |
     |                                                   |
     |  Platform   : WINDOWS X64 (SSE2)                  |
     |  OS Level   : Windows 8/10 & Srv 2016/2019 uam    |
     |  Compiler   : Intel Fortran XE 2019 MSVC++ 2019   |
     |  Hostname   : ENGR6321                            |
     |  Precision  : Single precision (I4R4)             |
     |                                                   |
     |  Unauthorized use infringes Ansys Inc. copyrights |
     |___________________________________________________|


 Command line options: i=C:\Users\rji229\Documents\fem_sph_circular\Test\1\Run.k 
                       ncpu=1 
                       memory=97m

 Input file: C:\Users\rji229\Documents\fem_sph_circular\Test\1\Run.k                         

 The native file format       : 32-bit small endian
 Memory size from command line:    97000000
 allocating  memory to    97000000 d         544 Initial allocation

 ********  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-2022                                    *
 * by Ansys, Inc.                                             *
 * All rights reserved                                        *
 *                                                            *
 ********  notice  ********  notice  ********  notice  ********

 Beginning of keyword reader                                   09/19/23 15:42:25
 expanding   memory to       20720 d         600  
 expanding   memory to       66080 d     2212719  
 expanding   memory to       98580 d     2212719  
 expanding   memory to     1578326 d     2272719  
 expanding   memory to     1606836 d     2272719  
 expanding   memory to     1629036 d     2280719  
 expanding   memory to     1744036 d     2280719  

                                                               09/19/23 15:42:25
 Open include file: SPH_Ball_15x15x15.k

                                                               09/19/23 15:42:25
 Open include file: SPH_Plate_90x90x30.k

                                                               09/19/23 15:42:25
 Open include file: MatFile_Y5_PA.k
 expanding   memory to     1765932 d     2280719  
 expanding   memory to     2012200 d     5465002  
 expanding   memory to     2256992 d     5465002  

 Memory required to process keyword     :    2256992  
 Additional dynamic memory required     :    8647293  

 contracting memory to           1 d     5465008 end of keyword reader

 LS-DYNA will perform a structural only analysis 

 Y5_PA_HV                                                                
                    ls-dyna smp s R13 R13.1.1_13_5480ec date 11/15/2022
 expanding   memory to        2001 d        1445 material, node and element data
 expanding   memory to        2071 d        1475 material properties
 expanding   memory to        2085 d     8324474 gradient of damage data
 expanding   memory to        2107 d     8324529 user control for contact interfaces
 contracting memory to        2087 d     8324531 thermal expansion 0
 expanding   memory to        2089 d     8324533 cross-sectional force definitions 0
 expanding   memory to      246880 d     8324547 temp index array for node sorting
 contracting memory to        2089 d     8569342 just because?
 expanding   memory to        4091 d    10283422 load curve data
 contracting memory to        2097 d    10283422 user_define_loading_set
 expanding   memory to        2297 d    10283422 regularize time history data
 contracting memory to        2101 d    10283422 *DEFINE_ELEMENT_DEATH
 expanding   memory to      492683 d    10528238 implicit contact 1
 contracting memory to        3101 d    10528344 contact segment data 1
 expanding   memory to     1227056 d    10528344 contact segment data 3
 contracting memory to        2101 d    14740307 contact definition
 expanding   memory to        2105 d    14740978 contact friction
 expanding   memory to     1715702 d    37752065 SPH elements
 expanding   memory to     1715802 d    38241723 rigid bodies for disp. termination
 expanding   memory to     1715872 d    38241723 save d3plot control data

 S t o r a g e   a l l o c a t i o n   

 Memory required to begin solution      :    2257K
 Additional dynamically allocated memory:      41M
                                   Total:      43M
       1 t 0.0000E+00 dt 1.10E-12 flush i/o buffers            09/19/23 15:42:33
       1 t 0.0000E+00 dt 1.10E-12 write d3plot file            09/19/23 15:42:33
      96 t 9.4113E-08 dt 9.41E-09 write d3plot file            09/19/23 15:43:22
     103 t 1.8341E-07 dt 1.83E-08 write d3plot file            09/19/23 15:43:25
     108 t 2.9539E-07 dt 2.95E-08 write d3plot file            09/19/23 15:43:28
 
 *** Warning 40042 (SOL+42)
     Part ID = 5000001 plasticity algorithm did not converge for MAT15
     time = 3.9317E-07, ncycle = 110
     111 t 3.9317E-07 dt 3.66E-08 write d3plot file            09/19/23 15:43:30
     113 t 4.6671E-07 dt 3.96E-08 write d3plot file            09/19/23 15:43:31
     116 t 5.9188E-07 dt 3.95E-08 write d3plot file            09/19/23 15:43:33
 
 *** Warning 40042 (SOL+42)
     Part ID = 5000001 plasticity algorithm did not converge for MAT15
     time = 6.3133E-07, ncycle = 116
     118 t 6.7157E-07 dt 4.11E-08 write d3plot file            09/19/23 15:43:34
     121 t 7.9410E-07 dt 4.17E-08 write d3plot file            09/19/23 15:43:37
     123 t 8.8038E-07 dt 4.54E-08 write d3plot file            09/19/23 15:43:38
 
 *** Warning 40042 (SOL+42)
     Part ID = 5000001 plasticity algorithm did not converge for MAT15
     time = 9.2573E-07, ncycle = 123
     125 t 9.6982E-07 dt 4.65E-08 write d3plot file            09/19/23 15:43:40
     127 t 1.0610E-06 dt 4.36E-08 write d3plot file            09/19/23 15:43:42
     130 t 1.1918E-06 dt 4.53E-08 write d3plot file            09/19/23 15:43:44
     132 t 1.2849E-06 dt 4.77E-08 write d3plot file            09/19/23 15:43:46
     134 t 1.3799E-06 dt 4.74E-08 write d3plot file            09/19/23 15:43:48
     136 t 1.4746E-06 dt 4.63E-08 write d3plot file            09/19/23 15:43:50
     138 t 1.5663E-06 dt 4.43E-08 write d3plot file            09/19/23 15:43:51
 
 *** Warning 40042 (SOL+42)
     Part ID = 5000001 plasticity algorithm did not converge for MAT15
     time = 1.6532E-06, ncycle = 139
     140 t 1.6532E-06 dt 4.69E-08 write d3plot file            09/19/23 15:43:53
     143 t 1.7969E-06 dt 4.83E-08 write d3plot file            09/19/23 15:43:56
     145 t 1.8806E-06 dt 3.90E-08 write d3plot file            09/19/23 15:43:57
     147 t 1.9626E-06 dt 4.46E-08 write d3plot file            09/19/23 15:43:59
     150 t 2.0788E-06 dt 3.60E-08 write d3plot file            09/19/23 15:44:02
     153 t 2.1979E-06 dt 4.44E-08 write d3plot file            09/19/23 15:44:05
     155 t 2.2817E-06 dt 3.75E-08 write d3plot file            09/19/23 15:44:06
     158 t 2.3922E-06 dt 3.80E-08 write d3plot file            09/19/23 15:44:09
     160 t 2.4720E-06 dt 4.60E-08 write d3plot file            09/19/23 15:44:11
 
 *** Warning 40042 (SOL+42)
     Part ID = 5000001 plasticity algorithm did not converge for MAT15
     time = 2.5625E-06, ncycle = 161
     162 t 2.5625E-06 dt 3.95E-08 write d3plot file            09/19/23 15:44:13
 
 *** Warning 40042 (SOL+42)
     Part ID = 5000001 plasticity algorithm did not converge for MAT15
     time = 2.6020E-06, ncycle = 162
 
 *** Warning 40042 (SOL+42)
     Part ID = 5000001 plasticity algorithm did not converge for MAT15
     time = 2.6867E-06, ncycle = 164
     165 t 2.6867E-06 dt 4.71E-08 write d3plot file            09/19/23 15:44:15
     167 t 2.7825E-06 dt 4.80E-08 write d3plot file            09/19/23 15:44:17
 
 *** Warning 40042 (SOL+42)
     Part ID = 5000001 plasticity algorithm did not converge for MAT15
     time = 2.8306E-06, ncycle = 167
     169 t 2.8726E-06 dt 3.94E-08 write d3plot file            09/19/23 15:44:19
     171 t 2.9541E-06 dt 4.63E-08 write d3plot file            09/19/23 15:44:21
     174 t 3.0939E-06 dt 4.62E-08 write d3plot file            09/19/23 15:44:24
     176 t 3.1808E-06 dt 4.38E-08 write d3plot file            09/19/23 15:44:26
     178 t 3.2728E-06 dt 4.89E-08 write d3plot file            09/19/23 15:44:28
     180 t 3.3658E-06 dt 4.24E-08 write d3plot file            09/19/23 15:44:30
     183 t 3.4967E-06 dt 4.81E-08 write d3plot file            09/19/23 15:44:33
 
 *** Warning 40042 (SOL+42)
     Part ID = 5000001 plasticity algorithm did not converge for MAT15
     time = 3.5926E-06, ncycle = 184
     185 t 3.5926E-06 dt 4.22E-08 write d3plot file            09/19/23 15:44:35
     187 t 3.6749E-06 dt 4.20E-08 write d3plot file            09/19/23 15:44:37
     189 t 3.7631E-06 dt 4.85E-08 write d3plot file            09/19/23 15:44:39
 
 *** Warning 40042 (SOL+42)
     Part ID = 5000001 plasticity algorithm did not converge for MAT15
     time = 3.8117E-06, ncycle = 189
     191 t 3.8592E-06 dt 4.25E-08 write d3plot file            09/19/23 15:44:41
     194 t 3.9923E-06 dt 4.82E-08 write d3plot file            09/19/23 15:44:43
     196 t 4.0878E-06 dt 4.72E-08 write d3plot file            09/19/23 15:44:45
     198 t 4.1826E-06 dt 4.38E-08 write d3plot file            09/19/23 15:44:47
     200 t 4.2747E-06 dt 4.86E-08 write d3plot file            09/19/23 15:44:49
     202 t 4.3722E-06 dt 4.90E-08 write d3plot file            09/19/23 15:44:51
     204 t 4.4688E-06 dt 4.47E-08 write d3plot file            09/19/23 15:44:53
     206 t 4.5568E-06 dt 4.54E-08 write d3plot file            09/19/23 15:44:55
     209 t 4.6991E-06 dt 4.69E-08 write d3plot file            09/19/23 15:44:58
     211 t 4.7946E-06 dt 4.44E-08 write d3plot file            09/19/23 15:45:00
     213 t 4.8792E-06 dt 4.17E-08 write d3plot file            09/19/23 15:45:01
     215 t 4.9669E-06 dt 4.79E-08 write d3plot file            09/19/23 15:45:04
     217 t 5.0629E-06 dt 4.76E-08 write d3plot file            09/19/23 15:45:06
     220 t 5.1926E-06 dt 4.49E-08 write d3plot file            09/19/23 15:45:09
     222 t 5.2864E-06 dt 4.41E-08 write d3plot file            09/19/23 15:45:11
     224 t 5.3770E-06 dt 4.62E-08 write d3plot file            09/19/23 15:45:13
     226 t 5.4674E-06 dt 4.83E-08 write d3plot file            09/19/23 15:45:15
     228 t 5.5637E-06 dt 4.82E-08 write d3plot file            09/19/23 15:45:17
     230 t 5.6554E-06 dt 4.69E-08 write d3plot file            09/19/23 15:45:19
     233 t 5.7966E-06 dt 4.41E-08 write d3plot file            09/19/23 15:45:22
     235 t 5.8892E-06 dt 4.64E-08 write d3plot file            09/19/23 15:45:24
     237 t 5.9818E-06 dt 4.82E-08 write d3plot file            09/19/23 15:45:25
     239 t 6.0788E-06 dt 4.88E-08 write d3plot file            09/19/23 15:45:27
     241 t 6.1725E-06 dt 4.63E-08 write d3plot file            09/19/23 15:45:28
     243 t 6.2657E-06 dt 4.86E-08 write d3plot file            09/19/23 15:45:30
     245 t 6.3630E-06 dt 4.79E-08 write d3plot file            09/19/23 15:45:31
     247 t 6.4585E-06 dt 4.50E-08 write d3plot file            09/19/23 15:45:33
     250 t 6.5981E-06 dt 4.74E-08 write d3plot file            09/19/23 15:45:35
     252 t 6.6944E-06 dt 4.91E-08 write d3plot file            09/19/23 15:45:36
     254 t 6.7913E-06 dt 4.90E-08 write d3plot file            09/19/23 15:45:38
     256 t 6.8895E-06 dt 4.78E-08 write d3plot file            09/19/23 15:45:40
     258 t 6.9860E-06 dt 4.74E-08 write d3plot file            09/19/23 15:45:41
     260 t 7.0814E-06 dt 4.89E-08 write d3plot file            09/19/23 15:45:43
     262 t 7.1792E-06 dt 4.88E-08 write d3plot file            09/19/23 15:45:44
     264 t 7.2768E-06 dt 4.78E-08 write d3plot file            09/19/23 15:45:46
     266 t 7.3730E-06 dt 4.87E-08 write d3plot file            09/19/23 15:45:47
     268 t 7.4703E-06 dt 4.86E-08 write d3plot file            09/19/23 15:45:49
     270 t 7.5675E-06 dt 4.86E-08 write d3plot file            09/19/23 15:45:50
     272 t 7.6649E-06 dt 4.66E-08 write d3plot file            09/19/23 15:45:52
     274 t 7.7580E-06 dt 4.66E-08 write d3plot file            09/19/23 15:45:53
     277 t 7.8993E-06 dt 4.65E-08 write d3plot file            09/19/23 15:45:56
     279 t 7.9947E-06 dt 4.85E-08 write d3plot file            09/19/23 15:45:57
     281 t 8.0890E-06 dt 4.70E-08 write d3plot file            09/19/23 15:45:59
 
 *** Warning 40042 (SOL+42)
     Part ID = 5000001 plasticity algorithm did not converge for MAT15
     time = 8.1844E-06, ncycle = 282
     283 t 8.1844E-06 dt 4.86E-08 write d3plot file            09/19/23 15:46:01
 
 *** Warning 40042 (SOL+42)
     Part ID = 5000001 plasticity algorithm did not converge for MAT15
     time = 8.2330E-06, ncycle = 283
     285 t 8.2819E-06 dt 4.88E-08 write d3plot file            09/19/23 15:46:02
     287 t 8.3787E-06 dt 4.88E-08 write d3plot file            09/19/23 15:46:04
     289 t 8.4710E-06 dt 4.78E-08 write d3plot file            09/19/23 15:46:05
     291 t 8.5661E-06 dt 4.87E-08 write d3plot file            09/19/23 15:46:07
     293 t 8.6636E-06 dt 4.89E-08 write d3plot file            09/19/23 15:46:08
     295 t 8.7594E-06 dt 4.79E-08 write d3plot file            09/19/23 15:46:10
     297 t 8.8559E-06 dt 4.85E-08 write d3plot file            09/19/23 15:46:11
     299 t 8.9530E-06 dt 4.84E-08 write d3plot file            09/19/23 15:46:13
     302 t 9.0990E-06 dt 4.87E-08 write d3plot file            09/19/23 15:46:15
     304 t 9.1965E-06 dt 4.56E-08 write d3plot file            09/19/23 15:46:17
     306 t 9.2881E-06 dt 4.69E-08 write d3plot file            09/19/23 15:46:18
     308 t 9.3827E-06 dt 4.90E-08 write d3plot file            09/19/23 15:46:20
     310 t 9.4806E-06 dt 4.82E-08 write d3plot file            09/19/23 15:46:21
     312 t 9.5729E-06 dt 4.67E-08 write d3plot file            09/19/23 15:46:23
     314 t 9.6638E-06 dt 4.66E-08 write d3plot file            09/19/23 15:46:24
     316 t 9.7579E-06 dt 4.85E-08 write d3plot file            09/19/23 15:46:26
 
 *** Warning 40042 (SOL+42)
     Part ID = 5000001 plasticity algorithm did not converge for MAT15
     time = 9.8512E-06, ncycle = 317
     319 t 9.8980E-06 dt 4.36E-08 write d3plot file            09/19/23 15:46:28
     321 t 9.9851E-06 dt 4.78E-08 write d3plot file            09/19/23 15:46:30
 
 *** termination time reached ***
     321 t 1.0033E-05 dt 4.78E-08 write d3dump01 file          09/19/23 15:46:30
     321 t 1.0033E-05 dt 4.78E-08 write d3plot file            09/19/23 15:46:30

 N o r m a l    t e r m i n a t i o n                          09/19/23 15:46:30

 S t o r a g e   a l l o c a t i o n   

 Memory required to complete solution   :    2257K
 Additional dynamically allocated memory:     150M
                                   Total:     152M

 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.0000E+00    0.41     4.9200E-01    0.20
    KW Reading ......... 1.0000E+00    0.41     3.3300E-01    0.14
    KW Writing ......... 0.0000E+00    0.00     3.1000E-02    0.01
  Initialization ....... 8.0000E+00    3.25     8.7360E+00    3.54
    Init Proc Phase 1 .. 7.0000E+00    2.85     7.0600E+00    2.86
    Init Proc Phase 2 .. 0.0000E+00    0.00     1.9000E-02    0.01
  Element processing ... 2.3200E+02   94.31     2.2922E+02   93.01
    SPH ................ 2.3200E+02   94.31     2.2922E+02   93.01
  Binary databases ..... 5.0000E+00    2.03     7.2970E+00    2.96
  ASCII database ....... 0.0000E+00    0.00     2.8000E-02    0.01
  Time step size ....... 0.0000E+00    0.00     1.5000E-02    0.01
  Others ............... 0.0000E+00    0.00     5.0000E-03    0.00
  Misc. 1 .............. 0.0000E+00    0.00     6.8000E-02    0.03
    Force to Accel ..... 0.0000E+00    0.00     6.6000E-02    0.03
  Misc. 2 .............. 0.0000E+00    0.00     3.7300E-01    0.15
  Misc. 3 .............. 0.0000E+00    0.00     1.7500E-01    0.07
  Misc. 4 .............. 0.0000E+00    0.00     4.3000E-02    0.02
    Timestep Init ...... 0.0000E+00    0.00     4.3000E-02    0.02
  ----------------------------------------------------------------
  T o t a l s            2.4600E+02  100.00     2.4645E+02  100.00

 Problem time       =    1.0033E-05
 Problem cycle      =       321
 Total CPU time     =       246 seconds (   0 hours  4 minutes  6 seconds)
 CPU time per zone cycle  =       3016.115 nanoseconds
 Clock time per zone cycle=       3018.953 nanoseconds

 Number of CPU's    1
 NLQ used/max     136/  136
 Start time   09/19/2023 15:42:25  
 End time     09/19/2023 15:46:30  
 Elapsed time     245 seconds for     321 cycles using  1 SMP thread 
             (      0 hour   4 minutes  5 seconds)

 N o r m a l    t e r m i n a t i o n                          09/19/23 15:46:30

 

Avatar
Rong Jin September 20, 2023
Error 30281

Hello, 

I am receiving this error when I try to run an model with SPH part

__________

 *** Error 30281 (INI+281)
     SPH initialization
     Material number 5
     material type 20 is not implemented

__________

I don't understand this error because material number 5 is not material type 20 in my model. 

I would appreciate every suggestion. 

Thanks, 

 

 

 

Avatar
kranthi balusu September 19, 2023
Negative Volume

Hello,

I have modeled a timber panel in LS-DYNA using *MAT_WOOD. The analysis goes well until a negative volume error terminates the analysis.  

I found some information and approaches that can help to overcome negative volumes in the following link:

https://www.dynasupport.com/howtos/material/negative-volume-in-soft-materials

In materials that undergo extremely large deformations, such as soft foams, an element may become so distorted that the volume of the element is calculated as negative. This may occur without the material reaching a failure criterion. There is an inherent limit to how much deformation a Lagrangian mesh can accommodate without some sort of mesh smoothing or remeshing taking place. A negative volume calculation in LS-DYNA will cause the calculation to terminate unless ERODE in *CONTROL_TIMESTEP is set to 1 and DTMIN in *CONTROL_TERMINATION is set to any nonzero value in which case the offending element is deleted and the calculation continues (in most cases). Even with ERODE and DTMIN set as described, a negative volume may cause an error termination.

I followed the approaches and set the ERODE=1 and DTMIN= nonzero, however, I understood that the analysis is very sensitive to the value of DTMIN which is a reduction factor for initial time step size to determine minimum time step. When I set it as a number equal to or greater than 1, all elements will be deleted at the beginning of the analysis. When I set it as a value between zero and 0.1 (e.g. DTMIN=0.08), all elements will be deleted at some point at the middle of the analysis. It is also noteworthy that when I set it as a number between 0.1 and one (e.g. DTMIN=0.8), error termination occurs because of out of range velocity of some nodes.

Therefore, I wonder how I can find the proper value for DTMIN so that I can get rid of such a negative volume error. I would appreciate it if you could give me some advice in this regard or any other approaches that can help to overcome negative volumes for my model.

Thanks in advance.

Regards,

Mehdi

Avatar
Mehdi Saloo September 19, 2023
Can Ansys Workbench LS-DYNA handle FSI problem? Ansys LS-DYNA features?

I am new to LS-DYNA and have a couple of questions about LS-DYNA and its relation to ANSYS. I want to do a fluid-structure interaction simulation with immersed boundary method. I am confused with the LS-DYNA available in ANSYS workbench. Based on an old user guide (2017) I found, it mentioned that workbench LS-DYNA cannot handle two-way FSI problem. Is this feature included in the new version of Ansys?

If not, is it true that the FSI problem can only be done via LS-RUN? Can I use Ansys product launcher to call LS-DYNA solver to run FSI problem directly?

I am also curious about whether you can use ANSYS-LS-DYNA to do preprocessing. I noticed that you can export k files from workbench setup for solids, can you do that for fluid as well? Can you export user-defined material from Ansys to k file?

Thanks!

Avatar
Eric Yu September 19, 2023
Finer mesh allows material to penetrate through

I have a simple ball impacting a plate (with some added erosion to simulate a perforation). At a medium sized mesh I get expected results

For the much finer mesh (below; mesh hidden since it is too fine). The ball seems to dimple the surface but then just ghosts right through the plate. Nothing was changed in the problem setup except for mesh size. 

 

I am using AUTOMATIC_SINGLE_SURFACE with SOFT = 1 but this happened with SOFT = 2, also. Presumably this is something to do with the contact stiffness, since all formulations SOFT=0,1,2 have contact stiffness dependent on either the mesh size or the timestep (via, in part, the mesh size). I would have thought this would be handled more robustly for such a simple problem then have to manually tweak contact scale factors... 

Perhaps there is an atlernative contact scheme that is more robust? If not, how do I know how much contact stiffness to add to get the same results as the coarse mesh. It is very counter intuitive that a coarse mesh produces better results in this case. 

Avatar
Jacob Ellis September 18, 2023
SOFT = 2 (Segment Based) Contact Downsides

Hi everyone,

I see a lot of discussion of the benefits of using SOFT = 2 over SOFT = 0 or 1 but outside of the added computational complexity are there any obvious downsides to SOFT = 2?  I feel like I want to move torwards using it as a default algorithm since my problems are often computationally small. Seems to be robust. Any thoughts on edge cases?

Avatar
Jacob Ellis September 18, 2023
Material model for non-newtonian fluid

Hello, I want to model fresh concrete as a Bingham/Herschel-Bulkley fluid using SPH method in Ansys-LS-DYNA. It'd be the best if I can input a loading curve of shear stress vs shear rate for the material model, or a combo of yield stress vs time and viscosity vs time. I was unable to find the exact material model compatible with SPH to use [I checked the material selector tool]. I was wondering if you can share some insights of alternative material models I can use? If fresh concrete specific models are unavailable, is there a generalized non-newtonian fluid model I can use? Am I able to customize a material model? Thank you!

Avatar
cs11 September 15, 2023
Scripting Command Language LsPrePost – Python

Hi,

I would like to create some parametric models in LsPrePost using Python. I have seen some examples on https://ftp.lstc.com/anonymous/outgoing/lsprepost/SCLexamples/lsppscripting.pdf 

However, I cannot find any documentation for the modules (DataCenter, LsPrePost), help, and how to start creating my models.

If you could please help me with finding more resources to achieve my goal. I'd appreciate.

Thank you.

Avatar
Stefano Del Rosso October 1, 2023
Drilling SPH particles

Hi all.

I also simulate rock drilling using SPH particles.

I use the increased density to speed up the solution to the problem. Increasing the density for SPH is not recommended, but I did such tests, I also do tests without increasing the density, they have not yet completed. But I have questions.

The task is similar to LSTC

https://www.lstc-cmmg.org/drilling

Only I use the rock and my feed and drilling speeds. I use an average speed of 25 m/s feed speed and 27 rad/s drilling speed. While solving the problem, I encountered several problems. 

1. Particles scatter wildly. And when the particles fly apart, they pull the solid ring behind them and it is greatly distorted:

2.I reduced the density a little, the solution became better and I got good stress. I also replaced the EROSION contact with AUTOMATIC:

But the particles still fly away strongly and do not break off in large pieces. Why? how to deal with this? It happens like this: the drill pierces the particles, they begin to fly away more and more until the end of the calculation time.

I use concrete material JHC