Dear Bart, What are you trying to achieve? As Jeremy explained, your geometry is non-conformal. Do you want to create a conformal mesh? If so you'll need to have a conformal geometry first.
Christophe > On 17 Aug 2020, at 10:16, Bart Deschoolmeester <[email protected]> > wrote: > > Hello, > Sorry to bother everybody again. > I’m a beginner and I seem to run in on one error after another. I’m using > GMSH 4.7.0 (I came across GMSH trough the onelab website) on Windows 10. > > See the script below (from a bigger 3D project where I try to narrow down the > errors I get). > I get a good 3D mesh but at the curve where the two objects meet, the mesh is > not aligned. > Using the coherence command I get a large list of errors of the type: > Error : could not find extruded vertex (0.08443279255020152, > -0.053826794978966, 0). > > Yet, if I change the point y-coordinates to a positive number: > Point(23) = {0, 0.1, 0}; > Point(24) = {0, 0.05, 0}; > Point(25) = {0, 0.05, WanD}; > Point(26) = {0, 0.1, WanD}; > > All works as expected. > > Thanks for the help. > Bart > > > ----------------------- > > KanR = DefineNumber[ 0.1, Name "Parameters/KanR" ]; > WanD = DefineNumber[ 0.2, Name "Parameters/WanD" ]; > KanMeshsize = DefineNumber[ 0.01, Name "Parameters/KanMeshsize" ]; > SetFactory("OpenCASCADE"); > Rectangle(1) = {-2.5, -2.5, 0, 5, 5, 0}; > Circle(5) = {0, 0, 0, KanR, 0, 2*Pi}; > Curve Loop(2) = {5}; > Plane Surface(2) = {2}; > BooleanDifference{ Surface{1}; Delete; }{ Surface{2}; Delete; } > Extrude {0, 0, WanD } { > Surface{1}; > } > > Point(23) = {0, -0.1, 0}; > Point(24) = {0, -0.05, 0}; > Point(25) = {0, -0.05, WanD}; > Point(26) = {0, -0.1, WanD}; > Line(32) = {23, 24}; > Line(33) = {24, 25}; > Line(34) = {25, 26}; > Line(35) = {26, 23}; > Curve Loop(30) = {34, 35, 32, 33}; > Plane Surface(14) = {30}; > Extrude {{0, 0, 1}, {0, 0, 0}, Pi} { > Surface{14}; Layers{25}; > } > > Field[1] = Distance; > Field[1].EdgesList = {5}; > Field[1].NNodesByEdge = 50; > Field[2] = Threshold; > Field[2].DistMax = 2.5; > Field[2].DistMin = 0; > Field[2].IField = 1; > Field[2].LcMax = 0.25; > Field[2].LcMin = KanMeshsize ; > > Background Field = 2; > > Mesh.CharacteristicLengthExtendFromBoundary = 0; > Mesh.CharacteristicLengthFromPoints = 0; > Mesh.CharacteristicLengthFromCurvature = 0; > _______________________________________________ > gmsh mailing list > [email protected] > http://onelab.info/mailman/listinfo/gmsh — Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine _______________________________________________ gmsh mailing list [email protected] http://onelab.info/mailman/listinfo/gmsh
