|
[Sponsors] |
August 6, 2008, 11:50 |
Description:
I tried to dis
|
#1 |
Guest
Posts: n/a
|
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 |
|
August 7, 2008, 02:38 |
I think I get the floating poi
|
#2 |
Guest
Posts: n/a
|
I think I get the floating point exception, because the solution diverged on the tetgen-mesh.
|
|
|
|