*heading
 Isotropic elasticity + Barlat plasticity
 Single-element uniaxial tensile tests with off-axis angle
 Check if the yield stress ratio (Sii/YIELDS) and r-value (PE22(global)/PE33(global)) agree with analytical solutions
*preprint, model=yes
*node
 111,  0.0, 0.0, 0.0
 112,  1.0, 0.0, 0.0
 113,  1.0, 1.0, 0.0
 114,  0.0, 1.0, 0.0
 115,  0.0, 0.0, 1.0
 116,  1.0, 0.0, 1.0
 117,  1.0, 1.0, 1.0
 118,  0.0, 1.0, 1.0
 171, 12.0, 0.0, 0.0
 172, 13.0, 0.0, 0.0
 173, 13.0, 1.0, 0.0
 174, 12.0, 1.0, 0.0
 175, 12.0, 0.0, 1.0
 176, 13.0, 0.0, 1.0
 177, 13.0, 1.0, 1.0
 178, 12.0, 1.0, 1.0
**
 311,  0.0, 8.0, 0.0
 312,  1.0, 8.0, 0.0
 313,  1.0, 9.0, 0.0
 314,  0.0, 9.0, 0.0
 315,  0.0, 8.0, 1.0
 316,  1.0, 8.0, 1.0
 317,  1.0, 9.0, 1.0
 318,  0.0, 9.0, 1.0
 371, 12.0, 8.0, 0.0
 372, 13.0, 8.0, 0.0
 373, 13.0, 9.0, 0.0
 374, 12.0, 9.0, 0.0
 375, 12.0, 8.0, 1.0
 376, 13.0, 8.0, 1.0
 377, 13.0, 9.0, 1.0
 378, 12.0, 9.0, 1.0
*ngen, nset=x0
 111, 171, 10
 112, 172, 10
 113, 173, 10
 114, 174, 10
*ngen, nset=x1
 115, 175, 10
 116, 176, 10
 117, 177, 10
 118, 178, 10
*ngen, nset=x3
 311, 371, 10
 312, 372, 10
 313, 373, 10
 314, 374, 10
*ngen, nset=x4
 315, 375, 10
 316, 376, 10
 317, 377, 10
 318, 378, 10
*nfill, nset=all
 x0, x3, 2, 100
 x1, x4, 2, 100
*nset, nset=left, generate
 111, 171, 10
 211, 271, 10
 311, 371, 10
 114, 174, 10
 214, 274, 10
 314, 374, 10
 115, 175, 10
 215, 275, 10
 315, 375, 10
 118, 178, 10
 218, 278, 10
 318, 378, 10
*nset, nset=fixdof2, generate
 111, 171, 10
 211, 271, 10
 311, 371, 10
 115, 175, 10
 215, 275, 10
 315, 375, 10
*nset, nset=fixdof3, generate
 111, 171, 10
 211, 271, 10
 311, 371, 10
*equation
 2,
 113 , 1, 1., 112 , 1, -1.
 2,
 116 , 1, 1., 112 , 1, -1.
 2,
 117 , 1, 1., 112 , 1, -1.
 2,
 123 , 1, 1., 122 , 1, -1.
 2,
 126 , 1, 1., 122 , 1, -1.
 2,
 127 , 1, 1., 122 , 1, -1.
 2,
 133 , 1, 1., 132 , 1, -1.
 2,
 136 , 1, 1., 132 , 1, -1.
 2,
 137 , 1, 1., 132 , 1, -1.
 2,
 143 , 1, 1., 142 , 1, -1.
 2,
 146 , 1, 1., 142 , 1, -1.
 2,
 147 , 1, 1., 142 , 1, -1.
 2,
 153 , 1, 1., 152 , 1, -1.
 2,
 156 , 1, 1., 152 , 1, -1.
 2,
 157 , 1, 1., 152 , 1, -1.
 2,
 163 , 1, 1., 162 , 1, -1.
 2,
 166 , 1, 1., 162 , 1, -1.
 2,
 167 , 1, 1., 162 , 1, -1.
 2,
 173 , 1, 1., 172 , 1, -1.
 2,
 176 , 1, 1., 172 , 1, -1.
 2,
 177 , 1, 1., 172 , 1, -1.
