Changing friction angle in triaxial test
Hi guys!
I am trying to do parametric study for triaxial test. I have played around the friction angle and alphaKr. However changing of these two parameters have no effect on my results. I am not sure where is my mistake . Would you please advice me.
Thanks heaps,
Here is my code:
from yade import pack
nRead=utils.
)
from yade.params import table
num_spheres=
targetPorosity = 0.8 #the porosity we want for the packing (3 specimens: (Ei,n) = (1,0.382), (2,0.387), (3,0.409) )
compFricDegree = table.compFricD
finalFricDegree = 30 # contact friction during the deviatoric loading
rate=0.001 # loading rate (strain rate)
damp=0.3 # damping coefficient
stabilityThresh
key='_kozoooooo
young=30e9 # contact stiffness k_n/Ds
mn,mx=Vector3(
thick = 0.01 # thickness of the plates
## create materials for spheres and plates
O.materials.
O.materials.
## create walls around the packing
walls=utils.
wallIds=
## use a SpherePack object to generate a random loose particles packing
sp=pack.
psdSizes=
#psdCumm=
psdCumm=
sp.makeCloud(
sp.toSimulation
triax=TriaxialS
thickness = thick,
stressMask = 7,
#the value of confining stress for the intitial (growth) phase
newton=
O.engines=[
useIncrement
always_
label=
),
GlobalStiffness
triax,
newton
]
### APPLYING CONFINING PRESSURE ###
#######
while 1:
O.run(1000, True)
#the global unbalanced force on dynamic bodies, thus excluding boundaries, which are not at equilibrium
unb=unbalance
#average stress
#note: triax.stress(k) returns a stress vector, so we need to keep only the normal component
meanS=
print 'unbalanced force:',unb,' mean stress: ',meanS
print 'void ratio='
print 'mean stress of engine', triax.meanStress
print '0.12,0.
if triax.porosity<0.7:
break
#######
### REACHING A SPECIFIED POROSITY PRECISELY ###
import sys #this is only for the flush() below
while triax.porosity>
# # we decrease friction value and apply it to all the bodies and contacts
print "\r Friction: ",compFricDegree," porosity:
#######
### DEVIATORIC LOADING ###
#######
#We move to deviatoric loading, let us turn internal compaction off to keep particles sizes constant
triax.internalC
# Change contact friction (remember that decreasing it would generate instantaneous instabilities)
#setContactFric
setContactFrict
#set stress control on x and z, we will impose strain rate on y (5)
triax.stressMask = 5
#now goal2 is the target strain rate
triax.goal2=-rate
#we assign constant stress to the other directions
triax.goal1=
triax.goal3=
##we can change damping here. What is the effect in your opinion?
#newton.damping=0.1
##Save temporary state in live memory. This state will be reloaded from the interface with the "reload" button.
O.saveTmp()
while triax.strain[1] < 0.25:
O.run(50); O.wait()
Question information
- Language:
- English Edit question
- Status:
- Answered
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Seti for more information if necessary.