yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #15936
Re: [Question #658909]: Question about FEM and DEM coupling
Question #658909 on Yade changed:
https://answers.launchpad.net/yade/+question/658909
Kun Zeng posted a new comment:
Hi Jan,
Sorry to post so many comments and questions above. I changed the
NameError part in my code into this:
rects = (movables + rect((p000,p100,p110,p010)) + rect((p200,p300,p310,p210)) +rect((p000,p010,p011,p001)) + rect((p300,p310,p311,p301)) + rect((p000,p100,p101,p001)) + rect((p100,p200,p201,p101)) + rect((p200,p300,p301,p201)) + rect((p010,p110,p111,p011)) + rect((p110,p210,p211,p111)) + rect((p210,p310,p311,p211)) + rect((p100,p200,p20b,p10b)) + rect((p110,p210,p21b,p11b)) + rect((p100,p110,p11b,p10b)) + rect((p200,p210,p21b,p20b)))
#rects = movables + sum((rect(r) for r in rects),[])
#config.permutations = [[1] + list(x) for x in config.permutations]
O.bodies.append(rects)
(1) This works and could run successfully. However I am not sure if it
is correct way?
(2) And I have a question to assign the material properties and boundary
conditions(like gravity) of the FEM element(I use the salome to generate
and mesh a cylinder). But I don't know how to apply the gravity to this
cylinder. I know the -unv2oofem.py should work for this purpose but how
can I change the example1-unv2oofem.py file to apply the gravity on the
cylinder? (Should I open an another question for this one?)
Thank you very much.
--
You received this question notification because your team yade-users is
an answer contact for Yade.