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
<[email protected]> 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, > [email protected] (par Internet) a écrit :
A solution could be the attached file.
Regards,
Benjamin
Le 19-03-2016 12:00, [email protected] a écrit :
Send gmsh mailing list submissions to
[email protected]
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
[email protected]
You can reach the person managing the list at
[email protected]
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"
<[email protected]>
To: [email protected]
Subject: [Gmsh] Question sur GMSH
Message-ID: <[email protected]>
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
[email protected]
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
[email protected]
http://onelab.info/mailman/listinfo/gmsh
//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} ; }
_______________________________________________
gmsh mailing list
[email protected]
http://onelab.info/mailman/listinfo/gmsh
