Error: Could not replace curve in Coherence
Hi!
I'm using the Extrude Rotate command to create a nozzle geometry. The geo file I'm using is as follows. The error says the original curve for the rotation could not be replaced in coherence. I cannot understand why. So I tried another simple geo file to verify my operations. The simple one is also listed as follows.
In fact, if I ignore such error and do the following operations to create my mesh. I get an unexpected mesh. The complete geo file is also attached. If possible, please have a look.
The Gmsh version is as follows,
Gmsh
version 2.14.0
Copyright (C) 1997-2016
Christophe Geuzaine and Jean-Francois Remacle
Credits and licensing information
Please send all questions and bug reports to the public mailing list gmsh@onelab.info
* Build OS: Linux
* Build date: 20161009
* Build host: debian
* Build options: Ann Bamg Bfgs Blas(Generic) Blossom Chaco DIntegration Dlopen Fltk Gmm Jpeg(Fltk) Kbipack Lapack(Generic) LinuxJoystick MathEx Med Mesh Metis Mmg3d Mpeg Netgen ONELAB ONELABMetamodel OpenCascade OpenGL OptHom PETSc Parser Plugins Png(Fltk) Post SLEPc Salome Solver Taucs TetGen/BR Tetgen1.5 Voro3D Zlib
* FLTK version: 1.3.4
* PETSc version: 3.5.2 (complex arithmetic)
* OCC version: 6.8.0
* MED version: 3.2.0
* Packaged by: geuzaineVisit http://gmsh.info for more information
Please help me! Thanks!
// The complex one leads to error.
Point(0)={-7.740000,3.000000,0};
Point(1)={-7.331390,3.000000,0};
Point(2)={-6.904915,2.996306,0};
Point(3)={-6.465840,2.954973,0};
Point(4)={-6.021626,2.871681,0};
Point(5)={-5.578380,2.753557,0};
Point(6)={-5.140597,2.608993,0};
Point(7)={-4.711511,2.446511,0};
Point(8)={-4.293491,2.274156,0};
Point(9)={-3.888378,2.099123,0};
Point(10)={-3.497692,1.927602,0};
Point(11)={-3.122748,1.764730,0};
Point(12)={-2.764733,1.614520,0};
Point(13)={-2.424742,1.479838,0};
Point(14)={-2.103754,1.362471,0};
Point(15)={-1.802524,1.263442,0};
Point(16)={-1.521763,1.182470,0};
Point(17)={-1.261841,1.118833,0};
Point(18)={-1.022890,1.071132,0};
Point(19)={-0.804744,1.037504,0};
Point(20)={-0.606929,1.015842,0};
Point(21)={-0.428671,1.004021,0};
Point(22)={-0.268948,1.000022,0};
Point(23)={-0.126539,0.999960,0};
Point(24)={0.000000,0.999960,0};
Spline(0)={0:24};
Point(25)={-12.000000,4.250000,0};
Point(26)={-11.072160,4.250000,0};
Point(27)={-10.195415,4.250000,0};
Point(28)={-9.356321,4.250000,0};
Point(29)={-8.542011,4.250000,0};
Point(30)={-7.740000,4.250000,0};
Point(31)={-6.807604,4.250000,0};
Point(32)={-5.591946,4.178856,0};
Point(33)={-4.340811,3.536167,0};
Point(34)={-3.120568,2.831660,0};
Point(35)={-2.066436,2.223057,0};
Point(36)={-1.263920,1.759724,0};
Point(37)={-0.709796,1.439801,0};
Point(38)={-0.352487,1.233509,0};
Point(39)={-0.132151,1.106297,0};
Point(40)={-0.000000,1.022424,0};
Point(41)={0.000000,0.999960,0};
Spline(1)={25:40};
Line(2)={40,41};
Point(42)={-12.000000,-50.000000,-50.000000,1.0};
Point(43)={-12.000000,-50.000000,50.000000,1.0};
Point(44)={-12.000000,50.000000,50.000000,1.0};
Point(45)={-12.000000,50.000000,-50.000000,1.0};
Point(50)={-7.740000,-0.500000,-0.500000,1.0};
Point(51)={-7.740000,-0.500000,0.500000,1.0};
Point(52)={-7.740000,0.500000,0.500000,1.0};
Point(53)={-7.740000,0.500000,-0.500000,1.0};
//+
Rotate {{1, 0, 0}, {0, 0, 0}, -Pi*0.75} {
Line{1}; Line{0}; Line{2};
}
//+
Line(3) = {50, 51};
Line(4) = {51, 52};
Line(5) = {52, 53};
Line(6) = {53, 50};
Line Loop(1) = {4, 5, 6, 3};
Plane Surface(1) = {1};
Extrude {7.74, 0, 0} {
Surface{1};
}
//+
Extrude {{1, 0, 0}, {0, 0, 0}, Pi*0.5} {
Line{0,1};
}
//+
Extrude {{1, 0, 0}, {0, 0, 0}, Pi/2} {
Line{29}; Line{33};
}
//+
Extrude {{1, 0, 0}, {0, 0, 0}, Pi/2} {
Line{41}; Line{37};
}
//+
Extrude {{1, 0, 0}, {0, 0, 0}, Pi/2} {
Line{45}; Line{49};
}
// The simple one has no error.
Point(0)={-2,2,0};
Point(1)={0,2,0};
Point(2)={0,1,0};
Point(3)={-1,1,0};
//+
Line(1) = {0, 1};
//+
Line(2) = {1, 2};
//+
Line(3) = {2, 3};
//+
Extrude {{1, 0, 0}, {0, 0, 0}, Pi/2} {
Line{1, 2, 3};
}
//+
Extrude {{1, 0, 0}, {0, 0, 0}, Pi/2} {
Line{4, 8, 12};
}
//+
Extrude {{1, 0, 0}, {0, 0, 0}, Pi/2} {
Line{24, 20, 16};
}
//+
Extrude {{1, 0, 0}, {0, 0, 0}, Pi/2} {
Line{28, 32, 36};
}
// The complete version
Point(0)={-7.740000,3.000000,0};
Point(1)={-7.331390,3.000000,0};
Point(2)={-6.904915,2.996306,0};
Point(3)={-6.465840,2.954973,0};
Point(4)={-6.021626,2.871681,0};
Point(5)={-5.578380,2.753557,0};
Point(6)={-5.140597,2.608993,0};
Point(7)={-4.711511,2.446511,0};
Point(8)={-4.293491,2.274156,0};
Point(9)={-3.888378,2.099123,0};
Point(10)={-3.497692,1.927602,0};
Point(11)={-3.122748,1.764730,0};
Point(12)={-2.764733,1.614520,0};
Point(13)={-2.424742,1.479838,0};
Point(14)={-2.103754,1.362471,0};
Point(15)={-1.802524,1.263442,0};
Point(16)={-1.521763,1.182470,0};
Point(17)={-1.261841,1.118833,0};
Point(18)={-1.022890,1.071132,0};
Point(19)={-0.804744,1.037504,0};
Point(20)={-0.606929,1.015842,0};
Point(21)={-0.428671,1.004021,0};
Point(22)={-0.268948,1.000022,0};
Point(23)={-0.126539,0.999960,0};
Point(24)={0.000000,0.999960,0};
Spline(0)={0:24};
Point(25)={-12.000000,4.250000,0};
Point(26)={-11.072160,4.250000,0};
Point(27)={-10.195415,4.250000,0};
Point(28)={-9.356321,4.250000,0};
Point(29)={-8.542011,4.250000,0};
Point(30)={-7.740000,4.250000,0};
Point(31)={-6.807604,4.250000,0};
Point(32)={-5.591946,4.178856,0};
Point(33)={-4.340811,3.536167,0};
Point(34)={-3.120568,2.831660,0};
Point(35)={-2.066436,2.223057,0};
Point(36)={-1.263920,1.759724,0};
Point(37)={-0.709796,1.439801,0};
Point(38)={-0.352487,1.233509,0};
Point(39)={-0.132151,1.106297,0};
Point(40)={-0.000000,1.022424,0};
Point(41)={0.000000,0.999960,0};
Spline(1)={25:40};
Line(2)={40,41};
Point(42)={-12.000000,-50.000000,-50.000000,1.0};
Point(43)={-12.000000,-50.000000,50.000000,1.0};
Point(44)={-12.000000,50.000000,50.000000,1.0};
Point(45)={-12.000000,50.000000,-50.000000,1.0};
Point(50)={-7.740000,-0.500000,-0.500000,1.0};
Point(51)={-7.740000,-0.500000,0.500000,1.0};
Point(52)={-7.740000,0.500000,0.500000,1.0};
Point(53)={-7.740000,0.500000,-0.500000,1.0};
//+
Rotate {{1, 0, 0}, {0, 0, 0}, -Pi*0.75} {
Line{1}; Line{0}; Line{2};
}
//+
Line(3) = {50, 51};
Line(4) = {51, 52};
Line(5) = {52, 53};
Line(6) = {53, 50};
Line Loop(1) = {4, 5, 6, 3};
Plane Surface(1) = {1};
Extrude {7.74, 0, 0} {
Surface{1};
}
//+
Extrude {{1, 0, 0}, {0, 0, 0}, Pi*0.5} {
Line{0,1};
}
//+
Extrude {{1, 0, 0}, {0, 0, 0}, Pi/2} {
Line{29}; Line{33};
}
//+
Extrude {{1, 0, 0}, {0, 0, 0}, Pi/2} {
Line{41}; Line{37};
}
//+
Extrude {{1, 0, 0}, {0, 0, 0}, Pi/2} {
Line{45}; Line{49};
}
//+
Line(60) = {88, 110};
//+
Line(61) = {139, 155};
//+
Line(62) = {196, 171};
//+
Line(63) = {63, 24};
//+
Line(64) = {54, 196};
//+
Line(65) = {139, 55};
//+
Line(66) = {59, 88};
//+
Line(67) = {44, 45};
//+
Line(68) = {45, 42};
//+
Line(69) = {42, 43};
//+
Line(70) = {43, 44};
//+
Line(71) = {140, 44};
//+
Line(72) = {156, 43};
//+
Line(73) = {25, 42};
//+
Line(74) = {95, 45};
//+
Line(75) = {51, 172};
//+
Line(76) = {50, 0};
//+
Line(77) = {53, 64};
//+
Line(78) = {52, 115};
//+ OutSqa Inlet
Line Loop(2) = {70, -71, 46, 72};
Plane Surface(61) = {2};
Line Loop(3) = {72, -69, -73, -54};
Plane Surface(62) = {3};
Line Loop(4) = {73, -68, -74, -34};
Plane Surface(63) = {4};
Line Loop(5) = {74, -67, -71, -42};
Plane Surface(64) = {5};
//+ NozIn to InSqa
Line Loop(6) = {10, 63, 31, -66};
Plane Surface(65) = {6};
Line Loop(7) = {11, 64, 59, -63};
Plane Surface(66) = {7};
Line Loop(8) = {8, -65, 51, -64};
Plane Surface(67) = {8};
Line Loop(9) = {65, 9, 66, 39};
Plane Surface(68) = {9};
//+ NozExit Wall
Line Loop(10) = {31, 60, -35, 2};
Plane Surface(69) = {10};
Line Loop(11) = {2, -59, 62, 55};
Plane Surface(70) = {11};
Line Loop(12) = {62, -47, -61, 51};
Plane Surface(71) = {12};
Line Loop(13) = {61, -43, -60, 39};
Plane Surface(72) = {13};
//+ Extrude OutSqa Inlet
Extrude {12, 0, 0} {
Surface{64}; Surface{61}; Surface{62}; Surface{63};
}
//+
Line(150) = {227, 40};
//+
Line(151) = {110, 197};
//+
Line(152) = {206, 155};
//+
Line(153) = {171, 217};
//+ NozExit Wall to MidLine
Line Loop(14) = {43, -152, 83, -151};
Plane Surface(167) = {14};
Line Loop(15) = {152, 47, 153, -104};
Plane Surface(168) = {15};
Line Loop(16) = {153, -127, 150, -55};
Plane Surface(169) = {16};
Line Loop(17) = {35, 151, 149, 150};
Plane Surface(170) = {17};
//+ Connect NozIn to InSqa
Line Loop(18) = {49, -64, -13, 75};
Plane Surface(171) = {18};
Line Loop(19) = {37, 65, -14, 78};
Plane Surface(172) = {19};
Line Loop(20) = {29, -66, -18, 77};
Plane Surface(173) = {20};
Line Loop(21) = {0, -63, -22, 76};
Plane Surface(174) = {21};
//+ Extrude MidSection to EndSection
Extrude {80, 0, 0} {
Surface{166}; Surface{100}; Surface{122}; Surface{144}; Surface{170}; Surface{169}; Surface{167}; Surface{168}; Surface{28}; Surface{65}; Surface{68}; Surface{67}; Surface{66}; Surface{69,70,71,72};
}
//+
Line Loop(549) = {50, -75, 4, 78};
//+
Plane Surface(550) = {549};
//+
Line Loop(551) = {78, -38, -77, -5};
//+
Plane Surface(552) = {551};
//+
Line Loop(553) = {77, -30, -76, -6};
//+
Plane Surface(554) = {553};
//+
Line Loop(555) = {3, 75, 58, -76};
//+
Plane Surface(556) = {555};
/* Delete { */
/* Line{0, 1}; */
/* } */
//+ NozIn Cylinder
Surface Loop(557) = {550, 52, 15, 67, 171, 172};
Volume(558) = {557};
Surface Loop(559) = {556, 60, 27, 66, 171, 174};
Volume(560) = {559};
Surface Loop(561) = {552, 40, 173, 172, 19, 68};
Volume(562) = {561};
Surface Loop(563) = {554, 32, 173, 174, 23, 65};
Volume(564) = {563};
//+
Line Loop(565) = {94, 152, -41};
//+
Plane Surface(566) = {565};
//+
Line Loop(567) = {85, -151, -33};
//+
Plane Surface(568) = {567};
//+
Line Loop(569) = {45, 153, -116};
//+
Plane Surface(570) = {569};
//+
Line Loop(571) = {138, 150, -1};
//+
Plane Surface(572) = {571};
/* Printf("%g,%g",newsl,newv); */
// NozOut to OutSqa
Surface Loop(573) = {117,168,48,566,570};
Volume(574) = {573};
Surface Loop(575) = {99,167,44,566,568};
Volume(576) = {575};
Surface Loop(577) = {165,170,36,568,572};
Volume(578) = {577};
Surface Loop(579) = {143,169,56,570,572};
Volume(580) = {579};
// InSqa 1: Volumn 1
// OutSqa 1: Volumn 2,3,4,5
// OutSqa 2: Volumn 6,7,8,9
// NozOut Cylinder 2: Volumn 10,11,12,13
// InSqa 2: Volumn 14
// InCylinder 2: Volumn 15,16,17,18
// NozWall 2: Volumn 19,20,21,22
//+ Inner Square and corresponding circles and outer square
Transfinite Line {3,4,5,6,8,9,10,11,352,353,354,355,
420,331,442,289,376,264,399,308,222,245,179,201,220,243,177,199,
51,47,59,55,31,35,39,43,104,127,149,83,102,125,147,81,
46,54,34,42,67,68,69,70,
30,38,50,58} = 4 Using Progression 1;
//+ Radial Inner
Transfinite Line {75,76,77,78,63,64,65,66,375,377,396,421} = 4 Using Progression 1;
//+ Radial Wall
Transfinite Line {2,60,61,388,463,465,486,508,
62} = 4 Using Progression 1;
//+ Radial Mid
Transfinite Line {150,151,152,153,265,267,286,309} = 4 Using Progression 1;
//+ Radial Outer
Transfinite Line {71,72,73,74,80,-82,105,-126,176,-178,-200,223} = 4 Using Progression 1.5;
//+ Streamwise: Outer Square 1
Transfinite Line {86,90,107,134,116,138,85,94,1,45,33,41} = 4 Using Progression 1;
//+ Streamwise: Outer Square 2
Transfinite Line {186,208,182,225,182,212,190,234,314,401,358,269,384,366,270,388,362,291,432,357,
181} = 4 Using Progression 1;
//+ Streamwise: Inner Square
Transfinite Line {0,22,49,13,37,14,29,18,
57} = 4 Using Progression 1;
//+
Transfinite Surface "*";
//+
Recombine Surface "*";
//+
Transfinite Volume "*";