[Gmsh] Question sur GMSH (PERRIN Pierre - CEREMA/DTerEst/Infra/DOA)
PERRIN Pierre - CEREMA/DTerEst/Infra/DOA
Pierre.Perrin at cerema.fr
Wed Mar 23 15:42:38 CET 2016
Bonjour,
Merci pour cette piste supplémentaire, qui me semble relativement
régulière, tout en restant pratique (cela évite le découpage en
plusieurs sous-face). J’obtiens ceci :
Avec la documentation de GMSH, il est vrai que malheureusement on a du
mal à savoir ce qui se cache derrière les nombreuses commandes
"Mesh.Option List" : des images seraient commodes, même si cela
alourdirait la documentation...
Toujours sur ma problématique, je ne connait pas bien le fonctionnement
des algorithmes de maillage de la littérature, mais je ne peux pas
m'empêcher de penser que l'on pourrait automatiser la détermination d'un
maillage uniquement composé de carrés de taille différentes.
Ne serait-ce pas intéressant de développer un tel algorithme pour
mailler (uniquement) les surfaces 2D rectangulaires, comprenant 1 ou
plusieurs impacts rectangulaires ?
Après je me doute que ce n'est pas nécessairement votre priorité.
Je profite de nos échanges pour vous féliciter et vous remercier au
sujet de GMSH : cet outils libre est vraiment très performant !
Cordialement
*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
Le 22/03/2016 17:48, > Christophe Geuzaine (par Internet) a écrit :
>> 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://onelab.info/pipermail/gmsh/attachments/20160323/dd553a7b/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ggibjfdc.png
Type: image/png
Size: 7325 bytes
Desc: not available
URL: <http://onelab.info/pipermail/gmsh/attachments/20160323/dd553a7b/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: logo Cerema
Type: image/png
Size: 6094 bytes
Desc: not available
URL: <http://onelab.info/pipermail/gmsh/attachments/20160323/dd553a7b/attachment-0003.png>
-------------- next part --------------
//PARAMETRES
h=0.2;
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};
//Transfinite Line{1,2,3,4,5,6,7,8} = 10
//SURFACES
Line Loop(1)={1,2,3,4};
Line Loop(2)={5,6,7,8};
Plane Surface(1)={1,2}; // la surface 1 est definie avec un trou consititué par la lineloop 2
//Transfinite Surface{1}; // le transfinite ne marche pas pas car la surface a 8 coins
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};
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;
Mesh.Algorithm = 8; // delquad
//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} ; }
More information about the gmsh
mailing list