*HEADING HYPERELASTIC MARLOW TRELOAR'S EXPERIMENTAL DATA BIAXIAL TEST DATA ONLY *RESTART,WRITE,FREQUENCY=5 *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,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 ** *NSET,NSET=FACE22 5,6,7 *EQUATION 2, FACE22,3,1, 8,3,-1 ** *ELEMENT,TYPE=C3D8RH,ELSET=ONE 1,1,2,3,4,5,6,7,8 *SOLID SECTION,ELSET=ONE,MATERIAL=TREL *MATERIAL,NAME=TREL *HYPERELASTIC,MARLOW *BIAXIAL TEST DATA 0.0, 0.0 0.9384, 0.0200 1.5900, 0.0600 2.4087, 0.1100 2.6220, 0.1400 3.3240, 0.2000 4.4278, 0.3100 5.1830, 0.4200 6.6024, 0.6800 7.7794, 0.9400 9.7857, 1.4900 12.6351, 2.0300 14.6804, 2.4300 17.4000, 2.7500 20.1058, 3.0700 22.4502, 3.2600 24.6530, 3.4500 *STEP,NLGEOM,INC=50 Step 1: Uniaxial Tension *STATIC,DIRECT .25,7. *BOUNDARY,OP=NEW FACE1,3 FACE3,2 FACE6,1 FACE4,1,1,7. *ENERGY PRINT *EL PRINT,FREQUENCY=5 S, E, *NODE PRINT,FREQUENCY=5 U,RF *NODE FILE,FREQUENCY=2 U,RF *OUTPUT,FIELD,FREQUENCY=0 *NODE OUTPUT,NSET=ALL U,RF *OUTPUT,HISTORY,FREQUENCY=1 *NODE OUTPUT,NSET=ALL U,RF *END STEP *STEP,NLGEOM,INC=20 Step 2: Unload *STATIC,DIRECT .5,7. *BOUNDARY,OP=MOD FACE4,1 *END STEP *STEP,NLGEOM,INC=50 Step 3: Planar Tension (Pure Shear) *STATIC,DIRECT .25,4. *BOUNDARY,OP=NEW FACE1,3 FACE3,2 FACE5,2 FACE6,1 FACE4,1,1,4. *END STEP *STEP,NLGEOM,INC=20 Step 4: Unload *STATIC,DIRECT .5,4. *BOUNDARY,OP=MOD FACE4,1 *END STEP *STEP,NLGEOM,INC=50 Step 5: Biaxial Tension *STATIC .1,4.,,0.25 *BOUNDARY,OP=NEW FACE1,3 FACE3,2 FACE6,1 FACE4,1,1,4. FACE5,2,2,4. *END STEP *STEP,NLGEOM,INC=20 Step 6: Unload *STATIC .2,4.,,0.5 *BOUNDARY,OP=MOD FACE4,1 FACE5,2 *END STEP