[Gmsh] [Getdp] Functions: symmetric gradients

[Ruth V. Sabariego]

You can't use Dof{eps[{d u}]}, that is not allowed. Sorry!
That is the reason why we split the displacement in components with a basis function per compoment: BF_NodeX, BF_NodeY and BF_NodeZ.


HTH,
Ruth

--
Dr. Ir. Ruth V. Sabariego
University of Liege, Electrical Engineering & Computer Science, 
Applied & Computational Electromagnetics (ACE),
phone: +32-4-3663737 - fax: +32-4-3662910 - http://ace.montefiore.ulg.ac.be/



On 25 Nov 2011, at 00:51, Saurabh Srivastava wrote:

> 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
> 
> 
> 
> 
> _______________________________________________
> getdp mailing list
> getdp at geuz.org
> http://www.geuz.org/mailman/listinfo/getdp

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20111125/8e9d09d9/attachment.html>