Question #688715 on Yade changed:
https://answers.launchpad.net/yade/+question/688715
Status: Answered => Open
mrhappy is still having a problem:
Hey,
Thank you for your reply. I followed [1] but it does not work. Cna you
check my code?
[1] https://yade-dem.org/doc/user.html#individual-interactions-on-demand
Thank you soo much for your help!
# MATERIAL PROPERTIES
intR=20 # allows near neighbour interaction (can be adjusted for every packing)
DENS=2500 # Density
YOUNG=1800
FRICT=7
ALPHA=0.1
TENS=1000000000000000
COH=1000000000000000000
iterMax = 1
O.materials.append(JCFpmMat(type=1,density=DENS,young=YOUNG,poisson=ALPHA,frictionAngle=radians(FRICT),tensileStrength=TENS,cohesion=COH,label='mat1'))
# PARTICLES
O.bodies.append([
# 0
sphere(center=(0,0,0),radius=0.5,fixed=True,material='mat1'),
# 1
sphere((0,1,0),radius=0.5,material='mat1')
])
nbSpheres = 2;
# BOUNDARY CONDITION
O.bodies[1].state.blockedDOFs = "xyzXYZ"
O.bodies[1].state.vel = (0,0.001,0)
O.bodies[1].state.angVel = (0,0,0)
# FUNCTIONAL COMPONENTS
# 1
def addPlotData1():
plot.addData(t=O.time,DX1 = O.bodies[1].state.pos[0],DY1 =
O.bodies[1].state.pos[1],DZ1 = O.bodies[1].state.pos[2], FX1 =
O.forces.f(1)[0], FY1 = O.forces.f(1)[1], FZ1 = O.forces.f(1)[2])
# 2
def addPlotData2():
plot.saveDataTxt('/home/nabid/Desktop/SCRIPTS/dataFile.txt',vars=('t','DX1','DY1','DZ1','FX1','FY1','FZ1'))
#PRINTING
from yade import plot
from pprint import pprint
i=createInteraction(0,1)
i.geom, i.phys
# SIMULATION LOOP
O.engines=[
PyRunner(command='addPlotData1()',iterPeriod=1),
PyRunner(command='addPlotData2()',iterPeriod=1),
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[ ]
),
GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=10,timestepSafetyCoefficient=0.5,
defaultDt=utils.PWaveTimeStep()),
NewtonIntegrator(damping=0.1),
]
#TIME STEP
O.dt=0.5e-4*PWaveTimeStep()
#PLOTTING
plot.plots={'DY1':('FX1','FY1')}
plot.plot()
O.run(iterMax)
# SAVE SIMULATION
--
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 : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
