[Jerome Berryhill]

I made a modification of the beamEndLoad case (bEL). The original code has a 30X1X1 beam, with one end fixed, and a shear traction applied in the Y direction applied to the opposite end. The beam is represented as a single block with (6,2,1) cell structure. I modified this, so that the last sixth of the beam is a separate block, and put a face on the -Y side of that block, with a traction in the +Y direction. That is, I replaced the shear on the end with a normal force from the bottom, in the same direction. I would expect the results to be nearly identical, but in fact they are very different.

Is this a bug? Or am I doing something wrong?
from blockMesh;
vertices
(

(0 0 0) // 0
(25 0 0) // 1
(25 1 0) // 2
(0 1 0) // 3
(0 0 1) // 4
(25 0 1) // 5
(25 1 1) // 6
(0 1 1) // 7
(30 0 0) // 8
(30 1 0) // 9
(30 0 1) // 10
(30 1 1) // 11

);

blocks
(
hex (0 1 2 3 4 5 6 7) (5 2 1) simpleGrading (1.0 1.0 1.0) // left part of beam
hex (1 8 9 2 5 10 11 6) (1 2 1) simpleGrading (1.0 1.0 1.0)
);

edges
(
);

boundary
(
pressureFace
{
type patch;
faces
(
(1 2 9 8)
);
}

fixedEnd
{
type patch;
faces
(
(0 4 7 3)
);
}

);

from 0/D

dimensions [ 0 1 0 0 0 0 0 ];

internalField uniform ( 0 0 0 );

boundaryField
{

pressureFace
{
type tractionDisplacement;
traction uniform (0 10000 0);
pressure uniform 0;
value uniform (0 0 0);
}

fixedEnd
{
type fixedValue;
value uniform (0 0 0);
}


defaultFaces
{
type empty;
}
}

[Jerome Berryhill]

Ooops. Never mind. I had x1, x2, x3 confused with x, y, z.