yade-users team mailing list archive

[ytang <question695925@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi All,

I use the following command to get the position of particles and then export it.
###########################################################
	xposition = [] ; yposition = [] ; zposition = [] ; ballid = [] ; radius = []
	for b in O.bodies:
		if isinstance(b.shape,Sphere):		
			xposition.append( b.state.pos[0] )  ;   yposition.append( b.state.pos[1] )
			zposition.append( b.state.pos[2] )  ;   radius.append( b.shape.radius )  ;	ballid.append( b.id )
	f = open('./global-position'+'-' +str(m)+'.txt',"w")
	f.write('x\ty\tz\tr\tbid\n')
	f.close
	j = 0
	while j < len(xposition):
		f.write(str(xposition[j])+'\t'+str(yposition[j])+'\t'+str(zposition[j])+'\t'+str(radius[j])+'\t'+str(ballid[j])+'\n')
		f.close
		j = j +1
#################################################################################
after that I use the TW module to get the micro strain filed. I can get the micro strain field.
##########################################################################
    O .bodies.append(ymport.text("./global"+"-"+"position"+"-"+str(i)+".txt"))
    TW = TesselationWrapper()
    TW.triangulate()
    TW.computeVolumes()
    TW.setState(0)
    O.reset()
    O.bodies.append(ymport.text("./global"+"-"+"position"+"-"+str(i+1)+".txt"))
    TW.setState(1)
    TW.defToVtk("./strain"+"-"+"control"+str(i)+".vtk")
##################################################################################
when i import the .vtk files into the PARAVIEW.
we can see two potions to coloring the strain. one is strain_deviator, another one is the strain_matrix.
here I use the strain_deviator.
undering the coloring mode. I choose the: rescale to data range. the colorbar gives me the range from 0 to 4.8e+00. which means the maximum value of the strain deviator is around 4.8. when we see the offcial website for the micro strain. the maximum value of the strain deviator is 0.478597 [1]. I also saw this post[2]. it mentioned that the strain devitor is the  second invariant. 
In my opinion, this value should be smaller than one. but for my case, it's larger than one. is there anyone who can give me some hints?

I post the .vtk files here [3].
best,
Yong.



references: [1]https://yade-dem.org/doc/user.html?highlight=micro%20strain
[2] https://answers.launchpad.net/yade/+question/266331
[3]https://www.dropbox.com/sh/lal17mtarct8u51/AAApMr6QgxUY53UopAHtinGEa?dl=0




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