Implicit
This LS-DYNA example shows the implicit modeling of Bolt Type A.
https://www.dynaexamples.com/show-cases/bolts/typea/implicit
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Implicit
This LS-DYNA example shows the implicit modeling of Bolt Type A.
*BOUNDARY_PRESCRIBED_MOTION_SET *BOUNDARY_SPC_SET *COMMENT *CONTROL_IMPLICIT_EIGENVALUE *CONSTRAINED_NODAL_RIGID_BODY *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_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
$-----------------------------------------------------------------------------
$
$ 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 a
$
$---+----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
$
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ CONSTRAINED cards (bolt head and nut)
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
*CONSTRAINED_NODAL_RIGID_BODY
$# pid cid nsid pnode iprt drflag rrflag
110 110 1
*CONSTRAINED_NODAL_RIGID_BODY
$# pid cid nsid pnode iprt drflag rrflag
120 120
$
$---+----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
*SET_PART_ADD
$# sid da1 da2 da3 da4 solver
10101
$# spid1 spid2 spid3 spid4 spid5 spid6 spid7 spid8
10001 10100
$
$---
$
*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