[Getdp] Functions: symmetric gradients

Saurabh Srivastava saurabh636 at gmail.com
Fri Nov 25 00:51:40 CET 2011


Hello,

I am interested in the elasticity example, as I realized a new operator was
D1/D2 was hard-coded to obtain symmetric gradient (which are unlike curl,
grad, or div) of displacements. I'm attempting to get around this using the
following approach ...essentially, I compute the exterior derivative of
displacement vector then pass it to a function which converts this 2nd rank
tensor into a vectorial (Voigt's) notation,

Thus I replace these lines ...

Formulation {
  { Name Mec2D_u ; Type FemEquation ;
    Quantity {
      { Name u ; Type Local ; NameOfSpace H_u_Mec2D ; }
    }
    Equation {
    //gu = Tensor[{d u}];
      Galerkin { [ C_m[] * Dof{D1 u}, {D1 u} ] ;
                 In Domain_Disp ; Jacobian Vol ; Integration GradGrad ; }
      Galerkin { [ -F[] , {u} ] ;
                 In Domain_Force; Jacobian SurLinVol; Integration GradGrad;
}
    }
  }
}

-----------------with these lines------------------------
Function {
  eps[] = Vector[CompXX[$1], CompYY[$1], CompXY[$1] ];
}

Formulation {
  { Name Mec2D_uxuy ; Type FemEquation ;
    Quantity {
      //{ Name ux; Type Local ; NameOfSpace H_uxuy_Mec2D ; }
      //{ Name uy; Type Local ; NameOfSpace H_uxuy_Mec2D ; }
      { Name u ; Type Local ; NameOfSpace H_u_Mec2D ; }
    }
    Equation {
      Galerkin { [ C_m[] * Dof{eps[{d u}]}, {eps[{d u}]} ]; //line 175
                 In Domain_Disp ; Jacobian Vol ; Integration GradGrad ; }
      Galerkin { [ -F[] , {u} ] ;
                 In Domain_Force; Jacobian SurLinVol; Integration GradGrad;
}
    }
  }

But I get the following error during pre-processing...

GetDP   : 'C:/Documents and Settings/xx/Desktop/getdp-2.1.1-Win32c/
elasticity2d2.pro', line 175 : syntax error ([)

Doesn't looks like its any bracketing issue as i fiddled around with that,
What am I doing wrong?

Any help is appreciated.

Thank You,
-ss
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/getdp/attachments/20111125/dc0dc3ba/attachment.html>


More information about the getdp mailing list