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