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

[問題] ansys workbench與comsol自然頻率解不同?

看板Cad_Cae (電腦輔助設計)作者 (阿泰斯)時間10年前 (2014/05/05 12:30), 編輯推噓6(6推 0噓 52→)
留言58則, 4人參與, 最新討論串1/1
題目是想分析一個正方體鐵塊，邊長1mm，密度7860kg/m^3，楊氏係數152*10^9Pa， 普松比0.27，在完全沒有邊界條件的情況下，分析其特徵頻率 comsol裡面的邊界條件中，完全沒設定的情況下，會自動設定整個模型的 "初始位移以及速度(都是0)"和"自由邊界條件" 而ansys我是選modal分析，材料常數自己新增一個，網格弄得跟comsol一模一樣 (4*4*4六面體網格)，邊界條件完全沒設定，然後直接分析其前六個特徵頻率 發現分析結果與comsol完全不同... 接著我測試其中一個面固定 發現comsol與ansys分析結果有點像，但是誤差還是很大...(比如:0.4567Hz和0.4667Hz) 照理講不該有這麼大的誤差@@，不然哪知道comsol與ansys哪個算出來的才是對的...? 請問一下ansys分析特徵頻率有甚麼特別需要注意的地方嗎? -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 175.180.80.71 ※ 文章網址: http://www.ptt.cc/bbs/Cad_Cae/M.1399264256.A.465.html
→
05/05 12:32, , 1F
補充: comsol 4.3b和workbench 13
05/05 12:32, 1F
→
05/05 12:36, , 2F
再補充: comsol有把鐵設定成"線性彈性材料"，ansys
05/05 12:36, 2F
→
05/05 12:36, , 3F
我自己新增材料的時候，只有輸入密度、楊氏係數、普松
05/05 12:36, 3F
→
05/05 12:36, , 4F
比，是這個原因嗎?
05/05 12:36, 4F
→
05/05 12:41, , 5F
補充，這是comsol算出來的特徵頻率:
05/05 12:41, 5F
→
05/05 12:42, , 6F
0.090931、0.107872、0.122555、0.13087、0.133412、
05/05 12:42, 6F
→
05/05 12:42, , 7F
0.143898 Hz。
05/05 12:42, 7F
→
05/05 12:43, , 8F
(無設定任何邊界條件下，除了comsol預設的)
05/05 12:43, 8F
推
05/05 18:42, , 9F
之前的筆記 有點相關 ANSYS用Fixed Support會正確跑出6個模態
05/05 18:42, 9F
→
05/05 18:42, , 10F
用Frictionless Support 前3個模態頻率幾乎是0 等同於剛體運動
05/05 18:42, 10F
→
05/05 18:42, , 11F
方塊本身無變形
05/05 18:42, 11F
推
05/06 00:04, , 12F
應該是前6個，含3個純旋轉的模態都應該是0
05/06 00:04, 12F
推
05/06 01:09, , 13F
是3個沒錯 XY平面被拘束住 Z無法移動 沿X軸.Y軸也無法旋轉
05/06 01:09, 13F
→
05/06 01:10, , 14F
沒Support才是6個0
05/06 01:10, 14F
→
05/06 12:34, , 15F
感謝!!! 我沒加任何邊界條件，comsol算出來前六個有解
05/06 12:34, 15F
→
05/06 12:35, , 16F
但是ansys算出來前六個都是0，而comsol跟ansys算出來
05/06 12:35, 16F
→
05/06 12:36, , 17F
第7個之後都是MHz起跳，數字都有點接近，但是誤差很大
05/06 12:36, 17F
→
05/06 12:37, , 18F
比如:1.53MHz跟1.54MHz，誤差非常大，網格已經弄到
05/06 12:37, 18F
→
05/06 12:37, , 19F
40*40*40的立方體網格了，請問為何會這樣@@?
05/06 12:37, 19F
→
05/06 12:38, , 20F
若是邊界條件問題，comsol裡面只設定"自由邊界條件"、
05/06 12:38, 20F
→
05/06 12:38, , 21F
"初始位移、速度"，可是ansys裡面沒有這些選樣阿!!!
05/06 12:38, 21F
→
05/06 13:31, , 22F
comsol裡面沒設定固定邊界、摩擦邊界...，前六個特徵
05/06 13:31, 22F
→
05/06 13:32, , 23F
頻率卻不是0...，怎會這麼奇怪...?
05/06 13:32, 23F
推
05/06 22:44, , 24F
1.53和1.54誤差有很大嗎@@，0.6%而已耶。
05/06 22:44, 24F
→
05/06 22:45, , 25F
comsol裡面可能有自動weak support的設定，使模型不會有剛體
05/06 22:45, 25F
→
05/06 22:45, , 26F
運動。
05/06 22:45, 26F
→
05/08 20:59, , 27F
感謝!!! 請問一下，您說的weak support是指什麼?
05/08 20:59, 27F
→
05/08 21:00, , 28F
因為comsol邊界條件有"固定邊界"這個選項，所以我沒設
05/08 21:00, 28F
→
05/08 21:00, , 29F
定"固定邊界"，就肯定沒有固定邊界這個條件。
05/08 21:00, 29F
→
05/08 21:01, , 30F
請問有哪些邊界是加了之後，前六個特徵頻率會不為0的?
05/08 21:01, 30F
→
05/08 21:02, , 31F
另外，您說的誤差只有0.6%，感覺好像很少，但是實際
05/08 21:02, 31F
→
05/08 21:03, , 32F
產品，是要看差幾Hz的@@
05/08 21:03, 32F
→
05/08 21:03, , 33F
應該說也是看%數，但是%數要求更精密
05/08 21:03, 33F
→
05/08 21:05, , 34F
比如說產品要求10MHz，我網格A求出10001000Hz
05/08 21:05, 34F
→
05/08 21:06, , 35F
網格B求出10002000Hz
05/08 21:06, 35F
→
05/08 21:07, , 36F
看誤差是相當的小...，9.9e-3%，但是差了1000Hz...
05/08 21:07, 36F
推
05/08 21:41, , 37F
但是這是數值分析，他有一定的誤差在，不會是你想的那樣完美
05/08 21:41, 37F
→
05/08 21:42, , 38F
而0.6%在FEM已經算是很小的誤差了，要再小你的網格可能得切
05/08 21:42, 38F
→
05/08 21:42, , 39F
的非常細才有可能。
05/08 21:42, 39F
→
05/08 21:43, , 40F
weak support就是程式會自動加上一個很小的彈簧支撐，使模型
05/08 21:43, 40F
→
05/08 21:44, , 41F
沒有剛體運動，而因為這個彈簧的k值很小，所以也不會有額外
05/08 21:44, 41F
→
05/08 21:44, , 42F
的反力發生。我不確定COMSOL有沒有，只是我的猜測而已。
05/08 21:44, 42F
→
05/08 22:00, , 43F
感謝!!!今天又對FEM更加了解了!!!
05/08 22:00, 43F
→
05/08 22:04, , 44F
另外，就上面這個例子，鐵塊是邊長1mm正方體，網格是
05/08 22:04, 44F
→
05/08 22:05, , 45F
每邊劃分40格的正方體網格，結果comsol吃60GB記憶體
05/08 22:05, 45F
→
05/08 22:05, , 46F
ansys只吃20幾GB的記憶體，差蠻多了
05/08 22:05, 46F
→
05/08 22:06, , 47F
只是comsol與ansys的前六個特徵頻率不一樣，以及之後的
05/08 22:06, 47F
→
05/08 22:06, , 48F
特徵頻率都有一些誤差，這樣而已。
05/08 22:06, 48F
→
05/08 22:07, , 49F
意思是說可能有兩種: 1.ansys吃比較少記憶體
05/08 22:07, 49F
→
05/08 22:07, , 50F
2.ansys邊界條件設定錯誤?
05/08 22:07, 50F
→
05/08 22:07, , 51F
還有一種可能就是comsol邊界條件設定錯誤...
05/08 22:07, 51F
推
05/08 22:25, , 52F
ansys和comsol管理記憶體的方式不同，使用比較少的記憶體是
05/08 22:25, 52F
→
05/08 22:26, , 53F
有可能的，這個沒辦法，程式就是這樣。
05/08 22:26, 53F
→
05/08 22:26, , 54F
至於前六個特徵頻率，ansys內如果沒任何BC，那前六個應該是
05/08 22:26, 54F
→
05/08 22:27, , 55F
0，或是接近0的數字，分別代表6個自由度的剛體運動。
05/08 22:27, 55F
→
05/08 22:27, , 56F
而comsol可能有很弱的拘束，所以計算出來的值會略大於0。
05/08 22:27, 56F
→
05/08 22:28, , 57F
這個沒什麼對或錯，單純看你要的是什麼，一般的做法就是前6
05/08 22:28, 57F
→
05/08 22:28, , 58F
個頻率不去看，就這樣。
05/08 22:28, 58F
‣ 更多 candy88257 的文章
文章代碼(AID): #1JPnG0Hb (Cad_Cae)
Cad_Cae 近期熱門文章
更多 近期熱門文章 >>
PTT數位生活區 即時熱門文章
更多 即時熱門文章 >>
‣ 更多 candy88257 的文章
文章代碼(AID): #1JPnG0Hb (Cad_Cae)