Shell2Solid

A dome has an impulsive pressure load. The dome contains shell and brick element joined with shell-brick interfaces.

*CONSTRAINED_SHELL_TO_SOLID
Impulsively Loaded Cap with Shells and Solids
LS-DYNA Manual Section: *CONSTRAINED_SHELL_TO_SOLID

Additional Sections:
*LOAD_SEGMENT

Example: Impulsively Loaded Cap with Shells and Solids

Filename: constrained.shell_solid.dome.k

Description:
A dome has an impulsive pressure load. The dome contains shell and brick
element joined with shell-brick interfaces.

Model:
Only 1/4 of the dome is modeled due to symmetry. The dome shells are
Hughes-Liu shell elements with three integration point through the thickness.
Four shell elements have a pressure load of 5,308 psi over 0.0017246
square inches. The termination time is 0.0004 seconds.

Input:
The model contains one shell-brick group that has 7 shell nodes tied to 5 brick
node (*CONSTRAINED_SHELL_TO_SOLID). The model contains four pressure surfaces
(*LOAD_SEGMENT). Five nodes are written to the time history ASCII database file
nodout (*DATABASE_HISTORY_NODE, *DATABASE_NODOUT).

Results:
The plots show the response of the dome.

Reference:
T. Littlewood

*CONSTRAINED_SHELL_TO_SOLID
*CONTROL_OUTPUT
*CONTROL_TERMINATION
*DATABASE_BINARY_D3PLOT
*DATABASE_BINARY_D3THDT
*DATABASE_ELOUT
*DATABASE_EXTENT_BINARY
*DATABASE_GLSTAT
*DATABASE_HISTORY_NODE
*DATABASE_HISTORY_SHELL
*DATABASE_HISTORY_SOLID
*DATABASE_NODOUT
*DEFINE_CURVE
*ELEMENT_SHELL
*ELEMENT_SOLID
*END
*KEYWORD
*LOAD_SEGMENT
*MAT_ELASTIC
*NODE
*PART
*SECTION_SHELL
*SECTION_SOLID
*SET_NODE_LIST
*TITLE

  *KEYWORD
  *TITLE
  Impulsively Loaded Cap with Shell-Brick Interfaces
  $
  $  LSTC Example
  $
  $  Last Modified: September 4, 1997
  $
  $   Units: lbf-s^2/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.0004
  $
  *CONTROL_OUTPUT
  $    npopt    neecho    nrefup    iaccop     opifs    ipnint    ikedit
           1         3
  $
  *DATABASE_BINARY_D3PLOT
  $       dt      lcdt
     0.00001
  $
  *DATABASE_BINARY_D3THDT
  $       dt      lcdt
   5.000E-07
  $
  *DATABASE_EXTENT_BINARY
  $    neiph     neips    maxint    strflg    sigflg    epsflg    rltflg    engflg
                                         1
  $
  $   cmpflg    ieverp    beamip
  
  $
  *DATABASE_ELOUT
  $       dt
   5.000E-07
  $
  *DATABASE_HISTORY_SHELL
  $      id1       id2       id3       id4       id5       id6       id7       id8
           1
  $
  *DATABASE_HISTORY_SOLID
  $      id1       id2       id3       id4       id5       id6       id7       id8
           1
  $
  *DATABASE_GLSTAT
  $       dt
   5.000E-07
  $
  *DATABASE_NODOUT
  $       dt
   5.000E-07
  $
  *DATABASE_HISTORY_NODE
  $      id1       id2       id3       id4       id5       id6       id7       id8
           1       116       284       361       326
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Constraints
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  $$$$  Constrain shell nodes on the inner radius of the dome shells 
  $$$$  to the outer radius of the dome solids
  $
  *CONSTRAINED_SHELL_TO_SOLID
  $      nid      nsid
         326         1
  $
  *SET_NODE_LIST
  $      sid       da1       da2       da3       da4
           1
  $     nid1      nid2      nid3      nid4      nid5      nid6      nid7      nid8
         116       158       200       242       284
  $
  $
  *CONSTRAINED_SHELL_TO_SOLID
  $      nid      nsid
         327         2
  $
  *SET_NODE_LIST
  $      sid       da1       da2       da3       da4
           2
  $     nid1      nid2      nid3      nid4      nid5      nid6      nid7      nid8
         117       159       201       243       285
  $
  $
  *CONSTRAINED_SHELL_TO_SOLID
  $      nid      nsid
         328         3
  $
  *SET_NODE_LIST
  $      sid       da1       da2       da3       da4
           3
  $     nid1      nid2      nid3      nid4      nid5      nid6      nid7      nid8
         118       160       202       244       286
  $
  $
  *CONSTRAINED_SHELL_TO_SOLID
  $      nid      nsid
         329         4
  $
  *SET_NODE_LIST
  $      sid       da1       da2       da3       da4
           4
  $     nid1      nid2      nid3      nid4      nid5      nid6      nid7      nid8
         119       161       203       245       287
  $
  $
  *CONSTRAINED_SHELL_TO_SOLID
  $      nid      nsid
         330         5
  $
  *SET_NODE_LIST
  $      sid       da1       da2       da3       da4
           5
  $     nid1      nid2      nid3      nid4      nid5      nid6      nid7      nid8
         120       162       204       246       288
  $
  $
  *CONSTRAINED_SHELL_TO_SOLID
  $      nid      nsid
         331         6
  $
  *SET_NODE_LIST
  $      sid       da1       da2       da3       da4
           6
  $     nid1      nid2      nid3      nid4      nid5      nid6      nid7      nid8
         121       163       205       247       289
  $
  $
  *CONSTRAINED_SHELL_TO_SOLID
  $      nid      nsid
         332         7
  *SET_NODE_LIST
  $      sid       da1       da2       da3       da4
           7
  $     nid1      nid2      nid3      nid4      nid5      nid6      nid7      nid8
         122       164       206       248       290
  $
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Define Loads
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Load 4 segments with pressure at 5,308 psi for 2.0E-04 seconds
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  *LOAD_SEGMENT
  $     lcid        sf        at        n1        n2        n3        n4
           1 1.000E+00 0.000E+00         1         2         4         3
           1 1.000E+00 0.000E+00         2         5         7         4
           1 1.000E+00 0.000E+00         3         4         8         6
           1 1.000E+00 0.000E+00         4         7         9         8
  $
  *DEFINE_CURVE
  $     lcid      sidr      scla      sclo      offa      offo
           1
  $
  $                  a                   o
             0.000E+00           5.308E+04
             2.000E-04           5.308E+04
             2.010E-04           0.000E+00
             1.000E+00           0.000E+00
  $
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Define Parts and Materials
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  *PART
  $      pid       sid       mid     eosid      hgid    adpopt
  Brick
           1         1         1         0         0         0
  Brick2
           2         1         1         0         0         0
  Shell
           3         2         1         0         0         0
  $
  $
  $$$$$$  Materials
  $
  *MAT_ELASTIC
  $      mid        ro         e        pr        da        db
           1   2.00e-4  29.00e+6     0.330
  $
  $
  $$$$$$  Sections
  $
  *SECTION_SOLID
  $      sid    elform
           1         0
  $
  *SECTION_SHELL
  $      sid    elform      shrf       nip     propt   qr/irid     icomp
           2 
  $
  $       t1        t2        t3        t4      nloc
   1.576E-02 1.576E-02 1.576E-02 1.576E-02
  $
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Define Nodes and Elements
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  *END