yade-users team mailing list archive

[procrastinator <question684615@xxxxxxxxxxxxxxxxxxxxx>]

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

hi all,

I want to get the particle stress at a specific region. but I get an error, can anyone give me some suggestions.

here is the MWE.



from yade import pack, plot, export
O.bodies.append(geom.facetBox((.5,.5,.5),(.5,.5,.5),wallMask=31))
sp=pack.SpherePack()
sp.makeCloud((0,0,0),(1,1,1),rMean=.08,rRelFuzz=.5)
sp.toSimulation()
O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
		[Ip2_FrictMat_FrictMat_FrictPhys()],
		[Law2_ScGeom_FrictPhys_CundallStrack()]
	),
	PyRunner(command='checkUnbalanced()',realPeriod=2),
	PyRunner(command='addPlotData()',iterPeriod=100),
	PyRunner(command='subbox()',iterPeriod=100),
	PyRunner(command='stress_export()',iterPeriod=100),
	NewtonIntegrator(gravity=(0,0,-9.81),damping=0.4),
]
O.dt=.5*PWaveTimeStep()
print(len(O.bodies))
def checkUnbalanced():
	if unbalancedForce()<.05:
		O.pause()
		#plot.saveDataTxt('bbb.txt')
def addPlotData():
	plot.addData(i=O.iter,unbalanced=unbalancedForce())
#plot.plots={'i':('unbalanced')}
#plot.plot()
######################## sub-region and coloring#############################
ball_list1 =[]
ball_list2 =[]
ball_list3 =[]
def subbox():
	for b in O.bodies:
		if 0< b.state.pos[2] <= 0.15:
			b.shape.color =Vector3(255,255,255) # white
			m = b.id
			ball_list1.append(m)
		elif 0.15 < b.state.pos[2] <= 0.3:
			b.shape.color =Vector3(0,100,000) # green
			n = b.id
			ball_list2.append(n)			
		elif    b.state.pos[2] > 0.3:
			b.shape.color =Vector3(00,00,255) # blue
			q = b.id
			ball_list3.append(q)
######################## sub-region and coloring#############################

######################################## get the per-perticle stess#####################

TW=TesselationWrapper()
TW.setState()
TW.computeVolumes()
vtk = export.VTKExporter("stress_field")
def stress_export():
	s=bodyStressTensors()
	for i in ball_list3:
		b = O.bodies(i)
		if isinstance(b.shape,Sphere):
			b.mystress = s[b.id]*4.*pi/3.*b.shape.radius**3/TW.volume(b.id)	
	vtk.exportSpheres(ids='ball_list3',what=[('pos','b.state.pos'),('stress','b.mystress')],useRef=False)
######################################## get the per-perticle stess#####################


when I run the code, here is the error:


TypeError                                 Traceback (most recent call last)
/usr/bin/yade in <module>()

/usr/bin/yade in stress_export()
     60         s=bodyStressTensors()
     61         for i in ball_list3:
---> 62                 b = O.bodies(i)
     63                 if isinstance(b.shape,Sphere):
     64                         b.mystress = s[b.id]*4.*pi/3.*b.shape.radius**3/TW.volume(b.id)

TypeError: 'BodyContainer' object is not callable


thanks,
yong








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