[Gmsh] Question sur GMSH (PERRIN Pierre - CEREMA/DTerEst/Infra/DOA)

Tue Mar 22 17:48:40 CET 2016

> On 22 Mar 2016, at 13:12, PERRIN Pierre - CEREMA/DTerEst/Infra/DOA <Pierre.Perrin at cerema.fr> wrote:
> 
> Bonjour,
> 
> Je vous remercie pour cette information. Finalement il est obligatoire de découper en plusieurs sous-surfaces si on souhaite obtenir un maillage quadrangle régulier.
> 

PS : pour un maillage non-structuré en quad "plus régulier", essayez également 

Mesh.Algorithm = 8; // delquad


> Sinon, en définissant les surfaces par soustraction, j'ai aussi trouvé la solution ci-jointe, mais qui n'est pas aussi régulière car transfinite surface ne s'applique pas si le nombre de coin d'une surface est supérieur à 4.
> Il me reste à tester si ces différentes solutions sont toutes bien acceptées par Code Aster pour résolution du problème EF.
> 
> 
> Merci de votre aide et de votre réactivité,
> 
> Cordialement
> 
> <hbefijgh.png>
> 
> Pierre PERRIN 
> Ingénieur d'études Ouvrages d'Art 
> Département Conception et Exploitation des Infrastructures 
> Division Ouvrages d'Art 
> Tél.: : +33(0)3 87 20 46 39 
> Mobilisé pour sauver le site de Metz
> <logo Cerema.png>
> Centre d’études et d’expertise sur les risques, l’environnement, la mobilité et l’aménagement   www.cerema.fr
> Direction territoriale Est 
> 1, boulevard Solidarité BP 85230 57076 METZ CEDEX 3 - Tél : +33 (0)3 87 20 43 00
> Siège social : Cité des Mobilités - 25, avenue François Mitterrand - CS 92 803 - F-69674 Bron Cedex - Tél : +33 (0)4 72 14 30 30 
> Le 19/03/2016 18:25, > ruard at artenum.com (par Internet) a écrit :
>> A solution could be the attached file. 
>> 
>> Regards, 
>> 
>> Benjamin 
>> 
>> Le 19-03-2016 12:00, gmsh-request at ace20.montefiore.ulg.ac.be a écrit : 
>>> Send gmsh mailing list submissions to 
>>>     gmsh at onelab.info 
>>> 
>>> To subscribe or unsubscribe via the World Wide Web, visit 
>>>     http://onelab.info/mailman/listinfo/gmsh 
>>> or, via email, send a message with subject or body 'help' to 
>>>     gmsh-request at onelab.info 
>>> 
>>> You can reach the person managing the list at 
>>>     gmsh-owner at onelab.info 
>>> 
>>> When replying, please edit your Subject line so it is more specific 
>>> than "Re: Contents of gmsh digest..." 
>>> 
>>> 
>>> Today's Topics: 
>>> 
>>>    1. Question sur GMSH (PERRIN Pierre - CEREMA/DTerEst/Infra/DOA) 
>>> 
>>> 
>>> ---------------------------------------------------------------------- 
>>> 
>>> Message: 1 
>>> Date: Fri, 18 Mar 2016 16:04:43 +0100 
>>> From: "PERRIN Pierre - CEREMA/DTerEst/Infra/DOA" 
>>>     <Pierre.Perrin at cerema.fr> 
>>> To: gmsh at geuz.org 
>>> Subject: [Gmsh] Question sur GMSH 
>>> Message-ID: <56EC190B.3050607 at cerema.fr> 
>>> Content-Type: text/plain; charset="utf-8"; Format="flowed" 
>>> 
>>> Bonjour, 
>>> 
>>> Je me pose une question apparemment simple sur GMSH, mais dont je ne 
>>> parviens pas ? trouver la r?ponse... 
>>> 
>>> Je souhaite mailler un grand rectangle, avec ? l'int?rieur un petit 
>>> rectangle. 
>>> Le petit rectangle ne repr?sente pas un trou, mais un impact de charge. 
>>> Je souhaite in fine avoir acc?s dans Code Aster ? un physical sur cet 
>>> impact pour pouvoir placer un chargement surfacique. 
>>> 
>>> 
>>> 
>>> Je vous joint mon fichier .geo, un peu pollu? par mes essais... Je 
>>> parviens ? mailler apparemment correctement, mais en r?alit? les 2 
>>> maillages sont distincts, ce que Code Aster n'appr?cie pas et me signale. 
>>> Vous serait-il possible de m'indiquer une m?thode qui me permette 
>>> d'arriver ? mes fins ? 
>>> 
>>> 
>>> En vous remerciant par avance, 
>>> 
>>> Cordialement 
>>> 
>>> Pierre PERRIN 
>>> 
>>> *Pour information * 
>>> Ce que je souhaite faire se fait tr?s facilement avec le module ?l?ments 
>>> finis de RDM6 (http://iut.univ-lemans.fr/ydlogi/rdm_version_6.html). 
>>> On d?fini 2 rectangles, puis un menu permet de passer les cot? de 
>>> l'impact en pointill?s, pour signifier simplement que l'on va adapter le 
>>> maillage ? ce contour, mais sans consid?rer de trou dans la plaque : 
>>> 
>>> puis lorsqu'on maille on obtient : 
>>> 
>>> Si on n'a pas d?fini les pointill?s, il y a par contre un trou : 
>>> 
>>> 
>>> -- 
>>> 
>>> *Pierre PERRIN * 
>>> *Ing?nieur d'?tudes Ouvrages d'Art 
>>> D?partement Conception et Exploitation des Infrastructures 
>>> Division Ouvrages d'Art * 
>>> *T?l.: : +33(0)3 87 20 46 39 
>>> * */Mobilis? pour sauver le site de Metz/* 
>>> 
>>> 
>>> Centre d??tudes et d?expertise sur les risques, l?environnement, la 
>>> mobilit? et l?am?nagement www.cerema.fr <http://www.cerema.fr> 
>>> Direction territoriale Est 
>>> 1, boulevard Solidarit? BP 85230 57076 METZ CEDEX 3 - T?l : +33 (0)3 87 
>>> 20 43 00 
>>> Si?ge social : Cit? des Mobilit?s - 25, avenue Fran?ois Mitterrand - CS 
>>> 92 803 - F-69674 Bron Cedex - T?l : +33 (0)4 72 14 30 30 
>>> 
>>> 
>>> 
>>> -------------- next part -------------- 
>>> An HTML attachment was scrubbed... 
>>> URL: 
>>> <http://onelab.info/pipermail/gmsh/attachments/20160318/1e8e28a6/attachment-0001.html> 
>>> -------------- next part -------------- 
>>> //PARAMETRES 
>>> h=0.1; 
>>> cote_x=20.0 ; 
>>> cote_y=10.0 ; 
>>> deb_impact_x=5.0; 
>>> deb_impact_y=5.0; 
>>> impact_x=0.6 ; 
>>> impact_y=0.4 ; 
>>> 
>>> 
>>> //POINTS 
>>> Point(1)={0,0,0,h}; 
>>> Point(2)={cote_x,0,0,h}; 
>>> Point(3)={cote_x,cote_y,0,h}; 
>>> Point(4)={0,cote_y,0,h}; 
>>> Point(5)={deb_impact_x,deb_impact_y,0,h}; 
>>> Point(6)={deb_impact_x+impact_y,deb_impact_y,0,h}; 
>>> Point(7)={deb_impact_x+impact_y,deb_impact_y+impact_y,0,h}; 
>>> Point(8)={deb_impact_x,deb_impact_y+impact_y,0,h}; 
>>> 
>>> //LINES 
>>> Line(1)={1,2}; 
>>> Line(2)={2,3}; 
>>> Line(3)={3,4}; 
>>> Line(4)={4,1}; 
>>> Line(5)={5,6}; 
>>> Line(6)={6,7}; 
>>> Line(7)={7,8}; 
>>> Line(8)={8,5}; 
>>> 
>>> //SURFACES 
>>> Line Loop(1)={1,2,3,4}; 
>>> Plane Surface(1)={1}; 
>>> Transfinite Surface{1}; 
>>> Recombine Surface{1}; 
>>> 
>>> // incorporation des lignes dans les surfaces : 
>>> // En GMSH ancien (v1.60) 
>>> //    ?? 
>>> // En GMSH moderne (v2.12) 
>>> //    Line {5,6,7,8} In Surface {1}; 
>>> 
>>> Line {5} In Surface {1}; 
>>> Line {6} In Surface {1}; 
>>> Line {7} In Surface {1}; 
>>> Line {8} In Surface {1}; 
>>> 
>>> 
>>> Line Loop(2)={5,6,7,8}; 
>>> Plane Surface(2)={2}; 
>>> Transfinite Surface{2}; 
>>> Recombine Surface{2}; 
>>> 
>>> 
>>> // Essais inutile au final : 
>>> //Compound Surface(3)={1,2}; 
>>> //Transfinite Surface{3}; 
>>> //Mesh.RemeshAlgorithm=1; 
>>> Coherence; 
>>> Mesh.Format=1; 
>>> 
>>> //PHYSICAL 
>>> Physical Line(1)={1, 3};                 //section d'encastrement 
>>> Physical Line(2)={2, 4};                 //section d'encastrement 
>>> Physical Surface(10000)={1}; 
>>> Physical Surface(20000)={2}; 
>>> Physical Surface(30000)={1,2}; 
>>> 
>>> //Color Yellow { Surface {1} ; } 
>>> //Color Red { Line {2} ; } 
>>> 
>>> 
>>> 
>>> 
>>> 
>>> ------------------------------ 
>>> 
>>> Subject: Digest Footer 
>>> 
>>> _______________________________________________ 
>>> gmsh mailing list 
>>> gmsh at onelab.info 
>>> http://onelab.info/mailman/listinfo/gmsh 
>>> 
>>> 
>>> ------------------------------ 
>>> 
>>> End of gmsh Digest, Vol 158, Issue 15 
>>> ************************************* 
>> 
>> 
>> 
> 
> <Charg_plaque_impact_PERRIN_ok_non_transfinite.geo>_______________________________________________
> gmsh mailing list
> gmsh at onelab.info
> http://onelab.info/mailman/listinfo/gmsh

-- 
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science 
http://www.montefiore.ulg.ac.be/~geuzaine

Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info