← Back to team overview

yade-users team mailing list archive

[Question #660540]: Minimum and Maximum values for parameters

 

New question #660540 on Yade:
https://answers.launchpad.net/yade/+question/660540

Hi guys!

I am trying to calibrate the model with the triaxial test.
I used CohFricMat as a material and Cohesionmoment for the Contact law.
The soil is sand with cohesion.
So I have to work with 5 parameters to match the result of the experimental test with the triaxial test.
These parameters are E, Ks/Kn, Phi, eTaRoll and alphaKr.
Hence, I have a general question about YADE which Are there any limitation for parameters? it means is there any minimum and maximum value for contact (i.e. CohFrictMat) parameters?

For instance, I used these value for my model:

from yade import pack,plot

############################################
### DEFINING VARIABLES AND MATERIALS ###
############################################

# The following 5 lines will be used later for batch execution
nRead=readParamsFromTable(
 num_spheres=5000,# number of spheres
 compFricDegree =30, # contact friction during the confining phase
 key='_triax_base_', # put you simulation's name here
 unknownOk=True
)
from yade.params import table

num_spheres=table.num_spheres# number of spheres
key=table.key
targetPorosity = 0.49 #the porosity we want for the packing
compFricDegree = 30.6# contact friction during the deviatoric loading
finalFricDegree = 10
rate=-0.1 # loading rate (strain rate)
damp=0.70 # damping coefficient
stabilityThreshold=0.01 # we test unbalancedForce against this value in different loops (see below)
young=1e6# contact stiffness
mn,mx=Vector3(0,0,0),Vector3(0.09,0.18,0.09) # corners of the initial packing

## create materials for spheres and plates
O.materials.append(CohFrictMat(alphaKr=40,young=young,poisson=0.2,frictionAngle=radians(30),normalCohesion=0,shearCohesion=0,momentRotationLaw=False,etaRoll=0.05,density=2600,isCohesive=False,label='spheres'))
O.materials.append(CohFrictMat(young=young,poisson=0.2,frictionAngle=radians(30),density=0,label='walls'))



Cheers
Hessam

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