[Yade-users] [Question #700830]: Is there any way to generate voronoi tessellation diagram?
New question #700830 on Yade: https://answers.launchpad.net/yade/+question/700830 Dear all, I'd like to generate voronoi tessellation diagram like the fig[1]. I tried the TW.defToVtk() and got Delaunay triangulation diagram but not voronoi tesselation diagram. It seems the voro++[2] is what I am looking for, but the result seems not good because I didn't see any spheres. Thanks for your any suggestion. [1]https://www.yade-dem.org/wiki/File:Cell_volume.png [2]http://math.lbl.gov/voro++/ -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #700754]: How are the monitoring points arranged
Question #700754 on Yade changed: https://answers.launchpad.net/yade/+question/700754 Status: Open => Answered Robert Caulk proposed the following answer: You decide what to monitor, how to monitor, when to monitor, and where to monitor. There are no predetermined "monitoring points." If you want, you can collect all data for all particles for all time steps (not recommended). -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #700813]: Batch run and run alone get different results
Question #700813 on Yade changed: https://answers.launchpad.net/yade/+question/700813 Status: Open => Needs information Robert Caulk requested more information: >However,when I run it through separate yadedaily,it stops when strain[2] is 0.0331 and peak strength is about 120MPa,when I run it through batch mode,it stops when strain[2] is 0.161 and peak strength is about 161MPa. How do these numbers change as you re-run them multiple times with exactly the same parameters? -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #700611]: TwoPhaseFlow - double free or corruption
Question #700611 on Yade changed: https://answers.launchpad.net/yade/+question/700611 Luis Barbosa gave more information on the question: Here is how I am running: + singularity run /beegfs/common/singularity/yade/yade_debian_bookwarm_1.0.sif yadedaily -j80 drain_importMatrixTriax7_10_c3.py Here is the output: double free or corruption (!prev) Interesting to note: 1. I have used tcmalloc, as suggested here: https://github.com/tensorflow/tensorflow/issues/6968 https://www.google.com/search?client=firefox-b-d&q=apt-get+install+google-perftools https://riptutorial.com/tensorflow/example/13427/use-the-tcmalloc-allocator But, it just gave me another output: "Core dumped" 2. >From the TwoPhaseFlowEngine I played around to see where exactly the error >happens and discovered that it is related to unsat.invasion(), Invasion 1 or >invasion 2 in the hpp file. In some scripts that I have double free or >corruption (!prev), when I switch .isPhaseTrapped=False, than it worked, but >in the script I provided above it still did not worked. 3. Sometimes, I am getting the message, when using Yadedaily, but not for Yade. <[36mNOFILTER[0m> :[93m750[32m void yade::CGT::Network::lineSolidPore(yade::CGT::Network::CellHandle, int) [with Tesselation = yade::CGT::_Tesselation >; yade::CGT::Network::CellHandle = CGAL::internal::CC_iterator, CGAL::Triangulation_cell_base_with_info_3, CGAL::Regular_triangulation_cell_base_3, CGAL::Triangulation_cell_base_3, CGAL::Triangulation_ds_cell_base_3, CGAL::Triangulation_vertex_base_with_info_3, CGAL::Regular_triangulation_vertex_base_3 > >, CGAL::Boolean_tag, CGAL::Boolean_tag >, CGAL::Alpha_shape_cell_base_3, CGAL::Triangulation_cell_base_with_info_3, CGAL::Regular_triangulation_cell_base_3 > >, CGAL::Boolean_tag, CGAL::Boolean_tag >, CGAL::Sequential_tag> > >, CGAL::Hidden_points_memory_policy >, std::__cxx11::list >, std::allocator > > > > >, CGAL::Boolean_tag, CGAL::Boolean_tag >, CGAL::Default, CGAL::Default, CGAL::Default>, false>][0m: Unhandled switch case:3 -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #700611]: TwoPhaseFlow - double free or corruption
Question #700611 on Yade changed:
https://answers.launchpad.net/yade/+question/700611
Status: Answered => Open
Luis Barbosa is still having a problem:
For this I have to provide txt files for the creation of the packs.
Sorry for the external links:
https://drive.google.com/file/d/16y5Ia4kUUm6oap8QlnSvQf9uCuMzJ7TN/view?usp=sharing
https://drive.google.com/file/d/1_rUoO-6b2j6e3LSRsMlh8C4vKRu8txoN/view?usp=sharing
https://drive.google.com/file/d/1oR9bHoru-
Us6wvHwn7UpkZwt9IX6GRvE/view?usp=sharing
https://drive.google.com/file/d/1cLyeXTjNuavGRxuFv27HyPxvtn-2EW5i/view?usp=sharing
=
from yade import pack
from yade import bodiesHandling
from yade import export
from yade import utils
from yade import ymport
import math
### DEFINING VARIABLES AND MATERIALS ###
# The following 5 lines will be used later for batch execution
nRead=readParamsFromTable(
num_spheres=3000,# number of spheres
compFricDegree = 1, # contact friction during the confining phase
key='_triax_base_', # put you simulation's name here
unknownOk=True
)
from yade.params import table
num_spheres=table.num_spheres# number of spheres
key=table.key
targetPorosity = 0.50 #the porosity we want for the packing
compFricDegree = table.compFricDegree # initial contact friction during the
confining phase (will be decreased during the REFD compaction process)
finalFricDegree = 30 # contact friction during the deviatoric loading
rate=0 # loading rate (strain rate)
damp=0.8 # damping coefficient
stabilityThreshold=0.01 # we test unbalancedForce against this value in
different loops (see below)
young=80e5 # contact stiffness200e4
young2=80e5
youngcoat=50e5
bondstr=1000#2e7
bondstr2=1000
bondstrcoat=10
## create materials for spheres and plates
mat=O.materials.append(JCFpmMat(type=1,young=young,poisson=0.3,frictionAngle=radians(compFricDegree),density=2000,tensileStrength=bondstr,cohesion=bondstr,jointNormalStiffness=0,jointShearStiffness=0,jointCohesion=bondstr,jointFrictionAngle=radians(0),jointDilationAngle=0.0,label='spheres'))
O.materials.append(JCFpmMat(type=1,young=20e7,poisson=0.3,frictionAngle=radians(0),density=2600,tensileStrength=0,cohesion=0,jointNormalStiffness=0,jointShearStiffness=0,jointCohesion=0,jointFrictionAngle=radians(0),jointDilationAngle=0.0,label='walls'))
O.materials.append(JCFpmMat(type=1,young=youngcoat,poisson=0.3,frictionAngle=radians(1),density=1500,tensileStrength=bondstrcoat,cohesion=bondstrcoat,jointNormalStiffness=0,jointShearStiffness=0,jointCohesion=bondstrcoat,jointFrictionAngle=radians(0),jointDilationAngle=0.0,label='spherescoat'))
## create walls around the packing
mn,mx=Vector3(0,0,0),Vector3(0.0015,0.0015,0.0015)
mnbox,mxbox=Vector3(0,0,0),Vector3(0.002,0.00195,0.002)
walls=aabbWalls([mnbox,mxbox],thickness=0,material='walls')
wallIds=O.bodies.append(walls)
O.bodies.append(ymport.textExt("matrix_vtest6.txt", format='x_y_z_r',
shift=Vector3(0,0.00025,0), scale=1.0,material='spheres',color=(0,1,1)))
O.bodies.append(ymport.textExt("coat_vtest5e5.txt", format='x_y_z_r',
shift=Vector3(0,0,0), scale=1.0,material='spherescoat',color=(0,1,1)))
Particle substitution by large
aggregate##
bodid=[]
for b in O.bodies:
if b and isinstance(b.shape,Sphere):
# print (b.shape.radius)
if b.state.pos[1]>0.00175:
bodid.append(b.id)
i=0
for p in bodid:
O.bodies.erase(bodid[i])
i=i+1
bodid=[]
a=[]
for b in O.bodies:# in sp:
if b and isinstance(b.shape,Sphere):
# print (b.shape.radius)
if b.shape.radius==0.0002:
bodid.append(b.id)
a.append(b.state.pos)
i=0
for p in bodid:
t=a[i]
f1=O.bodies.append(ymport.textExt("agg2e4_10e6.txt", format='x_y_z_r',
shift=t-Vector3(0,0,0.0002), scale=1.0,material='spheres',color=(0,1,1)))
O.bodies.erase(bodid[i])
i=i+1
bodiddd=[]
aaa=[]
for bbb in O.bodies:# in sp:
if bbb and isinstance(bbb.shape,Sphere):
# print (b.shape.radius)
if bbb.shape.radius==0.0005:
bodiddd.append(bbb.id)
aaa.append(bbb.state.pos)
iii=0
for ppp in bodiddd:
ttt=aaa[iii]
f3=O.bodies.append(ymport.textExt("agg5e5_18e6.txt", format='x_y_z_r',
shift=ttt-Vector3(0,0,0.0005), scale=1.0,material='spheres',color=(0,1,1)))
O.bodies.erase(bodiddd[iii])
iii=iii+1
##
#=
