*HEADING - INCOMPRESSIBLE HYPERELASTICITY (YEOH), UNIAXIAL TENSION *PREPRINT,MODEL=YES *NODE,NSET=ALL 1, 2,1. 3,1.,1., 4,0.,1., 5,0.,0.,1. 6,1.,0.,1. 7,1.,1.,1. 8,0.,1.,1. *NSET,NSET=FACE1 1,2,3,4 *NSET,NSET=FACE2 5,6,7,8 *NSET,NSET=FACE3 1,2,5,6 *NSET,NSET=FACE4 2 *NSET,NSET=FACE42 3,6,7 *NSET,NSET=FACE5 3 *NSET,NSET=FACE52 4,7,8 *NSET,NSET=FACE6 4,1,8,5 *EQUATION ** Since the S11 output is Cauchy or true stress, we need to ** determine the nominal stress for post-processing. ** Nodes 3,6,7 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) 2 FACE42,1,1, 2,1,-1 *EQUATION 2 FACE52,2,1,3,2,-1 *ELEMENT,TYPE=C3D8RH,ELSET=ONE 1,1,2,3,4,5,6,7,8 *SOLID SECTION,ELSET=ONE,MATERIAL=YEOH *MATERIAL,NAME=YEOH *HYPERELASTIC, YEOH 1.326,-0.326,0.1319,0.0,0.0,0.0 *MULLINS EFFECT,TEST DATA INPUT *UNIAXIAL TEST DATA 84.1596 , 2. 52.2701 , 1.9 32.8573 , 1.8 21.2217 , 1.7 14.1184 , 1.6 9.629 , 1.5 6.68741 , 1.4 4.70457 , 1.3 3.34586 , 1.2 2.41153 , 1.1 1.77447 , 1. 1.34805 , 900.E-03 1.06922 , 800.E-03 889.632E-03,700.E-03 770.943E-03,600.E-03 682.284E-03,500.E-03 598.896E-03,400.E-03 501.378E-03,300.E-03 375.081E-03,200.E-03 209.883E-03,100.E-03 0., 0. *UNIAXIAL TEST DATA 509.046 , 3. 309.315 , 2.85 190.291 , 2.7 120.689 , 2.55 79.0685 , 2.4 53.0791 , 2.25 36.1182 , 2.1 24.644 , 1.95 16.7047 , 1.8 11.1661 , 1.65 7.32415 , 1.5 4.71004 , 1.35 2.98974 , 1.2 1.91099 , 1.05 1.2786 , 900.E-03 934.216E-03, 750.E-03 751.632E-03, 600.E-03 631.26E-03 , 450.E-03 497.431E-03, 300.E-03 297.089E-03, 150.E-03 0., 0. *STEP,NLGEOM,INC=20 UNIAXIAL TENSION *STATIC,DIRECT 1.,20. *BOUNDARY,OP=NEW FACE1,3 FACE3,2 FACE6,1 FACE4,1,1,2. *EL PRINT,F=10 S, E, *NODE PRINT,F=10 U,RF *OUTPUT,HISTORY *NODE OUTPUT,NSET=FACE4 U,RF *ELEMENT OUTPUT,ELSET=ONE S, DMENER, SENER, *END STEP *STEP,NLGEOM,INC=20 *STATIC,DIRECT 1.,20. *BOUNDARY,OP=NEW FACE1,3 FACE3,2 FACE6,1 FACE4,1,1,0.0 *END STEP *STEP,NLGEOM,INC=20 *STATIC,DIRECT 1.,20. *BOUNDARY,OP=NEW FACE1,3 FACE3,2 FACE6,1 FACE4,1,1,3. *END STEP *STEP,NLGEOM,INC=20 *STATIC,DIRECT 1.,20. *BOUNDARY,OP=NEW FACE1,3 FACE3,2 FACE6,1 FACE4,1,1,0.0 *END STEP