*KEYWORD *TITLE template for metal forming $ $ $ use this input deck as a guide for setting up metal forming $ simulations with LS-DYNA $ $------------------------------------------------------------------------------- $ set the termination time, mass scaling time step size $ *CONTROL_TERMINATION $ ENDTIM x.x $ *CONTROL_TIMESTEP $ DTINIT TSSFAC ISDO TSLIMT DT2MS 0.0 0.0 0 0.0 x.x $ $------------------------------------------------------------------------------- $ The adaptive frequency should be small enough to adapt every 2-4 mm of punch $ travel. $ *CONTROL_ADAPTIVE $ ADPFREQ ADPTOL ADPOPT MAXLVL TBIRTH TDEATH LCADP IOFLAG x.x 4.000 2 3 0.0 0.0 0 0 $ ADPSIZE ADPASS IREFLG ADPENE 2.00 1 0 3.0000 $ *CONTROL_HOURGLASS $ IHQ QH 4 0.100 $ *CONTROL_SHELL $ WRPANG ITRIST IRNXX ISTUPD THEORY BWC MITER 20.000 2 -1 1 2 2 1 $ *CONTROL_CONTACT $ SLSFAC RWPNAL ISLCHK SHLTHK PENOPT THKCHG ORIEN 0.010 2 1 4 2 3 $ USRSTR USRFAC NSBCS INTERM XPENEN 0 0 0 0 0.0 $ $------------------------------------------------------------------------------- $ select an output frequency, additional databases. Add strain to plot file. $ *DATABASE_RCFORC $ dt x.x *DATABASE_BNDOUT $ dt x.x *DATABASE_BINARY_D3PLOT $ dt x.x *DATABASE_EXTENT_BINARY $ neiph neips maxint strflg sigflg epsflg rltflg engflg 0 0 0 1 0 0 0 0 $ cmpflg ieverp beamip dcomp shge stssz n3thdt 0 0 0 0 0 0 0 $ $------------------------------------------------------------------------------- $ $ Create the output file "dynain", including part #1 only. Add three $ nodes for constraining rigid body motion during springback. Code 7 = xyz, $ code 5 = yz, and code 3 = z-direction constrained. $ *INTERFACE_SPRINGBACK_DYNA3D $ psid 6 $ nid code xxx 7 xxx 5 xxx 3 $ *SET_PART_LIST 6 1 $ $------------------------------------------------------------------------------- $ enter material properties in parenthesis (xx) $ $ sheet *MAT_TRANSVERSELY_ANISOTROPIC_ELASTIC_PLASTIC 1 (rho) (young's) (pr) (sigy) 0.0000000 (r) 1 $ *DEFINE_CURVE 1 $ enter material stress vs PLASTIC strain x.x x.x x.x x.x x.x x.x $ $ $ die *MAT_RIGID 2 (rho) (young's) (pr) 1.0000000 7.0000000 7.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 $ $ punch *MAT_RIGID 3 (rho) (young's) (pr) 1.0000000 4.0000000 7.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 $ $ binder *MAT_RIGID 4 (rho) (young's) (pr) 1.0000000 4.0000000 7.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 $ $------------------------------------------------------------------------------- $ enter initial sheet thickness $ *SECTION_SHELL $ SECID ELFORM SHRF NIP PROPT QR/IRID ICOMP 1 16 0.0000000 7.0000000 0.0000000 0.0000000 0 (thick) (thick) (thick) (thick) $ *SECTION_SHELL $ SECID ELFORM SHRF NIP PROPT QR/IRID ICOMP 2 2 0.0000000 1.0000000 0.0000000 0.0000000 0 1.0000000 1.0000000 1.0000000 1.0000000 0.0000000 $ *SECTION_SHELL $ SECID ELFORM SHRF NIP PROPT QR/IRID ICOMP 3 2 0.0000000 1.0000000 0.0000000 0.0000000 0 1.0000000 1.0000000 1.0000000 1.0000000 0.0000000 $ *SECTION_SHELL $ SECID ELFORM SHRF NIP PROPT QR/IRID ICOMP 4 2 0.0000000 1.0000000 0.0000000 0.0000000 0 1.0000000 1.0000000 1.0000000 1.0000000 0.0000000 $------------------------------------------------------------------------------- $ *PART sheet $ PID SID MID EOSID HGID GRAV ADPOPT TMID 1 1 1 0 0 0 1 0 *PART die $ PID SID MID EOSID HGID GRAV ADPOPT TMID 2 2 2 0 0 0 0 0 *PART punch $ PID SID MID EOSID HGID GRAV ADPOPT TMID 3 3 3 0 0 0 0 0 *PART binder $ PID SID MID EOSID HGID GRAV ADPOPT TMID 4 4 4 0 0 0 0 0 $ $------------------------------------------------------------------------------- $ punch displacement $ *BOUNDARY_PRESCRIBED_MOTION_RIGID $ pid dir vad lcid scale 3 3 2 2 1.00 $ *DEFINE_CURVE 2 $ time punch disp x.x x.x x.x x.x x.x x.x $ $------------------------------------------------------------------------------- $ binder load $ *LOAD_RIGID_BODY $ pid dir lcid scale 4 3 3 1.00 $ *DEFINE_CURVE 3 $ time binder load x.x x.x x.x x.x x.x x.x $ $------------------------------------------------------------------------------- $ enter friction coefficient "fs" for each contact interface: $ $ sheet vs die $ *CONTACT_FORMING_ONE_WAY_SURFACE_TO_SURFACE $ ssid msid sstyp mstyp sboxid mboxid spr mpr 1 2 3 3 0 0 0 0 $ fs fd dc vc vdc penchk bt dt x.x 0.0 0.0 0.0 20.0 0 0.0 0.0 $ sfs sfm sst mst sfst sfmt fsf vsf 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 $ $ $ sheet vs punch $ *CONTACT_FORMING_ONE_WAY_SURFACE_TO_SURFACE $ ssid msid sstyp mstyp sboxid mboxid spr mpr 1 3 3 3 0 0 0 0 $ fs fd dc vc vdc penchk bt dt x.x 0.0 0.0 0.0 20.0 0 0.0 0.0 $ sfs sfm sst mst sfst sfmt fsf vsf 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 $ $ $ sheet vs binder $ *CONTACT_FORMING_ONE_WAY_SURFACE_TO_SURFACE $ ssid msid sstyp mstyp sboxid mboxid spr mpr 1 4 3 3 0 0 0 0 $ fs fd dc vc vdc penchk bt dt x.x 0.0 0.0 0.0 20.0 0 0.0 0.0 $ sfs sfm sst mst sfst sfmt fsf vsf 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 $ $------------------------------------------------------------------------------- $ *NODE $ $ enter or include node data $ $ *ELEMENT_SHELL $ $ enter or include element data $ use part numbering: 1=sheet, 2=die, 3=punch, 4=binder $ *END