CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Community Contributions > OpenFOAM CC Toolkits for Fluid-Structure Interaction

[solids4Foam] Solids Simulation: Stucked & FOAM Warning: From function eigenValues(const tensor&)

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 5, 2021, 06:52
Default Solids Simulation: Stucked & FOAM Warning: From function eigenValues(const tensor&)
  #1
New Member
 
Philipp Conen
Join Date: Jul 2021
Location: Germany, NRW
Posts: 22
Rep Power: 5
philippconen is on a distinguished road
Dear Foamers,

while doing a little benchmark with solids4Foam and ANSYS Mechanical I run some simulations of a Cantilever. The mesh is made in ANSYS and imported via fluentMeshToFoam.

Dealing with 150.000 cells made the model in OpenFOAM stuck. On my local computer I saw an increasing amount of used RAM until nothing more was available, so my computer crashed.
Also on the super computer of my university the simulation seems to stucks (without knowledge about the used resources - ?/64GB RAM).

My question would be, if the stucked simulation is related to the warning OpenFOAM gave me? What could be the reason for the stuck? How can I fix it?

Here are my pre logs:
Code:
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
| solids4Foam::convertCaseFormat start                                |
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

foam-extend-4.0 loaded: no changes made

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | foam-extend: Open Source CFD                    |
|  \\    /   O peration     | Version:     4.1-20200408                       |
|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
|    \\/     M anipulation  | For copyright notice see file Copyright         |
\*---------------------------------------------------------------------------*/
Build  : 4.1-20200408
Exec   : fluentMeshToFoam mech0_005.msh -writeSets
Date   : Oct 05 2021
Time   : 10:54:04
Host   : "cc-v01"
PID    : 12758
CtrlDict : "/home/ms/p/pc196535/Schreibtisch/MasterthesisDLR/Simulationen/solids4Foam/Sensitivitaetsanalyse/Balken_Sensitivitaetsanalyse/Balken_tipLoad/Balken_tipLoad_0_005/system/controlDict"
Case   : /home/ms/p/pc196535/Schreibtisch/MasterthesisDLR/Simulationen/solids4Foam/Sensitivitaetsanalyse/Balken_Sensitivitaetsanalyse/Balken_tipLoad/Balken_tipLoad_0_005
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Dimension of grid: 3
Number of points: 168861

number of faces: 468500
Number of cells: 150000
Reading points
Reading points
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Other readCellGroupData: 2 1 249f0 1 4
Reading uniform cells
Read zone1:1 name:interior-solid patchTypeID:interior
Reading zone data
Read zone1:2 name:solid patchTypeID:solid
Reading zone data
Read zone1:5 name:support patchTypeID:wall
Reading zone data
Read zone1:6 name:tip patchTypeID:wall
Reading zone data
Read zone1:7 name:top patchTypeID:wall
Reading zone data
Read zone1:8 name:bot patchTypeID:wall
Reading zone data
Read zone1:9 name:right patchTypeID:wall
Reading zone data
Read zone1:10 name:left patchTypeID:wall
Reading zone data


FINISHED LEXING


dimension of grid: 3
Creating shapes for 3-D cells
Building patch-less mesh...--> FOAM Warning : 
    From function polyMesh::polyMesh(... construct from shapes...)
    in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 627
    Found 37000 undefined faces in mesh; adding to default patch.
done.

Building boundary and internal patches.
maxZoneID: 10
Creating patch 0 for zone: 1 start: 1 end: 431500 type: interior name: interior-solid
Creating patch 1 for zone: 5 start: 431501 end: 432000 type: wall name: support
Creating patch 2 for zone: 6 start: 432001 end: 432500 type: wall name: tip
Creating patch 3 for zone: 7 start: 432501 end: 447500 type: wall name: top
Creating patch 4 for zone: 8 start: 447501 end: 462500 type: wall name: bot
Creating patch 5 for zone: 9 start: 462501 end: 465500 type: wall name: right
Creating patch 6 for zone: 10 start: 465501 end: 468500 type: wall name: left
Writing patch interior-solid of size 431500 to faceSet.
Patch interior-solid is internal to the mesh  and is not being added to the boundary.
Writing patch support of size 500 to faceSet.
Adding new patch support of type wall as patch 0
Writing patch tip of size 500 to faceSet.
Adding new patch tip of type wall as patch 1
Writing patch top of size 15000 to faceSet.
Adding new patch top of type wall as patch 2
Writing patch bot of size 15000 to faceSet.
Adding new patch bot of type wall as patch 3
Writing patch right of size 3000 to faceSet.
Adding new patch right of type wall as patch 4
Writing patch left of size 3000 to faceSet.
Adding new patch left of type wall as patch 5

