<div>Hi,<br><br></div><div>I need to visualize the stress results of my fem calculation</div><div>for a 4 node high-order element and I am using </div><div><br></div><div>$ElementNodeData</div><div>....</div><div><br></div>
<div>for plotting the stress field and other.</div><div><br></div><div><div>But, in this way I can depict only a linear function</div><div>between two adjacent nodes instead I would like to visualize the exact </div><div>
solution in term of stress that is quadratic. How can I do that?</div></div><div><br><br></div><div>Best Regards, </div><div> Giuseppe</div><div><br></div><div><br><div class="gmail_quote">2010/5/12 Giuseppe Zagari <span dir="ltr"><<a href="mailto:gzagari@unical.it">gzagari@unical.it</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div><br></div><div>Hi David, </div><div>thank you for the reply.</div><div><br></div><div>Something like this:</div><div>
<br></div><div>i#include <stdlib.h></div><div>#include <stdio.h></div><div>#include <sys/types.h></div>
<div>#include <signal.h></div><div>#include "MySolver.h"</div><div><br></div><div> int main(int argc, char* argv[])</div><div>{</div><div> int i, pid;</div><div><div><br></div><div> system("gmsh -pid -listen test.geo > /tmp/gmsh.pid &");</div>
<div> //system("gmsh -pid -listen -2 test.geo > /tmp/gmsh.pid &");</div><div> </div></div><div> for(i = 0; i < 5; i++)<br> { <br> system(" ??? "); // put here the command for refine</div>
<div> MySolver ms("test");</div><div> ms.analyze();</div><div> } </div><div><br></div><div> kill(pid, 9); </div><div>
} </div><br><div>I hope I understod right.</div><div><br></div><div>Best Regards, </div><div>GZ<br><div><br></div><div><br><div class="gmail_quote">2010/5/12 David Colignon <span dir="ltr"><<a href="mailto:David.Colignon@ulg.ac.be" target="_blank">David.Colignon@ulg.ac.be</a>></span><div>
<div></div><div class="h5"><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><br>
On 12/05/10 09:56, Giuseppe Zagari wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi,<br>
<br>
I am using gmsh gui for checking<br>
convergence of my finite elements following this procedure:<br>
open gmsh, open filename, go to mesh module, clicking 2D and refine<br>
many times as I want and run the analysis.<br>
</blockquote>
<br>
<br></div>
Hi Giuseppe,<br>
<br>
I suppose you mean that after refining, you save the mesh, perform your FE calculations, go back to the Gmsh GUI, click on refine, save, etc. ?<div><br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
But now I would like to use it in a script like this:<br>
<br>
gmsh -2 filename.geo<br>
analyze filename // it use both .geo and .msh files<br>
gmsh -refine filename.msh<br>
analyze filename.msh<br>
gmsh -refine filename.msh<br>
analyze filename.msh<br>
....<br>
<br>
but, used in this way, gmsh change the node numbering.<br>
</blockquote>
<br>
<br></div>
yes, if you save the mesh, quit Gmsh, relaunch Gmsh,reload the mesh, refine, save, quit, etc. like in your script above, the node numbering will probably change.<br>
<br>
One solution could be to launch Gmsh in "server" mode (launch Gmsh on the command line with "gmsh -listen") and interact with it through "sockets", like explained here (point 3):<br>
<br>
<a href="http://www.geuz.org/gmsh/doc/texinfo/gmsh.html#Solver-module-questions" target="_blank">http://www.geuz.org/gmsh/doc/texinfo/gmsh.html#Solver-module-questions</a><br>
<br>
or launch your FE code from within Gmsh :<br>
<br>
<a href="http://www.geuz.org/gmsh/doc/texinfo/gmsh.html#Solver" target="_blank">http://www.geuz.org/gmsh/doc/texinfo/gmsh.html#Solver</a><br>
<br>
<a href="http://www.geuz.org/gmsh/doc/texinfo/gmsh.html#Solver-module" target="_blank">http://www.geuz.org/gmsh/doc/texinfo/gmsh.html#Solver-module</a><br>
<br>
Regards,<br>
<br>
Dave<br>
<br>
-- <br>
David Colignon, Ph.D.<br>
Collaborateur Logistique du F.R.S.-FNRS<br>
CÉCI - Consortium des Équipements de Calcul Intensif<br>
ACE - Applied & Computational Electromagnetics<br>
Sart-Tilman B28<br>
Université de Liège<br>
4000 Liège - BELGIQUE<br>
Tél: +32 (0)4 366 37 32<br>
Fax: +32 (0)4 366 29 10<br>
WWW: <a href="http://hpc.montefiore.ulg.ac.be/" target="_blank">http://hpc.montefiore.ulg.ac.be/</a><br>
Agenda: <a href="http://www.google.com/calendar/embed?src=david.colignon%40gmail.com" target="_blank">http://www.google.com/calendar/embed?src=david.colignon%40gmail.com</a><br>
<br>
<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>
<br>
Please suggest me in wich way I can obtaine the same<br>
behavoiur of the gmsh-gui.<br>
<br>
<br>
Best Regards and Many Thanks for gmsh<br>
<br>
Giuseppe<br>
<br>
<br>
--<br>
Ing. Giuseppe Zagari, Ph. D.<br>
Dip. Modellistica per l'Ing. - UNICAL<br>
Via P. Bucci, Cubo 39/b 87036<br>
Arcavacata di Rende CS ITALY<br>
<br>
<br>
<br></div>
_______________________________________________<br>
gmsh mailing list<br>
<a href="mailto:gmsh@geuz.org" target="_blank">gmsh@geuz.org</a><br>
<a href="http://www.geuz.org/mailman/listinfo/gmsh" target="_blank">http://www.geuz.org/mailman/listinfo/gmsh</a><br>
</blockquote>
<br>
</blockquote></div></div></div><br><br clear="all"><br>-- <br><div class="im">Ing. Giuseppe Zagari, Ph. D.<br>Dip. Modellistica per l'Ing. - UNICAL<br>Via P. Bucci, Cubo 39/b 87036 <br>Arcavacata di Rende CS ITALY<br>
<br>
</div></div></div>
</blockquote></div><br><br clear="all"><br>-- <br>Ing. Giuseppe Zagari, Ph. D.<br>Dip. Modellistica per l'Ing. - UNICAL<br>Via P. Bucci, Cubo 39/b 87036 <br>Arcavacata di Rende CS ITALY<br><br>
</div>