yade-users team mailing list archive

[zengzhihong <question256715@xxxxxxxxxxxxxxxxxxxxx>]

New question #256715 on Yade:
https://answers.launchpad.net/yade/+question/256715

Hi all,
I'm a green hand of Yade users.I'm trying to build a mode of shearing of the concrete.
I consided to solve this basic on two simple examples which the one is the simpleshear and the other is uniax of concrete.
I mixed those two source, but still can not get the stimulate of the shearing concrete test
Here is part of the source:

from yade import plot
from yade.pack import *


O=Omega() 

length=0.1
height=0.02
width=0.04
thickness=0.001


dtSafety=.8
sphereRadius=3.5e-3
intRadius=1.5
strainRateTension=.05
strainRateCompression=.5


concreteId=O.materials.append(CpmMat(young=24e9,frictionAngle=atan(0.8),poisson=.2,density=4800,sigmaT=3.5e6,relDuctility=30,epsCrackOnset=1e-4,isoPrestress=0))

pred=inAlignedBox((0,0,-0.02),(0.01,0.02,0.02))

leftBox = box( center=(-thickness/2.0,(height)/2.0,0), extents=(thickness/2.0,5*(height/2.0+thickness),width/2.0) ,fixed=True,wire=True)

lowBox = box( center=(length/2.0,-thickness/2.0,0), extents=(length/2.0,thickness/2.0,width/2.0) ,fixed=True,wire=True)

rightBox = box( center=(length+thickness/2.0,height/2.0,0), extents=(thickness/2.0,5*(height/2.0+thickness),width/2.0) ,fixed=True,wire=True)

upBox = box( center=(length/2.0,height+thickness/2.0,0), extents=(length/2.0,thickness/2.0,width/2.0) ,fixed=True,wire=True)

behindBox = box( center=(length/2.0,height/2.0,-width/2.0-thickness/2.0), extents=(2.5*length/2.0,height/2.0+thickness,thickness/2.0), fixed=True,wire=True)

inFrontBox = box( center=(length/2.0,height/2.0,width/2.0+thickness/2.0), extents=(2.5*length/2.0,height/2.0+thickness,thickness/2.0), fixed=True,wire=True)

O.bodies.append([leftBox,lowBox,rightBox,upBox,behindBox,inFrontBox])

spheres=pack.randomDensePack(pack.inAlignedBox((0,0,-0.02),(0.1,0.02,0.02)),spheresInCell=2000,radius=3.5e-3,material=concreteId)
O.bodies.append(spheres)


O.engines=[
	ForceResetter
	InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=intRadius,label='is2aabb'),],verletDist=.05*sphereRadius),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=intRadius)],
		[Ip2_CpmMat_CpmMat_CpmPhys()],
		[Law2_ScGeom_CpmPhys_Cpm()],
	),
	NewtonIntegrator(damping=.4)

	]
O.dt=dtSafety*PWaveTimeStep()

O.engines = O.engines+[KinemCTDEngine(compSpeed=0.5,sigma_save=(),temoin_save=(),targetSigma=40000.0,LOG=False)]

ss2sc.interactionDetectionFactor=1.
is2aabb.aabbEnlargeFactor=-1.




First ,I remember that it would set UniaxialStrainer if we did a uniax test of concrete,but I don't know how to set the strain of the concrete when I wanna to do a shear test.

Second, I wonder if the KinemCTDEngine, KinemCNDEngine which make the strain by controlling walls can be work in the concrete material...?If they could ,what else parameters was needed? Otherwise, how can I put a shear-strain to the concrete material..?

Thanks for any suggestions,



-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.