Implicit — Welcome to LS-DYNA Examples

Home » Knowledge Base » Simple Examples » Examples Manual » Bolts » Bolt Type B » Implicit

This LS-DYNA example shows the implicit modeling of Bolt Type B.

Keywords

*BOUNDARY_PRESCRIBED_MOTION_SET
*BOUNDARY_SPC_SET
*COMMENT 
*CONTROL_IMPLICIT_EIGENVALUE
*CONTACT_AUTOMATIC_SINGLE_SURFACE_MORTAR_ID
*CONTACT_FORCE_TRANSDUCER_PENALTY_ID
*CONTROL_ACCURACY
*CONTROL_BULK_VISCOSITY
*CONTROL_CONTACT
*CONTROL_ENERGY
*CONTROL_HOURGLASS
*CONTROL_IMPLICIT_AUTO
*CONTROL_IMPLICIT_DYNAMICS
*CONTROL_IMPLICIT_GENERAL
*CONTROL_IMPLICIT_SOLUTION
*CONTROL_IMPLICIT_SOLVER
*CONTROL_MPP_IO_NODUMP
*CONTROL_OUTPUT
*CONTROL_RIGID
*CONTROL_SHELL
*CONTROL_SOLID
*CONTROL_SOLUTION
*CONTROL_TERMINATION
*CONTROL_TIMESTEP
*DEFINE_CURVE
*ELEMENT_BEAM
*ELEMENT_SHELL
*ELEMENT_SOLID
*END
*INCLUDE
*INITIAL_AXIAL_FORCE_BEAM
*KEYWORD
*MAT_NULL
*MAT_PLASTIC_KINEMATIC
*MAT_SPOTWELD
*NODE
*PARAMETER
*PARAMETER_EXPRESSION
*PART
*PART_CONTACT
*SECTION_BEAM
*SECTION_SHELL
*SECTION_SOLID
*SET_BEAM
*SET_NODE_LIST_TITLE
*SET_PART_ADD
*SET_PART_LIST
*TITLE

Reduced Input

$-----------------------------------------------------------------------------
$
$ Example provided by Dr. Nils Karajan (DYNAmore Corporation)
$
$ Date: Wed Apr 11 10:42:53 CET 2018
$
$ E-Mail: info@dynamore.de
$ Web: http://www.dynamore.de
$
$ Copyright, 2018 DYNAmore Corp.
$ Copying for non-commercial usage allowed if
$ copy bears this notice completely.
$
$-------------------------------------------------------------------------------
$
$ Run file as is. It was tested with LS-DYNA R9.2 with double precision.
$ 
$-------------------------------------------------------------------------------
$# Units: kg / mm / ms / kN / GPa / kN-mm
$-------------------------------------------------------------------------------
$
*KEYWORD
*PARAMETER
$#   prmr1      val1     prmr2      val2     prmr3      val3     prmr4      val4
$
$--- Simulation time 
R     tEnd       5.0
$
$--- Pre-force in beams / Pre-stress in solids (cross sect 74.9859 mm^2)
R bltForce      28.8
R bltStrss    0.3841
$
$--- Loading of the connection
R    distX       5.0
$
*PARAMETER_EXPRESSION
$#   prmr1                                                           expression1
$
$--- Plot intervals
R   dtPlot       TEnd/100.0
R  dtAscii       TEnd/1000.0
$
$--- Load application times
R dtPreStr       TEnd/5.0
R    tLoad       TEnd/2.0
$
$--- Implicit time integration
R      dt0       dtPreStr/20.
R    dtMin       dtPreStr/30.
R    dtMax       dtPreStr/10.
R   tDyDth       dtPreStr*2./3.
R   tDyBur       dtPreStr
R   tDyBir       tLoad
$
$--- Infinity time
R   tInfty       tEnd*1.01
$
$ =============
$ INCLUDE cards
$ =============
$#                                                                         title
Bolted Connection Type b
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ INITIAL cards
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
*INITIAL_AXIAL_FORCE_BEAM
$#    bsid      lcid     scale     kbend     
       100       100                   0
$
*SET_BEAM
$#     sid   
       100
$#      k1        k2        k3        k4        k5        k6        k7        k8
   1000000                                                                      
$
*DEFINE_CURVE
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp     lcint
       100           &dtPreStr &bltForce                                        
$#                a1                  o1  
                 0.0                 0.0
                 1.0                 1.0
$
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ BOUNDARY cards
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
*BOUNDARY_PRESCRIBED_MOTION_SET
$#    nsid       dof       vad      lcid        sf       vid     death     birth
        20         1         2        11                                        
*DEFINE_CURVE
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp     lcint
        11                                                                      
$#                a1                  o1  
                 0.0                 0.0
              &tLoad                 0.0
               &tEnd              &distX
             &tInfty              &distX
