yade-users team mailing list archive

[behzad <question258859@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi there,

It still have an issue with clump mass calculation with replaceByClumps(). 


id_Mat=O.materials.append(FrictMat(young=1e7,poisson=0.3,density=1000,frictionAngle=1))
Mat=O.materials[id_Mat]

sp=pack.SpherePack()
sp.makeCloud(minCorner=(-1.5,-1.5,.1),maxCorner=(1.5,1.5,2),rMean=.2,rRelFuzz=0.0,num=1,periodic=False)
O.bodies.append([sphere(c,r,material=Mat) for c,r in sp])

#Coke:
relRadList4=[0.000471583e3,0.00060213e3,0.000631454e3,0.000493144e3,0.000695336e3,0.000395284e3,0.000550567e3,0.00037674e3,0.000656561e3,
0.000218143e3,0.000223094e3,0.000484354e3,0.00036813e3,0.000514327e3,0.000395357e3,0.000377391e3,0.000454081e3,0.000355504e3,
0.00061375e3,0.00068747e3,0.000508243e3,0.000447592e3,0.000549008e3,0.000486356e3,0.000664425e3,0.000587608e3,0.00015184e3,
0.000695789e3,0.000484995e3,0.000443755e3]
relPosList4=[[0.000213504e3,0.000920219e3,5.17294e-002],[0.000181129e3,-0.000693148e3,-0.000301356e3],[0.000191565e3,0.000550929e3,0.000472726e3],[-6.9806e-002,0.000814473e3,0.000504795e3],[-0.00016687e3,-4.19741e-002,0.000273102e3],[-0.000632058e3,-0.000285906e3,0.000255655e3],[8.5103e-002,-2.87712e-003,0.000608655e3],[-0.000369852e3,6.12296e-003,-0.00056701e3],[7.38905e-002,0.000497426e3,2.97204e-002],[0.000135818e3,0.000562643e3,-0.000842692e3],[0.000748387e3,0.00077499e3,0.000784656e3],[7.88683e-002,0.000626452e3,-0.00032004e3],[0.000526651e3,-0.000553811e3,0.000109479e3],[0.000136113e3,0.000176874e3,-0.000254769e3],[-0.000209755e3,-0.000765094e3,-0.000832511e3],[-0.000945048e3,-0.00034182e3,-0.000150738e3],[0.000425963e3,-0.000381421e3,-0.000472185e3],[0.000544025e3,0.000492145e3,0.000918189e3],[-0.000262527e3,0.000319288e3,0.000244014e3],[-0.0002573e3,-0.000432994e3,-0.000293022e3],[0.000500309e3,-0.000741747e3,-0.000145635e3],[0.000354013e3,4.94286e-002,0.000772469e3],[7.29082e-002,-0.000364611e3,-0.000799481e3],[0.00015904e3,0.000954791e3,0.00058286e3],[0.000106567e3,-0.000252692e3,2.9281e-002],[-0.000616653e3,-0.000235914e3,-0.000163413e3],[-0.00114642e3,8.90673e-002,-0.000146111e3],[-0.000361004e3,-0.000311422e3,-0.000152367e3],[-0.000655366e3,-2.8257e-002,-0.000169021e3],[-0.000501808e3,-6.06848e-002,-0.000385601e3]]


templates= []
templates.append(utils.clumpTemplate(relRadii=relRadList4,relPositions=relPosList4))
O.bodies.replaceByClumps(templates,[1.0])


As it's written in the manual, the volume of the clump must be the same as the ball which is replaced by the clump. 

The initial ball here has the radius of 0.2 which gives the volume of 0.0335. The material has the density of 1000, so mass of the ball is 33.5. 

However, mass of the clump is 99.808

O.bodies[32].state.mass
 ->  [14]: 99.80883330235272

How this is calculated? and why 99.80? Why it's not 33.5? Can we control it? 
and I didn't find any clue about adaptClumpMasses(). 

Thanks

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