Shell 2 — Welcome to LS-DYNA Examples

Home » Knowledge Base » Simple Examples » Examples Manual » Control » Shell 2

A beam twisted 90 degrees about its length is constrained on on e edge and has a point load prescribed normal to the opposite end of the beam.

Introduction

*CONTROL_SHELL
Twisted Cantilever Beam
LS-DYNA Manual Section: *CONTROL_SHELL

Example: Twisted Cantilever Beam

Filename: control_shell.beam-twist.k

Description:
A beam twisted 90 degrees about its length is constrained on on e edge and has
a point load prescribed normal to the opposite end of the beam.

Model:
The beam measures 12.00 * 1.10 * 0.32 cubic inches. A concentrated load is
applied to one node on the end in the x-direction and the other node on the
end in the z-direction.

Input:
This model uses the Hughes-Liu five through the thickness integration points
(*CONTROL_SHELL, *SECTI ON_SHELL). The element has the shell normal update
calculation performed at each nodal fiber every cycle (*CONTROL_SHELL).
Note: This is another example that will not work correctly with the B-T shell
formulation (unless warping stiffness is added).

Results:
The beam oscillates about a neutral amplitude.

Reference:
Belytschko, Wang and Chiang.

Keywords

*CONTROL_OUTPUT
*CONTROL_SHELL
*CONTROL_TERMINATION
*DATABASE_BINARY_D3PLOT
*DATABASE_BINARY_D3THDT
*DATABASE_BNDOUT
*DATABASE_EXTENT_BINARY
*DATABASE_GLSTAT
*DATABASE_HISTORY_NODE
*DATABASE_NODOUT
*DEFINE_CURVE
*ELEMENT_SHELL
*END
*KEYWORD
*LOAD_NODE_POINT
*MAT_ELASTIC
*NODE
*PART
*SECTION_SHELL
*TITLE

Reduced Input

  *KEYWORD
  *TITLE
  Twisted Beam        
  $
  $  LSTC Example
  $
  $  Last Modified: September 15, 1997
  $
  $  Units: lbf-s2/in, in, s, lbf, psi, lbf-in
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Control Ouput
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  *CONTROL_TERMINATION
  $   endtim    endcyc     dtmin    endneg    endmas
       0.018
  $
  *CONTROL_OUTPUT
  $    npopt    neecho    nrefup    iaccop     opifs    ipnint    ikedit
           0         0         0         0                   2      1000
  $
  *CONTROL_SHELL
  $   wrpang    itrist     irnxx    istupd    theory       bwc     miter
                              -2                   1
  $
  $
  *DATABASE_EXTENT_BINARY
  $    neiph     neips    maxint    strflg    sigflg    epsflg    rltflg    engflg
                               4          
  $   cmpflg    ieverp    beamip
  
  $
  *DATABASE_BINARY_D3PLOT
  $       dt      lcdt
       0.001
  $
  *DATABASE_BINARY_D3THDT
  $       dt      lcdt
      0.0001
  $
  *DATABASE_BNDOUT
  $       dt
      0.0001
  $
  *DATABASE_GLSTAT
  $       dt
      0.0001
  $
  *DATABASE_NODOUT
  $       dt
      0.0001
  $
  *DATABASE_HISTORY_NODE
  $      id1       id2       id3       id4       id5       id6       id7       id8
          37
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Loading
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Load node 37, 38, 39 with 0.1667 lbs in both x and y direction.
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  *LOAD_NODE_POINT
  $      nid       dof      lcid        sf       cid        m1        m2        m3
          37         1         1 1.667E-01
          38         1         1 1.667E-01
          39         1         1 1.667E-01
          37         2         1 1.667E-01
          38         2         1 1.667E-01
          39         2         1 1.667E-01
  $
  *DEFINE_CURVE
  $     lcid      sidr      scla      sclo      offa      offo
           1
  $                  a                   o
             0.000E+00           1.000E+00
             1.000E+03           1.000E+00
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Define Parts and Materials
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  *PART
  $      pid       sid       mid     eosid      hgid    adpopt
  Twisted Beam    
           1         1         1
  $
  $
  *MAT_ELASTIC
  $      mid        ro         e        pr        da        db
           1   2.00e-4  29.00e+6     0.330
  $
  $
  *SECTION_SHELL
  $      sid    elform      shrf       nip     propt   qr/irid     icomp
           1         0                   4                   0
  $       t1        t2        t3        t4      nloc
       0.320     0.320     0.320     0.320
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Define Nodes and Elements
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Nodes 1, 2, 3 have fixed boundary conditions.
  $
  $
  *END

Figures

image01.gif

image02.gif

Animated Result

Result

Download