Question #707253 on Yade changed:
https://answers.launchpad.net/yade/+question/707253
Yipeng Xie posted a new comment:
Adding to the previous comment, I ran two sets of comparison tests. I
monitored the i.phys.moment_bending and the i.phys.normalForce at the
end of compaction and the end of shearing, and I found that the contacts
information is the same for both tests at the end of compaction, but
quite different at the end of shearing. From the output results, both
tests activate the rolling-resistance moment before the shearing and
assign the correct phys.maxRollPl to every existing and potential
contact. So I'm wondering what is the reason for the difference?
(Both for packing composed of 50 particles with radius = 0.1 and
CohesiveMat)
######### TEST No.1 #############
##set initial etaRoll=0.0;
##set 'always_use_moment_law' = True in compaction period;
##search for each sphere and assign the target etaRoll=0.3, and search for each
sphere and assign the target maxRollPl=0.03 before shearing start
from __future__ import print_function
from yade import pack,qt,plot
import matplotlib; matplotlib.rc('axes',grid=True)
import pylab
O.materials.append(CohFrictMat(young=1.e9,poisson=.8,alphaKr=1.0,frictionAngle=.0,momentRotationLaw=True,etaRoll=0.0,label='balls'))
sp = pack.SpherePack()
dmax=1.0
size=2*dmax
thickness=size/3.0
sp.makeCloud(minCorner=(0,0,thickness/2),maxCorner=(size,size,thickness/2),rMean=0.1,rRelFuzz=.0,num=50,periodic=True,seed=1)
sp.toSimulation()
O.cell.hSize = Matrix3(size,0,0, 0,size,0, 0,0,thickness) # RVE scale up
print(len(O.bodies))
for p in O.bodies:
p.state.blockedDOFs = 'zXY'
p.state.mass = 2650 * thickness * pi * p.shape.radius**2 # 0.1 = thickness
of cylindrical particle # RVE scale up
inertia = 0.5 * p.state.mass * p.shape.radius**2
p.state.inertia = (.5*inertia,.5*inertia,inertia)
O.dt = utils.PWaveTimeStep()
print(O.dt)
O.engines = [
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom6D()],
[Ip2_CohFrictMat_CohFrictMat_CohFrictPhys()],
[Law2_ScGeom6D_CohFrictPhys_CohesionMoment(always_use_moment_law=True,
useIncrementalForm=True)]
),
PeriTriaxController(
dynCell=True,
goal=(-4.e4,-4.e4,0),
stressMask=3,
relStressTol=.001,
maxUnbalanced=.001,
maxStrainRate=(.5,.5,.0),
doneHook='compactionFinished()',
label='biax'
),
NewtonIntegrator(damping=.1),
]
phase=0
def compactionFinished():
global phase
if phase==0:
#O.engines = O.engines[:3]+O.engines[4:]
for i in O.interactions:
print(i.phys.maxRollPl,i.phys.momentRotationLaw,i.phys.moment_bending,i.phys.normalForce)
O.materials[0].etaRoll=0.3
for i in O.interactions:
i.phys.maxRollPl=0.03
phase+=1
elif phase==1:
setContactFriction(radians(33))
# set the current cell configuration to be the reference one
O.cell.trsf=Matrix3.Identity
# change control type: keep constant confinement in x,y, 20%
compression in z
biax.goal=(-4.e4,-0.05,0)
biax.stressMask=1
# allow faster deformation along x,y to better maintain
stresses # RVE scale up, y value should be scale down a lot
biax.maxStrainRate=(1.0,0.05,0.0)
# next time, call triaxFinished instead of compactionFinished
biax.doneHook='biaxFinished()'
# do not wait for stabilization before calling triaxFinished
biax.maxUnbalanced=10
def biaxFinished():
for i in O.interactions:
print(i.phys.maxRollPl,i.phys.momentRotationLaw,i.phys.moment_bending,i.phys.normalForce)
print('Finished')
O.pause()
###print after compaction
0.0 True Vector3(0,0,0) Vector3(-4356.022418727271543,4190.262910390562865,0)
0.0 True Vector3(0,0,-0) Vector3(-946.2413718796913145,-3824.453444819895594,0)
0.0 True Vector3(0,0,-0) Vector3(2005.3174721931573,-1929.808863517112059,0)
0.0 True Vector3(0,0,0) Vector3(3190.850202655290104,-3070.077577822232342,0)
0.0 True Vector3(0,0,0) Vector3(2204.808551871029977,8902.688524042694553,0)
0.0 True Vector3(0,0,0) Vector3(4554.374440800664161,-4381.226637213368122,0)
0.0 True Vector3(0,0,-0) Vector3(1137.519578345387799,4598.20941843654964,0)
0.0 True Vector3(0,0,0) Vector3(-3344.408293120146936,-965.5064723078584166,0)
0.0 True Vector3(0,0,-0) Vector3(3144.948266091866572,907.7799091766298716,0)
0.0 True Vector3(0,0,0) Vector3(2247.276737117122593,9062.611694681330846,0)
0.0 True Vector3(0,0,0) Vector3(593.3228751609747178,2397.947073331891715,0)
0.0 True Vector3(0,0,-0) Vector3(-175.8118524235571272,-709.8813002115898598,0)
0.0 True Vector3(0,0,-0) Vector3(3130.931207949010513,-3013.271274804336372,0)
0.0 True Vector3(0,0,-0) Vector3(588.0179022407435241,2376.474661917997764,0)
0.0 True Vector3(0,0,-0) Vector3(1529.626901015044496,6182.86709493719809,0)
0.0 True Vector3(0,0,-0) Vector3(-4141.642333725674689,-1195.060840778490729,0)
0.0 True Vector3(0,0,-0) Vector3(-1792.775700468065679,1725.117708369524735,0)
0.0 True Vector3(0,0,0) Vector3(-67.53360237268792332,-271.5920647003094928,0)
0.0 True Vector3(0,0,-0) Vector3(4135.21852000556737,1194.014624050691964,0)
0.0 True Vector3(0,0,0) Vector3(-3810.761293476583887,-1099.517895306677246,0)
0.0 True Vector3(0,0,-0) Vector3(-2025.926394463087945,1949.234293640667602,0)
0.0 True Vector3(0,0,0) Vector3(-2675.563963173430238,2574.376774606120762,0)
0.0 True Vector3(0,0,-0) Vector3(5272.416890895456163,-5075.730797074852489,0)
0.0 True Vector3(0,0,0) Vector3(-85.79522611963324152,-216.9439744441962148,0)
0.0 True Vector3(0,0,-0) Vector3(-500.4888775940822256,480.3831482944814866,0)
0.0 True Vector3(0,0,0) Vector3(6790.002219449754193,1956.332687601257248,0)
0.0 True Vector3(0,0,0) Vector3(-4785.688251478742131,-1381.727736622080783,0)
0.0 True Vector3(0,0,0) Vector3(581.3424312033406522,2350.305985913877976,0)
0.0 True Vector3(0,0,-0) Vector3(-1741.833517726845912,1676.134517511269678,0)
0.0 True Vector3(0,0,0) Vector3(-5647.573462837440275,-1523.645580069464813,0)
0.0 True Vector3(0,0,-0) Vector3(-195.9442612484697008,-660.6786580392528094,0)
0.0 True Vector3(0,0,-0) Vector3(-482.6886363669209459,425.6844903374354203,0)
0.0 True Vector3(0,0,0) Vector3(5403.182807112851151,-5193.547357533195282,0)
0.0 True Vector3(0,0,0) Vector3(-4037.211573523102743,3885.827380731004723,0)
0.0 True Vector3(0,0,-0) Vector3(-5782.157319227181688,-1672.553876581613167,0)
0.0 True Vector3(0,0,-0) Vector3(-511.7428450696506275,537.7315348333218026,0)
0.0 True Vector3(0,0,-0) Vector3(-412.5180930583922532,-1667.067608065389095,0)
0.0 True Vector3(0,0,0) Vector3(-440.0569357205779397,462.5879847317094118,0)
0.0 True Vector3(0,0,-0) Vector3(-79.40905557252953884,-321.1603524112854302,0)
0.0 True Vector3(0,0,-0) Vector3(574.9178308474665755,-507.406089794769116,0)
0.0 True Vector3(0,0,-0) Vector3(6757.742281344556432,1950.531419672886159,0)
0.0 True Vector3(0,0,0) Vector3(3169.105186400296589,-3049.983867139284484,0)
0.0 True Vector3(0,0,0) Vector3(4934.170380698124063,1424.503266712679988,0)
0.0 True Vector3(0,0,0) Vector3(-3206.510905196595559,-925.5671278104202884,0)
0.0 True Vector3(0,0,0) Vector3(2278.564969719079272,9189.057761753021623,0)
0.0 True Vector3(0,0,-0) Vector3(-244.3011851037109068,-987.693651036061965,0)
0.0 True Vector3(0,0,0) Vector3(1613.0290283192885,-1551.795665594201864,0)
0.0 True Vector3(0,0,-0) Vector3(-2624.200635546494141,2525.699160412983474,0)
0.0 True Vector3(0,0,0) Vector3(-1894.742558239544678,-7654.642720827022458,0)
0.0 True Vector3(0,0,0) Vector3(-5453.070612922661894,-1574.730209720701396,0)
0.0 True Vector3(0,0,0) Vector3(550.7112834928792608,-530.2544805781109289,0)
0.0 True Vector3(0,0,0) Vector3(-6360.602085101942976,-1836.321129405905594,0)
0.0 True Vector3(0,0,0) Vector3(1278.409820090018457,5166.437438540249786,0)
0.0 True Vector3(0,0,0) Vector3(-501.3525301434042376,-2026.262393923267155,0)
0.0 True Vector3(0,0,-0) Vector3(-2271.811739081067572,-9191.198042164664002,0)
0.0 True Vector3(0,0,-0) Vector3(2118.342531565449917,-2038.576273405231632,0)
0.0 True Vector3(0,0,0) Vector3(6493.214469140548317,1876.119468843037794,0)
0.0 True Vector3(0,0,0) Vector3(-5400.851266028918872,5200.117448928389422,0)
0.0 True Vector3(0,0,0) Vector3(-2649.949579478481155,2549.48043116090048,0)
0.0 True Vector3(0,0,-0) Vector3(-5273.094760494164802,5068.296797843814602,0)
0.0 True Vector3(0,0,0) Vector3(642.908085010988998,2599.278453244339289,0)
0.0 True Vector3(0,0,0) Vector3(-5498.449301278576968,-1587.832342098636673,0)
0.0 True Vector3(0,0,-0) Vector3(-989.4557685983584179,-3999.091185152995877,0)
0.0 True Vector3(0,0,0) Vector3(-1629.28978454572939,-6586.428084369147655,0)
0.0 True Vector3(0,0,-0) Vector3(1019.938206925837108,4123.007007379715105,0)
0.0 True Vector3(0,0,0) Vector3(-4355.703135797158211,4189.954115430237835,0)
0.0 True Vector3(0,0,0) Vector3(2674.353110719725009,-2573.193619222766301,0)
0.0 True Vector3(0,0,-0) Vector3(1257.53382450168715,5082.085015788150486,0)
0.0 True Vector3(0,0,0) Vector3(5671.070048103227236,1741.359121719895256,0)
0.0 True Vector3(0,0,-0) Vector3(131.4846787252013769,662.5274659514100222,0)
0.0 True Vector3(0,0,-0) Vector3(-3257.224446709202766,3134.789185079573599,0)
0.0 True Vector3(0,0,-0) Vector3(-5416.749901594610492,-1564.545644680121768,0)
0.0 True Vector3(0,0,0) Vector3(-2204.759659924044172,-8902.875961359055509,0)
0.0 True Vector3(0,0,-0) Vector3(178.7324749728237521,718.8491860076645708,0)
0.0 True Vector3(0,0,0) Vector3(-584.7473032754612632,561.2564641809614159,0)
0.0 True Vector3(0,0,0) Vector3(-6.158419127361661083,26.0734180612055404,0)
0.0 True Vector3(0,0,0) Vector3(582.8097420617934858,-560.4734523253433736,0)
0.0 True Vector3(0,0,-0) Vector3(63.03895382092825628,-5.235699592228107768,0)
0.0 True Vector3(0,0,0) Vector3(-1398.024721625517486,-403.3393896789778523,0)
0.0 True Vector3(0,0,0) Vector3(4389.74447505059652,1266.925791252061799,0)
0.0 True Vector3(0,0,-0) Vector3(8.508941415349747572,34.414349564750502,0)
0.0 True Vector3(0,0,-0) Vector3(495.3726959074963929,-476.3910924743690316,0)
0.0 True Vector3(0,0,0) Vector3(-2971.403754700405898,2859.808144833682491,0)
0.0 True Vector3(0,0,0) Vector3(-547.5349918831619789,527.0987733962641641,0)
0.0 True Vector3(0,0,0) Vector3(-524.0691227176483835,504.4724037886144288,0)
0.0 True Vector3(0,0,-0) Vector3(24.54646675697048508,210.3727536274160457,0)
0.0 True Vector3(0,0,-0) Vector3(-1214.164619003100142,-4906.390546028451354,0)
0.0 True Vector3(0,0,0) Vector3(-260.111874640118117,-1051.591946580115064,0)
0.0 True Vector3(0,0,0) Vector3(2302.619237149820037,-2214.882467420897228,0)
0.0 True Vector3(0,0,0) Vector3(-4371.180785150242627,4204.733916816947385,0)
0.0 True Vector3(0,0,0) Vector3(-1098.412608630502746,-4439.050616960828847,0)
0.0 True Vector3(0,0,0) Vector3(6790.669141671738544,1963.812786986948595,0)
0.0 True Vector3(0,0,0) Vector3(2484.756067007789625,717.591842242276698,0)
0.0 True Vector3(0,0,-0) Vector3(-6143.255952945281024,-1773.664794304884481,0)
0.0 True Vector3(0,0,-0) Vector3(-2514.251945581449036,2418.579763056834508,0)
0.0 True Vector3(0,0,0) Vector3(-2548.34829534631217,2452.706293236090914,0)
0.0 True Vector3(0,0,-0) Vector3(4094.651600411214986,1181.542155068722195,0)
0.0 True Vector3(0,0,0) Vector3(-1968.954938907677843,-7954.994515831598619,0)
0.0 True Vector3(0,0,-0) Vector3(-6472.564400726898384,-1870.002264900543423,0)
0.0 True Vector3(0,0,0) Vector3(1283.919363400329985,5189.285403755833613,0)
0.0 True Vector3(0,0,0) Vector3(-1734.005746328204168,-7010.486056481552623,0)
0.0 True Vector3(0,0,0) Vector3(2673.178560915678645,771.9995470702701823,0)
0.0 True Vector3(0,0,-0) Vector3(-3349.020566802118537,-966.8133839347973435,0)
0.0 True Vector3(0,0,-0) Vector3(48.1053186181306387,40.59859085302196746,0)
0.0 True Vector3(0,0,0) Vector3(-2721.179722815830246,-785.4536521531675817,0)
0.0 True Vector3(0,0,-0) Vector3(-2241.433496616614775,-9067.109488036401672,0)
0.0 True Vector3(0,0,0) Vector3(6134.077892066618915,1769.510159572369503,0)
0.0 True Vector3(0,0,0) Vector3(1097.108672801957937,4433.645924170054059,0)
0.0 True Vector3(0,0,-0) Vector3(57.85017131547884617,-0.111204314725056902,0)
0.0 True Vector3(0,0,0) Vector3(-65.94555466807858579,-266.2428061520735127,0)
0.0 True Vector3(0,0,0) Vector3(-65.74496288266283273,-41.69951648649276876,0)
0.0 True Vector3(0,0,0) Vector3(3521.296799534754427,1016.086768376104146,0)
0.0 True Vector3(0,0,0) Vector3(3899.798315409473616,-3753.377884397875278,0)
0.0 True Vector3(0,0,-0) Vector3(4776.018541079056376,1378.922310522730641,0)
0.0 True Vector3(0,0,0) Vector3(802.838169049608382,-772.7767530636782567,0)
0.0 True Vector3(0,0,-0) Vector3(2843.329766835983719,-2735.623420421653464,0)
0.0 True Vector3(0,0,0) Vector3(-1600.714997317364578,-6469.500386790699849,0)
0.0 True Vector3(0,0,-0) Vector3(165.0939271744561836,667.629745238112605,0)
0.0 True Vector3(0,0,-0) Vector3(6527.522237594699618,1882.62701116793437,0)
0.0 True Vector3(0,0,0) Vector3(4385.541404957642953,1265.639190142398093,0)
0.0 True Vector3(0,0,-0) Vector3(1226.649577137013011,353.8899420795115134,0)
0.0 True Vector3(0,0,-0) Vector3(-3194.509206210168031,3073.64356591098749,0)
0.0 True Vector3(0,0,0) Vector3(2759.971945101356141,796.6365154151252455,0)
###print after shearing
0.03 True Vector3(0,0,-214.4743318728430665)
Vector3(-8816.246343126065767,-30293.31533904825483,0)
0.03 True Vector3(0,0,176.173187748141288)
Vector3(8089.379084490512469,-5721.041454786492068,0)
0.03 True Vector3(0,0,116.7907157318923197)
Vector3(8325.128956915803428,28982.16753363684984,0)
0.03 True Vector3(0,0,-96.18034735227443832)
Vector3(3947.581959152435047,41620.24546381638356,0)
0.03 True Vector3(0,0,258.172014379664688)
Vector3(8033.969204063756479,27432.31094291702175,0)
0.03 True Vector3(0,0,458.9136141964506805)
Vector3(-7905.467941164684817,-27521.60820142569719,0)
0.03 True Vector3(0,0,-500.9618617967993828)
Vector3(7892.961721389036029,27345.38354414847709,0)
0.03 True Vector3(0,0,192.3317880573678735)
Vector3(8451.970780218442087,29095.3757851041737,0)
0.03 True Vector3(0,0,-181.923089687621939)
Vector3(-25749.52581722119794,14633.4112628354651,0)
0.03 True Vector3(0,0,-466.059931156629716)
Vector3(-7783.585210045047461,-27514.92844738679196,0)
0.03 True Vector3(0,0,-343.9833186713068471)
Vector3(-6710.506162038463117,-39733.66002893296536,0)
0.03 True Vector3(0,0,-30.23573307304739899)
Vector3(-14513.55371828253192,10089.93999155055826,0)
0.03 True Vector3(0,0,-561.0428888761980488)
Vector3(-10903.78008043897353,-25417.82920288496825,0)
0.03 True Vector3(0,0,71.46986470965566696)
Vector3(-1920.847322427021027,1409.197001112073167,0)
0.03 True Vector3(0,0,-28.5325262470542107)
Vector3(4126.450099188157765,-2776.780579063361984,0)
0.03 True Vector3(0,0,-380.8760992530237672)
Vector3(-7962.684475832433236,-27585.51734013333771,0)
0.03 True Vector3(0,0,304.3048195140005419)
Vector3(-26451.29323578983167,15358.69761184588242,0)
0.03 True Vector3(0,0,394.5895764872377072)
Vector3(-9030.924345013181664,-26677.59201235684304,0)
0.03 True Vector3(0,0,-110.2349019237797307)
Vector3(7866.311713633559521,27060.91184042056193,0)
0.03 True Vector3(0,0,83.03516545246040437)
Vector3(-10222.45165068235474,-27852.59223357996234,0)
0.03 True Vector3(0,0,-443.451596141523396)
Vector3(-10655.3821168558643,-27896.63256981404629,0)
0.03 True Vector3(0,0,-566.2844015574634113)
Vector3(-22977.12999113381738,18257.50409055047567,0)
0.03 True Vector3(0,0,841.4209957266232323)
Vector3(26915.20944321317074,-13032.70401250191389,0)
0.03 True Vector3(0,0,428.0135621121204395)
Vector3(8748.939869277201069,27115.40367642032288,0)
0.03 True Vector3(0,0,-554.6589676924617152)
Vector3(24356.62644333421849,-18374.61978310428458,0)
0.03 True Vector3(0,0,-401.7804097563335972)
Vector3(10290.67374970011224,26505.15422501744251,0)
0.03 True Vector3(0,0,205.4580076550560364)
Vector3(-9904.461108112627699,-26653.06965819338802,0)
0.03 True Vector3(0,0,24.42100244068323178)
Vector3(-263.0057603094581964,-770.3754735679190162,0)
0.03 True Vector3(0,0,0.9669539127732386108)
Vector3(-9.735932918598127372,-30.72621607052947823,0)
0.03 True Vector3(0,0,-105.2765914955640767)
Vector3(-8265.101462530456047,5108.503524772123455,0)
0.03 True Vector3(0,0,530.9530353513898717)
Vector3(-8265.610490029659559,-28259.50785004990394,0)
0.03 True Vector3(0,0,-225.0316619171468631)
Vector3(-8681.275558859981174,-30162.49757713791041,0)
0.03 True Vector3(0,0,-1.719516418752231601)
Vector3(-49.22809264870047485,29.35741661190321139,0)
0.03 True Vector3(0,0,-65.11179702162021954)
Vector3(14884.96078048554409,-9365.304116045823321,0)
0.03 True Vector3(0,0,321.4110678513143284)
Vector3(8175.388552229303059,28139.62182729845517,0)
0.03 True Vector3(0,0,-835.7385371978712101)
Vector3(26855.22966644575354,-13003.2132782778117,0)
0.03 True Vector3(0,0,165.7794534400023281)
Vector3(6458.345817565497782,27075.78287987131262,0)
0.03 True Vector3(0,0,-301.6981303296558963)
Vector3(-8381.537062648527353,-28877.33491302212133,0)
0.03 True Vector3(0,0,-417.6020708769462999)
Vector3(7882.076726600218535,27292.79761297579171,0)
0.03 True Vector3(0,0,-36.31847331676575408)
Vector3(-3468.767805706687795,2596.320309882021775,0)
0.03 True Vector3(0,0,274.7339402960665211)
Vector3(-28874.62430638623846,27605.79806072095744,0)
0.03 True Vector3(0,0,19.00374933084872353)
Vector3(3936.536409070291484,-2289.747862948297097,0)
0.03 True Vector3(0,0,-84.21293669378326285)
Vector3(-4212.240358517611639,2737.219394434043352,0)
0.03 True Vector3(0,0,40.13541234301694516)
Vector3(9809.763000628168811,26753.96198014824404,0)
0.03 True Vector3(0,0,29.45752892877754547)
Vector3(2118.439481931869068,-1258.668498614760438,0)
0.03 True Vector3(0,0,-151.2177883890932719)
Vector3(3913.949266266366067,26393.5717130322555,0)
0.03 True Vector3(0,0,-201.6846214619720854)
Vector3(-10413.19326058576735,-27939.0153688458704,0)
0.03 True Vector3(0,0,0.07218472811154619251)
Vector3(2.060078199999796578,-1.243250668134362691,0)
0.03 True Vector3(0,0,-472.2020673523372238)
Vector3(9194.926533437201215,28093.86629321664077,0)
######### TEST No.2 #############
##set initial etaRoll=0.3;
##set 'always_use_moment_law' = False in compaction period;
##set 'always_use_moment_law' = True before shearing start
from __future__ import print_function
from yade import pack,qt,plot
import matplotlib; matplotlib.rc('axes',grid=True)
import pylab
O.materials.append(CohFrictMat(young=1.e9,poisson=.8,alphaKr=1.0,frictionAngle=.0,momentRotationLaw=True,etaRoll=0.3,label='balls'))
sp = pack.SpherePack()
dmax=1.0
size=2*dmax
thickness=size/3.0
sp.makeCloud(minCorner=(0,0,thickness/2),maxCorner=(size,size,thickness/2),rMean=0.1,rRelFuzz=.0,num=50,periodic=True,seed=1)
sp.toSimulation()
O.cell.hSize = Matrix3(size,0,0, 0,size,0, 0,0,thickness) # RVE scale up
print(len(O.bodies))
for p in O.bodies:
p.state.blockedDOFs = 'zXY'
p.state.mass = 2650 * thickness * pi * p.shape.radius**2 # 0.1 = thickness
of cylindrical particle # RVE scale up
inertia = 0.5 * p.state.mass * p.shape.radius**2
p.state.inertia = (.5*inertia,.5*inertia,inertia)
O.dt = utils.PWaveTimeStep()
print(O.dt)
O.engines = [
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom6D()],
[Ip2_CohFrictMat_CohFrictMat_CohFrictPhys()],
[Law2_ScGeom6D_CohFrictPhys_CohesionMoment(always_use_moment_law=False,
useIncrementalForm=True)]
),
PeriTriaxController(
dynCell=True,
goal=(-4.e4,-4.e4,0),
stressMask=3,
relStressTol=.001,
maxUnbalanced=.001,
maxStrainRate=(.5,.5,.0),
doneHook='compactionFinished()',
label='biax'
),
NewtonIntegrator(damping=.1),
]
phase=0
def compactionFinished():
global phase
if phase==0:
#O.engines = O.engines[:3]+O.engines[4:]
for i in O.interactions:
print(i.phys.maxRollPl,i.phys.momentRotationLaw,i.phys.moment_bending,i.phys.normalForce)
phase+=1
elif phase==1:
O.engines[2].lawDispatcher.functors[0].always_use_moment_law=True
setContactFriction(radians(33))
# set the current cell configuration to be the reference one
O.cell.trsf=Matrix3.Identity
# change control type: keep constant confinement in x,y, 20%
compression in z
biax.goal=(-4.e4,-0.05,0)
biax.stressMask=1
# allow faster deformation along x,y to better maintain
stresses # RVE scale up, y value should be scale down a lot
biax.maxStrainRate=(1.0,0.05,0.0)
# next time, call triaxFinished instead of compactionFinished
biax.doneHook='biaxFinished()'
# do not wait for stabilization before calling triaxFinished
biax.maxUnbalanced=10
def biaxFinished():
for i in O.interactions:
print(i.phys.maxRollPl,i.phys.momentRotationLaw,i.phys.moment_bending,i.phys.normalForce)
print('Finished')
O.pause()
###print after compaction
0.03 True Vector3(0,0,0) Vector3(-4356.022418727271543,4190.262910390562865,0)
0.03 True Vector3(0,0,0) Vector3(-946.2413718796913145,-3824.453444819895594,0)
0.03 True Vector3(0,0,0) Vector3(2005.3174721931573,-1929.808863517112059,0)
0.03 True Vector3(0,0,0) Vector3(3190.850202655290104,-3070.077577822232342,0)
0.03 True Vector3(0,0,0) Vector3(2204.808551871029977,8902.688524042694553,0)
0.03 True Vector3(0,0,0) Vector3(4554.374440800664161,-4381.226637213368122,0)
0.03 True Vector3(0,0,0) Vector3(1137.519578345387799,4598.20941843654964,0)
0.03 True Vector3(0,0,0) Vector3(-3344.408293120146936,-965.5064723078584166,0)
0.03 True Vector3(0,0,0) Vector3(3144.948266091866572,907.7799091766298716,0)
0.03 True Vector3(0,0,0) Vector3(2247.276737117122593,9062.611694681330846,0)
0.03 True Vector3(0,0,0) Vector3(593.3228751609747178,2397.947073331891715,0)
0.03 True Vector3(0,0,0) Vector3(-175.8118524235571272,-709.8813002115898598,0)
0.03 True Vector3(0,0,0) Vector3(3130.931207949010513,-3013.271274804336372,0)
0.03 True Vector3(0,0,0) Vector3(588.0179022407435241,2376.474661917997764,0)
0.03 True Vector3(0,0,0) Vector3(1529.626901015044496,6182.86709493719809,0)
0.03 True Vector3(0,0,0) Vector3(-4141.642333725674689,-1195.060840778490729,0)
0.03 True Vector3(0,0,0) Vector3(-1792.775700468065679,1725.117708369524735,0)
0.03 True Vector3(0,0,0) Vector3(-67.53360237268792332,-271.5920647003094928,0)
0.03 True Vector3(0,0,0) Vector3(4135.21852000556737,1194.014624050691964,0)
0.03 True Vector3(0,0,0) Vector3(-3810.761293476583887,-1099.517895306677246,0)
0.03 True Vector3(0,0,0) Vector3(-2025.926394463087945,1949.234293640667602,0)
0.03 True Vector3(0,0,0) Vector3(-2675.563963173430238,2574.376774606120762,0)
0.03 True Vector3(0,0,0) Vector3(5272.416890895456163,-5075.730797074852489,0)
0.03 True Vector3(0,0,0) Vector3(-85.79522611963324152,-216.9439744441962148,0)
0.03 True Vector3(0,0,0) Vector3(-500.4888775940822256,480.3831482944814866,0)
0.03 True Vector3(0,0,0) Vector3(6790.002219449754193,1956.332687601257248,0)
0.03 True Vector3(0,0,0) Vector3(-4785.688251478742131,-1381.727736622080783,0)
0.03 True Vector3(0,0,0) Vector3(581.3424312033406522,2350.305985913877976,0)
0.03 True Vector3(0,0,0) Vector3(-1741.833517726845912,1676.134517511269678,0)
0.03 True Vector3(0,0,0) Vector3(-5647.573462837440275,-1523.645580069464813,0)
0.03 True Vector3(0,0,0) Vector3(-195.9442612484697008,-660.6786580392528094,0)
0.03 True Vector3(0,0,0) Vector3(-482.6886363669209459,425.6844903374354203,0)
0.03 True Vector3(0,0,0) Vector3(5403.182807112851151,-5193.547357533195282,0)
0.03 True Vector3(0,0,0) Vector3(-4037.211573523102743,3885.827380731004723,0)
0.03 True Vector3(0,0,0) Vector3(-5782.157319227181688,-1672.553876581613167,0)
0.03 True Vector3(0,0,0) Vector3(-511.7428450696506275,537.7315348333218026,0)
0.03 True Vector3(0,0,0) Vector3(-412.5180930583922532,-1667.067608065389095,0)
0.03 True Vector3(0,0,0) Vector3(-440.0569357205779397,462.5879847317094118,0)
0.03 True Vector3(0,0,0) Vector3(-79.40905557252953884,-321.1603524112854302,0)
0.03 True Vector3(0,0,0) Vector3(574.9178308474665755,-507.406089794769116,0)
0.03 True Vector3(0,0,0) Vector3(6757.742281344556432,1950.531419672886159,0)
0.03 True Vector3(0,0,0) Vector3(3169.105186400296589,-3049.983867139284484,0)
0.03 True Vector3(0,0,0) Vector3(4934.170380698124063,1424.503266712679988,0)
0.03 True Vector3(0,0,0) Vector3(-3206.510905196595559,-925.5671278104202884,0)
0.03 True Vector3(0,0,0) Vector3(2278.564969719079272,9189.057761753021623,0)
0.03 True Vector3(0,0,0) Vector3(-244.3011851037109068,-987.693651036061965,0)
0.03 True Vector3(0,0,0) Vector3(1613.0290283192885,-1551.795665594201864,0)
0.03 True Vector3(0,0,0) Vector3(-2624.200635546494141,2525.699160412983474,0)
0.03 True Vector3(0,0,0) Vector3(-1894.742558239544678,-7654.642720827022458,0)
0.03 True Vector3(0,0,0) Vector3(-5453.070612922661894,-1574.730209720701396,0)
0.03 True Vector3(0,0,0) Vector3(550.7112834928792608,-530.2544805781109289,0)
0.03 True Vector3(0,0,0) Vector3(-6360.602085101942976,-1836.321129405905594,0)
0.03 True Vector3(0,0,0) Vector3(1278.409820090018457,5166.437438540249786,0)
0.03 True Vector3(0,0,0) Vector3(-501.3525301434042376,-2026.262393923267155,0)
0.03 True Vector3(0,0,0) Vector3(-2271.811739081067572,-9191.198042164664002,0)
0.03 True Vector3(0,0,0) Vector3(2118.342531565449917,-2038.576273405231632,0)
0.03 True Vector3(0,0,0) Vector3(6493.214469140548317,1876.119468843037794,0)
0.03 True Vector3(0,0,0) Vector3(-5400.851266028918872,5200.117448928389422,0)
0.03 True Vector3(0,0,0) Vector3(-2649.949579478481155,2549.48043116090048,0)
0.03 True Vector3(0,0,0) Vector3(-5273.094760494164802,5068.296797843814602,0)
0.03 True Vector3(0,0,0) Vector3(642.908085010988998,2599.278453244339289,0)
0.03 True Vector3(0,0,0) Vector3(-5498.449301278576968,-1587.832342098636673,0)
0.03 True Vector3(0,0,0) Vector3(-989.4557685983584179,-3999.091185152995877,0)
0.03 True Vector3(0,0,0) Vector3(-1629.28978454572939,-6586.428084369147655,0)
0.03 True Vector3(0,0,0) Vector3(1019.938206925837108,4123.007007379715105,0)
0.03 True Vector3(0,0,0) Vector3(-4355.703135797158211,4189.954115430237835,0)
0.03 True Vector3(0,0,0) Vector3(2674.353110719725009,-2573.193619222766301,0)
0.03 True Vector3(0,0,0) Vector3(1257.53382450168715,5082.085015788150486,0)
0.03 True Vector3(0,0,0) Vector3(5671.070048103227236,1741.359121719895256,0)
0.03 True Vector3(0,0,0) Vector3(131.4846787252013769,662.5274659514100222,0)
0.03 True Vector3(0,0,0) Vector3(-3257.224446709202766,3134.789185079573599,0)
0.03 True Vector3(0,0,0) Vector3(-5416.749901594610492,-1564.545644680121768,0)
0.03 True Vector3(0,0,0) Vector3(-2204.759659924044172,-8902.875961359055509,0)
0.03 True Vector3(0,0,0) Vector3(178.7324749728237521,718.8491860076645708,0)
0.03 True Vector3(0,0,0) Vector3(-584.7473032754612632,561.2564641809614159,0)
0.03 True Vector3(0,0,0) Vector3(-6.158419127361661083,26.0734180612055404,0)
0.03 True Vector3(0,0,0) Vector3(582.8097420617934858,-560.4734523253433736,0)
0.03 True Vector3(0,0,0) Vector3(63.03895382092825628,-5.235699592228107768,0)
0.03 True Vector3(0,0,0) Vector3(-1398.024721625517486,-403.3393896789778523,0)
0.03 True Vector3(0,0,0) Vector3(4389.74447505059652,1266.925791252061799,0)
0.03 True Vector3(0,0,0) Vector3(8.508941415349747572,34.414349564750502,0)
0.03 True Vector3(0,0,0) Vector3(495.3726959074963929,-476.3910924743690316,0)
0.03 True Vector3(0,0,0) Vector3(-2971.403754700405898,2859.808144833682491,0)
0.03 True Vector3(0,0,0) Vector3(-547.5349918831619789,527.0987733962641641,0)
0.03 True Vector3(0,0,0) Vector3(-524.0691227176483835,504.4724037886144288,0)
0.03 True Vector3(0,0,0) Vector3(24.54646675697048508,210.3727536274160457,0)
0.03 True Vector3(0,0,0) Vector3(-1214.164619003100142,-4906.390546028451354,0)
0.03 True Vector3(0,0,0) Vector3(-260.111874640118117,-1051.591946580115064,0)
0.03 True Vector3(0,0,0) Vector3(2302.619237149820037,-2214.882467420897228,0)
0.03 True Vector3(0,0,0) Vector3(-4371.180785150242627,4204.733916816947385,0)
0.03 True Vector3(0,0,0) Vector3(-1098.412608630502746,-4439.050616960828847,0)
0.03 True Vector3(0,0,0) Vector3(6790.669141671738544,1963.812786986948595,0)
0.03 True Vector3(0,0,0) Vector3(2484.756067007789625,717.591842242276698,0)
0.03 True Vector3(0,0,0) Vector3(-6143.255952945281024,-1773.664794304884481,0)
0.03 True Vector3(0,0,0) Vector3(-2514.251945581449036,2418.579763056834508,0)
0.03 True Vector3(0,0,0) Vector3(-2548.34829534631217,2452.706293236090914,0)
0.03 True Vector3(0,0,0) Vector3(4094.651600411214986,1181.542155068722195,0)
0.03 True Vector3(0,0,0) Vector3(-1968.954938907677843,-7954.994515831598619,0)
0.03 True Vector3(0,0,0) Vector3(-6472.564400726898384,-1870.002264900543423,0)
0.03 True Vector3(0,0,0) Vector3(1283.919363400329985,5189.285403755833613,0)
0.03 True Vector3(0,0,0) Vector3(-1734.005746328204168,-7010.486056481552623,0)
0.03 True Vector3(0,0,0) Vector3(2673.178560915678645,771.9995470702701823,0)
0.03 True Vector3(0,0,0) Vector3(-3349.020566802118537,-966.8133839347973435,0)
0.03 True Vector3(0,0,0) Vector3(48.1053186181306387,40.59859085302196746,0)
0.03 True Vector3(0,0,0) Vector3(-2721.179722815830246,-785.4536521531675817,0)
0.03 True Vector3(0,0,0) Vector3(-2241.433496616614775,-9067.109488036401672,0)
0.03 True Vector3(0,0,0) Vector3(6134.077892066618915,1769.510159572369503,0)
0.03 True Vector3(0,0,0) Vector3(1097.108672801957937,4433.645924170054059,0)
0.03 True Vector3(0,0,0) Vector3(57.85017131547884617,-0.111204314725056902,0)
0.03 True Vector3(0,0,0) Vector3(-65.94555466807858579,-266.2428061520735127,0)
0.03 True Vector3(0,0,0) Vector3(-65.74496288266283273,-41.69951648649276876,0)
0.03 True Vector3(0,0,0) Vector3(3521.296799534754427,1016.086768376104146,0)
0.03 True Vector3(0,0,0) Vector3(3899.798315409473616,-3753.377884397875278,0)
0.03 True Vector3(0,0,0) Vector3(4776.018541079056376,1378.922310522730641,0)
0.03 True Vector3(0,0,0) Vector3(802.838169049608382,-772.7767530636782567,0)
0.03 True Vector3(0,0,0) Vector3(2843.329766835983719,-2735.623420421653464,0)
0.03 True Vector3(0,0,0) Vector3(-1600.714997317364578,-6469.500386790699849,0)
0.03 True Vector3(0,0,0) Vector3(165.0939271744561836,667.629745238112605,0)
0.03 True Vector3(0,0,0) Vector3(6527.522237594699618,1882.62701116793437,0)
0.03 True Vector3(0,0,0) Vector3(4385.541404957642953,1265.639190142398093,0)
0.03 True Vector3(0,0,0) Vector3(1226.649577137013011,353.8899420795115134,0)
0.03 True Vector3(0,0,0) Vector3(-3194.509206210168031,3073.64356591098749,0)
0.03 True Vector3(0,0,0) Vector3(2759.971945101356141,796.6365154151252455,0)
###print after shearing
0.03 True Vector3(0,0,-456.156872891021294)
Vector3(-14661.1491645947026,-50551.1421652338322,0)
0.03 True Vector3(0,0,191.5794920813423516)
Vector3(5959.500559615104066,-4178.919790439877943,0)
0.03 True Vector3(0,0,102.9315762440743498)
Vector3(13440.36226713041287,46611.41919590990437,0)
0.03 True Vector3(0,0,-409.3002336293453709)
Vector3(2911.792492020297686,38166.21208162078983,0)
0.03 True Vector3(0,0,-40.75776403756982802)
Vector3(13728.4971468439926,47108.97444583152537,0)
0.03 True Vector3(0,0,535.1675585975327749)
Vector3(-13231.41591492013322,-45910.89679806610366,0)
0.03 True Vector3(0,0,-411.8793489654030964)
Vector3(13572.64497299293726,46977.68588109024131,0)
0.03 True Vector3(0,0,121.7188710382010441)
Vector3(13741.94947429019703,47428.8974333227161,0)
0.03 True Vector3(0,0,-19.69673690056714221)
Vector3(2994.560265308720773,-2024.662093290342,0)
0.03 True Vector3(0,0,-444.5079800344556133)
Vector3(-20181.98683050602631,11639.54026002786486,0)
0.03 True Vector3(0,0,-999.2785022751947963)
Vector3(-13727.54938945461436,-48402.97233861769928,0)
0.03 True Vector3(0,0,-588.8319275378147495)
Vector3(-5241.27097675931509,-35362.79578722034785,0)
0.03 True Vector3(0,0,-181.0174737834959728)
Vector3(-11885.23734660024638,8157.132473081320313,0)
0.03 True Vector3(0,0,-1132.219664084880151)
Vector3(-18461.3708385161226,-42988.51846401869261,0)
0.03 True Vector3(0,0,76.89194983684900819)
Vector3(-2402.352635714303233,1759.045763842013457,0)
0.03 True Vector3(0,0,497.463799184994798)
Vector3(-22159.4133038127693,13190.18129945255168,0)
0.03 True Vector3(0,0,165.8217587006903955)
Vector3(-15053.02153083200028,-44632.96708828026749,0)
0.03 True Vector3(0,0,352.4229778019040396)
Vector3(11517.80586767749992,34385.55631481562159,0)
0.03 True Vector3(0,0,-76.39416954965584239)
Vector3(13394.80315565880483,46230.41078894612292,0)
0.03 True Vector3(0,0,-591.6875206775177958)
Vector3(-12457.93446399561253,-33098.08408562492696,0)
0.03 True Vector3(0,0,-721.3944353030937009)
Vector3(-13083.41825527771653,-34462.06593824215815,0)
0.03 True Vector3(0,0,12.24238085242219753)
Vector3(-18299.06544585805386,14530.76968333611694,0)
0.03 True Vector3(0,0,654.8861013673189291)
Vector3(21385.88043460712652,-9618.768455212222761,0)
0.03 True Vector3(0,0,-499.3646622326259603)
Vector3(-13808.49562210596196,-47805.24246019552083,0)
0.03 True Vector3(0,0,-151.4631484986659586)
Vector3(11180.27620064112489,33988.99476485303603,0)
0.03 True Vector3(0,0,-8.197316924984093589)
Vector3(20791.44474459827325,-15517.37988793503428,0)
0.03 True Vector3(0,0,-730.9392014022250805)
Vector3(12284.60784758849513,31856.44124393725724,0)
0.03 True Vector3(0,0,657.8409039327362962)
Vector3(-12768.81440873680731,-33371.8086284124729,0)
0.03 True Vector3(0,0,60.33808067123651142)
Vector3(-6459.161063367578208,4081.558115505431033,0)
0.03 True Vector3(0,0,890.5063503400473337)
Vector3(-14145.01684309212578,-48605.44002568037831,0)
0.03 True Vector3(0,0,-8.873199333858828197)
Vector3(-14196.61767822329602,-49162.03084244739875,0)
0.03 True Vector3(0,0,-128.574332632859182)
Vector3(12153.34452567220615,-7849.496269127708729,0)
0.03 True Vector3(0,0,31.68875356030893897)
Vector3(13630.69107657153836,47057.2578627731782,0)
0.03 True Vector3(0,0,-587.2467513739522929)
Vector3(21642.50430361524195,-9734.634371212461701,0)
0.03 True Vector3(0,0,553.5985692974071526)
Vector3(10974.98989855506625,45860.9089808800054,0)
0.03 True Vector3(0,0,-284.7975541284719156)
Vector3(-13569.40980124090493,-46867.11448837029457,0)
0.03 True Vector3(0,0,-315.2273498622839725)
Vector3(13212.7461601718278,45791.13922822170571,0)
0.03 True Vector3(0,0,60.51517959206349673)
Vector3(-7271.516182375891731,5420.673308076224203,0)
0.03 True Vector3(0,0,250.3014945236063227)
Vector3(-24222.34083532287696,25158.10233777345638,0)
0.03 True Vector3(0,0,160.2837784547463968)
Vector3(7693.653192366281473,-4562.148502936391196,0)
0.03 True Vector3(0,0,-5.711262218682723102)
Vector3(-3247.391729161092826,2128.263358801883442,0)
0.03 True Vector3(0,0,-343.6041825359137647)
Vector3(12378.94883738586759,32724.2395769592149,0)
0.03 True Vector3(0,0,27.75163929978003097)
Vector3(2262.490320530358986,-1365.903227244946038,0)
0.03 True Vector3(0,0,-4.262773808438868883)
Vector3(-109.2012054117491431,90.91404734919984776,0)
0.03 True Vector3(0,0,15.63930588642283404)
Vector3(6702.522905384783371,44996.4658333945772,0)
0.03 True Vector3(0,0,-665.4798273646936195)
Vector3(-13100.62417867341173,-34089.82459094492515,0)
0.03 True Vector3(0,0,-18.50272966843086664)
Vector3(2782.543905807062401,-1736.117902883793477,0)
0.03 True Vector3(0,0,-37.94819380841906309)
Vector3(-3040.047935505311671,2157.095407734354467,0)
--
You received this question notification because your team yade-users is
an answer contact for Yade.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