Default patch type set to empty

Writing mesh... to "constant/polyMesh"  done.

Only one cell group: no set written

End

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | foam-extend: Open Source CFD                    |
|  \\    /   O peration     | Version:     4.1-20200408                       |
|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
|    \\/     M anipulation  | For copyright notice see file Copyright         |
\*---------------------------------------------------------------------------*/
Build  : 4.1-20200408
Exec   : setsToZones -noFlipMap
Date   : Oct 05 2021
Time   : 10:54:07
Host   : "cc-v01"
PID    : 12759
CtrlDict : "/home/ms/p/pc196535/Schreibtisch/MasterthesisDLR/Simulationen/solids4Foam/Sensitivitaetsanalyse/Balken_Sensitivitaetsanalyse/Balken_tipLoad/Balken_tipLoad_0_005/system/controlDict"
Case   : /home/ms/p/pc196535/Schreibtisch/MasterthesisDLR/Simulationen/solids4Foam/Sensitivitaetsanalyse/Balken_Sensitivitaetsanalyse/Balken_tipLoad/Balken_tipLoad_0_005
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Searched : "constant/polyMesh/sets"
Found    : 
7
(
interior-solid
tip
bot
right
top
left
support
)


Adding set right as a faceZone.
Adding set top as a faceZone.
Adding set left as a faceZone.
Adding set interior-solid as a faceZone.
Adding set tip as a faceZone.
Adding set bot as a faceZone.
Adding set support as a faceZone.
Writing mesh.

End
Here are my solver logs:
Code:
// ************************************************************************* //
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | foam-extend: Open Source CFD                    |
|  \\    /   O peration     | Version:     4.1-20200408                       |
|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
|    \\/     M anipulation  | For copyright notice see file Copyright         |
\*---------------------------------------------------------------------------*/
Build  : 4.1-20200408
Exec   : solids4Foam
Date   : Oct 05 2021
Time   : 10:54:08
Host   : "cc-v01"
PID    : 12762
CtrlDict : "/home/ms/p/pc196535/Schreibtisch/MasterthesisDLR/Simulationen/solids4Foam/Sensitivitaetsanalyse/Balken_Sensitivitaetsanalyse/Balken_tipLoad/Balken_tipLoad_0_005/system/controlDict"
Case   : /home/ms/p/pc196535/Schreibtisch/MasterthesisDLR/Simulationen/solids4Foam/Sensitivitaetsanalyse/Balken_Sensitivitaetsanalyse/Balken_tipLoad/Balken_tipLoad_0_005
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

/*---------------------------------------------------------------------------*\
|    For further information on the solids4Foam toolbox implementations,      |
|    please see the following publications:                                   |
|                                                                             |
|    P. Cardiff, A Karac, P. De Jaeger, H. Jasak, J. Nagy, A. Ivankovic,      |
|    Z. Tukovic: An open-source finite volume toolbox for solid mechanics and |
|    fluid-solid interaction simulations. arXiv:1808.10736v2, 2018, available |
|    at https://arxiv.org/abs/1808.10736.                                     |
|                                                                             |
|    Z. Tukovic, A. Karac, P. Cardiff, H. Jasak, A. Ivankovic: OpenFOAM       |
|    finite volume solver for fluid-solid interaction.  Transactions of       |
|    Famena, 42 (3), pp. 1-31, 2018, 10.21278/TOF.42301.                      |
\*---------------------------------------------------------------------------*/

Selecting physicsModel solid

