CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Bugs

SolidEqulibriumDisplacementFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   August 6, 2008, 11:50
Default Description: I tried to dis
  #1
karen
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
  Reply With Quote

Old   August 7, 2008, 02:38
Default I think I get the floating poi
  #2
karen
Guest
 
Posts: n/a
I think I get the floating point exception, because the solution diverged on the tetgen-mesh.
  Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On



All times are GMT -4. The time now is 03:33.