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

嵇康

于 2014-02-11 14:41:49 发布

阅读量4.4k

点赞数 1

CC 4.0 BY-SA版权

分类专栏：数值分析&FEM

本文链接：https://blog.youkuaiyun.com/conjimmy/article/details/19071871

数值分析&FEM 专栏收录该内容

26 篇文章

订阅专栏

本文介绍了一个使用DistMesh2D进行网格划分的例子，包括定义边界和内部区域的函数，以及生成网格的具体步骤。数据文件包含了节点坐标及材料属性等详细信息。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

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