yade-users team mailing list archive

[Ruidong LI <question707278@xxxxxxxxxxxxxxxxxxxxx>]

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.