Hello,
Your Surface is not "in" the volume but on its boundary, your surface is
counted twice.
Try BooleanFragments, this code should work:
SetFactory("OpenCASCADE");
Box(1) = {0, 0, 0, 100, 10, 10};
Rectangle(7) = {40,0,0,20,10,0};
BooleanFragments{Surface{7}; Delete;}{Volume{1}; Delete;}
Physical Volume("beam") = {1};
Physical Surface("left") = {8};
Physical Surface("right") = {9};
Physical Surface("load") = {7};
Best,
Le 13/02/2020 à 12:35, Vincent Libertiaux a écrit :
Hi, gmsh community,
I created a simple case (beam fixed at both extremities) with a
distributed load on a fraction of its upper surface. I reproduce the
short *.geo file:
SetFactory("OpenCASCADE");
Box(1) = {0, 0, 0, 100, 10, 10};
Rectangle(7) = {40,0,0,20,10,0};
Surface{7} In Volume{1};
surfaces[] = Boundary{Volume{1};};
Physical Volume("beam") = {1};
Physical Surface("left") = {surfaces[0]};
Physical Surface("right") = {surfaces[1]};
Physical Surface("load") = {7};
When I try to mesh it, I get the following error:
Error : PLC Error: A segment and a facet intersect at point
Error : Invalid boundary mesh (segment-facet intersection) on
surface 7, intersection (40,4.7255,6.14662e-11)
I think it comes from the fact that I am trying to impose a surface on
another, but I can't find an easy way (rather than defining the
geometry with Points, Lines, ... or determining the nodes on which I'd
like to apply the loading from the msh file).
Could someone please point me the right direction ?
Thanks,
kind regards,
Vincent
_______________________________________________
gmsh mailing list
[email protected]
http://onelab.info/mailman/listinfo/gmsh
--
Bertrand THIERRY
CNRS researcher
Laboratoire Jacques Louis Lions
Office 310, Corridor 15-25
+33 1 44 27 72 01
http://bthierry.pages.math.cnrs.fr/
_______________________________________________
gmsh mailing list
[email protected]
http://onelab.info/mailman/listinfo/gmsh
