Re: [Yade-users] [Question #683949]: particle position data
Question #683949 on Yade changed:
https://answers.launchpad.net/yade/+question/683949
Status: Needs information => Open
Jinny Kim gave more information on the question:
Hello, Jan Stránský,
I am sorry for this inconvenience. I have attached the simple version of
my simulation code below. It also have the same issue with mine It
looks like having weird data at iteration of 14, 25, 36, 47 (I am
using Yade 1.20.0.)
#
import math
import numpy as np
import sys
import os.path
def saveData():
global ids
global i
p_num=ids+1
pkdID=np.zeros((p_num,1))
p_radii=np.zeros((p_num,1))
x_pos=np.zeros((p_num,1))
y_pos=np.zeros((p_num,1))
z_pos=np.zeros((p_num,1))
x_vel=np.zeros((p_num,1))
y_vel=np.zeros((p_num,1))
z_vel=np.zeros((p_num,1))
x_ang=np.zeros((p_num,1))
y_ang=np.zeros((p_num,1))
z_ang=np.zeros((p_num,1))
p_mass=np.zeros((p_num,1))
eflag=np.zeros((p_num,1))
print(p_num)
for k in range(p_num):
pkdID[k][0]=O.bodies[k].id
x_pos[k][0]=O.bodies[k].state.pos[0]
y_pos[k][0]=O.bodies[k].state.pos[1]
z_pos[k][0]=O.bodies[k].state.pos[2]
x_vel[k][0]=O.bodies[k].state.vel[0]
y_vel[k][0]=O.bodies[k].state.vel[1]
z_vel[k][0]=O.bodies[k].state.vel[2]
x_ang[k][0]=O.bodies[k].state.angVel[0]
y_ang[k][0]=O.bodies[k].state.angVel[1]
z_ang[k][0]=O.bodies[k].state.angVel[2]
p_mass[k][0]=O.bodies[k].state.mass
p_radii[k][0]=O.bodies[k].shape.radius
f1=open(os.path.join("/home/ae_h1/yzk0056/Desktop/Semesters/2019 Summer
Semester/DEM
Project/DEMsimulation/Obtain_initial_position/Sim_regular_rotation_Apophis/constant_gravity_2/","Output_"+str(i)+".txt"),"w")
for p in range(p_num):
f1.write("%d\t%e\t%e\t%e\t%e\t%e\t%e\t%e\t%e\t%e\t%e\t%e\n" %
(pkdID[p][0],p_mass[p][0],p_radii[p][0],x_pos[p][0],y_pos[p][0],z_pos[p][0],x_vel[p][0],y_vel[p][0],z_vel[p][0],x_ang[p][0],y_ang[p][0],z_ang[p][0]))
f1.close()
i+=1
Density=2900
Young=3.3e7
poisson=0.3
c_r = 0.1
r_fc = 0.0 #1st case: low angle of repose
frictionAngle=radians(35)
sphereColor = (.8,.8,0.)#dirty yellow
i=0
from yade import pack, plot, ymport, export
## Import particles
iron_ore=FrictMat(young=Young,poisson=poisson,density=Density,frictionAngle=frictionAngle)
Mat=O.materials.append(iron_ore)
from yade import pack, plot, ymport, export
sp=pack.SpherePack()
ids=sp.makeCloud(minCorner=(-0.4,-0.2,0.01),maxCorner=(0.4,0.2,1),rMean=2.0e-2,rRelFuzz=0.,num=25,periodic=False)
sp.toSimulation(material=Mat)
for pp in O.bodies:
ids=pp.id
## Import wall's geometry
O.materials.append(FrictMat(young=Young,poisson=poisson,density=Density,frictionAngle=frictionAngle,label='frictionlessWalls'))
fctIds =
O.bodies.append(geom.facetBox((0.,0.,-0.01),(0.8,0.4,0.0),wallMask=0b110011,color=(1,0,0),wire=False))
g=6.5e-2
itr_num=3
## Engines
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_MindlinCapillaryPhys(krot=r_fc,en=c_r,es=c_r)],
[Law2_ScGeom_MindlinPhys_Mindlin(includeMoment=True)]
),
NewtonIntegrator(damping=0,gravity=(0,0,-g),label='NI'),
PyRunner(iterPeriod=itr_num,command='saveData()')
]
## time step (how to run the script)
O.dt=2.0e-4 # should be 1.0e-6
O.run(1440)
from yade import qt
qt.View()
#
Thank you for your help in advance.
Sincerely,
Jinny Kim
--
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
[Yade-users] [Question #683967]: GlobalStiffnessTimeStepper with Potential Blocks
New question #683967 on Yade: https://answers.launchpad.net/yade/+question/683967 Hi, I want to use the GlobalStiffnessTimeStepper in conjunction with the Potential Blocks code. Currently, this engine assumes that the stiffnesses "phys->kn" and "phys->ks" are in [N/m] units [1]. In my case, in the KnKsPBLaw.cpp the "phys->kn" and "phys->ks" are expressed in [N/m^3], while when we multiply them with the contact area, we get the stiffness parameters: "phys->Knormal_area" and "phys->Kshear_area", which are in [N/m] units [2, 3]. I would like your opinion whether I could modify the GlobalStiffnessTimeStepper to consider the "phys->Knormal_area" instead of the "phys->Kn" (and respectively for the shear), when the particle shape is "PotentialBlock" instead of "Sphere". Do you think this would slow down the engine or find any theoretical drawbacks? Cheers, Vasileios [1] https://gitlab.com/yade-dev/trunk/blob/master/pkg/dem/GlobalStiffnessTimeStepper.cpp#L63 [2] https://gitlab.com/yade-dev/trunk/blob/master/pkg/dem/KnKsPBLaw.cpp#L121 [3] https://gitlab.com/yade-dev/trunk/blob/master/pkg/dem/KnKsPBLaw.cpp#L186 -- 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 #683946]: How to update yade (stable) version
Question #683946 on Yade changed: https://answers.launchpad.net/yade/+question/683946 Status: Answered => Solved Leonard confirmed that the question is solved: Hi A Dinesh, Yaniv and Jan, Thank you guys, all your instruction solved my problem. Cheers Leonard -- 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 #683946]: How to update yade (stable) version
Question #683946 on Yade changed: https://answers.launchpad.net/yade/+question/683946 Jan Stránský proposed the following answer: If you use virtual systems, you can just install appropriate Ubuntu (19.04 as pointed by Yaniv) and install "normal" yade [1] cheers Jan [1] https://launchpad.net/ubuntu/disco/+source/yade -- 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 #683687]: Unable to save/reload FlowEngine scene
Question #683687 on Yade changed: https://answers.launchpad.net/yade/+question/683687 Robert Caulk posted a new comment: Not advised * -- 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 #683687]: Unable to save/reload FlowEngine scene
Question #683687 on Yade changed: https://answers.launchpad.net/yade/+question/683687 Status: Open => Solved Robert Caulk confirmed that the question is solved: Currently yade is unable to save and reload flowengine triangulations. Therefore it is i advised to use this method. -- 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 #683949]: particle position data
Question #683949 on Yade changed: https://answers.launchpad.net/yade/+question/683949 Status: Open => Needs information Jan Stránský requested more information: Hello, please provide a complete MWE [1], otherwise we can just guess.. thanks Jan [1] https://www.yade-dem.org/wiki/Howtoask -- 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 #683934]: Bond broken and non movement of cohesive particles
Question #683934 on Yade changed: https://answers.launchpad.net/yade/+question/683934 Status: Open => Needs information Jan Stránský requested more information: Hello, > https://drive.google.com/drive/folders/1bCS_w917aP0Tr_Jfy6YfN20c9W0cdJYM?usp=sharing > Note: Kindly download the particle file using the download link above to > impor the particles. Note: Kindly make your question self-contained [1] (also see below) Try to make a MWE [1], I think something like 8 + 8 predefined particles for lower and upper part should be enough to illustrate the problem. Often you can find the problem yourself creating such MWE. cheers Jan [1] https://www.yade-dem.org/wiki/Howtoask, section "Please, no external links!" -- 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 #683819]: Error during installation: libyade.so
Question #683819 on Yade changed: https://answers.launchpad.net/yade/+question/683819 Jan Stránský posted a new comment: > Actually, I tried to install 2019 at first, but I encountered the problem that I described above. Then I turned to a different version (2018.02b) and the same problem occurred. my note was about trying different version of **operating system**, not Yade. In your case, you have quite new Yade version running on quite old Ubuntu. There could be some mismatch with used libraries (new Yade expecting something which is not present yet on the system). The problem is, that nobody will happy to trying this case, not saying about fixing the issue.. cheers Jan -- 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
[Yade-users] [Question #683952]: Cleaning up Shediac a week after Dorian
New question #683952 on Yade: https://answers.launchpad.net/yade/+question/683952 Like many parts of New Brunswick, Shediac is on the mend from Dorian's visit.The wind wreaked havoc at the marina, leaving a tangle of boats piled into each other, after a couple of docks let go during the storm. It's a different scene at the club a week later. All but one boat is out of the water and club members met on Saturday to discuss what to do next. Just about all of the boats have been removed from the water at marina. There were about 90 boats in the water when Dorian hit. (Gary Moore/CBC) Gerry O'Brien, the manager at the Shediac Bay Yacht Club, said about 120 members of the club attended the meeting. They discussed what engineering would need to be done before the club could move to the next step.O'Brien said the board of directors gave presentations about what happened before and after the storm. -- 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
