Dear Sir.
I have a problem with a face visualization in Gmsh.
After extruding of non plane surface (cone surface) it is imposible to
show surfaces and surface's numbers in Gmsh.
But except that, model seems to be correct also meshing is correct.
Do you know this problem?
Is it impossible to extrude non plane (cone) surfaces?
Unfortunately I'm trying to prepare a complicated geometry, and I
divide the problem into several files.
The main file to open in Gmsh is: "t1.geo"
Probem with extrusion is in the file: "t1_2_palce_wsad.geo" (sorry for
such name) in line 30.
Other files contain:
t1_1_dno2.geo : begining of geometry,
t1.par : geometrical parameters,
t1.m : mesh parameters,
t1.v_s : physical volumes and surfaces definition
Gmsh version: 3.0.2, 2017 05 13
FLTK version: 1.3.4; PETSc version 3.7.5; OCC version: 7.1.0; MED
version: 3.2.0
Operating system: Windows 8.1 Pro 64 Polish language version
Sincerely,
and with request for help
Roman Przylucki
Include "t1.par";
Include "t1_2_palce_wsad.geo";
//nowy model starego tygla 15.05.2017
//2D mesh Mesh.Algorithm
//2D mesh algorithm (1=MeshAdapt, 2=Automatic, 5=Delaunay, 6=Frontal, 7=BAMG,
8=DelQuad)
//Mesh.Algorithm = 2;
//3D mesh algorithm (1=Delaunay, 2=New Delaunay, 4=Frontal, 5=Frontal Delaunay,
6=Frontal Hex, 7=MMG3D, 9=R-tree)
Mesh.Algorithm3D = 2; //nowy delaunay dziala
//Mesh.Algorithm3D = 4; //frontal
//+
Coherence;
//+
Coherence;
Include "t1.par";
//definicja siatek
//podzial dla punktow
ms = 4/s;
//mesh dla dna
//podzia³y po k¹cie
//szczelina powietrzna
md1a=4;
//scianka palca i pochodne
md1b=4;
//pozostale wnetrze dna i palca
md1c=15;
//podzialy wzdluzne
//od srodka do ...
md1d=15;
//zapas
md1e=8;
//naciecie
md1f=8;
//scianka palca
md1g=5;
//x otworu palca
md1h=10;
// szczelina wzbudnik wsad
md1i=4;
//scianka wzbudnika
md1j=4;
// otwor wzbudnika
md1k=10;
//zewnetrzna czesc dna
md1l=12;
//podzialy pionowe
//srodek
md12a=7;
//reszta
md12b=11;
//koniec siatek dla dna
// siatki dla wsadu i palca
//podzialy pionowe
//na wysokosci palca bez gornego dekla
mwpa=41;
//podzial wysokosci gornego dekla
mwpb=4;
//parametry geometryczne
//deklaracje parametrow geometrycznych
pi = 3.141592;
// skala
s = 1e-3;
//kat wycinka segmentu
alfa=2*pi/16;
//kat
//parametry dla dna
//po x czyli promieniu
//wymiar x dna tygla (ca³ego)
x_d = 44*s;
//g³êbokosc naciêcia dna
xnac_d=4*s;
//x zapasu dna
xzapas_d=10*s;
//x punktu dodatkowego (nieu¿ywane)
xpd_d=x_d/4;
//po z
//pol szczeliny miedzy palcami
//to samo w dnie
spp=0.5*s;
//koniec parametrow dna
//wymiary dla palca
//grubosc scianki palca rownomierna dla wszystkich wymiarow
gs_p=2*s;
//zewnetrzny wymiar x palca
xz_p=15*s;
//wymiar x wnetrza kanalu palca
xw_p=xz_p-2*gs_p;
//fi i x szczeliny powietrznej na xd
//d(alfa)=Arcsin(dz/d)
fi1_x_d=Asin(spp/(x_d));
//dx_fi1_d=2*x_d*Sin(fi1_x_d/2);
//fi grubosci scianki palca na xd+gs_p
fi2_x_d=Asin(gs_p/(x_d+gs_p));
//wymiar x calego dna = o 20 mm szersze od palcow
xdt_d=x_d+xz_p+20*s;
//wymiary dla wzbudnika i zwojow
//szczelina wzbudnik - palec
xsp_i=2*s;
//wewnetrzny promien wzbudnika
xw_i=x_d+xz_p+xsp_i;
//x wymiary duzego i ma³ego zwoju
xd_z=10*s;
xm_z=10*s;
//zewnetrzny promien wzbudnika du¿ego i ma³ego
xzd_i=xw_i+xd_z;
xzm_i=xw_i+xm_z;
//grubosc scianki zwoju wzbudnika
gs_i=2*s;
//definiowanie bocznika magnetycznego
//szczelina miêdzy bocznikiem a wzbudnikiem
xsi_b=3*s;
//poczatek bocznika
xp_b=xzd_i+xsi_b;
//wymiar x bocznika
x_b=15*s;
//wymiary otoczenia po x
x_ot= 2*xp_b;
//po y czyli wysokosci
//wysokosc zewnetrznej czesci dna
yz_d=10*s;
//glebokosc dna
yg_d=6*s;
//wysokosc srodka dna
ys_d=yz_d-yg_d;
//wysokosc tygla na promieniu / odleg³osci x
//wys = x*(yz_d-ys_d)/xd +ys_d;
//wysokosc na wsp naciecia
ynac_d=(x_d-xnac_d)*(yz_d-ys_d)/x_d+ys_d;
//wysokosc na wsp zaokraglenia dna
yzapas_d=(x_d-xzapas_d)*(yz_d-ys_d)/x_d+ys_d;
//wysokosc calego palca (razem z dnem)
h_p=140*s+ys_d;
//grubosc gornego dekla palca
gd_p=2*s;
//cala wysokosc palca nad dno (jego zewnêtrzna czesc)
hnd_p=h_p-yz_d;
//wysokosc palca nad dno (jego zewnêtrzna czesc) do poczatku dekla
hndpd_p=hnd_p-gd_p;
//wysokosc palca nad srodek dna
hnds_p=hnd_p+yg_d;
/*
//wysokosc wsadu nad zewnetrzna czesc gory palca
naddatek_wsadu=30*s;
//bezwzgledna wysokosc wsadu nad zewnetrzna czesc palca
bw_w=h_p+naddatek_wsadu;
//wysokosc translacji na promieniu naciecia
ty_xn=(yz_d-yzaok_d)+hnd_p;
//wysokosci zwojow i wzbudnika
//odleg³osc 1 zwoju od podstawy
ypi_i=10*s;
//polozenie y 1 zwoju
yd_i=yz_d+ypi_i;
//wymiar y malega i duzego zwoju
ym_z=14*s;
yd_z=20*s;
//szczelina miêdzy zwojami po y
yzz_i=4*s;
//translaca 1 zwoju (ma³ego)
t1z_y=ym_z+yzz_i;
//y poczatkowe 1 duzego zwoju
yddz_i=yd_i+2*ym_z+2*yzz_i;
//translacja 3 zwoju (du¿ego)
t3z_y=yd_z+yzz_i;
//calkowita bezwzgledna wysokosc wzbudnika
ytot_i=yddz_i+2*yd_z+yzz_i;
//translacja gory wzbudnika do wysokosci konca palcow
tg_i= h_p-ytot_i;
//wymiarowanie bocznika magnetycznego po y
// odleglosc bocznika od podstawy
yp_b=10*s;
//polozenie y spodu bocznika
yd_b=yz_d+yp_b;
//wymiar y bocznika domyslne = wysokosci wzbudnika
y_b=2*ym_z+2*yd_z+3*yzz_i;
//translacja gory bocznika do wysokosci konca palcow
tg_b=h_p-yd_b-y_b;
//wymiary otoczenia zewnetrznego
//ponizej dna
y_pod_dno_ot=h_p/2;
//powyzej naddatku wsadu
y_nad_wsad_ot= 1.5*(naddatek_wsadu+h_p);
//po z czyli przod tyl, kat
//kat dla szczeliny miedzy palcami na promieniu (x) naciecia x=xd-xnac_d
fi1=Asin(s12pp_z/(xd-xnac_d));
fi11=alfa-2*fi1;
//kat dla szczeliny miedzy palcami na promieniu (x) zaokraglenia x=xd-xzaok_d_d
fi2=Asin(s12pp_z/(xd-xzaok_d));
fi21=alfa-2*fi2;
//kat dla szczeliny miedzy palcami na promieniu (x) wnetrza dna x=xd
fi3=Asin(s12pp_z/(xd));
fi31=alfa-2*fi3;
//kat dla szczeliny miedzy palcami na promieniu (x) wnetrza
//blizsze wnetrze palca x=xd+gsp
fi4=Asin(s12pp_z/(xd+gsp));
fi41=alfa-2*fi4;
//kat dla szczeliny miedzy palcami na promieniu (x) wnetrza
//dalsze wnetrze palca x=xd+xz_p-gsp
fi5=Asin(s12pp_z/(xd+xz_p-gsp));
fi51=alfa-2*fi5;
//kat dla szczeliny miedzy palcami na promieniu (x) zewnetrze palca x=xd+xz_p
fi6=Asin(s12pp_z/(xd+xz_p));
fi61=alfa-2*fi6;
//kat dla szczeliny miedzy palcami na promieniu (x) zewnetrze dna x = xdt
fi7=Asin(s12pp_z/(xdt));
fi71=alfa-2*fi7;
//kat dla grubosci scianki palca (blizej tygla)
fi8=Asin(gsp/(xd+gsp));
fi81=alfa-(2*fi8+2*fi4);
//kat dla grubosci scianki palca (dalej tygla)
fi9=Asin(gsp/(xd+xz_p-gsp));
fi91=alfa-(2*fi9+2*fi5);
//wymiar z bocznika
z_b=15;
//stopien wypelnienia bocznikiem 0 do 1
stopien_wyp=0.5;
//kat wypelnienia bocznika
fi10= (alfa-alfa*(1-stopien_wyp));
//kat poczatku wzbudnika
fi101=(alfa-fi10)/2;
//gestosci siatek
//siatka po x (raczej)
s1=xd/10;
//s1=3*s;
*/
//numery objetosci i powierzchni
//t1_dno objetosci i obszary
//dno tygla bez palcow ale z otoczeniem dna (t1_dno)
//numery objetosci fizycznych
//miedz
dno_tygla=5001;
zewn_dna=5002; //poza palcami
//powietrze
szcz_dna=5003; //naciecie tygla wewnatrz tygla przod i tyl
szcz_dna_z=5004; //naciecie poza palcami tygla przod i tyl
//numery powierzchni bocznych dla dna bez palcow
s_dno_p = 6001;
s_dno_t = 6002;
s_dno_z_p = 6003;
s_dno_z_t = 6004;
//palce tygla w dnie
//objetosci
palce_dno=5005;
szcz_palcow_dno=5006;
palce_pow_dno=5007;
//powierzchnie palcow w dnie
s_palce_dno_p=6005;
s_palce_dno_t=6006;
//koniec t1_dno objetosci i obszary
//poczatek t1_palce objetosci i obszary
//objetosci
palce=5008;
szcz_palcow=5009;
palce_pow=5010;
//powierzchnie palcow
s_palce_p=6007;
s_palce_t=6008;
//koniec t1 palce objetosci i obszary
//wsad objetosci i obszary
//objetosci
//koniec t1 wsad objetosci i obszary
/*
wsad_w_tyglu=5001;
wsad_nad_tyglem=5002;
palec_tygla=5003;
bocznik=5004;
//zwoje od do³u
z1=5005;
z2=5006;
z3=5007;
z4=5008;
powietrze_w_palcu=5050;
powietrze_obok_bocznika=5051;
powietrze_wzbudnik=5052;
powietrze_wzb_palec=5053;
//reszta powietrza bez szczelin
powietrze_reszta=5054;
powietrze_szcz_przod=5055;
powietrze_szcz_tyl=5056;
//uwaga tu tez jest doliczony wsad nad tyglem
powietrze_cale=5060;
dno_tygla=5009;
dno_za_palcami=5010;
//=======================================
//wersja 12.02.2016
//po³¹czone objêtosci
//wsad jest tylko wy tyglu wsad poza tyglem to powietrze
//palec tygla i dno tygla sa razem jako tygiel
//po³aczone sa wszystkie powietrza
//=======================================
wsad=5101;
tygiel=5102;
bocznik_a=5103;
powietrze=5104;
z1_a=5105;
z2_a=5106;
z3_a=5107;
z4_a=5108;
//DEFINIOWANIE POWIERZCHNI
//powierzchnia przodu wsp z = 0
//to sp_
//powierzchnie tylu
//to st_
sp_wsad_tygiel=6001;
st_wsad_tygiel=6101;
sp_wsad_nad=6021;
st_wsad_nad=6121;
sp_z1=6002;
st_z1=6102;
sp_z2=6003;
st_z2=6103;
sp_z3=6004;
st_z3=6104;
sp_z4=6005;
st_z4=6105;
sp_tygiel=6006;
st_tygiel=6106;
//powietrze numeruje od dolu od osi modelu
//powietrze pod wsadem
sp_pow_1=6007;
st_pow_1=6107;
//powietrze szczelin
sp_pow_2=6008;
st_pow_2=6108;
//reszta pow do wysokosci tygla
sp_pow_3=6009;
st_pow_3=6109;
//powietrze nad tyglem
sp_pow_4=6010;
st_pow_4=6110;
//brzeg na zewnatrz
s_bok=6050;
//=======================================
//wersja 12.02.2016
//po³¹czone powierzchnie
//wsad jest tylko wy tyglu wsad poza tyglem to powietrze
//po³aczone sa wszystkie powietrza
//=======================================
//wsad
sp_wsad=6201;
st_wsad=6301;
//zwoje
sp_z1_a=6202;
st_z1_a=6302;
sp_z2_a=6203;
st_z2_a=6303;
sp_z3_a=6204;
st_z3_a=6304;
sp_z4_a=6205;
st_z4_a=6305;
//tygiel
sp_tygiel_a=6206;
st_tygiel_a=6306;
//powietrze
sp_pow= 6207;
st_pow= 6307;
//dirichlet
s_bok_a=6250;
*/
//plik 1
//definiowanie dna
Include "t1.par";
Include "t1.m";
Include "t1.v_s";
SetFactory("OpenCASCADE");
// punkty zewnetrznej linii podstawy dna
//v1d11
v1d11_p1=newp; Point(v1d11_p1)={0,0,0,ms};
v1d11_p2=newp; Point(v1d11_p2)={x_d-xzapas_d,0,0,ms};
v1d11_p3=newp; Point(v1d11_p3)={x_d-xnac_d,0,0,ms};
v1d11_p4=newp; Point(v1d11_p4)={x_d,0,0,ms};
v1d11_p5=newp; Point(v1d11_p5)={x_d+gs_p,0,0,ms};
v1d11_p6=newp; Point(v1d11_p6)={x_d+xz_p-gs_p,0,0,ms};
v1d11_p7=newp; Point(v1d11_p7)={x_d+xz_p,0,0,ms};
//to powyzej by³o z palcami
// to poni¿ej poza palcami czyli wzbudnik
v1d11_p8=newp; Point(v1d11_p8)={xw_i,0,0,ms};
v1d11_p9=newp; Point(v1d11_p9)={xw_i+gs_i,0,0,ms};
v1d11_p10=newp; Point(v1d11_p10)={xzm_i-gs_i,0,0,ms};
v1d11_p11=newp; Point(v1d11_p11)={xzm_i,0,0,ms};
v1d11_p12=newp; Point(v1d11_p12)={xdt_d,0,0,ms};
//linie w podstawie
v1d11_l1=newreg; Line(v1d11_l1) = {v1d11_p1,v1d11_p2};
v1d11_l2=newreg; Line(v1d11_l2) = {v1d11_p2,v1d11_p3};
v1d11_l3=newreg; Line(v1d11_l3) = {v1d11_p3,v1d11_p4};
v1d11_l4=newreg; Line(v1d11_l4) = {v1d11_p4,v1d11_p5};
v1d11_l5=newreg; Line(v1d11_l5) = {v1d11_p5,v1d11_p6};
v1d11_l6=newreg; Line(v1d11_l6) = {v1d11_p6,v1d11_p7};
//to powyzej to linie z palcami
// to poni¿ej to linie poza palcami czyli wzbudnik i reszta
v1d11_l7=newreg; Line(v1d11_l7) = {v1d11_p7,v1d11_p8};
v1d11_l8=newreg; Line(v1d11_l8) = {v1d11_p8,v1d11_p9};
v1d11_l9=newreg; Line(v1d11_l9) = {v1d11_p9,v1d11_p10};
v1d11_l10=newreg; Line(v1d11_l10) = {v1d11_p10,v1d11_p11};
v1d11_l11=newreg; Line(v1d11_l11) = {v1d11_p11,v1d11_p12};
//wyciaganie w pionie
v1d12_e1[]= Extrude {0,ys_d,0} {Point{v1d11_p1};};
v1d12_e2[]= Extrude {0,yzapas_d,0} {Point{v1d11_p2};};
v1d12_e3[]= Extrude {0,ynac_d,0} {Point{v1d11_p3};};
//x[]=Boundary{Line{v1d12_e1[0] };};
//Printf("punkty !!! %g , %g, %g, %g,%g , %g, %g, %g" ,
x[0],x[1],x[2],x[3],x[4],x[5],x[6],x[7]);
//wyciaganie w górê reszty dna o sta³ej wysokosci y
//mesh
//layers
v1d12_e4[]= Extrude {0,yz_d,0}
{Line{v1d11_l4, v1d11_l5, v1d11_l6, v1d11_l7, v1d11_l8,v1d11_l9,
v1d11_l10, v1d11_l11 };
Layers {md12b-1};Recombine;};
//linie po wyciagnieciu czyli przod (pierwsza warstwa)
// powyzej czyli (dno rzeczywstego tygla)
//czesc o nierownej wysokosci
p[]=Boundary{Line{v1d12_e1[0], v1d12_e2[0], v1d12_e3[0],v1d12_e4[2] };};
//Printf("punkty !!! %g , %g, %g, %g,%g , %g, %g, %g" ,
p[0],p[1],p[2],p[3],p[4],p[5],p[6],p[7]);
v1d21_l1=newreg; Line(v1d21_l1) = {p[1],p[3]};
v1d21_l2=newreg; Line(v1d21_l2) = {p[3],p[5]};
v1d21_l3=newreg; Line(v1d21_l3) = {p[5],p[7]};
//boczne pow. o nierównym wymiarze y
v1d1_L1=newreg; Line Loop (v1d1_L1) = { v1d11_l1, v1d12_e2[0], -v1d21_l1, -
v1d12_e1[0] };
v1d1_s1 = newreg; Plane Surface (v1d1_s1) = {v1d1_L1};
v1d1_L2=newreg; Line Loop (v1d1_L2) = { v1d11_l2, v1d12_e3[0], -v1d21_l2, -
v1d12_e2[0] };
v1d1_s2 = newreg; Plane Surface (v1d1_s2) = {v1d1_L2};
v1d1_L3=newreg; Line Loop (v1d1_L3) = { v1d11_l3, v1d12_e4[2] , -v1d21_l3, -
v1d12_e3[0] };
v1d1_s3 = newreg; Plane Surface (v1d1_s3) = {v1d1_L3};
v1d_e1[] = Extrude{ {0,1,0}, {0,0,0}, fi1_x_d}
{Surface{ v1d1_s1,
v1d1_s2, v1d1_s3, v1d12_e4[1], v1d12_e4[5],v1d12_e4[9],
v1d12_e4[13], v1d12_e4[17], v1d12_e4[21],v1d12_e4[25],
v1d12_e4[29] };
Layers{md1a}; Recombine; };
//grubosc scianki wzbudnika
//Printf("elementy %g , %g, %g, %g,%g , %g, %g, %g, %g" ,
v1d_e1[0],v1d_e1[1],v1d_e1[2],v1d_e1[3],v1d_e1[4],v1d_e1[5],v1d_e1[6],v1d_e1[7],v1d_e1[8]);
v1d_e2[] = Extrude{ {0,1,0}, {0,0,0}, fi2_x_d}
{Surface{ v1d_e1[0],
v1d_e1[5],v1d_e1[11],v1d_e1[17],v1d_e1[23],
v1d_e1[29],v1d_e1[35],v1d_e1[41],v1d_e1[47],v1d_e1[53],v1d_e1[59] };
Layers{md1b}; Recombine; };
v1d_e3[] = Extrude{ {0,1,0}, {0,0,0}, alfa-(fi2_x_d+fi1_x_d)}
{Surface{ v1d_e2[0],
v1d_e2[5],v1d_e2[11],v1d_e2[17],v1d_e2[23],
v1d_e2[29],v1d_e2[35],v1d_e2[41],v1d_e2[47],v1d_e2[53],v1d_e2[59] };
Layers{md1c}; Recombine; };
v1d_e4[] = Extrude{ {0,1,0}, {0,0,0}, fi2_x_d}
{Surface{ v1d_e3[0],
v1d_e3[5],v1d_e3[11],v1d_e3[17],v1d_e3[23],
v1d_e3[29],v1d_e3[35],v1d_e3[41],v1d_e3[47],v1d_e3[53],v1d_e3[59] };
Layers{md1b}; Recombine; };
v1d_e5[] = Extrude{ {0,1,0}, {0,0,0}, fi1_x_d}
{Surface{ v1d_e4[0],
v1d_e4[5],v1d_e4[11],v1d_e4[17],v1d_e4[23],
v1d_e4[29],v1d_e4[35],v1d_e4[41],v1d_e4[47],v1d_e4[53],v1d_e4[59] };
Layers{md1a}; Recombine; };
//zrobione ca³e dno
//======================
/*
//++++++++++++++++++++++
//definicja objêtosci i powierzchni
//objetosci
//miedz
Physical Volume(dno_tygla)={ v1d_e1[1], v1d_e2[1],v1d_e3[1], v1d_e4[1],
v1d_e5[1],
v1d_e1[6], v1d_e2[6],v1d_e3[6], v1d_e4[6],
v1d_e5[6],
v1d_e2[12],v1d_e3[12], v1d_e4[12],
v1d_e2[36],v1d_e3[36], v1d_e4[36],
v1d_e2[42],v1d_e3[42], v1d_e4[42],
v1d_e2[48],v1d_e3[48], v1d_e4[48],
v1d_e2[54],v1d_e3[54], v1d_e4[54],
v1d_e2[60],v1d_e3[60], v1d_e4[60] };
//powietrze
Physical Volume(szcz_dna)={ v1d_e1[12], v1d_e5[12]};
Physical Volume(szcz_dna_z)={v1d_e1[36], v1d_e5[36],
v1d_e1[42], v1d_e5[42],
v1d_e1[48], v1d_e5[48],
v1d_e1[54], v1d_e5[54],
v1d_e1[60], v1d_e5[60]};
//powierzchnie przednie z=0 i tylne dna
Physical Surface(s_dno_p)={v1d1_s1,v1d1_s2,v1d1_s3,
v1d12_e4[13],v1d12_e4[17],v1d12_e4[21],v1d12_e4[25],v1d12_e4[29] };
Physical Surface(s_dno_t)={v1d_e5[0],v1d_e5[5],v1d_e5[11],
v1d_e5[35],v1d_e5[41],v1d_e5[47],v1d_e5[53],v1d_e5[59]};
//objetosci palca w dnie
//miedz
Physical Volume(palce_dno)={ v1d_e2[18],v1d_e3[18], v1d_e4[18],
v1d_e2[24], v1d_e4[24],
v1d_e2[30],v1d_e3[30], v1d_e4[30]};
//powietrze
Physical Volume(szcz_palcow_dno)={ v1d_e1[18], v1d_e5[18],
v1d_e1[24], v1d_e5[24],
v1d_e1[30], v1d_e5[30]};
Physical Volume(palce_pow_dno)={v1d_e3[24]};
//powierzchnie przednie z=0 i tylne palcow w dnie
Physical Surface(s_palce_dno_p)={ v1d12_e4[1],v1d12_e4[5],v1d12_e4[9]};
Physical Surface(s_palce_dno_t)={ v1d_e5[17],v1d_e5[23],v1d_e5[29] };
//++++++++++++++++++++++
*/
//mesh
//po promieniu
//wnetrze
Transfinite Line{ v1d11_l1} = md1d;
Transfinite Line{ v1d21_l1} = md1d;
l[]=Boundary{Surface{v1d_e1[0],v1d_e2[0],v1d_e3[0],v1d_e4[0],v1d_e5[0]};};
For m In {0:18:4}
Transfinite Line{ l[m]} = md1d;
EndFor
For m In {2:18:4}
Transfinite Line{ -l[m]} = md1d;
EndFor
//zapas
Transfinite Line{ v1d11_l2} = md1e;
Transfinite Line{ v1d21_l2} = md1e;
l[]=Boundary{Surface{v1d_e1[5],v1d_e2[5],v1d_e3[5],v1d_e4[5],v1d_e5[5]};};
For m In {0:18:4}
Transfinite Line{ l[m]} = md1e;
EndFor
For m In {2:18:4}
Transfinite Line{ -l[m]} = md1e;
EndFor
//naciecie
Transfinite Line{ v1d11_l3} = md1f;
Transfinite Line{ v1d21_l3} = md1f;
l[]=Boundary{Surface{v1d_e1[11],v1d_e2[11],v1d_e3[11],v1d_e4[11],v1d_e5[11]};};
For m In {0:18:4}
Transfinite Line{ l[m]} = md1f;
EndFor
For m In {2:18:4}
Transfinite Line{ -l[m]} = md1f;
EndFor
//przod palca
Transfinite Line{ v1d11_l4} = md1g;
Transfinite Line{ v1d12_e4[0]} = md1g;
l[]=Boundary{Surface{v1d_e1[17],v1d_e2[17],v1d_e3[17],v1d_e4[17],v1d_e5[17]};};
For m In {1:19:4}
Transfinite Line{ l[m]} = md1g;
EndFor
For m In {3:19:4}
Transfinite Line{ -l[m]} = md1g;
EndFor
//srodek palca
Transfinite Line{ v1d11_l5} = md1h;
Transfinite Line{ v1d12_e4[4]} = md1h;
l[]=Boundary{Surface{v1d_e1[23],v1d_e2[23],v1d_e3[23],v1d_e4[23],v1d_e5[23]};};
For m In {1:19:4}
Transfinite Line{ l[m]} = md1h;
EndFor
For m In {3:19:4}
Transfinite Line{ -l[m]} = md1h;
EndFor
//tyl palca
Transfinite Line{ v1d11_l6} = md1g;
Transfinite Line{ v1d12_e4[8]} = md1g;
l[]=Boundary{Surface{v1d_e1[29],v1d_e2[29],v1d_e3[29],v1d_e4[29],v1d_e5[29]};};
For m In {1:19:4}
Transfinite Line{ l[m]} = md1g;
EndFor
For m In {3:19:4}
Transfinite Line{ -l[m]} = md1g;
EndFor
//szczelina palec wzbudnik
Transfinite Line{ v1d11_l7} = md1i;
Transfinite Line{ v1d12_e4[12]} = md1i;
l[]=Boundary{Surface{v1d_e1[35],v1d_e2[35],v1d_e3[35],v1d_e4[35],v1d_e5[35]};};
For m In {1:19:4}
Transfinite Line{ l[m]} = md1i;
EndFor
For m In {3:19:4}
Transfinite Line{ -l[m]} = md1i;
EndFor
//scianka wzbudnika wewnetrzna
Transfinite Line{ v1d11_l8} = md1j;
Transfinite Line{ v1d12_e4[16]} = md1j;
l[]=Boundary{Surface{v1d_e1[41],v1d_e2[41],v1d_e3[41],v1d_e4[41],v1d_e5[41]};};
For m In {1:19:4}
Transfinite Line{ l[m]} = md1j;
EndFor
For m In {3:19:4}
Transfinite Line{ -l[m]} = md1j;
EndFor
//otwor wzbudnika
Transfinite Line{ v1d11_l9} = md1k;
Transfinite Line{ v1d12_e4[20]} = md1k;
l[]=Boundary{Surface{v1d_e1[47],v1d_e2[47],v1d_e3[47],v1d_e4[47],v1d_e5[47]};};
For m In {1:19:4}
Transfinite Line{ l[m]} = md1k;
EndFor
For m In {3:19:4}
Transfinite Line{ -l[m]} = md1k;
EndFor
//scianka wzbudnika zewnetrzna
Transfinite Line{ v1d11_l10} = md1j;
Transfinite Line{ v1d12_e4[24]} = md1j;
l[]=Boundary{Surface{v1d_e1[53],v1d_e2[53],v1d_e3[53],v1d_e4[53],v1d_e5[53]};};
For m In {1:19:4}
Transfinite Line{ l[m]} = md1j;
EndFor
For m In {3:19:4}
Transfinite Line{ -l[m]} = md1j;
EndFor
//do konca dna
Transfinite Line{ v1d11_l11} = md1l;
Transfinite Line{ v1d12_e4[28]} = md1l;
l[]=Boundary{Surface{v1d_e1[59],v1d_e2[59],v1d_e3[59],v1d_e4[59],v1d_e5[59]};};
For m In {1:19:4}
// Printf("line %g ",l[m]);
Transfinite Line{ l[m]} = md1l;
EndFor
For m In {3:19:4}
// Printf("line %g ",l[m]);
Transfinite Line{ -l[m]} = md1l;
EndFor
//koniec podzialow po x
//po wysokosci
Transfinite Line{v1d12_e1[0]}=md12a;
Transfinite Line{v1d12_e2[0],v1d12_e3[0] } = md12b;
//(reszta jest zrobiona na pocz¹tku w extrude)
//prostokat od srodka dna
Recombine Surface{v1d1_s1};
// 2 nastepny
Transfinite Surface{v1d1_s2,v1d1_s3};
//pozostale s¹ wyci¹gniête po wysokosci
//uwaga wyci¹gniecie = zwykly podzial -1 !!!
//plik 2
//definiowanie palcow i wsadu
//wyciaganie ich zdna
Include "t1.par";
Include "t1.m";
Include "t1.v_s";
Include "t1_1_dno2.geo";
SetFactory("OpenCASCADE");
/*
For m In {4:15:5}
Printf("indeks, e1 %g, %g ", m, v1d_e1[m]);
//Printf("indeks, e2 %g, %g ", m, vp_e1[m]);
EndFor
*/
// Printf (" surface %g, %g, %g ",v1d_e1[22],v1d_e1[28],v1d_e1[34] );
//wyciaganie wsadu i dna na wysokosc palcow
//mesh
//layers
//!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
//problem
//!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
vp_e1[]= Extrude {0,hndpd_p,0} {Surface{
/*
below in first section is set of cone surfaces to extrude
which cause the problem
if you comment it is OK
*/
// /*
v1d_e1[4],v1d_e1[9],v1d_e1[15],
v1d_e2[4],v1d_e2[9],v1d_e2[15],
v1d_e3[4],v1d_e3[9],v1d_e3[15],
v1d_e4[4],v1d_e4[9],v1d_e4[15],
v1d_e5[4],v1d_e5[9],v1d_e5[15],
// */
v1d_e1[22],v1d_e1[28],v1d_e1[34],
v1d_e2[22],v1d_e2[28],v1d_e2[34],
v1d_e3[22],v1d_e3[28],v1d_e3[34],
v1d_e4[22],v1d_e4[28],v1d_e4[34],
v1d_e5[22],v1d_e5[28],v1d_e5[34]
};
Layers {mwpa-1};};
/*
For m In {0:40:1}
Printf("indeks, pow gorne %g, %g ", m, vp_e1[m]);
EndFor
//mesh
//layers
vp_e2[]= Extrude {0,gd_p,0} {Surface{ vp_e1[0], vp_e1[6], vp_e1[12],
vp_e1[18], vp_e1[24], vp_e1[30],
vp_e1[36], vp_e1[42], vp_e1[48],
vp_e1[54], vp_e1[60], vp_e1[66],
vp_e1[72], vp_e1[78], vp_e1[84] };};
//zrobione ca³e palce
//======================
For m In {1:84:6}
Printf("indeks, pow gorne %g, %g ", m, vp_e1[m]);
EndFor
//++++++++++++++++++++++
//definicja objêtosci i powierzchni dla palcow
//objetosci
//miedz
Physical Volume(palce)={ vp_e1[1], vp_e1[7], vp_e1[13],
vp_e1[19], vp_e1[25], vp_e1[31],
vp_e1[37], vp_e1[43], vp_e1[49],
vp_e1[55], vp_e1[61], vp_e1[67],
vp_e1[73], vp_e1[79], vp_e1[85] };
//powietrze
Physical Volume(szcz_dna)={ v1d_e1[12], v1d_e5[12]};
Physical Volume(szcz_dna_z)={v1d_e1[36], v1d_e5[36],
v1d_e1[42], v1d_e5[42],
v1d_e1[48], v1d_e5[48],
v1d_e1[54], v1d_e5[54],
v1d_e1[60], v1d_e5[60]};
//powierzchnie przednie z=0 i tylne dna
Physical Surface(s_dno_p)={v1d1_s1,v1d1_s2,v1d1_s3,
v1d12_e4[13],v1d12_e4[17],v1d12_e4[21],v1d12_e4[25],v1d12_e4[29] };
Physical Surface(s_dno_t)={v1d_e5[0],v1d_e5[5],v1d_e5[11],
v1d_e5[35],v1d_e5[41],v1d_e5[47],v1d_e5[53],v1d_e5[59]};
//objetosci palca w dnie
//miedz
Physical Volume(palce_dno)={ v1d_e2[18],v1d_e3[18], v1d_e4[18],
v1d_e2[24],v1d_e3[24], v1d_e4[24],
v1d_e2[30],v1d_e3[30], v1d_e4[30]};
//powietrze
Physical Volume(szcz_palcow_dno)={ v1d_e1[18], v1d_e5[18],
v1d_e1[24], v1d_e5[24],
v1d_e1[30], v1d_e5[30]};
//powierzchnie przednie z=0 i tylne palcow w dnie
Physical Surface(s_palce_dno_p)={ v1d12_e4[1],v1d12_e4[5],v1d12_e4[9]};
Physical Surface(s_palce_dno_t)={ v1d_e5[17],v1d_e5[23],v1d_e5[29] };
//++++++++++++++++++++++
//mesh
//po promieniu
//wnetrze
// punkty zewnetrznej linii podstawy dna
//v1d11
v1d11_p1=newp; Point(v1d11_p1)={0,0,0,ms};
v1d11_p2=newp; Point(v1d11_p2)={x_d-xzapas_d,0,0,ms};
v1d11_p3=newp; Point(v1d11_p3)={x_d-xnac_d,0,0,ms};
v1d11_p4=newp; Point(v1d11_p4)={x_d,0,0,ms};
v1d11_p5=newp; Point(v1d11_p5)={x_d+gs_p,0,0,ms};
v1d11_p6=newp; Point(v1d11_p6)={x_d+xz_p-gs_p,0,0,ms};
v1d11_p7=newp; Point(v1d11_p7)={x_d+xz_p,0,0,ms};
//to powyzej by³o z palcami
// to poni¿ej poza palcami czyli wzbudnik
v1d11_p8=newp; Point(v1d11_p8)={xw_i,0,0,ms};
v1d11_p9=newp; Point(v1d11_p9)={xw_i+gs_i,0,0,ms};
v1d11_p10=newp; Point(v1d11_p10)={xzm_i-gs_i,0,0,ms};
v1d11_p11=newp; Point(v1d11_p11)={xzm_i,0,0,ms};
v1d11_p12=newp; Point(v1d11_p12)={xdt_d,0,0,ms};
//linie w podstawie
v1d11_l1=newreg; Line(v1d11_l1) = {v1d11_p1,v1d11_p2};
v1d11_l2=newreg; Line(v1d11_l2) = {v1d11_p2,v1d11_p3};
v1d11_l3=newreg; Line(v1d11_l3) = {v1d11_p3,v1d11_p4};
v1d11_l4=newreg; Line(v1d11_l4) = {v1d11_p4,v1d11_p5};
v1d11_l5=newreg; Line(v1d11_l5) = {v1d11_p5,v1d11_p6};
v1d11_l6=newreg; Line(v1d11_l6) = {v1d11_p6,v1d11_p7};
//to powyzej to linie z palcami
// to poni¿ej to linie poza palcami czyli wzbudnik i reszta
v1d11_l7=newreg; Line(v1d11_l7) = {v1d11_p7,v1d11_p8};
v1d11_l8=newreg; Line(v1d11_l8) = {v1d11_p8,v1d11_p9};
v1d11_l9=newreg; Line(v1d11_l9) = {v1d11_p9,v1d11_p10};
v1d11_l10=newreg; Line(v1d11_l10) = {v1d11_p10,v1d11_p11};
v1d11_l11=newreg; Line(v1d11_l11) = {v1d11_p11,v1d11_p12};
//wyciaganie w pionie
v1d12_e1[]= Extrude {0,ys_d,0} {Point{v1d11_p1};};
v1d12_e2[]= Extrude {0,yzapas_d,0} {Point{v1d11_p2};};
v1d12_e3[]= Extrude {0,ynac_d,0} {Point{v1d11_p3};};
//x[]=Boundary{Line{v1d12_e1[0] };};
//Printf("punkty !!! %g , %g, %g, %g,%g , %g, %g, %g" ,
x[0],x[1],x[2],x[3],x[4],x[5],x[6],x[7]);
//wyciaganie w górê reszty dna o sta³ej wysokosci y
//mesh
//layers
v1d12_e4[]= Extrude {0,yz_d,0}
{Line{v1d11_l4, v1d11_l5, v1d11_l6, v1d11_l7, v1d11_l8,v1d11_l9,
v1d11_l10, v1d11_l11 };
Layers {md12b-1};Recombine;};
//linie po wyciagnieciu czyli przod (pierwsza warstwa)
// powyzej czyli (dno rzeczywstego tygla)
//czesc o nierownej wysokosci
p[]=Boundary{Line{v1d12_e1[0], v1d12_e2[0], v1d12_e3[0],v1d12_e4[2] };};
//Printf("punkty !!! %g , %g, %g, %g,%g , %g, %g, %g" ,
p[0],p[1],p[2],p[3],p[4],p[5],p[6],p[7]);
v1d21_l1=newreg; Line(v1d21_l1) = {p[1],p[3]};
v1d21_l2=newreg; Line(v1d21_l2) = {p[3],p[5]};
v1d21_l3=newreg; Line(v1d21_l3) = {p[5],p[7]};
//boczne pow. o nierównym wymiarze y
v1d1_L1=newreg; Line Loop (v1d1_L1) = { v1d11_l1, v1d12_e2[0], -v1d21_l1, -
v1d12_e1[0] };
v1d1_s1 = newreg; Plane Surface (v1d1_s1) = {v1d1_L1};
v1d1_L2=newreg; Line Loop (v1d1_L2) = { v1d11_l2, v1d12_e3[0], -v1d21_l2, -
v1d12_e2[0] };
v1d1_s2 = newreg; Plane Surface (v1d1_s2) = {v1d1_L2};
v1d1_L3=newreg; Line Loop (v1d1_L3) = { v1d11_l3, v1d12_e4[2] , -v1d21_l3, -
v1d12_e3[0] };
v1d1_s3 = newreg; Plane Surface (v1d1_s3) = {v1d1_L3};
v1d_e1[] = Extrude{ {0,1,0}, {0,0,0}, fi1_x_d}
{Surface{ v1d1_s1,
v1d1_s2, v1d1_s3, v1d12_e4[1], v1d12_e4[5],v1d12_e4[9],
v1d12_e4[13], v1d12_e4[17], v1d12_e4[21],v1d12_e4[25],
v1d12_e4[29] };
Layers{md1a}; Recombine; };
//grubosc scianki wzbudnika
//Printf("elementy %g , %g, %g, %g,%g , %g, %g, %g, %g" ,
v1d_e1[0],v1d_e1[1],v1d_e1[2],v1d_e1[3],v1d_e1[4],v1d_e1[5],v1d_e1[6],v1d_e1[7],v1d_e1[8]);
v1d_e2[] = Extrude{ {0,1,0}, {0,0,0}, fi2_x_d}
{Surface{ v1d_e1[0],
v1d_e1[5],v1d_e1[11],v1d_e1[17],v1d_e1[23],
v1d_e1[29],v1d_e1[35],v1d_e1[41],v1d_e1[47],v1d_e1[53],v1d_e1[59] };
Layers{md1b}; Recombine; };
v1d_e3[] = Extrude{ {0,1,0}, {0,0,0}, alfa-(fi2_x_d+fi1_x_d)}
{Surface{ v1d_e2[0],
v1d_e2[5],v1d_e2[11],v1d_e2[17],v1d_e2[23],
v1d_e2[29],v1d_e2[35],v1d_e2[41],v1d_e2[47],v1d_e2[53],v1d_e2[59] };
Layers{md1c}; Recombine; };
v1d_e4[] = Extrude{ {0,1,0}, {0,0,0}, fi2_x_d}
{Surface{ v1d_e3[0],
v1d_e3[5],v1d_e3[11],v1d_e3[17],v1d_e3[23],
v1d_e3[29],v1d_e3[35],v1d_e3[41],v1d_e3[47],v1d_e3[53],v1d_e3[59] };
Layers{md1b}; Recombine; };
v1d_e5[] = Extrude{ {0,1,0}, {0,0,0}, fi1_x_d}
{Surface{ v1d_e4[0],
v1d_e4[5],v1d_e4[11],v1d_e4[17],v1d_e4[23],
v1d_e4[29],v1d_e4[35],v1d_e4[41],v1d_e4[47],v1d_e4[53],v1d_e4[59] };
Layers{md1a}; Recombine; };
//zrobione ca³e dno
//======================
//++++++++++++++++++++++
//definicja objêtosci i powierzchni
//objetosci
//miedz
Physical Volume(dno_tygla)={ v1d_e1[1], v1d_e2[1],v1d_e3[1], v1d_e4[1],
v1d_e5[1],
v1d_e1[6], v1d_e2[6],v1d_e3[6], v1d_e4[6],
v1d_e5[6],
v1d_e2[12],v1d_e3[12], v1d_e4[12],
v1d_e2[36],v1d_e3[36], v1d_e4[36],
v1d_e2[42],v1d_e3[42], v1d_e4[42],
v1d_e2[48],v1d_e3[48], v1d_e4[48],
v1d_e2[54],v1d_e3[54], v1d_e4[54],
v1d_e2[60],v1d_e3[60], v1d_e4[60] };
//powietrze
Physical Volume(szcz_dna)={ v1d_e1[12], v1d_e5[12]};
Physical Volume(szcz_dna_z)={v1d_e1[36], v1d_e5[36],
v1d_e1[42], v1d_e5[42],
v1d_e1[48], v1d_e5[48],
v1d_e1[54], v1d_e5[54],
v1d_e1[60], v1d_e5[60]};
//powierzchnie przednie z=0 i tylne dna
Physical Surface(s_dno_p)={v1d1_s1,v1d1_s2,v1d1_s3,
v1d12_e4[13],v1d12_e4[17],v1d12_e4[21],v1d12_e4[25],v1d12_e4[29] };
Physical Surface(s_dno_t)={v1d_e5[0],v1d_e5[5],v1d_e5[11],
v1d_e5[35],v1d_e5[41],v1d_e5[47],v1d_e5[53],v1d_e5[59]};
//objetosci palca w dnie
//miedz
Physical Volume(palce_dno)={ v1d_e2[18],v1d_e3[18], v1d_e4[18],
v1d_e2[24],v1d_e3[24], v1d_e4[24],
v1d_e2[30],v1d_e3[30], v1d_e4[30]};
//powietrze
Physical Volume(szcz_palcow_dno)={ v1d_e1[18], v1d_e5[18],
v1d_e1[24], v1d_e5[24],
v1d_e1[30], v1d_e5[30]};
//powierzchnie przednie z=0 i tylne palcow w dnie
Physical Surface(s_palce_dno_p)={ v1d12_e4[1],v1d12_e4[5],v1d12_e4[9]};
Physical Surface(s_palce_dno_t)={ v1d_e5[17],v1d_e5[23],v1d_e5[29] };
//++++++++++++++++++++++
//mesh
//po promieniu
//wnetrze
Transfinite Line{ v1d11_l1} = md1d;
Transfinite Line{ v1d21_l1} = md1d;
l[]=Boundary{Surface{v1d_e1[0],v1d_e2[0],v1d_e3[0],v1d_e4[0],v1d_e5[0]};};
For m In {0:18:4}
Transfinite Line{ l[m]} = md1d;
EndFor
For m In {2:18:4}
Transfinite Line{ -l[m]} = md1d;
EndFor
//zapas
Transfinite Line{ v1d11_l2} = md1e;
Transfinite Line{ v1d21_l2} = md1e;
l[]=Boundary{Surface{v1d_e1[5],v1d_e2[5],v1d_e3[5],v1d_e4[5],v1d_e5[5]};};
For m In {0:18:4}
Transfinite Line{ l[m]} = md1e;
EndFor
For m In {2:18:4}
Transfinite Line{ -l[m]} = md1e;
EndFor
//naciecie
Transfinite Line{ v1d11_l3} = md1f;
Transfinite Line{ v1d21_l3} = md1f;
l[]=Boundary{Surface{v1d_e1[11],v1d_e2[11],v1d_e3[11],v1d_e4[11],v1d_e5[11]};};
For m In {0:18:4}
Transfinite Line{ l[m]} = md1f;
EndFor
For m In {2:18:4}
Transfinite Line{ -l[m]} = md1f;
EndFor
//przod palca
Transfinite Line{ v1d11_l4} = md1g;
Transfinite Line{ v1d12_e4[0]} = md1g;
l[]=Boundary{Surface{v1d_e1[17],v1d_e2[17],v1d_e3[17],v1d_e4[17],v1d_e5[17]};};
For m In {1:19:4}
Transfinite Line{ l[m]} = md1g;
EndFor
For m In {3:19:4}
Transfinite Line{ -l[m]} = md1g;
EndFor
//srodek palca
Transfinite Line{ v1d11_l5} = md1h;
Transfinite Line{ v1d12_e4[4]} = md1h;
l[]=Boundary{Surface{v1d_e1[23],v1d_e2[23],v1d_e3[23],v1d_e4[23],v1d_e5[23]};};
For m In {1:19:4}
Transfinite Line{ l[m]} = md1h;
EndFor
For m In {3:19:4}
Transfinite Line{ -l[m]} = md1h;
EndFor
//tyl palca
Transfinite Line{ v1d11_l6} = md1g;
Transfinite Line{ v1d12_e4[8]} = md1g;
l[]=Boundary{Surface{v1d_e1[29],v1d_e2[29],v1d_e3[29],v1d_e4[29],v1d_e5[29]};};
For m In {1:19:4}
Transfinite Line{ l[m]} = md1g;
EndFor
For m In {3:19:4}
Transfinite Line{ -l[m]} = md1g;
EndFor
//szczelina palec wzbudnik
Transfinite Line{ v1d11_l7} = md1i;
Transfinite Line{ v1d12_e4[12]} = md1i;
l[]=Boundary{Surface{v1d_e1[35],v1d_e2[35],v1d_e3[35],v1d_e4[35],v1d_e5[35]};};
For m In {1:19:4}
Transfinite Line{ l[m]} = md1i;
EndFor
For m In {3:19:4}
Transfinite Line{ -l[m]} = md1i;
EndFor
//scianka wzbudnika wewnetrzna
Transfinite Line{ v1d11_l8} = md1j;
Transfinite Line{ v1d12_e4[16]} = md1j;
l[]=Boundary{Surface{v1d_e1[41],v1d_e2[41],v1d_e3[41],v1d_e4[41],v1d_e5[41]};};
For m In {1:19:4}
Transfinite Line{ l[m]} = md1j;
EndFor
For m In {3:19:4}
Transfinite Line{ -l[m]} = md1j;
EndFor
//otwor wzbudnika
Transfinite Line{ v1d11_l9} = md1k;
Transfinite Line{ v1d12_e4[20]} = md1k;
l[]=Boundary{Surface{v1d_e1[47],v1d_e2[47],v1d_e3[47],v1d_e4[47],v1d_e5[47]};};
For m In {1:19:4}
Transfinite Line{ l[m]} = md1k;
EndFor
For m In {3:19:4}
Transfinite Line{ -l[m]} = md1k;
EndFor
//scianka wzbudnika zewnetrzna
Transfinite Line{ v1d11_l10} = md1j;
Transfinite Line{ v1d12_e4[24]} = md1j;
l[]=Boundary{Surface{v1d_e1[53],v1d_e2[53],v1d_e3[53],v1d_e4[53],v1d_e5[53]};};
For m In {1:19:4}
Transfinite Line{ l[m]} = md1j;
EndFor
For m In {3:19:4}
Transfinite Line{ -l[m]} = md1j;
EndFor
//do konca dna
Transfinite Line{ v1d11_l11} = md1l;
Transfinite Line{ v1d12_e4[28]} = md1l;
l[]=Boundary{Surface{v1d_e1[59],v1d_e2[59],v1d_e3[59],v1d_e4[59],v1d_e5[59]};};
For m In {1:19:4}
// Printf("line %g ",l[m]);
Transfinite Line{ l[m]} = md1l;
EndFor
For m In {3:19:4}
// Printf("line %g ",l[m]);
Transfinite Line{ -l[m]} = md1l;
EndFor
//koniec podzialow po x
//po wysokosci
Transfinite Line{v1d12_e1[0]}=md12a;
Transfinite Line{v1d12_e2[0],v1d12_e3[0] } = md12b;
//(reszta jest zrobiona na pocz¹tku w extrude)
//prostokat od srodka dna
Recombine Surface{v1d1_s1};
// 2 nastepny
Transfinite Surface{v1d1_s2,v1d1_s3};
//pozostale s¹ wyci¹gniête po wysokosci
//uwaga wyci¹gniecie = zwykly podzial -1 !!!
*/
_______________________________________________
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
