← Back to team overview

yade-users team mailing list archive

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