← Back to team overview

yade-users team mailing list archive

Re: [Question #229268]: Modelling a free-fall test of a concrete block

 

Question #229268 on Yade changed:
https://answers.launchpad.net/yade/+question/229268

    Status: Open => Answered

Jan Stránský proposed the following answer:
Hello,

welcome to Yade comunity :-)


> My first question is, basically, is Yade capable of simulate my main
> purpose? or, should I modify the models or code so that Yade works this
> test? I'm using Yade-daily version which, currently is 3+3602+47~precise1
>

In general, Yade is capable to simulate such scenario. However, the
concrete behavior is very complicated and Cpm model might not by accurate
in all possible situations, but it is under permanent development :-) see
[1] for more information.


idConcrete=O.materials.append(CpmMat(young=24e9,frictionAngle=atan(0.8),poisson=.2,density=4800,sigmaT=3.5e6,epsCrackOnset=1e-4,isoPrestress=0))
> # Concrete block material definition
>
> idGround=O.materials.append(CpmMat(young=30e30,density=1000,poisson=.2,frictionAngle=.5))
> # Concrete ground material definition
>
>
CpmMat has several parameters, which are not optional and have to be
defined. They are:
- epsCrackOnset
- relDuctility
so just add them as CpmMat arguments and the script will run OK.

If you have any other question, don't hesitate to ask
cheers
Jan


PS: some comments on the script (not necessary to read them imediatelly,
but before you proceed you should read them):
- it is always good to set manually the time step (e.g. O.dt =
0.5*utils.PWaveTimeStep() ) or to use GlobalTimeStepper, [2,3]
- Law2_Dem3DofGeom.. does not exist eny more, so you can replace body
InteractionLoop [4]:

  InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=1.5,label='bo1s'),Bo1_Wall_Aabb()])
        InteractionLoop(

[Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=1.5,label='ig2ss'),Ig2_Wall_Sphere_ScGeom()],
                [Ip2_CpmMat_CpmMat_CpmPhys()],
                [Law2_Dem3DofGeom_CpmPhys_Cpm(epsSoft=0) if not ScGeom else
Law2_ScGeom_CpmPhys_Cpm()],
        ),
- it is usual to run at least one step inside the script. In the case of
cohesive models, furthermore you should reset the enlargeFactors set
to Bo1_Sphere_Aabb and Ig2_Sphere_Sphere_Dem3DofGeom [4,5]. i.e. at the end
of the script:
  O.step()
  bo1s.aabbEnlargeFactor = ig2ss.interactionDetectionFactor = 1.0


[1] http://beta.arcig.cz/~eudoxos/smilauer2010-phd-thesis.pdf
[2] https://yade-dem.org/doc/yade.utils.html#yade._utils.PWaveTimeStep
[3]
https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.GlobalStiffnessTimeStepper
[4]
http://bazaar.launchpad.net/~yade-pkg/yade/git-trunk/view/head:/examples/concrete/uniax.py
[5] https://yade-dem.org/doc/user.html#creating-interactions (see the end
of example script)

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.