*DEFORMABLE_TO_RIGID
Interaction of Pendulums
LS-DYNA Manual Section: *DEFORMABLE_TO_RIGID
Additional Sections:
*BOUNDARY_SPC_NODE
*LOAD_BODY_Y *RIGID_DEFORMABLE_R2D
Example: Interaction of Pendulums
Filenames: deformable_to_ri gid.pendulum.k
deformable_to_rigid.pendulum.res
Execution lines:
ls940 i= deformable_to_rigid.pendulum.k
ls940 i= deformable_to_rigid.pendulum.res r=d3dump01
Description:
Two spheres are connected to wires to form two pendulums. One sphere is in a
horizontal position with gravitational acceleration, base acceleration and is
given an initial velocity in the vertical direction. The other sphere is in
the vertical direction. The spheres are treated as rigid bodies while no
contact or deformation occurs (i.e., when the horizontal pendulum swings down
towards the vertical pendulum). The spheres are switched to deformable through
a restart file so that they become flexible during contact.
Model:
Both spheres are modeled using shell elements. The pendulum wires are modeled
using elastic beams. Automatic single surface contact is used during the impact
phase.
Reference:
Reid, J.D.
deformable_to_rigid.pendulum.k:$ - uses *DEFORMABLE_TO_RIGID option to decrease execution time before impact
deformable_to_rigid.pendulum.k:*BOUNDARY_SPC_NODE
deformable_to_rigid.pendulum.k:*CONTACT_AUTOMATIC_SINGLE_SURFACE
deformable_to_rigid.pendulum.k:*CONTROL_CONTACT
deformable_to_rigid.pendulum.k:*CONTROL_ENERGY
deformable_to_rigid.pendulum.k:*CONTROL_OUTPUT
deformable_to_rigid.pendulum.k:*CONTROL_SHELL
deformable_to_rigid.pendulum.k:*CONTROL_TERMINATION
deformable_to_rigid.pendulum.k:*DATABASE_BINARY_D3PLOT
deformable_to_rigid.pendulum.k:*DATABASE_BINARY_D3THDT
deformable_to_rigid.pendulum.k:*DATABASE_EXTENT_BINARY
deformable_to_rigid.pendulum.k:*DATABASE_GLSTAT
deformable_to_rigid.pendulum.k:*DATABASE_HISTORY_NODE
deformable_to_rigid.pendulum.k:*DATABASE_MATSUM
deformable_to_rigid.pendulum.k:*DATABASE_NODOUT
deformable_to_rigid.pendulum.k:*DATABASE_RBDOUT
deformable_to_rigid.pendulum.k:*DEFINE_BOX
deformable_to_rigid.pendulum.k:*DEFINE_CURVE
deformable_to_rigid.pendulum.k:*DEFORMABLE_TO_RIGID
deformable_to_rigid.pendulum.k:*ELEMENT_BEAM
deformable_to_rigid.pendulum.k:*ELEMENT_SHELL
deformable_to_rigid.pendulum.k:*END
deformable_to_rigid.pendulum.k:*INITIAL_VELOCITY
deformable_to_rigid.pendulum.k:*KEYWORD
deformable_to_rigid.pendulum.k:*LOAD_BODY_Y
deformable_to_rigid.pendulum.k:*MAT_ELASTIC
deformable_to_rigid.pendulum.k:*MAT_PLASTIC_KINEMATIC
deformable_to_rigid.pendulum.k:*NODE
deformable_to_rigid.pendulum.k:*PART
deformable_to_rigid.pendulum.k:*SECTION_BEAM
deformable_to_rigid.pendulum.k:*SECTION_SHELL
deformable_to_rigid.pendulum.k:*TITLE
deformable_to_rigid.pendulum.res:*CONTROL_TERMINATION
deformable_to_rigid.pendulum.res:*DATABASE_BINARY_D3PLOT
deformable_to_rigid.pendulum.res:*END
deformable_to_rigid.pendulum.res:*KEYWORD
deformable_to_rigid.pendulum.res:*RIGID_DEFORMABLE_R2D
deformable_to_rigid.pendulum.res:*TITLE
*KEYWORD
*TITLE
Pendulum with 2 spheres colliding
$
$ LSTC Example
$
$ - uses *DEFORMABLE_TO_RIGID option to decrease execution time before impact
$
$ - one sphere is given an initial velocity (gravity alone just takes
$ too long for the pendulum to swing)
$
$
$ Last Modified: September 16, 1997
$
$ Units: mm, kg, ms, kN, GPa, kN-mm
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ Control Ouput
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
$
*CONTROL_TERMINATION
$ endtim endcyc dtmin endeng endmas
11.0 0 0.0 0.0 0.0
$
*CONTROL_CONTACT
$ slsfac rwpnal islchk shlchk penopt thkchg orien
2
$ usrstr usrfrc nsbcs interm xpene
$
*CONTROL_ENERGY
$ hgen rwen slnten rylen
2 2
$
*CONTROL_OUTPUT
$ npopt neecho nrefup iaccop opifs ipnint ikedit
1 3
$
*CONTROL_SHELL
$ wrpang itrist irnxx istupd theory bwc miter
1 2
$
$
*DATABASE_BINARY_D3PLOT
$ dt lcdt
1.00
$
*DATABASE_EXTENT_BINARY
$ neiph neips maxint strflg sigflg epsflg rltflg engflg
$ cmpflg ieverp beamip
1
$
*DATABASE_BINARY_D3THDT
$ dt lcdt
999999
$
$
*DATABASE_GLSTAT
$ dt
0.10
$
*DATABASE_MATSUM
$ dt
0.10
$
*DATABASE_NODOUT
$ dt
0.10
$
*DATABASE_HISTORY_NODE
$ id1 id2 id3
350 374 678 713
$
*DATABASE_RBDOUT
$ dt
0.10
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ Define Contacts - Sliding Interfaces
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
$
*CONTACT_AUTOMATIC_SINGLE_SURFACE
$ ssid msid sstyp mstyp
0
$ Equating ssid to zero means that all segments are included
$
$ fs fd
0.08 0.08
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ Gravity
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*LOAD_BODY_Y
$ lcid sf lciddr xc yc zc
1 0.00981
$
*DEFINE_CURVE
$ lcid sidr scla sclo offa offo
1
$ abscissa ordinate
$
0.00 1.000
10000.00 1.000
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ Boundary and Initial Conditions
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
$
$$$$ Constrain translation of end points of beams
$
*BOUNDARY_SPC_NODE
$ nid cid dofx dofy dofz dofrx dofry dofrz
45004 0 1 1 1 0 0 0
45005 0 1 1 1 0 0 0
45010 0 1 1 1 0 0 0
45011 0 1 1 1 0 0 0
$
$
$$$$ The nodes within box 5 are given an initial velocity.
$
*INITIAL_VELOCITY
$ nsid nsidex boxid
5
$
$ vx vy vz wx wy wz
0.0 -12.0 0.0
$
*DEFINE_BOX
$ boxid xmm xmx ymn ymx zmn zmx
5 -120.0 -80.0 80.0 120.0 -30.0 30.0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ Define Parts and Materials
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
$
$$$$$ SPHERES
$
*PART
$ pid sid mid eosid hgid adpopt
sphere1
1 1 1
sphere2
2 2 1
$
$
$$$$ Materials
$
$ Aluminum
$
*MAT_PLASTIC_KINEMATIC
$ mid ro e pr sigy etan beta
1 2.70e-6 68.9 0.330 0.286 0.00689
$ src srp fs
$
$$$$ Sections
$
$
*SECTION_SHELL
$
$ sid elform shrf nip propt qr/irid icomp
1 2
$ t1 t2 t3 t4 nloc
1.0 1.0 1.0 1.0
$
*SECTION_SHELL
$
$ sid elform shrf nip propt qr/irid icomp
2 2
$ t1 t2 t3 t4 nloc
1.0 1.0 1.0 1.0
$
$
$$$$$ PENDULUM WIRES - ELASTIC BEAMS
$
*PART
$ pid sid mid eosid hgid adpopt
Pendulum Wires - Elastic Beams
45 45 45
$
$
*MAT_ELASTIC
$ mid ro e pr da db k
45 7.86e-6 210.0 0.30
$
$
*SECTION_BEAM
$ sid elform shrf qr/irid cst
45 3 1.00000 1.0
$ res: a iss itt irr sa
10.0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ Deformable Switching
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
$
*DEFORMABLE_TO_RIGID
$ pid mrb
1
$
*DEFORMABLE_TO_RIGID
$ pid mrb
2
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ Define Nodes
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
$
$
*END