Hello,<br><br>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, <br>
<br>Thus I replace these lines ...<br><br>Formulation {<br> { Name Mec2D_u ; Type FemEquation ;<br> Quantity {<br> { Name u ; Type Local ; NameOfSpace H_u_Mec2D ; }<br> }<br> Equation {<br> //gu = Tensor[{d u}];<br>
Galerkin { [ C_m[] * Dof{D1 u}, {D1 u} ] ;<br> In Domain_Disp ; Jacobian Vol ; Integration GradGrad ; } <br> Galerkin { [ -F[] , {u} ] ; <br> In Domain_Force; Jacobian SurLinVol; Integration GradGrad; } <br>
}<br> }<br>}<br><br>-----------------with these lines------------------------<br>Function {<br> eps[] = Vector[CompXX[$1], CompYY[$1], CompXY[$1] ];<br>}<br><br>Formulation {<br> { Name Mec2D_uxuy ; Type FemEquation ;<br>
Quantity {<br> //{ Name ux; Type Local ; NameOfSpace H_uxuy_Mec2D ; }<br> //{ Name uy; Type Local ; NameOfSpace H_uxuy_Mec2D ; }<br> { Name u ; Type Local ; NameOfSpace H_u_Mec2D ; }<br> }<br> Equation {<br>
Galerkin { [ C_m[] * Dof{eps[{d u}]}, {eps[{d u}]} ]; //line 175<br> In Domain_Disp ; Jacobian Vol ; Integration GradGrad ; } <br> Galerkin { [ -F[] , {u} ] ; <br> In Domain_Force; Jacobian SurLinVol; Integration GradGrad; } <br>
}<br> }<br><br>But I get the following error during pre-processing...<br><br>GetDP : 'C:/Documents and Settings/xx/Desktop/getdp-2.1.1-Win32c/<a href="http://elasticity2d2.pro">elasticity2d2.pro</a>', line 175 : syntax error ([)<br>
<br>Doesn't looks like its any bracketing issue as i fiddled around with that, What am I doing wrong?<br><br>Any help is appreciated.<br><br>Thank You,<br>-ss<br><br><br><br><br>