how to get fabric tensor and micro information?
Hi,
I'm new with Yade and using it for my Ph.D. work in Geotechnics. I have coded a simple triaxial test but don't know how can I get micro-scale information such as fabric tensor, contact normal, branch vectors, interparticle forces, strains, etc.
My supervisor asked me to provide all of those in 2 days and I have no idea what should I do :-(
Would you please tell me how I can get that information? (and also if my code has any issue)
This is my code:
*********
# encoding: utf-8
# the script demonstrates a simple case of triaxial simulation using TriaxialCompres
from yade import pack
sp=pack.
## corners of the initial packing
mn, mx=Vector3(0,0,0), Vector3(10,10,10)
## box between mn and mx, average radius -+ 1/2(20%) , 2k spheres
sp.makeCloud(
## create material #0, which will be used as defult
O.materials.
O.materials.
## Assign section, copy spheres from the packing into the scene
## use defult material, don't care about that for now
O.bodies.
## create walls around the packing
walls=aabbWalls
wallIds=
triax=TriaxialC
wall_bottom_
wall_top_
wall_left_
wall_right_
wall_back_
wall_front_
internalCompac
## define the rest of triaxial parameters here, see in pkg/dem/
sigmaIsoCompac
sigmaLateralCo
max_vel=10,
strainRate=0.01,
label="triax",
)
O.engines=[
ForceResetter(), # reset forces on boundaries in each step
InsertionSortC
InteractionLoop(
[Ig2_Sphere_
[Ip2_FrictMat_
[Law2_
),
GlobalStiffnes
triax,
# you can add TriaxialStateRe
NewtonIntegrat
]
from yade import plot
O.engines=
def history():
plot.addData(
s11=-triax.
s22=-triax.
s33=-triax.
i=O.iter)
plot.plots=
O.saveTmp()
plot.plot()
*********
Thankyou
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 ehsan benabbas for more information if necessary.