Description:

I tried to discretize the beam-domain of the beamEndLoad-tutorial in 3D instead of blocks (about the blockMesh-utility) with tetraeders (about tetgen) and the solidEquilibriumDisplacementFoam-Solver doesn't work anymore.

Details:
I took the beamEndLoad-tutorial, write the mesh of the beam as 3 dimensional and refine it:

convertToMeters 1;

vertices
(

(0 0 0)
(30 0 0)
(30 1 0)
(0 1 0)
(0 0 1)
(30 0 1)
(30 1 1)
(0 1 1)

);

blocks
(
hex (0 1 2 3 4 5 6 7) (60 2 2) simpleGrading (1.0 1.0 1.0)
);

edges
(
);

patches
(
patch
topSurface
(
(3 7 6 2)
)

patch
bottomSurface
(
(4 0 1 5)
)

patch
fixedEnd
(
(0 4 7 3)
)

patch
tractionEnd
(
(1 2 6 5)
)
patch
defaultFaces
(
(0 3 2 1)
(5 6 7 4)
)
);

(Also I switch the planeStress-option in the mechanicalProperties-file)
If I applicate the solidEquilibriumDisplacementFoam I get a solution.

Now I take the same shape of the beam, but create a mesh with tetgen. For that I generate a tetgen-mesh with
tetgen -qfa0.125 beam.poly
where beam.poly is the following file:

# Part 1 - node list
# node count, 3 dim, no attribute, no boundary marker
8 3 0 0
# Node index, node coordinates
1 0.0 0.0 0.0
2 30 0.0 0.0
3 30 1.0 0.0
4 0.0 1.0 0.0
5 0.0 0.0 1.0
6 30 0.0 1.0
7 30 1.0 1.0
8 0.0 1.0 1.0

# Part 2 - facet list
# facet count, boundary marker
6 1
# facets
1 0 1 # 1 polygon, no hole, boundary marker
4 1 2 3 4 # front
1 0 1
4 5 6 7 8 # back
1 0 2
4 1 2 6 5 # bottom
1 0 3
4 2 3 7 6 # right
1 0 4
4 3 4 8 7 # top
1 0 5
4 4 1 5 8 # left

# Part 3 - hole list
0 # no hole

# Part 4 - region list
0 # no region

The obtained mesh (beam.1.node, beam.1.ele, beam.1.face) I converted with tetgenToFoam and applicate the SolidEquilibriumDisplacementFoam-Solver.
Now I get a floating point exception.

So I think I did nothing else than discretize the mesh instead of blocks with tetraeders (about tetgen) and the solidEquilibriumDisplacementFoam-solver didn't work anymore.

Is this a bug?


Application: solidEquilibriumDisplacementFoam

Source: $FOAM_SOLVERS/stressAnalysis/solidEquilibriumDisplacementFoam

Platform: Linux, Version 2.6.18, i686 GNU/Linux

Version: OpenFOAM 1.5

I think I get the floating point exception, because the solution diverged on the tetgen-mesh.