yade-users team mailing list archive

[liucheng83 <question204380@xxxxxxxxxxxxxxxxxxxxx>]

Question #204380 on Yade changed:
https://answers.launchpad.net/yade/+question/204380

    Status: Answered => Open

liucheng83 is still having a problem:
Hi, Jan
Thank you for you reply,  I thought it maybe there is two sim1 in description in multi.table, but just 9 line alone, the problem is still being here.
 And I try the yade-stable with the version "Welcome to Yade 0.80.1+2+10~precise1, from 2012-06-01"  in ub 12.04. The problem is still here, sometimes it is ok with one line in multi.table.

The multi.table is here 
-----------------------
!OMP_NUM_THREADS description nSpheres young
1	sim1	500	2.00E+006
1	sim2	500	4.00E+006
1	sim3	500	6.00E+006
1	sim4	500	8.00E+006
1	sim5	500	1.00E+007
1	sim6	500	1.20E+007
1	sim7	500	1.40E+007
1	sim8	500	1.60E+007
1	sim9	500	1.80E+007
1	sim1	600	2.00E+006
1	sim2	600	4.00E+006
1	sim3	600	6.00E+006
1	sim4	600	8.00E+006
1	sim5	600	1.00E+007
1	sim6	600	1.20E+007
1	sim7	600	1.40E+007
1	sim8	600	1.60E+007
1	sim9	600	1.80E+007
1	sim1	700	2.00E+006
1	sim2	700	4.00E+006
1	sim3	700	6.00E+006
1	sim4	700	8.00E+006
1	sim5	700	1.00E+007
1	sim6	700	1.20E+007
1	sim7	700	1.40E+007
1	sim8	700	1.60E+007
1	sim9	700	1.80E+007

----------------------
The  periodic-triax.py is here
-----------------------

# 2012-07-15 lc
#"Test and demonstrate use of PeriTriaxController."
from yade import *
#from yade import pack,qt,plot
from yade import pack,plot
O.periodic=True

O.load('./ggg/pack_ini'+'.xml')

############################################
###   DEFINING VARIABLES AND MATERIALS   ###
############################################

# The following 5 lines will be used later for batch execution
nRead=utils.readParamsFromTable(
	nSpheres=500,# number of spheres
	young=5e6,
	unknownOk=True
)
from yade.params import table

young=table.young
poisson=0.3
density=2700
nSpheres=table.nSpheres # number of spheres
compFricDegree = 2 # contact friction during the confining phase
finalFricDegree = 30  # contact friction during the deviatoric loading
compPressure= 1e6


# put you simulation's name here
key='E'+str(young) 
key_N='N'+str(nSpheres)+'E'+str(young) 
	    
O.materials[0].young=young


O.dt=0.5*utils.PWaveTimeStep()

O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb()]),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom()],
		[Ip2_FrictMat_FrictMat_FrictPhys()],
		[Law2_ScGeom_FrictPhys_CundallStrack()]
	),
	PeriTriaxController(dynCell=True,mass=mass,maxUnbalanced=0.01,relStressTol=0.001,goal=(-1*compPressure,-1*compPressure,-1*compPressure),stressMask=7,globUpdate=5,maxStrainRate=(.1,.1,.1),doneHook='triaxDone()',label='triax'),
	NewtonIntegrator(damping=.2),
]


def addPlotData():
	plot.addData(unbalanced=utils.unbalancedForce(),i=O.iter,
		sxx=-triax.stress[0],syy=-triax.stress[1],szz=-triax.stress[2],
		exx=-triax.strain[0],eyy=-triax.strain[1],ezz=-triax.strain[2],
		# save all available energy data
		Etot=O.energy.total(),**O.energy
	)

# enable energy tracking in the code
O.trackEnergy=True


# define another plot
plot.plots={'i':('exx','eyy','ezz',None,'sxx','syy','szz'),
     ' i':('ezz')   
}

phase=0
def triaxDone():
	global phase
	if phase==0:	
		O.cell.trsf=Matrix3.Identity
		O.engines=O.engines+[PyRunner(command='addPlotData()',iterPeriod=100)]
		triax.goal=(-1*compPressure,-1*compPressure,-0.2)
		triax.stressMask=3
	        # allow faster deformation along x,y to better maintain stresses
	        triax.maxStrainRate=(1.,1.,.01) 
		phase+=1
	elif phase==1:
		print 'Here we are: stress',triax.stress,'strain',triax.strain,'stiffness',triax.stiff
		print 'Done, pausing now.'
		O.pause()		
		plot.saveDataTxt('./ggg/results'+key+O.tags['description'])
		
O.run()
utils.waitIfBatch()
----------------------

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