'Vector3' object is not callable

Asked by Rinaldi Giuseppe on 2012-01-25

Hi all,

   I'm trying to modify the utils.normShearStressTensors in order to obtain the local stress in a periodic Cell during a simpleShear test.

I define a function LocalStress.
It has as input the coordinates of the 2 corner of the box where I want to evaluate the stress.

this is the function:

* Return the stress tensor decomposed in 2 contributions, from normal and shear forces.
The formulation follows the [Thornton2000]_ article
"Numerical simulations of deviatoric shear deformation of granular media", eq (3) and (4)
 */

def localStress(undCorner, upCorner): #inputs are the 2 corners of the box where we want to evaluate the stress

         for k in O.interactions:

            if (k.isReal): # check if is a real interaction

              mPt=k.geom.contactPoint # assign to mPt the coordinates of the contact point

              if mPt(0)<upCorner(0) and mPt(0)<undCorner(0): # check if the contact point is inside the box
                     if mPt(1)<upCorner(1) and mPt(1)<undCorner(1):
                            if mPt(2)<upCorner(2) and mPt(2)<undCorner(2):

                   n=k.geom.normal

                                 # if compression has positive sign, we need to change sign for both normal and shear force
                   # make copy to Fs since shearForce is used at multiple places (less efficient, but avoids confusion)

                   Fs=k.phys.shearForce

                   N=k.phys.normalForce.n

                                 T=Fs.norm()

                                 if(T>0): t=Fs/T

                   R=.5*(k.geom.refR1+k.geom.refR2)

              for i in range(0,3):
                 for j in range(i,3):
                    sigNorm[i,j]=R*N*n[i]*n[j]
             sigTan[i,j]=R*T*n[i]*t[j]

         vol=(upCorner(0)-undCorner(0))*(upCorner(1)-undCorner(1))*(upCorner(1)-undCorner(1)); # evaluate the volume of the box
         sigNorm=2/vol*sigNorm; # evaluate the stress tensor due to normal force
         sigTan=2/vol*sigTan; # evaluate the stress tensor due to tangential force

         #// fill terms under the diagonal
         sigNorm[1,0]=sigNorm[0,1]
         sigNorm[2,0]=sigNorm[0,2]
         sigNorm[2,1]=sigNorm[1,2]
         sigTan[1,0]=sigTan[0,1]
         sigTan[2,0]=sigTan[0,2]
         sigTan[2,1]=sigTan[1,2]

         return (sigNorm,sigTan)

then in the script:

LocStress=localStress([0.0001,0.0001,0.0001],[0.0002,0.0002,0.0002])

and this is what I get:

Yade [1]: ---------------------------------------------------------------------------
TypeError Traceback (most recent call last)

/usr/lib/yade-daily/py/yade/__init__.pyc in <module>()
----> 1
      2
      3
      4
      5

/usr/lib/yade-daily/py/yade/__init__.pyc in saveData()
    197 damping=contactLaw.normDampDissip # viscous damping
    198 # -------------------------------------------------------------------------

--> 199 LocStress=localStress((0.00001,0.00001,0.00001),(0.00002,0.00002,0.00002))
    200
    201 plot.addData(

/usr/lib/yade-daily/py/yade/__init__.pyc in localStress(undCorner, upCorner)
     55 mPt=k.geom.contactPoint # assign to mPt the coordinates of the contact point
     56
---> 57 if mPt(0)<upCorner(0) and mPt(0)<undCorner(0): # check if the contact point is inside the box
     58 if mPt(1)<upCorner(1) and mPt(1)<undCorner(1):
     59 if mPt(2)<upCorner(2) and mPt(2)<undCorner(2):

TypeError: 'Vector3' object is not callable

What's wrong?

Thanks for any suggestions,

Giuseppe

Question information

Language:
English Edit question
Status:
Solved
For:
Yade Edit question
Assignee:
No assignee Edit question
Solved by:
Rinaldi Giuseppe
Solved:
2012-01-27
Last query:
2012-01-27
Last reply:
2012-01-25

Use square brackets: corner[0] instead of corner(0).

Rinaldi Giuseppe (ciosrin) said : #2

thanks Bruno Chareyre , that solved my question.