yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #21692
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.