Drop Test of a Container (Explicit Dynamics)
6.1.3.1.Step 1: Set Preferences.
You will first set preferences in order to filter quantities that pertain to this discipline only.
Main Menu> Preferences
(check) “Individual discipline(s) to show in the GUI” = Structural
(check) “Discipline options” = LS-DYNA Explicit
[OK]
6.1.4.Input Geometry
6.1.4.1.Step 2: Read in geometry of the container.
You will begin by reading in a file that includes the model of the container.
Utility Menu> File> Read Input from ...
File name: container.inp
UNIX version:
/ansys_inc/v130/ansys/data/models/container.inp
PC version:
\Program Files\Ansys Inc\V130\ANSYS\data\models\container.inp
[OK]
6.1.5.Define Element Type, Real Constants, Material Model Properties
6.1.5.1.Step 3: Define element type.
Since the model is a container made of thin walls hitting the surface of a table, shell elements will be used to simulate the physical model. More precisely,
SHELL163, an explicit thin structural shell, is the element of choice. This 4-node element can take both in-plane and normal loads.
Main Menu> Preprocessor> Element Type> Add/Edit/Delete
[Add...]
“LS-DYNA Explicit” (left column)
“Thin Shell 163” (right column)
[OK]
Now, specify the S/R corotational Hughes-Liu element formulation to eliminate certain hourglass modes.
[Options...]
(drop down) “Element Formulation” = S/R corotation
[OK]
[Close]
6.1.5.2.Step 4: Define real constants.
Next, define the thickness of the shell elements by defining the appropriate real constant.
Main Menu> Preprocessor> Real Constants
[Add...]
[OK] to choose element Type 1 for
SHELL163.
[OK] to define real constant set No. 1.
“SHRF” = 1
“NIP” = 2 (For demonstration purposes in this tutorial, the default value of 2 is sufficient. However, for most nonlinear analyses, the NIP should be set higher than 2.)
“T1” = 0.1
“T2” = 0.1
“T3” = 0.1
“T4” = 0.1
“NLOC” = Mid surface
[OK]
[Close]
6.1.5.3.Step 5: Specify material models.
Now specify material models for the contact and target surfaces.
Main Menu> Preprocessor> Material Props> Material Models
(double-click) “LS-DYNA”, then “Rigid Material” to specify the table top, which is “Material Model Number 1”.
“DENS” = 7.3e-4
“EX” = 30e6
“NUXY” = 0.292
(drop down) “Translational Constraint Parameter” = All disps.
(drop down) “Rotational Constraint Parameter” = All rotations
[OK]
Material> New Model to specify the material for the container, which is “Material Model Number 2”.
“Define Material ID” = 2
[OK]
(double-click) “Nonlinear”, then “Inelastic”, then “Kinematic Hardening”, then “Bilinear Kinematic”
“DENS” = 2.5e-4
“EX” = 10.3e6
“NUXY” = 0.334
“Yield Stress” = 5000
“Tangent Modulus” = 20000
[OK]
Material> Exit
Toolbar: SAVE_DB
6.1.6.Generate Mesh
6.1.6.1.Step 6: Mesh the container.
In an explicit solver, such as ANSYS LS-DYNA, the analysis time is highly dependent on the smallest element in the model. Therefore, for efficient results, a uniform element size should be used.
For this example, use the default element size. This should give a 3x3 mesh on each face of the container.
Main Menu> Preprocessor> Meshing> MeshTool
(drop down) “Element Attributes” = Global; then [Set]
(drop down) “Material number” = 2
[OK]
(drop down) “Mesh” = Areas
(check) “Mapped”
[Mesh]
[Pick All]
Toolbar: SAVE_DB
6.1.6.2.Step 7: Generate table top elements.
Now mesh the table top by direct element generation. The table top will be defined as a rigid surface and therefore only one element is needed to model this surface. Create this model by defining its four corner nodes and then building the elements from those nodes.
Utility Menu> PlotCtrls> Pan, Zoom, Rotate
[Obliq]
[Close]
Utility Menu> PlotCtrls> Numbering
(check) “Node numbers” = On
[OK]
(drop down in MeshTool) “Element Attributes” = Global; then [Set]
(drop down) “Material number” = 1
[OK]
[Close] MeshTool
Main Menu> Preprocessor> Modeling> Create> Nodes> In Active CS
“Node number” = Leave blank so it defaults to the next available node number.
“X Y Z Location in active CS” = -50, -20, -50
[Apply] to generate the node in the far left corner.
“X Y Z Location in active CS” = -50, -20, 50
[Apply] to generate the node in the near left corner.
“X Y Z Location in active CS” = 50, -20, 50
[Apply] to generate the node in the near right corner.
“X Y Z Location in active CS” = 50, -20, -50
[OK] to generate the node in the far right corner.
Now create the element.
Main Menu> Preprocessor> Modeling> Create> Elements> Auto Numbered> Thru Nodes
In a counterclockwise direction, pick the four nodes just created.
[OK]
Utility Menu> Plot> Elements
6.1.6.3.Step 8: Create container component.
Most of the contact algorithms require contact parameters, which can be components, part IDs, or part assembly IDs. For this tutorial, create a component made of nodes from the container.
Utility Menu> Select> Entities
(first drop down) “Elements”
(second drop down) “By Attributes”
(check) “Material num”
“Min,Max,Inc.” = 2
[Apply]
(first drop down) “Nodes”
(second drop down) “Attached to”
(check) “Elements”
[OK]
Utility Menu> Select> Comp/Assembly> Create Component
“Component name” = BOX
(drop down) “Component is made of” = Nodes
[OK]
Utility Menu> Plot> Nodes
Utility Menu> Select> Everything
6.1.6.4.Step 9: Create table top component.
Now, create a component made of nodes from the table top.
Utility Menu> Select> Entities
(first drop down) “Elements”
(second drop down) “By Attributes”
(check) “Material num”
“Min,Max,Inc” = 1
[Apply]
(first drop down) “Nodes”
(second drop down) “Attached to”
(check) “Elements”
[OK]
Utility Menu> Select> Comp/Assembly> Create Component
“Component name” = TABLE
(drop down) “Component is made of” = Nodes
[OK]
Utility Menu> Plot> Nodes
Utility Menu> Select> Everything
Toolbar: SAVE_DB
6.1.6.5.Step 10: Specify contact parameters.
In an analysis in which the contact behavior is known, the appropriate contact algorithm can be chosen for the best results. However, due to the nature of this problem, contact conditions will be unpredictable. Thus, by choosing automatic surface to surface contact, the program will automatically adjust for the changes which occur during this simulation.
Main Menu> Preprocessor> LS-DYNA Options> Contact> Define Contact
“Contact Type” = “Surface to Surf” (left column); “Automatic (ASTS)” (right column)
[OK]
(drop down) “Contact Component or Part no.” = BOX
(drop down) “Target Component or Part no.” = TABLE
[OK]
Toolbar: SAVE_DB
6.1.7.Apply Loads
6.1.7.1.Step 11: Apply initial velocity to the container.
Main Menu> Solution> Initial Velocity> On Nodes> w/Nodal Rotate
(drop down) “Input velocity on component” = BOX
“Global (Translational) Y-component” = -200
[OK]
In a transient dynamic analysis, loads must be defined for the duration of the analysis. To do so, create an array containing this data.
Utility Menu> Parameters> Array Parameters> Define/Edit
[Add...]
“Parameter name” = TIME
[OK]
[Edit...]
“1” = 0
“2” = 1
File> Apply/Quit
[Add...]
“Parameter name” = ACCG
[OK]
(highlight) “ACCG”
[Edit...]
“1” = 386.4
“2” = 386.4
File> Apply/Quit
[Close]
6.1.7.2.Step 12: Apply acceleration to the container.
Main Menu> Solution> Loading Options> Specify Loads
“Load Labels” = ACLY
(drop down) “Component name or PART number:” = BOX
(drop down) “Parameter name for time values:” = TIME
(drop down) “Parameter name for data values:” = ACCG
[OK]
6.1.8.Obtain Solution
6.1.8.1.Step 13: Specify output controls.
Main Menu> Solution> Time Controls> Solution Time
“Terminate at Time:” = 1.0
[OK]
Main Menu> Solution> Output Controls> File Output Freq> Number of Steps
“Specify Results File Output Interval” = 50
“Specify Time History Output Interval” = 50
[OK]
Main Menu> Solution> Analysis Options> Energy Options
[OK] to turn on all energy options
Toolbar: SAVE_DB
6.1.8.2.Step 14: Solve.
Main Menu> Solution> Solve
Review the information in the status window, then choose:
File> Close (Windows) or
Close (X11/Motif), to close the window.
[OK] to initiate solution.
Note: A verification window will appear to notify you of a warning message. This warning was generated due to constraining the table top by defining it as a rigid element rather than using degrees of freedom constraints. Therefore the warning can be ignored.
[Yes]
[Close]
6.1.9.Review Results
6.1.9.1.Step 15: Animate stress contours.
Utility Menu> Plot> Elements
Utility Menu> PlotCtrls> Pan, Zoom, Rotate
[Front]
[Close]
Utility Menu> PlotCtrls> Numbering
(check) “Node numbers” = Off
[OK]
Main Menu> General Postproc> Read Results> First set
[Note for UNIX systems only: To capture the animation sequence in terminal segment memory, it is necessary to reduce the size of the Graphics Window or you will run out of terminal memory. Although this depends on your particular system, and the size of the animation file to be created, for this exercise it is suggested that the ratio of the original window size to the reduced window size should be about 3:1.]
Utility Menu> PlotCtrls> Animate> Over Results
(check) “Auto contour scaling” = On (Make sure that the box is checked.)
“Contour data for animation” = Stress (left column); von Mises SEQV (right column)
[OK]
Make choices in the Animation Controller (not shown), if necessary, then [Close].
6.1.9.2.Step 16: Animate deformed shape.
Utility Menu> PlotCtrls> Animate> Over Results
(check) “Auto contour scaling” = On
“Contour data for animation” = DOF Solution (left column); Deformed Shape (right column)
[OK]
Make choices in the Animation Controller (not shown), if necessary, then [Close].
6.1.9.3.Step 17: Exit the ANSYS program.
Toolbar: Quit
(check) “Quit - No Save!”
[OK]