How to do model bearing capacity of soil
I need to make model bearing capacity of soil .By wanting to know the weight that can be received . What to do model ?
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Jan Stránský
- Solved:
- Last query:
- Last reply:
Revision history for this message
|
#1 |
You could use Terzaghi's theory for bearing capacity :-)
[1] https:/
Revision history for this message
|
#3 |
Yes . I use terzaghi’s theory .Can. you building model bearing ?
Revision history for this message
|
#4 |
I am sorry, I do not understand what you are asking. Can you please rephrase your question?
Revision history for this message
|
#5 |
I want to make a model bearing capacity of soil . But I can’t resize the pressure wall. I want to resize the pressed wall, what to do?
This is my code#######
from yade import pack
sp=pack.
## corners of the initial packing
mn,mx=Vector3(
## box size, radius, number of ball
sp.makeCloud(
## create material #0, which will be used as default
O.materials.
O.materials.
## copy spheres from the packing into the scene
## use default material, don't care about that for now
O.bodies.
## create walls around the packing
walls=aabbWalls
wallIds=
# wall of each side
triax=TriaxialC
wall_bottom_
wall_top_
wall_left_
wall_right_
wall_back_
wall_front_
internalCompac
sigmaIsoCompac
sigmaLateralCo
max_vel=10,
strainRate=0.01,
label="triax"
)
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
GlobalStiffnes
triax,
NewtonIntegrat
]
from yade import plot
O.engines=
def history():
plot.addData(
s11=-
s22=-
s33=-
i=O.iter)
plot.plots={'i': ('e11',
O.saveTmp()
plot.plot()
Revision history for this message
|
#6 |
Revision history for this message
|
#7 |
>But I can’t resize the pressure wall. I want to resize the pressed wall, what to do?
What is the pressure/pressed wall? It would help to understand if you mention the methods you have already tried. (e.g. I tried to change _insert command here_ but I see that it doesn't change the size).
I think maybe you want to use facets [1] to model the application of weight, but I am still a bit unsure what you are asking.
Revision history for this message
|
#8 |
I want to make a model of soil weight bearing. Which the pressure that I use is needed as a square of various sizes.
but my code is not scalable. What should I do?
Revision history for this message
|
#9 |
Hello,
> What should I do?
help us to help you :-)
If you expect a help from us, you should try to formulate your problem both shortly but at the same time as specifically as possible. Also read [2].
It may happen that the formulation is not complete or too vague (like your description here). Then you should provide required information. E.g. answer:
- What is the pressure/pressed wall? (Robert's question)
- it would help to understand if you mention the methods you have already tried (Robert's info request)
- have you tried [1]? if yes, why is did not help?
Some other notes:
> I want to make a model of soil weight bearing.
yes, this is clear form the original post, you do not need to repeat this, rather focus on the information that are missing.
> But I can’t resize the pressure wall.
Yade wall is an infinite plate, is cannot be resized. Use a different type of a boundary (as already proposed by Robert), e.g. facet or box.
cheers
Jan
Revision history for this message
|
#10 |
ok sir, so i saw your video that you make" rod pressure soil" .I really want to make a pressure that you make . i try to use that code for my model but it is not working.i can not understand a code that you make and trying to run a code. it does not work too.can you please help me to make a pressure for my model. this is my model
from yade import pack, plot
# create rectangular box from facets
O.bodies.
# create empty sphere packing
# sphere packing is not equivalent to particles in simulation, it contains only the pure geometry
O.materials.
# sphere packing is not equivalent to particles in simulation, it contains only the pure geometry
sp=pack.
# generate randomly spheres with uniform radius distribution
O.materials.
sp.makeCloud(
# add the sphere pack to the simulation
sp.toSimulation()
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
# handle sphere+sphere and facet+sphere collisions
[Ig2_
[Ip2_
[Law2_
),
NewtonIntegrat
# call the checkUnbalanced function (defined below) every 2 seconds
]
O.dt=.5*
#
O.saveTmp()
Revision history for this message
|
#11 |
ok sir, so i saw your video that you make" rod pressure soil" .I really want to make a pressure that you make . i try to use that code for my model but it is not working.i can not understand a code that you make and trying to run a code. it does not work too.can you please help me to make a pressure for my model. this is my model
from yade import pack, plot
# create rectangular box from facets
O.bodies.
# create empty sphere packing
# sphere packing is not equivalent to particles in simulation, it contains only the pure geometry
O.materials.
# sphere packing is not equivalent to particles in simulation, it contains only the pure geometry
sp=pack.
# generate randomly spheres with uniform radius distribution
O.materials.
sp.makeCloud(
# add the sphere pack to the simulation
sp.toSimulation()
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
# handle sphere+sphere and facet+sphere collisions
[Ig2_
[Ip2_
[Law2_
),
NewtonIntegrat
# call the checkUnbalanced function (defined below) every 2 seconds
]
O.dt=.5*
#
O.saveTmp()
Revision history for this message
|
#12 |
if you don't understand that the meaning that i try to say please tell me . my English is not good . thank you
Revision history for this message
|
#13 |
> i saw your video that you make" rod pressure soil"
please provide a link
> but it is not working.
> it does not work too
please be more specific. If it gives errors, provide the errors. If the results are not as expected, describe actual and expected results....
> i can not understand a code that you make
please provide the code (if you want help to understand it)
cheers
Jan
Revision history for this message
|
#14 |
Dear jan,
this's link
https:/
i want to put the compressin like this video for my script :
from yade import pack, plot
# create rectangular box from facets
O.bodies.
# create empty sphere packing
# sphere packing is not equivalent to particles in simulation, it contains only the pure geometry
O.materials.
# sphere packing is not equivalent to particles in simulation, it contains only the pure geometry
sp=pack.
# generate randomly spheres with uniform radius distribution
O.materials.
sp.makeCloud(
# add the sphere pack to the simulation
sp.toSimulation()
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
# handle sphere+sphere and facet+sphere collisions
[Ig2_
[Ip2_
[Law2_
),
NewtonIntegrat
# call the checkUnbalanced function (defined below) every 2 seconds
]
O.dt=.5*
#
O.saveTmp()
thank for helping me
Revision history for this message
|
#15 |
Some notes:
- define material prior to creating bodies
- use box [1] for simple geometries instead of facets
a MWE:
###
from yade import pack
O.materials.
# create rectangular box from boxes (maybe there is some library function, but I did not find it quickly)
cx,cy,cz = .5,.5,.5 # center of the box
dx,dy,dz = .5,.5,.5 # half-dimensions of the box
t = 0.025 # half-thickness of walls
topx,topy = 0.2, 0.3 # top box half-dimensions
left = box((cx-
right = box((cx+
front = box((cx,
back = box((cx,
bottom = box((cx,
top = box((cx,
O.bodies.
sp=pack.
# generate randomly spheres with uniform radius distribution
sp.makeCloud(
# add the sphere pack to the simulation
sp.toSimulation()
O.engines=[
ForceResetter(),
InsertionSor
InteractionLoop(
[
[
[
),
NewtonIntegr
]
O.dt=.5*
top.state.vel = (0,0,-.5) # setting velocity is the simplest boundary condition, but you can also apply force or whatever...
###
cheers
Jan
Revision history for this message
|
#16 |
this is what i need thank you so much.
But i got a little bit problem. Some particles at the mid-bottom and a box of facet are got through the box. Is there anyway i can fix it?
here is video:
https:/
Revision history for this message
|
#17 |
You can:
- stop the simulation before it occurs
- use somewhat sophisticated simulation control (then just prescribe velocity and let it run)
cheers
Jan
Revision history for this message
|
#19 |
Thanks Jan Stránský, that solved my question.
Revision history for this message
|
#20 |
Dear jan ,
I appreciate the things you do so much indeed. So now the top box drop down with velocity that set up, I would like to know that can the top box drop down by force.and How to change it ?
Revision history for this message
|
#21 |
O.forces.
or for older Yade versions
O.forces.
cheers
Jan
[1] https:/