*HEADING STATIC INDENTATION OF ELASTOMERIC FOAM WITH A HEMISPHERICAL PUNCH AND CAX4 elements ******************************************************** ** Measure the punch reaction force during the following 2 steps: ** Step 1: Displace punch 200 mm downwards. ** Step 2: ALE ** Step 2: Displace punch 50mm further down. ** (Units: N, mm, sec) ******************************************************** *PREPRINT,MODEL=YES,HISTORY=NO *RESTART,WRITE *NODE 1, 0., 300. 19, 0., 0. 481, 300., 300. 499, 300., 0. 601, 600., 300. 619, 600., 0. *NSET,NSET=N1 1, *NSET,NSET=N19 19, *NSET,NSET=N481 481, *NSET,NSET=N499 499, *NFILL,NSET=TOP,BIAS=.85 N1,N481,12,40 *NGEN,NSET=TOP 481,601,40 *NFILL,NSET=BOT,BIAS=.85 N19,N499,12,40 *NGEN,NSET=BOT 499,619,40 *NFILL,NSET=ALLN,BIAS=.9 TOP,BOT,9,2 *NSET,NSET=CENTER,GENERATE 1,19 *ELEMENT,TYPE=CAX4 1,3,43,41,1 *ELGEN,ELSET=FOAM 1,15,40,10,9,2,1 *ELSET,ELSET=CENT,GENERATE 1,141,10 61,69 *ELSET,ELSET=ETOP,GENERATE 1,101,10 *ELSET,ELSET=UPPER,GEN 1,141,10 *NODE,NSET=PUNCH 1000,0.,510. *SOLID SECTION,ELSET=FOAM,MATERIAL=HYPER *MATERIAL, NAME=HYPER *HYPERELASTIC 80.0, 20.0, 0.01 *SOLID SECTION,ELSET=FOAM,MATERIAL=FOAM *ORIENTATION,NAME=RECT,SYSTEM=RECTANGULAR 1.,0.,0., 0.,1.,0. 3,0. **************************************************** *MATERIAL,NAME=FOAM *HYPERFOAM,N=2,TEST DATA INPUT,POISSON=0.0,MODULI=INSTANTANEOUS *UNIAXIAL TEST DATA -.0217, -.05 -.0317, -.10 -.0367, -.15 -.0402, -.20 -.0433, -.25 -.0467, -.30 -.0504, -.35 -.0542, -.40 -.0604, -.45 -.0668, -.50 -.0759, -.55 -.0909, -.60 -.1083, -.65 -.1410, -.70 -.1933, -.75 -.2896, -.80 *SIMPLE SHEAR TEST DATA .0140, .08, .0046 .0334, .16, .0166 .0533, .24, .0366 .0853, .32, .0573 .1280, .40, .0817 .1653, .48, .1098 .2080, .56, .1394 .2560, .64, .1666 .2987, .72, .1904 *VISCOELASTIC,TIME=RELAXATION TEST DATA *SHEAR TEST DATA,SHRINF=0.5000 1.0000, 0.0001 0.9695, 0.001 0.9417, 0.002 0.8722, 0.005 0.7913, 0.010 0.7043, 0.020 0.6233, 0.050 0.5736, 0.100 0.5271, 0.200 0.5013, 0.500 0.5000, 1.000 *DENSITY 1.E-11, **************************************************** *BOUNDARY BOT,1,2 CENTER,1 1000,1 1000,6 **************************************************** *SURFACE,TYPE=SEGMENTS,NAME=IMPACTOR START,141.42,641.42 CIRCL,-1.,300.,0.,500. *SURFACE,TYPE=ELEMENT,NAME=TARGET ETOP,S3 *RIGID BODY,REF NODE=1000,ANALYTICAL SURFACE=IMPACTOR *CONTACT PAIR,INTERACTION=ROUGH TARGET,IMPACTOR *SURFACE INTERACTION, NAME=ROUGH **************************************************** *STEP,NLGEOM,INC=200,AMPLITUDE=RAMP Displace punch 200 mm downwards. *STATIC .0015,1.,,.05 *BOUNDARY 1000,2,2,-200. *PRINT,CONTACT=YES *CONTACT CONTROLS,FRICTION ONSET=DELAY *CONTACT PRINT,SLAVE=TARGET **CONTACT FILE,SLAVE=TARGET,FREQUENCY=10 *ENERGY PRINT,FREQUENCY=5 *ENERGY PRINT,ELSET=FOAM,FREQ=5 *ENERGY PRINT,ELSET=CENT,FREQ=5 *EL PRINT,FREQUENCY=50,ELSET=CENT S, E, *NODE PRINT,FREQUENCY=25 U, RF, *OUTPUT,FIELD,FREQUENCY=4 *NODE OUTPUT U,V,A,RF *ELEMENT OUTPUT S,LE,NE *CONTACT OUTPUT,VARIABLE=PRESELECT,NSET=TOP *OUTPUT,HISTORY,FREQUENCY=1 *NODE OUTPUT,NSET=PUNCH U,V,A,RF *ELEMENT OUTPUT,ELSET=CENT S,LE,NE *ENERGY OUTPUT,VARIABLE=PRESELECT *CONTROLS, PARAMETERS=FIELD, FIELD=DISPLACEMENT 1e-04,1e-03 *END STEP **************************************************** *STEP,NLGEOM,INC=200,AMPLITUDE=RAMP ADAPTIVE MESHING *STATIC 1.0, 1.0 *PRINT,CONTACT=YES *CONTACT CONTROLS,FRICTION ONSET=DELAY *CONTACT PRINT,SLAVE=TARGET **CONTACT FILE,SLAVE=TARGET,FREQUENCY=10 *ENERGY PRINT,FREQUENCY=5 *ENERGY PRINT,ELSET=FOAM,FREQ=5 *ENERGY PRINT,ELSET=CENT,FREQ=5 *EL PRINT,FREQUENCY=50,ELSET=CENT S, E, *NODE PRINT,FREQUENCY=25 U, RF, *OUTPUT,FIELD,FREQUENCY=4 *NODE OUTPUT U,V,A,RF *ELEMENT OUTPUT S,LE,NE *CONTACT OUTPUT,VARIABLE=PRESELECT,NSET=TOP *OUTPUT,HISTORY,FREQUENCY=1 *NODE OUTPUT,NSET=PUNCH U,V,A,RF *ELEMENT OUTPUT,ELSET=CENT S,LE,NE *ENERGY OUTPUT,VARIABLE=PRESELECT *Adaptive Mesh, Elset=Foam, Frequency=1, Mesh Sweeps=6 *END STEP ********************************************************* *STEP,NLGEOM,INC=200,AMPLITUDE=RAMP Displace punch 50 mm further downwards. *STATIC 0.1, 1.0 *BOUNDARY 1000,2,2,-250. *PRINT,CONTACT=YES *CONTACT CONTROLS,FRICTION ONSET=DELAY *CONTACT PRINT,SLAVE=TARGET **CONTACT FILE,SLAVE=TARGET,FREQUENCY=10 *ENERGY PRINT,FREQUENCY=5 *ENERGY PRINT,ELSET=FOAM,FREQ=5 *ENERGY PRINT,ELSET=CENT,FREQ=5 *EL PRINT,FREQUENCY=50,ELSET=CENT S, E, *NODE PRINT,FREQUENCY=25 U, RF, *OUTPUT,FIELD,FREQUENCY=4 *NODE OUTPUT U,V,A,RF *ELEMENT OUTPUT S,LE,NE *CONTACT OUTPUT,VARIABLE=PRESELECT,NSET=TOP *OUTPUT,HISTORY,FREQUENCY=1 *NODE OUTPUT,NSET=PUNCH U,V,A,RF *ELEMENT OUTPUT,ELSET=CENT S,LE,NE *ENERGY OUTPUT,VARIABLE=PRESELECT *END STEP