Ansys Ls-Dyna Users Guide
Ansys Ls-Dyna Users Guide
Ansys Ls-Dyna Users Guide
Guide
ANSYS Release 10.0
002184
August 2005
ANSYS, Inc.
Southpointe
275 Technology Drive
Canonsburg, PA 15317
[email protected]
http://www.ansys.com
(T) 724-746-3304
(F) 724-514-9494
Disclaimer Notice
THIS ANSYS SOFTWARE PRODUCT AND PROGRAM DOCUMENTATION INCLUDE TRADE SECRETS AND ARE CONFIDENTIAL AND PROPRIETARY PRODUCTS
OF ANSYS, INC., ITS SUBSIDIARIES, OR LICENSORS. The software products and documentation are furnished by ANSYS, Inc., its subsidiaries, or affiliates
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. and ANSYS Europe, Ltd. are UL registered ISO 9001:2000 Companies.
Third-Party Software
See the online documentation in the product help files for the complete Legal Notice for ANSYS proprietary software and third-party software. The ANSYS
third-party software information is also available via download from the Customer Portal on the ANSYS web page. If you are unable to access the thirdparty legal notices, please contact ANSYS, Inc.
Published in the U.S.A.
Table of Contents
1. Introduction ........................................................................................................................................ 11
1.1. Overview of Steps in an Explicit Dynamic Analysis ......................................................................... 11
1.2. Commands Used in an Explicit Dynamic Analysis ........................................................................... 11
1.3. A Guide to Using this Document ................................................................................................... 13
1.4. Where to Find Explicit Dynamics Example Problems ...................................................................... 14
1.5. Additional Information ................................................................................................................. 14
2. Elements .............................................................................................................................................. 21
2.1. Solid and Shell Elements ............................................................................................................... 21
2.1.1. SOLID164 ............................................................................................................................. 21
2.1.2. SHELL163 ............................................................................................................................. 23
2.1.2.1. General Shell Formulations .......................................................................................... 23
2.1.2.2. Membrane Element Formulation ................................................................................. 24
2.1.2.3. Triangular Shell Formulations ...................................................................................... 24
2.1.3. PLANE162 ............................................................................................................................ 26
2.1.4. SOLID168 ............................................................................................................................. 28
2.2. Beam and Link Elements ............................................................................................................... 29
2.2.1. BEAM161 ............................................................................................................................. 29
2.2.2. LINK160 ............................................................................................................................. 210
2.2.3. LINK167 ............................................................................................................................. 210
2.3. Discrete Elements ....................................................................................................................... 210
2.3.1. COMBI165 Spring-Damper ................................................................................................. 210
2.3.2. MASS166 ........................................................................................................................... 211
2.4. General Element Capabilities ....................................................................................................... 211
3. Analysis Procedure .............................................................................................................................. 31
3.1. Build the Model ............................................................................................................................ 31
3.1.1. Define Element Types and Real Constants ............................................................................. 31
3.1.2. Specify Material Properties ................................................................................................... 32
3.1.3. Define the Model Geometry ................................................................................................. 32
3.1.4. Mesh the Model ................................................................................................................... 32
3.1.5. Define Contact Surfaces ....................................................................................................... 33
3.1.6. General Modeling Guidelines ............................................................................................... 34
3.2. Apply Loads and Obtain the Solution ............................................................................................ 34
3.2.1. Loads ................................................................................................................................... 34
3.2.2. Initial Velocities .................................................................................................................... 35
3.2.3. Constraints .......................................................................................................................... 36
3.2.4. DOF Coupling ...................................................................................................................... 36
3.2.5. Data Smoothing ................................................................................................................... 36
3.2.6. Specify Explicit Dynamics Controls ....................................................................................... 36
3.2.7. Save Database and Solve ...................................................................................................... 37
3.3. Review the Results ........................................................................................................................ 37
3.4. The Definition of Part .................................................................................................................... 37
3.4.1. Part Assemblies .................................................................................................................. 310
3.5. Adaptive Meshing ....................................................................................................................... 310
4. Loading ................................................................................................................................................ 41
4.1. General Loading Options .............................................................................................................. 41
4.1.1. Components ........................................................................................................................ 42
4.1.2. Array Parameters ................................................................................................................. 42
4.1.3. Applying Loads .................................................................................................................... 44
4.1.4. Data Curves ......................................................................................................................... 45
4.1.4.1. Using Data Curves with Material Models ...................................................................... 45
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
vi
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
vii
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
ix
List of Figures
2.1. Integration Points ............................................................................................................................... 26
4.1. Constrained Shell to Solid ................................................................................................................... 49
6.1. LS-DYNA Drawbead Representation .................................................................................................... 69
7.1. Surface of the Two-invariant Cap Model ............................................................................................ 722
9.1. Hourglass Deformations ..................................................................................................................... 91
11.1. Time Step Sizes Before and After Subcycling .................................................................................... 111
16.1. High Speed Impact of a Metal Bar .................................................................................................... 161
16.2. Lagrangian Impact Solution ............................................................................................................ 162
16.3. Eulerian Channel Flow Solution ....................................................................................................... 162
16.4. ALE Impact Solution ........................................................................................................................ 163
17.1. Two Views of the Target .................................................................................................................. 177
17.2. Drop Test Set-up Dialog Box - Basic Tab ......................................................................................... 1710
17.3. Drop Test Set-up Dialog Box - Velocity Tab .................................................................................... 1712
17.4. Drop Test Set-up Dialog Box - Target Tab ....................................................................................... 1713
17.5. Drop Test Set-up Dialog Box - Status Tab ....................................................................................... 1715
17.6. Graph and Time-History Variables Dialog Box ................................................................................ 1718
List of Tables
3.1. Loads Applicable in an Explicit Dynamics Analysis ............................................................................... 35
3.2. LS-DYNA Solution and Output Control Options ................................................................................... 36
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Chapter 1: Introduction
ANSYS LS-DYNA combines the LS-DYNA explicit finite element program with the powerful pre- and postprocessing
capabilities of the ANSYS program. The explicit method of solution used by LS-DYNA provides fast solutions for
short-time, large deformation dynamics, quasi-static problems with large deformations and multiple nonlinearites,
and complex contact/impact problems. Using this integrated product, you can model your structure in ANSYS,
obtain the explicit dynamic solution via LS-DYNA, and review results using the standard ANSYS postprocessing
tools.
You can also transfer geometry and results information between ANSYS and ANSYS LS-DYNA to perform sequential
implicit-explicit / explicit-implicit analyses, such as those required for droptest, springback and other applications.
The following ANSYS LS-DYNA general information topics are available:
1.1. Overview of Steps in an Explicit Dynamic Analysis
1.2. Commands Used in an Explicit Dynamic Analysis
1.3. A Guide to Using this Document
1.4. Where to Find Explicit Dynamics Example Problems
1.5. Additional Information
2.
Apply loads and obtain the solution (with the SOLUTION processor)
3.
This document describes procedures and concepts that are unique to an explicit dynamic analysis performed
using the ANSYS LS-DYNA product. It does not describe all details of the three steps listed above. If you are familiar with the ANSYS program, you already know how to perform these steps, and this document will provide additional information you need to perform an explicit dynamic analysis. If you have not used ANSYS before, you
should review the following two manuals (in addition to this document) to learn the basic analysis procedures:
When using ANSYS LS-DYNA, it is recommended that you use the default settings provided by the program.
Many times, these settings are appropriate for solving the problems.
EDALE
Chapter 1: Introduction
EDASMP
EDBOUND
EDBVIS
EDBX
EDCADAPT
EDCGEN
EDCLIST
EDCMORE
EDCNSTR
EDCONTACT
EDCPU
EDCRB
EDCSC
EDCTS
EDCURVE
EDDAMP
EDDBL
EDDC
EDDRELAX
EDDUMP
EDENERGY
EDFPLOT
EDGCALE
EDHGLS
EDHIST
EDHTIME
EDINT
EDIS
EDIPART
EDLCS
EDLOAD
Specifies loads
EDMP
EDNB
EDNDTSD
EDNROT
EDOPT
EDOUT
EDPART
EDPC
EDPL
EDPVEL
12
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
EDRD
EDREAD
EDRI
Defines inertia properties for a new rigid body that is created when a
deformable part is switched to rigid
EDRST
EDSHELL
EDSOLV
EDSP
EDSTART
EDTERM
EDTP
EDVEL
EDWELD
EDWRITE
PARTSEL
RIMPORT
REXPORT
UPGEOM
Adds displacements from a previous analysis and updates the geometry to the deformed configuration
For detailed alphabetized descriptions of the ANSYS commands (including specific menu path information for
each command), see the ANSYS Commands Reference.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
13
Chapter 1: Introduction
Other more advanced capabilities that are not required to complete an analysis, but that may be useful in some
situations, are presented in:
Chapter 9, Hourglassing
Chapter 10, Mass Scaling
Chapter 11, Subcycling
Chapter 13, Restarting
Chapter 14, Explicit-to-Implicit Sequential Solution
Chapter 15, Implicit-to-Explicit Sequential Solution
Chapter 16, Arbitrary Lagrangian-Eulerian Formulation
Chapter 17, Drop Test Module
Finally, the appendices contain reference information on the following topics:
Appendix A:, Comparison of Implicit and Explicit Methods
Appendix B:, Material Model Examples
Appendix C:, ANSYS LS-DYNA to LS-DYNA Command Mapping
14
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Chapter 2: Elements
The following elements can be used in an explicit dynamic analysis:
LINK160
BEAM161
PLANE162
SHELL163
SOLID164
COMBI165
MASS166
LINK167
SOLID168
All of the explicit dynamic elements listed above are 3-D except for PLANE162, and all have reduced integration
by default when applicable (for example, reduced integration is not the default for mass or link elements). Reduced
integration means that the number of points for numerical integration in the element formation process is less
than that necessary for exact integration. Therefore, one integration point is the default formulation for both
solid and shell elements. Full integration options are also available for the solid and shell elements. For additional
details, see Chapter 9, Hourglassing and also the LS-DYNA Theoretical Manual.
These elements assume a linear displacement function; higher order elements with a quadratic displacement
function are not available. Therefore, the explicit dynamic elements are not available with extra shape functions,
midside nodes, or p-elements. Explicit elements with linear displacement functions and one point integration
are best suited for nonlinear applications with large deformations and material failure.
It should be noted that the explicit elements are not linked directly to material properties. For example, the
SOLID164 element supports more than 20 different material models, including elastic, plastic, rubber, and foam
models. Unless stated otherwise (see Chapter 6, Contact Surfaces), the minimum material property requirements
for all elements are density, Poisson's ratio, and elastic modulus. See Chapter 7, Material Models of this document
for detailed information on the material properties that can be used in an explicit dynamic analysis. See also the
ANSYS Elements Reference for a complete description of each element, including details of input and output
quantities.
The following explicit element topics are available:
2.1. Solid and Shell Elements
2.2. Beam and Link Elements
2.3. Discrete Elements
2.4. General Element Capabilities
Chapter 2: Elements
the LS-DYNA Theoretical Manual. If hourglass phenomenon is a concern, such as with foams, the fully integrated
formulation may perform better because hourglass control is not required; but it is about four times more costly
in terms of CPU time.
KEYOPT(5) defines the element continuum treatment. Two different formulations are available: Lagrangian (default)
and Arbitrary Lagrangian-Eulerian (ALE). See Arbitrary Lagrangian-Eulerian Formulation in the ANSYS LS-DYNA
User's Guide for more information.
Wedge, pyramid, and tetrahedra shaped SOLID164 elements are simply degenerate bricks (i.e., some of the nodes
are repeated). These shapes are often too stiff in bending and cause problems in some situations. Therefore,
these degenerate shapes should be avoided. If a tetrahedron mesh is desired, it is recommended that SOLID168
elements be used instead of degenerated SOLID164 elements.
You can use any of the following material models for solid elements:
Isotropic Elastic
Orthotropic Elastic
Anisotropic Elastic
Bilinear Kinematic
Plastic Kinematic
Viscoelastic
Blatz-Ko Rubber
Bilinear Isotropic
Composite Damage
Concrete Damage
Geological Cap
Honeycomb
Mooney-Rivlin Rubber
Elastic-Plastic Hydrodynamic
Viscous Foam
Crushable Foam
Johnson-Cook Plasticity
Null
Zerilli-Armstrong
22
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Bamman
Steinberg
Elastic Fluid
2.1.2. SHELL163
Twelve different formulations are available for the SHELL163 element. Use KEYOPT(1) to specify the desired formulation. As with solid elements, the number of integration points per element directly impacts CPU time.
Therefore, for general analyses, the reduced integration shell formulations are recommended. The following is
an outline of the different formulations available for SHELL163. Coordinate systems for output of quantities associated with SHELL163 may vary for different formulations. See the ANSYS Elements Reference for more information.
Belytschko-Leviathan (KEYOPT(1) = 8)
Has four integration points in plane and does not need hourglass control.
Hughes-Liu (KEYOPT(1) = 1, 6, 7, 11) = > four different formulations which can offset the mid-plane of
element away from the nodes.
KEYOPT(1) = 1 General Hughes-Liu. Has one-point integration and is 2.5 times slower that BelytschkoTsay.
KEYOPT(1) = 11 Fast (corotational) Hughes-Liu. Has one-point integration and is 1.5 times slower than
Belytschko-Tsay.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
23
Chapter 2: Elements
KEYOPT(1) = 6 S/R Hughes-Liu. Has four integration points with no hourglassing but is 20 times slower
that Belytschko-Tsay.
KEYOPT(1) = 7 S/R corotational Hughes-Liu. Has four integration points with no hourglassing. 8.8 times
slower than Belytschko-Tsay. Recommended if hourglassing is a problem in an analysis.
No hourglassing
The description of SHELL163 in the ANSYS LS-DYNA User's Guide also has a complete list of available shell formulations.
Degenerate quadrilateral shell elements are prone to lock under transverse shear; therefore, C0 triangular shell
elements (based on work by Belytschko and coworkers) have been implemented. Triangular shells can be mixed
with quadrilateral shells within the same material property set, provided that the element sorting flag (the ITRST
field on the EDSHELL command) is set to 1 (which is the default).
You can use any of the following material models for shell elements:
Isotropic Elastic
Orthotropic Elastic
24
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Bilinear Kinematic
Plastic Kinematic
Blatz-Ko Rubber
Bilinear Isotropic
Composite Damage
Johnson-Cook Plasticity
Bamman
Note When the Mooney-Rivlin Rubber material model is used with SHELL163 elements, the LS-DYNA
code will automatically use a total Lagrangian modification of the Belytschko-Tsay formulation instead
of using the formulation you specify via KEYOPT(1). This program-chosen formulation is required to address
the special needs of the hyperelastic material.
All shell element formulations can have an arbitrary number of integration points through the thickness. Typically,
2 integration points (default) are required through the thickness for elastic behavior, while 3 or more integration
points are required for plastic behavior. The number of integration points through the thickness is controlled
using the second real constant:
R,NSET,,R2 where R2 = number of integration points (NIP)
3-D plane stress constitutive subroutines are implemented for the shell elements; these update the stress tensor
such that the stress component normal to the shell mid surface is zero. The integration points are stacked vertically
at the centroid of the element, as shown in Figure 2.1: Integration Points.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
25
Chapter 2: Elements
#
" !
Through-thickness directions at each node are initially normal to the element surface but rotate with the nodes.
The through-the-thickness integration is needed to calculate bending moments and in-plane forces. The distribution of strain is always linear; whereas the stress is more complicated in nature and depends on the material
law.
Two integration points are sufficient for linear elastic materials, while more points are required for nonlinear
materials. Stress output is given at the outermost integration points, not at the surfaces (despite the nomenclature
of postprocessors, which refer to top and bottom surfaces), so care is needed in interpretation of results. For
elastic materials, stresses can be extrapolated to the surfaces. For nonlinear materials the usual procedure is to
choose four or five integration points through the thickness and to ignore the inaccuracy (i.e., the difference in
stress between the surface and the outermost integration point). The locations of the outermost integration
points of Gauss quadrature are given in the following table:
Mid Plane
0
2 points
0.7746
0.8611
5 points
Outer Surface
3 points
4 points
Outermost Point
0.5774
0.9062
1.0000
Note Keep in mind that this integration could be made exactly beforehand when linear elastic material is used, but it is not done in ANSYS LS-DYNA because nonlinear behavior is generally being modeled.
Also, for elements using full integration, the output stress is the averaged stress values from the 2x2 integration points of the same layer.
2.1.3. PLANE162
The PLANE162 element is a 2-D, 4-node solid element that can be used either as a planar (X-Y plane) or as an
axisymmetric (Y-axis of symmetry) element. KEYOPT(3) allows you to specify a plane stress, axisymmetric, or
plane strain option for the element. For the axisymmetric element formulation, you can specify either an area
or volume weighted option using KEYOPT(2).
26
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Isotropic Elastic
Orthotropic Elastic
Elastic Fluid
Viscoelastic
Bilinear Isotropic
Bilinear Kinematic
Plastic Kinematic
Composite Damage
Johnson-Cook Plasticity
Bamman
For the plane stress option (KEYOPT(3) = 0), you can also choose the following materials:
For the axisymmetric and plane strain options (KEYOPT(3) = 1 or 2), you can also choose the following materials:
Blatz-Ko Rubber
Mooney-Rivlin Rubber
Elastic-Plastic Hydrodynamic
Crushable Foam
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
27
Chapter 2: Elements
Honeycomb
Null
Zerilli-Armstrong
Steinberg
2.1.4. SOLID168
The SOLID168 element is a 10-node tetrahedron element. It uses five point integration; see the element description
for SOLID168 in the ANSYS Elements Reference.
You can use any of the following material models for solid elements:
Isotropic Elastic
Isotropic Elastic
Orthotropic Elastic
Anisotropic Elastic
Bilinear Kinematic
Plastic Kinematic
Viscoelastic
Blatz-Ko Rubber
Bilinear Isotropic
Composite Damage
Concrete Damage
Geological Cap
Honeycomb
Mooney-Rivlin Rubber
Elastic-Plastic Hydrodynamic
Viscous Foam
Crushable Foam
Johnson-Cook Plasticity
Null
Zerilli-Armstrong
28
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Bamman
Steinberg
Elastic Fluid
SOLID168 elements are well suited to modeling irregular meshes such as those produced from various CAD/CAM
systems and the ANSYS Workbench. It is important to note, however, that models made up entirely of SOLID168
elements may not be as accurate as hexahedral (SOLID164) models. For this reason, we recommend using a
combination of SOLID168 and SOLID164 elements within a model as a good modeling practice. The SOLID168
elements could be used for the irregular portions of the model while the SOLID164 elements could be used to
mesh the more uniform features.
Although SOLID164 and SOLID168 elements can be used in the same model, you must be careful to develop an
adequate mesh transition between these two element types. This is due to the fact that SOLID168 elements
conatin mid-side nodes, which are not allowed in SOLID164 elements.
Isotropic Elastic
Bilinear Kinematic
Plastic Kinematic
Viscoelastic
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
29
Chapter 2: Elements
2.2.2. LINK160
The LINK160 spar (truss) element is similar to the Belytschko-Schwer beam, but carries axial loads only. This element
assumes a straight bar, axially loaded at its ends, and of uniform properties from end to end. The materials you
can use for this element are isotropic elastic, plastic kinematic (rate dependent), and bilinear kinematic.
2.2.3. LINK167
The LINK167 element is a tension only spar that can be used to model cables. It is similar to a spring element in
that the relationships between force and deformation are input directly by the user. This element requires that
the cable option be specified with the EDMP command (see the EDMP command description).
Elastoplastic Spring
When using COMBI165, be sure to specify a unique set of real constants, element type, and material properties
(R, ET, and TB commands, respectively) for each part to ensure that parts are uniquely defined.
210
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
2.3.2. MASS166
The mass element is defined by a single node and a mass value (Force*Time2/Length). Mass elements can often
be used in a model to realistically characterize the mass of a structure without incorporating a large number of
solid or shell elements. In an automobile impact analysis, for example, mass elements can model engine components whose deformation behavior is not of primary interest. Such a use of mass elements will reduce the
number of elements required in an analysis, hence reducing the total computation time required to obtain a
solution.
You can also use the MASS166 element to define a lumped rotary inertia at a node. To use this option, set KEYOPT(1) = 1 in the MASS166 element definition and input six moment of inertia values (IXX, IXY, IXZ, IYY, IYZ, IZZ)
via the element real constants. You can not input a mass value for this option; therefore, you must define a second
mass element at the same node (with KEYOPT(1) = 0) to account for mass.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
211
212
2.
3.
If you have never used an ANSYS product before, you should review the ANSYS Basic Analysis Guide and the ANSYS
Modeling and Meshing Guide to learn the general procedures used to build a model in ANSYS.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
32
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Avoid degenerate shell and solid element shapes (such as a triangular shell or a tetrahedral solid). They
are generally too stiff and less accurate than quadrilateral and hexahedral shapes.
Try to achieve a mesh with uniform element sizes (i.e., avoid areas with relatively small elements). A large
difference in element sizes can cause a small minimum time step size and, therefore, a long run time. If a
relatively few number of small elements are required to mesh a particular geometry, mass scaling (see
Chapter 10, Mass Scaling) can be used to increase the minimum time step.
Do not use the SmartSizing method of element size control (SMRTSIZE command), which can create a
large variation in element sizes within the mesh. Instead, use the ESIZE and related commands to control
element sizes.
Do not make the mesh too coarse when reduced element formulations are used, or elements may experience hourglassing.
If hourglassing is a problem, try to use fully integrated elements in a part of the model or in the entire
model.
After specifying the desired mesh controls, mesh the model using typical meshing commands (AMESH, VMESH,
etc.). In the GUI, meshing functions are accessed by picking Main Menu> Preprocessor> Meshing> Mesh or
via the Mesh Tool (Main Menu> Preprocessor> Meshing> Mesh Tool).
Determine the type of contact which best defines the physical model (EDCGEN command)
Identify contact surfaces (CM, or EDPART, or EDASMP commands, with the EDCGEN command)
Specify additional input required for the chosen contact type (EDCGEN and EDCONTACT commands)
If automatic contact is not being used, you must define contact surfaces by grouping the nodes on each surface
into a component with the CM command. Once these components are created, use the EDCGEN command to
specify contact between the desired surfaces (i.e., node components). You must also identify the type of contact
with the EDCGEN command. For the single surface contact algorithms, the outer surfaces of the contacting
bodies are determined by ANSYS LS-DYNA.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
33
Use rigid bodies whenever possible to represent relatively stiff, unyielding objects in the model. Using rigid
bodies simplifies the solution and results in shorter run times. (See Chapter 8, Rigid Bodies for information
on defining rigid bodies.)
Use realistic values for material properties. For example, do not use an unrealistically high elastic modulus
to represent a rigid body, and do not use unrealistic thickness values for shell elements.
Consider using damping (EDDAMP command) to prevent unrealistic oscillations in your models structural
response. Refer to the EDDAMP command description in the ANSYS Commands Reference for detailed
information.
If you have performed 2-D dynamic analyses with the regular ANSYS program, consider extruding these
models to 3-D and analyzing them with ANSYS LS-DYNA. You may achieve more accurate results in a
shorter run time.
Note that the submodeling and substructuring features of the ANSYS program cannot be used in ANSYS
LS-DYNA.
3.2.1. Loads
In an explicit dynamic analysis, all loads must be specified over time using component logic or part IDs, array
parameters, and the EDLOAD command. (In the GUI, load options are found under Main Menu> Solution>
Loading Options> Specify Loads.) Basic input for this command is a component name or part number and two
array parameter names or a load curve ID number (LCID). The component specified must contain the nodes or
elements on which the load is being applied. The array parameters specified must contain time varying load
data (one array for time values and one array for the corresponding load values; the two arrays must be the same
length). As an alternative to inputting the array parameters on the EDLOAD command, you can define the load
curve using the EDCURVE command (Main Menu> Solution> Loading Options> Curve Options) and input
the load curve ID on EDLOAD. Valid loads are shown in the table below.
34
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Label
Displacements
UX, UY, UZ
Rotations
Forces
FX, FY, FZ
Moments
MX, MY, MZ
Velocities
VX, VY, VZ
AX, AY, AZ
Base Accelerations
Angular Velocities
Temperature
TEMP
PRESS
Note We recommend that you specify velocity time histories instead of displacement time histories.
A piecewise linear displacement time history may lead to discontinuous velocities and infinite accelerations.
The load symbol will appear automatically on the active graphics window. The load symbol is erased automatically
when you replot. To turn the display of this symbol on or off, issue EDFPLOT (or pick Main Menu> Solution>
Loading Options> Show Forces in the GUI).
To visualize the applied load curve, use the EDPL command (or pick Main Menu> Solution> Loading Options>
Plot Load Curve in the GUI).
For more detailed information on applying loads with the EDLOAD command, see Section 4.1: General Loading
Options.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
35
3.2.3. Constraints
In addition to loads and initial velocities, constraints can also be applied to the model. Using the D command,
you can apply constraints only to the displacement (UX, UY, UZ) and rotation (ROTX, ROTY, ROTZ) degrees of
freedom; the constraint value must be zero. (In the GUI, pick Main Menu> Solution> Constraints> Apply.)
Other specialized types of constraints are also available in an explicit dynamic analysis (such as symmetry
boundary planes and welded constraints). See Section 4.2.1: Constraints and Section 4.2.2: Welds for more information.
Command
GUI Path
Termination Time
TIME
SHELL/BEAM Output
EDINT
EDRST
EDHTIME
Energy Options
The above table does not show all controls available in an explicit dynamic analysis. However, most of the default
settings for the LS-DYNA control options (output controls, file controls, damping options, etc.) are sufficient for
most explicit dynamic analyses and need not be modified.
A brief description of those options that are recommended follows:
36
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
POST1 is used to review results over the entire model at specific time-points.
POST26 is used to track specific nodal and element result items over a more detailed load history.
The following are basic guidelines for postprocessing explicit dynamic results.
The database must contain the same model for which the solution was calculated.
All stresses and strains output from LS-DYNA are in the global Cartesian coordinate system. However, if
you are using composite materials, stresses can be in a local (element) coordinate system.
For more information on postprocessing explicit dynamic results, see Chapter 12, Postprocessing later in this
document. For a complete description of all postprocessing functions, see Chapter 4, An Overview of Postprocessing in the ANSYS Basic Analysis Guide.
37
MAT TYP
REL
ESY
SEC
NODES
11
12
13
14
15
16
17
18
21
22
23
24
25
26
27
28
31
32
33
34
35
36
37
38
41
42
43
44
45
46
47
48
51
52
53
54
55
56
57
58
61
62
63
64
65
66
67
68
71
72
73
74
75
76
77
78
81
82
83
84
85
86
87
88
11
101
102
103
104
105
106
107
108
12
111
112
113
114
115
116
117
118
13
121
122
123
124
125
126
127
128
38
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
MAT TYP
REL
ESY
SEC
NODES
14
131
132
133
134
135
136
137
138
15
141
142
143
144
145
146
147
148
To define a user-specified PART ID, first select elements 3 and 5 and create an element component E35. EDPART,ADD,10,E35 specifies PART ID 10 for elements 3 and 5. PART ID 10 (for the component E35) will not change
if you issue EDPART,CREATE or EDPART,UPDATE after reselecting all of the elements in the model. Consider the
definition of a second user-specified PART ID 60 for element 14. Selecting element 14 and issuing EDPART,ADD,60
changes the PART ID for element 14 from 6 to 60. After selecting all of the elements, EDPART,UPDATE shows
the following part list.
PART
MAT
TYP
REAL
USED
10
2 - user defined
60
1 - user defined
In the above list, elements 1,2 and 4 have MAT = TYPE = REAL = 1 and form PART 1. Elements 3 and 5 also have
MAT = TYPE = REAL = 1 but have the user-specified PART ID 10. Element 14 with MAT = 1, TYPE = 3 and REAL =
2 has the user-specified PART ID 60. The PART IDs for the rest of elements were generated by the ANSYS LS-DYNA
program. Element 10 had MAT = TYPE = REAL = 2, but was deleted after the part list was created initially.
Therefore, USED for PART ID 4 is zero. The USED field in the part list indicates the number of elements having
the PART ID. USED is calculated from all currently selected elements. If the element selection changes, USED will
change. To get the correct number of USED, select the correct number of elements and issue EDPART,UPDATE.
Note For this example, issuing EDPART,UPDATE will not change the PART IDs in the part list. However,
issuing EDPART,CREATE will create a new part list. PART 5 will become the new PART 4 and PART 6 will
become the new PART 5. The user-specified PART IDs 10 and 60 will not change. Because the automatically
generated PART IDs can change as a result of EDPART,CREATE, this can invalidate previously defined
part-based loading, contact specifications, etc.
The following procedure is recommended:
a.
Build the model as usual until you have to specify a PART ID number on a part-related command.
b.
Create the part list (EDPART,CREATE) or specify PART IDs for some groups of elements (EDPART,ADD);
then list the parts (EDPART,LIST).
c.
Use the appropriate PART ID number(s) from the list for the command(s).
d.
Continue modeling.
e.
Update the part list (EDPART,UPDATE) if the elements or their attributes have changed, which includes
the addition of new elements.
f.
List the updated PART ID numbers and make use of them on further part-related command(s).
Issuing ELIST with the part listing key on (PTKEY = 1) will show the PART ID for each elements as follows.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
39
REL
ESY
SEC
NODES
11
12
13
14
15
16
17
18
10
21
22
23
24
25
26
27
28
31
32
33
34
35
36
37
38
10
41
42
43
44
45
46
47
48
51
52
53
54
55
56
57
58
61
62
63
64
65
66
67
68
71
72
73
74
75
76
77
78
81
82
83
84
85
86
87
88
11
101
102
103
104
105
106
107
108
12
111
112
113
114
115
116
117
118
13
121
122
123
124
125
126
127
128
14
60
131
132
133
134
135
136
137
138
15
141
142
143
144
145
146
147
148
In addition to listing PART IDs in the element list, PART IDs can also be plotted with the EPLOT command by
turning on the part display key (/PNUM,PART,1). Element selection using the ESEL command can also be performed
based on PART IDs. These two capabilities replace the obsolete PARTSEL command.
The CDWRITE command writes the part information to the Jobname.CDB file. This information can then be
automatically read into ANSYS with the CDREAD command. However, if more than one Jobname.CDB file is
read, the part list from the last Jobname.CDB file overwrites the existing part list, if any. In general, you must issue
the EDPART,CREATE command to recreate the part list of the total model. This will affect all part-related commands
contained in the Jobname.CDB file. Therefore, you can join models, but not part-related inputs, which must be
modified using the newly-created PART ID numbers.
In limited cases, an update of the part list (EDPART,UPDATE) is possible. This requires that no used combination
of MAT/TYPE/REAL appears more than once in the list. However, partial changes to the part-related commands
may be necessary.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Note Adaptive meshing is only valid for parts consisting of SHELL163 elements.
When adaptivity is turned on for a part, the mesh for that part will be automatically regenerated during an analysis to ensure adequate element aspect ratios are maintained throughout the deformation process. Adaptive
meshing is most commonly used in the analysis of large deformation processes such as metal forming (adaptive
meshing is typically used for the blank). To activate adaptive meshing for more than one part within a model,
you must reissue the EDADAPT command for each different PART ID. By default, adaptive meshing is off for all
parts within a model.
After specifying which parts will be adaptively meshed, you must also define adaptive meshing parameters using
the EDCADAPT command. The EDCADAPT command globally sets the control options for all PART IDs that are
to be adaptively meshed as defined by the EDADAPT command. The parameters controlled by the EDCADAPT
command are listed below.
Frequency (FREQ) - The time interval between which adaptive mesh refinements will take place. For example, if FREQ is set to 0.01, elements will be adaptively remeshed every .01 seconds (assuming the time
unit used is seconds) if they are above the specified angle tolerance. Since the default value for the FREQ
field is 0.0, this field must be specified when using adaptive meshing in an analysis.
Angle Tolerance (TOL) - The adaptive angle tolerance (in degrees) for which adaptive meshing will occur
(default = 1e31). The TOL field controls the aspect ratio between elements and is very important for ensuring accurate results. If the relative angle between elements exceeds the specified TOL value, the elements
will be automatically refined.
Adaptivity Option (OPT) - There are two different angle options available for mesh adaptivity. For OPT =
1, the angle change that is compared to the specified TOL value is computed based on the original mesh
configuration. For OPT = 2, the angle change that is compared to the specified TOL value is computed
based on the previously refined mesh.
Mesh Refinement Levels (MAXLVL) - The MAXLVL field controls the number of times an element can be
remeshed during the entire analysis. For a single original element, a MAXLVL of 1 would allow the creation
of one additional element, a MAXLVL of 2 would allow the creation of up to four additional elements, and
a MAXLVL of 3 would allow the creation of up to 16 additional elements. High values of MAXLVL will yield
more accurate results, but can dramatically increase the model size.
Remeshing Birth and Death Times (BTIME and DTIME) - The remeshing birth and death times control
when adaptive meshing is turned on and off in an analysis. For example, if you set BTIME = .01 and DTIME
= .1, adaptive meshing will only take place in the analysis between .01 and .1 seconds (assuming the time
unit used is seconds).
Interval of Remeshing Curve (LCID) - This data curve (specified on the EDCURVE command) defines the
interval of remeshing as a function of time. The abscissa of the data curve is time, and the ordinate is the
varied adaptive time interval. If this option is nonzero, the adaptive frequency (FREQ) will be replaced by
this data curve. Note, however, that a nonzero FREQ value is still required to initiate the first adaptive loop.
Minimum Element Size (ADPSIZE) - Minimum element size to be adapted based on element edge length.
If this parameter is undefined, the edge length limit is ignored.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
311
One or Two Pass Option (ADPASS) - If ADPASS = 0, two pass adaptivity is used, and the calculation is repeated
after adaptive remeshing (this is the default). If ADPASS = 1, one pass adaptivity is used, and the calculation
is not repeated after adaptive remeshing. See Figures 30.9(a) and 30.9(b) in the ANSYS/LS-DYNA Theoretical
Manual for graphical representations of these two options.
Uniform Refinement Level Flag (IREFLG) - Values of 1, 2, 3, etc. allow 4, 16, 64, etc. elements, respectively,
to be created uniformly for each original element.
Penetration Flag (ADPENE) - Depending on whether the value of ADPENE is positive (approach) or negative
(penetrate), the program will adapt the mesh when the contact surfaces approach or penetrate the tooling
surface by this value. The tooling adaptive refinement is based on the curvature of the tooling. If ADPENE
is positive, the refinement generally occurs before contact takes place; consequently, it is possible to use
one pass adaptivity (ADPASS = 1).
Shell Thickness Level (ADPTH) - Absolute shell thickness level below which adaptivity should begin. This
option works only if the adaptive angle tolerance (TOL) is nonzero. If thickness based adaptive remeshing
is desired without angle change, set TOL to a large angle. (If ADPTH = 0.0, this option is not used.)
Maximum Element Limit (MAXEL) - Maximum number of elements at which adaptivity will be terminated.
Adaptivity is stopped if this number of elements is exceeded.
For most problems, you should do an initial analysis without adaptive meshing. If the results of the initial analysis
show a highly distorted mesh and if the solution appears to be incorrect, then you can rerun the analysis with
adaptive meshing turned on. Adaptive meshing may also be used when an analysis terminates in LS-DYNA with
a negative volume element error.
When adaptive meshing is turned on, the number of elements within the model will change during the solution.
After each adaptive loop, the mesh will be updated and a new results file will be written with the extension RSnn,
where nn is the adaptive mesh level. (An adaptive loop will occur at every time increment specified by FREQ or
at the time intervals specified by LCID.) For example, a model that is adaptively remeshed two times will produce
two results files, Jobname.RS01 and Jobname.RS02. For details on postprocessing adaptively meshed results,
see Section 12.2.3: Postprocessing after Adaptive Meshing.
Note A new results file will be created for each loop, even if the mesh does not change.
A new time history file having the extension HInn is also written each loop; however, these files may not be as
useful for postprocessing as the RSnn files. In addition, LS-DYNA creates a set of adaptive meshing files with the
root name "adapt". Therefore, you should not use "adapt" as your jobname when adaptive meshing is activated.
312
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Chapter 4: Loading
After building the model, the next step is to apply loads to the structure in preparation for solution. In order to
properly model your structures behavior, it is necessary to apply loads with respect to a specified time interval.
This chapter discusses the following topics related to loads:
Temperature loading
Dynamic relaxation
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Chapter 4: Loading
loads in ANSYS LS-DYNA are applied using a pair of array parameters, one corresponding to the time and the
other corresponding to the loading condition.
Note Although nodal accelerations (AX, AY, AZ) and nodal velocities (VX, VY, VX) appear as degrees
of freedom, they are not physical DOFs and cannot be constrained with the D command. To apply these
nodal loads, use the EDLOAD command.
In ANSYS LS-DYNA, all loads are applied in one load step. This is much different from an implicit analysis where
loads are often applied in multiple steps. In ANSYS LS-DYNA, for certain kinds of loads, you can also specify when
a load is imposed on a body (birth time), and when the load is removed (death time) using the EDLOAD command.
Please refer to Birth Time, Death Time, and CID Support in the EDLOAD command to check the applicability of
the birth/death time.
To apply a load to your model, you will need to follow these steps:
Designate portions of the model that will receive the load as components (or parts, for rigid bodies)
Define array parameters containing time intervals and load data values
Define the load direction using the EDLCS command if the load is not acting in the global coordinate
system
4.1.1. Components
With the exception of loads on rigid bodies, all loads are applied to components in an explicit analysis. Therefore,
the first step in applying loads is to gather portions of your model into nodal or element groups called components.
Each component should consist of those portions of the model that will receive the same loads and that are related
to each other through such things as material properties, position in the model, expected behavior, etc.
For example, if you want to analyze a baseball being hit into a wall, you would probably define the nodes of the
ball as one component, the nodes on the bat as another component, and the nodes on the wall as a third component.
You may define as many components as you wish, and then apply loads to each of them. The components must
consist of either nodes or elements. (Components consisting of elements are only used when applying pressure
loads.) To define a component, you must first select only those entities that you want to be included in the
component. Then use the CM command or the following menu path in the GUI to define the component: Utility
Menu> Select> Comp/Assembly> Create Component.
See Selecting and Components in the ANSYS Basic Analysis Guide for more information on components.
For the case of rigid bodies, loads are applied to a specific part number rather than to components. This is due
to the fact that rigid bodies already consist of a specific set of nodes and elements as defined by the command
EDMP,RIGID,MAT.
42
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
!
!
!
!
!
!
select
define
select
define
select
define
all
the
the
the
the
the
Now select the nodes at the base of the bat handle (bathand) for which a displacement versus time loading curve
will be applied.
Define an array named time to store the time values. Remember to use time values that are consistent with all
loads, dimensions, and material properties in your model. Next, define an array named xdisp to contain the corresponding X displacements of the component bathand. Likewise, define arrays ydisp and zdisp to contain the
corresponding Y and Z displacements of bathand.
nsel,s,node,...
! select the nodes at the bat handle base
cm,bathand,node
! define the component bathand
nsel,all
*dim,time,,4
! dimension the array parameter time
*dim,xdisp,,4
! dimension the array parameter xdisp
*dim,ydisp,,4
! dimension the array parameter ydisp
*dim,zdisp,,4
! dimension the array parameter zdisp
time(1)=0,0.25,0.5,0.75,1 ! times at specified displacements
xdisp(1)=0,-1,-2,-1,3
! X displacement of the bat handle
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
43
Chapter 4: Loading
ydisp(1)=0,1,2,3,4
zdisp(1)=0,3,6,8,9
The example given is actually a simple version of a more complex phenomenon. In a more accurate simulation,
additional displacement locations (and corresponding load curves) should be defined to better simulate the
true motion of the bat handle. Further, the initial velocities of the nodes comprising the bat are all different. Finally,
the ball is a composite of several different materials and material models.
44
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
! X displacement of bathand
! Y displacement of bathand
! Z displacement of bathand
Alternatively, the very stiff bat could have been modeled as a rigid body to simplify the required input and to
reduce the CPU time. For this case, the appropriate rigid body loads would be applied to valid part numbers
(and not to nodal components).
edload,add,rbux,,2,time,xdisp
edload,add,rbuy,,2,time,ydisp
edload,add,rbuz,,2,time,zdisp
If the pressure that the bat imparts on the baseball was known as a function of time, the loading could be accomplished without defining any contact surfaces. Instead, the element component containing the elements on the
surface of the ball would be loaded with the PRESS label of the EDLOAD command.
edload,add,press,1,cover,battime,batload
In this last scenario, the elements contained in the element component cover are loaded on face number one
(face number is input in the KEY field) with the load curve produced by the battime and batload array parameters.
It is important to note that in the examples listed above, the pressure load uses the KEY field of the EDLOAD
command to specify a load key. Load keys (1,2,3, etc.) are associated with surface loads and are listed under
"surface loads" in the input data tables for each element type in the ANSYS Elements Reference. For most load
types other than pressure loads, you can use the KEY field to specify a coordinate system identification number,
CID, from the EDLCS command. The load will act in a direction you define on EDLCS, or in the global coordinate
system direction if CID is not specified. See Section 4.1.5: Defining Loads in a Local Coordinate System for more
information.
Note To avoid timing problems on some platforms, it is a good practice to always add a very small
time value (such as 1.0 10-6) to the value in the final item in the time array. For example, instead of the
value 3.0, such an array might contain the following value for the last item:
timeint(1)=0,1,2,3.00001
The addition of this very small padding factor does not affect the accuracy of the results.
In addition to adding loads, you can also list and delete loads with the commands EDLOAD,LIST and EDLOAD,DELE. You can display load symbols on the element plot with the EDFPLOT command, and you can plot
load curves with the EDPL command.
45
Chapter 4: Loading
Similar to placing loads on components, data curves are grouped into array parameters and then associated
with a particular curve reference number (LCID) which can be used by a specified material model (PLAW, HONEY,
etc.) or contact type (RNTR, ROTR) and drawbead. Use the following steps to specify data curves:
1.
Define an array parameter that contains the abscissa values for the material or friction force behavior
(e.g., effective plastic strain, effective strain rate, displacement, etc.).
2.
Define a second array parameter that contains the ordinate values for the material behavior or friction
force (e.g., initial yield stress, elastic modulus, force, etc.).
3.
Specify the data curve that will be used to identify this data [EDCURVE]. During this step you will choose
a data curve ID number that will be used to associate this data with a particular material behavior when
setting up the data table [TBDATA].
After the parameters have been specified, use one of the following menu paths in the GUI to define a data curve:
Main Menu> Preprocessor> Material Props> Curve Options
Main Menu> Solution> Loading Options> Curve Options
Defined data curves can be listed using EDCURVE,LIST. They can be plotted using EDCURVE,PLOT and can be
deleted using EDCURVE,DELETE.
The following example shows how data curves can be used to define a piecewise linear plasticity curve
(TB,PLAW,,,,8) for steel:
! "3" was arbitrarily chosen as the material reference (MAT) number.
mp, ex,3,30.0e6
! elastic (Youngs) modulus (psi)
mp,dens,3,7.33e-4
! mass density (lbf-sec^2/in^4)
mp,nuxy,3,0.30
! Poissons ratio (unitless)
!Note: First convert engineering stress versus engineering strain data
!
into true stress versus true (hencky) strain data. Then subtract
!
off the elastic true strain from the total true strain to find
!
the plastic true strain, which is used with the total true stress
!
in LS-DYNA *MAT_PIECEWISE_LINEAR_PLASTICITY material model #24.
!-----------------------------------------------------------------------! Stress-Strain Data used with Piecewise Linear Plasticity (Power Law 8):
!-----------------------------------------------------------------------!
Total
Total
Total
Total
Elastic
Plastic
! Stress/
Eng.
Eng.
True
True
True
True
! Strain
Stress
Strain
Stress
Strain
Strain
Strain
! Point
(psi)
(in/in)
(psi)
(in/in)
(in/in)
(in/in)
!-----------------------------------------------------------------------!
1
0
0.0000
0
0.0000
0.0000
0.0000
!
2
60,000
0.0020
60,120
0.0020
0.0020
0.0000
!
3
77,500
0.0325
80,020
0.0320
0.0027
0.0293
!
4
83,300
0.0835
90,260
0.0802
0.0030
0.0772
!
5
98,000
0.1735
115,000
0.1600
0.0038
0.1562
!
6
98,300
0.2710
124,940
0.2398
0.0042
0.2356
!
7
76,400
1.2255
170,030
0.8000
0.0057
0.7943
!-----------------------------------------------------------------------!Note: The first point on the stress/strain curve is NOT entered.
!
Start with the second point (where ordinate = yield stress).
!
Also, please follow the limits imposed by the *SET command.
*dim,strn,,6 ! define array for effective plastic true strain data
*dim,strs,,6 ! define array for effective total true stress data
strn(1)= 0.0, 0.0293, 0.0772, 0.1562, 0.2356, 0.7943 ! strain (in/in)
strs(1)= 60120., 80020., 90260., 115000., 124940., 170030. ! stress (psi)
edcurve,add,1,strn,strs ! curve #1: abscissa=strain & ordinate=stress
tb,plaw,3,,,8
! specify power law #8 for material (MAT) #3
tbdata,6,1
! use load curve #1 for stress/strain data
!Note: If desired, a plastic failure strain can be defined. Further,
!
strain rate effects can be included by specifying the necessary
!
strain rate parameters or the load curve defining the strain rate
!
scaling effect on the yield stress. Please refer to Chapter 7 of
!
this guide for a complete description of this material model.
46
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
As shown above, the LCID is referenced in the ninth field of the EDLOAD command. It is important to note that
if you use an LCID with the EDLOAD command, array parameters should not also be specified to define the
load. Using an LCID instead of array parameters on the EDLOAD command is particularly useful when a specific
load curve is to be applied to multiple components or load labels.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
47
Chapter 4: Loading
4.2.1. Constraints
Unlike the ANSYS implicit program, ANSYS LS-DYNA does distinguish between zero and nonzero constraints.
Nonzero constraints are handled as loads (together with a load curve; see earlier discussion in this chapter). Only
zero constraints can be applied with the D command, i.e., the value specified must always be 0 (zero). No other
values are valid, as the D constraints are only used to fix certain portions of the model. You may use symmetry/antisymmetry boundary conditions to apply these zero constraints.
You can use the EDNROT command to apply zero constraints in terms of a rotated nodal coordinate system.
You must first define the local coordinate system with the EDLCS command.
When modeling small symmetric sections of a geometry, you may need to define sliding or cyclic symmetry. You
can use the EDBOUND command to define symmetry boundary planes for sliding or cyclic symmetry. You use
nodal components to identify the boundaries and a direction vector to define a normal (sliding symmetry) or
axis of rotation (cyclic symmetry).
Various other types of constraints can be modeled in ANSYS LS-DYNA using the EDCNSTR command. The
available constraint types are extra node set (ENS), nodal rigid body (NRB), shell edge to solid (STS), and rivet
(RIVET). In the GUI, you access these constraints by picking:
Main Menu> Solution> Constraints> Apply> Additional Nodal
The extra node set constraint type (EDCNSTR,ADD,ENS) allows the addition of nodes (via a nodal component)
to an existing rigid body that was defined with the EDMP command. The nodal component that is added must
not be attached to any other rigid body. The extra nodes that are added to a rigid body may be located anywhere
in the model, and may have coordinates outside those of the original rigid body. The ENS option has many potential applications, including placing nodes where joints will be attached between rigid bodies, defining nodes
where point loads will be applied, and defining a lumped mass at a specific location.
Unlike typical rigid bodies that are defined with the EDMP command, nodal rigid bodies defined with the EDCNSTR,ADD,NRB command are not associated with a part number, but are based on a node component. The NRB
option can be advantageous when modeling rigid (welded) joints in a model. For a rigid joint, portions of different
flexible components (having different MAT IDs) act together as a rigid body. It is difficult to define this type of
rigid body with a unique MAT ID (and corresponding part number). However, the rigid joint can be easily defined
using a nodal rigid body. Because nodal rigid bodies are not associated with a part number, other options that
use rigid bodies (such as loads applied with the EDLOAD command ) cannot be used with a nodal rigid body.
The shell to solid edge option (EDCNSTR,ADD,STS) ties regions of solid elements to regions of shell elements.
As shown in Figure 4.1: Constrained Shell to Solid, a single shell node may be tied to up to nine solid nodes
that define a 'fiber' vector. The solid element nodes that define the fiber remain linear (straight) throughout the
analysis but can move relative to each other in the fiber direction. The shell node must be coincident with one
of the solid element nodes along the fiber.
48
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
& $21$0 !()%&
# '
4
3
5
7
6
8
9$
Similar to spotwelds defined with the EDWELD command, the RIVET option (EDCNSTR,ADD,RIVET) defines a
massless rigid constraint between two nodes that are noncoincident. Unlike a spotweld, however, failure cannot
be specified for a rivet. When a rivet is defined, the distance between the nodes is kept constant throughout any
motion that occurs during a simulation. Nodes connected by a rivet cannot be part of any other constraints
specified in the model.
4.2.2. Welds
In an explicit dynamic analysis, it is common to model components that are physically welded together. This is
especially true in automotive applications where parts are often assembled using welds. In such cases, welded
constraints can be modeled in ANSYS LS-DYNA using the EDWELD command. Two different types of welds can
be modeled: massless spotwelds and generalized welds. Nodes connected by the EDWELD command cannot
be constrained in any other way.
For a massless spotweld, you must specify two noncoincident nodes. You can define failure within the spotweld
by inputting failure parameters on the EDWELD command. This failure is based on the following relationship:
fn
S
n
exp n
fs
+
S
s
exp s
In the GUI, use the following menu path to define a massless spotweld: Main Menu> Preprocessor> LS-DYNA
Options> Spotweld> Massless Spotwld
A generalized weld is typically used to model a long-welded section of two parts. For a generalized weld, you
must specify a valid nodal component. Coincident nodes are permitted. However, if you use coincident nodes,
you must specify a local coordinate system (EDLCS) which will be used for output data. Failure within the generalized weld is also defined with the failure parameters on the EDWELD command using the same relationship
as shown in the equation above.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
49
Chapter 4: Loading
In the GUI, use the following menu path to define a generalized weld: Main Menu> Preprocessor> LS-DYNA
Options> Spotweld> Genrlizd Spotwld
Define the entity to which you want to apply the initial velocity. This may be a single node, a nodal
component [CM], a part [EDPART], or a part assembly [EDASMP].
2.
Decide whether the VELO or VGEN option is appropriate for your application.
3.
Specify the initial velocity [EDVEL or EDPVEL] by using one of the following menu paths in the GUI:
Main Menu> Solution> Initial Velocity> On Nodes> w/Nodal Rotate (VELO option)
Main Menu> Solution> Initial Velocity> On Nodes> w/Axial Rotate (VGEN option)
Main Menu> Solution> Initial Velocity> On Parts> w/Nodal Rotate (VELO option)
Main Menu> Solution> Initial Velocity> On Parts> w/Axial Rotate (VGEN option)
For either the VELO or the VGEN option, input the translational velocities relative to the global Cartesian coordinate
system (VX, VY, and VZ fields on the EDVEL or EDPVEL command). For the VELO option, input the nodal rotations
in terms of the global Cartesian coordinate system (OMEGAX, OMEGAY, and OMEGAZ fields). For the VGEN option,
input the magnitude of the angular velocity (OMEGAX), the coordinate axis for rotation (XC, YC, and ZC fields) and
directional angles (ANGX, ANGY, and ANGZ fields) relative to the global X, Y, and Z axes.
If you do not specify initial velocities using EDVEL or EDPVEL, all initial velocities in the model will be zero.
Likewise, if you only specify the Cname field on the EDVEL command (e.g., EDVEL,VGEN,Cname) or the PID field
on the EDPVEL command, (e.g., EDPVEL,VGEN,PID), zero initial velocities will be applied since the remaining
fields on these two commands default to zero.
If you wish to change an initial velocity previously specified using EDVEL, just respecify a new velocity with the
same component name or node number. This new value will overwrite the old data for the component or node.
To list or delete initial velocities previously applied to node components or nodes, use the EDVEL,LIST and
EDVEL,DELE commands (in the GUI, use Main Menu> Solution> Initial Velocity> On Nodes> List and Main
Menu> Solution> Initial Velocity> On Nodes> Delete).
410
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Or, assuming the component ball was part number 1 in the model, you could specify the velocity of the baseball
as:
edpvel,velo,1,,,-1600.0 ! Part number 1 initially moving
! in -Z direction at 1600 in/sec
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
411
Chapter 4: Loading
dilatational and shear activation flags. For example, the following command defines a nonreflecting boundary
on the component ground:
EDBOUND,ADD,GROUND,1,1
To list and delete defined nonreflecting boundaries, use the commands EDNB,LIST and EDNB,DELE, respectively.
Temperature results from an ANSYS thermal analysis applied as non-uniform temperature loads (that do
not vary with time) in a subsequent explicit dynamic analysis (LDREAD; also requires a sequential solution)
The first method utilizes the EDLOAD,,TEMP command and the general loading procedure to apply a timevarying temperature to a specific nodal component. You must define two array parameters to represent the
load; the first parameter contains the time values and the second parameter contains the temperature values.
You may choose to define a load curve with these parameters or to input them directly on the EDLOAD command.
You can also scale the temperature values using the SCALE parameter on EDLOAD. See Section 4.1: General
Loading Options for a complete description of how to apply this type of load.
The second method allows you to apply a uniform constant temperature to all nodes in the model. This method
may be used to model a structure subjected to steady-state thermal loading. You can use the TUNIF command
or the BFUNIF,TEMP command to apply this type of temperature load.
The third method allows you to apply the temperatures calculated in an ANSYS thermal analysis as loads in an
explicit dynamic analysis. This method is useful for modeling temperature-dependent phenomena such as forging.
To use this method, you must perform an implicit-to-explicit sequential solution. In the explicit phase, you use
the LDREAD command to read the temperature data from the thermal (implicit) analysis results file (Jobname.RTH)
and apply them as loads to the nodes in the model. You can only transfer the temperatures from one specified
time point in the thermal analysis. Refer to Chapter 15, Implicit-to-Explicit Sequential Solution of this guide for
a detailed description of this loading procedure.
For all three methods of temperature loading, a reference temperature can be input via the TREF command.
The thermal loading is defined as the difference between the applied temperature and the reference temperature.
If the reference temperature is not specifically defined, it defaults to zero.
In order for any of these thermal loads to take effect, you must use the temperature dependent bilinear isotropic
material model or the elastic viscoplastic thermal material model. The temperature dependent bilinear isotropic
model can be used to represent a thermoelastic material (if needed) by omitting the yield strength and tangent
modulus. Please refer to Chapter 7, Material Models of this guide for complete descriptions of these material
models.
Note It is important to note that the EDLOAD method of temperature loading cannot be used together
with the LDREAD or TUNIF / BFUNIF methods of temperature loading. In addition, EDLOAD cannot be
used to list or delete temperatures applied with LDREAD, TUNIF, or BFUNIF.
412
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
There are essentially five different types of analyses that can be conducted with ANSYS LS-DYNA concerning
dynamic relaxation. They are discussed below:
1.
Transient Dynamic Analysis Only (default - EDDRELAX,OFF): In this case, the PHASE parameter of the
EDLOAD command is set to zero and dynamic relaxation is not used. This is the default setting.
2.
Stress Initialization Only with No Transient Analysis (EDDRELAX,DYNA): This case basically uses the ANSYS
LS-DYNA explicit solver to approximate the solution of a static analysis that is really best handled by an
implicit solver such as ANSYS. The PHASE parameter on the EDLOAD command is set to one. The termination time (TIME command) must be set to zero to prevent unloading of the static load. True dynamic
relaxation is used.
3.
Stress Initialization with a No-Load Transient Analysis (EDDRELAX,DYNA): This case is the same as the
previous case, except that the termination time is set to the desired value. At time zero, the structure is
instantly unloaded and free to vibrate.
4.
Stress Initialization with a Loaded Transient Analysis (EDDRELAX,DYNA): This is similar to the previous
case, except that the load curve specified by the EDLOAD command (via PHASE = 2) will be used in both
pseudo time for the dynamic relaxation to obtain the preload and in real time for the transient analysis.
If a ramped load curve is used instead of a constant load curve, then the structure will unload at time
zero and reload again. Therefore, if a ramped load curve is needed for the explicit solver to accurately
converge on the static solution, it would be better to use two EDLOAD commands, the first (with PHASE
= 1) to apply the preload using a ramped load curve and the second (with PHASE = 0) to continue the
load with a load curve that is not ramped.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
413
Chapter 4: Loading
Note If you issue EDLOAD multiple times for the same component (or Part ID) and the same
load label, the values associated with the last EDLOAD command are used. Consequently, you
cannot use multiple values of the PHASE parameter for a given component (or Part ID) and load
label. To work around this limitation, you must use duplicate nodal components for the stress
initialization and transient curves. (This workaround is applicable only if you use components to
define the load.)
5.
414
Implicit-to-Explicit Sequential Solutions (EDDRELAX,ANSYS): This case uses the displacement results
from the ANSYS implicit solver to apply the preload for the ANSYS LS-DYNA explicit solver via stress initialization to a prescribed geometry. The PHASE parameter on the EDLOAD command must be set to
zero. This scenario is explained in detail in Chapter 15, Implicit-to-Explicit Sequential Solution of this
guide. Fields 2-6 of the EDDRELAX command are ignored for this case.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
The header records, including geometric quantities (such as nodes and elements), are written on the
two results files Jobname.RST and Jobname.HIS. (At this point, the ANSYS LS-DYNA database should
contain all relevant information. Before issuing the SOLVE command, you should issue a SAVE command
to write all model information to the file Jobname.DB.)
2.
The input file for the LS-DYNA program, Jobname.K, is written using all information that has been entered
so far.
3.
Control is transferred from the ANSYS LS-DYNA program to the LS-DYNA program. The solutions produced
by the LS-DYNA solver are written to the results files Jobname.RST and Jobname.HIS. If the EDOPT,ADD,,BOTH option was specified before issuing the SOLVE command, results files for the LS-POST
postprocessor (files d3plot and d3thdt) will also be written.
When the solution is complete without errors or warnings, the ANSYS LS-DYNA GUI notifies the user that the
solution is done, and control is transferred back to the ANSYS LS-DYNA program. The results can be viewed using
the POST1 and POST26 processors of the ANSYS LS-DYNA program. If errors or warnings are produced, messages
pop-up and are also displayed in the ANSYS Output Window. The messages state the number of errors and
warnings that were produced, and make reference to the LS-DYNA message file, where details of the errors and
warnings are written. These same details are also written to the LS-DYNA d3hsp file.
Termination time - Use the TIME command to specify an end-time for the analysis. The calculation will
stop when the accumulation of time steps reaches that end-time.
CPU time limit - You can use the EDCPU command to specify a CPU time limit (in seconds). The calculation
will stop when that time limit is reached.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Termination criteria - You can use the EDTERM command to specify that the solution stop when a specific
node or rigid body reaches a certain position, or when a specific node comes into contact with another
surface. You may set up multiple criteria using this command. (See the EDTERM command for details.)
You should always specify an analysis end-time using the TIME command. The other termination controls are
optional. The solution will terminate when any one of the specified termination criteria is met.
When you are performing a small restart (EDSTART,2) or a full restart (EDSTART,3) analysis, if the previous analysis terminated due to a satisfied termination criteria that was set by EDTERM, you must change or delete that
criteria so that it will not cause the restart to terminate immediately.
Specify the number of processors to be used. You can set this value via the config100.ans file, via the
/CONFIG command, or via the SETNPROC macro. For more information on all of these methods, please
refer to Section 15.1: Activating Parallel Processing in the ANSYS Advanced Analysis Techniques Guide. If
you do not set the number of processors to be greater than one, you will not be able to see the appropriate menu selections in the GUI to specify SMP, nor will the EDRUN command settings have any effect.
2.
When you are using shared memory parallel processing, the calculations may be executed in different order,
depending on CPU availability and the workload on each CPU. Because of this, you may see slight differences in
the results when you run the same job multiple times. To avoid these differences, you can specify that consistency
be maintained (EDRUN,SMP,1). Maintaining consistency can result in an increase of up to 15% in CPU time.
If you are using all CPUs on a platform (ANSYS recommends always running one fewer CPU than is available),
you should close all other applications on those CPUs. Any calculations or machine resources that compete with
the ANSYS LS-DYNA application running in SMP mode will reduce the SMP performance significantly.
The SMP capability is expected to scale linearly when used with up to 8 processors. When used with more than
8 processors, any additional scalability is minimal; the larger the problem, the higher the potential for scalability.)
52
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
0
0
0
0
mins)
mins)
mins)
mins)
At this point, you can enter one of the following four switches:
sw1
sw2
ANSYS LS-DYNA responds with time and cycle numbers and continues to run. This allows
you to see how far the solution has actually progressed.
sw3
sw4
The first estimate of CPU time reported in the console window (before the issuance of CTRL-C) is usually inaccurate
(see estimated total cpu time in the output example shown above). You can use CTRL-C to interrupt the execution
of the LS-DYNA solver, then type the sense switch sw2 to get a better estimate of the execution time and cycle
numbers. The ANSYS LS-DYNA program writes all important messages (errors, warnings, failed elements, contact
problems, etc.) to the ANSYS LS-DYNA console window and to the LS-DYNA ASCII output file d3hsp.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
53
Option 1 - Plot the elements with the smallest time step sizes. Each element is shaded red or yellow based
on its time step size. The number of elements plotted is based on the VALUE1 field (default = 100) of the
EDTP command.
Option 2 - Plot and list the elements with the smallest time step sizes. This option produces the same plot
as Option 1, and also provides a listing of the smallest elements.
Option 3 - This option produces a plot similar to Option 1, except that all selected elements are plotted.
Elements beyond the number specified by VALUE1 are blue and translucent, with the degree of translucency
based on the VALUE2 field of the EDTP command.
In the plot produced by EDTP the elements are shaded red (smallest), yellow (intermediate), or blue (largest)
based on their time step sizes.
Note Care should be taken when using the EDTP command for large models because the time step
size calculations may take a significant amount of CPU time.
54
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Although these capabilities cannot be directly accessed, users who are familiar with the LS-DYNA program can
still use these features indirectly in the ANSYS LS-DYNA program. The procedure for including additional features
is outlined below.
Note Only the LS-DYNA features that are described in the ANSYS LS-DYNA documentation and available
through the ANSYS LS-DYNA interface are supported features. Although you can access other features
as described below, features that you implement by editing the LS-DYNA input file are not supported.
After modeling has been completed, you can issue the EDWRITE command instead of issuing the SOLVE command. (In the GUI, pick Main Menu> Solution> Write Jobname.k.) EDWRITE will create the LS-DYNA input file
Jobname.K. (Note that if the analysis is a small restart, the input file will be named Jobname.R, and if it is a full
restart, the input file will be named Jobname_nn.K.)
The difference between the EDWRITE command and the SOLVE command is that step number 3 discussed
under Solution Process at the beginning of this chapter will not be executed. That is, Jobname.K will be written
(along with headers to the ANSYS results files), but the solution process of the LS-DYNA program will not be
started. After issuing EDWRITE, use one of the following two methods to include additional functionality in the
explicit dynamic analysis.
5.7.1. Method A
1.
2.
Edit the LS-DYNA input file Jobname.K to incorporate the items that are not supported directly through
the ANSYS LS-DYNA program.
3.
Execute the LS-DYNA executable separately in the same directory where the two results files, Jobname.RST and Jobname.HIS, reside. The solution produced by the LS-DYNA executable will be appended
to these results files.
4.
After completion of the LS-DYNA execution, enter the ANSYS LS-DYNA program and view the results
using the postprocessors (POST1 and POST26).
5.7.2. Method B
1.
Using another window, go to the directory where the ANSYS LS-DYNA program was started. Do not exit
the ANSYS LS-DYNA program in the original window.
2.
Edit the LS-DYNA input file Jobname.K to incorporate the items that are not supported directly through
the ANSYS LS-DYNA program.
3.
Execute the LS-DYNA executable separately in the same directory (where the two results files, Jobname.RST and Jobname.HIS, reside). The solution produced by the LS-DYNA executable will be appended
to these results files. (Remember, you have not exited the ANSYS LS-DYNA program in the original window.)
4.
After completion of the LS-DYNA execution, go back to the original window and view the results using
the ANSYS LS-DYNA postprocessors.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
55
For example, if you have ANSYS Mechanical with LS-DYNA at ANSYS Release 10.0 and your input file is crashtest.k,
you would type:
lsdyna100 i=crashtest.k pr=ansysds
Note The lsdyna100 script is only available on UNIX platforms. On the PC, use the ls970.exe (or
ls970_DP.exe) executable instead to run LS-DYNA.
You may need additional command line arguments, depending on what type of analysis you are running. For
an implicit-to-explicit sequential solution (REXPORT and EDDRELAX commands; see Chapter 15, Implicit-toExplicit Sequential Solution), the argument m=drelax is also required. For a restart analysis, (EDSTART command;
see Chapter 13, Restarting), the argument r=d3dumpnn is also required, where nn is the number of the dump
file (01, 02, ... 99) to be used as a starting point. For any type of analysis, you can specify the memory to be used
by including the command line argument memory=n, where n is the desired amount of memory in words.
When using these methods to modify the file Jobname.K, nodes and elements must not be changed. Also, note
that the ANSYS LS-DYNA database is not updated with the changes you make in file Jobname.K. Thus, if an attempt
is made to view the modified items in the PREP7, POST1 or POST26 processors after the solution is complete,
the original version, and not the modified version, will be shown. (For these reasons, editing the LS-DYNA keyword
input file, Jobname.K, is NOT a supported feature of ANSYS LS-DYNA.) In all cases, the complete set of results
can be postprocessed in the LS-POST postprocessor (assuming you requested LS-DYNA results files via the EDWRITE command).
56
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Automatic Single Surface (ASSC): Contact is established when a surface of one body contacts itself or the
surface of another body. This type is easy to use because no contact or target surface definitions are required. It is efficient for self-contacting problems or large deformation problems where general areas of
contact are not known beforehand.
Automatic General (AG): This contact type is similar to Automatic Single Surface in that contact is established
when a surface of one body contacts itself or the surface of another body. Automatic General contact is
also easy to define since you do not need to specify contact and target surfaces. The main advantage of
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Node-to-Surface (NTS): Contact is established when a contacting node penetrates a target surface. This
type is commonly used for general contact between two surfaces and is most efficient when a smaller
surface comes into contact with a larger surface, such as a thin rod impacting a flat plate.
Surface-to-Surface: (STS) Contact is established when a surface of one body penetrates the surface of another body. This type is commonly used for arbitrary bodies that have large contact areas and is very efficient for bodies that experience large amounts of relative sliding with friction, such as a block sliding on
a plane.
Alternatively, contact surfaces can be specified using part numbers or part assembly numbers that are currently
defined within the finite element model. Part assembly numbers are defined with the EDASMP command.
The input line below demonstrates how the EDCGEN command is used to define contact between two different
parts / part assemblies in a model:
EDCGEN,STS,1,2,.25,.23 ! Generate surface-to-surface contact
! between PARTS 1 and 2
Additionally, contact can be defined between contact and target surfaces using a combination of PART / part
assembly and component definitions, as illustrated below:
EDCGEN,NTS, N1,2,.3,.28 ! Generate node-to-surface contact
! between component N1 and PART 2
EDCGEN,ESTS,1,N2,.15,.15 ! Generate eroding surface-to-surface
! contact between PART 1 and component N2
EDCGEN,STS,1,1,.1,.1 ! Generate surface-to-surface contact
! between PART 1 and itself
You can also use the EDCGEN command for defining contact between part assemblies as shown in the following
sample input:
EDCGEN,STS,5,6,.3,.28 ! Generate surface-to-surface contact
! between part assemblies 5 and 6
For the specialized case of single surface contact (ASSC, AG, ESS, and SS), no contact or target surface definitions
are required. As explained later in this section, single surface contact is the most general type of contact definition
as all external surfaces within a model can be in contact at any point during an analysis. The program will ignore
62
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Note--When defining contact entities in an explicit analysis, no initial penetrations are allowed. Therefore, use
great care when defining contact components.
Flim = VC i A cont
where Acont is the area of the segment contacted by the node in contact. The suggested value for VC is to use
the yield stress in shear:
VC =
o
3
VDC
i crit
100
VDC is the viscous damping coefficient (input on EDCGEN command as percent of critical damping; e.g., VDC =
20 indicates 20%). crit is determined in the following fashion by ANSYS LS-DYNA:
crit = 2m
where
m = min(mslave , mmaster )
= ki
mslave + mmaster
; k is interface stiffness
mslave i mmaster
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
63
You can use the EDPC command to select and plot contact entities. The plot will consist of nodes or elements,
depending on the method that was used to define the contact surfaces (that is, components or parts). You can
specify a minimum contact entity number, maximum contact entity number, and contact entity number increment
by using the MIN, MAX, and INC fields of the EDPC command. Hence, by issuing the command EDPC,1,2,1 for
the above listing, the defined entities for the STS and ASSC contact would be selected and plotted. It is important
to note that for single surface contact definitions, all external surfaces in the model will be selected and plotted.
64
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Nodes to surface
Surface to surface
General (Normal)
SS
NTS
STS, OSTS
Automatic
ANTS
ASTS
Rigid
RNTR
ROTR
Tied
TDNS
TDSS, TSES
TNTS
TSTS
ENTS
ESTS
Eroding
ESS
Edge
SE
Drawbead
DRAWBEAD
Forming
FNTS
FSTS, FOSS
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
65
Node-to-Surface Contact (NTS, ANTS, RNTR, TDNS, TNTS, ENTS, DRAWBEAD, FNTS)
Node-to-surface contact is a contact type which is established when a contacting node penetrates a
target surface. This type of contact is commonly used for general contact between two surfaces. Use the
same rules as in ANSYS implicit to determine which surfaces are target or contact:
The flat or concave surface is the target. The convex surface is the contact surface.
The coarser mesh is the target surface. The finer mesh is the contact surface.
In the case of Drawbead contact, the bead is always the nodal contact surface and the blank is always
the target surface.
3.
Surface-to-Surface Contact (STS, OSTS, ASTS, ROTR, TDSS, TSTS, ESTS, FSTS, FOSS, TSES)
Surface-to-surface contact is established when a surface of one body penetrates the surface of another
body. Surface-to-surface contact is the most general type of contact as it is commonly used for bodies
that have arbitrary shapes with relatively large contact areas. This type of contact is most efficient for
bodies that experience large amounts of relative sliding, such as a block sliding on a plane or a sphere
sliding within a groove.
2.
66
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
4.
5.
6.
fn
fn,fail
m2
f
+ s
fs,fail
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
67
7.
8.
68
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
9.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
69
BOXID
XMIN
Minimum x coordinate
XMAX
Maximum x coordinate
YMIN
Minimum y coordinate
YMAX
Maximum y coordinate
ZMIN
Minimum z coordinate
ZMAX
Maximum z coordinate
Once defined, a BOXID can be used in the BOXID1 and BOXID2 fields of the EDCGEN command. BOXID1 corresponds to the contact box and BOXID2 corresponds to the target contact box.
EDCGEN,NTS,Cont,Targ, ... DTIME,BOXID1,BOXID2
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
area ]
solid elements:
volume
contact depth = min
, 0 .5
area
area
k=
fs Area2 K
Volume
fs Area K
Minimum Diagonal
611
612
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Isotropic (#1)
Orthotropic (#2)
Anisotropic (#2)
Elastic Fluid (# 1)
Viscoelastic (#6)
Foam Models
Honeycomb (# 26)
Null (# 9)
Zerilli-Armstrong (# 65)
Steinberg (# 11)
Elastoplastic Spring
Cable (# 71)
Rigid Model
Rigid (# 20)
72
Choose menu path Main Menu> Preprocessor> Material Props> Material Models. The Define Material Model Behavior dialog box appears.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
In the Material Models Available window on the right, double-click on LS-DYNA, then on one of the
material model categories: Linear, Nonlinear, Equation of State, Discrete Element Properties, or Rigid
Material.
3.
Double-click on a material subcategory. For example, under the Nonlinear category, the subcategories
are: Elastic, Inelastic, and Foam Material Models.
4.
Continue to double-click on any further material subcategories until a data input dialog box appears.
The options that appear in this dialog box will be all of the individual material property options that are
valid for the material model you have selected.
5.
Enter the desired values, and click on OK. The model that you defined is listed by model type and number
in the Material Models Defined window on the left.
You can later edit the values by double-clicking on the material model listed on the left side in the Material
Models Defined window until the associated data input dialog box appears. There, you can edit the data and
click on OK again.
You can copy the contents of an existing material model by choosing Edit> Copy and specifying the new model
number. The copied material model is listed on the left side in the Material Models Defined window with the
new number. Its contents are identical to the model whose contents you copied.
You can delete a material model by clicking on the model number to highlight it, then choosing Edit> Delete.
For further information on defining materials using the GUI, see Section 1.1.4.4: Material Model Interface in the
ANSYS Basic Analysis Guide. Also, see Section 4.1.1.10: Using Tree Structures in the ANSYS Operations Guide for
information specifically on the tree structure layout of the material model interface.
If you define, edit, copy, or delete material models via the GUI, ANSYS will automatically issue the correct commands and write them to the log file.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
73
For an example input listing, see Section B.2.2: Orthotropic Elastic Example: Aluminum Oxide.
%#
$#
#
"
!
)10(&
)'
74
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
K=
E
3(1 2 )
If VAL1 (on EDMP), EX, and NUXY are all specified, VAL1 will be used as the bulk modulus.
1
Sij = G Cij V 1 2 ij
V
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
75
1
W = C10 (I1 3) + C01(I 2 3) + C 1 + D(I 3 1)2
I2
c
C
C = 10 + C01
2
C (5 2) + C01(11 5)
D = 10
2(1 2 )
I1, I2, and I3 are invariants of the right Cauchy-Green Tensor.
Input Poissons ratio () and density with the MP command. (For Poisson's ratio, a value greater than .49 is recommended; smaller values may not work.) Input the Mooney-Rivlin constants with the TB and TBDATA commands.
Data at only one temperature is permitted and must be specified in locations 1 and 2 for the data table:
TB,MOONEY,,,,0
TBDATA, 1, C10
TBDATA, 2, C01
As an alternative to inputting C10 and C01 directly, you can set these constants to zero and supply tabulated
uniaxial data via a load curve. The program will calculate the constants based on the experimental data input in
locations 3-6 of the TBDATA command. To use this input method, you must set the TBOPT field of the TB command
to 2:
TB,MOONEY,,,,2
TBDATA, 1, C10 (set to zero to use experimental data)
TBDATA, 2, C01 (set to zero to use experimental data)
TBDATA, 3, C3 (specimen gauge length L0)
TBDATA, 4, C4 (specimen gauge width)
TBDATA, 5, C5 (specimen thickness)
TBDATA, 6, C6 (load curve ID)
The load curve definition that provides the uniaxial data should give the change in gauge length L versus the
corresponding force. In compression, both the force and the change in gauge length must be specified as negative values. In tension, the force and change in gauge length should be input as positive values. The principal
stretch ratio in the uniaxial direction, 1, is then given by:
76
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
1 =
L o + L
Lo
ij ( )
t
ij = 2o ( t )
d
( t ) = G + (G 0 - G ) e - t
In the model, elastic bulk behavior is assumed when calculating the incrementally integrated pressure from the
volume, V: p = K ln V. The parameters G , G0 , K (Bulk modulus) and are required to define the linear viscoelastic
material model. Input these values with TB,EVISC and locations 46, 47, 48, and 61 of the TBDATA command:
TB,EVISC
TBDATA, 46, G0
TBDATA, 47, G
TBDATA, 48, K
TBDATA, 61, 1/
Note For this material option, you must specify density (DENS) with the MP command.
For an example input listing, see Section B.2.6: Viscoelastic Example: Glass.
77
78
2R
2R + 1 2
1122 + 2
12
R +1
R +1
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
P
22
P
33
The stress-strain behavior can be specified at only one temperature. Input the elastic modulus (EX), density
(DENS), and Poissons ratio (NUXY) with the MP command. Input the yield stress, tangent modulus, anisotropic
hardening parameter, and a load curve ID for effective yield stress versus effective plastic strain with TB,PLAW,,,,7
and locations 1 - 4 of the TBDATA command:
TB,PLAW,,,,7
TBDATA, 1, Y (yield stress)
TBDATA, 2, Etan (tangent modulus)
TBDATA, 3, R (anisotropic hardening parameter)
TBDATA, 4, LCID (load curve ID for yield stress vs. plastic strain)
For an example input listing, see Section B.2.8: Transversely Anisotropic Elastic Plastic Example: 1010 Steel.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
79
This plasticity model is only available for shell and 2-D elements. The model directly follows the plasticity theory
introduced in the Transversely Anisotropic Elastic Plastic model described earlier in this section. You can refer
to that model for the theoretical basis.
To utilize the Transversely Anisotropic FLD model, you must input the density (DENS), elastic modulus (EX), and
Poisson's ratio (NUXY) with the MP command. As shown below, the additional input parameters are specified
using TB,PLAW,,,,10 and locations 1-5 of the TBDATA command.
TB,PLAW,,,,10
TBDATA,1,y (yield stress)
TBDATA,2,Etan (tangent modulus)
TBDATA,3,R (anisotropic hardening parameter)
TBDATA,4,LCID1 (load curve for defining effective stress vs. plastic strain)
TBDATA,5,LCID2 (load curve for defining FLD)
For an example input listing, see Section B.2.9: Transversely Anisotropic FLD Example: Stainless Steel.
710
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
1 + P ( + E eff )
Y =
p p
C 0
where 0 is the initial yield stress, is the strain rate, C and P are the Cowper-Symonds strain rate parameters,
peff is the effective plastic strain, and Ep is the plastic hardening modulus which is given by
Ep =
Etan E
E Etan
The stress-strain behavior can only be specified at one temperature. Input the elastic modulus (EX), density
(DENS), and Poissons ratio (NUXY) with the MP command. Input the yield stress, tangent slope, hardening
parameter, strain rate parameters C and P, and the failure strain with TB,PLAW,,,,1 and locations 1 - 6 of the TBDATA
command:
TB,PLAW,,,,1
TBDATA, 1, Y (yield stress)
TBDATA, 2, Etan (tangent modulus)
TBDATA, 3, (hardening parameter)
TBDATA, 4, C (strain rate parameter)
TBDATA, 5, P (strain rate parameter)
TBDATA, 6, f (failure strain)
For an example input listing, see Section B.2.11: Plastic Kinematic Example: 1018 Steel.
2( Y )m = a K1 + K 2 m+ a K1 K 2 m + c 2K 2 m
where Y is the yield stress, a and c are anisotropic material constants, m is Barlat exponent, and K1 and K2 are
defined by:
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
711
K1 =
xx + h yy
2
2
xx h yy
2 2
K2 =
+ p xy
2
where h and p are additional anisotropic material constants. For the exponential hardening option, the material
yield strength is given by:
y = k (0 + p )n
where k is the strength coefficient, 0 is the initial strain at yield, p is the plastic strain, and n is the hardening
coefficient. All of the anisotropic material constants, excluding p which is determined implicitly, are determined
from Barlat and Lian width to thickness strain ratio (R) values as shown:
a=2 2
R00
R90
1 + R00 1 + R90
c=2a
R00 1 + R90
1 + R00 R90
h=
The width to thickness strain ratio for any angle can be calculated from:
R =
2mm
Y
xx yy
Above, is the uniaxial tension in the direction. The stress-strain behavior can be specified at only one temperature. Input the elastic modulus (EX), density (DENS), and Poissons ratio (NUXY) with the MP command. Input
the hardening rule type, HR (equal to 1 for linear or 2 for exponential), tangent modulus (for HR = 1) or strength
coefficient (for HR = 2), yield stress (for HR = 1) or hardening coefficient (for HR = 2), the Barlat exponent, m, the
width to thickness strain ratio values, R00, R45, and R90, and the orthotropic material axes with TB,PLAW,,,,3 and
locations 1 - 8 of the TBDATA command:
TB,PLAW,,,,3
TBDATA, 1, HR (hardening rule type)
TBDATA, 2, Etan or k (tangent modulus or strength coefficient)
TBDATA, 3, Y or n (yield stress or hardening coefficient)
TBDATA, 4, m (Barlat exponent)
TBDATA, 5, R00
TBDATA, 6, R45
TBDATA, 7, R90
TBDATA,8,CSID (defines orthotropic material axes)
712
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
If CSID = 2, material axes are determined by a local coordinate system specified with the EDLCS command (see
the command description for the details of how the axes are oriented). Before specifying the material property,
you must define the local coordinate system using EDLCS, and then issue EDMP,ORTHO,,VAL1, where VAL1 is
the coordinate system ID from EDLCS.
For an example input listing, see Section B.2.12: 3 Parameter Barlat Example: Aluminum 5182.
= S1 S2 m + S2 S3 m + S3 S1 m
where m is the flow potential exponent and Si are the principal values of the symmetric matrix Sij.
S xx = 1/ 3 [c ( xx yy ) b (zz xx )]
S yy = 1/ 3 [a ( yy zz ) c ( xx yy )]
Szz = 1/ 3 [b (zz xx ) a ( yy zz )]
S yz = f yz
Szx = g zx
S xy = f xy
where a, b, c, f, g, and h represent the anisotropic material constants. When a=b=c=f=g=h=1, isotropic material
behavior is modeled and the yield surface reduces to the Tresca surface for m = 1 and the von Mises surface for
m = 2 or 4. For this material option, the yield strength is given by:
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
713
y = k (P + 0 ) n
where k is the strength coefficient, P is the plastic strain, 0 is the initial strain at yield, and n is the hardening
coefficient. The stress-strain behavior can be specified at only one temperature.
Input the elastic modulus (EX), density (DENS), and Poissons ratio (NUXY) with the MP command. Input the
strength coefficient, the initial strain at yield, the hardening coefficient, the flow potential exponent, and the
Barlat anisotropic constants a-h with TB,PLAW,,,,6 and locations 1 - 10 of the TBDATA command.
TB,PLAW,,,,6
TBDATA, 1, k (strength coefficient)
TBDATA, 2, 0 (initial strain)
TBDATA, 3, n (hardening coefficient)
TBDATA, 4, m (flow potential (Barlat) exponent)
TBDATA, 5, a
TBDATA, 6, b
TBDATA, 7, c
TBDATA, 8, f
TBDATA, 9, g
TBDATA, 10, h
For an example input listing, see Section B.2.13: Barlat Anisotropic Plasticity Example: 2008-T4 Aluminum.
yy = k mn
where is the strain, is the strain rate, k is the material constant, m is the hardening coefficient, and n is the
strain rate sensitivity coefficient. The stress-strain behavior can be specified at only one temperature. Input the
elastic modulus (EX), density (DENS), and Poissons ratio (NUXY) with the MP command. Input the material constant,
hardening coefficient, strain rate sensitivity coefficient, and initial strain rate with TB,PLAW,,,,4 and locations 1 4 of the TBDATA command:
TB,PLAW,,,,4
TBDATA, 1, k (material constant)
TBDATA, 2, m (hardening coefficient)
TBDATA, 3, n (strain rate sensitivity coefficient)
TBDATA, 4, (initial strain rate)
0
For an example input listing, see Section B.2.14: Rate Sensitive Powerlaw Plasticity Example: A356 Aluminum.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
where 0 is the initial yield strength, is the effective strain rate, (p)eff is the effective plastic strain, and Eh is
given by:
Eh =
EEtan
E Etan
The stress-strain behavior can be specified at only one temperature. Input the elastic modulus (EX), density
(DENS), and Poissons ratio (NUXY) using the MP command. Input the load curve ID for defining the initial yield
stress vs. effective strain rate, the tangent modulus, the load curve ID for defining the elastic modulus vs. effective
strain rate (optional), the load curve ID for defining the tangent modulus vs. effective strain rate (optional), and
the load curve ID for defining the von Mises stress at failure vs. effective strain rate using TB,PLAW,,,,5 and locations
1 - 5 of the TBDATA command. For shell elements, you have the option of specifying Mn Time in location 6, instead
of LCID4 in location 5, to define failure of the material. Mn Time is the minimum time step size for automatic
element deletion.
TB,PLAW,,,,5
TBDATA, 1, LCID1 (load curve ID for defining the initial yield stress vs. effective strain rate)
TBDATA, 2, Etan (tangent (plastic hardening) modulus)
TBDATA, 3, LCID2 (the load curve ID for defining the elastic modulus vs. effective strain rate)
TBDATA, 4, LCID3 (the load curve ID for defining the tangent modulus vs. effective strain rate)
TBDATA, 5, LCID4 (the load curve ID for defining von Mises stress at failure vs. effective strain rate)
TBDATA, 6, Mn Time (minimum time step size for automatic element deletion, for shell elements only)
For an example input listing, see Section B.2.15: Strain Rate Dependent Plasticity Example: 4140 Steel.
P
P
P
P
y (eff , eff ) = y (eff ) 1 + ( eff )P
(P )
where is the effective plastic strain rate, C and P are strain rate parameters, and y eff is the yield stress
without considering strain rate. This yield stress quantity can be input as either y and Etan (TBDATA entries 1
and 2) or as a total true stress vs. effective plastic strain curve, LCID1 (TBDATA entry 6).
The second method uses the LCID2 curve (TBDATA entry 7) to scale the yield stress with respect to strain rate.
If LCID2 is input, it is used instead of the Cowper-Symonds scaling term defined by the strain rate parameters C
and P.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
715
716
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Input the strain rate parameters C and P (TBDATA entries 4 and 5) to use the Cowper-Symonds model.
Input the LCID2 curve (TBDATA entry 7) to scale the yield stress with respect to strain.
Define separate stress vs. strain curves for different strain rates (TBDATA entries 11 - 30).
Input the elastic modulus (EX), density (DENS), and Poissons ratio (NUXY) with the MP command. Input the remaining parameters with TB,PLAW,,,,11 and locations 1 - 30 of the TBDATA command:
TB, PLAW,,,,11
TBDATA,1, Y (yield stress)
TBDATA,2, Etan (tangent modulus)
TBDATA,3, F (effective plastic strain at failure)
TBDATA,4, C (strain rate parameter)
TBDATA,5, P (strain rate parameter)
TBDATA,6, LCID1 (load curve defining total true stress vs. effective plastic strain)
TBDATA,7, LCID2 (load curve that scales LCID1 with respect to strain rate)
TBDATA,8, EPSTHIN (thinning plastic strain at failure)
TBDATA,9, EPSMAJ (major in-plane strain at failure)
TBDATA,10, NUMINT (number of through-thickness integration points that must fail for element failure)
TBDATA,11, Rate1 (strain rate value 1)
TBDATA,12, LCID3 (stress vs. effective plastic strain curve for Rate1)
TBDATA,13, Rate2 (strain rate value 2)
TBDATA,14, LCID4 (stress vs. effective plastic strain curve for Rate2)
TBDATA,15, Rate3 (strain rate value 3)
TBDATA,16, LCID5 (stress vs. effective plastic strain curve for Rate3)
TBDATA,17, Rate4 (strain rate value 4)
TBDATA,18, LCID6 (stress vs. effective plastic strain curve for Rate4)
TBDATA,19, Rate5 (strain rate value 5)
TBDATA,20, LCID7 (stress vs. effective plastic strain curve for Rate5)
TBDATA,21, Rate6 (strain rate value 6)
TBDATA,22, LCID8 (stress vs. effective plastic strain curve for Rate6)
TBDATA,23, Rate7 (strain rate value 7)
TBDATA,24, LCID9 (stress vs. effective plastic strain curve for Rate7)
TBDATA,25, Rate8 (strain rate value 8)
TBDATA,26, LCID10 (stress vs. effective plastic strain curve for Rate8)
TBDATA,27, Rate9 (strain rate value 9)
TBDATA,28, LCID11 (stress vs. effective plastic strain curve for Rate9)
TBDATA,29, Rate10 (strain rate value 10)
TBDATA,30, LCID12 (stress vs. effective plastic strain curve for Rate10)
If load curve LCID1 (total true stress vs. effective plastic strain) is input, the yield stress (Y) and tangent modulus
(Etan) values will be overwritten. If LCID2 (the scaling curve) is input, strain rate parameters C and P will be overwritten. If C and P, LCID2, and Rate1 (first strain rate value) are all set to zero, strain rate dependency is ignored.
For an example input listing, see Section B.2.17: Modified Piecewise Linear Plasticity Example: PVC.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
717
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
P
Y = 1 + k[e + P ] n
eff
C
where is the strain rate, C and P are the Cowper-Symonds strain rate parameters, e is the elastic strain to yield,
eff
( p )
is the effective plastic strain, k is the strength coefficient, and n is the hardening coefficient. The stressstrain behavior can be specified at only one temperature. Input the elastic modulus (EX), density (DENS), and
Poissons ratio (NUXY) with the MP command. Input the strength coefficient, hardening coefficient, and strain
rate parameters C and P with TB,PLAW,,,,2 and locations 1 - 4 of the TBDATA command:
TB, PLAW,,,,2
TBDATA, 1, k (strength coefficient)
TBDATA, 2, n (hardening coefficient)
TBDATA, 3, C (strain rate parameter)
TBDATA, 4, P (strain rate parameter)
For an example input listing, see Section B.2.18: Powerlaw Plasticity Example: Aluminum 1100.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
719
Vm
where 0 is the initial yield stress, peff is the effective plastic strain, and the other parameters are as described in
the TBDATA input shown below.
Strain rate is accounted for using the Cowper and Symonds model (where C and P are the Cowper-Symonds
strain rate parameters), which scales the yield stress with the factor:
p
1 + eff
C
To convert from these constants, set the viscoelastic constants, Vk and Vm, to the following values:
1
1
Vk = p
C
Vm =
1
p
You must input density (DENS) using the MP command. You must also input the elastic modulus (EX), Poisson's
ratio (NUXY), and the coefficient of thermal expansion (ALPX) by either using the MP command or by inputting
load curves in locations 13, 14, and 18 of the TBDATA command. You can use the load curve input to represent
each quantity as a function of temperature. Input the remaining parameters (described below) with the
TB,PLAW,,,,12 command and locations 1 - 18 of the TBDATA command.
TB,PLAW,,,,12
TBDATA,1, 0 (initial yield stress)
TBDATA,2, QR1 (isotropic hardening parameter)
TBDATA,3, CR1 (isotropic hardening parameter)
TBDATA,4, QR2 (isotropic hardening parameter)
TBDATA,5, CR2 (isotropic hardening parameter)
TBDATA,6, QX1 (kinematic hardening parameter)
TBDATA,7, CX1 (kinematic hardening parameter)
TBDATA,8, QX2 (kinematic hardening parameter)
TBDATA,9, CX2 (kinematic hardening parameter)
TBDATA,10, Vk (viscous parameter)
TBDATA,11, Vm (viscous parameter)
TBDATA,12, LCID1 (load curve ID for total true stress vs. effective plastic strain (if specified, TBDATA 2-9 are
ignored))
TBDATA,13, LCID2 (load curve ID for elastic modulus vs. temperature)
TBDATA,14, LCID3 (load curve ID for Poisson's ratio vs. temperature)
TBDATA,15, LCID4 (load curve ID for initial yield stress vs. temperature; if input, 0 is ignored)
720
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
y = 0 + Eh p
You can calculate the plastic hardening modulus, Eh, in terms of the Youngs modulus and tangent modulus:
Eh =
E tE
E Et
If the effective true plastic strain and stress values are specified, the stress-strain behavior is defined by the data
points along the effective true stress vs. true plastic strain curve. Up to 16 data points may be specified. Linear
extrapolation will be used if the strain values exceed the maximum input values; therefore, you should input
values to cover the full range of strains expected in the analysis. Input the density (DENS), elastic modulus (EX),
and shear modulus (GXY) with the MP command. Input the parameters described below with TB,PLAW,,,,9 and
locations 1-45 of the TBDATA command:
TB,PLAW,,,,9
TBDATA,1, 0 (initial yield stress)
TBDATA,2,Eh (hardening modulus)
TBDATA,3,PC (pressure cutoff value)
TBDATA,4, f (failure strain)
TBDATA,5-20,1-16 (effective strain data curve values)
TBDATA,21-36,1-16 (effective stress data curve values)
TBDATA,37,Co (linear polynomial equation of state constant)
TBDATA,38,C1 (linear polynomial equation of state constant)
TBDATA,39,C2 (linear polynomial equation of state constant)
TBDATA,40,C3 (linear polynomial equation of state constant)
TBDATA,41,C4 (linear polynomial equation of state constant)
TBDATA,42,C5 (linear polynomial equation of state constant)
TBDATA,43,C6 (linear polynomial equation of state constant)
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
721
1
Sij Sij
2
and is the objective scalar measure of the distortional or shearing stress. The first invariant of the stress, J1, is the
trace of the stress tensor.
$ "
%#!
"
2D
1 space, as shown in Figure 7.1: Surface of the
The cap model consists of three surfaces in pressure
Two-invariant Cap Model. Surface f1 is the failure envelope, f2 is the cap surface, and f3 is the tension cutoff. The
functional form of f1 is
f1 =
where Fe is given by
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Tmises X ( n ) L ( n )
J
2D
1 space, and therefore does not harden unless kinematic
This failure envelope surface is fixed in
hardening is present. Next, there is a cap surface f2 in the figure, with f2 given by
f2 = J2D Fc (J1 , )
where Fc is defined by
Fc (J1 , )
1
R
[X ( ) L ( )]2 J1 L ( ) 2
X( ) = + RFe ( )
and L() is defined by
if > 0
L ( )
0 if 0
The hardening parameter is related to the plastic volume change through the hardening law
p
p = W {1 exp[ D ( X( ) X 0 )]}
Geometrically, is seen in the figure as the J1 coordinate of the intersection of the cap surface and the failure
surface. Finally, there is the tension cutoff surface, denoted f3 in the figure. The function f3 is given by
f3 T J1
where T is the input material parameter which specifies the maximum hydrostatic tension sustainable by the
2D
1 space is then bounded by the failure envelope surface above, the
material. The elastic domain in
tension cutoff surface on the left, and the cap surface on the right.
Input the density (DENS) and shear modulus (GXY) with the MP command. Input the parameters described below
with the TB,GCAP command and locations 113 of the TBDATA command.
TB,GCAP
TBDATA,1,K (bulk modulus)
TBDATA,2, (failure envelope parameter)
TBDATA,3, (failure envelope linear coefficient)
TBDATA,4, (failure envelope exponential coefficient)
TBDATA,5, (failure envelope exponent)
TBDATA,6,R (cap surface axis ratio)
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
723
where (ij)sk is the skeletal stress, po is the initial foam pressure, is the ratio of foam to polymer density, ij is the
Kronecker delta, and is the volumetric strain which is defined by:
= V 1 + 0
where V is the relative volume and o is the initial volumetric strain. The yield condition is applied to the principal
trial stresses and is defined by:
Y = a + b (1 + c )
where a, b, and c are user defined input constants. The stress-strain behavior can be specified at only one temperature. Input the elastic modulus (EX) and density (DENS) with the MP command. Poissons ratio for this
model is assumed to be zero. Input the yield stress constants a, b, and c, the initial foam pressure, the ratio of
foam to polymer density, and the initial volumetric strain with TB,FOAM,,,,1 and locations 1 - 6 of the TBDATA
command:
TB, FOAM,,,,1
TBDATA, 1, a
TBDATA, 2, b
TBDATA, 3, c
TBDATA, 4, po (initial foam pressure)
TBDATA, 5, (ratio of foam to polymer density)
TBDATA, 6, o (initial volumetric strain)
724
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
absorbs energy. The elastic stiffness, E , and initial viscous coefficient, V , are both nonlinear functions of the
relative volume, V:
E1 = E1V n1
V2 = V2 1 V n 2
where E1 is the initial elastic stiffness, V2 is the initial viscous coefficient, and n1 and n2 are the Power Laws for
the elastic stiffness and viscous coefficient, respectively. The stress-strain behavior can be specified at only one
temperature. Input the initial elastic stiffness (EX), Poissons ratio (NUXY), and density (DENS) with the MP command. Input the Power Law for the elastic stiffness, the initial viscous coefficient, the elastic stiffness for viscosity
(required to prevent time step problems), and the pearly for viscosity with TB,FOAM,,,,3 and locations 1 - 4 of
the TBDATA command:
TB, FOAM,,,,3
TBDATA, 1, n1 (Power Law for the elastic stiffness)
TBDATA, 2, V2 (initial viscous coefficient)
TBDATA, 3, E1 (elastic stiffness for viscosity)
TBDATA, 4, n2 (Power Law for the viscous coefficient)
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
725
n +1 = n + E ij
ij
ij
tn +1/ 2
where ij is the strain rate, E is the elastic modulus, and t is time. The model includes a tensile stress cutoff value
which defines failure under tensile loads. For stresses below the tensile cutoff value, the model predicts similar
response between tensile and compressive loading. It is important to have a nonzero value for the cutoff stress
to prevent deterioration of the material under small tensile loads. Input the elastic modulus (EX), density (DENS),
and Poissons ratio (NUXY) with the MP command. Input the load curve ID for defining stress vs. volumetric strain,
the tension cutoff value, and the viscous damping coefficient with TB,FOAM,,,,4 and locations 1 - 3 of the TBDATA
command:
TB, FOAM,,,,4
TBDATA, 1, LCID (load curve ID for defining yield stress vs. volumetric strain)
TBDATA, 2, TC (tension cutoff value)
TBDATA, 3, DAMP (viscous damping coefficient, values between 0.05 and 0.5 are recommended)
G=
E
= Elastic shear modulus for fully compacted honeycomb material
o
2(1 + )
and
1 V
= max min
,1, 0
1 Vf
726
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
727
< 0, =
1
0
Input the required constants using the TBDATA command. The starting location (LOC) for the constants depends
on which material model you are using. (See the description for the specific material model to find out what
starting location value should be used.)
TBDATA,LOC,C0
TBDATA,LOC+1,C1
TBDATA,LOC+2,C2
TBDATA,LOC+3,C3
TBDATA,LOC+4,C4
TBDATA,LOC+5,C5
TBDATA,LOC+6,C6
TBDATA,LOC+7,E0 (initial internal energy)
TBDATA,LOC+8,V0 (initial relative volume)
p=
a
0C2 1 + 1 0 2
2
2
2
3
S3
1 (S1 1) S2
+1
( + 1)2
+ ( 0 + a )E
p = 0C2 + ( 0 + a )E
where C is the intercept of the vs-vp curve; S1, S2, and S3 are the coefficients of the slope of the vs-vp curve, 0
=
is the Gruneisen gamma, a is the first order volume correction to 0, and
1
0
Input the required constants using the TBDATA command. The starting location (LOC) for the constants depends
on which material model you are using. (See the description for the specific material model to find out what
starting location value should be used.)
TBDATA,LOC,C
TBDATA,LOC+1,S1
728
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
P = Ci ( vi ) + Ti ( vi )E
Input the required constants using the TBDATA command. The starting location (LOC) for the constants depends
on which material model you are using. (See the description for the specific material model to find out what
starting location value should be used.)
TBDATA,LOC,GAMA (temperature constant)
TBDATA,LOC+1,E0 (initial internal energy)
TBDATA,LOC+2,V0 (initial relative volume)
TBDATA,LOC+3 to LOC+12,v1-v10 (volumetric strain data values 1-10, natural log of the relative volume)
TBDATA,LOC+13 to LOC+22,C1-C10 (volumetric pressure values at vi)
TBDATA,LOC+23 to LOC+32,T1-T10 (temperature values at vi)
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
729
y = ( A + B p )(1 + c ln *)(1 T *m )
where:
A, B, C, n, and m = material constant
p = effective plastic strain
p
*=
0
T * = homologous temperature =
T Troom
Tmelt Troom
f = D1 + D2 exp D3 * 1 + D4 ln * 1 + D5 T *
*=
p
eff
D=
p
f
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
(T )
y = C1 + C2 (p )0.5 e( C3 + C4 In (*)) T + C5
(293)
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
731
(T )
y = C1 + C2 e( C3 + C4 In (*)) T + C5 (p )n + C6
(293)
where:
(T )
2
= B1 + B2 T + B3 T
(293)
and the relationship between heat capacity and temperature is given by:
Cp = G1 + G2T + G3T2 + G4T3
An equation of state must also be defined with the Zerilli-Armstrong model.
Input the density (DENS), elastic modulus (EX) and Poissons ratio (NUXY) with the MP command. Input the
parameters described below with TB,EOS and locations 1-19 of the TBDATA command:
TB,EOS,,,,3,EOSOPT
TBDATA,1,o (initial strain)
TBDATA,2,N (strain exponent for BCC metals)
TBDATA,3,Troom (room temperature)
TBDATA,4,PC (pressure cutoff value)
TBDATA,5,SPALL (Spall type. Set SPALL = 1.0 to use minimum pressure limit; SPALL = 2.0 to use maximum
principal stress; SPALL = 3.0 to use minimum pressure cutoff.)
TBDATA,6,C1 (flow stress coefficient)
TBDATA,7,C2 (flow stress coefficient)
TBDATA,8,C3 (flow stress coefficient)
TBDATA,9,C4 (flow stress coefficient)
TBDATA,10,C5 (flow stress coefficient)
TBDATA,11,C6 (flow stress coefficient)
TBDATA,12,FAIL (failure strain)
TBDATA,13,B1 (temperature coefficient)
TBDATA,14,B2 (temperature coefficient)
TBDATA,15,B3 (temperature coefficient)
TBDATA,16,G1 (heat capacity coefficient)
TBDATA,17,G2 (heat capacity coefficient)
TBDATA,18,G3 (heat capacity coefficient)
TBDATA,19,G4 (heat capacity coefficient)
When using the Zerilli-Armstrong Model, you can define one of three types of equations of state: Linear Polynomial (EOSOPT = 1), Gruneisen (EOSOPT = 2), and Tabulated (EOSOPT = 3). Each has its own set of required constants
that are input starting with position (LOC) 20 of the TBDATA command. See the descriptions of the three
equation of state models in Section 7.2.6: Equation of State Models.
732
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
E - Ec
fE /E E
G = G o 1 + bpV 1/3 - h i
- 300 e i m i
3R
where Go, b, h, and f are experimentally determined material constants, p is the pressure, V is the relative volume,
and Ec is the cold compression energy:
x
E c (x) =
pdx
2 ( o a 1 / 2)
x =1V
and Em is the melting energy:
E m (x) = E c ( x ) + 3RTm ( x )
which is in terms of the melting temperature Tm(x):
Tm (x) =
2 ( o a 1/ 3)
R =
R
A
where R is the gas constant and A is the atomic weight. If R is not defined, LS-DYNA computes it with R in the
cm-gram-microsecond system of units.
The yield strength is given by:
Ei - E c
fE / E E
i
m
i
- 300 e
3R
y = 1 + b pV 1/3 - h
o
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
733
is given by:
y = 1 + ( i + p )
o
where i is the initial plastic strain. Whenever o exceeds m, o is set equal to m. After the material melts,
y and G are set to one half their initial values.
If the coefficients EC0-EC9 are not defined in the input, LS-DYNA will fit the cold compression energy to a ten
term polynomial expansion either as a function of or , depending on the input variable FLAG:
E c (i ) =
E c ( i ) =
EC i
i
i=0
9
EC i
i
i=0
1
o
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
735
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
E area
(initial length offset )
The area and offset are defined by means of the real constants for LINK167. For a slack cable, the offset should
be a negative length. For an initial tensile force, the offset should be positive. If a load curve is specified, Youngs
modulus is ignored and the load curve is used instead. The points on the load curve are defined as engineering
stress versus engineering strain, i.e., the change in length over the initial length.
Use the MP and EDMP commands to input the required values:
MP,DENS
MP,EX
EDMP,CABLE,MAT,Load Curve ID
For an example input listing, see Section B.2.26: Cable Material Example: Steel.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
737
No constraints (default)
Constrained X displacement
Constrained Y displacement
Constrained Z displacement
No constraints (default)
Constrained X rotation
Constrained Y rotation
Constrained Z rotation
For example, the command EDMP,RIGID,2,7,7 would constrain the rigid body elements of material 2 in all degrees
of freedom.
738
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
739
740
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
8.3. Loading
As described in Chapter 4, Loading, displacements and velocities are applied to rigid bodies using the EDLOAD
command. It is important to note that all rigid body displacements and velocities are applied to a PART ID rather
than to a component of nodes. For example, a typical EDLOAD command may look like this:
EDLOAD,ADD,RBUX, ,2,PAR1,PAR2
This command defines a UX displacement load on the rigid body identified as PART number 2. The PART number
is input in the Cname field of the EDLOAD command. (Note that for other load types, a component name is input
in this field. However, for the rigid body load labels, a PART number must be input instead of a component name.)
The PART number must correspond to a rigid body that has been defined with the EDMP command.
82
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
83
84
Chapter 9: Hourglassing
Despite being robust for large deformations and saving extensive amounts of computer time, the one-point
(reduced) integration solid and shell elements used in ANSYS LS-DYNA are prone to zero-energy modes. These
modes, commonly referred to as hourglassing modes, are oscillatory in nature and tend to have periods that are
much shorter than those of the overall structural response (i.e., they result in mathematical states that are not
physically possible). They typically have no stiffness and give a zigzag appearance to a mesh (see Figure 9.1: Hourglass Deformations) known as hourglass deformations. The occurrence of hourglass deformations in an analysis
can invalidate the results and should always be minimized.
Chapter 9: Hourglassing
Hourglassing deformations can also be resisted by adding elastic stiffness to a model. Hourglassing can be a
problem with small displacement situations, particularly when dynamic relaxation is used. In these cases, it is
often beneficial to add elastic stiffness to the model instead of using bulk viscosity methods. This can be done
by increasing the hourglassing coefficient (HGCO) of the EDHGLS command. However, use care when increasing
these coefficients because they may over-stiffen the model's response in large deformation problems and cause
instabilities when HGCO exceeds 0.15.
The hourglass control methods discussed so far have been for the entire model. The final method of hourglass
control is to locally increase a portion of the model's stiffness by using the command EDMP,HGLS. For this
command, the material number, the hourglass control type (viscous or stiffness), the hourglass coefficient, and
bulk viscosity coefficients must be specified. (The defaults for hourglass coefficient and bulk viscosities should
be sufficient.) Using this method, hourglass control is specified for a given material and not for the entire model.
This allows resisting hourglass deformations in high risk areas of the model without dramatically changing the
stiffness characteristics of the entire model.
When performing an explicit dynamics analysis with reduced integration elements, it is always important to
determine whether hourglassing effects have significantly degraded the results. As a general guideline, the
hourglassing energy should not exceed 10% of the internal energy. The hourglass energy can be compared to
the internal energy by reviewing the ASCII files GLSTAT and MATSUM and can be plotted in POST26 (see
Chapter 12, Postprocessing). To make sure that hourglass energy results are reported in these files, the HGEN
field of the EDENERGY command must be set to 1.
92
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
l
l
tmin = min = 2
c
c
c=
E
(1 2 ) i
= Poissons ratio
= specific mass density
E = Youngs Modulus
In the ANSYS LS-DYNA program, you can control the minimum time step size by including mass scaling in the
analysis. It may be necessary to use mass scaling if the program calculated time step is too small. When mass
scaling is requested, element density is adjusted to achieve a user specified time step size:
2
t specified
(1 2 ) i i
li
E
i =
( t specified )2 i E
l2 i (1 2 )
i
in element i
The EDCTS command is used to specify mass scaling in ANSYS LS-DYNA. Use this command to apply mass scaling
to the program in one of two ways, depending on the value of the specified time step size DTMS:
If DTMS is positive, the same time step size is used for all elements, and mass scaling is applied to all elements.
If DTMS is negative, mass scaling is applied only to elements whose calculated time step size is smaller
that DTMS.
Of the two methods listed above, the second is typically more efficient and therefore recommended. Although
proper use of mass scaling will add a small amount of mass to the model and slightly change a structure's center
of mass, the CPU reduction it achieves far outweighs the minor errors introduced. For example, when using mass
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
102
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
2.
3.
Solid elements
Beam elements
Shell elements
Discrete elements are intentionally excluded since these elements generally contribute insignificantly to the
calculational costs.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
112
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Element and nodal components should be limited to conserve disk space. For details on creating components,
see the ANSYS Basic Analysis Guide.
Note to LS-DYNA Users--If you are familiar with other ASCII files created by the LS-DYNA program, you may want
to issue EDOUT,ALL to write all possible ASCII files.
122
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
This command method reads in all results data and animates the plot. For faster plotting, use the INRES command
to select specific types of solution data before issuing the above commands:
INRES,Item
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
123
With SOLID164, stress and strain solutions are saved only at the element centroid, regardless of whether
you use 1-point or 8-point integration (KEYOPT(1)).
With BEAM161 and SHELL163 elements, stress and strain solutions are saved only at the element center
(for each layer for shell elements), regardless of whether you use full or reduced integration.
When using BEAM161 and SHELL163 elements, you must specify both the number of integration points
for which you want ANSYS LS-DYNA to calculate the results (using real constants) and the number of these
points for which you want results stored (using the EDINT command). By default, the number of integration
points (NIP) is set to 2 for shell elements; you can specify up to 5 points (layers) for the Gauss integration
rule. For more than 5 layers, you need to use either the trapezoidal integration rule or the user-defined
integration rule. (Note--The trapezoidal rule is not recommended for less than twenty integration points,
especially in bending.) Also by default, the number of layers you specify using EDINT,SHELLIP is set to 3
for shell elements. If NIP = 2, the results at the bottom correspond to integration point 1, the results at
the top correspond to integration point 2, and the results at the middle are an average of the top and
bottom results. For beam elements, results are saved for 4 integration points by default (BEAMIP = 4 on
the EDINT command). For the resultant beam formulation (KEYOPT(1) = 2), there is no stress output, regardless of the BEAMIP setting. If you set BEAMIP = 0, no stress output is written for any of the beam elements. In this case, the beams will not appear in any POST1 plots because the program assumes they are
failed elements.
With SHELL163 elements use the LAYER,NUM command to specify the layer for which you want results.
Layers are numbered starting from the bottom layer and moving up for stress data. However, for strain
data, layer 1 is the bottom layer and layer 2 is the top layer, no matter how many layers exist. Strain information is saved only at these two layers. Further, data is only available at the middle of the layer, not at the
top or bottom surface of a layer (or the element, for that matter). To get results close to the element surface,
specify a large number of integration points through the thickness of the shell. However, saving data for
all of these layers can become expensive.
SHELL163 is not affected by the SHELL,LOC command. By default, output data for the top layer is displayed
when plotting data (PLNSOL, PLESOL), but data for the top and bottom layers appears on printed results.
124
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Caution: The element attributes are not saved on the RSnn file. Therefore, care must be taken when
using select logic during postprocessing with these files. Also, you should not attempt to perform a
solution using only the information contained in these files, nor should you attempt to save a new
database from these files.
You should not use the standard method for animating results (ANDATA command) after adaptive meshing
because it is only valid for animating a single results file (Jobname.RST). Instead, you must use the ANMRES
command to animate results when multiple results files are present (Jobname.RS01, Jobname.RS02, etc.).
Before starting the animation, you must use the /FILNAME command to specify the jobname that the results
will be animated across. Because results are animated across multiple files with the same jobname, simply specifying the FILE command will not allow multiple results animations.
The following sample input demonstrates a typical animation session using ANMRES. In this example, the jobname
is 'impact' and the animation runs from results file impact.rs01 to impact.rs10. The results file increment is two
and the time delay is 0.5 seconds.
/FILNAME,impact
/POST1
FILE,impact,rs01
SET,last
PLNSOL,s,eqv
ANMRES,.5,1,10,2
In the GUI, the ANMRES command is accessed from the utility menu: Utility Menu> PlotCtrls> Animate> Over
Results. Choose the options to animate over multiple results files and then use the push-button animation
controls in the dialog box to start, stop, resume, and otherwise control the animation.
Miscellaneous output data, such as different energies, work, and reaction forces, stored in ASCII files
(GLSTAT, RCFORC, SLEOUT, MATSUM, SPCFORC, NODOUT, and RBDOUT)
Remember when reviewing results in POST26 that you must use the Jobname.HIS file to have adequate results
for a useful time-history record. The Jobname.HIS file must be loaded into POST26 using the FILE command
(Main Menu> TimeHist Postpro> Settings> File):
FILE,Jobname,HIS
If you do not load the Jobname.HIS file into POST26, the results will be read from the Jobname.RST file, which
may provide inadequate data for time-history postprocessing. (This is because the EDHTIME setting is usually
much larger than the EDRST setting.)
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
125
ESOL
STORE
PLVAR
PRVAR
Element and nodal solution results can be obtained from both the Jobname.HIS and Jobname.RST. Care must
be taken when reading from the .HIS file, however, such that only nodes and elements included in the component
specified by the EDHIST,NCOMP are specified.
For storing and displaying solutions at different integration points for the BEAM161 and SHELL163 elements,
use the LAYERP26 (Main Menu> TimeHist Postpro> Define Variables) and the STORE command (Main Menu>
TimeHist Postpro> Store Data) for each layer, as shown below for SHELL163:
LAYERP26,num1
ESOL
NSOL
STORE,MERGE
LAYERP26,num2
ESOL
NSOL
STORE,MERGE
.
.
.
LAYERP26,numn
ESOL
NSOL
STORE,MERGE
You must specify STORE,MERGE for each layer, or your final plot will contain data from the last layer only.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
127
128
To complete an analysis that was terminated or that requires more than the user-defined CPU limit.
To perform an analysis in stages and monitor results at the end of each stage.
The ability to perform an analysis in stages provides great flexibility for explicit dynamic applications. After the
completion of each stage, a restart "dump" file is written. This file contains all information required to continue
the analysis. Results can be checked at each stage by postprocessing the output in the normal way. You may
then modify the model and continue the analysis. For example, you may decide to delete excessively distorted
elements, materials that are no longer important, or contact surfaces that are no longer needed. You may also
choose to change loading or consider other materials that were not present in the previous analysis. You can
also change output frequencies of the various results files. By continuously monitoring and fine-tuning the
analysis, you can increase efficiency and reduce the chance of wasting computer time on an incorrect analysis.
Restarting can also be used to diagnose an analysis that encountered problems. You can restart at a point before
the problem (numerical difficulties or error message) occurred, and request more frequent outputs to the results
files. By viewing the development of the error, you may be able to identify where the first symptoms appeared
and what caused them.
The following restart topics are available:
13.1. The Restart Dump File
13.2. The EDSTART Command
13.3. Effect on Output Files
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Only the commands mentioned above can be used in a small restart analysis (refer to the ANSYS Commands
Reference for details on their usage). Because some of these commands are applicable to a restart as well as a
new analysis, it is important that you issue EDSTART,2 first so that subsequent commands are processed correctly
for the restart.
In a small restart, you should generally extend the time of the calculation (TIME command). If the previous analysis (new or restart) finished at the specified end-time and no new time is input for the subsequent restart, the
restart analysis will stop immediately with only one step. You may also need to modify termination criteria that
were set in the previous analysis using the EDTERM command. If the previous analysis terminated due to one
132
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
2.
3.
4.
5.
6.
7.
133
If you anticipate the need to delete certain elements from the model in a future full restart analysis, use
a different element type number, material number, or real constant number for those elements in the
original analysis, even if they share the same attributes with other elements in the model. This will cause
a unique part number to be assigned to those elements so that they may be deleted later without altering
the makeup of any other parts in the model.
If you need to add more elements to the model during the full restart, use a different element type number,
material number, or real constant number for those elements, even if they share the same attributes with
other elements in the full restart analysis. Again, this will cause a unique part number to be assigned to
the new elements, and previously existing parts will not be altered.
If you do not follow the above guidelines, you may inadvertently create parts in the restart analysis that do not
match parts from the previous analysis. In this case, stress initialization may fail for those parts.
For discrete elements (COMBI165), the initialization is all or nothing. If you initialize any discrete elements in
the full restart, all discrete elements will be initialized.
Although you can change almost any aspect of the analysis in a full restart, a few features are not supported or
are supported in a limited way, as described below.
Contact specifications: You cannot add or delete any contact specifications (EDCGEN and EDDC) in the
full restart. However, you can list contact specifications (EDCLIST) that were defined in the previous analysis.
Initial velocities: You cannot change initial velocities (EDVEL and EDPVEL) in the full restart. For any portions
of the model that are carried over from the previous analysis, the velocities at the beginning of the full
restart will be the same as the velocities at the end of the previous analysis. You cannot define initial velocities for new nodes or parts added in the full restart; the initial velocity of new model entities are assumed
134
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Adaptive meshing: Adaptive meshing (EDADAPT and EDCADAPT) is not supported in a full restart. In
addition, a full restart is not possible if adaptive meshing was used in the previous analysis.
Mass scaling: As in a small restart, mass scaling (EDCTS) is only supported in a full restart if it was active
in the previous analysis.
You can perform multiple full restart analyses and you can mix full restart analyses with other restart analysis
types (simple or small). Just issue EDSTART,3 with a different d3dump file at the beginning of each full restart.
Each time EDSTART,3 is issued, the jobname will be changed from the current jobname to Jobname_nn (nn =
01, 02,) automatically.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
135
136
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Run the explicit analysis as described earlier, using Jobname1. You must use SHELL163 and/or SOLID164
to model the working piece in order to analyze the springback effect in a subsequent ANSYS implicit
analysis. In addition, you must use one of the following element formulations for the SHELL163 elements:
KEYOPT(1) = 1, 2, 6, 7, 8, 9, 10, 11, or 12 (10 and 12 are recommended). Solve and finish the analysis.
You should always check your explicit analysis solution from ANSYS LS-DYNA carefully before proceeding
with the ANSYS implicit analysis. Specifically, check whether there is any undesirable dynamic effect left
in the structure at the end of the explicit run (using POST26).
2.
3.
Change to Jobname2 to prevent the explicit results files from being overwritten.
Command(s): /FILNAME,Jobname2
GUI: Utility Menu> File> Change Jobname
4.
5.
Convert explicit element types to corresponding companion implicit element types. (Note that the 2-D
explicit element, PLANE162, and composite shell element SHELL163 with KEYOPT(3) = 1, cannot be used
in a sequential solution.) The companion explicit-implicit element type pairs are:
Explicit Element Type
LINK160
LINK8
BEAM161
BEAM4
SHELL163
SHELL181
SOLID164
SOLID185
COMBI165
COMBIN14
MASS166
MASS21
LINK167
LINK10
Although all explicit element types are converted, only SHELL163 and SOLID164 data (stresses and
thicknesses for SHELL163, and stresses for SOLID164) are transferred to SHELL181 and SOLID185 (via the
RIMPORT command; see step 12).
Command(s): ETCHG,ETI
GUI: Main Menu> Preprocessor> Elem Type> Switch Elem Type
6.
142
Redefine the key options, real constants, material properties, boundary conditions, and loading values
on any implicit elements that are converted from explicit element types. (For SHELL163 elements that
were converted to SHELL181, you do NOT need to redefine the real constants, but you do need to redefine
the other values. For SOLID164 elements that were converted to SOLID185 elements, you need to specify
the uniform reduced integration option by setting KEYOPT(2) = 1.) The TYPE, REAL, and MAT numbers
from the explicit elements are retained, but the actual key option and real constant values are reset to
zero or the default settings.
Command(s): KEYOPT, R, MP, etc.
GUI: Main Menu> Preprocessor> Element Type/Real Constants/Material Properties/Loads
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Turn off shape checking because elements may have undergone considerable deformation during the
explicit analysis.
Command(s): SHPP,OFF
GUI: Main Menu> Preprocessor> Checking Ctrls> Shape Checking
8.
9.
Unselect or delete any unnecessary elements (mainly those making up any rigid bodies from the explicit
analysis), or convert them to null elements. Any explicit elements that are not either unselected, deleted,
converted to null elements, or converted to implicit will remain active in ANSYS, which will produce an
error and terminate the analysis. Also, if the rigid bodies in the explicit analysis were made up of SHELL163
or SOLID164 elements, these elements must be unselected, deleted, or converted to NULL elements
before importing stresses and thicknesses (from SHELL163 to SHELL181, or SOLID164 to SOLID185) by
the RIMPORT command (See Step 12); otherwise, the implicit analysis will be terminated.
Command(s): ESEL, EDELE
GUI: Utility Menu> Select> Entities or
Main Menu> Preprocessor> Delete> Elements
143
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
nsle,s,1
shpp,off
finish
/solution
nlgeom,on
! Turn on nonlinear geometry
rimport,dyna,,,,,stamp1,rst ! Import stresses and thicknesses from
! explicit analysis
save
solve
finish
Dynamic effects due to a high punch speed used in the explicit (LS-DYNA) run.
Use a lower punch speed in the explicit run. You can do this by redefining the tool speed in the explicit
analysis (EDLOAD command).
Slenderness of the structure; the ratio of part size to its thickness is large.
A more refined mesh may reduce this slenderness effect.
The model contains local instabilities: the part, being very thin, tends to wrinkle along edges or sharp
curvatures.
Use a refined mesh to accurately model curvatures of the formed part.
If you have checked all of the above conditions and the analysis is still not converging, try the following methods
(in the order listed):
144
Use a linear springback analysis (if such an approximation is acceptable) with the command NLGEOM,OFF.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
145
146
In an implicit-to-explicit sequential solution, you must first run an ANSYS implicit structural analysis to apply a
preload to the structure being analyzed. In this implicit analysis, completely constrain all of the nodes of any
elements that will only be used in the explicit analysis (e.g., the bird in a bird-strike problem). The nodal displacements and rotations from the ANSYS implicit solution are written to the ANSYS LS-DYNA dynamic relaxation file
drelax.
Note Temperatures from the ANSYS implicit structural solution are also written to the drelax file, but
are not used by LS-DYNA. See Section 15.1.1: Special Considerations for Thermal Loading for more information on how temperature loads are handled.
After defining additional loads, initial velocities, different material models (e.g., adding plasticity), etc., the explicit
dynamic analysis can be conducted. The first part of this analysis uses the displacement results stored in the
drelax file to do a stress initialization to a prescribed geometry. This preload is applied in pseudo time over 101
time steps to damp out any kinetic energy. The transient portion of the analysis then begins at time zero with a
stable preloaded structure.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Run the implicit analysis as described earlier, using Jobname1. Keep in mind that this analysis must be
small strain with linear material behavior. The only element types that can be used for an implicit-toexplicit sequential solution are:
LINK8
BEAM4
SHELL181
SOLID185
COMBIN14
MASS21
LINK10
2.
Define any additional nodes and elements that are necessary to complete the explicit solution (for example, the bird in a bird-strike simulation, or a rigid surface that a phone would impact in a droptest).
These additional nodes and elements may not be part of the implicit analysis, but they need to be defined
here nonetheless. These additional nodes must be constrained (using D,ALL,ALL,0).
Command(s): N, E
GUI: Main Menu> Preprocessor> Modeling> Create> Nodes or Elements
3.
4.
5.
Change to Jobname2 to prevent the implicit results files from being overwritten.
Command(s): /FILNAME,Jobname2
GUI: Utility Menu> File> Change Jobname
6.
7.
Convert implicit element types to corresponding companion explicit element types. Note that the 2-D
explicit element, PLANE162, cannot be used in this type of sequential solution. (PLANE162 is allowed in
a thermal implicit-to-explicit sequential solution; see Section 15.2: Thermal Implicit-to-Explicit Solution
for details.) The corresponding companion implicit-explicit element type pairs are:
Implicit Element Type
LINK8
LINK160
BEAM4
BEAM161
SHELL181
SHELL163
SOLID185
SOLID164
COMBIN14
COMBI165
MASS21
152
MASS166
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
LINK10
LINK167
Command(s): ETCHG,ITE
GUI: Main Menu> Preprocessor> Element Type> Switch Elem Type
Implicit elements not listed above can also be used, as long as they are defined by the same number of
nodes, but they will not automatically be converted to explicit elements when ETCHG is issued. These
elements must be converted manually using EMODIF. Higher-order implicit elements can also be used,
but must also be converted manually using EMODIF with the corner nodes only. Do NOT delete or unselect
the midside nodes - these nodes must be written to the LS-DYNA input file. The drelax file contains
solutions for these nodes, but the ANSYS LS-DYNA explicit elements do not use these nodes in their
definition.
Command(s): EMODIF
GUI: Main Menu> Preprocessor> Modeling> Move/Modify> Nodes
Note Element types LINK8 and LINK10 lack a third node; however, their corresponding companion explicit element types, LINK160 and LINK167, require a third (orientation) node. If you
are using element types LINK8 or LINK10, you must first convert the element type using ETCHG,ITE,
and then manually define the third node of LINK160 or LINK167 elements using N and EMODIF.
Also, if you are converting BEAM4 to BEAM161, you may need to manually define the third node
of BEAM161 elements as well. However, BEAM4 allows you to define a third, optional node. If
you have defined this third node on BEAM4, then the conversion to BEAM161 will be completed
automatically when you issue ETCHG,ITE. If you did not define the third node on BEAM4, then
you must manually define it on BEAM161 using N and EMODIF.
8.
Redefine the key options, real constants, boundary conditions, and loading values on the explicit elements.
The TYPE, REAL, and MAT numbers from the implicit elements are retained, but the actual key option
and real constant values are reset to zero or the default settings.
Command(s): KEYOPT, R, MP, etc.
GUI: Main Menu> Preprocessor> Element Type, Real Constants, Material Props, or LS-DYNA
Options
9.
Remove constraints from the additional nodes or elements defined in Step 2, above.
Command(s): DDELE
GUI: Main Menu> Preprocessor> LS-DYNA Options> Constraints> Delete
153
!
!
!
!
!
!
!
!
!
!
!
!
!
!
d,all,uy,0
finish
/solution
rexport,dyna,,,,,drop1,rst ! Create DRELAX file from implicit
! results
eddrelax,ansys
! Specify stress initialization by prescribed
! geometry
edpart,create
! Create parts for loading
edpart,list
! List parts
edvel,
! Apply initial velocities
edload,add,rbvx,,2,time,load,0 ! Apply phase = 0 loads
save
solve
finish
154
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Run the implicit thermal analysis using Jobname1. The only element types that can be used for a thermal
implicit-to-explicit sequential solution are:
PLANE55
SOLID70
SHELL57
Keep in mind that you cannot combine 2-D and 3-D elements in the subsequent explicit analysis.
Therefore, you should use only 2-D or only 3-D elements in the thermal analysis phase.
2.
You may define any additional nodes and elements that are necessary to complete the explicit solution.
These additional nodes and elements may not be part of the thermal analysis, but they can be defined
here if desired. These additional nodes must be constrained to a constant temperature (using
D,NODE,TEMP,VALUE).
Command(s): N, E
GUI: Main Menu> Preprocessor> Modeling> Create> Nodes or Elements
3.
4.
5.
Change to Jobname2 to prevent the thermal results files from being overwritten.
Command(s): /FILNAME,Jobname2
GUI: Utility Menu> File> Change Jobname
6.
155
Convert thermal element types to corresponding companion explicit element types. The corresponding
companion thermal-explicit element type pairs are:
Thermal (Implicit) Element Type
PLANE55
PLANE162
SHELL57
SHELL163
SOLID70
SOLID164
Command(s): ETCHG,TTE
GUI: Main Menu> Preprocessor> Element Type> Switch Elem Type
8.
Redefine the key options, real constants, and material property values on the explicit elements. The TYPE,
REAL, and MAT numbers from the implicit elements are retained, but the actual key option and real
constant values are reset to zero or the default settings. Keep in mind that for temperature loading to
take effect, you must use the Temperature Dependent Bilinear Isotropic material model in the explicit
analysis.
Command(s): KEYOPT, R, MP, etc.
GUI: Main Menu> Preprocessor> Element Type> Real Constants> Material Properties
9.
Define additional nodes and elements, if needed. You may have already completed this task in the
thermal analysis portion of the sequential solution (see Step 2). However, if desired, you could wait until
the explicit analysis to complete this step.
Command(s): N, E
GUI: Main Menu> Preprocessor> Modeling> Create> Nodes or Elements
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
nsel,s,loc,z,0
d,all,temp,300.0
nsel,s,loc,z,11
cm,ntip,node
esln
cm,etip,elem
bf,all,hgen,20.0e3
esel,all
nsel,s,ext
sf,all,conv,6.0,0.0
nsel,all
save
solve
fini
/post1
set,last
plnsol,temp
/wait,3
fini
/filnam,explicit
/prep7
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
157
etchg,tte
mp,dens,1,8030.0
edmp,hgls,1,5
158
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
159
1510
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
In a Eulerian based formulation, the mesh is stationary and the material flows through the mesh. The Eulerian
approach originated in the fluid dynamics field and is best suited for very large deformation flow problems, such
as the channel flow depicted in Figure 16.3: Eulerian Channel Flow Solution. In the Eulerian method, new
(damaged) surfaces are automatically created. The greatest disadvantage of the Eulerian approach is that a fine
mesh is required to capture the material response, making the method very computationally expensive. This is
particularly true for problems that contain regions where the structural response is desired and the strains are
relatively small.
The Arbitrary Lagrangian-Eulerian (ALE) approach is a very effective alternative for simulating large deformation
problems. In its most basic sense, the ALE method defines that the mesh motion is independent of the motion
of the material being analyzed. Although the mesh motion may be arbitrary, it typically deforms with the material in 'near Lagrangian' flow fields. The greatest advantage of the ALE method is that it allows smoothing of a
162
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
For more information about the ALE formulation, refer to Livermore Software Technology Corporations LS-DYNA
Theoretical Manual.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
163
for finite difference or finite element calculations by using the solutions of Laplace equations as the mesh
lines. The same method can be used to smooth selected points in an unstructured 3-D mesh provided that
it is at least locally structured.
EFAC - Equilibrium smoothing weight factor (available only for PLANE162). For this method, artificial springs
are attached to each ALE element node. The springs are used to adjust the position of each node from the
equilibrium solution. This approach can overcome possible calculation instabilities found in the other
smoothing methods.
The EDALE command also allows you to define start and end times for ALE smoothing.
You can define two additional ALE options with the EDGCALE command:
ADV - Number of cycles between advection
METH - Advection method (donor cell or Van Leer)
In general, it is not worthwhile to advect an element unless at least twenty percent of its volume will be transported, because the gain in the time step size will not offset the cost of the advection calculations. It is best to begin
an ALE analysis with a Van Leer advection technique (METH = 1).
If you are working in the GUI, all ALE related options are accessed by picking Main Menu> Solution> Analysis
Options> ALE Options.
Note In order to activate the ALE formulation for a given set of elements, you must specify a smoothing
weight factor (EDALE command), and you must specify the cycles between advections (NADV argument
on the EDGCALE command). If you do not supply both of these inputs, the default Lagrangian formulation
will be used for all elements, including those for which KEYOPT(5) = 1.
The remap step maps the solution from a distorted Lagrangian mesh onto the new mesh. The underlying assumptions of the remap step are 1) the topology of the mesh is fixed (that is, the element nodal connectivity
remains unchanged), and 2) the mesh motion during a step is less than the characteristic lengths of the surrounding elements. The algorithms for performing the remap step are taken from the computational fluid dynamics
community and are referred to as advection algorithms.
The donor cell algorithm is a first order Godunov method applied to the advection equation. Aside from its first
order accuracy, it is a good advection algorithm; it is stable, monotonic and simple. The Van Leer algorithm is a
higher order Godunov method that improves the estimates of the initial values of left and right states for the
Riemann problem at the nodes. The donor cell algorithm assumes that the distribution of the initial value function
is constant over an element. Van Leer replaces the piecewise constant distribution with a higher order interpolation function that is subject to an element level conservation constraint.
164
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
If you are starting ANSYS from a command line (UNIX) or DOS prompt (Windows), include -dtm in the line
as well as a productname that represents an ANSYS LS-DYNA product. For example, to start the ANSYS
LS-DYNA product with the DTM in interactive, graphics mode, type the following in a command line:
ansys100 -p dyna -dtm -g
If you are starting ANSYS using the ANSYS Launcher, on the Launch tab choose an ANSYS LS-DYNA license
and select the Drop Test Module for ANSYS LS-DYNA add-on.
Refer to Chapter 3, Running the ANSYS Program in the ANSYS Operations Guide for more information on starting
ANSYS.
You can also have ANSYS automatically include the DTM option by setting the following environment variable:
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
See the ANSYS Installation and Configuration Guide for your platform for more details.
We recommend that a certain degree of damping be applied to the model to reduce the oscillatory response.
You can apply both alpha (mass weighted) and beta (stiffness weighted) damping using the EDDAMP command.
Note After completing the model and before entering the DTM, it is good practice to save the existing
database under a unique name (pick Utility Menu> File> Save As).
172
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
173
2.
Click on the Status tab since this displays your most updated input data.
3.
View the object and the target in the ANSYS Graphics window.
If you choose the Status tab, ANSYS attempts to automatically create the target (unless you have already created
it by choosing the Target tab) and displays a statement in black print at the end of the Status tab page indicating
that all the set-up data is acceptable and you can proceed to solve the analysis, or a statement in yellow print
indicating that you should proceed with caution because the data produced error messages, or a statement
displayed in red indicating that the set-up data is unacceptable.
If you receive either the second or third statement regarding error messages, check them in the ANSYS Output
window. See Section 17.4.5: Status Tab of the Drop Test Set-up Dialog Box in this user's guide for more information
on using the Status tab. Make changes if necessary, then choose OK. At this time, all of the standard ANSYS
menus become usable again. You then initiate the LS-DYNA solution by choosing Main Menu> Drop Test>
Solve.
Note A database that you save from the current version of ANSYS LS-DYNA DTM is not compatible
with a previous version of ANSYS LS-DYNA DTM.
In the DTM, screen coordinates are always defined so that the screen Y axis is in the direction opposite to
the acceleration due to gravity, (g).
174
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
The screen Z axis is along the normal from the screen to the viewer, or "out".
During initialization, a keyword is created that is related to the filtering of menu picks in the DTM. The
keyword is DTFILT. You should avoid using this word as a keyword or a parameter name, as it could interfere
with the filtering of menu options.
Within the DTM, a component containing all the nodes in the object is automatically created and given
the name _dtdrpob. This component name is used in issuing commands associated with LS-DYNA. You
should avoid using this name.
Parameters _dtcgnum and _dtlownum are used for time-history postprocessing plots and lists. You should
avoid using these names.
The DTM assumes that there are no constraints imposed on the object being dropped. However, coupling
of components within an object is allowed.
Some parameter names used internally by the Drop Test Module begin with the underscore character (_),
so you should avoid using parameter names that begin with an underscore. (This is true for any ANSYS
analysis.)
If you do not choose to orient the object, the default assumption is the current screen Y axis view, which
is directed opposite to the acceleration due to gravity (g).
You can modify any input set up parameter any number of times before entering the solution phase. All
parameters in the Drop Test Module are subject to redefinition after being initially defined. You will see
that whenever an input parameter that affects the initial target/object orientation is changed, the target
is recreated automatically to reflect the change when you click on the Target tab or the Status tab.
After you specify a drop height, you will notice that, by default, the target is located nearly in contact with
the object. This is due to the default option to begin the analysis near impact time. The target is created
with just a small separation between the object and target. You may want to adjust the Solution time
option to begin the analysis at drop time. After doing so, you will observe that the redefined target is
at a position consistent with the specified drop height. However, while this option provides full animation
of the drop scenario, it will require more computation time.
To use the option to start the analysis near impact time, the initial angular velocity must be zero.
The Drop Test Module does not alter the model of the object in any way. The DTM does create a model
of a target after you specify a new orientation for the object.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
175
176
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
It is expected that the default target position, size, orientation, material properties, and contact parameters will
be acceptable for most analyses. But if necessary, you can change the characteristics of the target by following
the procedures outlined in the following sections.
177
The material number (MAT) used is one greater than the highest material number used by the model of the object
being dropped. The default property values used for the target depend on the material properties used for the
object. The Drop Test Module searches through all materials belonging to the object, finds the one with the
lowest elastic modulus, and stores the corresponding material number. The modulus, density, and Poissons ratio
corresponding to this material number are also specified as the modulus, density, and Poissons ratio of the target
material.
If you know the actual target material properties, it is preferred that you input them directly. To specify the
properties, choose Main Menu> Drop Test> Set Up, choose the Target tab, then, under the Material Properties
heading, enter the Young's modulus, density, and Poisson's ratio. See Section 17.4.4: Target Tab of the Drop Test
Set-up Dialog Box in this user's guide for more information.
178
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
179
zero (default)
386.4
32.2
9.81
9810
Any value that you choose must have units that are consistent with the units used in the model of the object
being dropped. You only need to input the magnitude of gravity. The direction is assumed to be in the negative
Y direction, referenced to the screen coordinates. The object's coordinate components of the gravity loading
are imposed on the object internally by the DTM, using three separate EDLOAD commands just before the
problem is solved.
Drop Height -- Specify values for the drop height and the height reference point.
Height field: (default = 0) You type the value of the drop height in units consistent with the units of length
used in creating the object being dropped. The height is always measured along the Y axis (referenced
to the screen coordinates) from the center of the top face of the target to some height reference point in
the object.
Reference drop down list: You choose the reference point from the following options:
1710
Lowest Object Point (default): Reference point is the point on the object with the minimum screen
Y coordinate value. If this option is used, then the node on the object with the minimum Y coordinate
value is the reference point. If there is more than one node at this Y coordinate, then one of them is
chosen arbitrarily to be the reference point.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Object CG: Reference point is the object's center of gravity as calculated by ANSYS. If the center of
gravity has been previously calculated for some other purpose, then the information is stored, and
the calculation is not repeated.
If the analysis is set to begin near impact time (default), then you will notice that the separation between the
center of the target and the reference point, measured along the screen Y axis, is not equal to the specified drop
height. This is because the target is created so that there is only a small separation, which is the case just before
impact. However, if the analysis is set to begin at drop time, the separation between the target and the reference
point equals the drop height.
By default, the target is oriented such that its top face is normal to the screen Y axis. If the orientation of the target
is modified from the default orientation (see Section 17.3.2.3: Target Orientation in this User's Guide), the location
of the center of the top face of the target at the time of the drop is not altered. Therefore, the definition of the
drop height outlined above is not changed.
Note You should always define the drop height so that all nodes on the object are above the plane of
the top face of the target prior to the drop.
Set Orientation -- Specify the orientation of the model with respect to the gravitational field.
Rotate button: You click this button to access the DTM Rotate dialog box, which allows you to dynamically
rotate the object on the screen. This dialog box contains a Help button that you can click for more information.
Pick Nodes button: You click this button to access a picking dialog box, which allows you to define a
vector either by picking two nodes, or by picking one node and having the object's center of gravity be
the second node. The up direction is defined as the direction going from the first node that you pick to
the second node that you pick. This dialog box contains a Help button that you can click for more information.
Solution Time -- Specify when the drop test analysis will start and how long it will run.
Start analysis near impact time radio button (default): You click this option to save computation time
that would be expended as the object drops through the air. To use this option, the angular velocity must
be zero. The velocity at the start of the analysis is calculated automatically assuming the object is a rigid
body and there is no resistance due to the air. If an initial translational velocity is applied, it is taken into
account in calculating the velocity at the start of the analysis.
Start analysis at drop time radio button: You click this option to include the entire time it takes for the
object to drop from the specified height in the analysis. Although extra computation time is expended
that is not directly involved with the drop test impact, this option accounts for rigid body rotation and is
useful if you want to produce animations that include frames of the entire drop test, starting from when
the object drops.
Run time after impact field: You type the time in seconds that you want the analysis to run after the
object hits the target.
Number of Results Output -- Specify the number of times that the output is written to the results file and to
the time-history file.
On Results file field: (default = 100) You type the number of times that the output is written to results
file Jobname.RST. See the EDRST command for details.
On time history file field: (default = 1000) You type the number of times that the output is written to
time-history file Jobname.HIS. See the EDHTIME command for details. Note here that a component containing the lowest point and the node closest to the object's center of gravity is automatically created for
time-history postprocessing. You may define additional nodes. Refer to the EDHIST command for details.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
1711
OK: Applies the settings on the Drop Test Set-up dialog box to the ANSYS database and closes the dialog
box.
Cancel: Cancels any changes that you have made on the Drop Test Set-up dialog box since the last time
that you clicked the OK button.
Help: Displays this help topic.
Related Information -- For further information about the options on the Basic tab, see the following sections
in this user's guide:
Use the Velocity tab to specify the object velocities, as described below.
Translational Velocity -- Specify the initial translational velocity of the object being dropped at the beginning
of the drop time, even if the analysis is set to begin near impact. If you apply a translational velocity other than
zero, and choose to start the analysis near impact time, the initial angular velocity must be zero. If these conditions
are not met, ANSYS will automatically start the analysis at drop time, and it will automatically recreate the target.
Velocity relative to drop down list: You choose either object coordinate system or screen coordinate
system (default). If you enter the velocity relative to the screen coordinate system, the components are
resolved into the object coordinate system internally by the DTM.
X Component field: (default = 0) The X component of the translational velocity.
Y Component field: (default = 0) The Y component of the translational velocity.
Z Component field: (default = 0) The Z component of the translational velocity.
Angular Velocity -- Specify the initial angular velocity of the object being dropped. The angular velocity is about
the object centroid in units of radians/(unit time). Applying an angular velocity other than zero is only valid when
1712
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Velocity relative to drop down list: You choose either object coordinate system (default) or screen coordinate system. If you enter the velocity relative to the screen coordinate system, the components are
resolved into the object coordinate system internally by the DTM before being input on the EDVEL
command.
X Component field (default = 0): The X component of the angular velocity.
Y Component field (default = 0): The Y component of the angular velocity.
Z Component field (default = 0): The Z component of the angular velocity.
Action Buttons
OK: Applies the settings on the Drop Test Set-up dialog box to the ANSYS database and closes the dialog
box.
Cancel: Cancels any changes that you have made on the Drop Test Set-up dialog box since the last time
that you clicked the OK button.
Help: Displays this help topic.
Related Information -- For further information about the options on the Velocity tab, see the following sections
in this user's guide:
When you choose the Target tab, the target is automatically created.
Also, use the Target tab to specify the options described below.
Dimensions -- Specify the target's position and size.
Center target using drop down list. You can change the position of the target by changing the location
of the point on the center of the top face of the target in the screen X-Z plane. The point's Y coordinate
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
1713
Lowest Object Point: (default) The lowest object point is the point on the object with the minimum
screen Y coordinate value. If only a single node exists on the object at this point, then the screen X
coordinate value of this node is used as the point's X coordinate, and the screen Z value is used as the
point's Z value. If there is more than one node at the minimum Y point, then the X coordinate of the
point is found by first selecting only those nodes located at the minimum Y point. Then, the point's X
value is equal to the sum of the maximum and minimum screen X values in the selected set of nodes,
divided by 2. Similarly, the point's Z value is equal to the sum of the maximum and minimum screen
Z values in the selected set of nodes, divided by 2.
Object CG: The screen X value at the object center of gravity is used as the point's X value, and the
screen Z value of the cg is used as the point's Z value.
Length scale factor field: (default = 1.0) You type a scale factor that changes the length of the target from
the default length. The default length is determined by ANSYS, based on the overall dimensions of the
object being dropped.
Thickness scale factor field: (default = 1.0) You type a scale factor that changes the thickness of the target
from the default thickness. The default thickness is equal to 0.1 times the default length.
Orientation angle -- Specify the rotation of the target about the screen Z axis (normal to the screen and pointing
toward you).
Rotation about screen Z field: (default = 0) You type the rotation angle in degrees about the screen Z axis.
The drop height loses its meaning when the normal to the targets large faces is perpendicular to the screen Y
axis. Therefore, if you modify the orientation so that the angle is greater than 89o or less than -89o with the screen
Y axis, then the input is considered to be invalid. A warning is printed to the screen, and the target is returned
to its previous orientation.
Material Properties -- Specify the material properties of the target. Material properties are used for a contact
stiffness calculation and a time step calculation.
Young's Modulus field: (default = object's lowest Young's modulus) You type the target's Young's modulus in units consistent with the analysis.
Density field: (default = object's highest density) You type the target's density in units consistent with
the analysis.
Poisson's Ratio field: (default = object's Poisson's ratio) You type the target's Poisson's ratio.
Static friction coeff field: (default = 0.0) You type the static friction coefficient, which is FS as described
in the Notes section of the EDCGEN command.
Dynamic friction coeff field: (default = 0.0) You type the dynamic friction coefficient, which is FD as described in the Notes section of the EDCGEN command.
Exponential decay coeff field: (default = 0.0) You type the exponential decay coefficient, which is DC as
described in the Notes section of the EDCGEN command.
Action Buttons
1714
OK: Applies the settings on the Drop Test Set-up dialog box to the ANSYS database and closes the dialog
box.
Cancel: Cancels any changes that you have made on the Drop Test Set-up dialog box since the last time
that you clicked the OK button.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Related Information -- For further information about the options on the Target tab, see the following sections
in this user's guide, except as noted:
Use the Status tab to check the values currently set for acceleration due to gravity, initial translational and rotational velocity, and the total time at the end of solution. One of the following statements is also displayed:
Successful set-up. Ready to drop. [black print] This statement indicates that all the current set-up data
is acceptable and you can click OK to commit the settings.
Proceed with caution. See error messages. [yellow print] This statement indicates that some of the
current set-up data could be problematic in enabling a solution, or is otherwise non-standard (for example,
a zero entry detected for the value of gravity). You should examine the error messages in the ANSYS
Output window to determine what is causing them to occur. Then, if necessary, change settings as appropriate and recheck the Status tab until the first message in black appears for problem situations. Click OK
to commit the settings.
Unsuccessful set-up. See error messages. [red print] This statement indicates that the current set-up
data has produced errors that are serious enough to prevent a solution. The target will not be created in
this situation. You should examine the error messages in the ANSYS Output window to determine the
causes of these problems. Then change settings as appropriate and recheck the Status tab until the first
message in black appears. Click OK to commit the settings.
It is recommended that you check this information in the Status tab along with a view of the object and the
target in the in the ANSYS Graphics window before you click OK in the Drop Test Set-up dialog box as this
commits all the set up data to the ANSYS database.
Action Buttons
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
1715
OK: Applies the settings on the Drop Test Set-up dialog box to the ANSYS database and closes the dialog
box.
Cancel: Cancels any changes that you have made on the Drop Test Set-up dialog box since the last time
that you clicked the OK button.
Help: Displays this help topic.
Related Information -- For further information about the options on the Status tab, see the following section
in this user's guide:
Pick Two Nodes option: To specify the model's orientation by picking two nodes:
1.
2.
3.
4.
The object's orientation changes according to the vector defined by the two nodes that you picked. The
object is aligned such that the vector you created will be parallel to the Y direction of the screen coordinate
system.
Pick One Node and Use CG option: To specify the model's orientation by picking one node and using
the object's center of gravity (calculated by ANSYS):
1.
2.
3.
Choose OK again in the Pick Nodes dialog box because the second node is the object's center of
gravity by default.
The object's orientation changes according to the vector defined by the node that you picked and the
object's center of gravity.
Related Information -- For further information related to the Pick Nodes dialog box, see the following sections:
Section 17.4.2: Basic Tab of the Drop Test Set-up Dialog Box in this user's guide
Chapter 5, Graphical Picking in the ANSYS Operations Guide
1716
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Range Minimum, Maximum text fields -- Specify the first and last data sets over which the animation should
span. If the Minimum field is left blank or equals 0, the default is the first data set. If the Maximum field is left
blank or equals 0, the default is the last data set.
Increment result set text field -- Specify the increment between result data set (default = 1).
Include last SBST (substep) for each LDST (load step) check box -- Specify whether or not to force the last
result data set in a load step to be animated. The default is unchecked, meaning that the last data set is not included in the animation.
Animation time delay (sec) text field -- Specify the time delay between frames during the animation (default
= 0.5 seconds).
Display Type two-column selection list -- Specify what parameter is to be animated in the left column (for example,
Stress, Stress-total, Stress-plastic), and the specific subset of the chosen parameter in the right column (for example,
1st Principal Stress, von Mises Stress).
Action Buttons
OK: Applies the settings on the Animate Over Results dialog box to the ANSYS database and closes the
dialog box.
Cancel: Cancels any changes that you have made on the Animate Over Results dialog box since the last
time that you clicked the OK button.
Help: Displays information on the ANDATA command from the ANSYS Commands Reference.
After you choose OK, the animation will be created and will be displayed. An Animation Controller dialog box
will also be displayed that includes controls for starting and stopping the animation, as well as changing its direction. This dialog box includes its own Help button where more detailed information is available.
Related Information -- For further information related to the Animate Over Results dialog box, see the following
sections:
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
1717
Note The contents of the Graph Time-History Variables dialog box and the List Time-History
Variables dialog box are the same so this help topic describes both of them. Only the Graph TimeHistory Variables dialog box is shown above.
Use the applicable dialog box to display the time-history results of a drop test analysis as described below:
Nodes to graph/list drop down list: You select one of the following choices as the node(s) to graph or list:
CG and Lowest Pt.: (default) Both the center of gravity and the lowest point of the object, as initially
specified will be graphed or listed.
Center of grav.: Only the center of gravity of the object will be graphed or listed.
Initial Low Pt.: Only the lowest point of the object, as initially specified will be graphed or listed.
Item and Comp to be graphed/listed two-column selection list: You first select Displacement (default), Velocity,
or Acceleration as the item to be graphed or listed, then select the x, y, or z component of the chosen item referenced to the object or screen coordinate system. The default item and component is displacement in the Y
direction referenced to the object coordinate system.
Action Buttons
OK: Applies the settings on the Graph/List Time-History Variables dialog box to the ANSYS database
and closes the dialog box.
Cancel: Cancels any changes that you have made on the Graph/List Time-History Variables dialog box
since the last time that you clicked the OK button.
Help: Displays this help topic.
Related Information -- For further information related to the Graph/List Time-History Variables dialog box,
see the following sections in this user's guide:
1718
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
{ut + t } = [K ]1 {Fta t }
+
For linear problems, the solution is unconditionally stable when the stiffness matrix [K] is linear, and large time
steps can be taken.
For nonlinear problems:
Convergence tools are provided, but convergence is not guaranteed for highly nonlinear problems.
{Ftext }
is the applied external and body force vector
{Ftint }
is the internal force vector which is given by:
Fhg is the hourglass resistance force, and Fcont is the contact force.
The velocities and displacements are then evaluated:
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
t t + t / 2 = .5 ( t t + t t + t )
and
t t t / 2 = .5 ( t t t t + t )
The geometry is updated by adding the displacement increments to the initial geometry {X0}:
{ x t + t } = { xo } + {ut + t }
For nonlinear problems:
The equations become uncoupled and can be solved for directly (explicitly).
No inversion of the stiffness matrix is required. All nonlinearities (including contact) are included in the
internal force vector.
Very small time steps are required to maintain the stability limit.
t t crit =
2
max
A2
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
max = 2
c
l
where
c=
t =
l
c
t is the time needed for the wave to propagate through the rod of length l. Note that the critical time step size
for explicit time integration depends on element length and material properties (sonic speed).
t = 0 . 9
l
c
The characteristic length, l, and the wave propagation velocity, c, are dependent on element type.
For beam elements:
l = length of the element, and
c=
I=
A
max (I1,I 2,I 3,I 4 )
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
A3
I=
2A
max (I1,I 2,I 3 )
c=
A4
E
(1 2 )
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
LS-DYNA Command
LS-DYNA
MAT #
Example
Isotropic Elastic
*MAT_ELASTIC
Yes
Orthotropic Elastic
*MAT_ORTHOTROPIC_ELASTIC
Yes
Anisotropic Elastic
*MAT_ANISOTROPIC_ELASTIC
Yes
Elastic Fluid
*MAT_ELASTIC_FLUID
No
Blatz-Ko Rubber
*MAT_BLATZ-KO_RUBBER
Yes
27
Yes
Viscoelastic
*MAT_VISCOELASTIC
Yes
Bilinear Isotropic
Plasticity
*MAT_PLASTIC_KINEMATIC
Yes
Temperature De*MAT_ELASTIC_PLASTIC_THERMAL
pendent Bilinear Isotropic
No
37
Yes
39
Yes
Bilinear Kinematic
*MAT_PLASTIC_KINEMATIC
Yes
Plastic Kinematic
*MAT_PLASTIC_KINEMATIC
Yes
3 Parameter Barlat
Plasticity
*MAT_3-PARAMETER_BARLAT
36
Yes
Barlat Anisotropic
Plasticity
*MAT_BARLAT_ANISOTROPIC _PLASTICITY
33
Yes
Rate Sensitive
*MAT_RATE_SENSITIVE_POWERLAW _PLASTICITY
Powerlaw Plasticity
64
Yes
19
Yes
Piecewise Linear
Plasticity
24
Yes
*MAT_PIECEWISE_LINEAR_PLASTICITY
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
LS-DYNA Command
LS-DYNA
MAT #
Example
123
Yes
22
No
Concrete Damage
72
No
18
Yes
Elastic Viscoplastic
Thermal
106
Yes
10
No
Geological Cap
*MAT_GEOLOGICAL_CAP_MODEL
25
Yes
*MAT_CLOSED_CELL_FOAM
53
No
Viscous Foam
*MAT_VISCOUS_FOAM
62
No
*MAT_LOW_DENSITY_FOAM
57
No
Crushable Foam
*MAT_CRUSHABLE_FOAM
63
No
Honeycomb Foam
*MAT_HONEYCOMB
26
No
Tabulated EOS
*EOS_TABULATED
Bamman
*MAT_BAMMAN
*MAT_CONCRETE_DAMAGE
*MAT_ELASTIC_VISCOPLASTIC _THERMAL
No
51
No
15
Yes
Johnson-Cook
Gruneisen EOS
15
Yes
Yes
Yes
65
No
11
Yes
N/A
No
General Nonlinear
Spring
*MAT_SPRING_GENERAL_NONLINEAR
N/A
No
Nonlinear Elastic
Spring
*MAT_SPRING_NONLINEAR_ELASTIC
N/A
No
N/A
No
N/A
No
Maxwell Viscosity
Spring
*MAT_SPRING_MAXWELL
N/A
No
Linear Viscosity
Damper
*MAT_DAMPER_VISCOUS
N/A
No
*MAT_JOHNSON_COOK
*EOS_GRUNEISEN
*EOS_GRUNEISEN
Zerilli-Armstrong
B2
*MAT_ZERILLI_ARMSTRONG
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
LS-DYNA Command
LS-DYNA
MAT #
Example
Nonlinear Viscosity
Damper
*MAT_DAMPER_NONLINEAR_VISCOUS
N/A
No
Cable
*MAT_CABLE_DISCRETE_BEAM
71
Yes
Rigid
*MAT_RIGID
20
Yes
! Pa
MP,nuxy,1,.29
! No units
MP,dens,1,7850
! kg/m3
! Pa
MP,ey,1,358.1e9
! Pa
MP,ez,1,358.1e9
! Pa
MP,gxy,126.9e9
! Pa
MP,gxz,126.9e9
! Pa
MP,gyz,126.9e9
! Pa
MP,nuxy,1,.20
! No units
MP,nuxz,1,.20
! No units
MP,nuyz,1,.20
! No units
MP,dens,1,3750
! kg/m3
! kg/m3
TB,ANEL,1
TBDATA,1,121.0e9
! C11 (Pa)
TBDATA,2,48.1e9
! C12 (Pa)
TBDATA,3,121.0e9
! C22 (Pa)
TBDATA,4,44.2e9
! C13 (Pa)
TBDATA,5,44.2e9
! C23 (Pa)
TBDATA,6,51.3e9
! C33 (Pa)
TBDATA,10,18.5
! C44 (Pa)
TBDATA,15,18.5
! C55 (Pa)
TBDATA,21,24.2
! C66 (Pa)
! kg/m3
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
B3
! Pa
! lb/in3
MP,nuxy,1,.499
! No units
TB,MOONEY,1
TBDATA,1,80
! C10 (psi)
TBDATA,2,20
! C01 (psi)
! kg/m3
TB,EVISC,1
TBDATA,46,27.4e9
TBDATA,47,0.0
! Go (Pa)
! G (Pa)
TBDATA,48,60.5e9
TBDATA,61,.53
! 1/
! Pa
MP,nuxy,1,.31
! No units
MP,dens,1,8490
! kg/m3
TB,BISO,1
TBDATA,1,900e6
TBDATA,2,445e6
! Pa
MP,nuxy,1,.29
! No units
MP,dens,1,7845
! kg/m3
TB,PLAW,,,,7
TBDATA,1,128.5e6
TBDATA,2,202e5
TBDATA,3,1.41
! r-value
TBDATA,4,1
*DIM,STRAIN,,5
*DIM,YLDSTRES,,5
Strain(1) = 0,.05,.1,.15,.2
YldStres(1)=207e6,210e6,214e6,218e6,220e6 ! yield stress
B4
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
! Pa
MP,nuxy,1,.29
! No units
MP,dens,1,.00285
! kg/m3
TB,PLAW,1,,,10
TBDATA,1,20e3
TBDATA,2,5000
TBDATA,3,.2
! Hardening parameter
TBDATA,5,1
*DIM,mnstrn,,6
*DIM,mjstrn,,6
mnstrn(1) = -30,-10,0,20,40,50
mjstrn(1) = 80,40,29,39,45,44
EDCURVE,ADD,1,mnstrn (1),mjstrn(1)
! Pa
MP,nuxy,1,.36
! No units
MP,dens,1,4650
! kg/m3
TB,BKIN,1
TBDATA,1,70e6
TBDATA,2,112e6
! Pa
MP,nuxy,1,.27
! No units
MP,dens,1,7865
! kg/m3
TB,PLAW,,,,1
TBDATA,1,310e6
TBDATA,2,763e6
TBDATA,4,40.0
! C (s-1)
TBDATA,5,5.0
!P
TBDATA,6,.75
! Failure strain
! Pa
MP,nuxy,1,.34
! No units
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
B5
! kg/m3
TB,PLAW,,,,3
TBDATA,1,1
TBDATA,2,25e6
TBDATA,3,145e6
TBDATA,4,0.170
! Barlat exponent, m
TBDATA,5, .73
! R00
TBDATA,6,.68
! R45
TBDATA,7,.65
! R90
TBDATA,8,0
! CSID
! Pa
MP,nuxy,1,.34
! No units
MP,dens,1,2720
! kg/m3
TB,PLAW,,,,6
TBDATA,1,1.04
! k (MPa)
TBDATA,2,.65
TBDATA,3,.254
!n
TBDATA,4,11
! Barlat exponent, m
TBDATA,5, 1.017
!a
TBDATA,6,1.023
!b
TBDATA,7,.9761
!c
TBDATA,8,.9861
!f
TBDATA,9,.9861
!g
TBDATA,9,.8875
!h
! Pa
MP,nuxy,1,.33
! No units
MP,dens,1,2750
! kg/m3
TB,PLAW,,,,4
TBDATA,1,1.002
! k (MPa)
TBDATA,2,.7
!m
TBDATA,3,.32
!n
TBDATA,4,5.0
! Pa
MP,nuxy,1,.29
! No units
MP,dens,1,7850
! kg/m3
B6
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
! LCID yield stress vs. strain rate (see first EDCURVE command
below)
TBDATA,2,22e5
TBDATA,3,2
*DIM,StrnRate,,5
*DIM,YldStres,,5
*DIM,ElasMod,,5
StrnRate(1) = 0,.08,.16,.4,1.0
YldStres(1) = 207e6,250e6,275e6,290e6,300e6
ElasMod(1) = 209e9,211e9,212e9,215e9,218e9
EDCURVE,ADD,1,StrnRate(1),YldStres(1)
EDCURVE,ADD,2,StrnRate(1),ElasMod(1)
! Pa
MP,nuxy,1,.30
! No units
MP,dens,1,7830
! kg/m3
TB,PLAW,,,,8
TBDATA,1,207e6
TBDATA,3,.75
! Failure strain
TBDATA,4,40.0
TBDATA,5,5.0
TBDATA,6,1
! LCID for true stress vs. true strain (see EDCURVE below)
*DIM,TruStran,,5
*DIM,TruStres,,5
TruStran(1)=0,.08,.16,.4,.75
TruStres(1)=207e6,250e6,275e6,290e6,3000e6
EDCURVE,ADD,1,TruStran (1),TruStres(1)
! psi
MP,DENS,1,.03
! lb/in3
MP,NUXY,1,.45
! No units
TB,PLAW,1,,,11
! Modified PLP
TBDATA,1,1000
! Yield strength
TBDATA,2,50000
! Tangent modulus
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
B7
TBDATA,11,1
! Rate (1/s)
TBDATA,12,1
TBDATA,13,4
! Rate (4/s)
TBDATA,14,2
TBDATA,15,45
! Rate 45/s
TBDATA,16,3
! Pa
MP,nuxy,1,.33
! No units
MP,dens,1,2710
! kg/m3
TB,PLAW,,,,2
TBDATA,1,0.598
!k
TBDATA,2,0.216
!n
TBDATA,3,6500.0
! C (s-1)
TBDATA,4,4.0
!P
! kg/m3
MP,alpx,1,10e-6
! m/m-C
MP,ex,1,3.06e9
! Pa
TB,PLAW,1,,,12
TBDATA,1,3.447e7
TBDATA,12,1
TBDATA,13,2
B8
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
TBDATA,15,4
TBDATA,18,5
! Stress-strain curve
*DIM,strn1,,10
*DIM,strs1,,10
strn1(1)=0,.004,.006,.008,.01,.014,.018,.022,.026,.03
strs1(1)=0,1.224e7,1.836e7,2.455e7,3.06e7,3.54e7,3.75e7,3.90e7,4.05e7,4.12e7
EDCURVE,ADD,1,strn1(1),strs1(1)
! Temperature dependent curves
*DIM,temp1,,4
*DIM,nuxy1,,4
*DIM,ex1,,4
*DIM,ystrs1,,4
*DIM,alpha1,,4
temp1(1)=70,120,250,500
ex1(1)=3.06e9,1.85e9,1.12e9,6.79e8
nuxy1(1)=.445,.451,.454,.456
ystrs1(1)=3.447e7,1.85e7,1.042e7,5.627e6
alpha1(1)=1e-5,8.4e-6,7.6e-6,7.2e-6
EDCURVE,ADD,2,temp1(1),ex1(1)
EDCURVE,ADD,3,temp1(1),nuxy1(1)
EDCURVE,ADD,4,temp1(1),ystrs1(1)
EDCURVE,ADD,5,temp1(1),alpha1(1)
! Note: Be sure to define a temperature range that exceeds the actual resulting temperatures. In this particular case, LS-DYNA will expect all temperatures to be between 0
and 500 degrees C. This requirement is on the LS-DYNA side.
! ksi
MP,dens,1,2.226e-7
! kip-sec2/in4
TB,GCAP,1
TBDATA,1,2100
! K (ksi)
TBDATA,2,.7
! (ksi)
TBDATA,3,.1
TBDATA,4,.2
! (ksi)
TBDATA,5,1.473
! (ksi-1)
TBDATA,6,10.8
!R
TBDATA,7,.00154
! D (ksi-1)
TBDATA,8,.884
!W
TBDATA,9,18
! X0 (ksi)
TBDATA,10,0
!C
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
B9
!N
TBDATA,12,1
! Ftype
TBDATA,13,-0.3
! Toff
! Pa
MP,nuxy,1,.30
! No units
MP,dens,1,7850
! kg/m3
TB,EOS,1,,,1,1
TBDATA,1,350.25e6
! A (Pa)
TBDATA,2,275e6
! B (Pa)
TBDATA,3,.36
!n
TBDATA,4,.022
!c
TBDATA,5,1.0
!m
TBDATA,6,1400
TBDATA,7,30
TBDATA,8,10
TBDATA,9,4500
! Specific heat
TBDATA,10,240e6
! Failure stress
TBDATA,11,-.8
! Failure value D1
TBDATA,12,2.1
! Failure value D2
TBDATA,13,-.5
! Failure value D3
TBDATA,14,.0002
! Failure value D4
TBDATA,15,.61
! Failure value D5
TBDATA,17,140e9
! Pa
MP,nuxy,1,.35
! No units
MP,dens,1,8330
! kg/m3
TB,EOS,1,,,1,2
TBDATA,1,89.63e6
! A (Pa)
TBDATA,2,291.64e6
! B (Pa)
TBDATA,3,.31
!n
TBDATA,4,.025
!c
TBDATA,5,1.09
!m
TBDATA,6,1200
TBDATA,7,30
TBDATA,8,10
TBDATA,9,4400
! Specific heat
TBDATA,10,240e6
! Failure stress
TBDATA,11,-.54
! Failure value D1
B10
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
! Failure value D2
TBDATA,13,-3.03
! Failure value D3
TBDATA,14,.014
! Failure value D4
TBDATA,15,1.12
! Failure value D5
TBDATA,16,.394
!C
TBDATA,17,1.489
! S1
TBDATA,18,0.0
! S2
TBDATA,19,0.0
! S3
TBDATA,20,2.02
! 0
TBDATA,21,.47
!A
! Pa
MP,nuxy,1,.3
! No units
MP,dens,1,7500
! kg/m3
TB,EOS,1,,,2,1
TBDATA,1,0.0
! Pressure cut-off
TBDATA,3,1.5
TBDATA,4,.7
TBDATA,17,16e5
! Pa
MP,nuxy,1,.34
! No units
MP,dens,1,2500
! kg/m3
TB,EOS,1,,,2,2
TBDATA,1,-10000
TBDATA,3,2.0
TBDATA,4,.5
TBDATA,16,.5386
!C
TBDATA,17,1.339
! S1
TBDATA,18,0.0
! S2
TBDATA,19,0.0
! S3
TBDATA,20,1.97
! 0
TBDATA,21,.48
!A
! (Pa)
MP,dens,1,.285
! (kg/m3)
TB,EOS,1,,,5,2
TBDATA,1,49.3e3
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
B11
! Hardening coefficient
TBDATA,3,.35
!n
TBDATA,5,.36e6
TBDATA,10,32
! Atomic weight
TBDATA,11,2380
TBDATA,15,2
! Spall type
TBDATA,16,1
TBDATA,17,-50
TBDATA,18,200
TBDATA,29,.457
!C
TBDATA,30,1.49
! S1
TBDATA,31,0.0
! S2
TBDATA,32,0.0
! S3
TBDATA,33,1.93
! 0
TBDATA,34,1.4
!A
! Pa
MP,nuxy,1,.3
! No units
EDMP,cable,1,1
*DIM,EngStran,,4
*DIM,EngStres,,4
EngStran(1) = .02,.04,.06,.08
EngStres(1) = 207e6,210e6,215e6,220e6
EDCURVE,ADD,1,EngStran (1),EngStres(1)
! Pa
MP,nuxy,1,.3
! No units
MP,dens,1,7580
! kg/m3
EDMP,rigid,1,7,7
B12
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
EDADAPT
*PART
EDALE
*CONTROL_ALE
EDASMP
*SET_PART_LIST
EDBOUND
*BOUNDARY_SLIDING_PLANE *BOUNDARY_CYCLIC
EDBVIS
*CONTROL_BULK_VISCOSITY
EDBX
*DEFINE_BOX
EDCADAPT
*CONTROL_ADAPTIVE
EDCGEN
*CONTACT
EDCMORE
*CONTACT
EDCNSTR
*CONSTRAINED_EXTRA_NODES_SET
*CONSTRAINED_NODAL_RIGID_BODY
*CONSTRAINED_SHELL_TO_SOLID
*CONSTRAINED_RIVET
EDCONTACT
*CONTROL_CONTACT
EDCPU
*CONTROL_CPU
EDCRB
*CONSTRAINED_RIGID_BODIES
EDCSC
*CONTROL_SUBCYCLE
EDCTS
*CONTROL_TIMESTEP
EDCURVE
*DEFINE_CURVE
EDDAMP
*DAMPING_PART_MASS *DAMPING_PART_STIFFNESS
EDDC
*DELETE_CONTACT
EDDRELAX
*CONTROL_DYNAMIC_RELAXATION
EDDUMP
*DATABASE_BINARY_D3DUMP
EDENERGY
*CONTROL_ENERGY
EDGCALE
*CONTROL_ALE
EDHGLS
*CONTROL_HOURGLASS
EDHIST
*DATABASE_HISTORY_NODE *DATABASE_HISTORY
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
EDHTIME
*DATABASE
*DATABASE_BINARY_D3THDT
EDINT
*INTEGRATION_BEAM
*INTEGRATION_SHELL
EDIPART
*PART_INERTIA
EDIS
*STRESS_INITIALIZATION
*STRESS_INITIALIZATION_DISCRETE
EDLCS
*DEFINE_COORDINATE_VECTOR
*DEFINE_COORDINATE_SYSTEM
*DEFINE_VECTOR
*LOAD_NODE_SET
*LOAD_THERMAL_VARIABLE
EDLOAD, , PRESS
*LOAD_SEGMENT
*LOAD_SEGMENT_SET
*LOAD_SHELL_SET
*HOURGLASS
EDNB
*BOUNDARY_NON_REFLECTING
EDNROT
*BOUNDARY_SPC_SET
EDOPT
*DATABASE_FORMAT
EDOUT
*DATABASE_OPTION
EDPART
*PART
EDPVEL
*SET_NODE
*INITIAL_VELOCITY
*INITIAL_VELOCITY_GENERATION
*CHANGE_VELOCITY
*CHANGE_VELOCITY_ZERO
C2
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
EDRC
*RIGID_DEFORMABLE_CONTROL
EDRD
*DEFORMABLE_TO_RIGID
*RIGID_DEFORMABLE_D2R
*RIGID_DEFORMABLE_R2D
EDRI
*DEFORMABLE_TO_RIGID_INERTIA
EDRST
*DATABASE_BINARY_D3PLOT
EDSHELL
*CONTROL_SHELL
EDSP
*CHANGE_SMALL_PENETRATION
EDSTART
EDTERM
*TERMINATION_NODE
*TERMINATION_BODY
EDVEL
*SET_NODE
*INITIAL_VELOCITY
*INITIAL_VELOCITY_GENERATION
*CHANGE_VELOCITY
*CHANGE_VELOCITY_ZERO
EDWELD
*CONSTRAINED_SPOTWELD
*CONSTRAINED_GENERALIZED_WELD_SPOT
BFUNIF or TUNIF
*LOAD_THERMAL_CONSTANT
LDREAD
*LOAD_THERMAL_CONSTANT_NODE
REXPORT
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
C3
C4
! kg/m^3
! W/m-C
! J/kg-C
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
!!!etchg,tts
!!!keyopt,1,2,1
et,2,SOLID185
emodif,all,type,2
etdele,1
numcmp,type
keyopt,1,2,1
!
!
!
!
!
mplist,all,all
fini
/solu
antype,static
D2
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
time,1.0
nsubst,2,1000,2
outres,all,last
ldread,temp,last,,,,thermal,rth
tref,0.0
bflist,all,temp
nsel,s,loc,y,1
cm,nbeamtop,node
esln
cm,ebeamtop,elem
sf,all,press,1.1e5
nsel,s,loc,z,0
d,all,uz,0.0
nsel,s,loc,z,0
nsel,r,loc,y,0.5
d,all,uy,0
nsel,s,loc,z,0
nsel,r,loc,x,0.5
d,all,ux,0
nsel,all
esel,all
/title, ANSYS Implicit Static Preload Analysis
save
solve
fini
/post1
set,last
/dscale,,10
plnsol,s,eqv,2
/wait,3
! Note:
!
!
! exaggerate displacement
save
fini
! = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =
/filnam,dynamic
/psf,pres,norm,2
/dscale,,1.0
/prep7
sfdele,all,all
etchg,ite
!
!
!
!
!
!
!
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
D3
500 degrees C.
edmp,hgls,1,5
edpart,create
eddamp,all,,0.10
eddamp,1,,1.0e-6
mplist,all,all
tblist,all,all
fini
/solu
time,0.2
rexport,dyna,,,,,static,rst
! Note:
!
!
!
!
!
bflist,all,temp
eddrelax,ansys
edrst,100
edhtime,200
edhist,nbeamtop
edhist,ebeamtop
edenergy,1,1,1,1
edout,glstat
edopt,add,,both
edwrite,both,,k
nsel,all
esel,all
/title, ANSYS LS-DYNA Dynamic Relaxation and Transient Analyses
save
solve
fini
/post1
set,first
/dscale,,10
! exaggerate displacement
D4
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
D5
D6
Bibliography
[1.] F. Barlat and J. Lian. "Plastic Behavior and Stretchability of Sheet Metals. Part I: A Yield Function for Orthotropic
Sheets Under Plane Stress Conditions". Int. Journal of Plasticity, Vol. 5. pg. 51-66. 1989.
[2.] F. Barlat, D. J. Lege, and J. C. Brem. "A Six-Component Yield Function for Anistropic Materials". Int. Journal of
Plasticity, Vol. 7. pg. 693-712. 1991.
[3.] R. Hill. "A Theory of the Yielding and Plastic Flow of Anisotropic Metals". Proceedings of the Royal Society of London,
Series A., Vol. 193. pg. 281197. 1948.
[4.] F. K. Chang and K. Y. Chang. "A Progressive Damage Model for Laminated Composites Containing Stress Concentration. Journal of Composite Materials, 21. pg. 834-855. 1987.
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Index
A
Action buttons
Drop Test Animate Over Results dialog box, 1717
Drop Test Graph Time-History dialog box, 1718
Drop Test List Time-History dialog box, 1718
Drop Test Set-up dialog box, 1712, 1713, 1714,
1715
Status tab, 1715
Target tab, 1714
Velocity tab, 1713
Adaptive meshing
in explicit dynamic analysis, 310
Animating Results, Postprocessing, 122
ANSYS LS-DYNA overview, 11
Applying loads
explicit dynamic analysis, 34
Array parameters, 42
B
Beam Elements, 29
BFUNIF command, 412
C
Commands used in explicit dynamic analyses, 11
Components, 42
Constraints, 48
explicit dynamic analysis, 36
Constraints and initial conditions, 47
Contact Surface
Definitions, 61
Contact Surfaces , 61
Contact Options , 65
Definition of Contact Types, 65
Coupling and constraint equations, 411
Creating Components, POST26, 121
D
Data curves, 45
Data Smoothing, 126
defining material models, 72
Discrete Elements, 210
Drop Test Graph Time-History Variables dialog box,
1718
Drop Test List Time-History Variables dialog box, 1718
Drop Test Pick Nodes dialog box, 1716
Drop Test Set-up dialog box
Basic tab, 1710, 1710, 1710, 1711, 1711, 1711
Drop Height, 1710
Gravity, 1710
E
EDHTIME command, 37
EDLOAD command, 412
EDNDTSD command, 36
EDRST command, 37
EDSTART Command, 131
elastic material models, 71
Element Output Data, 123
Element Solution, Postprocessing, 126
Elements
Beam and Link Elements, 29
BEAM161, 29
COMBI165 Spring-Damper, 210
Discrete Elements, 210
explicit dynamics, 21
explicit solid and shell elements, 21
General Element Capabilities, 211
LINK160, 210
LINK167, 210
MASS166, 211
PLANE162, 26
SHELL163, 23
SOLID164, 21
SOLID168, 28
Triangular Shell Formulations, 24
equation of state material models, 72
Example problems, explicit dynamics, 14
Explicit dynamic analysis
adaptive meshing, 310
applying loads, 34
commands, 11
constraints, 36
data smoothing, 36
defining contact, 33
defining parts, 37
displaying load symbols, 35
DOF coupling, 36
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Index
example problems, 14
initial velocities, 35
modeling, 31
overview, 11
viewing load curves, 35
explicit dynamic analysis
output control options, 36
DOF coupling
drelax file, 151
dynamic relaxation, 151
explicit dynamic analysis, 36, 36
explicit-implicit element type pairs, 141
implicit-explicit element type pairs, 151
importing stresses, 141
thermal-explicit element type pairs, 155
explicit-to-implicit sequential solution, 141
F
foam material models, 72
H
Hourglassing , 91
I
implicit-to-explicit sequential solution, 151
Initial velocity, 410
explicit dynamic analysis, 35
L
LDREAD command, 412
Link Elements, 29
Load curves
displaying in explicit dynamics, 35
Load symbols
displaying in explicit dynamics, 35
Load types
explicit dynamic analysis, 35
Loading, 41
array parameters, 42
components, 42
constraints, 48
constraints and initial conditions, 47
coupling and constraint equations, 411
data curves, 45
dynamic relaxation, 413
general loading, 44
general loading options, 41
initial velocity, 410
nonreflecting boundary, 411
temperature, 412
Loads
explicit dynamic analysis, 34
Index2
M
mass scaling, 101
material models, 71
3-parameter Barlat, 711
anisotropic elastic, 74
Bamman plasticity model , 729
Barlat anisotropic plasticity, 713
bilinear isotropic, 77
bilinear kinematic (strain rate independent), 710
Blatz-Ko rubber, 75
cable, 737
closed cell foam, 724
composite damage, 718
concrete damage, 718
crushable foam, 726
defining, 72
descriptions, 73
elastic, 71, 73
elastic fluid, 75
elastic viscoplastic thermal, 719
elastic-plastic hydrodynamic, 721
elastoplastic spring, 736
equation of state, 72
foam, 72
general nonlinear spring, 735
Gruneisen equation of state, 728
honeycomb, 726
inelastic tension- or compression-only spring, 736
Johnson-Cook plasticity, 729
linear elastic spring, 735
linear polynomial equation of state, 727
linear viscous damper, 737
low density foam, 725
maxwell viscoelastic spring, 736
modified piecewise linear plasticity, 716
Mooney-Rivlin rubber, 76
nonlinear elastic, 71
nonlinear elastic spring, 736
nonlinear inelastic, 71
nonlinear viscous damper, 737
null, 731
orthotropic elastic, 74
others, 72
piecewise linear plasticity, 715
plastic kinematic, 711
Power Law plasticity, 719
pressure dependent plasticity, 71
rate sensitive Power Law plasticity, 714
rigid, 738
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Index
spring-damper (discrete), 72
Steinberg model, 733
strain rate dependent plasticity, 714
tabulated equation of state, 729
temperature dependent bilinear isotropic, 78
transversely anisotropic FLD, 79
transversely anisotropic hardening, 78
viscoelastic, 77
viscous foam, 725
Zerilli-Armstrong model, 731
Modeling
in explicit dynamic analysis, 31
Monitoring an LS-DYNA solution , 53
N
Nodal Solution, Postprocessing, 126
nonlinear elastic material models, 71
nonlinear inelastic material models, 71
Nonreflecting boundary, 411
O
other material models, 72
output control options
explicit dynamic analysis, 36
Output Controls, Postprocessing, 121
Output Data, 123
Output Files, POST26, 121
R
restarting, 131
EDSTART Command, 131
effect on Output Files, 135
Restart Dump File, 131
Rigid bodies, 81
inertia properties, 81
loading, 82
switching (D2R and R2D), 82
S
scaling, 101
Shell elements
explicit dynamics, 21
formulations, 23
Shell Formulations, Triangular, 24
SMP, 52
Solid elements
explicit dynamics, 21
Solution
editing the LS-DYNA input file, 54
solution control and monitoring , 53
Solution process, 51
Spring-Damper, 210
spring-damper (discrete) material models, 72
Subcycling, 111
P
Parallel processing in ANSYS LS-DYNA, 52
Parameters,array, 42
Parts
in explicit dynamic analysis, 37
POST26, Creating Components, 121
POST26, Output Files, 121
POST26, use, 125
Postprocessing, 121
Animating Results, 122
Creating Components for POST26, 121
Data Smoothing, 126
Element Output Data, 123
Finding Additional Information, 127
in explicit dynamic analysis, 122
Nodal and Element Solutions, 126
Output Controls, 121
Reading ASCII Files for Miscellaneous Output Data,
126
Results vs. History, 121
Using POST26 with ANSYS LS-DYNA, 125
Writing the Output Files for POST26, 121
pressure dependent plasticity material models, 71
ANSYS LS-DYNA User's Guide . ANSYS Release 10.0 . 002184 . SAS IP, Inc.
Index3