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. _______________________________________________ Mailing list: https://launchpad.net/~yade-users Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
