Dear gmsh list:<br><br>Hi, It is my first time to post here. Now I am trying to build a vessel model with some stenosis in the inner side.<br><br>following is my geo code<br><br>---<br>// Gmsh project created on Tue May 17 04:19:43 2011<br>
lc = 0.1;<br>Point(1) = {0, 0, 0, lc};<br>Point(2) = {1, 0, 0, lc};<br>Point(3) = {0, 1, 0, lc};<br>Point(4) = {-1, 0, 0, lc};<br>Point(5) = {0, -1, 0, lc};<br><br>Point(6) = {1.1, 0, 0, lc};<br>Point(7) = {0, 1.1, 0, lc};<br>
Point(8) = {-1.1, 0, 0, lc};<br>Point(9) = {0, -1.1, 0, lc};<br><br>Circle(1) = {5, 1, 4};<br>Circle(2) = {4, 1, 3};<br>Circle(3) = {3, 1, 2};<br>Circle(4) = {2, 1, 5};<br>Circle(5) = {9, 1, 8};<br>Circle(6) = {8, 1, 7};<br>
Circle(7) = {7, 1, 6};<br>Circle(8) = {6, 1, 9};<br><br>Line Loop(1) = {1, 2, 3, 4};<br>Plane Surface(1) = {1};<br><br>Line Loop(2) = {5, 6, 7, 8};<br>Plane Surface(2) = {2, 1};<br><br>vol1[] = Extrude {0, 0, 10} {<br>  Surface{1};<br>
};<br>Printf("vol1 = %g", vol1[1]);<br><br>vol2[] = Extrude {0, 0, 10} {<br>  Surface{2};<br>};<br>Printf("vol2 = %g", vol2[1]);<br><br>r = 0.2;<br>p0 = newp; Point(p0) = {0, 1.1, 5, lc};<br>p1 = newp; Point(p1) = {-r, 1.1, 5, lc};<br>
p2 = newp; Point(p2) = {r, 1.1, 5, lc};<br>p3 = newp; Point(p3) = {0, 1.1, 5+r, lc};<br>c1 = newreg; Circle(c1) = {p2, p0, p3};<br>c2 = newreg; Circle(c2) = {p3, p0, p1};<br>l1 = newreg; Line(l1) = {p1, p2};<br>lp0 = newreg;<br>
Line Loop(lp0) = {c1, c2, l1};<br>s0 = newreg;<br>Plane Surface(s0) = {lp0};<br><br>vol3[] = Extrude {{1, 0, 0}, {1.1, 1.1, 5}, Pi} {<br>  Surface{77};<br>};<br>Printf("vol3 = %g", vol3[1]);<br><br>Compound Volume(4) = {vol2[1], vol3[1]}; <br>
---<br><br>you may notice i use vol3 to simulate the stenosis, vol1 is the tunnel to let blood go through, vol2 is the wall of the vessel.<br><u>but I don't know how to compound vol3 with vol2, can anyone help me out? appreciate!!</u><br>
<br><u>another question: If I try to rotate a surface 360 degree to form a volume by using extrude command, why I can not see the volume in the window?</u><br><br>thank you for your help.<br><br>Alfred<br>