$
$---
$
*BOUNDARY_SPC_SET
$#    nsid       cid      dofx      dofy      dofz     dofrx     dofry     dofrz
        10                   1         1         1         1         1         1
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ PART / SECTION / MATERIALS
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$--- Bolt
$
*PART_CONTACT
$#                                                                         title
Bolt
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
       100       100       100                                                  
$#      fs        fd        dc        vc      optt       sft       ssf   
                                              10.0                    
*SECTION_BEAM
$#   secid    elform      shrf   qr/irid       cst     scoor       nsm   
       100         9       1.0         2         1                    
$#     ts1       ts2       tt1       tt2     print     
      10.0      10.0                              
$
*MAT_SPOTWELD
$#     mid        ro         e        pr      sigy        eh        dt     tfail
       100   7.85e-6     210.0       0.3      0.64       9.5                    
$#   efail       nrr       nrs       nrt       mrr       mss       mtt        nf
                                                                                
$
$--- Bolt Head and Nut
$
*PART_CONTACT
$#                                                                         title
Bolt Head                                                                       
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
       110       110       110                                                  
$#      fs        fd        dc        vc      optt       sft       ssf   
                                               6.4                    
*PART_CONTACT
$#                                                                         title
Bolt Nut                                                                        
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
       120       110       110                                                  
$#      fs        fd        dc        vc      optt       sft       ssf   
                                               6.4                    
$
*SECTION_SHELL
$#   secid    elform      shrf       nip     propt   qr/irid     icomp     setyp
       110        16       1.0         5       3.0                              
$#      t1        t2        t3        t4      nloc     marea      idof    edgset
       6.4
*MAT_PLASTIC_KINEMATIC
$#     mid        ro         e        pr      sigy      etan      beta
       110   7.85e-6     210.0       0.3      0.64       9.5          
$#     src       srp        fs        vp
                                     1.0
$
$--- Contact Beams
$
*PART_CONTACT
$#                                                                         title
Contact Beam Bolt (imp only)
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
       101       101       310                                                  
$#      fs        fd        dc        vc      optt       sft       ssf   
                                              10.0                    
*SECTION_BEAM
$
$#   secid    elform      shrf   qr/irid       cst     scoor       nsm   
       101         1       1.0         2         1                    
$#     ts1       ts2       tt1       tt2     nsloc     ntloc     
      10.0      10.0                                        
$
*PART_CONTACT
$#                                                                         title
Contact Beams Plates
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
       310       310       310                                                  
$#      fs        fd        dc        vc      optt       sft       ssf   
                                               0.5                    
*SECTION_BEAM
$
$#   secid    elform      shrf   qr/irid       cst     scoor       nsm   
       310         1       1.0         2         1                    
$#     ts1       ts2       tt1       tt2     nsloc     ntloc     
       0.5       0.5                                        
$
*MAT_NULL
$#     mid        ro        pc        mu     terod     cerod        ym        pr
       310   7.85e-6                                             210.0       0.3
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ CONTROL cards for implicit
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
*CONTROL_IMPLICIT_GENERAL
$#  imflag       dt0    imform      nsbs       igs     cnstn      form    zero_v
         1      &dt0         2         1         2                              
$
$--- Switch on/off dynamic effects over time
$
*CONTROL_IMPLICIT_DYNAMICS
$#   imass     gamma      beta    tdybir    tdydth    tdybur     irate     alpha
       -42      0.60   0.38000
*DEFINE_CURVE
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp     lcint
        42                                                                      
$#                a1                  o1  
                 0.0                 1.0
             &tDyDth                 1.0
             &tDyBur                 0.0
        &tDyBir*0.99                 0.0
             &tDyBir                 1.0
             &tInfty                 1.0
$
$--- Set Key points in time to reach exactly
$
*CONTROL_IMPLICIT_AUTO
$#   iauto    iteopt    itewin     dtmin     dtmax     dtexp     kfail    kcycle
         1        40        10    &dtMin       -24                        
*DEFINE_CURVE
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp     lcint
        24                                                                      
$#                a1                  o1  
             &tDyDth              &dtMax
             &tDyBur              &dtMax
             &tDyBir              &dtMax
$
$---
$
*CONTROL_IMPLICIT_SOLUTION
$#  nsolvr    ilimit    maxref     dctol     ectol     rctol     lstol    abstol
        12         6        12    1.0e-3    1.0e-2                       1.0e-20
$#   dnorm    diverg     istif   nlprint    nlnorm   d3itctl     cpchk
         1                             3         4         1          
$#  arcctl    arcdir    arclen    arcmth    arcdmp    arcpsi    arcalf    arctim
                                                                                
$#   lsmtd     lsdir      irad      srad      awgt      sred
         4                                                   
*CONTROL_IMPLICIT_SOLVER
$#  lsolvr    lprint     negev     order      drcm    drcprm   autospc   autotol
                   1                   2                             2          
