Hi all,
Can someone please help me in changing the .inp file below to allow for explicit analysis. I want to directly compare the results from the two simulations. I realise some of the commands in this .inp file are old and almost obsolete, but all I want to change is the step type. What I have tried both of the following:
1. Changed the lines
*STEP, NAME=close_roller, nlgeom=YES, inc=500 *STATIC 0.05, 1.000000, 1e-06, 1.000000
to
*Step, name=close_roller, nlgeom=YES *Dynamic, Explicit , 0.400080 *Bulk Viscosity 0.06, 1.2 ** Mass Scaling: Semi-Automatic ** sheet *Fixed Mass Scaling, elset=E_BAR_0, factor=2000.
2. Changed the lines
*STEP, NAME=rotate_roller, nlgeom=YES, inc=2000 *STATIC 0.05, 0.400080, 1e-06, 0.400080
to
*Step, name=rotate_roller, nlgeom=YES *Dynamic, Explicit , 0.400080
but I keep getting the error ``KEYWORD CARDS FOR STEP DEPENDENT INPUT MUST APPEAR AFTER THE FIRST *STEP CARD." The .inp file runs fine in the implicit case
** Roll Set-up Specifications ** **note: delh is gauge half-reduction(m), r: reduction (%), U: roller surface velocity, h0 is half-thickness(m) **Nx=180; **Ny=18; **Nz=0; **h0=0.05; **h1=0.04; **r=0.2; **delh_tot=0.01; **delh_top=0.010000000000000002; **delh_bot=0.010000000000000002; **R_top=0.13; **R_bot=0.13; **fr_top=0.3; **fr_bot=0.3; **U_top=1.24974997499; **U_bot=1.24974997499; **ms=2; **width=0; **length=1.0; **l_rollbite=0.05; **delta=1.0; **refNode=150000000; **materialproperties= ** elasticity=(1e+11, 0.3) ** plasticity= ** (173205080.757, 0) ** density=2700.0 **.INP File Starts Here ** *HEADING Plate Rolling: Plane strain, ABAQUS/Standard (Analytical rigid surfaces) SECTION CONTROLS USED (HOURGLASS=STIFFNESS)) Units - N, m ** Workpiece Geometry - 2*h0 thickness** *NODE 1, 0.05, -0.05 3259, 0.05, 0.05 181, -0.95, -0.05 3439, -0.95, 0.05 *NGEN, NSET=N_SURF_0 1, 181, 1 *NGEN, NSET=N_SURF_1 3259, 3439, 1 *NFILL,NSET=N_BAR_0 N_SURF_0, N_SURF_1, 18, 181 *ELEMENT, TYPE=CPE4R, ELSET=E_BAR_0 1, 1, 182, 183, 2 *ELGEN, ELSET=E_BAR_0 1, 180, 1, 1, 18, 181, 180 *ELSET,ELSET=E_SURF_0, GEN 1, 180, 1 *ELSET,ELSET=E_SURF_1, GEN 3061, 3240, 1 ** Define bar material behaviour ** *SOLID SECTION, ELSET=E_BAR_0, MAT=WeakPerfectPlasticAprox, CONTROLS=SECT 1., *SECTION CONTROLS, NAME=SECT, HOURGLASS=STIFFNESS *MATERIAL,NAME=WeakPerfectPlasticAprox *ELASTIC 1e+11, 0.3 *PLASTIC 300000000.0,0 *DENSITY 2700.0, ** Top Roll ** *NODE,NSET=REFTOP 150000000, 0.0 , 0.181 *SURFACE, TYPE=ELEMENT, NAME=SURFTOP, REGION TYPE=SLIDING E_SURF_1,S2 *SURFACE, TYPE=SEGMENTS, NAME=RIGIDTOP START,0.0,0.051 CIRCL,-0.13,0.181,0.0,0.181 CIRCL,0.0,0.311,0.0,0.181 CIRCL,0.13,0.181,0.0,0.181 CIRCL,0.0,0.051,0.0,0.181 *RIGID BODY, REFNODE=150000000, ANALYTICAL SURFACE=RIGIDTOP ** Bottom Roll ** *NODE,NSET=REFBOT 150000001, 0.0 , -0.181 *SURFACE, TYPE=ELEMENT, NAME=SURFBOT, REGION TYPE=SLIDING E_SURF_0,S4 *SURFACE, TYPE=SEGMENTS, NAME=RIGIDBOT START,0.0,-0.311 CIRCL,-0.13,-0.181,0.0,-0.181 CIRCL,0.0,-0.051,0.0,-0.181 CIRCL,0.13,-0.181,0.0,-0.181 CIRCL,0.0,-0.311,0.0,-0.181 *RIGID BODY, REFNODE=150000001, ANALYTICAL SURFACE=RIGIDBOT ** RAMPING DEFINITION ** *AMPLITUDE, NAME=apply_squeeze_ramp, DEFINITION=SMOOTH STEP 0,0,0.8,1 *AMPLITUDE, NAME=apply_roll_ramp, DEFINITION=SMOOTH STEP 0,0,1,1 ** Surface Interactions ** *SURFACE INTERACTION,NAME=FRICTTOP *FRICTION 0.3, *SURFACE BEHAVIOR, pressure-overclosure=EXPONENTIAL 0.0001, 100. *SURFACE INTERACTION,NAME=FRICTBOT *FRICTION 0.3, *SURFACE BEHAVIOR, pressure-overclosure=EXPONENTIAL 0.0001, 100. *CONTACT PAIR,INTERACTION=FRICTTOP, type=SURFACE TO SURFACE SURFTOP,RIGIDTOP *CONTACT PAIR,INTERACTION=FRICTBOT, type=SURFACE TO SURFACE SURFBOT,RIGIDBOT ** STEPPING INFORMATION ** *STEP, NAME=close_roller, nlgeom=YES, inc=500 *STATIC 0.05, 1.000000, 1e-06, 1.000000 ** Top Roll Boundary Conditions ** *BOUNDARY 150000000,1,1 150000000,6,6 *BOUNDARY,TYPE=DISPLACEMENT, AMPLITUDE=apply_squeeze_ramp 150000000,2,2,-0.011 ** Bottom Roll Boundary Conditions** *BOUNDARY 150000001,1,1 150000001,6,6 *BOUNDARY,TYPE=DISPLACEMENT, AMPLITUDE=apply_squeeze_ramp 150000001,2,2,0.011 ** Request output data ** *OUTPUT, FIELD, variable=PRESELECT, NUMBER INTERVAL=3 CSTRESS *OUTPUT, FIELD, NUMBER INTERVAL=3 *ELEMENT OUTPUT S,PEEQ, *NODE OUTPUT U,RF, V, *OUTPUT,HISTORY, FREQUENCY=4, NAME=REFTOP *NODE OUTPUT,NSET=REFTOP RF2,RM3 *OUTPUT,HISTORY, FREQUENCY=4, NAME=REFBOT *NODE OUTPUT,NSET=REFBOT RF2,RM3 *OUTPUT,HISTORY, FREQUENCY=4 SSPEEQ,SSTORQ,SSFORC *INCREMENTATION OUTPUT *OUTPUT, HISTORY, VAR=PRESELECT *ENERGY OUTPUT, VAR=PRESELECT *ENDSTEP *STEP, NAME=rotate_roller, nlgeom=YES, inc=2000 *STATIC 0.05, 0.400080, 1e-06, 0.400080 ** Top Roll Boundary Conditions ** *BOUNDARY 150000000,1,1 *BOUNDARY,TYPE=DISPLACEMENT 150000000,2,2,-0.011 150000000,3,5 *BOUNDARY,TYPE=VELOCITY, AMPLITUDE=apply_roll_ramp 150000000,6,6,9.613461346115374 ** Bottom Roll Boundary Conditions** *BOUNDARY 150000001,1,1 *BOUNDARY,TYPE=DISPLACEMENT 150000001,2,2,0.011 150000001,3,5 *BOUNDARY,TYPE=VELOCITY, AMPLITUDE=apply_roll_ramp 150000001,6,6,-9.613461346115374 ** Request output data ** *OUTPUT, FIELD, variable=PRESELECT, NUMBER INTERVAL=20 CSTRESS *OUTPUT, FIELD, NUMBER INTERVAL=20 *ELEMENT OUTPUT S,PEEQ, *NODE OUTPUT U,RF, V, *OUTPUT,HISTORY, FREQUENCY=4, NAME=REFTOP *NODE OUTPUT,NSET=REFTOP RF2,RM3 *OUTPUT,HISTORY, FREQUENCY=4, NAME=REFBOT *NODE OUTPUT,NSET=REFBOT RF2,RM3 *OUTPUT,HISTORY, FREQUENCY=4 SSPEEQ,SSTORQ,SSFORC *INCREMENTATION OUTPUT *OUTPUT, HISTORY, VAR=PRESELECT *ENERGY OUTPUT, VAR=PRESELECT *ENDSTEP
Abaqus