*Heading Mullins Effect, Compressible Van der Waals Model, 2D Element Cyclic Uniaxial Tension ** Job name: Job-1 Model name: x_mmecdo2cut_vdw.inp ** ** PARTS ** *Part, name=PART-1 *End Part ** ** ASSEMBLY ** *Assembly, name=Assembly ** *Instance, name=PART-1-1, part=PART-1 *Node 1, 0., 0. 2, 0., 1. 3, 1., 0. 4, 1., 1. *Element, type=CPE4R 1, 1, 3, 4, 2 *Elset, elset=ONE 1, ** Region: (Section-1-ONE:ONE) *Elset, elset=_I1, internal 1, ** Section: Section-1-ONE *Solid Section, elset=_I1, material=TREL 1., *End Instance *Nset, nset=ALL, instance=PART-1-1, generate 1, 4, 1 *Nset, nset=BOT, instance=PART-1-1 1, 3 *Nset, nset=LEFT, instance=PART-1-1 1, 2 *Nset, nset=RIGHT, instance=PART-1-1 3, 4 *Nset, nset=TOP, instance=PART-1-1 3, 4 *Nset, nset=SET, instance=PART-1-1 3 *Nset, nset=FOUR, instance=PART-1-1 4 *EQUATION ** Since the S11 output is Cauchy or true stress, we need to ** determine the nominal stress for post-processing. ** Node 3 is tied to node 4 in dof 1 so that: ** Nominal stress (dof 1) = RF1 (@ node 4) / Original area ** (w/c is 1 x 1 = 1) 2 SET,1,1, FOUR,1,-1 *End Assembly ** ** MATERIALS ** *Material, name=TREL *Hyperelastic, van der waals 200., 10., 0.1, 0., 0.001 *Mullins effect 3., 100., 0.1 *Density 1. ** ---------------------------------------------------------------- ** ** STEP: Step-1 ** *Step, name=Step-1 UNIAXIAL TENSION *Dynamic, Explicit , 20., ** ** BOUNDARY CONDITIONS ** *Boundary BOT, 2, 2 LEFT, 1, 1 *Boundary, type=VELOCITY FOUR, 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=FOUR U,RF *End Step ** ---------------------------------------------------------------- ** ** STEP: Step-2 ** *Step, name=Step-2 *Dynamic, Explicit , 20., ** ** BOUNDARY CONDITIONS ** *Boundary, type=VELOCITY FOUR, 1, 1, -0.1 ** *End Step ** ---------------------------------------------------------------- ** ** STEP: Step-3 ** *Step, name=Step-3 *Dynamic, Explicit , 20., ** ** BOUNDARY CONDITIONS ** *Boundary, type=VELOCITY FOUR, 1, 1, 0.15 ** *End Step ** ---------------------------------------------------------------- ** ** STEP: Step-4 ** *Step, name=Step-4 *Dynamic, Explicit , 20., ** ** BOUNDARY CONDITIONS ** *Boundary, type=VELOCITY FOUR, 1, 1, -0.15 *End Step