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