Re: [Yade-users] [Question #697615]: displacement vs friction graph plotting

2021-06-28 Thread Jan Stránský
Question #697615 on Yade changed:
https://answers.launchpad.net/yade/+question/697615

Status: Open => Answered

Jan Stránský proposed the following answer:
So "friction at interface" is friction force?

In your O.engines, you have no collider, therefore no collision detection and 
no forces.
Have a look at [3].

Cheers
Jan

[3] https://yade-dem.org/doc/user.html#base-engines

-- 
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 #697615]: displacement vs friction graph plotting

2021-06-23 Thread Pawan
Question #697615 on Yade changed:
https://answers.launchpad.net/yade/+question/697615

Status: Needs information => Open

Pawan gave more information on the question:
Sorry for the inconvenience, I am new to this platform.
I have sorted my code till here

from yade import plot
idSteel=O.materials.append(FrictMat(young=210e9,poisson=.25,frictionAngle=.8,label="steel"))
 #material
sp=sphere((0,0,0),1,material=idSteel)  #sphere
plate=box((0,0,-1),(10,100,0),material=idSteel)  #plate
plateID=O.bodies.append(plate) #plateID=0
spID=O.bodies.append(sp) #spID=1
sp.state.blockedDOFs = "xyXYZ" #sphere movement z dir allowed
O.forces.setPermF(1,(0,0,-1))  # 1 mag. force on sphere in -ve z direction
plate.state.blockedDOFs = "xyzXYZ"  #plate movement blocked
plate.state.vel =(0,.02,0)  #plate given const. velocity in y direction

#engine
O.engines=[ForceResetter(),
InteractionLoop(
 [Ig2_Box_Sphere_ScGeom()],
 [Ip2_FrictMat_FrictMat_FrictPhys()],
 [Law2_ScGeom_FrictPhys_CundallStrack()]
),
NewtonIntegrator(),
PyRunner(iterPeriod=500,command="addPlotData()")
]

def addPlotData():
  Dz=plate.state.displ()
  Fz=O.forces.f(plateID)
  plot.addData(i=O.iter,Fz=Fz,Dz=Dz)
plot.plots={'Dz':('Fz')}
plot.plot()
O.dt=0.5*PWaveTimeStep()
O.saveTmp()

Here what I am trying to do is fixing a sphere particle on a moving plate. A 
normal load is applied on the sphere the same thing as we do in the shear 
test ..but here with a single particle. 
Now I want to record the friction at the interface with the displacement of the 
plate.

The problem that I feel is the normal force applied is not getting
transferred to the interface and hence no friction is recorded.

For the friction force, I have taken force on the plateas the only
force on the plate in a moving direction will be friction.

Kindly help me with this .

-- 
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 #697615]: displacement vs friction graph plotting

2021-06-20 Thread Jan Stránský
Question #697615 on Yade changed:
https://answers.launchpad.net/yade/+question/697615

Status: Open => Needs information

Jan Stránský requested more information:
Hello,

> Can you please suggest the code for plotting graph of "Displacement of plate 
> v/s friction at interface" .
> How can we plot "Displacement of plate versus mag. of Friction at interface"

Can you please suggest what is "friction at interface" to be plotted [1]?
Friction coefficient? Tangential force? Some inelastic part of the tangential 
force? ... ?


How to track and plot values, have a look at [1].
"The plot module provides simple interface and storage for tracking various 
data. Although originally conceived for plotting only, it is widely used for 
tracking variables in general."
You can then plot data directly in yade or save the data to be plotted by a 
different program.

Cheers
Jan

[1] https://www.yade-dem.org/wiki/Howtoask
[2] https://yade-dem.org/doc/user.html#tracking-variables

-- 
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 #697615]: displacement vs friction graph plotting

2021-06-20 Thread Pawan
Question #697615 on Yade changed:
https://answers.launchpad.net/yade/+question/697615

Pawan gave more information on the question:
#creating sphere
sp=sphere((0,0,0),1) 
#sp.state.blockedDOFs = "xyzXYZ"

#creating plate
plate=box((0,0,-1.1),(10,10,.1))

#adding sphere and plate
O.bodies.append(sp)
O.bodies.append(plate)

#providing plate velocity
plate.state.vel =(0,.002,0)

#normal force on sphere
sp.setPermF(5,(0,0,-1))

#tribology
mat = FrictMat(frictionAngle=0.5)
O.materials.append(mat)


O.dt=dtSafety*PWaveTimeStep()
print('Timestep',O.dt)

O.engines=[
ForceResetter(),
InteractionLoop(
[Ig2_Box_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]
),
]

We have completed code till here. How can we  plot "Displacement of
plate versus mag. of Friction at interface"

-- 
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 #697615]: displacement vs friction graph plotting

2021-06-18 Thread Pawan
New question #697615 on Yade:
https://answers.launchpad.net/yade/+question/697615

We have written a code for a sphere on a horizontal plate with some contact law 
between them and a normal load is applied on the sphere. 

sp=sphere((0,0,0),1)
plate=box((0,0,-1.1),(10,10,.1))
O.bodies.append(sp)
O.bodies.append(plate)
plate.state.vel =(0,.002,0)
sp.setPermF(5,(0,0,-1))
mat = FrictMat(frictionAngle=0.5)
O.materials.append(mat)
InteractionLoop(
   [Ig2_Box_Sphere_ScGeom()],
   [Ip2_FrictMat_FrictMat_FrictPhys()],
   [Law2_ScGeom_FrictPhys_CundallStrack()]
)

Can you please suggest the code for plotting graph of "Displacement of plate 
v/s friction at interface" . 
Thanks in Advance

-- 
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