[Gmsh] How to write Plugin crack command for post processing
LU Xiaoxin
Xiaoxin.Lu at u-pem.fr
Thu Sep 29 17:42:33 CEST 2016
Hello,
I am a PhD student using GMSH to generate crack faces. However, I could
not find any example for this kind of post processing.
My .geo file is attached below, I put 'plugin crack' command inside but
it did not work, I was wondering how to generate a post process file for
this command.
Thank you for your help.
best regards,
Xiaoxin
//start of the file
lc1 = 3;
lc2 = 3;
// We proceed by defining some elementary entities describing a cube:
Point(1) = {0,0,0,lc1};
Point(2) = {40,0,0,lc1};
Point(3) = {40,40,0,lc1};
Point(4) = {0,40,0,lc1};
Point(5) = {0,0,40,lc1};
Point(6) = {40,0,40,lc1};
Point(7) = {40,40,40,lc1};
Point(8) = {0,40,40,lc1};
Line(30)={1,2};
Line(31)={2,3};
Line(32)={3,4};
Line(33)={4,1};
Line(34)={5,6};
Line(35)={6,7};
Line(36)={7,8};
Line(37)={8,5};
Line(38)={1,5};
Line(39)={2,6};
Line(40)={3,7};
Line(41)={4,8};
Line Loop(42)={30,31,32,33};
Plane Surface(14)={42};
Line Loop(43)={-37,-36,-35,-34};
Plane Surface(16)={43};
Line Loop(44)={38,34,-39,-30};
Plane Surface(18)={44};
Line Loop(45)={-32,40,36,-41};
Plane Surface(20)={45};
Line Loop(46)={-33,41,37,-38};
Plane Surface(22)={46};
Line Loop(47)={39,35,-40,-31};
Plane Surface(24)={47};
vol1[]=Extrude{0,0,40}{Surface{14};};
Delete {Volume{vol1[1]};}
vol2[]=Extrude{0,40,0}{Surface{18};};
Delete {Volume{vol2[1]};}
vol3[]=Extrude{40,0,0}{Surface{22};};
Delete {Volume{vol3[1]};}
Physical Surface(1)={14};
Physical Surface(2)={vol1[0]};
Physical Surface(3)={18};
Physical Surface(4)={vol2[0]};
Physical Surface(5)={22};
Physical Surface(6)={vol3[0]};
Surface Loop(1)={14,vol2[0],vol3[0],18,22,vol1[0]};
Volume(1)={1};
Physical Volume(1)={1};
//n=1; add a crack surface.
Point(5000)={1.645000e+01,7.466800e+00,5.890600e+00,lc2};
Point(5001)={1.728000e+01,2.112300e+01,1.204100e+01,lc2};
Point(5002)={2.823700e+01,1.636900e+01,2.111600e+01,lc2};
Point(5003)={2.740700e+01,2.713500e+00,1.496600e+01,lc2};
Line(5004)={5000,5001};
Line(5005)={5001,5002};
Line(5006)={5002,5003};
Line(5007)={5003,5000};
Line Loop(5008)={5004,5005,5006,5007};
Plane Surface(5009)={5008};
Physical Surface(11)={5009};
Surface{5009} In Volume{1};
Plugin(Crack).Dimension = 2 ;
Plugin(Crack).PhysicalGroup = 11 ;
Plugin(Crack).OpenBoundaryPhysicalGroup = 0 ;
Plugin(Crack).Run ;
//end of the file
More information about the gmsh
mailing list