Reduced Input
*Node and *Element Cards are removed.
*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