yade-users team mailing list archive

[Roxana Saghafian Larijani <question704035@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Needs information => Open

Roxana Saghafian Larijani gave more information on the question:

Sorry, eta is an attribute that I have added to ViscElCapMat in the source code, so please consider the following script  without eta:

from yade import ymport
from yade import pack
from yade import utils, plot,wrapper
from yade import export


fr = 0.38
rho = 3700
En = 0.22
Et = 0.22
poi=0.3
yoM=200e6
r = 0.0005

Gamma = 0.072
Theta = 0
vB = 6.03 * 1e-9
CapType="Rabinovich"
Tc=0.00045

##defining material

mat=O.materials.append(
        ViscElCapMat(frictionAngle=fr, density=rho, poisson=poi, young= yoM, Vb=vB, gamma=Gamma, theta=Theta, Capillar=True, CapillarType=CapType,en=En, et=Et,tc=Tc)
)
mat2=O.materials.append(
        ViscElCapMat(frictionAngle=fr, density=rho, poisson=poi, young= yoM, Vb=vB, gamma=Gamma, theta=90, Capillar=True, CapillarType=CapType, en=En, et=Et,tc=Tc)
)

#defining the spheres

sp=pack.SpherePack()
sp.makeCloud((-0.035,-0.035,-0.01),(0.035,0.035,0.01),psdSizes=[0.001,0.002,0.003,0.004], psdCumm=[0.0,0.8,0.95,1.0])
sp.toSimulation(material=mat)
Nprtcl=len(O.bodies)
print(Nprtcl)
Tt= utils.PWaveTimeStep()
O.dt = 0.5*Tt
print(O.dt)

#liquidMigration
VV=0.03
Vmin=0.0
for s in O.bodies:
        if not type(s.shape)==wrapper.Sphere:
                continue
        s.state.Vf=VV * (4/3) * 3.14*(s.shape.radius)**3
        s.state.Vmin=Vmin
##

Drum=geom.facetCylinder(material=mat2,center=(0.0,0.0,0.0),radius=0.05,height=0.03,orientation=Quaternion(Vector3(0,0,1),(pi/2.0)))
walls = O.bodies.append(Drum)

##engine
O.engines = [
        ForceResetter(),
        InsertionSortCollider([Bo1_Sphere_Aabb(), Bo1_Facet_Aabb()]),
        InteractionLoop(
                [Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
                [Ip2_ViscElCapMat_ViscElCapMat_ViscElCapPhys()],
                [Law2_ScGeom_ViscElCapPhys_Basic()],

        ),

        NewtonIntegrator(gravity=[0, -9.8, 0]),
        RotationEngine(ids=walls,rotationAxis=[0,0,1],rotateAroundZero=True, zeroPoint=[0.0,0.0,0.0], angularVelocity=2.0),
        LiqControl(label='lqc'),
        VTKRecorder(iterPeriod=1000,recorders=['spheres','facets','colors'],fileName='test-')]

O.run( 500000,True)

O.save('test.bz2')

Regards,
Roxana Saghafian

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