**
 2,
 213 , 1, 1., 212 , 1, -1.
 2,
 216 , 1, 1., 212 , 1, -1.
 2,
 217 , 1, 1., 212 , 1, -1.
 2,
 223 , 1, 1., 222 , 1, -1.
 2,
 226 , 1, 1., 222 , 1, -1.
 2,
 227 , 1, 1., 222 , 1, -1.
 2,
 233 , 1, 1., 232 , 1, -1.
 2,
 236 , 1, 1., 232 , 1, -1.
 2,
 237 , 1, 1., 232 , 1, -1.
 2,
 243 , 1, 1., 242 , 1, -1.
 2,
 246 , 1, 1., 242 , 1, -1.
 2,
 247 , 1, 1., 242 , 1, -1.
 2,
 253 , 1, 1., 252 , 1, -1.
 2,
 256 , 1, 1., 252 , 1, -1.
 2,
 257 , 1, 1., 252 , 1, -1.
 2,
 263 , 1, 1., 262 , 1, -1.
 2,
 266 , 1, 1., 262 , 1, -1.
 2,
 267 , 1, 1., 262 , 1, -1.
 2,
 273 , 1, 1., 272 , 1, -1.
 2,
 276 , 1, 1., 272 , 1, -1.
 2,
 277 , 1, 1., 272 , 1, -1.
**
 2,
 313 , 1, 1., 312 , 1, -1.
 2,
 316 , 1, 1., 312 , 1, -1.
 2,
 317 , 1, 1., 312 , 1, -1.
 2,
 323 , 1, 1., 322 , 1, -1.
 2,
 326 , 1, 1., 322 , 1, -1.
 2,
 327 , 1, 1., 322 , 1, -1.
 2,
 333 , 1, 1., 332 , 1, -1.
 2,
 336 , 1, 1., 332 , 1, -1.
 2,
 337 , 1, 1., 332 , 1, -1.
 2,
 343 , 1, 1., 342 , 1, -1.
 2,
 346 , 1, 1., 342 , 1, -1.
 2,
 347 , 1, 1., 342 , 1, -1.
 2,
 353 , 1, 1., 352 , 1, -1.
 2,
 356 , 1, 1., 352 , 1, -1.
 2,
 357 , 1, 1., 352 , 1, -1.
 2,
 363 , 1, 1., 362 , 1, -1.
 2,
 366 , 1, 1., 362 , 1, -1.
 2,
 367 , 1, 1., 362 , 1, -1.
 2,
 373 , 1, 1., 372 , 1, -1.
 2,
 376 , 1, 1., 372 , 1, -1.
 2,
 377 , 1, 1., 372 , 1, -1.
*nset, nset=right, generate
 112, 172, 10
 212, 272, 10
 312, 372, 10
**
*element, type=c3d8r
 11, 111, 112, 113, 114, 115, 116, 117, 118
*elgen, elset=all
 11, 7, 10, 1, 3, 100, 10
*elset, elset=xy_00
 11
*elset, elset=xy_15
 12
*elset, elset=xy_30
 13
*elset, elset=xy_45
 14
*elset, elset=xy_60
 15
*elset, elset=xy_75
 16
*elset, elset=xy_90
 17
*elset, elset=yz_00
 21           
*elset, elset=yz_15
 22           
*elset, elset=yz_30
 23           
*elset, elset=yz_45
 24           
*elset, elset=yz_60
 25           
*elset, elset=yz_75
 26           
*elset, elset=yz_90
 27
*elset, elset=zx_00
 31           
*elset, elset=zx_15
 32           
*elset, elset=zx_30
 33           
*elset, elset=zx_45
 34           
*elset, elset=zx_60
 35           
*elset, elset=zx_75
 36           
*elset, elset=zx_90
 37
** --------------------------------------------------------------------
*orientation, name=global
 1., 0., 0., 0., 1., 0.
 3, 0.
*orientation, name=ori_123_0
 1., 0., 0., 0., 1., 0.
 3, 0.
*orientation, name=ori_123_15
 1., 0., 0., 0., 1., 0.
 3, -15.
*orientation, name=ori_123_30
 1., 0., 0., 0., 1., 0.
 3, -30.
*orientation, name=ori_123_45
 1., 0., 0., 0., 1., 0.
 3, -45.
*orientation, name=ori_123_60
 1., 0., 0., 0., 1., 0.
 3, -60.
*orientation, name=ori_123_75
 1., 0., 0., 0., 1., 0.
 3, -75.
*orientation, name=ori_123_90
 1., 0., 0., 0., 1., 0.
 3, -90.
*orientation, name=ori_231_0
 0., 0., 1., 1., 0., 0.
 1, 0.
*orientation, name=ori_231_15
 0., 0., 1., 1., 0., 0.
 1, -15.
*orientation, name=ori_231_30
 0., 0., 1., 1., 0., 0.
 1, -30.
