Question #707278 on Yade changed:
https://answers.launchpad.net/yade/+question/707278
Status: Answered => Open
Ruidong LI is still having a problem:
Hi! Jan
Thanks for your reply. Your answer is rigorous and helpful.
> area = 2* pi*r_cyl*h_cyl/len(nodes) # cylinder surface / number of
nodes ? The question mark in my comment is meant seriously. For regular
grid it is probably OK. Also consider that the "boundary" nodes (bottom
and top) maybe should have half force then "inner" nodes?
Yes. That's a really critical problem that needed to be clarified. What
confused me currently is that I found the number of nodes is not consistent
with what I expected. For example, in the beginning, we defined the number of
facets by the following code [1]:
nw=40 # number of facets along the perimeter
nh=25 # number of facets along the height
nodes = [b for b in O.bodies if isinstance(b.shape, GridNode)]
print(len(nodes))
If we print the number of nodes, we will get 8000 where my initial
suppose is 40*26 = 1040. That indicates that the ID of nodes forming
Pfacet are independent of each other. In other words, more than one node
can be created at the same position. So, I am wondering if maybe it is
correct to use ' confining Stress * cylinder surface / number of nodes'
to represent the average force applied on each node? We don't need to
classify the top and bottom nodes anymore as the force applied on nodes
is the same. It is identical to obtaining the node force by calculating
the force applied on each Pfacet and distributing it into three nodes.
[1] https://answers.launchpad.net/yade/+question/681410
--
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
