← Back to team overview

yade-users team mailing list archive

Re: [Question #235138]: how to define friction at the boundaries in triaxial test

 

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

    Status: Open => Answered

Christian Jakob proposed the following answer:
in your script i see some inconsistencies, that can lead to unexpected
results:

##### 1. problem:

idSand=O.materials.append(FrictMat(young=young,poisson=pr,frictionAngle=radians(fa),label='idSand'))
#frictionAngle=fa !!!

##### 2. problem:

O.bodies.append(geom.facetBox((0.035,0.035,0.035),(0.035,0.035,0.035),wallMask=15, material = BoundaryMat, dynamic=1))
#this is ok

O.bodies.append(geom.facetBox((0.035,0.035,0.035),(0.035,0.035,0.035),wallMask=15, dynamic=1))
#missing material = BoundaryMat !!!

O.bodies.append(utils.geom.facetBox((0.035,0.035,0.035),(0.035,0.035,0.035),wallMask=16,dynamic=0))
#missing material = BoundaryMat !!!

O.bodies.append(utils.geom.facetBox((0.035,0.035,0.105),(0.035,0.035,0.035),wallMask=15))
#missing material = BoundaryMat

##### 3.problem:

sp=pack.SpherePack()
sp.makeCloud(mincorner,maxcorner,psdSizes=psdSizes,num=number,psdCumm=psdCumm,distributeMass=1)

sp.toSimulation()

you did not set material for spheres. but i do not know how to do this with sp.toSimulation()
i prefer to use:

O.bodies.append([sphere(c,r,material=idSand) for c,r in sp])
#instead of sp.toSimulation()

regards,

christian

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.