Selecting solidModel coupledUnsLinearGeometryLinearElastic
Selecting dynamicFvMesh staticFvMesh
Patch tip    Traction boundary field 
Patch top    Traction boundary field 
Patch bot    Traction boundary field 
Patch right    Traction boundary field 
Patch left    Traction boundary field 
    under-relaxation method: fixed

/*---------------------------------------------------------------------------*\
|    For further information on this class, please see the following          |
|    publication:                                                             |
|                                                                             |
|    P. Cardiff, Z. Tukovic, H, Jasak, A. Ivankovic: A block-coupled finite   |
|    volume methodology for linear elasticity and unstructured meshes.        |
|    Computers and Structures, 2016, 175 100-122,                             |
|    DOI: 10.1016/-j.compstruc.2016.07.004.                                   |
\*---------------------------------------------------------------------------*/

Creating the mechanicalModel
Selecting mechanical law linearElastic
Time = 1

Evolving solid solver
Solving the momentum equation for D
--> FOAM Warning : 
    From function eigenValues(const tensor&)
    in file primitives/Tensor/tensor/tensor.C at line 180
    complex eigenvalues detected for tensor: (4.6875e-06 -1.5625e-06 1.5625e-06 -1.5625e-06 4.6875e-06 1.5625e-06 1.5625e-06 1.5625e-06 4.6875e-06)
--> FOAM Warning : 
    From function eigenValues(const tensor&)
    in file primitives/Tensor/tensor/tensor.C at line 180
    complex eigenvalues detected for tensor: (4.6875e-06 -1.5625e-06 1.5625e-06 -1.5625e-06 4.6875e-06 1.5625e-06 1.5625e-06 1.5625e-06 4.6875e-06)
--> FOAM Warning : 
    From function eigenValues(const tensor&)
    in file primitives/Tensor/tensor/tensor.C at line 180
    complex eigenvalues detected for tensor: (4.6875e-06 -1.5625e-06 -1.5625e-06 -1.5625e-06 4.6875e-06 -1.5625e-06 -1.5625e-06 -1.5625e-06 4.6875e-06)
--> FOAM Warning : 
    From function eigenValues(const tensor&)
    in file primitives/Tensor/tensor/tensor.C at line 180
    complex eigenvalues detected for tensor: (4.6875e-06 -1.5625e-06 -1.5625e-06 -1.5625e-06 4.6875e-06 -1.5625e-06 -1.5625e-06 -1.5625e-06 4.6875e-06)
--> FOAM Warning : 
    From function eigenValues(const tensor&)
    in file primitives/Tensor/tensor/tensor.C at line 180
    complex eigenvalues detected for tensor: (4.6875e-06 -1.5625e-06 -1.5625e-06 -1.5625e-06 4.6875e-06 -1.5625e-06 -1.5625e-06 -1.5625e-06 4.6875e-06)
--> FOAM Warning : 
    From function eigenValues(const tensor&)
    in file primitives/Tensor/tensor/tensor.C at line 180
    complex eigenvalues detected for tensor: (4.6875e-06 -1.5625e-06 -1.5625e-06 -1.5625e-06 4.6875e-06 -1.5625e-06 -1.5625e-06 -1.5625e-06 4.6875e-06)
EigenSolver : writeMatlabFiles is false
EigenSparseLU Eigen linear solver
The results of the benchmark I can access until now are available in the attachments.

Greetings!
Attached Images
File Type: png cantileverTime.png (9.4 KB, 7 views)
File Type: png cantileverDisp.png (13.7 KB, 4 views)
philippconen is offline   Reply With Quote

Reply

Tags
aborting, ansys, cantiliver, solids4foam, warning

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[mesh manipulation] refineMesh Error mohsen.boojari OpenFOAM Meshing & Mesh Conversion 3 March 1, 2018 23:07
is internalField(U) equivalent to zeroGradient? immortality OpenFOAM Running, Solving & CFD 7 March 29, 2013 02:27
channelFoam for a 3D pipe AlmostSurelyRob OpenFOAM 3 June 24, 2011 14:06
latest OpenFOAM-1.6.x from git failed to compile phsieh2005 OpenFOAM Bugs 25 February 9, 2010 05:37
Warning 097- AB Siemens 6 November 15, 2004 05:41


All times are GMT -4. The time now is 05:59.