作者candy88257 (阿泰斯)
看板Cad_Cae
标题[问题] CAE"四面体"与"角柱形"网格之品质控制?
时间Thu Apr 24 14:07:08 2014
我是用学校版的Comsol，有很多个版本，4.0~4.3都有
在网格的时候，可以用六面体网格、四面体网格或角柱形网格
六面体网格比较好控制网格品质，只要每个网格形状都接近立方体的话，网格品质就会
比较好
四面体网格则不知道如何控制网格品质，因为许多模型没办法用六面体网格，所以会改
用四面体网格
而角柱形网格是五面体网格，就是有一面是三角形，然後再拉伸，变成角柱形的网格，
这种网格不知道计算出来的结果会不会准确?
想问一下四面体网格如何控制网格品质?
因为很奇怪的是，我用六面体网格，网格品质是1，网格数大约10万，记忆体要吃到快
100GB，而我用四面体网格，网格数约20万，记忆体吃不到50GB，所以依照
"硬体的消耗量"，势必要研究一下四面体网格的网格品质如何控制
恳求大家回应!!!
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 175.180.205.96
※ 文章网址: http://webptt.com/cn.aspx?n=bbs/Cad_Cae/M.1398319631.A.CBC.html
1F：推 fgj:记忆体100GB=16GB*8? 元素本身同类型? 是否有中间节点? 04/24 23:08
2F：→ candy88257:没有节点，每个元素都是正方体，大概10万个。 04/24 23:14
3F：→ candy88257:物理量是压电元件，频域分析。 04/24 23:15
4F：→ candy88257:请问一下，正方体网格数10万上下，记忆体吃到100GB很不 04/24 23:16
5F：→ candy88257:正常吗? 怎会有此疑问? 04/24 23:16
6F：推 gamer:10万个网格吃到100g记忆体很不合理..... 04/24 23:21
7F：→ gamer:我觉得一定有什麽地方搞错了。 04/24 23:22
8F：→ candy88257:零件就很简单，但是很小，约3mm*1mm*0.5mm，是否是零件 04/24 23:27
9F：→ candy88257:大小的关系呢? 零件是长方体。 04/24 23:27
10F：→ candy88257:是在"线性求解器"的进度超过30%的时候，记忆体才会爆冲 04/24 23:34
11F：→ candy88257:到100GB左右，很奇怪@@ 04/24 23:34
12F：→ candy88257:"线性求解器"的进度不是整体进度，"参数求解器"的进度 04/24 23:37
13F：→ candy88257:才是整体进度。 04/24 23:37
14F：→ candy88257:大大们是否有尝试用压电元件物理量，频域分析，测试我 04/24 23:38
15F：→ candy88257:上面说的尺寸呢? 3mm*1mm*0.5mm，用映射+扫描，把厚度 04/24 23:39
16F：→ candy88257:弄成6格就好，网格数量就会很多了，这时候分析的话， 04/24 23:39
17F：→ candy88257:记忆体会吃到将近100GB 04/24 23:39
18F：→ candy88257:我是说"线性求解器"超过30%的时候。 04/24 23:40
19F：→ candy88257:抱歉，上面尺寸说错，是3mm*1mm*0.05mm 04/24 23:42
20F：→ candy88257:每个元素都是正六面体喔! 正确尺寸是:3mm*1mm*0.05mm 04/24 23:43
21F：→ candy88257:抱歉!!!尺寸又弄错，是1.5mm*1mm*0.05mm 04/24 23:47
23F：推 gamer:0.05mm弄6层代表每一层的厚度只有0.008mm，你又做到1:1，这 04/25 00:49
24F：→ gamer:样mesh会变超密，我是建议不需要用到这麽密，可以先试着把 04/25 00:50
25F：→ gamer:长宽比拉大一点(5:1)，然後慢慢缩小看收敛性如何。 04/25 00:50
26F：→ gamer:但是以10万个element的K矩阵大小，100GB还是太夸张.... 04/25 00:51
27F：→ candy88257:感谢建议!!! 但是我是说真的!!! 六面体网格的情况下， 04/25 11:22
28F：→ candy88257:我上面说的那个模型会吃到100GB!!! 04/25 11:23
29F：→ candy88257:设定方面，真的没有什麽好错或不错的，就是模型画好， 04/25 11:23
30F：→ candy88257:给材料，给电极，给网格，给频域，然後分析... 04/25 11:24
31F：→ candy88257:已经测试很久一段时间了，要是记忆体没吃这麽多的话， 04/25 11:25
32F：→ candy88257:早就测试出一个满意的结果了。 04/25 11:25
33F：→ candy88257:六面体网格跟四面体网格真的差很多，四面体网格数20万 04/25 11:26
34F：→ candy88257:，记忆体吃大概10GB... 04/25 11:26
35F：→ candy88257:跟六面体网格真的是天朗之别 04/25 11:26
36F：→ candy88257:对了! 请问各位大大的收敛性都怎麽看? 是看特定某个值 04/25 16:12
37F：→ candy88257:随着网格变化的变化吗? 04/25 16:12
38F：→ candy88257:但是在频域分析里面，频率间隔不同，收敛性也会不同耶 04/25 16:13
39F：→ candy88257:能用特徵频率去看收敛性吗? 04/25 16:13
40F：推 largesperm:你单位有看错吗? 要插满 100GB 的记忆体也不容易 04/25 16:43
41F：→ candy88257:没看错，CAE是公尺，模型是mm，故我输入1.5e-3 04/25 17:20
42F：→ candy88257:、1e-3、0.05e-3 04/25 17:21
43F：→ candy88257:至於电脑是学校资产 04/25 17:24
44F：→ candy88257:ptt有常用的免空吗？我把存档放上来给大家看@@ 04/25 17:42
45F：推 largesperm:我指的单位是记忆体 04/25 17:44
47F：→ candy88257:看工作管理员也是吃大概这麽多 04/25 17:50
48F：推 largesperm:能截一下工作管理员的图吗 04/25 18:14
49F：推 largesperm:我觉得应该不会到 100GB 04/25 18:16
50F：→ candy88257:抱歉，请稍等一下，我正在计算其他的东西@@ 04/25 18:17
51F：→ candy88257:抱歉@@ 可能明天才会PO图了，目前一直在测试四面体网格 04/25 19:42
52F：→ candy88257:的收敛性。 04/25 19:42
53F：→ gamer:是mesh完就100GB?还是计算到某个点变100GB?如果是後者，我认 04/25 23:06
54F：→ gamer:为应该解发散了吧。这种状况不正常，建议找原厂来看看。 04/25 23:08
55F：→ gamer:收敛性可以做收敛测试，同样的条件算频域的值在多密的网格下 04/25 23:09
56F：→ gamer:解的差异才会趋近最小值，该模型就是已收敛，再改密网格也不 04/25 23:11
57F：→ gamer:会对解的精确性有明显的帮助。 04/25 23:11
58F：→ candy88257:是後者，是线性求解器求解超过30%的时候，记忆体会爆冲 04/25 23:44
59F：→ candy88257:任何模型分析都会这样@@，比如说我分析一个简单的模型 04/25 23:45
60F：→ candy88257:线性求解器进度低於30%时，只吃2GB，但是超过30%後， 04/25 23:45
61F：→ candy88257:记忆体可以吃到5GB或更高。 04/25 23:46
62F：→ candy88257:线性求解器的进度不是总进度，是参数求解器的进度才是 04/25 23:46
63F：→ candy88257:总进度。 04/25 23:46
64F：→ candy88257:补充，刚刚说线性求解器进度超过30%後，记忆体会爆冲， 04/25 23:47
65F：→ candy88257:但是线性求解器进度冲到100%後，进度会归零，这时候 04/25 23:47
66F：→ candy88257:记忆体的使用量又会降低到原本。 04/25 23:48
67F：→ candy88257:上面是指频域分析的时候，求解器进度有分参数与线性。 04/25 23:49
68F：推 gamer:记忆体量突然大伏增加应该是开始在解k矩阵的关系。 04/26 01:40
69F：推 gamer:此时因为需要多次叠代解逆矩阵，所以会使用大量记忆体， 04/26 01:42
70F：→ gamer:但是在线性运算时，此动作需要的资源和node的自由度有比较 04/26 01:43
71F：→ gamer:明显的关系，10万个HEX. elements，如果都是一阶，顶多80万 04/26 01:44
72F：→ gamer:个node，不应该吃这麽多资源才是。 04/26 01:45
73F：→ candy88257:刚看了一下，有4524052个自由度量... 04/26 02:12
75F：→ candy88257:自由度是由甚麽来控制的@@? 04/26 02:25
76F：→ candy88257:对了！请问一下您说的“一阶”是什麽？研究好久研究不 04/26 05:12
77F：→ candy88257:出来。 04/26 05:12
78F：→ candy88257:请问一下，网格的"阶数"是去那里调整? 网路上有找到 04/26 13:07
79F：→ candy88257:阶数与自由度的关系，但是没找到调整阶数的地方@@ 04/26 13:08
80F：→ candy88257:已测试"几何形状阶数"，完全无关。 04/26 13:17
81F：→ candy88257:刚刚用"固体力学"下去测试同样模型，自由度约360万， 04/26 13:51
82F：→ candy88257:也就是说不同物理量，自由度不同，请问您是否有测试过 04/26 13:52
83F：→ candy88257:用压电元件测试1.5*1*0.05mm这麽模型?用六面体网格，Y 04/26 13:52
84F：→ candy88257:厚度切六格，网格数约13万1:1网格。 04/26 13:53
86F：→ candy88257:上面那个档案是Comsol 4.3b存档。 04/26 18:54
87F：推 gamer:我没用过comsol，所以档案的部分可能要请其他有经验的板友 04/26 22:37
88F：→ gamer:帮忙，其他的问题我回在一篇文章好了。 04/26 22:37
89F：→ candy88257:感谢!!! 04/26 23:27
90F：推 henkey:看了您的档案，单单一个压电板，网格不需要切的那麽密吧! 04/28 20:32
91F：→ candy88257:感谢回应!!! 但是切很粗的话，频率没有收敛阿@@ 04/28 21:26