The nsteps and the computing time of Peri3dController

Asked by liukeqi

Hello,
    I want to couple DEM and FEM in order to implement the multiscale analysis just like the method in Ref[1]. I choose the engine Peri3dController to do it, because this engine has many advantages, such as it can use strain to get stress directly and it can save the former strain state by using p3d.strain=strain. In the FEM, I use the central difference method to compute, and it has stability condition that the loading step is very small. The result is that the compute times are very big and each strain increment is very small.
    In the DEM, I called the engine Peri3dController each compute time and I found if I use the constant nSteps=4000, it will be a very big computation cost. So I want to find an appropriate nSteps that can guarantee the result is relative steady and meanwhile can be small to reduce the computation cost. I also find if I called the Peri3dController too many times, the compute time will increase along with each calling. I do not whether it is related to that I use saveTmp() and loadTmp().
    I accord to the file example/test/peri3dController_triaxialCompression.py, the file use two engine, can I use one engine and only change the goal of the Peri3dController to get the same result?
    Can you give me some suggestions? Thank you.

Keqi

[1] http://www.engmech.cz/2012/proceedings/pdf/018_Stransky_J-FT.pdf

Question information

Language:
English Edit question
Status:
Solved
For:
Yade Edit question
Assignee:
No assignee Edit question
Solved by:
Jan Stránský
Solved:
Last query:
Last reply:

This question was reopened

Revision history for this message
Bruno Chareyre (bruno-chareyre) said :
#1

Hi,
I don't quite understand the question overall, I'm sorry, but did you know a multscale DEMxFEM coupling was already available?
Check https://yade-dem.org/doc/FEMxDEM.html
Bruno

Revision history for this message
liukeqi (891979456-h) said :
#3

Thanks Bruno. I knew the project of https://yade-dem.org/doc/FEMxDEM.html. But I want to couple YADE with another FEM software. I mean that the memory occupation is more and more along with the times I called Peri3dController continuously in one *.py. Are there some method to release the memory occupation and guarantee the function strain=p3d.strain of Peri3dController available?

Revision history for this message
Best Jan Stránský (honzik) said :
#4

Hello,

have a look at Bruno's FEMxDEM link.

So I want to find an appropriate nSteps that can guarantee the result is
> relative steady

It strongly depends on the definition of "relative steady", material and
other factors. I think trial-and-error approach is the best one to find
such nSteps value

cheers
Jan

2016-08-01 6:08 GMT+02:00 liukeqi <email address hidden>:

> Question #313433 on Yade changed:
> https://answers.launchpad.net/yade/+question/313433
>
> Status: Solved => Open
>
> liukeqi is still having a problem:
> Thanks Bruno. I knew the project of https://yade-
> dem.org/doc/FEMxDEM.html. But I want to couple YADE with another FEM
> software. I mean that the memory occupation is more and more along with
> the times I called Peri3dController continuously in one *.py. Are there
> some method to release the memory occupation and guarantee the function
> strain=p3d.strain of Peri3dController available?
>
> --
> You received this question notification because your team yade-users is
> an answer contact for Yade.
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-users
> Post to : <email address hidden>
> Unsubscribe : https://launchpad.net/~yade-users
> More help : https://help.launchpad.net/ListHelp
>

Revision history for this message
liukeqi (891979456-h) said :
#5

Thanks Jan Stránský, that solved my question.