Question #403169 on Yade changed:
https://answers.launchpad.net/yade/+question/403169
Status: Open => Answered
Chareyre proposed the following answer:
Hello. This question is to vague hence unlikely to receive a response. What
is the meaning of "dramatically altered"? Is it certain that the script
generates a sensible mechanical system for every size? Etc
Le 19 oct. 2016 06:08, "liukeqi" <question403...@answers.launchpad.net> a
écrit :
> New question #403169 on Yade:
> https://answers.launchpad.net/yade/+question/403169
>
> Hi, everyone
> I want to make a simple example to simulate the behaviour of concrete
> under pressure load. I used the periodic boundary and I suffered a
> question. I found that the stress-strain relationship will dramatically
> altered as I changed the size of the periodic boundary(that is the number
> of the bodies), but I need the relative stability result to convince
> others. I want to know how many number of bodies will lead to a steady
> result that will not dramatically alter as the number of bodies continue
> increasing. The following is my code, and you can alter the initSize to
> change the size of the periodic boundary. Can you give me some suggestions?
> Thank you.
>
> #################################################
> #!/usr/bin/python
> import string
> from yade import plot,qt
> from yade.pack import *
> from yade import pack, plot
>
> O.materials.append(CpmMat(young=24e9,density=4800,frictionAngle=atan(0.8),
> poisson=0.2,sigmaT=600e6,epsCrackOnset=1e-4,relDuctility=35))
>
> initSize=0.025
>
> sp=pack.randomPeriPack(radius=.000625,initSize=Vector3(
> initSize,initSize,initSize))
> O.periodic=True
> sp.toSimulation()
>
>
> O.dt=1e-7
>
>
> plot.plots={'ey':('sy',)}
> def plotAddData():
> plot.addData(
> sy=p3d.stress[1],
> ey=p3d.strain[1],
> ex=p3d.strain[0]
> )
>
> EnlargeFactor=1.5
> EnlargeFactor=1.0
> O.engines=[
> ForceResetter(),
> InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=
> EnlargeFactor,label='bo1s')]),
> InteractionLoop(
> [Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=
> EnlargeFactor,label='ig2ss')],
> [Ip2_CpmMat_CpmMat_CpmPhys()],
> [Law2_ScGeom_CpmPhys_Cpm()]),
> NewtonIntegrator(),
> Peri3dController(
>
> nSteps=4000,
> # how many time steps the simulation will last
> # after reaching
> nSteps do doneHook action
> doneHook='print
> "Simulation with Peri3dController finished."; O.pause()',
>
> # the prescribed
> path (step,value of stress/strain) can be defined in absolute values
>
> # or in relative
> values
>
> # if the goal
> value is 0, the absolute stress/strain values are always considered (step
> values remain relative)plot.plots={'ex':('sx',)}
>
> # if ##Path is not
> explicitly defined, it is considered as linear function between (0,0) and
> (nSteps,goal)
> # as in yzPath and
> xyPath
> # the relative
> values are really relative (zxPath gives the same - except of the sign from
> goal value - result as yyPath)
>
>
> # variables used
> in the first step
> label='p3d'
> ),
> PyRunner(command='plotAddData()',iterPeriod=1),
> ]
>
>
>
> p3d.stressMask=0b100001 # prescribed ex,ey,sz,syz,ezx,sxy;
> e..strain; s..stress
> p3d.goal=(0,-0.008,0,0,0,0) #xx, yy, zz, yz, zx, xy
> O.step()
> bo1s.aabbEnlargeFactor=ig2ss.interactionDetectionFactor=1.
> O.run(); O.wait()
> plot.plot(subPlots=False)
>
>
>
>
>
>
>
> --
> 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 : yade-users@lists.launchpad.net
> Unsubscribe : https://launchpad.net/~yade-users
> More help : https://help.launchpad.net/ListHelp
>
>
--
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 : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
