A rigid sphere drops onto an elastic plate. The sphere contains shell elements automatically generated with a “Geometric Contact Entity” spherical surface.
*CONTACT_ENTITY Rigid Sphere Impacts Plate LS-DYNA Manual Section: *CONTACT_ENTITY Additional Sections: *BOUNDARY_PRESCRIBED_MOTION_RIGID Example: Rigid Sphere Impacts Plate Filename: contact_entity.sphere.k Description: A rigid sphere drops onto an elastic plate. The sphere contains shell elements automatically generated with a "Geometric Contact Entity" spherical surface. Model: The plate of elastic material measures 40 * 40 * 2 mm3 and contains 64 Belytschko-Tsay shell elements. The sphere has a radius of 6.0 mm and the distance from the center of the cube to the plate is 8.5 mm. The inertia properties of the sphere are defined by the properties of the rigid brick element. A geometric contact entity defines the spherical contact surface. The sphere moves toward the plate with a uniform motion. The termination time is 0.0005 seconds. Input: The Geometric Contact Entity defines the outer master surface on the rigid sphere (*CONTACT_ENTITY). The nodes on the plate are slave nodes (*SET_NODE_LIST), and are in the "Geometric Entity". A load curve definition defines the movement of the sphere (* BOUNDARY_PRESCRIBED_MOTION_RIGID, *DEFINE_CURVE). The displacement condition for rigid bodies is input by part number, not by listing the nodes included in the definition. Reference: Schweizerhof, K. and Weimar, K.
$ * Part 1 - plate $ * Part 2 - sphere - Contact Entity *BOUNDARY_PRESCRIBED_MOTION_RIGID *CONTACT_ENTITY *CONTROL_HOURGLASS *CONTROL_TERMINATION *CONTROL_TIMESTEP *DATABASE_BINARY_D3PLOT *DATABASE_GCEOUT *DATABASE_RBDOUT *DEFINE_CURVE *ELEMENT_SHELL *END *KEYWORD *MAT_ELASTIC *MAT_RIGID *NODE *PART *SECTION_SHELL *SECTION_SOLID *SET_NODE_LIST *TITLE
*KEYWORD
*TITLE
Geometric Contact Entity
$
$ LSTC Example
$
$ * Part 1 - plate
$ Shells (2.0 mm thick)
$ elastic material
$ translational constraints on all four edges (z-dir only)
$
$ * Part 2 - sphere - Contact Entity
$ Defined as shells and rigid material but really there are no part 2
$ elements defined explicitly. The contact entity is really part 2.
$ center (x,y,z) = (20,20,9)
$ radius = 6 mm
$
$
$ ===> Due to the course mesh of the plate, there is considerable amount
$ of penetration of the sphere into the plate.
$
$ Last Modified: April 10, 1997
$
$ Units: mm, ton, s, N, MPa, N-mm
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ Control/Ouput
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
$
*CONTROL_TERMINATION
$ endtim endcyc dtmin endneg endmas
.5000E-3
$
*CONTROL_TIMESTEP
$ dtinit scft isdo tslimt dtms lctm erode ms1st
0.1
$
*CONTROL_HOURGLASS
$ ihq qh
4
$
$
*DATABASE_BINARY_D3PLOT
$ dt/cycl lcdt
0.0200E-3
$
*DATABASE_RBDOUT
$ dt/cycl lcdt
0.005e-3
$
*DATABASE_GCEOUT
$ dt/cycl lcdt
0.005e-3
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ Contact Entity
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
$
*CONTACT_ENTITY
$ pid geotyp sid styp sf df cf intord
2 2 1 0 1.0
$
$ bt dt so
$
$ xc yc zc ax ay az
0.00 0.00 .00 1.00 0.00 0.00
$
$ bx by bz
0.00 1.00 0.00
$
$ inout g1 g2 g3 g4 g5 g6 g7
0 20.00 20.00 9.00 6.00
$
$
*SET_NODE_LIST
$ sid da1 da2 da3 da4
1
$
$ nid1 nid2 nid3 nid4 nid5 nid6 nid7 nid8
1 2 3 4 5 6 7 8
9 10 11 12 13 14 15 16
17 18 19 20 21 22 23 24
25 26 27 28 29 30 31 32
33 34 35 36 37 38 39 40
41 42 43 44 45 46 47 48
49 50 51 52 53 54 55 56
57 58 59 60 61 62 63 64
65 66 67 68 69 70 71 72
73 74 75 76 77 78 79 80
81
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ Define Parts and Materials
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
$
$$$ Part 1 shell: plate - elastic material
$
$$$ Part 2 solid: sphere - rigid material ==> contact entity
$
*PART
plate
$ pid sid mid eosid hgid adpopt
1 1 1 0 0 0
$
*PART
sphere
$ pid sid mid eosid hgid adpopt
2 2 2 0 0 0
$
$
$$$$ Materials
$
*MAT_ELASTIC
$ mid ro e pr da db k
1 2.00e-08 100000. 0.300
$
*MAT_RIGID
$ mid ro e pr n couple m alias
2 2.00e-08 100000. 0.300 4
$
$ cmo con1 con2
$
$ lco/a1 a2 a3 v1 v2 v3
$
$
$$$$ Sections
$
*SECTION_SHELL
$ sid elform shrf nip propt qr/irid icomp
1 0.83333 2.0 3.0
$
$ t1 t2 t3 t4 nloc
2.0 2.0 2.0 2.0
$
*SECTION_SOLID
$ sid elform
2
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ Boundary and Initial Conditions
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
$
*BOUNDARY_PRESCRIBED_MOTION_RIGID
$ mid dof vad lcid sf vid
2 3 2 1 1.000E+00
$
*DEFINE_CURVE
$ lcid sidr scla sclo offa offo
1
$
$ abscissa ordinate
0.0 0.0
0.0050 -150.0
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ Define Nodes and Elements
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
$
$
*END
