CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Modifying stressedFoam for geometric non-linearity

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

LinkBack Thread Tools Search this Thread Display Modes
Old   May 25, 2020, 12:20
Smile Modifying stressedFoam for geometric non-linearity
New Member
Craig Thomas
Join Date: Dec 2019
Location: Leeds, UK
Posts: 12
Rep Power: 5
mn14cat is on a distinguished road
Hello fellow foamers,

My name is Craig and this is my first time posting to the forum, so I hope that I am following the Forum rules correctly.

To provide background and context to my post, I am currently working on a research project which I am using foam-extend v4.0 to model cartilage and conduct a stress analysis. I am developing both my model and solver in an incremental fashion, starting with infinitesimal strain moving to finite strain.

I am trying to modify the stressedFoam solver to account for finite strains. I am aware that foam-extend already offers solvers which can handle this task, however I am required to develop my own solver and therefore have started with a solver I feel comfortable using and can take the steps of modifying.

My model is a simple column shape; 250 mm in height and with a 400 mm^2 cross sectional area.

I have fixed one end (the bottom) and applied a load to other (the top). the walls are free to deform in a uniaxial direction.

I have used the traction boundary condition on the walls and top of the column and everything works as expected for the linear-elastic case.

When I try to now implement non linearity to the model I am receiving an error which I do not understand which I have included below. I have also attached both the .C of the solver and traction displacement boundary for those who wish to see how I have implemented the nonlinearity.

| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: |
| \\/ M anipulation | For copyright notice see file Copyright |
Build : 4.0-268bb07d15d8
Exec : NLStressedFoam
Date : May 25 2020
Time : 15:41:49
Host : mn14cat-XPS-13-7390
PID : 17826
CtrlDict : "/home/mn14cat/foam/mn14cat-4.0/run/exercise/NonLinStressedFoam/NonlinColumn/system/controlDict"
Case : /home/mn14cat/foam/mn14cat-4.0/run/exercise/NonLinStressedFoam/NonlinColumn
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

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

Create mesh for time = 0

Reading mechanical properties

Normalising E : E/rho

Calculating Lame's coefficients

Plane Strain

mu = 500000 Pa/rho
lambda = 4.66668e+06 Pa/rho
threeK = 1.5e+07 Pa/rho
Reading thermal properties

Reading field U

Calculating displacement field

Iteration: 1


request for volTensorField grad(U).T() & grad(U) from objectRegistry region0 failed
available objects of type volTensorField are


From function objectRegistry::lookupObject<Type>(const word&) const
in file /home/mn14cat/foam/foam-extend-4.0/src/foam/lnInclude/objectRegistryTemplates.C at line 139.

FOAM aborting

Aborted (core dumped)

I would be extremely grateful for any help understanding where I have gone wrong in modifying the solver, and hope that I have provided enough detail for anyone who is willing to help.

Thanks in advance,
Attached Files
File Type: c NLTractionDisplacementFvPatchVectorField.C (6.4 KB, 1 views)
File Type: c NLStressedFoam.C (3.8 KB, 1 views)
mn14cat is offline   Reply With Quote


foam-extend 4.0, modify solver, non-linear problems

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
Linearity and Non linearity sampathevs FLUENT 1 February 6, 2014 08:17
stressedFoam with wall motion looses symmetry and does not converge gtg258f OpenFOAM 1 June 6, 2012 13:28

All times are GMT -4. The time now is 16:26.