*COMMENT *CONTROL_IMPLICIT_EIGENVALUE
$#    neig    center     lflag    lftend     rflag    rhtend    eigmth    shfscl
        20
$#   solid     ibeam    ishell   itshell    mstres    evdump   mstrscl
                                                                  0.01
$
$---
$
*CONTROL_ACCURACY
$#     osu       inn    pidosu      iacc
         1         4                   1
$
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ CONTROL cards - mostly recommendations for explicit crash applications
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
*CONTROL_SOLUTION
$#    soln       nlq     isnan     lcint     
                             1      1001
*CONTROL_ENERGY
$#    hgen      rwen    slnten     rylen     
         2         2         2         2
*CONTROL_HOURGLASS
$#     ihq        qh  
         8       0.1
*CONTROL_BULK_VISCOSITY
$#      q1        q2      type     btype     
       1.5      0.06        -1          
*CONTROL_CONTACT
$#  slsfac    rwpnal    islchk    shlthk    penopt    thkchg     orien    enmass
      0.80       1.0         2                   1                              
$#  usrstr    usrfrc     nsbcs    interm     xpene     ssthk      ecdt   tiedprj
                                                                                
$#   sfric     dfric       edc       vfc        th     th_sf    pen_sf      
                                                                      
$#  ignore    frceng   skiprwg    outseg   spotstp   spotdel   spothin       
         2                             1         2         1          
$#    isym    nserod    rwgaps    rwgdth     rwksf      icov    swradf    ithoff
                             1                 1.0                              
$#  shledg    pstiff    ithcnt    tdcnof     ftall    unused    shltrw      
                                                                      
*CONTROL_RIGID
$#     lmf      jntf    orthmd     partm    sparse    metalf    plotel     rbsms
                                                                                
$#  norbic   gjadstf   gjadvsc   tjadstf   tjadvsc
                    
*CONTROL_SHELL
$#  wrpang     esort     irnxx    istupd    theory       bwc     miter      proj
      30.0         1                   1                                       1
$# rotascl    intgrd    lamsht    cstyp6    thshel      
                                                  
$# psstupd   sidt4tu     cntco    itsflg    irquad    w-mode   stretch      icrq
                                                                                
$#  nfail1    nfail4   psnfail    keepcs     delfr   drcpsid    drcprm   intperr
                   1                                               1.0          
*CONTROL_SOLID
$#   esort   fmatrix   niptets    swlocl    psfail   t10jtol    icohed    tet13k
         1                             1                                        
$#   pm1     pm2     pm3     pm4     pm5     pm6     pm7     pm8     pm9    pm10
                                                                                
*CONTROL_OUTPUT
$#   npopt    neecho    nrefup    iaccop     opifs    ipnint    ikedit    iflush
         1         3                   1                                        
$#   iprtf    ierode   tet10s8    msgmax    ipcurv      gmdt   ip1dblt      eocs
                   1               10000                                        
$#   tolev    newleg    frfreq     minfo    solsig    msgflg    cdetol      
                                                                      
*CONTROL_MPP_IO_NODUMP
$
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ CONTACT cards
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
*CONTACT_AUTOMATIC_SINGLE_SURFACE_MORTAR_ID
$#     cid                                                                 title
     10001 Global Contact
$
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
     10101                   2                                                  
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
      0.15      0.10       1.0                                              
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
                                                                                
$#    soft    sofscl    lcidab    maxpar     sbopt     depth     bsort    frcfrq
                                                                               1
$#  penmax    thkopt    shlthk     snlog      isym     i2d3d    sldthk    sldstf
                                                                                
$#    igap    ignodprfac/mpadtstif/mpar2   unused     unused    flangl   cid_rcf
         3         2                                                            
$#   q2tri    dtpchk     sfnbr    fnlscl   dnlscl       tcso    tiedid    shledg
                                                                                
$#  sharec    cparm8    ipback     srnde    fricsf      icor     ftorq    region
                                                                                
$#  pstiff   ignroff   beam-cs       
                              
$
*CONTACT_FORCE_TRANSDUCER_PENALTY_ID
$#     cid                                                                 title
     10002 Contact Force Between Sheets
$
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
        10        20         3         3                                        
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
                                                                                
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
                                                                                
$ =============
$ CONTROL cards
$ =============
$
*CONTROL_TERMINATION
$#  endtim    endcyc     dtmin    endeng    endmas     nosol     
     &tEnd                                                        
*CONTROL_TIMESTEP
$#  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st
                                                                                
$#  dt2msf   dt2mslc     imscl    unused    unused     rmscl     
                                                                               
*END

Figures

Bolt_B_Impl_NoFringe

Bolt_B_Impl_Fringe

Animated Result No Fringe

Bolt_B_Impl_NoFringe_Movie

Animated Result Fringe

Bolt_B_Impl_Fringe_Movie

Download