← Back to team overview

yade-users team mailing list archive

Re: [Question #687871]: Outputs of the Triaxial code by Bruno Chareyre

 

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

    Status: Answered => Solved

ehsan benabbas confirmed that the question is solved:
Jan and Robert, thanks for your feedback.

When I want to use basic units, it implicitly means that I should make
them consistent in my mind to a specific unit. For example, make them
consistent based on "mm". Exactly like what I did. Accordingly, I made
changes for "m" and "N" as follows:

Sima iso = 300 Kpa = 3e5 N/m2
Porosity = 0.4
Nu = 0.5
Phi = 30
Density = 2000 Kg/m3
E = 5e7 pa

Which needs to make these changes in the code:

num_spheres=20000     
targetPorosity = 0.4
young=5e7
mn,mx=Vector3(0,0,0),Vector3(0.02,0.02,0.02)
O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=radians(compFricDegree),density=2000,label='spheres'))
O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=0,density=0,label='walls'))
sp.makeCloud(mn,mx,0.0002,0.5,num_spheres,False, 0.95,seed=1)
triax.goal1=triax.goal2=triax.goal3=-300000
triax.goal1=-300000
triax.goal3=-300000

In this case, it takes about 2 mins to get the results. But plots are
not conventional. Would you please make these changes and run the
modified [1] and let me know your feedback?

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.