*HEADING - INCOMPRESSIBLE HYPERELASTICITY (YEOH), PLANAR 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,R=1.25 *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. *BIAXIAL TEST DATA 483.638 , 2. 307.653 , 1.9 197.176 , 1.8 129.372 , 1.7 87.2224 , 1.6 60.1221 , 1.5 42.0717 , 1.4 29.5985 , 1.3 20.7535 , 1.2 14.3889 , 1.1 9.79593 , 1. 6.51097 , 900.E-03 4.21196 , 800.E-03 2.66104 , 700.E-03 1.67218 , 600.E-03 1.09216 , 500.E-03 788.463E-03, 400.E-03 640.653E-03, 300.E-03 533.528E-03, 200.E-03 354.571E-03, 100.E-03 0.,0. *BIAXIAL TEST DATA 2.52169E+03, 3. 1.58161E+03, 2.85 1.00121E+03, 2.7 650.501 , 2.55 435.203 , 2.4 298.236 , 2.25 207.459 , 2.1 145.069 , 1.95 100.914 , 1.8 69.3758 , 1.65 46.6858 , 1.5 30.4757 , 1.35 19.1023 , 1.2 11.363 , 1.05 6.33716 , 900.E-03 3.29594 , 750.E-03 1.64788 , 600.E-03 902.022E-03,450.E-03 635.266E-03,300.E-03 456.6E-03 ,150.E-03 0.,0. *PLANAR TEST DATA 98.9568 , 2. 62.4102 , 1.9 39.7517 , 1.8 25.958 , 1.7 17.4286 , 1.6 11.9816 , 1.5 8.37799 , 1.4 5.92302 , 1.3 4.21898 , 1.2 3.02851 , 1.1 2.2017 , 1. 1.63739 , 900.E-03 1.26259 , 800.E-03 1.02129 , 700.E-03 868.284E-03, 600.E-03 765.72E-03 , 500.E-03 681.104E-03, 400.E-03 586.018E-03, 300.E-03 455.248E-03, 200.E-03 266.316E-03, 100.E-03 0.,0. *PLANAR TEST DATA 561.053 , 3. 345.872 , 2.85 215.537 , 2.7 138.188 , 2.55 91.3716 , 2.4 61.8989 , 2.25 42.5327 , 2.1 29.3416 , 1.95 20.1378 , 1.8 13.6458 , 1.65 9.07634 , 1.5 5.90754 , 1.35 3.77131 , 1.2 2.3926 , 1.05 1.55638 , 900.E-03 1.08902 , 750.E-03 847.405E-03, 600.E-03 711.932E-03, 450.E-03 581.221E-03, 300.E-03 368.084E-03, 150.E-03 0.,0. *STEP,NLGEOM,INC=20 PLANAR TENSION *STATIC,DIRECT 1.,20. *BOUNDARY,OP=NEW FACE1,3 FACE3,2 FACE6,1 FACE4,1,1,2. FACE5,2,2,0. *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 FACE5,2,2,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. FACE5,2,2,0. *END STEP *STEP,NLGEOM,INC=20 *STATIC,DIRECT 1.,20. *BOUNDARY,OP=NEW FACE1,3 FACE3,2 FACE6,1 FACE4,1,1,0.0 FACE5,2,2,0.0 *END STEP