DEFINE_ELEMENT_SET

Output
Attention: This command is in the beta stage and the format may change over time.
*DEFINE_ELEMENT_SET
"Optional title"
coid
entype, enid, fid, padding
Parameter definition
VariableDescription
coid Command ID
entype Entity type
options: P, PS
enid Entity identification number
fid Function ID
padding Padding
default: no padding
Description

This command is used to define a set of elements (SET_ELEMENT) at termination time. A typical application is to identify and mark critical regions of a model, that need to be refined in a subsequent simulation.

The set is written to the file element_set_X.k, where X=coid. The FUNCTION is evaluated for each element in the part or part set. An element will be included in the part set if the function returns a positive value.

The parameter padding is a positive integer. It is used to include one or several layers of elements surrounding elements satisfying the inclusion criterion.

Example
Change element type in impact zone

In this example a cylindrical impactor punches a hole in a plate. The model is first run with only linear elements (Step 1). The purpose of this step is to identify the region undergoing large deformations (here defined as elements with effective plastic strains larger than 0.5). The region (with or without padding) is defined as an element set. In subsequent simulations (Step 2) elements in the set are converted to 3rd order hexahedra (for an increased accuracy).


Step 1:

*UNIT_SYSTEM
SI
#
# --- TIME ---
#
*TIME
2.0e-5
*OUTPUT
1.0e-5
#
# --- MESH ---
#
*COMPONENT_BOX
"Plate"
1, 1, 150, 150, 3
-0.075, -0.075, -0.003, 0.075, 0.075, 0.0
*COMPONENT_CYLINDER
"Punch"
2, 2, 10, 4
0, 0, 0, 0, 0, 0.03, 0.02
*CHANGE_P-ORDER
P, 2, 3
*SMOOTH_MESH
P, 2, 45.0
#
# --- MATERIAL ---
#
*MAT_METAL
"Plate"
1, 7800.0, 210.0e9, 0.3, 1
1
*FUNCTION
1
4.0e8 + 4.0e8*epsp^0.3
*PROP_DAMAGE_CL
1, 4
8.0e8
#
*MAT_METAL
"Punch"
2, 7800.0, 210.0e9, 0.3
2
*FUNCTION
2
1.0e9 + 1.0e9*epsp^0.3
#
# --- PART ---
#
*PART
"Plate"
1, 1, 0, 0, 0, 0, 3.0
"Punch"
2, 2
#
# --- INITIAL VELOCITY ---
#
*INITIAL_VELOCITY
P, 2, 0, 0, -600.0
#
# --- CONTACT ---
#
*CONTACT
"All to All"
1
ALL, 0, ALL, 0
#
# --- ELEMENT SET ---
#
*DEFINE_ELEMENT_SET
"Without padding"
1
P, 1, 123
#
*DEFINE_ELEMENT_SET
"With padding"
2
P, 1, 123, 2
#
*FUNCTION
123
epsp - 0.5
*END

Step 2 (without padding):

*INCLUDE
element_set_1.k
*COMPONENT_BOX
"Plate"
1, 1, 150, 150, 3
-0.075, -0.075, -0.003, 0.075, 0.075, 0.0
*CHANGE-P-ORDER
ES, 1, 3

Step 2 (with padding):

*INCLUDE
element_set_2.k
*COMPONENT_BOX
"Plate"
1, 1, 150, 150, 3
-0.075, -0.075, -0.003, 0.075, 0.075, 0.0
*CHANGE-P-ORDER
ES, 2, 3