Re: [Yade-users] [Question #680892]: O.engines[].dead=True
Question #680892 on Yade changed: https://answers.launchpad.net/yade/+question/680892 Status: Answered => Solved Francesco confirmed that the question is solved: Thank you Robert! Your idea looks to be good. I used this option for the translation of other bodies and it actually works. I'll do what you suggested me. Cheers, Francesco -- 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 #680892]: O.engines[].dead=True
Question #680892 on Yade changed: https://answers.launchpad.net/yade/+question/680892 Status: Answered => Open Francesco is still having a problem: Hi, I attached the code. I commented some functions related to data storing that could create problems. Thank you so much. Francesco Code: from numpy import linspace from yade import pack, qt import numpy as np import os import glob import re import smtplib from email.MIMEMultipart import MIMEMultipart from email.MIMEText import MIMEText from email.MIMEBase import MIMEBase from email import encoders ## # Parameters ## mi=10**-3 # Millimetres transformation mu=10**-6 # Micron transformation x1=10*mi# Box half side z1=10*mi# Box heigth dist=100*mi # Nozzle target distance h=z1+(dist) # Nozzle heigth from z=0 rIn=6*mi# Inner radius rOut=12*mi # Cylinder radius eps=0.1*mi # Tollerance itereps=5 rPack=rOut-eps # Pack radius rd=rOut-rIn # Difference between radius l2=80*mi# Cylinder length l_cloud=50*mi # Particles cloud length c1=(4*l_cloud)/9# External particles pack range c2=(1*l_cloud)/9# Middle particles pack range l3=1*mi # Distance from the nozle indentation v=57# Particles initial velocity vi_y=-v*cos(pi*angle/180) vi_z=-v*sin(pi*angle/180) vf1=57 # Fluid velocity 1 vf2=220 # Fluid velocity 2 rMean=40*mu # Mean particles radius Rho=3.61# Air density 1 Cd=0.47 # Sphere drag coefficient wt=10 # Box rotation angle=22.5 # Nozzle angle box_angle=45 vx=1 vz=1 ## # Nozzle creation ## delta=((pi*0.5)/16) poly=((rIn,h), (rOut-(rd*cos(delta)),h+(rIn*sin(delta))), (rOut-(rd*cos(2*delta)),h+(rIn*sin(2*delta))), (rOut-(rd*cos(3*delta)),h+(rIn*sin(3*delta))), (rOut-(rd*cos(4*delta)),h+(rIn*sin(4*delta))), (rOut-(rd*cos(5*delta)),h+(rIn*sin(5*delta))), (rOut-(rd*cos(6*delta)),h+(rIn*sin(6*delta))), (rOut-(rd*cos(7*delta)),h+(rIn*sin(7*delta))), (rOut-(rd*cos(8*delta)),h+(rIn*sin(8*delta))), (rOut-(rd*cos(9*delta)),h+(rIn*sin(9*delta))), (rOut-(rd*cos(10*delta)),h+(rIn*sin(10*delta))), (rOut-(rd*cos(11*delta)),h+(rIn*sin(11*delta))), (rOut-(rd*cos(12*delta)),h+(rIn*sin(12*delta))), (rOut-(rd*cos(13*delta)),h+(rIn*sin(13*delta))), (rOut-(rd*cos(14*delta)),h+(rIn*sin(14*delta))), (rOut-(rd*cos(15*delta)),h+(rIn*sin(15*delta))), (rOut,h+rIn),(rOut,h+rIn+(l2/8)),(rOut,h+rIn+(l2/4)), (rOut,h+rIn+((3*l2)/8)),(rOut,h+rIn+((4*l2)/8)),(rOut,h+rIn+((5*l2)/8)), (rOut,h+rIn+((6*l2)/8)),(rOut,h+rIn+((7*l2)/8)),(rOut,h+rIn+(l2)) ) thetas=linspace(0,2*pi,num=32,endpoint=True) Nozzle=O.materials.append(FrictMat(density=7.8*10**3,young=210e9,poisson=0.3,label="Nozzle")) pts=pack.revolutionSurfaceMeridians([[(pt[0],pt[1]) for pt in poly] for theta in thetas],thetas) surf=pack.sweptPolylines2gtsSurface(pts) nozzle=O.bodies.append(pack.gtsSurface2Facets(surf)) ## # Creating Target ## Aluminum=O.materials.append(FrictMat(density=2.81*10**3,young=71e9,poisson=0.33,label="Aluminum")) p=utils.box((x1,0,0),(x1,x1,z1),fixed=True,material="Aluminum") box=O.bodies.append(p) l=len(O.bodies) t_id=l-1 p1=utils.box((-x1,0,0),(x1,x1,z1),fixed=True,material="Aluminum") box1=O.bodies.append(p1) l_1=len(O.bodies) t_id1=l_1-1 p2=utils.box((2*x1,0,0),(2*x1,x1,z1),fixed=True,material="Aluminum") box2=O.bodies.append(p2) l_2=len(O.bodies) t_id2=l_2-1 p3=utils.box((-2*x1,0,0),(2*x1,x1,z1),fixed=True,material="Aluminum") box3=O.bodies.append(p3) l_3=len(O.bodies) t_id3=l_3-1 p4=utils.box((x1,-2*x1,0),(x1,x1,z1),fixed=True,material="Aluminum") box4=O.bodies.append(p4) l_4=len(O.bodies) t_id4=l_4-1 ## # Creating shots ###
[Yade-users] [Question #680892]: O.engines[].dead=True
New question #680892 on Yade: https://answers.launchpad.net/yade/+question/680892 Hi yade users, I'm using the rotation and translation engines to move bodies and pack of particles. The simulation starts with the engines set to dead= false. After I use Pyrunners in order to stop the movements. For example these are the engines: RotationEngine(angularVelocity=-wt,dead=False,ids=nozzle,rotateAroundZero=True,rotationAxis=(1,0,0),zeroPoint=[0,y_rotation,0]), RotationEngine(angularVelocity=-wt,dead=False,ids=pack_rotation,rotateAroundZero=True,rotationAxis=(1,0,0),zeroPoint=[0,y_rotation,0]), This is the function of the pyrunner: def motion_5(): O.engines[11].dead=True O.engines[12].dead=True And this is the pyrunner: PyRunner(command='motion_5()',iterPeriod=1,firstIterRun=iterrotation_2,nDo=2) I do not understand why the rotations don't stop. Does somebody find the same problem? Thank you so much, Francesco -- 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 #669519]: Drag force
Question #669519 on Yade changed: https://answers.launchpad.net/yade/+question/669519 Status: Answered => Solved Francesco confirmed that the question is solved: Thank you Raphael, I will follow your advices. Best regards, Francesco -- 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 #669519]: Drag force
Question #669519 on Yade changed: https://answers.launchpad.net/yade/+question/669519 Status: Answered => Open Francesco is still having a problem: Firstly thank you for the quick response, I will investigate the potentiality of the function HydroForceEngine. Essentially, I'm not interested to precisely describe the fluid motion. I only want to impose the drag force to all the particles in the 3D tube considering a flow presenting a constant velocity profile. I know it's an approximation but to start I consider it acceptable. So I would like to ask you how much the engine would increase the computational costs of the simulation. I'm working with 6/10 particles inside a 50 mm length tube. Thank you again, Francesco -- 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 #669519]: Drag force
New question #669519 on Yade: https://answers.launchpad.net/yade/+question/669519 Hello, I'm new in Yade and I need to simulate the effect of an air flow inside a tube. I initially modeled it by imposing a given velocity to the particles but i would like to improve my work. I read about the dragEngine that simulates the drag effect on particles that move at a specific velocity. I thought that I should need to specify the air flow velocity in the dragEngine. In this way I would obtain a drag force that essentially would represent my air flow. Unfortunately i observed that this option has not been already developed. Before than apply it manually I would like to ask you if there is any alternative function able to solve my problem. Thank you so much. Best regards, Francesco -- 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 #668594]: Cylindrical loose sphere pack
Question #668594 on Yade changed: https://answers.launchpad.net/yade/+question/668594 Status: Answered => Solved Francesco confirmed that the question is solved: Hello Robert, Thank you, your solution is exactly what I need! I'm a beginner as regards Yade and I didn't understand the potentiality of pack.filterSpherePack. Best regards, Francesco -- 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 #668594]: Cylindrical loose sphere pack
New question #668594 on Yade: https://answers.launchpad.net/yade/+question/668594 Hello, I need to generate a loose pack of spheres with a specific Particle Size Distribution (psd). I identified the .makeCloud function as the most suitable for my need. Unfortunately, I need to fill a cylindrical body and the .makeCloud function allows only to create a spheres pack contained into a box. The other options as the function .randomDensePack combined with pack.inCylinder allow to create a cylindrical spheres pack but do not permit to specify the PSD. Essentially I would need the properties of the function pack.inCylinder for the shape of my pack, but the possibility to specify the PSD guaranteed by the .makeCloud function. So I would like to ask you if it is possible to obtain a cylindrical pack of spheres imposing a specific Particle Size Distribution. Thank you. Best regards, Francesco -- 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
