Can I add properties and color into export text file?
Hi,
I am trying to save the position, properties, and color. Can I export it into one text file for future importing uses for example I will import the spherer packing from text when I implement gravity deposition?
##My Code##
import random
import math
from yade import geom, pack, utils, plot, ymport, export
import numpy as np
# Materials Properties
gravel = CohFrictMat(young = 1e7, poisson = 0.25, density = 2700, label = 'gravel')
asphalt_binder = CohFrictMat(young = 1e7, poisson = 0.25, density = 1060, frictionAngle = radians(40), normalCohesion = 2e5, shearCohesion = 2e5, label = 'asphalt_binder')
# add properties
O.materials.
O.materials.
# Define cylinder with funnel parameters
center = (0, 0, 0)
diameter = 0.102
height = 0.18
# create cylindrical body with radius 0.102 m and height 0.064 m
cylinder = geom.facetCylin
# add cylinder to simulation
O.bodies.
# Define cylinder with funnel parameters
center1 = (0,0,height/2)
dBunker = 0.4
dOutput = 0.102
hBunker = 0
hOutput = 0.15
hPipe = 0
# create funnel as a bunker with diameter 0.102 m, height 0.064 m
funnel = geom.facetBunke
# add funnel to simulation
O.bodies.
def makeRandomBlock
bunkerLimit = diameter/2
#
# Make grid blocks inside a bunker
#
grid = np.arange(
inGrid = []
for yi in grid:
for xi in grid:
c = [[xi, yi], [xi+gridsize, yi], [xi+gridsize, yi+gridsize], [xi, yi+gridsize]]
#
# Check if the block corners are outside of the bunker or not
#
out = False
for p in c:
r = (p[0]**2 + p[1]**2)**0.5
if r > bunkerLimit:
#
# If the block is inside the bunker, keep it
#
if not out:
layer = math.ceil(numblocks / len(inGrid))
blocks = []
for l in range(layer):
for g in inGrid:
zi = minz + l*gridsize
p1 = g[0].copy()
p2 = g[2].copy()
random.
return blocks
minZ = 2.35
dbunker = 0.4
gridSize = dbunker/8
numblocks = 51
blockList = makeRandomBlock
for i in range(numblocks):
print("Making cloud block", i+1, "/", numblocks)
corner = blockList.pop()
sp = pack.SpherePack()
# 15.75 mm 13 particles
if (i < 13):
n1 = sp.makeCloud(
# 11 mm 51 particles
n2 = sp.makeCloud(
# 7.125 mm 51x11 = 561 particles
n3 = sp.makeCloud(
# 3.555 mm 51x100 = 5,100 particles
n4 = sp.makeCloud(
# 1.77 mm 51x360 = 18,360 particles
n5 = sp.makeCloud(
# 0.6 mm 51x19000 = 969,000 particles
n6 = sp.makeCloud(
# 0.3 mm 51x78510 = 4,004,010 particles
n7 = sp.makeCloud(
sp.
for body in O.bodies:
if not isinstance(
continue
if body.shape.radius == 0.01575/2 :
if body.shape.radius == 0.011/2:
if body.shape.radius == 0.007125/2:
if body.shape.radius == 0.003555/2:
if body.shape.radius == 0.00177/2 :
if body.shape.radius == 0.0006/2 :
if body.shape.radius == 0.0003/2 :
# give properties to sphere
for i in range(n1,n6):
O.bodies[i].mat = gravel
for i in range(n6, n7):
O.bodies[i].mat = asphalt_binder
# Save body IDs and materials to a text file
with open("body_
for b in O.bodies:
export.
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Huan
- Solved:
- Last query:
- Last reply: