CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Solid Large deformation for icoFsiFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 12, 2008, 07:59
Default Hi everyone, Lovely day for a
  #1
Member
 
feng wang
Join Date: Mar 2009
Posts: 67
Rep Power: 8
fw407 is on a distinguished road
Hi everyone,
Lovely day for all of you!

firstly, I think it is necessary to explain what I am trying to do :I am trying to add the auto mesh motion to stressedFoam , and add some material models to stressedFoam so that a large deformation solver for solid mechanics can be created. when this is achieved, I want to merge this large deformation solver with icoFsiFoam to simulate the interaction of blood and blood vessel in an artery .


Unfortunately I run some problems at the starting point. the solver combing icoDyFoam and stressedFoam does not solve the mesh motion, and even though a "motionU" file is created for every time step, but no data is written into it and it remains the same in the "0" time step.


I am still a new foamer(just have runned icoFsiFoam for several weeks),and I have read some threads about mesh motion and icoDyFoam, but still have no idea how to solve this problem (maybe very easy for most people), could anyone help me?


Thanks in advance!


feng
fw407 is offline   Reply With Quote

Old   June 12, 2008, 19:12
Default Hi there, after a whole da
  #2
Member
 
feng wang
Join Date: Mar 2009
Posts: 67
Rep Power: 8
fw407 is on a distinguished road
Hi there,

after a whole day, I haven't succeed making the mesh move for my solver. Maybe something is wrong with it, the follows are the codes,if you find something stupid, please correct me, I will really appreciate your help!


#include "fvCFD.H"
#include "dynamicFvMesh.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

int main(int argc, char *argv[])
{

# include "setRootCase.H"

# include "createTime.H"
# include "createDynamicFvMesh.H"
# include "readMechanicalProperties.H"
# include "readThermalProperties.H"
# include "createFields.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Info<< "\nCalculating displacement field\n" << endl;

for (runTime++; !runTime.end(); runTime++)
{
Info<< "Iteration: " << runTime.timeName() << nl << endl;

# include "readStressedFoamControls.H"

int iCorr = 0;
scalar initialResidual = 0;

do
{
mesh.update();

volTensorField gradU = fvc::grad(U);

if (thermalStress)
{
solve
(
fvm::ddt(T) == fvm::laplacian(DT, T)
);
}

fvVectorMatrix UEqn
(
fvm::d2dt2(U)
==
fvm::laplacian(2*mu + lambda, U, "laplacian(DU,U)")

+ fvc::div
(
mu*gradU.T() + lambda*(I*tr(gradU)) - (mu + lambda)*gradU,
"div(sigma)"
)
);

if (thermalStress)
{
UEqn += threeKalpha*fvc::grad(T);
}

//UEqn.setComponentReference(1, 0, vector::X, 0);
//UEqn.setComponentReference(1, 0, vector::Z, 0);

initialResidual = UEqn.solve().initialResidual();

} while (initialResidual > convergenceTolerance && ++iCorr < nCorr);

# include "calculateStress.H"

Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
<< " ClockTime = " << runTime.elapsedClockTime() << " s"
<< nl << endl;
}
runTime.write();
deleteDemandDrivenData(Tptr);
Info<< "End\n" << endl;

return(0);
}

Kind regards

feng
fw407 is offline   Reply With Quote

Old   June 13, 2008, 10:58
Default Hi there, I solve this prob
  #3
Member
 
feng wang
Join Date: Mar 2009
Posts: 67
Rep Power: 8
fw407 is on a distinguished road
Hi there,

I solve this problem by using the method employed in icoFsiFoam to move the solid mesh(I should have seen it).

Now I want to model the deformation of the blood vessel which should be polymeric material. I find the development version of OF has added new solver for stress analysis like newStressedFoam and also something related to rheology and material models.If these stuff are the ones I need to model the blood vessel material(such as hyperelastic), Could some one give me some hints how to do that?

Thanks in advance!
fw407 is offline   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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Icofsifoam varun OpenFOAM Running, Solving & CFD 8 April 27, 2011 06:10
IcoFsiFoam vinz OpenFOAM Running, Solving & CFD 14 November 3, 2010 07:20
Solver for large deformation fw407 OpenFOAM Running, Solving & CFD 3 June 10, 2008 19:06
FSI with "solid deformation" siavashghassemi CFX 3 September 4, 2007 16:35
Large solid Ernie FLUENT 0 July 12, 2006 02:41


All times are GMT -4. The time now is 02:36.