yade-users team mailing list archive

[Yuri Bezmenov <question707210@xxxxxxxxxxxxxxxxxxxxx>]

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

Yuri Bezmenov gave more information on the question:
Further I tried with mask option but it export all the bodies in
simulation

MWE:

from builtins import range
from yade import pack, geom, plot, export
from yade.gridpfacet import *
rate=-0.5 #1.66e-4
compFricDegree=30
finalFricDegree = 40
stabilityThreshold = 0.1
young = 5e6
mn, mx = Vector3(0, 0, 0), Vector3(0.9, 0.6, 0.6)
O.engines = [
        ForceResetter(),
        InsertionSortCollider([
                Bo1_Sphere_Aabb(),
                Bo1_GridConnection_Aabb(),
		Bo1_Facet_Aabb(), Bo1_Wall_Aabb()
        ]),
        InteractionLoop(
                [       Ig2_Sphere_Sphere_ScGeom(),
			Ig2_Facet_Sphere_ScGeom(), Ig2_Wall_Sphere_ScGeom(), 
                        Ig2_Sphere_GridConnection_ScGridCoGeom(),
                        Ig2_GridNode_GridNode_GridNodeGeom6D(),
                        Ig2_GridConnection_GridConnection_GridCoGridCoGeom()], 
		[Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True, setCohesionOnNewContacts=True), Ip2_FrictMat_FrictMat_FrictPhys()], 
		[Law2_ScGeom_FrictPhys_CundallStrack(), Law2_ScGridCoGeom_FrictPhys_CundallStrack(), Law2_GridCoGridCoGeom_FrictPhys_CundallStrack(), 
		Law2_ScGeom6D_CohFrictPhys_CohesionMoment(),]
        ),
        GlobalStiffnessTimeStepper(timestepSafetyCoefficient=0.1, label='ts'),
        NewtonIntegrator(gravity=(0,-10,0), damping=.5, label='newton'),
	PyRunner(command='addPlotData()', iterPeriod=500),
	#PyRunner(command='expo()', iterPeriod=500),
	VTKRecorder(fileName='/files/Geogrid/paraview/pullout', recorders=['all'], iterPeriod=500, mask=10),
	VTKRecorder(fileName='/files/Geogrid/paraview1/pullout', recorders=['all'], iterPeriod=500, mask=11),

]
def addPlotData():
	plot.addData(
		unbalanced=unbalancedForce(), 
		i=O.iter
	)
plot.plots = {
        'i':('unbalanced'),
}
plot.plot()
O.materials.append(FrictMat(young=young, poisson=0.5, frictionAngle=radians(compFricDegree), density=2600, label='spheres'))
O.materials.append(FrictMat(young=young, poisson=0.5, frictionAngle=0, density=0, label='walls'))
O.materials.append(CohFrictMat(young=1.325e8, poisson=0.40, density=900, frictionAngle=radians(10), normalCohesion=1.1e8, shearCohesion=1.1e8, momentRotationLaw=True, label='spheremat'))
L = 0.608  #length [m]
l = 0.304  #width [m]
nbL = 16+1  #number of nodes for the length	[#]
nbl = 8+1  #number of nodes for the width	[#]
r = 0.008/2	# L / 100.  #radius
color = [255. / 255., 102. / 255., 0. / 255.]
nodesIds = []
for i in range(0, nbL):
	for j in range(0, nbl):
		nodesIds.append(O.bodies.append(gridNode([i*L/nbL, 0.3, 0.15+j*l/nbl], r, wire=False, fixed=False, material='spheremat', color=color)))
a=[]
b=[]
#Create connection between the nodes
for i in range(0, len(nodesIds)):
	for j in range(i + 1, len(nodesIds)):
		dist = (O.bodies[i].state.pos-O.bodies[j].state.pos).norm()
		if (dist <= (L/nbL)*1.05):
			O.bodies.append(gridConnection(i, j, r, color=color))
			a=a+[(i,O.bodies[-1].id), (O.bodies[-1].id,j)]
			b=b+[(i,j)]
O.bodies.append(geom.facetBox((.45, .5, .3), (.45, .5, .3), wallMask=55, material='walls'))
sp = pack.SpherePack()
sp.makeCloud((0, 0, 0), (0.9, 0.2, 0.6), rMean=.05, rRelFuzz=.3)
sp.makeCloud((0, 0.4, 0), (0.9, 0.6, 0.6), rMean=.05, rRelFuzz=.3)
O.bodies.append([sphere(center, rad, material='spheres') for center, rad in sp])
for b in O.bodies:
	if b.id<153:
		b.groupMask=10
	else:
		if b.id<433:
			b.groupMask=11
		else:
			b.groupMask=12

#vtkExporter = export.VTKExporter('/home/harshal/files/Geogrid/paraview1/geogrid')
#def expo():
	#vtkExporter.exportSpheres(what=dict(dist='b.state.pos.norm()',pos='b.state.pos'), ids=range(0,433)) 
	#vtkExporter.exportInteractions(what=dict(kn='i.phys.kn', kr='i.phys.kr', ks='i.phys.ks'), ids=b)   ###Here I replace a and b 
	# vtkExporter.exportContactPoints(what={'nn': 'i.geom.normal'})
	# vtkExporter.exportPolyhedra(what=dict(n='b.state.pos.norm()'), ids=plate.id)
O.run(10000,True)
for i in range (0,nbl):
	O.bodies[i].state.blockedDOFs='x'
	O.bodies[i].state.vel = (rate,0,0)
O.run(15000,True)

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