Particle Size distribution for fine-grained soil
Hi every body
I am trying to simulate triaxial condition on an fine-grained type of soil (CL-ML), first of all I have to match the particle size distribution according the real condition of the soil. In below I mention my code with real particle distribution in terms of psdSizes and psdCumm.
I wonder, if there is any way to match this PSD size?
I also changed the bin number but there is no improvement.
#######
from mpl_toolkits.
import numpy as np
import matplotlib.pyplot as plt
from yade import pack,qt
utils.readParam
from yade.params import table
seed=table.seed
num_spheres=
compFricDegree = table.compFricD
confiningS=-1e5
## creat a packing with a specific particle side distribution (PSD)
psdSizes,
sp=pack.
mn,mx=Vector3(
sp.makeCloud(
#sp.psd(
## create material #0, which will be used as default
O.materials.
O.materials.
## create walls around the packing
walls=aabbWalls
wallIds=
O.bodies.
triax=TriaxialS
internalCompac
goal1=confiningS,
goal2=confiningS,
goal3=confiningS,
max_vel=10,
label="triax"
)
newton=
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
#Finally, two different contact laws for sphere-box and sphere-sphere
[Law2_
useIncrement
always_
label=
),
GlobalStiffnes
triax,
newton
]
#from yade import qt
#qt.View()
file=open(
Diameter=[]
idsphere=[]
for center,rad in sp:
Diameter.
for i in range(0,
idsphere.
file.write(
file.write(
file.write(
file.close()
plt.figure()
plt.subplot(221)
plt.plot(
plt.semilogx(
plt.xlabel(
plt.ylabel(
plt.grid(True)
plt.legend(
plt.subplot(222)
plt.ylabel('Sphere Diameter(mm)')
plt.xlabel('Sphere id')
plt.grid(True)
plt.plot(
plt.show()
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Hossein
- Solved:
- Last query:
- Last reply: