yade-users team mailing list archive

[liukeqi <question404012@xxxxxxxxxxxxxxxxxxxxx>]

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

Description changed to:
Hi,
  I am attempting to use paraview to do some post processing of a model using cpm material. 
  
  I have some confusion as following.

[1] I want to output some results that can indicate the damage condition
of the model. So which information should I be concerned, such as
‘cmpDamage’,  ‘cmpStress’ or ‘forceN’? What the specific meaning of
cpmDamage and cmpStress are?

[2] I used the paraview to display the data set ‘spheres‘ and chosen cpmDamage following the 'Coloring' button. But the paraview appear error information such as ‘Warning: Range [1,10] invalid for log scaling.  Changing to [0,0]. The cpmDamage information did not save, why? 
If I chosen cpmStress, it appears 'Warning: Range [0,9.31449e-05] invalid for log scaling.  Changing to [9.31449e-06,9.31449e-05].'. Can not paraview display the negative values?

[3] Can the paraview output the orientation distribution of sphere? If
it can, How to do it?

[4] Can you give me some post processing example of cpm material?

Thank you.
Liu keqi

The following is my code.
**********************************************************
#!/usr/bin/python           # This is server.py file
import string
from yade import plot,qt
from yade.pack import *
from yade import pack, plot


O.materials.append(CpmMat(young=160e9,density=4800,frictionAngle=atan(0.8),poisson=0.1,sigmaT=12e8,epsCrackOnset=0.4e-4,relDuctility=500))

initSize=0.4

sp=pack.randomPeriPack(radius=0.0175,initSize=Vector3(initSize,initSize,initSize),rRelFuzz=0.7)
O.periodic=True
sp.toSimulation()



O.dt=1e-7


plot.plots={'ey':('sy',)}
def plotAddData():
	plot.addData(
		sy=p3d.stress[1],
		ey=p3d.strain[1],
		ex=p3d.strain[0]
	)

EnlargeFactor=1.5
EnlargeFactor=1.0
O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=EnlargeFactor,label='bo1s')]),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=EnlargeFactor,label='ig2ss')],
		[Ip2_CpmMat_CpmMat_CpmPhys()],[Law2_ScGeom_CpmPhys_Cpm()]),
	NewtonIntegrator(),
	Peri3dController(
					
							nSteps=4000, 		
							doneHook='print "Simulation with Peri3dController finished."; O.pause()',
							label='p3d'
							),
	PyRunner(command='plotAddData()',iterPeriod=1),
	VTKRecorder(fileName='3d-vtk-',recorders=['all','cpm'],iterPeriod=1000),
]


p3d.stressMask=0b100001       
p3d.goal=(0,-0.008,0,0,0,0)  
O.step()
bo1s.aabbEnlargeFactor=ig2ss.interactionDetectionFactor=1.
O.run(); O.wait()

plot.plot(subPlots=False)

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