Skip to content
Snippets Groups Projects
Commit c110dc69 authored by Christophe Geuzaine's avatar Christophe Geuzaine
Browse files

extrusion benchmarks

parent e2c43cbd
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
benchmarks/3d/stator-newxtr-rec.geo 0 → 100644
+ 249
− 0
View file @ c110dc69
lc =0.001;
s=1;
h0 = 2.0e-3; //3.e-3;
t1 = 1.0e-3; //3.9e-3;
t2 = 1.0e-3;
Wd = 6.e-3;
R = 27.5e-3;
r = 14.0e-3;
nArc = 15;
nSeg1 =2 ; // Side of Stator and Tooth (Vertical)
nSeg2 =2 ; // Piezo
nSeg3 =6;//5 ; // Long radius
nSeg4 =4;//3 ; // Wd
Point(1) = {0., 0., 0., lc} ;
Point(2) = {r, 0., 0., lc} ;
Point(3) = {(R-Wd/2), 0., 0., lc} ;
Point(4) = {(R+Wd/2), 0., 0., lc} ;
Point(5) = {0., 0., t1, lc} ;
Point(6) = {r, 0., t1, lc} ;
Point(7) = {(R-Wd/2), 0., t1, lc} ;
Point(8) = {(R+Wd/2), 0., t1, lc} ;
Point(9) = {0., 0., (t1+h0), lc};
Point(10)={(R-Wd/2), 0., (t1+h0), lc};
Point(11)={(R+Wd/2), 0., (t1+h0), lc};
Point(12) = {0., 0., -t2, lc} ;
Point(13) = {r, 0., -t2, lc} ;
Point(14) = {(R-Wd/2), 0., -t2, lc} ;
Point(15) = {(R+Wd/2), 0., -t2, lc} ;
Line(1) = {2,3} ;
Line(2) = {3,4} ;
Line(3) = {6,7} ;
Line(4) = {7,8} ;
Line(5) = {10,11} ;
Line(6) = {13,14} ;
Line(7) = {14,15} ;
Line(8) = {2,6} ;
Line(9) = {3,7} ;
Line(10) = {4,8} ;
Line(11) = {7,10} ;
Line(12) = {8,11} ;
Line(13) = {13,2} ;
Line(14) = {14,3} ;
Line(15) = {15,4} ;
Line Loop(1) = {3,-9,-1,8}; //surface centrale
Plane Surface(2) = {1};
Line Loop(3) = {4,-10,-2,9};
Plane Surface(4) = {3};
Line Loop(5) = {5,-12,-4,11};
Plane Surface(6) = {5};
Line Loop(7) = {1,-14,-6,13};
Plane Surface(8) = {7};
Line Loop(9) = {2,-15,-7,14};
Plane Surface(10) = {9};
// 1st gap
Extrude Surface {2, {0,0,1}, {0,0,0}, Pi/180}{
Recombine ; Layers { {1}, {9000+100}, {1} } ;
};
Extrude Surface {4, {0,0,1}, {0,0,0}, Pi/180}{
Recombine ; Layers { {1}, {9004+100}, {1} } ;
};
Extrude Surface {6, {0,0,1}, {0,0,0}, Pi/180}{
Recombine ; Layers { {1}, {9008+100}, {1} } ;
};
Extrude Surface {8, {0,0,1}, {0,0,0}, Pi/180}{
Recombine ; Layers { {1}, {9012+100}, {1} } ;
};
Extrude Surface {10, {0,0,1}, {0,0,0}, Pi/180}{
Recombine ; Layers { {1}, {9016+100}, {1} } ;
};
// 1st piezo
Extrude Surface {37, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Recombine ; Layers { {1,1,1,nArc,1,1,1}, {9000,9000,9000,9000,9000,9000,9000 },{0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {59, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Recombine ; Layers { {1,1,1,nArc,1,1,1}, {9004,9004,9004,9004,9004,9004,9004 },{0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1 } } ;
};
Extrude Surface {81, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Recombine ; Layers { {1,1,1,nArc,1,1,1}, {9008,9008,9008,9008,9008,9008,9008},{0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1 } } ;
};
Extrude Surface {103, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Recombine ; Layers { {1,1,1,nArc,1,1,1}, {9012,9012,9012,9012,9012,9012,9012 },{0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1 } } ;
};
Extrude Surface {125, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Recombine ; Layers { {1,1,1,nArc,1,1,1}, {9016,9016,9016,9016,9016,9016,9016 },{0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1 } } ;
};
// 2nd gap
Extrude Surface {147, {0,0,1}, {0,0,0}, Pi/180}{
Recombine ; Layers { {1}, {9001+200}, {1} } ;
};
Extrude Surface {169, {0,0,1}, {0,0,0}, Pi/180}{
Recombine ; Layers { {1}, {9005+200}, {1} } ;
};
Extrude Surface {191, {0,0,1}, {0,0,0}, Pi/180}{
Recombine ; Layers { {1}, {9009+200}, {1} } ;
};
Extrude Surface {213, {0,0,1}, {0,0,0}, Pi/180}{
Recombine ; Layers { {1}, {9013+200}, {1} } ;
};
Extrude Surface {235, {0,0,1}, {0,0,0}, Pi/180}{
Recombine ; Layers { {1}, {9017+200}, {1} } ;
};
// 2nd pizeo
Extrude Surface {257, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Recombine ; Layers { {1,1,1,nArc,1,1,1}, {9001,9001,9001,9001,9001,9001,9001,9001}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {279, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Recombine ; Layers { {1,1,1,nArc,1,1,1}, {9005, 9005,9005,9005, 9005,9005,9005}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {301, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Recombine ; Layers { {1,1,1,nArc,1,1,1}, {9009, 9009,9009,9009, 9009,9009,9009}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {323, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Recombine ; Layers { {1,1,1,nArc,1,1,1}, {9013, 9013,9013,9013, 9013,9013,9013}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {345, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Recombine ; Layers { {1,1,1,nArc,1,1,1}, {9017, 9017,9017,9017, 9017,9017,9017}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
// 3rd gap
Extrude Surface {257+110, {0,0,1}, {0,0,0}, Pi/180}{
Recombine ; Layers { {1}, {9002+300}, {1} } ;
};
Extrude Surface {279+110, {0,0,1}, {0,0,0}, Pi/180}{
Recombine ; Layers { {1}, {9006+300}, {1} } ;
};
Extrude Surface {301+110, {0,0,1}, {0,0,0}, Pi/180}{
Recombine ; Layers { {1}, {9010+300}, {1} } ;
};
Extrude Surface {323+110, {0,0,1}, {0,0,0}, Pi/180}{
Recombine ; Layers { {1}, {9014+300}, {1} } ;
};
Extrude Surface {345+110, {0,0,1}, {0,0,0}, Pi/180}{
Recombine ; Layers { {1}, {9018+300}, {1} } ;
};
// 3rd piezo
Extrude Surface {257+220, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Recombine ; Layers { {1,1,1,nArc,1,1,1}, {9002, 9002,9002,9002, 9002,9002,9002 }, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {279+220, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Recombine ; Layers { {1,1,1,nArc,1,1,1}, {9006, 9006,9006,9006, 9006,9006,9006}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {301+220, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Recombine ; Layers { {1,1,1,nArc,1,1,1}, {9010, 9010,9010,9010, 9010,9010,9010}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {323+220, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Recombine ; Layers { {1,1,1,nArc,1,1,1}, {9014, 9014,9014,9014, 9014,9014,9014}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {345+220, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Recombine ; Layers { {1,1,1,nArc,1,1,1}, {9018, 9018,9018,9018, 9018,9018,9018}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
// 4th gap
Extrude Surface {257+330, {0,0,1}, {0,0,0}, Pi/180}{
Recombine ; Layers { {1}, {9003+400}, {1} } ;
};
Extrude Surface {279+330, {0,0,1}, {0,0,0}, Pi/180}{
Recombine ; Layers { {1}, {9007+400}, {1} } ;
};
Extrude Surface {301+330, {0,0,1}, {0,0,0}, Pi/180}{
Recombine ; Layers { {1}, {9011+400}, {1} } ;
};
Extrude Surface {323+330, {0,0,1}, {0,0,0}, Pi/180}{
Recombine ; Layers { {1}, {9015+400}, {1} } ;
};
Extrude Surface {345+330, {0,0,1}, {0,0,0}, Pi/180}{
Recombine ; Layers { {1}, {9019+400}, {1} } ;
};
// 4th pizeo
Extrude Surface {257+440, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Recombine ; Layers { {1,1,1,nArc,1,1,1 }, {9003,9003,9003,9003,9003,9003,9003},{0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {279+440, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Recombine ; Layers { {1,1,1,nArc,1,1,1}, {9007,9007,9007,9007,9007,9007,9007}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {301+440, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Recombine ; Layers { {1,1,1,nArc,1,1,1}, {9011,9011,9011,9011,9011,9011,9011}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1}} ;
};
Extrude Surface {323+440, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Recombine ; Layers { {1,1,1,nArc,1,1,1}, {9015,9015,9015,9015,9015,9015,9015}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {345+440, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Recombine ; Layers { {1,1,1,nArc,1,1,1}, {9019,9019,9019,9019,9019,9019,9019}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
PZTA = 1000;
PZTB = 1500;
Stator = 2000;
Electrode1 = 1001;
Electrode2 = 1002;
Electrode3 = 1003;
Electrode4 = 1004;
Electrode_Upper1 = 1005;
Electrode_Upper2 = 1006;
Electrode_Upper3 = 1007;
Electrode_Upper4 = 1008;
Inner_Cylinder = 1009;
Point99 = 1010;
Point101 = 1011;
Cover = 1012;
Physical Volume (PZTA) = {9012,9013,9016,9017 , 9012+100,9013+200,9016+100,9017+200 } ;
Physical Volume (PZTB) = {9014,9015,9018,9019 , 9014+300,9015+400,9018+300,9019+400 } ;
Physical Volume (Stator) = {9000,9001,9002,9003, 9004,9005,9006,9007, 9008,9009,9010,9011,
9000+100,9001+200,9002+300,9003+400,
9004+100,9005+200,9006+300,9007+400,
9008+100,9009+200,9010+300,9011+400 } ;
Physical Surface (Electrode1) = {208,230} ;
Physical Surface (Electrode2) = {428,450} ;
Physical Surface (Electrode3) = {648,670} ;
Physical Surface (Electrode4) = {861,882} ;
Physical Surface (Electrode_Upper1) = {142,164} ;
Physical Surface (Electrode_Upper2) = {362,384} ;
Physical Surface (Electrode_Upper3) = {582,604} ;
Physical Surface (Electrode_Upper4) = {802,823} ;
Physical Surface (Inner_Cylinder) = {36,146,256,366, 476,586,696,806, 102,212,322,432,542,652,762,865} ;
Physical Point (Point99) = {201};
Physical Point (Point101) = {203};
benchmarks/3d/stator-newxtr.geo 0 → 100644
+ 249
− 0
View file @ c110dc69
lc =0.001;
s=1;
h0 = 2.0e-3; //3.e-3;
t1 = 1.0e-3; //3.9e-3;
t2 = 1.0e-3;
Wd = 6.e-3;
R = 27.5e-3;
r = 14.0e-3;
nArc = 15;
nSeg1 =2 ; // Side of Stator and Tooth (Vertical)
nSeg2 =2 ; // Piezo
nSeg3 =6;//5 ; // Long radius
nSeg4 =4;//3 ; // Wd
Point(1) = {0., 0., 0., lc} ;
Point(2) = {r, 0., 0., lc} ;
Point(3) = {(R-Wd/2), 0., 0., lc} ;
Point(4) = {(R+Wd/2), 0., 0., lc} ;
Point(5) = {0., 0., t1, lc} ;
Point(6) = {r, 0., t1, lc} ;
Point(7) = {(R-Wd/2), 0., t1, lc} ;
Point(8) = {(R+Wd/2), 0., t1, lc} ;
Point(9) = {0., 0., (t1+h0), lc};
Point(10)={(R-Wd/2), 0., (t1+h0), lc};
Point(11)={(R+Wd/2), 0., (t1+h0), lc};
Point(12) = {0., 0., -t2, lc} ;
Point(13) = {r, 0., -t2, lc} ;
Point(14) = {(R-Wd/2), 0., -t2, lc} ;
Point(15) = {(R+Wd/2), 0., -t2, lc} ;
Line(1) = {2,3} ;
Line(2) = {3,4} ;
Line(3) = {6,7} ;
Line(4) = {7,8} ;
Line(5) = {10,11} ;
Line(6) = {13,14} ;
Line(7) = {14,15} ;
Line(8) = {2,6} ;
Line(9) = {3,7} ;
Line(10) = {4,8} ;
Line(11) = {7,10} ;
Line(12) = {8,11} ;
Line(13) = {13,2} ;
Line(14) = {14,3} ;
Line(15) = {15,4} ;
Line Loop(1) = {3,-9,-1,8}; //surface centrale
Plane Surface(2) = {1};
Line Loop(3) = {4,-10,-2,9};
Plane Surface(4) = {3};
Line Loop(5) = {5,-12,-4,11};
Plane Surface(6) = {5};
Line Loop(7) = {1,-14,-6,13};
Plane Surface(8) = {7};
Line Loop(9) = {2,-15,-7,14};
Plane Surface(10) = {9};
// 1st gap
Extrude Surface {2, {0,0,1}, {0,0,0}, Pi/180}{
Layers { {1}, {9000+100}, {1} } ;
};
Extrude Surface {4, {0,0,1}, {0,0,0}, Pi/180}{
Layers { {1}, {9004+100}, {1} } ;
};
Extrude Surface {6, {0,0,1}, {0,0,0}, Pi/180}{
Layers { {1}, {9008+100}, {1} } ;
};
Extrude Surface {8, {0,0,1}, {0,0,0}, Pi/180}{
Layers { {1}, {9012+100}, {1} } ;
};
Extrude Surface {10, {0,0,1}, {0,0,0}, Pi/180}{
Layers { {1}, {9016+100}, {1} } ;
};
// 1st piezo
Extrude Surface {37, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Layers { {1,1,1,nArc,1,1,1}, {9000,9000,9000,9000,9000,9000,9000 },{0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {59, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Layers { {1,1,1,nArc,1,1,1}, {9004,9004,9004,9004,9004,9004,9004 },{0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1 } } ;
};
Extrude Surface {81, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Layers { {1,1,1,nArc,1,1,1}, {9008,9008,9008,9008,9008,9008,9008},{0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1 } } ;
};
Extrude Surface {103, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Layers { {1,1,1,nArc,1,1,1}, {9012,9012,9012,9012,9012,9012,9012 },{0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1 } } ;
};
Extrude Surface {125, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Layers { {1,1,1,nArc,1,1,1}, {9016,9016,9016,9016,9016,9016,9016 },{0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1 } } ;
};
// 2nd gap
Extrude Surface {147, {0,0,1}, {0,0,0}, Pi/180}{
Layers { {1}, {9001+200}, {1} } ;
};
Extrude Surface {169, {0,0,1}, {0,0,0}, Pi/180}{
Layers { {1}, {9005+200}, {1} } ;
};
Extrude Surface {191, {0,0,1}, {0,0,0}, Pi/180}{
Layers { {1}, {9009+200}, {1} } ;
};
Extrude Surface {213, {0,0,1}, {0,0,0}, Pi/180}{
Layers { {1}, {9013+200}, {1} } ;
};
Extrude Surface {235, {0,0,1}, {0,0,0}, Pi/180}{
Layers { {1}, {9017+200}, {1} } ;
};
// 2nd pizeo
Extrude Surface {257, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Layers { {1,1,1,nArc,1,1,1}, {9001,9001,9001,9001,9001,9001,9001,9001}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {279, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Layers { {1,1,1,nArc,1,1,1}, {9005, 9005,9005,9005, 9005,9005,9005}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {301, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Layers { {1,1,1,nArc,1,1,1}, {9009, 9009,9009,9009, 9009,9009,9009}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {323, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Layers { {1,1,1,nArc,1,1,1}, {9013, 9013,9013,9013, 9013,9013,9013}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {345, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Layers { {1,1,1,nArc,1,1,1}, {9017, 9017,9017,9017, 9017,9017,9017}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
// 3rd gap
Extrude Surface {257+110, {0,0,1}, {0,0,0}, Pi/180}{
Layers { {1}, {9002+300}, {1} } ;
};
Extrude Surface {279+110, {0,0,1}, {0,0,0}, Pi/180}{
Layers { {1}, {9006+300}, {1} } ;
};
Extrude Surface {301+110, {0,0,1}, {0,0,0}, Pi/180}{
Layers { {1}, {9010+300}, {1} } ;
};
Extrude Surface {323+110, {0,0,1}, {0,0,0}, Pi/180}{
Layers { {1}, {9014+300}, {1} } ;
};
Extrude Surface {345+110, {0,0,1}, {0,0,0}, Pi/180}{
Layers { {1}, {9018+300}, {1} } ;
};
// 3rd piezo
Extrude Surface {257+220, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Layers { {1,1,1,nArc,1,1,1}, {9002, 9002,9002,9002, 9002,9002,9002 }, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {279+220, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Layers { {1,1,1,nArc,1,1,1}, {9006, 9006,9006,9006, 9006,9006,9006}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {301+220, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Layers { {1,1,1,nArc,1,1,1}, {9010, 9010,9010,9010, 9010,9010,9010}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {323+220, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Layers { {1,1,1,nArc,1,1,1}, {9014, 9014,9014,9014, 9014,9014,9014}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {345+220, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Layers { {1,1,1,nArc,1,1,1}, {9018, 9018,9018,9018, 9018,9018,9018}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
// 4th gap
Extrude Surface {257+330, {0,0,1}, {0,0,0}, Pi/180}{
Layers { {1}, {9003+400}, {1} } ;
};
Extrude Surface {279+330, {0,0,1}, {0,0,0}, Pi/180}{
Layers { {1}, {9007+400}, {1} } ;
};
Extrude Surface {301+330, {0,0,1}, {0,0,0}, Pi/180}{
Layers { {1}, {9011+400}, {1} } ;
};
Extrude Surface {323+330, {0,0,1}, {0,0,0}, Pi/180}{
Layers { {1}, {9015+400}, {1} } ;
};
Extrude Surface {345+330, {0,0,1}, {0,0,0}, Pi/180}{
Layers { {1}, {9019+400}, {1} } ;
};
// 4th pizeo
Extrude Surface {257+440, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Layers { {1,1,1,nArc,1,1,1 }, {9003,9003,9003,9003,9003,9003,9003},{0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {279+440, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Layers { {1,1,1,nArc,1,1,1}, {9007,9007,9007,9007,9007,9007,9007}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {301+440, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Layers { {1,1,1,nArc,1,1,1}, {9011,9011,9011,9011,9011,9011,9011}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1}} ;
};
Extrude Surface {323+440, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Layers { {1,1,1,nArc,1,1,1}, {9015,9015,9015,9015,9015,9015,9015}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
Extrude Surface {345+440, {0,0,1}, {0,0,0}, Pi/2-Pi/180}{
Layers { {1,1,1,nArc,1,1,1}, {9019,9019,9019,9019,9019,9019,9019}, {0.02, 0.05, 0.09, 0.91, 0.95, 0.98, 1} } ;
};
PZTA = 1000;
PZTB = 1500;
Stator = 2000;
Electrode1 = 1001;
Electrode2 = 1002;
Electrode3 = 1003;
Electrode4 = 1004;
Electrode_Upper1 = 1005;
Electrode_Upper2 = 1006;
Electrode_Upper3 = 1007;
Electrode_Upper4 = 1008;
Inner_Cylinder = 1009;
Point99 = 1010;
Point101 = 1011;
Cover = 1012;
Physical Volume (PZTA) = {9012,9013,9016,9017 , 9012+100,9013+200,9016+100,9017+200 } ;
Physical Volume (PZTB) = {9014,9015,9018,9019 , 9014+300,9015+400,9018+300,9019+400 } ;
Physical Volume (Stator) = {9000,9001,9002,9003, 9004,9005,9006,9007, 9008,9009,9010,9011,
9000+100,9001+200,9002+300,9003+400,
9004+100,9005+200,9006+300,9007+400,
9008+100,9009+200,9010+300,9011+400 } ;
Physical Surface (Electrode1) = {208,230} ;
Physical Surface (Electrode2) = {428,450} ;
Physical Surface (Electrode3) = {648,670} ;
Physical Surface (Electrode4) = {861,882} ;
Physical Surface (Electrode_Upper1) = {142,164} ;
Physical Surface (Electrode_Upper2) = {362,384} ;
Physical Surface (Electrode_Upper3) = {582,604} ;
Physical Surface (Electrode_Upper4) = {802,823} ;
Physical Surface (Inner_Cylinder) = {36,146,256,366, 476,586,696,806, 102,212,322,432,542,652,762,865} ;
Physical Point (Point99) = {201};
Physical Point (Point101) = {203};
benchmarks/3d/stator-oldxtr.geo 0 → 100644
+ 238
− 0
View file @ c110dc69
r = 7.0e-3;
Wd = 2.e-3;
r0 = 3.0e-3;
nbr_tooth = 64;
teta_tooth = 0.0520881;
nbr_gap = 4;
teta_gap = Pi/180/2;
lcext =r*teta_tooth;
lcint =0.0008;
Point(1) = {0, 0, 0, lcext };
p=0; lt=0;
For i In {0:nbr_tooth-1}
t = i*2*Pi/nbr_tooth + 2*Pi/nbr_tooth/2 ;
p1 = newp;
Point(p1) = {r*Cos(t-teta_tooth/2) , r*Sin(t-teta_tooth/2) , 0 , lcext };
p2 = newp; points[p] = p2; p++;
Point(p2) = {(r+Wd)*Cos(t-teta_tooth/2) , (r+Wd)*Sin(t-teta_tooth/2) , 0 , lcext };
p3 = newp;
Point(p3) = {r*Cos(t+teta_tooth/2) , r*Sin(t+teta_tooth/2) , 0 , lcext };
p4 = newp; points[p] = p4; p++;
Point(p4) = {(r+Wd)*Cos(t+teta_tooth/2) , (r+Wd)*Sin(t+teta_tooth/2) , 0 , lcext };
l1 = newreg; Line(l1) = {p2,p1};
l2 = newreg; Line(l2) = {p3,p4};
l3 = newreg; Circle(l3) = {p1,1,p3};
l4 = newreg; Circle(l4) = {p2,1,p4};
s1 = newreg; Line Loop(s1) = {-l1,-l2,-l3,l4};
Plane Surface(s1+1) = {s1};
surftooth[i] = s1+1;
looptooth[lt] = l1; lt++;
looptooth[lt] = l2; lt++;
looptooth[lt] = l3; lt++;
EndFor
points[p] = points[0];
g=0; g0=0;
For i In {0:nbr_gap-1}
t = i*2*Pi/nbr_gap;
p1 = newp; gap[g] = p1; g++;
Point(p1) = {(r+Wd)*Cos(t-teta_gap/2) , (r+Wd)*Sin(t-teta_gap/2) , 0 , lcext };
p2 = newp; gap[g] = p2; g++;
Point(p2) = {(r+Wd)*Cos(t+teta_gap/2) , (r+Wd)*Sin(t+teta_gap/2) , 0 , lcext };
p3 = newp; gap0[g0] = p3; g0++;
Point(p3) = {r0*Cos(t-teta_gap/2) , r0*Sin(t-teta_gap/2) , 0 , lcint };
p4 = newp; gap0[g0] = p4; g0++;
Point(p4) = {r0*Cos(t+teta_gap/2) , r0*Sin(t+teta_gap/2) , 0 , lcint };
EndFor
gap0[g0]=gap0[0];
gap[g]=gap[0];
j=1; g=0; le=0; lge=0;
For i In {1:p-2:2}
If( Modulo(j,nbr_tooth/nbr_gap) )
loopext[le] = newreg;
Circle(loopext[le]) = {points[i],1,points[i+1]};
le++;
EndIf
If( !Modulo(j,nbr_tooth/nbr_gap) )
loopgapext[lge] = newreg;
Circle(loopgapext[lge]) = {gap[g],1,gap[g+1]};
lge++;
Printf("ii=%g", i);
loopext[le] = newreg;
Circle(loopext[le]) = {gap[g+1], 1, points[i-(2*nbr_tooth/nbr_gap-1)]};
le++;
loopext[le] = newreg;
Circle(loopext[le]) = {points[i], 1, gap[g+2]};
le++;
g+=2;
EndIf
j++;
EndFor
loopext[le] = loopext[0];
li=0; lgi=0; lg=0;
For i In {0:g0-2:2}
loopgapint[lgi] = newreg;
Circle(loopgapint[lgi]) = {gap0[i],1,gap0[i+1]};
lgi++;
loopgap[lg] = newreg;
Line(loopgap[lg]) = {gap[i], gap0[i]};
lg++;
loopgap[lg] = newreg;
Line(loopgap[lg]) = {gap0[i+1] , gap[i+1]};
lg++;
loopint[li] = newreg;
Circle(loopint[li]) = {gap0[i+1],1,gap0[i+2]};
li++;
EndFor
loopgap[lg] = loopgap[0];
i_looptooth = 0;
i_loopext = 0;
For i In {1:nbr_gap}
l1 = loopgap[2*i-1];
l2 = loopgap[2*i];
For j In {i_looptooth : i_looptooth+nbr_tooth*3/nbr_gap-1}
Printf("tooth=%g", looptooth[j]);
EndFor
For j In {i_loopext : i_loopext+nbr_tooth/nbr_gap+1-1}
Printf("ext=%g", loopext[j]);
EndFor
Printf("int=%g", loopint[i-1]);
Printf("gap=%g,%g", l1, l2);
ll = newreg;
Line Loop(ll) = {
looptooth[{i_looptooth : i_looptooth+nbr_tooth*3/nbr_gap-1}] ,
loopext[{i_loopext : i_loopext+nbr_tooth/nbr_gap+1-1}] ,
-loopint[i-1],
l1, l2
};
Plane Surface(ll+1) = ll;
surfair[i-1] = ll+1;
ll = newreg;
Line Loop(ll) = {
-loopgap[2*i-2],
-loopgap[2*i-1],
-loopgapint[i-1],
loopgapext[i-1]
};
Plane Surface(ll+1) = ll;
surfgap[i-1] = ll+1;
i_looptooth += nbr_tooth*3/nbr_gap;
i_loopext += nbr_tooth/nbr_gap+1;
EndFor
// Old extrusion mesh generator: gmsh -extrude < newstator.par
//
//
// All geometrical entities are automatically numbered:
//
// volumes: V_ + layer * N_ + surf->num
// New XY surfaces: S_ + layer * N_ + surf->num
// perp. surfaces: L_ + layer * N_ + curve->num
// perp. lines: P_ + layer * N_ + point->Num
L_ = 100e6 ;
S_ = 200e6 ;
V_ = 300e6 ;
P_ = 400e6 ;
N_ = 1e6 ;
// Volumes:
For i In {0:nbr_tooth-1}
toothvol1[i] = V_ + 1*N_ + surftooth[i] ;
toothvol2[i] = V_ + 2*N_ + surftooth[i] ;
toothvol3[i] = V_ + 3*N_ + surftooth[i] ;
EndFor
For i In {0:nbr_gap-1}
corpsvol1[i] = V_ + 1*N_ + surfair[i];
corpsvol2[i] = V_ + 2*N_ + surfair[i];
corpsvol3[i] = V_ + 3*N_ + surfair[i];
gapvol1[i] = V_ + 1*N_ + surfgap[i];
gapvol2[i] = V_ + 2*N_ + surfgap[i];
gapvol3[i] = V_ + 3*N_ + surfgap[i];
EndFor
// 1, 2, 3, ... -> volmumes piezo
ii = 0;
For i In {1:nbr_gap}
Physical Volume(i) = { toothvol1[{ii:(ii+nbr_tooth/nbr_gap-1)}], corpsvol1[i-1], gapvol1[i-1] };
ii += nbr_tooth/nbr_gap;
EndFor
// 100 -> volume metal
Physical Volume(100) = {
toothvol2[{0:nbr_tooth-1}],
toothvol3[{0:nbr_tooth-1}],
corpsvol2[{0:nbr_gap-1}],
gapvol2[{0:nbr_gap-1}]
};
// Surfaces
For i In {0:nbr_tooth-1}
toothsur0[i] = surftooth[i] ;
toothsur1[i] = S_ + 1*N_ + surftooth[i] ;
EndFor
For i In {0:nbr_gap-1}
corpssur0[i] = surfair[i];
corpssur1[i] = S_ + 1*N_ + surfair[i];
gapsur1[i] = S_ + 1*N_ + surfgap[i];
EndFor
// 1000 -> surface plan de masse
Physical Surface(1000) = { toothsur1[] , corpssur1[] , gapsur1[] };
// 1001, 1002, ... -> surfaces electrodes
ii = 0;
For i In {1:nbr_gap}
Physical Surface(1000+i) = {
toothsur0[{ii:(ii+nbr_tooth/nbr_gap-1)}],
corpssur0[i-1]
};
ii += nbr_tooth/nbr_gap;
EndFor
j=0;
For i In {0:nbr_gap-1}
cylint1[j] = L_ + 1*N_ + loopint[i];
cylint1[j+1] = L_ + 1*N_ + loopgapint[i];
cylint2[j] = L_ + 2*N_ + loopint[i];
cylint2[j+1] = L_ + 2*N_ + loopgapint[i];
j+=2;
EndFor
// 1100 -> surface cylindre interieur
Physical Surface(1100) = {
cylint1[],
cylint2[]
};
benchmarks/3d/stator-oldxtr.par 0 → 100644
+ 10
− 0
View file @ c110dc69
3
1
0.001
1
1
0.001
1
1
0.001
1
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment