Commit 138e652f authored by Felipe Bordeu's avatar Felipe Bordeu
Browse files

(MaterialHelp.py) axisymetric HookeIso work

parent 79e985a8

src/BasicTools/FE/MaterialHelp.py

+9 −9
Original line number Diff line number Diff line
@@ -76,12 +76,6 @@ class HookeLaw() : 
        if dim == 1:

            return np.array([[E],])



        if dim == 2:

            if planeStress:

              return ((E/(1.-nu**2.))*

              np.array([[1. , nu, 0     ],

                        [nu, 1. , 0     ],

                        [0 , 0 , (1.-nu)/2.]]));

        if axisymetric:

            k = nu/(1-nu)

            s = (1-2*nu)/(2*(1-nu) )
@@ -91,6 +85,12 @@ class HookeLaw() : 
                        [k , k , 1, 0],

                        [0 , 0 , 0, s]]));



        if dim == 2:

            if planeStress:

              return ((E/(1.-nu**2.))*

              np.array([[1. , nu, 0     ],

                        [nu, 1. , 0     ],

                        [0 , 0 , (1.-nu)/2.]]));

            else:

              return  ((E/((1.+nu)*(1-2*nu)))*

              np.array([[1.-nu, nu  , 0 ],