distmesh划分网格之正方形内含椭圆

>>fh=@(p) 0.05+0.3*dellipse(p,[3,1]);
>>fd=@(p) ddiff(drectangle(p,-10,10,-10,10),dellipse(p,[3,1]));

>>[p,t]=distmesh2d(fd,fh,0.05,[-10,-10;10,10],[-10,-10;-10,10;10,-10;10,10]);

 

 

 

数据文件 exam4_2_ellipse_crack.dat3如下：

1283
1 -10.000000 -10.000000
2 -10.000000 10.000000
3 10.000000 -10.000000
4 10.000000 10.000000
5 -10.000000 -1.328691
6 -10.000000 -3.764229
7 -10.000000 -6.512389
8 -6.035232 -10.000000
9 -10.000000 1.010601
10 -10.000000 3.444287
11 -7.392375 -7.690254
12 -2.708587 -10.000000

。。。

5
1   10E6      0.3   1.0  1.0e3
2   2.60e9    0.20  1.0  23.0e3
3   2.85e10   0.20  1.0  25.0e3
4   1.85e10   0.20  1.0  23.0e3
5   2.85e10   0.20  1.0  22.0e3
8
2     2 0
4     2 0
690   2 0
318   2 0
1278  2 0
199   2 0
13    2 0
62    2 0