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

Hi, everyone. I have a problem with ‘JCFpmMat’.

What’s the difference between ‘JCFpmMat’ and ‘FrictMat’ if we set the 
tensileStrength and cohesion to 0?

For example, there are two scripts and the difference between them (materials 
and contact InteractionLoop) is as follows:

Script1:

O.materials.append(FrictMat(young = WYoung , poisson = WPoisson , frictionAngle 
= radians(WFrictionAngle) , density = WDensity,label='walls'))
O.materials.append(FrictMat(young = Young , poisson = Poisson , frictionAngle = 
radians(compFricDegree),density = Density,label='spheres'))

O.engines=[
        ForceResetter(),
        InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb()]),
        InteractionLoop(
                [Ig2_Sphere_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()],
                [Ip2_FrictMat_FrictMat_FrictPhys()],
                [Law2_ScGeom_FrictPhys_CundallStrack()]
        ),
        
GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.8),
        triax,
        newton
]

#
#
#

Script 2:

O.materials.append(JCFpmMat(type=0,young = WYoung , poisson = WPoisson , 
frictionAngle = radians(WFrictionAngle) , density = WDensity,label='walls'))
O.materials.append(JCFpmMat(type=1,young = Young , poisson = Poisson , 
frictionAngle = radians(compFricDegree),density = 
Density,tensileStrength=TENS,cohesion=COH,label='spheres'))


triax.goal1=triax.goal2=triax.goal3=-100000
setContactFriction(radians(finalFricDegree))
O.engines=[
        ForceResetter(),
        
InsertionSortCollider([Bo1_Box_Aabb(),Bo1_Sphere_Aabb(aabbEnlargeFactor=intR,label='Saabb')]),
        InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=intR,label='SSgeom'),Ig2_Box_Sphere_ScGeom()],
[Ip2_JCFpmMat_JCFpmMat_JCFpmPhys(cohesiveTresholdIteration=1,label='interactionPhys')],
[Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM(recordCracks=True,Key=OUT,label='interactionLaw')]),
        
GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.8),
        triax,
        newton
]
#
#
#
The two scripts are same except the above content. However, the results of the 
two triaxal tests differed greatly.

-- 
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

Reply via email to