*Heading
Mullins effect, Compressible OGDEN model
3D element, cyclic uniaxial tension
Comparison of nominal stresses
** Job name: Job-1 Model name: x_mmecdo3cut_ogden.inp
**
** PARTS
**
*Part, name=PART-1
*End Part
**
** ASSEMBLY
**
*Assembly, name=Assembly
**
*Instance, name=PART-1-1, part=PART-1
*Node
*Element, type=C3D8R
*Elset, elset=ONE
** Region: (Section-1-ONE:ONE)
*Elset, elset=_I1, internal
** Section: Section-1-ONE
*Solid Section, elset=_I1, material=OGDEN
1.,
*End Instance
*Nset, nset=ALL, instance=PART-1-1, generate
*Nset, nset=FACE1, instance=PART-1-1, generate
*Nset, nset=FACE2, instance=PART-1-1, generate
*Nset, nset=FACE3, instance=PART-1-1
*Nset, nset=FACE4, instance=PART-1-1
*Nset, nset=FACE5, instance=PART-1-1
*Nset, nset=FACE6, instance=PART-1-1
*Nset, nset=SET, instance=PART-1-1
*NSET, NSET=TWO, instance=PART-1-1
*EQUATION
** Since the S11 output is Cauchy or true stress, we need to
** determine the nominal stress for post-processing.
** Nodes 6,7,3 are tied to node 2 in dof 1 so that:
** Nominal stress (dof 1) = RF1 (@ node 2) / Original area
** (w/c is 1 x 1 = 1)
SET,1,1, PART-1-1.2,1,-1
*End Assembly
**
** MATERIALS
**
*Material, name=OGDEN
*HYPERELASTIC,N=2,OGDEN
160.,2.,40.,-2.,0.001
*Mullins effect
*viscoelastic,time=prony
0.5,,3.
*Density
1.
*Amplitude, name=loading
0.,0.,1.,1.
*Amplitude, name=unloading
0.,1.,1.,0.
** ----------------------------------------------------------------
**
** STEP: Step-1
**
*Step, name=Step-1
UNIAXIAL TENSION
*Dynamic, Explicit
, 20.,
**
** BOUNDARY CONDITIONS
**
*Boundary
FACE1, 3, 3
FACE3, 2, 2
FACE6, 1, 1
*Boundary, type=VELOCITY
TWO, 1, 1, 0.1
**
** HISTORY OUTPUT
**
*Output, history, time interval=5
**
*Energy output
allke, allse
*Element Output, elset=PART-1-1.ONE
SENER, S, DMENER
*Node Output, nset=TWO
U,RF,
*End Step
** ----------------------------------------------------------------
**
** STEP: Step-2
**
*Step, name=Step-2
*Dynamic, Explicit
, 20.,
**
** BOUNDARY CONDITIONS
**
*Boundary, type=VELOCITY
TWO, 1, 1, -0.1
*End Step
** ----------------------------------------------------------------
**
** STEP: Step-3
**
*Step, name=Step-3
*Dynamic, Explicit
, 20.,
**
** BOUNDARY CONDITIONS
**
*Boundary,type=VELOCITY
TWO, 1, 1, 0.15
*End Step
** ----------------------------------------------------------------
**
** STEP: Step-4
**
*Step, name=Step-4
*Dynamic, Explicit
, 20.,
**
** BOUNDARY CONDITIONS
**
*Boundary, type=VELOCITY
TWO, 1, 1, -0.15
**
*End Step