*HEADING HYPERELASTIC TEST DATA INPUT TRELOAR'S EXPERIMENTAL DATA *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 *ELEMENT,TYPE=C3D8RH,ELSET=ONE 1,1,2,3,4,5,6,7,8 *SOLID SECTION,ELSET=ONE,MATERIAL=TREL *MATERIAL,NAME=TREL *HYPERELASTIC,MOONEY-RIVLIN,TEST DATA INPUT *UNIAXIAL TEST DATA 1.5506, 0.1338 2.4367, 0.2675 3.1013, 0.3567 4.2089, 0.6242 5.3165, 0.8917 5.9810, 1.1592 6.8671, 1.4268 8.8608, 2.0510 10.6329, 2.5860 12.4051, 3.0318 16.1709, 3.7898 19.9367, 4.3694 23.4810, 4.8153 27.4684, 5.1720 31.0127, 5.4395 34.5570, 5.7070 38.3228, 5.9299 42.0886, 6.0637 45.6329, 6.1975 49.3987, 6.3312 53.1646, 6.4650 56.9304, 6.5541 64.2405, 6.6433 ** To fit uniaxial test data only, remove the *BIAXIAL and *PLANAR TEST ** DATA cards and the stress-strain data below. *BIAXIAL TEST DATA 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 *PLANAR TEST DATA ** *PLANAR TEST DATA is the new name for *SHEAR TEST DATA (Pure Shear). 0.6000, 0.0690 1.6000, 0.1034 2.4000, 0.1724 3.3600, 0.2828 4.2000, 0.4276 6.0000, 0.8483 7.8000, 1.3862 9.6000, 2.0000 11.1200, 2.4897 12.9600, 3.0345 14.8800, 3.4483 16.5800, 3.7793 18.2000, 4.0621 *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: Biaxial Tension *STATIC,DIRECT .25,4. *BOUNDARY,OP=NEW FACE1,3 FACE3,2 FACE6,1 FACE4,1,1,4. FACE5,2,2,4. *END STEP *STEP,NLGEOM,INC=20 Step 4: Unload *STATIC,DIRECT .5,4. *BOUNDARY,OP=MOD FACE4,1 FACE5,2 *END STEP *STEP,NLGEOM,INC=50 Step 5: 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