Consolidation method in Direct Shear Test test
Dear All,
I'm in a lot of trouble
I am doing an analysis on a direct(simple) shear test.
I have created a program like below.
So I have one question. It's about consolidation.
I was able to generate shear , but not consolidation.
So is there a way to make it consolidation?
My idea was to apply a constant stress/
However, just as I can create a plate or facet, I cannot push a plate or facet, with a specified force or stress.
Can you give me a program to apply a constant stress/force to the tabletop?
Thank you very much.
#######
from __future__ import print_function
from yade import pack,plot,
from yade import export,qt
from yade.gridpfacet import *
import math
import numpy as np
import os
# load parameters from file if run in batch
# default values are used if not run from batch
# the "if 0:" block will be never executed, therefore the "else:" block will be
# to use cloud instead of regular packing, change to "if 1:" or something similar
frictangle = 0.436
density = 2400.0
young = 3e7
gravity = (0.0, 0.0, -9.81)
velocity_cyl=0.10
# create cloud of spheres and insert them into the simulation
# we give corners, mean radius, radius variation
#######
#Create sphere
#######
mat_sp = CohFrictMat(young = young, poisson = 0.35,frictionAn
O.materials.
pred=pack.
sp0=pack.
O.bodies.
#######
#Boundary creation
#######
idsCyl1 = O.bodies.
idsCyl2 = O.bodies.
#idsCyl3 = O.bodies.
radius1 = .06
radius2 = .15
fixBoxIds=[]
moveBoxIds=[]
i=0
for i in range(0,365,5):
r = math.radians(i)
x1 = radius1 * math.cos(r)
y1 = radius1 * math.sin(r)
x2 = radius2 * math.cos(
y2 = radius2 * math.sin(
x3 = radius2 * math.cos(
y3 = radius2 * math.sin(
x4 = radius1 * math.cos(
y4 = radius1 * math.sin(
x5 = radius1 * math.cos(
y5 = radius1 * math.sin(
x6 = radius2 * math.cos(r)
y6 = radius2 * math.sin(r)
fixBoxIds.
fixBoxIds.
moveBoxIds.
moveBoxIds.
#######
#engines
#######
O.engines = [
# interaction loop
),
]
#######
#output
#######
sp_num = len([b for b in O.bodies if isinstance(b.shape, Sphere)]) # 粒子の総数をカウント
print("number of spheres = ",sp_num)
O.dt=0.
def plotPlotData():
pipe_force_x = sum([O.
pipe_force_y = sum([O.
plot.addData(
plot.saveDataT
plot.plots=
plot.plot()
#######
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- 内山康太郎
- Solved:
- Last query:
- Last reply: