yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #06169
[Question #201766]: Problems with law-test.py
New question #201766 on Yade:
https://answers.launchpad.net/yade/+question/201766
Hi,
I am not sure how to calculate kn,ks from E, v form frictMat when I use Law2_ScGeom_FrictPhys_CundallStrack(),
In script law-test.py, E=1e7, v=0.3, D1=0.2, D2=0.4, so kn=0.133E7, if ks=v*kn, then ks= 0.0399e7, but when I test
the script with this path:
-----
#LawTester(ids=[0,1],path=[
# (-1e-5,0,0),(-.1,0,0),(-.1,.1,0),(-1e-5,.1,0), # towards, shear, back to intial normal distance
# (-.02,.1,.1),(-.02,-.1,.1),(-.02,-.1,-.1),(-.02,.1,-.1),(-.02,.1,.1), # go in square in the shear plane without changing normal deformation
# (-1e-4,0,0) # back to the origin, but keep some overlap to not delete the interaction
# ],pathSteps=[100],doneHook='tester.dead=True; O.pause()',label='tester',rotWeight=.2,idWeight=.2),
---------
I got Fn/un=0.133E7=kn, but I got Fs/us1=4.44e5 in the third step from (-.1,0,0) to (-.1,.1,0), why?
Another problem is when path go from (-.1,.1,0) to (-1e-5,.1,0), Fs will go down because of max Fs from Fn, but
the point Fs go down is not equal to the Fn*tan(friction angle), why? Thank you!
--
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.