>>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