PTT數位生活區 / Cad_Cae (電腦輔助設計)

[問題] CAE"四面體"與"角柱形"網格之品質控制?

看板Cad_Cae (電腦輔助設計)作者 (阿泰斯)時間10年前 (2014/04/24 14:07), 編輯推噓11(11推 0噓 80→)
留言91則, 5人參與, 最新討論串1/1
我是用學校版的Comsol，有很多個版本，4.0~4.3都有 在網格的時候，可以用六面體網格、四面體網格或角柱形網格 六面體網格比較好控制網格品質，只要每個網格形狀都接近立方體的話，網格品質就會 比較好 四面體網格則不知道如何控制網格品質，因為許多模型沒辦法用六面體網格，所以會改 用四面體網格 而角柱形網格是五面體網格，就是有一面是三角形，然後再拉伸，變成角柱形的網格， 這種網格不知道計算出來的結果會不會準確? 想問一下四面體網格如何控制網格品質? 因為很奇怪的是，我用六面體網格，網格品質是1，網格數大約10萬，記憶體要吃到快 100GB，而我用四面體網格，網格數約20萬，記憶體吃不到50GB，所以依照 "硬體的消耗量"，勢必要研究一下四面體網格的網格品質如何控制 懇求大家回應!!! -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 175.180.205.96 ※ 文章網址: http://www.ptt.cc/bbs/Cad_Cae/M.1398319631.A.CBC.html
推
04/24 23:08, , 1F
記憶體100GB=16GB*8? 元素本身同類型? 是否有中間節點?
04/24 23:08, 1F
→
04/24 23:14, , 2F
沒有節點，每個元素都是正方體，大概10萬個。
04/24 23:14, 2F
→
04/24 23:15, , 3F
物理量是壓電元件，頻域分析。
04/24 23:15, 3F
→
04/24 23:16, , 4F
請問一下，正方體網格數10萬上下，記憶體吃到100GB很不
04/24 23:16, 4F
→
04/24 23:16, , 5F
正常嗎? 怎會有此疑問?
04/24 23:16, 5F
推
04/24 23:21, , 6F
10萬個網格吃到100g記憶體很不合理.....
04/24 23:21, 6F
→
04/24 23:22, , 7F
我覺得一定有什麼地方搞錯了。
04/24 23:22, 7F
→
04/24 23:27, , 8F
零件就很簡單，但是很小，約3mm*1mm*0.5mm，是否是零件
04/24 23:27, 8F
→
04/24 23:27, , 9F
大小的關係呢? 零件是長方體。
04/24 23:27, 9F
→
04/24 23:34, , 10F
是在"線性求解器"的進度超過30%的時候，記憶體才會爆衝
04/24 23:34, 10F
→
04/24 23:34, , 11F
到100GB左右，很奇怪@@
04/24 23:34, 11F
→
04/24 23:37, , 12F
"線性求解器"的進度不是整體進度，"參數求解器"的進度
04/24 23:37, 12F
→
04/24 23:37, , 13F
才是整體進度。
04/24 23:37, 13F
→
04/24 23:38, , 14F
大大們是否有嘗試用壓電元件物理量，頻域分析，測試我
04/24 23:38, 14F
→
04/24 23:39, , 15F
上面說的尺寸呢? 3mm*1mm*0.5mm，用映射+掃描，把厚度
04/24 23:39, 15F
→
04/24 23:39, , 16F
弄成6格就好，網格數量就會很多了，這時候分析的話，
04/24 23:39, 16F
→
04/24 23:39, , 17F
記憶體會吃到將近100GB
04/24 23:39, 17F
→
04/24 23:40, , 18F
我是說"線性求解器"超過30%的時候。
04/24 23:40, 18F
→
04/24 23:42, , 19F
抱歉，上面尺寸說錯，是3mm*1mm*0.05mm
04/24 23:42, 19F
→
04/24 23:43, , 20F
每個元素都是正六面體喔! 正確尺寸是:3mm*1mm*0.05mm
04/24 23:43, 20F
→
04/24 23:47, , 21F
抱歉!!!尺寸又弄錯，是1.5mm*1mm*0.05mm
04/24 23:47, 21F
→
04/24 23:52, , 22F
這是吃的記憶體的截圖: http://ppt.cc/wotp
04/24 23:52, 22F
推
04/25 00:49, , 23F
0.05mm弄6層代表每一層的厚度只有0.008mm，你又做到1:1，這
04/25 00:49, 23F
→
04/25 00:50, , 24F
樣mesh會變超密，我是建議不需要用到這麼密，可以先試著把
04/25 00:50, 24F
→
04/25 00:50, , 25F
長寬比拉大一點(5:1)，然後慢慢縮小看收斂性如何。
04/25 00:50, 25F
→
04/25 00:51, , 26F
但是以10萬個element的K矩陣大小，100GB還是太誇張....
04/25 00:51, 26F
→
04/25 11:22, , 27F
感謝建議!!! 但是我是說真的!!! 六面體網格的情況下，
04/25 11:22, 27F
→
04/25 11:23, , 28F
我上面說的那個模型會吃到100GB!!!
04/25 11:23, 28F
→
04/25 11:23, , 29F
設定方面，真的沒有什麼好錯或不錯的，就是模型畫好，
04/25 11:23, 29F
→
04/25 11:24, , 30F
給材料，給電極，給網格，給頻域，然後分析...
04/25 11:24, 30F
→
04/25 11:25, , 31F
已經測試很久一段時間了，要是記憶體沒吃這麼多的話，
04/25 11:25, 31F
→
04/25 11:25, , 32F
早就測試出一個滿意的結果了。
04/25 11:25, 32F
→
04/25 11:26, , 33F
六面體網格跟四面體網格真的差很多，四面體網格數20萬
04/25 11:26, 33F
→
04/25 11:26, , 34F
，記憶體吃大概10GB...
04/25 11:26, 34F
→
04/25 11:26, , 35F
跟六面體網格真的是天朗之別
04/25 11:26, 35F
→
04/25 16:12, , 36F
對了! 請問各位大大的收斂性都怎麼看? 是看特定某個值
04/25 16:12, 36F
→
04/25 16:12, , 37F
隨著網格變化的變化嗎?
04/25 16:12, 37F
→
04/25 16:13, , 38F
但是在頻域分析裡面，頻率間隔不同，收斂性也會不同耶
04/25 16:13, 38F
→
04/25 16:13, , 39F
能用特徵頻率去看收斂性嗎?
04/25 16:13, 39F
推
04/25 16:43, , 40F
你單位有看錯嗎? 要插滿 100GB 的記憶體也不容易
04/25 16:43, 40F
→
04/25 17:20, , 41F
沒看錯，CAE是公尺，模型是mm，故我輸入1.5e-3
04/25 17:20, 41F
→
04/25 17:21, , 42F
、1e-3、0.05e-3
04/25 17:21, 42F
→
04/25 17:24, , 43F
至於電腦是學校資產
04/25 17:24, 43F
→
04/25 17:42, , 44F
ptt有常用的免空嗎？我把存檔放上來給大家看@@
04/25 17:42, 44F
推
04/25 17:44, , 45F
我指的單位是記憶體
04/25 17:44, 45F
→
04/25 17:50, , 46F
這是計算時的截圖@@ http://ppt.cc/wotp
04/25 17:50, 46F
→
04/25 17:50, , 47F
看工作管理員也是吃大概這麼多
04/25 17:50, 47F
推
04/25 18:14, , 48F
能截一下工作管理員的圖嗎
04/25 18:14, 48F
推
04/25 18:16, , 49F
我覺得應該不會到 100GB
04/25 18:16, 49F
→
04/25 18:17, , 50F
抱歉，請稍等一下，我正在計算其他的東西@@
04/25 18:17, 50F
→
04/25 19:42, , 51F
抱歉@@ 可能明天才會PO圖了，目前一直在測試四面體網格
04/25 19:42, 51F
→
04/25 19:42, , 52F
的收斂性。
04/25 19:42, 52F
→
04/25 23:06, , 53F
是mesh完就100GB?還是計算到某個點變100GB?如果是後者，我認
04/25 23:06, 53F
→
04/25 23:08, , 54F
為應該解發散了吧。這種狀況不正常，建議找原廠來看看。
04/25 23:08, 54F
→
04/25 23:09, , 55F
收斂性可以做收斂測試，同樣的條件算頻域的值在多密的網格下
04/25 23:09, 55F
→
04/25 23:11, , 56F
解的差異才會趨近最小值，該模型就是已收斂，再改密網格也不
04/25 23:11, 56F
→
04/25 23:11, , 57F
會對解的精確性有明顯的幫助。
04/25 23:11, 57F
→
04/25 23:44, , 58F
是後者，是線性求解器求解超過30%的時候，記憶體會爆衝
04/25 23:44, 58F
→
04/25 23:45, , 59F
任何模型分析都會這樣@@，比如說我分析一個簡單的模型
04/25 23:45, 59F
→
04/25 23:45, , 60F
線性求解器進度低於30%時，只吃2GB，但是超過30%後，
04/25 23:45, 60F
→
04/25 23:46, , 61F
記憶體可以吃到5GB或更高。
04/25 23:46, 61F
→
04/25 23:46, , 62F
線性求解器的進度不是總進度，是參數求解器的進度才是
04/25 23:46, 62F
→
04/25 23:46, , 63F
總進度。
04/25 23:46, 63F
→
04/25 23:47, , 64F
補充，剛剛說線性求解器進度超過30%後，記憶體會爆衝，
04/25 23:47, 64F
→
04/25 23:47, , 65F
但是線性求解器進度衝到100%後，進度會歸零，這時候
04/25 23:47, 65F
→
04/25 23:48, , 66F
記憶體的使用量又會降低到原本。
04/25 23:48, 66F
→
04/25 23:49, , 67F
上面是指頻域分析的時候，求解器進度有分參數與線性。
04/25 23:49, 67F
推
04/26 01:40, , 68F
記憶體量突然大伏增加應該是開始在解k矩陣的關係。
04/26 01:40, 68F
推
04/26 01:42, , 69F
此時因為需要多次疊代解逆矩陣，所以會使用大量記憶體，
04/26 01:42, 69F
→
04/26 01:43, , 70F
但是在線性運算時，此動作需要的資源和node的自由度有比較
04/26 01:43, 70F
→
04/26 01:44, , 71F
明顯的關係，10萬個HEX. elements，如果都是一階，頂多80萬
04/26 01:44, 71F
→
04/26 01:45, , 72F
個node，不應該吃這麼多資源才是。
04/26 01:45, 72F
→
04/26 02:12, , 73F
剛看了一下，有4524052個自由度量...
04/26 02:12, 73F
→
04/26 02:15, , 74F
04/26 02:15, 74F
→
04/26 02:25, , 75F
自由度是由甚麼來控制的@@?
04/26 02:25, 75F
→
04/26 05:12, , 76F
對了！請問一下您說的“一階”是什麼？研究好久研究不
04/26 05:12, 76F
→
04/26 05:12, , 77F
出來。
04/26 05:12, 77F
→
04/26 13:07, , 78F
請問一下，網格的"階數"是去那裡調整? 網路上有找到
04/26 13:07, 78F
→
04/26 13:08, , 79F
階數與自由度的關係，但是沒找到調整階數的地方@@
04/26 13:08, 79F
→
04/26 13:17, , 80F
已測試"幾何形狀階數"，完全無關。
04/26 13:17, 80F
→
04/26 13:51, , 81F
剛剛用"固體力學"下去測試同樣模型，自由度約360萬，
04/26 13:51, 81F
→
04/26 13:52, , 82F
也就是說不同物理量，自由度不同，請問您是否有測試過
04/26 13:52, 82F
→
04/26 13:52, , 83F
用壓電元件測試1.5*1*0.05mm這麼模型?用六面體網格，Y
04/26 13:52, 83F
→
04/26 13:53, , 84F
厚度切六格，網格數約13萬1:1網格。
04/26 13:53, 84F
→
04/26 18:53, , 85F
補上檔案: http://ppt.cc/-jE1
04/26 18:53, 85F
→
04/26 18:54, , 86F
上面那個檔案是Comsol 4.3b存檔。
04/26 18:54, 86F
推
04/26 22:37, , 87F
我沒用過comsol，所以檔案的部分可能要請其他有經驗的板友
04/26 22:37, 87F
→
04/26 22:37, , 88F
幫忙，其他的問題我回在一篇文章好了。
04/26 22:37, 88F
→
04/26 23:27, , 89F
感謝!!!
04/26 23:27, 89F
推
04/28 20:32, , 90F
看了您的檔案，單單一個壓電板，網格不需要切的那麼密吧!
04/28 20:32, 90F
→
04/28 21:26, , 91F
感謝回應!!! 但是切很粗的話，頻率沒有收斂阿@@
04/28 21:26, 91F
‣ 更多 candy88257 的文章
文章代碼(AID): #1JMAeFoy (Cad_Cae)
Cad_Cae 近期熱門文章
更多 近期熱門文章 >>
PTT數位生活區 即時熱門文章
更多 即時熱門文章 >>
‣ 更多 candy88257 的文章
文章代碼(AID): #1JMAeFoy (Cad_Cae)