yade-users team mailing list archive

[weijie <question693836@xxxxxxxxxxxxxxxxxxxxx>]

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

weijie gave more information on the question:
Hi Vasileios, and thank you again.

I consider which k, r, R value is more suitable for particles from the
perspective of modeling.This is my MWE[1].I found that the particles
created with this are very different from the particles of Potential
Blocks with the same parameters, and it seems that the volume and
Principal inertia cannot be directly calculated with b.shape.volume,
b.shape.inertia.

Best regards,
Jie

[1]
distanceToCentre= 0.025
b=Body()
b.mask=1
b.aspherical=True
wire=False
color=Vector3(random.random(),random.random(),random.random())
highlight=False
aa= [-0.7295950392572228, 0.46966065591120043, 0.978755366731517, 0.4159312753165438, -0.7851930775036051, 0.32622712056941044, 0.2487629499110137, -0.32802297814039366, 0.22757387395233566, -0.34064170979743924, 0.6535371560905189, -0.534030672328059, -0.9483691932498606]
bb= [-0.5797330029669138, -0.8677701055770507, 0.0908583947065325, -0.4698271127991969, 0.37706716796452805, 0.04463765109291187, 0.9516154740960145, -0.9017571077715157, 0.6925524728306885, -0.10267768395748879, 0.3532767701547774, 0.8086775588116711, -0.0540657356796434]
cc= [-0.3627681407761932, -0.1624620329672607, -0.18380066432307973, 0.7786293458972144, -0.49121500575982796, -0.9442369119611344, 0.18040228438849618, 0.28148720112193015, -0.6845299148104668, 0.9345696971138125, 0.6693912975817099, 0.24668167116191958, 0.3125264301143626]
dd= [0.01697034586028, 	0.01248528317756, 	0.01173016874627, 	0.01854063787194, 	0.01405551559515, 	0.01349450583878 ,	0.01047944143112, 	0.01105292067252 ,	0.02109918761652, 	0.01609025523102 ,	0.01787444301714,0.01542370609022 ,	0.01180116564426   ]
r=min(dd)/2 
b.shape=PotentialParticle(k=0.1, r=r, R=0.5*R, a=aa, b=bb, c=cc, d=dd, isBoundary=False, color=color, wire=wire, highlight=highlight, minAabb=sqrt(3)*Vector3(distanceToCentre,distanceToCentre,distanceToCentre), maxAabb=sqrt(3)*Vector3(distanceToCentre,distanceToCentre,distanceToCentre), maxAabbRotated=1.02*Vector3(distanceToCentre,distanceToCentre,distanceToCentre), minAabbRotated=1.02*Vector3(distanceToCentre,distanceToCentre,distanceToCentre), AabbMinMax=True, id=count)
V=(2*distanceToCentre)**3 
geomInert=(1./6.)*V*(2*distanceToCentre)**2 
utils._commonBodySetup(b, V, Vector3(geomInert,geomInert,geomInert), material='frictionless', pos=[0,0,0], fixed=False)
b.state.pos =[0,0,0]
b.state.ori = Quaternion((random.random(),random.random(),random.random()),random.random()) 
O.bodies.append(b)

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