*Heading
Mullins effect, Compressible OGDEN model with viscoelasticity
3D element, cyclic uniaxial tension
Comparison of nominal stresses
** Job name: Job-1 Model name: mmecdo3cut_ogden
**
** PARTS
**
*Part, name=PART-1
*End Part
**
** ASSEMBLY
**
*Assembly, name=Assembly
**
*Instance, name=PART-1-1, part=PART-1
*Node
*Element, type=C3D8RH
*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,,10.0
** ----------------------------------------------------------------
**
** STEP: Step-1
**
*Step, name=Step-1, nlgeom, inc=20
UNIAXIAL TENSION
*visco
1., 20.,
**
** BOUNDARY CONDITIONS
**
*Boundary
FACE1, 3, 3
FACE3, 2, 2
FACE6, 1, 1
TWO, 1, 1, 2.
**
** HISTORY OUTPUT
**
*Output, history, frequency=10
***Output, history, frequency=1
**
*Element Output, elset=PART-1-1.ONE
SENER, S, DMENER
*Node Output, nset=TWO
U,RF,
*El Print, freq=999999
*Node Print, freq=999999
*End Step
** ----------------------------------------------------------------
**
** STEP: Step-2
**
*Step, name=Step-2, nlgeom
*visco
1., 20.,
**
** BOUNDARY CONDITIONS
**
*Boundary, op=NEW
FACE1, 3, 3
FACE3, 2, 2
FACE6, 1, 1
TWO, 1, 1
*End Step
** ----------------------------------------------------------------
**
** STEP: Step-3
**
*Step, name=Step-3, nlgeom
*visco
1., 20.,
**
** BOUNDARY CONDITIONS
**
*Boundary, op=NEW
FACE3, 2, 2
FACE6, 1, 1
TWO, 1, 1, 3.
FACE1, 3, 3
*End Step
** ----------------------------------------------------------------
**
** STEP: Step-4
**
*Step, name=Step-4, nlgeom
*visco
1., 20.,
**
** BOUNDARY CONDITIONS
**
*Boundary, op=NEW
FACE1, 3, 3
FACE3, 2, 2
FACE6, 1, 1
TWO, 1, 1
**
*End Step