*orientation, name=ori_231_45
 0., 0., 1., 1., 0., 0.
 1, -45.
*orientation, name=ori_231_60
 0., 0., 1., 1., 0., 0.
 1, -60.
*orientation, name=ori_231_75
 0., 0., 1., 1., 0., 0.
 1, -75.
*orientation, name=ori_231_90
 0., 0., 1., 1., 0., 0.
 1, -90.
*orientation, name=ori_312_0
 0., 1., 0., 0., 0., 1.
 2, 0.
*orientation, name=ori_312_15
 0., 1., 0., 0., 0., 1.
 2, -15.
*orientation, name=ori_312_30
 0., 1., 0., 0., 0., 1.
 2, -30.
*orientation, name=ori_312_45
 0., 1., 0., 0., 0., 1.
 2, -45.
*orientation, name=ori_312_60
 0., 1., 0., 0., 0., 1.
 2, -60.
*orientation, name=ori_312_75
 0., 1., 0., 0., 0., 1.
 2, -75.
*orientation, name=ori_312_90
 0., 1., 0., 0., 0., 1.
 2, -90.
**
*solid section, elset=xy_00, material=alloy_2090t3, orientation=ori_123_0
*solid section, elset=xy_15, material=alloy_2090t3, orientation=ori_123_15
*solid section, elset=xy_30, material=alloy_2090t3, orientation=ori_123_30
*solid section, elset=xy_45, material=alloy_2090t3, orientation=ori_123_45
*solid section, elset=xy_60, material=alloy_2090t3, orientation=ori_123_60
*solid section, elset=xy_75, material=alloy_2090t3, orientation=ori_123_75
*solid section, elset=xy_90, material=alloy_2090t3, orientation=ori_123_90
*solid section, elset=yz_00, material=alloy_2090t3, orientation=ori_231_0
*solid section, elset=yz_15, material=alloy_2090t3, orientation=ori_231_15
*solid section, elset=yz_30, material=alloy_2090t3, orientation=ori_231_30
*solid section, elset=yz_45, material=alloy_2090t3, orientation=ori_231_45
*solid section, elset=yz_60, material=alloy_2090t3, orientation=ori_231_60
*solid section, elset=yz_75, material=alloy_2090t3, orientation=ori_231_75
*solid section, elset=yz_90, material=alloy_2090t3, orientation=ori_231_90
*solid section, elset=zx_00, material=alloy_2090t3, orientation=ori_312_0
*solid section, elset=zx_15, material=alloy_2090t3, orientation=ori_312_15
*solid section, elset=zx_30, material=alloy_2090t3, orientation=ori_312_30
*solid section, elset=zx_45, material=alloy_2090t3, orientation=ori_312_45
*solid section, elset=zx_60, material=alloy_2090t3, orientation=ori_312_60
*solid section, elset=zx_75, material=alloy_2090t3, orientation=ori_312_75
*solid section, elset=zx_90, material=alloy_2090t3, orientation=ori_312_90
**
*elset, elset=xy
 xy_00, xy_15, xy_30, xy_45, xy_60, xy_75, xy_90
*elset, elset=yz
 yz_00, yz_15, yz_30, yz_45, yz_60, yz_75, yz_90
*elset, elset=zx
 zx_00, zx_15, zx_30, zx_45, zx_60, zx_75, zx_90
** --------------------------------------------------------------------
*material, name=alloy_2090t3
*density
 2.7e-9
*elastic
 70e3, 0.33
*plastic
 150, 0.
 265, 0.247
*potential, type=barlat, power=8.
 -0.069888, 0.936408, 0.079143, 1.003060, 0.524741, 1.363180,  0.477161, 0.534530
  0.511885, 0.981171, 0.476741, 0.575316, 0.866827, 1.145010, -0.079294, 0.702310
  0.573550, 0.525830
** --------------------------------------------------------------------
*boundary
 left, 1, 1
 fixdof2, 2, 2
 fixdof3, 3, 3
** --------------------------------------------------------------------
*step, name=loading, nlgeom=yes
*dynamic, explicit
 , 1.
*variable mass scaling, type=below min, dt=1e-4, freq=10000
*boundary, type=velocity
 right, 1, 1,  5e-2
*output, field, number interval=100, time marks=yes
*node output
 U, RF
*element output
 LE, PE, S, PEEQ, YIELDPOT, YIELDS
*output, history, freq=1
*element output, elset=all
 PEEQ, YIELDPOT, YIELDS, PE
*node output, nset=right
 U1, RF1 
*energy output, variable=all
*end step