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

removed all examples -> redistributed in benchmarks

parent 8fb1d9ed
No related branches found
No related tags found
No related merge requests found
Point(101) = {0,0,0,1} ;
Point(102) = {0,1,1,1} ;
Point(100) = {1,1,0,1} ;
Point(103) = {1,0,1,1} ;
Line(2001) = {100, 101} ;
Line(2000) = {102, 103} ;
Point(1) = {0.0,0.0,0.0,1.0};
Point(2) = {1,0.0,0.0,1.0};
Point(3) = {1,0.0,1,1.0};
Point(4) = {1,2,1,1.0};
Point(5) = {0,2,0,1.0};
Point(6) = {1,2,0,1.0};
Line(1) = {4,6};
Line(2) = {6,5};
Line(3) = {5,4};
Line(4) = {2,1};
Line(5) = {1,3};
Line(6) = {3,2};
Line(7) = {2,6};
Line(8) = {5,1};
Line(9) = {3,4};
Line Loop(10) = {9,-3,8,5};
Plane Surface(11) = {10};
Line Loop(12) = {-7,-6,9,1};
Plane Surface(13) = {12};
Line Loop(14) = {-4,7,2,8};
Plane Surface(15) = {14};
Line Loop(16) = {6,4,5};
Plane Surface(17) = {16};
Line Loop(18) = {1,2,3};
Plane Surface(19) = {18};
Surface Loop(20) = {13,15,17,-11,-19};
Complex Volume(21) = {20};
Transfinite Line {1,2,3,4,5,6,7,8,9} = 3;
Transfinite Surface {17} = {1,3,2};
Transfinite Surface {15} = {1,2,6,5};
Transfinite Surface {11} = {1,3,4,5};
Transfinite Surface {13} = {3,2,6,4};
Transfinite Surface {19} = {5,4,6};
Transfinite Volume {21} = {1,3,2,5,4,6};
/*
Recombine Surface {17,19};
*/
Recombine Surface {11,13,15};
l = .3;
Point(1) = {-2,-1,0,l};
Point(2) = {-2,1,0,l};
Point(3) = {0,-1,0,l};
Point(4) = {0,1,0,l};
Point(5) = {-3,-2,0,l};
Point(6) = {-3,2,0,l};
Point(7) = {1,-2,0,l};
Point(8) = {1,2,0,l};
Line(1) = {6,8};
Line(2) = {8,7};
Line(3) = {7,5};
Line(4) = {5,6};
Line(5) = {2,1};
Line(6) = {1,3};
Line(7) = {3,4};
Line(8) = {4,2};
Line Loop(9) = {5,6,7,8};
Plane Surface(10) = {9};
Line Loop(11) = {2,3,4,1};
Plane Surface(12) = {11,9};
/* this will try to produce quadrangles for surfaces 10 and 12 */
Recombine Surface {10,12} ;
/* For surface 10, we force a structured mesh (transfinite meshing) */
/* We force the number of points on the lines 5,6,7,8 to be 10 and 5 */
Transfinite Line {5,7} = 10 ;
Transfinite Line {6,8} = 5 ;
/* we define transfinite interpolation on surface 10, determined by the
ordering of the for ***vertices*** 1,3,2,4 */
Transfinite Surface {10} = {1,3,4,2} ;
/*
result :
Surface 10 should be made entirely of quads
Surface 12 should be mixed tri/quad
comments :
Transfinite surfaces can only have 3 or 4 sides
to have a specified repartion of the nodes for transfinite lines :
Transfinite Line {X} = Y Using Power Z ;
Transfinite Line {X} = Y Using Bump Z ;
*/
Point(1) = {0,0,0,.25};
Point(2) = {0,1,0,.25};
Point(3) = {1,0,0,.25};
Point(4) = {1,1,0,.25};
Line(1) = {1,2};
Line(2) = {2,4};
Line(3) = {4,3};
Line(4) = {3,1};
Line Loop(5) = {4,1,2,3};
Plane Surface(6) = 5;
Extrude( 6, { 1.,0.,0.}, {0,-1,0.}, 1.57);
Extrude(24, {0,0,-1}, {-1,-1,0}, 1.57);
Coherence;
1
1
0.04
// Mudanca de algumas coordenadas 22/05/2001 - Liège.
/* Desenho do servo motor em 3D - 20/07/2000 */
/* Analise de um servomotor de imas permanentes - GRUCAD - julho de 2000 */
rR1 = 10.5e-03;
rR2 = 23.243e-03;
rR3 = 23.862e-03;
rR4 = 23.9e-03;
rR5 = 25.6e-03;
rS1 = 26.02e-03;
rS2 = 26.62e-03;
rS3 = 26.96e-03;
rS4 = 38.16e-03;
rS5 = 38.27e-03;
rS6 = 40.02e-03;
rS7 = 46.0e-03;
Gap = rS1-rR5;
rB1=rR5+Gap/3; rB1b=rB1;
//rB1=rR5+Gap/3.7; rB1b=rR5+Gap/2.3; //To see both different surfaces
rB2=rR5+Gap*2/3;
NbrDiv = 16;
A0 = 45.0;
/* Rotor position(deg):use 0 for the mesh and other angles for drawing*/
//A1 = 4.0001*45.0/NbrDiv;
A1 = 90. ;
deg2rad = Pi/180;
/* Mesh parameters */
s = 1.;
pR1 =(rR2-rR1)/10.*s;
pR2 =(rR2-rR1)/12.*s;
pS1 =(rS7-rS1)/7.*s;
pS2 =(rS7-rS1)/12.*s;
pS3 =(rS6-rS3)/10.*s;
// Uma das faces do servo motor
/* Rotor */
Point(1)={0,0,0,pR1};
A = 0.+A1; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(2)={rR1*cosA,rR1*sinA,0,pR1};
Point(3)={rR4*cosA,rR4*sinA,0,pR1};
Point(4)={rR5*cosA,rR5*sinA,0,pR1};
Point(5)={rB1*cosA,rB1*sinA,0,pR2};
A = 4.455+A1; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(6)={rR2*cosA,rR2*sinA,0,pR1};
A = 4.938+A1; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(7)={rR3*cosA,rR3*sinA,0,pR1};
A = 6.165+A1; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(8)={rR5*cosA,rR5*sinA,0,pR1};
A = 38.835+A1; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(9)={rR5*cosA,rR5*sinA,0,pR1};
A = 40.062+A1; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(10)={rR3*cosA,rR3*sinA,0,pR1};
A = 40.545+A1; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(11)={rR2*cosA,rR2*sinA,0,pR1};
A = 45.0+A1; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(12)={rR1*cosA,rR1*sinA,0,pR1};
Point(13)={rR4*cosA,rR4*sinA,0,pR1};
Point(14)={rR5*cosA,rR5*sinA,0,pR1};
Point(15)={rB1*cosA,rB1*sinA,0,pR2};
/* Gap */
Point(16)={rB1b,0,0,pS1};
Point(17)={rB2,0,0,pS1};
A = 45.0; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(18)={rB1b*cosA,rB1b*sinA,0,pS1};
Point(19)={rB2*cosA,rB2*sinA,0,pS1};
/* Stator */
Point(20)={rS1,0,0,pS1};
Point(21)={rS7,0,0,pS2};
A = 5.71+0; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(61)={rS7*cosA,rS7*sinA,0,pS3};
A = 5.71+15; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(62)={rS7*cosA,rS7*sinA,0,pS3};
A = 5.71+30; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(63)={rS7*cosA,rS7*sinA,0,pS3};
//A = 22.5; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
//Point(22)={rS7*cosA,rS7*sinA,0,pS2};
A = 45.0; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(23)={rS1*cosA,rS1*sinA,0,pS1};
Point(24)={rS7*cosA,rS7*sinA,0,pS2};
/* Pontos da ranhura 01 */
A2 = 0.0;
A = 2.77+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(25)={rS5*cosA,rS5*sinA,0,pS3};
A = 4.0+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(26)={rS3*cosA,rS3*sinA,0,pS3};
A = 5.52+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(27)={rS1*cosA,rS1*sinA,0,pS3};
A = 5.56+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(28)={rS2*cosA,rS2*sinA,0,pS3};
A = 5.56+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(29)={rS4*cosA,rS4*sinA,0,pS3};
A = 5.65+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(30)={rS6*cosA,rS6*sinA,0,pS3};
A = 9.35+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(31)={rS6*cosA,rS6*sinA,0,pS3};
A = 9.44+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(32)={rS4*cosA,rS4*sinA,0,pS3};
A = 9.44+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(33)={rS2*cosA,rS2*sinA,0,pS3};
A = 9.48+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(34)={rS1*cosA,rS1*sinA,0,pS3};
A = 11.0+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(35)={rS3*cosA,rS3*sinA,0,pS3};
A = 12.23+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(36)={rS5*cosA,rS5*sinA,0,pS3};
/* Pontos da ranhura 02 */
A2 = 15.0;
A = 2.77+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(37)={rS5*cosA,rS5*sinA,0,pS3};
A = 4.0+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(38)={rS3*cosA,rS3*sinA,0,pS3};
A = 5.52+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(39)={rS1*cosA,rS1*sinA,0,pS3};
A = 5.56+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(40)={rS2*cosA,rS2*sinA,0,pS3};
A = 5.56+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(41)={rS4*cosA,rS4*sinA,0,pS3};
A = 5.65+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(42)={rS6*cosA,rS6*sinA,0,pS3};
A = 9.35+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(43)={rS6*cosA,rS6*sinA,0,pS3};
A = 9.44+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(44)={rS4*cosA,rS4*sinA,0,pS3};
A = 9.44+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(45)={rS2*cosA,rS2*sinA,0,pS3};
A = 9.48+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(46)={rS1*cosA,rS1*sinA,0,pS3};
A = 11.0+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(47)={rS3*cosA,rS3*sinA,0,pS3};
A = 12.23+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(48)={rS5*cosA,rS5*sinA,0,pS3};
/* Pontos da ranhura 03 */
A2 = 30.0;
A = 2.77+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(49)={rS5*cosA,rS5*sinA,0,pS3};
A = 4.0+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(50)={rS3*cosA,rS3*sinA,0,pS3};
A = 5.52+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(51)={rS1*cosA,rS1*sinA,0,pS3};
A = 5.56+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(52)={rS2*cosA,rS2*sinA,0,pS3};
A = 5.56+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(53)={rS4*cosA,rS4*sinA,0,pS3};
A = 5.65+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(54)={rS6*cosA,rS6*sinA,0,pS3};
A = 9.35+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(55)={rS6*cosA,rS6*sinA,0,pS3};
A = 9.44+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(56)={rS4*cosA,rS4*sinA,0,pS3};
A = 9.44+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(57)={rS2*cosA,rS2*sinA,0,pS3};
A = 9.48+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(58)={rS1*cosA,rS1*sinA,0,pS3};
A = 11.0+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(59)={rS3*cosA,rS3*sinA,0,pS3};
A = 12.23+A2; sinA=Sin(A*deg2rad);cosA=Cos(A*deg2rad);
Point(60)={rS5*cosA,rS5*sinA,0,pS3};
// linhas e circulos do rotor
Line(1)={1,2};
Line(2)={2,3};
Line(3)={3,4};
Line(4)={4,5};
Line(5)={3,7};
Line(6)={6,7};
Line(7)={7,8};
Line(8)={6,11};
Line(9)={11,10};
Line(10)={10,9};
Line(11)={1,12};
Line(12)={12,13};
Line(13)={13,14};
Line(14)={14,15};
Line(15)={13,10};
Circle(16)={2,1,12};
Circle(17)={4,1,8};
Circle(18)={8,1,9};
Circle(19)={9,1,14};
Circle(20)={5,1,15};
/* linhas e circulos do entreferro */
Line(21)={16,17};
Line(22)={17,20};
Line(23)={18,19};
Line(24)={19,23};
Circle(25)={16,1,18};
Circle(26)={17,1,19};
/* linhas e circulos do estator */
Line(27)={20,21};
Line(28)={23,24};
Circle(29)={20,1,27};
Circle(30)={34,1,39};
Circle(31)={46,1,51};
Circle(32)={58,1,23};
Circle(33)={21,1,61};
Circle(34)={61,1,62};
Circle(89)={62,1,63};
Circle(90)={63,1,24};
Line(91)={30,61};
Line(92)={42,62};
Line(93)={54,63};
/* linhas e circulos da ranhura 01 */
Line(35)={27,28};
Line(36)={28,26};
Line(37)={26,25};
Circle(38)={25,29,30};
Line(39)={30,31};
Circle(40)={31,32,36};
Line(41)={36,35};
Line(42)={35,33};
Line(43)={33,34};
Line(44)={26,35};
/* linhas e circulos da ranhura 02 */
Line(45)={39,40};
Line(46)={40,38};
Line(47)={38,37};
Circle(48)={37,41,42};
Line(49)={42,43};
Circle(50)={43,44,48};
Line(51)={48,47};
Line(52)={47,45};
Line(53)={45,46};
Line(54)={38,47};
/* linhas e circulos da ranhura 03 */
Line(55)={51,52};
Line(56)={52,50};
Line(57)={50,49};
Circle(58)={49,53,54};
Line(59)={54,55};
Circle(60)={55,56,60};
Line(61)={60,59};
Line(62)={59,57};
Line(63)={57,58};
Line(64)={50,59};
/* superficies do rotor */
Line Loop(65) = {-11,1,16};
Plane Surface(66) = {65};
Line Loop(67) = {12,15,-9,-8,6,-5,-2,16};
Plane Surface(68) = {67};
Line Loop(69) = {7,-17,-3,5};
Plane Surface(70) = {69};
Line Loop(71) = {-18,-7,-6,8,9,10};
Plane Surface(72) = {71};
Line Loop(73) = {-19,-10,-15,13};
Plane Surface(74) = {73};
Line Loop(75) = {18,19,14,-20,-4,17};
Plane Surface(76) = {75};
/* superficies do entreferro (banda de movimento) */
Transfinite Line{20,25,26}=NbrDiv+1 Using Power 1.;
Transfinite Line{21,23}=2;
Line Loop(77) = {-23,-25,21,26};
Plane Surface(78) = {77};
Transfinite Surface{78}={16,17,19,18};
Recombine Surface{78}; // For Hexahedra in the Moving Band
/* superficies do estator */
Line Loop(79) = {24,-32,-63,-62,-64,-56,-55,-31,-53,-52,-54,-46,-45,-30,-43,-42,-44,-36,-35,-29,-22,26};
Plane Surface(80) = {79};
Line Loop(94) = {37,38,91,-33,-27,29,35,36};
Plane Surface(95) = {94};
Line Loop(96) = {42,43,30,45,46,47,48,92,-34,-91,39,40,41};
Plane Surface(97) = {96};
Line Loop(98) = {52,53,31,55,56,57,58,93,-89,-92,49,50,51};
Plane Surface(99) = {98};
Line Loop(100) = {62,63,32,28,-90,-93,59,60,61};
Plane Surface(101) = {100};
/*
Line Loop(81) = {37,38,39,40,41,42,43,30,45,46,47,48,49,50,51,52,53,31,55,56,57,58,59,60,61,62,63,32,28,-90,-89,-34,-33,-27,29,35,36};
Plane Surface(82) = {81};
*/
/* superficie da ranhura 01 */
Line Loop(83) = {-41,-40,-39,-38,-37,44};
Plane Surface(84) = {83};
/* superficie da ranhura 02 */
Line Loop(85) = {-51,-50,-49,-48,-47,54};
Plane Surface(86) = {85};
/* superficie da ranhura 03 */
Line Loop(87) = {-61,-60,-59,-58,-57,64};
Plane Surface(88) = {87};
// Continue with 101 ...
// ----------------------------------------------------
// ----------------------------------------------------
L_ = 100e6 ;
S_ = 200e6 ;
V_ = 300e6 ;
P_ = 400e6 ;
N_ = 1e6 ;
Physical Volume(1301) = {V_+N_+68}; // Rotor
//Physical Volume(1302) = {V_+N_+95, V_+N_+97, V_+N_+99, V_+N_+101}; // Stator
Physical Volume(1320) = {V_+N_+95}; // StatorSup0
Physical Volume(1321) = {V_+N_+97}; // StatorSup1
Physical Volume(1322) = {V_+N_+99}; // StatorSup2
Physical Volume(1323) = {V_+N_+101}; // StatorSup3
Physical Volume(1300) = {V_+N_+66}; // AirInt
Physical Volume(1303) = {V_+N_+76}; // AirRotor
Physical Volume(1304) = {V_+N_+80}; // AirStator
Physical Volume(1305) = {V_+N_+78}; // AirBM
Physical Volume(1310) = {V_+N_+72}; // Magnet
Physical Volume(1311) = {V_+N_+70}; // Cale1
Physical Volume(1312) = {V_+N_+74}; // Cale2
Physical Volume(1401) = {V_+N_+84}; // Ind1
Physical Volume(1402) = {V_+N_+86}; // Ind2
Physical Volume(1403) = {V_+N_+88}; // Ind3
Physical Surface(2001) = {L_+N_+16}; // SurfInt
Physical Surface(2002) = {L_+N_+33,L_+N_+34,L_+N_+89,L_+N_+90}; // SurfExt
Physical Surface(1011) = {L_+N_+ 1, L_+N_+2, L_+N_+3, L_+N_+4, L_+N_+21,L_+N_+22, L_+N_+27}; // SurfA
Physical Surface(1012) = {L_+N_+11,L_+N_+12,L_+N_+13,L_+N_+14, L_+N_+23,L_+N_+24, L_+N_+28}; // SurfB
Physical Surface(1501) = {L_+N_+20}; // SurfR
Physical Surface(1502) = {L_+N_+25}; // SurfBM
Physical Line(2021) = {P_+N_+15}; // PointR
Physical Line(2022) = {P_+N_+18}; // PointBM
// SurfR (PointR) & SurfBM (PointBM) are lines (points) located at the same
// position, but with distinct nodes
// For h-formulation (cut ...)
Physical Surface(2210) = {L_+N_+37,L_+N_+38,/*L_+N_+39,*/L_+N_+40,L_+N_+41,L_+N_+44}; // SkinInd1
Physical Surface(2211) = {L_+N_+91}; // CutInd1
Physical Surface(2212) = {L_+N_+39}; // SkinInd1_2
Physical Surface(2220) = {L_+N_+47,L_+N_+48,/*L_+N_+49,*/L_+N_+50,L_+N_+51,L_+N_+54}; // SkinInd2
Physical Surface(2221) = {L_+N_+92}; // CutInd2
Physical Surface(2222) = {L_+N_+49}; // SkinInd2_2
Physical Surface(2230) = {L_+N_+57,L_+N_+58,/*L_+N_+59,*/L_+N_+60,L_+N_+61,L_+N_+64}; // SkinInd3
Physical Surface(2231) = {L_+N_+93}; // CutInd3
Physical Surface(2232) = {L_+N_+59}; // SkinInd3_2
Physical Point(2015) = {2}; // PointRef
i = 0*S_+0*N_ ;
Physical Surface(2300) =
{ i+68, i+95,i+97,i+99,i+101, i+66,i+76,i+80,i+78, i+72, i+70,i+74, i+84,i+86,i+88 }; // SurfZ0
i = S_+1*N_ ;
Physical Surface(2301) =
{ i+68, i+95,i+97,i+99,i+101, i+66,i+76,i+80,i+78, i+72, i+70,i+74, i+84,i+86,i+88 }; // SurfZ1
lc = 0.01;
dext = .05;
dint = .055;
Point(1) = {0,0,0,lc};
Point(2) = {dint,0,0,lc};
Point(3) = {0,dint,0,lc};
Point(4) = {0,0,dint,lc};
Point(5) = {dext,0,0,lc};
Point(6) = {0,dext,0,lc};
Point(7) = {0,0,dext,lc};
Circle(1) = {7,1,5};
Circle(2) = {4,1,2};
Circle(3) = {5,1,6};
Circle(4) = {2,1,3};
Circle(5) = {7,1,6};
Circle(6) = {4,1,3};
Line(7) = {1,5};
Line(8) = {5,2};
Line(9) = {1,6};
Line(10) = {6,3};
Line(11) = {1,7};
Line(12) = {7,4};
Line Loop(17) = {2,-8,-1,12};
Plane Surface(18) = {17};
Line Loop(19) = {-4,-8,3,10};
Plane Surface(20) = {19};
Line Loop(21) = {-6,-12,5,10};
Plane Surface(22) = {21};
Line Loop(23) = {-6,2,4};
Ruled Surface(24) = {23};
Line Loop(25) = {-5,1,3};
Ruled Surface(26) = {25};
Surface Loop(27) = {24,-22,-18,20,-26};
Complex Volume(28) = {27};
Transfinite Line {1,2,3,4,5,6} = 10 ;
Transfinite Line {-8,-10,-12} = 5 Using Power 1.6 ;
Transfinite Surface {22} = {4,7,6,3};
Transfinite Surface {20} = {3,2,5,6};
Transfinite Surface {18} = {4,2,5,7};
Transfinite Surface {24} = {3,4,2};
Transfinite Surface {26} = {6,7,5};
Recombine Surface {18,20,22};
/*
Recombine Surface {24,26};
*/
Transfinite Volume {28} = {3,4,2,6,7,5};
Include "../tutorial/view1.pos" ;
View[0].Name = "gloups.prout" ;
General.Trackball = 0 ;
For num In {1:2}
General.RotationX += 10 ;
Print StrCat( StrPrefix(Sprintf(View[0].Name)) , Sprintf("-%g.jpg", num) ) ;
Draw ;
EndFor
l= .3;
Point(1) = {-1.00000E+00,-1.00000E-00, 0.00000E+00,l};
Point(2) = {-1.00000E+00, 1.00000E-00, 0.00000E+00,l};
Point(3) = { 1.00000E+00, 1.00000E-00, 0.00000E+00,l};
Point(4) = { 1.00000E+00,-1.00000E-00, 0.00000E+00,l};
Point(5) = {-0.00000E+00,-0.00000E-00, 0.00000E+00,l};
Point(6) = {-0.00000E+00, 1.00000E-00, 0.00000E+00,l};
Point(7) = { 1.00000E+00, 0.00000E-00, 0.00000E+00,l};
Line(1) = Liste[2,6];
Line(2) = Liste[6,5];
Line(3) = Liste[5,7];
Line(4) = Liste[7,4];
Line(5) = Liste[4,1];
Line(6) = Liste[1,2];
Line Loop(7) = Liste[3,4,5,6,1,2];
Plane Surface(8) = Liste[7];
Extrude ( 8, { 0.00000E+00, 0.00000E+00, 1.00000E+00} );
Line(35) = Liste[6,3];
Line(36) = Liste[3,7];
Line Loop(37) = Liste[3,-36,-35,2];
Plane Surface(38) = Liste[37];
Extrude ( 38, { 0.00000E+00, 0.00000E+00,-1.00000E+00} );
Extrude ( 8, { 0.00000E+00, 0.00000E+00,-1.00000E+00} );
Surface Loop(80) = Liste[56,-79,-68,22,25,-71,-74,28,31,-76,-53,-38,19,34,33,-50];
Complex Volume(81) = Liste[80];
dT1 = 10.e-3 ;
dT2 = 10.e-3 ;
dT3 = 10.e-3 ;
lT1 = 10.e-3 ;
lT2 = 10.e-3 ;
lT3 = 10.e-3 ;
e = 0.65e-3/2. ;
pp = 4.99e-3 ;
p0 = pp ;
pT1 = pp ; pT2 = pp ; pT3 = pp ;
pT13 = pp ; pT23 = pp ;
Point(1) = { -dT3/2.,0,0,p0};
Point(2) = { -dT3/2.-lT1, 0, 0, pT1};
Point(3) = { -dT3/2.-lT1, dT1, 0, pT1};
Point(4) = { -dT3/2., dT1 , 0, pT13};
Point(5) = { -dT3/2., dT1+lT3, 0, pT3};
Point(6) = { dT3/2., dT1+lT3, 0, pT3};
Point(7) = { dT3/2., dT2 , 0, pT23};
Point(8) = { dT3/2.+lT2, dT2, 0, pT2};
Point(9) = { dT3/2.+lT2, 0, 0, pT2};
Point(10) = { dT3/2.,0,0,p0};
Line(1) = {1,2};
Line(2) = {2,3};
Line(3) = {3,4};
Line(4) = {4,1};
Line(5) = {1,10};
Line(6) = {10,7};
Line(7) = {7,4};
Line(8) = {4,5};
Line(9) = {5,6};
Line(10) = {6,7};
Line(11) = {7,8};
Line(12) = {8,9};
Line(13) = {9,10};
Line Loop(14) = {1,2,3,4};
Plane Surface(15) = {14};
Line Loop(16) = {11,12,13,6};
Plane Surface(17) = {16};
Line Loop(18) = {8,9,10,7};
Plane Surface(19) = {18};
Line Loop(20) = {5,6,7,4};
Plane Surface(21) = {20};
Translate({0,0,e}) {
Duplicata { Surface(15); }
}
Coherence;
Translate({0,0,e}) {
Duplicata { Surface(21); }
}
Coherence;
Translate({0,0,e}) {
Duplicata { Surface(17); }
}
Coherence;
Translate({0,0,e}) {
Duplicata { Surface(19); }
}
Coherence;
Line(46) = {2,12};
Line(47) = {3,14};
Line(48) = {1,11};
Line(49) = {4,16};
Line(50) = {5,26};
Line(51) = {6,28};
Line(52) = {7,20};
Line(53) = {10,18};
Line(54) = {9,24};
Line(55) = {8,22};
Line Loop(56) = {23,-47,-2,46};
Plane Surface(57) = {56};
Line Loop(58) = {22,-46,-1,48};
Plane Surface(59) = {58};
Line Loop(60) = {-49,-3,47,24};
Plane Surface(61) = {60};
Line Loop(62) = {-25,-49,4,48};
Plane Surface(63) = {62};
Line Loop(64) = {40,-50,-8,49};
Plane Surface(65) = {64};
Line Loop(66) = {41,-51,-9,50};
Plane Surface(67) = {66};
Line Loop(68) = {-42,-51,10,52};
Plane Surface(69) = {68};
Line Loop(70) = {-30,-52,7,49};
Plane Surface(71) = {70};
Line Loop(72) = {28,-53,-5,48};
Plane Surface(73) = {72};
Line Loop(74) = {29,-52,-6,53};
Plane Surface(75) = {74};
Line Loop(76) = {34,-55,-11,52};
Plane Surface(77) = {76};
Line Loop(78) = {36,-53,-13,54};
Plane Surface(79) = {78};
Line Loop(80) = {-35,-55,12,54};
Plane Surface(81) = {80};
Surface Loop(82) = {63,27,-59,-57,-61,-15};
Complex Volume(83) = {82};
Surface Loop(84) = {75,-39,77,-81,17,79};
Complex Volume(85) = {84};
Surface Loop(86) = {71,45,-65,-67,69,-19};
Complex Volume(87) = {86};
Surface Loop(88) = {71,33,-73,-75,-21,63};
Complex Volume(89) = {88};
/* selon e */
Transfinite Line {-46,-48,-53,-54,-47,-49,-52,-55,-50,-51}
= 8 Using Power 1.2 ;
/* vers les angles */
Transfinite Line {-3,-24,1,22, 11,34,-13,-36, 8,40,-10,-42}
= 11 Using Power 1.4 ;
/* selon y pour f1 et f2 */
Transfinite Line {-2,-23,4,25,-6,-29,12,35}
= 11 Using Bump -1.4 ;
/* selon x pour f3 */
Transfinite Line {9,41,7,30,5,28}
= 11 Using Bump -1.4 ;
Transfinite Surface {15} = {1,2,3,4} ;
Transfinite Surface {27} = {11,12,14,16} ;
Transfinite Surface {17} = {9,10,7,8} ;
Transfinite Surface {39} = {24,18,20,22} ;
Transfinite Surface {19} = {7,4,5,6} ;
Transfinite Surface {45} = {20,16,26,28} ;
Transfinite Surface {21} = {10,1,4,7} ;
Transfinite Surface {33} = {18,11,16,20} ;
Transfinite Surface {57} = {2,12,14,3} ;
Transfinite Surface {63} = {1,11,16,4} ;
Transfinite Surface {75} = {10,18,20,7} ;
Transfinite Surface {81} = {9,24,22,8} ;
Transfinite Surface {65} = {4,16,26,5} ;
Transfinite Surface {69} = {7,20,28,6} ;
Transfinite Surface {59} = {12,11,1,2} ;
Transfinite Surface {61} = {14,16,4,3} ;
Transfinite Surface {73} = {11,18,10,1} ;
Transfinite Surface {71} = {16,20,7,4} ;
Transfinite Surface {67} = {26,28,6,5} ;
Transfinite Surface {79} = {18,24,9,10} ;
Transfinite Surface {77} = {20,22,8,7} ;
Recombine Surface {15,27,17,39,19,45,21,33,57,63,75,81,59,61,73,71,67,79,77,65,69} ;
Transfinite Volume {83} = {1,2,12,11,4,3,14,16} ;
Transfinite Volume {85} = {9,10,18,24,8,7,20,22} ;
Transfinite Volume {87} = {7,4,16,20,6,5,26,28} ;
Transfinite Volume {89} = {10,1,11,18,7,4,16,20} ;
Physical Volume(101) = {83,85,87,89} ; /* Tooth */
Physical Line(2001) = { 2,47,23,46} ; /* LineFlux1 */
Physical Line(2002) = {55,12,54,35} ; /* LineFlux2 */
Physical Line(2003) = { 9,51,41,50} ; /* LineFlux3 */
Physical Surface(1001) = {57} ; /* Flux1 */
Physical Surface(1002) = {81} ; /* Flux2 */
Physical Surface(1003) = {67} ; /* Flux3 */
Physical Surface(1100) = {27,33,39, 59,73,79, 61,65,45,69,77} ; /* SkinTooth */
dT1 = 10.e-3 ;
dT2 = 10.e-3 ;
dT3 = 10.e-3 ;
lT1 = 10.e-3 ;
lT2 = 10.e-3 ;
lT3 = 10.e-3 ;
e = 0.65e-2/2. ;
pp = 4.99e-3 ;
p0 = pp ;
pT1 = pp ; pT2 = pp ; pT3 = pp ;
pT13 = pp ; pT23 = pp ;
Point(1) = { -dT3/2.,0,0,p0};
Point(2) = { -dT3/2.-lT1, 0, 0, pT1};
Point(3) = { -dT3/2.-lT1, dT1, 0, pT1};
Point(4) = { -dT3/2., dT1 , 0, pT13};
Point(5) = { -dT3/2., dT1+lT3, 0, pT3};
Point(6) = { dT3/2., dT1+lT3, 0, pT3};
Point(7) = { dT3/2., dT2 , 0, pT23};
Point(8) = { dT3/2.+lT2, dT2, 0, pT2};
Point(9) = { dT3/2.+lT2, 0, 0, pT2};
Point(10) = { dT3/2.,0,0,p0};
Line(1) = {1,2};
Line(2) = {2,3};
Line(3) = {3,4};
Line(4) = {4,1};
Line(5) = {1,10};
Line(6) = {10,7};
Line(7) = {7,4};
Line(8) = {4,5};
Line(9) = {5,6};
Line(10) = {6,7};
Line(11) = {7,8};
Line(12) = {8,9};
Line(13) = {9,10};
Line Loop(14) = {1,2,3,4};
Plane Surface(15) = {14};
Line Loop(16) = {11,12,13,6};
Plane Surface(17) = {16};
Line Loop(18) = {8,9,10,7};
Plane Surface(19) = {18};
Line Loop(20) = {5,6,7,4};
Plane Surface(21) = {20};
Translate({0,0,e}) {
Duplicata { Surface(15); }
}
Coherence;
Translate({0,0,e}) {
Duplicata { Surface(21); }
}
Coherence;
Translate({0,0,e}) {
Duplicata { Surface(17); }
}
Coherence;
Translate({0,0,e}) {
Duplicata { Surface(19); }
}
Coherence;
Line(46) = {2,12};
Line(47) = {3,14};
Line(48) = {1,11};
Line(49) = {4,16};
Line(50) = {5,26};
Line(51) = {6,28};
Line(52) = {7,20};
Line(53) = {10,18};
Line(54) = {9,24};
Line(55) = {8,22};
Line Loop(56) = {23,-47,-2,46};
Plane Surface(57) = {56};
Line Loop(58) = {22,-46,-1,48};
Plane Surface(59) = {58};
Line Loop(60) = {-49,-3,47,24};
Plane Surface(61) = {60};
Line Loop(62) = {-25,-49,4,48};
Plane Surface(63) = {62};
Line Loop(64) = {40,-50,-8,49};
Plane Surface(65) = {64};
Line Loop(66) = {41,-51,-9,50};
Plane Surface(67) = {66};
Line Loop(68) = {-42,-51,10,52};
Plane Surface(69) = {68};
Line Loop(70) = {-30,-52,7,49};
Plane Surface(71) = {70};
Line Loop(72) = {28,-53,-5,48};
Plane Surface(73) = {72};
Line Loop(74) = {29,-52,-6,53};
Plane Surface(75) = {74};
Line Loop(76) = {34,-55,-11,52};
Plane Surface(77) = {76};
Line Loop(78) = {36,-53,-13,54};
Plane Surface(79) = {78};
Line Loop(80) = {-35,-55,12,54};
Plane Surface(81) = {80};
Surface Loop(82) = {63,27,-59,-57,-61,-15};
Complex Volume(83) = {82};
Surface Loop(84) = {75,-39,77,-81,17,79};
Complex Volume(85) = {84};
Surface Loop(86) = {71,45,-65,-67,69,-19};
Complex Volume(87) = {86};
Surface Loop(88) = {71,33,-73,-75,-21,63};
Complex Volume(89) = {88};
nbpts1 = 1.0;
Transfinite Line{46,47,49,48,53,52,50,51,55,54} = nbpts1 Using Power 1.0;
nbpts2 = 10;
Transfinite Line{9,41,7,30,3,24,11,34,13,36,5,28,1,22} = nbpts2 Using Power 1.0;
nbpts3 = 15;
Transfinite Line{2,23,4,25,8,40,10,42,6,29,12,35} = nbpts3 Using Power 1.0;
Transfinite Surface{27} = {12,11,16,14};
Transfinite Surface{33} = {11,18,20,16};
Transfinite Surface{39} = {18,24,22,20};
Transfinite Surface{15} = {2,1,4,3};
Transfinite Surface{21} = {1,10,7,4};
Transfinite Surface{17} = {10,9,8,7};
Transfinite Surface{45} = {16,20,28,26};
Transfinite Surface{19} = {4,7,6,5};
Point(1) = {0,0,0,.2};
Point(2) = {0,1,0,.2};
Point(3) = {1,0,0,.2};
Point(4) = {0,0,.5,.2};
Line(7) = {3,4};
Circle(8) = {3,1,2};
Line(9) = {1,3};
Line(10) = {1,4};
Line(11) = {1,2};
Line(12) = {4,2};
Line Loop(13) = {9,8,-11};
Plane Surface(14) = {13};
Line Loop(15) = {10,12,-11};
Plane Surface(16) = {15};
Line Loop(17) = {-8,7,12};
Ruled Surface(18) = {17};
Line Loop(19) = {-10,9,7};
Plane Surface(20) = {19};
l = .3;
Point(1) = {-2,-1,0,l};
Point(2) = {-2,1,0,l};
Point(3) = {0,-1,0,l};
Point(4) = {0,1,0,l};
Point(5) = {-3,-2,0,l};
Point(6) = {-3,2,0,l};
Point(7) = {1,-2,0,l};
Point(8) = {1,2,0,l};
Line(1) = {6,8};
Line(2) = {8,7};
Line(3) = {7,5};
Line(4) = {5,6};
Line(5) = {2,1};
Line(6) = {1,3};
Line(7) = {3,4};
Line(8) = {4,2};
Loop(9) = {5,6,7,8};
Plane Surface(10) = {9};
Loop(11) = {2,3,4,1};
Plane Surface(12) = {11,9};
Transfinite Line 5 = 10 Using Power 1;
Transfinite Line 7 = 10 Using Power 1;
Transfinite Line 6 = 10 Using Power 1;
Transfinite Line 8 = 10 Using Power 1;
Recombine Surface {10} = 30;
l = 1;
Point(1) = {-2,-1,0,l};
Point(2) = {-2,1,0,l};
Point(3) = {0,-1,0,l};
Point(4) = {0,1,0,l};
Point(5) = {-3,-2,0,l};
Point(6) = {-3,2,0,l};
Point(7) = {1,-2,0,l};
Point(8) = {1,2,0,l};
Line(1) = {6,8};
Line(2) = {8,7};
Line(3) = {7,5};
Line(4) = {5,6};
Line(5) = {2,1};
Line(6) = {1,3};
Line(7) = {3,4};
Line(8) = {4,2};
/*
Loop(9) = {5,6,7,8};
Plane Surface(10) = {9};
Loop(11) = {2,3,4,1};
Plane Surface(12) = {11,9};
Transfinite Line {5,6,7,8} = 3;
Physical Point(9999) = {2};
Physical Line(2222) = {6,7};
Physical Line(3333) = {-7,-6};
Physical Association(111111)= {3333,2222};
*/
Line Loop(9) = {5,6,7,8};
Plane Surface(10) = {9};
Physical Surface(1111) = 10 ;
Point(1) = {0.0,0.0,0.0,1.0};
Point(2) = {1,0.0,0.0,1.0};
Point(3) = {.5,.1,0.0,1.0};
Point(4) = {.5,.2,0.0,1.0};
Point(5) = {0,.5,0.0,1.0};
Point(6) = {1,.55,0.0,1.0};
Point(7) = {0,0,1,1.0};
Line(1) = {2,6};
Line(2) = {5,1};
Spline(3) = {1,3,2};
Spline(4) = {5,4,6};
Line Loop(5) = {-1,-3,-2,4};
Plane Surface(6) = {5};
Point(8) = {0,.5,1,1.0};
Point(9) = {1,.5,1,1.0};
Point(10) = {1,0,1,1.0};
Line(7) = {6,9};
Line(8) = {9,10};
Line(9) = {10,2};
Line(10) = {10,7};
Line(11) = {7,1};
Line(12) = {5,8};
Line(13) = {8,7};
Line(14) = {8,9};
Line Loop(15) = {14,-7,-4,12};
Ruled Surface(16) = {15};
Line Loop(17) = {7,8,9,1};
Ruled Surface(18) = {17};
Line Loop(19) = {-10,-8,-14,13};
Ruled Surface(20) = {19};
Line Loop(21) = {-11,-13,-12,2};
Ruled Surface(22) = {21};
Line Loop(23) = {-9,10,11,3};
Ruled Surface(24) = {23};
Surface Loop(25) = {16,20,24,18,6,22};
Complex Volume(26) = {25};
Transfinite Line {7,9,11,12} = 5;
Transfinite Line {3,10,14} = 10;
Transfinite Line {4} = 10 Using Bump 1.5;
Transfinite Line {-1,2} = 20 Using Power 1.6;
Transfinite Line {8,13} = 20 ;
Transfinite Surface {6} = {1,2,6,5};
Transfinite Surface {20} = {7,10,9,8};
Transfinite Surface {18} = {10,2,6,9};
Transfinite Surface {22} = {7,1,5,8};
Transfinite Surface {16} = {8,9,6,5};
Transfinite Surface {24} = {7,10,2,1};
Transfinite Volume {26} = {7,10,2,1,8,9,6,5};
Recombine Surface {6,20,18,22,16,24};
Point(1) = {0.0,0.0,0.0,1.0};
Point(2) = {1,0.0,0.0,1.0};
Point(3) = {.5,.1,0.0,1.0};
Point(4) = {.5,.2,0.0,1.0};
Point(5) = {0,.5,0.0,1.0};
Point(6) = {1,.55,0.0,1.0};
Point(7) = {0,0,1,1.0};
Line(1) = {2,6};
Line(2) = {5,1};
Spline(3) = {1,3,2};
Spline(4) = {5,4,6};
Point(8) = {0,.5,1,1.0};
Point(9) = {1,.5,1,1.0};
Point(10) = {1,0,1,1.0};
Line(7) = {6,9};
Line(8) = {9,10};
Line(9) = {10,2};
Line(10) = {10,7};
Line(11) = {7,1};
Line(12) = {5,8};
Line(13) = {8,7};
Line(14) = {8,9};
This diff is collapsed.
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment